Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “semiclassical”

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

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

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

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Interfacial spin–orbit torques

Spin–orbit torques offer a promising mechanism for electrically controlling magnetization dynamics in nanoscale heterostructures. While spin–orbit torques occur predominately at interfaces, the physical mechanisms underlying these torques can originate in both the bulk layers and at interfaces. Classifying spin–orbit torques based on the region that they originate in provides clues as to how to optimize the effect. While most bulk spin–orbit torque contributions are well studied, many of the interfacial contributions allowed by symmetry have yet to be fully explored theoretically and experimentally. To facilitate progress, we review interfacial spin–orbit torques from a semiclassical viewpoint and relate these contributions to recent experimental results. In this study, we show the relationship between different interface transport parameters within the same model. For charges and spins flowing perpendicular to the interface, interfacial spin–orbit coupling both modifies the mixing conductance of the magnetoelectronic circuit theory and gives rise to spin memory loss. For in-plane electric fields, interfacial spin–orbit coupling gives rise to torques described by spin–orbit filtering, spin swapping, and precession. In addition, these same interfacial processes generate spin currents that flow into the non-magnetic layer. For in-plane electric fields in trilayer structures, the spin currents generated at the interface between one ferromagnetic layer and the non-magnetic spacer layer can propagate through the non-magnetic layer to produce novel torques on the other ferromagnetic layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High electron mobility of Al x Ga1− x N evaluated by unfolding the DFT band structure

We calculate the alloy-disorder-limited electron mobility of AlxGa1−xN from first principles. AlxGa1−xN is a technologically important ultra-wide-bandgap alloy with promise in light emitting diodes and high-power transistors. Alloying introduces statistical disorder, which causes electrons to scatter between different crystal-momentum states, leading to a reduction in mobility for intermediate alloy compositions. The corresponding lifetime, which appears as an energy broadening in the band structure, can be evaluated by unfolding the band structure from the supercell basis to the primitive-cell basis. We fit the first-principles band broadening with a model scattering potential and evaluate the low-field electron mobility using the semiclassical Boltzmann transport equation in the relaxation-time approximation. Our calculated mobility is in agreement with experimental values. We also find the lowest alloy-scattering electron mobility (total electron mobility) across the entire composition range to be 186 cm2/V s (136 cm2/V s), which is comparable to the highest electron mobility predicted in the competitor system, β-(AlxGa1−x)2O3. Our results elucidate the intrinsic limits imposed by alloy disorder on electron transport in AlxGa1−xN.

Pant, Nick (ORCID:0000000254506045)↗
DOE OSTI2021

Direct coherent switching with decay of mixing for intersystem crossing dynamics of thioformaldehyde: The effect of decoherence

We evaluate the effect of electronic decoherence on intersystem crossing in the photodynamics of thioformaldehyde. First, we show that the state-averaged complete-active-space self-consistent field electronic structure calculations with a properly chosen active space of 12 active electrons in 10 active orbitals can predict the potential energy surfaces and the singlet–triplet spin–orbit couplings quite well for CH 2 S, and we use this method for direct dynamics by coherent switching with decay of mixing (CSDM). We obtain similar dynamical results with CSDM or by adding energy-based decoherence to trajectory surface hopping, with the population of triplet states tending to a small steady-state value over 500 fs. Without decoherence, the state populations calculated by the conventional trajectory surface hopping method or the semiclassical Ehrenfest method gradually increase. This difference shows that decoherence changes the nature of the results not just quantitatively but qualitatively.

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

CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)—an opensource software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi’s golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.

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

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al . J. Am. Chem. Soc. 142, 10459–10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. Here, in this work, we explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck–Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

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

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

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

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗
DOE OSTI2025

Can increasing the size and flexibility of a molecule reduce decoherence and prolong charge migration?

Coherent superposition of electronic states, created by ionizing a molecule, can initiate ultrafast dynamics of the electron density. Correlation between nuclear and electron motions, however, typically dissipates the electronic coherence in only a few femtoseconds, especially in larger and more flexible molecules. We, therefore, use ab initio semiclassical dynamics to study decoherence in a sequence of analogous organic molecules of increasing size and find, surprisingly, that extending the carbon skeleton in propynal analogs slows down decoherence and prolongs charge migration. To elucidate this observation, we decompose the overall decoherence into contributions from individual vibrational modes and show that: 1) The initial decay of electronic coherence is caused by high- and intermediate-frequency vibrations via momentum separation of nuclear wavepackets evolving on different electronic surfaces. 2) At later times, the coherence disappears completely due to the increasing position separation in the low-frequency modes. 3) In agreement with another study, we observe that only normal modes that preserve the symmetry of the molecule induce decoherence. All together, we justify the enhanced charge migration by a combination of increased hole-mixing and the disappearance of decoherence contributions from specific vibrational modes—CO stretching in butynal and various H rockings in pentynal.

Science & Technology - Other Topics↗
DOE OSTI2020

Pattern of perturbations from a coherent quantum inflationary horizon

It is proposed that if quantum states of space-time are coherent on null surfaces, holographic Planck-scale fluctuations of inflationary horizons dominate the formation of primordial scalar curvature perturbations. It is shown that the reduction of quantum states on nearly-spherical emergent horizon surfaces around each observer creates a distinctive pattern whose correlations in the angular domain differ from the standard quantum theory of inflation. Causal constraints are used in a semiclassical model to formulate candidate directional symmetries. It is suggested that this hypothesis could provide a physical explanation for several well known anomalies measured in CMB anisotropy. New exact symmetries are predicted, such as a vanishing temperature correlation function at 90 degrees angular separation, that can be tested with current data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Attosecond synchronization of extreme ultraviolet high harmonics from crystals

The interaction of strong near-infrared (NIR) laser pulses with wide-bandgap dielectrics produces high harmonics in the extreme ultraviolet (XUV) wavelength range. These observations have opened up the possibility of attosecond metrology in solids, which would benefit from a precise measurement of the emission times of individual harmonics with respect to the NIR laser field. In this work, we show that, when high-harmonics are detected from the input surface of a magnesium oxide crystal, a bichromatic probing of the XUV emission shows a clear synchronization largely consistent with a semiclassical model of electron–hole recollisions in bulk solids. On the other hand, the bichromatic spectrogram of harmonics originating from the exit surface of the 200 μm-thick crystal is strongly modified, indicating the influence of laser field distortions during propagation. Our tracking of sub-cycle electron and hole re-collisions at XUV energies is relevant to the development of solid-state sources of attosecond pulses.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Laser-assisted radiative recombination in a cold hydrogen plasma

Abstract We study the process of laser-assisted radiative recombination of an electron with a proton in a cold hydrogen plasma employing the semiclassical Kramers’ approach which involves calculation of classical trajectories in combined laser and Coulomb fields and the use of the correspondence principle. Due to the Coulomb focusing effect, recombination is the most effective when the initial electron momentum is parallel to the laser polarization. Orders of magnitude enhancement of the cross section, as compared to the laser-free case, is observed in this case. With increasing angle between the electron momentum and polarization, the recombination cross section drops. However, even after averaging over Maxwellian velocity distribution we obtain a substantial enhancement of the recombination rate constant, as compared to the zero-field case. For the field intensities in the range 30–350 MW cm −2 , the enhancement occurs in the region of the radiation wavelength from 5 to 20 µ m and for the plasma temperature from 20 to 300 K.

Fabrikant, I. I. (ORCID:0000000293849454)↗
DOE OSTI2024

Effective action approach for preheating

Abstract We present a semiclassical non-perturbative approach for calculating the preheating process at the end of inflation. Our method involves integrating out the decayed particles within the path integral framework and subsequently determining world-line instanton solutions in the effective action. This enables us to obtain the effective action of the inflaton, with its imaginary part linked to the phenomenon of particle creation driven by coherent inflaton field oscillations. Additionally, we utilize the Bogoliubov transformation to investigate the evolution of particle density within the medium after multiple inflaton oscillations. We apply our approach to various final state particles, including scalar fields, tachyonic fields, and gauge fields. The non-perturbative approach provides analytical results for preheating that are in accord with previous methods.

Astronomy & Astrophysics↗
DOE OSTI2025

Next-generation tunnel FETs: exploring material perspectives and areal tunneling configurations

The end of Dennard scaling, which facilitated proportional increases in computing power without added energy costs until the mid-2000s, has underscored the urgent need for innovative semiconductor devices that can enhance energy efficiency. Tunnel field-effect transistors (TFETs) have emerged as promising candidates to surpass the energy efficiency of conventional metal oxide semiconductor field-effect transistors (MOSFETs). Unlike MOSFETs, which rely on thermionic emission to overcome the source-channel potential barrier, TFETs operate through quantum tunneling, potentially enabling sub-60 mV dec −1 subthreshold swing (SS) for low-voltage operation. However, lateral TFETs have faced challenges in achieving adequate on-state current (I ON ) and a broad SS operation window, limiting their practical utility. This review article advocates for areal TFETs, which utilize face-to-face tunnel junctions that ideally offer step-function current turn-on characteristics and allow I ON to scale with device area rather than width. We highlight recent advancements in integrating 2D materials into tunneling structures, which could facilitate efficient band-to-band tunneling through atomically thin layers, while addressing challenges of gate field screening. We then discuss the nearer-term prospects of epitaxial areal TFETs comprising III–V compound semiconductors and group-IV semiconductors based on recent experimental progress. The review examines both quantum mechanical and semiclassical modeling approaches for TFETs, including techniques to reduce the computational complexity. The article delves into ongoing challenges in material synthesis, interface engineering, device fabrication, and integration pathways, concluding with recommendations for future research directions to overcome the fundamental power density limitations of conventional transistor technology.

2D materials↗
DOE OSTI2026

Single pion production and pion propagation in A chilles

We extend the applicability of A chilles (A CHIcagoLand Lepton Event Simulator) by incorporating the single pion production mechanism in a fully exclusive fashion. The electroweak interaction vertex is modeled by combining the state-of-the-art dynamical coupled-channels approach with realistic hole spectral functions, which account for correlations in both the initial target state and the residual spectator system. Final-state interactions are treated using a semiclassical intranuclear cascade that leverages nuclear configurations to determine the correlated spatial distribution of protons and neutrons. The meson-baryon scattering amplitudes used in the cascade are computed within the dynamical coupled-channels framework, consistent with the electroweak vertex. To model pion absorption, we employ the optical potential approach of Oset and Salcedo [Nucl. Phys. A484, 557 (1988).]. As an alternative approach, we explicitly model the production and propagation of resonances which mediate pion-nucleon scattering and pion absorption. We validate our approach against pion-nucleon and pion-nucleus scattering data, and present comparisons with electron- and neutrino-nucleus measurements from e⁢4⁢𝜈, T2K, MINER⁢𝜈⁢A, and MicroBooNE.

Isaacson, Joshua [Michigan State Univ., East Lansi↗
DOE OSTI2025

Intrinsic Thermal Hall Effect in Mott Insulators

Here, in light of recent experimental data indicating a substantial thermal Hall effect in square lattice antiferromagnetic Mott insulators, we investigate whether a simple Mott insulator can sustain a finite thermal Hall effect. We verify that the answer is “no” if one performs calculations within a spin-only low-energy effective spin model with noninteracting magnons. However, by performing determinant quantum Monte Carlo simulations, we show the single-band 𝑡−𝑡′−𝑈 Hubbard model coupled to an orbital magnetic field does support a finite thermal Hall effect when 𝑡′ ≠0 and 𝐵 ≠0 in the Mott insulating phase. We argue that the (carrier agnostic) necessary conditions for observing a finite thermal Hall effect are time-reversal and particle-hole symmetry breaking. By considering magnon-magnon scattering using a semiclassical Boltzmann analysis, we illustrate a physical mechanism by which finite transverse thermal conductivity may arise, consistent with our symmetry argument and numerical results. Our results indicate that square and triangular lattices with SU(2) symmetry can support a finite thermal Hall effect and call for a critical reexamination of thermal Hall effect data in insulating magnets, as the magnon contribution should not be excluded a priori .

Mott insulators↗