Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular energetics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP 6 ) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP 6 and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Site-specific electronic structure of covalently linked bimetallic dyads from nitrogen K-edge x-ray absorption spectroscopy

A nitrogen K-edge x-ray absorption near-edge structure (XANES) survey is presented for tetrapyrido[3,2-a:',3'-c:3'',2''-h:2''',3'''-j]phenazine (tpphz)-bridged bimetallic assemblies that couple chromophore and catalyst transition metal complexes for light driven catalysis, as well as their individual molecular constituents. We demonstrate the high N site sensitivity of the N pre-edge XANES features, which are energetically well-separated for the phenazine bridge N atoms and for the individual metal-bound N atoms of the inner coordination sphere ligands. By comparison with the time-dependent density functional theory calculated spectra, we determine the origins of these distinguishable spectral features. We find that metal coordination generates large shifts toward higher energy for the metal-bound N atoms, with increasing shift for 3d < 4d < 5d metal bonding. This is attributed to increasing ligand-to-metal σ donation that increases the effective charge of the bound N atoms and stabilizes the N 1s core electrons. In contrast, the phenazine bridge N pre-edge peak is found at a lower energy due to stabilization of the low energy electron accepting orbital localized on the phenazine motif. While no sensitivity to ground state electronic coupling between the individual molecular subunits was observed, the spectra are sensitive to structural distortions of the tpphz bridge. Here these results demonstrate N K-edge XANES as a local probe of electronic structure in large bridging ligand motifs, able to distinctly investigate the ligand-centered orbitals involved in metal-to-ligand and ligand-to-ligand electron transfer following light absorption.

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

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State

The COVID-19 pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the “D-form” to the “G-form”) that carried an amino acid substitution D614G in its “Spike” protein. The G-form is more infectious in vitro and associated with increased viral loads. To gain insight into the molecular-level underpinnings of these characteristics, we employed microsecond all-atom simulations. Here we show that changes in the protein energetics favor a higher population of infection-capable (open) states through release of hydrogen bonds of an asymmetry present in the D-form but not the G-form. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are significant for vaccine design. The Molecular Dynamics Simulations datasets generated in this study, namely that of the soluble form SARS-CoV-2 Spike protein are made available here. Details of the trajectory and supporting files are given in the README.txt file.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether “staples” into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with pchloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔT m ' = +27 °C) and temperature-induced heme loss (ΔT 50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

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

Experimental and Theoretical Studies of Ultrafast Vibrational Energy Transfer Dynamics in Energetic Materials

Energy transfer through anharmonically-coupled vibrations influences the earliest chemical steps in shockwave-induced detonation in energetic materials. A mechanistic description of vibrational energy transfer is therefore necessary to develop predictive models of energetic material behavior. We performed transient broadband infrared spectroscopy on hundreds of femtoseconds to hundreds of picosecond timescales as well as density functional theory and molecular dynamics simulations to investigate the evolution of vibrational energy distribution in thin film samples of pentaerythritol tetranitrate (PETN) , 1,3,5 - trinitroperhydro - 1,3,5 - triazine (RDX) , and 2,4,6 - triamino 1,3,5 - trinitrobenzene (TATB). Experimental results show dynamics on multiple timescales, providing strong evidence for coupled vibrations in these systems, as well as material-dependent evolution on tens to hundreds of picosecond timescales. Theoretical results also reveal pathways and distinct timescales for energy transfer through coupled vibrations in the three investigated materials, providing further insight into the mechanistic underpinnings of energy transfer dynamics in energetic material sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

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

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗
DOE OSTI2021

From inert gas to fertilizer, fuel and fine chemicals: N 2 reduction and fixation

The 100th anniversary of a leading nitrogen fixation technology developer like CASALE SA is a reason to reflect over the 20th century successful solution of the problem of world food supply, and to look out for solutions for sustainable developments with respect to ammonia production. Here, we review the role of nitrogen as essential chemical constituent in photosynthesis and biology, and component of ammonia as it is used as fertilizer for primary production by photosynthesis for farming and food supply and its future role as energy carrier. While novel synthesis methods and very advanced synchrotron based x-ray analytical techniques are being developed, we feel it is important to refer to the historical and economical context of nitrogen. The breaking of the N≡N triple bond remains a fundamental chemical and energetic problem in this context. We review the electrochemical ammonia synthesis as an energetically and environmentally benign method. Two relatively novel X-ray spectroscopy methods, which are relevant for the molecular understanding of the catalysts and biocatalysts, i.e. soft X-ray absorption spectroscopy and nuclear resonant vibration spectroscopy are presented. We illustrate the perceived reality in fertilizer usage on the field, and fertilizer production in the factory complex with photos and thus provide a contrast to the academic view of the molecular process of ammonia function and production.

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

Engineering Electrical Conductivity in Stable Zirconium-Based PCN-222 MOFs with Permanent Mesoporosity

Electrically conductive metal–organic frameworks (MOFs) featuring zirconium-based nodes are of great interest for electrochemical and optoelectronic applications due to their exceptional thermal and chemical stability, although the number of such MOFs remains limited. Here in this paper, we demonstrate that electron deficient molecules such as transition metal(IV) bis(dicarbollide) (TM IV (C 2 B 9 H 11 ) 2 ; TM IV CB for short, where TM = Ni, Pd, and Pt) or C 60 , spatially infiltrated in the microporous channels of the zirconium-based porphyrinic MOFs PCN-222 and PCN-222-Zn, generate highly stable and electrically conductive frameworks in which the mesopores remain accessible to other guests. Solid-state density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations indicate incorporation of Ni IV CB and C 60 in the MOF micropores is energetically and structurally feasible. Interestingly, in contrast to the free species, strong donor–acceptor interactions between the MOF and Ni IV CB restrain it to a single conformation. Calculated electronic structures and charge-hopping conduction probabilities illustrate that efficient charge-transfer (CT) from photoexcited linkers to the guest molecules facilitates charge hopping in the framework, making the MOFs electrically conductive. The donor–acceptor conjugates also enhance exciton dissociation at their heterojunctions, fostering the formation of long-lived electron-trapped states with potential utility for photo- and electrochemical devices.

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

Toward Understanding and Controlling Organic Reactions on Metal Oxide Catalysts

Metal oxides have structurally complex surfaces on which a variety of adsorption site types can occur, including cation sites, anion sites, oxygen vacancy sites, and Brønsted acid sites. These sites can catalyze the catalytic transformation of organic molecules via diverse routes, thus enabling H abstraction, O abstraction, C–C bond formation, and other reactions. This Perspective provides an update on recent advances and future directions for various organic reactions on metal oxide catalyst surfaces, particularly for C–H activation of alkanes and for C–C bond formation with organic oxygenate reactants. Here, we put emphasis on the molecular scale details, on the active site structures required to enable the formation of kinetically relevant transition states, energetic descriptors, as well as contemporary ideas to enable low activation energies. This progress has been enabled by specialized experiments and the increased capabilities of modern electronic structure calculations.

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

Microscopic and Macroscopic Characterization of Grain Boundary Energy and Strength in Silicon Carbide via Machine-Learning Techniques

Predicting the properties of grain boundaries poses a challenge because of the complex relationships between structural and chemical attributes both at the atomic and continuum scales. Grain boundary systems are typically characterized by parameters used to classify local atomic arrangements in order to extract features such as grain boundary energy or grain boundary strength. The present work utilizes a combination of high-throughput atomistic simulations, macroscopic and microscopic descriptors, and machine-learning techniques to characterize the energy and strength of silicon carbide grain boundaries. Additionally, a diverse data set of symmetric tilt and twist grain boundaries are described using macroscopic metrics such as misorientation, the alignment of critical low-index planes, and the Schmid factor, but also in terms of microscopic metrics, by quantifying the local atomic structure and chemistry at the interface. These descriptors are used to create random-forest regression models, allowing for their relative importance to the grain boundary energy and decohesion stress to be better understood. Results show that while the energetics of the grain boundary were best described using the microscopic descriptors, the ability of the macroscopic descriptors to reasonably predict grain boundaries with low energy suggests a link between the crystallographic orientation and the resultant atomic structure that forms at the grain boundary within this regime. For grain boundary strength, neither microscopic nor macroscopic descriptors were able to fully capture the response individually. However, when both descriptor sets were utilized, the decohesion stress of the grain boundary could be accurately predicted. These results highlight the importance of considering both macroscopic and microscopic factors when constructing constitutive models for grain boundary systems, which has significant implications for both understanding the fundamental mechanisms at work and the ability to bridge length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel

In this LDRD “Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel” a peer-reviewed lab-scale hypersonic wind tunnel (WT) was planned to be assembled in-house, along with time and space resolved and related flow and materials diagnostics to derive high quality measurements and the science that is able to inform the development of our modeling and simulation codes for the hypersonic regime. The baseline capability to achieve Mach 3 flow was scheduled to be completed in Q4 2020, in the first year of the LDRD. However due to the COVID-19 delays and travel restrictions starting in March, 2020 this task could not be completed to carry out the materials response studies in-house, or those planned with external partners that were engaged during this LDRD (NASA Langley, NASA Ames, Texas A&M university, and Virginia Tech). Still, using our hypersonic computational fluid simulation codes, many of the hypersonic components have been designed (for Mach 3, 5, and 7) and procured (or in procurement phase), and a new dedicated hypersonic laboratory was facilitized, 2 WT test sections have been manufactured, and a Mach 3 nozzle was 3D printed for prototyping. In parallel, unique materials diagnostics have been tested under static conditions and fully specified and awaiting validation in our partner facilities. Due to the new hypersonic initiative announced in January, 2020 we proposed a related scope increase in hypersonics research requiring energy interactions and new diagnostics well beyond what could be supported by an LDRD ER. A decision was made to end the current LDRD ER and formally propose a new LDRD SI that leverages the current effort, which will be starting in FY21. As a result of the uniqueness of this LDRD work, our early efforts have been encouraged and well received in the hypersonic research community, and is now part of multiple unsolicited internal and external project proposals expected to be funded in the coming years, which is encouraging to make LLNL a leader in experimentally probing extreme physics of local energetics critical to hypersonics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

In situ molecular imaging of ion clusters reveals the acid gas capture capacity and mechanism of water-lean ionic liquids

Water-lean solvents are a promising technology for capturing acid gases like carbon dioxide (CO 2 ). In situ liquid time-of-flight secondary ionization mass spectroscopy (ToF-SIMS) is used to study a representative solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (2-EEMPA) with different CO 2 loadings to reveal the complex solvent structure upon CO 2 capture. Characteristic peaks of 2-EEMPA, such as m/z – 215 C 11 H 23 N 2 O 2 – (deprotonated 2-EEMPA) and m/z + 217 C 11 H 25 N 2 O 2 + (protonated 2-EEMPA), are detected due to acid gas uptake. Also, solvent molecules and carboxylate ion pairs, such as m/z – 259 C 12 H 23 N 2 O 4 – [(deprotonated 2-EEMPA∙∙∙CO 2 )] and m/z + 261 C 12 H 25 N 2 O 4 + (protonated 2-EEMPA∙∙∙CO 2 ), are observed. Interestingly, more than one CO 2 molecule can be captured per each solvent molecule as evidenced in SIMS mass spectra, for example, m/z – 321 C 13 H 25 N 2 O 7 – [(deprotonated 2-EEMPA)∙∙∙2CO 2 ∙∙∙H 2 O], m/z + 305 C 13 H 25 N 2 O 4 + [(protonated 2-EEMPA)∙∙∙2CO 2 ], m/z – 389 C 17 H 29 N 2 O 8 – [(deprotonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙3CH 2 ], and m/z + 373 C 16 H 25 N 2 O 8 + [(protonated 2-EEMPA)∙∙∙3CO 2 ∙∙∙2C]. However, the monomer of 2-EEMPA and CO 2 seems to be most prevalent. Furthermore, solvent clusters are detected in loaded solvents, for instance m/z + 433 C 22 H 49 N 4 O 4 + [(2-EEMPA)2∙∙∙H] and m/z + 646 [(2-EEMPA) 3 -2H], while capturing CO 2 at different amounts. Relative abundance of cluster ions provides a semi-qualitative venue to assess the free energies of gas capture energetics, indicating the relative stability trend within the same solvent system, previously impossible. These observed ion clusters are verified with molecular modeling, where dimer, trimer, and cluster ions are validated for their presence either due to weak molecular interactions or hydrogen bonds. In situ molecular imaging of ionic liquids and molecular modeling reveals that the acid gas capture mechanism by ionic liquids includes both physical adsorption and chemical bonding with multiple reaction pathways, engaging cluster formation and alteration of solvent structures.

Acid gas capture↗
DOE OSTI2024

Ni and Co Incorporation in Forsterite: A Density Functional Theory Study with Hubbard Correction

Ni and Co are critical elements for the world economy and modern technologies. Mafic and ultramafic deposits represent low-grade yet abundant alternatives to traditional Ni and Co ores. Here, in this work, density functional theory (DFT) with the Hubbard U correction (DFT+U) was used to simulate the incorporation of Ni and Co in forsterite (Mg 2 SiO 4 ), the Mg endmember of olivine, a common mineral in mafic and ultramafic rocks. Hubbard U terms for Ni and Co were parametrized using a series of oxide, hydroxide, carbonate, silicate, and sulfide minerals relevant to extraction and recovery of Ni and Co from mafic and ultramafic deposits. Electronic, energetic, magnetic, and structural properties were considered in the parametrization. For each of Ni and Co, an effective Hubbard correction (U eff ) value that optimized agreement with either experimental data or a hybrid exchange-correlation functional for all of the minerals considered is reported. DFT+U ab initio molecular dynamics (AIMD) simulations of Ni and Co incorporated into the M1 and M2 octahedral sites of forsterite were then performed. Ni and Co substitution in the M1 site was more energetically favorable than substitution in the M2 site, in agreement with published partition coefficients. AIMD trajectories were used to compute extended X-ray absorption fine structure (EXAFS) spectra of Ni in the M1 and M2 sites for direct fitting to a published experimental spectrum of Ni in a natural San Carlos olivine sample. The results of the fit indicated that ordering of Ni in the M1 site was not as strong at the low Ni concentrations relevant to mafic and ultramafic silicate minerals as that at the higher concentrations of the Ni-Mg olivine solid solutions studied to date.

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

First-Principles Study of Molecular Hydrogen Activation by Defects in Boron Nitride

Here, we used density functional theory simulations in combination with ab initio thermodynamics to determine the H 2 partial pressure (p H2 )-dependent energetics associated with H 2 activation and recovery at various defect sites in hexagonal boron nitride (h-BN). We found that some defects are very reactive with hydrogen, thereby definitely trapping hydrogen in defective h-BN. However, depending on hydrogen partial pressure, less reactive defect sites can be populated. Because of the lower binding capability of these sites, they would allow hydrogen to be recycled and recovered. For small defect sizes, we found that hydrogen preferentially binds to nitrogen sites by forming N–H bonds, and if no N sites are available then boron sites would be the next to bind hydrogen. Hydrogen dissociation via frustrated Lewis pair is found to be more favorable than forming only N–H bonds but only if the defect size is large enough to accommodate steric effects. For specific conditions such as T = 400 K, p H2 = 1 bar, and only considering one molecular H 2 per defect, three defects, namely, the N monovacancy, 3V(1B2N), and hexagonal 6V(3B3N) could play a role in both the activation and recycling of H 2 as they would be reacting enough to allow a favorable splitting of H 2 while not binding too strongly to allow its recovery. More broadly, a range of p H2 and hydrogen loading conditions were investigated for different types of defects and the finding suggests that p H2 could be used to fine-tune the Gibbs free energy of hydrogenation, thereby allowing several types of defects at different hydrogen loading contents to play a role in the activation/recovery process of H 2 in defective h-BN.

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

Infrared bands of neutral gas-phase carbon clusters in a broad spectral range

The identification of species in the interstellar medium requires precise and molecule-specific spectroscopic information in the laboratory framework, in broad spectral ranges and under conditions relevant to interstellar environments. In this work, we measure the gas-phase infrared spectra of neutral carbon clusters, C N (N = 6–11), in a molecular beam. The C N distribution is formed by photofragmentation of C 60 molecules, concurrently showing a top-down formation mechanism. A broad spectral range in the infrared between 500–3200 cm -1 (20–3.125 μm) is investigated. We observe strong bands between 5 and 6 μm, in conjunction with novel features in the 3 μm region. Density functional theory calculations reveal that these short wavelength modes correspond to combination bands with significant infrared intensity. Moreover, we identify the N ≤ 10 clusters as linear, while C 11 adopts a ring configuration, placing the linear-to-ring transition at N = 11 under our molecular beam conditions. The linearity of C 10 is discussed based on the formation pathway from larger clusters in energetic conditions. Given the vast and very precise infrared information already been released from the James Webb Space Telescope mission, this infrared spectroscopic data set in conjunction with information on formation mechanisms is of major relevance for identifying neutral carbon clusters in astronomical environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗