Kramers' unrestricted Hartree–Fock and second-order Møller–Plesset perturbation methods using relativistic effective core potentials with spin–orbit operators: Test calculations for HI and CH3I

The Kramers' restricted Hartree–Fock (KRHF) and second-order Møller–Plesset perturbation (KRMP2) methods using relativistic effective core potentials (RECP) with spin–orbit operators and two-component spinors are extended to the unrestricted forms, KUHF and KUMP2. As in the conventional unrestricted methods, the KUHF and KUMP2 methods are capable of qualitatively describing the bond breaking for a single bond. As a result, it is possible to estimate spin–orbit effects along the dissociation curve at the HF and MP2 levels of theory as is demonstrated by the test calculations on the ground states of HI and CH₃I. Since the energy lowering due to spin–orbit interactions is larger for the I atom than for the closed-shell molecules, dissociation energies are reduced and bond lengths are slightly elongated by the inclusion of the spin–orbit interactions. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 91–98, 1998     

Spin-orbit effects on structures of closed-shell polyatomic molecules containing heavy atoms calculated by two-component Hartree–Fock method

We have implemented geometry optimization using an analytic gradient to a two-component Kramers' restricted Hartree–Fock (KRHF) method for polyatomic molecules with closed-shell configurations. The KRHF method is a Hartree–Fock method based on relativistic effective core potentials with effective spin-orbit operators. The derivatives of spin-orbit integrals are obtained by numerical differentiation. Geometries for the various forms of polyatomic hydrides containing row 6 p-block elements are optimized with and without spin-orbit interactions. The structural changes due to spin-orbit interactions are small, but show definite trends, which correlate well with the p_1/2 spinor population. Atomization energies are reduced significantly by incorporating spin-orbit interactions for all molecules considered. The KRHF calculations of several methylhalides demonstrate that the spinor energies from the KRHF method can be useful for the interpretation of experimental photoelectron spectra of molecules exhibiting spin-orbit splittings. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1526–1533, 1998     

A comparison of Hartree—Fock,MP2, and DFT results for the HCN dimer and crystal

A number of hydrogen-bond related quantities—geometries, interaction energies, dipole moments, dipole moment derivatives, and harmonic vibrational frequencies—were calculated at the Hartree—Fock, MP2, and different DFT levels for the HCN dimer and the periodic HCN crystal. The crystal calculations were performed with the Hartree—Fock program CRYSTAL92, which routinely allows an a posteriori electron-correlation correction of the Hartree—Fock obtained lattice energy using different correlation-only functionals. Here, we have gone beyond this procedure by also calculating the electron-correlation energy correction during the structure optimization, i.e., after each CRYSTAL92 Hartree—Fock energy evaluation, the a posteriori density functional scheme was applied. In a similar manner, we optimized the crystal structure at the MP2 level, i.e., for each Hartree—Fock CRYSTAL92 energy evaluation, an MP2 correction was performed by summing the MP2 pair contributions from all HCN molecules within a specified cutoff distance. The crystal cell parameters are best reproduced at the Hartree—Fock and the nongradient-corrected HF + LDA and HF + VWN levels. The BSSE-corrected MP2 method and the HF + P91, HF + LDA, and HF + VWN methods give lattice energies in close agreement with the ZPE-corrected experimental lattice energy. The (HCN)₂ dimer properties are best reproduced at the MP2 level, at the gradient-corrected DFT levels, and with the B3LYP and BHHLYP methods. © 1996 John Wiley & Sons, Inc.     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (II): The restricted open‐shell approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent two‐component relativistic calculations, termed the Kramers‐restricted open‐shell HF (KROHF) method, is developed. The present KROHF method is defined as a relativistic analogue of ROHF using time‐reversal symmetry and quaternion algebra, based on the Kramers‐unrestricted HF (KUHF) theory reported in our previous study (Int. J. Quantum Chem., doi: 10.1002/qua.25356 ). As seen in the nonrelativistic ROHF theory, the ambiguity of the KROHF Fock operator gives physically meaningless spinor energies. To avoid this problem, the canonical parametrization of KROHF to satisfy Koopmans' theorem is also discussed based on the procedure proposed by Plakhutin et al. (J. Chem. Phys. 2006 , 125, 204110). Numerical assessments confirmed that KROHF using Plakhutin's canonicalization procedure correctly gives physical spinor energies within the frozen‐orbital approximation under spin–orbit interactions.     

Atomic Xα calculations based on EHFS(α) = Eexp

The total energies and one-electron energies for first- and second-row atoms were calculated by using the Hartree–Fock and the Hartree–Fock-Slater Hamiltonian with X_α orbitals, u_i(α_exp); α was parametrized from E^HFS (α_exp) = E^exp. The E^HF (α_exp) total energies are always higher than the Hartree–Fock energies for the atoms. The relation of the calculated ionization potential to the experimental ionization potential depends on the α used to define u_i(α), α_exp, or α_HF.     

Ab initio calculation of the Dzyaloshinskii–Moriya parameters: Spin–orbit GSO‐HF,DFT, and CI approaches

We present ab initio methods to determine the Dzyaloshinskii–Moriya (DM) parameter, which provides the anisotropic effects of noncollinear spin systems. For this purpose, we explore various general spin orbital (GSO) approaches, such as Hartree–Fock (HF), density functional theory (DFT), and configuration interaction (CI), with one‐electron spin–orbit coupling (SOC1). As examples, two simple D_3h‐symmetric models, H₃ and B(CH₂)₃, are examined. Implications of the computational results are discussed in relation to as isotropic and anisotropic interactions of molecular‐based magnets. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Supplementary d and f functions in molecular wave functions at large and small internuclear separations

The CO, CO₂, CS, CIF, and SO₂ molecules were used to test the dependence of supplementary d and f function exponents to changes in bond lengths and bond angles in MO calculations utilizing Gaussian basis sets in Hartree–Fock and Moller–Plesset calculations. Using Dunning–Hay double zeta basis sets, optimizations were performed at internuclear separations from 100–200 pm and beyond. The energy cost of not reoptimizing d function exponents when bonds are stretched or compressed is much smaller for correlated calculations than for those at the Hartree–Fock level and is greatest at the lower end of the range of internuclear distances. The problem is much less serious at all levels when multiple sets of d functions are used. © 1993 John Wiley & Sons, Inc.     

Ab initio Hartree–Fock and multiple-scattering Xα calculation of the g factors for CuF2

Ab initio Hartree–Fock and multiple-scattering wave functions are calculated for linear CuF₂. These wave functions are used to calculate the spin–orbit coupling in a new way where the neglect of two- and many-center terms is avoided and where experimental or calculated spin–orbit coupling constants for the atomic ions are used. The calculated value of g is too small by the MS Xα method and too large by the ab initio method, indicating too much 3d–L interaction in the MS Xα case and too little in the ab initio case.     

HF and MP2 calculations on CN−, N2, AlF,SiO, PN,SC, ClB,and P2 using correlated molecular wave functions

Contracted basis sets of double zeta valence quality plus polarization functions (DZP) and augmented DZP basis sets, which were recently constructed for the first‐ and second‐row atoms, are applied to study the electronic ground states of the diatomic molecules CN^?, N₂, AlF, SiO, PN, SC, ClB, and P₂. At the Hartree–Fock (HF) and/or Møller–Plesset second‐order (MP2) levels, total and molecular orbital energies, dissociation energies, bond lengths, harmonic vibrational frequencies, and dipole moments are calculated and compared with available experimental data and with the results obtained from correlation consistent polarized valence basis sets of Dunning's group. For N₂, calculations of polarizabilities at the HF and MP2 levels with the sets presented above are also done and compared with results reported in the literature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Nonrelativistic and quasirelativistic model potential calculations on AgH and Ag2

The major relativistic effects are included into the model potential (MP) method of Bonifacic and Huzinaga. The effects are incorporated on the level of Cowan and Griffin's relativistic Hartree–Fock (RHF) method. The model potential parameters are determined using the results of nonrelativistic Hartree–Fock (NHF) and RHF calculations. A new scheme of selection of the basis functions for use in atomic and molecular MP calculations is proposed. To obtain agreement with the Hartree–Fock calculations on AgH and Ag₂, the 4p shell has to be included explicitly in the MP calculations. The explicit treatment of the 4p electrons and the resulting reduction of the core size are necessary in order to overcome difficulties with approximate representation of the large 4p–4d core-valence interactions on the MP level.     

Spin–other–orbit and spin–orbit interactions in ƒ2 and ƒ3 electron configurations

Expressions of the matrix elements of the spin–other–orbit and spin–orbit interactions for the various multiplets of all the states of ?²- and ?³-electron configurations are reported and used to evaluate the Hartree–Fock values of these interactions in the neutral atoms Ce(4?²), Pr(4?³), Ho(4?¹¹) and Er(4?¹²). The required values of the spin–spin parameters M, and the spin-orbit parameter ζ for these atoms were obtained using numerical Hartree–Fock wave functions.     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (I): The unrestricted approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent, two‐component relativistic calculations, termed the Kramers‐unrestricted HF (KUHF) method, is developed. The present KUHF method, which is formulated as a relativistic counterpart of nonrelativistic UHF, is based on quaternion algebra and partly uses time‐reversal symmetry. The fundamental characteristics of KUHF are discussed in this study. From numerical assessments, it was revealed that KUHF gives a corresponding solution to nonrelativistic UHF; furthermore, KUHF properly describes spin‐orbit interactions. In addition, KUHF can improve the self‐consistent field convergence behavior in spin‐dependent calculations, for example, for f‐block elements.     

Multiplicity,instability, and SCF convergence problems in Hartree–Fock solutions

We present a study of the instability and convergence of Hartree–Fock (HF) ab initio solutions for the diatomic systems H₂, LiH, CH, C₂, and N₂. In our study, we consider real molecular orbitals (MOs) and analyze the classes of single‐determinant functions associated to Hartree–Fock–Roothaan (HFR) and Hartree–Fock–Pople–Nesbet (HFPN) equations. To determine the multiple HF solutions, we used either an SCF iterative procedure with aufbau and non‐aufbau ordering rules or the algebraic method (AM). Stability conditions were determined using TICS and ASDW stability matrices, derived from the maximum and minimum method of functions (MMF). We examined the relationship between pure SCF convergence criterion with the aufbau ordering rule, and the classification of the HF solution as an extremum point in its respective class of functions. Our results show that (i) in a pure converged SCF calculation, with the aufbau ordering rule, the solutions are not necessarily classified as a minimum of the HF functional with respect to the TICS or ASDW classes of solutions, and (ii) for all studied systems, we obtained local minimum points associated only with the aufbau rule and the solutions of lower energies. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 600–610, 2000     

Theoretical investigations on thallium halides: Relativistic and electron correlation effects in T1X and T1X3 compounds (XF,C1, Br,and I)

Relativistic and electron correlation effects in thallium halides TlX and TlX₃ (X?F, Cl, Br, and I) are investigated by extensive ab initio configuration interaction calculations. Spin–orbit coupling is included at the Hartree–Fock level for the diatomic TlBr and TlI. At the best level of treatment of electron correlation (quadratic configuration interaction), the calculated molecular properties are in good agreement with experimental results, i.e., for the diatomic thallium halides deviations from experimental values are <0.06 Å for bond distances, <0.14 mdyn/Å for force constants, <35 kJ/mol for dissociation energies, and <0.3 D for dipole moments. The convergence of the Møller–Plesset series up to the fourth order is discussed. Two alternative structures of TlI₃ are compared. At the Møller–Plesset level of theory, the trigonal planar structure with thallium in the oxidation state + 3 is the preferred gas phase arrangement compared with the bent arrangement containing a linear I unit and thallium in the oxidation state + 1, the difference being ca. 95 kJ/mol. Vibrational frequencies are predicted for all trigonal planar thallium(III) halides. © 1993 John Wiley & Sons, Inc.     

Ab initio investigation of water clusters (H2O)n (n = 2–34)

Various properties of typical structures of water clusters in the n = 2–34 size regime with the change of cluster size have been systematically explored. Full optimizations are carried out for the structures presented in this article at the Hartree–Fock (HF) level using the 6‐31G(d) basis set by taking into account the positions of all atoms within the cluster. The influence of the HF level on the results has been reflected by the comparison between the binding energies of (H₂O)_n (n = 2–6, 8, 11, 13, 20) calculated at the HF level and those obtained from high‐level ab initio calculations at the second‐order Møller–Plesset (MP2) perturbation theory and the coupled cluster method including singles and doubles with perturbative triples (CCSD(T)) levels. HF is inaccurate when compared with MP2 and CCSD(T), but it is more practical and allows us to study larger systems. The computed properties characterizing water clusters (H₂O)_n (n = 2–34) include optimal structures, structural parameters, binding energies, hydrogen bonds, charge distributions, dipole moments, and so on. When the cluster size increases, trends of the above various properties have been presented to provide important reference for understanding and describing the nature of the hydrogen bond. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

A microsolvation approach to the prediction of the relative enthalpies and free energies of hydration for ammonium ions

Hartree–Fock (HF) and second-order Møller–Plesset (MP2) calculations were used to investigate the structures and thermochemistry of methylammonium–water clusters (Me_4-m NH _m ⁺ (H₂O) _n , m=1–4, n=1–4). Water molecules were treated ab initio and with effective fragment potentials (EFP). In addition to a thorough phase-space search, the importance of basis set, electron correlation, and thermodynamic effects was systematically examined. Cluster structures resulted from hydrogen bond formation between the ammonium group and water molecules; upon saturation of the hydrogen bonding sites of the ammonium group, water molecules entered the second hydration shell. With only four water molecules, the experimental relative enthalpies of hydration were well reproduced at the HF level, while the MP2 relative free energies were in best agreement with experiment. Absolute energies of hydration were calculated using an empirical correction. These results strongly suggest that a HF-based microsolvation approach employing a small number of water molecules can be used to compute relative enthalpies of hydration.     

Accurate Quantum‐Chemical Prediction of Enthalpies of Formation of Small Molecules in the Gas Phase

The coupled‐cluster approach, including single and double excitations and perturbative corrections for triple excitations, is capable of predicting molecular electronic energies and enthalpies of formation of small molecules in the gas phase with very high accuracy (specifically, with error bars less than 5 kJ mol^?1), provided that the electronic wavefunction is dominated by the Hartree–Fock configuration. This capability is illustrated by calculations on molecules containing O–H and O–F bonds, namely OH, FO, H₂O, HOF, and F₂O. To achieve this very high accuracy, it is imperative to account for electron‐correlation effects in a quantitative manner, either by using explicitly correlated two‐particle basis functions (R12 functions) or by extrapolating to the limit of a complete basis. Besides taking into account harmonic zero‐point vibrational energies, it is also necessary to account for anharmonic corrections to the zero‐point vibrational energies, to include the core orbitals into the coupled‐cluster calculations, and to account for spin–orbit corrections and scalar relativistic effects. These additional corrections constitute small but significant contributions in the range of 1–4 kJ mol^?1 to the enthalpies of formation of the aforementioned molecules. The highly accurate coupled‐cluster results, obtained by employing R12 functions and by including various corrections, are compared with standard Kohn–Sham density‐functional calculations as well as with the Gaussian‐2 and complete‐basis‐set model chemistries.     

Geometry,basis set,and correlation energy dependence of computed protonation energies of imino bases

The nitrogen protonation energies of the imino bases HN?CHR, where R is H, CH₃, NH₂, OH, and F, have been evaluated to determine the dependence of absolute and relative protonation energies on geometry, basis set, and correlation effects. Reliable absolute protonation energies require a basis set larger than a split-valence plus polarization basis, the inclusion of correlation, and optimized geometries of at least Hartree–Fock 4-31G quality. Consistent relative protonation energies can be obtained at the Hartree–Fock level with smaller basis sets. Extending the split-valence basis set by the addition of polarization functions on all atoms decreases the computed absolute Hartree–Fock nitrogen protonation energies of the imino bases HN?CHR except when R is F, but increases the oxygen protonation energies of the carbonyl bases O?CHR.