Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PATH INTEGRALS”

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

Chirped laser dispersion spectroscopy for spectroscopic chemical sensing with simultaneous range detection

Spectroscopic chemical detection requires knowledge or determination of an optical path for accurate quantification of path-integrated concentration of species. Continuous- wave-laser-based spectroscopic systems operating in an open integrated-path remote sensing configuration are usually not equipped for optical path determination. Here we demonstrate a measurement technique capable of simultaneous spectroscopic chemical quantification and range finding. Here, the range-finding functionality is implemented with chirped laser dispersion spectroscopy. The methodology is potentially useful for remote chemical sensing in a hard-target LIDAR configuration and for automatic calibration of gas cells with unknown or varying lengths.

03 NATURAL GAS↗
DOE OSTI2014

Solvation of molecules in superfluid helium enhances the “interaction induced localization” effect

Atomic nuclei become delocalized at low temperatures as a result of quantum effects, whereas they are point-like in the high temperature (classical) limit. For non-interacting nuclei, the delocalization upon lowering the temperature is quantitatively described in terms of the thermal de Broglie wavelength of free particles. Clearly, light non-interacting nuclei – the proton being a prominent one – are much more delocalized at low temperatures compared to heavy nuclei, such as non-interacting oxygen having water in mind. However, strong interactions due to chemical bonding in conjunction with ultra-low temperatures characteristic to superfluid helium nanodroplets change this common picture substantially for nuclei in molecules or clusters. It turns out that protons shared in hydrogen bonds undergo an extreme “interaction induced localization” at temperatures on the order of 1 K, which compresses the protonic spatial distributions to the size of the much heavier donor or acceptor atoms, such as O or Cl nuclei, corresponding to about 0.1% of the volume occupied by a non-interacting proton at the same temperature. Moreover, applying our recently developed hybrid ab initio path integral molecular dynamics/bosonic path integral Monte Carlo quantum simulation technique to a HCl/water cluster, HCl(H{sub 2}O){sub 4}, we find that helium solvation has a significant additional localizing effect of up to about 30% in volume. In particular, the solvent-induced excess localization is the stronger the lesser the given nucleus is already localized in the gas phase reference situation.

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

Ab initio molecular dynamics with nuclear quantum effects at classical cost: Ring polymer contraction for density functional theory

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ring polymer contraction scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, this approach is as fast as a typical classical ab initio molecular dynamics simulation and 35× faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

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

Can the ring polymer molecular dynamics method be interpreted as real time quantum dynamics?

The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The present quantum dynamics formulation also serves as the basis for developing new quantum dynamical methods that utilize the cyclic nature of the imaginary time path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Summing over trajectories of stochastic dynamics with multiplicative noise

We demonstrate that previous path integral formulations for the general stochastic interpretation generate incomplete results exemplified by the geometric Brownian motion. We thus develop a novel path integral formulation for the overdamped Langevin equation with multiplicative noise. The present path integral leads to the corresponding Fokker-Planck equation, and naturally generates a normalized transition probability in examples. Our result solves the inconsistency of the previous path integral formulations for the general stochastic interpretation, and can have wide applications in chemical and physical stochastic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Spin structures and baby universes

We extend a 2d topological model of the gravitational path integral to include sums over spin structure, corresponding to Neveu-Schwarz (NS) or Ramond (R) boundary conditions for fermions. This path integral corresponds to a correlator of boundary creation operators on a non-trivial baby universe Hilbert space, and vanishes when the number of R boundaries is odd. This vanishing implies a non-factorization of the correlator, which necessitates a dual interpretation of the bulk path integral in terms of a product of partition functions (associated to NS boundaries) and Witten indices (associated to R boundaries), averaged over an ensemble of theories with varying Hilbert space dimension and different numbers of bosonic and fermionic states. We also consider a model with End-of-the-World (EOW) branes, for which the dual ensemble then includes a sum over randomly chosen fermionic and bosonic states. We propose two modifications of the bulk path integral which restore an interpretation in a single dual theory: (i) a geometric prescription where we add extra boundaries with a sum over their spin structures, and (ii) an algebraic prescription involving “spacetime D-branes”. We extend our ideas to Jackiw-Teitelboim gravity, and propose a dual description of a single unitary theory with spin structure in a system with eigenbranes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Stationary bubbles and their tunneling channels toward trivial geometry

In the path integral approach, one has to sum over all histories that start from the same initial condition in order to obtain the final condition as a superposition of histories. Applying this into black hole dynamics, we consider stable and unstable stationary bubbles as a reasonable and regular initial condition. We find examples where the bubble can either form a black hole or tunnel toward a trivial geometry, i.e., with no singularity nor event horizon. We investigate the dynamics and tunneling channels of true vacuum bubbles for various tensions. In particular, in line with the idea of superposition of geometries, we build a classically stable stationary thin-shell solution in a Minkowski background where its fate is probabilistically given by non-perturbative effects. Since there exists a tunneling channel toward a trivial geometry in the entire path integral, the entire information is encoded in the wave function. This demonstrates that the unitarity is preserved and there is no loss of information when viewed from the entire wave function of the universe, whereas a semi-classical observer, who can see only a definitive geometry, would find an effective loss of information. Ultimately, this may provide a resolution to the information loss dilemma.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Computing real time correlation functions on a hybrid classical/quantum computer

Quantum devices may overcome limitations of classical computers in studies of nuclear structure functions and parton Wigner distributions of protons and nuclei. In this talk, we discuss a worldline approach to compute nuclear structure functions in the high energy Regge limit of QCD using a hybrid quantum computer, by expressing the fermion determinant in the QCD path integral as a quantum mechanical path integral over 0 + 1-dimensional fermionic and bosonic world-lines in background gauge fields. Our simplest example of computing the well-known dipole model result for the structure function F 2 in the high energy Regge limit is feasible with NISQ era technology using few qubits and shallow circuits. We note this example can be scaled up in complexity and extended in scope to compute structure functions, scattering amplitudes and other real-time correlation functions in QCD, relevant for example to describe non-equilibrium transport of quarks and gluons in a Quark-Gluon-Plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Real-time spin systems from lattice field theory

We construct a lattice field theory method for computing the real-time dynamics of spin systems in a thermal bath. This is done by building on previous work of Takano with Schwinger-Keldysh and functional differentiation techniques. We derive a Schwinger-Keldysh path integral for generic spin Hamiltonians, then demonstrate the method on a simple system. Our path integral has a sign problem, which generally requires exponential run time in the system size, but requires only linear storage. The latter may place this method at an advantage over exact diagonalization, which is exponential in both. Our path integral is amenable to contour deformations, a technique for reducing sign problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The wavefunction of a quantum S 1 × S 2 universe

We study quantum gravity corrections to the no-boundary wavefunction describing a universe with spatial topology S 1 × S 2 . It has been suggested that quantum effects become increasingly important when the size of the circle is large relative to the sphere. In this paper, we confirm this claim by an explicit four-dimensional one-loop calculation of the gravitational path integral preparing such a state. In the process, we clarify some aspects of the gravitational path integral on complex spacetimes. These quantum corrections play a crucial role in ensuring that the norm of the wavefunction is naturally expressed in terms of a path integral over S 2 × S 2 at the classical level. We extend some of the analysis to more general spatial topologies, as well as to the inclusion of fermions.

cosmological models↗
DOE OSTI2016

Worldline approach for numerical computation of electromagnetic Casimir energies: Scalar field coupled to magnetodielectric media

Here, we present a worldline method for the calculation of Casimir energies for scalar fields coupled to magnetodielectric media. The scalar model we consider may be applied in arbitrary geometries, and it corresponds exactly to one polarization of the electromagnetic field in planar layered media. Starting from the field theory for electromagnetism, we work with the two decoupled polarizations in planar media and develop worldline path integrals, which represent the two polarizations separately, for computing both Casimir and Casimir-Polder potentials. We then show analytically that the path integrals for the transverse-electric polarization coupled to a dielectric medium converge to the proper solutions in certain special cases, including the Casimir-Polder potential of an atom near a planar interface, and the Casimir energy due to two planar interfaces. We also evaluate the path integrals numerically via Monte Carlo path-averaging for these cases, studying the convergence and performance of the resulting computational techniques. Lastly, while these scalar methods are only exact in particular geometries, they may serve as an approximation for Casimir energies for the vector electromagnetic field in other geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

On exact-WKB analysis, resurgent structure, and quantization conditions

There are two well-known approaches to studying nonperturbative aspects of quantum mechanical systems: saddle point analysis of the partition functions in Euclidean path integral formulation and the exact-WKB analysis based on the wave functions in the Schrödinger equation. In this work, based on the quantization conditions obtained from the exact-WKB method, we determine the relations between the two formalism and in particular show how the two Stokes phenomena are connected to each other: the Stokes phenomenon leading to the ambiguous contribution of different sectors of the path integral formulation corresponds to the change of the “topology” of the Stoke curves in the exact-WKB analysis. We also clarify the equivalence of different quantization conditions including Bohr-Sommerfeld, path integral and Gutzwiller’s ones. In particular, by reorganizing the exact quantization condition, we improve Gutzwiller’s analysis in a crucial way by bion contributions (incorporating complex periodic paths) and turn it into an exact result. Furthermore, we argue the novel meaning of quasi-moduli integral and provide a relation between the Maslov index and the intersection number of Lefschetz thimbles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

On observers in holographic maps

A straightforward gravitational path integral calculation implies that closed universes are trivial, described by a one dimensional Hilbert space. Two recent papers by Harlow-Usatyuk-Zhao and Abdalla-Antonini-Iliesiu-Levine have sought to ameliorate this issue by defining special rules to incorporate observers into the path integral. However, the proposed rules are different, leading to differing results for the Hilbert space dimension. Moreover, the former work offers a holographic map realized using a non-isometric code construction to complement their path integral result and clarify its physics. In this work, we propose a non-isometric code that implements the second construction, allowing thorough comparison. Our prescription may be thought of as simply removing the portion of the map that acts on the observer, while preserving the rest, creating an effective holographic boundary at the observer-environment interface. This proposal can be directly applied to general holographic maps for both open and closed universes of any dimension.

AdS-CFT correspondence↗
DOE OSTI2019

Lifting of states in 2-dimensional N = 4 supersymmetric CFTs

We consider states of the D1-D5 CFT where only the left-moving sector is excited. As we deform away from the orbifold point, some of these states will remain BPS while others can `lift'. The lifting can be computed by a path integral containing two twist deformations; however, the relevant 4-point amplitude cannot be computed explicitly in many cases. We analyze an older proposal by Gava and Narain where the lift can be computed in terms of a finite number of 3-point functions. A direct Hamiltonian decomposition of the path integral involves an infinite number of 3-point functions, as well the first order correction to the starting state. We note that these corrections to the state account for the infinite number of 3-point functions arising from higher energy states, and one can indeed express the path-integral result in terms of a finite number of 3-point functions involving only the leading order states that are degenerate. The first order correction to the supercharge $\bar G^{(1)}$ gets replaced by a projection $\bar G^{(P)}$; this projected operator can also be used to group the states into multiplets whose members have the same lifting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Equilibrium clumped-isotope effects in doubly substituted isotopologues of ethane

We combine path-integral Monte Carlo methods with a new intramolecular potential energy surface to quantify the equilibrium enrichment of doubly substituted ethane isotopologues due to clumped-isotope effects. Ethane represents the simplest molecule to simultaneously exhibit 13C–13C, 13C–D, and D–D clumped-isotope effects, and the analysis of corresponding signatures may provide useful geochemical and biogeochemical proxies of formation temperatures or reaction pathways. Utilizing path-integral statistical mechanics, we predict equilibrium fractionation factors that fully incorporate nuclear quantum effects, such as anharmonicity and rotational-vibrational coupling which are typically neglected by the widely used Urey model. The magnitude of the calculated fractionation factors for the doubly substituted ethane isotopologues indicates that isotopic clumping can be observed if rare-isotope substitutions are separated by up to three chemical bonds, but the diminishing strength of these effects suggests that enrichment at further separations will be negligible. The Urey model systematically underestimates enrichment due to 13C–D and D–D clumped-isotope effects in ethane, leading to small relative errors in the apparent equilibrium temperature, ranging from 5 K at 273.15 K to 30 K at 873.15 K. We additionally note that the rotameric dependence of isotopologue enrichment must be carefully considered when using the Urey model, whereas the path-integral calculations automatically account for such effects due to configurational sampling. These findings are of direct relevance to future clumped-isotope studies of ethane, as well as studies of 13C–13C, 13C–D, and D–D clumped-isotope effects in other hydrocarbons.

Geochemistry & Geophysics↗
DOE OSTI2015

Quantum-classical correspondence in steady states of nonadiabatic systems

We first present nonadiabatic path integral which is exact formulation of quantum dynamics in nonadiabatic systems. Then, by applying the stationary phase approximations to the nonadiabatic path integral, a semiclassical quantization condition, i.e., quantum-classical correspondence, for steady states of nonadiabatic systems is presented as a nonadiabatic trace formula. The present quantum-classical correspondence indicates that a set of primitive hopping periodic orbits, which are invariant under time evolution in the phase space of the slow degree of freedom, should be quantized. The semiclassical quantization is then applied to a simple nonadiabatic model and accurately reproduces exact quantum energy levels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗