View metadata, citation and similar papers at core.ac.uk brought to you by CORE
provided by DSpace at Waseda University
Efficient Two-Component Relativistic Method to Obtain Electronic and Molecular Structures of Heavy-Element Systems
ͷిࢠঢ়ଶ͓ΑͼࢠߏͷͨΊͷܥॏݩૉ ޮతͳ 2 ૬ର๏
February 2017
Waseda University Graduate School of Advanced Science and Engineering Department of Chemistry and Biochemistry, Research on Electronic State Theory
Yuya NAKAJIMA தౢɹ༟
Contents
Chapter 1 General introduction 1
References ...... 6
Chapter 2 Theoretical background 7
2.1 Dirac Hamiltonian ...... 7
2.2 The IODKH method ...... 9
2.3 The LUT scheme ...... 12
2.4 Spin-free and spin-dependent formalisms ...... 14
2.5 The FCP method ...... 16
2.6 The GHF method ...... 19
2.7 The Analytical energy derivative for the GHF method . . . 21
References ...... 24
Chapter 3 Analytical energy gradient for spin-free infinite-order
Douglas–Kroll–Hess method with local unitary transfor-
mation 25
3.1 Introduction ...... 25
3.2 Theory and implementation ...... 27
3.2.1 Energy gradient for IODKH ...... 27
3.2.2 Analytical derivative for the space transformation
matrices ...... 30
3.2.3 Energy gradient for LUT-IODKH ...... 32
3.2.4 Implementation ...... 33
i 3.3 Numerical assessments ...... 34
3.3.1 Computational details ...... 34
3.3.2 Numerical gradient values ...... 37
3.3.3 Accuracies of IODKH/C and LUT-IODKH/C methods 38
3.3.4 Computational cost of the LUT scheme ...... 45
3.3.5 Metal complexes ...... 47
3.3.6 Heavier analogues of ethylene ...... 48
3.3.7 Harmonic frequencies of diatomic molecules . . . . 49
3.4 Conclusion ...... 50
References ...... 52
Chapter 4 Implementation of spin-dependent relativistic analytical
energy gradient 57
4.1 Introduction ...... 57
4.2 Implementation ...... 60
4.3 Numerical assessments ...... 61
4.3.1 Computational details ...... 61
4.3.2 Accuracy of SD-IODKH and LUT-SD-IODKH . . . . 63
4.3.3 Computational cost of (LUT-)SD-IODKH method . . 67 − 4.3.4 Application in fac Ir(ppy3)...... 70 4.4 Conclusion ...... 72
References ...... 74
Chapter 5 Implementation of LUT-IODKH in GAMESS program 79
5.1 Introduction ...... 79
5.2 Implementation ...... 79
5.2.1 Relativistic correction in GAMESS program . . . . . 79
5.2.2 Threshold setting in the IOTC method ...... 80
ii 5.2.3 Integral evaluation at quadruple precision ...... 81
5.2.4 Combination with DC ...... 82
5.2.5 Input options for LUT-IOTC ...... 83
5.2.6 Major capabilities of LUT-IOTC ...... 83
5.3 Numerical assessment ...... 84
5.3.1 Computational details ...... 84
5.3.2 Total energies of heavy atoms and molecules . . . . 85
5.3.3 Threshold dependence of IOTC and LUT-IOTC . . . 86
5.3.4 Computational cost of LUT-IOTC ...... 87
5.4 Conclusion ...... 89
References ...... 91
Chapter 6 Relaxation of core orbitals in the frozen core potential treat-
ment 93
6.1 Introduction ...... 93
6.2 Theory and implementation ...... 94
6.2.1 FCP with relaxation of core electrons ...... 94
6.2.2 Implementation ...... 96
6.3 Numerical assessments ...... 98
6.3.1 Computational details ...... 98
6.3.2 Computational cost of FCP-CR ...... 99
6.3.3 Accuracy of FCP-CR ...... 103
6.3.4 Core ionization energy and core level shift ...... 104
6.3.5 Accuracy of an iterative procedure between valence
and core calculations ...... 107
6.4 Conclusion ...... 108
Appendix 6.A Dependence of core ionization potential energies 111
iii References ...... 113
Chapter 7 Relativistic effect on enthalpy of formation for transition
metal complexes 117
7.1 Introduction ...... 117
7.2 Computational details ...... 119
7.3 Results and discussion ...... 121
7.3.1 Accuracy of WFT ...... 121
7.3.2 Functional dependence ...... 123
7.3.3 Geometry difference between PP and AE methods . 124
7.3.4 Effect of the levels of relativistic Hamiltonians . . . . 124
7.3.5 Contribution of frozen core orbitals ...... 126
7.4 Conclusion ...... 128
References ...... 130
Chapter 8 General Conclusion 133
Acknowledgments 137
List of Achievements 139
iv List of abbreviations
2c two-component
4c four-component
AE all-electron
AO atomic orbital
BSS Barysz-Sadlej-Snijders
CA composite approach
CCSD coupled cluster singles and doubles
CCSD(T) coupled cluster singles, doubles, and perturbative triples
CIE core ionization energy
CLS core level shift
CPU central processing unit
DC divide-and-conquer
DFT density functional theory
DLU local approximation to the unitary decoupling transformation dmpe 1,2-bis(dimethylphosphino)ethane dOEI derivative of one-electron integral dppe 1,2-bis((pentafluorophenyl)phosphino)ethane dTEI derivative of two-electron integral
ECP effective core potential
ESC elimination of small components
F-dppe 1,2-bis-((pentafluorophenyl)phosphino)ethane
FCP frozen core potential
v FCP-CR frozen-core potential with relaxation of core electrons
FW Fouldy-Wouthuysen
GAMESS General Atomic and Molecular Electronic Structure System
GHF general Hartree–Fock
HF Hartree–Fock
HFR Hartree–Fock–Roothaan
IODK infinite-order Douglas–Kroll–Hess
IOTC infinite-order two-component
LUT local unitary tranformation
MAE mean absolute error
MaxE maximum error
MCP model core potential
MP model potential
MP2 second-order Møller–Plesset perturbation
NESC normalized elimination of the small component
NMR nuclear magnetic resonance
NR non-relativistic
OEI one-electron integral
PP pseudo-potential ppy 2–phenylpyridine
RA regular approximation
RECP relativistic effective core potential
RI resolution of identity
SAC-CI symmetry-adapted cluster configuration interaction
SCF self-consistent field
SD spin-dependent
SF spin-free
vi TEI two-electron integral
UT unitary transformation
WFT wave function theory
X2C exact two-component
vii
Chapter 1 General introduction
Quantum chemistry plays an essential role in the qualitative and/or quantitative prediction and analysis of energetics and molecular properties such as geometries, spectra, and reactivities. To date, the basic, commonly used equation of quantum chemistry has been the NR Schrodinger¨ equation. In significant cases involving relativistic effects, such as orbital contractions and splittings in heavy-element systems, intersystem crossing, and core-electron related properties, relativistic effects are mostly accounted for through corrections to the NR treatment. On the other hand, using the Dirac equation, as a basic equation that satisfies the Lorentz invariance for electron motion, can account for relativistic effects. In 2002, Barysz et al. proposed a rigorous 2c method, referred to as the IODKH method1, that treats only electronic states. Seino and Nakai proposed the LUT technique that is able to perform efficient calculations without the loss of accuracy, using the original IODKH method2–6. The construction and extension of a theory based on both IODKH and LUT methods would enable the efficient treatment of relativistic effects using basic equations.
Another efficient method for the treatment of heavy elements is the ECP method, which reduces the number of electrons that are treated explicitly7. The
ECP method replaces the effect of core orbitals on the valence electrons with a potential. Constructing this potential to include relativistic effects provides a rel- ativistic treatment in a convenient manner. In 2014, Seino et al. proposed the FCP
1 method8, which describes core potentials by utilizing explicit core-orbital infor- mation obtained through atomic calculations. However, the effects of the chemical environment are not taken into account in the FCP method for molecules. This thesis extends the LUT-IODKH and FCP methods from a theoretical perspective, and applies the extended method to the determination of relativistic effects in significant systems. In order to perform geometry optimizations and frequency calculations, the analytical energy gradient for LUT-IODKH was developed and extended to the spin-dependent method9,10. This gradient method and the LUT scheme were then implemented in the GAMESS quantum chemical package11.
Furthermore, FCP-CR, which relaxes the core orbitals that were treated as frozen orbitals by FCP, was also developed12. These methods were utilized for enthalpy of formation calculations, in the gas phase, of transition metal complexes.
This thesis includes seven chapters, in addition this general introduction chapter (Chapter 1).
Chapter 2 summarizes the theoretical background of the 2c relativistic method, as well as the IODKH, LUT, and FCP methods. Moreover, the energy and gradient expressions for the GHF method are provided.
Chapter 3 extends the IODKH method to include analytical energy gradi- ents, and gradients combined with the LUT scheme. The energy gradient of the
IODKH Hamiltonian, with respect to nuclear coordinates, is analytically derived.
Numerical assessments of equilibrium bond lengths in diatomic molecules con- taining heavy elements indicate that the accuracy of this method is close to that of 4c. When compared to electron-correlation methods such as MP2 and DFT, good agreements with equilibrium bond lengths are observed. Especially for heavy-element systems, the simultaneous consideration of both electron corre- lation and relativistic effects was found to be required for accurate calculations.
For harmonic frequency and force constant calculations, relativistic effects were
2 found to contribute to 10% of the obtained values, even for first-row transition metals. In order to treat larger-sized molecules, the LUT technique is applied to the gradient calculation. The resulting CPU time required for relativistic unitary transformations in 1d- to 3d-silver clusters was found to be similar to those of NR methods, and scaled quasi-linearly with respect to the number of atoms.
Following on from the previous chapter, Chapter 4 extends and implements the SD energy gradient. The NR and SF relativistic methods are only described within real number space because no spin operator, represented by Pauli spin matrices, are included. To treat the spin operator explicitly, the GHF method, and its energy gradient with respect to nuclear coordinates, is derived from the viewpoint of its implementation within complex number space. In addition, the IODKH and LUT-IODKH methods are extended to the SD treatment. These enable efficient, self-consistent, relativistic geometry optimizations that include
SD relativistic effects. The atomic forces obtained by the analytical SD-IODKH gradient method are enough close to the corresponding numerical gradients in atomic units. Equilibrium bond lengths of diatomic molecules containing fifth- and sixth-row elements are close to those calculated using 4c methods. A com- parison of the SF and SD relativistic effects reveals that the contributions of these effects depend on the type of orbital partitioning within the bond. In addition, the CPU time for the relativistic transformation scales quasi-linearly with respect to the number of atoms.
Chapter 5 implements the LUT-IODKH method into the GAMESS quantum chemical package and assesses its performance numerically. GAMESS includes relativistic correction functionalities based on 2c methods such as NESC, RESC,
first- to third-order DKH, and IOTC. One of them, IOTC, is equivalent to IODKH from a theoretical point of view. Based on the original code for IOTC, the LUT technique, and its gradient extension, has been implemented in order to perform
3 LUT-IOTC calculations, such as excited-state calculations, in combination with
GAMESS functionalities. In order to verify the integrity of this implementation, the total energies of heavy atoms and molecules are calculated. Several of the
IOTC program-default settings were altered to obtain numerically reliable results.
Furthermore, when combined with the DC method, which has been implemented since 2009 in GAMESS, overall computational linear-scaling is achieved, from relativistic transformations to electron correlation method such as CCSD(T).
Chapter 6 describes studies into theoretically improving FCP for describing the core orbitals of molecules. Conventional FCP only describes valence orbitals of molecules accurately. In this chapter, FCP is extended by relaxing the core orbitals that have been treated as frozen orbitals, by means of valence-orbital potentials, for the accurate descriptions of molecular core orbitals. This enables accurate AE calculations using core potentials. Numerical assessments of gold clusters, composed of up to 86 atoms, indicate that FCP-CR is faster than AE. The orbital energies and total energies of coinage metal dimers are as accurate as those obtained by AE. Furthermore, because FCP-CR possesses core-orbital wavefunc- tions, it can provide core-electron related properties. Using these wavefunctions the CIEs of the 1s orbitals of C in vinyl acetate, and of the 4f orbitals of W in tungsten complexes, are qualitatively reproduced.
Chapter 7 discusses the evaluation of the gas-phase enthalpies of formation of transition metal complexes using the methods developed in the previous chapters.
These complexes include early and late transition metals, light and heavy met- als, low- and high-valency complexes, soft and hard ligands, and organometallic complexes. Contributions of relativistic effects to the enthalpies of formation, cal- culated using NR, ECP and FCP methods, are assessed and compared. In addition, the correlation between experimental and DFT-calculated data is examined.
Finally, future prospects are discussed after the results of the individual
4 chapters are briefly summarized.
5 References
(1) Barysz, M.; Sadlej, A. J. J. Chem. Phys. 2002, 116, 2696.
(2) Seino, J.; Nakai, H. J. Chem. Phys. 2012, 136, 244102.
(3) Seino, J.; Nakai, H. J. Chem. Phys. 2012, 137, 144101.
(4) Seino, J.; Nakai, H. J. Chem. Phys. 2013, 139, 034109.
(5) Seino, J.; Nakai, H. J. Comput. Chem. Jpn. 2014, 13,1.
(6) Seino, J.; Nakai, H. Int. J. Quantum Chem. 2015, 115, 253.
(7) Schwerdtfeger, P. Chem. Phys. Chem 2011, 12, 3143.
(8) Seino, J.; Tarumi, M.; Nakai, H. Chem. Phys. Lett. 2014, 592, 341.
(9) Nakajima, Y.; Seino, J.; Nakai, H. J. Chem. Phys. 2013, 139, 244107.
(10) Nakajima, Y.; Seino, J.; Nakai, H. J. Chem. Theory Comput. 2016, 15, 68.
(11) Schmidt, M. W.; Baldridge, K. K.; Boatz, J. A.; Elbert, S. T.; Gordon, M. S.;
Jensen, J. H.; Koseki, S.; Matsunaga, N.; Nguyen, K. A.; Su, S.; Windus, T. L.;
Dupuis, M.; Montgomery Jr., J. A. J. Comput. Chem. 1993, 14, 1347.
(12) Nakajima, Y.; Seino, J.; Nakai, H. Chem. Phys. Lett. 2016, 663, 97.
6 Chapter 2 Theoretical background
Relativistic quantum chemistry is based on the Dirac Hamiltonian, which con- tains electronic and positronic states. In the 2c relativistic treatment, the electronic state (or large component thereof), which is essential for chemistry, is only treated explicitly after several Hamiltonian manipulations, so as to ESCs and apply UTs that are directly related to the accuracy and efficiency of the approximations.
Moreover, the relativistic 2c Hamiltonian includes spin operators, or Pauli spin matrices, which require a generalized treatment of SCF-calculated wavefunctions in order to achieve highly accurate results. In this chapter, brief reviews of theoret- ical backgrounds are provided and include the IODKH Hamiltonian1, for accurate relativistic transformations, and its derivatives; the LUT2 and FCP3 schemes for efficient relativistic transformations without the loss of accuracy; and the GHF method, and its derivatives, for explicit use with spin operators. In addition, methods for splitting the SF and SD terms are also presented.
2.1 Dirac Hamiltonian
Relativistic quantum chemistry is based on the Dirac equation in order to include relativistic effects. The Dirac Hamiltonian is a one-particle Hamiltonian, which is
7 expressed as follows: ⎛ ⎞ ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎜ψ ⎟ ⎜ψ ⎟ ⎜ 1⎟ ⎜ 1⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ψ ⎟ ⎜ψ ⎟ ⎜ 2⎟ = ⎜ 2⎟ h4 ⎜ ⎟ E ⎜ ⎟ (2.1) ⎜ ⎟ ⎜ ⎟ ⎜ψ ⎟ ⎜ψ ⎟ ⎜ 3⎟ ⎜ 3⎟ ⎜ ⎟ ⎜ ⎟ ⎝ψ ⎠ ⎝ψ ⎠ 4 4 with ⎛ ⎞ ⎜ ⎟ ⎜ − ⎟ ⎜ V 0 cpz c px ipy ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ 0 Vc(p + ip ) −cp ⎟ = ⎜ x y z ⎟ , h4 ⎜ ⎟ (2.2) ⎜ ⎟ ⎜ cp c p − ip V − 2c2 0 ⎟ ⎜ z x y ⎟ ⎜ ⎟ ⎝ + − − 2 ⎠ c(px ipy) cpz 0 V 2c where V is a scalar potential, px, py, and pz are momentum operators along with x-, y-, and z-axis, respectively, and c is the speed of light. Note that atomic units ψ are used hereafter. The upper two wavefunction components in Eq. (2.1), 1 and ψ ψ 2, are referred to as the large components, while the lower two components, 3 ψ and 4, are referred to as the small components. These two types of wavefunction are interpreted to be the electronic and positronic wavefunctions, respectively.
Individual wavefunction components correspond to α- and β-spin, respectively.
Furthermore, using Pauli spin matrices, the Dirac Hamiltonian can be rewritten as follows: ⎛ ⎞ ⎜ σ · p ⎟ ⎜ V12 c ⎟ = ⎜ ⎟ h4 ⎜ ⎟ (2.3) ⎝ σ · p − 2 ⎠ c V12 2c 12 with σ = σ ,σ ,σ , x y z (2.4) ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎜01⎟ σ = ⎜ ⎟ , (2.5) x ⎝⎜ ⎠⎟ 10 ⎛ ⎞ ⎜ ⎟ ⎜ − ⎟ ⎜0 i⎟ σ = ⎜ ⎟ , (2.6) y ⎝⎜ ⎠⎟ i 0
8 and ⎛ ⎞ ⎜ ⎟ ⎜ ⎟ ⎜10⎟ σ = ⎜ ⎟ , (2.7) z ⎝⎜ ⎠⎟ 0 −1
× where 1n is the n n identity matrix.
2.2 The IODKH method
As mentioned above, the Dirac equation includes both electronic and positronic information. In general chemistry however, more attention is given to the elec- tronic information, which, conceptually, leads to a decoupling of the electronic and positronic states. The Dirac Hamiltonian is block-diagonalized using a uni- tary transformation U, which is expressed as: ⎛ ⎞ ⎜ + ⎟ ⎜h 0 ⎟ = † = ⎜ 2 ⎟ , h2 U h4U ⎜ ⎟ (2.8) ⎝⎜ −⎠⎟ 0h2
+ − where h2 and h2 are the electronic and positronic Hamiltonians, respectively. Practical calculations use only the electronic Hamiltonian. In this thesis, the
IODKH method, which is one of the most accurate relativistic transformation schemes, is used. Many electronic, two-component, relativistic Hamiltonians are generally expressed by:
= + + ++ , , H2 h2 (i) g2 (i j) (2.9) i i>j
+ ++ where h2 and g2 are one- and two-particle Hamiltonians, respectively. For the + one-particle Hamiltonian, h2 , the IODKH Hamiltonian was adopted, and is given by: + = Ω† Ω, h2 G (2.10) where the second-step transformation, Ω, is defined by
Ω = + † −1/2, (12 Y Y) (2.11)
9 and the Hamiltonian transformed in Eq. (2.10), G, is expressed as
p2 = 2 + † + † + Π† σ · p σ · p Π, G p b12 Y − 12 Y M VM ( V ) (2.12) 1 ep with = + , M K(12 bpY) (2.13)
Π = α − −1 , K( b12 p Y) (2.14)
+ 1/2 ep 1 K = , (2.15) 2ep
= 1 , b + (2.16) ep 1 1/2 = + α2 2 , ep 1 p (2.17) and
α = 1/c. (2.18)
Here, p is the momentum operator and V is the sum of the nucleus-electron
= attraction operators with respect to an atom A, i.e., V A VA. The operator Y is obtained by the decoupling conditions: + = α3 − −1 σ · p σ · p epY Yep pKbVK p K V bK − − + α2 p 1Kσ · pVσ · pp 1KY − YKVK (2.19) + α4 pKbVbKpY − YKbσ · pVσ · pbK − + α3Y Kbσ · pVσ · pKp 1 − KVKbp Y.
The present thesis focuses on the relativistic effect on the one-particle interac-
++ tion. Hence, the two-particle Hamiltonian g2 is replaced as the NR Coulomb interaction: 1 g(i, j) = . (2.20) rij
10 This treatment is denoted as IODKH/C. It should be noted that relativistic two- particle interactions can be included using the Breit-Pauli approximation,4 as well as the FW5 and IODKH transformations.6
In the pragmatic evaluation of Eq. (2.10),the matrix transformation method and the RI technique proposed by Hess are adopted.7 The first step in the technique of Hess is the determination of the eigenvectors, {k},ofp2, which are expanded by the linear combination of the primitive Gaussian functions {φ},
| = |φ , k r Crk (2.21) r ffi where the coe cient Crk is determined by diagonalizing a matrix whose element φ 2 φ { } is r p s . The resulting (pseudo) eigenvectors, k , satisfy: 2 = ω δ , k p k k kk (2.22)
ω { } where k is a (pseudo) eigenvalue. The explicit matrix expressions, in k -space, for Eqs. (2.10)–(2.19) can be written as: + † = Ω Ω , h Gk k k k (2.23) 2 kk kk k ,k − / † 1 2 Ω = + , kk k 12 Y Y k (2.24) 2 2 † p = + Gkk k p b12 k Ykk k − 12 k Yk k 1 ep k ,k (2.25) † † + | | +Π | σ · p σ · p | Π , Mkk k V k Mk k kk k V k k k k ,k with = + , Mkk Kk 1 bkpkYkk (2.26) − Π = α − 1 , kk Kk bk pk Ykk (2.27)
11 and ⎡ ⎢ 1 ⎢ − = ⎢α3 | | − 1 | σ · p σ · p | Ykk ⎢ pkKkbk k V k Kk p Kk k V k bk Kk + ⎣ k ek ek − − + α2 1 | σ · p σ · p | 1 − | | pk Kk k K V k pk Kk Yk k Ykk Kk k V k Kk k + α4 | | pkKkbk k V k bk Kk pk Yk k k − | σ · p σ · p | Ykk Kk bk k V k bk Kk ) − + α3 | σ · p σ · p | 1 Ykk Kk bk k V k Kk pk k ,k − | | . Kk k V k Kk bk pk Yk k
(2.28)
Here, the approximate RI relationship is inserted between the individual opera- tors, which is expressed as
|k k|≈1. (2.29) k The matrix elements, which are the functions of p2, f p2 , such as the operators
K, b, and p operators in Eqs. (2.15)–(2.17) become: 2 = ω δ . k f (p ) k f ( k) kk (2.30)
2.3 The LUT scheme
The IODKH method achieves high accuracy compared to the 4c method although the computational cost associated with the relativistic transformation is consider- ably large. On the basis of the nature of unitary transformations and relativistic effects, approximated unitary transformations can be constructed in order to im- prove the efficiency. In this section, the concept of the LUT scheme and the
LUT-IODKH Hamiltonian are provided.
The LUT scheme arises from the localization of the relativistic effect. The
12 entire one-particle unitary transformation U can be approximated in a block- diagonal matrix form consisting of the subsystem contributions:
UTotal ≈ UA ⊕ UB ⊕ UC ··· , (2.31) where {A, B, C, ···}denotes the subsystems. Here, atomic partitioning is adopted.
The relativistic effect on the kinetic energy is dominant at each atom while the effect on the nuclear attraction is dominant not only at each atom, but also on its interactions within a small distances τ. The matrix representation of the one-particle Hamiltonian with basis function {χ} is expressed as χA + χB ≈ χA LUT χB μ h2 ν μ h2 ν ⎧ ⎪ ⎪ A + + NR B ⎪ χμ T + V + V χν (A = B) ⎪ A C ⎪ CA ⎪ ⎨⎪ (2.32) χA + + + + NR + NR χB , ≤ τ = ⎪ μ V V T V ν (A B RAB ) ⎪ A B C ⎪ CA,B ⎪ ⎪ ⎪ A NR NR B ⎪ χμ T + V χν (A B, R >τ) ⎩ C AB C
NR NR where T and VC denote the standard NR kinetic and nuclear-attraction oper- ators: 1 TNR = p21 , (2.33) 2 2 and NR = . VC VC (2.34)
+ + T and VC are the relativistic operators given by: p2 + = Ω† 2 + † Ω, T p b12 Y − 12 Y (2.35) 1 ep and + = Ω† † + Π† σ · p σ · p Π Ω. VC M VCM ( VC ) (2.36)
13 Note that the operators Y, Ω, M, and Π are determined within individual sub- systems because U contains only single-subsystem information. Hence, the com- putational cost of the relativistic transformation using the LUT scheme scales linearly with respect to the system size.
2.4 Spin-free and spin-dependent formalisms
SF and SD terms are considered separately in standard relativistic quantum chem- istry. One of the reasons for this is that the SF formalism is a straightforward extension to NR theory. This subsection explains how the SF and SD formalisms are obtained.
In order to obtain the SF formalism, the spin-related terms containing Pauli spin matrices are separated by the Dirac relationship8,
σ · pVσ · p = pV · p + iσ · (pV × p) (2.37) where i denotes the imaginary unit. Since the dot products and cross products can be evaluated in {χ}-space, the explicit expression containing these terms is given by: χ | p · p | χ = χ · + · + · χ V pxV px pyV py pzV pz (2.38) and ⎛ ⎞ ⎜ ⎟ ⎜ 0 χ (pV × p) χ ⎟ χ | σ · p × p | χ = ⎜ y ⎟ i ( V ) ⎜ ⎟ ⎝ χ p × p χ ⎠ ( V )y 0 ⎛ ⎞ (2.39) ⎜ ⎟ ⎜ χ (pV × p) χ χ (pV × p) χ ⎟ + ⎜ z x ⎟ i ⎜ ⎟ ⎝ χ p × p χ − χ p × p χ ⎠ ( V )x ( V )z with χ p × p χ = χ · − · χ , ( V )x pyV pz pzV py (2.40)
χ p × p χ = χ · − · χ , ( V )y pxV pz pzV px (2.41)
14 and χ p × p χ = χ · − · χ . ( V )z pxV py pyV px (2.42)
The momentum operator p is performed on the basis set {χ}. Partial integrations of Eqs. (2.40)–(2.42) leads to: ∂χ ∂χ ∂χ ∂χ χ (pV × p) χ = V − V , (2.43) x ∂y ∂z ∂z ∂y ∂χ ∂χ ∂χ ∂χ χ (pV × p) χ = V − V , (2.44) y ∂x ∂z ∂z ∂x and ∂χ ∂χ ∂χ ∂χ χ (pV × p) χ = V − V . (2.45) z ∂x ∂y ∂y ∂x
In these equations, x, y, and z refer to electronic coordinates of electrons. In the SF formalism, the IODKH and LUT-IODKH Hamiltonians are expressed as follows: + † = Ω Ω , h Gk k k k (2.46) 2 kk kk k ,k − / † 1 2 Ω = + , kk k 12 Y Y k (2.47) 2 2 † p = + Gkk k p b12 k Ykk k − 12 k Yk k 1 ep k ,k (2.48) † † + | | +Π | p · p | Π , Mkk k V k Mk k kk k V k k k k ,k with = + , Mkk Kk 1 bkpkYkk (2.49) − Π = α − 1 , kk Kk bk pk Ykk (2.50)
15 and ⎡ ⎢ 1 ⎢ − = ⎢α3 | | − 1 | p · p | Ykk ⎢ pkKkbk k V k Kk p Kk k V k bk Kk + ⎣ k ek ek − − + α2 1 | p · p | 1 − | | pk Kk k V k pk Kk Yk k Ykk Kk k V k Kk k + α4 | | pkKkbk k V k bk Kk pk Yk k k (2.51) − | p · p | Ykk Kk bk k V k bk Kk ) − + α3 | p · p | 1 Ykk Kk bk k V k Kk pk k ,k − | | , Kk k V k Kk bk pk Yk k and ⎧ ⎪ ⎪ A + + NR B ⎪ χμ T + V + V χν (A = B) ⎪ A C ⎪ CA ⎪ ⎨⎪ A LUT B χA + + + + NR + NR χB , ≤ τ χμ h χν = ⎪ μ V V T V ν (A B RAB ) 2 ⎪ A B C ⎪ CA,B ⎪ ⎪ ⎪ A NR NR B ⎪ χμ T + V χν (A B, R >τ) ⎩ C AB C (2.52) with + = Ω† † + Π† p · p Π Ω. VC M VCM ( VC ) (2.53)
2.5 The FCP method
The IODKH method, when combined with the LUT scheme, extends the appli- cability of relativistic quantum chemistry to the calculation of relatively large molecules without loss of accuracy. However, all-electron calculations still have problems that include large basis set requirements for heavy elements, prob- lematic convergence during SCF calculations, and the demand for high memory and disk capacities. For the further extension of relativistic quantum chemistry, decreasing the number of electrons treated explicitly offers a solution to these
16 problems. In this subsection, the FCP method, that provides a solution to the problems mentioned above, is explained.
All-electron calculations are based on the HFR equation:
F(D)C = SCE (2.54) with
Fμν(D) = Hμν + Gμν(D), (2.55) 1 Gμν(D) = Dλρ χμχρ|g|χνχλ − χμχρ|g|χλχν , (2.56) λρ 2 and occ = , Dμν 2 CμiCνi (2.57) i where F is the Fock matrix; D is the density matrix, C is the orbital coefficient matrix, S is the overlap matrix, E is the orbital energy matrix, H is the one-electron integral, G is the two-electron integral, g is the two-electron Hamiltonian, χ are basis functions, while “occ” denotes occupied orbitals, and C is the complex conjugate of C.
The density matrix is divided into core and valence contributions as follows:
D = DC + DV (2.58) with core C = Dμν 2 CμcCνc (2.59) c and valence V = . Dμν 2 CμvCνv (2.60) v In the case of basis sets that segmented and contracted, the contributions from core and valence AOs to the core-density matrix, are large (μ and ν in core AOs) and small (μ or ν in valence AOs), respectively. The core AOs contribute moderately
17 to the valence-density matrix. In the case of generally contracted basis sets, the core and valence density matrices are negligibly affected by valence and core
AOs, respectively, since the individual core and valence regions are sufficiently described by core and valence AOs. When generally contracted basis sets are adopted, the entire density matrix can be approximately represented in a block- diagonalized form, as follows:
D ≈ dC ⊕ dV, (2.61) where dC and dV are the core- and valence-contributing density matrices, respec- tively, which include only core and valence AO indices. Since the individual atomic core-density matrices contribute mainly to the molecular core-density ma- trix, the introduction of the direct sum approximation to the core-density matrix sufficiently reproduces the original molecular core-density matrix:
C ≈ ˜ C = ˜ C ⊕ ˜ C ⊕··· D D DA DB (2.62) and C ≈ ˜ C = ˜ C ⊕ ˜ C ⊕···, d d dA dB (2.63)
C where {A, B, ···}are individual atoms. The atomic core-density matrices, D˜ and C d˜ , correspond to the core-density matrices obtained from the individual atomic calculations, which are frozen for molecular calculations.
These two approximations lead to the Huzinaga-Cantu equation9, which is expressed as:
fV(dV)cV = sVcVεV (2.64) with V V V V C C = μν + − + − + μν, , fμν(d ) H 2Jμν Kμν 2Jμν,A Kμν,A P A (2.65) A A
18 valence V V V V V V V 1 V V V V 2Jμν − Kμν = dλρ χμ χρ g χν χλ − χμ χρ g χλ χν , (2.66) λ,ρ 2 core C C C V C V C 1 V C C V 2Jμν, − Kμν, = d˜λρ χμ χρ g χν χλ − χμ χρ g χλ χν , (2.67) A A A 2 λ,ρ∈A and
core = − ε χV ϕ ϕ χV . Pμν,A 2 c,A μ c,A c,A ν (2.68) c∈A ϕ ε C C Here, c,A is the core orbital, c,A is the core orbital energy, J and K are Coulomb and exchange terms for the core region obtained by atomic all-electron calcula- tions, JV and KV are Coulomb and exchange terms for the valence region, and P represents projection/shift terms. Finally, the total, core, and valence energies are expressed as follows:
ETotal = EC + EV (2.69) with C − C + C 1 nA ZB nB nB ZA EC = tr DC H + F DC + (2.70) 2 RAB A