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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

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

Ultrafast infrared transient absorption spectroscopy of gas-phase Ni(CO) 4 photodissociation at 261 nm

Here, we employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin–orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

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

Electronic heat generation in semiconductors: Non-equilibrium excitation and evolution of zone-edge phonons via electron–phonon scattering in photo-excited germanium

We investigate experimentally and using first-principles theory the generation of phonons and the relaxation of carriers on picosecond timescales across the Brillouin zone of photo-excited Ge by inter-valley electron–phonon scattering. The phonons generated are typical of those generated in semiconductor devices, contributing to the accumulation of heat within the material. We simulate the time-evolution of phonon populations, based on first-principles band structure and electron–phonon and phonon–phonon matrix elements, and compare them to data from time-resolved x-ray diffuse scattering experiments, performed at the Linac Coherent Light Source x-ray free-electron laser facility, following photo-excitation by a 50 fs near-infrared optical pulse. We show that the intensity of the non-thermal x-ray diffuse scattering signal, which is observed to grow substantially near the L-point of the Brillouin zone over 3–5 ps, is due to phonons generated by scattering of carriers between the Δ and L valleys. These phonons have low group velocities, resulting in a heat bottleneck. With the inclusion of phonon decay through 3-phonon processes, the simulations also account for other non-thermal features observed in the x-ray diffuse scattering intensity, which are due to anharmonic phonon–phonon scattering of the phonons initially generated by electron–phonon scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Theoretical upper limits of the thermal conductivity of Si3N4

Silicon nitride (Si3N4) is a promising substrate for high-power electronics due to its superior mechanical properties and potential outstanding thermal conductivity (κ). As experiments keep pushing the upper limit of κ of Si3N4, it is believed that it can reach 450 W/mK, similar to SiC, based on classical models and molecular dynamics simulations. In this work, we reveal from first principles that the theoretical κ upper limits of β-Si3N4 are only 169 and 57 W/mK along the c and a axes at room temperature, respectively. Those of α-Si3N4 are about 116 and 87 W/mK, respectively. The predicted temperature-dependent κ matches well with the highest available experimental data, which supports the accuracy of our calculations, and suggests that the κ upper limit of Si3N4 has already been reached in the experiment. Compared to other promising semiconductors (e.g., SiC, AlN, and GaN), Si3N4 has a much lower κ than expected even though the chemical bonding and mechanical strengths are close or even stronger. We find the underlying reason is that Si3N4 has much lower phonon lifetimes and mean free paths (<0.5 μm) due to the larger three-phonon scattering phase space and stronger anharmonicity. Interestingly, we find that the larger unit cell (with more basis atoms) that leads to a smaller fraction of acoustic phonons is not the reason for lower κ. Grain size-dependent κ indicates that the grain boundary scattering plays a negligible role in most experimental samples. This work clarifies the theoretical κ upper limits of Si3N4 and can guide experimental research.

Physics↗
DOE OSTI2023

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
DOE OSTI2023

Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations

While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute–solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.

Chemistry↗
DOE OSTI2024

Acoustoelastic characterization of plates using zero group velocity Lamb modes

Acoustoelasticity, a characteristic of material anharmonicity, gives rise to a link between wave propagation velocity and the stress state in materials. Ultrasonic techniques to monitor this coupling, particularly with high sensitivity and in a noncontact manner, can have widespread application both in the quantification of applied and residual stress and in the characterization of nonlinear material behavior through measurement of higher order elastic constants. Here, we use a laser ultrasonic technique to excite and detect zero group velocity (ZGV) Lamb wave resonances in aluminum plates under uniaxial loading. A laser line source is used to excite these resonances at different orientations with respect to the applied load, and the signals are detected using an interferometer. The effects of stress and source orientation on ZGV resonance frequencies are validated using the theory of acoustoelastic Lamb wave propagation. In addition, a model-based inversion technique is used to extract Murnaghan's third-order elastic constants from measurements of the stress dependence of the first two ZGV modes generated parallel and perpendicular to the applied load. Laser generation and detection of ZGV resonances is shown to be an effective and powerful approach for the noncontact and nondestructive acoustoelastic characterization of elastic waveguides.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Infrared signature of the hydroperoxyalkyl intermediate (·QOOH) in cyclohexane oxidation: An isomer-resolved spectroscopic study

Infrared (IR) action spectroscopy is utilized to characterize carbon-centered hydroperoxy-cyclohexyl radicals (·QOOH) transiently formed in cyclohexane oxidation. The oxidation pathway leads to three nearly degenerate ·QOOH isomers, β-, γ-, and δ-QOOH, which are generated in the laboratory by H-atom abstraction from the corresponding ring sites of the cyclohexyl hydroperoxide (CHHP) precursor. The IR spectral features of jet-cooled and stabilized ·QOOH radicals are observed from 3590 to 7010 cm −1 (∼10–20 kcal mol −1 ) at energies in the vicinity of the transition state (TS) barrier leading to OH radicals that are detected by ultraviolet laser-induced fluorescence. The experimental approach affords selective detection of β-QOOH, arising from its significantly lower TS barrier to OH products compared to γ and δ isomers, which results in rapid unimolecular decay and near unity branching to OH products. The observed IR spectrum of β-QOOH includes fundamental and overtone OH stretch transitions, overtone CH stretch transitions, and combination bands involving OH or CH stretch with lower frequency modes. The assignment of β-QOOH spectral features is guided by anharmonic frequencies and intensities computed using second-order vibrational perturbation theory. The overtone OH stretch (2ν OH ) of β-QOOH is shifted only a few wavenumbers from that observed for the CHHP precursor, yet they are readily distinguished by their prompt vs slow dissociation rates to OH products.

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

Coherent spin wave excitation with radio-frequency spin–orbit torque

Spin waves, collective perturbations of magnetic moments, are both fundamental probes for magnetic physics and promising candidates for energy-efficient signal processing and computation. Traditionally, coherent propagating spin waves have been generated by radio frequency (RF) inductive Oersted fields from current-carrying electrodes. An alternative mechanism, spin–orbit torque (SOT), offers more localized excitation through interfacial spin accumulation but has been mostly limited to DC to kHz frequencies. SOT driven by RF currents, with potentially enhanced pumping efficiency and unique spin dynamics, remains largely unexplored, especially in magnetic insulators. Here, we conduct a comprehensive theoretical and computational investigation into the generation of coherent spin waves via RF-SOT in the prototypical yttrium iron garnet. We characterize the excitation of forward volume, backward volume, and surface modes in both linear and nonlinear regimes, employing single and interdigitated electrode configurations. We reveal and explain several unique and surprising features of RF-SOT compared to inductive excitation, including higher efficiency, distinct mode selectivity, and directional symmetry, a ~ $3π/4$ phase offset, reduced anharmonic distortion in the nonlinear regime, and the absence of second harmonic generation. These insights position RF-SOT as a promising new mechanism for future magnonic and spintronic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Isolating solvent–solute hydrogen bonding interactions via 2D IR solvation shell spectroscopy

The solvation shell around a solute is a fundamental feature of liquid-phase solutions, determining the behavior and properties of both the solute and the overall solution. Direct experimental measurements of the solvation shell properties are challenging due to the strong signals generated from the bulk solvent, which overwhelm the small contribution of the solvation shell. Here, we use ultrafast two dimensional infrared (2D IR) spectroscopy and intermolecular cross-peaks to isolate the IR absorption spectrum of methanol molecules in the solvation shell surrounding the solute N-methylacetamide. We demonstrate that the intermolecular coupling between the solvent and solute vibrations is indirectly mediated by a low-frequency hydrogen-bonding mode, suggesting an important mechanism for anharmonic coupling induced by hydrogen bonds. From the relative frequency shifts and cross-peak anisotropy, we find that methanol molecules surrounding N-methylacetamide form stronger and distinctly oriented hydrogen bonds than those in the bulk solvent. Here, we also compare these results with the solvent spectra of the solute N,N-dimethylacetamide to investigate how solute structural changes alter the solvation shell and the contribution of N–H hydrogen bond donation. Our results are supported by molecular dynamics simulations, which provide detailed insights into the hydrogen-bonding distributions. Through these results, we demonstrate 2D solvation shell spectroscopy to be a valuable method for investigating solvation structures and dynamics without interference from the bulk solvent.

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

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. I. Conformer- and isomer-resolved infrared spectra

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by infrared (IR) action spectroscopy under jet-cooled conditions with selective detection of hydroxyl (OH) radical products. Two distinct •QOOH isomers with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group are identified by comparison with theoretically calculated IR absorption features for multiple conformers, including C–O, ethyl, and O–O torsions of the two isomers. Here, a master-equation analysis is developed and utilized to explore the conformational cooling process and the resultant conformer population distributions of the two isomers. Most of the IR features observed, including those in the strong fundamental and first overtone OH stretch regions as well as weaker combination bands involving OH stretch with torsion or OOH bend, are in good accord with computed anharmonic frequencies for the most populated conformers of the •QOOH_Et isomer. Relatively weak IR features most evident in the fundamental and first overtone OH stretch regions are ascribed to multiple conformers of the less stable •QOOH_Me isomer, along with a weak feature that is uniquely attributed to a combination band involving asymmetric CH 2 stretch and HCH bend of •QOOH_Me.

Chemical reaction jet spectroscopy↗
DOE OSTI2020

Association of Cl with C 2 H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Significance A key issue in chemical kinetics is advancing the theoretical framework to handle reactions beyond the domain of textbook transition-state theory by including––for example––anharmonicity, barrierless transition states, transition states in series, and the effect of conformational flexibility on equilibrium constants. Other key issues are validating affordable electronic structure methods for direct dynamics and the direct calculation of high-pressure limiting rate constants, which are often obtainable experimentally only by extrapolation. The present article addresses all these issues for a prototype radical–molecule association reaction and thereby demonstrates how to combine improved theoretical methods to provide rate constants in cases where experimental data are uncertain or missing.

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

Nonadiabatic coupling of the dynamical structure to the superconductivity in YSr 2 Cu 2.75 Mo 0.25 O 7.54 and Sr 2 CuO 3.3

A crucial issue in cuprates is the extent and mechanism of the coupling of the lattice to the electrons and the superconductivity. Here we report Cu K edge extended X-ray absorption fine structure measurements elucidating the internal quantum tunneling polaron (iqtp) component of the dynamical structure in two heavily overdoped superconducting cuprate compounds, tetragonal YSr 2 Cu 2.75 Mo 0.25 O 7.54 with superconducting critical temperature, T c = 84 K and hole density p = 0.3 to 0.5 per planar Cu, and the tetragonal phase of Sr 2 CuO 3.3 with T c = 95 K and p = 0.6. In YSr 2 Cu 2.75 Mo 0.25 O 7.54 changes in the Cu-apical O two-site distribution reflect a sequential renormalization of the double-well potential of this site beginning at T c , with the energy difference between the two minima increasing by ~6 meV between T c and 52 K. Sr 2 CuO 3.3 undergoes a radically larger transformation at T c , >1-Å displacements of the apical O atoms. The principal feature of the dynamical structure underlying these transformations is the strongly anharmonic oscillation of the apical O atoms in a double-well potential that results in the observation of two distinct O sites whose Cu–O distances indicate different bonding modes and valence-charge distributions. The coupling of the superconductivity to the iqtp that originates in this nonadiabatic coupling between the electrons and lattice demonstrates an important role for the dynamical structure whereby pairing occurs even in a system where displacements of the atoms that are part of the transition are sufficiently large to alter the Fermi surface. The synchronization and dynamic coherence of the iqtps resulting from the strong interactions within a crystal would be expected to influence this process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Vibrational relaxation dynamics in layered perovskite quantum wells

Significance Halide perovskites, especially layered perovskites, offer a number of advantages to creating bright and efficient light-emitting devices and other optoelectronic applications. The organic–inorganic hybrid layered perovskite features complex lattice dynamics due to the ionic character of the crystal and the softness arising from noncovalent bonds between molecular moieties and the inorganic network. Especially, the packing geometry of the organic barriers leads to a structural deformation of inorganic octahedral, which strongly affects the properties that are crucial for device applications. In this work, we use high-resolution resonant impulsive stimulated Raman spectroscopy of a variety of ligand substitutions in layered perovskites. We find the composition of organic ligands can substantially change the dephasing rate of optical phonons and their temperature dependence, due to varying degrees of anharmonicity in the lattice and dynamic structural disorder.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Direct visualization of polaron formation in the thermoelectric SnSe

SnSe is a layered material that currently holds the record for bulk thermoelectric efficiency. The primary determinant of this high efficiency is thought to be the anomalously low thermal conductivity resulting from strong anharmonic coupling within the phonon system. Here we show that the nature of the carrier system in SnSe is also determined by strong coupling to phonons by directly visualizing polaron formation in the material. We employ ultrafast electron diffraction and diffuse scattering to track the response of phonons in both momentum and time to the photodoping of free carriers across the bandgap, observing the bimodal and anisotropic lattice distortions that drive carrier localization. Relatively large (18.7 Å), quasi-one-dimensional (1D) polarons are formed on the 300-fs timescale with smaller (4.2 Å) 3D polarons taking an order of magnitude longer (4 ps) to form. This difference appears to be a consequence of the profoundly anisotropic electron–phonon coupling in SnSe, with strong Fröhlich coupling only to zone-center polar optical phonons. These results demonstrate a high density of polarons in SnSe at optimal doping levels. Strong electron-phonon coupling is critical to the thermoelectric performance of this benchmark material and, potentially, high performance thermoelectrics more generally.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High temperature superconductivity in the candidate phases of solid hydrogen

As the simplest element in nature, unraveling the phase diagram of hydrogen is a primary task for condensed matter physics. As conjectured many decades ago, in the low-temperature and high-pressure part of the phase diagram, solid hydrogen is expected to become metallic with a high superconducting transition temperature. The metallization may occur via band gap closure in the molecular solid or via a transition to the atomic solid. Recently, a few experimental studies pushed the achievable pressures into the 400–500 GPa range. There are strong indications that at some pressure in this range metallization via either of these mechanisms occurs, although there are disagreements between experimental reports. Furthermore, there are multiple good candidate crystal phases that have emerged from recent computational and experimental studies which may be realized in upcoming experiments. Therefore, it is crucial to determine the superconducting properties of these candidate phases. In a recent study, we reported the superconducting properties of the C2/c-24 phase, which we believe to be a strong candidate for metallization via band gap closure. Here, we report the superconducting properties of the Cmca-12, Cmca-4 and I4 1 /amd-2 phases including the anharmonic effects using a Wannier function-based dense k-point and q-point sampling. Here we find that the Cmca-12 phase has a superconducting transition temperature that rises from 86 K at 400 GPa to 212 K at 500 GPa, whereas the Cmca-4 and I4 1 /amd-2 phases show a less pressure-dependent behavior with their T c in the 74–94 K and 307–343 K ranges, respectively. These properties can be used to distinguish between crystal phases in future experiments. Understanding superconductivity in pure hydrogen is also important in the study of high-T c hydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The effect of temperature oscillations on energy storage rectification in harmonic systems

Abstract Rectification, the preferential transport of a current in one direction through a system, has garnered significant attention in molecules because of its importance for controlling thermal and electronic currents at the nanoscale. Here, we report the presence of energy storage rectification effects in a molecular chain. This phenomenon is generated by subjecting a harmonic molecular chain to an oscillating temperature gradient and showing that the energy absorption rate of the system depends on the direction of the gradient. We examine how the energy storage rectification ratios in the chain are affected by the oscillating gradient, asymmetry in the chain, and the system parameters. We find that energy storage rectification can be observed in harmonic lattice structures with time-dependent temperatures and that, correspondingly, anharmonicity is not required to generate this rectification mechanism in such systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Theoretical investigation of decoherence channels in athermal phonon sensors

The creation and evolution of nonequilibrium phonons is central in applications ranging from cosmological particle searches to decoherence processes in qubits. However, the fundamental understanding of decoherence pathways for athermal phonon distributions in solid-state systems remains an open question. Using first-principles calculations, we investigate the primary decay channels of athermal phonons in two technologically relevant semiconductors—Si and GaAs. We quantify the contributions of anharmonic, isotopic, and interfacial scattering in these materials. From this, we construct a model to estimate the thermal power in a readout scheme as a function of time. We discuss the implication of our results on noise limitations in current phonon sensor designs and strategies for improving coherence in next-generation phonon sensors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗