Search

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H2Z2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces. [ABSTRACT FROM AUTHOR]

The superposition of atomic potentials (SAP) approach has recently been shown to be a simple and efficient way to initialize electronic structure calculations [S. Lehtola, J. Chem. Theory Comput. 15, 1593–1604 (2019)]. Here, we study the differences between effective potentials from fully numerical density functional and optimized effective potential calculations for fixed configurations. We find that the differences are small, overall, and choose exchange-only potentials at the local density approximation level of theory computed on top of Hartree–Fock densities as a good compromise. The differences between potentials arising from different atomic configurations are also found to be small at this level of theory. Furthermore, we discuss the efficient Gaussian-basis implementation of SAP via error function fits to fully numerical atomic radial potentials. The guess obtained from the fitted potentials can be easily implemented in any Gaussian-basis quantum chemistry code in terms of two-electron integrals. Fits covering the whole periodic table from H to Og are reported for non-relativistic as well as fully relativistic four-component calculations that have been carried out with fully numerical approaches. [ABSTRACT FROM AUTHOR]

We present a formulation and implementation of the calculation of (orbital-unrelaxed) expectation values at the 4-component relativistic coupled cluster level with spin-orbit coupling included from the start. The Lagrangian-based analytical energy derivative technique constitutes the basic theoretical framework of this work. The key algorithms for single reference relativistic coupled cluster have been implemented using routines for general tensor contractions of up to rank-2 tensors in which the direct product decomposition scheme is employed to benefit from double group symmetry. As a sample application, we study the electric field gradient at the bismuth nucleus in the BiX (X = N, P) series of molecules, where the effect of spin-orbit coupling is substantial. Our results clearly indicate that the current reference value for the nuclear quadrupole moment of 209Bi needs revision. We also have applied our method to the calculation of the parity violating energy shift of chiral molecules. The latter property is strictly zero in the absence of spin-orbit coupling. For the H2X2 (X = O,S,Se,Te) series of molecules the effect of correlation is found to be quite small. [ABSTRACT FROM AUTHOR]

Localized molecular orbitals are often used for the analysis of chemical bonds, but they can also serve to efficiently and comprehensibly compute linear response properties. While conventional canonical molecular orbitals provide an adequate basis for the treatment of excited states, a chemically meaningful identification of the different excited-state processes is difficult within such a delocalized orbital basis. In this work, starting from an initial set of supermolecular canonical molecular orbitals, we provide a simple one-step top-down embedding procedure for generating a set of orbitals, which are localized in terms of the supermolecule but delocalized over each subsystem composing the supermolecule. Using an orbital partitioning scheme based on such sets of localized orbitals, we further present a procedure for the construction of local excitations and charge-transfer states within the linear response framework of time-dependent density functional theory (TDDFT). This procedure provides direct access to approximate diabatic excitation energies and, under the Tamm–Dancoff approximation, also their corresponding electronic couplings—quantities that are of primary importance in modeling energy transfer processes in complex biological systems. Our approach is compared with a recently developed diabatization procedure based on subsystem TDDFT using projection operators, which leads to a similar set of working equations. Although both of these methods differ in the general localization strategies adopted and the type of basis functions (Slaters vs Gaussians) employed, an overall decent agreement is obtained. [ABSTRACT FROM AUTHOR]

Recently, a general framework suitable for general frozen-density embedding (FDE) methods was published [S. Höfener, A. S. P. Gomes, and L. Visscher, J. Chem. Phys. 136, 044104 (2012)]. In the present article, we report the fragmentation of a supermolecule while treating all subsystems with coupled-cluster theory and the interaction of the subsystems with density-functional theory. This variant is denoted wave-function theory in wave-function theory FDE, or coupled-cluster theory in coupled-cluster theory FDE. Main target of this approach is not the embedding of a single molecule in large solvation shells, but rather the possibility to divide a complex system consisting of several molecules when all subsystems are to be treated with, e.g., coupled-cluster methods to provide a balanced and unbiased description. We present numerical results for hydrogen-bonded complexes which exhibit rather strong interactions. Cases with weakly interacting subsystems are expected to exhibit even higher accuracy. This facilitates the study of properties of larger complexes such as DNA base pairs with coupled-cluster methods. [ABSTRACT FROM AUTHOR]

The triiodide ion I3- exhibits a complex photodissociation behavior, the dynamics of which are not yet fully understood. As a first step toward determining the full potential energy surfaces of this species for subsequent simulations of its dissociation processes, we investigate the performance of different electronic structure methods [time-dependent density functional theory, complete active space perturbation theory to second order (CASPT2), Fock-space coupled cluster and multireference configuration interaction] in describing the ground and excited states of the triiodide ion along the symmetrical dissociation path. All methods apart from CASPT2 include scalar relativity and spin-orbit coupling in the orbital optimization, providing useful benchmark data for the more common two-step approaches in which spin-orbit coupling is introduced in the configuration interaction. Time-dependent density functional theory with the statistical averaging of model orbital potential functional is off the mark for this system. Another choice of functional may improve performance with respect to vertical excitation energies and spectroscopic constants, but all functionals are likely to face instability problems away from the equilibrium region. The Fock-space coupled cluster method was shown to perform clearly best in regions not too far from equilibrium but is plagued by convergence problems toward the dissociation limit due to intruder states. CASPT2 shows good performance at significantly lower computational cost, but is quite sensitive to symmetry breaking. We furthermore observe spikes in the CASPT2 potential curves away from equilibrium, signaling intruder state problems that we were unable to curb through the use of level shifts. Multireference configuration interaction is, in principle, a viable option, but its computational cost in the present case prohibits use other than for benchmarking purposes. [ABSTRACT FROM AUTHOR]

A new approach for relativistic correlated electron structure calculations is proposed by which a transformation to a two-spinor basis is carried out after solving the four-component relativistic Hartree–Fock equations. The method is shown to be more accurate than approaches that apply an a priori transformation to a two-spinor basis. We also demonstrate how the two-component relativistic calculations with properly transformed two-electron interaction can be simulated at the four-component level by projection techniques, thus allowing an assessment of errors introduced by more approximate schemes. [ABSTRACT FROM AUTHOR]

We present an extension of the frozen-density embedding (FDE) scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] that can be applied to subsystems connected by covalent bonds, as well as a practical implementation of such an extended FDE scheme. We show how the proposed scheme can be employed for quantum chemical calculations of proteins by treating each constituting amino acid as a separate subsystem. To assess the accuracy of the extended FDE scheme, we present calculations for several dipeptides and for the protein ubiquitin. [ABSTRACT FROM AUTHOR]

Relativistic coupled cluster theory is used to determine accurate electric field gradients in order to provide a theoretical value for the nuclear quadrupole moment of 139La. Here we used the diatomic lanthanum monohalides LaF, LaCl, LaBr, and LaI as accurate nuclear quadrupole coupling constants are available from rotational spectroscopy by Rubinoff et al. [J. Mol. Spectrosc. 218, 169 (2003)]. The resulting nuclear quadrupole moment for 139La (0.200±0.006 barn) is in excellent agreement with earlier work using atomic hyperfine spectroscopy [0.20(1) barn]. [ABSTRACT FROM AUTHOR]

We have extended the frozen-density embedding (FDE) scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] to include external magnetic fields and applied this extension to the nonrelativistic calculation of nuclear magnetic resonance (NMR) shieldings. This leads to a formulation in which the electron density and the induced current are calculated separately for the individual subsystems. If the current dependence of the exchange-correlation functional and of the nonadditive kinetic-energy functional are neglected, the induced currents in the subsystems are not coupled and each of them can be determined without knowledge of the induced current in the other subsystem. This allows the calculation of the NMR shielding as a sum of contributions of the individual subsystems. As a test application, we have calculated the solvent shifts of the nitrogen shielding of acetonitrile for different solvents using small geometry-optimized clusters consisting of acetonitrile and one solvent molecule. By comparing to the solvent shifts obtained from supermolecular calculations we assess the accuracy of the solvent shifts obtained from FDE calculations. We find a good agreement between supermolecular and FDE calculations for different solvents. In most cases it is possible to neglect the contribution of the induced current in the solvent subsystem to the NMR shielding, but it has to be considered for aromatic solvents. We demonstrate that FDE can describe the effect of induced currents in the environment accurately. [ABSTRACT FROM AUTHOR]

Improvements on current one-component extraction procedures of spin-orbit pseudopotentials are investigated for high accuracy computation of spin-orbit coupling energies. By means of the perturbation-theory formalism we first show that spin-orbit pseudopotentials, extracted at the one-component self-consistent-field level from a reference all-electron Dirac-Coulomb or Dirac-Coulomb-Breit calculation, include valence spin-orbit polarization and relaxation effects. As a consequence the use of these pseudopotentials in uncontracted spin-orbit configuration interaction (CI) with singles from the reference ground-state configuration gives rise to double counting of these spin-orbit effects. Two new methods that avoid such double counting have been investigated. The first, so-called “explicit” method, calculates explicitly, by means of a four-component spin-orbit CI, the double-counted spin-orbit effects and removes them from the pseudopotentials. Due to the nonadditivity of the core and valence spin-orbit effects as well as the so-called “pseudovariational collapse,” this method is shown to be cumbersome. In the second “implicit” method the spin-orbit pseudopotential is extracted at the spin-orbit polarized and relaxed level by means of a single-excitation spin-orbit CI calculation. Atomic tests on iodine demonstrate the ability of the latter method to solve the double-counting problem. [ABSTRACT FROM AUTHOR]

We present calculations at the relativistic coupled cluster theory that predict the 1Σ0+ ground state of CUO to lie 58.2 kJ/mol below the first excited state, 3[uppercase_phi_synonym]2. This can be contrasted with the outcome of earlier density functional theory and complete active space second order perturbation theory (CASPT2) calculations that both predicted a 3[uppercase_phi_synonym]2 ground state upon inclusion of spin-orbit coupling in the calculations. Our result gives further justification to the interpretation of the measured frequency shifts of this species in various noble gas matrices as being caused by significant interaction between the uranium and the heavier noble gas atoms. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The nuclear quadrupole moment of [sup 115]In has been determined by combining the experimental nuclear quadrupole coupling constants and electric field gradients, calculated at the four-component CCSD(T) level of theory, of four indium halides. Our recommended value for the nuclear quadrupole moment of [sup 115]In is 770(8) millibarn. A basis set study at the Dirac-Hartree-Fock level shows a slow convergence of the electric field gradient with respect to higher angular momentum functions. [ABSTRACT FROM AUTHOR]

An implementation of the relativistic multireference Fock-space coupled cluster method is presented which allows simultaneous calculation of potential surfaces for different oxidation states and electronic levels of a molecule, yielding values for spectroscopic constants and transition energies. The method is tested in pilot calculations on the I[sub 2] and HgH molecules, and is shown to give a good and balanced description of various electronic states and energies. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Resonance Raman spectroscopy has long been established as one of the most sensitive techniques for detection, structure characterization, and probing the excited-state dynamics of biochemical systems. However, the analysis of resonance Raman spectra is much facilitated when measurements are accompanied by Density Functional Theory (DFT) calculations that are expensive for large biomolecules. In this work, resonance Raman spectra are therefore computed with the Density Functional Tight-Binding (DFTB) method in the time-dependent excited-state gradient approximation. To test the accuracy of the tight-binding approximations, this method is first applied to typical resonance Raman benchmark molecules, such as β-carotene, and compared to results obtained with pure and range-separated exchange–correlation functionals. We then demonstrate the efficiency of the approach by considering a computationally challenging heme variation. Overall, we find that the vibrational frequencies and excited-state properties (energies and gradients) that are needed to simulate the spectra are reasonably accurate and suitable for interpretation of experiments. We can therefore recommend DFTB as a fast computational method to interpret resonance Raman spectra. [ABSTRACT FROM AUTHOR]

Examines molecular relativistic effects on calculated electric field gradient (EFG) of the hydrogen halides using relativistic four-component method. Importance of spin-orbit coupling and picture change effect; Analysis of nuclear magnetic resonance of nuclei; Observation on the coupling of the nuclear electric quadruple moment with EFG.

The formalism for a relativistic open-shell CCSD(T) method is presented and implemented in a computer program, RELCCSD. The code can be used for calculations with 2- or 4-component relativistic reference wave functions and allows a full inclusion of the spin–orbit coupling. The code is interfaced to the MOLFDIR program system. We illustrate its use with ab initio calculations of the fine structure splittings of Cl, FO, ClO, O+2, and O-2. The triples correction is found to make a large contribution to the Cl atom splitting, which is within 23 cm-1, of the experimental value. The molecular results are within 4 cm-1 of the experimental values where these are available. The value for FO is predicted to be -195±4 cm-1, in good agreement with experiment. © 1996 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The NL-eEDM collaboration is building an experimental setup to search for the permanent electric dipole moment of the electron in a slow beam of cold barium fluoride molecules [NL-eEDM Collaboration, Eur. Phys. J. D 72, 197 (2018)]. Knowledge of the molecular properties of BaF is thus needed to plan the measurements and, in particular, to determine the optimal laser-cooling scheme. Accurate and reliable theoretical predictions of these properties require the incorporation of both high-order correlation and relativistic effects in the calculations. In this work, theoretical investigations of the ground and lowest excited states of BaF and its lighter homologs, CaF and SrF, are carried out in the framework of the relativistic Fock-space coupled cluster and multireference configuration interaction methods. Using the calculated molecular properties, we determine the Franck-Condon factors (FCFs) for the A 2 Π 1 / 2 → X 2 Σ 1 / 2 + transition, which was successfully used for cooling CaF and SrF and is now considered for BaF. For all three species, the FCFs are found to be highly diagonal. Calculations are also performed for the B 2 Σ 1 / 2 + → X 2 Σ 1 / 2 + transition recently exploited for laser-cooling of CaF; it is shown that this transition is not suitable for laser-cooling of BaF, due to the nondiagonal nature of the FCFs in this system. Special attention is given to the properties of the A′2Δ state, which in the case of BaF causes a leak channel, in contrast to CaF and SrF species where this state is energetically above the excited states used in laser-cooling. We also present the dipole moments of the ground and excited states of the three molecules and the transition dipole moments (TDMs) between the different states. Finally, using the calculated FCFs and TDMs, we determine that the A 2 Π 1 / 2 → X 2 Σ 1 / 2 + transition is suitable for transverse cooling in BaF. [ABSTRACT FROM AUTHOR]

We report in this paper an implementation of a 4-component relativistic Hamiltonian based Equation-of-Motion Coupled-Cluster with singles and doubles (EOM-CCSD) theory for the calculation of ionization potential, electron affinity, and excitation energy. In this work, we utilize the previously developed double group symmetry-based generalized tensor contraction scheme and also extend it in order to carry out tensor contractions involving non-totally symmetric and odd-ranked tensors. Several approximated spin-free and two-component Hamiltonians can also be accessed in this implementation. We have applied this method to the halogen monoxide (XO, X = Cl, Br, I, At, Ts) species, in order to assess the quality of a few other recent EOM-CCSD implementations, where spin-orbit coupling contribution has been approximated in different degrees. Besides, we have also studied various excited states of CH2IBr, CH2I2, and I 3 − (as well as single electron attachment and detachment electronic states of the same species) where comparison has been made with a closely related multi-reference coupled-cluster method, namely, Intermediate Hamiltonian Fock Space Coupled-Cluster singles and doubles theory. [ABSTRACT FROM AUTHOR]

We present a formulation of Laplace-transformed atomic orbital-based second-order Møller-Plesset perturbation theory (MP2) energies for two-component Hamiltonians in the Kramers-restricted formalism. This low-order scaling technique can be used to enable correlated relativistic calculations for large molecular systems. We show that the working equations to compute the relativistic MP2 energy differ by merely a change of algebra (quaternion instead of real) from their non-relativistic counterparts. With a proof-of-principle implementation we study the effect of the nuclear charge on the magnitude of half-transformed integrals and show that for light elements spin-free and spin-orbit MP2 energies are almost identical. Furthermore, we investigate the effect of separation of charge distributions on the Coulomb and exchange energy contributions, which show the same long-range decay with the inter-electronic/atomic distance as for nonrelativistic MP2. A linearly scaling implementation is possible if the proper distance behavior is introduced to the quaternion Schwarz-type estimates as for non-relativistic MP2. [ABSTRACT FROM AUTHOR]

In this work we investigate different approaches for calculating electric field gradients in order to provide accurate theoretical values for the nuclear quadrupole moments (NQM) for aluminum and gallium. Electron correlation is included in a fully four-component framework at the CCSD(T) level. The resulting NQM for [sup 27]Al (146.0±0.4 mb) is in good agreement with earlier work, while the value for [sup 69]Ga (171±2 mb) is higher than suggested on basis of previous molecular calculations. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

New implementations of the Lévy-Leblond, zeroth-order regular approach (ZORA) and spin-free Dirac equation are presented within the framework of the four-component relativistic program system DIRAC. This implementation allows systematic incorporation of relativistic effects at different levels of theory and corresponding computational cost. One of the possibilities of the new code is to neglect the effect of spin-orbit coupling in the orbital optimization process and introduce it in a later stage of the calculation. This method is shown to be unstable despite the boundedness of the spin-orbit operator itself. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Relativistic four-component random phase approximation (RPA) calculations of indirect nuclear spin-spin coupling constants in MH[sub 4] (M=C, Si, Ge, Sn, Pb) and Pb(CH[sub 3])[sub 3]H are presented. The need for tight s-functions also in relativistic four-component calculations is verified and explained, and the effect of omission of (SS-LL) and (SS-SS) two-electron integrals is investigated. Already in GeH[sub 4] we see a relativistic increase in the coupling constant by 12%, and for PbH[sub 4] the effect is a 156% increase for the one-bond coupling. Large relativistic effects are also computed for the two-bonds couplings. We find that the relativistic effects on the one-bond couplings are mainly due to scalar relativistic factors rather than spin-orbit corrections. © 2000 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In the new field of quantum plasmonics, plasmonic excitations of silver and gold nanoparticles are utilized to manipulate and control light–matter interactions at the nanoscale. While quantum plasmons can be described with atomistic detail using Time-Dependent Density Functional Theory (DFT), such studies are computationally challenging due to the size of the nanoparticles. An efficient alternative is to employ DFT without approximations only for the relatively fast ground state calculations and use tight-binding approximations in the demanding linear response calculations. In this work, we use this approach to investigate the nature of plasmonic excitations under the variation of the separation distance between two nanoparticles. We thereby provide complementary characterizations of these excitations in terms of Kohn–Sham single–orbital transitions, intrinsic localized molecular fragment orbitals, scaling of the electron–electron interactions, and probability of electron tunneling between monomers. [ABSTRACT FROM AUTHOR]

We present a method for a subsystem-based calculation of indirect nuclear spin-spin coupling tensors within the framework of current-spin-density-functional theory. Our approach is based on the frozen-density embedding scheme within density-functional theory and extends a previously reported subsystem-based approach for the calculation of nuclear magnetic resonance shielding tensors to magnetic fields which couple not only to orbital but also spin degrees of freedom. This leads to a formulation in which the electron density, the induced paramagnetic current, and the induced spin-magnetization density are calculated separately for the individual subsystems. This is particularly useful for the inclusion of environmental effects in the calculation of nuclear spin-spin coupling constants. Neglecting the induced paramagnetic current and spin-magnetization density in the environment due to the magnetic moments of the coupled nuclei leads to a very efficient method in which the computationally expensive response calculation has to be performed only for the subsystem of interest. We show that this approach leads to very good results for the calculation of solvent-induced shifts of nuclear spin-spin coupling constants in hydrogen-bonded systems. Also for systems with stronger interactions, frozen-density embedding performs remarkably well, given the approximate nature of currently available functionals for the non-additive kinetic energy. As an example we show results for methylmercury halides which exhibit an exceptionally large shift of the one-bond coupling constants between 199 Hg and 13C upon coordination of dimethylsulfoxide solvent molecules. [ABSTRACT FROM AUTHOR]

Building on the framework recently reported for determining general response properties for frozen-density embedding [S. Höfener, A. S. P. Gomes, and L. Visscher, J. Chem. Phys. 136, 044104 (2012)], in this work we report a first implementation of an embedded coupled-cluster in density-functional theory (CC-in-DFT) scheme for electronic excitations, where only the response of the active subsystem is taken into account. The formalism is applied to the calculation of coupled-cluster excitation energies of water and uracil in aqueous solution. We find that the CC-in-DFT results are in good agreement with reference calculations and experimental results. The accuracy of calculations is mainly sensitive to factors influencing the correlation treatment (basis set quality, truncation of the cluster operator) and to the embedding treatment of the ground-state (choice of density functionals). This allows for efficient approximations at the excited state calculation step without compromising the accuracy. This approximate scheme makes it possible to use a first principles approach to investigate environment effects with specific interactions at coupled-cluster level of theory at a cost comparable to that of calculations of the individual subsystems in vacuum. [ABSTRACT FROM AUTHOR]

Quantum-mechanical methods that are both computationally fast and accurate are not yet available for electronic excitations having charge transfer character. In this work, we present a significant step forward towards this goal for those charge transfer excitations that take place between non-covalently bound molecules. In particular, we present a method that scales linearly with the number of non-covalently bound molecules in the system and is based on a two-pronged approach: The molecular electronic structure of broken-symmetry charge-localized states is obtained with the frozen density embedding formulation of subsystem density-functional theory; subsequently, in a post-SCF calculation, the full-electron Hamiltonian and overlap matrix elements among the charge-localized states are evaluated with an algorithm which takes full advantage of the subsystem DFT density partitioning technique. The method is benchmarked against coupled-cluster calculations and achieves chemical accuracy for the systems considered for intermolecular separations ranging from hydrogen-bond distances to tens of Ångstroms. Numerical examples are provided for molecular clusters comprised of up to 56 non-covalently bound molecules. [ABSTRACT FROM AUTHOR]

In this article, we present a consistent derivation of a density functional theory (DFT) based embedding method which encompasses wave-function theory-in-DFT (WFT-in-DFT) and the DFT-based subsystem formulation of response theory (DFT-in-DFT) by Neugebauer [J. Neugebauer, J. Chem. Phys. 131, 084104 (2009)] as special cases. This formulation, which is based on the time-averaged quasi-energy formalism, makes use of the variation Lagrangian techniques to allow the use of non-variational (in particular: coupled cluster) wave-function-based methods. We show how, in the time-independent limit, we naturally obtain expressions for the ground-state DFT-in-DFT and WFT-in-DFT embedding via a local potential. We furthermore provide working equations for the special case in which coupled cluster theory is used to obtain the density and excitation energies of the active subsystem. A sample application is given to demonstrate the method. [ABSTRACT FROM AUTHOR]

We report the development of a general order relativistic coupled-cluster (CC) code. Our implementation is based on Kramers-paired molecular spinors, utilizes double group symmetry, and is applicable with the full Dirac-Coulomb and several approximate relativistic Hamiltonians. The available methods include iterative and perturbative single-reference CC approaches with arbitrary excitations as well as a state-selective multi-reference CC ansatz. To illustrate the performance of the new code, benchmark calculations have been performed for the total energies, bond lengths, and vibrational frequencies of the monoxides of Group IVa elements. The trends due to the simultaneous inclusion of relativity as well as higher-order electron correlation effects are analyzed. The newly developed code significantly widens the scope of the ab initio relativistic calculations, for both molecules and atoms alike, surpassing the accuracy and reliability of the currently available implementations in the literature. [ABSTRACT FROM AUTHOR]

We have investigated the functional derivative of the nonadditive kinetic-energy bifunctional, which appears in the embedding potential that is used in the frozen-density embedding formalism, in the limit that the separation of the subsystems is large. We have derived an exact expression for this kinetic-energy component of the embedding potential and have applied this expression to deduce its exact form in this limit. Comparing to the approximations currently in use, we find that while these approximations are correct at the nonfrozen subsystem, they fail completely at the frozen subsystem. Using test calculations on two model systems, a H2O...Li+ complex and a cluster of aminocoumarin C151 surrounded by 30 water molecules, we show that this failure leads to a wrong description of unoccupied orbitals, which can lead to convergence problems caused by too low-lying unoccupied orbitals and which can further have serious consequences for the calculation of response properties. Based on our results, a simple correction is proposed, and we show that this correction is able to fix the observed problems for the model systems studied. [ABSTRACT FROM AUTHOR]

Electric field gradients at the nuclei of halogen atoms are calculated using a finite field approach. The four-component Dirac-Coulomb Hamiltonian serves as the framework, all electrons are correlated by the relativistic Fock-space coupled cluster method with single and double excitations, and the Gaunt term, the main part of the Breit interaction, is included. Large basis sets (e.g., 28s24p21d9f4g2h Gaussian-type functions for I) are used. Combined with experimental nuclear quadrupole coupling constants, accurate estimates of the nuclear quadrupole moments are obtained. The calculated values are Q(35Cl)=-81.1(1.2) mb, Q(79Br)=302(5) mb, and Q(127I)=-680(10) mb. Currently accepted reference values [Pyykkö, Mol. Phys. 99, 1617 (2001)] are -81.65(80), 313(3), and -710(10) mb, respectively. Our values are lower for the heavier halogens, corroborating the recent work of van Stralen and Visscher [Mol. Phys. 101, 2115 (2003)], who obtained Q(127I)=-696(12) mb in a series of molecular calculations. [ABSTRACT FROM AUTHOR]

In this paper we explore the use of the novel relativistic intermediate Hamiltonian Fock-space coupled-cluster method in the calculation of the electronic spectrum for small actinyl ions (NpO2+, NpO22+, and PuO22+). It is established that the method, in combination with uncontracted double-zeta quality basis sets, yields excitation energies in good agreement with experimental values, and better than those obtained previously with other theoretical methods. We propose the reassignment of some of the peaks that were observed experimentally, and confirm other assignments. [ABSTRACT FROM AUTHOR]

Relativistic Dirac-Coulomb (DC) Hartree-Fock calculations are employed to obtain the analytic electric field gradient (EFG) on the antimony nucleus in the SbN, SbP, SbF, and SbCl molecules. The electronic correlation contribution to the EFGs is included with the DC-CCSD(T) and DC-CCSD-T approaches, also in the four-component framework, using a finite-difference method. The total EFG results, along with the experimental nuclear quadrupole coupling constants from microwave spectroscopy, allow to derive the nuclear quadrupole moments of 121Sb and 123Sb, respectively, as -543(11) and -692(14) mb. [ABSTRACT FROM AUTHOR]

The orbital-free frozen-density embedding scheme within density-functional theory [T. A. Wesolowski and A. Warshel, J. Phys. Chem. 97, 8050 (1993)] is applied to the calculation of induced dipole moments of the van der Waals complexes CO2...X (X=He, Ne, Ar, Kr, Xe, Hg). The accuracy of the embedding scheme is investigated by comparing to the results of supermolecule Kohn-Sham density-functional theory calculations. The influence of the basis set and the consequences of using orbital-dependent approximations to the exchange-correlation potential in embedding calculations are examined. It is found that in supermolecular Kohn-Sham density-functional calculations, different common approximations to the exchange-correlation potential are not able to describe the induced dipole moments correctly and the reasons for this failure are analyzed. It is shown that the orbital-free embedding scheme is a useful tool for applying different approximations to the exchange-correlation potential in different subsystems and that a physically guided choice of approximations for the different subsystems improves the calculated dipole moments significantly. [ABSTRACT FROM AUTHOR]

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1–0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal–d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The extension of a relativistic double group configuration interaction (CI) formalism to the use of 2- and 4-spinors is presented. We first elucidate the theoretical aspects of the formalism that is needed to work with spinors that are optimized with a Hartree-Fock scheme that includes spin-orbit coupling. We then describe a new general implementation for the computation of sigma vectors and n-particle density matrices that occur in direct CI algorithms. Sample calculations of the spin-orbit splitting in atoms with one particle or hole in an otherwise closed shell configuration ℓ[SUP1] (ℓ = 1,2,3) and molecules containing such atoms illustrate the advantage of treating this effect in a basis of true spinors rather than in a basis of scalar relativistic orbitals as is conventionally done. [ABSTRACT FROM AUTHOR]

An implementation of a two-component all-electron (jj/ωω) treatment of both scalar and spin–orbit relativistic effects in the MOLFDIR program suite is presented. Relativity is accounted for by Douglas–Kroll transformed one-electron operators: scalar (spin-free) and so called mean-field spin–orbit terms. The interelectronic interaction is represented by the nonrelativistic Coulomb operator. High-level correlated calculations of properties of several systems (FO, ClO, Cl, O[sub 2][sup +], O[sub 2][sup -], Tl, and TlH) where spin–orbit effects play a dominant role are presented and compared with other data. Agreement with Dirac–Coulomb(–Gaunt) reference values is in general very good. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree–Fock, Kohn–Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model. [ABSTRACT FROM AUTHOR]

We report a time-dependent density functional based tight-binding (TD-DFTB) scheme for the calculation of UV/Vis spectra, explicitly taking into account the excitation of nuclear vibrations via the adiabatic Hessian Franck-Condon method with a harmonic approximation for the nuclear wavefunction. The theory of vibrationally resolved UV/Vis spectroscopy is first summarized from the viewpoint of TD-DFTB. The method is benchmarked against time-dependent density functional theory (TD-DFT) calculations for strongly dipole allowed excitations in various aromatic and polar molecules. Using the recent 3ob:freq parameter set of Elstner's group, very good agreement with TD-DFT calculations using local functionals was achieved. [ABSTRACT FROM AUTHOR]

A counterintuitive gauge-dependence of NMR shieldings for rare-gas dimers is presented and analyzed. It is shown that common belief about the existence of a natural gauge-origin for spherical atoms with respect to NMR shielding calculations is wrong. [ABSTRACT FROM AUTHOR]

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB. [ABSTRACT FROM AUTHOR]

The potential energy curve for the ground-state of radium dimer (Ra2) is provided by means of atomic and molecular relativistic coupled cluster calculations. The short-range part of this curve is defined by an equilibrium bond length of 5.324 Å, a dissociation energy of 897 cm-1, and a harmonic vibrational frequency of 20.5 cm-1. The asymptotic behavior at large interatomic distances is characterized by the van der Waals coefficients C6 = 5.090 × 10³, C8 = 6.978 × 105, and C10 = 8.786 × 107 atomic units. The two regions are matched in an analytical potential to provide a convenient representation for use in further calculations, for instance, to model cold collisions between radium atoms. This might become relevant in future experiments on ultracold, optically trapped, radioactive radium atoms that are used to search for a permanent electric dipole moment. [ABSTRACT FROM AUTHOR]

We present a study of the electronic structure of the [UO2]+, [UO2]2 +, [UO2]3+, NUO, [NUO]+, [NUO]2+, [NUN]-, NUN, and [NUN]+ molecules with the intermediate Hamiltonian Fock-space coupled cluster method. The accuracy of mean-field approaches based on the eXact-2-Component Hamiltonian to incorporate spin-orbit coupling and Gaunt interactions are compared to results obtained with the Dirac-Coulomb Hamiltonian. Furthermore, we assess the reliability of calculations employing approximate density functionals in describing electronic spectra and quantities useful in rationalizing Uranium (VI) species reactivity (hardness, electronegativity, and electrophilicity). [ABSTRACT FROM AUTHOR]

We employed an explicit expression for the dispersion (D) energy in conjunction with Kohn-Sham (KS) density functional theory and frozen-density embedding (FDE) to calculate interaction energies between DNA base pairs and a selected set of amino acid pairs in the hydrophobic core of a small protein Rubredoxin. We use this data to assess the accuracy of an FDE-D approach for the calculation of intermolecular interactions. To better analyze the calculated interaction energies we furthermore propose a new energy decomposition scheme that is similar to the well-known KS bond formation analysis [F. M. Bickelhaupt and E. J. Baerends, Rev. Comput. Chem. 15, 1 (2000)], but differs in the electron densities used to define the bond energy. The individual subsystem electron densities of the FDE approach sum to the total electron density which makes it possible to define bond energies in terms of promotion energies and an explicit interaction energy. We show that for the systems considered only a few freeze-and-thaw cycles suffice to reach convergence in these individual bond energy components, illustrating the potential of FDE-D as an efficient method to calculate intermolecular interactions. [ABSTRACT FROM AUTHOR]

The electronic spectrum of the CUO molecule was investigated with the IHFSCC-SD (intermediate Hamiltonian Fock-space coupled cluster with singles and doubles) method and with TD-DFT (time-dependent density functional theory) employing the PBE and PBE0 exchange-correlation functionals. The importance of both spin-orbit coupling and correlation effects on the low-lying excited-states of this molecule are analyzed and discussed. Noble gas matrix effects on the energy ordering and vibrational frequencies of the lowest electronic states of the CUO molecule were investigated with density functional theory (DFT) and TD-DFT in a supermolecular as well as a frozen density embedding (FDE) subsystem approach. This data is used to test the suitability of the FDE approach to model the influence of different matrices on the vertical electronic transitions of this molecule. The most suitable potential was chosen to perform relativistic wave function theory in density functional theory calculations to study the vertical electronic spectra of the CUO and CUONg4 with the IHFSCC-SD method. [ABSTRACT FROM AUTHOR]

The frozen-density embedding (FDE) scheme [Wesolowski and Warshel, J. Phys. Chem. 97, 8050 (1993)] relies on the use of approximations for the kinetic-energy component vT[ρ1,ρ2] of the embedding potential. While with approximations derived from generalized-gradient approximation kinetic-energy density functional weak interactions between subsystems such as hydrogen bonds can be described rather accurately, these approximations break down for bonds with a covalent character. Thus, to be able to directly apply the FDE scheme to subsystems connected by covalent bonds, improved approximations to vT are needed. As a first step toward this goal, we have implemented a method for the numerical calculation of accurate references for vT. We present accurate embedding potentials for a selected set of model systems, in which the subsystems are connected by hydrogen bonds of various strength (water dimer and F–H–F-), a coordination bond (ammonia borane), and a prototypical covalent bond (ethane). These accurate potentials are analyzed and compared to those obtained from popular kinetic-energy density functionals. [ABSTRACT FROM AUTHOR]

We have assessed the accuracy of a representative set of currently available approximate kinetic-energy functionals used within the frozen-density embedding scheme for the NgAuF (Ng=Ar, Kr, Xe) molecules, which we partitioned into a Ng and a AuF subsystem. Although it is weak, there is a covalent interaction between these subsystems which represents a challenge for this subsystem density functional theory approach. We analyzed the effective-embedding potentials and resulting electron density distributions and provide a quantitative analysis of the latter from dipole moment differences and root-mean-square errors in the density with respect to the supermolecular Kohn–Sham density functional theory reference calculation. Our results lead to the conclusion that none of the tested approximate kinetic-energy functionals performs well enough to describe the bond between the noble gas and gold adequately. This observation contributes to the growing evidence that the current procedure to obtain approximate kinetic-energy functionals by reparametrizing functionals obtained via the “conjointness” hypothesis of Lee, Lee, and Parr [Phys. Rev. A 44, 768 (1991)] is insufficient to treat metal-ligand interactions with covalent character. [ABSTRACT FROM AUTHOR]

A valence-only (V) dipole moment surface (DMS) has been computed for water at the internally contracted multireference configuration interaction level using the extended atom-centered correlation-consistent Gaussian basis set aug-cc-pV6Z. Small corrections to these dipole values, resulting from core correlation (C) and relativistic (R) effects, have also been computed and added to the V surface. The resulting DMS surface is hence called CVR. Interestingly, the C and R corrections cancel out each other almost completely over the whole grid of points investigated. The ground-state CVR dipole of H216O is 1.8676 D. This value compares well with the best ab initio one determined in this study, 1.8539±0.0013 D, which in turn agrees well with the measured ground-state dipole moment of water, 1.8546(6) D. Line intensities computed with the help of the CVR DMS shows that the present DMS is highly similar to though slightly more accurate than the best previous DMS of water determined by Schwenke and Partridge [J. Chem. Phys. 113, 16 (2000)]. The influence of the precision of the rovibrational wave functions computed using different potential energy surfaces (PESs) has been investigated and proved to be small, due mostly to the small discrepancies between the best ab initio and empirical PESs of water. Several different measures to test the DMS of water are advanced. The seemingly most sensitive measure is the comparison between the ab initio line intensities and those measured by ultralong pathlength methods which are sensitive to very weak transitions. [ABSTRACT FROM AUTHOR]

The ground and excited states of the UO2 molecule have been studied using a Dirac-Coulomb intermediate Hamiltonian Fock-space coupled cluster approach (DC-IHFSCC). This method is unique in describing dynamic and nondynamic correlation energies at relatively low computational cost. Spin-orbit coupling effects have been fully included by utilizing the four-component Dirac-Coulomb Hamiltonian from the outset. Complementary calculations on the ionized systems UO2+ and UO22+ as well as on the ions U4+ and U5+ were performed to assess the accuracy of this method. The latter calculations improve upon previously published theoretical work. Our calculations confirm the assignment of the ground state of the UO2 molecule as a 3[uppercase_phi_synonym]2u state that arises from the 5f17s1 configuration. The first state from the 5f2 configuration is found above 10 000 cm-1, whereas the first state from the 5f16d1 configuration is found at 5 047 cm-1. [ABSTRACT FROM AUTHOR]

The nuclear quadrupole moment for 197Au has been determined on the base of the state-of-art relativistic molecular calculations. The experimental shifts in the nuclear coupling constants in the series of molecules AuF, XeAuF, KrAuF, ArAuF, (OC)AuF, and AuH have been combined with highly accurate determinations of the electric field gradient (EFG) at the gold nucleus, obtained by molecular relativistic Dirac-Coulomb-Gaunt Hartree-Fock calculations. The electronic correlation contribution to the EFG is included with the CCSD(T) and CCSD-T approaches, also in the four-component framework, using a finite-difference method. In order to estimate the accuracy of their approach the authors have thoroughly investigated the convergence of the results with respect to the basis set employed and the size of the correlated orbital space. The effect of the full Breit electron-electron interaction on the nuclear quadrupole moment of gold has also been considered explicitly for the AuF molecule. They obtain for 197Au a nuclear quadrupole moment of 510±15 mb, which deviates by about 7% from the currently accepted muonic value. [ABSTRACT FROM AUTHOR]