1
|
Hoyer CE, Liao C, Shumilov KD, Zhang T, Li X. State Interaction for Relativistic Four-Component Methods: Choose the Right Zeroth-Order Hamiltonian for Late-Row Elements. J Chem Theory Comput 2024. [PMID: 39257190 DOI: 10.1021/acs.jctc.4c00797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
We present several schemes based on the spin-separation of the Dirac-Coulomb-Breit Hamiltonian for the perturbative treatment of relativistic four-component Hamiltonians within the state interaction framework. While state interaction approaches traditionally use zeroth-order scalar-relativistic states, we develop augmented zeroth-order Hamiltonians with increasing accuracy and investigate convergence to the variational limit as a function of the choice of zeroth-order Hamiltonian. The state interaction schemes developed in this work are benchmarked using ground-state fine-structure splitting of late-row atoms and diatomic hyrides. Although the scalar-relativistic zeroth-order Hamiltonian exhibits significant errors in ground-state fine-structure splitting, the predictive accuracy can be improved by augmenting the zeroth-order Hamiltonian with one- and two-electron vector-relativistic operators (e.g., spin-orbit, spin-spin, orbit-orbit). This work lays the theoretical foundation for the development of low-scaling, high-accuracy perturbative relativistic methods suitable for late-row elements.
Collapse
Affiliation(s)
- Chad E Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Can Liao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Kirill D Shumilov
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Tianyuan Zhang
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
2
|
Majumder R, Sokolov AY. Consistent Second-Order Treatment of Spin-Orbit Coupling and Dynamic Correlation in Quasidegenerate N-Electron Valence Perturbation Theory. J Chem Theory Comput 2024; 20:4676-4688. [PMID: 38795071 DOI: 10.1021/acs.jctc.4c00458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2024]
Abstract
We present a formulation and implementation of second-order quasidegenerate N-electron valence perturbation theory (QDNEVPT2) that provides a balanced and accurate description of spin-orbit coupling and dynamic correlation effects in multiconfigurational electronic states. In our approach, the energies and wave functions of electronic states are computed by treating electron repulsion and spin-orbit coupling operators as equal perturbations to the nonrelativistic complete active-space wave functions, and their contributions are incorporated fully up to the second order. The spin-orbit effects are described using the Breit-Pauli (BP) or exact two-component Douglas-Kroll-Hess (DKH) Hamiltonians within spin-orbit mean-field approximation. The resulting second-order methods (BP2- and DKH2-QDNEVPT2) are capable of treating spin-orbit coupling effects in nearly degenerate electronic states by diagonalizing an effective Hamiltonian expanded in a compact non-relativistic basis. For a variety of atoms and small molecules across the entire periodic table, we demonstrate that DKH2-QDNEVPT2 is competitive in accuracy with variational two-component relativistic theories. BP2-QDNEVPT2 shows high accuracy for the second- and third-period elements, but its performance deteriorates for heavier atoms and molecules. We also consider the first-order spin-orbit QDNEVPT2 approximations (BP1- and DKH1-QDNEVPT2), among which DKH1-QDNEVPT2 is reliable but less accurate than DKH2-QDNEVPT2. Both DKH1- and DKH2-QDNEVPT2 hold promise as efficient and accurate electronic structure methods for treating electron correlation and spin-orbit coupling in a variety of applications.
Collapse
Affiliation(s)
- Rajat Majumder
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
3
|
Liao C, Hoyer CE, Banerjee Ghosh R, Jenkins AJ, Knecht S, Frisch MJ, Li X. Comparison of Variational and Perturbative Spin-Orbit Coupling within Two-Component CASSCF. J Phys Chem A 2024. [PMID: 38489510 DOI: 10.1021/acs.jpca.3c08031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2024]
Abstract
The modeling of spin-orbit coupling (SOC) remains a challenge in computational chemistry due to the high computational cost. With the rising popularity of spin-driven processes and f-block metals in chemistry and materials science, it is incumbent on the community to develop accurate multiconfigurational SOC methods that scale to large systems and understand the limits of different treatments of SOC. Herein, we introduce an implementation of perturbative SOC in scalar-relativistic two-component CASSCF (srX2C-CASSCF-SO). Perspectives on the limitations and accuracy of srX2C-CASSCF-SO are presented via benchmark calculations.
Collapse
Affiliation(s)
- Can Liao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Chad E Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Rahoul Banerjee Ghosh
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Stefan Knecht
- Algorithmiq Ltd, Kanavakatu 3C, FI-00160 Helsinki, Finland
- ETH Zürich, Department of Chemistry and Applied Life Sciences, Vladimir-Prelog-Weg 1-5/10, CH-8093 Zürich, Switzerland
| | - Michael J Frisch
- Gaussian Inc., 340 Quinnipiac Street, Bldg 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
4
|
Zhang C, Zheng X, Liu J, Asthana A, Cheng L. Analytic gradients for relativistic exact-two-component equation-of-motion coupled-cluster singles and doubles method. J Chem Phys 2023; 159:244113. [PMID: 38153147 DOI: 10.1063/5.0175041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 11/20/2023] [Indexed: 12/29/2023] Open
Abstract
A first implementation of analytic gradients for spinor-based relativistic equation-of-motion coupled-cluster singles and doubles method using an exact two-component Hamiltonian augmented with atomic mean-field spin-orbit integrals is reported. To demonstrate its applicability, we present calculations of equilibrium structures and harmonic vibrational frequencies for the electronic ground and excited states of the radium mono-amide molecule (RaNH2) and the radium mono-methoxide molecule (RaOCH3). Spin-orbit coupling is shown to quench Jahn-Teller effects in the first excited state of RaOCH3, resulting in a C3v equilibrium structure. The calculations also show that the radium atoms in these molecules serve as efficient optical cycling centers.
Collapse
Affiliation(s)
- Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Ayush Asthana
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
5
|
Liao C, Kasper JM, Jenkins AJ, Yang P, Batista ER, Frisch MJ, Li X. State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives. JACS AU 2023; 3:358-367. [PMID: 36873704 PMCID: PMC9975852 DOI: 10.1021/jacsau.2c00659] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/04/2022] [Revised: 01/17/2023] [Accepted: 01/18/2023] [Indexed: 06/18/2023]
Abstract
Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au25(SR)18 - is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.
Collapse
Affiliation(s)
- Can Liao
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| | - Joseph M. Kasper
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Andrew J. Jenkins
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| | - Ping Yang
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Enrique R. Batista
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Michael J. Frisch
- Gaussian
Inc., 340 Quinnipiac Street, Bldg 40, Wallingford, Connecticut06492, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| |
Collapse
|
6
|
Majumder R, Sokolov AY. Simulating Spin-Orbit Coupling with Quasidegenerate N-Electron Valence Perturbation Theory. J Phys Chem A 2023; 127:546-559. [PMID: 36599072 DOI: 10.1021/acs.jpca.2c07952] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We present the first implementation of spin-orbit coupling effects in fully internally contracted second-order quasidegenerate N-electron valence perturbation theory (SO-QDNEVPT2). The SO-QDNEVPT2 approach enables the computations of ground- and excited-state energies and oscillator strengths combining the description of static electron correlation with an efficient treatment of dynamic correlation and spin-orbit coupling. In addition to SO-QDNEVPT2 with the full description of one- and two-body spin-orbit interactions at the level of two-component Breit-Pauli Hamiltonian, our implementation also features a simplified approach that takes advantage of spin-orbit mean-field approximation (SOMF-QDNEVPT2). The accuracy of these methods is tested for the group 14 and 16 hydrides, 3d and 4d transition metal ions, and two actinide dioxides (neptunyl and plutonyl dications). The zero-field splittings of group 14 and 16 molecules computed using SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the available experimental data. For the 3d transition metal ions, the SO-QDNEVPT2 method is significantly more accurate than SOMF-QDNEVPT2, while no substantial difference in the performance of two methods is observed for the 4d ions. Finally, we demonstrate that for the actinide dioxides the results of SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the data from previous theoretical studies of these systems. Overall, our results demonstrate that SO-QDNEVPT2 and SOMF-QDNEVPT2 are promising multireference methods for treating spin-orbit coupling with a relatively low computational cost.
Collapse
Affiliation(s)
- Rajat Majumder
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| |
Collapse
|
7
|
Meitei OR, Mayhall NJ. Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules. J Chem Theory Comput 2021; 17:2906-2916. [PMID: 33861603 DOI: 10.1021/acs.jctc.1c00121] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.
Collapse
Affiliation(s)
- Oinam Romesh Meitei
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
8
|
Rask AE, Zimmerman PM. Toward Full Configuration Interaction for Transition-Metal Complexes. J Phys Chem A 2021; 125:1598-1609. [DOI: 10.1021/acs.jpca.0c07624] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Alan E. Rask
- Department of Chemistry, University of Michigan, 930N. University Avenue, Ann Arbor 48109, Michigan, United States
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, 930N. University Avenue, Ann Arbor 48109, Michigan, United States
| |
Collapse
|
9
|
Carreras A, Jiang H, Pokhilko P, Krylov AI, Zimmerman PM, Casanova D. Calculation of spin–orbit couplings using RASCI spinless one-particle density matrices: Theory and applications. J Chem Phys 2020; 153:214107. [DOI: 10.1063/5.0029146] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Abel Carreras
- Donostia International Physics Center (DIPC), Manuel de Lardizabal Pasalekua 4, 20018 Donostia, Euskadi, Spain
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Pavel Pokhilko
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - David Casanova
- Donostia International Physics Center (DIPC), Manuel de Lardizabal Pasalekua 4, 20018 Donostia, Euskadi, Spain
| |
Collapse
|
10
|
Lin HH, Maschio L, Kats D, Usvyat D, Heine T. Fragment-Based Restricted Active Space Configuration Interaction with Second-Order Corrections Embedded in Periodic Hartree–Fock Wave Function. J Chem Theory Comput 2020; 16:7100-7108. [DOI: 10.1021/acs.jctc.0c00576] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Hung-Hsuan Lin
- Theoretische Chemie, Technische Universität Dresden, Dresden, Germany
| | - Lorenzo Maschio
- Dipartimento di Chimica, Università di Torino, Torino, Italy
| | - Daniel Kats
- Max-Planck-Institut für Festkörperforschung, Stuttgart, Germany
| | - Denis Usvyat
- Institut für Chemie, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Thomas Heine
- Theoretische Chemie, Technische Universität Dresden, Dresden, Germany
| |
Collapse
|