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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Vibrational Spectroscopy of Ions and Radicals Present in the Interstellar Medium and in Planetary Atmospheres: A Theoretical Study

Anharmonic vibrational frequencies and intensities are calculated for OH(H2O)n and H(H2O)n radicals (that form on icy particles of the interstellar medium), HCO radical (the main intermediate in the synthesis of organic molecules in space), NH2(-) and C2H(-) anions, H5(+) cation, and other systems relevant to interstellar chemistry. In addition to pure ions and radicals, their complexes with water are studied to assess the effects of water environment on infrared spectra. The calculations are performed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. Fundamental, overtone, and combination excitations are computed. The results are in good agreement with available experimental data and provide reliable predictions for vibrational excitations not yet measured in laboratory experiments. The data should be useful for interpretation of astronomically observed spectra and identification of ions and radicals present in the interstellar medium and in planetary atmospheres.

Chaban, Galina M.↗

NASA NTRS1996

Effects of H2O Vapor on Vibrational Relaxation in Expanding and Contracting Flows

As opposed to previous explanations based on the effects of anharmonicity of simple diatomic molecules, traces of water vapor are suggested to be the most likely cause of the anomalously fast vibrational relaxation of such gases observed in supersonic and hypersonic nozzles. The mechanism is the strong V-VR coupling with H2O molecules that dramatically facilitates the collisional transfer of vibrational energy. Slight moisture content is thus a real world aspect of gas dynamics that must be considered in characterizations of shock tubes, reflected shock tunnels, and expansion tubes.

Meador, Willard E.↗

NASA NTRS2004

Quantitative Analysis of Spectral Interference of Spontaneous Raman Scattering in High-Pressure Fuel-Rich H2-Air Combustion

We present a theoretical study of the spectral interferences in the spontaneous Raman scattering spectra of major combustion products in 30-atm fuel-rich hydrogen-air flames. An effective methodology is introduced to choose an appropriate line-shape model for simulating Raman spectra in high-pressure combustion environments. The Voigt profile with the additive approximation assumption was found to provide a reasonable model of the spectral line shape for the present analysis. The rotational/vibrational Raman spectra of H2, N2, and H2O were calculated using an anharmonic-oscillator model using the latest collisional broadening coefficients. The calculated spectra were validated with data obtained in a 10-atm fuel-rich H2-air flame and showed excellent agreement. Our quantitative spectral analysis for equivalence ratios ranging from 1.5 to 5.0 revealed substantial amounts of spectral cross-talk between the rotational H2 lines and the N2 O-/Q-branch; and between the vibrational H2O(0,3) line and the vibrational H2O spectrum. We also address the temperature dependence of the spectral cross-talk and extend our analysis to include a cross-talk compensation technique that removes the nterference arising from the H2 Raman spectra onto the N2, or H2O spectra.

Kojima, Jun↗

NASA NTRS2014

Vibrational Frequencies and Spectroscopic Constants for 1(sup 3)A' HNC and 1(sup 3)A' HOC+ from High-Accuracy Quartic Force Fields

The spectroscopic constants and vibrational frequencies for the 1(sup 3)A' states of HNC, DNC, HOC+, and DOC+ are computed and discussed in this work. The reliable CcCR quartic force field based on high-level coupled cluster ab initio quantum chemical computations is exclusively utilized to provide the anharmonic potential. Then, second order vibrational perturbation theory and vibrational configuration interaction methods are employed to treat the nuclear Schroedinger equation. Second-order perturbation theory is also employed to provide spectroscopic data for all molecules examined. The relationship between these molecules and the corresponding 1(sup 3)A' HCN and HCO+ isomers is further developed here. These data are applicable to laboratory studies involving formation of HNC and HOC+ as well as astronomical observations of chemically active astrophysical environments.

Astrochemistry↗

NASA NTRS2017

Thermal Expansion from Stochastic Nonlinear Acoustic Fields

The methods of stochastic mechanics are applied to nonlinear acoustic fields generated by random acoustic radiation sources associated with oscillating lattice sites in nonlinear lattices. The assumption of stochastically independent zero-point and temperature-dependent acoustic fields leads to expressions of the thermodynamic internal and Helmholtz free energies per unit mass in terms of modal energies per unit mass that account for the nonlinearity of the propagation modes. The thermodynamic state functions canonically transform to the familiar results of a system of quantized, simple harmonic oscillators in the linear field limit. The thermal expansion coefficient, derived from the Helmholtz free energy per unit mass, is obtained as a sum of the zero-point modal nonlinearity parameters, weighted by temperature-dependent, modal heat capacities. The relationship is fully anharmonic, in contrast to commonly used quasi-harmonic models, and predicts with excellent agreement the experimentally observed null thermal expansion of vitreous silica at the ‘cross-over’ temperature corresponding to the balance between long wavelength modes that contribute negative nonlinearity parameters to the modal sum and short wavelength modes that contribute positive nonlinearity parameters.

Cantrell, John H., Jr.↗

NASA NTRS2019

Thermodynamic Analysis of Boundary Layer Species During Ablation

The focus of this work is on carbon clusters, hydrocarbons and other carbonaceous species which are constituents of pyrolysis gases injected into the boundary layer. Some of these molecules are observed to have absorption spectra in the VUV and UV region that match the emission spectra of atomic nitrogen and oxygen. Hence, they can potentially absorb the radiation impinging on the heat shield. Accurate thermodynamic data is not available for many of these species over the relevant temperature range for reentry. The objective of this work is to determine the thermochemical properties of potential radiation absorbing molecules using accurate ab initio quantum chemistry calculations. We have compiled a database of formation enthalpies, anharmonic vibrational frequencies and rotational constants for the ground and low-lying excited electronic states of these carbonaceous species. Using these thermochemical properties, we compute equilibrium mole fractions and observe that the mole fractions of some species of interest are significant, thus they could affect the radiative heat flux. Furthermore, comparisons of the formation enthalpies and mole fractions are made with data from literature.

heatshields↗

NASA NTRS2020

A new model of monodeuterated ethane (C2H5D) spectrum: Enabling sensitive constraints on the D/H in ethane emission in comets

Similar to ethane (C2H6), monodeuterated ethane (C2H5D) has a bright series of Q-branches in the CH stretch region, between about 2960 - 3020 cm (3.38 - 3.31 µm). As the progressions in the different isotopes are slightly offset from one another, the new C2H5D model enables sensitive constrains of the D/H ratio in spectra from comets, in organics across the Solar System and beyond. Specifically, the D/H ratio of organic molecules in comets is a critical cosmogenic indicator, which not only provides insights into the formation processes of our Solar System, but also allows for a test of our understanding of solid-state astrochemistry. High D/H ratios are predicted in a number of astrochemical models, which can be benchmarked using the new spectroscopic information of C2H5D and C2H6. We present a quantum band model for C2H5D in the mid-infrared. Our model is based on rotationally resolved spectra of C2H5D and C2H6 obtained at 85K between 2050 - 3050 per cm (4.88 - 3.28 µm) using a Bruker IFS-125HR equipped with a cryogenic Herriott cell at JPL. High-level quantum chemical calculations were carried out at the CCSD(T)/ANO1 level of theory to aid vibrational assignments, with the anharmonic frequencies and vibrational corrections determined from second-order vibrational perturbation theory (VPT2). As a first demonstration of the applicability of the model, we obtained a stringent upper limit to the organic D/H ratio (18.2 VSMOW) in comet C/2007 W1 (Boattini) from high resolution fluorescence spectra obtained with Keck/NIRSPEC.

Kirstin D. Doney↗

NASA NTRSDate not supplied

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Recent rapid progresses in physics theory and computational power have made it possible to predict the martensitic transition temperatures (MTTs) in shape memory alloys (SMAs) from first principles [1-3]. In particular, rigorous while time-consuming thermodynamic integration has been employed to compute the anharmonic phonon free energies, which play a crucial role in determining martensitic phase transitions in SMAs. However, this approach has only been applied to simple binaries, and its accuracy is unsatisfying for certain SMAs such as the most commonly used NiTi. In this work, we report on several new developments to our method that bring first-principles theory and experiment much closer into agreement including the MTT of NiTi, and that improve the computational efficiency significantly. We have applied our refined approach to investigate the Ni0.5Ti0.5-xHfx and PdxNi0.5-xTi0.5 ternaries, and the predicted MTT for each composition is within 100K compared with experiment. We will address various techniques to overcome the difficulty encountered in studying ternaries. Our theoretical approach is expected to be a broadly applicable and predictive theory for designing complex SMAs with desirable properties. [1] J.B. Haskins, A.E. Thompson,and J.W. Lawson, Phys. Rev B 94, 214110 (2016). [2] J.B. HaskinsandJ.W. Lawson, J. App. Phys. 121, 205103 (2017). [3] J.B. Haskins, H. Malmir, S. J. Honrao, L. A. Sandoval, and J.W. Lawson, Acta Materialia 212, 116872 (2017).

Zhigang Wu↗

NASA NTRSDate not supplied

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Zhigang Wu↗

NASA NTRS2026

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly↗

DOE OSTI2020

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗

DOE OSTI2021

Excitation dynamics in inductively coupled fluxonium circuits

We propose a near-term quantum simulator based on the fluxonium qubits inductively coupled to form a chain. This system provides long coherence time, large anharmonicity, and strong coupling, making it suitable to study Ising spin models. At the half-flux quantum sweet spot, the system is described by the transverse field Ising model (TFIM). We evaluate the propagation of qubit excitations through the system. As disorder increases, the excitations become localized. A single qubit measurement using the circuit QED methods is sufficient to identify localization transition without introducing tunable couplers. We argue that inductively coupled fluxoniums provide opportunities to study localization and many-body effects in highly coherent quantum systems.

Barış Özgüler, A.↗

DOE OSTI2025

Lattice-Renormalized Tunneling Models for Superconducting Qubit Materials

We present a lattice-renormalized formalism for configurational tunneling two-level systems (TLS) that overcomes limitations of minimum-energy-path and light-particle models. Derived from the nuclear Hamiltonian, our formulation introduces composite phonon coordinates to capture lattice distortions between degenerate potential wells. This approach resolves deficiencies in prior models and enables accurate computation of tunnel splittings and excitation spectra for hydrogen-based TLS in bcc Nb. Our results bound experimental tunnel splittings and reveal strong anharmonic couplings between tunneling atoms and lattice phonons, establishing a direct link between TLS dynamics and phonon-mediated strain interactions. The formalism further generalizes to multi-level systems (MLS), providing insight into defect-induced decoherence in superconducting qubits and guiding strategies for materials design to suppress TLS-related loss.

Pritchard, P. G. [Northwestern U.]↗

DOE OSTI2026

Neural-Network Inverse Design of SRF Cavities and Transmons for Bosonic Quantum Computation

Three-dimensional superconducting radio-frequency (SRF) cavities provide exceptionally long-lived electromagnetic modes and, when coupled to nonlinear elements such as transmon qubits, become promising architectures for bosonic quantum information processing. The inverse design of such systems, i.e., recovering device geometries that produce specified electromagnetic and coupling targets, is generally a one-to-many problem. The qubit-cavity coupling strength depends sensitively on both the transmon geometry and its position within the cavity's electromagnetic field. As these systems scale up and their design parameter spaces grow, the cost of conventional iterative simulation becomes prohibitive. We present two deep neural network (DNN) approaches that address this inverse-design problem at complementary levels of the design stack. The first proposes SRF cavity geometries that produce target cavity observables. The second proposes transmon qubit designs that produce target qubit-cavity parameters - the coupling rate, qubit frequency, and anharmonicity $(g, ν_q, α)$. The recovered candidate designs match the targets to within ~5% (cavity) and ~2% (transmon), confirmed by end-to-end re-simulation. Both approaches map desired device behavior directly to candidate designs, a fast alternative to the iterative simulation studies usually required.

Yaker, Joseph [Fermilab; Northwestern U.]↗

DOE OSTI2026

Component-Level Inverse Design of Transmon Qubits Using Neural Networks

Designing a superconducting qubit to realize specific Hamiltonian parameters typically requires iterating through a time and compute-intensive forward loop in which the designer chooses a layout geometry, simulates it, extracts circuit parameters such as capacitances, and refines the geometry. We study the inverse version of this task using a neural-network workflow that maps target Hamiltonian parameters directly to component-level layout parameters, which we subsequently demonstrate on a planar transmon layout. During training, we pair the inverse model with a frozen forward surrogate model and evaluate the loss in Hamiltonian space rather than in layout-parameter space. In validation against a conventional EM solver, 97% of generated designs produce usable geometries, and the inverse-plus-surrogate pipeline reaches mean percent errors of 0.73% for qubit frequency and 1.58% for anharmonicity, comparable to or below the fabrication and simulation-to-measurement uncertainty expected for academic-process transmon devices of this type. A single pipeline query takes ~60 ms on CPU, versus ~2 min for a conventional EM capacitance extraction on the same hardware, a speedup of approximately 2,000x. Batching minimizes the AI model inference overhead, reducing the runtime to 3.1 microseconds per sample on CPU and 2.6 microseconds per sample on GPU at a batch size of 2048, resulting in speedups of 3.9 x 10^7 and 4.6 x 10^7, respectively, relative to a single conventional CPU EM extraction. Our results indicate that component-level inverse design usefully extends and complements conventional EM simulation, including for small datasets on the order of 1,000 samples.

Seidel, Olivia [Fermilab; Texas U., Arlington]↗

DOE OSTI2026

Transition-state lattice modes and the breakdown of adiabatic tunneling for hydrogen and deuterium in bcc Nb

Interstitial hydrogen and deuterium in body-centered-cubic metals constitute archetypal quantum tunneling systems. Their relevance has been renewed by the connection between hydrogenic tunneling in Nb and defect-induced decoherence in superconducting qubits, motivating a predictive microscopic theory. Existing theoretical treatments invoke an adiabatic separation between the light interstitial and the host lattice, an assumption whose validity has not been rigorously established for hydrogenic species. Here, we show that the experimentally measured tunnel splittings of O-trapped H and D in bcc Nb are quantitatively reproduced only within a five-dimensional (5D) Lattice-Renormalized Born-Oppenheimer (LRBO) framework. This approach treats three interstitial modes and two judiciously selected lattice modes, which includes a transition-state mode, on equal quantum footing. By recasting nested Born-Oppenheimer hierarchies within this same formalism and benchmarking against modern \textit{ab initio} potential energy surfaces, we show that adiabatic separation of the light particle from lattice dynamics is satisfied only in the positive-muon ($μ^{+}$) mass limit. In contrast, tunneling for H and D is fundamentally a collective, nonadiabatic process mediated by anharmonic lattice couplings. Finally, we show that the breakdown of adiabaticity can be anticipated from simple energy estimates involving the ground-state light-particle energy evaluated at a small number of fixed lattice configurations, providing a practical criterion for assessing the validity of adiabatic tunneling theories in other systems.

Pritchard, P. Graham [Northwestern U.]↗

DOE OSTI2021

TUMME: Tsinghua University Minnesota Master Equation program

We report that TUMME is a program for assembling and solving master equations for gas-phase chemical kinetics based on chemically significant eigenmodes. TUMME has interfaces to the Gaussian, Polyrate, and/or MSTor output files that allow the master equation code to obtain the microcanonical flux coefficients needed for the coefficient matrix of the master equation. The flux coefficients for reactions with barriers can be calculated by multi-structural variational transition state theory with small-curvature tunneling (MS-VTST/SCT) or by simpler approximations to this such as conventional transition state theory without tunneling (also called RRKM theory). The flux coefficients for barrierless reactions are provided by a hard-sphere model. TUMME is written in double precision with Python 3; quadruple and octuple precision are also available for some subtasks in C++. The Python code can run in serial or parallel (MP or MPI), and the C++ code can run on a single processor or on multiple processors with OpenMP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Response of vibrational properties and thermal conductivity of perovskites to pressure

We study the response of SrTiO 3 and KTaO 3 's vibrational properties and thermal conductivity to pressurization. Our goal is to improve the understanding of the relationship between crystal structure, vibrational dynamics, and thermal conductivity in perovskites. We measure the thermal conductivity of SrTiO 3 and KTaO 3 up to 28 GPa by time-domain thermoreflectance. We also perform Raman scattering and stimulated Brillouin scattering measurements of SrTiO 3 and KTaO 3 to characterize changes in vibrational dynamics with pressure. The thermal conductivity of SrTiO 3 increases under pressure with a slope comparable to that of other perovskites whose thermal conductivity has been measured versus pressure. Alternatively, the thermal conductivity of KTaO 3 has a stronger pressure dependence than that of other materials with similar crystal structure. We correlate pressure-induced changes in Raman and Brillouin spectra with pressure-induced changes in thermal conductivity. We show that pressure-induced changes in phonon lifetimes dominate the pressure dependence of thermal conductivity. This study provides benchmark knowledge of why depends on pressure and improves understanding of structure/thermal-property relationships.

36 MATERIALS SCIENCE↗