0000000000215241
AUTHOR
Thomas Bondo Pedersen
Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using cholesky decompositions
A new implementation of the approximate coupled cluster singles and doubles CC2 linear response model using Cholesky decomposition of the two-electron integrals is presented. Significantly reducing storage demands and computational effort without sacrificing accuracy compared to the conventional model, the algorithm is well suited for large-scale applications. Extensive basis set convergence studies are presented for the static and frequency-dependent electric dipole polarizability of benzene and C60, and for the optical rotation of CNOFH2 and (−)-trans-cyclooctene (TCO). The origin-dependence of the optical rotation is calculated and shown to persist for CC2 even at basis set convergence. …
The CCSD(T) model with Cholesky decomposition of orbital energy denominators
A new implementation of the coupled cluster singles and doubles with approximate triples correction method [CCSD(T)] using Cholesky decomposition of the orbital energy denominators is described. The new algorithm reduces the scaling of CCSD(T) from N-7 to N-6, where N is the number of orbitals. The Cholesky decomposition is carried out using simple analytical expressions that allow us to evaluate a priori the order in which the decomposition should be carried out and to obtain the relevant parts of the vectors whenever needed in the calculation. Several benchmarks have been carried out comparing the performance of the conventional and Cholesky CCSD(T) implementations. The Cholesky implement…
Carbon nanorings: A challenge to theoretical chemistry
High-level quantum-chemical methods show that the binding in the inclusion complex of hexamethylbenzene (HMB) in 6-cycloparaphenilacetylene (6-CPPA) cannot be explained only in terms of electrostatic interactions - caused by the polarization associated to curved Ï-conjugated systems - and the inclusion of dispersion forces is definitely needed. The theoretical description of van der Waals interactions is notoriously complicated and in fact some DFT methods cannot even predict the existence of the relatively small supramolecular nanoring studied here. However, ab initio MP2 calculations agree with experimental data and show that, in the considered complex, the HMB fragment is placed at the …
Origin invariant calculation of optical rotation without recourse to London orbitals
We propose an inherently origin invariant formulation of specific optical rotation obtained from the conventional velocity gauge expression by subtracting the static limit from the electric dipole-magnetic dipole polarizability: ããpâ,LâããÏâããpâ,LâããÏ-ããpâ,Lâãã0. London (gauge including) atomic orbitals are therefore not needed to ensure origin invariance and, consequently, this modified velocity gauge formulation may be employed in conjunction with variational and non-variational quantum chemical methods alike. Sample calculations on S-propylene oxide and 3,4- methylenedioxymethamphetamine using self-consistent field and coupled cluster wave functions…
The Dalton quantum chemistry program system
Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, MOller-Plesset, confi ...
Study of the benzene⋅N2 intermolecular potential-energy surface
The intermolecular potential-energy surface pertaining to the interaction between benzene and N2 is investigated theoretically and experimentally. Accurate intermolecular interaction energies are evaluated for the benzene–N2 van der Waals complex using the coupled cluster singles and doubles including connected triples [CCSD(T)] method and the aug-cc-pVDZ basis set extended with a set of 3s3p2d1f1g midbond functions. After fitting the energies to an analytic function, the intermolecular Schrödinger equation is solved to yield energies, rotational constants, and Raman-scattering coefficients for the lowest intermolecular levels of several benzene–N2 isotopomers. Experimentally, intermolecula…
Calibration of Cholesky Auxiliary Basis Sets for Multiconfigurational Perturbation Theory Calculations of Excitation Energies
The accuracy of auxiliary basis sets derived from Cholesky decomposition of two-electron integrals is assessed for excitation energies calculated at the state-average complete active space self-consistent field (CASSCF) and multiconfigurational second order perturbation theory (CASPT2) levels of theory using segmented as well as generally contracted atomic orbital basis sets. Based on 196 valence excitations in 26 organic molecules and 72 Rydberg excitations in 3 organic molecules, the results show that Cholesky auxiliary basis sets can be used without compromising the accuracy of the multiconfigurational methods. Specifically, with a decomposition threshold of 10(-4) au, the mean error due…
Fast noniterative orbital localization for large molecules
We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implem…
Cholesky decomposition techniques in electronic structure theory
We review recently developed methods to efficiently utilize the Cholesky decomposition technique in electronic structure calculations. The review starts with a brief introduction to the basics of the Cholesky decomposition technique. Subsequently, examples of applications of the technique to ab inito procedures are presented. The technique is demonstrated to be a special type of a resolution-of-identity or density-fitting scheme. This is followed by explicit examples of the Cholesky techniques used in orbital localization, computation of the exchange contribution to the Fock matrix, in MP2, gradient calculations, and so-called method specific Cholesky decomposition. Subsequently, examples o…
Computational and experimental investigation of intermolecular states and forces in the benzene-helium van der Waals complex
A study of the intermolecular potential-energy surface (IPS) and the intermolecular states of the perprotonated and perdeuterated benzene–He complex is reported. From a fit to ab initio data computed within the coupled cluster singles and doubles including connected triples model for 280 interaction geometries, an analytic IPS including two- to four-body atom–atom terms is obtained. This IPS, and two other Lennard-Jones atom–atom surfaces from the literature, are each employed in dynamically exact (within the rigid-monomer approximation) calculations of J = 0 intermolecular states of the isotopomers. Rotational constants and Raman-scattering coefficients for intermolecular vibrational trans…
Variation of polarizability in the [4n+2] annulene series: from [22]- to [66]-annulene.
Using correlated ab initio methods, the polarizability of large [4n + 2]-annulenes is determined, showing that there exists an almost linear relation between the exaltation of magnetic susceptibility (a measure of aromaticity) and an equivalent enlargement of polarizability.
Reduced scaling in electronic structure calculations using Cholesky decompositions
The small numerical rank of the two-electron integral matrix for large molecular systems and large basis sets was demonstrated. Though, the current implementation still requires some improvements on the calculations done in the inner most loop of the decomposition do not exploit the parsity in the Cholesky vectors. With respect to the practical applicability of the presented method an efficient approach to geometrical derivatives was imperative. Such an approach was obtained including certain derivative product functions and decomposing an expanded integral matrix.
Rovibrational structure of the Ar–CO complex based on a novel three-dimensional ab initio potential
The first three-dimensional ab initio intermolecular potential energy surface of the Ar–CO van der Waals complex is calculated using the coupled cluster singles and doubles including connected triples model and the augmented correlation-consistent polarized valence quadruple zeta (aug-cc-pVQZ) basis set extended with a (3s3p2d1f1g) set of midbond functions. The three-dimensional surface is averaged over the three lowest vibrational states of CO. Rovibrational energies are calculated up to 50 cm−1 above the ground state, thus enabling comprehensive comparison between theory and available experimental data as well as providing detailed guidance for future spectroscopic investigations of highe…
Accurate ab initio density fitting for multiconfigurational self-consistent field methods
Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra- μ -acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxi…
Theoretical absorption spectrum of the Ar–CO van der Waals complex
The three-dimensional intermolecular electric dipole moment surface of Ar–CO is calculated at the coupled cluster singles and doubles level of theory with the aug-cc-pVTZ basis set extended with a 3s3p2d1f1g set of midbond functions. Using the rovibrational energies and wave functions of our recent study [J. Chem. Phys. 117, 6562 (2002)], temperature-dependent spectral intensities are evaluated and compared to available experimental data. Based on the theoretical spectrum, alternative assignments of the experimentally observed lines in the fundamental band of CO around 2160 and 2166 cm−1 are suggested. Thomas.Bondo@uv.es
Coupled cluster calculations of interaction energies in benzene–fluorobenzene van der Waals complexes
Benzene-fluorobenzene complexes are used as model systems to simulate the interactions of the SBB-HCAII protein-ligand complex. Using the second-order Moller-Plesset [MP2] and the coupled cluster singles and doubles including connected triple excitations models recently implemented with Cholesky decompositions we evaluate accurate interaction energies for several benzene-fluorobenzene van der Waals complexes. We consider edge-to-face interactions and compare the results to those from a recent MP2 study and to experimental findings. In contrast to experimental trends, we find that the interaction tends to decrease with increasing fluorination and conclude that benzene-fluorobenzene complexes…