First-principles Molecular Dynamics Simulations
|
|
- Anis Collins
- 6 years ago
- Views:
Transcription
1 First-principles Molecular Dynamics Simulations François Gygi University of California, Davis MICCoM Computational School, Jul 17-19,
2 Outline Molecular dynamics simulations Electronic structure calculations First-Principles Molecular Dynamics (FPMD) Applications 2
3 Molecular Dynamics An atomic-scale simulation method Compute the trajectories of all atoms extract statistical information from the trajectories Atoms move according to Newton s law: m i!! Ri = F i
4 Molecular dynamics: general principles Integrate Newton s equations of motion for N atoms m!! i Ri (t) = F i (R 1,, R N ) i =1,, N F i (R 1,, R N ) = i E(R 1,, R N ) Compute statistical averages from time averages (ergodicity hypothesis) A = Ω dr 3N dp 3N A(r, p)e βh (r,p) 1 T A(t)dt Examples of A(t): potential energy, pressure, T 0
5 Simple energy model Model of the hydrogen molecule (H 2 ): harmonic oscillator E(R 1, R 2 ) = E( R 1 R 2 ) = α( R 1 R 2 d 0 ) 2 This model does not describe intermolecular interactions
6 Simple energy model Model of the hydrogen molecule including both intra- and intermolecular interactions: E(R 1,, R N ) = E intra ( R i R j ) + E inter ( R i R j ) {i, j} M i M j M ' This model does not describe adequately changes in chemical bonding
7 Simple energy model Description of the reaction H 2 +H H + H 2 The model fails!
8 What is a good energy model?
9 Atomistic simulation of complex structures Complex structures Nanoparticles Assemblies of nanoparticles Embedded nanoparticles Liquid/solid interfaces 9
10 A difficult case: Structural phase transitions in CO 2 Molecular phases polymeric phase
11 The energy is determined by quantum mechanical properties First-Principles Molecular Dynamics: Derive interatomic forces from quantum mechanics Ni-tris(2-aminoethylamine)
12 First-Principles Molecular Dynamics Monte Carlo Molecular Dynamics Statistical Mechanics FPMD Electronic Structure Theory Quantum Chemistry Density Functional Theory R. Car and M. Parrinello (1985) 12
13 Electronic structure calculations Problem: determine the electronic properties of an assembly of atoms using the laws of quantum mechanics. Solution: solve the Schroedinger equation!
14 The Schroedinger equation for N electrons A partial differential equation for the wave function ψ: r i R 3, ψ L 2 (R 3N ) i! t ψ(r 1,, r N,t) = H(r 1,, r N,t) ψ(r 1,, r N,t) H is the Hamiltonian operator: H 2! 2 r1,, rn, t) = i + V ( r,, rn, t) 2m ( 1 i
15 The time-independent Schroedinger equation If the Hamiltonian is time-independent, we have ψ(r 1,, r N,t) = ψ(r 1,, r N ) e iet/! where ψ(r) is the solution of the timeindependent Schroedinger equation: H( r,, rn ) ψ ( r1,, rn ) = Eψ ( r1,, r 1 N ) energy
16 Solving the Schroedinger equation The time-independent Schroedinger equation can have many solutions: H(r 1,, r N )ψ n (r 1,, r N ) = E n ψ n (r 1,, r N ) n = 0,1, 2 The ground state wave function ψ 0 describes the state of lowest energy Ε 0 Excited states are described by ψ 1, ψ 2,.. and have energies Ε 1, Ε 2,.. > Ε 0
17 Hamiltonian operator for N electrons and M nuclei Approximation: treat nuclei as classical particles Nuclei are located at positions R i, electrons at r i H(r 1,, r N, R 1,, R M ) =!2 2m e N 2 i i=1 N M i=1 j=1 Z j e 2 r i R j + M Z + i Z j e M M! 2 i R i R i R j 2 i< j i=1 N i< j e 2 r i r j
18 The adiabatic approximation The Hamiltonian describing an assembly of atoms is time-dependent because atoms move: ) ( )) ( ( ), ( ), ( 2 ), ( e-e ion 2 2 r r r r V t R r V t V t V m t H j j i i + = + =! time-dependence through ionic positions
19 The adiabatic approximation If ions move sufficiently slowly, we can assume that electrons remain in the electronic ground state at all times ψ ( r, t) = ψ ( r) H ( r,{ Ri ( t)}) ψ 0( r) = 0 E ψ ( r) 0 0 Ground state energy Ground state wave function
20 Mean-field approximation The problem of solving the N-electron Schroedinger equation is formidable (N! complexity). Wave functions must be antisymmetric (Pauli principle) Assuming that electrons are independent (i.e. feel the same potential) reduces this complexity dramatically. The potential is approximated by an average effective potential ),,,,,, ( ),,,,,, ( 1 1 N i j N j i r r r r r r r r ψ ψ = exchanged ),, ( ),, ( ),, ( N n n N n N E H r r r r r r ψ ψ =
21 Independent particles, solutions are Slater determinants A Slater determinant is a simple form of antisymmetric wave function : ψ ( r 1,, r N ) = det{ ϕ i ( r j )} The one-particle wave functions ϕ i satisfy the one-particle Schroedinger equation: h( r) ϕ ( r) h( r) i = ε ϕ ( r) 2! = 2m i 2 i + V eff ( r) Note: effective potential
22 Electron-electron interaction H(r 1,, r N, R 1,, R M ) =!2 2m e N 2 i i=1 N M i=1 j=1 Z j e 2 r i R j + M Z + i Z j e M M! 2 i R i R i R j 2 i< j i=1 N i< j e 2 r i r j
23 Density Functional Theory Introduced by Hohenberg & Kohn (1964) Chemistry Nobel prize to W.Kohn (1999) The electronic density is the fundamental quantity from which all electronic properties can be derived E = E ρ [ ] E [ ρ] = T [ ρ] + V(r)ρ(r)dr + E xc ρ [ ] Problem: the functional E[ρ] is unknown!
24 The Local Density Approximation Kohn & Sham (1965) E ρ [ ] = T ρ [ ] + V Approximations: (r)ρ(r)dr + E xc ρ [ ] The kinetic energy is that of a non-interacting electron gas of same density The exchange-correlation energy density depends locally on the electronic density E xc = E xc [ ρ(r) ] = ε xc (ρ(r))ρ(r)dr
25 The Local Density Approximation V e-e = ρ( r!) r r! d r! +V (ρ(r)) XC The mean-field approximation is sometimes not accurate, in particular for strongly correlated electrons excited state properties
26 The Kohn-Sham equations Coupled, non-linear, integro-differential equations: ( Δφ i +V(ρ, r)φ i = ε i φ i i =1 N el * * ρ( r #) V(ρ, r) = V ion (r)+ r r# d r # +V XC (ρ(r), ρ(r)) * ) N el * ρ(r) = φ i (r) 2 * i=1 * + * φ i (r)φ j (r)dr = δ ij
27 Numerical methods Basis sets: solutions are expanded on a basis of N orthogonal functions φ i (r) = N c ij ϕ j (r) j=1 ϕ j (r)ϕ k (r) = δ jk Ω R 3 Ω The solution of the Schroedinger equation reduces to a linear algebra problem
28 Numerical methods: choice of basis Gaussian basis (non-orthogonal) ϕ i (r) = e α i r R 2 Plane wave basis (orthogonal) ϕ q (r) = e iq R Other representations of solutions: values on a grid finite element basis
29 Numerical methods: choice of basis Hamiltonian matrix: H ij = ϕ i H ϕ j = ϕ i (r) Hϕ j (r) Ω d 3 r Schroedinger equation: an algebraic eigenvalue problem Hc n = ε n c n c n C N
30 Numerical methods: choice of basis Non-orthogonal basis sets lead to generalized eigenvalue problems S ij = φi φ j = Ω i φ ( r) φ ( r) j d 3 r δ ij Hc n = ε n Sc n c n C N
31 Solving large eigenvalue problems The size of the matrix H often exceeds Direct diagonalization methods cannot be used Iterative methods: Lanczos type methods subspace iteration methods Many algorithms focus on one (or a few) eigenpairs Electronic structure calculations involve many eigenpairs (~ # of electrons) robust methods are necessary
32 Solving the Kohn-Sham equations: fixed-point iterations The Hamiltonian depends on the electronic density ( Δφ i +V(ρ, r)φ i = ε i φ i i =1 N el * * ρ( r #) V(ρ, r) = V ion (r)+ r r# d r # +V XC (ρ(r), ρ(r)) * ) N el * ρ(r) = φ i (r) 2 * i=1 * + * φ i (r)φ j (r)dr = δ ij
33 Self-consistent iterations For k=1,2, Compute the density ρ k Solve the Kohn-Sham equations The iteration may converge to a fixed point
34 Simplifying the electron-ion interactions: Pseudopotentials The electron-ion interaction is singular V e-ion (r) = Ze2 r R Only valence electrons play an important role in chemical bonding Valence electrons core electrons
35 Simplifying the electron-ion interactions: Pseudopotentials The electron-ion potential can be replaced by a smooth function near the atomic core " $ V e-ion (r) = Ze2 # r R $ % f ( r R ) r R > r cut r R < r cut Core electrons are not included in the calculation (they are assumed to be "frozen")
36 Pseudopotentials: Silicon Solutions of the Schroedinger equation for Si including all electrons (core+valence): ψ 3s ψ 3p Potential = -Z/r Core Valence wavefunctions
37 Pseudopotentials: Silicon Solutions of the Schroedinger equation for Si including all electrons (zoom on core region): rψ 1s rψ 2s rψ 2p Potential = -Z/r Core region
38 Pseudopotentials: Silicon The electron-ion potential can be replaced by a smooth function near the atomic core ψ 3s ψ 3p pseudopotentials -Z/r Core Valence wavefunctions
39 Summary: First-principles electronic structure Time-independent Schroedinger equation Mean-field approximation Simplified electron-electron interaction: Density Functional Theory, Local Density Approximation Simplified electron-ion interaction: Pseudopotentials
40 Molecular dynamics: Computation of ionic forces Hamiltonian: H(λ) Hellman-Feynman theorem: if ψ 0 (λ) is the electronic ground state of H(λ) E λ λ0 = λ ψ 0 (λ) H(λ) ψ 0 (λ) = ψ 0 (λ 0 ) H(λ) λ λ0 ψ 0 (λ 0 ) For ionic forces: λ=r i (ionic positions) F i = E R i = ψ 0 H ψ 0 = ψ 0 V e-ion (r R j ) ψ 0 R i R i j
41 Integrating the equations of motion: the Verlet algorithm The equations of motion are coupled, second order ordinary differential equations Any ODE integration method can be used Time-reversible integrators are preferred The Verlet algorithm (or leapfrog method) is time-reversible x(t + Δt) = 2x(t) x(t + Δt) + Δt 2 m F(x(t))
42 Integrating the equations of motion: the Verlet algorithm Derivation of the Verlet algorithm: Taylor expansion of x(t) x(t + Δt) = x(t) + Δt dx dt + Δt 2 2 x(t Δt) = x(t) Δt dx dt + Δt 2 2 d 2 x dt + Δt d 2 x dt Δt d 3 x dt +O(Δt 4 ) 3 d 3 x dt +O(Δt 4 ) 3 Add the two Taylor expansions: x(t + Δt) + x(t Δt) = 2x(t) + Δt 2 d 2 x dt 2 +O(Δt 4 )
43 Integrating the equations of motion: the Verlet algorithm use Newton s law m d 2 x dt 2 = f (x(t)) x(t + Δt) + x(t Δt) = 2x(t) + Δt 2 d 2 x dt 2 +O(Δt 4 ) x(t + Δt) = 2x(t) x(t Δt) + Δt 2 m F(x(t)) +O(Δt 4 )
44 First-Principles Molecular Dynamics Molecular Dynamics Density Functional Theory m i d 2 dt 2 R i = F i FPMD ( Δ +V eff )ϕ i (x) = ε i ϕ i (x) ( ) = ϕ i (x) 2 ρ x n i=1 Newton equations Kohn-Sham equations R. Car and M. Parrinello (1985) 44
45 FPMD: the Recipe Choose a starting geometry: atomic positions Choose an exchange-correlation functional Choose appropriate pseudopotentials Run! Publish!! 45
46 FPMD: the Recipe Choose a starting geometry: atomic positions Choose an exchange-correlation functional Choose appropriate pseudopotentials Run! Publish!! Test! Test sensitivity to starting geometry, finite size effects Test sensitivity to duration of the simulation Test accuracy of the basis set Test choice of exchange-correlation functionals Test accuracy of pseudopotentials 46
47 First-Principles Molecular Dynamics applications Solid state physics Surface physics Nanotechnology High pressure physics Chemical Physics Biochemistry Mechanisms of drug action Solvation processes The absence of empirical parameters makes this approach widely applicable and predictive.
48 Nanoparticles Exploration of multiple locally stable structures Electronic properties at finite temperature Cd 34 Se 34 48
49 Embedded nanoparticles, Liquids and Liquid-Solid Interfaces assemblies of nanoparticles Annealing of structures at finite temperature Calculation of band gaps and band alignments Si/ZnS S. Wippermann, M. Vörös, A. Gali, F. Gygi, G. Zimanyi, G.Galli, Phys. Rev. Lett. 112, (2014).
50 Liquids and Liquid-Solid Interfaces Structure of water at the interface Electronic structure band alignment of bulk solid and liquid Spectroscopy IR and Raman spectra H 2 O/Si(100)H 50
51 Liquid-solid interfaces Water on oxide surfaces H 2 O/WO 3 Simulation of surface relaxation and dynamics structure of defects electronic structure spectroscopic signature of water at the interface H 2 O/WO 3 51
52 Electronic properties: Polarization The electronic polarization (per unit cell) of an infinite periodic system is ill-defined P = 1 % Ω e Z R ( ' l l + rρ(r) dr* & ) l P depends on the choice of origin The change in polarization caused by a small perturbation is well defined The electric current caused by a perturbation (e.g. a deformation) can be computed R. Resta, Rev. Mod. Phys. 66, 899 (1994). 52
53 Electronic properties: Polarization The electronic polarization (per unit cell) of an infinite periodic system is ill-defined P = 1 % Ω e Z R ( ' l l + rρ(r) dr* & ) l P depends on the choice of origin The change in polarization caused by a small perturbation is well defined The electric current caused by a perturbation (e.g. a deformation) can be computed R. Resta, Rev. Mod. Phys. 66, 899 (1994). 53
54 Wannier functions A set of localized orbitals that span the same subspace as the Kohn-Sham eigenvectors minimize the spread σ 2 = φ ( x φ x φ ) 2 φ Wannier centers: centers of charge of each Wannier function Polarization can be expressed in terms of the centers P = 1 % Ω e Z ( ' lr l + e r w n (r) dr* & ) l n N. Marzari, A. Mostofi, J. Yates, I. Souza and D. Vanderbilt, Rev. Mod. Phys. 84, 1419 (2012). 54
55 Time-dependent polarization of nanoparticles PBE DFT MD 300K dt=1.9 fs Cd 34 Se Debye 55
56 IR Spectroscopy IR spectra during MD simulations Autocorrelation function of P(t) ( ) = 2πω 2 β 3cVn ( ω) α ω e iωt P µ ( 0) P ν ( t) µν dt 56
57 Raman Spectroscopy Compute the polarizability at each MD step Use Density Functional Perturbation Theory (Baroni, Giannozzi, Testa, 1987) Use a finite-difference formula with P(t) and finite field 57
58 On-the-fly Computation of Raman spectra (D 2 O) 64 Position of O-D stretching band: PBE functional yields a red shifted peak, compared to expt. Low frequency bands: satisfactory agreement with expt. Peak Intensities in good agreement with expt. Q. Wan, L. Spanu, G. Galli, F. Gygi, JCTC 9, 4124 (2013)
59 Solving the Kohn-Sham equations in a finite electric field In finite systems: add a linear potential H KS = p2 2m +V(r) eex The spectrum is not bounded below (no "ground state") In periodic systems: define the electric enthalpy: F φ [ ] = E KS φ [ ] ΩP φ [ ] E I. Souza, J. Iniguez, D. Vanderbilt, Phys. Rev. Lett. 89, (2002). 59
60 Si(100):H-H2O interface DFT MD of the Si/H2O interface under finite field Si(100)-(3x3):H-H2O, canted dihydride surface termination, 116 water molecules Analysis of time-dependent polarization Comparison with IR spectra L. Yang, F. Niu, S. Tecklenburg, M. Pander, S. Nayak, A. Erbe, S. Wippermann, F. Gygi, G. Galli 60
61 Validation of DFT: PBE vs PBE0 vs Oxygen-oxygen pair correlation function in (H 2 O) 32 PBE (439±29K) PBE (367±25K) PBE0 (438±29K) PBE0 (374±27K) Exp (300K) C.Zhang, D.Donadio, F.Gygi, G.Galli, JCTC 7, 1443 (2011) 61
62 Is my simulation reproducible? D 2 O Power spectrum of ionic velocities (32 x 10 ps runs) intensity (arb. units) frequency(cm -1 ) 62
63 Validating/comparing levels of theory Need for (quantitative) statistical analysis compute confidence intervals An accurate determination of structural and electronic properties requires multiple uncorrelated simulations Autocorrelation times may vary for different quantities Example: the PBE400 dataset First-principles MD simulations of water 63
64 Summary Basic features of FPMD Approximations of electronic structure calculations Extensions: polarization, finite electric field Applications Next FPMD steps: Today 1:45 pm: Qbox tutorial Tomorrow 1:30 pm: Qbox hands-on exercises 64
Compression Algorithms for Electronic Structure Computations
Compression Algorithms for Electronic Structure Computations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org IPAM Workshop, Computation
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 informationHigh pressure core structures of Si nanoparticles for solar energy conversion
High pressure core structures of Si nanoparticles for solar energy conversion S. Wippermann, M. Vörös, D. Rocca, A. Gali, G. Zimanyi, G. Galli [Phys. Rev. Lett. 11, 4684 (213)] NSF/Solar DMR-135468 NISE-project
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 informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
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 informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
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 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 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 informationMultiple Exciton Generation in Si and Ge Nanoparticles with High Pressure Core Structures
Multiple Exciton Generation in Si and Ge Nanoparticles with High Pressure Core Structures S. Wippermann, M. Vörös, D. Rocca, A. Gali, G. Zimanyi, G. Galli NanoMatFutur DPG-214, 4/3/214 Multiple Exciton
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 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 informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More 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 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 informationSearch for materials to harvest light
Solar nanocomposites with complementary charge extraction pathways for electrons and holes: Si embedded in ZnS S. Wippermann, M. Vörös, F. Gygi, A. Gali, G. Zimanyi, G. Galli NanoMatFutur DPG-2014, 04/03/2014
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 informationWhy use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationDensity Functional Theory (DFT)
Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid
More informationElectrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation
Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationIntro to ab initio methods
Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio
More informationBefore we start: Important setup of your Computer
Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann
More informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
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 informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
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 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 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 informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
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 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 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 informationSpring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor
2145-25 Spring College on Computational Nanoscience 17-28 May 2010 Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor Stefano BARONI SISSA & CNR-IOM DEMOCRITOS Simulation Center
More informationTHE UNIVERSITY OF CHICAGO FIRST PRINCIPLES SIMULATIONS OF VIBRATIONAL SPECTRA OF AQUEOUS SYSTEMS A DISSERTATION SUBMITTED TO
THE UNIVERSITY OF CHICAGO FIRST PRINCIPLES SIMULATIONS OF VIBRATIONAL SPECTRA OF AQUEOUS SYSTEMS A DISSERTATION SUBMITTED TO THE FACULTY OF THE INSTITUTE FOR MOLECULAR ENGINEERING IN CANDIDACY FOR THE
More informationThe electronic structure of materials 1
Quantum mechanics 2 - Lecture 9 December 18, 2013 1 An overview 2 Literature Contents 1 An overview 2 Literature Electronic ground state Ground state cohesive energy equilibrium crystal structure phase
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
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 informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
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 informationElectronic Structure Calculations, Density Functional Theory and its Modern Implementations
Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011 Outline 1 of Atomistic calculations
More informationSi-nanoparticles embedded in solid matrices for solar energy conversion: electronic and optical properties from first principles
Si-nanoparticles embedded in solid matrices for solar energy conversion: electronic and optical properties from first principles S. Wippermann, M. Vörös, D. Rocca, T. Li, A. Gali, G. Zimanyi, F. Gygi,
More informationElectronic structure calculations: fundamentals George C. Schatz Northwestern University
Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
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 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 informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationChemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase
Supplementary Information Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Jia Chen, Ye-Fei Li, Patrick Sit, and Annabella Selloni Department of Chemistry,
More informationYuan Ping 1,2,3*, Robert J. Nielsen 1,2, William A. Goddard III 1,2*
Supporting Information for the Reaction Mechanism with Free Energy Barriers at Constant Potentials for the Oxygen Evolution Reaction at the IrO2 (110) Surface Yuan Ping 1,2,3*, Robert J. Nielsen 1,2, William
More informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More informationAdditional background material on the Nobel Prize in Chemistry 1998
Additional background material on the Nobel Prize in Chemistry 1998 The Royal Swedish Academy of Sciences has decided to award the 1998 Nobel Prize in Chemistry with one half to Professor WALTER KOHN,
More informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute for Microstructure Physics OUTLINE LECTURE I Phenomena to be described by TDDFT Some generalities on functional theories LECTURE
More informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
More information5.61 Physical Chemistry Final Exam 12/16/09. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry
5.6 Physical Chemistry Final Exam 2/6/09 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.6 Physical Chemistry Final Examination () PRINT your name on the cover page. (2) It
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 information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
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 informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
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 informationIntroduction to Vibrational Spectroscopy
Introduction to Vibrational Spectroscopy Harmonic oscillators The classical harmonic oscillator The uantum mechanical harmonic oscillator Harmonic approximations in molecular vibrations Vibrational spectroscopy
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
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 informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationLinear response time-dependent density functional theory
Linear response time-dependent density functional theory Emmanuel Fromager Laboratoire de Chimie Quantique, Université de Strasbourg, France fromagere@unistra.fr Emmanuel Fromager (UdS) Cours RFCT, Strasbourg,
More informationNotes on Density Functional Theory
Notes on Density Functional Theory 1 Basic Theorems The energy, E, of a system with a given Hamiltonian H is a functional of the (normalized, many-particle) wave function Ψ. We write this functional as
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 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 informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationTime-dependent density functional theory (TDDFT)
04/05/16 Hands-on workshop and Humboldt-Kolleg: Density-Functional Theory and Beyond - Basic Principles and Modern Insights Isfahan University of Technology, Isfahan, Iran, May 2 to 13, 2016 Time-dependent
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 informationLarge scale GW calculations
M. Govoni, UoC & ANL T.A. Pham, LLNL Nguyen H.-Viet, Vietnam Acad. of Sci.&Technology D. Rocca, U.Lorraine, France H. Wilson RMIT Univ., Australia I.Hamada, P. Scherpelz, NIMS, UoC Japan & UoC Large scale
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 informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationAdvanced TDDFT II. Double Excitations in TDDFT
Advanced TDDFT II. Double Excitations in TDDFT f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York First, quick recall of how we get excitations in TDDFT: Linear
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
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 informationAb-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle
Ab-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle Alain Delgado (a,b), Stefano Corni (b), Carlo Andrea Rozzi (b) Stefano Pittalis
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 informationExercises : Questions
Exercises 18.05.2017: Questions Problem 1 where Calculate the following commutators: a) [ Ĥ, ˆp ], b) [ Ĥ, ˆr ], Ĥ = 1 2m ˆp2 + V ˆr), 1) ˆp 2 = ˆp 2 x + ˆp 2 y + ˆp 2 z and V ˆr) = V ˆx, ŷ, ẑ) is an arbitrary
More information2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density
1. SIESTA DFT In a Nutshell Introduction to SIESTA Atomic Orbitals Capabilities Resources 2. TranSIESTA Transport in the Nanoscale - motivation Boundary Conditions: Open systems Greens functions and charge
More informationTopological Physics in Band Insulators IV
Topological Physics in Band Insulators IV Gene Mele University of Pennsylvania Wannier representation and band projectors Modern view: Gapped electronic states are equivalent Kohn (1964): insulator is
More informationQuantum 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 information5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.61 Physical Chemistry Exam III (1) PRINT your name on the cover page. (2) It is suggested that you READ THE ENTIRE EXAM before
More informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
More informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
More information