Ab initio Molecular Dynamics
|
|
- Lester Hensley
- 5 years ago
- Views:
Transcription
1 Ab initio Molecular Dynamics Jürg Hutter Department of Chemistry, University of Zurich
2 Overview Equations of motion and Integrators General Lagrangian for AIMD Car Parrinello MD and Born Oppenheimer MD ab initio Langevin MD Stability and efficiency Thermostats External Drivers
3 Equations of Motion (EOM) Newton s EOM for a set of classical point particles in a potential. dv (R) M I RI = dr I These EOM generate for a given number of particles N in a volume V the micro canonical ensemble (NVE ensemble). The total energy E is a constant of motion!
4 Total Energy Total Energy = Kinetic energy + Potential energy Kinetic energy = T (Ṙ) = N Potential energy = V (R) I=1 M I 2 Ṙ2 We will use the total energy as an indicator for the numerical accuracy of simulations.
5 Lagrange Equation with d L dt Ṙ = L R L(R, Ṙ) = T (Ṙ) V (R) Equivalent to Newton s EOM in Cartesian coordinates, but is more general and flexible. Extended systems Constraints
6 Integration of EOM Discretization of time R(t) R(t + τ) R(t + 2τ) R(t + mτ) V (t) V (t + τ) V (t + 2τ) V (t + mτ) Efficiency: minimal number of force evaluations, minimal number of stored quantities Stability: minimal drift in constant of motion (energy) Accuracy: minimal distance to exact trajectory
7 Sources of Errors Type of integrator predictor-corrector, time-reversible, symplectic Time step τ short time accuracy measured as O(τ n ) Consistency of forces and energy e.g. cutoffs leading to non-smooth energy surfaces Accuracy of forces e.g. convergence of iterative force calculations (SCF, constraints)
8 Velocity Verlet Integrator R(t + τ) = R(t) + τv (t) + τ 2 V (t + τ) = V (t) + 2M f (t) τ [f (t) + f (t + τ) 2M Efficiency: 1 force evaluation, 3 storage vectors Stability: time reversible Accuracy: O(τ 2 ) Simple adaptation for constraints (shake, rattle, roll) Simple adaptation for multiple time steps and thermostats
9 Test on Required Accuracy of Forces Classical Force Field Calculations, 64 molecules, 330 K TIP3P (flexible), SPME (α = 0.44, GMAX = 25),
10 Stability: Accuracy of Forces Stdev. f Stdev. Energy Drift Hartree/Bohr µhartree µhartree/ns Kelvin/ns
11 Born Oppenheimer MD: The Easy Way dv (R) M I RI = dr I EOM V (R) = min KS({Φ(r)}; R) + const.] Φ Kohn Sham BO potential Forces f KS (R) = d min Φ E KS ({Φ(r)}; R) dr = E KS R + const. R + i (E KS + const.) Φ i } Φ {{ i R } =0
12 Computational Details System 64 water molecules density 1gcm 3 Temperature 330K Timestep 0.5fs DFT Calculations GPW, TZV2P basis (2560 bsf), PBE functional Cutoff 280 Rydberg, ɛ default = OT-DIIS, Preconditioner FULL_SINGLE_INVERSE Reference trajectory (1ps), ɛ SCF = 10 10
13 Stability in BOMD Unbiased initial guess; Φ(t) = Φ 0 (R(t)) ɛ SCF MAE E KS MAE f Drift Hartree Hartree/Bohr Kelvin/ns Consistent with results from classical MD Note accuracy of forces!
14 Efficiency: Initial Guess of Wavefunction 4th order Gear predictor (PS extrapolation in CP2K) Method ɛ SCF Iterations Drift (Kelvin/ns) Guess Gear(4) Gear(4) Gear(4) What is the problem? Time reversibility has been broken!
15 Generalized Lagrangian L(q, q, x, ẋ) = 1 2 M q µẋ2 E(q, y) + kµg( x y ) y = F (q, x) G( x y ) wavefunction optimization wavefunction retention potential Equations of motion M q = E q E y µẍ = E y F x + kµ F q + kµ G F y q [ G x + G y F x Extension of Niklasson Lagrangian, PRL (2008) ]
16 Dynamical System
17 Car Parrinello Molecular Dynamics y = x G( x y ) = 0 Lagrangian L(q, q, x, ẋ) = 1 2 M q µẋ2 E(q, x) Equations of motion M q = E q µẍ = E x
18 Properties of CPMD Accuracy: Medium Distance from BO surface controlled by mass µ Requires renormalization of dynamic quantities (e.g. vibrational spectra) Stability: Excellent All forces can be calculated to machine precision easily Efficiency: Good Efficiency is strongly system dependent (electronic gap) Requires many nuclear gradient calculations Not implemented in CP2K!
19 BOMD y = Min x E(q, x) and µ = 0 Lagrangian L(q, q) = 1 2 M q 2 + E(q, y) L(x, ẋ) = 1 2ẋ2 + kg( x y ) M q ẍ Equations of motion = E q decoupled equations = k G x
20 BOMD with Incomplete SCF Convergence y = F(q, x) Min x E(q, x) and µ = 0 Equations of motion M q = E q E F y q ẍ = k G x G F y x SCF Error: Neglect of force terms E F y q G F and y x. EOM are coupled through terms neglected!
21 ASPC Integrator Integration of electronic DOF (x) has to be accurate: good wavefunction guess gives improved efficiency stable: do not destroy time-reversibility of nuclear trajectory ASPC: Always Stable Predictor Corrector J. Kolafa, J. Comput Chem. 25: (2004) ASPC(k): time-reversible to order 2k + 1
22 Orthogonality Constraint Wavefunction extrapolation for non-orthogonal basis sets C init = K B j+1 C t jτ C T t jτ S t jτ }{{} PS j=0 C t τ Coefficients B i are given by ASPC algorithm x = C init y[t] = C t S overlap matrix J. VandeVondele et al. Comp. Phys. Comm. 167: (2005)
23 Importance of Time-Reversibility Method ɛ SCF Iterations Drift (Kelvin/ns) Guess ASPC(3) ASPC(3) Gear(4) Gear(4) Gear(4)
24 Efficiency and Drift Method ɛ SCF Iterations Drift (Kelvin/ns) Guess ASPC(4) ASPC(4) ASPC(4) ASPC(4)
25 Efficiency and Drift Method ɛ SCF Iterations Drift (Kelvin/ns) ASPC(4) ASPC(5) ASPC(6) ASPC(7) ASPC(8)
26 Langevin BOMD Starting point: BOMD with ASPC(k) extrapolation Analysis of forces f BO (R) = f HF (R) + f Pulay (R) +f nsc (R) }{{} f (R) f BO : correct BO force f HF : Hellmann Feynman force f Pulay : Pulay force f nsc : non-self consistency error force f : approximate BO force
27 Forces in Approximate BOMD Approximate f nsc by fnsc = [( Vxc (ρ i ) ] ) dr ρ i ρ + V H ( ρ) I ρ i with ρ = ρ o ρ i, ρ o final (output) density, ρ i initial (predicted) density. Now assume f (R) + f nsc (R) = f BO (R) γ D Ṙ where γ D is a constant friction parameter.
28 Langevin EOM M R = f BO (R) (γ D + γ L )Ṙ + Θ with Θ a Gaussian random noise term and Θ(0)Θ(t) = 6(γ D + γ L )Mk B T δ(t) Given temperature T = k B 3 MṘ2 and an arbitrary γ L this determines γ D
29 Langevin BOMD 2nd Generation Car Parrinello (SGCP) T. D. Kühne et al., Phys. Rev. Lett (2007) single SCF step plus force correction needed extremely efficient for systems with slow SCF convergence γ D is small: correct statistics and dynamics difficult to stabilize in complex systems
30 SGCP: Examples T. Musso et al. Eur. Phys. J. B, (2018) Analysis of Force Errors Liquid Silicon at 3000 K
31 SGCP: Nanomesh h-bn on Rh(111)
32 SGCP: Force Errors
33 SGCP: h-bn on Rh(111) Comparision of MD Results
34 SGCP: h-bn on Rh(111) Performance
35
36 Micro-Canonical Ensemble System variables: N, V, E Microstate definition: Γ = {r N, p N } Energy H(Γ) = i p 2 i 2m i + U(r N ) All states have the same weight δ(h(γ) E) Partitionfunction Ω(N, V, E) = Γ δ(h(γ) E) = 1 N 1 h 3N dr N dp N δ(h(γ) E)
37 Molecular Dynamics Fixed number of particles: N Fixed volume: V Periodic boundary condition to avoid surface effects Total energy E(t) = i m i 2 ṙ 2 i (t) + U(r(t))
38 Constant of Motion: E de(t) dt = i F i = du(r) dr i Ė = i m i 2 2ṙ i r i + i = m i r i m i ṙ i F i m i + i du(r) ṙ i dr i ( F i )ṙ i = i (ṙ i F i ṙ i F i ) = 0 The energy is conserved. MD samples correct microstates.
39 Definition of Temperature Instantaneous temperature T (t) = 2 1 3Nk B 2 i m i ṙ 2 i (t) Temperature T = T (t)
40 Equipartition Theorem The equipartition theorem is a general formula that relates the temperature of a system with its average energies. For a system described by the Hamiltonian (energy function) H(r, ṙ) r i H r j = p i H p j = δ ij k B T For an ideal gas H = H kin = 3 2 k BT
41 Canonical Ensemble System variables: N, V, T Closed system exchanging energy with a heat bath at temperature T. How to simulate on a computer? We can directly manipulate velocities (temperature). No extended heat bath needed. Molecular dynamics with periodic boundary conditions, but manipulate ṙ i during the integration.
42 Control of velocities/temperature Differential control: the temperature is fixed to the predescribed value and no fluctuations occur. Proportional control: velocities are corrected at each integration step through a coupling constant towards the prescribed temperature. The coupling constant determines the strength of the fluctuations around T. Integral control: the Hamiltonian is extended and variables are introduced wich reflect the effect of an external system which fix the temperature. Time evolution is derived from the extended Hamiltonian. Stochastic control: position and velocities are propagated according to modified equations of motion with friction and stochastic forces. The final equations give the correct mean value of the temperature.
43 NVE H(Γ) = Γ = {r N, p N } N i=1 1 2m i p i + U(r N ) microstate All accessible microstates are equally probable const. W NVE (Γ) = δ(h(γ) E) statistical weight Distribution of microstates Ω NVE = Γ δ(h(γ) E) 1 N! 1 h 3N dr N p N δ(h(γ) E) 5 / 1
44 Canonical Ensemble (NVT) Microstates follow the Boltzmann s distribution [ W NVT (Γ) = exp H(rN, p N ] ) k B T Canonical partition function: describes the statistical properties Q NVT = Γ Probability distribution [ exp H(Γ) ] 1 k B T N!h 3N P(Γ) = W NVT(Γ) Q NVT dr N dp N [ exp H(rN, p N ] ) k B T and the expectation value of the total energy is (β 1/k B T ) { E = ln Q NVT = 1 [ H(Γ) exp H(Γ) ] } β Q NVT k B T Γ
45 Velocity Rescaling At t and temperature T (t) all velocities are multiplied by λ T = 1 2 N i=1 To control the temperature 2 m i (λv i ) Nk B 2 N i=1 λ = T 0 /T (t) 2 m i vi 2 3 Nk B = (λ 2 1)T (t) Easy to implement Sampling does not correspond to any ensemble Not recommended for production MD Good for equilibration Not time-reversible, not deterministic 24 / 1
46 Perform one (or more) MD steps Compute the instantanous kinetic energy Rescale the velocities by a factor 25 / 1
47 Jumps by Rescaling 20 15!=1 K 10 5 K t 20!= t 26 / 1
48 Berendsen Thermostat Weak coupling to a heat bath: temperature correction dt (t) dt = 1 τ (T bath T (t)) τ coupling parameter: larger τ weaker coupling ) Exponential decay of the system towards the desired temperature T = δt τ (T bath T (t)) τ, thermostat is inactive too small τ, unrealistically low fluctuations τ = δt, standard rescaling Same problems as velocity rescaling T Time 27 / 1
49 Like a friction Continuous formulation of the velocity rescaling (1 ζδt) = ṗ i = f i ζp i 1 δt τ p 2 i i m i Nk B T 1 Friction positive when K > Nk B T /2, cooling. 1/ !=1 K t 20 15!=10 K t
50 Langevin Thermostat Random forces due to impacts exerted by molecules in the environment: irregular motion and friction. At each step all v i are corrected by a random force and a friction term m r i = U r i mγṙ i + W i (t) Frictional force and random motion are related, i.e. fluctuation dissipation theorem is satisfied W i (t), W j (t ) = δ ij δ(t t )6mΓk B T Gaussian distribution with infinitely short correlation time (great number of collisions) Transition probability P t = i ( p i P + U P + m i D 2 P m i r i r i p i p 2 + Γ p ) ip i p i
51 Stochastic Behavior The canonical distribution is then recovered as ( ) lim P = C exp 1 t k B T { mp 2 i 2 + U(rN )} i and friction and diffusion are related D = Γ m k BT Local thermostat: noise term depends on the particle Samples from canonical ensemble and ergodic. Allows larger time steps compared to non-stochastic thermostats. Momentum transfer is destroyed, i.e., no diffusion coefficients. 34 / 1
52 Nose Thermostat Heat bath integral part of the system: s, ṡ, and Q (coupling) L = i ms 2 ṙ 2 i 2 U(r N ) Qṡ2 gk B T ln s The conjugate momenta in the extended system are p i L ṙ i = m i s 2 ṙ i p s L ṡ = Qṡ where the auxiliary variable is a time-scaling parameter dt = sdt p i = p/s real system momentum Microcanonical ensemble in the extended system: for g = 3N + 1 canonical ensemble for the real system 38 / 1
53 Nose Hoover Using only differentiation in real time t : 1 s H Nose = i d dt ṙ i = p i m i ṗ i = U ξp i r i ξ = sp s /Q ṡ s ξ = 1 Q = d ln s dt ( i = ξ p 2 i m i g β ) p i + U(r N ) + ξ2 Q 2m i 2 + g ln s β conserved 39 / 1
54 Friction Coefficient Define the relaxation time as ν T 3Nk B T /Q dynamics of friction coefficient by (feed-back mechanism) p 2 [ ] i ξ = νt 2 i m i 3Nk B T 1 T = νt 2 T 1 Too large Q: canonical distribution after long simulation Too small Q: high-frequency T oscillations
55 Ergodicity Problems Microcanonic Andersen Nose-Hoover v r r r NVE: closed loop, constant energy shell Andersen: points corresponding to a Gaussian velocity distribution Nose Hoover: band trajectory determined by the initial conditions Not sufficiently chaotic to sample all phase-space; trapped in a subspace 41 / 1
56 Nose Hoover Chain Ergodicity is improved by thermostatting the thermostat variable: M Nose Hoover thermostats ṙ i = p i m i ṗ i = U ξ 1 p i r i ξ 1 = 1 Q 1 ( i p 2 i g m i β ) ξ i ξ 2 ξ j = 1 Q j ( Qj 1 ξ 2 j 1 k B T ) ξ j ξ j+1 ξ M = 1 ( QM 1 ξm 1 2 k B T ) Q M H NHC = i p i 2m i + U(r N ) + j Q j ξ 2 j 2 M + gk B Ts 1 + k B Ts j j=2 42 / 1
57 Thermostats Comparison 600 velocity velocity rescaling position velocity Berendsen position velocity Nose-Hoover chain position 2 4 Time [ps] / 1
58 Thermostats Table tune cont. L/G correct ergodic cons. q. determ. Velocity rescaling G? X Andersen X L X X Berendsen X X G? X Nosé-Hoover X X L/G X X X Nosé-Hoover chains X X L/G X X X X Langevin X X L X X X Stochastic velocity rescaling X X L/G X X X 44 / 1
59 Other Ensembles N, V, E microcanonical inner energy U N, V, T canonical Helmholtz free energy F F = U TS N, P, H isobaric-isoenthalpic enthalpy H H = U + pv N, P, T isobaric-isothermal Gibbs free energy G G = U TS + pv = H TS
60 How to simulate an isobaric ensemble? We can manipulate the volume of the system. This emulates an external piston used to keep the pressure constant. Instantaneous pressure p(t) fluctuates, forces on the container Averaging p(t) gives the observable internal pressure.
61 Pressure: Virial Theorem From equipartition theorem pi 2 H kin = = 1 U(r) r i 2m i 2 r i i The total potential U(r) is coming from in internal potential V int (r) and and a wall potential W (r). W (r) r i = p df r = p dv (divr) = 3pV r i i i The virial theorem is then 3pV = 2 H kin i V (r) r i r i
62 Isobaric-Isoenthalpic Ensemble The system exchange work with an external surface: Ω is a variable Cubic box Ω h = 0 Ω Ω 1 3 Scaled coordinates r i = Ω 1/3 s i Lagrangian in terms of scaled coordinates External applied pressure P L = 1 m i Ω 2 3 ṡ2 2 i U(s N, Ω) i P Ω work on the system to balance P: the internal pressure fluctuates 46 / 1
63 Barostat Correct the internal pressure by scaling the voluem, i.e. the inter-particle distances The system is coupled to a barostat Isobaric-Isoenthalpic ensemble NPH P(r N, p N, Ω) δ(c K(p N ) U(r N, Ω) ΠΩ) where Π is the internal pressure Isobaric-Isothermal ensemble NPT ( P(r N, p N, Ω) exp K(pN ) + U(r N ), Ω) + ΠΩ k B T 47 / 1
64 Berendsen Barostat Deviation from desired pressure dp(t) dt = P 0 P(t) τ P τ P relaxation time Rescale volume by η the volume, η 1/3 coordinates γ is the compressibility. η(t) = 1 δt τ P γ(p 0 P(t)) Asotropic or anisotropic, τ P strength of coupling, not correct ensemble distribution
65 Andersen Barostat Associate a kinetic energy to the variable volume L = 1 m i Ω 2 3 ṡ2 2 i W Ω 2 U(s N, Ω) PΩ i The coupling mimics the action of a piston of mass W (strength of coupling) Equations of motion ( W Ω = 1 Ω 2 3 3Ω s i = 1 U 2 m i Ω 2 3 s i 3 m i ṡ 2 i i i Ω Ω s i ) U s i P = P int P s i Feedback loop to adjust the simulation volume: internal pressure fluctuates around the external pressure 49 / 1
66 Correct Distribution Conserved Hamiltonian of the 6N + 2 dimensional system H = 1 2 i π 2 i m i Ω Π2 2W + U(sN, Ω) + PΩ low W will result in rapid box size oscillations large W slow adjustmen of the volume Partition function (N, P, H) = dω dπ dπ N ) s N δ (H + Π2 2W H trajectories consistent with the NPH ensemble The combination with one of the constant temperature methods allows the simulation in the NPT ensemble. 50 / 1
67 Experiment with barostat and thermostat Cons. [kj/mol] V [A 3 ] P [katm] Berendsen Thermostat and Barostat τt =0.1 τp =1, 2, 4, time (ps) Cons. [kj/mol] V [A 3 ] P [katm] Nose--Hoover Thermostat and Barostat τt =0.1 τp =1, 2, 4, time (ps) 51 / 1
68 Multiple Time Step Method (r-respa) M.E. Tuckerman et al. JCP, (1992) Equation of motion ṙ = p m ṗ = F fast + F slow il = Liouville operator ( ) p m r + F fast + F slow p p + il ref + il slow Integrator x(δt) = exp(ilδt)x(0) = exp(il slow δt/2) [exp(il ref δt)] n exp(il slow δt/2) + O(δt 3 )
69 MTS in CP2K Two FORCE_EVAL sections defining methods Example DFT with hybrid functionals F ref = F GGA F slow = F hybrid
70 PLUMED: MD Driver and Enhanced Sampling Definition of many Collective Variables (CV) Analysis features for MD and MetaDynamics Bias methods for enhanced sampling MD Logarithmic Mean Force Dynamics Method Experiment Directed Simulation Methods
71 i-pi: External MD Driver Kapil et al., Comp. Phys. Comm. 236, (2018) Driver software communicating with force engines (e.g. CP2K) through socket interfaces
72 i-pi: Features MD and Path-Integral MD (PIMD) for NVE, NVT, NPT Ring Polymer Contraction (RPC) MD and Centroid MD Generalized Langevin Equation (GLE) thermostats PI+GLE and PIGLET Ring Polymer Instantons Thermodynamic Integration, Geomerty optimization, Saddle point search Harmonic vibrations Multiple Time Step algorithms Metadynamics (interface to PLUMED) Replica Exchage MD Second Generation CP like integration
73 i-pi: Multiple Time Steps
74 Summary: BOMD in CP2K Born Oppenheimer MD with ASPC(3) is default in CP2K SCF convergence criteria depends on system is a reasonable starting guess Best used together with OT and the FULL_SINGLE_INVERSE preconditioner, for large system an iterative update of the preconditioner should be used PRECOND_SOLVER INVERSE_UPDATE Langevin dynamics is an option needs some special care, analysis of forces
MD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
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 informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationGear methods I + 1/18
Gear methods I + 1/18 Predictor-corrector type: knowledge of history is used to predict an approximate solution, which is made more accurate in the following step we do not want (otherwise good) methods
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 informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More informationAb initio Molecular Dynamics Born Oppenheimer and beyond
Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator
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 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 informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationAb Initio Molecular Dynamic
Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master / RFCT Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD)
More informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationAb Initio Molecular Dynamic
Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD) 1 - Which
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
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 informationClassical Molecular Dynamics
Classical Molecular Dynamics Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Motivation n Types of
More informationLangevin Dynamics in Constant Pressure Extended Systems
Langevin Dynamics in Constant Pressure Extended Systems D. Quigley and M.I.J. Probert CMMP 2004 1 Talk Outline Molecular dynamics and ensembles. Existing methods for sampling at NPT. Langevin dynamics
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 informationChapter 4 Molecular Dynamics and Other Dynamics
Chapter 4 Molecular Dynamics and Other Dynamics Molecular dynamics is a method in which the motion of each individual atom or molecule is computed according to Newton s second law. It usually involves
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationNanoscale simulation lectures Statistical Mechanics
Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationCHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued
CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
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 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 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 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 informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
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 informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
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 informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationConstant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems.
Constant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems. David Quigley A thesis submitted for the degree of Doctor of Philosophy University
More informationFree 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 informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationUNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES
UNIVERSITY OF OSLO FCULTY OF MTHEMTICS ND NTURL SCIENCES Exam in: FYS430, Statistical Mechanics Day of exam: Jun.6. 203 Problem :. The relative fluctuations in an extensive quantity, like the energy, depends
More informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationComputational Physics
Computational Physics Molecular Dynamics Simulations E. Carlon, M. Laleman and S. Nomidis Academic year 015/016 Contents 1 Introduction 3 Integration schemes 4.1 On the symplectic nature of the Velocity
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationEnhanced sampling via molecular dynamics II: Unbiased approaches
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
More informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
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 informationAb initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory Ondrej Marsalek 1 1, a) and Thomas E. Markland Department of Chemistry,
More informationModeling of Nanostructures and Materials Jacek A. Majewski Summer Semester 2013 Lecture 6 Lecture Modeling of Nanostructures Molecular Dynamics
Chair of Condensed Matter Physics nstitute of Theoretical Physics Faculty of Physics, Universityof Warsaw Summer Semester 013 Lecture Modeling of Nanostructures and Materials Jacek A. Majewski E-mail:
More informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More informationA Brief Introduction to Ab Initio Molecular Dynamics
A Brief Introduction to Ab Initio Molecular Dynamics Matt Probert Condensed Matter Dynamics Group, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Overview of Talk In this talk
More informationAnalysis of MD Results Using Statistical Mechanics Methods. Molecular Modeling
Analysis of MD Results Using Statistical Mechanics Methods Ioan Kosztin eckman Institute University of Illinois at Urbana-Champaign Molecular Modeling. Model building. Molecular Dynamics Simulation 3.
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationTHERMOSTATS AND THERMAL TRANSPORT IN SOLIDS
Hands-On Tutorial Workshop, July 19 th 2011 THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS Christian Carbogno FRITZ-HABER-INSTITUT DER MAX-PLANCK-GESELLSCHAFT, BERLIN - GERMANY Hands-On Tutorial Workshop,
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationMolecular Dynamics at Constant Pressure: Allowing the System to Control Volume Fluctuations via a Shell Particle
Entropy 2013, 15, 3941-3969; doi:10.3390/e15093941 Review OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume
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 informationCitation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n.
University of Groningen Stochastic concepts in molecular simulation Hess, Berk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check
More informationDeclaration. I, Derrick Beckedahl, declare that
Abstract In this dissertation we reformulate the configurational temperature Nosè-Hoover thermostat proposed by Braga and Travis (C. Braga and K.P. Travis, 005), using the antisymmetric matrix formalism
More informationarxiv: v1 [cond-mat.stat-mech] 28 Mar 2008
Canonical sampling through velocity-rescaling Giovanni Bussi, Davide Donadio, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials
More informationMolecular dynamics simulations
Chapter 8 Molecular dynamics simulations Most physical systems are collections of interacting objects. For example, a drop of water contains more than 10 21 water molecules, and a galaxy is a collection
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationHeat Transport in Glass-Forming Liquids
Heat Transport in Glass-Forming Liquids by VINAY VAIBHAV The Institute of Mathematical Sciences CIT Campus, Taramani, Chennai 600 113, India. A thesis submitted in partial fulfillment of requirements for
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
More informationFrom quantum to classical statistical mechanics. Polyatomic ideal gas.
From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
More informationSub -T g Relaxation in Thin Glass
Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition
More informationIV. Classical Molecular Dynamics
IV. Classical Molecular Dynamics Basic Assumptions: 1. Born-Oppenheimer Approximation 2. Classical mechanical nuclear motion Unavoidable Additional Approximations: 1. Approximate potential energy surface
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
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 informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationA (short) practical introduction to kinetic theory and thermodynamic properties of gases through molecular dynamics
A (short) practical introduction to kinetic theory and thermodynamic properties of gases through molecular dynamics Miguel A. Caro mcaroba@gmail.com March 28, 2018 Contents 1 Preface 3 2 Review of thermodynamics
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationNext generation extended Lagrangian first principles molecular dynamics
Next generation extended Lagrangian first principles molecular dynamics Anders M. N. Niklasson Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (Dated: June 1, 017) LA-UR-17-3007
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 informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationStructure and Dynamics of Hydrogen-Bonded Systems October 2009
2054-7 Structure and Dynamics of Hydrogen-Bonded Systems 26-27 October 2009 Aqueous Water and Its Liquid-Vapor Surface from First-Principles by a Car-Parrinello-like Approach Thomas D. KUEHNE Harvard University
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 informationEfficient use of CP2K
Efficient use of CP2K Fawzi Mohamed 1 1 Department of Chemisty Humboldt University of Berlin June 10, 2009 cp2k Fortran 95 opensource (GPL) general framework for different methods DFT no official release
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
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 informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez Part I: ( background ) This Lecture Post-Modern Classical Mechanics Hamiltonian Mechanics Numerical Considerations Thermodynamics,
More informationSupporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationCalculating Vibrational Spectra from Molecular Dynamics
Calculating Vibrational Spectra from Molecular Dynamics A Simulating a Trajectory with Wannier Centers To calculate IR spectra from Molecular Dynamics, it is necessary to have dipole information for the
More informationQuantum (Path Integral) Molecular Dynamics
Quantum (Path Integral) Molecular Dynamics Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Background
More informationStatistical Mechanics. Atomistic view of Materials
Statistical Mechanics Atomistic view of Materials What is statistical mechanics? Microscopic (atoms, electrons, etc.) Statistical mechanics Macroscopic (Thermodynamics) Sample with constrains Fixed thermodynamics
More informationCP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
More informationThe micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction
ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,
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 informationStatistical thermodynamics (mechanics)
Statistical thermodynamics mechanics) 1/15 Macroskopic quantities are a consequence of averaged behavior of many particles [tchem/simplyn.sh] 2/15 Pressure of ideal gas from kinetic theory I Molecule =
More informationChapter 5 - Systems under pressure 62
Chapter 5 - Systems under pressure 62 CHAPTER 5 - SYSTEMS UNDER PRESSURE 5.1 Ideal gas law The quantitative study of gases goes back more than three centuries. In 1662, Robert Boyle showed that at a fixed
More information