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

Dynamical stability of a many-body Kapitza pendulum

We consider a many-body generalization of the Kapitza pendulum: the periodically-driven sine–Gordon model. We show that this interacting system is dynamically stable to periodic drives with finite frequency and amplitude. This finding is in contrast to the common belief that periodically-driven unbounded interacting systems should always tend to an absorbing infinite-temperature state. The transition to an unstable absorbing state is described by a change in the sign of the kinetic term in the Floquet Hamiltonian and controlled by the short-wavelength degrees of freedom. We investigate the stability phase diagram through an analytic high-frequency expansion, a self-consistent variational approach, and a numeric semiclassical calculation. Classical and quantum experiments are proposed to verify the validity of our results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Oscillatory reaction cross sections caused by normal mode sampling in quasiclassical trajectory calculations

The quasiclassical trajectory (QCT) method is an efficient and important tool for studying the dynamics of bimolecular reactions. In this method, the motion of the atoms is simulated classically, and the only quantum effect considered is that the initial vibrational states of reactant molecules are semiclassically quantized. A sensible expectation is that the initial ensemble of classical molecular states generated this way should be stationary, similarly to the quantum state it is supposed to represent. The most widely used method for sampling the vibrational phase space of polyatomic molecules is based on the normal mode approximation. In the present work, it is demonstrated that normal mode sampling provides a nonstationary ensemble even for a simple molecule like methane, because real potential energy surfaces are anharmonic in the reactant domain. The consequences were investigated for reaction CH{sub 4} + H → CH{sub 3} + H{sub 2} and its various isotopologs and were found to be dramatic. Reaction probabilities and cross sections obtained from QCT calculations oscillate periodically as a function of the initial distance of the colliding partners and the excitation functions are erratic. The reason is that in the nonstationary ensemble of initial states, the mean bond length of the breaking C–H bond oscillates in time with the frequency of the symmetric stretch mode. We propose a simple method, one-period averaging, in which reactivity parameters are calculated by averaging over an entire period of the mean C–H bond length oscillation, which removes the observed artifacts and provides the physically most reasonable reaction probabilities and cross sections when the initial conditions for QCT calculations are generated by normal mode sampling.

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

Hartree potential dependent exchange functional

We introduce a novel non-local ingredient for the construction of exchange density functionals: the reduced Hartree parameter, which is invariant under the uniform scaling of the density and represents the exact exchange enhancement factor for one- and two-electron systems. The reduced Hartree parameter is used together with the conventional meta-generalized gradient approximation (meta-GGA) semilocal ingredients (i.e., the electron density, its gradient, and the kinetic energy density) to construct a new generation exchange functional, termed u-meta-GGA. This u-meta-GGA functional is exact for the exchange of any one- and two-electron systems, is size-consistent and non-empirical, satisfies the uniform density scaling relation, and recovers the modified gradient expansion derived from the semiclassical atom theory. For atoms, ions, jellium spheres, and molecules, it shows a good accuracy, being often better than meta-GGA exchange functionals. Our construction validates the use of the reduced Hartree ingredient in exchange-correlation functional development, opening the way to an additional rung in the Jacob’s ladder classification of non-empirical density functionals.

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

Fluctuations in quantum mechanics and field theories from a new version of semiclassical theory. II.

This is the second paper on the semiclassical approach based on the density matrix given by the Euclidean time path integral with fixed coinciding end points. The classical path, interpolating between this point and the classical vacuum (called a “flucton”), as well as systematic one- and two-loop corrections were calculated in the first paper [M. A. Escobar-Ruiz, E. Shuryak, and A. V. Turbiner, Phys. Rev. D 93, 105039 (2016)] for a double-well potential. Here, we extend them for a number of quantum-mechanical problems, such as an anharmonic oscillator and the sine-Gordon potential. The method is based on a systematic expansion in Feynman diagrams and thus can be extended to quantum field theories (QFTs). We show that the loop expansion in quantum mechanics resembles the leading-log approximations in QFT. In this sequel, we present a complete set of results obtained using this method in a unified way. Alternatively, starting from the Schrödinger equation we derive a generalized Bloch equation whose semiclassical-like, iterative solution generates the loop expansion. We rederive the two-loop expansions for all three of the above potentials and extend them to three loops, which has not yet been done via Feynman diagrams. All results for both methods are fully consistent with each other. An asymmetric (tilted) double-well potential (nondegenerate minima) is also studied using the second method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods for the simulation of nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. As a result, the remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.

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

Bion non-perturbative contributions versus infrared renormalons in two-dimensional $\mathbb{C}$P N - 1 models

We derive the semiclassical contributions from the real and complex bions in the two-dimensional $\mathbb{C}$ P N 1 sigma model on $\mathbb{R}$× S 1 with a twisted boundary condition. The bion configurations are saddle points of the complexified Euclidean action, which can be viewed as bound states of a pair of fractional instantons with opposite topological charges. We first derive the bion solutions by solving the equation of motion in the model with a potential which simulates an interaction induced by fermions in the $\mathbb{C}$ P N - 1 quantum mechanics. The bion solutions have quasi-moduli parameters corresponding to the relative distance and phase between the constituent fractional instantons. By summing over the Kaluza-Klein modes of the quantum fluctuations around the bion backgrounds, we find that the effective action for the quasi-moduli parameters is renormalized and becomes a function of the dynamical scale (or the renormalized coupling constant). Based on the renormalized effective action, we obtain the semiclassical bion contribution in a weak coupling limit by making use of the Lefschetz thimble method. We find in the supersymmetric case that the bion contribution vanishes as expected from supersymmetry. In non-supersymmetric cases, the non-perturbative contribution has an imaginary ambiguity which is consistent with the expected infrared renormalon ambiguity. Our results explicitly demonstrate that the complex bion can explain the infrared renormalon.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Locality of correlation in density functional theory

The Hohenberg-Kohn density functional was long ago shown to reduce to the Thomas-Fermi (TF) approximation in the non-relativistic semiclassical (or large-Z) limit for all matter, i.e., the kinetic energy becomes local. Exchange also becomes local in this limit. Numerical data on the correlation energy of atoms support the conjecture that this is also true for correlation, but much less relevant to atoms. We illustrate how expansions around a large particle number are equivalent to local density approximations and their strong relevance to density functional approximations. Analyzing highly accurate atomic correlation energies, we show that E{sub C} → −A{sub C} ZlnZ + B{sub C}Z as Z → ∞, where Z is the atomic number, A{sub C} is known, and we estimate B{sub C} to be about 37 mhartree. The local density approximation yields A{sub C} exactly, but a very incorrect value for B{sub C}, showing that the local approximation is less relevant for the correlation alone. This limit is a benchmark for the non-empirical construction of density functional approximations. We conjecture that, beyond atoms, the leading correction to the local density approximation in the large-Z limit generally takes this form, but with B{sub C} a functional of the TF density for the system. The implications for the construction of approximate density functionals are discussed.

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

Semiclassical initial value representation for the quantum propagator in the Heisenberg interaction representation

One of the challenges facing on-the-fly ab initio semiclassical time evolution is the large expense needed to converge the computation. In this paper, we suggest that a significant saving in computational effort may be achieved by employing a semiclassical initial value representation (SCIVR) of the quantum propagator based on the Heisenberg interaction representation. We formulate and test numerically a modification and simplification of the previous semiclassical interaction representation of Shao and Makri [J. Chem. Phys. 113, 3681 (2000)]. The formulation is based on the wavefunction form of the semiclassical propagation instead of the operator form, and so is simpler and cheaper to implement. The semiclassical interaction representation has the advantage that the phase and prefactor vary relatively slowly as compared to the “standard” SCIVR methods. This improves its convergence properties significantly. Using a one-dimensional model system, the approximation is compared with Herman-Kluk’s frozen Gaussian and Heller’s thawed Gaussian approximations. The convergence properties of the interaction representation approach are shown to be favorable and indicate that the interaction representation is a viable way of incorporating on-the-fly force field information within a semiclassical framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2024

Monte Carlo evaluation of the semiclassical multi-step direct reaction series

We show how the multi-step direct reaction series can be evaluated using Monte Carlo methods. Nucleon-nucleon collisions occur according to the random selection of a nucleon’s attenuation factor along its classical trajectory. The particles and hole excited in a collision are selected from a local Fermi distribution. We assume that the particles continue to propagate and possibly collide again before leaving the nucleus. We assume that holes collide in place, to possibly produce other particles and holes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1979

Close-coupling calculations with an INDOX/1s static potential, semiclassical exchange, and a semi-empirical polarisation potential for electron-CO2 elastic scattering and rotational excitation

A semiempirical molecular-orbital method for modelling the effective potential for electron-molecule scattering is applied to elastic scattering and rotational excitation of CO2 at 20 eV impact energy. Agreement with experiment is reasonably good. The calculated rotationally summed integral cross section is 67.8 a sub 0 squared.

Onda, K.↗
DOE OSTI2021

Shell-structure and asymmetry effects in level densities

Level density [Formula: see text] is derived for a nuclear system with a given energy [Formula: see text], neutron [Formula: see text], and proton [Formula: see text] particle numbers, within the semiclassical extended Thomas–Fermi and periodic-orbit theory beyond the Fermi-gas saddle-point method. We obtain [Formula: see text], where [Formula: see text] is the modified Bessel function of the entropy [Formula: see text], and [Formula: see text] is related to the number of integrals of motion, except for the energy [Formula: see text]. For small shell structure contribution one obtains within the micro–macroscopic approximation (MMA) the value of [Formula: see text] for [Formula: see text]. In the opposite case of much larger shell structure contributions one finds a larger value of [Formula: see text]. The MMA level density [Formula: see text] reaches the well-known Fermi gas asymptote for large excitation energies, and the finite micro-canonical limit for low excitation energies. Fitting the MMA [Formula: see text] to experimental data on a long isotope chain for low excitation energies, due mainly to the shell effects, one obtains results for the inverse level density parameter [Formula: see text], which differs significantly from that of neutron resonances.

Physics↗
NASA NTRS1978

Experimental and finite element investigation of the buckling characteristics of a beaded skin panel for a hypersonic aircraft

As part of NASA's continuing research into hypersonics and 85 square foot hypersonic wing test section of a proposed hypersonic research airplane was laboratory tested. The project reported on in this paper has carried the hypersonic wing test structure project one step further by testing a single beaded panel to failure. The primary interest was focused upon the buckling characteristics of the panel under pure compression with boundary conditions similar to those found in a wing mounted condition. Three primary phases of analysis are included in the report. These phases include: experimental testing of the beaded panel to failure; finite element structural analysis of the beaded panel with the computer program NASTRAN; a summary of the semiclassical buckling equations for the beaded panel under purely compressive loads. Comparisons between each of the analysis methods are also included.

Siegel, W. H.↗
NASA NTRS2001

Semiclassical Models for Virtual Antiparticle Pairs, the Unit of Charge e, and the QCD Coupling alpha(sub s)

New semiclassical models of virtual antiparticle pairs are used to compute the pair lifetimes, and good agreement with the Heisenberg lifetimes from quantum field theory (QFT) is found. The modeling method applies to both the electromagnetic and color forces. Evaluation of the action integral of potential field fluctuation for each interaction potential yields approximately Planck's constant/2 for both electromagnetic and color fluctuations, in agreement with QFT. Thus each model is a quantized semiclassical representation for such virtual antiparticle pairs, to good approximation. When the results of the new models and QFT are combined, formulae for e and alpha(sub s)(q) are derived in terms of only Planck's constant and c.

Batchelor, David↗
DOE OSTI2019

The numerical investigation of colliding optical solitons as an all-optical-gate using the method of Lines

The bidirectional propagation and collision of counter-propagating optical solitons in a two-level resonant medium is analyzed. For sufficiently short duration of the pulses and sufficiently high power levels, self-induced transparency allows the pulses to propagate without losses at anomalously low velocities. Here, we use an extended semiclassical approach, called Maxwell–Bloch model. The quantization of the electromagnetic field is neglected, however, a full quantum mechanical formulation is used to describe the behavior of the matter. The numerical simulations are performed by using the method of lines. For improved accuracy, spline interpolation of the fields is introduced. Results are presented to demonstrate the use of soliton collision for all-optical switching, in particular, as OR- and XOR-gates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Nonautonomous complex wave solutions for the (2+1)-dimensional Heisenberg ferromagnetic spin chain equation with variable coefficients

This paper studies the (2+1)-dimensional Heisenberg ferromagnetic spin chain equation with variable coefficients (2D-VcHFSCE), which describes the Heisenberg ferromagnetic spin chains’s nonlinear spin dynamics with bilinear and anisotropic interactions in the semiclassical limit. Nonautonomous complex wave solutions are obtained such as similarity, bright and dark soliton solutions. Furthermore, analytic solutions are also presented based on the G{sup ′}/G-expansion method. The propagation properties of the obtained solutions are shown by some 3D- and contour plots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Open quantum systems and Schwinger-Keldysh holograms

We initiate the study of open quantum field theories using holographic methods. Specifically, we consider a quantum field theory (the system) coupled to a holographic field theory at finite temperature (the environment). We investigate the effects of integrating out the holographic environment with an aim of obtaining an effective dynamics for the resulting open quantum field theory. The influence functionals which enter this open effective action are determined by the real-time (Schwinger-Keldysh) correlation functions of the holographic thermal environment. To evaluate the latter, we exploit recent developments, wherein the semiclassical gravitational Schwinger-Keldysh saddle geometries were identified as complexified black hole spacetimes. We compute real-time correlation functions using holographic methods in these geometries, and argue that they lead to a sensible open effective quantum dynamics for the system in question, a question that hitherto had been left unanswered. In addition to shedding light on open quantum systems coupled to strongly correlated thermal environments, our results also provide a principled computation of Schwinger-Keldysh observables in gravity and holography. In particular, these influence functionals we compute capture both the dissipative physics of black hole quasi- normal modes, as well as that of the fluctuations encoded in outgoing Hawking quanta, and interactions between them. We obtain results for these observables at leading order in a low frequency and momentum expansion in general dimensions, in addition to determining explicit results for two dimensional holographic CFT environments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Full-dimensional multi-state simulation of the photodissociation of thioanisole

The photodissociation of thioanisole is very interesting because the experiments of Lim and Kim provide evidence for mode-specific effects on the product distribution. They showed that, with a specific S–CH 3 stretching mode being excited as the reagent is excited to the S 1 electronic state, there is a sharp increase in the proportion of the ground-state product to the excited-state product. In this work, we report 78 011 full-dimensional semiclassical multi-state trajectories of the photodissociation process using the coherent switching with decay of mixing dynamics method. The potential surfaces and couplings are based on electronic structure calculations that include dynamic correlation through second order perturbation theory. We report results for four sets of initial conditions, one corresponding roughly to 0–0 excitation and three corresponding to exciting one vibrational mode, to look for mode-specific effects. The simulations show no significant mode-specific effect on the product energy distributions, but they do show an effect on the distribution of minimum-energy gaps in the trajectories and on the lifetime for dissociation. In particular, excitation of the S–CH 3 stretching mode leads to trajectories passing closer to the S 1 -S 2 conical intersection and to shorter lifetimes. This provides a possible explanation of why experimental results are different for excitation of this vibration

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗