Ab initio thermochemistry involving heavy atoms: an investigation of the reactions Hg + IX (X = I, Br, Cl, O)

Accurate 0 K enthalpies have been calculated for reactions of mercury with a series of small iodine-containing molecules (I2, IBr, ICl, and IO). The calculations have been carried out with the coupled cluster singles and doubles method with a perturbative correction for connected triple excitations [CCSD(T)] using sequences of correlation consistent basis sets and accurate relativistic pseudopotentials. Corrections have been included to account for core-valence correlation, spin-orbit coupling, scalar relativity, and the Lamb shift. In a few cases coupled cluster calculations with iterative triple (CCSDT) and quadruple (CCSDTQ) excitations have been carried out to estimate the effects of higher order electron correlation. The pseudopotential calculations have also been compared to all electron calculations using second- and third-order Douglas-Kroll-Hess Hamiltonians. In addition to the reaction enthalpies, heats of formation, bond lengths, and harmonic vibrational frequencies have been calculated for the stable triatomic products HgI2, HgIBr, HgICl, and HgIO. Accurate dissociation energies, equilibrium bond lengths, and harmonic vibrational frequencies have also been calculated for each of the diatomic molecules involved in this study (HgI, HgBr, HgCl, HgO, I2, IBr, ICl, and IO). The reported enthalpies are expected to have accuracies of 1 kcal/mol or better.     

On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions

A coupled cluster composite approach has been used to accurately determine the spectroscopic constants, bond dissociation energies, and heats of formation for the X1(2)II(3/2) states of the halogen oxides ClO, BrO, and IO, as well as their negative ions ClO-, BrO-, and IO-. After determining the frozen core, complete basis set (CBS) limit CCSD(T) values, corrections were added for core-valence correlation, relativistic effects (scalar and spin-orbit), the pseudopotential approximation (BrO and IO), iterative connected triple excitations (CCSDT), and iterative quadruples (CCSDTQ). The final ab initio equilibrium bond lengths and harmonic frequencies for ClO and BrO differ from their accurate experimental values by an average of just 0.0005 A and 0.8 cm-1, respectively. The bond length of IO is overestimated by 0.0047 A, presumably due to an underestimation of molecular spin-orbit coupling effects. Spectroscopic constants for the spin-orbit excited X2(2)III(1/2) states are also reported for each species. The predicted bond lengths and harmonic frequencies for the closed-shell anions are expected to be accurate to within about 0.001 A and 2 cm-1, respectively. The dissociation energies of the radicals have been determined by both direct calculation and through use of negative ion thermochemical cycles, which made use of a small amount of accurate experimental data. The resulting values of D0, 63.5, 55.8, and 54.2 kcal/mol for ClO, BrO, and IO, respectively, are the most accurate ab initio values to date, and those for ClO and BrO differ from their experimental values by just 0.1 kcal/mol. These dissociation energies lead to heats of formation, DeltaH(f) (298 K), of 24.2 +/- 0.3, 29.6 +/- 0.4, and 29.9 +/- 0.6 kcal/mol for ClO, BrO, and IO, respectively. Also, the final calculated electron affinities are all within 0.2 kcal/mol of their experimental values. Improved pseudopotential parameters for the iodine atom are also reported, together with revised correlation consistent basis sets for this atom.     

Gradients for two-component quasirelativistic methods. Application to dihalogenides of element 116

The authors report the implementation of geometry gradients for quasirelativistic two-component Hartree-Fock and density functional methods using either the zero-order regular approximation Hamiltonian or spin-dependent effective core potentials. The computational effort of the resulting program is comparable to that of corresponding nonrelativistic calculations, as it is dominated by the evaluation of derivative two-electron integrals, which is the same for both types of calculations. Besides the implementation of derivatives of matrix elements of the one-particle Hamiltonian with respect to nuclear displacements, the calculation of the derivative exchange-correlation energy for the open shell case involves complicated expressions because of the noncollinear approach chosen to define the spin density. A pilot application to dihalogenides of element 116 shows how spin-orbit coupling strongly affects the chemistry of the superheavy p-block elements. While these molecules are bent at a scalar-relativistic level, spin-orbit coupling is so strong that only the 7p3/2 atomic orbitals of element 116 are involved in bonding, which favors linear molecular geometries for dihalogenides with heavy terminal halogen atoms.     

等电子三原子铀化物OUO2+、NUN和NUO+的结构和谐振频率的CCSD(T)计算研究

用小核相对论有效势和CCSD(T)方法计算了三原子铀化物OUO²⁺, NUN和NUO⁺的平衡键长和谐振频率. 计算结果显示U原子内层5s5p5d 电子相关能对这些化合物性质的影响非常小. 除NUN的弯曲振动频率,旋轨耦合效应对这些化合物的结构和频率的影响并不明显. 本文的计算结果与其他研究组的计算结果以及已有的实验值相比符合较好, 这表明作为单参考态方法, CCSD(T)能够对这些体系的键长和频率给出较精确的计算结果. 与此前密度泛函理论(DFT)的计算结果相比, CCSD(T)方法与PBE0泛函的结果吻合最好. 本文的工作有助于在用密度泛函方法研究这些体系时选择合适的交换相关泛函, 也为今后的实验研究提供了新的理论数据.     

Comparative study of relativistic density functional methods applied to actinide species AcO(2)(2+) and AcF(6) for Ac = U,Np

A two-component relativistic density functional method based on the Douglas-Kroll-Hess transformation has been applied to the actinyls and hexafluorides of U and Np. All-electron scalar relativistic calculations as well as calculations including spin-orbit interaction have been compared to results obtained with a pseudopotential approach. In addition, several exchange-correlation potentials have been applied to examine their performance for the bond lengths and vibrational frequencies of the title compounds. The calculations confirm the well-known accuracy of the LDA approach for geometries and frequencies, which is corroborated for the hexafluorides where gas phase experimental data are available. Comparison with results of accurate wave function based methods provides further confirmation of this finding. Gradient-corrected functionals tend to overestimate bond lengths and underestimate frequencies also for actinide compounds. The results obtained with Stoll-Preuss (small core) effective core potentials agree very well with those of all-electron calculations, while calculations with Hay-Martin large core pseudopotentials are somewhat less accurate. For all molecules and properties considered, spin-orbit effects have been found negligible concomitant with the closed-shell electronic structure of the U(VI) compounds and the open-shell situation of the Np(VI) compounds with a single valence f electron.     

Spin-orbit relativistic long-range corrected time-dependent density functional theory for investigating spin-forbidden transitions in photochemical reactions

A long-range corrected (LC) time-dependent density functional theory (TDDFT) incorporating relativistic effects with spin-orbit couplings is presented. The relativistic effects are based on the two-component zeroth-order regular approximation Hamiltonian. Before calculating the electronic excitations, we calculated the ionization potentials (IPs) of alkaline metal, alkaline-earth metal, group 12 transition metal, and rare gas atoms as the minus orbital (spinor) energies on the basis of Koopmans' theorem. We found that both long-range exchange and spin-orbit coupling effects are required to obtain Koopmans' IPs, i.e., the orbital (spinor) energies, quantitatively in DFT calculations even for first-row transition metals and systems containing large short-range exchange effects. We then calculated the valence excitations of group 12 transition metal atoms and the Rydberg excitations of rare gas atoms using spin-orbit relativistic LC-TDDFT. We found that the long-range exchange and spin-orbit coupling effects significantly contribute to the electronic spectra of even light atoms if the atoms have low-lying excitations between orbital spinors of quite different electron distributions.     

Relativistic all-electron two-component self-consistent density functional calculations including one-electron scalar and spin-orbit effects

We have implemented a Gaussian basis-set two-component self-consistent field method based on the fourth-order nuclear-only Douglas-Kroll-Hess approximation. Two-electron spin-orbit effects are included using Boettger's screened-nuclear spin-orbit approximation. In our two-component approach, the spin-orbit interaction is taken into account in a variational fashion employing a generalized Kohm-Sham scheme which allows one to work with hybrid density functionals. For open-shell systems we adopt the noncollinear spin-density approximation. Results are presented for equilibrium bond lengths, harmonic vibrational frequencies, and bond dissociation energies with local spin-density, generalized gradient approximation, and hybrid functionals in a set of benchmark molecules.     

pCCSD: parameterized coupled-cluster theory with single and double excitations

The primary characteristics of single reference coupled-cluster (CC) theory are size-extensivity and size-consistency, invariance under orbital rotations of the occupied or virtual space, the exactness of CC theory for N electron systems when the cluster operator is truncated to N-tuple excitations, and the relative insensitivity of CC theory to the choice of the reference determinant. In this work, we propose a continuous class of methods which display the desirable features of the coupled-cluster approach with single and double excitations (CCSD). These methods are closely related to the CCSD method itself and are inspired by the coupled electron pair approximation (CEPA). It is demonstrated that one can systematically improve upon CCSD and obtain geometries, harmonic vibrational frequencies, and total energies from a parameterized version of CCSD or pCCSD(α,β) by selecting a specific member from this continuous family of approaches. In particular, one finds that one such approach, the pCCSD(-1,1) method, is a significant improvement over CCSD for the calculation of equilibrium structures and harmonic frequencies. Moreover, this method behaves surprisingly well in the calculation of potential energy surfaces for single bond dissociation. It appears that this methodology has significant promise for chemical applications and may be particularly useful in applications to larger molecules within the framework of a high accuracy local correlation approach.     

Calculation of current densities using gauge-including atomic orbitals

A method for calculating the various components of the magnetically induced current-density tensor using gauge-including atomic orbitals is described. The method is formulated in the framework of analytical derivative theory, thus enabling implementation at the Hartree-Fock self-consistent-field (HF-SCF) as well as at electron-correlated levels. First-order induced current densities have been computed up to the coupled-cluster singles and doubles level (CCSD) augmented by a perturbative treatment of triple excitations [CCSD(T)] for carbon dioxide and benzene and up to the full coupled-cluster singles, doubles, and triples (CCSDT) level in the case of ozone. The applicability of the gauge including magnetically induced current method to larger molecules is demonstrated by computing first-order current densities for porphin and hexabenzocoronene at the HF-SCF and density-functional theory level. Furthermore, a scheme for obtaining quantitative values for the induced currents in a molecule via numerical integration over the current flow is presented. For benzene, a perpendicular magnetic field induces a (field dependent) ring current of 12.8 nA T(-1) at the HF-SCF level using a triple-zeta basis set augmented with polarization functions (TZP). At the CCSD(T)/TZP level the induced current was found to be 11.4 nA T(-1). Gauge invariance and its relation to charge-current conservation is discussed.     

Structure and bonding of the MCN molecules, M=Cu,Ag,Au,Rg

High-precision calculations are reported for the title series with M=Cu, Ag, Au, using CCSD(T) with the latest pseudopotentials and basis sets up to cc-pVQZ. The bond lengths for M=Cu, Ag, Au agree with experiment within better than 1 pm. The role of deep-core excitations is studied. The second-order spin-orbit effects are evaluated at the density functional theory level, including M=Rg. A qualitative bonding analysis suggests multiple M-C bonding. The calculated vibrational frequencies are expected to be more accurate than the available experimental estimates. The M-C bond lengths obey Cu相似文献   

Two-component relativistic methods for the heaviest elements

Different generalized Douglas-Kroll transformed Hamiltonians (DKn, n=1, 2,...,5) proposed recently by Hess et al. are investigated with respect to their performance in calculations of the spin-orbit splittings. The results are compared with those obtained in the exact infinite-order two-component (IOTC) formalism which is fully equivalent to the four-component Dirac approach. This is a comprehensive investigation of the ability of approximate DKn methods to correctly predict the spin-orbit splittings. On comparing the DKn results with the IOTC (Dirac) data one finds that the calculated spin-orbit splittings are systematically improved with the increasing order of the DK approximation. However, even the highest-order approximate two-component DK5 scheme shows certain deficiencies with respect to the treatment of the spin-orbit coupling terms in very heavy systems. The meaning of the removal of the spin-dependent terms in the so-called spin-free (scalar) relativistic methods for many-electron systems is discussed and a computational investigation of the performance of the spin-free DKn and IOTC methods for many-electron Hamiltonians is carried out. It is argued that the spin-free IOTC rather than the Dirac-Coulomb results give the appropriate reference for other spin-free schemes which are based on approximate two-component Hamiltonians. This is illustrated by calculations of spin-free DKn and IOTC total energies, r(-1) expectation values, ionization potentials, and electron affinities of heavy atomic systems.     

Analyzing molecular properties calculated with two-component relativistic methods using spin-free natural bond orbitals: NMR spin-spin coupling constants

An analysis method for static linear response properties employing two-component (spin-orbit) relativistic density functional theory along with scalar relativistic "natural localized molecular orbitals" (NLMOs) and "natural bond orbitals" (NBOs) has been developed. The spin-orbit NLMO/NBO analysis has been applied to study the indirect spin-spin coupling (J-coupling) constants in Tl-I, PbH(4), and a dinuclear Pt-Tl bonded complex with a very large Pt-Tl coupling constant (expt.: 146.8 kHz). For Tl-I it is shown that the analysis scheme based on scalar relativistic NLMOs is applicable even if spin-orbit coupling is responsible for most of the coupling's magnitude. For PbH(4) it is shown that electron delocalization plays a much larger role for the Pb-H coupling than it is the case for the C-H coupling in methane. For the Pt-Tl complex the analysis clearly demonstrates the strong influence of the ligands on the Pt-Tl coupling constant and quantifies the effect of the delocalization of the Pt-Tl bond on the Pt-Tl coupling constant.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

New coupled cluster approaches based on the unrestricted Hartree-Fock reference for treating molecules with multireference character

Here we review the basic formalism, implementation details, and performance of two newly developed coupled cluster (CC) methods based on the unrestricted Hartree-Fock (UHF) reference for treating molecules with multireference character. These two approaches can be considered to be approximations to the CC singles, doubles, and triples (CCSDT) method. The key concept of these two approaches is the corresponding orbitals, which are unitary transformations of canonical UHF molecular orbitals so that all spin orbitals are grouped into unique orbital pairs. In one approach called CCSDT(5P), a subset of triple excitations involving up to five-pair indices is included. In another approach called CCSD(T)-h, the contribution of connected triple excitations is treated in a hybrid way. With the concept of active corresponding orbitals, triple excitations can be automatically partitioned into two subsets, and the amplitudes of these two subsets are determined via solving different equations. Both CCSD(T)-h and CCSDT(5P) computationally scale as the seventh power of the system size. A survey of a number of applications demonstrates that CCSD(T)-h is an excellent approximation to the full CCSDT method, and CCSDT(5P) provides a good approximation to CCSDT for single-bond breaking processes. The overall performance of CCSDT(5P) is less accurate than that of CCSD(T)-h, but significantly better than that of the widely used CCSD(T).     

Analyzing NMR shielding tensors calculated with two-component relativistic methods using spin-free localized molecular orbitals

A recently developed analysis method [J. Chem. Phys. 127, 124106 (2007)] for NMR spin-spin coupling constants employing two-component (spin-orbit) relativistic density functional theory along with scalar relativistic natural localized molecular orbitals (NLMOs) and natural bond orbitals (NBOs) has been extended for analyzing NMR shielding tensors. Contributions from a field-dependent basis set (gauge-including atomic orbitals) have been included in the formalism. The spin-orbit NLMO/NBO nuclear magnetic shielding analysis has been applied to methane, plumbane, hydrogen iodide, tetracholoplatinate(II), and hexachloroplatinate(IV).     

Accurate computational thermochemistry from explicitly correlated coupled-cluster theory

Explicitly correlated coupled-cluster theory has developed into a valuable computational tool for the calculation of electronic energies close to the limit of a complete basis set of atomic orbitals. In particular at the level of coupled-cluster theory with single and double excitations (CCSD), the space of double excitations is quickly extended towards a complete basis when Slater-type geminals are added to the wave function expansion. The purpose of the present article is to demonstrate the accuracy and efficiency that can be obtained in computational thermochemistry by a CCSD model that uses such Slater-type geminals. This model is denoted as CCSD(F12), where the acronym F12 highlights the fact that the Slater-type geminals are functions f(r ₁₂) of the interelectronic distances r ₁₂ in the system. The performance of explicitly correlated CCSD(F12) coupled-cluster theory is demonstrated by computing the atomization energies of 73 molecules (containing H, C, N, O, and F) with an estimated root-mean-square deviation from the values compiled in the Active Thermochemical Tables of σ = 0.10 kJ/mol per valence electron. To reach this accuracy, not only the frozen-core CCSD basis-set limit but also high-order excitations (connected triple and quadruple excitations), core–valence correlation effects, anharmonic vibrational zero-point energies, and scalar and spin–orbit relativistic effects must be taken into account.