Local unitary transformation method toward practical electron correlation calculations with scalar relativistic effect in large-scale molecules

Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory

Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.

Range Separation Method for Density Functional Theory Based on Two-Electron Infinite-Order Two-Component Hamiltonian

The range separation method for density functional theory (DFT) was extended to a two-component relativistic theory based on the unitary transformation of one- and two-electron operators and a density operator. In the framework of the spin-free infinite-order two-component Hamiltonian, we implemented several types of two-electron integrals of range-separated two-electron interactions arising from the unitary transformation. Numerical assessments were performed using long-range-corrected (LC)-DFT, which utilizes the range separation of an exchange functional. The present method successfully reproduced the reference values obtained by the four-component LC-DFT calculations when the whole unitary transformations of one-electron, full-range, and range-separated two-electron operators and a density operator were considered. An efficient scheme for the unitary transformation, which is termed the local unitary transformation (LUT), was also applied to the range-separated two-electron term and other operators. The LUT method reduced the computational costs of the LC-DFT calculations significantly without any loss of accuracy.

Toward a correct treatment of core properties with local hybrid functionals

In local hybrid functionals (LHs), a local mixing function (LMF) determines the position-dependent exact-exchange admixture. We report new LHs that focus on an improvement of the LMF in the core region while retaining or partly improving upon the high accuracy in the valence region exhibited by the LH20t functional. The suggested new pt-LMFs are based on a Padé form and modify the previously used ratio between von Weizsäcker and Kohn-Sham local kinetic energies by different powers of the density to enable flexibly improved approximations to the correct high-density and iso-orbital limits relevant for the innermost core region. Using TDDFT calculations for a set of K-shell core excitations of second- and third-period systems including accurate state-of-the-art relativistic orbital corrections, the core part of the LMF is optimized, while the valence part is optimized as previously reported for test sets of atomization energies and reaction barriers (Haasler et al., J Chem Theory Comput 2020, 16, 5645). The LHs are completed by a calibration function that minimizes spurious nondynamical correlation effects caused by the gauge ambiguities of exchange-energy densities, as well as by B95c meta-GGA correlation. The resulting LH23pt functional relates to the previous LH20t functional but specifically improves upon the core region.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Relativistic Segmented Correlation Consistent Basis Sets for the 5p and 6p Elements

The relativistic all-electron triple-ζ (TZ) and quadruple-ζ (QZ) correlation-consistent basis sets for the 5p and 6p elements were reoptimized to have a segmented contraction scheme. Properties computed with the segmented basis sets using the coupled-cluster level of theory with single, double, and perturbative triple excitations closely match the results obtained using the generally contracted analogues. Deviations in the ionization potentials and electron affinities computed using the segmented basis sets compared with those computed with the generally contracted basis sets are within 1 kcal mol-1. The relative deviation in the computed bond lengths for a selection of diatomic molecules is within 0.02 Å, and the computed harmonic vibrational frequencies differ by less than 3%. The mean absolute deviations (MADs) in the computed bond dissociation energies are 1.01 and 1.31 kcal mol-1 for the TZ and QZ basis sets, respectively, when only the valence electrons are correlated and 0.14 and 0.10 kcal mol-1 for the TZ and QZ basis sets, respectively, when the valence and outer-core electrons are correlated. The segmented basis sets also retain the systematically convergent behavior intrinsic to the correlation-consistent family of basis sets while affording speedups for the average time required to form the Fock matrix from 32.8 to 82.9× compared to the time required with the generally contracted versions when used with integral algorithms employed by many computational chemistry software packages.

RAQET: Large-scale two-component relativistic quantum chemistry program package

The Relativistic And Quantum Electronic Theory (RAQET) program is a new software package, which is designed for large-scale two-component relativistic quantum chemical (QC) calculations. The package includes several efficient schemes and algorithms for calculations involving large molecules which contain heavy elements in accurate relativistic formalisms. These calculations can be carried out in terms of the two-component relativistic Hamiltonian, wavefunction theory, density functional theory, core potential scheme, and evaluation of electron repulsion integrals. Furthermore, several techniques, which have frequently been used in non-relativistic QC calculations, have been customized for relativistic calculations. This article introduces the brief theories and capabilities of RAQET with several calculation examples. © 2018 Wiley Periodicals, Inc.

Computerized implementation of higher-order electron-correlation methods and their linear-scaling divide-and-conquer extensions

We have implemented a linear-scaling divide-and-conquer (DC)-based higher-order coupled-cluster (CC) and Møller-Plesset perturbation theories (MPPT) as well as their combinations automatically by means of the tensor contraction engine, which is a computerized symbolic algebra system. The DC-based energy expressions of the standard CC and MPPT methods and the CC methods augmented with a perturbation correction were proposed for up to high excitation orders [e.g., CCSDTQ, MP4, and CCSD(2)_TQ ]. The numerical assessment for hydrogen halide chains, polyene chains, and first coordination sphere (C1) model of photoactive yellow protein has revealed that the DC-based correlation methods provide reliable correlation energies with significantly less computational cost than that of the conventional implementations. © 2017 Wiley Periodicals, Inc.

Implementation of Analytical Energy Gradient of Spin-Dependent General Hartree-Fock Method Based on the Infinite-Order Douglas-Kroll-Hess Relativistic Hamiltonian with Local Unitary Transformation

An analytical energy gradient for the spin-dependent general Hartree-Fock method based on the infinite-order Douglas-Kroll-Hess (IODKH) method was developed. To treat realistic systems, the local unitary transformation (LUT) scheme was employed both in energy and energy gradient calculations. The present energy gradient method was numerically assessed to investigate the accuracy in several diatomic molecules containing fifth- and sixth-period elements and to examine the efficiency in one-, two-, and three-dimensional silver clusters. To arrive at a practical calculation, we also determined the geometrical parameters of fac-tris(2-phenylpyridine)iridium and investigated the efficiency. The numerical results confirmed that the present method describes a highly accurate relativistic effect with high efficiency. The present method can be a powerful scheme for determining geometries of large molecules, including heavy-element atoms.

Extension of accompanying coordinate expansion and recurrence relation method for general-contraction basis sets

An algorithm of the accompanying coordinate expansion and recurrence relation (ACE-RR), which is used for the rapid evaluation of the electron repulsion integral (ERI), has been extended to the general-contraction (GC) scheme. The present algorithm, denoted by GC-ACE-RR, is designed for molecular calculations including heavy elements, whose orbitals consist of many primitive functions with and without higher angular momentum such as d- and f-orbitals. The performance of GC-ACE-RR was assessed for (ss|ss)-, (pp|pp)-, (dd|dd)-, and (ff|ff)-type ERIs in terms of contraction length and the number of GC orbitals. The present algorithm was found to reduce the central processing unit time compared with the ACE-RR algorithm, especially for higher angular momentum and highly contracted orbitals. Compared with HONDOPLUS and GAMESS program packages, GC-ACE-RR computations for ERIs of three-dimensional gold clusters Aun (n = 1, 2, …, 10, 15, 20, and 25) are more than 10 times faster.