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]

On a comprehensive database with 1,644 datapoints, covering several aspects of main‐group as well as of transition metal chemistry, we assess the performance of 60 density functional approximations (DFA), among them 36 double hybrids (DH). All calculations are performed using a Slater type orbital (STO) basis set of triple‐ζ (TZ) quality and the highly efficient pair atomic resolution of the identity approach for the exchange‐ and Coulomb‐term of the KS matrix (PARI‐K and PARI‐J, respectively) and for the evaluation of the MP2 energy correction (PARI‐MP2). Employing the quadratic scaling SOS‐AO‐PARI‐MP2 algorithm, DHs based on the spin‐opposite‐scaled (SOS) MP2 approximation are benchmarked against a database of large molecules. We evaluate the accuracy of STO/PARI calculations for B3LYP as well as for the DH B2GP‐PLYP and show that the combined basis set and PARI‐error is comparable to the one obtained using the well‐known def2‐TZVPP Gaussian‐type basis set in conjunction with global density fitting. While quadruple‐ζ (QZ) calculations are currently not feasible for PARI‐MP2 due to numerical issues, we show that, on the TZ level, Jacob's ladder for classifying DFAs is reproduced. However, while the best DHs are more accurate than the best hybrids, the improvements are less pronounced than the ones commonly found on the QZ level. For conformers of organic molecules and noncovalent interactions where very high accuracy is required for qualitatively correct results, DHs provide only small improvements over hybrids, while they still excel in thermochemistry, kinetics, transition metal chemistry and the description of strained organic systems. [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]

We present a robust and non-heuristic algorithm that finds all extremum points of the error distribution function of numerically Laplace-transformed orbital energy denominators. The extremum point search is one of the two key steps for finding the minimax approximation. If pre-tabulation of initial guesses is supposed to be avoided, strategies for a sufficiently robust algorithm have not been discussed so far. We compare our non-heuristic approach with a bracketing and bisection algorithm and demonstrate that 3 times less function evaluations are required altogether when applying it to typical non-relativistic and relativistic quantum chemical systems. [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]

Abstract: The first implementation of analytical first-order one-electron molecular properties at the Dirac–Coulomb MP2 level of theory is described. The formalism presented allows for the use of inactive spinors. The method has been applied to the calculation of parity-violation energies in H2X2 (X=O, S, Se, and Te). Electron correlation does not play an important role for this property in this kind of systems, but the effect is shown to be very geometry dependent. [Copyright &y& Elsevier]

Abstract: In this work, we investigate the influence of correlating core electrons on the relativistic corrections to the equilibrium geometries, harmonic frequencies and dissociation energies of hydrogen astatide, using moderately large basis sets. Our results suggest that relativistic effects tend to increase bond length and decrease harmonic frequencies and dissociation energies, in line with previous studies involving valence correlation. In all cases the driving force to the relativistic effects is the spin–orbit coupling, with scalar effects acting cooperatively except for the bond lengths. Core correlation counteracts the relativistic increase of bond lengths, but not strongly enough to reverse any trends. [Copyright &y& Elsevier]

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]

Relativistic ab initio methods are used to compute the electric field gradient at the iodine nucleus in nine different closed-shell diatomic molecules. Combining these theoretical electric field gradients with experimental nuclear quadrupole coupling constants gives a consistent value of the nuclear quadrupole moment of 127 I of -696(12) millibarn. We argue that this value is more precise than the current standard value of the nuclear quadrupole moment of 127 I and recommend adjusting the reference value accordingly. The precision of this determination is still determined by technical limitations in the theoretical work, in particular the neglect of the two-electron Gaunt interaction in the Hamiltonian and correlation contributions beyond those described at the CCSD(T) level of theory, but the errors are reduced relative to the theoretical work that underlies the current standard value of this nuclear quadrupole moment. As a secondary study we also considered the calculation of the small electric field gradient at the gold nucleus in the AuI molecule and conclude that this computation remains a challenge for theoreticians. [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]

Frozen density embedding (FDE) theory is one of the major techniques aiming to bring modeling of extended chemical systems into the realm of high accuracy calculations. To improve its accuracy it is of interest to develop kinetic energy density functional approximations specifically for FDE applications. In the study reported here we focused on optimizing parameters of a generalized gradient approximation‐like kinetic energy functional with the purpose of better describing electron excitation energies. We found that our optimized parametrizations, named excPBE and excPBE‐3 (as these are derived from a Perdew‐Burke‐Ernzerhof‐like parametrization), could not yield improvements over available functionals when applied on a test set of systems designed to probe solvatochromic shifts. Moreover, as several different functionals yielded very similar errors to the simple local‐density approximation (LDA), it is questionable whether it is worthwhile to go beyond the LDA in this context. [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]

Graphical abstract Abstract In the present article, we introduce the relativistic Cholesky-decomposed density (CDD) matrix second-order Møller-Plesset perturbation theory (MP2) energies. The working equations are formulated in terms of the usual intermediates of MP2 when employing the resolution-of-the-identity approximation (RI) for two-electron integrals. Those intermediates are obtained by substituting the occupied and virtual quaternion pseudo-density matrices of our previously proposed two-component (2C) atomic orbital-based MP2 (Helmich-Paris et al., 2016) by the corresponding pivoted quaternion Cholesky factors. While working within the Kramers-restricted formalism, we obtain a formal spin-orbit overhead of 16 and 28 for the Coulomb and exchange contribution to the 2C MP2 correlation energy, respectively, compared to a non-relativistic (NR) spin-free CDD-MP2 implementation. This compact quaternion formulation could also be easily explored in any other algorithm to compute the 2C MP2 energy. The quaternion Cholesky factors become sparse for large molecules and, with a block-wise screening, block sparse-matrix multiplication algorithm, we observed an effective quadratic scaling of the total wall time for heavy-element containing linear molecules with increasing system size. The total run time for both NR and 2C calculations was dominated by the contraction to the exchange energy. We have also investigated a bulky Te-containing supramolecular complex. For such bulky, three-dimensionally extended molecules the present screening scheme has a much larger prefactor and is less effective. [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]

In this article, we document a new implementation of the fuzzy cells scheme for numerical integration in polyatomic systems [Becke, J. Chem. Phys. 1998, 88, 2547] and compare its efficiency and accuracy with respect to an integration scheme based on the Voronoi space partitioning. We show that the accuracy of both methods is comparable, but that the fuzzy cells scheme is better suited for geometry optimization. For this method, we also introduce the locally dense grid concept and present a proof-of-concept application. © 2013 Wiley Periodicals, Inc. [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]

Photocatalytic water oxidation remains the bottleneck in many artificial photosynthesis devices. The efficiency of this challenging process is inherently linked to the thermodynamic and electronic properties of the chromophore and the water oxidation catalyst (WOC). Computational investigations can facilitate the search for favorable chromophore‐catalyst combinations. However, this remains a demanding task due to the requirements on the computational method that should be able to correctly describe different spin and oxidation states of the transition metal, the influence of solvation and the different rates of the charge transfer and water oxidation processes. To determine a suitable method with favorable cost/accuracy ratios, the full catalytic cycle of a molecular ruthenium based WOC is investigated using different computational methods, including density functional theory (DFT) with different functionals (GGA, Hybrid, Double Hybrid) as well as the semi‐empirical tight binding approach GFN‐xTB. A workflow with low computational cost is proposed that combines GFN‐xTB and DFT and provides reliable results. GFN‐xTB geometries and frequencies combined with single‐point DFT energies give free energy changes along the catalytic cycle that closely follow the full DFT results and show satisfactory agreement with experiment, while significantly decreasing the computational cost. This workflow allows for cost efficient determination of energetic, thermodynamic and dynamic properties of WOCs. [ABSTRACT FROM AUTHOR]

Relativistic basis sets of double-zeta, triple-zeta, and quadruple-zeta quality have been optimized for the lanthanide elements La-Lu. The basis sets include SCF exponents for the occupied spinors and for the 6p shell, exponents of correlating functions for the valence shells (4f, 5d and 6s) and the outer core shells (4d, 5s and 5p), and diffuse functions, including functions for dipole polarization of the 4f shell. A finite nuclear size was used in all optimizations. The basis sets are illustrated by calculations on YbF. Prescriptions are given for constructing contracted basis sets. The basis sets are available as an internet archive and from the Dirac program web site, . [ABSTRACT FROM AUTHOR]