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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2020

Developing new and understanding old approximations in TDDFT

When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange–correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn–Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange–correlation potential from a coupling-constant path integration.

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

Topological field theory of dynamical systems

Here, it is shown that the path-integral representation of any stochastic or deterministic continuous-time dynamical model is a cohomological or Witten-type topological field theory, i.e., a model with global topological supersymmetry (Q-symmetry). As many other supersymmetries, Q-symmetry must be perturbatively stable due to what is generically known as non-renormalization theorems. As a result, all (equilibrium) dynamical models are divided into three major categories: Markovian models with unbroken Q-symmetry, chaotic models with Q-symmetry spontaneously broken on the mean-field level by, e.g., fractal invariant sets (e.g., strange attractors), and intermittent or self-organized critical (SOC) models with Q-symmetry dynamically broken by the condensation of instanton-antiinstanton configurations (earthquakes, avalanches, etc.) SOC is a full-dimensional phase separating chaos and Markovian dynamics. In the deterministic limit, however, antiinstantons disappear and SOC collapses into the 'edge of chaos.' Goldstone theorem stands behind spatio-temporal self-similarity of Q-broken phases known under such names as algebraic statistics of avalanches, 1/f noise, sensitivity to initial conditions, etc. Other fundamental differences of Q-broken phases is that they can be effectively viewed as quantum dynamics and that they must also have time-reversal symmetry spontaneously broken. Q-symmetry breaking in non-equilibrium situations (quenches, Barkhausen effect, etc.) is also briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Rayleigh-Taylor-induced magnetic fields in laser-irradiated plastic foils

Experimental observations of magnetic fields generated by Rayleigh-Taylor growth in laser-irradiated planar foils are presented. X-ray and monoenergetic proton radiographic techniques were used to probe plastic foils with seeded surface perturbations at different times during the evolution. Protons deflected by fields in the target cause modulations in proton fluence at the seed wavelength of 120 {mu}m. Path-integrated magnetic-field strengths were inferred from modulations in proton fluence using a discrete-Fourier-transform analysis technique and found to increase from 10 to 100 T-{mu}m during linear growth. Electron thermal conduction was shown to be unaffected by Rayleigh-Taylor-induced magnetic fields during the linear growth phase.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Quantum mechanics in fractional and other anomalous spacetimes

We formulate quantum mechanics in spacetimes with real-order fractional geometry and more general factorizable measures. In spacetimes where coordinates and momenta span the whole real line, Heisenberg's principle is proven and the wave-functions minimizing the uncertainty are found. In spite of the fact that ordinary time and spatial translations are broken and the dynamics is not unitary, the theory is in one-to-one correspondence with a unitary one, thus allowing us to employ standard tools of analysis. These features are illustrated in the examples of the free particle and the harmonic oscillator. While fractional (and the more general anomalous-spacetime) free models are formally indistinguishable from ordinary ones at the classical level, at the quantum level they differ both in the Hilbert space and for a topological term fixing the classical action in the path integral formulation. Thus, all non-unitarity in fractional quantum dynamics is encoded in a contribution depending only on the initial and final states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Criteria for the accuracy of small polaron quantum master equation in simulating excitation energy transfer dynamics

The small polaron quantum master equation (SPQME) proposed by Jang et al. [J. Chem. Phys. 129, 101104 (2008)] is a promising approach to describe coherent excitation energy transfer dynamics in complex molecular systems. To determine the applicable regime of the SPQME approach, we perform a comprehensive investigation of its accuracy by comparing its simulated population dynamics with numerically exact quasi-adiabatic path integral calculations. We demonstrate that the SPQME method yields accurate dynamics in a wide parameter range. Furthermore, our results show that the accuracy of polaron theory depends strongly upon the degree of exciton delocalization and timescale of polaron formation. Finally, we propose a simple criterion to assess the applicability of the SPQME theory that ensures the reliability of practical simulations of energy transfer dynamics with SPQME in light-harvesting systems.

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

Refinement of the experimental dynamic structure factor for liquid para-hydrogen and ortho-deuterium using semi-classical quantum simulation

The dynamic structure factor of liquid para-hydrogen and ortho-deuterium in corresponding thermodynamic states (T = 20.0 K, n = 21.24 nm{sup −3}) and (T = 23.0 K, n = 24.61 nm{sup −3}), respectively, has been computed by both the Feynman-Kleinert linearized path-integral (FK-LPI) and Ring-Polymer Molecular Dynamics (RPMD) methods and compared with Inelastic X Ray Scattering spectra. The combined use of computational and experimental methods enabled us to reduce experimental uncertainties in the determination of the true sample spectrum. Furthermore, the refined experimental spectrum of para-hydrogen and ortho-deuterium is consistently reproduced by both FK-LPI and RPMD results at momentum transfers lower than 12.8 nm{sup −1}. At larger momentum transfers the FK-LPI results agree with experiment much better for ortho-deuterium than for para-hydrogen. More specifically we found that for k ∼ 20.0 nm{sup −1} para-hydrogen provides a test case for improved approximations to quantum dynamics.

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

Polymer escape from a confining potential

The rate of escape of polymers from a two-dimensionally confining potential well has been evaluated using self-avoiding as well as ideal chain representations of varying length, up to 80 beads. Long timescale Langevin trajectories were calculated using the path integral hyperdynamics method to evaluate the escape rate. A minimum is found in the rate for self-avoiding polymers of intermediate length while the escape rate decreases monotonically with polymer length for ideal polymers. The increase in the rate for long, self-avoiding polymers is ascribed to crowding in the potential well which reduces the free energy escape barrier. An effective potential curve obtained using the centroid as an independent variable was evaluated by thermodynamic averaging and Kramers rate theory then applied to estimate the escape rate. While the qualitative features are well reproduced by this approach, it significantly overestimates the rate, especially for the longer polymers. The reason for this is illustrated by constructing a two-dimensional effective energy surface using the radius of gyration as well as the centroid as controlled variables. This shows that the description of a transition state dividing surface using only the centroid fails to confine the system to the region corresponding to the free energy barrier and this problem becomes more pronounced the longer the polymer is. A proper definition of a transition state for polymer escape needs to take into account the shape as well as the location of the polymer.

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

Water dimer equilibrium constant calculation: A quantum formulation including metastable states

We present a full quantum evaluation of the water second virial coefficient B(T) based on the Takahashi-Imada second order approximation. As the associated trace Tr[e{sup −βH{sub A}{sub B}}−e{sup −βH{sub A}{sub B}{sup o}}] is performed in the coordinate representation, it does also include contribution from the whole continuum, i.e., resonances and collision pairs of monomers. This approach is compared to a Path Integral Monte Carlo evaluation of this coefficient by Schenter [J. Chem. Phys. 117, 6573 (2002)] for the TIP4P potential and shown to give extremely close results in the low temperature range (250–450 K) reported. Using a recent ab initio flexible potential for the water dimer, this new formulation leads to very good agreement with experimental values over the whole range of temperatures available. The virial coefficient is then used in the well known relation K{sub p}(T) = −(B(T) − b{sub M})/RT where the excluded volume b{sub M} is assimilated to the second virial coefficient of pure water monomer vapor and approximated from the inner repulsive part of the interaction potential. This definition, which renders b{sub M} temperature dependent, allows us to retrieve the 38 cm{sup 3} mol{sup −1} value commonly used, at room temperature. The resulting values for K{sub p}(T) are in agreement with available experimental data obtained from infrared absorption spectra of water vapor.

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

Mixed quantum classical calculation of proton transfer reaction rates: From deep tunneling to over the barrier regimes

We present mixed quantum classical calculations of the proton transfer (PT) reaction rates represented by a double well system coupled to a dissipative bath. The rate constants are calculated within the so called nontraditional view of the PT reaction, where the proton motion is quantized and the solvent polarization is used as the reaction coordinate. Quantization of the proton degree of freedom results in a problem of non-adiabatic dynamics. By employing the reactive flux formulation of the rate constant, the initial sampling starts from the transition state defined using the collective reaction coordinate. Dynamics of the collective reaction coordinate is treated classically as over damped diffusive motion, for which the equation of motion can be derived using the path integral, or the mixed quantum classical Liouville equation methods. The calculated mixed quantum classical rate constants agree well with the results from the numerically exact hierarchical equation of motion approach for a broad range of model parameters. Moreover, we are able to obtain contributions from each vibrational state to the total reaction rate, which helps to understand the reaction mechanism from the deep tunneling to over the barrier regimes. The numerical results are also compared with those from existing approximate theories based on calculations of the non-adiabatic transmission coefficients. It is found that the two-surface Landau-Zener formula works well in calculating the transmission coefficients in the deep tunneling regime, where the crossing point between the two lowest vibrational states dominates the total reaction rate. When multiple vibrational levels are involved, including additional crossing points on the free energy surfaces is important to obtain the correct reaction rate using the Landau-Zener formula.

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

The Nosé–Hoover looped chain thermostat for low temperature thawed Gaussian wave-packet dynamics

We have used a generalised coherent state resolution of the identity to map the quantum canonical statistical average for a general system onto a phase-space average over the centre and width parameters of a thawed Gaussian wave packet. We also propose an artificial phase-space density that has the same behaviour as the canonical phase-space density in the low-temperature limit, and have constructed a novel Nosé–Hoover looped chain thermostat that generates this density in conjunction with variational thawed Gaussian wave-packet dynamics. This forms a new platform for evaluating statistical properties of quantum condensed-phase systems that has an explicit connection to the time-dependent Schrödinger equation, whilst retaining many of the appealing features of path-integral molecular dynamics.

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

How to remove the spurious resonances from ring polymer molecular dynamics

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

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

On the room-temperature phase diagram of high pressure hydrogen: An ab initio molecular dynamics perspective and a diffusion Monte Carlo study

The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results.

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

Adding flexibility to the “particles-on-a-sphere” model for large-amplitude motion: POSflex force field for protonated methane

The so-called “particles-on-a-sphere” (POS) model has been introduced a while ago in order to describe in simple terms large-amplitude motion of polyatomic hydrides, XH{sub n}. The POS model of protonated methane, CH{sub 5}{sup +}, has been shown to capture well the essence of the fluxional nature of this enigmatic floppy molecule. Here, we extend this model to the POSflex force field by adding flexibility to the C–H bonds, which are constrained to a common fixed bond length in the original model. This makes the present model extremely efficient for computer simulation, including path integral molecular dynamics in order to assess the crucial quantum effects on nuclear motion at low temperatures. Moreover, the POSflex force field can be conveniently used to study microsolvation effects upon combining it with intermolecular pair potentials to account for solute-solvent interactions. Upon computing static properties as well as thermal and quantum fluctuation effects at ambient and low temperatures, respectively, it is shown that the POSflex model is very well suited to describe the structural properties of bare CH{sub 5}{sup +}, including hydrogen scrambling and thus fluxionality in the first place. The far- to mid-infrared spectrum up to the bending band is roughly described, whereas the model fails to account for the well-structured stretching band by construction.

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

Entropic penalties in circular DNA assembly

The thermodynamic properties of DNA circular molecules are investigated by a new path integral computational method which treats in the real space the fundamental forces stabilizing the molecule. The base pair and stacking contributions to the classical action are evaluated separately by simulating a broad ensemble of twisted conformations. We obtain, for two short sequences, a free energy landscape with multiple wells corresponding to the most convenient values of helical repeat. Our results point to a intrinsic flexibility of the circular structures in which the base pair fluctuations move the system from one well to the next thus causing the local unwinding of the helix. The latter is more pronounced in the shorter sequence whose cyclization causes a higher bending stress. The entropic reductions associated to the formation of the ordered helicoidal structure are estimated.

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

Ice and water droplets on graphite: A comparison of quantum and classical simulations

Ice and water droplets on graphite have been studied by quantum path integral and classical molecular dynamics simulations. The point-charge q-TIP4P/F potential was used to model the interaction between flexible water molecules, while the water-graphite interaction was described by a Lennard-Jones potential previously used to reproduce the macroscopic contact angle of water droplets on graphite. Several energetic and structural properties of water droplets with sizes between 10{sup 2} and 10{sup 3} molecules were analyzed in a temperature interval of 50–350 K. The vibrational density of states of crystalline and amorphous ice drops was correlated to the one of ice Ih to assess the influence of the droplet interface and molecular disorder on the vibrational properties. The average distance of covalent OH bonds is found 0.01 Å larger in the quantum limit than in the classical one. The OO distances are elongated by 0.03 Å in the quantum simulations at 50 K. Bond distance fluctuations are large as a consequence of the zero-point vibrations. The analysis of the H-bond network shows that the liquid droplet is more structured in the classical limit than in the quantum case. The average kinetic and potential energy of the ice and water droplets on graphite has been compared with the values of ice Ih and liquid water as a function of temperature. The droplet kinetic energy shows a temperature dependence similar to the one of liquid water, without apparent discontinuity at temperatures where the droplet is solid. However, the droplet potential energy becomes significantly larger than the one of ice or water at the same temperature. In the quantum limit, the ice droplet is more expanded than in a classical description. Liquid droplets display identical density profiles and liquid-vapor interfaces in the quantum and classical limits. The value of the contact angle is not influenced by quantum effects. Contact angles of droplets decrease as the size of the water droplet increases which implies a positive sign of the line tension of the droplet.

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

Simultaneous streak and frame interferometry for electron density measurements of laser produced plasmas

A system of two collinear probe beams with different wavelengths and pulse durations was used to capture simultaneously snapshot interferograms and streaked interferograms of laser produced plasmas. The snapshots measured the two dimensional, path-integrated, electron density on a charge-coupled device while the radial temporal evolution of a one dimensional plasma slice was recorded by a streak camera. This dual-probe combination allowed us to select plasmas that were uniform and axisymmetric along the laser direction suitable for retrieving the continuous evolution of the radial electron density of homogeneous plasmas. Demonstration of this double probe system was done by measuring rapidly evolving plasmas on time scales less than 1 ns produced by the interaction of femtosecond, high intensity, laser pulses with argon gas clusters. Experiments aimed at studying homogeneous plasmas from high intensity laser-gas or laser-cluster interaction could benefit from the use of this probing scheme.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Role of temperature on static correlational properties in a spin-polarized electron gas

We have studied the effect of temperature on the static correlational properties of a spin-polarized three-dimensional electron gas (3DEG) over a wide coupling and temperature regime. This problem has been very recently studied by Brown et al. using the restricted path-integral Monte Carlo (RPIMC) technique in the warm-dense regime. To this endeavor, we have used the finite temperature version of the dynamical mean-field theory of Singwi et al, the so-called quantum STLS (qSTLS) approach. The static density structure factor and the static pair-correlation function are calculated, and compared with the RPIMC simulation data. We find an excellent agreement with the simulation at high temperature over a wide coupling range. However, the agreement is seen to somewhat deteriorate with decreasing temperature. The pair-correlation function is found to become small negative for small electron separation. This may be attributed to the inadequacy of the mean-field theory in dealing with the like spin electron correlations in the strong-coupling domain. A nice agreement with RPIMC data at high temperature seems to arise due to weakening of both the exchange and coulomb correlations with rising temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗