Search NASA⌕ Search

SEARCH · Search NASA

Results for “anharmonic”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 1: Path-integral calculations with CCSD(T) quality potentials

Isotopic compositions of alkanes are typically assumed to be kinetically controlled, but recently is has been proposed that alkanes can isotopically equilibrate for both C and H isotopes during natural gas generation. Evaluation of this requires knowledge of the isotopic equilibrium between alkanes and other common hydrogen and carbon bearing species. Here, in this study, we calculate isotopic equilibria within and between gaseous dihydrogen (H 2 ), water (H 2 O), methane (CH 4 ), ethane (C 2 H 6 ) and propane (C 3 H 8 ), including isotope fractionation among molecules, clumped isotope effects, as well as among sites of propane (i.e., the site-specific isotope effects) from 0°C to 500°C using a path-integral method paired with high-level descriptions of molecular potentials and the diagonal correction to the Born Oppenheimer approximation. While path-integral calculations with high- level CCSD(T) potentials are available for the isotopic equilibria involving methane, the path-integral calculations for ethane and propane have only been performed based on lower-level descriptions of the molecular potentials. We analyze the relative importance of various approximations that are commonly employed when isotopic equilibria are evaluated. We find that clumped isotope effects can be calculated to the same accuracy using computationally inexpensive combination of the Bigeleisen-Mayer-Urey model with the molecular potential from density functional theory. In contrast, fractionation and site preferences of both deuterium and carbon-13 benefit from the use of the higher level CCSD(T) potentials and accounting for anharmonic effects. Additionally, for fractionation and site preference of deuterium corrections to Born-Oppenheimer approximation can also be important.

03 NATURAL GAS↗
DOE OSTI2024

Structure and disorder resulting in ultralow thermal conductivity in the defect chalcopyrite AgInSnSe 4

Quaternary chalcogenides are attractive for a variety of technological fields of interest due to the diverse physical properties these materials possess. Herein, the structural and thermal properties of the disordered chalcopyrite AgInSnSe 4 are reported revealing a distinct relation between its structure and physical properties. Modeling of the experimental temperature-dependent thermal properties revealed an ultralow thermal conductivity (0.47 Wm -1 K -1 at room temperature) due to lattice anharmonicity, a low speed of sound and a low Debye temperature. Moreover, first principles electronic structure calculations revealed that the Ag-Se tetrahedra within the crystal structure possess relatively weak bonds due to occupied antibonding states from p-d orbital hybridization, which suppress the thermal conductivity. Here, in order to further quantify our findings, we extended our analyses to include comparisons with other ternary and quaternary adamantine materials, all of which possess a fourfold tetrahedral coordination of atoms, in revealing the origin of the thermal properties in AgInSnSe 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

A synergistic experimental and theoretical study elucidating the electronic and thermal properties in spinel CuInSnS 4

Quaternary chalcogenides of different structure types continue to be of interest due to the novel physical properties they exhibit and for applications ranging from optoelectronics to energy-related technologies. Herein we report on the electronic and thermal properties of spinel CuInSnS 4 . UV–vis–NIR spectroscopy indicates strong absorption in the visible-light region with an indirect bandgap of 1.52 eV. First-principles computations corroborate our experimental data and reveal that the band-edge electronic properties are due to strong p – d ands – p hybridization and antibonding interactions. Analyses of temperature-dependent thermal properties together with first principles simulations reveal strong anharmonicity and low-lying optical branches that hybridize with the acoustic modes effectively suppressing thermal conductivity. This intrinsically very low thermal conductivity is much lower than that of related chalcogenides. Finally, our findings are discussed in the light of ongoing interest in quaternary metal chalcogenide materials and may aid in the development of this and similar chalcogenides for applications of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Systematic exploration of the thermochemistry for a set of peroxy hydroperoxy-alkyl radicals

Here, the thermochemistry of peroxy hydroperoxy-alkyl ($\overline{O}$OQOOH) radicals has a significant influence on the reactivity of fuels and on the formation of highly oxygenated molecules (HOMs) in the atmosphere. Theoretical characterization of these radicals can be arduous due to their molecular size and complex fundamental interactions, such as hydrogen-bonding and torsional anharmonicity, and difficult to validate in the absence of any direct experimental thermochemical data. In this work, we systematically explore the thermochemistry of a set of $\overline{O}$OQOOH radicals with five increasingly affordable approaches with considerations for these interactions. In doing so, we present a novel conformer selection approach suited to predict properties at combustion temperatures. As a corollary, we also highlight the shortcomings in the standard choice of the ground conformer as the reference. The set of molecules is comprised of 149 C 2 -C 8 $\overline{O}$OQOOH radicals, selected to encompass a wide variety of branching and substitution patterns. Comparisons amongst the approaches help quantify the errors arising from various simplifying assumptions. For the largest of these radicals, and with the most affordable approach, final predictions of Gibbs energies are assigned a 2 σ uncertainty of 4 kcal mol -1 in the negative temperature coefficient region.

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

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Leveraging Polymorphism in YbCuBi to Map Transport and Elastic Properties

AMX Zintl compounds with the hexagonal ZrBeSi structure have gained significant attention for their remarkable vacancy tolerance and low thermal conductivity. Their 2D honeycomb sublattice, composed of M–X covalent bonds, is believed to contribute to high anharmonicity and unusual thermal transport properties. In this study, we explore the temperature-dependent polymorphism of YbCuBi as a model system to investigate the relationship between the structure and elastic and thermal transport properties in AMX Zintls. YbCuBi undergoes a structural transition from the “flat” Cu–Bi layers in the ZrBeSi structure to corrugated layers in the LiGaGe structure below 410 K, resulting in a distortion of its centrosymmetric structure. To probe the effects of this crystallographic transition, we employ inelastic neutron scattering and temperature-dependent resonant ultrasound spectroscopy. These experimental findings, coupled with first-principles calculations and thermal conductivity measurements, allow us to elucidate a direct relationship between corrugation of the honeycomb lattice and the observed changes in elastic and thermal transport properties. These insights can be extended to other Zintl phases with similar structure types, providing a platform for the rational design of functional materials with tailored thermal properties.

Chemical structure↗
DOE OSTI2025

Interplay between Cation “Coloring” and Stereochemically Active Lone Pairs in AgBiS 2 Thin Films

Solution-processed AgBiS 2 thin films were fabricated using novel thiol− amine precursor inks to investigate the stereochemical activity of Bi 3+ 6s 2 lone pairs and their impact on the structure. A dual-space analysis combining Bragg diffraction and hard X-ray photoelectron spectroscopy (HAXPES) revealed a rock salt-like average structure with local distortions linked to cation coloring. Density functional theory (DFT) and crystal orbital Hamilton population (COHP) analyses confirmed that local Bi-rich and Ag-rich nanodomains amplify stereochemical activity, whereas more mixed and cation-order nanodomains are less stereochemically active. This local, nanoscopic mixing of segregated and ordered domains would indeed explain an average Fm3̅m structure that is rock salt-like and that does not manifest the full anharmonicity and noncentrosymmetry evidenced in canonical structures with stereochemical expression. These findings provide insights into the local structural and electronic complexities governing the optoelectronic properties of AgBiS 2 thin films.

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

Ab Initio-Based Bond Order Potential for Arsenene Polymorphs Developed via Hierarchical Reinforcement Learning

Arsenene, a less-explored two-dimensional material, holds the potential for applications in wearable electronics, memory devices, and quantum systems. This study introduces a bond-order potential model with Tersoff formalism, the ML-Tersoff, which leverages multireward hierarchical reinforcement learning (RL), trained on an ab initio data set. This data set covers a spectrum of properties for arsenene polymorphs, enhancing our understanding of its mechanical and thermal behaviors without the complexities of traditional models requiring multiple parameter sets. Our RL strategy utilizes decision trees coupled with a hierarchical reward strategy to accelerate convergence in high-dimensional continuous search spaces. Unlike the Stillinger-Weber approach, which demands separate formalisms for buckled and puckered forms, the ML-Tersoff model concurrently captures multiple properties of the two polymorphs by effectively representing the local environment, thereby avoiding the need for different atomic types. Here, we apply the ML model to understand the mechanical and thermal properties of the arsenene polymorphs and nanostructures. We observe an inverse relationship between the critical strain and temperature in arsenene. Thermal conductivity calculations in nanosheets show good agreement with ab initio data, reflecting a decrease in thermal conductivity attributable to increased anharmonic effects at higher temperatures. We also apply the model to predict the thermal behavior of arsenene nanotubes.

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

Infrared Fingerprint and Unimolecular Decay Dynamics of the Hydroperoxyalkyl Intermediate (•QOOH) in Cyclopentane Oxidation

A transient carbon-centered hydroperoxyalkyl intermediate (•QOOH) in the oxidation of cyclopentane is identified by IR action spectroscopy with time-resolved unimolecular decay to hydroxyl (OH) radical products that are detected by UV laser-induced fluorescence. Two nearly degenerate •QOOH isomers, β- and γ-QOOH, are generated by H atom abstraction of the cyclopentyl hydroperoxide precursor. Fundamental and first overtone OH stretch transitions and combination bands of •QOOH are observed and compared with anharmonic frequencies computed by second-order vibrational perturbation theory. An OH stretch transition is also observed for a conformer arising from torsion about a low-energy CCOO barrier. Definitive identification of the β-QOOH isomer relies on its significantly lower transition state (TS) barrier to OH products, which results in rapid unimolecular decay and near unity branching to OH products. A benchmarking approach is utilized to compute high-accuracy stationary point energies, most importantly TS barriers, for cyclopentane oxidation (C 5 H 9 O 2 ), building on higher level reference calculations for ethane oxidation (C 2 H 5 O 2 ). The experimental OH product appearance rates are compared with computed statistical microcanonical rates using RRKM theory, including heavy-atom tunneling, thereby validating the computed TS barrier. The results are extended to thermal unimolecular decay rate constants at temperatures and pressures relevant to cyclopentane combustion via master-equation modeling. Furthermore, the various torsional and ring puckering states of the wells and transition states are explicitly considered in these calculations.

Chemical calculations↗
DOE OSTI2025

Discrepancies between Theory and Experiment in Determining the Ionization Energy of NF 3

High-accuracy ab initio thermochemical predictions for the ionization energy of NF 3 , the barrier height (to inversion) of NF 3 + , and the dissociative ionization threshold of NF 3 to NF 2 + + F are presented and incorporated into Active Thermochemical Tables. The adiabatic ionization energy of the first ionization band of NF 3 , calculated at 12.647 ± 0.010 eV, is at odds with previous experimental interpretations by nearly 0.36 eV due to unfavorable Franck-Condon factors associated with this transition. The barrier (to inversion) height is calculated to be about 0.6 eV lower in energy than the prior interpretation, which instigates a discussion of the supposed vibrational structure of the first ionization band of NF 3 . Updated assignments of the photoelectron spectrum are proposed, and the loss in vibrational spacing on the high-energy side of the experimental ionization band is discussed. Rudimentary anharmonic Franck-Condon simulations qualitatively reproduce the broad spectral features observed in experiment.

Active Thermochemical Tables↗
DOE OSTI2025

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2025

Characterization of Two Positional Isomers of the Cs + Gly Complex Using Two-Color, IR–IR Photobleaching of the Cryogenically Cooled Ions

Metal ion binding to amino acid residues is an important interaction motif that controls the tertiary structures of oligopeptides. Analyses of the vibrational band patterns displayed by the amino acid scaffolds are commonly used to characterize the local docking motifs. Here we carry out two-color, IR-IR photobleaching measurements to obtain isomer-selective vibrational spectra of the Cs + Gly ion-molecule complex isolated in a cryogenically cooled, radiofrequency ion trap. The distinct band patterns of two non-interconverting isomers are observed and traced to different bidentate binding motifs between Cs + and the glycine scaffold. In one isomer, the ion attaches to the oxygen atoms of the carboxyl group whereas in the other it docks to the amino nitrogen and the carbonyl oxygen. Attachment to the acid head group yields a very diffuse absorption associated the OH group engaged in a strong intramolecular H-bond that closes a 5 membered ring. Furthermore, the band assignments, rearrangement pathways and electrostatic distortion of the electron density distributions in the glycine scaffold by the proximal ion are explored with electronic structure calculations and anharmonic theory.

Infrared spectroscopy↗
DOE OSTI2026

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2026

Signatures of Antisymmetric Vibrations in the Ultrafast Dynamics of Quadrupolar Dyes

Antisymmetric molecular vibrations are central to ultrafast, nonadiabatic photophysical and photochemical processes such as conical intersection dynamics, Herzberg–Teller couplings and, potentially, singlet fission. Direct spectroscopic identification of such vibrations is, however, challenging, since they are typically Raman inactive and affect optical transitions only weakly. Here, we report experimental signatures of vibronic couplings to a high-frequency antisymmetric vibration in the excited state dynamics of a quasi-quadrupolar molecule by ultrafast two-dimensional electronic spectroscopy (2DES). The early time, sub-50 fs 2DES maps reveal an asymmetric peak pattern with characteristic low-energy cross-peaks. We show that these peaks arise from stimulated emission transitions from an anharmonic, double-minimum excited state potential energy surface formed by vibronic coupling to a high-frequency antisymmetric mode. Simulations based on a phenomenological essential state model support the results. Our findings offer a new approach for identifying antisymmetric vibrations in ultrafast 2DES and track excited state wavepacket motion before relaxation washes out the spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Vibrational Fermi Resonance in Atomically Thin Black Phosphorus

Fermi resonance is a phenomenon involving the hybridization of two coincidentally quasi-degenerate states that is observed in the vibrational or electronic spectra of molecules. Despite numerous examples in molecular systems, vibrational Fermi resonances in dispersive semiconducting systems remain largely unexplored due to the rarity of occurrence. Here we report a vibrational Fermi resonance in atomically thin black phosphorus. The Fermi resonance arises via anharmonic mixing of a fundamental Raman mode and a Davydov component of an infrared mode, leading to a doublet with mixed character. The extent of Fermi coupling can be modulated by the application of external biaxial strain. The consequences of Fermi hybridization are revealed by electronic resonance effects in the thickness-dependent and excitation-wavelength-dependent Raman spectrum, which is predicted by ab initio hybrid functional simulations including excitonic interactions. This work reveals new insight into electron–phonon coupling in black phosphorus and demonstrates a novel method for modulating Fermi resonances in 2D semiconductors.

74 ATOMIC AND MOLECULAR PHYSICS↗