Enhanced sampling via molecular dynamics II: Unbiased approaches
|
|
- Lester Kelly
- 6 years ago
- Views:
Transcription
1 Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY 10003
2 Spatial-Warping Transformations Z. Zhu, MET, S. O. Samuelson, G. J. Martyna, Phys. Rev. Lett (2002) P. Minary, MET, G. J. Martyna SIAM J. Sci. Comput. 30, 2055 (2007) Science, December 22, 2006
3 Molecular dynamics in a double well. V(x) p x = m dv p = = F( x) + heat bath dx x time Crossing the barrier is a rare event!
4 Spatial-warping transformations Z. Zhu, MET, S. O. Samuelson, G. J. Martyna, Phys. Rev. Lett (2002) V(x) Hamiltonian: 2 p H( x, p) = + V( x) 2m Canonical partition function: 2 p x= a x= a Q = dp dx exp β + V ( x) 2m Change the integration variables: u = f ( x) = a + dy e β a Transformed partition function: x V ( y) 1 βv ( f ( u)) dx = e du 2 p Q = dp du exp β + V f ( u) V f ( u) 2m ( 1 ) ( 1 )
5 V ( x ) Spatial-warping transformations (cont d) Choosing the transformation: V ( x) V( x) a< x< a = 0 Otherwise x = a x= a V(x) V ( x) V ( x ) V ( x ) 1 1 ( ( )) ( ( )) V f u V f u MD: p d u = p = Fu ( ) + heat bath F( u) = V f ( u) V f ( u) m du Reference potential spatial warping (REPSWA) 1 1 ( ) ( )
6 How is the space warped? Barrier: kt B Barrier: 5kT B Barrier: 10kT B
7 V No Transformation Transformation 5kT 10kT 10kT
8 How to do biomolecules and polymers P. Minary, MET, G. J. Martyna SIAM J. Sci. Comput. 30, 2055 (2007) Sampling hindered by barriers in the dihedral angles: V ( φ) Ramachandran dihedral angles (alanine dipeptide): φ ψ ψ φ
9 Dihedral angle transformations Transformed angle: φ = + d e u min φ φ φ χ min βv ( χ) V ( φ) φ min
10 P. Minary, G. J. Martyna and MET SIAM J. Sci. Comp. 30, 2055 (2008) r 6 V ( φ,{ r}) = V ( φ) S ( φ) + α V ( r ( φ,{ r}) r ) S ( r ( φ,{ r}) ri ) tors 1 n.b. 4 i 2 4 i nn
11
12 A 400-mer alkane chain No nonbonded interactions 397 collective variables! Nonbonded interactions off
13 No Transformation REPSWA Transformations RIS Model value: 10
14 β-barrel Proteins
15 Honeycutt and Thirumalai PNAS 87, 3526 (1990)
16 MD REMC Transformations PT replicas = 16
17 Equations of motion: Integrators in molecular dynamics p q =, p = Fq ( ) m Time evolution: q( t) (0) (0) il t q q = e = φ t p( t) p(0) p(0) where il p = + F( q) m q p t p t = + 2 p m q 2 p 3 exp( ilt) exp Fq ( ) exp t exp Fq ( ) O( t) Pseudocode: p p+ 0.5* t* F; q q+ t* p; GetForce( q); p p+ 0.5* t* F;
18 Multiple time scale (r-respa) integration MET, G. J. Martyna and B. J. Berne, J. Chem. Phys. 97, 1990 (1992) p q = m p = F + F fast slow p il = + F + F = il + il m q p p fast slow ref slow
19 Multiple Time Scale Integration (cont d) Trotter factorized propagator: 3 exp( il t) = exp( ilslow t / 2)exp( ilref t)exp( ilslow t / 2) + O( t ) = exp( ilslow t / 2) exp( ilref δt) exp( ilslow t / 2) n t exp( il t) = exp Fslow 2 p δt p δt exp Ffast exp δt exp Ffast 2 p m q 2 p t exp F 2 slow p n
20 From: G. J. Martyna, MET, D. Tobias, and M. L. Klein Mol. Phys. 87, 1117 (1996).
21 Illustration of resonance A. Sandu and T. S. Schlick, J. Comput. Phys. 151, 74 (1999) F = ω x F = Ω x fast 2 2 slow t 2 p = p(0) Ω x(0) 2 p ' x( t) = x(0)cos( ω t) + sin( ω t) ω p = p(0)cos( ω t) ωx(0)sin( ω t) t p t = p Ω x t 2 2 ( ) ( ) x( t) x(0) = A( ω, Ω, t) p( t) p(0)
22 Illustration of resonance (cont d) 2 Ω t 1 cos( ω t) sin( ω t) sin( ω t) 2ω ω A( ω, Ω, t) = t Ω 2 Ω t ω sin( ω t) tω cos( ω t) cos( ω t) sin( ω t) 4ω 2ω Note: det(a) = 1 Depending on Δt, eigenvalues of A are either complex conjugate pairs 2 < Tr( A) < 2 or eigenvalues are both real Tr( A) 2 Leads to resonances ( Tr(A) 2) at Δt = nπ/ω
23 Isokinetic dynamics Constraint the kinetic energy of a system: 2 i 3N 1 p = kt 2m 2 Introduce constraint via a Lagrange multiplier: i i p r = p = F λp i i i i i m i Derivative of constraint yields multiplier: Fp i i / mi pi pi i p i = [ Fi λpi] = 0 λ = m m p p / m Partition function generated: i i i i i i i i p 3N 1 Ω= d p δ kt d r e i 2mi 2 Resonance-free but is a poor thermostat. 2 N i N βu( r)
24 Consider the equations of motion: Resonance-free stochastic isokinetic thermostat P. Minary, MET, G. J. Martyna PRL (2004); B. Leimkuhler, D. Margul, MET Mol. Phys. (in press) dq = vdt F dv = λv dt m dv = λv dt v v dt, k = 1,..., L 1, k 1, k 2, k 1, k Gv ( ) dv = dt γv dt + σdw 1, k 2, k 2, k k Q2 G( u) Isokinetic constraint: Multiplier: = = Qu β k= 1 λ L mv Q v v v v L L , k = Λ (, 1,1,..., 1, L) = β L + 1 k= 1 L L vf Q v v L + 1 Λ( vv,,..., v ) 1,1 1, L 2 1 1, k 2, k
25 Multiple time-step integrator Assume forces have fast and slow components: F = Ff + Fs il q = v q il = il + il + il + il + il ( f ) ( s) q v v N OU il F v v L ( f ) f ( f ) ( f ) v = λf F 1, k m λ v k = 1 v1, k F il v v L ( s) s ( s) ( s) v = λf λf 1, k m v k = 1 v1, k Gv ( ) il v v v v L L L 1, k N = λn λn 1, k 2, k 1, k + v k= 1 v1, k k= 1 v1, k k= 1 Q2 v2, k L L Qv v ( f ) vff ( s) vf s L + 1 k= 1 λf =, λf =, λ N = Λ Λ Λ 2 1 1, k 2, k
26 Multiple time-step integrator ( s) ( f ) /2 /2 ( f ) n ( s) il t iln t/2 ilv t/2 ilv δt/2 ilqδt il il OUδt qδt ilv δt/2 ilv t/2 iln t/2 e = e e e e e e e e e Solution of the Ornstein-Uhlenbeck process: e v = v (0) e + σ Rt ( ) ilouδ t γδ t 2, k 2, k 1 e 2γ 2γδ t
27 Numerical illustration of resonance Harmonic oscillator with quartic perturbation L 9 Uq ( ) = q q Nbins t 1 = 3, δt = ς () t = [ Pq ( i;) t Pexact ( qi) ] 100 N i= 1 bins P(q) q
28 Flexible SPC/E water U() r = U () r + U () r = U () r + U () r + U () r intra nb intra short long qq i jerfc( αr) ij Ushort () r = ULJ ( rij ) + ij rij 1 4π 2 2 α U ( r) = e S( k) q long 2 V kk, < k k S( k) = RESPA schemes i qe i ikr F. Paesani et al. J. Chem. Phys. 125, (2006) max k /4α 2 2 i π i qq U () r = U (; r α, k ) + 1 δ θ( r r ) ( ) ij (1) i j ij long recip res ij cut i> j rij ( ) ij erf ( αr) qq U () r = U (; r α, k ) 1 δ 1 θ( r r ) erf ( αr) (2) i j ij long recip max ij cut i> j rij
29 δ t = 0.5 fs t = 3 fs T =? RESPA 1: L = 4, r cut = 6 Å RESPA 2: L = 4, r cut = 8 Å
30 RESPA 1: L = 4, r cut = 6 Å L = 4 (colors) L = 1 (dashed) L = 2 (circles)
31 Parallel Tempering Swendsen and Wang, Phys. Rev. Lett. 57, 2607 (1986) T M TM 1 T 1 f r ( K ) K ( ) = ( K ) βk U ( r ) e Z( NVT,, ) K F M (1) ( M) ( K) ( r,..., r ) f( r ) = K = 1 ( K+ 1) ( K) ( K+ 1) ( K) ( K) ( K+ 1) fk( r ) fk+ 1( r ) A( r, r r, r ) = min 1, = min 1, e fk( r ) fk+ 1( r ) K, K+ 1 ( K) ( K+ 1) ( β β ) KK, + 1 = K K+ 1 U r ( K) ( K+ 1) ( ) U( ) r
32 Comparison for a 50-mer chain using CHARMM22 REMC/PT replicas = 10; 300 < T < 1000, 10 replicas 20% acceptance Nonbonded interactions off
33 Comparison for 50-mer using CHARMM22 all interactions
34 Comparison for 50-mer with all interactions PT replicas = 10; REPSWA α = 0.8; Every 10 th dihedral not transformed
35 Replica exchange with solute tempering P. Liu, B. Kim, R. A. Friesner, B. J. Berne PNAS (2005); L. Wang, R. A. Friesner, B. J. Berne J. Phys. Chem. B (2011) A downside of replica-exchange is that U ~ N, therefore, the energy differences between neighboring replica need to shrink as the system size grows, which means more and more replica are needed in order for the acceptance probability to remain approximately fixed. Replica exchange with solute tempering (REST) was developed to allow just part of the system (the solute ) to be subject to attempted exchanges. It relies on the idea that different replicas can be subject to different potentials. f K () r = βkuk( r) e Z( NVT,, ) ( ( K ) ( K 1) ) ( ( K ( ) ( ) ( 1) ) ( ( K U r U r + β U r + U r ) )) = β + KK, + 1 K K K K+ 1 K+ 1 K+ 1 If the potential energy of a solute-solvent system is: U () r = U () r + U () r + U () r K s sv v The latest version of REST (called REST2) uses the following ladder of potentials: U () β β r = U () r + U () r + U () r K K K s sv v β0 β0 K
36 The ladder of potentials Replica exchange with solute tempering allows each replica to be run at a single temperature T 0 so that the solute probability e = e β ( β / β ) U ( r) β U ( r) 0 K 0 s K s For full potential ladder: βu () r = βu () r + ββ U () r + βu () r 0 K K s 0 K sv 0 v Acceptance probability determined by β K, K + 1 = K K + 1 s s + sv βk + βk+ 1 ( K) ( K+ 1) ( ) ( 1) ( ) ( ) ( ) 0 K K+ β β U r U r ( U ( r ) U( r ))
37 Solute potential energy Convergence of alanine dipeptide in solution
Free energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationNewtonian and extended Lagrangian dynamics
Newtonian and extended Lagrangian dynamics Gianni Cardini and Riccardo Chelli Dipartimento di Chimica Ugo Schiff Università di Firenze, Via della Lastruccia 3, 50019 Sesto Fno, Firenze giannicardini@unifiit,
More informationMolecular Dynamics Lecture 3
Molecular Dynamics Lecture 3 Ben Leimkuhler the problem of the timestep in MD constraints - SHAKE/RATTLE multiple timestepping stochastic methods for holonomic constraints stochastic multiple timestepping
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationMethods of Computer Simulation. Molecular Dynamics and Monte Carlo
Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?
More informationLow-temperature isomers of the water hexamer
Low-temperature isomers of the water hexamer Volodymyr Babin and Francesco Paesani Department of Chemistry & Biochemistry, University of California, San Diego, La Jolla, CA 92093-0314 Hundreds of years
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationSampling the conformation space of complex systems, such as
Replica exchange with solute tempering: A method for sampling biological systems in explicit water Pu Liu*, Byungchan Kim*, Richard A. Friesner, and B. J. Berne Department of Chemistry and Center for Biomolecular
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationEXTENDING THE BOUNDARIES OF STRUCTURAL MODELING EPIGENETIC EFFECTS ON NUCLEOSOME POSITIONING
EXTENDING THE BOUNDARIES OF STRUCTURAL MODELING EPIGENETIC EFFECTS ON NUCLEOSOME POSITIONING Peter Minary Computational Structural Biology Group Stanford University Stanford, CA 94305 ~ 0.3 nm TO MODEL
More informationSet the initial conditions r i. Update neighborlist. r i. Get new forces F i
Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling
More informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
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 informationRare Event Simulations
Rare Event Simulations Homogeneous nucleation is a rare event (e.g. Liquid Solid) Crystallization requires supercooling (µ solid < µ liquid ) Crystal nucleus 2r Free - energy gain r 3 4! 3! GBulk = ñr!
More informationEfficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers
Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics
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 informationApproach to Thermal Equilibrium in Biomolecular
Approach to Thermal Equilibrium in Biomolecular Simulation Eric Barth 1, Ben Leimkuhler 2, and Chris Sweet 2 1 Department of Mathematics Kalamazoo College Kalamazoo, Michigan, USA 49006 2 Centre for Mathematical
More informationLaser Induced Control of Condensed Phase Electron Transfer
Laser Induced Control of Condensed Phase Electron Transfer Rob D. Coalson, Dept. of Chemistry, Univ. of Pittsburgh Yuri Dakhnovskii, Dept. of Physics, Univ. of Wyoming Deborah G. Evans, Dept. of Chemistry,
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationMolecular Dynamics Lecture 2
Molecular Dynamics Lecture 2 Ben Leimkuhler the problem of the timestep in md splitting methods in Langevin dynamics constraints multiple timestepping Spring School, WIAS, Berlin 2017 fruit fly alanine
More informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
More informationOpen quantum systems
Open quantum systems Wikipedia: An open quantum system is a quantum system which is found to be in interaction with an external quantum system, the environment. The open quantum system can be viewed as
More informationAb initio molecular dynamics
Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution
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 informationA Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space
Proceedings of Institute of Mathematics of NAS of Ukraine 00, Vol. 43, Part, 645 65 A Symmetric Treatment of Damped Harmonic Oscillator in Extended Phase Space S. NASIRI and H. SAFARI Institute for Advanced
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 informationM e t ir c S p a c es
A G M A A q D q O I q 4 78 q q G q 3 q v- q A G q M A G M 3 5 4 A D O I A 4 78 / 3 v D OI A G M 3 4 78 / 3 54 D D v M q D M 3 v A G M 3 v M 3 5 A 4 M W q x - - - v Z M * A D q q q v W q q q q D q q W q
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationElectronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani
Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations
More informationOvercoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 5 1 AUGUST 1998 Overcoming stability limitations in biomolecular dynamics. I. Combining force splitting via extrapolation with Langevin dynamics in LN Eric
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationLuttinger Liquid at the Edge of a Graphene Vacuum
Luttinger Liquid at the Edge of a Graphene Vacuum H.A. Fertig, Indiana University Luis Brey, CSIC, Madrid I. Introduction: Graphene Edge States (Non-Interacting) II. III. Quantum Hall Ferromagnetism and
More informationSupporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Dipak
More informationQuantum Physics III (8.06) Spring 2008 Final Exam Solutions
Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,
More informationCompeting sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit
Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit Paul Dupuis Division of Applied Mathematics Brown University IPAM (J. Doll, M. Snarski,
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationAspects of nonautonomous molecular dynamics
Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski
More informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationReplica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)
pubs.acs.org/jpcb Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) Lingle Wang, Richard A. Friesner, and B. J. Berne* Department of Chemistry,
More informationSpeeding up path integral simulations
Speeding up path integral simulations Thomas Markland and David Manolopoulos Department of Chemistry University of Oxford Funded by the US Office of Naval Research and the UK EPSRC Outline 1. Ring polymer
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 0 Simple Biasing Methods David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo simulation Markov chain to generate
More informationLecture 6: 2D FET Electrostatics
Lecture 6: 2D FET Electrostatics 2016-02-01 Lecture 6, High Speed Devices 2014 1 Lecture 6: III-V FET DC I - MESFETs Reading Guide: Liu: 323-337 (he mainly focuses on the single heterostructure FET) Jena:
More informationThe Geometry of Euler s equation. Introduction
The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated
More information3. Riley 12.9: The equation sin x dy +2ycos x 1 dx can be reduced to a quadrature by the standard integrating factor,» Z x f(x) exp 2 dt cos t exp (2
PHYS 725 HW #4. Due 15 November 21 1. Riley 12.3: R dq dt + q C V (t); The solution is obtained with the integrating factor exp (t/rc), giving q(t) e t/rc 1 R dsv (s) e s/rc + q() With q() and V (t) V
More informationPractice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009.
Practice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009. Solutions (1) Let Γ be a discrete group acting on a manifold M. (a) Define what it means for Γ to act freely. Solution: Γ acts
More informationLecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The Sine-Gordon Equation 8. Wednesday, August 28
Lecture 2 Wednesday, August 28 Contents 1 Fermi s Method 2 2 Lattice Oscillators 3 3 The Sine-Gordon Equation 8 1 1 Fermi s Method Feynman s Quantum Electrodynamics refers on the first page of the first
More informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationNosé-Hoover Thermostats
Nosé- Nosé- Texas A & M, UT Austin, St. Olaf College October 29, 2013 Physical Model Nosé- q Figure: Simple Oscillator p Physical Model Continued Heat Bath T p Nosé- q Figure: Simple Oscillator in a heat
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationFluctuation theorem in systems in contact with different heath baths: theory and experiments.
Fluctuation theorem in systems in contact with different heath baths: theory and experiments. Alberto Imparato Institut for Fysik og Astronomi Aarhus Universitet Denmark Workshop Advances in Nonequilibrium
More informationSupplemental Material for Global Langevin model of multidimensional biomolecular dynamics
Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics,
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationLAMMPS Performance Benchmark on VSC-1 and VSC-2
LAMMPS Performance Benchmark on VSC-1 and VSC-2 Daniel Tunega and Roland Šolc Institute of Soil Research, University of Natural Resources and Life Sciences VSC meeting, Neusiedl am See, February 27-28,
More information4.3 Lecture 18: Quantum Mechanics
CHAPTER 4. QUANTUM SYSTEMS 73 4.3 Lecture 18: Quantum Mechanics 4.3.1 Basics Now that we have mathematical tools of linear algebra we are ready to develop a framework of quantum mechanics. The framework
More informationarxiv:hep-th/ v3 16 May 1996
BNL-63106 An Exact Solution for Quantum Tunneling in a Dissipative System arxiv:hep-th/9605081v3 16 May 1996 Li Hua Yu National Synchrotron Light Source, Brookhaven National Laboratory, N.Y.11973 Abstract
More informationSupporting information
Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationPopulation annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice
Population annealing study of the frustrated Ising antiferromagnet on the stacked triangular lattice Michal Borovský Department of Theoretical Physics and Astrophysics, University of P. J. Šafárik in Košice,
More informationIonic Liquids simulations : obtention of structural and transport properties from molecular dynamics. C. J. F. Solano, D. Beljonne, R.
Ionic Liquids simulations : obtention of structural and transport properties from molecular dynamics C. J. F. Solano, D. Beljonne, R. Lazzaroni Ionic Liquids simulations : obtention of structural and transport
More information5.4 Given the basis e 1, e 2 write the matrices that represent the unitary transformations corresponding to the following changes of basis:
5 Representations 5.3 Given a three-dimensional Hilbert space, consider the two observables ξ and η that, with respect to the basis 1, 2, 3, arerepresentedby the matrices: ξ ξ 1 0 0 0 ξ 1 0 0 0 ξ 3, ξ
More informationSerial Generalized Ensemble Simulations of Biomolecules with Self-Consistent Determination of Weights
pubs.acs.org/jctc Serial Generalized Ensemble Simulations of Biomolecules with Self-Consistent Determination of Weights Riccardo Chelli*,, and Giorgio F. Signorini Dipartimento di Chimica, Universita di
More informationEffective Temperatures in Driven Systems near Jamming
Effective Temperatures in Driven Systems near Jamming Andrea J. Liu Department of Physics & Astronomy University of Pennsylvania Tom Haxton Yair Shokef Tal Danino Ian Ono Corey S. O Hern Douglas Durian
More information1 The Quantum Anharmonic Oscillator
1 The Quantum Anharmonic Oscillator Perturbation theory based on Feynman diagrams can be used to calculate observables in Quantum Electrodynamics, like the anomalous magnetic moment of the electron, and
More informationTwo recent works on molecular systems out of equilibrium
Two recent works on molecular systems out of equilibrium Frédéric Legoll ENPC and INRIA joint work with M. Dobson, T. Lelièvre, G. Stoltz (ENPC and INRIA), A. Iacobucci and S. Olla (Dauphine). CECAM workshop:
More informationLecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review
Lecture contents A few concepts from Quantum Mechanics Particle in a well Two wells: QM perturbation theory Many wells (atoms) BAND formation Tight-binding model Solid state physics review Approximations
More informationarxiv: v1 [physics.chem-ph] 24 Apr 2018
arxiv:1804.08913v1 [physics.chem-ph] 24 Apr 2018 Fast-Forward Langevin Dynamics with Momentum Flips Mahdi Hijazi, 1 David M. Wilkins, 1, a) and Michele Ceriotti 1 Laboratory of Computational Science and
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationComputational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationSerial Replica Exchange
1416 J. Phys. Chem. B 2007, 111, 1416-1423 Serial Replica Exchange Morten Hagen, Byungchan Kim, Pu Liu, Richard A. Friesner, and B. J. Berne* Department of Chemistry and Center for Biomolecular Simulation,
More informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
More informationEC Control Engineering Quiz II IIT Madras
EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ
More informationProlongation structure for nonlinear integrable couplings of a KdV soliton hierarchy
Prolongation structure for nonlinear integrable couplings of a KdV soliton hierarchy Yu Fa-Jun School of Mathematics and Systematic Sciences, Shenyang Normal University, Shenyang 110034, China Received
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationLINEAR ALGEBRA (PMTH213) Tutorial Questions
Tutorial Questions The tutorial exercises range from revision and routine practice, through lling in details in the notes, to applications of the theory. While the tutorial problems are not compulsory,
More informationMolecular dynamics with multiple time scales: The selection of efficient reference system propagators
Molecular dynamics with multiple time scales: The selection of efficient reference system propagators Steven J. Stuart, Ruhong Zhou, and B. J. Berne Department of Chemistry and the Center for Biomolecular
More informationPhys460.nb Back to our example. on the same quantum state. i.e., if we have initial condition (5.241) ψ(t = 0) = χ n (t = 0)
Phys46.nb 89 on the same quantum state. i.e., if we have initial condition ψ(t ) χ n (t ) (5.41) then at later time ψ(t) e i ϕ(t) χ n (t) (5.4) This phase ϕ contains two parts ϕ(t) - E n(t) t + ϕ B (t)
More informationPhonons (Classical theory)
Phonons (Classical theory) (Read Kittel ch. 4) Classical theory. Consider propagation of elastic waves in cubic crystal, along [00], [0], or [] directions. Entire plane vibrates in phase in these directions
More informationTime-dependent DMRG:
The time-dependent DMRG and its applications Adrian Feiguin Time-dependent DMRG: ^ ^ ih Ψ( t) = 0 t t [ H ( t) E ] Ψ( )... In a truncated basis: t=3 τ t=4 τ t=5τ t=2 τ t= τ t=0 Hilbert space S.R.White
More informationof Long-Term Transport
Atomistic Modeling and Simulation of Long-Term Transport Phenomena in Nanomaterials K.G. Wang, M. and M. Ortiz Virginia Polytechnic Institute and State University Universidad de Sevilla California Institute
More informationG : Statistical Mechanics
G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the
More informationOptimized statistical ensembles for slowly equilibrating classical and quantum systems
Optimized statistical ensembles for slowly equilibrating classical and quantum systems IPAM, January 2009 Simon Trebst Microsoft Station Q University of California, Santa Barbara Collaborators: David Huse,
More information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
More informationnd A L T O SOLO LOWELL. MICHIGAN. THURSDAY. APRIL Spring Activities Heads Up and Forward (Editorial By " T h e Committee'')
- 6 7 8 9 3-6 7 8 9 3 G UDY 3 93 VU XXXV U XY F K FD D j V K D V FY G F D Y X K X DD Y j \ V F \ VD GD D U Y 78 K U D Y U Y 484?35 V 93 7 4 U x K 77 - D :3 F K > 6 D x F 5 - - x - G 7 43 8 D $35 K F $5
More informationLangevin Dynamics of a Single Particle
Overview Langevin Dynamics of a Single Particle We consider a spherical particle of radius r immersed in a viscous fluid and suppose that the dynamics of the particle depend on two forces: A drag force
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationConsider a system of n ODEs. parameter t, periodic with period T, Let Φ t to be the fundamental matrix of this system, satisfying the following:
Consider a system of n ODEs d ψ t = A t ψ t, dt ψ: R M n 1 C where A: R M n n C is a continuous, (n n) matrix-valued function of real parameter t, periodic with period T, t R k Z A t + kt = A t. Let Φ
More informationEran Rabani, S. A. Egorov, a) and B. J. Berne Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 15 15 OCTOBER 1998 A comparison of exact quantum mechanical and various semiclassical treatments for the vibronic absorption spectrum: The case of fast vibrational
More informationBurgers equation in the complex plane. Govind Menon Division of Applied Mathematics Brown University
Burgers equation in the complex plane Govind Menon Division of Applied Mathematics Brown University What this talk contains Interesting instances of the appearance of Burgers equation in the complex plane
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
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 information