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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

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

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

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

On the temperature and density dependence of dislocation drag from phonon wind

At extreme strain rates, where fast moving dislocations govern plastic deformation, anharmonic phonon scattering imparts a drag force on the dislocations. In this paper, we present calculations of the dislocation drag coefficients of aluminum and copper as functions of temperature and density. We discuss the sensitivity of the drag coefficients to changes in the third-order elastic constants with temperature and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2022

Semi-empirical quantum optics for mid-infrared molecular nanophotonics

Nanoscale infrared (IR) resonators with sub-diffraction limited mode volumes and open geometries have emerged as new platforms for implementing cavity quantum electrodynamics at room temperature. The use of IR nanoantennas and tip nanoprobes to study strong light–matter coupling of molecular vibrations with the vacuum field can be exploited for IR quantum control with nanometer spatial and femtosecond temporal resolution. In order to advance the development of molecule-based quantum nanophotonics in the mid-IR, we propose a generally applicable semi-empirical methodology based on quantum optics to describe light–matter interaction in systems driven by mid-IR femtosecond laser pulses. The theory is shown to reproduce recent experiments on the acceleration of the vibrational relaxation rate in infrared nanostructures. It also provides physical insights on the implementation of coherent phase rotations of the near-field using broadband nanotips. We then apply the quantum framework to develop general tip-design rules for the experimental manipulation of vibrational strong coupling and Fano interference effects in open infrared resonators. We finally propose the possibility of transferring the natural anharmonicity of molecular vibrational levels to the resonator near-field in the weak coupling regime to implement intensity-dependent phase shifts of the coupled system response with strong pulses and develop a vibrational chirping model to understand the effect. The semi-empirical quantum theory is equivalent to first-principles techniques based on Maxwell’s equations, but its lower computational cost suggests its use as a rapid design tool for the development of strongly coupled infrared nanophotonic hardware for applications ranging from quantum control of materials to quantum information processing.

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

Infrared spectroscopic signature of a hydroperoxyalkyl radical (•QOOH)

Infrared (IR) action spectroscopy is utilized to characterize a prototypical carbon-centered hydroperoxyalkyl radical (•QOOH) transiently formed in the oxidation of volatile organic compounds. The •QOOH radical formed in isobutane oxidation, 2-hydroperoxy-2-methylprop-1-yl, •CH 2 (CH 3 ) 2 COOH, is generated in the laboratory by H-atom abstraction from tert-butyl hydroperoxide (TBHP). IR spectral features of jet-cooled and stabilized •QOOH radicals are observed from 2950 to 7050 cm –1 at energies that lie below and above the transition state barrier leading to OH radical and cyclic ether products. The observed •QOOH features include overtone OH and CH stretch transitions, combination bands involving OH or CH stretch and a lower frequency mode, and fundamental OH and CH stretch transitions. Most features arise from a single vibrational transition with band contours well simulated at a rotational temperature of 10 K. In each case, the OH products resulting from unimolecular decay of vibrationally activated •QOOH are detected by UV laser-induced fluorescence. Assignments of observed •QOOH IR transitions are guided by anharmonic frequencies computed using second order vibrational perturbation theory, a 2 + 1 model that focuses on the coupling of the OH stretch with two low-frequency torsions, as well as recently predicted statistical •QOOH unimolecular decay rates that include heavy-atom tunneling. Most of the observed vibrational transitions of •QOOH are readily distinguished from those of the TBHP precursor. Furthermore, the distinctive IR transitions of •QOOH, including the strong fundamental OH stretch, provide a general means for detection of •QOOH under controlled laboratory and real-world conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗