Junji Seino

Local unitary transformation method for large-scale two-component relativistic calculations: Case for a one-electron Dirac Hamiltonian

An accurate and efficient scheme for two-component relativistic calculations at the spin-free infinite-order Douglas-Kroll-Hess (IODKH) level is presented. The present scheme, termed local unitary transformation (LUT), is based on the locality of the relativistic effect. Numerical assessments of the LUT scheme were performed in diatomic molecules such as HX and X(2) (X = F, Cl, Br, I, and At) and hydrogen halide clusters, (HX)(n) (X = F, Cl, Br, and I). Total energies obtained by the LUT method agree well with conventional IODKH results. The computational costs of the LUT method are drastically lower than those of conventional methods since in the former there is linear-scaling with respect to the system size and a small prefactor.

Unveiling a New Aspect of Simple Arylboronic Esters: Long-Lived Room-Temperature Phosphorescence from Heavy-Atom-Free Molecules

Arylboronic esters can be used as versatile reagents in organic synthesis, as represented by Suzuki-Miyaura cross-coupling. Here we report a serendipitous finding that simple arylboronic esters are phosphorescent in the solid state at room temperature with a lifetime on the order of several seconds. The phosphorescence properties of arylboronic esters are remarkable in light of the general notion that phosphorescent organic molecules require heavy atoms and/or carbonyl groups for the efficient generation of a triplet excited state. Theoretical calculations on phenylboronic acid pinacol ester indicated that this molecule undergoes an out-of-plane distortion at the (pinacol)B-Cipso moiety in the T1 excited state, which is responsible for its phosphorescence. A compound survey with 19 arylboron compounds suggested that the phosphorescence properties might be determined by solid-state molecular packing rather than by the patterns and numbers of boron substituents on the aryl units. The present finding may update the general notion of phosphorescent organic molecules.

Local unitary transformation method for large-scale two-component relativistic calculations. II. Extension to two-electron Coulomb interaction

The local unitary transformation (LUT) scheme at the spin-free infinite-order Douglas-Kroll-Hess (IODKH) level [J. Seino and H. Nakai, J. Chem. Phys. 136, 244102 (2012)], which is based on the locality of relativistic effects, has been extended to a four-component Dirac-Coulomb Hamiltonian. In the previous study, the LUT scheme was applied only to a one-particle IODKH Hamiltonian with non-relativistic two-electron Coulomb interaction, termed IODKH/C. The current study extends the LUT scheme to a two-particle IODKH Hamiltonian as well as one-particle one, termed IODKH/IODKH, which has been a real bottleneck in numerical calculation. The LUT scheme with the IODKH/IODKH Hamiltonian was numerically assessed in the diatomic molecules HX and X(2) and hydrogen halide molecules, (HX)(n) (X = F, Cl, Br, and I). The total Hartree-Fock energies calculated by the LUT method agree well with conventional IODKH/IODKH results. The computational cost of the LUT method is reduced drastically compared with that of the conventional method. In addition, the LUT method achieves linear-scaling with respect to the system size and a small prefactor.

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

In order to perform practical electron correlation calculations, the local unitary transformation (LUT) scheme at the spin-free infinite-order Douglas-Kroll-Hess (IODKH) level [J. Seino and H. Nakai, J. Chem. Phys. 136, 244102 (2012); and ibid. 137, 144101 (2012)], which is based on the locality of relativistic effects, has been combined with the linear-scaling divide-and-conquer (DC)-based Hartree-Fock (HF) and electron correlation methods, such as the second-order Mo̸ller-Plesset (MP2) and the coupled cluster theories with single and double excitations (CCSD). Numerical applications in hydrogen halide molecules, (HX)n (X = F, Cl, Br, and I), coinage metal chain systems, Mn (M = Cu and Ag), and platinum-terminated polyynediyl chain, trans,trans-{(p-CH3C6H4)3P}2(C6H5)Pt(C≡C)4Pt(C6H5){(p-CH3C6H4)3P}2, clarified that the present methods, namely DC-HF, MP2, and CCSD with the LUT-IODKH Hamiltonian, reproduce the results obtained using conventional methods with small computational costs. The combination of both LUT and DC techniques could be the first approach that achieves overall quasi-linear-scaling with a small prefactor for relativistic electron correlation calculations.

Analytical energy gradient based on spin-free infinite-order Douglas-Kroll-Hess method with local unitary transformation

In this study, the analytical energy gradient for the spin-free infinite-order Douglas-Kroll-Hess (IODKH) method at the levels of the Hartree-Fock (HF), density functional theory (DFT), and second-order Møller-Plesset perturbation theory (MP2) is developed. Furthermore, adopting the local unitary transformation (LUT) scheme for the IODKH method improves the efficiency in computation of the analytical energy gradient. Numerical assessments of the present gradient method are performed at the HF, DFT, and MP2 levels for the IODKH with and without the LUT scheme. The accuracies are examined for diatomic molecules such as hydrogen halides, halogen dimers, coinage metal (Cu, Ag, and Au) halides, and coinage metal dimers, and 20 metal complexes, including the fourth-sixth row transition metals. In addition, the efficiencies are investigated for one-, two-, and three-dimensional silver clusters. The numerical results confirm the accuracy and efficiency of the present method.

Semi-local machine-learned kinetic energy density functional with third-order gradients of electron density

A semi-local kinetic energy density functional (KEDF) was constructed based on machine learning (ML). The present scheme adopts electron densities and their gradients up to third-order as the explanatory variables for ML and the Kohn-Sham (KS) kinetic energy density as the response variable in atoms and molecules. Numerical assessments of the present scheme were performed in atomic and molecular systems, including first- and second-period elements. The results of 37 conventional KEDFs with explicit formulae were also compared with those of the ML KEDF with an implicit formula. The inclusion of the higher order gradients reduces the deviation of the total kinetic energies from the KS calculations in a stepwise manner. Furthermore, our scheme with the third-order gradient resulted in the closest kinetic energies to the KS calculations out of the presented functionals.

An ab initio study of nuclear volume effects for isotope fractionations using two‐component relativistic methods

We investigate the accuracy of two‐component Douglas–Kroll–Hess (DKH) methods in calculations of the nuclear volume term (≡ lnKnv) in the isotope fractionation coefficient. lnKnv is a main term in the chemical equilibrium constant for isotope exchange reactions in heavy element. Previous work based on the four‐component method reasonably reproduced experimental lnKnv values of uranium isotope exchange. In this work, we compared uranium reaction lnKnv values obtained from the two‐component and four‐component methods. We find that both higher‐order relativistic interactions and spin‐orbit interactions are essential for quantitative description of lnKnv. The best alternative is the infinite‐order Douglas–Kroll–Hess method with infinite‐order spin‐orbit interactions for the one‐electron term and atomic‐mean‐field spin‐same‐orbit interaction for the two‐electron term (IODKH‐IOSO‐MFSO). This approach provides almost equivalent results for the four‐component method, while being 30 times faster. The IODKH‐IOSO‐MFSO methodology should pave the way toward computing larger and more general molecules beyond the four‐component method limits.

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.

Accompanying coordinate expansion and recurrence relation method using a transfer relation scheme for electron repulsion integrals with high angular momenta and long contractions

An efficient algorithm for the rapid evaluation of electron repulsion integrals is proposed. The present method, denoted by accompanying coordinate expansion and transferred recurrence relation (ACE-TRR), is constructed using a transfer relation scheme based on the accompanying coordinate expansion and recurrence relation method. Furthermore, the ACE-TRR algorithm is extended for the general-contraction basis sets. Numerical assessments clarify the efficiency of the ACE-TRR method for the systems including heavy elements, whose orbitals have long contractions and high angular momenta, such as f- and g-orbitals.

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.