Implementation of the diagonalization‐free algorithm in the self‐consistent field procedure within the four‐component relativistic scheme

A recently developed Thouless‐expansion‐based diagonalization‐free approach for improving the efficiency of self‐consistent field (SCF) methods (Noga and Šimunek, J. Chem. Theory Comput. 2010, 6, 2706) has been adapted to the four‐component relativistic scheme and implemented within the program package ReSpect. In addition to the implementation, the method has been thoroughly analyzed, particularly with respect to cases for which it is difficult or computationally expensive to find a good initial guess. Based on this analysis, several modifications of the original algorithm, refining its stability and efficiency, are proposed. To demonstrate the robustness and efficiency of the improved algorithm, we present the results of four‐component diagonalization‐free SCF calculations on several heavy‐metal complexes, the largest of which contains more than 80 atoms (about 6000 4‐spinor basis functions). The diagonalization‐free procedure is about twice as fast as the corresponding diagonalization. © 2014 Wiley Periodicals, Inc.

[1]  Leonardo Belpassi,et al.  RSC_CP_C1CP20569B 1..27 , 2011 .

[2]  F. Jensen.  Erratum: “Polarization consistent basis sets: Principles” [J. Chem. Phys. 115, 9113 (2001)] , 2002 .

[3]  Marco Häser,et al.  Improvements on the direct SCF method , 1989 .

[4]  P. Dirac.  The quantum theory of the electron , 1928 .

[5]  Frank Jensen,et al.  Polarization consistent basis sets: Principles , 2001 .

[6]  Florian Weigend,et al.  A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency , 2002 .

[7]  Kimihiko Hirao,et al.  A new computational scheme for the Dirac–Hartree–Fock method employing an efficient integral algorithm , 2001 .

[8]  J. Leszczynski,et al.  QM/MM Studies of Light-responsive Biological Systems , 2021, Challenges and Advances in Computational Chemistry and Physics.

[9]  W. Kutzelnigg.  Basis set expansion of the dirac operator without variational collapse , 1984 .

[10]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[11]  J. Noga,et al.  Solving the Independent-Particle Model via Nonunitary Transformation Based on Variational Coupled Cluster Singles. , 2010, Journal of chemical theory and computation.

[12]  Florian Weigend,et al.  Acceleration of self‐consistent‐field convergence by combining conventional diagonalization and a diagonalization‐free procedure , 2011, J. Comput. Chem..

[13]  W. Kutzelnigg.  Solved and unsolved problems in relativistic quantum chemistry , 2012 .

[14]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[15]  K. Dyall,et al.  Revised relativistic basis sets for the 5d elements Hf–Hg , 2009 .

[16]  Lorenzo Pisani,et al.  Relativistic Dirac–Fock calculations for closed‐shell molecules , 1994, J. Comput. Chem..

[17]  H. Quiney,et al.  Electron density fitting for the Coulomb problem in relativistic density-functional theory. , 2006, The Journal of chemical physics.

[18]  T. Shiozaki.  Communication: An efficient algorithm for evaluating the Breit and spin-spin coupling integrals. , 2013, The Journal of chemical physics.

[19]  T. Shiozaki,et al.  Large-scale Dirac-Fock-Breit method using density fitting and 2-spinor basis functions. , 2013, The Journal of chemical physics.

[20]  K. Hirao,et al.  Pseudospectral approach to relativistic molecular theory. , 2004, The Journal of chemical physics.

[21]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[22]  Pekka Pyykkö,et al.  Relativistic effects in structural chemistry , 1988 .

[23]  H. Quiney,et al.  Relativistic density functional theory using Gaussian basis sets , 2002 .

[24]  Jörg Kussmann,et al.  Linear‐scaling self‐consistent field methods for large molecules , 2013 .

[25]  E. Engel,et al.  Systematic study of the Ib diatomic molecules Cu 2 , Ag 2 , and Au 2 using advanced relativistic density functionals , 1999 .

[26]  D. Thouless.  Stability conditions and nuclear rotations in the Hartree-Fock theory , 1960 .

[27]  Trygve Helgaker,et al.  The augmented Roothaan-Hall method for optimizing Hartree-Fock and Kohn-Sham density matrices. , 2008, The Journal of chemical physics.

[28]  Trygve Helgaker,et al.  Principles of direct 4-component relativistic SCF: application to caesium auride , 1997 .

[29]  Trygve Helgaker,et al.  Four‐component relativistic Kohn–Sham theory , 2002, J. Comput. Chem..

[30]  With contributions from , 2007 .

[31]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[32]  M. Kaupp,et al.  A fully relativistic method for calculation of nuclear magnetic shielding tensors with a restricted magnetically balanced basis in the framework of the matrix Dirac-Kohn-Sham equation. , 2008, The Journal of chemical physics.

[33]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[34]  O. Malkina,et al.  Restricted magnetically balanced basis applied for relativistic calculations of indirect nuclear spin-spin coupling tensors in the matrix Dirac-Kohn-Sham framework , 2009 .

[35]  P. Dirac.  The Quantum Theory of the Electron. Part II , 1928 .

[36]  H. Quiney,et al.  Poisson-transformed density fitting in relativistic four-component Dirac-Kohn-Sham theory. , 2008, The Journal of chemical physics.

[37]  Trond Saue,et al.  Relativistic Hamiltonians for chemistry: a primer. , 2011, Chemphyschem : a European journal of chemical physics and physical chemistry.

[38]  T. Fleig.  Invited review: Relativistic wave-function based electron correlation methods , 2012 .

[39]  K. Dyall.  Relativistic Quadruple-Zeta and Revised Triple-Zeta and Double-Zeta Basis Sets for the 4p, 5p, and 6p Elements , 2006 .

[40]  Wenjian Liu.  Perspectives of relativistic quantum chemistry: the negative energy cat smiles. , 2012, Physical chemistry chemical physics : PCCP.

[41]  Michael Dolg,et al.  The Beijing four-component density functional program package (BDF) and its application to EuO, EuS, YbO and YbS , 1997 .

[42]  Lucas Visscher,et al.  RELATIVISTIC QUANTUM-CHEMISTRY - THE MOLFDIR PROGRAM PACKAGE , 1994 .

[43]  M. Head‐Gordon,et al.  An exact reformulation of the diagonalization step in electronic structure calculations as a set of second order nonlinear equations. , 2004, The Journal of chemical physics.

[44]  Lucas Visscher,et al.  Approximate molecular relativistic Dirac-Coulomb calculations using a simple Coulombic correction , 1997 .

[45]  P. Taylor,et al.  All-electron molecular Dirac-Hartree-Fock calculations - The group IV tetrahydrides CH4, SiH4, GeH4, SnH4, and PbH4 , 1991 .

[46]  Kimihiko Hirao,et al.  A new implementation of four-component relativistic density functional method for heavy-atom polyatomic systems , 2001 .

[47]  Harry M. Quiney,et al.  Ab initio relativistic quantum chemistry: four-components good, two-components bad!* , 1998 .

[48]  Richard E. Stanton,et al.  Kinetic balance: A partial solution to the problem of variational safety in Dirac calculations , 1984 .

[49]  Lucas Visscher,et al.  The Dirac equation in quantum chemistry: Strategies to overcome the current computational problems , 2002, J. Comput. Chem..

[50]  Kimihiko Hirao,et al.  A highly efficient algorithm for electron repulsion integrals over relativistic four-component Gaussian-type spinors , 2002 .

[51]  Sergio Rampino,et al.  Efficient Parallel All-Electron Four-Component Dirac-Kohn-Sham Program Using a Distributed Matrix Approach II. , 2010, Journal of chemical theory and computation.