Monte Carlo strategies for first-principles simulations of elemental systems
|
|
- Stella Griffith
- 6 years ago
- Views:
Transcription
1 Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 1 / 28
2 Introduction Atomistic simulations Require a potential, which gives the energy of a system as a function of atomic coordinates. Classical (empirical) V (r N ) = i V 1 (r i ) + i>j V 2 ( rij ) + i>j>k V 3 ( rij, r jk, r ik ) +... Quantum-mechanical (semi-empirical or first-principles ) Wavefunction-based methods H ψ = E ψ Density Functional Theory (DFT) E = F[ρ] only real choice for condensed-phase systems. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 2 / 28
3 Introduction First-principles simulations of phase behavior First-principles energies are desirable when chemical interactions are significant, rather than physical ones: high-pressure melting of Si 1, Fe 2, Cu 3, MgO 4, Mo 5, C 6 others... High-pressure solid-solid transitions (minerology) Phosphorous LLE 7 Metallic hydrogen by QMC 8 Water (by Monte Carlo 9 and CPMD 10 ) 1 Sugino and Car, PRL 74 (1995) Alfè et al., Nature 401 (1999) Vočadlo et al., JCP 120 (2004) Alfè, PRL 94 (2005) Belonoshko et al., PRL 92 (2004) Wang et al., PRL 95 (2005) Ghiringhelli and Meijer, JCP 122 (2005) Delaney et al., PRL 97 (2006) McGrath et al., Comp. Phys. Comm 169 (2005) 289; Kuo et al. JCPB 108 (2004) 12990; McGrath et al., ChemPhysChem 6 (2005) Many. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 3 / 28
4 Introduction Monte Carlo basics Starting at some initial state i... 1 Generate a new trial state j. For instance: Displace one atom Insert one atom Change one atom to a different element Many others! 2 Evaluate E = E j E i 3 Accept or reject the new state according to π j /π i = e β E, etc. 4 Repeat 1-3 many times. For few-body empirical potentials, evaluating E is an O(N) operation. For many-body potentials, its O(N 2 ) or greater. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 4 / 28
5 Introduction Molecular dynamics vs Monte Carlo Most phase coexistence calculations with first-principles potentials are done using molecular dynamics. Direct simulation of a two-phase system Thermodynamic integration and reference systems Most phase coexistence calculations with empirical potentials are done using Monte Carlo. Open system simulations - direct access to µ Multiple-cell (Gibbs Ensemble) simulations Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 5 / 28
6 Introduction Molecular dynamics vs Monte Carlo Most phase coexistence calculations with first-principles potentials are done using molecular dynamics. Direct simulation of a two-phase system Thermodynamic integration and reference systems Most phase coexistence calculations with empirical potentials are done using Monte Carlo. Open system simulations - direct access to µ Multiple-cell (Gibbs Ensemble) simulations In molecular dynamics, all the atoms are moved at once. In Monte Carlo, atoms are (usually) moved a few at a time. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 5 / 28
7 Introduction Molecular dynamics vs Monte Carlo Most phase coexistence calculations with first-principles potentials are done using molecular dynamics. Direct simulation of a two-phase system Thermodynamic integration and reference systems Most phase coexistence calculations with empirical potentials are done using Monte Carlo. Open system simulations - direct access to µ Multiple-cell (Gibbs Ensemble) simulations In molecular dynamics, all the atoms are moved at once. In Monte Carlo, atoms are (usually) moved a few at a time. For many-body potentials, this is extremely expensive. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 5 / 28
8 Introduction Benchmark DFT calculations Scale of the problem Determining liquid-state properties and/or phase coexistence requires: 1 N particles (atoms or molecules) (as few as 30 atoms in some cases.) 2 Simulating for a sufficient duration: Molecular dynamics: 10 ps up to many ns. ( steps) Monte Carlo: O(10 4 ) trial moves per particle. These are serial processes: because equilibration may take significant time, many very short trajectories are not equivalent to one long one. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 6 / 28
9 Introduction Benchmark DFT calculations Cubic scaling with system size Most of our calculations are performed using NWChem, a well-parallelized and widely-used code from PNL 11. Benchmark calculations: startup SCF cycles, revpbe XC functional with TM pseudopotentials, Γ-point only. Calculations run on crystalline Aluminum with random vacancies. 11 Valiev et al., Comp. Phys. Comm. 181 (2010) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 7 / 28
10 Introduction Benchmark DFT calculations Parallel efficiency on Ranger Good scaling can be obtained up to 512 processors, but only for sufficiently large problems. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 8 / 28
11 Introduction Benchmark DFT calculations Time to self-consistent solution The time required to converge the DFT energy to a specified tolerance can vary significant from one configuration to the next. This seriously impacts efficiency of any method where multiple energy evaluations are performed in parallel. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 9 / 28
12 Introduction Benchmark DFT calculations Time to self-consistent solution (2) If the number of atoms fluctuates (as in an open ensemble) the variation is larger: Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 10 / 28
13 Introduction Benchmark DFT calculations Benchmark summary Computational cost scales as N 3. Good parallel scaling obtained, but only for relatively large N/P The total wall-clock time per calculation will not be reduced much below O(1 min). Time to SCF solution varies substantially even between configurations with similar density and N. Multiple simultaneous evaluations present a load-balancing problem. (Not shown) Computational cost can vary dramatically between different XC functionals. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 11 / 28
14 Presampling Code and early work Pre-sampling with an approximate potential 12 Use an approximate potential, E a, to generate large moves with a high probability of acceptance by the correct potential E c. One cycle: 1 For I = 1, M 1 Generate displacement (or particle exchange, etc...) 2 Accept or reject based on e β Ea (for canonical ensemble) 2 Accept or reject the final state after M moves on: e β( Ec Ea) ; generally, α ij = min ( π c j πi a πi c πj a If M N = N particles, then this requires only one E c evaluation per cycle, instead of N., 1 ) 12 Akhmatskaya et al., CPL 267 (1997) 105; Liu and Chen J. Am. Stat. Assoc., 93 (1998) 1032; Iftimie et al., JCP 113 (2000) 4852; Liu, Monte Carlo Strategies for Scientific Computing (2001) Springer; Hetényi et al., JCP 117 (2002) 8203; Gelb, JCP 118 (2003) 7747; Zwin and Luo, JCP 123 (2005) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 12 / 28
15 Presampling Code and early work Code design and implementation Capabilities: 1 NVT, µvt, NPT, NhT, NVE, NPH, GEMC, binary-mixture GEMC, NPH-GEMC, and three-cell NVE-GEMC simulations 2 Parallel tempering 3 Use of external electronic structure codes unmodified 4 Approximate-potential optimization Architecture choices: Written in Python and pypar Extensive use of Numpy, Scipy, Multiprocessing modules Expensive routines in C/C++ (through Weave) Built-in architecture aware processor scheduler Interface with external codes through standard I/O and text files Current code: NWChem and CPMD integrated as drivers, and all basic functionality met in lines. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 13 / 28
16 Presampling Code and early work Some early work: Lithium NPT Consistently high densities Large XC functional dependence Some cutoff dependence Little system-size dependence Poor parallel tempering efficiency Gelb and Carnahan, CPL 417 (2006) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 14 / 28
17 Presampling Code and early work Some early work: Lithium GEMC Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 15 / 28
18 Presampling Efficiency Pre-sampling efficiency analysis For cycles of M moves: Cost of N moves with the correct potential: Nw c Cost of N moves using presampling: (N/M)w c + Nw a, Speedup is then: S = 1 w a /w c + 1/M M However, this doesn t address efficiency! If E c E a grows linearly with M, S = 1 w a /w c + 1/M e σm Acceptance probability drops exponentially with presampling chain length. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 16 / 28
19 Presampling Efficiency Proxy system: embedded-atom potentials Generic form (Daw and Baskes, 1983): V = G i (ρ i ) + ( ) V 2 rij i i>j ρ i = j ρ a j (r ij) G i is the embedding energy for placing atom i into the electron density provided by all the other atoms. ρ a j (r ij) is the spherically-averaged electron density from atom j at position r i. Sutton-Chen potential (Sutton and Chen, 1990): V /ɛ = i j>i ( a r ij ) n c ρi i ρ i = j i ( a ) m r ij Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 17 / 28
20 Presampling Efficiency Pre-sampling performance with cycle length Test system: qsc simulations of liquid copper, with different potential truncations ( correct is nm.) Speedups as predicted in this instance. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 18 / 28
21 Presampling Three-stage presampling Three-stage sampling Introduce a third, middle -level energy evaluation: E a. Modified algorithm: 1 Every cycle of M moves, accept/reject the final state by: E 1 = E a E a 2 Every M cycles of M moves, accept/reject the final state by: E 2 = E c E a Application: E a : empirical potential E a : first-principles potential, low-quality E c : first-principles potential, high-quality Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 19 / 28
22 Presampling Three-stage presampling Low-cutoff DFT as a presampler Test: 50 random (but nonoverlapping) configurations of Al with N = 30. Densities between 0.01 Å 3 and 0.05 Å 3. Take correct potential as E = 20.0 H cutoff. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 20 / 28
23 Presampling Three-stage presampling Low-cutoff DFT as a presampler Compare differences between consecutive configurations. Absolute performance is less critical for E c E a. E c will be much smaller in practice (Configurations will be strongly correlated.) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 21 / 28
24 Presampling Three-stage presampling Three-stage pre-sampling scaling analysis In one overall cycle of M M moves, do: 1 E c evaluation M E a evaluations M M E a evaluations. S = If the cost w a is again negligible, we get [ 1 M M + 1 w a + M ] 1 w a M w c M w c [ ] 1 S = M w a /w c + 1/M [ ] wc M w a Overall sampling quality improved by a factor of w c /w a, ignoring acceptance rate issues. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 22 / 28
25 Presampling Three-stage presampling Example application Liquid Li at 1500 K, 1 bar; LDA/TM theory. Total simulation time for each run: 140 hours. Acceptance frequency at top stage with M = 10: 95%. The three-stage sampling is twice as efficient! Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 23 / 28
26 Coarse-grained parallelism Tree sampling Generate X new configurations in parallel, and try them one at a time; if one is accepted, discard the rest. The probability of accepting a move is Inefficiency (states never considered): I = X j=1 P tree = 1 (1 P acc) X X(X j)p acc(1 P acc) j 1 Constant P acc Constant X P acc X P tree I/X P acc X P tree I/X For 4 work, get 3 the throughput (moves accepted per wallclock time.) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 24 / 28
27 Coarse-grained parallelism Tree sampling performance This works as advertised, provided that: We use a three-stage scheme, where the middle-level potential is DFT run for a fixed number of SCF iterations, rather than to a specified convergence. Test case: 32 Al atoms at 3000 K and 1 kbar. E a : quantum-corrected Sutton-Chen potential E a : revpbe, 10.0 H, max 50 SCF cycles. E c : revpbe, 10.0 H, converged M % accepted (P acc ) 5 63% 10 57% 20 40% 40 22% Also, > 95% acceptance on E a /E c (M = 5). counts run time (seconds) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 25 / 28
28 Coarse-grained parallelism Tree sampling performance (2) There are now many factors influencing simulation efficiency: M, low-level cycle-length Dynamically adjustable for target P acc. M, mid-level cycle-length Number of SCF iterations at middle level; determines P acc. X, number of tree branches Trade-off with more processors per DFT evaluation. 2nd-generation trees? Target higher P acc, because need to accept either y 1 or y 2 for good efficiency! For P acc = 0.5, accept one move 75% of the time, and two 56% of the time. y 2 y 21 x y 1 y 11 y 12 y 22 Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 26 / 28
29 Conclusions Summary and prognosis Quality of the empirical approximation potential is critical. On-the-fly refinement of the approximate potential? Mid-level low-quality first-principles potentials are a time-saver, but do not substantially improve sampling. (Not shown) OBMC sampling on the whole system is inefficient at least for small numbers of replicas. Tree sampling on the whole system can be very effective. Requires fixed-iteration-count middle-level potential. Multi-generation trees possible, but not yet tested. Alleviates poor scaling of DFT codes for large P. Support from NSF (# ) Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 27 / 28
30 Coarse-grained parallelism Exploiting parallelism: whole-system OBMC (Also called Multiple-Try Monte Carlo (Liu 2001)) Generate X new configurations {y i } in parallel, each by an approximate chain of M moves, and choose one according to P(k) = π c (y k )/ j πc (y j ). Then generate X 1 backwards chains {x j } to calculate the final acceptance probability: W (y k ) W (x 0 ) = i πc (y i ) j πc (x j ) πa (x 0 ) π a (y k ) Sequence requires wall-clock time for two consecutive first-principles calculations. Performance: so far, not very good. For X 4, no significant improvement over regular sampling. Issue: since forward move often uphill in energy, the backwards chains tend to be very uphill. Lev Gelb (UT Dallas) Monte Carlo Strategies XSEDE12 28 / 28
Accelerated 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 informationComplementary approaches to high T- high p crystal structure stability and melting!
Complementary approaches to high T- high p crystal structure stability and melting! Dario ALFÈ Department of Earth Sciences & Department of Physics and Astronomy, Thomas Young Centre@UCL & London Centre
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 informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
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 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 informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
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 informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
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 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 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 informationMD 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 informationAll-Electron Path Integral Monte Carlo (PIMC) Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
All-Electron Path Integral Monte Carlo (PIMC) Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas Kevin Driver and Burkhard Militzer Department of Earth and Planetary Science University
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 informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationSystematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids
Systematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids Cameron F. Abrams Department of Chemical and Biological Engineering Drexel University Philadelphia, PA USA 9 June
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 informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationATOMISTIC MODELING OF DIFFUSION IN ALUMINUM
ATOMISTIC MODELING OF DIFFUSION IN ALUMINUM S. GRABOWSKI, K. KADAU and P. ENTEL Theoretische Physik, Gerhard-Mercator-Universität Duisburg, 47048 Duisburg, Germany (Received...) Abstract We present molecular-dynamics
More informationStructure-Property Relationships of Porous Materials for Carbon Dioxide Separation and Capture
Supporting Information Structure-Property Relationships of Porous Materials for Carbon Dioxide Separation and Capture Christopher E. Wilmer, 1 Omar K. Farha, 2 Youn-Sang Bae, 3,a Joseph T. Hupp, 2 and
More informationLecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications
Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
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 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 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 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 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 Simulation Background
Molecular Simulation Background Why Simulation? 1. Predicting properties of (new) materials 2. Understanding phenomena on a molecular scale 3. Simulating known phenomena? Example: computing the melting
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationOPENATOM for GW calculations
OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW
More informationReferences in the Supporting Information:
Identification of the Selective Sites for Electrochemical Reduction of CO to C2+ Products on Copper Nanoparticles by Combining Reactive Force Fields, Density Functional Theory, and Machine Learning Supporting
More informationQuantum computing with superconducting qubits Towards useful applications
Quantum computing with superconducting qubits Towards useful applications Stefan Filipp IBM Research Zurich Switzerland Forum Teratec 2018 June 20, 2018 Palaiseau, France Why Quantum Computing? Why now?
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationab initio Electronic Structure Calculations
ab initio Electronic Structure Calculations New scalability frontiers using the BG/L Supercomputer C. Bekas, A. Curioni and W. Andreoni IBM, Zurich Research Laboratory Rueschlikon 8803, Switzerland ab
More informationModule 16. Diffusion in solids II. Lecture 16. Diffusion in solids II
Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,
More informationQMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
More informationOn the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study
On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study Wynand Dednam and André E. Botha Department of Physics, University of South Africa, P.O.
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More informationSupporting Information for. Dynamics Study"
Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationCCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation
CCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation Biswajit Santra 1, Angelos Michaelides 1,2, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG, Berlin,
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 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 informationMOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM
MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM MOLECULAR DYNAMICS SIMULATION STUDY OF SOLID-LIQUID INTERFACE PROPERTIES OF HCP MAGNESIUM By YUNFEI BAI, B. Eng.
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 informationBottom-up modelling of charge transport in organic semiconductors
Bottom-up modelling of charge transport in organic semiconductors David L. Cheung Department of Chemistry & Centre for Scientific Computing University of Warwick LCOPV 2010 Boulder 7th-10th August 2010
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
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 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 informationMolecular Dynamics of Covalent Crystals
Molecular Dynamics of Covalent Crystals J. Hahn and H.-R. Trebin Institut für Theoretische und Angewandte Physik, Universität Stuttgart, D-70550 Stuttgart, Germany Abstract. A molecular mechanics-like
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
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 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 informationCHAPTER 16 A MACROSCOPIC DESCRIPTION OF MATTER
CHAPTER 16 A MACROSCOPIC DESCRIPTION OF MATTER This brief chapter provides an introduction to thermodynamics. The goal is to use phenomenological descriptions of the microscopic details of matter in order
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationAb initio molecular dynamics simulation on temperature-dependent properties of Al Si liquid alloy
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 16 (4) 57 514 PII: S953-8984(4)7691-8 Ab initio molecular dynamics simulation on temperature-dependent properties
More informationAb initio statistical mechanics
Ab initio statistical mechanics Luca M. Ghiringhelli Fritz-Haber-Institut der MPG, Berlin Hands-on Workshop Density-Functional Theory and Beyond: Accuracy, Efficiency and Reproducibility in Computational
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
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 informationThree examples of a Practical Exact Markov Chain Sampling
Three examples of a Practical Exact Markov Chain Sampling Zdravko Botev November 2007 Abstract We present three examples of exact sampling from complex multidimensional densities using Markov Chain theory
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 informationPath Integral Monte Carlo Simulations on the Blue Waters System. Burkhard Militzer University of California, Berkeley
Path Integral Monte Carlo Simulations on the Blue Waters System Burkhard Militzer University of California, Berkeley http://militzer.berkeley.edu Outline 1. Path integral Monte Carlo simulation method
More informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
More informationUniform-acceptance force-biased Monte Carlo: A cheap way to boost MD
Dresden Talk, March 2012 Uniform-acceptance force-biased Monte Carlo: A cheap way to boost MD Barend Thijsse Department of Materials Science and Engineering Delft University of Technology, The Netherlands
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
More informationBasic introduction of NWChem software
Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationQuantum Monte Carlo methods
Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground
More informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
More informationDepartment of Chemical Engineering University of California, Santa Barbara Spring Exercise 3. Due: Thursday, 5/3/12
Department of Chemical Engineering ChE 210D University of California, Santa Barbara Spring 2012 Exercise 3 Due: Thursday, 5/3/12 Objective: To learn how to write & compile Fortran libraries for Python,
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 informationMicroporous Carbon adsorbents with high CO 2 capacities for industrial applications
Microporous Carbon adsorbents with high CO 2 capacities for industrial applications Santiago Builes, a,b Thomas Roussel,* b Camelia Matei Ghimbeu, c Julien Parmentier, c Roger Gadiou, c Cathie Vix-Guterl
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationLarge-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
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 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 informationHigh Temperature High Pressure Properties of Silica From Quantum Monte Carlo
High Temperature High Pressure Properties of Silica From Quantum Monte Carlo K.P. Driver, R.E. Cohen, Z. Wu, B. Militzer, P. Lopez Rios, M. Towler, R. Needs, and J.W. Wilkins Funding: NSF, DOE; Computation:
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
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 informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More informationWinmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17
Winmostar tutorial LAMMPS Melting point V7.025 X-Ability Co,. Ltd. question@winmostar.com Contents Configure I. Build solid phase II. Equilibration of solid phase III. Equilibration of liquid phase IV.
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationMonte Carlo Study of the Crystalline and Amorphous NaK Alloy
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 08C (207) 25 22 International Conference on Computational Science, ICCS 207, 2-4 June 207, Zurich, Switzerland Monte Carlo
More information1. Thermodynamics 1.1. A macroscopic view of matter
1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.
More informationPredictive Computing for Solids and Liquids
Predictive Computing for Solids and Liquids So Hirata Department of Chemistry May 214 Blue Waters Symposium 1 Schrödinger equation for a water molecule 1-particle, 3-dimensional partial differential equation
More informationRenormalization of Tensor- Network States Tao Xiang
Renormalization of Tensor- Network States Tao Xiang Institute of Physics/Institute of Theoretical Physics Chinese Academy of Sciences txiang@iphy.ac.cn Physical Background: characteristic energy scales
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 informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
More informationEquilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach
THE JOURNAL OF CHEMICAL PHYSICS 123, 124912 2005 Equilibrium sampling of self-associating polymer solutions: A parallel selective tempering approach Chakravarthy Ayyagari, a Dmitry Bedrov, and Grant D.
More information