Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 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↗
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 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

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↗
DOE OSTI2025

Systematic study of the validity of the eikonal model including uncertainties

Nuclear reactions at intermediate beam energies are often interpreted using the eikonal model. In the analysis of complex reaction probes, where few-body reaction methods are needed, the eikonal method may be used as an efficient way for describing the fragment-target reaction process. In this work, we perform a systematic study to test the validity of the eikonal approximation for nucleon-nucleus reactions. We also quantify uncertainties due to the nucleon optical potential on reaction observables. We inspect the validity of the eikonal model and its semiclassical correction by comparing it to exact solutions (obtained from solving the optical-model equation with a finite-differences method) for a wide range of reactions. We also study the effect of relativistic corrections, both kinematic and dynamic, by effectively incorporating the relativistic effects at intermediate energies. The uncertainties from a Bayesian global optical potential (KDUQ) are propagated to the observables of interest. Our study includes neutron and proton reactions on 27 Al , 40 Ca , 90 Zr , and 208 Pb , for a wide range of energies 𝐸 lab = 0–400 MeV. We calculate neutron-total cross sections (elastic and reactions) as well as proton-absorption cross sections as a function of beam energy, using the eikonal model, the eikonal model with a semiclassical correction, and the exact solution. Here, we also compute angular distributions for the methods above. Our results show that for the proton-absorption cross section, the eikonal model can be used down to around 60 MeV and the semiclassical correction extends its use to 30 MeV. However, the validity of the eikonal model for the neutron-total cross section only goes down to ≈120 MeV, a range extended to ≈ 50 MeV when using the semiclassical correction. We find the semiclassical correction to the eikonal model to be less effective in describing the angular distributions. The 1⁢𝜎 uncertainty intervals on the observables we studied is less than 5% for most of the energies considered, but increases rapidly for higher energies, namely energies outside the range of KDUQ (𝐸 lab > 200MeV).

Cluster models↗
DOE OSTI2017

Fluctuations of the gluon distribution from the small- x effective action

The computation of observables in high-energy QCD involves an average over stochastic semiclassical small-x gluon fields. The weight of various configurations is determined by the effective action. We introduce a method to study fluctuations of observables, functionals of the small-x fields, which does not explicitly involve dipoles. We integrate out those fluctuations of the semiclassical gluon field under which a given observable is invariant. Thereby we obtain the effective potential for that observable describing its fluctuations about the average. Here, we determine explicitly the effective potential for the covariant gauge gluon distribution both for the McLerran-Venugopalan (MV) model and for a (nonlocal) Gaussian approximation for the small-x effective action. This provides insight into the correlation of fluctuations of the number of hard gluons versus their typical transverse momentum. We find that the spectral shape of the fluctuations of the gluon distribution is fundamentally different in the MV model, where there is a pileup of gluons near the saturation scale, versus the solution of the small-x JIMWLK renormalization group, which generates essentially scale-invariant fluctuations above the absorptive boundary set by the saturation scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system’s complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Chemistry↗
NASA NTRS1979

Solution of Boltzmann equation for highly nonequilibrium diatomic gases rotational translational energy relaxation

The direct simulation Monte Carlo method is applied to solve the Boltzmann equation for collisions between internally excited diatomic gases in highly nonequilibrium states. The semiclassical transition probability is incorporated in the simulation for energy exchange between rotational and translational energy. The results provide details on the fundamental mechanisms of gas kinetics where analytical methods are impractical. The validity of the local Maxwellian assumption and relaxation time, rotational-translational energy transition, and a velocity analysis of the inelastic collision are discussed in detail.

Yoshikawa, K. K.↗
DOE OSTI2022

Microscopic-macroscopic level densities for low excitation energies

Level density ρ(E,Q) is derived within the micro-macroscopic approximation (MMA) for a system of strongly interacting Fermi particles with the energy E and additional integrals of motion Q , in line with several topics of the universal and fruitful activity of A. S. Davydov. Within the extended Thomas Fermi and semiclassical periodic orbit theory beyond the Fermi-gas saddle-point method, we obtain ρ ∝ I ν (S)/S ν , where I ν (S) is the modified Bessel function of the entropy S . For small shell-structure contribution, one finds ν = κ/2 + 1, where κ is the number of additional integrals of motion. This integer number is a dimension of Q, Q = { N, Z , …} for the case of two-component atomic nuclei, where N and Z are the numbers of neutrons and protons, respectively. For much larger shell structure contributions, one obtains ν = κ /2 + 2. The MMA level density ρ reaches the well-known Fermi gas asymptote for large excitation energies and the finite micro-canonical combinatoric limit for low excitation energies. Further, the additional integrals of motion can also be the projection of the angular momentum of a nuclear system for nuclear rotations of deformed nuclei, number of excitons for collective dynamics, and so on. Fitting the MMA total level density ρ( E , Q) for a set of the integrals of motion Q = { N, Z }, to experimental data on a long nuclear isotope chain for low excitation energies, one obtains the results for the inverse level-density parameter K , which differs significantly from those of neutron resonances due to shell, isotopic asymmetry, and pairing effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Direct and indirect spin current generation and spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers

Spin-orbit torques in ferromagnet/nonmagnet/ferromagnet trilayers are studied using a combination of symmetry analysis, circuit theory, semiclassical simulations, and first-principles calculations using the nonequilibrium Green's function method with supercell disorder averaging. Here, we focus on unconventional processes involving the interplay between the two ferromagnetic layers, which are classified into direct and indirect mechanisms. The direct mechanism involves spin current generation by one ferromagnetic layer and its subsequent absorption by the other. In the indirect mechanism, the in-plane spin-polarized current from one ferromagnetic layer “leaks” into the other layer, where it is converted into an out-of-plane spin current and reabsorbed by the original layer. The direct mechanism results in a predominantly dampinglike torque, which damps the magnetization towards a certain direction 𝐬 𝑑 . The indirect mechanism results in a predominantly fieldlike torque with respect to a generally different direction 𝐬 𝑓 . Similarly to the current-in-plane giant magnetoresistance, the indirect mechanism is only active if the thickness of the nonmagnetic spacer is smaller than or comparable to the mean free path. Numerical calculations for a semiclassical model based on the Boltzmann equation confirm the presence of both direct and indirect mechanisms of spin current generation. First-principles calculations reveal sizable unconventional spin-orbit torques in Co/Cu/Co, Py/Cu/Py, and Co/Pt/Co trilayers and provide strong evidence of indirect spin current generation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Shells, orbit bifurcations, and symmetry restorations in Fermi systems

The periodic-orbit theory based on the improved stationary-phase method within the phase-space path integral approach is presented for the semiclassical description of the nuclear shell structure, concerning themain topics of the fruitful activity ofV.G. Soloviev. We apply this theory to study bifurcations and symmetry breaking phenomena in a radial power-law potential which is close to the realistic Woods–Saxon one up to about the Fermi energy. Using the realistic parametrization of nuclear shapes we explain the origin of the double-humped fission barrier and the asymmetry in the fission isomer shapes by the bifurcations of periodic orbits. The semiclassical origin of the oblate–prolate shape asymmetry and tetrahedral shapes is also suggested within the improved periodic-orbit approach. The enhancement of shell structures at some surface diffuseness and deformation parameters of such shapes are explained by existence of the simple local bifurcations and new non-local bridge-orbit bifurcations in integrable and partially integrable Fermi-systems. We obtained good agreement between the semiclassical and quantum shell-structure components of the level density and energy for several surface diffuseness and deformation parameters of the potentials, including their symmetry breaking and bifurcation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗