Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anharmonicity”

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

Water Catalysis of the Reaction of Methanol with OH Radical in the Atmosphere is Negligible

Abstract Faced with the contradictory results of two recent experimental studies [Jara‐Toro et al., Angew. Chem. Int. Ed . 2017 , 56 , 2166 and Chao et al., Angew. Chem. Int. Ed . 2019 , 58 , 5013] of the possible catalytic effect of water vapor on CH 3 OH + OH reaction, we report calculations that corroborate the conclusion made by Chao et al. and extend the rate constant evaluation down to 200 K. The rate constants of the CH 3 OH + OH reaction catalyzed by a water molecule are computed as functions of temperature and relative humidity using high‐level electronic structure and kinetics calculations. The Wuhan–Minnesota Scaling (WMS) method is used to provide accurate energetics to benchmark a density functional for direct dynamics. Both high‐frequency and low‐frequency anharmonicities are included. Variational and tunneling effects are treated by canonical variational transition state theory with multidimensional small‐curvature tunneling. And, most significantly, we include multistructural effects in the rate constant calculations. Our calculations show that the catalytic effect of water vapor is not observable at 200–400 K.

Wu, Junjun↗
DOE OSTI2023

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

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

Exploring the Tl 2H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Abstract Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet TlH potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 MH (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Chemistry↗
DOE OSTI2024

Out-of-time-ordered-correlators for the pure inverted quartic oscillator: classical chaos meets quantum stability

Out-of-time-ordered-correlators (OTOCs) have been suggested as a means to diagnose chaotic behavior in quantum mechanical systems. Recently, it was found that OTOCs display exponential growth for the inverted quantum harmonic oscillator, mirroring the fact that this system is classically and quantum mechanically unstable. In this work, I study OTOCs for the inverted anharmonic (pure quartic) oscillator in quantum mechanics, finding only oscillatory behavior despite the classically unstable nature of the system. For higher temperature, OTOCs seem to exhibit saturation consistent with a value of –2< x 2 > T < p 2 > T at late times. I provide analytic evidence from the spectral zeta-function and the WKB method as well as direct numerical solutions of the Schrödinger equation that the inverted quartic oscillator possesses a real and positive energy eigenspectrum, and normalizable wave-functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗
DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗
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 OSTI2023

Systematically derived thermodynamic properties for alkane oxidation

Key combustion properties, such as ignition delay time, show strong sensitivity to the thermochemistry of the main species in the standard radical oxidation pathway at low temperatures (600-1000K). Significant uncertainties persist in current estimates of thermodynamic properties, particularly for the larger species of relevance to common practical fuels. In this article, we use advanced computational schemes to evaluate thermodynamic properties for the fuel, fuel radical, peroxy, hydroperoxy-alkyl radical, and hydroperoxide species for a set of 17 fuels containing up to 9 carbon atoms, and with various degrees of branching in the alkane backbone. The procedure, termed STAR-1D, combines conformer sampling to find the minimum geometry, B2PLYP-D3/cc-pVTZ harmonic frequency evaluations, and omega B97X-D/cc-pVTZ one-dimensional torsional mapping. It includes two physically-based scaling routines: a frequency-dependent scaling to B2PLYP-D3/cc-pVTZ anharmonic frequencies and a scaling of one-dimensional omega B97X-D/ccpVTZ torsional profiles to reproduce the product of the B2PLYP-D3/cc-pVTZ frequencies in the harmonic limit. Substantive comparisons with existing experimental databases, together with careful examinations of key theoretical assumptions, are used to explore the accuracy of the predictions. These computationally intensive explorations of 195 species were facilitated by automated thermochemistry software. Further, high accuracy 0 K heats of formation from a separate study are used with the STAR-1D computations to generate NASA polynomial representations. In addition to their intrinsic value, the present results also provide a reliable database for the optimization of group additivity or machine learning schemes for scaling to larger combustion systems. Towards this end a complementary extensive conformational analysis is carried out for the medium sized species and the thermodynamic properties of the lowest energy hydrogen-bonded and non-hydrogen bonded conformers are contrasted for larger hydroperoxy-alkyl radical species.

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

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

Neutron scattering study of tantalum monohydride and monodeuteride

Powder samples of TaH 0.89 and TaD0.96 are synthesized under a hydrogen (deuterium) pressure of 2.8 GPa and a temperature of 250 °C, then quenched to the liquid nitrogen temperature, recovered to ambient pressure and studied by neutron diffraction (ND) and inelastic neutron scattering (INS). The ND study shows that both hydrogen and deuterium atoms occupy tetrahedral interstitial sites in a distorted body centered cubic (bcc) crystal structure of metal atoms, while the ordering scenarios in TaH 0.89 and TaD 0.96 are different. Hydrogen and deuterium atoms are ordered in a layered fashion, forming long period superstructures with space groups P$\bar{4}$ and P222, respectively, so that the unit cells of the crystal structures of TaH 0.89 and TaD 0.96 are $\sqrt{2}$ × $\sqrt{2}$ × 7 and $\sqrt{2}$ × $\sqrt{2}$ × 8 supercells of the initial cubic unit cell. The INS study demonstrates a pronounced “soft” (trumpet-like) anharmonicity of the potential well for H and D atoms.

08 HYDROGEN↗
DOE OSTI2023

Thermalization of the Kerr index of refraction in acetone and methanol using femtosecond pump-probe scattering spectroscopy

Femtosecond time-resolved pump-probe scattering is used to temporally investigate the index thermalization Δn(t) of the molecular motions associated with the Kerr nonlinear index n 2 . This is done using an intense 300 fs 1034 nm pump pulse and temporally probed by a 517 nm and supercontinuum pulse in acetone and methanol. The optically pumped molecular states change the index of refraction and through the processes of Raman and Rayleigh scattering are shown to engage in a 4-step energy relaxation process. Here, the pump first initiates the Kerr effect as well as populate resonant vibrational molecular motions, which act as a “mother” energy state. The electronic component of the Kerr effect is shown to couple with these vibrational modes through the process of Born-Oppenheimer coupling, which can be measured through the Raman scattering of the probe. The “mother” coupled state then decays into “daughter” anharmonic non-resonant states, which then further decay into “granddaughter” and thermal bath states. The population of these states alter the index of refraction in time, Δn(t), and can be measured through the changing in the probe beam’s Rayleigh scattering. This thermalization process takes ~5.5 ps in acetone and ~3.7 ps in methanol. The scattering signals from each Δn(t) thermalization stage in the localized region, (“mother”, to “daughter”, to “granddaughter”, to the bath states, to the ground state), are fitted to theoretical decay equations that show the rise and fall times of the energy decay routes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Synergistic role of Sn-doping on the thermal and electrical properties of sinnerite Cu 6 As 4 S 9

Sinnerite Cu 6 As 4 S 9 has been identified as a promising material for optoelectronic applications, with potential for other energy-related applications; however, knowledge of the electrical and thermal properties as well as the mechanisms underlying transport in sinnerite is lacking. Here, we present an investigation of the synthesis, structural, thermal and electrical transport properties of stoichiometric and Sn-doped Cu 6 As 4 S 9 . Sinnerite has a triclinic lattice structure, with highly distorted local atomic coordination environments that, in part, results in its complex structure and bonding. Our results and analyses indicate As 4s 2 lone pair-induced distortions and strong lattice anharmonicity that leads to a relatively short phonon mean free path, resulting in intrinsically very low thermal conductivity. The electrical resistivity for both compositions, Cu 6 As 4-x Sn x S 9 (x = 0, 0.2), are relatively high and varies little with temperature, typical of degenerate semiconductors. An increase in mobility and electrical conductivity was obtained by Sn doping in Cu 6 As 3.8 Sn 0.2 S 9 . This work demonstrates an effective route to synthesize bulk sinnerite as well as advances the knowledge of the properties of sinnerite, as this and other ternary chalcogenides continue to be of interest for potential technologically significant applications.

36 MATERIALS SCIENCE↗
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 OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis of several high-resolution infrared bands of perdeutero-spiropentane, C5D8

Perdeutero-spiropentane (C5D8) has been synthesized and infrared spectra were taken from 400 to 4000 cm-1 at low sample pressures in a 22 meter multipass cell. The spectral resolution (0.0015 cm-1) reveals thousands of individual rovibrational transitions for several fundamental bands occurring in the spectral region of 500 to 1100 cm-1. C5D8 ground state constants reported in a previous study of a ?15 parallel band at 1054 cm-1 [1] were improved significantly in the present work by adding combination-differences derived from ?16 parallel and ?24 perpendicular band analyses. The larger data set at higher J values also allowed determination of the ground state K? = 2 splitting constant (?2? = 0.76(2) ×10-9 cm-1). For ?16, the band origin is 847.16933(5) cm-1 and changes in upper state rotational constants are ?A = –0.0013411(5) and ?B = –0.0002021(1) (in units of cm-1). Overlapping, and in Fermi resonance with ?16, is a second parallel band which, from Gaussian calculations at the anharmonic level, is assigned as the b2 combination ?8+?12, with band origin of 840.5434(1) and ?A = –0.0003500(7) and ?B = –0.0003502(1) cm-1. For ?24 analogous values are 556.41453(2) and ?A = –0.00032295(9) and ?B = –0.00016782(4) cm-1. The analyses also yields accurate values for the changes in both quartic and sextic centrifugal distortion constants and, for ?24, values for the l-doubling constants q+ and q-. Theoretical calculations of the rovibrational parameters aided in the analysis and the numerical values are in good agreement with experimental results.

Nibler, Joseph W.↗