-
Times Square sampling: an adaptive algorithm for free energy estimation
Authors:
Cristian Predescu,
Michael Snarski,
Avi Robinson-Mosher,
Duluxan Sritharan,
Tamas Szalay,
David E. Shaw
Abstract:
Estimating free energy differences, an important problem in computational drug discovery and in a wide range of other application areas, commonly involves a computationally intensive process of sampling a family of high-dimensional probability distributions and a procedure for computing estimates based on those samples. The variance of the free energy estimate of interest typically depends strongl…
▽ More
Estimating free energy differences, an important problem in computational drug discovery and in a wide range of other application areas, commonly involves a computationally intensive process of sampling a family of high-dimensional probability distributions and a procedure for computing estimates based on those samples. The variance of the free energy estimate of interest typically depends strongly on how the total computational resources available for sampling are divided among the distributions, but determining an efficient allocation is difficult without sampling the distributions. Here we introduce the Times Square sampling algorithm, a novel on-the-fly estimation method that dynamically allocates resources in such a way as to significantly accelerate the estimation of free energies and other observables, while providing rigorous convergence guarantees for the estimators. We also show that it is possible, surprisingly, for on-the-fly free energy estimation to achieve lower asymptotic variance than the maximum-likelihood estimator MBAR, raising the prospect that on-the-fly estimation could reduce variance in a variety of other statistical applications.
△ Less
Submitted 9 December, 2021;
originally announced December 2021.
-
The $\textit{u}$-series: A separable decomposition for electrostatics computation with improved accuracy
Authors:
Cristian Predescu,
Adam K. Lerer,
Ross A. Lippert,
Brian Towles,
J. P. Grossman,
Robert M. Dirks,
David E. Shaw
Abstract:
The evaluation of electrostatic energy for a set of point charges in a periodic lattice is a computationally expensive part of molecular dynamics simulations (and other applications) because of the long-range nature of the Coulomb interaction. A standard approach is to decompose the Coulomb potential into a near part, typically evaluated by direct summation up to a cutoff radius, and a far part, t…
▽ More
The evaluation of electrostatic energy for a set of point charges in a periodic lattice is a computationally expensive part of molecular dynamics simulations (and other applications) because of the long-range nature of the Coulomb interaction. A standard approach is to decompose the Coulomb potential into a near part, typically evaluated by direct summation up to a cutoff radius, and a far part, typically evaluated in Fourier space. In practice, all decomposition approaches involve approximations---such as cutting off the near-part direct sum---but it may be possible to find new decompositions with improved tradeoffs between accuracy and performance. Here we present the $\textit{u-series}$, a new decomposition of the Coulomb potential that is more accurate than the standard (Ewald) decomposition for a given amount of computational effort, and achieves the same accuracy as the Ewald decomposition with approximately half the computational effort. These improvements, which we demonstrate numerically using a lipid membrane system, arise because the $\textit{u}$-series is smooth on the entire real axis and exact up to the cutoff radius. Additional performance improvements over the Ewald decomposition may be possible in certain situations because the far part of the $\textit{u}$-series is a sum of Gaussians, and can thus be evaluated using algorithms that require a separable convolution kernel; we describe one such algorithm that reduces communication latency at the expense of communication bandwidth and computation, a tradeoff that may be advantageous on modern massively parallel supercomputers.
△ Less
Submitted 4 November, 2019;
originally announced November 2019.
-
Entropic effects in large-scale Monte Carlo simulations
Authors:
Cristian Predescu
Abstract:
The efficiency of Monte Carlo samplers is dictated not only by energetic effects, such as large barriers, but also by entropic effects that are due to the sheer volume that is sampled. The latter effects appear in the form of an entropic mismatch or divergence between the direct and reverse trial moves. We provide lower and upper bounds for the average acceptance probability in terms of the Reny…
▽ More
The efficiency of Monte Carlo samplers is dictated not only by energetic effects, such as large barriers, but also by entropic effects that are due to the sheer volume that is sampled. The latter effects appear in the form of an entropic mismatch or divergence between the direct and reverse trial moves. We provide lower and upper bounds for the average acceptance probability in terms of the Renyi divergence of order 1/2. We show that the asymptotic finitude of the entropic divergence is the necessary and sufficient condition for non-vanishing acceptance probabilities in the limit of large dimensions. Furthermore, we demonstrate that the upper bound is reasonably tight by showing that the exponent is asymptotically exact for systems made up of a large number of independent and identically distributed subsystems. For the last statement, we provide an alternative proof that relies on the reformulation of the acceptance probability as a large deviation problem. The reformulation also leads to a class of low-variance estimators for strongly asymmetric distributions. We show that the entropy divergence causes a decay in the average displacements with the number of dimensions n that are simultaneously updated. For systems that have a well-defined thermodynamic limit, the decay is demonstrated to be n^{-1/2} for random-walk Monte Carlo and n^{-1/6} for Smart Monte Carlo (SMC). Numerical simulations of the LJ_38 cluster show that SMC is virtually as efficient as the Markov chain implementation of the Gibbs sampler, which is normally utilized for Lennard-Jones clusters. An application of the entropic inequalities to the parallel tempering method demonstrates that the number of replicas increases as the square root of the heat capacity of the system.
△ Less
Submitted 27 June, 2007; v1 submitted 26 March, 2007;
originally announced March 2007.
-
Sampling diffusive transition paths
Authors:
Thomas F. Miller III,
Cristian Predescu
Abstract:
We address the problem of sampling double-ended diffusive paths. The ensemble of paths is expressed using a symmetric version of the Onsager-Machlup formula, which only requires evaluation of the force field and which, upon direct time discretization, gives rise to a symmetric integrator that is accurate to second order. Efficiently sampling this ensemble requires avoiding the well-known stiffne…
▽ More
We address the problem of sampling double-ended diffusive paths. The ensemble of paths is expressed using a symmetric version of the Onsager-Machlup formula, which only requires evaluation of the force field and which, upon direct time discretization, gives rise to a symmetric integrator that is accurate to second order. Efficiently sampling this ensemble requires avoiding the well-known stiffness problem associated with sampling infinitesimal Brownian increments of the path, as well as a different type of stiffness associated with sampling the coarse features of long paths. The fine-feature sampling stiffness is eliminated with the use of the fast sampling algorithm (FSA), and the coarse-feature sampling stiffness is avoided by introducing the sliding and sampling (S&S) algorithm. A key feature of the S&S algorithm is that it enables massively parallel computers to sample diffusive trajectories that are long in time. We use the algorithm to sample the transition path ensemble for the structural interconversion of the 38-atom Lennard-Jones cluster at low temperature.
△ Less
Submitted 17 February, 2007; v1 submitted 21 October, 2006;
originally announced October 2006.
-
Generalized moments of spectral functions from short-time correlation functions
Authors:
Cristian Predescu
Abstract:
We present an integral transformation capable of extracting moments of arbitrary Paley-Wiener entire functions against a given spectral distribution based solely on short-time values of the correlation function in a small open disk about the origin. The integral is proven to converge absolutely to the expected result for those correlation functions that can be extended analytically to the entire…
▽ More
We present an integral transformation capable of extracting moments of arbitrary Paley-Wiener entire functions against a given spectral distribution based solely on short-time values of the correlation function in a small open disk about the origin. The integral is proven to converge absolutely to the expected result for those correlation functions that can be extended analytically to the entire complex plane, with the possible exception of two branch cuts on the imaginary axis. It is only the existence of an analytic continuation that is required and not the actual values away from the small disk about the origin. If the analytic continuation exists only for a strip |Im(z)| < τ_0, then the integral transformation remains valid for all Paley-Wiener functions obtained by Fourier-Laplace transforming a compactly supported distribution, with the support included in the interval (-2τ_0, 2τ_0). Finally, if the support of the distribution is contained in the interval $(-τ_0, τ_0)$, then the generalized moment can be evaluated from the short-time values of the correlation function exponentially fast
△ Less
Submitted 26 February, 2006; v1 submitted 12 February, 2006;
originally announced February 2006.
-
Highly optimized fourth-order short-time approximation for path integrals
Authors:
Cristian Predescu
Abstract:
This paper is withdrawn due to copyright restrictions. The final version will become available at this url: http://pubs.acs.org/journals/jpcbfh/
This paper is withdrawn due to copyright restrictions. The final version will become available at this url: http://pubs.acs.org/journals/jpcbfh/
△ Less
Submitted 3 December, 2005; v1 submitted 10 September, 2005;
originally announced September 2005.
-
Thermodynamics and equilibrium structure of Ne_38 cluster: Quantum Mechanics versus Classical
Authors:
Cristian Predescu,
Pavel A. Frantsuzov,
Vladimir A. Mandelshtam
Abstract:
The equilibrium properties of classical LJ_38 versus quantum Ne_38 Lennard-Jones clusters are investigated. The quantum simulations use both the Path-Integral Monte-Carlo (PIMC) and the recently developed Variational-Gaussian-Wavepacket Monte-Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The VGW method is used to locate and…
▽ More
The equilibrium properties of classical LJ_38 versus quantum Ne_38 Lennard-Jones clusters are investigated. The quantum simulations use both the Path-Integral Monte-Carlo (PIMC) and the recently developed Variational-Gaussian-Wavepacket Monte-Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The VGW method is used to locate and characterize the low energy states of Ne_38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne_38 is a liquid-like structure. Among the several liquid-like states with energies below the two symmetric states (O_h and C_5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne_38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained.
△ Less
Submitted 9 August, 2005;
originally announced August 2005.
-
Design of high-order short-time approximations as a problem of matching the covariance of a Brownian motion
Authors:
Cristian Predescu
Abstract:
One of the outstanding problems in the numerical discretization of the Feynman-Kac formula calls for the design of arbitrary-order short-time approximations that are constructed in a stable way, yet only require knowledge of the potential function. In essence, the problem asks for the development of a functional analogue to the Gauss quadrature technique for one-dimensional functions. In PRE 69,…
▽ More
One of the outstanding problems in the numerical discretization of the Feynman-Kac formula calls for the design of arbitrary-order short-time approximations that are constructed in a stable way, yet only require knowledge of the potential function. In essence, the problem asks for the development of a functional analogue to the Gauss quadrature technique for one-dimensional functions. In PRE 69, 056701 (2004), it has been argued that the problem of designing an approximation of order νis equivalent to the problem of constructing discrete-time Gaussian processes that are supported on finite-dimensional probability spaces and match certain generalized moments of the Brownian motion. Since Gaussian processes are uniquely determined by their covariance matrix, it is tempting to reformulate the moment-matching problem in terms of the covariance matrix alone. Here, we show how this can be accomplished.
△ Less
Submitted 28 August, 2005; v1 submitted 8 August, 2005;
originally announced August 2005.
-
Moments of spectral functions: Monte Carlo evaluation and verification
Authors:
Cristian Predescu
Abstract:
The subject of the present study is the Monte Carlo path-integral evaluation of the moments of spectral functions. Such moments can be computed by formal differentiation of certain estimating functionals that are infinitely-differentiable against time whenever the potential function is arbitrarily smooth. Here, I demonstrate that the numerical differentiation of the estimating functionals can be…
▽ More
The subject of the present study is the Monte Carlo path-integral evaluation of the moments of spectral functions. Such moments can be computed by formal differentiation of certain estimating functionals that are infinitely-differentiable against time whenever the potential function is arbitrarily smooth. Here, I demonstrate that the numerical differentiation of the estimating functionals can be more successfully implemented by means of pseudospectral methods (e.g., exact differentiation of a Chebyshev polynomial interpolant), which utilize information from the entire interval $(-β\hbar / 2, β\hbar/2)$. The algorithmic detail that leads to robust numerical approximations is the fact that the path integral action and not the actual estimating functional are interpolated. Although the resulting approximation to the estimating functional is non-linear, the derivatives can be computed from it in a fast and stable way by contour integration in the complex plane, with the help of the Cauchy integral formula (e.g., by Lyness' method). An interesting aspect of the present development is that Hamburger's conditions for a finite sequence of numbers to be a moment sequence provide the necessary and sufficient criteria for the computed data to be compatible with the existence of an inversion algorithm. Finally, the issue of appearance of the sign problem in the computation of moments, albeit in a milder form than for other quantities, is addressed.
△ Less
Submitted 22 July, 2005;
originally announced July 2005.
-
Structure of Si(114) determined by global optimization methods
Authors:
F. C. Chuang,
C. V. Ciobanu,
C. Predescu,
C. Z. Wang,
K. M. Ho
Abstract:
In this article we report the results of global structural optimization of the Si(114) surface, which is a stable high-index orientation of silicon. We use two independent procedures recently developed for the determination of surface reconstructions, the parallel-tempering Monte Carlo method and the genetic algorithm. These procedures, coupled with the use of a highly-optimized interatomic pote…
▽ More
In this article we report the results of global structural optimization of the Si(114) surface, which is a stable high-index orientation of silicon. We use two independent procedures recently developed for the determination of surface reconstructions, the parallel-tempering Monte Carlo method and the genetic algorithm. These procedures, coupled with the use of a highly-optimized interatomic potential for silicon, lead to finding a set of possible models for Si(114), whose energies are recalculated with ab-initio density functional methods. The most stable structure obtained here without experimental input coincides with the structure determined from scanning tunneling microscopy experiments and density functional calculations by Erwin, Baski and Whitman [Phys. Rev. Lett. 77, 687 (1996)].
△ Less
Submitted 22 November, 2004;
originally announced November 2004.
-
On the efficient Monte Carlo implementation of path integrals
Authors:
Cristian Predescu
Abstract:
We demonstrate that the Levy-Ciesielski implementation of Lie-Trotter products enjoys several properties that make it extremely suitable for path-integral Monte Carlo simulations: fast computation of paths, fast Monte Carlo sampling, and the ability to use different numbers of time slices for the different degrees of freedom, commensurate with the quantum effects. It is demonstrated that a Monte…
▽ More
We demonstrate that the Levy-Ciesielski implementation of Lie-Trotter products enjoys several properties that make it extremely suitable for path-integral Monte Carlo simulations: fast computation of paths, fast Monte Carlo sampling, and the ability to use different numbers of time slices for the different degrees of freedom, commensurate with the quantum effects. It is demonstrated that a Monte Carlo simulation for which particles or small groups of variables are updated in a sequential fashion has a statistical efficiency that is always comparable to or better than that of an all-particle or all-variable update sampler. The sequential sampler results in significant computational savings if updating a variable costs only a fraction of the cost for updating all variables simultaneously or if the variables are independent. In the Levy-Ciesielski representation, the path variables are grouped in a small number of layers, with the variables from the same layer being statistically independent. The superior performance of the fast sampling algorithm is shown to be a consequence of these observations. Both mathematical arguments and numerical simulations are employed in order to quantify the computational advantages of the sequential sampler, the Levy-Ciesielski implementation of path integrals, and the fast sampling algorithm.
△ Less
Submitted 15 November, 2004;
originally announced November 2004.
-
The fast sampling algorithm for Lie-Trotter products
Authors:
Cristian Predescu
Abstract:
A fast algorithm for path sampling in path integral Monte Carlo simulations is proposed. The algorithm utilizes the Levy-Ciesielski implementation of Lie-Trotter products to achieve a mathematically proven computational cost of n*log_2(n) with the number of time slices n, despite the fact that each path variable is updated separately, for reasons of optimality. In this respect, we demonstrate th…
▽ More
A fast algorithm for path sampling in path integral Monte Carlo simulations is proposed. The algorithm utilizes the Levy-Ciesielski implementation of Lie-Trotter products to achieve a mathematically proven computational cost of n*log_2(n) with the number of time slices n, despite the fact that each path variable is updated separately, for reasons of optimality. In this respect, we demonstrate that updating a group of random variables simultaneously results in loss of efficiency.
△ Less
Submitted 12 November, 2004; v1 submitted 1 November, 2004;
originally announced November 2004.
-
On the efficiency of exchange in parallel tempering Monte Carlo simulations
Authors:
Cristian Predescu,
Mihaela Predescu,
Cristian V. Ciobanu
Abstract:
If accepted, the paper will be available at http://pubs.acs.org/journals/jpcbfh/
If accepted, the paper will be available at http://pubs.acs.org/journals/jpcbfh/
△ Less
Submitted 5 November, 2004; v1 submitted 28 October, 2004;
originally announced October 2004.
-
Optimal choice of dividing surface for the computation of quantum reaction rates
Authors:
Cristian Predescu,
William H. Miller
Abstract:
See http://pubs.acs.org/journals/jpcafh/ for the latest revision
See http://pubs.acs.org/journals/jpcafh/ for the latest revision
△ Less
Submitted 4 November, 2004; v1 submitted 18 August, 2004;
originally announced August 2004.
-
Reconstruction of thermally-symmetrized quantum autocorrelation functions from imaginary-time data
Authors:
Cristian Predescu
Abstract:
In this paper, I propose a technique for recovering quantum dynamical information from imaginary-time data via the resolution of a one-dimensional Hamburger moment problem. It is shown that the quantum autocorrelation functions are uniquely determined by and can be reconstructed from their sequence of derivatives at origin. A general class of reconstruction algorithms is then identified, accordi…
▽ More
In this paper, I propose a technique for recovering quantum dynamical information from imaginary-time data via the resolution of a one-dimensional Hamburger moment problem. It is shown that the quantum autocorrelation functions are uniquely determined by and can be reconstructed from their sequence of derivatives at origin. A general class of reconstruction algorithms is then identified, according to Theorem 3. The technique is advocated as especially effective for a certain class of quantum problems in continuum space, for which only a few moments are necessary. For such problems, it is argued that the derivatives at origin can be evaluated by Monte Carlo simulations via estimators of finite variances in the limit of an infinite number of path variables. Finally, a maximum entropy inversion algorithm for the Hamburger moment problem is utilized to compute the quantum rate of reaction for a one-dimensional symmetric Eckart barrier.
△ Less
Submitted 6 October, 2004; v1 submitted 11 May, 2004;
originally announced May 2004.
-
Phase changes in selected Lennard-Jones X_{13-n}Y_n clusters
Authors:
Dubravko Sabo,
Cristian Predescu,
J. D. Doll,
David L. Freeman
Abstract:
Detailed studies of the thermodynamic properties of selected binary Lennard-Jones clusters of the type X_{13-n}Y_n (where n=1,2,3) are presented. The total energy, heat capacity and first derivative of the heat capacity as a function of temperature are calculated by using the classical and path integral Monte Carlo methods combined with the parallel tempering technique. A modification in the pha…
▽ More
Detailed studies of the thermodynamic properties of selected binary Lennard-Jones clusters of the type X_{13-n}Y_n (where n=1,2,3) are presented. The total energy, heat capacity and first derivative of the heat capacity as a function of temperature are calculated by using the classical and path integral Monte Carlo methods combined with the parallel tempering technique. A modification in the phase change phenomena from the presence of impurity atoms and quantum effects is investigated.
△ Less
Submitted 22 March, 2004;
originally announced March 2004.
-
Reconstruction of silicon surfaces: a stochastic optimization problem
Authors:
Cristian V. Ciobanu,
Cristian Predescu
Abstract:
Over the last two decades, scanning tunnelling microscopy (STM) has become one of the most important ways to investigate the structure of crystal surfaces. STM has helped achieve remarkable successes in surface science such as finding the atomic structure of Si(111) and Si(001). For high-index Si surfaces the information about the local density of states obtained by scanning does not translate d…
▽ More
Over the last two decades, scanning tunnelling microscopy (STM) has become one of the most important ways to investigate the structure of crystal surfaces. STM has helped achieve remarkable successes in surface science such as finding the atomic structure of Si(111) and Si(001). For high-index Si surfaces the information about the local density of states obtained by scanning does not translate directly into knowledge about the positions of atoms at the surface. A commonly accepted strategy for identifying the atomic structure is to propose several possible models and analyze their corresponding {\em simulated} STM images for a match with the experimental ones. However, the number of good candidates for the lowest-energy structure is very large for high-index surfaces, and heuristic approaches are not likely to cover all the relevant structural models. In this article, we take the view that finding the atomic structure of a surface is a problem of stochastic optimization, and we address it as such. We design a general technique for predicting the reconstruction of silicon surfaces with arbitrary orientation, which is based on parallel-tempering Monte Carlo simulations combined with an exponential cooling. The advantages of the method are illustrated using the Si(105) surface as example, with two main results: (a) the correct single-step rebonded structure [e.g., Fujikawa {\em et al.}, Phys. Rev. Lett. 88, 176101 (2002)] is obtained even when starting from the paired-dimer model [Mo {\em et al.}, Phys. Rev. Lett. 65, 1020 (1990)] that was assumed to be correct for many years, and (b) we have found several double-step reconstructions that have lower surface energies than any previously proposed double-step models.
△ Less
Submitted 8 February, 2004;
originally announced February 2004.
-
The incomplete beta function law for parallel tempering sampling of classical canonical systems
Authors:
Cristian Predescu,
Mihaela Predescu,
Cristian V. Ciobanu
Abstract:
We show that the acceptance probability for swaps in the parallel tempering Monte Carlo method for classical canonical systems is given by a universal function that depends on the average statistical fluctuations of the potential and on the ratio of the temperatures. The law, called the incomplete beta function law, is valid in the limit that the two temperatures involved in swaps are close to o…
▽ More
We show that the acceptance probability for swaps in the parallel tempering Monte Carlo method for classical canonical systems is given by a universal function that depends on the average statistical fluctuations of the potential and on the ratio of the temperatures. The law, called the incomplete beta function law, is valid in the limit that the two temperatures involved in swaps are close to one another. An empirical version of the law, which involves the heat capacity of the system, is developed and tested on a Lennard-Jones cluster. We argue that the best initial guess for the distribution of intermediate temperatures for parallel tempering is a geometric progression and we also propose a technique for the computation of optimal temperature schedules. Finally, we demonstrate that the swap efficiency of the parallel tempering method for condensed-phase systems decreases naturally to zero at least as fast as the inverse square root of the dimensionality of the physical system.
△ Less
Submitted 2 December, 2003; v1 submitted 21 October, 2003;
originally announced October 2003.
-
Heat capacity estimators for random series path-integral methods by finite-difference schemes
Authors:
Cristian Predescu,
Dubravko Sabo,
J. D. Doll,
David L. Freeman
Abstract:
Previous heat capacity estimators used in path integral simulations either have large variances that grow to infinity with the number of path variables or require the evaluation of first and second order derivatives of the potential. In the present paper, we show that the evaluation of the total energy by the T-method estimator and of the heat capacity by the TT-method estimator can be implement…
▽ More
Previous heat capacity estimators used in path integral simulations either have large variances that grow to infinity with the number of path variables or require the evaluation of first and second order derivatives of the potential. In the present paper, we show that the evaluation of the total energy by the T-method estimator and of the heat capacity by the TT-method estimator can be implemented by a finite difference scheme in a stable fashion. As such, the variances of the resulting estimators are finite and the evaluation of the estimators requires the potential function only. By comparison with the task of computing the partition function, the evaluation of the estimators requires k + 1 times more calls to the potential, where k is the order of the difference scheme employed. Quantum Monte Carlo simulations for the Ne_13 cluster demonstrate that a second order central-difference scheme should suffice for most applications.
△ Less
Submitted 4 November, 2003; v1 submitted 31 July, 2003;
originally announced July 2003.
-
Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods
Authors:
Cristian Predescu
Abstract:
In this article, I provide significant mathematical evidence in support of the existence of short-time approximations of any polynomial order for the computation of density matrices of physical systems described by arbitrarily smooth and bounded from below potentials. While for Theorem 2, which is ``experimental'', I only provide a ``physicist's'' proof, I believe the present development is math…
▽ More
In this article, I provide significant mathematical evidence in support of the existence of short-time approximations of any polynomial order for the computation of density matrices of physical systems described by arbitrarily smooth and bounded from below potentials. While for Theorem 2, which is ``experimental'', I only provide a ``physicist's'' proof, I believe the present development is mathematically sound. As a verification, I explicitly construct two short-time approximations to the density matrix having convergence orders 3 and 4, respectively. Furthermore, in the Appendix, I derive the convergence constant for the trapezoidal Trotter path integral technique. The convergence orders and constants are then verified by numerical simulations. While the two short-time approximations constructed are of sure interest to physicists and chemists involved in Monte Carlo path integral simulations, the present article is also aimed at the mathematical community, who might find the results interesting and worth exploring. I conclude the paper by discussing the implications of the present findings with respect to the solvability of the dynamical sign problem appearing in real-time Feynman path integral simulations.
△ Less
Submitted 22 December, 2003; v1 submitted 3 June, 2003;
originally announced June 2003.
-
Energy estimators for random series path-integral methods
Authors:
Cristian Predescu,
Dubravko Sabo,
J. D. Doll,
David L. Freeman
Abstract:
We perform a thorough analysis on the choice of estimators for random series path integral methods. In particular, we show that both the thermodynamic (T-method) and the direct (H-method) energy estimators have finite variances and are straightforward to implement. It is demonstrated that the agreement between the T-method and the H-method estimators provides an important consistency check on th…
▽ More
We perform a thorough analysis on the choice of estimators for random series path integral methods. In particular, we show that both the thermodynamic (T-method) and the direct (H-method) energy estimators have finite variances and are straightforward to implement. It is demonstrated that the agreement between the T-method and the H-method estimators provides an important consistency check on the quality of the path integral simulations. We illustrate the behavior of the various estimators by computing the total, kinetic, and potential energies of a molecular hydrogen cluster using three different path integral techniques. Statistical tests are employed to validate the sampling strategy adopted as well as to measure the performance of the parallel random number generator utilized in the Monte Carlo simulation. Some issues raised by previous simulations of the hydrogen cluster are clarified.
△ Less
Submitted 4 November, 2003; v1 submitted 27 May, 2003;
originally announced May 2003.
-
Numerical implementation of some reweighted path integral methods
Authors:
Cristian Predescu,
Dubravko Sabo,
J. D. Doll
Abstract:
The reweighted random series techniques provide finite-dimensional approximations to the quantum density matrix of a physical system that have fast asymptotic convergence. We study two special reweighted techniques that are based upon the Levy-Ciesielski and Wiener-Fourier series, respectively. In agreement with the theoretical predictions, we demonstrate by numerical examples that the asymptoti…
▽ More
The reweighted random series techniques provide finite-dimensional approximations to the quantum density matrix of a physical system that have fast asymptotic convergence. We study two special reweighted techniques that are based upon the Levy-Ciesielski and Wiener-Fourier series, respectively. In agreement with the theoretical predictions, we demonstrate by numerical examples that the asymptotic convergence of the two reweighted methods is cubic for smooth enough potentials. For each reweighted technique, we propose some minimalist quadrature techniques for the computation of the path averages. These quadrature techniques are designed to preserve the asymptotic convergence of the original methods.
△ Less
Submitted 19 May, 2003;
originally announced May 2003.
-
Reweighted techniques: definition and asymptotic convergence
Authors:
Cristian Predescu
Abstract:
I define and characterize the reweighted methods, which are techniques used in conjunction with the random series implementation of the Feynman-Kac formula. I prove several convergence results valid for all series representations and then I specialize the results for the Levy-Ciesielski and Wiener-Fourier series. As opposed to the partial averaging method on which they are based, the reweighted…
▽ More
I define and characterize the reweighted methods, which are techniques used in conjunction with the random series implementation of the Feynman-Kac formula. I prove several convergence results valid for all series representations and then I specialize the results for the Levy-Ciesielski and Wiener-Fourier series. As opposed to the partial averaging method on which they are based, the reweighted techniques do not involve any modification of the physical potential. Rather, the underlying idea is to develop some specialized constructions of the Brownian bridge that enters the Feynman-Kac formula, so as to simulate the partial averaging effect. For the Levy-Ciesielski series representation, I develop a reweighted technique which has o(1/n^2) convergence for potentials having first order Sobolev derivatives. It is suggested that the asymptotic convergence may reach O(1/n^3) for potentials having second order Sobolev derivatives. The method preserves the favorable log(n) scaling for the time necessary to compute a path at a given discretization point. For the Wiener-Fourier series representation, the particular reweighted method designed in the present article is shown to have O(1/n^3) convergence if the potential has second order Sobolev derivatives. The convergence constant has superior dependence with the inverse temperature as compared to the partial averaging method for the same series. Because the expression of the convergence constant does not actually involve the second order derivatives of the potential, it is conjectured that the O(1/n^3) convergence extends to the potentials having first order Sobolev derivatives only.
△ Less
Submitted 9 February, 2003;
originally announced February 2003.
-
Asymptotic convergence of the partial averaging technique
Authors:
Cristian Predescu,
J. D. Doll,
David L. Freeman
Abstract:
We study the asymptotic convergence of the partial averaging method, a technique used in conjunction with the random series implementation of the Feynman-Kac formula. We prove asymptotic bounds valid for most series representations in the case when the potential has first order Sobolev derivatives. If the potential has also second order Sobolev derivatives, we prove a sharper theorem which gives…
▽ More
We study the asymptotic convergence of the partial averaging method, a technique used in conjunction with the random series implementation of the Feynman-Kac formula. We prove asymptotic bounds valid for most series representations in the case when the potential has first order Sobolev derivatives. If the potential has also second order Sobolev derivatives, we prove a sharper theorem which gives the exact asymptotic behavior of the density matrices. The results are then specialized for the Wiener-Fourier series representation. It is found that the asymptotic behavior is o(1/n^2) if the potential has first order Sobolev derivatives. If the potential has second order Sobolev derivatives, the convergence is shown to be O(1/n^3) and we give the exact expressions for the convergence constants of the density matrices.
△ Less
Submitted 27 January, 2003;
originally announced January 2003.
-
Random Series and Discrete Path Integral methods: The Levy-Ciesielski implementation
Authors:
Cristian Predescu,
J. D. Doll
Abstract:
We perform a thorough analysis of the relationship between discrete and series representation path integral methods, which are the main numerical techniques used in connection with the Feynman-Kac formula. First, a new interpretation of the so-called standard discrete path integral methods is derived by direct discretization of the Feynman-Kac formula. Second, we consider a particular random ser…
▽ More
We perform a thorough analysis of the relationship between discrete and series representation path integral methods, which are the main numerical techniques used in connection with the Feynman-Kac formula. First, a new interpretation of the so-called standard discrete path integral methods is derived by direct discretization of the Feynman-Kac formula. Second, we consider a particular random series technique based upon the Levy-Ciesielski representation of the Brownian bridge and analyze its main implementations, namely the primitive, the partial averaging, and the reweighted versions. It is shown that the n=2^k-1 subsequence of each of these methods can also be interpreted as a discrete path integral method with appropriate short-time approximations. We therefore establish a direct connection between the discrete and the random series approaches. In the end, we give sharp estimates on the rates of convergence of the partial averaging and the reweighted Levy-Ciesielski random series approach for sufficiently smooth potentials. The asymptotic rates of convergence are found to be O(1/n^2), in agreement with the rates of convergence of the best standard discrete path integral techniques.
△ Less
Submitted 17 May, 2004; v1 submitted 10 December, 2002;
originally announced December 2002.
-
Optimal series representations for numerical path integral simulations
Authors:
Cristian Predescu,
J. D. Doll
Abstract:
By means of the Ito-Nisio theorem, we introduce and discuss a general approach to series representations of path integrals. We then argue that the optimal basis for both ``primitive'' and partial averaged approaches is the Wiener sine-Fourier basis. The present analysis also suggests a new approach to improving the convergence of primitive path integral methods. Current work indicates that this…
▽ More
By means of the Ito-Nisio theorem, we introduce and discuss a general approach to series representations of path integrals. We then argue that the optimal basis for both ``primitive'' and partial averaged approaches is the Wiener sine-Fourier basis. The present analysis also suggests a new approach to improving the convergence of primitive path integral methods. Current work indicates that this new technique, the ``reweighted'' method, converges as the cube of the number of path variables for ``smooth'' potentials. The technique is based on a special way of approximating the Brownian bridge which enters the Feynman-Kac formula and it does not require the Gaussian transform of the potential for its implementation.
△ Less
Submitted 6 November, 2002; v1 submitted 19 September, 2002;
originally announced September 2002.
-
The Partial Averaging method
Authors:
Cristian Predescu
Abstract:
The partial averaging technique is defined and used in conjunction with the random series implementation of the Feynman-Kac formula. It enjoys certain properties such as good rates of convergence and convergence for potentials with coulombic singularities. In this work, I introduce the reader to the technique and I analyze the basic mathematical properties of the method. I show that the method i…
▽ More
The partial averaging technique is defined and used in conjunction with the random series implementation of the Feynman-Kac formula. It enjoys certain properties such as good rates of convergence and convergence for potentials with coulombic singularities. In this work, I introduce the reader to the technique and I analyze the basic mathematical properties of the method. I show that the method is convergent for all Kato-class potentials that have finite Gaussian transform.
△ Less
Submitted 4 June, 2003; v1 submitted 27 September, 2002;
originally announced September 2002.
-
Local Variational Principle
Authors:
Cristian Predescu
Abstract:
A generalization of the Gibbs-Bogoliubov-Feynman inequality for spinless particles is proven and then illustrated for the simple model of a symmetric double-well quartic potential. The method gives a pointwise lower bound for the finite-temperature density matrix and it can be systematically improved by the Trotter composition rule. It is also shown to produce groundstate energies better than th…
▽ More
A generalization of the Gibbs-Bogoliubov-Feynman inequality for spinless particles is proven and then illustrated for the simple model of a symmetric double-well quartic potential. The method gives a pointwise lower bound for the finite-temperature density matrix and it can be systematically improved by the Trotter composition rule. It is also shown to produce groundstate energies better than the ones given by the Rayleigh-Ritz principle as applied to the groundstate eigenfunctions of the reference potentials. Based on this observation, it is argued that the Local Variational Principle performs better than the equivalent methods based on the centroid path idea and on the Gibbs-Bogoliubov-Feynman variational principle, especially in the range of low temperatures.
△ Less
Submitted 4 June, 2003; v1 submitted 5 September, 2002;
originally announced September 2002.