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 163 records · Page 9

Global Estimates of Errors in Quantum Computation by the Feynman–Vernon Formalism

The operation of a quantum computer is considered as a general quantum operation on a mixed state on many qubits followed by a measurement. The general quantum operation is further represented as a Feynman–Vernon double path integral over the histories of the qubits and of an environment, and afterward tracing out the environment. The qubit histories are taken to be paths on the two-sphere S{sup 2} as in Klauder’s coherent-state path integral of spin, and the environment is assumed to consist of harmonic oscillators initially in thermal equilibrium, and linearly coupled to to qubit operators Ŝ{sub z}. The environment can then be integrated out to give a Feynman–Vernon influence action coupling the forward and backward histories of the qubits. This representation allows to derive in a simple way estimates that the total error of operation of a quantum computer without error correction scales linearly with the number of qubits and the time of operation. It also allows to discuss Kitaev’s toric code interacting with an environment in the same manner.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sign Learning Kink-based (SiLK) Quantum Monte Carlo for molecular systems

The Sign Learning Kink (SiLK) based Quantum Monte Carlo (QMC) method is used to calculate the ab initio ground state energies for multiple geometries of the H{sub 2}O, N{sub 2}, and F{sub 2} molecules. The method is based on Feynman’s path integral formulation of quantum mechanics and has two stages. The first stage is called the learning stage and reduces the well-known QMC minus sign problem by optimizing the linear combinations of Slater determinants which are used in the second stage, a conventional QMC simulation. The method is tested using different vector spaces and compared to the results of other quantum chemical methods and to exact diagonalization. Our findings demonstrate that the SiLK method is accurate and reduces or eliminates the minus sign problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Going with the flow: a new solution to the sign problem

Here, we discuss a novel solution for the sign problem which prevents first principle Monte-Carlo computations of QCD at finite chemical potential (especially important for both the search for the critical point and neutron star physics) as well as real time quantities such as transport coefficients. The solution is based on deforming the region of integration in the path integral into a complex manifold where the sign problem can be mitigated substantially. We explain the new Monte-Carlo algorithm based on this idea and give examples of interacting quantum field theories (bosonic and fermionic) with nonzero chemical potential as well as real time dynamics where this method successfully solves the sign problem. Furthermore this approach generalizes the ”Lefschetz thimble” method that received much attention lately. We also compare/contrast with the complex Langevin method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Can quantum transition state theory be defined as an exact t = 0 + limit?

Here, the definition of the classical transition state theory (TST) as a t → 0 + limit of the flux-side time correlation function relies on the assumption that simultaneous measurement of population and flux is a well defined physical process. However, the noncommutativity of the two measurements in quantum mechanics makes the extension of such a concept to the quantum regime impossible. For this reason, quantum TST (QTST) has been generally accepted as any kind of quantum rate theory reproducing the TST in the classical limit, and there has been a broad consensus that no unique QTST retaining all the properties of TST can be defined. Contrary to this widely held view, Hele and Althorpe (HA) recently suggested that a true QTST can be defined as the exact t → 0 + limit of a certain kind of quantum flux-side time correlation function and that it is equivalent to the ring polymer molecular dynamics (RPMD) TST. This work seeks to question and clarify certain assumptions underlying these suggestions and their implications. First, the time correlation function used by HA as a starting expression is not related to the kinetic rate constant by virtue of linear response theory, which is the first important step in relating a t = 0 + limit to a physically measurable rate. Second, a theoretical analysis calls into question a key step in HA’s proof which appears not to rely on an exact quantum mechanical identity. The correction of this makes the true t = 0 + limit of HA’s QTST different from the RPMD-TST rate expression, but rather equal to the well-known path integral quantum transition state theory rate expression for the case of centroid dividing surface. An alternative quantum rate expression is then formulated starting from the linear response theory and by applying a recently developed formalism of real time dynamics of imaginary time path integrals. In conclusion, it is shown that the t → 0 + limit of the new rate expression vanishes in the exact quantum limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Compressible fluids with Maxwell-type equations, the minimal coupling with electromagnetic field and the Stefan–Boltzmann law

In this work we have obtained a higher-derivative Lagrangian for a charged fluid coupled with the electromagnetic fluid and the Dirac’s constraints analysis was discussed. A set of first-class constraints fixed by noncovariant gauge condition were obtained. The path integral formalism was used to obtain the partition function for the corresponding higher-derivative Hamiltonian and the Faddeev–Popov ansatz was used to construct an effective Lagrangian. Through the partition function, a Stefan–Boltzmann type law was obtained. - Highlights: • Higher-derivative Lagrangian for a charged fluid. • Electromagnetic coupling and Dirac’s constraint analysis. • Partition function through path integral formalism. • Stefan–Boltzmann-kind law through the partition function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Entanglement between two gravitating universes

Abstract We study two disjoint universes in an entangled pure state. When only one universe contains gravity, the path integral for the n th Rényi entropy includes a wormhole between the n copies of the gravitating universe, leading to a standard ‘island formula’ for entanglement entropy consistent with unitarity of quantum information. When both universes contain gravity, gravitational corrections to this configuration lead to a violation of unitarity. However, the path integral is now dominated by a novel wormhole with 2 n boundaries connecting replica copies of both universes. The analytic continuation of this contribution involves a quotient by Z n replica symmetry, giving a cylinder connecting the two universes. When entanglement is large, this configuration has an effective description as a ‘swap wormhole’, a geometry in which the boundaries of the two universes are glued together by a ‘swaperator’. This description allows precise computation of a generalized entropy-like formula for entanglement entropy. The quantum extremal surface computing the entropy lives on the Lorentzian continuation of the cylinder/swap wormhole, which has a connected Cauchy slice stretching between the universes—a realization of the ER = EPR idea. The new wormhole restores unitarity of quantum information.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating quantum instanton calculations of the kinetic isotope effects

Path integral implementation of the quantum instanton approximation currently belongs among the most accurate methods for computing quantum rate constants and kinetic isotope effects, but its use has been limited due to the rather high computational cost. Here, we demonstrate that the efficiency of quantum instanton calculations of the kinetic isotope effects can be increased by orders of magnitude by combining two approaches: The convergence to the quantum limit is accelerated by employing high-order path integral factorizations of the Boltzmann operator, while the statistical convergence is improved by implementing virial estimators for relevant quantities. After deriving several new virial estimators for the high-order factorization and evaluating the resulting increase in efficiency, using ⋅H{sub α} + H{sub β}H{sub γ} → H{sub α}H{sub β} + ⋅ H{sub γ} reaction as an example, we apply the proposed method to obtain several kinetic isotope effects on CH{sub 4} + ⋅ H ⇌ ⋅ CH{sub 3} + H{sub 2} forward and backward reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coarse-grained fixed-point tensor networks and holographic reflected entropy in 3D gravity

We use the framework of fixed-point BCFT tensor networks to present a microscopic CFT derivation of the correspondence between reflected entropy (RE) and entanglement wedge cross section (EW) in AdS 3 /CFT 2 , for both bipartite and multipartite settings. These fixed-point tensor networks, obtained by triangulating Euclidean CFT path integrals, allow us to explicitly construct the canonical purification via cutting-and-gluing CFT path integrals. Employing modular flow in the large-c limit, we demonstrate that these intrinsic CFT manipulations reproduce bulk geometric prescriptions, without assuming the AdS/CFT dictionary. The emergence of bulk geometry is traced to coarse-graining over heavy states in the large-c limit. Universal coarse-grained BCFT data for compact 2D CFTs, through the relation to Liouville theory with ZZ boundary conditions, yields hyperbolic geometry on the Cauchy slice. The corresponding averaged replica partition functions reproduce all candidate EWs, arising from different averaging patterns, with the dominant one providing the correct RE and EW. In this way, many heuristic tensor-network intuitions in toy models are made precise and established directly from intrinsic CFT data.

AdS-CFT correspondence↗

Application of a New Ensemble Conserving Quantum Dynamics Simulation Algorithm to Liquid para-Hydrogen and ortho-Deuterium

Here, we apply the Feynman-Kleinert Quasi-Classical Wigner (FK-QCW) method developed in our previous work [Smith et al., J. Chem. Phys. 142, 244112 (2015)] for the determination of the dynamic structure factor of liquid para-hydrogen and ortho-deuterium at state points of (T = 20.0 K, n = 21.24 nm -3 ) and (T = 23.0 K, n = 24.61 nm -3 ), respectively. When applied to this challenging system, it is shown that this new FK-QCW method consistently reproduces the experimental dynamic structure factor reported by Smith et al. [J. Chem. Phys. 140, 034501 (2014)] for all momentum transfers considered. Moreover, this shows that FK-QCW provides a substantial improvement over the Feynman-Kleinert linearized path-integral method, in which purely classical dynamics are used. Furthermore, for small momentum transfers, it is shown that FK-QCW provides nearly the same results as ring-polymer molecular dynamics (RPMD), thus suggesting that FK-QCW provides a potentially more appealing algorithm than RPMD since it is not formally limited to correlation functions involving linear operators.

36 MATERIALS SCIENCE↗

Holography: Quantum Mechanical Aspects of Black Holes (Final Report)

Current theories of fundamental physics are based on quantum field theory, unifying quantum mechanics and special relativity. A similar unification of quantum field theory with gravity has been elusive. String theory is a consistent model of quantum gravity and it is important to identify general lessons that can be applied to the real world. The main technique that emerged from such studies has been the gravitational path integral. During the period of this grant, the PI has explored the implication of the gravitational path integral in novel setups to uncover new features of rotating black holes in four dimensions, as well as increasing our understanding of some specific black holes in string theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring the impact of ions on oxygen K-edge X-ray absorption spectroscopy in NaCl solution using the GW-Bethe-Salpeter-equation approach

X-ray absorption spectroscopy (XAS) is a powerful experimental tool to probe the local structure in materials with the core hole excitations. Here, the oxygen K-edge XAS spectra of the NaCl solution and pure water are computed by using a recently developed GW-Bethe-Salpeter equation approach, based on configurations modeled by path-integral molecular dynamics with the deep-learning technique. The neural network is trained on ab initio data obtained with strongly constrained and appropriately normed density functional theory. The observed changes in the XAS features of the NaCl solution, compared to those of pure water, are in good agreement between experimental and theoretical results. We provided detailed explanations for these spectral changes that occur when NaCl is solvated in pure water. Specifically, the presence of solvating ion pairs leads to localization of electron-hole excitons. As a result, our theoretical XAS results support the theory that the effects of the solvating ions on the H-bond network are mainly confined within the first hydration shell of ions, however beyond the shell the arrangement of water molecules remains to be comparable to that observed in pure water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of first-principles properties of deuterium–tritium on inertial confinement fusion target designs

A comprehensive knowledge of the properties of high-energy-density plasmas is crucial to understanding and designing low-adiabat, inertial confinement fusion (ICF) implosions through hydrodynamic simulations. Warm-dense-matter (WDM) conditions are routinely accessed by low-adiabat ICF implosions, in which strong coupling and electron degeneracy often play an important role in determining the properties of warm dense plasmas. The WDM properties of deuterium–tritium (DT) mixtures and ablator materials, such as the equation of state, thermal conductivity, opacity, and stopping power, were usually estimated by models in hydro-codes used for ICF simulations. In these models, many-body and quantum effects were only approximately taken into account in the WMD regime. Moreover, the self-consistency among these models was often missing. To examine the accuracy of these models, we have systematically calculated the static, transport, and optical properties of warm dense DT plasmas, using first-principles (FP) methods over a wide range of densities and temperatures that cover the ICF “path” to ignition. These FP methods include the path-integral Monte Carlo (PIMC) and quantum-molecular dynamics (QMD) simulations, which treat electrons with many-body quantum theory. The first-principles equation-of-state table, thermal conductivities (κ{sub QMD}), and first principles opacity table of DT have been self-consistently derived from the combined PIMC and QMD calculations. They have been compared with the typical models, and their effects to ICF simulations have been separately examined in previous publications. In this paper, we focus on their combined effects to ICF implosions through hydro-simulations using these FP-based properties of DT in comparison with the usual model simulations. We found that the predictions of ICF neutron yield could change by up to a factor of ∼2.5; the lower the adiabat of DT capsules, the more variations in hydro-simulations. The FP-based properties of DT are essential for designing ICF ignition targets. Future work on first-principles studies of ICF ablator materials is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flat entanglement spectra in fixed-area states of quantum gravity

We use the Einstein-Hilbert gravitational path integral to investigate gravitational entanglement at leading order O(1/G). We argue that semiclassical states prepared by a Euclidean path integral have the property that projecting them onto a subspace in which the Ryu-Takayanagi or Hubeny-Rangamani-Takayanagi surface has definite area gives a state with a flat entanglement spectrum at this order in gravitational perturbation theory. This means that the reduced density matrix can be approximated as proportional to the identity to the extent that its Renyi entropies Sn are independent of n at this order. The n-dependence of Sn in more general states then arises from sums over the RT/HRT-area, which are generally dominated by different values of this area for each n. This provides a simple picture of gravitational entanglement, bolsters the connection between holographic systems and tensor network models, clarifies the bulk interpretation of algebraic centers which arise in the quantum errorcorrecting description of holography, and strengthens the connection between bulk and boundary modular Hamiltonians described by Jafferis, Lewkowycz, Maldacena, and Suh.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spin-textured neutron beams with orbital angular momentum

In this work, we present a rigorous theoretical framework underpinning the technique of spin-echo modulated small-angle neutron scattering (SEMSANS), and show how the technique can be extended in order to generate spin-textured neutron beams with orbital angular momentum (OAM) via birefringent neutron spin-polarization devices known as magnetic Wollaston prisms. Neutron OAM beams are mathematically characterized by a “cork-screw” phase singularity e iℓΦ about the propagation axis where ℓ is the OAM quantum number. To understand the precise relationship between the emergent OAM state and the variety of spin textures realized by various setups, we have developed a path-integral approach that in the interferometric limit makes a judicious use of magnetic Snell's law. We show that our proposed technique produces a complex two-dimensional pattern of spin-OAM entangled states which may be useful as a probe of quantum magnetic materials. We compare our path-integral approach to the well-known single-path Larmor precession model and present a pedagogical derivation of magnetic Snell's law of refraction for both massive and massless particles based on Maupertuis's action principle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Gluing. Part I. Integrals and symmetries

We review some aspects of the cutting and gluing law in local quantum field theory (QFT) and study it from a new point of view. In particular, we emphasize the description of gluing by a path integral over a space of polarized boundary conditions, which are given by leaves of some Lagrangian foliation in the phase space. We think of this path integral as a non-local (d – 1)-dimensional gluing theory associated to the parent local d-dimensional QFT. This is a novel point of view paving the way for applications of the standard QFT techniques (that do not rely on locality) to the gluing theory. We describe various properties of this procedure and spell out conditions under which symmetries of the parent theory lead to symmetries of the gluing theory. The purpose of this paper is to set up a playground for the companion paper where these techniques are applied to obtain new results in supersymmetric theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗

Renormalization of excitonic properties by polar phonons

Here we employ quasiparticle path integral molecular dynamics to study how the excitonic properties of model semiconductors are altered by electron–phonon coupling. We describe ways within a path integral representation of the system to evaluate the renormalized mass, binding energy, and radiative recombination rate of excitons in the presence of a fluctuating lattice. To illustrate this approach, we consider Fröhlich-type electron–phonon interactions and employ an imaginary time influence functional to incorporate phonon-induced effects nonperturbatively. The effective mass and binding energies are compared with perturbative and variational approaches, which provide qualitatively consistent trends. We evaluate electron-hole recombination rates as mediated through both trap-assisted and bimolecular processes, developing a consistent statistical mechanical approach valid in the reaction limited regime. These calculations demonstrate how phonons screen electron–hole interactions, generically reducing exciton binding energies and increasing their radiative lifetimes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗