Multiscale Simulations
|
|
- Priscilla Hudson
- 5 years ago
- Views:
Transcription
1 Multiscale Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia nstitute of Technology Atlanta, GA 30332, U.S.A.
2 DFT-MD coupling Topics First-Principles Molecular Dynamics DFT-KMC coupling on-the-fly KMC simulation MD-FEM coupling Quasicontinuum method Coarse-Grained Molecular Dynamics
3 Modeling & Simulation at Multiple Scales s ms μs ns pico-s femto-s Time Scales Finite Element Analysis Dislocation Dynamics Kinetic Monte Carlo Molecular Dynamics / Force Field Tight Binding Density Functional Theory Self-Consistent Field (Hartree-Fock) Quantum Monte Carlo Length Scales nm μm mm m Various methods used to simulate at different length and time scales
4 Zoo of Multiscale Simulation Methods First-principles MD (quantumatomistic coupling) Ehrenfest MD Born-Oppenheimer MD Car-Parrinello MD on-the-fly KMC (DFT-KMC coupling) QM/MM coupling Mathematical Homogenization Heterogeneous Multiscale Method Multiscale FEM quasi-continuum coarse-grained molecular dynamics variational multiscale method concurrent coupling coupled atomistic/discretedislocation adaptive multiscale modeling, bridging scale method bridging domain method DD/FEM coupling TB/MD/FEM coupling <new species born each year>
5 First-Principles Molecular Dynamics The major idea is to replace the predefined potentials in classical molecular dynamics (MD) by first-principles electronic structure calculation on-the-fly (i.e. keep electronic variables as active degrees of freedom in MD). The Algorithm: 1. solve the electronic structure problem for a set of ionic coordinates 2. evaluate forces 3. move atoms 4. repeat
6 Classical MD Based on Newtonian dynamics where approx M R () t = F () t = V { R ()} t e is a few-body-additive-interaction approximation of the true potential energy surface. The electrons follow adiabatically the classical nuclear motion and can be integrated out so that the nuclei evolve on a single global Born-Oppenheimer potential energy surface. A priori construction of the global potential energy surface suffers from the dimensionality bottleneck. Multiscale Systems Engineering Research Group ( ) N approx e J = 1 < J ( ) () 1 ( ) ( 2) { R ( )} = R ( ) + ( R ( ), R ( )) V t V t V t t ( R (), R (), R ()) + + ( 3) V t t t J K < J< K
7 Time-Dependent Schrödinger Equation where 2 H = 2M i Φ r R t = Φ r R t i i t 2 ({ },{ }; ) H ({ },{ }; ) ZZe e J i i 2me i< j r r < J R R J, i r R i j i = ({ },{ }) + V r R i n e i 2M i 2m e 2 2 = + H ({ r},{ R }) e i 2M Total wavefunction can be decomposed as t ( ) ( ) ( ) i Φ { r},{ R }; t Ψ { r}; t χ { R }; t exp dt' E ( t') i i e t0 with a phase factor E * * = d r d RΨ χ H Ψ χ e e Multiscale Systems Engineering Research Group Ze
8 Ehrenfest Molecular Dynamics Simultaneously solve E M R () t = V { ()} t e ( R ) = Ψ H Ψ e * = drψ H Ψ e i Ψ= t H e Ψ 2 2 = Ψ+ V ({ r},{ R }) Ψ i n e i i 2me
9 Born-Oppenheimer Molecular Dynamics n ground state BOMD, the time-independent electronic structure problem is solved self-consistently from each time for a given configuration of nuclei M R () t = min Ψ H Ψ Ψ 0 e 0 H Ψ = E Ψ e Electrons are explicitly set to be fully relaxed for a given configuration of nuclei, in contrast to Ehrenfest MD where electron relaxation is implicit by solving the timedependent Schrödinger equation.
10 BOMD with HF-SCF Assuming single Slater determinant Ψ = det{ ψ 0 j } M R () t = min Ψ H Ψ { ψ } 0 e 0 j Constrained (orthonormal orbitals) minimization min Ψ H Ψ st.. ψ ψ = δ { ψ } 0 e 0 i j ij j Define Lagrangian L = Ψ H Ψ + ( ) 0 e 0 Λ ψ ψ δ ij, ij i j ij The necessary condition of optimality L ψ j = 0 leads to Hartree-Fock equations HF H ψ = Λ ψ e j ij j i Then the new equations of motion are M R () t = min Ψ H Ψ { ψ } 0 0 j e H ψ = Λ ψ HF e j i ij j Multiscale Systems Engineering Research Group
11 BOMD with DFT H Based on the Hellmann-Feynman theorem, MD force is He F = Ψ H Ψ = Ψ H Ψ = Ψ Ψ 0 e 0 0 e R R Recall ( ) 2 e J { },{ } = r R + + e i i i 2 m e i< j r r < J R R J, i r R i j i The force is computed by DFT as ZZe He F = Ψ Ψ = 0 0 R ZZe Ze = + Ψ Ψ Ze 2 2 J R (,..., ) (,..., ) d... d 0 1 M 0 1 M 1 M J r r r r r r R R i r R J i 2 2 ZZe Ze J ρ() d J r r R R r R J = +
12 BOMD with DFT 1. Fix positions of nuclei {R 1,,R N }, solve DFT equations self-consistently; Ze ρ( ') ' r dr + + Vxc () r ψi () r = εψ i i () r 2 r R r r' ρ() r = f ψ () r 2. Find electrostatic force on each atom; ZZe Ze F r dr 2 2 J = () + ρ J R R r R J i 3. Perform a time step and find new positions of nuclei; i i 2 4. Repeat; M R () t = F
13 Drawbacks of BOMD The need to fully relax electronic subsystem while moving the atoms makes it computationally expensive. Full self-consistency at each MD step may not be necessary, especially when system is far from its equilibrium, since one simply needs a rough idea of the force field for a given atomic configuration
14 Car-Parrinello Molecular Dynamics R. Car and M. Parrinello (1985) Unified approach for molecular dynamics and density-functional theory. Phys. Rev. Lett. 55: Publication and citation analysis: : number of publications which appeared up to the year n that contain the keyword ab initio molecular dynamics (or synonyma first-principles MD, Car- Parrinello simulations etc.) in title, abstract or keyword list. : number of publications which appeared up to the year n that cite the 1985 paper by Car and Parrinello
15 Car-Parrinello Lagrangian L Car and Parrinello (1985) postulated the Lagrangian 1 1 = M R + μ ψ ψ Ψ H Ψ + Λ ( ψ ψ δ ) i 0 0 i ij 2 2 CP i i e ij i i ij Kinetic Energy Potential Energy Constraints where μ i s are fictitious masses for the dynamics of orbitals ψ i s.
16 Car-Parrinello Equations of Motion From Euler-Lagrange differential equations in classical mechanics (to ensure δ L ( R, R, ψ, ψ ; tdt ) = 0 ) CP i i d L L CP CP = dt R R d L L CP CP = dt ψ ψ Car-Parrinello equations of motion is derived as i M ( t) R = Ψ H Ψ + constraints({ ψ },{ }) R R i { R } 0 e 0 i δ δ μψ () t = Ψ H Ψ + constraints({ },{ }) δψ δψ { ψ R } i i 0 e 0 i i i
17 CPMD with DFT Car-Parrinello equations of motion M R E [{ ψ },{ R }] DFT i = R 2 Ze = ρ(,) r t dr r R ( t) δe [{ ψ },{ R }] DFT i μψ (,) r t = + ε ψ (,) r t i i ij j δψ (,) r t i j = H ψ (,) r t + ε ψ (,) r t DFT i ij j j
18 Verlet algorithm in CPMD Orbital Dynamics First, a verlet step ignoring orthogonality constraint Δt H ψ DFT i ψ ( t +Δ t) = 2ψ ( t) ψ ( t Δ t) + i i i μ Then, restore orthogonality ψ ( t +Δ t) = ψ ( t +Δ t) + ξ ψ ( t) i i i 2 ( ) ( ) Computationally dynamics is applied to c i s in the reciprocal space with Kohn-Sham orbitals 1 ψ (, rk) = c ( Gk, )exp i( + ) i i G k r Ω G i t
19 CPMD Code The CPMD code is a planewave implementation of DFT for first-principles molecular dynamics First version by Jürg Hutter at BM Zurich Research Lab with dozens of other contributors The code is copyrighted jointly by BM Corp and by Max Planck nstitute, Stuttgart t is distributed free of charge to non-profit organizations (
20 CPMD capabilities Wavefunction optimization: direct minimization and diagonalization Geometry optimization: local optimization and simulated annealing Molecular dynamics: NVE, NVT, NPT Path integral MD Response functions Excited states Time-dependent DFT (excitations, MD in excited states) Coarse-grained non-markovian metadynamics Wannier, EPR, Vibrational analysis QM/MM See on-line manual at:
21 On-the-fly KMC Searching saddle points on the potential energy surface (PES) on-the-fly while performing KMC simulation Search Saddle Points on PES Find Activation Energy Calculate Rate Constants (TST/hTST) Simulate Phase Transition by KMC
22 On-the-fly KMC 1. Start from a minimum configuration; 2. Randomly generate a set of configurations around the minimum and search the saddle points by the Dimer method; 3. Locate the saddle points connect to the minimum; 4. nsert the new events and propensities in the event table in KMC and simulate one step; 5. Repeat;
23 Diffusion of adatom on Al(100) Multiscale Systems Engineering Research Group
24 ΔR V V 1 V 0 2 F 1 F F R 2 Dimer energy: E = V1+ V2 Curvature along Dimer: ˆN C Dimer Method [Henkelman & Jónsson 1999] = E (2 V ) ( ΔR) Estimate: F1 = V1 F2 = V2 E ΔR V ˆ 0 = + ( F1 F2) in 2 4 F R = ( F 1 + F 2) 2 2. Dimer rotates to find the lowest curvature mode of PES, i.e. minimize Dimer energy 3. Translate Dimer towards uphill according to 4. Repeat = 1 2 = F F F F ( F F )/2 F F if C > 0 = FR 2F if C < 0
25 Quasicontinuum (QC) Method Based on the full atomistic model, use mesh to reduce the degrees of freedom representative atoms (repatoms)
26 QC Method Displacement of atom i: Displacement of N atoms: {u 1,u 2,,,u N } Empirically the total energy is the sum of site energy of each atom the Stillinger-Weber (SW) type site energy for atom i is with two-body potential and three-body potential where
27 QC Method The total potential energy of the system (atoms + external loads) is where f i u i is the potential energy of the applied load f i on atom i The goal of the static QC method is to find the atomic displacements that minimize the total potential such that the number of degrees of freedom is substantially reduced from 3N; the computation of the total potential is accurately approximated without the need to explicitly compute the site energy of all the atoms; the critical regions can evolve with the deformation by addition/removal of repatoms.
28 QC Method Reduce DoF Any atom not chosen as a repatom is constrained to move according to the interpolated displacements This first approximation of the QC then, it to replace the energy E tot by E tot,h : with continuum displacement field where S α is the interpolation function with local support
29 QC Method Local Energy The Cauchy-Born rule assumes that a uniform deformation gradient at the macro-scale can be mapped directly to the same uniform deformation on the micro-scale. Thus, every atom in a region subject to a uniform deformation gradient will be energetically equivalent. The energy within an element in crystals can be estimated by computing the energy of one atom in the deformed state and multiplying by the number of atoms in the element. Energy density for each element is where Ω 0 is the unit cell volume and E 0 is the energy of the unit cell when its lattice vectors are distorted according to F The total energy of an element is Multiscale Systems Engineering Research Group
30 QC Method Nonlocal Energy Local energy does not approximate well where deformation of crystal is non-uniform (e.g. surfaces and interfaces) and shorter than the cutoff radius of inter-atomic potential. Energy-based formulation: Nonlocal energy is weighted sum of those of repatoms as where Force-based formulation: The force on repatom β is determined by its neighborhood C β (α,β, for repatom) Atomic-level force in neighborhood c of repatom α weights of the atomic forces Multiscale Systems Engineering Research Group
31 QC Method detailed issues Local-nonlocal energy coupling where N loc +N nl =N rep Local/nonlocal criterion: whether a repatom should be local or nonlocal? whether there is significant variantion of the deformation gradient Effects of local-nonlocal interface: Polycrystals Elastic/plastic deformation decomposition
32 QC Applications nanoindentation (Smith et al. 2001) Silicone Phase transformation and dislocation nucleation observed Different phases observed
33 QC Applications dislocation behavior (Rodney & Phillips 1999) Dislocation junction under shear stress
34 Coarse Grained (CG) Molecular Dynamics displacement of mesh node j (does not have to coincident with an atom) u = j μ fj μ u μ, is a weighted average of displacements of atoms μ s. j Displacement field u( x) = N ( x) u j j j μ Coarse grained energy is the average of the canonical ensemble of the atomistic Hamiltonian on the constrained phase space {(x,p)} ( ) u, MD E u, u = H = d x d p H e Δ/ Z where β=1/(kt), Z d d e βh MD = x p Δ is partition function, and μ μ Δ= δ ( u f ) ( f / m ) j j u δ u μ μ jμ j p μ μ jμ μ enforces constraints. The atomistic Hamiltonian MD k k MD u μ μ MD k k βh p p H = + E + D T μ μ coh 1 T u u μ μ μν ν μ 2m μ μν, 2 Cohesive energy of atoms
35 CGMD Partition function Z = Z kin Z pot Potential part of partition function is CG potential energy then is w/ stiffness Computationally, the full CG energy is where T ( N N ) u K 2 2 Z ( u, β) = CC β e β pot k atom node j jk k 3 1 T E ( u ) = log Z = ( N N ) kt + u K u β 2 2 u pot k pot atom node j jk k ( ) int Thermal 1 E u, u = U + ( M u u + u K u ) T T k k jk j k j jk k 2 j, k coh atom atom node U = N E + 3( N N ) kt int K = ( f D f ) 1 1 jk μν, jμ μν kν Harmonic is internal energy and x x K f m f m N N 1 1 = ( ) = ( ) ( ) jk μ jμ μ kμ μ μ j μ k μ is CG mass matrix
36 CGMD Applications MEMS/NEMS NEMS silicon microresonator The coarse grained (CG) region comprises most of the volume The molecular dynamics (MD) region contains most of the simulated degrees of freedom the CG mesh is refined to the atomic scale where it joins with the MD lattice.
37 Summary First-Principles MD On-the-fly KMC Quasicontinuum method Coarse-Grain Molecular Dynamics
38 Further Readings First-Principles Molecular Dynamics Marx D. and Hutter J. (2000) Ab nitio Moleular Dynamics: Theory and mplementation. n: J. Grotendorst (Ed.) Modern Methods and Algorithms of Quantum Chemistry (Jülich: John von Neumann nstitute for Computing, SBN ), Vol.3, pp On-the-fly KMC Mei, D., Ge, Q., Neurock, M., Kieken, L., and Lerou, J. (2004) First-principles-based kinetic Monte Carlo simulation of nitric oxide decomposition over Pt and Rh surfaces under lean-burn conditions. Molecular Physics, 102(4): Kratzer P. and Scheffler M. (2002) Reaction-limited island nucleation in molecular beam epitaxy of compound semiconductors. Physical Review Letters, 88(3): (1-4) Battaile C.C., Srolovitz D.J., Oleinik.., Pettifor D.G., Sutton A.P., Harris S.J., and Butler J.E. (1999) Etching effects during the chemical vapor deposition of (100) diamond. Journal of Chemical Physics, 111(9): Henkelman, G. and Jónsson, H. (2001) Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table. Journal of Chemical Physics, 115(21): Trushin, O., Karim, A., Kara, A. and Rahman, T.S. (2005) Self-learning kinetic Monte Carlo method: Application to Cu(111). Physical Review B, 72: (1-9). Xu, L. and Henkelman, G. (2008) Adaptive kinetic Monte Carlo for first-principles accelerated dynamics. Journal Chemical Physics, 129: (1-9). Mei D., Xu L., and Henkelman G. (2009) Potential energy surface of Methanol decomposition on Cu(110). Journal of Physical Chemistry C, 113:
39 Further Readings Quasicontinuum Tadmor E.B., Ortiz M., and Phillips R. (1996) Quansicontinuum analysis of defects in solids. Philosophical Magazine A, 73(6): Shenoy V.B., Miller R., Tadmor E.B., Rodney D., Phillips R., and Ortiz M. (1998) An adaptive methodology for atomic scale mechanics: the quasicontinuum method. J. Mech. Phys. Sol., 47: Knap J. and Ortiz M. (2001) Ana analysis of the quasicontinuum method. J. Mech. Phys. Sol., 49: Smith G.S., Tadmor E.B., Bernstein N. and Kaxiras E. (2001) Multiscale simulations of Silicon Nanoindentation. Acta Mater., 49: Rodney D. and Phillips R. (1999) Structure and strength of dislocation juncitons: An atomic level analysis. Physical Review Letters, 82(8): Coarse-grained molecular dynamics Rudd R.E. and Broughton J.Q. (1998) Coarse-grained molecular dynamics and the atomic limit of finite elements. Physical Review B, 58(10):R5893-R5896 Rudd R.E. and Broughton J.Q. (2005) Coarse-grained molecular dynamics: Nonlinear finite elements and finite temperature. Physical Review B, 72(14): Variational multiscale method Hughes T.J.R., Feijóo G.R., Mazzei L., and Quincy J.-B. (1998) The variational multiscale method a paradigm for computational mechanics. Computational Methods in Applied Mechanics & Engineering, 166(1-2):3-24
40 Further Readings Concurrent coupling Broughton J.Q., Abraham F.F., Bernstein N., Kaxiras E. (1999) Concurrent coupling of length scales: methodology and application. Physical Review B, 60(4): Bridging domain method Xiao S.P. and Belytschko T. (2004) A bridging domain method for coupling continua with molecular dynamics. Computational Methods in Applied Mechanics & Engineering, 193(17-20): Bridging scale method Liu W.K., Park H.S., Qian D., Karpov E.G., Kadowaki H., Wagner G.J. (2006) Bridging scale methods for nanomechanics and materials. Computational Methods in Applied Mechanics & Engineering, 195(13-16): Wagner G.J., Liu W.K. (2003) Coupling of atomistic and continuum simulations using a bridging scale decomposition. Journal of Computational Physics, 190(1): Park H.S., Karpov E.G., Liu W.K., Klein P.A. (2005) The bridging scale for twodimensional atomistic/continuum coupling. Philosophical Magazine, 85(1):
Quantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
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 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 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 informationA Quasicontinuum for Complex Crystals
A Quasicontinuum for Complex Crystals Ellad B. Tadmor Aerospace Engineering and Mechanics University of Minnesota Collaborators: Previous Work: Current Work: U. V. Waghmare, G. S. Smith, N. Bernstein,
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 informationCar-Parrinello Molecular Dynamics
Car-Parrinello Molecular Dynamics Eric J. Bylaska HPCC Group Moral: A man dreams of a miracle and wakes up with loaves of bread Erich Maria Remarque Molecular Dynamics Loop (1) Compute Forces on atoms,
More informationError Control for Molecular Statics Problems
International Journal for Multiscale Computational Engineering, 4(5&6)647 662(26) Error Control for Molecular Statics Problems Serge Prudhomme, Paul T. Bauman & J. Tinsley Oden Institute for Computational
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 28 Jun 2001
arxiv:cond-mat/665v [cond-mat.mtrl-sci] 28 Jun 2 Matching Conditions in -Continuum Modeling of Materials Weinan E and Zhongyi Huang 2 Department of Mathematics and PACM, Princeton University and School
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 informationError Control for Molecular Statics Problems
Error Control for Molecular Statics Problems Serge Prudhomme 1, Paul T. Bauman 2, and J. Tinsley Oden 3 Institute for Computational Engineering and Sciences The University of Texas at Austin Austin, Texas
More informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
More informationElectronic-structure calculations at macroscopic scales
Electronic-structure calculations at macroscopic scales M. Ortiz California Institute of Technology In collaboration with: K. Bhattacharya, V. Gavini (Caltech), J. Knap (LLNL) BAMC, Bristol, March, 2007
More informationAn Atomistic-based Cohesive Zone Model for Quasi-continua
An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationAb Initio Molecular Dynamics: Theory and Implementation
John von Neumann nstitute for Computing Ab nitio Molecular Dynamics: Theory and mplementation Dominik Marx and Jürg Hutter published in Modern Methods and Algorithms of Quantum Chemistry, Proceedings,
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 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 informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
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 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 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 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 informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
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 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 informationA Coupled Discrete/Continuous Method for Computing Lattices. Application to a Masonry-Like Structure
Author manuscript, published in "International Journal of Solids and Structures 48, 21 (2011) 3091-3098" DOI : 10.1016/j.ijsolstr.2011.07.002 A Coupled Discrete/Continuous Method for Computing Lattices.
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 informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationFirst-Principles Quantum Simulations of Nanoindentation
Chapter 10 First-Principles Quantum Simulations of Nanoindentation Qing Peng Additional information is available at the end of the chapter http://dx.doi.org/10.5772/48190 1. Introduction Modern nanotechnology
More informationConical Intersections. Spiridoula Matsika
Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The
More informationAtomistic aspects: the role of DFT calculations
Atomistic aspects: the role of DFT calculations Peter Kratzer Fritz-Haber-Institut der MPG, Berlin, Germany How to obtain information about the energetics at the atomic level From statistical physics to
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationMultiscale Modeling with Extended Bridging Domain Method
Multiscale Modeling with Extended Bridging Domain Method Mei Xu, Robert Gracie, and Ted Belytschko Northwestern University 2145 Sherman Ave, Evanston, IL, 60208 1 Introduction This Chapter describes the
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationMolecular Mechanics: The Ab Initio Foundation
Molecular Mechanics: The Ab Initio Foundation Ju Li GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA 2 Outline Why are electrons quantum? Born-Oppenheimer
More informationto appear in International Journal for Multiscale Computational Engineering
to appear in International Journal for Multiscale Computational Engineering GOAL-ORIENTED ATOMISTIC-CONTINUUM ADAPTIVITY FOR THE QUASICONTINUUM APPROXIMATION MARCEL ARNDT AND MITCHELL LUSKIN Abstract.
More informationSet the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
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 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 informationSummary of the new Modelling Vocabulary
Summary of the new Modelling Vocabulary These two pages attempts to summarise in a concise manner the Modelling Vocabulary. What are Models? What are Simulations? Materials Models consist of Physics or
More informationNuggets on coarse- graining and mul,scale computa,onal schemes Maria Fyta Ins,tut für Computerphysik, Universität Stu<gart Stu<gart, Germany
Nuggets on coarse- graining and mul,scale computa,onal schemes Maria Fyta Ins,tut für Computerphysik, Universität Stu
More informationHigh Performance Multiscale Simulation for Crack Propagation
1 High Performance Multiscale Simulation for Crack Propagation Guillaume Anciaux, Olivier Coulaud and Jean Roman ScAlApplix Project HPSEC 2006-18th August 2 Outline 1. Introduction Motivations State of
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 informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Molecular Mechanics by Ju Li. Given August 9, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use
More informationThe electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More informationKinetic Monte Carlo: from transition probabilities to transition rates
Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot
More informationMOLECULAR MODELLING OF STRESSES AND DEFORMATIONS IN NANOSTRUCTURED MATERIALS
Int. J. Appl. Math. Comput. Sci., 2004, Vol. 14, No. 4, 541 548 MOLECULAR MODELLING OF STRESSES AND DEFORMATIONS IN NANOSTRUCTURED MATERIALS GWIDON SZEFER Institute of Structural Mechanics, AGH University
More informationA Multilattice Quasicontinuum for Phase Transforming Materials: Cascading Cauchy Born Kinematics
Journal of Computer-Aided Materials Design manuscript No. (will be inserted by the editor) A Multilattice Quasicontinuum for Phase Transforming Materials: Cascading Cauchy Born Kinematics M. Dobson R.
More informationReactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)
Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li
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 informationSimulation Methods II
Simulation Methods II Maria Fyta Institute for Computational Physics Universität Stuttgart Summer Term 2018 SM II - contents First principles methods Hartree-Fock and beyond Density-funtional-theory Ab
More informationFirst-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design
First-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design Peter Kratzer 5/2/2018 Peter Kratzer Abeokuta School 5/2/2018 1 / 34 Outline
More informationKinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys
Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys, Scott Dunham Department of Electrical Engineering Multiscale Modeling Hierarchy Configuration energies and transition
More informationQUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology http://www.tut.fi/semiphys CONTENTS MOTIVATION PATH INTEGRAL APPROACH
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
More informationAb initio phonon calculations in mixed systems
Ab initio phonon calculations in mixed systems Andrei Postnikov apostnik@uos.de Outline: Experiment vs. ab initio theory Ways of theory: linear response and frozen phonon approaches Applications: Be x
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
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 informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationIndentation of Silicon: Phase Transition? Indentation of Silicon: Phase Transition?
Indentation of Silicon: Phase Transition? Kallman et al., Phys. Rev. B. 47, 7705 (1993). Smith et al., Acta. Mater. 49, 4089 (2001). MD: NVT 350,000 atoms Cut away of crystalline Si indented with a tetrahedral
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
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 informationQuantum Molecular Dynamics Basics
Quantum Molecular Dynamics Basics Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More information5 Multiscale Modeling and Simulation Methods
5 Multiscale Modeling and Simulation Methods This chapter is dedicated to a discussion of multiscale modeling and simulation methods. These approaches are aimed to provide a seamless bridge between atomistic
More informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationAb initio molecular dynamics : BO
School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal Why ab initio MD Free of parametrization and only fundamental constants required. Bond
More informationMulti-paradigm modeling of fracture of a silicon single crystal under mode II shear loading
Journal of Algorithms & Computational Technology Vol. 2 No. 2 203 Multi-paradigm modeling of fracture of a silicon single crystal under mode II shear loading Markus J. Buehler 1, *, Alan Cohen 2, and Dipanjan
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 informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationTemperature-related Cauchy Born rule for multiscale modeling of crystalline solids
Computational Materials Science 37 (2006) 374 379 www.elsevier.com/locate/commatsci Temperature-related Cauchy Born rule for multiscale modeling of crystalline solids Shaoping Xiao a, *, Weixuan Yang b
More informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
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 informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationAdaptive modelling in atomistic-to-continuum multiscale methods
Journal of the Serbian Society for Computational Mechanics / Vol. 6 / No. 1, 2012 / pp. 169-198 (UDC: 530.145.6) Adaptive modelling in atomistic-to-continuum multiscale methods E. Marenic 1,2, J. Soric
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 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 informationMultiscale modeling of nano/micro systems by a multiscale continuum field theory
DOI 10.1007/s00466-010-0538-5 ORIGINAL PAPER Multiscale modeling of nano/micro systems by a multiscale continuum field theory Xiaowei Zeng Xianqiao Wang James D. Lee Yajie Lei Received: 8 March 2010 /
More informationModel operator approach to calculations of the Lamb shifts in relativistic many-electron atoms
Model operator approach to calculations of the Lamb shifts in relativistic many-electron atoms Vladimir Shabaev a in collaboration with Ilya Tupitsyn a and Vladimir Yerokhin b a St. Petersburg State University
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationDepartment of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an
The statistical characteristics of static friction J. Wang, G. F. Wang*, and W. K. Yuan Department of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an 710049, China * E-mail: wanggf@mail.xjtu.edu.cn
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
More informationJournal of the Mechanics and Physics of Solids
Journal of the Mechanics and Physics of Solids 59 (2011) 775 786 Contents lists available at ScienceDirect Journal of the Mechanics and Physics of Solids journal homepage: www.elsevier.com/locate/jmps?
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationPure and zero-point vibrational corrections to molecular properties
Pure and zero-point vibrational corrections to molecular properties Kenneth Ruud UiT The Arctic University of Norway I UNIVERSITETET TROMSØ June 30 2015 Outline Why consider vibrational effects? General
More informationCondensed matter physics FKA091
Condensed matter physics FKA091 Ermin Malic Department of Physics Chalmers University of Technology Henrik Johannesson Department of Physics University of Gothenburg Teaching assistants: Roland Jago &
More informationModeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques
Loughborough University Institutional Repository Modeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques This item was submitted to Loughborough University's
More informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationCoarse-graining molecular dynamics in crystalline solids
in crystalline solids Xiantao Li Department of Mathematics Penn State University. July 19, 2009. General framework Outline 1. Motivation and Problem setup 2. General methodology 3. The derivation of the
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationModeling and Simulation of Wave Propagation Based on Atomistic Field Theory
Xianqiao Wang James D. Lee Department of Mechanical and Aerospace Engineering, George Washington University, Washington, DC 20052 Qian Deng Department of Mechanical and Aerospace Engineering, University
More informationSTOCHASTIC SIMULATIONS WITH APPLICATIONS IN MATERIAL SCIENCES
STOCHASTIC SIMULATIONS WITH APPLICATIONS IN MATERIAL SCIENCES Deák Robert Ph.D. thesis In joint Ph.D. supervision between: Department of Physics Babeş Bolyai University, Cluj-Napoca, Romania Scientific
More informationPrinciples of Molecular Spectroscopy
Principles of Molecular Spectroscopy What variables do we need to characterize a molecule? Nuclear and electronic configurations: What is the structure of the molecule? What are the bond lengths? How strong
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 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 information