1
|
Liu W. Comment on "Theoretical examination of QED Hamiltonian in relativistic molecular orbital theory" [J. Chem. Phys. 159, 054105 (2023)]. J Chem Phys 2024; 160:187101. [PMID: 38739035 DOI: 10.1063/5.0174011] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Accepted: 01/02/2024] [Indexed: 05/14/2024] Open
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
2
|
Liu W. Unified construction of relativistic Hamiltonians. J Chem Phys 2024; 160:084111. [PMID: 38415836 DOI: 10.1063/5.0188794] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 02/09/2024] [Indexed: 02/29/2024] Open
Abstract
It is shown that the four-component (4C), quasi-four-component (Q4C), and exact two-component (X2C) relativistic Hartree-Fock equations can be implemented in a unified manner by making use of the atomic nature of the small components of molecular 4-spinors. A model density matrix approximation can first be invoked for the small-component charge/current density functions, which gives rise to a static, pre-molecular mean field to be combined with the one-electron term. As a result, only the nonrelativistic-like two-electron term of the 4C/Q4C/X2C Fock matrix needs to be updated during the iterations. A "one-center small-component" approximation can then be invoked in the evaluation of relativistic integrals, that is, all atom-centered small-component basis functions are regarded as extremely localized near the position of the atom to which they belong such that they have vanishing overlaps with all small- or large-component functions centered at other nuclei. Under these approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians share precisely the same structure and accuracy. Beyond these is the effective quantum electrodynamics Hamiltonian that can be constructed in the same way. Such approximations lead to errors that are orders of magnitude smaller than other sources of errors (e.g., truncation errors in the one- and many-particle bases as well as uncertainties of experimental measurements) and are, hence, safe to use for whatever purposes. The quaternion forms of the 4C, Q4C, and X2C equations are also presented in the most general way, based on which the corresponding Kramers-restricted open-shell variants are formulated for "high-spin" open-shell systems.
Collapse
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
3
|
Liao C, Lambros E, Sun Q, Dyall KG, Li X. Exploring Locality in Molecular Dirac-Coulomb-Breit Calculations: A Perspective. J Chem Theory Comput 2023; 19:9009-9017. [PMID: 38090757 DOI: 10.1021/acs.jctc.3c01012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
The Dirac-Coulomb-Breit (DCB) operator is widely recognized for its ability to accurately capture relativistic effects and spin-physics in molecular calculations. However, due to its high computational cost, there is a need to develop low-scaling approximations without compromising accuracy. To tackle this challenge, it becomes essential to gain a deeper understanding of the DCB operator's behavior. This work aims to explore local integral approximations, shedding light on the locality of the parts of the charge-current distribution due to the small component. In particular, we propose an atomic Breit approximation that leverages an analysis of the behavior observed in a series of gold chains. Through benchmark studies of metal complexes, we evaluated the accuracy and performance of the proposed atomic Breit approximation. This work provides a comprehensive understanding of the behavior of the charge-current distribution in terms of its contributions from its AO basis constituents, facilitating the development of low-scaling methods that strike a balance between computational efficiency and accuracy.
Collapse
Affiliation(s)
- Can Liao
- Department of Chemistry, University of Washington, Seattle, Washington, 98195 United States
| | - Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington, 98195 United States
| | - Qiming Sun
- AxiomQuant Investment Management LLC, Shanghai, 200120 China
| | | | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington, 98195 United States
| |
Collapse
|
4
|
Audinet T, Toulouse J. Effective quantum electrodynamics: One-dimensional model of the relativistic hydrogen-like atom. J Chem Phys 2023; 158:244108. [PMID: 37352423 DOI: 10.1063/5.0152956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Accepted: 06/06/2023] [Indexed: 06/25/2023] Open
Abstract
We consider a one-dimensional effective quantum electrodynamics (QED) model of the relativistic hydrogen-like atom using delta-potential interactions. We discuss the general exact theory and the Hartree-Fock approximation. The present one-dimensional effective QED model shares the essential physical feature of the three-dimensional theory: the nuclear charge polarizes the vacuum state (creation of electron-positron pairs), which results in a QED Lamb-type shift of the bound-state energy. Yet, this 1D effective QED model eliminates some of the most serious technical difficulties of the three-dimensional theory coming from renormalization. We show how to calculate the vacuum-polarization density at zeroth order in the two-particle interaction and the QED Lamb-type shift of the bound-state energy at first order in the two-particle interaction. The present work may be considered a step toward the development of a quantum-chemistry effective QED theory of atoms and molecules.
Collapse
Affiliation(s)
- Timothée Audinet
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| | - Julien Toulouse
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
- Institut Universitaire de France, F-75005 Paris, France
| |
Collapse
|
5
|
Hoyer CE, Lu L, Hu H, Shumilov KD, Sun S, Knecht S, Li X. Correlated Dirac-Coulomb-Breit multiconfigurational self-consistent-field methods. J Chem Phys 2023; 158:044101. [PMID: 36725503 DOI: 10.1063/5.0133741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The fully correlated frequency-independent Dirac-Coulomb-Breit Hamiltonian provides the most accurate description of electron-electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac-Coulomb-Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac-Coulomb-Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. We also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.
Collapse
Affiliation(s)
- Chad E Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Hang Hu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kirill D Shumilov
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Shichao Sun
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Stefan Knecht
- Algorithmiq Ltd., Kanavakatu 3C, FI-00160 Helsinki, Finland
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
6
|
Liu W. Perspective: Simultaneous treatment of relativity, correlation, and
QED. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science Shandong University Qingdao Shandong China
| |
Collapse
|
7
|
Sun S, Ehrman JN, Sun Q, Li X. Efficient Evaluation of the Breit Operator in the Pauli Spinor Basis. J Chem Phys 2022; 157:064112. [DOI: 10.1063/5.0098828] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The frequency-independent Coulomb-Breit operator gives rise to the most accurate treatment of two-electron interaction in the non-quantum-electrodynamics regime. The Breit interaction in the Coulomb gauge consists of magnetic and gauge contributions. The high computational cost of the gauge term limits the application of the Breit interaction in relativistic molecular calculations. In this work, we apply the Pauli component integral-density matrix contraction scheme for gauge interaction with a maximum spin- and component separation scheme. We also present two different computational algorithms for evaluating gauge integrals. One is the generalized Obara-Saika algorithm, where the Laplace transformation is used to transform the gauge operator into Gaussian functions and the Obara-Saika recursion is used for reducing the angular momentum. The other algorithm is the second derivative of inverse Coulomb interaction evaluated with Rys-quadrature. This work improves the efficiency of performing Dirac-Hartree-Fock with variational treatment of Breit interaction for molecular systems. We use this formalism to examine relativistic trends in the periodic table, and analyze the relativistic two-electron interaction contributions in heavy-element complexes.
Collapse
Affiliation(s)
- Shichao Sun
- Chemistry, University of Washington, United States of America
| | | | - Qiming Sun
- Anxian Investment Management Co. Ltd, China
| | - Xiaosong Li
- Department of Chemistry, University of Washington, United States of America
| |
Collapse
|
8
|
Zhang N, Xiao Y, Liu W. SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:224007. [PMID: 35287124 DOI: 10.1088/1361-648x/ac5db4] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 03/14/2022] [Indexed: 06/14/2023]
Abstract
The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.
Collapse
Affiliation(s)
- Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, People's Republic of China
| |
Collapse
|
9
|
Ferenc D, Jeszenszki P, Mátyus E. On the Breit interaction in an explicitly correlated variational Dirac–Coulomb framework. J Chem Phys 2022; 156:084110. [DOI: 10.1063/5.0075097] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Breit interaction is implemented in the no-pair variational Dirac–Coulomb (DC) framework using an explicitly correlated Gaussian basis reported in the previous paper [P. Jeszenszki, D. Ferenc, and E. Mátyus, J. Chem. Phys. 156, 084111 (2022)]. Both a perturbative and a fully variational inclusion of the Breit term are considered. The no-pair DC plus perturbative Breit and the no-pair DC–Breit energies are compared with perturbation theory results including the Breit–Pauli Hamiltonian and leading-order non-radiative quantum electrodynamics corrections for low Z values. Possible reasons for the observed deviations are discussed.
Collapse
Affiliation(s)
- Dávid Ferenc
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest, H-1117, Hungary
| | - Péter Jeszenszki
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest, H-1117, Hungary
| | - Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest, H-1117, Hungary
| |
Collapse
|
10
|
Jeszenszki P, Ferenc D, Mátyus E. Variational Dirac–Coulomb explicitly correlated computations for atoms and molecules. J Chem Phys 2022; 156:084111. [DOI: 10.1063/5.0075096] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Dirac–Coulomb equation with positive-energy projection is solved using explicitly correlated Gaussian functions. The algorithm and computational procedure aims for a parts-per-billion convergence of the energy to provide a starting point for further comparison and further developments in relation with high-resolution atomic and molecular spectroscopy. Besides a detailed discussion of the implementation of the fundamental spinor structure, permutation, and point-group symmetries, various options for the positive-energy projection procedure are presented. The no-pair Dirac–Coulomb energy converged to a parts-per-billion precision is compared with perturbative results for atomic and molecular systems with small nuclear charge numbers. Paper II [D. Ferenc, P. Jeszenszki, and E. Mátyus, J. Chem. Phys. 156, 084110 (2022).] describes the implementation of the Breit interaction in this framework.
Collapse
Affiliation(s)
- Péter Jeszenszki
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| | - Dávid Ferenc
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| | - Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| |
Collapse
|
11
|
Song Y, Guo Y, Lei Y, Zhang N, Liu W. The Static-Dynamic-Static Family of Methods for Strongly Correlated Electrons: Methodology and Benchmarking. Top Curr Chem (Cham) 2021; 379:43. [PMID: 34724123 DOI: 10.1007/s41061-021-00351-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Accepted: 09/15/2021] [Indexed: 11/28/2022]
Abstract
A series of methods (SDSCI, SDSPT2, iCI, iCIPT2, iCISCF(2), iVI, and iCAS) is introduced to accurately describe strongly correlated systems of electrons. Born from the (restricted) static-dynamic-static (SDS) framework for designing many-electron wave functions, SDSCI is a minimal multireference (MR) configuration interaction (CI) approach that constructs and diagonalizes a [Formula: see text] matrix for [Formula: see text] states, regardless of the numbers of orbitals and electrons to be correlated. If the full molecular Hamiltonian H in the QHQ block (which describes couplings between functions of the first-order interaction space Q) of the SDSCI CI matrix is replaced with a zeroth-order Hamiltonian [Formula: see text] before the diagonalization is taken, we obtain SDSPT2, a CI-like second-order perturbation theory (PT2). Unlike most variants of MRPT2, SDSPT2 treats single and multiple states in the same way and is particularly advantageous in the presence of near degeneracy. On the other hand, if the SDSCI procedure is repeated until convergence, we will have iterative CI (iCI), which can converge quickly from the above to the exact solutions (full CI) even when starting with a poor guess. When further combined with the selection of important configurations followed by a PT2 treatment of dynamic correlation, iCI becomes iCIPT2, which is a near-exact theory for medium-sized systems. The microiterations of iCI for relaxing the coefficients of contracted many-electron functions can be generalized to an iterative vector interaction (iVI) approach for finding exterior or interior roots of a given matrix, in which the dimension of the search subspace is fixed by either the number of target roots or the user-specified energy window. Naturally, iCIPT2 can be employed as the active space solver of the complete active space (CAS) self-consistent field, leading to iCISCF(2), which can further be combined with iCAS for automated selection of active orbitals and assurance of the same CAS for all states and all geometries. The methods are calibrated by taking the Thiel set of benchmark systems as examples. Results for the corresponding cations, a new set of benchmark systems, are also reported.
Collapse
Affiliation(s)
- Yangyang Song
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, 266237, Shandong, China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, 266237, Shandong, China
| | - Yibo Lei
- Key Laboratory of Synthetic and Natural Functional Molecule of the Ministry of Education, College of Chemistry and Materials Science, Shaanxi key Laboratory of Physico-Inorganic Chemistry, Northwest University, Xi'an, 710127, Shaanxi, China
| | - Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, 266237, Shandong, China.
| |
Collapse
|
12
|
Paquier J, Giner E, Toulouse J. Relativistic short-range exchange energy functionals beyond the local-density approximation. J Chem Phys 2020; 152:214106. [DOI: 10.1063/5.0004926] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Julien Paquier
- Laboratoire de Chimie Théorique (LCT), Sorbonne Université and CNRS, F-75005 Paris, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique (LCT), Sorbonne Université and CNRS, F-75005 Paris, France
| | - Julien Toulouse
- Laboratoire de Chimie Théorique (LCT), Sorbonne Université and CNRS, F-75005 Paris, France
- Institut Universitaire de France, F-75005 Paris, France
| |
Collapse
|
13
|
Affiliation(s)
- Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| |
Collapse
|
14
|
Zhang Y, Suo B, Wang Z, Zhang N, Li Z, Lei Y, Zou W, Gao J, Peng D, Pu Z, Xiao Y, Sun Q, Wang F, Ma Y, Wang X, Guo Y, Liu W. BDF: A relativistic electronic structure program package. J Chem Phys 2020; 152:064113. [DOI: 10.1063/1.5143173] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yong Zhang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Bingbing Suo
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Zikuan Wang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Ning Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, People’s Republic of China
| | - Yibo Lei
- Key Laboratory of Synthetic and Natural Functional Molecule Chemistry of Ministry of Education, Shaanxi Key Laboratory of Physico-Inorganic Chemistry, College of Chemistry and Materials Science, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Wenli Zou
- Shaanxi Key Laboratory for Theoretical Physics Frontiers, Institute of Modern Physics, Northwest University, Xi’an, Shaanxi 710127, People’s Republic of China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, Hubei 430070, People’s Republic of China
| | - Daoling Peng
- College of Chemistry, South China Normal University, Guangzhou, Guangdong 510006, People’s Republic of China
| | - Zhichen Pu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing 100871, People’s Republic of China
| | - Qiming Sun
- Tencent America LLC, Palo Alto, California 94306, USA
| | - Fan Wang
- Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Institute of Atomic and Molecular Physics, Sichuan University, Chengdu, Sichuan 610065, People’s Republic of China
| | - Yongtao Ma
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Xiaopeng Wang
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, People’s Republic of China
| |
Collapse
|
15
|
De Santis M, Rampino S, Quiney HM, Belpassi L, Storchi L. Charge-Displacement Analysis via Natural Orbitals for Chemical Valence in the Four-Component Relativistic Framework. J Chem Theory Comput 2018; 14:1286-1296. [DOI: 10.1021/acs.jctc.7b01077] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Matteo De Santis
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Sergio Rampino
- Istituto di Scienze e Tecnologie Molecolari, Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Harry M. Quiney
- ARC Centre of Excellence for Advanced Molecular Imaging, School of Physics, The University of Melbourne, 3010 Victoria, Australia
| | - Leonardo Belpassi
- Istituto di Scienze e Tecnologie Molecolari, Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Consortium for Computational Molecular and Materials Sciences (CMS)2, Via Elce di Sotto, 8, 06123 Perugia, Italy
| | - Loriano Storchi
- Istituto di Scienze e Tecnologie Molecolari, Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento di Farmacia, Università degli Studi ‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| |
Collapse
|
16
|
Abstract
The foundations, formalisms, technicalities, and practicalities of relativistic time-dependent density functional theories (R-TD-DFT) for spinor excited states of molecular systems containing heavy elements are critically reviewed.
Collapse
Affiliation(s)
- Wenjian Liu
- Beijing National Center for Molecular Sciences
- Institute of Theoretical and Computational Chemistry
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
| | - Yunlong Xiao
- Beijing National Center for Molecular Sciences
- Institute of Theoretical and Computational Chemistry
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
| |
Collapse
|
17
|
Sasmal S, Talukdar K, Nayak MK, Vaval N, Pal S. Electron–nucleus scalar–pseudoscalar interaction in PbF: Z-vector study in the relativistic coupled-cluster framework. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1332396] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Sudip Sasmal
- Electronic Structure Theory Group, Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune, India
| | - Kaushik Talukdar
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| | - Malaya K. Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, India
| | - Nayana Vaval
- Electronic Structure Theory Group, Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune, India
| | - Sourav Pal
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
18
|
Niskanen J, Jänkälä K, Huttula M, Föhlisch A. QED effects in 1s and 2s single and double ionization potentials of the noble gases. J Chem Phys 2017; 146:144312. [PMID: 28411594 DOI: 10.1063/1.4979991] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present calculations on the quantum electrodynamics (QED) effects in 1s and 2s single and double ionization potentials of noble gases from Ne to Rn as perturbations on relativistic four-component Dirac-Fock wavefunctions. The most dominant effect originates from the self-energy of the core-electron that yields corrections of similar order as the transverse interaction. For 1s ionization potentials, a match within few eV against the known experimental values is obtained, and our work reveals considerable QED effects in the photoelectron binding energies across the periodic table-most strikingly even for Ne. We perform power-law fits for the corrections as a function of Z and interpolate the QED correction of ∼-0.55 eV for S1s. Due to this, the K-edge electron spectra of the third row and below need QED for a match in the absolute energy when using state-of-the-art instrumentation.
Collapse
Affiliation(s)
- J Niskanen
- Helmholtz Zentrum Berlin für Materialien und Energie, Institute for Methods and Instrumentation for Synchrotron Radiation Research, Albert-Einstein-St. 15, D-12489 Berlin, Germany
| | - K Jänkälä
- Nano and Molecular Systems Research Unit, University of Oulu, P.O. Box 3000, 90014 Oulu, Finland
| | - M Huttula
- Nano and Molecular Systems Research Unit, University of Oulu, P.O. Box 3000, 90014 Oulu, Finland
| | - A Föhlisch
- Helmholtz Zentrum Berlin für Materialien und Energie, Institute for Methods and Instrumentation for Synchrotron Radiation Research, Albert-Einstein-St. 15, D-12489 Berlin, Germany
| |
Collapse
|
19
|
Pašteka LF, Eliav E, Borschevsky A, Kaldor U, Schwerdtfeger P. Relativistic Coupled Cluster Calculations with Variational Quantum Electrodynamics Resolve the Discrepancy between Experiment and Theory Concerning the Electron Affinity and Ionization Potential of Gold. PHYSICAL REVIEW LETTERS 2017; 118:023002. [PMID: 28128629 DOI: 10.1103/physrevlett.118.023002] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2016] [Indexed: 06/06/2023]
Abstract
The first ionization potential (IP) and electron affinity (EA) of the gold atom have been determined to an unprecedented accuracy using relativistic coupled cluster calculations up to the pentuple excitation level including the Breit and QED contributions. We reach meV accuracy (with respect to the experimental values) by carefully accounting for all individual contributions beyond the standard relativistic coupled cluster approach. Thus, we are able to resolve the long-standing discrepancy between experimental and theoretical IP and EA of gold.
Collapse
Affiliation(s)
- L F Pašteka
- Centre for Theoretical Chemistry and Physics, New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0745 Auckland, New Zealand
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University, Mlynská dolina, 84104 Bratislava, Slovakia
| | - E Eliav
- School of Chemistry, Tel Aviv University, 6997801 Tel Aviv, Israel
| | - A Borschevsky
- Van Swinderen Institute for Particle Physics and Gravity, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - U Kaldor
- School of Chemistry, Tel Aviv University, 6997801 Tel Aviv, Israel
| | - P Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, New Zealand Institute for Advanced Study, Massey University Auckland, Private Bag 102904, 0745 Auckland, New Zealand
| |
Collapse
|
20
|
Almoukhalalati A, Knecht S, Jensen HJA, Dyall KG, Saue T. Electron correlation within the relativistic no-pair approximation. J Chem Phys 2016; 145:074104. [DOI: 10.1063/1.4959452] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
|
21
|
Liu W, Lindgren I. Erratum: “Going beyond ‘no-pair relativistic quantum chemistry’” [J. Chem. Phys. 139, 014108 (2013)]. J Chem Phys 2016; 144:049901. [DOI: 10.1063/1.4940777] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wenjian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Ingvar Lindgren
- Department of Physics, University of Gothenburg, Gothenburg, Sweden
| |
Collapse
|
22
|
Abstract
Abstract
Any quantum mechanical calculation on electronic structure ought to choose first an appropriate Hamiltonian H and then an Ansatz for parameterizing the wave function Ψ, from which the desired energy/property E(λ) can finally be calculated. Therefore, the very first question is: what is the most accurate many-electron Hamiltonian H? It is shown that such a Hamiltonian i.e. effective quantum electrodynamics (eQED) Hamiltonian, can be obtained naturally by incorporating properly the charge conjugation symmetry when normal ordering the second quantized fermion operators. Taking this eQED Hamiltonian as the basis, various approximate relativistic many-electron Hamiltonians can be obtained based entirely on physical arguments. All these Hamiltonians together form a complete and continuous ‘Hamiltonian ladder’, from which one can pick up the right one according to the target physics and accuracy. As for the many-electron wave function Ψ, the most intriguing questions are as follows. (i) How to do relativistic explicit correlation? (ii) How to handle strong correlation? Both general principles and practical strategies are outlined here to handle these issues. Among the electronic properties E(λ) that sample the electronic wave function nearby the nuclear region, nuclear magnetic resonance (NMR) shielding and nuclear spin-rotation (NSR) coupling constant are especially challenging: they require body-fixed molecular Hamiltonians that treat both the electrons and nuclei as relativistic quantum particles. Nevertheless, they have been formulated rigorously. In particular, a very robust ‘relativistic mapping’ between the two properties has been established, which can translate experimentally measured NSR coupling constants to very accurate absolute NMR shielding scales that otherwise cannot be obtained experimentally. Since the most general and fundamental issues pertinent to all the three components of the quantum mechanical equation HΨ = EΨ (i.e. Hamiltonian H, wave function Ψ, and energy/property E(λ)) have fully been understood, the big picture of relativistic molecular quantum mechanics can now be regarded as established.
Collapse
Affiliation(s)
- Wenjian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
23
|
Bhattacharyya S, Opalka D, Poluyanov LV, Domcke W. The (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller Hamiltonian Including Spin–Orbit Coupling for Trigonal Systems. J Phys Chem A 2014; 118:11962-70. [DOI: 10.1021/jp506793z] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Daniel Opalka
- Department
of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Leonid V. Poluyanov
- Institute
of Chemical Physics, Academy of Sciences, Chernogolovka, Moscow 142432, Russia
| | - Wolfgang Domcke
- Department
of Chemistry, Technische Universität München, D-85748 Garching, Germany
| |
Collapse
|
24
|
Li Z, Xiao Y, Liu W. On the spin separation of algebraic two-component relativistic Hamiltonians: Molecular properties. J Chem Phys 2014; 141:054111. [DOI: 10.1063/1.4891567] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Zhendong Li
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Yunlong Xiao
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People's Republic of China
| | - Wenjian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, State Key Laboratory of Rare Earth Materials Chemistry and Applications, College of Chemistry and Molecular Engineering, and Center for Computational Science and Engineering, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
25
|
Liu Y, Wu J. Structure and thermodynamic properties of relativistic electron gases. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:012141. [PMID: 25122285 DOI: 10.1103/physreve.90.012141] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/25/2014] [Indexed: 06/03/2023]
Abstract
Relativistic effect is important in many quantum systems but theoretically complicated from both fundamental and practical perspectives. Herein we introduce an efficient computational procedure to predict the structure and energetic properties of relativistic quantum systems by mapping the Pauli principle into an effective pairwise-additive potential such that the properties of relativistic nonquantum systems can be readily predicted from conventional liquid-state methods. We applied our theoretical procedure to relativistic uniform electron gases and compared the pair correlation functions with those for systems of nonrelativistic electrons. A simple analytical expression has been developed to correlate the exchange-correlation free energy of relativistic uniform electron systems.
Collapse
Affiliation(s)
- Yu Liu
- Department of Chemical and Environmental Engineering and Department of Mathematics, University of California, Riverside, California 92521, USA
| | - Jianzhong Wu
- Department of Chemical and Environmental Engineering and Department of Mathematics, University of California, Riverside, California 92521, USA
| |
Collapse
|
26
|
SEINO J, NAKAI H. Large-Scale and Highly Accurate Relativistic Quantum-Chemical Scheme:toward Establishment ofTheoretical Foundation for Element Strategy. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2014. [DOI: 10.2477/jccj.2013-0010] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
27
|
Tang LY, Tang YB, Shi TY, Mitroy J. Solution of the Dirac Coulomb equation for helium-like ions in the Poet-Temkin model. J Chem Phys 2013; 139:134112. [PMID: 24116557 DOI: 10.1063/1.4823486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The Dirac-Coulomb equation for the helium atom is studied under the restrictions of the Poet-Temkin model which replaces the 1/r12 interaction by the simplified 1/r> form. The effective reduction in the dimensionality made it possible to obtain binding energies for the singlet and triplet states in this model problem with a relative precision from 10(-8) to 10(-10). The energies for the singlet state were consistent with a previous configuration interaction calculation [H. Tatewaki and Y. Watanabe, Chem. Phys. 389, 58 (2011)]. Manifestations of Brown-Ravenhall disease were noted at higher values of nuclear charge and ultimately limited the accuracy of the Poet-Temkin model energy. The energies from a no-pair configuration interaction (CI) calculation (the negative-energy states for the appropriate hydrogen-like ion were excluded from the CI expansion) were found to be different from the unrestricted B-spline calculation.
Collapse
Affiliation(s)
- Li-Yan Tang
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, People's Republic of China
| | | | | | | |
Collapse
|
28
|
Analytical Nuclear Gradients of Density-Fitted Dirac–Fock Theory with a 2-Spinor Basis. J Chem Theory Comput 2013; 9:4300-3. [DOI: 10.1021/ct400719d] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|