Hartree-Fock Algorithm

Total Page:16

File Type:pdf, Size:1020Kb

Hartree-Fock Algorithm Hartree-Fock algorithm I. Pre-computed quantities/notations 1. Overlap Matrix S Definition: ∗ 3 푆휇푣 = ∫ 푋휇 (푟⃗)푋푣(푟⃗)푑 푟 푤ℎ푒푟푒 푋휇(푟⃗)푖푠 푡ℎ푒 퐴푂 푏푎푠푖푠 푠푒푡 휇 = 1, … , 푁 푏푎푠푖푠 Properties: i) S is a NXN matrix where N is the dimension of basis. 1, 푖푓 휇 = 푣 ii) For orthonormal basis set 푆 = 훿 = { 휇푣 휇푣 0, 푖푓 휇 ≠ 푣 iii) S is a Hermitian Matrix and all elements in S are real: ∗ 푆휇푣 = 푆푣휇 = 푆푣휇 2. Core Hamiltonian h Definition: ∗ 3 ℎ휇푣 = ∫ 푋휇(푟⃗) ℎ̂푋푣(푟⃗)푑 푟 1 푤ℎ푒푟푒 ℎ̂ = − ∇2 + 푉̂ 푖푠 푡ℎ푒 퐾푖푛푒푡푖푐 푒푛푒푟푔푦 푎푛푑 푝표푡푒푛푡푖푎푙 푒푛푒푟푔푦 2 푖푛 푡ℎ푒 푓푖푒푙푑 표푓 푡ℎ푒 푛푢푐푙푒푖 푓표푟 푒푙푒푐푡푟표푛푠 Properties: h is a Hermitian matrix and all elements in h are real: ∗ ℎ휇푣 = ℎ푣휇 = ℎ푣휇 3. A critical quantity in the matrix for electron-electron interaction: Definition: The electron-electron interaction is described using the two-electron Integrals: ∗ ∗ 1 3 3 ⟨휇푣|휎휏⟩ = ∫ 푋휇(푟⃗⃗⃗1⃗)푋푣(푟⃗⃗⃗2⃗) 푋휎(푟⃗1)푋휏(푟⃗⃗⃗2⃗)푋 푑 푟1푑 푟2 푟12 Properties: In practice the AO orbitals we use are real. Therefore the integrals are the same under the permutations(휇 ↔ 휎 푎푛푑 푣 ↔ 휏 ). Moreover we can always permute coordinates 1 and 2 in the above integral which means we interchange (휇휎) ↔ (푣휏) 1 Overall the spatial integrals have an 8-fold permutation symmetry: ⟨휇푣|휎휏⟩ = ⟨휎푣|휇휏⟩ = ⟨휇휏|휎푣⟩ = ⟨휎휏|휇푣⟩ = ⟨푣휇|휏휎⟩ = ⟨휏휇|푣휎⟩ = ⟨푣휎|휏휇⟩ = ⟨휏휎|푣휇⟩ The antisymmetric integrals are defined as: ⟨휇푣||휎휏⟩ = ⟨휇푣|휎휏⟩ − ⟨휇푣|휏휎⟩ They have even more symmetry in the sense that interchanging two labels in the bra or in the ket just introduces a minum sign: We have: ⟨휇푣||휎휏⟩ = ⟨푣휇||휏휎⟩ = ⟨휎휏||휇푣⟩ = ⟨휇휏||휎푣⟩ = ⟨휏휎||푣휇⟩ = ⟨휏휇||푣휎⟩ And ⟨휇푣||휎휏⟩ = −⟨푣휇||휏휎⟩ = −⟨휎휏||푣휇⟩ = −⟨휇휏||푣휎⟩ = −⟨휏휎||휇푣⟩ = −⟨휏휇||휎푣⟩ The σ,μ,ν,τ labels for the orbitals should be re-interpreted as two spin-orbitals.(Due to its anti-symmetry property) II. HF equations 푭푪 = 푺푪푬 푪Ɨ푺푪 = 푰 푭흁풗 = 풉흁풗 + ∑⟨휇휎||푣휏⟩푫흈흉 휎,휏 1, 퐹표푟 표푐푐푢푝푖푒푑 표푟푏푖푡푎푙 푫 = 푪n푪햋 푛: 표푐푐푢푝푎푡푖표푛 푛푢푚푏푒푟 푛 = { 0, 퐹표푟 푢푛표푐푐푢푝푖푒푑 표푟푏푖푡푎푙 1 퐸 = 푻풓(풉 + 푭) ∗ 푫 퐻퐹 2 III. Deviation of Hartree-Fock equations 1. Choose a basis set and the wave function can be regarded as a linear combination of the basis set: For one occupied orbital |휑푎⟩ = ∑ |푋휆⟩퐶휆푎 휆 The total wave function: 휑푎(푟1), 휑푎(푟2), … 휑푎(푟푁) 휑 (푟 ), 휑 (푟 ), … 휑 (푟 ) 1 | 푏 1 푏 2 푏 푁 | |훹 ⟩ = … ( 휑 , … , 휑 푎푟푒 표푐푐푢푝푖푒푑 표푟푏푖푡푎푙푠) 퐻퐹 푁! | | 푎 푧 휑푧(푟1), 휑푧(푟2), … 휑푧(푟푁) 2 Different orbitals are orthonormal to each other: 1, 푎 = 푏 ⟨휑 |휑 ⟩ = ∫ 휑∗(푟⃗) ∗ 휑 (푟⃗)푑푟 = 훿 = { 푎 푏 푎 푏 푎푏 0, 푎 ≠ 푏 2. Calculate the total energy: 1 퐸 = ⟨훹 |Ĥ|훹 ⟩ = ∑⟨푎|ℎ|푎⟩ + ∑⟨푎푏||푎푏⟩ 퐻퐹 퐻퐹 퐻퐹 2 푎∈푎 푎,푏 1 = ∑ 퐶∗ 퐶 ⟨휇|ℎ|푣⟩ + ∑ 퐶∗ 퐶 퐶∗ 퐶 ⟨휇푣||휎휏⟩ 휇푣 푣푎 2 휇푎 휎푎 푣푏 휏푏 푎,휇,푣 푎,푏 3. Density Matrix: ∗ ∗ ∑ 퐶휇푎퐶푣푎 = ∑ 퐶푣푎퐶휇푎 = 푫풗흁 푎 푎 푫 = ∑ 푪 (푪Ɨ) = ∑ 푪 푛 푪햋 = 푪퐧푪햋 풗흁 풗풂 풂흁 풗풑 푝 풑흁 풂 푷:풐풗풆풓 풂풍풍 풐풓풃풊풕풂풍풔 1 표푐푐푢푝푖푒푑 푛 = { 푝 0 푢푛표푐푐푢푝푖푒푑 4. Substitute the density matrix into the expression of 퐸퐻퐹: ∑⟨푎|ℎ|푎⟩ = ∑⟨휇|ℎ|푣⟩퐷휇푣 = ∑ ℎ휇푣퐷푣휇 = 푻풓(푯푫) 푎 휇.푣 휇,푣 1 1 ∑⟨푎푏||푎푏⟩ = ∑ ⟨휇푣||휎휏⟩퐷 퐷 2 2 휎휇 휏푣 푎 휇,푣,휎,휏 ∗ 푫흈흁 = ∑ 퐶휎푎퐶휇푎 푎 ∗ 푫흉풗 = ∑ 퐶휏푏퐶푣푏 푏 1 퐸 = ∑ 풉 푫 + ∑ ⟨휇푣||휎휏⟩푫 푫 퐻퐹 흁풗 흁풗 2 흈흁 흉풗 휇,푣 휇,푣,휎,휏 5. The orthonormality of the orbitals in the AO basis is expressed by: ∗ ∗ Ɨ ⟨푎|푏⟩ = 훿푎푏 ⇒ ∑ 퐶휇푎⟨휇|푣⟩퐶푣푏 = 퐶푎휇푆휇푣퐶푣푏 = 푪 푺푪 = 푰 휇,푣 3 6. Derivation of Fock equation: FC=SCE (briefly): Lagrange’s method to optimize the parameters by minimizing the total energy as a function of C: 휁(푪) = 퐸 (푫) + (푪Ɨ푺푪 − 푰 ) ∗ 휦 퐻퐹 풂풃 푎푏 풂풃 휦풂풃 = 휀푎훿푎푏 (퐿푎푔푟푎푛푔푒 푀푢푙푡푖푝푙푖푒푟푠) IV. Special feature of Density Matrix: invariant under rotation of occupied orbitals 푫(풂) = ∑ 퐶 퐶∗ = ∑ 푪 (푪햋) 흁풗 휇푎 푣푎 흁풂 풂풗 푎 풂 Define the rotation operator: U U is a unitary operator: |푖⟩ = ∑ |푎⟩푈푎푖 푎,푏 Dirac notation: Bra: |푖⟩ = ∑|푎⟩푈푎푖 = ∑|휇⟩퐶휇푎푈푎푖 푎 휇,푎 Ket: ∗ ∗ ∗ ⟨푖| = ∑ 푈푖푎⟨푎| = ∑ 푈푖푎퐶푎휇⟨휇| 푎 휇,푎 ⟨푖|푗⟩ = ∑ 푈픣 ⟨푎|푏⟩ 푈 = (푼햋푼) = 푰 푎푖 푏푗 풊풋 풊풋 푎,푏 Under this rotational transformation the Density Matrix is invariant: 푫′ = 푼Ɨ푫푼 = 푼Ɨ(푪푛푪Ɨ)푼 ′ 햋 ∗ Ɨ Ɨ Ɨ Ɨ 푫흁풗 = 푼 (∑(퐶휇푖퐶푖푣)) 푼 = ∑ 퐶휇푎(푈푎푖푈푖푏)퐶푏푣 = ∑ 푪흁풂휹풂풃푪풃풗 = ∑ 푪흁풃푪풃풗 = 푫흁풗 푖 푎,푖,푏 풂,풃 푏 A deduce: Hence we have shown that molecular orbitals are not unique. One can define a different set of orbitals that are related to the original set by a unitary transformation (or a "rotation"). The density matrix is invariant under such a change of orbitals. Take water molecular as an example: As we learned in general chemistry, the orbital of water molecular can be expressed as 푠푝3 hybride. 4 One can express the orbitals in two ways: The first way is separate the orbitals into 4 localized orbitals : The second way is to express it to symmetry and anti-symmetry orbitals: 5 V. Algorithm for Hartree-Fock method 1. Read integrals from file: H,V,S,set D=0 initially 2. Cal_X: Calculate X from S: 푺 = 푼Ɨ푺푼(푑푖푎푔표푛푎푙푖푠푒 푆 푚푎푡푟푖푥) ퟏ 1 − 햋 − 푿 = 푺 ퟐ = 푼 푠 2푼 ퟏ − 퐗 = 퐒 ퟐ (퐗햋푺푿 = 푰) Where small s is the eigenvalues of Matrix S 3. Procedure: Calculate density matrix from fock matrix( initially Fμv = ℎ휇푣): 푭푪 = 푺푪휺 푫 = 푪풏푪햋 In practice: 푭′ = 푿햋푭푿 푭′푪′ = 푪′훆 푪 = 푿푪′ 푫 = 푪푛푪햋 4. Procedure: Calculate fock matrix from density matrix For closed shells: 푭흁풗 = 풉흁풗 + ∑(2⟨휇휎|푣휏⟩ − ⟨휇휎|푣휏⟩)푫흈흉 휎,휏 5. Self-consistency check__converagence: 풊풏 풐풖풕 if max (abs(퐃훍퐯 − 푫흁풗 )) < 푡ℎ푟푒푠ℎ표푙푑 where threshold is a constant (1e-7 e.g.) 6. Calculate energy 1 퐸 = 푻풓(풉 + 푭) ∗ 푫 2 6 The main program should be similar to this: Main {read_integrals(n_occ,n_total, S, h, V); Cal_X; D_in=0; ##Iterate: While(iter<=maximum and density matrix converged) {Calc_Fock(F,D_in); Calc_energy(F,D_in); Calc_density_form_Fock(D_out); Check_convergance;}; Return 0; }; VI. Unrestricted Hartree-Fock Problems: When solving open-shelled molecular by Hartree-Fock approximation, one need to take spin- orbitals into account. The way to deal with it is to divide the Fock Matrix into two part__α electron part and β electron part and calculate them individually. 7 Take oxygen molecule as an example, as we know that the MO for O2 molecule is as fellow: CHEM 209 course note, SCOTT HOPKINS One separate the filled orbital into spin-up orbitals ( α orbitals) and spin-down orbitals( β orbitals): α orbitals β orbitals 8 Compared to α orbitals, the energy level of β orbitals are slight higher. This is due to the difference of the exchange terms in their Fock matrix: i. Density Matrix: (휶) ∗ 푫흁풗 = ∑ 퐶휇푖퐶푣푖 푖∈훼−푠푝푖푛 표푐푐 (휷) ∗ 푫흁풗 = ∑ 퐶휇푖퐶푣푖 푖∈훽 표푐푐 In this way one can separate the Density Matrix into α-spin orbital part and β-spin orbital part. ii. Fock Matrix: 푭흁풗 = 풉흁풗 + ∑⟨휇휎|푣휏⟩푫흉흈 − ∑⟨휇휎|휏푣⟩푫흉흈 휎,휏 휎,휏 I II Coulomb interaction Exchange term One can separate the interaction term of Fock matrix into Coulomb interaction part and Exchange term. The first term can be understood as coulomb potential classically, while the second term is due to the anti-symmetry property of electrons and has no classical counterpart. iii. Coulomb interaction: (휶) (휷) 푱흁풗 = ∑⟨휇휎|푣휏⟩푫흉흈 = ∑ ⟨휇휎훼|푣휏훼⟩푫흉흈 + ∑ ⟨휇휎훽|푣휏훽⟩푫흉흈 휎,휏 휎,휏 휎,휏 훼−푠푝푖푛 훽−푠푝푖푛 (휶) (휷) (휶) ⟨휇휎|푣휏⟩ [푫흉흈 + 푫흉흈 ] 푱 흁풗 = ∑ 풔풑풂풕풊풂풍 휎,휏 푡표푡푎푙 푫흉흈 푆푝푎푡푖푎푙 (풔풑풂풕풊풂풍) (풔풖풎) 푱흁풗 = ∑ ⟨휇휎|푣휏⟩푫흉흈 퐷푖푟푒푐푡 퐶표푢푙표푚푏 푃표푡푒푛푡푖푎푙 휎,휏 푆푝푎푡푖푎푙 9 iv. Exchange interactions (휶) (휶) 푫흉흈 → ∑ ⟨휇훼휎훼|휏훼푣훼⟩푫흉흈 휎,휏 훼−푠푝푖푛 (휶) (휶) → 푲흁풗 = ∑ ⟨휇휎|휏푣⟩푫흉흈 휎.휏 푠푝푎푡푖푎푙 (휷) (휷) 푫흉흈 → ∑ ⟨휇훽휎훽|휏훽푣훽⟩푫흉흈 휎,휏 훽−푠푝푖푛 (휷) (휷) → 푲흁풗 = ∑ ⟨휇휎|휏푣⟩푫흉흈 휎.휏 푠푝푎푡푖푎푙 V. General HF algorithm 1. If 푫(휶) = 푫(휷), then it is a closed shell molecule. And 푭(휶) = 푭(휷), one can only need to calculate Fock Matrix for either alpha or beta orbital once. For closed shell: (휶) 푭흁풗 = 풉흁풗 + ∑(2⟨휇휎|푣휏⟩ − ⟨휇휎|휏푣⟩푫흉흈 휎,휏 2. For unrestricted HF calculations, one need to treat different orbitals correspond to alpha spins and beta spins separately. The total energy should be the addition of the two orbitals: 1 (흈) (흈) (흈) 퐸 = ∑ (풉 + 푭 ) 푫 퐻퐹 휎=(훼,훽) 2 흁풗 흁풗 흁풗 The UHF can give the lowest energy of the molecule, but the spin-symmetry was broken. The wavefunction optimized in this method is not the eigenfunction of spin: ퟐ 푺 |훹퐻퐹⟩ ≠ 푠(푠 + 1)|훹퐻퐹⟩ 10 .
Recommended publications
  • Arxiv:1809.04476V2 [Physics.Chem-Ph] 18 Oct 2018 (Sub)States
    Quantum System Partitioning at the Single-Particle Level Quantum System Partitioning at the Single-Particle Level Adrian H. M¨uhlbach and Markus Reihera) ETH Z¨urich, Laboratorium f¨ur Physikalische Chemie, Vladimir-Prelog-Weg 2, CH-8093 Z¨urich,Switzerland (Dated: 17 October 2018) We discuss the partitioning of a quantum system by subsystem separation through unitary block- diagonalization (SSUB) applied to a Fock operator. For a one-particle Hilbert space, this separation can be formulated in a very general way. Therefore, it can be applied to very different partitionings ranging from those driven by features in the molecular structure (such as a solute surrounded by solvent molecules or an active site in an enzyme) to those that aim at an orbital separation (such as core-valence separation). Our framework embraces recent developments of Manby and Miller as well as older ones of Huzinaga and Cantu. Projector-based embedding is simplified and accelerated by SSUB. Moreover, it directly relates to decoupling approaches for relativistic four-component many-electron theory. For a Fock operator based on the Dirac one-electron Hamiltonian, one would like to separate the so-called positronic (negative-energy) states from the electronic bound and continuum states. The exact two-component (X2C) approach developed for this purpose becomes a special case of the general SSUB framework and may therefore be viewed as a system- environment decoupling approach. Moreover, for SSUB there exists no restriction with respect to the number of subsystems that are generated | in the limit, decoupling of all single-particle states is recovered, which represents exact diagonalization of the problem.
    [Show full text]
  • A New Scalable Parallel Algorithm for Fock Matrix Construction
    A New Scalable Parallel Algorithm for Fock Matrix Construction Xing Liu Aftab Patel Edmond Chow School of Computational Science and Engineering College of Computing, Georgia Institute of Technology Atlanta, Georgia, 30332, USA [email protected], [email protected], [email protected] Abstract—Hartree-Fock (HF) or self-consistent field (SCF) ber of independent tasks and using a centralized dynamic calculations are widely used in quantum chemistry, and are scheduler for scheduling them on the available processes. the starting point for accurate electronic correlation methods. Communication costs are reduced by defining tasks that Existing algorithms and software, however, may fail to scale for large numbers of cores of a distributed machine, particularly in increase data reuse, which reduces communication volume, the simulation of moderately-sized molecules. In existing codes, and consequently communication cost. This approach suffers HF calculations are divided into tasks. Fine-grained tasks are from a few problems. First, because scheduling is completely better for load balance, but coarse-grained tasks require less dynamic, the mapping of tasks to processes is not known communication. In this paper, we present a new parallelization in advance, and data reuse suffers. Second, the centralized of HF calculations that addresses this trade-off: we use fine- grained tasks to balance the computation among large numbers dynamic scheduler itself could become a bottleneck when of cores, but we also use a scheme to assign tasks to processes to scaling up to a large system [14]. Third, coarse tasks reduce communication. We specifically focus on the distributed which reduce communication volume may lead to poor load construction of the Fock matrix arising in the HF algorithm, balance when the ratio of tasks to processes is close to one.
    [Show full text]
  • Domain Decomposition and Electronic Structure Computations: a Promising Approach
    Domain Decomposition and Electronic Structure Computations: A Promising Approach G. Bencteux1,4, M. Barrault1, E. Canc`es2,4, W. W. Hager3, and C. Le Bris2,4 1 EDF R&D, 1 avenue du G´en´eral de Gaulle, 92141 Clamart Cedex, France {guy.bencteux,maxime.barrault}@edf.fr 2 CERMICS, Ecole´ Nationale des Ponts et Chauss´ees, 6 & 8, avenue Blaise Pascal, Cit´eDescartes, 77455 Marne-La-Vall´ee Cedex 2, France, {cances,lebris}@cermics.enpc.fr 3 Department of Mathematics, University of Florida, Gainesville, FL 32611-8105, USA, [email protected] 4 INRIA Rocquencourt, MICMAC project, Domaine de Voluceau, B.P. 105, 78153 Le Chesnay Cedex, FRANCE Summary. We describe a domain decomposition approach applied to the spe- cific context of electronic structure calculations. The approach has been introduced in [BCH06]. We survey here the computational context, and explain the peculiari- ties of the approach as compared to problems of seemingly the same type in other engineering sciences. Improvements of the original approach presented in [BCH06], including algorithmic refinements and effective parallel implementation, are included here. Test cases supporting the interest of the method are also reported. It is our pleasure and an honor to dedicate this contribution to Olivier Pironneau, on the occasion of his sixtieth birthday. With admiration, respect and friendship. 1 Introduction and Motivation 1.1 General Context Numerical simulation is nowadays an ubiquituous tool in materials science, chemistry and biology. Design of new materials, irradiation induced damage, drug design, protein folding are instances of applications of numerical sim- ulation. For convenience we now briefly present the context of the specific computational problem under consideration in the present article.
    [Show full text]
  • Matrix Algebra for Quantum Chemistry
    Matrix Algebra for Quantum Chemistry EMANUEL H. RUBENSSON Doctoral Thesis in Theoretical Chemistry Stockholm, Sweden 2008 Matrix Algebra for Quantum Chemistry Doctoral Thesis c Emanuel Härold Rubensson, 2008 TRITA-BIO-Report 2008:23 ISBN 978-91-7415-160-2 ISSN 1654-2312 Printed by Universitetsservice US AB, Stockholm, Sweden 2008 Typeset in LATEX by the author. Abstract This thesis concerns methods of reduced complexity for electronic structure calculations. When quantum chemistry methods are applied to large systems, it is important to optimally use computer resources and only store data and perform operations that contribute to the overall accuracy. At the same time, precarious approximations could jeopardize the reliability of the whole calcu- lation. In this thesis, the selfconsistent eld method is seen as a sequence of rotations of the occupied subspace. Errors coming from computational ap- proximations are characterized as erroneous rotations of this subspace. This viewpoint is optimal in the sense that the occupied subspace uniquely denes the electron density. Errors should be measured by their impact on the over- all accuracy instead of by their constituent parts. With this point of view, a mathematical framework for control of errors in HartreeFock/KohnSham calculations is proposed. A unifying framework is of particular importance when computational approximations are introduced to eciently handle large systems. An important operation in HartreeFock/KohnSham calculations is the calculation of the density matrix for a given Fock/KohnSham matrix. In this thesis, density matrix purication is used to compute the density matrix with time and memory usage increasing only linearly with system size. The forward error of purication is analyzed and schemes to control the forward error are proposed.
    [Show full text]
  • Problem Set #1 Solutions
    PROBLEM SET #2 SOLUTIONS by Robert A. DiStasio Jr. Q1. 1-particle density matrices and idempotency. (a) A matrix M is said to be idempotent if M 2 = M . Show from the basic definition that the HF density matrix is idempotent when expressed in an orthonormal basis. An element of the HF density matrix is given as (neglecting the factor of two for the restricted closed-shell HF density matrix): * PCCμν = ∑ μi ν i . (1) i In matrix form, (1) is simply equivalent to † PCC= occ occ (2) where Cocc is the occupied block of the MO coefficient matrix. To reformulate the nonorthogonal Roothaan-Hall equations FC= SCε (3) into the preferred eigenvalue problem ~~ ~ FC= Cε (4) the AO basis must be orthogonalized. There are several ways to orthogonalize the AO basis, but I will only discuss symmetric orthogonalization, as this is arguably the most often used procedure in computational quantum chemistry. In the symmetric orthogonalization procedure, the MO coefficient matrix in (3) is recast as ~ + 1 − 1 ~ CSCCSC=2 ⇔ = 2 (5) which can be substituted into (3) yielding the eigenvalue form of the Roothaan-Hall equations in (4) after the following algebraic manipulation: FC = SCε − 1 ~ − 1 ~ FSCSSC( 2 )(= 2 )ε − 1 − 1 ~ −1 − 1 ~ SFSCSSSC2 ( 2 )(= 2 2 )ε . (6) −1 − 1 ~ −1 − 1 ~ ()S2 FS2 C= () S2 SS2 Cε ~~ ~ FC = Cε ~ ~ −1 − 1 −1 − 1 In (6), the dressed Fock matrix, F , was clearly defined as F≡ S2 FS 2 and S2 SS2 = S 0 =1 was used. In an orthonormal basis, therefore, the density matrix of (2) takes on the following form: ~ ~† PCC= occ occ .
    [Show full text]
  • Relativistic Cholesky-Decomposed Density Matrix
    Relativistic Cholesky-decomposed density matrix MP2 Benjamin Helmich-Paris,1, 2, a) Michal Repisky,3 and Lucas Visscher2 1)Max-Planck-Institut f¨urKohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 M¨ulheiman der Ruhr 2)Section of Theoretical Chemistry, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands 3)Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromø Norway (Dated: 15 October 2018) In the present article, we introduce the relativistic Cholesky-decomposed density (CDD) matrix second-order Møller–Plesset perturbation theory (MP2) energies. The working equations are formulated in terms of the usual intermediates of MP2 when employing the resolution-of-the-identity approximation (RI) for two-electron inte- grals. Those intermediates are obtained by substituting the occupied and vir- tual quaternion pseudo-density matrices of our previously proposed two-component atomic orbital-based MP2 (J. Chem. Phys. 145, 014107 (2016)) by the corresponding pivoted quaternion Cholesky factors. While working within the Kramers-restricted formalism, we obtain a formal spin-orbit overhead of 16 and 28 for the Coulomb and exchange contribution to the 2C MP2 correlation energy, respectively, compared to a non-relativistic (NR) spin-free CDD-MP2 implementation. This compact quaternion formulation could also be easily explored in any other algorithm to compute the 2C MP2 energy. The quaternion Cholesky factors become sparse for large molecules and, with a block-wise screening, block sparse-matrix multiplication algorithm, we observed an effective quadratic scaling of the total wall time for heavy-element con- taining linear molecules with increasing system size.
    [Show full text]
  • Lecture 16 Hartree Fock Implementation
    Quantum Chemistry Lecture 16 Implementation of Hartree-Fock Theory NC State University The linear combination of atomic orbitals (LCAO) Calculations of the energy and properties of molecules requires hydrogen-like wave functions on each of the nuclei. The Hartree-Fock method begins with assumption That molecular orbitals can be formed as a linear combination of atomic orbitals. The basis functions fm are hydrogen-like atomic orbitals that have been optimized by a variational procedure. The HF procedure is a variational procedure to minimize the coefficients Cmi. Note that we use the index m for atomic orbitals and i or j for molecular orbitals. Common types of atomic orbitals Slater-type orbitals (STOs) The STOs are like hydrogen atom wave functions. The problem with STOs arises in multicenter integrals. The Coulomb and exchange integrals involve electrons on different nuclei and so the distance r has a different origin. Gaussian-type orbitals (GTOs) Gaussian orbitals can be used to mimic the shape of exponentials, i.e. the form of the solutions for the hydrogen atom. Multicenter Gaussian integrals can be solved analytically. STOs vs GTOs GTOs are mathematically easy to work with STOs vs GTOs GTOs are mathematically easy to work with, but the shape of a Gaussian is not that similar to that of an exponential. STOs vs GTOs Therefore, linear combinations of Gaussians are used to imitate the shape of an exponential. Shown is a representation of the 3-Gaussian model of a STO. Double-zeta basis sets Since the remaining atoms have a different exponential dependence than hydrogen it is often convenient to include more parameters.
    [Show full text]
  • Fock Matrix Construction for Large Systems
    Fock Matrix Construction for Large Systems Elias Rudberg Theoretical Chemistry Royal Institute of Technology Stockholm 2006 Fock Matrix Construction for Large Systems Licentiate thesis c Elias Rudberg, 2006 ISBN 91-7178-535-3 ISBN 978-91-7178-535-0 Printed by Universitetsservice US AB, Stockholm, Sweden, 2006 Typeset in LATEX by the author Abstract This licentiate thesis deals with quantum chemistry methods for large systems. In particu- lar, the thesis focuses on the efficient construction of the Coulomb and exchange matrices which are important parts of the Fock matrix in Hartree–Fock calculations. The meth- ods described are also applicable in Kohn–Sham Density Functional Theory calculations, where the Coulomb and exchange matrices are parts of the Kohn–Sham matrix. Screening techniques for reducing the computational complexity of both Coulomb and exchange com- putations are discussed, as well as the fast multipole method, used for efficient computation of the Coulomb matrix. The thesis also discusses how sparsity in the matrices occurring in Hartree–Fock and Kohn– Sham Density Functional Theory calculations can be used to achieve more efficient storage of matrices as well as more efficient operations on them. As an example of a possible type of application, the thesis includes a theoretical study of Heisenberg exchange constants, using unrestricted Kohn–Sham Density Functional Theory calculations. iii Preface The work presented in this thesis has been carried out at the Department of Theoretical Chemistry, Royal Institute of Technology, Stockholm, Sweden. List of papers included in the thesis Paper 1 Efficient implementation of the fast multipole method, Elias Rudberg and Pawe l Sa lek,J.
    [Show full text]
  • Theoretical Studies of Charge Recombination Reactions
    Theoretical studies of charge recombination reactions Sifiso Musa Nkambule Akademisk avhandling f¨oravl¨aggandeav licentiatexamen vid Stockholms universitet, Fysikum. January, 2015. Abstract This thesis is based on theoretical studies that have been done for low-energy reactions involving small molecular systems. It is mainly focusing on mutual neutralization of oppositely charged ions and dissociative recombination of molecular ions. Both reactions involve highly excited electronic states that are coupled to each other. Employing electronic structure meth- ods, the electronic states relevant to the reactions are computed. It is necessary to go beyond the Born-Oppenheimer approximation and include non-adiabatic effects. This non-adiabatic behaviour plays a crucial role in driving the reactions. Schemes on how to go about and include these coupling elements are discussed in the thesis, which include computations of non-adiabatic and electronic couplings of the resonant states to the ionisation continuum. The nuclear dynamics are studied either semi-classically, using the Landau-Zener method or quantum mechanically, employing the time-independent and time-dependant Schr¨odinger + equations. Reactions studied here are mutual neutralization in collision of H + H−, where both the semi-classical and quantum mechanical methods are employed. Total and differential cross sections are computed for all hydrogen isotopes. We also perform quantum mechanical + studies of mutual neutralization in the collision of He + H−. For this system not only the non-adiabatic couplings among the neutral states have to be considered, but also the electronic coupling to the ionization continuum. Worth noting in the two systems is that reactants and the final products belong to different electronic states, which a coupled together, hence for the reaction to take place there must be non-adiabatic transition.
    [Show full text]
  • Quaternion Symmetry in Relativistic Molecular Calculations: the Dirac–Hartree–Fock Method T
    Quaternion symmetry in relativistic molecular calculations: The Dirac–Hartree–Fock method T. Saue, H. Jensen To cite this version: T. Saue, H. Jensen. Quaternion symmetry in relativistic molecular calculations: The Dirac–Hartree– Fock method. Journal of Chemical Physics, American Institute of Physics, 1999, 111 (14), pp.6211- 6222. 10.1063/1.479958. hal-03327446 HAL Id: hal-03327446 https://hal.archives-ouvertes.fr/hal-03327446 Submitted on 27 Aug 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 14 8 OCTOBER 1999 Quaternion symmetry in relativistic molecular calculations: The Dirac–Hartree–Fock method T. Sauea) and H. J. Aa Jensen Department of Chemistry, University of Southern Denmark—Main campus: Odense University, DK-5230 Odense M, Denmark ͑Received 5 April 1999; accepted 9 July 1999͒ A symmetry scheme based on the irreducible corepresentations of the full symmetry group of a molecular system is presented for use in relativistic calculations. Consideration of time-reversal symmetry leads to a reformulation of the Dirac–Hartree–Fock equations in terms of quaternion algebra. Further symmetry reductions due to molecular point group symmetry are then manifested by a descent to complex or real algebra.
    [Show full text]
  • Hartree-Fock Theory
    Hartree-Fock theory Morten Hjorth-Jensen1 National Superconducting Cyclotron Laboratory and Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA1 May 16-20 2016 c 2013-2016, Morten Hjorth-Jensen. Released under CC Attribution-NonCommercial 4.0 license Why Hartree-Fock? Hartree-Fock (HF) theory is an algorithm for finding an approximative expression for the ground state of a given Hamiltonian. The basic ingredients are I Define a single-particle basis f αg so that HF h^ α = "α α with the Hartree-Fock Hamiltonian defined as HF HF h^ = t^+u ^ext +u ^ HF I The term u^ is a single-particle potential to be determined by the HF algorithm. HF I The HF algorithm means to choose u^ in order to have HF hH^i = E = hΦ0jH^jΦ0i that is to find a local minimum with a Slater determinant Φ0 being the ansatz for the ground state. HF I The variational principle ensures that E ≥ E0, with E0 the exact ground state energy. Why Hartree-Fock? We will show that the Hartree-Fock Hamiltonian h^HF equals our definition of the operator f^ discussed in connection with the new definition of the normal-ordered Hamiltonian (see later lectures), that is we have, for a specific matrix element HF X hpjh^ jqi = hpjf^jqi = hpjt^+u ^extjqi + hpijV^jqiiAS ; i≤F meaning that HF X hpju^ jqi = hpijV^jqiiAS : i≤F The so-called Hartree-Fock potential u^HF brings an explicit medium dependence due to the summation over all single-particle states below the Fermi level F .
    [Show full text]
  • Restricted Closed Shell Hartree Fock Roothaan Matrix Method Applied to Helium Atom Using Mathematica Introduction the Hartree Fo
    European J of Physics Education Volume 5 Issue 1 2014 Acosta, Tapia, Cab Restricted Closed Shell Hartree Fock Roothaan Matrix Method Applied to Helium Atom using Mathematica César R. Acosta*1 J. Alejandro Tapia1 César Cab1 1Av. Industrias no Contaminantes por anillo periférico Norte S/N Apdo Postal 150 Cordemex Physics Engineering Department Universidad Autónoma de Yucatán *[email protected] (Received: 20.12. 2013, Accepted: 13.02.2014) Abstract Slater type orbitals were used to construct the overlap and the Hamiltonian core matrices; we also found the values of the bi-electron repulsion integrals. The Hartree Fock Roothaan approximation process starts with setting an initial guess value for the elements of the density matrix; with these matrices we constructed the initial Fock matrix. The Mathematica software was used to program the matrix diagonalization process from the overlap and Hamiltonian core matrices and to make the recursion loop of the density matrix. Finally we obtained a value for the ground state energy of helium. Keywords: STO, Hartree Fock Roothaan approximation, spin orbital, ground state energy. Introduction The Hartree Fock method (Hartree, 1957) also known as Self Consistent Field (SCF) could be made with two types of spin orbital functions, the Slater Type Orbitals (STO) and the Gaussian Type Orbital (GTO), the election is made in function of the integrals of spin orbitals implicated. We used the STO to construct the wave function associated to the helium atom electrons. The application of the quantum mechanical theory was performed using a numerical method known as the finite element method, this gives the possibility to use the computer for programming the calculations quickly and efficiently (Becker, 1988; Brenner, 2008; Thomee, 2003; Reddy, 2004).
    [Show full text]