Search NASA⌕ Search

SEARCH · Search NASA

Results for “covalent functionalization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Permanent Magnets Featuring Heavy Main Group Elements for Magnetic Anisotropy

Permanent magnets are the functional component of electric motors and generators found in numerous renewable energy applications. To improve energy conversion in such applications, we require fundamentally new magnets that generate higher magnetic flux per volume while retaining the properties conferred by rare-earth elements incorporated into current technologies. We hypothesize that by engendering a covalent interaction between two elements, we can access a new regime of magnetic materials where the two components of a magnetic moment—spin and orbital angular momentum—come from two separate atoms to form a complete magnetic moment. Our previous research utilized high-pressure conditions to discover two new candidate materials ideal for assessing this hypothesis. The first material, FeBi 2 , enables the study of an unprecedented solid state metal-metal bonding interaction. The second, MnBi 2 represents the second member of the promising Mn–Bi family known for its magnetic properties. Importantly, MnBi 2 is isostructural to FeBi 2 . Together these chemically simple but magnetically rich materials provide an elegant platform for elucidating fundamental design principals of magnetic anisotropy while inspiring the synthesis of new magnetic materials. More generally, solid-state chemistry remains a synthetic black box. To fully harness the potential of such new materials, it is vital to create a window into both the structure and the properties of our high-pressure materials. One missing area of importance to energy science is high-pressure magnetometry. To this end, we assessed the magnetic structure of these materials through two objectives: (1) we performed magnetometry measurements at high-pressures, including performing high-pressure X-ray Magnetic Circular Dichroism (XMCD) experiments on these materials that are only synthesized at high-pressure and cannot be recovered to ambient conditions and (2) we targeted the recovery of these materials using dynamic compression approaches, starting with laser shock techniques. Future directions include expanding our focus to ternary phase space, employing high-throughput computational approaches to direct our search for magnetic materials across the Periodic Table. This award was originated at Northwestern, and reopened as a new grant at MIT. This is the close out for Northwestern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Selective Covalent Basal Plane Modification as a Probe of Hydroxide Ion Conduction Pathways in Magnesium Aluminum Layered Double Hydroxides

Understanding ionic conduction in layered double hydroxides (LDHs) is a crucial step towards utilizing them as solid, hydroxide ion‐conducting electrolytes in energy conversion applications. We selectively modified the interlayer and external surfaces of MgAl LDHs with tris(hydroxymethyl)aminomethane (TRIS) ligands. By adjusting the concentration of the TRIS surface modifier, the LDH basal plane surfaces could be functionalized everywhere (internally and externally) or only externally. External modification resulted in loss of OH‐conductivity compared to pristine LDHs, confirming that external platelet surfaces are the primary ion conduction pathway.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Vacancy and interstitial interactions with crystal/amorphous, metal/covalent interfaces

In this work, we use atomistic simulations to investigate the interaction of vacancies and interstitials with interfaces between a crystalline metal and an amorphous, covalently-bonded solid. We select the gold (Au)/silicon (Si) binary system as a model material and construct interface models along two different facets of crystalline Au and with amorphous Si (a-Si) created at three different quench rates. We compute formation energies of vacancies, self-interstitials, and interstitial impurities as a function of position relative to the interface and find that they have markedly lower values near the interface than in the interior of the adjoining phases. We conclude that crystal/amorphous, metal/covalent interfaces may be as effective at removing radiation-induced point defects as interfaces in polycrystalline metals composites. Moreover, irrespective of interface character, the average formation energies of all point defects at all the Au/a-Si interfaces we investigated are comparable. Thus, unlike in polycrystalline metals, where an interface’s crystallographic character has a marked effect on its interactions with point defects, all interface types in crystal/amorphous, metal/covalent composites may be equally effective at absorbing all radiation-induced defects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topochemical Deintercalation of Li from Layered LiNiB: toward 2D MBene

The pursuit of two-dimensional (2D) borides, MBenes, has proven to be challenging, not the least because of the lack of a suitable precursor prone to the deintercalation. Here, we studied room-temperature topochemical deintercalation of lithium from the layered polymorphs of the LiNiB compound with a considerable amount of Li stored in between [NiB] layers (33 at. % Li). Deintercalation of Li leads to novel metastable borides (Li~ 0.5 NiB) with unique crystal structures. Partial removal of Li is accomplished by exposing the parent phases to air, water, or dilute HCl under ambient conditions. Scanning transmission electron microscopy and solid-state 7 Li and 11 B NMR spectroscopy, combined with X-ray pair distribution function (PDF) analysis and DFT calculations, were utilized to elucidate the novel structures of Li~ 0.5 NiB and the mechanism of Li-deintercalation. We have shown that the deintercalation of Li proceeds via a “zip-lock” mechanism, leading to the condensation of single [NiB] layers into double or triple layers bound via covalent bonds, resulting in structural fragments with Li[NiB] 2 and Li[NiB] 3 compositions. The crystal structure of Li~ 0.5 NiB is best described as an intergrowth of the ordered single [NiB], double [NiB] 2 , or triple [NiB] 3 layers alternating with single Li layers; this explains its structural complexity. Here, the formation of double or triple [NiB] layers induces a change in the magnetic behavior from temperatureindependent paramagnets in the parent LiNiB compounds to the spin-glassiness in the deintercalated Li~ 0.5 NiB counterparts. LiNiB compounds showcase the potential to access a plethora of unique materials, including 2D MBenes (NiB).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Crystal structures of an E1–E2–ubiquitin thioester mimetic reveal molecular mechanisms of transthioesterification

E1 enzymes function as gatekeepers of ubiquitin (Ub) signaling by catalyzing activation and transfer of Ub to tens of cognate E2 conjugating enzymes in a process called E1–E2 transthioesterification. The molecular mechanisms of transthioesterification and the overall architecture of the E1–E2–Ub complex during catalysis are unknown. Here, we determine the structure of a covalently trapped E1–E2–ubiquitin thioester mimetic. Two distinct architectures of the complex are observed, one in which the Ub thioester (Ub(t)) contacts E1 in an open conformation and another in which Ub(t) instead contacts E2 in a drastically different, closed conformation. Altogether our structural and biochemical data suggest that these two conformational states represent snapshots of the E1–E2–Ub complex pre- and post-thioester transfer, and are consistent with a model in which catalysis is enhanced by a Ub(t)-mediated affinity switch that drives the reaction forward by promoting productive complex formation or product release depending on the conformational state.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Thermodynamic non-ideality and disorder heterogeneity in actinide silicate solid solutions

Non-ideal thermodynamics of solid solutions can greatly impact materials degradation behavior. We have investigated an actinide silicate solid solution system (USiO 4 –ThSiO 4 ), demonstrating that thermodynamic non-ideality follows a distinctive, atomic-scale disordering process, which is usually considered as a random distribution. Neutron total scattering implemented by pair distribution function analysis confirmed a random distribution model for U and Th in first three coordination shells; however, a machine-learning algorithm suggested heterogeneous U and Th clusters at nanoscale (~2 nm). The local disorder and nanosized heterogeneous is an example of the non-ideality of mixing that has an electronic origin. Partial covalency from the U/Th 5f–O 2p hybridization promotes electron transfer during mixing and leads to local polyhedral distortions. The electronic origin accounts for the strong non-ideality in thermodynamic parameters that extends the stability field of the actinide silicates in nature and under typical nuclear waste repository conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Molecule-surface interaction from van der Waals-corrected semilocal density functionals: The example of thiophene on transition-metal surfaces

Semilocal density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semilocal functionals require an additive nonlocal correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of aromatic thiophene (C 4 H 4 S) bound to transition-metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well-balanced short- and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy. Here, our methods indicate that the correct interplay is not present in any combination of recent meta-GGA's and rVV10. The simple short-range damping of the vdW correction scheme that is practically successful in the combination of GGA's and vdW approximations is less transferable in SCAN+rVV10 or in the revised version, revSCAN+rVV10. In addition, we present accurate random-phase-approximation-quality adsorption energies from a model based on the one of Zaremba and Kohn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Manipulating the Geometric and Electronic Structures of Manganese Nitrides for Ammonia Synthesis

Abstract Manganese (Mn) nitrides are important nitrogen (N) carriers for small‐scale intermittent ammonia (NH 3 ) synthesis. However, only 3∼8 % of lattice N are converted into NH 3 . In this study, the geometric and electronic structures of well‐defined Mn 4 N and Mn 2 N lattices were altered using transition metal heteroatoms (Cr, Fe, Co, Ni, Mo) to understand the driving force behind lattice N diffusion and extraction. Density Functional Theory (DFT) revealed that the binding of early hydrogenation product (NH) follows a linear relationship with lattice N over a wide range of close‐packed surfaces. But, the binding of NH 2 and NH 3 are more sensitive to the geometric and electronic structures. Further, the chemical bonding can be quantitatively characterized with the covalency derived from Crystal Orbital Hamiltonian Population (COHP). In the Eley Rideal‐Mars van Krevelan pathway, the overall NH 3 formation free energy ( ) and lattice N diffusion barrier ( ) are the respective thermodynamic and kinetic determining factors. Aided by a rate‐determining step (RDS) model, Mn 4 N modified with Fe, Co, Ni single‐atom dopants all show enhanced rates for NH 3 formation.

Shan, Nannan↗
DOE OSTI2022

Surface Immobilization of a Re(I) Tricarbonyl Phenanthroline Complex to Si(111) through Sonochemical Hydrosilylation

A sonochemical-based hydrosilylation method was employed to covalently attach a rhenium tricarbonyl phenanthroline complex to silicon(111). fac-Re(5-(p-Styrene)-phen)(CO) 3 Cl (5-(p-styrene)-phen = 5-(4-vinylphenyl)-1,10-phenanthroline) was reacted with hydrogen-terminated silicon(111) in an ultrasonic bath to generate a hybrid photoelectrode. Subsequent reaction with 1-hexene enabled functionalization of remaining atop Si sites. Attenuated total reflectance–Fourier transform infrared spectroscopy confirms attachment of the organometallic complex to silicon without degradation of the organometallic core, supporting hydrosilylation as a strategy for installing coordination complexes that retain their molecular integrity. Detection of Re(I) and nitrogen by X-ray photoelectron spectroscopy (XPS) further support immobilization of fac-Re(5-(p-styrene)-phen)(CO) 3 Cl. Cyclic voltammetry and electrochemical impedance spectroscopy under white light illumination indicate that fac-Re(5-(p-styrene)-phen)(CO) 3 Cl undergoes two electron reductions. Mott–Schottky analysis indicates that the flat band potential is 239 mV more positive for p-Si(111) co-functionalized with both fac-Re(5-(p-styrene)-phen)(CO) 3 Cl and 1-hexene than when functionalized with 1-hexene alone. XPS, ultraviolet photoelectron spectroscopy, and Mott–Schottky analysis show that functionalization with fac-Re(5-(p-styrene)-phen)(CO) 3 Cl and 1-hexene introduces a negative interfacial dipole, facilitating reductive photoelectrochemistry.

X-ray photoelectron spectroscopy↗
DOE OSTI2025

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗
DOE OSTI2023

Improved Fuel Cell Chemical Durability of an Heteropoly Acid Functionalized Perfluorinated Terpolymer-Perfluorosulfonic Acid Composite Membrane

Commercial proton exchange membrane heavy-duty fuel cell vehicles will require a five-fold increase in durability compared to current state-of-the art light-duty fuel cell vehicles. We describe a new composite membrane that incorporates silicotungstic heteroply acid (HPA), α -K 8 SiW 11 O 40 ▪13H 2 O, a radical decomposition catalyst and when acid-exchanged can potentially conduct protons. The HPA was covalently bound to a terpolymer of tetrafluoroethylene, vinylidene fluoride, and sulfonyl fluoride containing monomer (1,1,2,2,3,3,4,4-octafluoro-4-((1,2,2-trifluorovinyl)oxy)butane-1-sulfonyl fluoride) by dehydrofluorination followed by addition of diethyl (4-hydroxyphenyl) phosphonate, giving a perfluorosulfonic acid-vinylidene fluoride-heteropoly acid (PFSA-VDF-HPA). A composite membrane was fabricated using a blend of the PFSA-VDF-HPA and the 800EW 3M perfluoro sulfonic acid polymer. The bottom liner-side of the membrane tended to have a higher proportion of HPA moieties compared to the air-side as gravity caused the higher mass density PFSA-VDF-HPA to settle. The composite membrane was shown to have less swelling, more hydrophobic properties, and higher crystallinity than the pure PFSA membrane. The proton conductivity of the membrane was 0.130 ± 0.03 S cm −1 at 80 °C and 95% RH. Impressively, when the membrane with HPA-rich side was facing the anode, the membrane survived more than 800 h under accelerated stress test conditions of open-circuit voltage, 90 °C and 30% RH.

08 HYDROGEN↗
DOE OSTI2025

Critical Role of the Steric Factor in the Viscoelasticity of Vitrimers

Dynamic covalent networks (DCNs) are a promising solution to mitigate plastic-waste-related issues through improved recyclability enabled by dynamic bonds. However, our understanding of the mechanisms controlling their viscoelasticity, especially in vitrimers where bond exchange relies on associative reactions, remains limited. Here, in this study, we investigate the dynamics in model DCNs with boric ester functionalities, and the analysis of the temperature dependence of their terminal relaxation times revealed a puzzling result: extremely large Arrhenius prefactors associated with their dynamic bond rearrangement times. We ascribe this observation to the often-overlooked chemical steric factor that slows down chemical reactions, therefore decreasing the vitrimers’ bond exchange rate by many orders. The estimated steric factor of the bond exchange in our model DCNs is comparable to those observed in boronic ester exchange reactions between small molecules. Additional analysis of literature data revealed an overall low steric factor also for imine bond exchange, thus highlighting the role of this parameter in tremendously slowing down bond exchange in DCNs despite low activation energy barriers. We propose a general approach for designing vitrimers with desired viscoelastic and creep properties considering the critical role of the steric factor in bond rearrangement mechanisms, in addition to the traditionally considered activation energy and matrix properties.

Carden, Gregory P. [Univ. of Tennessee, Knoxville,↗
DOE OSTI2024

Can the Rate of a Catalytic Turnover Be Altered by Ligands in the Absence of Direct Binding Interactions?

Second sphere coordination effects ubiquitous in enzymatic catalysis occur through direct interactions, either covalent or non-covalent, with reaction intermediates and transition states. Here, we present herein evidence of indirect second sphere coordination effects in which ligation of water/alkanols far removed from the primary coordination sphere of the active site nevertheless alter energetic landscapes within catalytic redox cycles in the absence of direct physicochemical interactions with surface species mediating catalytic turnovers. Density functional theory, in situ X-ray absorption and infrared spectroscopy, and a wide array of steady-state and transient CO oxidation rate data suggest that the presence of peripheral water renders oxidation half-cycles within two-electron redox cycles over μ m -oxo-bridged trimers in MIL-100(M) more kinetically demanding. Communication between ligated water and the active site appears to occur through the Fe-O-Fe backbone, as inferred from spin density variations on the central μ m -oxygen 'junction'. Evidence is provided for the generality of these second sphere effects in that they influence different types of redox half-cycles or metals, and can be amplified or attenuated through choice of coordinating ligand. Specifically in the case of MIL-100(M) materials, the Cr isostructure can be made to kinetically mimic the Fe variant by disproportionately hindering oxidation half-cycles relative to the reduction half-cycles. Kinetic and spectroscopic inferences presented here significantly expand both the conceptual definition of second sphere effects as well as the palette of synthetic levers available for tuning catalytic redox performance through chemical ligation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Effect of Methyl Groups on Formation of Ordered or Layered Graphitic Materials from Aromatic Molecules

Abstract Developing functionally complex carbon materials from small aromatic molecules requires an understanding of how the chemistry and structure of its constituent molecules evolve and crosslink, to achieve a tailorable set of functional properties. Here, molecular dynamics (MD) simulations are used to isolate the effect of methyl groups on condensation reactions during the oxidative process and evaluate the impact on elastic modulus by considering three monodisperse pyrene‐based systems with increasing methyl group fraction. A parameter to quantify the reaction progression is designed by computing the number of new covalent bonds formed. Utilizing the previously developed MD framework, it is found that increasing methylation leads to an almost doubling of bond formation, a larger fraction of the new bonds oriented in the direction of tensile stress, and a higher basal plane alignment of the precursor molecules along the direction of tensile stress, resulting in enhanced tensile modulus. Additionally, via experiments, it is demonstrated that precursors with a higher fraction of methyl groups result in a higher alignment of molecules. Moreover, increased methylation results in the lower spread of single molecule alignment which may lead to smaller variations in tensile modulus and more consistent properties in carbon materials derived from methyl‐rich precursors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Electron Counting and High-Pressure Phase Transformations in Metal Hexaborides

Pressure-induced structural transitions of the alkaline earth hexaborides, CaB 6 , SrB 6 , and BaB 6 , are studied theoretically using electron counting rules and density functional theory calculations. Here, we demonstrate the applicability of gas-phase borane electron counting methods to solid-state metal borides under pressure and validate the assumptions of the rules by density functional theory (DFT) calculations. All three compounds share ambient-pressure and high-pressure structures, but BaB 6 differs from CaB 6 and SrB 6 at intermediate pressures. The unique BaB 6 phase is shown to break electron counting rules, while all other phases obey them. This anomaly is resolved by DFT, which reveals B–Ba covalency and unusual B–B π bonding under pressure. The relationships between structure and bonding can help us to understand the exotic behavior of lanthanide hexaborides and design new borides with desirable properties. Developing electron counting procedures for solids will enhance materials discovery efforts with chemical intuition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ion Pairing Mediates Molecular Organization Across Liquid/Liquid Interfaces

Liquid/liquid interfaces play a central role in scientific fields ranging from nanomaterial synthesis and soft matter electronics to nuclear waste remediation and chemical separations. This diversity of functions arises from an interface’s ability to respond to changing conditions in its neighboring bulk phases. Understanding what drives this interfacial flexibility can provide novel avenues for designing new functional interfaces. However, limiting this progress is an inadequate understanding of the subtle intermolecular and interphase interactions taking place at the molecular level. Here, we use surface-specific vibrational sum frequency generation spectroscopy combined with atomistic molecular dynamics simulations to investigate the self-assembly and structure of model ionic oligomers consisting of an oligodimethylsiloxane (ODMS) tail covalently attached to a positively charged methyl imidazolium (MIM + ) head group at buried oil/aqueous interfaces. We show how the presence of seemingly innocuous salts can impart dramatic changes to the ODMS tail conformations in the oil phase via specific ion effects and ion-pairing interactions taking place in the aqueous phase. These specific ion interactions are shown to drive enhanced amphiphile adsorption, induce morphological changes, and disrupt emergent hydrogen-bonding structures at the interface. Tuning these interactions allows for independent control over the oligomer structure in the oil phase versus interfacial population changes and represents key mechanistic insight that is needed to control chemical reactions at liquid/liquid interfaces.

36 MATERIALS SCIENCE↗