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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Lead-free Zr-doped ceria ceramics with low permittivity displaying giant electrostriction

Electrostrictors, materials developing mechanical strain proportional to the square of the applied electric field, present many advantages for mechanical actuation as they convert electrical energy into mechanical, but not vice versa. Both high relative permittivity and reliance on Pb as the key component in commercial electrostrictors pose serious practical and health problems. Here we describe a low relative permittivity (<250) ceramic, Zr x Ce 1-x O 2 (x < 0.2), that displays electromechanical properties rivaling those of the best performing electrostrictors: longitudinal electrostriction strain coefficient ~10 -16 m 2 /V 2 ; relaxation frequency ≈ a few kHz; and strain ≥0.02%. Combining X-ray absorption spectroscopy, atomic-level modeling and electromechanical measurements, here we show that electrostriction in Zr x Ce 1-x O 2 is enabled by elastic dipoles produced by anharmonic motion of the smaller isovalent dopant (Zr). Unlike the elastic dipoles in aliovalent doped ceria, which are present even in the absence of an applied elastic or electric field, the elastic dipoles in Zr x Ce 1-x O 2 are formed only under applied anisotropic field. The local descriptors of electrostrictive strain, namely, the cation size mismatch and dynamic anharmonicity, are sufficiently versatile to guide future searches in other polycrystalline solids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photoinduced anisotropic lattice dynamic response and domain formation in thermoelectric SnSe

Identifying and understanding the mechanisms behind strong phonon–phonon scattering in condensed matter systems is critical to maximizing the efficiency of thermoelectric devices. To date, the leading method to address this has been to meticulously survey the full phonon dispersion of the material in order to isolate modes with anomalously large linewidth and temperature-dependence. Here we combine quantitative MeV ultrafast electron diffraction (UED) analysis with Monte Carlo based dynamic diffraction simulation and first-principles calculations to directly unveil the soft, anharmonic lattice distortions of model thermoelectric material SnSe. A small single-crystal sample is photoexcited with ultrafast optical pulses and the soft, anharmonic lattice distortions are isolated using MeV-UED as those associated with long relaxation time and large displacements. We reveal that these modes have interlayer shear strain character, induced mainly by c -axis atomic displacements, resulting in domain formation in the transient state. These findings provide an innovative approach to identify mechanisms for ultralow and anisotropic thermal conductivity and a promising route to optimizing thermoelectric devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Magnetically driven phonon instability enables the metal–insulator transition in h-FeS

Hexagonal iron sulfide exhibits a fascinating coexistence of metal–insulator, structural and magnetic transitions, reflecting an intimate interplay of its spin, phonon and charge degrees of freedom. Here, we show how a subtle competition of energetic and entropic free-energy components governs its thermodynamics and the sequence of phase transitions it undergoes upon cooling. By means of comprehensive neutron and X-ray scattering measurements, and supported by first-principles electronic structure simulations, we identify the critical role of the coupling between antiferromagnetic ordering and instabilities of anharmonic phonons in the metallic phase in driving the metal–insulator transition. The antiferromagnetic ordering enables the emergence of two zone-boundary soft phonons, whose coupling to a zone-centre mode drives the lattice distortion opening the electronic bandgap. Simultaneously, spin–lattice coupling opens a gap in the magnon spectrum that controls the entropy component of the metal–insulator transition free energy. Furthermore, these results reveal the importance of spin–phonon coupling to tune anharmonic effects, thus opening new avenues to design novel technologically important materials harbouring the metal–insulator transition and magnetoelectric behaviours.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Size distribution of polycyclic aromatic hydrocarbons in space: an old new light on the 11.2/3.3 μm intensity ratio

The intensity ratio of the 11.2/3.3 μm emission bands is considered to be a reliable tracer of the size distribution of polycyclic aromatic hydrocarbons (PAHs) in the interstellar medium (ISM). This paper describes the validation of the calculated intrinsic infrared (IR) spectra of PAHs that underlie the interpretation of the observed ratio. The comparison of harmonic calculations from the NASA Ames PAH IR spectroscopic database to gas-phase experimental absorption IR spectra reveals a consistent underestimation of the 11.2/3.3 μm intensity ratio by 34%. IR spectra based on higher level anharmonic calculations, on the other hand, are in very good agreement with the experiments. While there are indications that the 11.2/3.3 μm ratio increases systematically for PAHs in the relevant size range when using a larger basis set, it is unfortunately not yet possible to reliably calculate anharmonic spectra for large PAHs. Based on these considerations, we have adjusted the intrinsic ratio of these modes and incorporated this in an interstellar PAH emission model. This corrected model implies that typical PAH sizes in reflection nebulae such as NGC 7023 – previously inferred to be in the range of 50 to 70 carbon atoms per PAH are actually in the range of 40 to 55 carbon atoms. The higher limit of this range is close to the size of the C 60 fullerene (also detected in reflection nebulae), which would be in line with the hypothesis that, under appropriate conditions, large PAHs are converted into the more stable fullerenes in the ISM.

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

Low-cost calculation and analysis of 2D IR spectra of model diiron trinitrosyl complexes in the NO stretch region with vibrational perturbation theory

Two-dimensional infrared spectroscopy offers unique capabilities for probing vibrational coupling in complex metal–ligand systems. In this paper, we combine two-dimensional infrared spectroscopy with vibrational perturbation theory to investigate vibrational coupling in a diiron trinitrosyl complex across three stable redox states. Although these systems are challenging for electronic structure methods, we demonstrate that key features of experimental 2D IR spectra can be accurately reproduced using reduced-dimensional anharmonic calculations with a small harmonic frequency scaling. Analysis reveals that N–O stretching modes maintain high locality across all redox states, with coupling patterns that directly reflect variations in Fe–N bond strength. Using curvilinear coordinate analysis, we demonstrate these differences result from systematic changes in cubic anharmonic force constants rather than mode delocalization. Our results establish N–O stretches as sensitive probes of metal–ligand bonding strength, expanding the toolkit for studying biologically relevant nitrosyl complexes.

Moran, Hayden A. [Texas A&M University, College St↗
DOE OSTI2020

Vibrational spectroscopy of shock compressed condensed phase nitrous oxide: Frequency shifts in the ν 1 mode

Over many years, we have published vibrational spectra of molecular materials at extreme pressure and temperature obtained using shock compression and coherent Raman spectroscopy methods. For diatomic molecules, we were able to extract vibrational temperatures from the intensities of spectrally resolved vibrational hot bands. Larger molecules in the condensed phase suffer band broadening effects that obscure the vibrational hot bands as the anharmonicities are typically smaller than the widths of the bands. This inability to resolve the hot bands inhibits the extraction of vibrational temperature and the measurement of the vibrational frequency of the fundamental. Here, we use a hot band model based on gas-phase anharmonic coupling coefficients to fit coherent anti-Stokes Raman spectra of the ν 1 vibrational mode of shock compressed condensed phase N 2 O with shock pressures and temperatures estimated from literature equations of state and compare to fits from a model using a single Gaussian peak. We report the resulting vibrational frequency shifts with shock pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental and computational studies of the kinetics of the reaction of hydrogen peroxide with the amidogen radical

In this work, the pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH 2 with H 2 O 2 to yield a second-order rate constant of (2.42 ± 0.55) × 10 –14 cm 3 molecule –1 s –1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10 –15 (T/298 K) 4.02 exp(–513 K/T) cm 3 molecule –1 s –1 over the temperature range of 298–3000 K.

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

Helium droplet infrared spectroscopy of the butyl radicals

Butyl radicals (n-, s-, i-, and tert-butyl) are formed from the pyrolysis of stable precursors (1-pentyl nitrite, 2-methyl-1-butyl nitrite, isopentyl nitrite, and azo-tert-butane, respectively). The radicals are doped into a beam of liquid helium droplets and probed with infrared action spectroscopy from 2700 to 3125 cm -1 , allowing for a low temperature measurement of the CH stretching region. The presence of anharmonic resonance polyads in the 2800–3000 cm -1 region complicates its interpretation. To facilitate spectral assignment, the anharmonic resonances are modeled with two model Hamiltonian approaches that explicitly couple CH stretch fundamentals to HCH bend overtones and combinations: a VPT2+K normal mode model based on coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] quartic force fields and a semi-empirical local mode model. Both of these computational methods provide generally good agreement with the experimental spectra.

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

Pitfalls in the n -mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Chemistry↗
DOE OSTI2024

Linear response of molecular polaritons

In this article, we show that the collective light–matter strong coupling regime where N molecular emitters couple to the photon mode of an optical cavity can be mapped to a quantum impurity model where the photon is the impurity that is coupled to a bath of anharmonic transitions. In the thermodynamic limit where N ≫ 1, we argue that the bath can be replaced with an effective harmonic bath, leading to a dramatic simplification of the problem into one of the coupled harmonic oscillators. We derive simple analytical expressions for linear optical spectra (transmission, reflection, and absorption) where the only molecular input required is the molecular linear susceptibility. This formalism is applied to a series of illustrative examples, showing the role of temperature, disorder, vibronic coupling, and optical saturation of the molecular ensemble, explaining that it is useful even when describing an important class of nonlinear optical experiments. For completeness, we provide Appendixes A–C that include a self-contained derivation of the relevant spectroscopic observables for arbitrary anharmonic systems (for both large and small N) within the rotating-wave approximation. While some of the presented results herein have already been reported in the literature, we provide a unified presentation of the results as well as new interpretations that connect powerful concepts in open quantum systems and linear response theory with molecular polaritonics.

Chemistry↗
DOE OSTI2023

Stress-dependent activation entropy in thermally activated cross-slip of dislocations

Cross-slip of screw dislocations in crystalline solids is a stress-driven thermally activated process essential to many phenomena during plastic deformation, including dislocation pattern formation, strain hardening, and dynamic recovery. Molecular dynamics (MD) simulation has played an important role in determining the microscopic mechanisms of cross-slip. However, due to its limited timescale, MD can only predict cross-slip rates in high-stress or high-temperature conditions. The transition state theory can predict the cross-slip rate over a broad range of stress and temperature conditions, but its predictions have been found to be several orders of magnitude too low in comparison to MD results. This discrepancy can be expressed as an anomalously large activation entropy whose physical origin remains unclear. Here, we resolve this discrepancy by showing that the large activation entropy results from anharmonic effects, including thermal softening, thermal expansion, and soft vibrational modes of the dislocation. We expect these anharmonic effects to be significant in a wide range of stress-driven thermally activated processes in solids.

Science & Technology - Other Topics↗
DOE OSTI2025

Origin of Intrinsically Low Thermal Conductivity in a Garnet-Type Solid Electrolyte: Linking Lattice and Ionic Dynamics with Thermal Transport

Understanding thermal transport in solid electrolytes is essential for improving the performance, reliability, and safety of all-solid-state batteries. Garnet-type lithium-ion conductors are promising candidates for solid electrolytes, yet their thermal-transport mechanisms remain poorly understood. Here, we connect the lattice and ion dynamics of single-crystal garnet-type Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 to its intrinsically low thermal conductivity. Our study reveals that the single crystals grown by the floating-zone method exhibit remarkably low glasslike thermal conductivity. Using first-principles calculations and inelastic-neutron-scattering measurements, we identify both the acoustic and numerous optical phonon modes, which stem from the complex crystal structure of the material. Notably, a low-energy optical branch exhibits an avoided crossing with acoustic phonons near 7 meV. These optical modes can enhance the scattering of heat-carrying acoustic phonons and reduce thermal conductivity. Furthermore, the calculated Grüneisen parameters are large, especially for the vibrational modes around 6 meV, indicating strong anharmonicity, with a noticeable contribution from lithium-ion vibrations. A two-channel thermal-transport model is employed to describe the weak temperature dependence of the thermal conductivity, which can be attributed to the substantial contribution of diffuson transport facilitated by the abundance of optical phonons and intrinsic anharmonicity. These results offer valuable insights into the thermal transport in a broad class of ionic conductors of interest for energy conversion and storage applications.

Ab initio calculations↗
DOE OSTI2022

One- and two-qubit gate infidelities due to motional errors in trapped ions and electrons

In this work, we derive analytic formulas that determine the effect of error mechanisms on one- and two-qubit gates in trapped ions and electrons. First, we analyze and derive expressions for the effect of driving field inhomogeneities on one-qubit gate fidelities. Second, we derive expressions for two-qubit gate errors, including static motional frequency shifts, trap anharmonicities, field inhomogeneities, heating, and motional dephasing. We show that, for small errors, each of our expressions for infidelity converges to its respective numerical simulation; this shows that our formulas are sufficient for determining error budgets for high-fidelity gates, obviating numerical simulations in future projects. All of the derivations are general to any internal qubit state, and any mixed state of the ion crystal's motion that is diagonal in the Fock state basis. Our treatment of static motional frequency shifts, trap anharmonicities, heating, and motional dephasing apply to both laser-based and laser-free gates, while our treatment of field inhomogeneities applies to laser-free systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum stabilization and flat hydrogen-based bands of nitrogen-doped lutetium hydride

We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.

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

Order-disorder versus displacive transitions in Jahn-Teller active layered materials

Large anharmonic vibrations often play a crucial role in dynamically stabilizing crystalline phases whose structures are unstable at low temperature. Although the average structure of such phases can be measured through diffraction experiments, their local structure remains a challenge to characterize and understand. Dynamically stabilized phases are often classified as order/disorder or displacive based on the qualitative nature of their local structure. A robust understanding of how chemistry determines this distinction in behavior, however, is lacking. This article presents a parametric study of an anharmonic vibrational model that describes the transition from a cooperative to a noncooperative Jahn-Teller distortion in layered oxides—a class of materials widely used in Li-ion and Na-ion batteries. The results illustrate how the shape of the energy landscape determines the extent to which the high-temperature phase has order/disorder vs displacive character. In this work, we find that the nature of the high-temperature phase is determined by a competition between the strength of the elastic coupling between Jahn-Teller distortions at nearby sites and the energy scale driving the Jahn-Teller distortion. A comparison of the model to energy landscapes calculated with density-functional theory suggests that the high-temperature phases of the Jahn-Teller active layered compounds LiNiO 2 , NaNiO 2 , LiMnO 2 , and NaMnO 2 exhibit order/disorder character.

36 MATERIALS SCIENCE↗