1
|
Uhlířová T, Cianchino D, Nottoli T, Lipparini F, Gauss J. Cholesky Decomposition in Spin-Free Dirac-Coulomb Coupled-Cluster Calculations. J Phys Chem A 2024; 128:8292-8303. [PMID: 39268870 DOI: 10.1021/acs.jpca.4c04353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
We present an implementation for the use of Cholesky decomposition (CD) of two-electron integrals within the spin-free Dirac-Coulomb (SFDC) scheme that enables to perform high-accuracy coupled-cluster (CC) calculations at costs almost comparable to those of their nonrelativistic counterparts. While for nonrelativistic CC calculations, atomic-orbital (AO)-based algorithms, due to their significantly reduced disk-space requirements, are the key to efficient large-scale computations, such algorithms are less advantageous in the SFDC case due to their increased computational cost in that case. Here, molecular-orbital (MO)-based algorithms exploiting the CD of the two-electron integrals allow us to reduce disk-space requirements and lead to computational cost in the CC step that is more or less the same as in the nonrelativistic case. The only remaining overhead in a CD-SFDC-CC calculation is due to the need to compute additional two-electron integrals, the somewhat higher cost of the Hartree-Fock calculation in the SFDC case, and additional cost in the transformation of the Cholesky vectors from the AO to the MO representation. However, these additional costs typically amount to less than 5-15% of the total wall time and are thus acceptable. We illustrate the efficiency of our CD scheme for SFDC-CC calculations on a series of illustrative calculations for the X(CO)4 molecules with X = Ni, Pd, Pt.
Collapse
Affiliation(s)
- Tereza Uhlířová
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz D-55128, Germany
| | - Davide Cianchino
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via G. Moruzzi 13, Pisa I-56124, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via G. Moruzzi 13, Pisa I-56124, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Universitá di Pisa, Via G. Moruzzi 13, Pisa I-56124, Italy
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz D-55128, Germany
| |
Collapse
|
2
|
Desmarais JK, Boccuni A, Flament JP, Kirtman B, Erba A. Perturbation Theory Treatment of Spin-Orbit Coupling. III: Coupled Perturbed Method for Solids. J Chem Theory Comput 2023; 19:1853-1863. [PMID: 36917759 DOI: 10.1021/acs.jctc.3c00088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A previously proposed noncanonical coupled-perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment for spin-orbit coupling is here generalized to infinite periodic systems. The scalar-relativistic periodic KS-DFT/HF solution, obtained with a relativistic effective core potential, is taken as the zeroth-order approximation. Explicit expressions are given for the total energy through third-order, which satisfy the 2N + 1 rule (i.e., requiring only the first-order perturbed wave function for determining the energy through third-order). Expressions for additional second-order corrections to the perturbed wave function (as well as related one-electron properties) are worked out at the uncoupled-perturbed level of theory. The approach is implemented in the Crystal program and validated with calculations of the total energy, electronic band structure, and density variables of spin-current DFT on the tungsten dichalcogenide hexagonal bilayer series (i.e., WSe2, WTe2, WPo2, WLv2), including 6p and 7p elements as a stress test. The computed properties through second- or third-order match well with those from reference two-component self-consistent field (2c-SCF) calculations. For total energies, E(3) was found to consistently improve the agreement against the 2c-SCF reference values. For electronic band structures, visible differences w.r.t. 2c-SCF remained through second-order in only the single-most difficult case of WLv2. As for density variables of spin-current DFT, the perturbed electron density, being vanishing in first-order, is the most challenging for the perturbation theory approach. The visible differences in the electron densities are, however, largest close to the core region of atoms and smaller in the valence region. Perturbed spin-current densities, on the other hand, are well reproduced in all tested cases.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alberto Boccuni
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523 ─ PhLAM ─ Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
3
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling, Part I: Double Perturbation Theory Based on a Single-Reference Initial Approximation. J Chem Theory Comput 2021; 17:4697-4711. [PMID: 34288690 DOI: 10.1021/acs.jctc.1c00343] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We develop a perturbation theory for solving the many-body Dirac equation within a given relativistic effective-core potential approximation. Starting from a scalar-relativistic unrestricted Hartree-Fock (SR UHF) solution, we carry out a double perturbation expansion in terms of spin-orbit coupling (SOC) and the electron fluctuation potential. Computationally convenient energy expressions are derived through fourth order in SOC, second order in the electron fluctuation potential, and a total of third order in the coupling between the two. Illustrative calculations on the halogen series of neutral and singly positive diatomic molecules show that the perturbation expansion is well-converged by taking into account only the leading (nonvanishing) term at each order of the electron fluctuation potential. Our perturbation theory approach provides a computationally attractive alternative to a two-component self-consistent field treatment of SOC. In addition, it includes coupling with the fluctuation potential through third order and can be extended (in principle) to multireference calculations, when necessary for both closed- and open-shell cases, using quasi-degenerate perturbation theory.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,IPREM, E2S UPPA, CNRS, Université de Pau et des Pays de l'Adour, 64053 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, CNRS, Université de Lille, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
4
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method. J Chem Theory Comput 2021; 17:4712-4732. [PMID: 34286577 DOI: 10.1021/acs.jctc.1c00460] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At2-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,Université de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, 64000 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
5
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
6
|
Stopkowicz S, Gauss J. A one-electron variant of direct perturbation theory for the treatment of scalar-relativistic effects. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1536812] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Stella Stopkowicz
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| |
Collapse
|
7
|
Cheng L, Wang F, Stanton JF, Gauss J. Perturbative treatment of spin-orbit-coupling within spin-free exact two-component theory using equation-of-motion coupled-cluster methods. J Chem Phys 2018; 148:044108. [DOI: 10.1063/1.5012041] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, China
| | - John F. Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| |
Collapse
|
8
|
Epifanovsky E, Klein K, Stopkowicz S, Gauss J, Krylov AI. Spin-orbit couplings within the equation-of-motion coupled-cluster framework: Theory, implementation, and benchmark calculations. J Chem Phys 2015; 143:064102. [DOI: 10.1063/1.4927785] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Evgeny Epifanovsky
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Kerstin Klein
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Stella Stopkowicz
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, University of Oslo, N-0315 Oslo, Norway
| | - Jürgen Gauss
- Institut für Physikalische Chemie, Universität Mainz, D-55099 Mainz, Germany
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
9
|
Cheng L, Stopkowicz S, Gauss J. Spin-free Dirac-Coulomb calculations augmented with a perturbative treatment of spin-orbit effects at the Hartree-Fock level. J Chem Phys 2013; 139:214114. [DOI: 10.1063/1.4832739] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
10
|
Kelley MS, Shiozaki T. Large-scale Dirac–Fock–Breit method using density fitting and 2-spinor basis functions. J Chem Phys 2013; 138:204113. [DOI: 10.1063/1.4807612] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
11
|
Mück LA, Gauss J. Communication: Spin-orbit splittings in degenerate open-shell states via Mukherjee's multireference coupled-cluster theory: A measure for the coupling contribution. J Chem Phys 2012; 136:111103. [DOI: 10.1063/1.3694132] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Helgaker T, Coriani S, Jørgensen P, Kristensen K, Olsen J, Ruud K. Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations. Chem Rev 2012; 112:543-631. [DOI: 10.1021/cr2002239] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trygve Helgaker
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy
| | - Poul Jørgensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|
13
|
Cheng L, Gauss J. Analytic second derivatives for the spin-free exact two-component theory. J Chem Phys 2011; 135:244104. [DOI: 10.1063/1.3667202] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
14
|
Schwalbach W, Stopkowicz S, Cheng L, Gauss J. Direct perturbation theory in terms of energy derivatives: Scalar-relativistic treatment up to sixth order. J Chem Phys 2011; 135:194114. [DOI: 10.1063/1.3659316] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
|
15
|
Saue T. Relativistic Hamiltonians for chemistry: a primer. Chemphyschem 2011; 12:3077-94. [PMID: 22076930 DOI: 10.1002/cphc.201100682] [Citation(s) in RCA: 312] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2011] [Indexed: 11/06/2022]
Affiliation(s)
- Trond Saue
- Laboratoire de Chimie et Physique Quantique (UMR 5626), CNRS/Université de Toulouse 3 (Paul Sabatier), 118 route de Narbonne, 31062 Toulouse, France.
| |
Collapse
|
16
|
Cheng L, Gauss J. Analytic energy gradients for the spin-free exact two-component theory using an exact block diagonalization for the one-electron Dirac Hamiltonian. J Chem Phys 2011; 135:084114. [DOI: 10.1063/1.3624397] [Citation(s) in RCA: 124] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
17
|
Cheng L, Gauss J. Analytical evaluation of first-order electrical properties based on the spin-free Dirac-Coulomb Hamiltonian. J Chem Phys 2011; 134:244112. [DOI: 10.1063/1.3601056] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
18
|
Stopkowicz S, Gauss J. Fourth-order relativistic corrections to electrical first-order properties using direct perturbation theory. J Chem Phys 2011; 134:204106. [DOI: 10.1063/1.3587633] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
|
19
|
Puzzarini C, Cazzoli G, López JC, Alonso JL, Baldacci A, Baldan A, Stopkowicz S, Cheng L, Gauss J. Spectroscopic investigation of fluoroiodomethane, CH2FI: Fourier-transform microwave and millimeter-/submillimeter-wave spectroscopy and quantum-chemical calculations. J Chem Phys 2011; 134:174312. [DOI: 10.1063/1.3583498] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|