DFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
|
|
- Timothy Allen
- 5 years ago
- Views:
Transcription
1 DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE
2 Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics, molecules and clusters versus periodic systems - scalar relativity - solution of the eigenvalue problem
3 Equations to be solved Hohenberg-Kohn theorems Electron density Ψ,,, ; { } [ ], = Ψ,,, ; { } Ψ,,, ; { } Variational principle: = 0,,?
4 Equations to be solved = Ψ Ψ = Ψ,, Ψ,, Ψ,,, Ψ,, + + Ψ,,, Ψ,, = =, +, with one- and two-particle density matrices:, = Ψ,, Ψ,,, = Ψ,, Ψ,,, Problem: kinetic energy and electron-electron interaction energy as a functional of ( ) --?
5 Equations to be solved Kohn-Sham approach: + = + ( ) + [ ] where functions represent the density: = = 0 1 = 0 (Generalized) Kohn-Sham equations: = + + +, [{ }], =
6 = Solving Kohn-Sham equations + + +, [{ }], =
7 Self-consistent field method Initial guess: e.g. density Update potential n 0 V j+1 Calculate potential V 0 Construct Hamiltonian n j+1 = n j +f(δn) Update density / density matrix (mixer, preconditioner) No Obtain new eigenvalues, eigenvectors, and density n j Converged? Yes
8 Convergence When is the calculation finished? Answer: When property does not change Very loose Convergence threshold Very strict Typical properties to converge: Change of the density D n Change of the sum of eigenvalues Change of the total energy Options differ between codes
9 Self-consistent field method (SCF) Initial guess: e.g. density Update potential n 0 V j+1 Calculate potential V 0 Construct Hamiltonian n j+1 = n j +f(δn) Update density / density matrix (mixer, preconditioner) No Obtain new eigenvalues, eigenvectors, and density n j Converged? Yes
10 Initial guess Importance of the initial guess Electronic structure exhibits multiple minima Qualitatively different behavior depending on guess Common issue for magnetic solids 1 O 2 3 O 2
11 Initial guess Importance of the initial guess Electronic structure exhibits multiple minima Qualitatively different behavior depending on guess Common issue for magnetic solids How do we find the right 1 O 2 initial guess? 3 O 2
12 Initial guess Random initialization Traditional approach (for bandstructure codes)
13 Initial guess Superposition of spheric atomic densities Most common (for bandstructure codes)
14 Initial guess Extended Hückel Model [1] Semi-empirical method: - Linear combination of atomic orbitals - Hamiltonian parameterized: - Density from orbitals: [1] R. Hoffmann, J Chem. Phys (1963), 1397; [2] R. S. Mulliken, J. Chem. Phys. (1946) 497; [3] M. Wolfsberg and L. Helmholtz, J. Chem. Phys (1952), 837
15 Self-consistent field method (SCF) Initial guess: e.g. density Update potential n 0 V j+1 Calculate potential V 0 Construct Hamiltonian n j+1 = n j +f(δn) Update density / density matrix (mixer, preconditioner) No Obtain new eigenvalues, eigenvectors, and density n j Converged? Yes
16 Density Update Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j D n Naive Mixing Example: H 2, d=1.5å, projected in 1 dimension
17 Density Update Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j D n Naive Mixing Example: H 2, d=1.5å, projected in 1 dimension Typically does not work - Best case: Osciallating, non-converging results - Worst case: Bistable, apparently converged solution
18 Density Update Linear Mixing Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j ad n Example: H 2, d=1.5å, projected in 1 dimension
19 Density Update Linear Mixing Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j ad n Example: H 2, d=1.5å, projected in 1 dimension
20 Density Update Linear Mixing Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j ad n Example: H 2, d=1.5å, projected in 1 dimension
21 Density Update Linear Mixing Calculate new density from Kohn-Sham orbitals Replace old density by new density n j 1 = n j ad n Example: H 2, d=1.5å, projected in 1 dimension system dependent No clear recipe to choose ideal α Guaranteed convergence
22 Density Update Pulay Mixing [1] (Direct Inversion of Iterative Subspace) Define residual: Change between input and output density Predict residual of next step from previous steps Idea: obtain b minimizing norm Pitfall: Only use limited number of previous steps [1] P. Pulay, Chem. Phys. Lett. 73, 393 (1980).
23 Density Update Pulay Mixing [1] (Direct Inversion of Iterative Subspace) Construct optimal density minimizing residual norm, form new density Learning from previous steps Ideal case: convergence independent of a Rules of thumb: - a = for insulators - a = 0.05 for metals [1] P. Pulay, Chem. Phys. Lett. 73, 393 (1980).
24 Density Update Pulay Mixing [1] (Direct Inversion of Iterative Subspace) Construct optimal density minimizing residual norm, form new density Learning from previous steps Ideal case: convergence independent of a Rules ofmuch thumb: faster than linear mixing - a = No 0.2 guarantee for insulators for congerence - a = 0.05 for metals [1] P. Pulay, Chem. Phys. Lett. 73, 393 (1980).
25 Density Update Kerker (charge) preconditioning Problem: Charge Sloshing Charge bounces back and forth (Lack of) balance between long-range and short-range charge re-distribution Typical for metals Often encountered in defective systems Localized defects (vacancies, interstitials) Surfaces, Slabs, Thin Films
26 Density Update Charge preconditioning Solution: Make α depend on ( ) Concept: Treat density change in Fourier space: = + Δ = + Δ = 1 Ω Δ Damp long-range oscillations: q 0 : Screening parameter = + q >> q 0 : G a Normal mixing q << q 0 : G aq 2 /q 0 2 Damping
27 Density Update Charge preconditioning = + q 0 : Screening parameter 8 layer Al slab, PBE calculation no a priori known ideal choice for q 0 Reasonable guess from Thomas-Fermi:
28 Density Update A special problem: Metallic systems Metallic properties determined by Fermi surface Fermi-surface: Collection of k-points at which e = e Fermi At e Fermi occupation changes from 1 to 0 Problem: Frequent level switching during SCF
29 Density Update A special problem: Metallic systems Solution: Replace step function by a more smooth function Fermi [1] Gaussian [2] Methfessel-Paxton [3] [1] N. Mermin, Phys. Rev.137, A1441 (1965).[2] C.-L. Fu, K.-H. Ho, Phys. Rev. B 28, 5480 (1983). [3] M. Methfessel, A. Paxton, Phys. Rev. B 40, 3616 (1989).
30 Density Update A special problem: Metallic systems Numerical Implications Level switching is damped Much faster convergence Total energy now depends on σ, no longer variational Back-extrapolation possible Pitfall: Adsorption of molecules Physical implications: States become fractionally occupied Cu (111), 5 layer slab, Gaussian smearing, PBE calculation, 12x12x1 k- points In DFT, well defined as ensemble average Broadening is not electron temperature (except: Fermi)
31 Summary SCF Ground state electron density determined iteratively Initial guess needs initial thought, can change results qualiatively Density update by Linear mixing (slow) Pulay mixing Convergence acceleration by Preconditioner Broadening of states
32 Representing KS states: Basis sets
33 The tool for this workshop: FHI-aims (NAOs)
34 Our choice: Numeric atom-centered basis functions
35 Our choice: Numeric atom-centered basis functions
36 Constructing a basis-set library
37 Iterative selection of NAO basis functions
38 Result: Hierarchical basis set library for all elements
39 Accuracy: (H 2 O) 2 hydrogen bond energy
40 Basis-set superposition error (BSSE) By construction, our NAO basis sets are close to completeness for very small BSSE for atomic dimers
41 Basis-set superposition error (BSSE)
42 Using numeric atom-centered basis functions
43 Numeric atom-centered basis functions: Integration
44 Overlapping atom-centered grids: Partitioning of unity
45 Integration in practice: Large systems - small errors!
46 Hartree potential (electrostatics): Overlapping multipoles
47 Electrostatics: convergence with
48 Periodic systems
49 Structure optimization
50 Global structure search Born-Oppenheimer energy surface can contain several minima Constitution isomery Configuration isomery Conformation isomery Polymorphism System in equilibrium is given by ensemble average over all minima Often dominated by global minimum (but watch out for tautomers) 50
51 Global structure search Methods to find the global minimum: o Basin Hopping o Molecular dynamics: o Simulated annealing o Minima Hopping o Metadynamics o Cluster expansion o Genetic algorithm o Diffusion methods
52 Local structure optimization Once we have a reasonable guess, find closest minimum Different approaches possible: Mapping of the potential energy surface Requires a lot of calculations Typically only for dynamics a CCO d CC 52
53 Local structure optimization Once we have a reasonable guess, find closest minimum Different approaches possible: Mapping of the whole potential energy surface Gradient free method: e.g., simplex Gradient-based methods Calculate gradient (a.k.a. forces ) =
54 Total energy gradient Search for minimum by following the gradient = Ψ Ψ + Ψ + Ψ = + = + First term vanishes Pulay forces Second term survives for atom-centered basis functions (vanish when basis set approaches completeness)
55 Total energy gradient Search for minimum by following the gradient = Ψ Ψ + Ψ + Ψ Additional contributions from atom-centered approximations Multipole expansion Relativistic corrections (Integration grids) All straightforward but lengthy Once we have the Force, how do we find the minimum?
56 Geometry update Steepest descent Follow negative gradient to find minimum Step-length α variable Guaranteed but slow convergence Oscillates near minimum Not suitable for saddle points Improved versions exists Line minimization (optimal step length) Conjugated gradient [1] [1] M. Hestenes, E. Stiefel (1952). "Methods of Conjugate Gradients for Solving Linear Systems"
57 Geometry update (Quasi)Newton Methods Approximate PES by quadratic function Δ + Δ + Δ Δ... Hessian Find minimum: Newton: calculate exact H Expensive! Cheaper method needed [1] J. Nocedal and S. J. Wright, Numerical optimization (Springer, 2006)
58 Geometry update (Quasi)Newton Methods Approximate PES by quadratic function Δ + Δ + Δ Δ... Hessian Find minimum: Newton: calculate exact H Quasi-Newton: approximate H Update as search progresses [1] [1] J. Nocedal and S. J. Wright, Numerical optimization (Springer, 2006)
59 Geometry update (Quasi)Newton Methods Guess initial Hessian Simple choice: Scaled unit matrix Chemically motivated choice (e.g.: Lindh [1]) streching bending torsion k parameterized [1] Update as search progresses [2] [1] R. Lindh et al., Chem. Phys. Lett. 241, 423 (1995). [2] J. Nocedal and S. J. Wright, Numerical optimization (Springer, 2006)
60 Geometry update (Quasi)Newton Methods Initial Hessian is critical for performance May even influence result (combination with convergence thresholds) Diagonal Hessian Model Hessian (Fischer) Courtesy of Elisabeth Wruss Lower in energy
61 Challenges of Quasi-Newton Methods Soft degrees of freedom can cause large ΔR Step control needed: Line search method: If new point is worse than old, interpolate Trust radius method Enforce upper limit for ΔR Evaluate quality of quadratic model Adjust ΔR max based on q
62 Conclusions (Global optimization: PES feature-rich, methods to find global minima exist) Local geometry optimization: Follow gradient Hellman-Feynman from moving potentials Pulay from moving basis functions + additional terms Quasi-Newton method de-facto standard Requires approximation and update of Hessian Step control by line search or trust radius method
63 Vibrations in the harmonic approximation
64 Why calculating vibrations? Vibrations give important information about the system: Classification of stationary points (minimum / saddle point) If saddle-point: Provides search direction Thermodynamic data Zero-point energy Partition sum Finite-temperature effects Connection to experiment: Infrared intensities: derivative of dipole moment Raman intensities: derivative of polarizabilty
65 How good is the harmonic approximatino? Morse-potential: Harmonic oscillator: Harmonic osc. Morse:
66 How good is the harmonic approximation? Morse-potential: Harmonic oscillator: Harmonic osc. Harmonic osc. Reasonable model for small displacements (~10% of bond length) Maximum displacement: Morse: About 0.1Å at room temperature
67 Vibrations Expand Energy in Taylor series: Hessian H
68 Calculating Vibrations Solve Newton s equations of motion: Exponential ansatz: Leads to generalized eigenvalue problem: Interpretation of results: Single negative ω: Transition state ~ /
69 From vibrations to thermodynamics Partition function Free energy for finite temperature Hessian from geometry optimization is not sufficent Analytic second derivative using perturbation theory [1] Numerical differentiation Computationally very expensive Contains parameter (displacement), needs to be checked carefully Often single displacement is not sensible to sample all vibrations [1] S. Baroni et al., Rev. Mod. Phys. 73, 515 (2001).
70 Conclusions Often calculated in harmonic approximation Yield information about stability of geometry Required for temperature effects Anharmonic effects via molecular dynamics
71 Relativity
72 Simple approximation to scalar relativity
73 Simple approximation to scalar relativity
74 Fixing ZORA
75 Atomic ZORA and scaled ZORA in practice
76 Computational scaling: two subproblems
77 CPU time versus system size
78 CPU time versus system size
79 Parallel eigenvalue solvers - the problem
80 Parallel eigenvalue solvers - the problem
81 A massively-parallel dense eigensolver: ELPA
82 The key improvement: two-step trigonalization
83 The key improvement: Two-step trigonalization
84 ELPA: Performance
85 ELPA: Performance (Cray XC30, 2013)
86 Generalized Kohn-Sham DFT: Exact exchange
87 Resolution of identity (RI)
88 Hybrid functionals -- scaling with system size Periodic GaAs, HSE06 hybrid functional
89 Hybrid functionals -- CPU scaling Periodic GaAs, HSE06 hybrid functional
90 Summary
91 Our approach to all-electron simulations: FHI-aims Igor Ying Zhang (FHI) Karsten Reuter (TU Munich) Volker Blum (Duke Univ.) Patrick Rinke (Aalto Univ., Finland) Ville Havu (Aalto Univ., Finland)
Making Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
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 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 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 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 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 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 informationTUTORIAL 6: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION
Hands-On Tutorial Workshop, July 29 th 2014 TUTORIAL 6: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION Christian Carbogno & Manuel Schöttler Fritz-Haber-Institut der Max-Planck-Gesellschaft,
More informationLecture 16: DFT for Metallic Systems
The Nuts and Bolts of First-Principles Simulation Lecture 16: DFT for Metallic Systems Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview of talk What
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
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 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 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 informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
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 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 informationIntroduction to Geometry Optimization. Computational Chemistry lab 2009
Introduction to Geometry Optimization Computational Chemistry lab 9 Determination of the molecule configuration H H Diatomic molecule determine the interatomic distance H O H Triatomic molecule determine
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 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 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 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 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 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 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 informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
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 informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationAccuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation
Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Multi-wavelets: High accuracy
More informationNon-linear optics, k p perturbation theory, and the Sternheimer equation
Non-linear optics, k p perturbation theory, and the Sternheimer equation David A. Strubbe Department of Materials Science and Engineering Massachusetts Institute of Technology, Cambridge, MA Formerly Department
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
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 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 informationWavelets for density functional calculations: Four families and three. applications
Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating
More informationElectronic energy optimisation in ONETEP
Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density
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 informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
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 informationHigher-Order Methods
Higher-Order Methods Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. PCMI, July 2016 Stephen Wright (UW-Madison) Higher-Order Methods PCMI, July 2016 1 / 25 Smooth
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
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 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 informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More 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 informationELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers Victor Yu and the ELSI team Department of Mechanical Engineering & Materials Science Duke University Kohn-Sham Density-Functional
More informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
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 informationTutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra
Tutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra Iurii Timrov SISSA Scuola Internazionale Superiore di Studi Avanzati, Trieste Italy itimrov@sissa.it Computer modelling of materials
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationComparison of methods for finding saddle points without knowledge of the final states
JOURNAL OF CHEMICAL PHYSICS VOLUME 121, NUMBER 20 22 NOVEMBER 2004 Comparison of methods for finding saddle points without knowledge of the final states R. A. Olsen and G. J. Kroes Leiden Institute of
More informationTable of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures
Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state
More informationLarge-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
More 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 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 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 informationThe Linearized Augmented Planewave (LAPW) Method
The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable
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 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 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 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 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 informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one
More informationChapter 8 Gradient Methods
Chapter 8 Gradient Methods An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Introduction Recall that a level set of a function is the set of points satisfying for some constant. Thus, a point
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 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 informationTools for QM studies of large systems
Tools for QM studies of large systems Automated, hessian-free saddle point search & characterization QM/MM implementation for zeolites Shaama Mallikarjun Sharada Advisors: Prof. Alexis T Bell, Prof. Martin
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 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 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 informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More 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 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 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 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 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 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 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 informationQUANTUM MECHANICS. Franz Schwabl. Translated by Ronald Kates. ff Springer
Franz Schwabl QUANTUM MECHANICS Translated by Ronald Kates Second Revised Edition With 122Figures, 16Tables, Numerous Worked Examples, and 126 Problems ff Springer Contents 1. Historical and Experimental
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 informationParameterization of a reactive force field using a Monte Carlo algorithm
Parameterization of a reactive force field using a Monte Carlo algorithm Eldhose Iype (e.iype@tue.nl) November 19, 2015 Where innovation starts Thermochemical energy storage 2/1 MgSO 4.xH 2 O+Q MgSO 4
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 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 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 informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationCHEM6085: Density Functional Theory Lecture 10
CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationOrbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
More informationThe Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii
The Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii Fundamental Points About Which to Be Thinking It s clear the PES is useful, so how can I construct it for an arbitrary
More informationA dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 15 15 OCTOBER 1999 A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives Graeme Henkelman and Hannes
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 informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationTUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION
TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION 1 INVESTIGATED SYSTEM: Silicon, diamond structure Electronic and 0K properties see W. Huhn, Tutorial 2, Wednesday August 2 2 THE HARMONIC
More informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More information