Reinhold Schneider, MATHEON TU Berlin. Direct Minimization for effective single particle models (DFT)
|
|
- Elmer Fox
- 5 years ago
- Views:
Transcription
1 Direct Minimization for effective single particle models (DFT)
2 Partners joint work with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry Deutsch, Co-ordinator, Atomistic Simulation group at CEA-Grenoble, Prof. Stefan Goedecker, Computational Physics Group at the Basel University, Switzerland; Prof. Xavier Gonze, Uni. de Phy.-Ch.et de Phys. des Mat., Univ. Catholique de Louvain; Joint project with H.J. Flad (Matheon - TUB) and W. Hackbusch (MPI Leipzig) DFG Priority Program: Modern and Universal First-Principle-Methods for Many-Electron Systems in Chemistry and Physics
3 Basic model - electronic Schrödinger equation Electronic Schrödinger equation N nonrelativistic electrons + Born Oppenheimer approximation HΨ = EΨ The Hamilton operator H = 1 2 i i i K ν=1 Z ν x i a ν x i x j i,j acts on anti-symmetric wave functions Ψ H 1 (R 3 {± 1 2 })N, Ψ(x 1, s 1,..., x N, s N ), x i = (x i, s i ) R 3 {± 1 2 }.
4 Hartree-Fock- (HF) Approximation Ground state energy E 0 = min{ Hψ, ψ : ψ, ψ = 1} approximation of ψ by a single Slater determinant Ψ SL Φ (x 1, s 1,..., x N, s N ) := 1 N! det(φ i(x j, s j )) Closed Shell Restricted HF (RHF): N := N 2 electron pairs minimization of the functional J HF (Φ) D E Φ J HF (Φ) := HΨ SL Φ, ΨSL Φ = + NX Z j=1 R 3 " φj (y) 2 φ i (x) 2 x y w.r.t. orthogonality constraints NX Z ` φi (x) 2 + 2V core(x) φ i (x) 2 + i=1 # 1 φ i (x)φ i (y) φ j (x)φ j (y) dy dx 2 x y Φ = (φ i ) N i=1 H1 (R 3 ) N and φ i, φ j = δi,j
5 Kohn-Sham model Theorem (Kohn-Hohenberg) The ground state energy E 0 is a functional of the electron density n. 1 2 E 0 E KS = inf{j KS (Φ) : φ i, φ j = δ ij } minimization of the Kohn Sham energy functional J KS (Φ) J KS (Φ) = ( Z 1 2 NX Z φ i 2 + nv core + 1 Z Z 2 i=1 n(x)n(y) x y dx dy E xc(n) φ i H 1 (R 3 ), electron density n(x) := N i=1 φ i(x) 2 E xc (n) exchange-correlation-energy (not known explicitly) e.g. LDA (local density approximation) n E xc (n) : R R )
6 General setting Φ := (φ 1,..., φ N ) (H 1 (R 3 )) N = V N Gelfand triple V := H 1 (R 3 ) L 2 (R 3 ) H 1 (R 3 ) = V Φ T Ψ := ( φ i, ψ j )i,j R N N scalar product Φ, Ψ := tr Φ T Ψ = N i=1 φ i, ψ i R AΦ := (Aφ 1,..., Aφ N ), A : V V Simplified Problem: minimize J SCF A (Φ) := N i=1 Aφ i, φ i = tr Φ T AΦ = Φ, AΦ w.r.t. to orthogonality constraints Φ T Φ = I.
7 Remarks J = J HF, J KS, J SCF are invariant under unitary transformations U U N J (Φ) = J (ΦU), ΦU := ( φ i u ij ) j=1,...,n i=1 gradient J (Φ) = A [Φ] Φ, A [Φ] : V V The density matrix operator D Φ := N φ i, L 2 φ i i=1 projects onto span Φ := span{φ i : 1 i N}, D 2 = D, tr D = N, D T = D
8 Stiefel and Grassmann manifolds Definition Stiefel manifold S V,N := S := {Φ = (φ i ) N i=1 φ i V, Φ T Φ I N N = 0 R N N } Grassmann manifold is a quotient manifold G V,N := G := S V,N /, Φ Φ Φ = ΦU, U U(N) (identify ONB spanning the same subspace span Φ) see [Edelman, Arias, Smith] for V = R n. There is a one-to-one correspondence between [Φ] G D Φ (density matrix operator)
9 Tangent space Proposition (Edelman, Arias, Smith (98); Blauert, Neelov, Rohwedder, S. (08)) tangent space T [Φ] G = {δψ V N (δψ) T Φ = 0 R N N } (I D Φ ) : V N T [Φ] G, is an orthogonal projection onto the tangent space T [Φ] G tangent space T Φ S = T [Φ] G + {ΦA : A T = A} = {Θ V N : Θ T Φ = Φ T Θ }
10 Existence and uniqueness known results Problem Minimize J : G R, J = J HF, J KS, J SCF A Theorem (Lieb-Simon (78), P.L. Lions (89)) There exists a minimizer Φ (H 1 (R 3 )) N for J HF provided that N Z ν. Uniqueness is not know yet? J SCF A has unique minimizer on Grassmann manifold G if there is a gap for the lowest eigenvalues of A: λ N < λ N+1 no uniqueness on Stiefel manifold S.
11 1st order optimality conditions For [Φ] G there hold J (Φ) = A [Φ] Φ (V ) N where A SCF [Φ] and A HF [Φ], AKS [Φ] are the Hamilton Fock operators = A Definition E.g. the Kohn-Sham Hamilton Fock operator is given by ( A KS [Φ] := V core + n 1 ) + v xc (n) = V (n) Nec. cond.: If [Ψ] = argmin {J (Φ) : [Φ] G} V N (V N h ) then A [Ψ] Ψ, δφ = 0 δφ T [Ψ] G V N (V N h ) (I D Ψ )A [Ψ] Ψ, δφ = 0 δφ V N (V N h )
12 Canonical HF - KS equations Lagrangian L(Φ, Λ) := J (Φ) trλ( Φ T Φ I) At a (local) minimum Ψ V, i.e. at stationary points (Ψ, λ) of L there holds Λ = A [Ψ] ˆψ i, ˆψ j. There exists a representation Ψ = ˆΨU of [ ˆΨ] G, s. t. λ i,j = λ i δ ij and λ 1 λ 2... s.t. There hold the canonical (HF or KS ) equations A [Ψ] ψ i = λ i ψ i, λ 1 λ 2... λ N < 0. Self consistent iteration: A [Ψ (n) ] ψ(n+1) i (Cances-LeBris) = λ (n+1) i ψ (n+1) i
13 Self consistent iteration - SCF We define A D n such that A [Ψ] = A DΨ for an appropriate D n P := {D : D = D, trd = N, D 2 D}. Algorithm (Self consistent iteration) 1 Ψn+1 := argmin {JA SCF (Φ) : Φ G} = D n argmin { A D nφ, Φ : Φ G} 2 D n+1 := µd Ψ n + (1 µ)d Ψn+1 P for some µ [0, 1] Roothaan algorithm: A [Ψ ( n)] ψ(n+1) i = λ (n+1) i ψ (n+1) i fails often Optimal damping algorithm (Cances-LeBris 02) converges D n+1 := µd Ψ n + (1 µ)d Ψn+1 P with optimal choice of µ [0, 1]
14 Gradient flow The minimizer Ψ is spanning an invariant subspace V Ψ := span{ψ 1,..., ψ N } of A [Ψ] corresponding to the N lowest eigenvalues. A local minimizer Ψ is a fix point of the gradient flow [Φ](t), Ψ = lim t Φ(t) or d dt Φ A [Φ]Φ(t), θ = 0 θ T [Φ] G V N (V N h ) d dt Φ(t) (I D Φ)A [Φ] Φ, θ = 0 θ V N (V N h )
15 Direct Minimization Algorithm (direct minimization) 1 A := J (Φ (n) ) = A [Φ (n) ] : V N V N 2 Λ (n) := Φ T AΦ = ( Aφ i, φ j )i,j=1,...,n R N N 3 ˆΦ (n+1) := Φ (n) αb 1 (AΦ (n) Φ (n) Λ (n) ), preconditioner B 4 ˆΦ (n+1) Φ (n+1), Φ (n+1) := P ˆΦ (n+1) by projection onto G, e.g. orthogonalization (AΦ (n) Φ (n) Λ (n) ) = [A, D Φ (n)]φ (n) = (I D Φ (n))aφ (n) B, ˆB : V V, B : V N V N, ˆB := (I D)B(I D) + DBD, D = D Ψ
16 Projection onto the Stiefel manifold Φ = Φ (n+1) = ( ˆφ i ) N i=1, and P Φ =: Φ V resp. G s.t. span { ˆφ i : i = 1,..., N} = span {φ i : i = 1,..., N} Projection P : V N S,resp. G Gram Schmidt orthogonalization Löwdin transformation Φ = L 1 Φ where LL T = Φ T Φ Diagonalization of Λ (n+1) = Φ T A [Φ] Φ, yields the first N eigenvalues λ (n) 1... λ (n) N of A [Φ].
17 Comment on direct minimization improvement by subspace acceleration: e.g. line search and DIIS Φ (n+1) = Function(ˆΦ (n+1), Φ (n),..., Φ (0) ) gradient directed convergence with Armijo line search B : V V, φ 2 B = φ, φ B := Bφ, φ φ 2 H 1 e.g.: B C, e.g. multigrid or convolution by FFT also valid if V := V h is a finite dimensional subspace for the simplified problem cf. e.g. Knyazev et al.
18 Convergence results Definition measure of the error between subspaces spanned by Φ, Ψ G (I D Ψ )Φ 2ˆB := N φ i D Ψ φ i 2ˆB i=1 Theorem ( Blauert, Neelov, Rohwedder, S. (08)) If Φ (0) U δ (Ψ), and (J (Ψ) Λ)Φ, Φ γ Φ 2 H 1N Φ T [Ψ] G, then there exists χ < 1 such that for all (I D Ψ )Φ (n+1) ˆB χ (I D Ψ )Φ (n) ˆB
19 Energy error Constrained optimization problem (revisited): u = argmin{j(v) : G(v) = 0} Lagrangian L(x) := L(u, Λ) = J(u) ΛG(u) (x = (u, Λ) X) Theorem (Rannacher et al. ) If L (x)y = 0 y X and L (x n )x n = 0, L(x) L(x n ) = 1 2 L (x n )(x x n ) + O( x x n 3 X ) J(u) J(u n ) = 1 J 2[ (u n )(u u n ) Λ n G (u n )(u u n ) (Λ Λ n )G(u n ) ] Here L(Φ, Λ) = J (Φ) + trλ( Φ T, Φ I)
20 Energy error Theorem ( Blauert, Neelov, Rohwedder, S. (08)) Let J C 2 (V N, R), then J (Φ (n) ) J (Ψ) + R 2 = 2 (Ψ Φ (n) ), (A [Φ (n) ] Φ(n) Φ (n) Λ (n) ) (I D Ψ )Φ (n) 2 V N where R 2 = O( (I D Ψ )Φ (n) 2 ). V N Corresponding results for a priori and a posterriori estimates for the Galerkin solution together with S. Schwinger and W.Hackbusch (08).
21 Pseudo-potentials pseudo-potentials (nuclei + core electrons) K V core (x) = Z j f Rj (x) + P j f Rj (x) = j=1 1 x R j if x R j > r j consider only N - valence electrons (pairs), 1 Φ(x) 1 P VL Φ(x) 1 (I P VL )Φ(x) (I P VL )Φ(x) x x x + non-local pseudo-potentials solutions are smooth orbitals! This allows e.g. relativistic corrections! But modeling error!!!
22 Discretization Let ψ = (ψ λ ) λ I be basis in V h V, dimv h = N <. Φ (n) = (ϕ n 1,..., ϕn N ) be given ϕ(n) i = λ I c(n) i,λ ψ λ = Ψ T c (n) i We build the Galerkin matrices V (n) H T = ( 1 2 ψ λ, ψ λ ), V c = ( V core (x)ψ λ, ψ λ ), (n) = ( V (x)ψ λ, ψ λ ), W (n) = ( ( W (n) + V (n) xc )ψ λ, ψ λ ) H with the Hartree potential V H and the exchange energy terms Hamilton Fock matrix F (n) = F Φ (n) = T + V c + V (n) H + W(n)
23 Basis functions Local Atomic Orbitals (LAO), Gaussian type orbitals (GTO): ψ λ (x) = e a λ R j x 2 R j x 2n Y l,m (x) adapted to atoms; 3,..., 20 basis functions per electron plane waves, trigonometric polynomials, ψ λ = e ik x C for crystals, periodic setting, requires (non-local) pseudo potentials (only valence electrons) 100,... 10, 000 basis functions per atom. local basis functions, finite elements, splines,wavelets, grid methods etc., without pseudo potentials , with pseudo potentials basis functions per atom
24 BigDFT EU NEST Project for developing a linear scaling wavelet code for density functional (DFT) Sim/BigDFT/ Code is incorporated in ABINIT: ( an Open Source plane wave DFT code (molecules and crystals),
25 Figure: BigDFT orth. wavelets representing φ i, interpolets for the density n, convergence rate for E 0 about O(h 14 )! (non-local) pseudo-potentials, semi-adaptivity (2 level of resolution),
26 BigDFT approximate Galerkin approximate V local G 1 G 2 G 3 diag(v(x j k )) GT 1 GT 2 GT 3 (magic filter G = (g[i j])) by Kronecker tensor products. Fast evaluation of nonlocal pseudo potentials Preconditioner C, C ( µi) 1 based on Multilevel preconditioning. Hartree potential V H (x) = R 3 n(y) x y dy for updating the Hamilton Fock operator is a global procedure: Fast Poisson Solver Newton potential based on Kronecker tensor products
27 BigDFT direct minimization, linear scaling (in progress),
28 Numerical results Figure: example: C 19 H 22 N 2 O (N = 55)
29 Numerical results Convergence history for the direct minimization scheme and together with DIIS acceleration 10 1 steepest descent 10 1 history size for DIIS = hgrid = 0.3 hgrid = 0.45 hgrid = hgrid = 0.3 hgrid = 0.45 hgrid = 0.7 absolute error absolute error iter, iter Figure: convergence
30 Numerical results Required resources 5000 memory 8000 runtime DIIS, hist. size = 6 steepest descent 7000 DIIS, hist. size = 6 steepest descent MByte sec hgrid, hgrid Figure: memory and computing time
31 Numerical results first N eigenvalues 0.8 least eigenvalues, hgrid = 0.3 least eigenvalues initial guess HOMO LUMO gap for initial guess final eigenvalues k, initial guess, hgrid = 0.3 HOMO LUMO gap for initial guess,hgrid = 0.3 final eigenvalues, hgrid = 0.3 initial guess, hgrid = 0.7 final eigenvalues, hgrid = 0.7 HOMO LUMO gap for initial guess,hgrid = k Figure: memory and computing time
32 Further comments: Linear scaling Representative of [Φ] G with minimal support Usual complexity scaling: O(N 2 dimv h ) = O(N 3 )! Often there exist orbitals Ψ = (ψ i ) N i=1 with [Ψ] = [Φ], ψ i (x) = e α x R i domain decomposition R 3 = M j=1 Ω j Ω, φ i 0 on Ω. linear scaling w. r. t. the number of particles (for fixed accuracy per atom) BigDFT and Bencteux, Hager, Cances & Lebris... (research in progress)
Direct 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 informationI ) Direct minimization in Hartree Fock and DFT - II) Converence o
I ) Direct minimization in Hartree Fock and DFT - II) Converence of Coupled Cluster Computational methods in many electron systems : MP2 Workshop Göteborg 2009 Ab initio computation Reliable computation
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
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 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 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 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
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 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 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 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 informationElectronic Structure Calculation - the backbone of computational material science
Electronic Structure Calculation - the backbone of computational material science Reinhold Schneider, MATHEON TU Berlin JvNeumann colloquium TU Munich 2012 Aim: Advanced Materials Numerical simulation
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 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 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 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 informationTensorization of Electronic Schrödinger Equation and Second Quantization
Tensorization of Electronic Schrödinger Equation and Second Quantization R. Schneider (TUB Matheon) Mariapfarr, 2014 Introduction Electronic Schrödinger Equation Quantum mechanics Goal: Calculation of
More informationELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS
FROM RESEARCH TO INDUSTRY 32 ème forum ORAP 10 octobre 2013 Maison de la Simulation, Saclay, France ELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS APPLICATION ON HPC, BLOCKING POINTS, Marc
More informationMathematical representations of quantum states
Mathematical representations of quantum states Eric CANCES Ecole des Ponts and INRIA, Paris, France IPAM, Los Angeles, Septembre 13-16, 2016 Outline of the tutorial 2 1. Quantum states of many-particle
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 informationDaubechies wavelets as a basis set for density functional pseudopotential calculations
Daubechies wavelets as a basis set for density functional pseudopotential calculations Luigi Genovese, Alexey Neelov, 2 Stefan Goedecker, 2 Thierry Deutsch, Seyed Alireza Ghasemi, 2 Alexander Willand,
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 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 informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
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 matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationSloppy Nuclear Energy Density Functionals: effective model optimisation. T. Nikšić and D. Vretenar
Sloppy Nuclear Energy Density Functionals: effective model optimisation T. Nikšić and D. Vretenar Energy Density Functionals the nuclear many-body problem is effectively mapped onto a one-body problem
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 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 informationMini-tutorial on Mathematical Aspects of Computational Quantum Chemistry SIAM MS 18, July 2018
Mini-tutorial on Mathematical Aspects of Computational Quantum Chemistry SIAM MS 18, July 2018 Benjamin Stamm Center for Computational Engineering Science, RWTH Aachen University, Aachen, Germany Outline
More informationFast Eigenvalue Solutions
Fast Eigenvalue Solutions Techniques! Steepest Descent/Conjugate Gradient! Davidson/Lanczos! Carr-Parrinello PDF Files will be available Where HF/DFT calculations spend time Guess ρ Form H Diagonalize
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 informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
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 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 informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
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 informationBigDFT tutorial. BigDFT tutorial FIRST YARMOUK SCHOOL. Thierry Deutsch. November, L_Sim - CEA Grenoble. Introduction DFT.
Big tutorial FIRST YARMOUK SCHOOL Big tutorial Thierry Deutsch L_Sim - CEA Grenoble November, 2. 2010 Outline 1 2 Density Functional Theory (quick view) 3 4 A basis for nanosciences: the Big project STREP
More informationLocalized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy
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 Joint work
More informationAdaptive Finite Element Approximations for Density Functional Models
HPPC Methods and Algorithms for Materials Defects, IMS Singapore, February 12, 2015 Adaptive Finite Element Approximations for Density Functional Models Aihui Zhou LSEC, Institute of Computational Mathematics
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 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 informationHandbook of Computational Quantum Chemistry
Handbook of Computational Quantum Chemistry David B. Cook Dept. of Chemistry University of Sheffield DOVER PUBLICATIONS, INC. Mineola, New York F Contents 1 Mechanics and molecules 1 1.1 1.2 1.3 1.4 1.5
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
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 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 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 informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
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 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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
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 informationThe Abinit project. Coding is based on modern software engineering principles
The Abinit project Abinit is a robust, full-featured electronic-structure code based on density functional theory, plane waves, and pseudopotentials. Abinit is copyrighted and distributed under the GNU
More informationNumerical tensor methods and their applications
Numerical tensor methods and their applications 8 May 2013 All lectures 4 lectures, 2 May, 08:00-10:00: Introduction: ideas, matrix results, history. 7 May, 08:00-10:00: Novel tensor formats (TT, HT, QTT).
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationRepresentation of operators in MADNESS (Multiresolution ADaptive Numerical Environment for Scientific Simulation)
Representation of operators in MADNESS (Multiresolution ADaptive Numerical Environment for Scientific Simulation) Gregory Beylkin University of Colorado at Boulder IPAM: Big Data Meets Computation workshop
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 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 Hartree-Fock approximation
Contents The Born-Oppenheimer approximation Literature Quantum mechanics 2 - Lecture 7 November 21, 2012 Contents The Born-Oppenheimer approximation Literature 1 The Born-Oppenheimer approximation 2 3
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 informationDFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
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 informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
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 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 informationTensor-Product Representation of Operators and Functions (7 introductory lectures) Boris N. Khoromskij
1 Everything should be made as simple as possible, but not simpler. A. Einstein (1879-1955) Tensor-Product Representation of Operators and Functions (7 introductory lectures) Boris N. Khoromskij University
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 informationMinimization 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 informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationFEM and sparse linear system solving
FEM & sparse linear system solving, Lecture 9, Nov 19, 2017 1/36 Lecture 9, Nov 17, 2017: Krylov space methods http://people.inf.ethz.ch/arbenz/fem17 Peter Arbenz Computer Science Department, ETH Zürich
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
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 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 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 informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
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 informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationDynamical low rank approximation in hierarchical tensor format
Dynamical low rank approximation in hierarchical tensor format R. Schneider (TUB Matheon) John von Neumann Lecture TU Munich, 2012 Motivation Equations describing complex systems with multi-variate solution
More informationQuantum Chemistry for Large Systems. Elias Rudberg
Quantum Chemistry for Large Systems Elias Rudberg Theoretical Chemistry School of Biotechnology Royal Institute of Technology Stockholm 2007 Quantum Chemistry for Large Systems Doctoral thesis c Elias
More informationGround-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices
Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices Heinz Siedentop Chennai, August 16, 2010 I. Introduction Atomic Schrödinger operator H N,Z := N n=1 self-adjointly realized
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
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 informationThe spectral action for Dirac operators with torsion
The spectral action for Dirac operators with torsion Christoph A. Stephan joint work with Florian Hanisch & Frank Pfäffle Institut für athematik Universität Potsdam Tours, ai 2011 1 Torsion Geometry, Einstein-Cartan-Theory
More informationKrylov Space Solvers
Seminar for Applied Mathematics ETH Zurich International Symposium on Frontiers of Computational Science Nagoya, 12/13 Dec. 2005 Sparse Matrices Large sparse linear systems of equations or large sparse
More informationImage restoration: numerical optimisation
Image restoration: numerical optimisation Short and partial presentation Jean-François Giovannelli Groupe Signal Image Laboratoire de l Intégration du Matériau au Système Univ. Bordeaux CNRS BINP / 6 Context
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 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 informationDensity Matrix Renormalization Group Algorithm - Optimization in TT-format
Density Matrix Renormalization Group Algorithm - Optimization in TT-format R Schneider (TU Berlin MPI Leipzig 2010 Density Matrix Renormalisation Group (DMRG DMRG algorithm is promininent algorithm to
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 informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationInner products. Theorem (basic properties): Given vectors u, v, w in an inner product space V, and a scalar k, the following properties hold:
Inner products Definition: An inner product on a real vector space V is an operation (function) that assigns to each pair of vectors ( u, v) in V a scalar u, v satisfying the following axioms: 1. u, v
More information