Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids
|
|
- Randall Price
- 6 years ago
- Views:
Transcription
1 Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids Kevin P. Driver, Ph.D. Defense, February 1, 011 Outline 1) Introduction: Quantum mechanics describes matter but is difficult to solve. ) Computational Methods: Density functional theory(dft) and quantum Monte Carlo (QMC) provide accurate and computationally affordable solutions to large quantum mechanical systems. 3) Accuracy: Expensive QMC and hybrid DFT provide benchmarks for less expensive DFT methods: Tests of DFT's accuracy provides a route to deal with DFT weaknesses. 4) Applications: Silicon self-interstitial defects and silica are difficult and important real world problems that simultaneously exhibit important technical tests for ab initio methods. 5) Conclusions: In addition to QMC, hybrid DFT(HSE) is a benchmarking standard for DFT. Results establish the feasibility of QMC and HSE to study and benchmark complex materials significantly beyond the scope of previous calculations.
2 Quantum Mechanics Describes Matter But No Exact Solution Quantum Mechanics is a non-deterministic theory that relates wave properties of matter to its energy: A particle is described by a wavefunction (prob amp) r, t in states of discrete energy, E with a probability of being found at point r, at time t of r, t Mathematical description is given by Schrödinger's Wave Equation: N where H is the full Hamiltonian operator: = H 1 mi =E H i i=1 is a N particle (many body) wavefunction. kinetic energy N Zi Z j i j r i r j potential Coulomb energy Coulomb term describes each electron interacts with all others in a system: many-body problem The coulomb term is inseparable, which means simple analytic solutions are not possible. Additionally, each electron has 3 dimensions; total number of dimensions for solids is large! mn >> me: allows separation of electron and nuclear degrees of freedom (Born-Oppenheimer)
3 Many-body electron interactions required for accuracy 1) Electron Exchange Interactions (Fermi Correlation/Hund #1) The interaction of electrons via the Pauli exclusion principle Largely responsible for the shape/volume of matter Fermion's obey Pauli r 1, r = r, r 1 ) Electron Correlation Interactions (Coulomb correlation/hund #) Coulomb interactions cause electrons to stay out of each others way. e e In solids, there are Avogadro's number of electrons interacting within essentially continuous bands of quantum states. Is there a clever way to proceed?
4 Computational Methods: Numerical (Approximate) Solutions to Schrödinger's Equation
5 Hartree Fock: Explicit Exchange, Averaged Correlation Method: consider each electron separately interacting with the average repulsion (mean field potential) of all others; allows Schrödinger's equation to be broken down into a set of one electron equations. Pros: Slater determinant wavefunction form obeys Pauli exclusion: exchange is exactly computed. Cons: Correlation is only accounted for by the average effect of repulsion; also includes self interaction. Mean field potential: U = e dr ' el n r ' r r ' 1 r 1 r 1 electron example: r, r = 1 n r = e i r 1 r r i i r are one electron wavefunctions = 1 r 1 r r 1 1 r includes permutation operator, P1 The electron potential terms in the Hamiltonian become orbital dependent operators: V HF =[ J K ] j r 1 = dr i r r r 1 * j r 1 dr i r i r 1 r r 1 j r Coulomb operator Hartree term V H Guess form of i r construct Self consistent Field solution converges VHF exchange operator V X =T V V H ion HF solve new i r H =E
6 Density Functional Theory: Approximate Exchange and Correlation =E Tot =E unique [ n r ] Hohenberg Kohn: GS H GS GS GS There is a 1 to 1 mapping from psi to n; ngs(r) uniquely determines the GS psi and E; it's the true ngs. Real system Electrons fully interact density, n Energy, E e Kohn Sham Exact mapping e e Fictitious model system Electrons treated independently; interact via Veff same density, if Veff is exact. same energy, E, but functional of n(r) 0 e R, V E TOT [ n]=t [ n] E ion [ n] E H [ n] E XC [ n] known form real unknown e e e e r, V i non interacting KE eff EXC includes all the information that is not known: exchange correlation energy and T T0 Minimizing ETOT[n] wrt density variation gives effective one particle Kohn Sham eqns: 1 H i =[ V eff r ] i r = i r where V eff r =V ion V H V XC DFT is exact for universal V, but it is unknown; VXC must be approximated. XC Guess i r construct n r = e i r Self consistent solution converges n i construct V eff new i r solve H i r = i r
7 Quantum Monte Carlo (QMC): Explicit Exchange and Correlation What is Monte Carlo? An efficient way of solving many dimensional integrals (mean value theorem). Evaluation: Randomly sample the integrand and average the sampled values. Why use Monte Carlo? Conventional integration methods (e.g. Simpson's rule) use a mesh of points and error in the result 1 falls off increasingly slow with the dimension of the problem. Error ~ Statistical error from Monte Carlo is independent of dimension. N What is Quantum Monte Carlo? Solution to the Schrödinger equation (SE) using Monte Carlo integration with a many body wavefunction. An alternative to DFT when accuracy is paramount. Explicit many body method including correlation and exchange from the outset. QMC comes in two forms: 1) Variational MC and ) more accurate Diffusion MC Basic outline of a DMC calculation: DFT Trial wavefunction Variational QMC Optimize wavefunction = Jastrow factor Orbital Determinant Correlation Exchange Diffusion QMC Project out Ground State
8 VMC: provides best solution for a fixed wavefunction form. Compute energy via variational principle evaluated using Monte Carlo integration. Sample configurations {R} according to the probability density function (Metropolis algorithm). Evaluate E for each sampled configuration and average the values. L E vmc = R PDF [ ] H R R dr 1 N MC local energy E L N MC E L Ri ± i=1 1 N MC E L E L The VMC Slater Jastrow wavefunction form: R =e J R D R where J(R) is the Jastrow correlation factor and D(R) is a Slater determinant J(R) is a polynomial expansion in particle separation: body (e n, e e) and 3 body (e e n) terms D(R) is made up of single particle electron orbitals (from DFT calculations). Jastrow parameters are optimized by minimizing the energy or variance of the energy. Iterate VMC and VMC optimization until energy change is below desired statistical error. equilibration loop statistics accumulation loop VMC Algorithm random configuration accept trial move reject update e positions accept trial move reject update e positions compute EL optimize Jastrow
9 DMC: projects best solution for fixed wavefunction nodes V(x) DMC transforms Schrödinger's equation to imaginary time diffusion/rate equation ( it): i t t =H =[ 1 V R ] Recast in terms of a Green's function propagator with time step that projects out ground state relative to damped excited states. N 0 x x E DMC = 1 N step N config set up configurations N step N config j=1 k=1 GS = lim N exp H initial 0 j=1 is represented statistically by a population of e configurations which, in each time step, undergo a random diffusion and E L R k, j potential dependent rate(birth/death) processes. DMC the Green's function is the PDF; Possible to sample negative regions. Nodes are fixed to ensure positive PDF. Allow e 's to diffuse for small time steps obeying fixed node and weighted Metropolis. compute rate process probability accumulate averages over all configs increment time step
10 Accuracy Thesis broadly asks two general questions: 1) How accurate are our methods? Efficiency accuracy trade-off. ) Can we solve some important, practical problems?
11 Strengths, Weaknesses, and Approximations in DFT Strength: DFT works great for many materials (thousands of calculations published). Weaknesses: In general, DFT cannot do: 1) Excited states ) Band Gaps 3) van der Waals attractions 4) Energy differences between structures with vastly different interatomic bonding. Source of DFT's weaknesses: el 1) Self interaction error in mean field approach: U = e dr ' n r ' r r ' ) Approximations to exchange-correlation functional (LDA, GGA, meta-gga, Hybrids) Functional performance is sometimes unreliable; There is skepticism amongst the zoo. Jacob's ladder classification systematically includes more and more physical information. Higher rungs of the ladder are not always better. Must find ways to determine the accuracy of these methods. Benchmark methods are needed to ultimately determine the appropriate functional for the job. 3) Pseudopotentials: Construct from atomic calculations; compare tests with all electron before use. ik r Basis sets: Planewaves are easy to converge for solids, but not molecules n r = c nk e k Gaussians are troublesome to converge, but efficient for molecules (and hybrid functionals).
12 Jacob's Ladder Approximation of XC functionals E XC [ n]= dr n r xc [ n r ] where xc is the exchange correlation energy per particle nonlocal nonlocal n, semilocal orbitals semilocal local Hybrid functionals are times more expensive than lower-rung (local, semilocal) functionals. HSE hybrid functional has performed very well for band gaps and properties of simple solids. Performance of this ladder of functionals has not been studied for large, complex solids. QMC benchmarks help determine which functional is most accurate for a given system/property.
13 Strengths, Weaknesses, and approximations in QMC Strengths: Accurate and reliable because XC is computed explicitly. Weaknesses: Computational expense requires approximations and xDFT CPU time 1) Statistical (VMC & DMC) error ~1/sqrt(N) Tunable approximations Increase number of samples (MC steps). J R ) Fixed VMC wavefunction form limits ability to estimate E of GS. R =e D R Expand J in more terms/parameters until energy converges; converge DFT orbitals. 3) Interpolation of PW orbitals (VMC & DMC) PW orbitals are grid interpolated for factor N speed up. Decrease grid size for better accuracy. 4) DMC time-step (DMC) determines accuracy at which wavefunction is sampled Reduce until energy converges 5) Finite-size (VMC & DMC) periodic simulation cell causes systematic error in energy 1 k-point/cell in QMC: use DFT k-point correction or larger supercells in QMC Self-interaction caused by Ewald summation. Fixed by restoring e-e interaction in cell (MPC). 6) Finite population size (DMC) causes bias in energy Increase configuration sample size 1) Pseudopotential necessary for solids Non-tunable approximations Compare results with different pseudopotentials;use DFT all-electron correction ) Pseudopotential locality calculate nonlocal channels with variational lattice DMC 3) Fixed-node apply backflow transformation to electron coordinates, or try different orbitals.
14 Applications 1) Silicon defects ) QMC study of high pressure silica phases 3) Hybrid functional study of silica Semiconductors: Silicon Defects Geophysics: High pressure silica
15 Application 1: Silicon defects affect semiconductor device performance Silicon chips require photographic and chemical processing. Ion implantation induced interstitials precipitate silicon self interstitial defects. Annealing causes defects to coalesce into extended defects. Extended defects broaden dopant profiles, limiting device size. Direct experimental detection is not possible; Watch diffusion of radio isotope profiles to estimate activation energy. Need accurate simulations to understand Si defect formation and migration energies. Defect Evolution (Tight Binding MD/DFT) Single Interstitials Tri Interstitials Anneal Di Interstitials Anneal
16 DFT and QMC Study Silicon Self Interstitial Defect Formation Energy 16 atom Supercell at L point Defect Type LDA lies roughly ev below QMC; GGA 1.5 ev below QMC. Likely due to self-interaction plus inability of LDA and GGA to compute energy differences between distorted structures. X and H defect formation energies are degenerate and T has highest by about 0.5 ev in 64-atom cell. Observe systematic improvement of functionals, consistent with Jacob's ladder. QMC results agree with less accurate results of Leung et al.
17 DFT and QMC Self Diffusion Barrier (Migration Energy) QMC and GGA results for self-diffusion from X to H to T to X. Lowest energy barrier from X to H is similar in both QMC and GGA (~300 mev). X to T barrier is twice as large in QMC as in GGA. Activation energy for diffusion is formation + migration energy Experiments estimate ev QMC predicts 5.(1) ev GGA predicts 4.15 ev
18 Di and Tri interstitials LDA lies roughly 3 ev below QMC; GGA ev below QMC. Observe systematic improvement of functionals, consistent with Jacob's ladder. A I is most stable di-interstitial and I3A is most stable tri-interstitital. GGA improves for less distorted defect.
19 Tunable QMC approximation tests DMC time step convergence Finite size convergence Bulk Si, L point X defect, L point DMC time step converged by 0.1 Ha-1 Finite size errors converged by 64-atom cell VMC energy converged by Jastrow polynomial order of 550. (order 5 in e-n expansion, order 5 in e-e expansion, no e-e-n expansion) VMC, X defect, 16 atom cell
20 Non Tunable QMC approximation tests DMC Formation energy does not vary with choice of DFT orbitals. DMC and VMC energies are equivalent with HF and LDA pseudopotentials by 3-atom supercell
21 Application Compression of silica Application: QMC and:dft Study High Pressure Silica Simplest of Earth's silicates; ubiquitous component of Earth; Complex phase transitions with pressure. Quartz to stishovite is a four fold to six fold coordination change. nd Stishovite undergoes a order transition to a CaCl structured phase. CaCl structure transforms to a PbO structure, which is stable to the core. Goals of this work: Explore feasibility of using QMC for high pressure/temperature properties of silica. Compute thermal equations of state and phase diagram (Combine QMC with DFT phonons). Compute thermodynamic properties of silica.
22 Compute static energies of silica phases at several volumes using QMC QMC Static Energy vs Volume Static energy is dominant crystal energy QMC is most accurate method for solids DFT 1)Choose XC functional and pseudopotential. )Relax crystal structures. 3)Compute energy and wavefunction. DFT QMC 1)Use DFT's relaxed crystal structures. )Optimize DFT wavefunction. 3)Compute energy stochastically, fit E(V) T=0 K Transition Pressure (GPa) Quartz Stishovite CaCl apbo Experiment QMC DFT(WC) 6 to 7 ~90 6.3(0.14) 88(8).1 86
23 Compute energy due to thermal vibrations (phonons) using DFT Quartz Phonon Dispersion (P=0 GPa) Frequency (cm 1) Compute vibrational free energy with DFT Currently too costly for QMC Vibrational energy is small. Typically well described in DFT Phonons computed within Quasi harmonic approx. = V Solve equations of motion for frequencies using small displacements. Dispersion data from Burkel, et al. Physica B, 63 64, pp (1999). Phonons computed from Helmholtz vibrational energy: F =E static V F vibration V, T = ktln Z 1 modes Z = exp i / kt i Implicit T dependence because of V dependence modes ℏ i V kt ln 1 exp ℏ i V / kt i=1 i=1 static energy F QHA = E 0 V zero point E 1 thermal vib E is force constant matrix; sum times phase factor gives dynamical matrix K F QHA=E 0 K u Phonon frequencies are eigenvalues of dynamical matrix DFT End of the day: F, T =F QMC T =0 F, T Static vib
24 Thermal Equations of State Fit F(V) with Vinet EoS then compute P= F V T Thermal EoS determines fundamental thermodynamic parameters and phase relations. QMC improves agreement with experiment in each phase: quartz, stishovite, PbO. Only small number of measurements for PbO, making QMC most accurate available. QMC gives internal estimate of error. Gray shading indicates one standard deviation of statistical error.
25 Silica Phase Boundaries G=F PV G 1 P T, V 1 =G P T,V Gibbs FE equal at equilibrium F 1 P T P T V 1 P T =F P T P T V P T Equilibrium phase boundaries computed from Gibbs free energies. Metastable quartz stishovite transition measured with thermocal or shock. QMC agrees well. QMC CaCl PbO transition lies between two measurements. No free silica at base of mantle. DFT(WCGGA) boundary 4 GPa too low for quartz stishovite and within statistical error of QMC for CaCl PbO
26 Thermal Expansivity Thermal Expansivity (10 5K 1) = 1 V V T P QMC and DFT(WCGGA) temperature and pressure dependence of. QMC shows best agreement for quartz. Experimentally, quartz appears at 846 K we only consider quartz. QMC and DFT show good agreement with stishovite measurements. Temperature (K) PbO curves are the best available.
27 Heat Capacity C p= H T P QMC and DFT(WCGGA) temperature and pressure dependence of Cp. Cp/R QMC and DFT results are nearly identical. Good agreement with experiment for quartz and stishovite. Temperature (K) QMC PbO curves are best available.
28 ... And many more thermodynamic properties Room temperature comparison to experiment for all computed thermodynamic quantities: Agreement with QMC is very good for all computed quantities
29 First QMC Calculation of elastic constants by straining the lattice Strain energy density relation: Stishovite DFT(WC) Strain the lattice: R ' =[ I ] R For a volume conserving strain: E 1 = c ijkl ij kl V 1500 processor hours c ijkl = 1 E V Elastic constants obtained from curvature of energy strain curve Double well at 80 Bohr3 indicates elastic instability of stishovite CaCl becoming more stable that stishovite under pressure Feasibility of Elastic Constants in QMC Elasticity is a tough problem for QMC: energy differences ~ ev Extremely expensive to get accurate error bars for large (100 atom) systems Through parallel computation on large supercomputers, it's possible to succeed.
30 QMC Energy vs. Strain Curves: The Brute Force Method 1) Take optimized input structures from DFT(WC) (we can't do forces in QMC yet) ) Run QMC on thousands of processors for many days until error bars are sufficiently small for each structure. VMC (500,000 hrs) DMC (additional 1.3 million hrs) QMC at this accuracy level is 100 times more expensive than DFT. QMC error bars must be made much smaller than the strain energy differences. VMC error bars decrease twice as fast as DMC error bars. Highest pressure curves are most difficult to fit and require smallest error bars (work on high pressure curves is still in progress).
31 c11 c1 (GPa) Stishovite Shear Modulus c ij = 1 E V Shear constant: 3 million CPU hours (1 summer). Total QMC silica project: 10 million CPU hours. Stishovite to CaCl transition is driven by instability in the elastic shear modulus. VMC modulus computed at several pressures and DMC checks at endpoints. Shear modulus computed from strain energy relation (brute force 1000 CPU cost of DFT) All methods roughly agree, with the shear modulus going unstable around 50 GPa.
32 Application 3: Hybrid DFT Study of Silica Quartz Stishovite Transition Hybrid functionals are potentially as accurate as QMC (particularly HSE) Hybrid functionals have not been well tested for large, complex systems due to expense associated with exact exchange. HSE improves efficiency by screening exchange. Hybrid functionals typically require localized (Gaussian) basis sets; difficult to converge
33 QMC and hybrid comparison of zero pressure volume Procedure Compute E at several volumes at zero temperature. Use Vinet EoS to fit E(V) for quartz and Stishovite All following properties are then computed from E(V) Compare dependence on codes, pseudopotentials, basis set, and XC functional. CRYSTAL: all-electron, Gaussian basis set ABINIT: norm-conserving, nonlocal pseudopotential (based on atomic calc), planewaves VASP: PAW pseudopotential (based on all-electron solid calc), planewaves GGA's: PBE, PW91, PBEsol, WC Hybrids: B3LYP, PBE0, HSE LDA better than GGAs for structural properties; HSE performs as well as LDA. No Jacob's Ladder trend in performance
34 QMC and hybrid comparison of zero bulk modulus K = V P V T =V F V LDA better than GGAs for elastic properties HSE performs better for stishovite than quartz. No Jacob's Ladder trend in performance
35 QMC and hybrid comparison of quartz stishovite transition pressures GGAs better than LDA for energy differences for two systems with very different structures. PBE, PW91, HSE, and PBE0 all predict the correct transition pressure. Overall conclusion: LDA is better than GGA for structural and elastic properties. GGA is better than LDA for energy differences HSE is good for all properties and the best compromise of accuracy and speed for benchmarks.
36 Conclusions Quantum mechanics describes matter, but its equations are difficult to solve for matter. Mean-field-based approaches (HF, DFT) provide manageable single-particle equations. QMC offers stochastic solution to the full many-body Schrödinger equation.-vmc uses variational principle for fixed many-body wavefunction form -DMC projects the ground state using a diffusing/branching, statistical wavefunction DFT is efficient, but it's accuracy is limited mainly by the XC functional. QMC is expensive, but computes XC explicitly at a high accuracy level. A benchmark method that compromises accuracy and speed is needed to validate XC functionals Jacob's ladder of functionals indicates hybrid functionals should be most accurate QMC benchmarks of silicon interstitial defects indicate LDA and GGA are not adequate, but HSE performs well. QMC benchmark calculations of silica use DFT phonons to produce QMC-based thermal equations of state and thermodynamic properties and the first elastic constant. LDA, PBE, WC are not adequate to study high pressure phases. QMC provides reliable geophysical implications. Hybrid studies of the quartz-stishovite silica transition indicates the HSE functional performs very well for all silica properties. HSE is a less expensive compromise for benchmark calculations of complex materials.
37 Backup Slides
38 Harmonic approximation to free energy F =E TS= ktln Z Helmholtz Free Energy: Partition function: Z = exp i / kt i 1 3 Assume microstates are harmonic: i = ℏ i, ℏ i, 5 ℏ i,... Partition function for harmonic quantum states: modes Z = exp i=1 F =E 0 static internal 1 1 ℏ i / kt modes s=0 modes i=1 modes exp s ℏ i / kt = i=1 modes ℏ i kt zero point i=1 exp 1 ℏ i / kt 1 exp ℏ i / kt ln 1 exp ℏ i / kt harmonic vibrational contribution (thermal)
39 Quasi harmonic adds anharmonicity the easy way Assume = V is a function of volume only F QHA=E 0 V 1 modes i=1 modes ℏ i V kt i=1 ln 1 exp ℏ i V / kt Notice the thermal term now has both T and V dependence T dependence comes from the partition function/definition of F V dependence comes from frequency By assuming a frequency volume dependence, we add in low order anharmonicity terms Since Fth=Fth(T,V), T and V dependence of all thermodynamic functions can be computed. (Thermal expansion and gruneisen are no longer zero) In DFT: Compute phonons for small displacements at various volumes of a structure Phonons are computed at T=0 (frequencies are T independent) When will this breakdown?
40 Statistical Error Propagation 1) Propagation of error equation from Taylor expanding a function about the mean values of its parmeters: F = i, j F F Ci C j Cov [i, j ] P = i ) Propagation with Monte Carlo Pi Ei E i actual data set Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Fit Vinet Equation to new set of energies and compute property Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Allow random Gaussian fluctuations on QMC energies Ei with stdv Ei Fit Vinet Equation to new set of energies and compute property Fit Vinet Equation to new set of energies and compute property standard deviation of synthetic data sets gives uncertainty in property......
41
42 QMC Locality Approximation Pseudopotentials separate core and valence electrons; effective potential reproduces core. The pseudopotential is made separately for each angular momentum state; must be nonlocal. DMC requires pseudopotentials for efficiency (Z6.5 scaling without, Z3 scaling using Pps) In DMC, nonlocal pseudopotentials create an addition sign problem if left untreated: H eff =K V loc V nloc result of Vnloc on trial wfn Schrodinger's equation becomes f R, = 1 m v f f d The locality approximation assumes H T T E T f V DMC DMC V T T f DMC = T in the last term. Without this approximation, the elements V R can cause G (R,R') to change sign. R e B Going beyond the locality approximation Use variational, lattice regularized DMC (Casula) Silica QMC calculations used this correction scheme. nloc
High 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 informationExploring deep Earth minerals with accurate theory
Exploring deep Earth minerals with accurate theory 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: NCAR, TeraGrid, NCSA,
More informationKevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2.
Quantum Monte Carlo Simulations of a Single Iron Impurity in MgO Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2 1 Department of Earth & Planetary Science University of California, Berkeley
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 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 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 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 informationDFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
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 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 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 informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
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 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 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 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 informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
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 informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationThe Schrödinger equation for many-electron systems
The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3
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 informationMODULE 2: QUANTUM MECHANICS. Principles and Theory
MODULE 2: QUANTUM MECHANICS Principles and Theory You are here http://www.lbl.gov/cs/html/exascale4energy/nuclear.html 2 Short Review of Quantum Mechanics Why do we need quantum mechanics? Bonding and
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
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 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 informationBasics of density-functional theory and fast guide to actual calculations Matthias Scheffler
Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler http://www.fhi-berlin.mpg.de/th/th.html I. From the many-particle problem to the Kohn-Sham functional II. From
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
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 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 informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
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 informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
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 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 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
More informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
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 informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationPhonon calculations with SCAN
Workshop on the SCAN density functional: Fundamentals, practices, and extensions Temple university, Philadelphia May 18th, 2017 Hands-on tutorial 3 Phonon calculations with SCAN Yubo Zhang and Jianwei
More informationComparison of various abinitio codes used in periodic calculations
Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
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 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 informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
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 informationPhonon wavefunctions and electron phonon interactions in semiconductors
Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
More informationDFT: Exchange-Correlation
DFT: Exchange-Correlation Local functionals, exact exchange and other post-dft methods Paul Tulip Centre for Materials Physics Department of Physics University of Durham Outline Introduction What is exchange
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationLecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review
Lecture contents A few concepts from Quantum Mechanics Particle in a well Two wells: QM perturbation theory Many wells (atoms) BAND formation Tight-binding model Solid state physics review Approximations
More informationIntroduction to Density Functional Theory
1 Introduction to Density Functional Theory 21 February 2011; V172 P.Ravindran, FME-course on Ab initio Modelling of solar cell Materials 21 February 2011 Introduction to DFT 2 3 4 Ab initio Computational
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 informationQuantum Monte Carlo wave functions and their optimization for quantum chemistry
Quantum Monte Carlo wave functions and their optimization for quantum chemistry Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France CEA Saclay, SPhN Orme des Merisiers April 2015 Outline
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 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 information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
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 informationQMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters
QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters Idoia G. de Gurtubay TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay. Quantum
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 informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationSelf-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,
More informationQuantum Monte Carlo backflow calculations of benzene dimers
Quantum Monte Carlo backflow calculations of benzene dimers Kathleen Schwarz*, Cyrus Umrigar**, Richard Hennig*** *Cornell University Department of Chemistry, **Cornell University Department of Physics,
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationdynamics computer code, whiccan simulate the motion of atoms at a surface and We coembininewton's equation for the nuclei,
DUI FILE COpy T DTIC ANNUAL LETTER REPORT BY 0. F. ELEC33 20KE AD-A224 069 P.1. FOR N00014-85-K-0442 U Theory of electronic states and formations energies of.defect comp exes, interstitial defects, and
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
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 informationComputational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen
! Computational methods: Coupled Electron Ion Monte Carlo Path Integral Monte Carlo Examples: Electron gas Hydrogen WHO DID THE WORK? Miguel Morales, Livermore Carlo Pierleoni: L Aquila, Italy Jeff McMahon
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
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 informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
More information3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional
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 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 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 informationStructural Calculations phase stability, surfaces, interfaces etc
Structural Calculations phase stability, surfaces, interfaces etc Keith Refson STFC Rutherford Appleton Laboratory September 19, 2007 Phase Equilibrium 2 Energy-Volume curves..................................................................
More informationGeometry Optimisation
Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background
More informationAl 2 O 3. Al Barrier Layer. Pores
Overview of DFT Calculations for Nanoscience and Nanotechnology Oğuz Gülseren Al 2 O 3 Al Barrier Layer Pores Jules Verne A Trip to the Moon, 1865 Nanostructures: a. Contain a countable number of atoms
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 informationAb initio Berechungen für Datenbanken
J Ab initio Berechungen für Datenbanken Jörg Neugebauer University of Paderborn Lehrstuhl Computational Materials Science Computational Materials Science Group CMS Group Scaling Problem in Modeling length
More informationMinimal Update of Solid State Physics
Minimal Update of Solid State Physics It is expected that participants are acquainted with basics of solid state physics. Therefore here we will refresh only those aspects, which are absolutely necessary
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationNew applications of Diffusion Quantum Monte Carlo
New applications of Diffusion Quantum Monte Carlo Paul R. C. Kent (ORNL) Graphite: P. Ganesh, J. Kim, C. Park, M. Yoon, F. A. Reboredo (ORNL) Platinum: W. Parker, A. Benali, N. Romero (ANL), J. Greeley
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationIndependent electrons in an effective potential
ABC of DFT Adiabatic approximation Independent electrons in an effective potential Hartree Fock Density Functional Theory MBPT - GW Density Functional Theory in a nutshell Every observable quantity of
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 informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
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 informationIntegrated Computational Materials Engineering Education
Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &
More informationThermodynamics of finite quantum systems
Thermodynamics of finite quantum systems Shell structure in finite quantum systems Erice summer school, July 25-30, 2010 Klavs Hansen Department of Physics, University of Gothenburg Topics Metal clusters
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 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 informationFixed-Node quantum Monte Carlo for Chemistry
Fixed-Node quantum Monte Carlo for Chemistry Michel Caffarel Lab. Physique et Chimie Quantiques, CNRS-IRSAMC, Université de Toulouse e-mail : caffarel@irsamc.ups-tlse.fr. p.1/29 The N-body problem of Chemistry
More informationab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 30, 2009 This afternoon s plan
ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 3, 29 This afternoon s plan introductory talk Phonons: harmonic vibrations for solids Phonons:
More information