ADAPTIVE TWO-STAGE INTEGRATORS FOR SAMPLING ALGORITHMS BASED ON HAMILTONIAN DYNAMICS
|
|
- Leslie Barton
- 5 years ago
- Views:
Transcription
1 ADAPTIVE TWO-STAGE INTEGRATORS FOR SAMPLING ALGORITHMS BASED ON HAMILTONIAN DYNAMICS E. Akhmatskaya a,c, M. Fernández-Pendás a, T. Radivojević a, J. M. Sanz-Serna b a Basque Center for Applied Mathematics (BCAM), Bilbao, Spain b Departamento de Matemáticas, Universidad Carlos III de Madrid, Madrid, Spain c IKERBASQUE, Basque Foundation for Science, Bilbao, Spain
2 SAMPLING WITH HAMILTONIAN DYNAMICS Let H(θ,p)= 1 2 pt M 1 p+u(θ ) is Hamiltonian θ position, p momentum, U potential energy, M mass matrix, (θ,p) R 2D, D - dimension Equations of motion associated with H: dθ dt = M 1 p; dp dt =-U θ (θ ) The algorithms sampling with Hamiltonian dynamics: ª Can be viewed as special cases of Generalized Hybrid Monte Carlo (GHMC) [Kennedy, Pendleton, 21] ª Numerically solve Hamiltonian equations
3 GHMC AT GLANCE Initialization: positions θ, momenta p (drawn from Gaussian distribution), step size h, length of trajectories L, number of samples N, φ Proposal: numerical integration of Hamiltonian equations with the symplectic method Ψ h over L steps and step-size h: Ψ T (θ,p) = ( θ!, p!),t = Lh Sampling: Metropolis accept / reject test # % FΨ T (θ,p) with probability min(1,exp(-βδh)) % ( θ!, p!) = $ % (θ, p) otherwise % &% Momentum flip F : (θ,p) = (θ,-p) ΔH = H(Ψ T (θ,p)) H(θ,p) = H( FΨ T (θ,p)) H(θ,p) Momentum refreshment: " p! $ % " # u! & ' = $ # cosϕ sinϕ sinϕ cosϕ % " ' F $ & # p u % ', < ϕ π / 2, u ~ N(,β 1 M) &
4 METHODS & APPLICATIONS Molecular Simulations Computational Statistics GHMC GSHMC MMHMC Sample with Shadow Hamiltonians HMC Hybrid MC HMC Hamiltonian φ=π/2
5 METHODS & APPLICATIONS Molecular Simulations Computational Statistics GHMC GSHMC MMHMC Sample with Shadow Hamiltonians Shadow Hamiltonians SHMC / S2HMC HMC Hybrid MC Adaptive parameters NUTS HMC Hamiltonian φ=π/2 Delayed rejections LAHMC / ECGHMC
6 IN FOCUS Molecular simulations: ª Hybrid Monte Carlo (HMC)[Duane, et. al, 1987] ª Generalized shadow hybrid Monte Carlo (GSHMC) [Akhmatskaya, Reich, 28] HMC: GHMC with a complete momentum update GSHMC: Importance sampling GHMC
7 IN FOCUS Computational statistics: ª Hamiltonian Monte Carlo (HMC)[Neil, 1993] ª Mixed & Match Hamiltonian Monte Carlo (MMHMC) [Radivojevic, Akhmatskaya, 216] HMC: Hybrid Monte Carlo formulated for statistics MMHMC: GSHMC adjusted for statistics ª For the target density π(θ) of position vector θ, momenta p conjugate to θ, with a mass matrix M (a preconditioner) Potential function: Hamiltonian: U(θ) = log(π (θ)) H(θ,p)= 1 2 pt M 1 p+u(θ )
8 WHAT INTEGRATOR TO CHOOSE? Other suggestions?
9 NUMERICAL INTEGRATORS FOR SIMULATING HAMILTONIAN DYNAMICS Wish list: ª Reversible ª Symplectic (preserve the symplectic structure of the Hamiltonian dynamics) => volume preserving ª Low integration error ª Computationally efficient ª Stable
10 VERLET INTEGRATOR Verlet or Verlet / Stormer or leapfrog: a golden standard for Hamiltonian dynamics based simulation methods ª ª ª Reversible & symplectic Second order of accuracy (global error of O (h 2 ), h time step) One force evaluation per step ª Stability interval: (,2) ª ª Easy to implement Verlet as a splitting integrator: ψ h = ϕ B h/2 ϕ A B h ϕ h/2 ϕ t X - t-flow of X={A, B} Can we beat it?
11 SPLITTING INTEGRATORS Potential competitors to Verlet are members of one-parameter family of 2-stage splitting integrators of Hamiltonian system with H(θ,p)= 1 2 pt M 1 p+u(θ ) A+B b is a parameter of the family; b=.25 for Verlet 2-stage splitting integrators: ψ h = ϕ B bh ϕ A B h/2 ϕ (1 2b)h ϕ A B h/2 ϕ bh ª Capable to provide higher accuracy than Verlet (substeps h/2) ª Can achieve bigger time steps (stability interval up to: (,4)) ª Two force evaluations per step but ª Fair comparison with Verlet: equal numbers of force evaluations h Verlet = h 2 stage / 2; L Verlet = 2L 2 stage
12 2-STAGE INTEGRATORS: SPECIAL CASE 1 Minimum-error (ME) integrator by McLachlan (1995): ª Criterion: to minimize the error εper step ε Ch 3 + O(h 5 ); h ª Resulting parameter: b.1932 ª Performance: ü outperforms Verlet at small time steps ü performance degrades for large h: stability interval (, 2.55)
13 2-STAGE INTEGRATORS: SPECIAL CASE 2 BCSS integrator by Blanes, Casas and Sanz-Serna (214): ª Criterion: to minimize expectation of the energy error,, for < h < 2 Δ = H(ψ h,l (θ,p)) H(θ,p) E(Δ) ρ(h,b) h 4 (2b 2 (1/ 2 b)h 2 + 4b 2 6b +1) 2 ρ(h,b) = 8(2 bh 2 )(2 (1/ 2 b)h 2 )(1 b(1/ 2 b)h 2 ) E(Δ) ª Resulting parameter: b.2113 ª Performance: ü Shows the best performance around h=2 ü Performance may drop for smaller or larger step sizes Any better values of b?
14 ADAPTIVE INTEGRATION APPROACH AIA Idea: to present a parameter b as a function of internal properties (frequencies) of a simulated system and a simulation time step, i.e., in such a way that: b = b(δt,ω) Given a simulation problem (ω)and a simulation time (Δt), b defines the unique integrator which guarantees the best conservation of energy for harmonic forces. ª The resulting integrator: a problem specific two-stage integrator ª Criterion: to minimize max <h<h ρ(h,b), h h(δt,ω) ª AIA uses the upper bound suggested in BCSS method
15 AIA: ALGORITHM Given a simulation problem (ω)and a simulation time (Δt): 1. Compute the non-dimensional quantity Molecular Simulation: S = 2 " D % 6 6 dh / 2 ω $ i ' # i=1 & Statistics : S = 2 ω, D Δt - highest frequencies and dimension of the system, - average energy error obtained from warm-up simulation with 2. If h>4 abort the integration 3. Find the optimal value of the parameter b as b = arg min max ρ(h,b) b (,.25) h (,h ) h = S ωδt (Schlick, 22; Skeel, 1999) dh Δt
16 AIA: REMARKS ª No computational overheads in simulations ª Extended to constrained dynamics: ü Presented the two-stage integrator as two concatenated velocity Verlet steps and combined with RATTLE
17 AIA FOR MODIFIED HAMILTONIAN METHODS ª MAIA: AIA for methods sampling with modified Hamiltonian ª The same ideas as in AIA but the expected value of modified energy Δ = H(ψ h,l (θ,p)) H(θ,p) w.r.t. to the modified density is minimised ª MAIA requires: ü Modified Hamiltonians for 2-stage splitting integrators ü Upper bound for expectation of modified energies
18 MODIFIED HAMILTONIANS FOR SPLITTING INTEGRATORS 4-th order modified Hamiltonians for 2-stage splitting integrators were derived in terms of quantities available during simulation - numerical time derivative of - parameter of a 2-stage integrator
19 MAIA: UNDERLYING ANALYSIS Model problem: Standard harmonic oscillator in non-dimensional variables (standard univariate Gaussian). The solution flow is given by Defining depends on a integrator. obtain: ü Numerical solution is exact solution of modified Hamiltonian ü Numerical solution stays on ellipse ü Maximum Δ: numerical solution starts at ellipse vertices on major semiaxis
20 UPPER BOUND FOR EXPECTATION OF MODIFIED ENERGY ERROR D-variate Gaussian target distribution (d coupled linear oscillators): D j=1 E(Δ) ρ(ω j h) ω j - the angular frequencies of the oscillators
21 INTEGRATOR S PARAMETER VS. TIME STEP.27 Sampling with Hamiltonians Sampling with shadow Hamiltonians b b h AIA BCSS ME h MAIA M-BCSS M-ME
22 EXPECTED ERROR VS. TIME STEP Expected error Sampling with Hamiltonians AIA BCSS ME Expected error Sampling with shadow Hamiltonians MAIA M-BCSS M-ME h h
23 EXPECTED ERROR VS. TIME STEP: ZOOM Expected error AIA BCSS ME MAIA M-BCSS M-ME h
24 TIME STEP vs. EXPECTED ERROR 4 Sampling with Hamiltonians 4 Sampling with shadow Hamiltonians 3 3 h 2 h 2 1 AIA BCSS ME Expected error 1 MAIA M-BCSS M-ME Expected error
25 TIME STEP vs. EXPECTED ERROR: ZOOM 4 3 h 2 1 AIA BCSS ME MAIA M-BCSS M-ME Expected error
26 (M)AIA IMPLEMENTATION ª Can be implemented in any MD/HMC software with no overheads ª Two major implementation steps: 1. Introduce AIA in a pre-processing unit to be run once before a simulation ü find b and pass to a simulator 2. Present 2-stage scheme in a kick/drift factorization form in a simulation module ª Current implementation ü MS: in multihmc-gromacs BCAM in-house modified version of GROMACS ü CS: in BCAM software package Hamiltonians in Computational Statistics (HaiCS)
27 A TESTING (M)AIA 18 o Molecular Simulation B Spider toxin in membrane/water environment: Coarse grained unconstrained system with 781 particles Villiin: Atomistic system with 9389 atoms. Hydrogens are constrained using the SHAKE/RATTLE algorithm Computational Statistics Gaussian distribution: D=1, 2 Metrics ESS time normalized effective sample size Rel. ESS relative ESS w.r.t. ESS achieved with Verlet IACF integrated autocorrelation function
28 ACCEPTANCE RATES: TOXIN & VILLIN AR (%) AR (%) AIA BCSS HMC Δt (fs) Δt (fs) 4 5 AR (%) AR (%) MAIA M-BCSS GSHMC Δt (fs) Δt (fs) 4 5 Toxin Villin
29 IACF: TOXIN IACF x time (h) HMC AIA BCSS IACF x time (h) GSHMC MAIA M-BCSS Δt (fs) Δt (fs)
30 TOXIN-BILAYER DISTANCE 3 HMC 3 HMC Distance (nm) BCSS AIA Target distance True distribution BCSS AIA Time (ps) Distance (nm)
31 ACCEPTANCE RATES: GAUSSIAN 1 HMC MMHMC D=1 75 AR (%) 5 25 AIA BCSS ME MAIA M-BCSS M-ME HMC t MMHMC t 75 D=2 AR (%) t t
32 RELATIVE ESS: GAUSSIAN, D=1 rel. min ESS rel. medess AIA BCSS ME HMC t MAIA M-BCSS M-ME MMHMC t
33 RELATIVE ESS: GAUSSIAN, D=2 rel. min ESS rel. medess AIA BCSS ME HMC t MAIA M-BCSS M-ME MMHMC t
34 ESS: GAUSSIAN min ESS AIA BCSS ME HMC HMC h MAIA M-BCSS M-ME MMHMC MMHMC h D=1 D=2 min ESS h h
35 CONCLUSIONS ª We present an alternative to the standard Verlet integrator, (M)AIA, which offers, for any chosen simulation problem and step size, a system-specific two-stage splitting integrator to provide the best conservation of energy for harmonic forces ª The method was derived for sampling with Hamiltonians and modified Hamiltonians and extended to constrained dynamics ª Testing M(AIA) in molecular simulation and statistical applications demonstrated its superiority over velocity Verlet, BCSS and ME integrators ü never worse than, BCSS and ME by design ü leads to higher acceptance rates and better sampling ü allows for longer time steps ü no computational overheads in simulations
36 WORK IN PROGRESS ª M(AIA) for n-stage splitting integrators (n>2) ª M(AIA) for multiple-time stepping methods ª Predicting a range of optimal time steps using the (M)AIA underlying analysis
37 Guggenheim, Bilbao Modeling & Simulation in Life & Materials Sciences: MSLMS BCAM
Mix & Match Hamiltonian Monte Carlo
Mix & Match Hamiltonian Monte Carlo Elena Akhmatskaya,2 and Tijana Radivojević BCAM - Basque Center for Applied Mathematics, Bilbao, Spain 2 IKERBASQUE, Basque Foundation for Science, Bilbao, Spain MCMSki
More informationEnabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package
Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package Mario Fernández Pendás MSBMS Group Supervised by Bruno Escribano and Elena Akhmatskaya BCAM 18 October
More informationGSHMC: An efficient Markov chain Monte Carlo sampling method. Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe)
GSHMC: An efficient Markov chain Monte Carlo sampling method Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe) 1. Motivation In the first lecture, we started from a
More informationConstant pressure hybrid Monte Carlo simulations in GROMACS
Constant pressure hybrid Monte Carlo simulations in GROMACS Mario Fernández-Pendás 1, Bruno Escribano 1, Tijana Radivojević 1, and Elena Akhmatskaya 1,2 1 Basque Center for Applied Mathematics, E-48009
More informationNew Hybrid Monte Carlo Methods for Efficient Sampling: from Physics to Biology and Statistics
TECHNICAL MATERIAL New Hybrid Monte Carlo Methods for Efficient Sampling: from Physics to Biology and Statistics Elena AKHMATSKAYA 1-3* and Sebastian REICH 4 1 Fujitsu Laboratories of Europe Ltd, Hayes
More informationDownloaded 07/28/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 36, No. 4, pp. A1556 A1580 c 2014 Society for Industrial and Applied Mathematics NUMERICAL INTEGRATORS FOR THE HYBRID MONTE CARLO METHOD SERGIO BLANES, FERNANDO CASAS, AND J.
More informationMathematical Modeling of Chemical Reactions at Basque Center for Applied Mathematics. Simone Rusconi. December 11, 2015
Mathematical Modeling of Chemical Reactions at Basque Center for Applied Mathematics Simone Rusconi December 11, 2015 Basque Center for Applied Mathematics (BCAM) BCAM - Basque Center for Applied Mathematics,
More informationFinite Element Model Updating Using the Separable Shadow Hybrid Monte Carlo Technique
Finite Element Model Updating Using the Separable Shadow Hybrid Monte Carlo Technique I. Boulkaibet a, L. Mthembu a, T. Marwala a, M. I. Friswell b, S. Adhikari b a The Centre For Intelligent System Modelling
More informationCombining stochastic and deterministic approaches within high efficiency molecular simulations
Cent. Eur. J. Math. 1-16 Author version Central European Journal of Mathematics Combining stochastic and deterministic approaches within high efficiency molecular simulations Research Bruno Escribano 1,
More informationDynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June
Dynamic modelling of morphology development in multiphase latex particle Elena Akhmatskaya (BCAM) and Jose M. Asua (POLYMAT) June 7 2012 Publications (background and output) J.M. Asua, E. Akhmatskaya Dynamic
More informationHamiltonian Monte Carlo with Fewer Momentum Reversals
Hamiltonian Monte Carlo with ewer Momentum Reversals Jascha Sohl-Dickstein December 6, 2 Hamiltonian dynamics with partial momentum refreshment, in the style of Horowitz, Phys. ett. B, 99, explore the
More informationCombining stochastic and deterministic approaches within high efficiency molecular simulations
Cent. Eur. J. Math. 11(4) 2013 787-799 DOI: 10.2478/s11533-012-0164-x Central European Journal of Mathematics Combining stochastic and deterministic approaches within high efficiency molecular simulations
More informationHamiltonian Monte Carlo Without Detailed Balance
Jascha Sohl-Dickstein Stanford University, Palo Alto. Khan Academy, Mountain View Mayur Mudigonda Redwood Institute for Theoretical Neuroscience, University of California at Berkeley Michael R. DeWeese
More informationMonte Carlo in Bayesian Statistics
Monte Carlo in Bayesian Statistics Matthew Thomas SAMBa - University of Bath m.l.thomas@bath.ac.uk December 4, 2014 Matthew Thomas (SAMBa) Monte Carlo in Bayesian Statistics December 4, 2014 1 / 16 Overview
More informationHamiltonian Monte Carlo
Hamiltonian Monte Carlo within Stan Daniel Lee Columbia University, Statistics Department bearlee@alum.mit.edu BayesComp mc-stan.org Why MCMC? Have data. Have a rich statistical model. No analytic solution.
More informationConservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt =
Conservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt = d (mv) /dt where p =mv is linear momentum of particle
More informationGradient-based Monte Carlo sampling methods
Gradient-based Monte Carlo sampling methods Johannes von Lindheim 31. May 016 Abstract Notes for a 90-minute presentation on gradient-based Monte Carlo sampling methods for the Uncertainty Quantification
More informationIntroduction to Hamiltonian Monte Carlo Method
Introduction to Hamiltonian Monte Carlo Method Mingwei Tang Department of Statistics University of Washington mingwt@uw.edu November 14, 2017 1 Hamiltonian System Notation: q R d : position vector, p R
More informationRiemann Manifold Methods in Bayesian Statistics
Ricardo Ehlers ehlers@icmc.usp.br Applied Maths and Stats University of São Paulo, Brazil Working Group in Statistical Learning University College Dublin September 2015 Bayesian inference is based on Bayes
More informationTesting a Fourier Accelerated Hybrid Monte Carlo Algorithm
Syracuse University SURFACE Physics College of Arts and Sciences 12-17-2001 Testing a Fourier Accelerated Hybrid Monte Carlo Algorithm Simon Catterall Syracuse University Sergey Karamov Syracuse University
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More information1 Geometry of high dimensional probability distributions
Hamiltonian Monte Carlo October 20, 2018 Debdeep Pati References: Neal, Radford M. MCMC using Hamiltonian dynamics. Handbook of Markov Chain Monte Carlo 2.11 (2011): 2. Betancourt, Michael. A conceptual
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 7 Beyond Atoms: Simulating Molecules David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu Review Fundamentals units, properties, statistical
More informationExponential Integration for Hamiltonian Monte Carlo
Wei-Lun Chao WEILUNC@USC.EDU Justin Solomon 2 JUSTIN.SOLOMON@STANFORD.EDU Dominik L. Michels 2 MICHELS@CS.STANFORD.EDU Fei Sha FEISHA@USC.EDU Department of Computer Science, University of Southern California,
More informationarxiv: v1 [stat.me] 6 Apr 2013
Generalizing the No-U-Turn Sampler to Riemannian Manifolds Michael Betancourt Applied Statistics Center, Columbia University, New York, NY 127, USA Hamiltonian Monte Carlo provides efficient Markov transitions
More informationDiscontinuous Hamiltonian Monte Carlo for discrete parameters and discontinuous likelihoods
Discontinuous Hamiltonian Monte Carlo for discrete parameters and discontinuous likelihoods arxiv:1705.08510v3 [stat.co] 7 Sep 2018 Akihiko Nishimura Department of Biomathematics, University of California
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture Molecular Dynamics Simulation David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu MD of hard disks intuitive Review and Preview collision
More informationGeometry & Dynamics for Markov Chain Monte Carlo
arxiv:1705.02891v1 [stat.co] 8 May 2017 Annu. Rev. Stat. Appl. 2018. 5:1 26 This article s doi: 10.1146/please add article doi)) Copyright c 2018 by Annual Reviews. All rights reserved Geometry & Dynamics
More informationPolynomial Filtered Hybrid Monte Carlo
Polynomial Filtered Hybrid Monte Carlo Waseem Kamleh and Michael J. Peardon CSSM & U NIVERSITY OF A DELAIDE QCDNA VII, July 4th 6th, 2012 Introduction Generating large, light quark dynamical gauge field
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More information19 : Slice Sampling and HMC
10-708: Probabilistic Graphical Models 10-708, Spring 2018 19 : Slice Sampling and HMC Lecturer: Kayhan Batmanghelich Scribes: Boxiang Lyu 1 MCMC (Auxiliary Variables Methods) In inference, we are often
More informationPostdoctoral Fellow, BCAM Basque Center for Applied Mathematics Research Group Modelling and Simulation in Life and Materials Sciences
Tijana Radivojević Curriculum Vitae Alameda Mazarredo 14, 48009 Bilbao, Spain (+34) 946 567 842 tradivojevic@bcamath.org bcamath.org/en/people/radivojevic Current position since November 2016 Postdoctoral
More informationChapter 29. Quantum Chaos
Chapter 29 Quantum Chaos What happens to a Hamiltonian system that for classical mechanics is chaotic when we include a nonzero h? There is no problem in principle to answering this question: given a classical
More informationBoltzmann-Gibbs Preserving Langevin Integrators
Boltzmann-Gibbs Preserving Langevin Integrators Nawaf Bou-Rabee Applied and Comp. Math., Caltech INI Workshop on Markov-Chain Monte-Carlo Methods March 28, 2008 4000 atom cluster simulation Governing Equations
More informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
More informationarxiv: v1 [math.pr] 14 Nov 2017
Acta Numerica (2018), pp. 1 92 c Cambridge University Press, 2018 doi:10.1017/s09624929 Printed in the United Kingdom arxiv:1711.05337v1 [math.pr] 14 Nov 2017 Geometric integrators and the Hamiltonian
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationIMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS
IMPLICIT MASS-MATRIX PENALIZATION OF HAMILTONIAN DYNAMICS WITH APPLICATION TO EXACT SAMPLING OF STIFF SYSTEMS PETR PLECHÁČ AND MATHIAS ROUSSET Résumé. An implicit mass-matrix penalization (IMMP) of Hamiltonian
More informationErnst Hairer 1, Robert I. McLachlan 2 and Robert D. Skeel 3
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ON ENERGY CONSERVATION OF THE SIMPLIFIED TAKAHASHI IMADA METHOD, Ernst Hairer 1,
More information17 : Optimization and Monte Carlo Methods
10-708: Probabilistic Graphical Models Spring 2017 17 : Optimization and Monte Carlo Methods Lecturer: Avinava Dubey Scribes: Neil Spencer, YJ Choe 1 Recap 1.1 Monte Carlo Monte Carlo methods such as rejection
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationMonte Carlo Algorithms
Last Modified Wednesday, 15 June 011 Monte Carlo Algorithms Tait Institute, School of Physics & Astronomy, The University of Edinburgh Tuesday, 14 June 011 STRONGnet Summer School 011, Bielefeld Outline
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationSecond quantization: where quantization and particles come from?
110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian
More informationApproximate inference in Energy-Based Models
CSC 2535: 2013 Lecture 3b Approximate inference in Energy-Based Models Geoffrey Hinton Two types of density model Stochastic generative model using directed acyclic graph (e.g. Bayes Net) Energy-based
More informationHamiltonian Monte Carlo for Scalable Deep Learning
Hamiltonian Monte Carlo for Scalable Deep Learning Isaac Robson Department of Statistics and Operations Research, University of North Carolina at Chapel Hill isrobson@email.unc.edu BIOS 740 May 4, 2018
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 17, 23) I. INTRODUCTION One frequently obtains detailed dynamical information
More informationClassical Equations of Motion
3 Classical Equations of Motion Several formulations are in use Newtonian Lagrangian Hamiltonian Advantages of non-newtonian formulations more general, no need for fictitious forces better suited for multiparticle
More informationChapter 2 Sampling Constrained Probability Distributions Using Spherical Augmentation
Chapter 2 Sampling Constrained Probability Distributions Using Spherical Augmentation Shiwei Lan and Babak Shahbaba Abstract Statistical models with constrained probability distributions are abundant in
More informationEach problem is worth 34 points. 1. Harmonic Oscillator Consider the Hamiltonian for a simple harmonic oscillator. 2ml 2 0. d 2
Physics 443 Prelim # with solutions March 7, 8 Each problem is worth 34 points.. Harmonic Oscillator Consider the Hamiltonian for a simple harmonic oscillator H p m + mω x (a Use dimensional analysis to
More informationA Noisy Monte Carlo Algorithm
A Noisy Monte Carlo Algorithm UK/99 06 May 1999 hep-lat/9905033 arxiv:hep-lat/9905033v2 12 May 2000 L. Lin, K.F. Liu, and J. Sloan Dept. of Physics and Astronomy, Univ. of Kentucky, Lexington, KY 40506
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationReinforcement Learning via Policy Optimization
Reinforcement Learning via Policy Optimization Hanxiao Liu November 22, 2017 1 / 27 Reinforcement Learning Policy a π(s) 2 / 27 Example - Mario 3 / 27 Example - ChatBot 4 / 27 Applications - Video Games
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationHigh-order actions and their applications to honor our friend and
High-order actions and their applications to honor our friend and collaborator Siu A. Chin Universitat Politecnica de Catalunya, Barcelona, Spain Lecture: Magnus Expansion and Suzuki s Method Outline of
More informationSmart walking: A new method for Boltzmann sampling of protein conformations
Smart walking: A new method for Boltzmann sampling of protein conformations Ruhong Zhou and B. J. Berne Department of Chemistry and the Center for Biomolecular Simulation, Columbia University, New York,
More informationBridging the Gap between Stochastic Gradient MCMC and Stochastic Optimization
Bridging the Gap between Stochastic Gradient MCMC and Stochastic Optimization Changyou Chen, David Carlson, Zhe Gan, Chunyuan Li, Lawrence Carin May 2, 2016 1 Changyou Chen Bridging the Gap between Stochastic
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationSimulations with MM Force Fields. Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi
Simulations with MM Force Fields Monte Carlo (MC) and Molecular Dynamics (MD) Video II.vi Some slides taken with permission from Howard R. Mayne Department of Chemistry University of New Hampshire Walking
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationANGULAR DEPENDENCE OF ELECTRON VELOCITY DISTRIBUTIONS IN LOW-PRESSURE INDUCTIVELY COUPLED PLASMAS 1
ANGULAR DEPENDENCE OF ELECTRON VELOCITY DISTRIBUTIONS IN LOW-PRESSURE INDUCTIVELY COUPLED PLASMAS 1 Alex V. Vasenkov 2, and Mark J. Kushner Department of Electrical and Computer Engineering Urbana, IL
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationE = φ 1 A The dynamics of a particle with mass m and charge q is determined by the Hamiltonian
Lecture 9 Relevant sections in text: 2.6 Charged particle in an electromagnetic field We now turn to another extremely important example of quantum dynamics. Let us describe a non-relativistic particle
More informationPaul Karapanagiotidis ECO4060
Paul Karapanagiotidis ECO4060 The way forward 1) Motivate why Markov-Chain Monte Carlo (MCMC) is useful for econometric modeling 2) Introduce Markov-Chain Monte Carlo (MCMC) - Metropolis-Hastings (MH)
More informationarxiv: v1 [stat.co] 17 Sep 2013
Spherical Hamiltonian Monte Carlo for Constrained Target Distributions Shiwei Lan, Bo Zhou, Babak Shahbaba September 18, 013 arxiv:1309.489v1 [stat.co] 17 Sep 013 Abstract We propose a new Markov Chain
More informationVan Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006)
Van Gunsteren et al. Angew. Chem. Int. Ed. 45, 4064 (2006) Martini Workshop 2015 Coarse Graining Basics Alex de Vries Every word or concept, clear as it may seem to be, has only a limited range of applicability
More informationChapter 1 Symplectic Integrator and Beam Dynamics Simulations
Chapter 1 and Beam Accelerator Physics Group, Journal Club November 9, 2010 and Beam NSLS-II Brookhaven National Laboratory 1.1 (70) Big Picture for Numerical Accelerator Physics and Beam For single particle
More informationAdvanced Molecular Dynamics
Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration
More information2 Canonical quantization
Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.
More informationMIT /30 Gelman, Carpenter, Hoffman, Guo, Goodrich, Lee,... Stan for Bayesian data analysis
MIT 1985 1/30 Stan: a program for Bayesian data analysis with complex models Andrew Gelman, Bob Carpenter, and Matt Hoffman, Jiqiang Guo, Ben Goodrich, and Daniel Lee Department of Statistics, Columbia
More informationA Lattice Path Integral for Supersymmetric Quantum Mechanics
Syracuse University SURFACE Physics College of Arts and Sciences 8-3-2000 A Lattice Path Integral for Supersymmetric Quantum Mechanics Simon Catterall Syracuse University Eric B. Gregory Zhongshan University
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationLearning Mixtures of Truncated Basis Functions from Data
Learning Mixtures of Truncated Basis Functions from Data Helge Langseth, Thomas D. Nielsen, and Antonio Salmerón PGM This work is supported by an Abel grant from Iceland, Liechtenstein, and Norway through
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationLagrangian Dynamical Monte Carlo
Lagrangian Dynamical Monte Carlo Shiwei Lan, Vassilios Stathopoulos, Babak Shahbaba, Mark Girolami November, arxiv:.3759v [stat.co] Nov Abstract Hamiltonian Monte Carlo (HMC) improves the computational
More informationObject-based High Contrast Traveltime Tomography
1 Object-based High Contrast Traveltime Tomography Yenting Lin Member, IEEE, Antonio Ortega, Fellow, IEEE arxiv:1303.3052v2 [physics.data-an] 6 Jul 2017 Abstract We consider a traveltime tomography problem
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises
More informationPreface On August 15, 2012, I received the following email message from Ryan Budney (ryan.budney@gmail.com). Hi Dick, Here s an MO post that s right up your alley. :) http://mathoverflow.net/questions/104750/about-a-letter-by-richard-palais-of-1965
More informationThe distance of the object from the equilibrium position is m.
Answers, Even-Numbered Problems, Chapter..4.6.8.0..4.6.8 (a) A = 0.0 m (b).60 s (c) 0.65 Hz Whenever the object is released from rest, its initial displacement equals the amplitude of its SHM. (a) so 0.065
More informationDraft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov
TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL Nikolai Shegunov, Ivan Hristov March 9, Seminar in IMI,BAS,Soa 1/25 Mathematical model We consider the Hamiltonian :
More information129 Lecture Notes More on Dirac Equation
19 Lecture Notes More on Dirac Equation 1 Ultra-relativistic Limit We have solved the Diraction in the Lecture Notes on Relativistic Quantum Mechanics, and saw that the upper lower two components are large
More informationBayesian Sampling Using Stochastic Gradient Thermostats
Bayesian Sampling Using Stochastic Gradient Thermostats Nan Ding Google Inc. dingnan@google.com Youhan Fang Purdue University yfang@cs.purdue.edu Ryan Babbush Google Inc. babbush@google.com Changyou Chen
More informationarxiv: v1 [physics.plasm-ph] 15 Sep 2013
Symplectic Integration of Magnetic Systems arxiv:309.373v [physics.plasm-ph] 5 Sep 03 Abstract Stephen D. Webb Tech-X Corporation, 56 Arapahoe Ave., Boulder, CO 80303 Dependable numerical results from
More informationPY 351 Modern Physics - Lecture notes, 3
PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence
More informationQuantum Statistical Calculations and Symplectic Corrector Algorithms
Quantum Statistical Calculations and Symplectic Corrector Algorithms Siu A. Chin Department of Physics, Texas A&M University, College Station, TX 77843, USA arxiv:cond-mat/03202v2 [cond-mat.stat-mech]
More informationManifold Monte Carlo Methods
Manifold Monte Carlo Methods Mark Girolami Department of Statistical Science University College London Joint work with Ben Calderhead Research Section Ordinary Meeting The Royal Statistical Society October
More informationLow Emittance Machines
CERN Accelerator School Advanced Accelerator Physics Course Trondheim, Norway, August 2013 Low Emittance Machines Part 1: Beam Dynamics with Synchrotron Radiation Andy Wolski The Cockcroft Institute, and
More informationReinforcement Learning
Reinforcement Learning Policy gradients Daniel Hennes 26.06.2017 University Stuttgart - IPVS - Machine Learning & Robotics 1 Policy based reinforcement learning So far we approximated the action value
More information2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
More informationLecture 7 and 8: Markov Chain Monte Carlo
Lecture 7 and 8: Markov Chain Monte Carlo 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering University of Cambridge http://mlg.eng.cam.ac.uk/teaching/4f13/ Ghahramani
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More information( r) = 1 Z. e Zr/a 0. + n +1δ n', n+1 ). dt ' e i ( ε n ε i )t'/! a n ( t) = n ψ t = 1 i! e iε n t/! n' x n = Physics 624, Quantum II -- Exam 1
Physics 624, Quantum II -- Exam 1 Please show all your work on the separate sheets provided (and be sure to include your name) You are graded on your work on those pages, with partial credit where it is
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationarxiv: v1 [stat.ml] 23 Jan 2019
arxiv:1901.08045v1 [stat.ml] 3 Jan 019 Abstract We consider the problem of sampling from posterior distributions for Bayesian models where some parameters are restricted to be orthogonal matrices. Such
More informationDimension-Independent likelihood-informed (DILI) MCMC
Dimension-Independent likelihood-informed (DILI) MCMC Tiangang Cui, Kody Law 2, Youssef Marzouk Massachusetts Institute of Technology 2 Oak Ridge National Laboratory 2 August 25 TC, KL, YM DILI MCMC USC
More information