Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Semiclassical methods”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Calculating time-resolved differential absorbance spectra for ultrafast pump-probe experiments with surface hopping trajectories

We report a surface hopping approach for modeling the full time- and frequency-resolved differential absorbance spectra (beyond the inhomogenous limit) obtained in ultrafast pump-probe experiments. In our approach, we combine dynamical information obtained from ensembles of classical trajectories propagated on the ground and on the excited potential energy surfaces to directly calculate optical response functions and hence spectral lineshapes. We demonstrate that our method is exact for the model problem of two shifted harmonic potentials with identical harmonic frequencies in the absence of electronic relaxation. We then consider a model three state system with electronic relaxation and show that our method is able to capture the effects of nonadiabatic excited state dynamics on the time-dependent differential absorbance spectra. Furthermore, by comparing our spectra against those spectra calculated with either an (1) inhomogenous expression, (2) ground-state Kubo theory, or (3) excited-state Kubo theory, we show that including dynamical information from both the ground and excited potential energy surfaces significantly improves the reliability of the semiclassical approximations. As such, our surface hopping method should find immediate use in modeling the time-dependent differential abosrbance spectra of ultrafast pump-probe experiments.

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

Efficient algorithms for semiclassical instanton calculations based on discretized path integrals

Path integral instanton method is a promising way to calculate the tunneling splitting of energies for degenerated two state systems. In order to calculate the tunneling splitting, we need to take the zero temperature limit, or the limit of infinite imaginary time duration. In the method developed by Richardson and Althorpe [J. Chem. Phys. 134, 054109 (2011)], the limit is simply replaced by the sufficiently long imaginary time. In the present study, we have developed a new formula of the tunneling splitting based on the discretized path integrals to take the limit analytically. We have applied our new formula to model systems, and found that this approach can significantly reduce the computational cost and gain the numerical accuracy. We then developed the method combined with the electronic structure calculations to obtain the accurate interatomic potential on the fly. We present an application of our ab initio instanton method to the ammonia umbrella flip motion.

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

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗
DOE OSTI2017

Semiclassical approach to dynamics of interacting fermions

Understanding the behaviour of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of this. Already in equilibrium, fermions are notoriously hard to handle due to the sign problem. Out of equilibrium, an important outstanding problem is the efficient numerical simulation of the dynamics of these systems. In this work we develop a new semiclassical phase-space approach (a.k.a. the truncated Wigner approximation) for simulating the dynamics of interacting fermions in arbitrary dimensions. As fermions are essentially non-classical objects, a phase-space is constructed out of all fermionic bilinears. Classical phase-space is thus comprised of highly non-local (hidden) variables representing these bilinears, and the cost of the method is that it scales quadratic rather than linear with system size. We demonstrate the strength of the method by comparing the results to the exact quantum dynamics of fermion expansion in the Hubbard model and quantum thermalization in the Sachdev–Ye–Kitaev (SYK) model for small systems, where the semiclassics nearly perfectly reproduces correct results. We furthermore analyse fermion expansion in a larger, intractable by exact methods, 2D Hubbard model, which is directly relevant to recent cold atom experiments. - Highlights: • We propose a semiclassical method to simulate dynamics of interacting fermions. • A phase-space is constructed out of all fermionic bilinears. • We study quantum thermalization in the Sachdev–Ye–Kitaev model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

SCF treatment of charge polarization effects in intermediate-energy electron scattering calculations with applications to N2

Converged rotational close coupling calculations of differential, integral, and momentum-transfer cross sections are obtained for seven model potentials for electron-N2 scattering at an impact energy of 30 eV. The model potentials involve a static potential calculated by the INDO/1s or INDOXI/1s method, and exchange potential calculated by the semiclassical exchange approximation from the INDO/1s or INDOXI/1s unperturbed electron density and a polarization potential. The polarization potentials used include the Buckley-Burke semiempirical one and various modifications of the INDOXI and INDO SCF (self-consistent-field) adiabatic polarization potentials. Without adjustment of parameters, good agreement is obtained with the angle dependence of the experimentally measured sum of the elastic and rotational excitation differential cross sections; however, the absolute value of the calculated cross sections is about 20-30% larger than the measured values in the best case, perhaps indicating that the model potentials are too strong or should have a nonzero imaginary part.

Onda, K.↗
DOE OSTI2014

Communication: Overcoming the root search problem in complex quantum trajectory calculations

Three new developments are presented regarding the semiclassical coherent state propagator. First, we present a conceptually different derivation of Huber and Heller's method for identifying complex root trajectories and their equations of motion [D. Huber and E. J. Heller, J. Chem. Phys. 87, 5302 (1987)]. Our method proceeds directly from the time-dependent Schrödinger equation and therefore allows various generalizations of the formalism. Second, we obtain an analytic expression for the semiclassical coherent state propagator. We show that the prefactor can be expressed in a form that requires solving significantly fewer equations of motion than in alternative expressions. Third, the semiclassical coherent state propagator is used to formulate a final value representation of the time-dependent wavefunction that avoids the root search, eliminates problems with caustics and automatically includes interference. We present numerical results for the 1D Morse oscillator showing that the method may become an attractive alternative to existing semiclassical approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Quantum Bound to Chaos and the Semiclassical Limit

We discuss the quantum bound on chaos in the context of the free propagation of a particle in an arbitrarily curved surface at low temperatures. The semiclassical calculation of the Lyapunov exponent can be performed in much the same way as the corresponding one for the ‘Loschmidt echo’. The bound appears here as the impossibility to scatter a wave, by effect of the curvature, over characteristic lengths smaller than the deBroglie wavelength.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Experimental investigation of terahertz quantum cascade laser with variable barrier heights

We report an experimental study of terahertz quantum cascade lasers with variable barrier heights based on the Al{sub x}Ga{sub 1–x}As/GaAs material system. Two new designs are developed based on semiclassical ensemble Monte Carlo simulations using state-of-the-art Al{sub 0.15}Ga{sub 0.85}As/GaAs three-quantum-well resonant phonon depopulation active region design as a reference. The new designs achieved maximum lasing temperatures of 188 K and 172 K, as compared to the maximum lasing temperature of 191 K for the reference structure. These results demonstrate that terahertz quantum cascade laser designs with variable barrier heights provide a viable alternative to the traditional active region designs with fixed barrier composition. Additional design space offered by using variable barriers may lead to future improvements in the terahertz quantum cascade laser performance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

Calculation of rotational transition probabilities in molecular collisions - Application to N2 + N2

A computational method is proposed to obtain rotational transition probabilities in collisions between two diatomic molecules. The potential method of Rabitz and an exponential approximation are used to solve the semiclassical coupled equations without invoking any perturbational technique. The collision trajectory is determined in the classical modified-wave-number approximation. The method can treat systems involving strong interactions and provide probabilities for transitions even with a multiquantum jump. A simultaneous transition in the rotational states of both molecules, i.e., the rotational-rotational energy transfer, is taken into account. An application to the system N2 + N2 is presented.

Itikawa, Y.↗
DOE OSTI2022

Effect of nuclear charge on laser-induced fusion enhancement in advanced fusion fuels

Here, based on the preliminary work done in a previous paper [Phys. Rev. C 105, 054001 (2022)], we investigate the effects of laser-induced fusion rate enhancements of different fusion fuels. In the aforementioned work, which considered $_{1}^{2}$H-$_{1}^{3}$H (D-T) and DHe 3 fusion, it was observed that a larger product of charge numbers of the fusion reactants leads to an increased laser-induced enhancement to the fusion cross section for static external electric fields. We investigate whether this trend persists for DT, DHe 3 , and pB 11 fusion for dynamical electric fields, using the semiclassical approaches of the Wentzel-Kramers-Brillouin method and the imaginary-time method, as well as the Kramers- Henneberger method, and the Volkoff-state approximation. We find that the fusion cross section of pB 11 indeed exhibits the largest relative enhancement for all laser parameters considered and may even surpass the cross section of DT and DHe 3 fusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Accelerating ring-polymer molecular dynamics with parallel-replica dynamics

Nuclear quantum effects are important for systems containing light elements, and the effects are more prominent in the low temperature regime where the dynamics also becomes sluggish. We show that parallel replica (ParRep) dynamics, an accelerated molecular dynamics approach for infrequent-event systems, can be effectively combined with ring-polymer molecular dynamics, a semiclassical trajectory approach that gives a good approximation to zero-point and tunneling effects in activated escape processes. The resulting RP-ParRep method is a powerful tool for reaching long time scales in complex infrequent-event systems where quantum dynamics are important. Two illustrative examples, symmetric Eckart barrier crossing and interstitial helium diffusion in Fe and Fe–Cr alloy, are presented to demonstrate the accuracy and long-time scale capability of this approach.

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

Low-temperature phonoemissive tunneling rates in single molecule magnets

Tunneling between the two lowest energy levels of single molecule magnets with Ising type anisotropy, accompanied by the emission or absorption of phonons, is considered. Quantitatively accurate calculations of the rates for such tunneling are performed for a model Hamiltonian especially relevant to the best studied example, Fe{sub 8}. Two different methods are used: high-order perturbation theory in the spin–phonon interaction and the non-Ising-symmetric parts of the spin Hamiltonian, and a novel semiclassical approach based on spin-coherent-state-path-integral instantons. The methods are found to be in good quantitative agreement with other, and consistent with previous approaches to the problem. The implications of these results for magnetization of molecular solids of these molecules are discussed briefly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

How to remove the spurious resonances from ring polymer molecular dynamics

Two of the most successful methods that are presently available for simulating the quantum dynamics of condensed phase systems are centroid molecular dynamics (CMD) and ring polymer molecular dynamics (RPMD). Despite their conceptual differences, practical implementations of these methods differ in just two respects: the choice of the Parrinello-Rahman mass matrix and whether or not a thermostat is applied to the internal modes of the ring polymer during the dynamics. Here, we explore a method which is halfway between the two approximations: we keep the path integral bead masses equal to the physical particle masses but attach a Langevin thermostat to the internal modes of the ring polymer during the dynamics. We justify this by showing analytically that the inclusion of an internal mode thermostat does not affect any of the established features of RPMD: thermostatted RPMD is equally valid with respect to everything that has actually been proven about the method as RPMD itself. In particular, because of the choice of bead masses, the resulting method is still optimum in the short-time limit, and the transition state approximation to its reaction rate theory remains closely related to the semiclassical instanton approximation in the deep quantum tunneling regime. In effect, there is a continuous family of methods with these properties, parameterised by the strength of the Langevin friction. Here, we explore numerically how the approximation to quantum dynamics depends on this friction, with a particular emphasis on vibrational spectroscopy. We find that a broad range of frictions approaching optimal damping give similar results, and that these results are immune to both the resonance problem of RPMD and the curvature problem of CMD.

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

Realistic magnetic thermodynamics by local quantization of a semiclassical Heisenberg model

Abstract Classical Monte Carlo simulation of the Heisenberg model poorly describes many thermodynamic phenomena due to its neglect of the quantum nature of spins. Alternatively, we discuss how to semiclassically approach the quantum problem and demonstrate a simple method for introducing a locally approximate form of spin quantization. While the procedure underestimates magnetic short-range order, our results suggest a simple correction for recovering realistic spin–spin correlations above the critical temperature. Moreover, ensemble fluctuations are found to provide reasonably accurate thermodynamics, largely reproducing quantum mechanically calculated heat capacities and experimental magnetometry for ferromagnetic Fe and antiferromagnetic RbMnF 3 . Extensions of the method are proposed to address remaining inaccuracies.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The microscopic Einstein-de Haas effect

The Einstein-de Haas (EdH) effect, where the spin angular momentum of electrons is transferred to the mechanical angular momentum of atoms, was established experimentally in 1915. While a semiclassical explanation of the effect exists, modern electronic structure methods have not yet been applied to model the phenomenon. In this paper, we investigate its microscopic origins by means of a noncollinear tight-binding model of an O2 dimer, which includes the effects of spin-orbit coupling, coupling to an external magnetic field, and vector Stoner exchange. By varying an external magnetic field in the presence of spin-orbit coupling, a torque can be generated on the dimer, validating the presence of the EdH effect. The avoided energy level crossings and the rate of change of magnetic field determine the evolution of the spin. We also find that the torque exerted on the nuclei by the electrons in a time-varying B field is not only due to the EdH effect. In conclusion, the other contributions arise from field-induced changes in the electronic orbital angular momentum and from the direct action of the Faraday electric field associated with the time-varying magnetic field.

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

Self-Broadening and Carbon-Dioxide Broadening of Lines of the H{sub 2}S Molecule

Carbon-dioxide-broadening coefficients and self-broadening coefficients of lines of the main isotopic modification of Н{sub 2}S are estimated on the basis of literature data. The J′-dependences of the above line-profile parameters of the hydrogen-sulfide molecule are examined. In the case of CO{sub 2} broadening, the half-widths of lines are calculated by a semiempirical method based on a parametric modification of the impact semiclassical model; the model parameters were determined from the fit to experimental data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

On the derivation of the semiclassical approximation to the quantum propagator

In order to rigorously derive the amplitude factor of the semiclassical approximation to the quantum propagator, we extend an existing method originally devised to evaluate Gaussian path-integral expressions. Using a result which relates the determinant of symmetric block-tridiagonal matrices to the determinants of their blocks, two difference equations are obtained. The first one allows to establish the connection of the amplitude factor to Jacobi’s accessory equations in the continuous-time limit, while the second one leads to an additional factor which, however, contributes to the final result only in exceptional cases. In order to demonstrate the wide applicability of these difference equations, we treat explicitly the case where the time-sliced Lagrangian is written in generalized coordinates, for which a general derivation has so far been unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Ring-polymer instanton theory of electron transfer in the nonadiabatic limit

We take the golden-rule instanton method derived in the previous paper [J. O. Richardson, R. Bauer, and M. Thoss, J. Chem. Phys. 143, 134115 (2015)] and reformulate it using a ring-polymer instanton approach. This gives equations which can be used to compute the rates of electron-transfer reactions in the nonadiabatic (golden-rule) limit numerically within a semiclassical approximation. The multidimensional ring-polymer instanton trajectories are obtained efficiently by minimization of the action. In this form, comparison with Wolynes’ quantum instanton method [P. G. Wolynes, J. Chem. Phys. 87, 6559 (1987)] is possible and we show that our semiclassical approach is the steepest-descent limit of this method. We discuss advantages and disadvantages of both methods and give examples of where the new approach is more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗