Search NASA⌕ Search

SEARCH · Search NASA

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 OSTI2019

Path-accelerated stochastic molecular dynamics: Parallel-in-time integration using path integrals

Massively parallel computer architectures create new opportunities for the performance of long-time scale molecular dynamics (MD) simulations. Here, we introduce the path-accelerated molecular dynamics method that takes advantage of distributed computing to reduce the wall-clock time of MD simulation via parallelization with respect to stochastic MD time steps. The marginal distribution for the time evolution of a system is expressed in terms of a path integral, enabling the use of path sampling techniques to numerically integrate MD trajectories. By parallelizing the evaluation of the path action with respect to time and by initializing the path configurations from a nonequilibrium distribution, the algorithm enables significant speedups in terms of the length of MD trajectories that can be integrated in a given amount of wall-clock time. The method is demonstrated for Brownian dynamics, although it is generalizable to other stochastic equations of motion including open systems. We apply the method to two simple systems, a harmonic oscillator and a Lennard-Jones liquid, and we show that in comparison to the conventional Euler integration scheme for Brownian dynamics, the new method can reduce the wall-clock time for integrating trajectories of a given length by more than three orders of magnitude in the former system and more than two in the latter. This new method for parallelizing MD in the dimension of time can be trivially combined with algorithms for parallelizing the MD force evaluation to achieve further speedup.

Rosa-Raíces, Jorge L.↗
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↗
NASA NTRS1993

Prediction of stable tearing of 2024-T3 aluminum alloy using the crack-tip opening angle approach

In this study, the crack-tip opening angle (CTOA) approach was incorporated into a damage growth finite element program, MADGIC (Micromechanics Analysis and Damage Growth in Composites), and was used to predict stable tearing in a middle-crack tension 2024-T3 aluminum alloy specimen. The MADGIC code is a displacement based finite element program implemented with an incremental elastic-plastic algorithm used to model elastic-plastic behavior and a nodal splitting and nodal force relaxation algorithm used to generate crack surfaces. Predictions of the applied stress as a function of crack extension and applied stress as a function of load-line displacement were in good agreement with experiments and with similar predictions made using an existing finite element program, ZIP2D. In addition, path integrals, namely, the J-integral and T*-integral, were also evaluated and compared with the CTOA approach. There appears to be a weak relationship between the CTOA and the T*-integral evaluated on a specific integration path during crack extension beyond maximum applied stress. This study further verifies that the CTOA can be used as an effective elastic-plastic fracture mechanics parameter to predict crack growth.

Bakuckas, J. G., Jr.↗
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↗
NASA NTRS2019

Physical Evaluation of GPM DPR Single- and Dual-Wavelength Algorithms

A physical evaluation of the rain profiling retrieval algorithms for the Dual-frequency Precipitation Radar (DPR) aboard the Global Precipitation Measurement (GPM) core satellite is carried out by applying them to the hydrometeor profiles generated from measured raindrop size distributions (DSD). The DSD-simulated radar profiles are used as input to the algorithms, and their estimates of hydrometeors' parameters are compared with the same quantities derived directly from the DSD data (or truth). The retrieval accuracy is assessed by the degree to which the estimates agree with the truth. To check the validity and robustness of the retrievals, the profiles are constructed for cases ranging from fully correlated (or uniform) to totally uncorrelated DSDs along the columns. Investigation into the sensitivity of the retrieval results to the model assumptions is made to characterize retrieval uncertainties and identify error sources. Comparisons between the single- and dual-wavelength algorithm performance are carried out with either a single- or dual-wavelength constraint of the path integral or differential path integral attenuation. The results suggest that the DPR dual-wavelength algorithm generally provides accurate range-profiled estimates of rainfall rate and mass-weighted diameter with the dual-wavelength estimates superior in accuracy to those from the single-wavelength retrievals.

Liang Liao↗
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↗