Fundamentals of Quantum Mechanics for Chemistry
G´eraldMONARD
Equipe de Chimie et Biochimie Th´eoriques UMR 7565 CNRS - Universit´eHenri Poincar´e Facult´edes Sciences - B.P. 239 54506 Vandœuvre-les-Nancy Cedex - FRANCE
http://www.monard.info/ Outline . . . 1. Fundamentals of Quantum Mechanics for Chemistry Hartree-Fock methods Density Functional Theory The QM scaling problem Semiempirical methods Molecular Mechanics 2. Fundamentals of QM/MM methods Partionning QM/MM interactions Cutting covalent bonds ONIOM Some available software 3. Selected QM/MM applications Solvent effects Spectroscopy Biochemistry . . . Outline 4. Fundamentals of Linear Scaling methods QM Bottlenecks General ideas and solutions Some available software 5. Focus on some Linear Scaling methods CG-DMS Mozyme Divide & Conquer 6. Selected Linear Scaling applications Energy Decomposition; Charge Transfer & Polarization Born-Oppenheimer Molecular Dynamics 7. Parallelization of QM/MM and Linear Scaling methods Fundamentals of Quantum Mechanics for Chemistry (1)
Some approximations And there was the Schr¨odinger equation. . .
H0Ψ0 = E0Ψ0
where:
H0 is an Hamiltonian operator that describes a molecular system
Ψ0 is a wavefunction (solution of the Schr¨odingerequation) that describe a state of the system
E0 the energy associated to Ψ0
1 equation + 2 unknowns (given H0) = +∞ solutions! From now on: ground state closed shell non-relativistic Fundamentals of Quantum Mechanics for Chemistry (2)
Born-Oppenheimer approximation nuclei are fixed point charges only electrons are represented by a wavefunction Ψ
HΨ = Eelec Ψ (1)
H = Te + VeN + Vee
1 −Z 1 K (2) = − ∑∆i + ∑∑ + ∑∑ 2 i i K riK i i>j rij | {z } | {z } | {z } kinetic energy e−-nuclei inter. e−–e− inter.
E = Eelec + Enuclei Z Z = hΨ|H|Ψi + ∑ ∑ K L K L>K RKL Fundamentals of Quantum Mechanics for Chemistry (3)
Orbital approximation each electron is described by a mono-electronic wavefunction: the Molecular Orbital (MO)
Ψ(1,2,...,n) = ψ1(1)ψ2(2)...ψn(n)
all MOs are combined in a Slater determinant (Pauli principle)
Linear Combination of Atomic Orbitals (LCAO)
Each MO ψi is developed on a basis set of functions φµ : the Atomic Orbitals (AO) AOs ψi = ∑ cµi φµ µ
the real coefficients ciµ are the unknown of the problem Fundamentals of Quantum Mechanics for Chemistry (4)
Variational Principle
The eletronic energy Eelec corresponds to a minimum with respect to each MO ψi ∂E ∀i elec = 0 ∂ψi Hartree-Fock equations
∀i Fψi = εi ψi
where F is the mono-electronic Fock operator The Hartree-Fock method (1)
The Fock operator
c F(1) = H (1) + ∑[Jj (1) − Kj (1)] j
Hc(1) is the one-electron core Hamiltonian
c 1 ZK H (1) = − ∆1 − ∑ 2 K R1K
Jj (1) is the Coulomb operator, Z ∗ 1 Jj (1) = ψj (2) ψj (2)dτ2 (2) r12
Kj (1) is the exchange operator Z ∗ 1 Kj (1)ψi (1) = ψj (1) ψj (2) ψi (2)dτ2 (2) r12 The Hartree-Fock method (2)
The Roothan-Hall equations N electrons, φµ AO basis set, closed shell, ground state occupied MO = 2 electrons; virtual MO = 0 electron
AO ψi = ∑ cµi φµ µ
the Hartree-Fock equations can be re-written:
FC = SCε
with
C the matrix of cµi coefficients ε the diagonal energy matrix S the overlap matrix:
Sµν =< φµ |φν > The Hartree-Fock method (3) The Density Matrix
From the MO coefficient cµi , it is possible to build a density matrix whose elements are:
MO Pµν = ∑ nj cµj cνj with nj = 0 or 2 (occupation number) j The Fock matrix
AO AO 1 F = Hc + P (µν|λη) − (µη|λν) µν µν ∑∑ λη 2 λ η with: Z Z 1 (µν|λη) = φµ (1)φν (1) φλ (2)φη (2)dr1dr2 (1) (2) r12 The Hartree-Fock energy
1 Eelec = ∑∑Pµν [Hµν + Fµν ] 2 µ ν The Hartree-Fock method (4)
The Hartree-Fock algorithm 1. Compute mono- and bielectronic integrals 2. Build core hamiltonian (invariant) Hc 3. Guess an initial density matrix 4. Build the Fock matrix F 5. Orthogonal transformation using S1/2
F0C0 = εC0
6. Diagonalization of the Fock matrix F0 The C0 coefficients are obtained 7. Inverse transformation C0 → C 8. Build the new density matrix Back to 4. unless convergence The Hartree-Fock method (5)
The Hartree-Fock method is an ab initio method No (empirical) parameters ¡ ab initio method Orbital approximation ¡ No electronic correlation
Other ab initio methods (post-Hartree-Fock methods) Møller-Plesset Perturbation Theory (MP2, MP4, etc) Configuration Interaction (CI) Coupled Cluster (e.g. CCSD(T)) MultiConfigurational Self-Consistent-Field (MCSCF) Density Functional Theory (1)
Background DFT (Density Functional Theory) methods are (almost) ab initio methods which include electronic correlation at a cost similar to a Hartree-Fock calculation. In most cases, a DFT calculation is even less costly than a HF calculation. DFT methods relie on the Hohenberg-Kohn theorem (1964) which states that the ground state energy E of a system is a functional of the electronic density of this system, ρ(~r). Any electronic density ρ0(~r) other than the real electronic density will necessary lead to a higher energy. (variational principle) Density Functional Theory (2)
A different approach To the opposite of ab initio methods, DFT methods try to find a simple 3-dimensional ρ(~r) function and not a complex 3N-dimensional wave function.
ρ : R3 −→ R ~r 7−→ ρ(~r) From the Hohenberg-Kohn theorem, the energy E depends on the electronic density. It is said that E is a functional of the electronic density: