Localized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy
|
|
- Clementine Stanley
- 5 years ago
- Views:
Transcription
1 Localized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 16 December 2009
2 Joint work with Michael Overton NYU Juan Meza LBNL Chao Yang LBNL February NC State, Raleigh NC
3
4 In the past Every attempt to employ mathematical methods in the study of chemical questions must be considered profoundly irrational and contrary to the spirit of chemistry. If mathematical analysis should ever hold a prominent place in chemistry - an aberration which is happily almost impossible - it would occasion a rapid and widespread degeneration of that science. Auguste Comte, 1830
5 In the future Whereas computational methods nowadays mainly supplement experimentally obtained information, they are expected increasingly to supersede this information. Stephen Goedecker, 1999
6
7 CdSe Solar Cell CMOS electrons > 1,000,000 electrons > 50,000 electrons
8
9 Many-body Schrödinger Equation H Ψ ( r 1, r 2,..., r ne ) = λψ ( r 1, r 2,..., r ne ) If the n u nuclei are fixed in space at positions ˆr j, n e H = - 1 Δ ri 2 1 n u j=1 n e i=1 z j r i ˆr j + 1 r 1 i,j n i r j e
10 Many-body Schrödinger Equation H Ψ ( r 1, r 2,..., r ne ) = λψ ( r 1, r 2,..., r ne ) If the n u nuclei are fixed in space at positions ˆr j, n e H = - 1 Δ ri 2 1 n u j=1 n e i=1 z j r i ˆr j + 1 r 1 i,j n i r j e Ψ contains all the information in the system and satisfies orthogonality and antisymmetry properties
11 Many-body Schrödinger Equation H Ψ ( r 1, r 2,..., r ne ) = λψ ( r 1, r 2,..., r ne ) If the n u nuclei are fixed in space at positions ˆr j, n e H = - 1 Δ ri 2 1 n u j=1 n e i=1 z j r i ˆr j + 1 r 1 i,j n i r j e Ψ contains all the information in the system and satisfies orthogonality and antisymmetry properties Ψ ( r 1, r 2,..., r ne ) 2 dr 1... dr ne gives probability density
12 Many-body Schrödinger Equation H Ψ ( r 1, r 2,..., r ne ) = λψ ( r 1, r 2,..., r ne ) If the n u nuclei are fixed in space at positions ˆr j, n e H = - 1 Δ ri 2 1 n u j=1 n e i=1 z j r i ˆr j + 1 r 1 i,j n i r j e Ψ contains all the information in the system and satisfies orthogonality and antisymmetry properties Ψ ( r 1, r 2,..., r ne ) 2 dr 1... dr ne gives probability density λ represents the quantized energy
13 Many-body Schrödinger Equation H Ψ ( r 1, r 2,..., r ne ) = λψ ( r 1, r 2,..., r ne ) If the n u nuclei are fixed in space at positions ˆr j, n e H = - 1 Δ ri 2 1 n u j=1 n e i=1 z j r i ˆr j + 1 r 1 i,j n i r j e Ψ contains all the information in the system and satisfies orthogonality and antisymmetry properties Ψ ( r 1, r 2,..., r ne ) 2 dr 1... dr ne gives probability density λ represents the quantized energy If we attempt to solve on a 32 x 32 x 32 grid, and n e = 5, H has dimension
14
15 Density Functional Theory Hohenberg and Kohn 1964: Proved that at ground state, total energy is function of charge density ρ ( r) = n e Ψ* ( r, r 2,..., r ne ) Ψ ( r, r 2,..., r ne ) dr 2 dr 3... dr ne
16 Density Functional Theory Hohenberg and Kohn 1964: Proved that at ground state, total energy is function of charge density ρ ( r) = n e Ψ* ( r, r 2,..., r ne ) Ψ ( r, r 2,..., r ne ) dr 2 dr 3... dr ne Kohn and Sham 1965: Proposed a practical formulation, using n e single-particle, orthogonal wavefunctions ψ j that do not interact. Electron-electron interactions modeled by exchange-correlation energy.
17 The Kohn-Sham Equations Charge Density: n e ρ (r) = ψ* i (r ) ψ i ( r ) 1
18 The Kohn-Sham Equations Charge Density: n e ρ (r) = ψ* i (r ) ψ i ( r ) Exchange-correlation energy: E xc (ρ) = ρ (r ) ε xc [ρ ( r )]dr 1
19 The Kohn-Sham Equations Charge Density: n e ρ (r) = ψ* i (r ) ψ i ( r ) Exchange-correlation energy: E xc (ρ) = ρ (r ) ε xc [ρ ( r )]dr 1 Kohn-Sham energy: n e E KS [{ψ j }] = ψ j ( Δψ j ) dr + V ion ( r ) ρ ( r ) dr ρ ( r ) ρ ( r ) r r drdr + E xc ( ρ ).
20 Kohn-Sham Equations continued min E KS ( ψ j ) s.t. ψi * ψ j = δ i,j for i = 1, 2,..., n e
21 Kohn-Sham Equations continued KKT Conditions: min E KS ( ψ j ) s.t. ψi * ψ j = δ i,j for i = 1, 2,..., n e H KS ( ρ ) ψ i = λ i ψ i, ψ j * ψ j = δ i,j, for i = 1, 2,..., n e where H KS = Δ + V ion ( r ) +ρ (r ) * 1 r + V xc ( ρ) = Δ + V ( ρ (r) )
22 Non-orthogonal Kohn-Sham Equations
23 Non-orthogonal Kohn-Sham Equations Overlap matrix: S jk = ψ* j ψ k dr Charge-density: ρ (r) = j ψ* j (r ) (S -1 ) jk ψ k ( r) k
24 Non-orthogonal Kohn-Sham Equations Overlap matrix: S jk = ψ* j ψ k dr Charge-density: ρ (r) = j ψ* j (r ) (S -1 ) jk ψ k ( r) k Kohn-Sham Energy: E KS ( {ψ j }) = ( S -1 ) jk j k ψ j ( Δψ k ) dr + V ion ( r ) ρ ( r ) dr ρ ( r ) ρ ( r ) r r drdr + E xc ( ρ)
25 Non-orthogonal Kohn-Sham Equations Overlap matrix: S jk = ψ* j ψ k dr Charge-density: ρ (r) = j ψ* j (r ) (S -1 ) jk ψ k ( r) k Kohn-Sham Energy: E KS ( {ψ j }) = ( S -1 ) jk j k ψ j ( Δψ k ) dr + V ion ( r ) ρ ( r ) dr ρ ( r ) ρ ( r ) r r drdr + E xc ( ρ) This energy is independent of basis
26
27 Numerical Methods Step 1: Choose a discretization
28 Numerical Methods Step 1: Choose a discretization Plane-waves basis: ψ ( r) = k α k e ig k r 1 Advantages: ρ ( r ) * r Matrix Multiplication (FFT) Disadvantages: Inherently nonlocal
29 Numerical Methods Step 1: Choose a discretization Plane-waves basis: ψ ( r) = k α k e ig k r 1 Advantages: ρ ( r ) * r Matrix Multiplication (FFT) Disadvantages: Inherently nonlocal Finite-differences: ψ ( r j ) = 1 h [ψ ( r j + h ) ψ ( r j h)] Advantages: Localized, highly sparse 1 Disadvantages:ρ ( r ) * r requires solution to a Poisson problem
30 Discretized Kohn-Sham Equations E KS (X) = E kineticionic ( X ) +E Hartree ( X ) +E xc ( X) ( ( E kineticionic = trace (X*X) -1 X* 1 ) ) 2 L + V ion X E Hartree = 1 2 ρt L ρ E xc = ρ T ε xc [ρ] [ ρ (X) = diag X ( X*X ) -1 ] X*
31 Discretized Kohn-Sham Equations E KS (X) = E kineticionic ( X ) +E Hartree ( X ) +E xc ( X) ( ( E kineticionic = trace (X*X) -1 X* 1 ) ) 2 L + V ion X E Hartree = 1 2 ρt L ρ E xc = ρ T ε xc [ρ] [ ρ (X) = diag X ( X*X ) -1 ] X* E (XG) = E (X)
32 Simpli ed Kohn-Sham Equations (The Model Problem) ( ( E (X) = trace (X*X ) -1 1 X* 2 L + α ) ) 1 2 V X + α 2 4 ρt L -1 ρ [ ρ = diag X ( X*X ) -1 ] X*
33 Simpli ed Kohn-Sham Equations (The Model Problem) ( ( E (X) = trace (X*X ) -1 1 X* 2 L + α ) ) 1 2 V X + α 2 4 ρt L -1 ρ Look Familiar? [ ρ = diag X ( X*X ) -1 ] X*
34 Simpli ed Kohn-Sham Equations (The Model Problem) ( ( E (X) = trace (X*X ) -1 1 X* 2 L + α ) ) 1 2 V X + α 2 4 ρt L -1 ρ [ ρ = diag X ( X*X ) -1 ] X* Look Familiar? x*ax ( x*x trace (X*X ) -1 ) X*AX
35 Numerical Methods Step 1: Choose a discretization Step 2: Compute a solution
36 Numerical Methods Step 1: Choose a discretization Step 2: Compute a solution SCF: Solve the nonlinear Kohn-Sham equations
37 Numerical Methods Step 1: Choose a discretization Step 2: Compute a solution SCF: Solve the nonlinear Kohn-Sham equations CG: Directly minimize the Kohn-Sham Energy
38 Self-Consistent Field Iteration [ ½Δ + V(ρ(r))]ψ = Eψ i i i ρ(r) = ψ ² i V(ρ(r))
39 Self-Consistent Field Iteration [ ½Δ + V(ρ(r))]ψ = Eψ i i i O ( n 3 e ) ρ(r) = ψ ² i V(ρ(r))
40 Self-Consistent Field Iteration [ ½Δ + V(ρ(r))]ψ = Eψ ρ(r) = ψ ² i i i i O ( n 3 e ) Convergence slow, non-monotonic V(ρ(r))
41 Self-Consistent Field Iteration [ ½Δ + V(ρ(r))]ψ = Eψ ρ(r) = ψ ² i i i i O ( n 3 e ) Convergence slow, non-monotonic Charge-mixing to accelerate convergence V(ρ(r))
42 CG (Fletcher-Reeves) 1. RM = - X E ( X 0 ) 2. PM = RM 3. XM = X 0 4. repeat until convergence criteria are satisfied 5. Choose αm from a line search on E ( XM + αpm) 6. X k+1 = XM + αmpm 7. R k+1 = - X E ( X k+1 ) 8. β k+1 = (rk+1 ) T r k+1 (r k ) T r k 9. P k+1 = R k+1 + β k+1 ( PM) 10. k = k end loop
43
44 Two Equal Solutions Orthogonal Solution A Localized Solution Student Version of MATLAB Student Version of MATLAB
45 Simpli ed Kohn-Sham Equations (The Model Problem) ( ( E (X) = trace (X*X ) -1 1 X* 2 L + α ) ) 1 2 V X + α 2 4 ρt L -1 ρ [ ρ = diag X ( X*X ) -1 ] X*
46 Localization and Sparsity nz = nz = 150
47
48 Maintaining Localization/sparsity An iterative algorithm need not maintain the sparsity pattern.
49 Maintaining Localization/sparsity An iterative algorithm need not maintain the sparsity pattern. Naive approach: Truncate, i.e. set X k+1 = (X k+1 ) T
50 Maintaining Localization/sparsity An iterative algorithm need not maintain the sparsity pattern. Naive approach: Truncate, i.e. set X k+1 = (X k+1 ) T Before truncation After truncation Student Version of MATLAB Student Version of MATLAB
51 Maintaining Localization/sparsity An iterative algorithm need not maintain the sparsity pattern. Naive approach: Truncate, i.e. set X k+1 = (X k+1 ) T E and Gao approach: Minimize truncation error (localize)
52 Maintaining Localization/sparsity An iterative algorithm need not maintain the sparsity pattern. Naive approach: Truncate, i.e. set X k+1 = (X k+1 ) T E and Gao approach: Minimize truncation error (localize) Compute invertible G which minimizes XG (XG) T F Remember E (XG) = E (X) for any invertible G
53 Sparse Subspace E((XMG T )) E(XMG) E(XMG) = E(XM) (XMG) T XM
54 CG + Localization (CG L by E and Gao) 1. RM = - X E ( X 0 ) 2. PM = RM 3. XM = X 0 4. repeat until convergence criteria are satisfied 5. Choose αm from a line search on E ( XM + αpm) 6. X k+1 = XM + αmpm 7. G = arg min G X k+1 G (X k+1 G ) T 8. X k+1 = (X k+1 G) T 9. R k+1 = - X E ( X k+1 ) 10. β k+1 = (rk+1 ) T r k+1 (r k ) T r k 11. P k+1 = R k+1 + β k+1 ( PMG) T 12. k = k end loop
55 Avoiding Local Minima: CG L versus truncation n = 100, n e = localized truncated
56 CG + Localization (CG L by E and Gao) 1. RM = - X E ( X 0 ) 2. PM = RM 3. XM = X 0 4. repeat until convergence criteria are satisfied 5. Choose αm from a line search on E ( XM + αpm) 6. X k+1 = XM + αmpm 7. G = arg min G X k+1 G (X k+1 G ) T 8. X k+1 = (X k+1 G) T 9. R k+1 = - X E ( X k+1 ) 10. β k+1 = (rk+1 ) T r k+1 (r k ) T r k 11. P k+1 = R k+1 + β k+1 ( PMG) T 12. k = k end loop
57 The Plan Develop an optimization algorithm that
58 The Plan Develop an optimization algorithm that maintains the property of avoiding local minima
59 The Plan Develop an optimization algorithm that maintains the property of avoiding local minima evaluates the energy only at sparse iterates.
60 LocOpt: Localize the search direction { Sparse Dense { Pk { Sparse + Pk Localization (( + Pk)Gk) T
61 LocOptCG 1. RM = - X E ( X 0 ) 2. PM = RM 3. XM = X 0 4. repeat until convergence criteria are satisfied 5. G = arg min G ( XM + PM ) G (( XM + PM ) G ) T 6. PM = (( XM + PM ) G ) T XM 7. Choose αm from a line search on E ( XM + α PM) 8. X k+1 = XM + αm PM 9. R k+1 = - X E ( X k+1 ) 10. β k+1 = (rk+1 ) T r k+1 (r k ) T r k 11. P k+1 = R k+1 + β k+1 PM 12. k = k end loop
62 Avoiding Local Minima: LocOptCG vs: CG L!0.625! E and Gao o - Our method!0.635!0.64!0.645!0.65! Tolerance = 10 ⁷; n e = 5, n = 100
63 Monotonic Convergence !1 - E and Gao - Our method log(e - E*) 10!2 10!3 10!4 10!5 10!6 10!7 10!8 10! n e = 5, n = 100
64 A real problem: Methane (CH₄)
65 A real problem: Methane (CH₄) 8 Electrons in 4 Orbitals
66 A real problem: Methane (CH₄) 8 Electrons in 4 Orbitals Hamiltonian dimension 5832x5832
67 A real problem: Methane (CH₄) 8 Electrons in 4 Orbitals Hamiltonian dimension 5832x5832 Known solution from : Ry
68 Methane (CH₄) Solution
69 Methane (CH₄) Solution nz = 2762
70 Methane (CH₄) Solution CG LMBFGS nz = 2762 Energy: Ry
71 Methane (CH₄) Solution (cont) LocOpt: Rho -RSDFT: Rho
72 A Variation for { Sparse Dense { + PM PM + PM)GM { Sparse Localization (( ) T (( + PM)GM) T E T
73 Something magic happens 10 2 log(em (-17.62))) LMBFGS + old CG + new LMBFGS + new
74
75 Scale & Conditioning E (XG) = E (X)
76 Scale & Conditioning E (XG) = E (X) In particular, what if G = Diag[α 1, α 2,..., α ne ]?
77 Scale & Conditioning E (XG) = E (X) In particular, what if G = Diag[α 1, α 2,..., α ne ]? Can we normalize columns?
78 Scale & Conditioning E (XG) = E (X) In particular, what if G = Diag[α 1, α 2,..., α ne ]? Can we normalize columns? How does this affect search direction scale?
79 Localization & Descent min s.t. ( XM + PM ) G (( XM + PM ) G ) T F G is invertible
80 Localization & Descent min s.t. ( XM + PM ) G (( XM + PM ) G ) T F G is invertible ( trace E T P ) < 0
81 Localization & Descent min s.t. ( XM + PM ) G (( XM + PM ) G ) T F G is invertible ( trace E T P ) ( trace ( E ) T ( )) T ( XM + PM ) G < 0.
82 Localization & Descent min s.t. ( XM + PM ) G (( XM + PM ) G ) T F G is invertible ( trace ( E ) T ( )) T ( XM + PM ) G < 0
83 The Gradient E KS (X) = E kineticionic ( X ) +E Hartree ( X ) +E xc ( X) E kineticionic -XS -1 X* ( L + V ) XS -1 + (L + V ) XS -1 E Hartree Diag ( L -1 ρ ) XS -1 XS -1 X* Diag ( L -1 ρ ) XS -1 E xc... ρ (X) = diag [ X ( X*X ) -1 ] X*
84 The Overlap Matrix x S S -1 Student Version of MATLAB Student Version of MATLAB
85
86 Integration with BFGS (LMBFGS) 10 2 CG !2 BFGS 10!4 10!6 10!8 10!10 10!12 10!14 10! n e = 5, n = 100. Bad Search Directions: CG (0), BFGS (7)
87 Solve (more) realistic problems Based on by Chelikowsky, Troullier, and Saad
88 Solve (more) realistic problems Based on by Chelikowsky, Troullier, and Saad Open-source Matlab code
89 Solve (more) realistic problems Based on by Chelikowsky, Troullier, and Saad Open-source Matlab code In real-space, using finite-differences
90 Solve (more) realistic problems Based on by Chelikowsky, Troullier, and Saad Open-source Matlab code In real-space, using finite-differences Limited in problem size and mesh length
91 Solve (more) realistic problems Based on by Chelikowsky, Troullier, and Saad Open-source Matlab code In real-space, using finite-differences Limited in problem size and mesh length Just started designing and writing a object-oriented C/C++ implementation to run on machines
92
Minimization of the Kohn-Sham Energy with a Localized, Projected Search Direction
Minimization of the Kohn-Sham Energy with a Localized, Projected Search Direction Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 29 October 2009 Joint work with Michael
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 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 informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part II Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationA new future for finite-element methods in Quantum Chemistry?
A new future for finite-element methods in Quantum Chemistry? Kenneth Ruud Department of Chemistry University of Tromsø 9037 Tromsø Norway Geilo, 29/1 2007 A new future for finite-element methods in Quantum
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 informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
More 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 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 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 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 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 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 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 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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More informationIntroduction to Computational Chemistry Computational (chemistry education) and/or. (Computational chemistry) education
Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools to help increase student understanding of material
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
More informationSolutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid
Solutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid These problems have been solved as a part of the Independent study class with prof. Neepa Maitra. Compiled
More information7.3 The Jacobi and Gauss-Siedel Iterative Techniques. Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP.
7.3 The Jacobi and Gauss-Siedel Iterative Techniques Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP. 7.3 The Jacobi and Gauss-Siedel Iterative Techniques
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 informationVASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria The Many-Body Schrödinger equation 0 @ 1 2 X i i + X i Ĥ (r 1,...,r
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
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 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 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 informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More information2.5 Time dependent density functional theory
.5 Time dependent density functional theory The main theorems of Kohn-Sham DFT state that: 1. Every observable is a functional of the density (Hohenger-Kohn theorem).. The density can be obtained within
More informationAdaptively Compressed Polarizability Operator
Adaptively Compressed Polarizability Operator For Accelerating Large Scale ab initio Phonon Calculations Ze Xu 1 Lin Lin 2 Lexing Ying 3 1,2 UC Berkeley 3 ICME, Stanford University BASCD, December 3rd
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 informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
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 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 informationElectronic Structure Calculations, Density Functional Theory and its Modern Implementations
Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011 Outline 1 of Atomistic calculations
More 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 informationCOMPUTATIONAL CHEMISTRY. Introduction to the Theory and Applications of Molecular and Quantum Mechanics
COMPUTATIONAL CHEMISTRY Introduction to the Theory and Applications of Molecular and Quantum Mechanics COMPUTATIONAL CHEMISTRY Introduction to the Theory and Applications of Molecular and Quantum Mechanics
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 informationApplications of Matrix Functions Part III: Quantum Chemistry
Applications of Matrix Functions Part III: Quantum Chemistry Emory University Department of Mathematics and Computer Science Atlanta, GA 30322, USA Prologue The main purpose of this lecture is to present
More informationCompressed representation of Kohn-Sham orbitals via selected columns of the density matrix
Lin Lin Compressed Kohn-Sham Orbitals 1 Compressed representation of Kohn-Sham orbitals via selected columns of the density matrix Lin Lin Department of Mathematics, UC Berkeley; Computational Research
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More 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 informationIntroduction to Parallelism in CASTEP
to ism in CASTEP Stewart Clark Band University of Durham 21 September 2012 Solve for all the bands/electrons (Band-) Band CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei:
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More 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 informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationPhani Motamarri and Vikram Gavini Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109, USA
A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization Phani Motamarri and Vikram Gavini Department
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 informationReinhold Schneider, MATHEON TU Berlin. Direct Minimization for effective single particle models (DFT)
Direct Minimization for effective single particle models (DFT) Partners joint work with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
More informationSolving the Generalized Poisson Equation Using the Finite-Difference Method (FDM)
Solving the Generalized Poisson Equation Using the Finite-Difference Method (FDM) James R. Nagel September 30, 2009 1 Introduction Numerical simulation is an extremely valuable tool for those who wish
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 informationThe Conjugate Gradient Method
The Conjugate Gradient Method Jason E. Hicken Aerospace Design Lab Department of Aeronautics & Astronautics Stanford University 14 July 2011 Lecture Objectives describe when CG can be used to solve Ax
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 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 informationarxiv: v1 [math.na] 18 Mar 2012
DECAY PROPERTIES OF SPECTRAL PROJECTORS WITH APPLICATIONS TO ELECTRONIC STRUCTURE MICHELE BENZI, PAOLA BOITO, AND NADER RAZOUK arxiv:1203.3953v1 [math.na] 18 Mar 2012 Abstract. Motivated by applications
More informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute for Microstructure Physics OUTLINE LECTURE I Phenomena to be described by TDDFT Some generalities on functional theories LECTURE
More informationIntroduction to Scientific Computing
(Lecture 5: Linear system of equations / Matrix Splitting) Bojana Rosić, Thilo Moshagen Institute of Scientific Computing Motivation Let us resolve the problem scheme by using Kirchhoff s laws: the algebraic
More informationSteepest Descent. Juan C. Meza 1. Lawrence Berkeley National Laboratory Berkeley, California 94720
Steepest Descent Juan C. Meza Lawrence Berkeley National Laboratory Berkeley, California 94720 Abstract The steepest descent method has a rich history and is one of the simplest and best known methods
More informationCHEM6085: Density Functional Theory
Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron
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 informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-3: Unconstrained optimization II Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428
More informationElectron-Proton Correlation, Theory, and Tunneling Splittings. Sharon Hammes-Schiffer Penn State University
Nuclear-Electronic Orbital Approach: Electron-Proton Correlation, Multicomponent Density Functional Theory, and Tunneling Splittings Sharon Hammes-Schiffer Penn State University Nuclear Quantum Effects
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 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 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 informationApplied Machine Learning for Biomedical Engineering. Enrico Grisan
Applied Machine Learning for Biomedical Engineering Enrico Grisan enrico.grisan@dei.unipd.it Data representation To find a representation that approximates elements of a signal class with a linear combination
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 informationSpring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor
2145-25 Spring College on Computational Nanoscience 17-28 May 2010 Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor Stefano BARONI SISSA & CNR-IOM DEMOCRITOS Simulation Center
More informationChapter 4. Unconstrained optimization
Chapter 4. Unconstrained optimization Version: 28-10-2012 Material: (for details see) Chapter 11 in [FKS] (pp.251-276) A reference e.g. L.11.2 refers to the corresponding Lemma in the book [FKS] PDF-file
More informationChemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer. Lecture 25, November 5, 2003
Chemistry 3502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2003 Christopher J. Cramer Lecture 25, November 5, 2003 (Some material in this lecture has been adapted from Cramer, C.
More informationConsider the following example of a linear system:
LINEAR SYSTEMS Consider the following example of a linear system: Its unique solution is x + 2x 2 + 3x 3 = 5 x + x 3 = 3 3x + x 2 + 3x 3 = 3 x =, x 2 = 0, x 3 = 2 In general we want to solve n equations
More informationNotes on PCG for Sparse Linear Systems
Notes on PCG for Sparse Linear Systems Luca Bergamaschi Department of Civil Environmental and Architectural Engineering University of Padova e-mail luca.bergamaschi@unipd.it webpage www.dmsa.unipd.it/
More informationMethods that avoid calculating the Hessian. Nonlinear Optimization; Steepest Descent, Quasi-Newton. Steepest Descent
Nonlinear Optimization Steepest Descent and Niclas Börlin Department of Computing Science Umeå University niclas.borlin@cs.umu.se A disadvantage with the Newton method is that the Hessian has to be derived
More informationNotes on Density Functional Theory
Notes on Density Functional Theory Rocco Martinazzo E-mail: rocco.martinazzo@unimi.it Contents 1 Introduction 1 Density Functional Theory 7 3 The Kohn-Sham approach 11 1 Introduction We consider here a
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationContinuous Optimisation, Chpt 6: Solution methods for Constrained Optimisation
Continuous Optimisation, Chpt 6: Solution methods for Constrained Optimisation Peter J.C. Dickinson DMMP, University of Twente p.j.c.dickinson@utwente.nl http://dickinson.website/teaching/2017co.html version:
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 informationDensity-functional theory in quantum chemistry. Trygve Helgaker. From Quarks to the Nuclear Many-Body Problem
1 Density-functional theory in quantum chemistry Trygve Helgaker Centre for Theoretical and Computational Chemistry, University of Oslo, Norway From Quarks to the Nuclear Many-Body Problem A conference
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 informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
More 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 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 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 informationAdaptive Application of Green s Functions:
Adaptive Application of Green s Functions: Fast algorithms for integral transforms, with a bit of QM Fernando Pérez With: Gregory Beylkin (Colorado) and Martin Mohlenkamp (Ohio University) Department of
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 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 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 informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More 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 informationDecay Properties of Spectral Projectors with Applications to Electronic Structure
SIAM REVIEW Vol. 55, No. 1, pp. 3 64 c 2013 Society for Industrial and Applied Mathematics Decay Properties of Spectral Projectors with Applications to Electronic Structure Michele Benzi Paola Boito Nader
More informationDirect Minimization in Density Functional Theory
Direct Minimization in Density Functional Theory FOCM Hongkong 2008 Partners joint paper with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
More informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
More 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 informationPreconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
More informationNumerical Methods for Electronic Structure Calculations of Materials
SIAM REVIEW Vol. 52, No. 1, pp. 3 54 c 2010 Society for Industrial and Applied Mathematics Numerical Methods for Electronic Structure Calculations of Materials Yousef Saad James R. Chelikowsky Suzanne
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 information