Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions

Cholesky Decomposition and the Second-Derivative Two-Electron Integrals Required for the Computation of Magnetizabilities Using Gauge-Including Atomic Orbitals

The computation of magnetizability tensors using gauge-including atomic orbitals is discussed in the context of Cholesky decomposition (CD) for the two-electron repulsion integrals with a focus on the involved doubly differentiated integrals. Three schemes for their handling are suggested: the first exploits the density fitting (DF) aspect of Cholesky decomposition, the second uses expressions obtained by differentiating the CD expression for the unperturbed two-electron integrals, while the third addresses the issue that the first two schemes are not able to represent the doubly differentiated integrals with arbitrary accuracy. This scheme uses a separate Cholesky decomposition for the cross terms in the doubly differentiated two-electron integrals. Test calculations reveal that all three schemes are able to represent the integrals with similar accuracy and yield indistinguishable results for the values of the computed magnetizability tensor elements. Thus, we recommend our first scheme which has the lowest computational cost for routine computations. The applicability of our CD schemes is further shown in large-scale Hartree-Fock calculations of the magnetizability tensor of coronene (C24H12) with a doubly polarized triple-ζ basis consisting of 684 basis functions.

Linear Response pCCD-Based Methods: LR-pCCD and LR-pCCD+S Approaches for the Efficient and Reliable Modeling of Excited State Properties

In this work, we derive working equations for the linear response pair coupled cluster doubles (LR-pCCD) ansatz and its extension to singles (S), LR-pCCD+S. These methods allow us to compute electronic excitation energies and transition dipole moments based on a pCCD reference function. We benchmark the LR-pCCD+S model against the linear response coupled-cluster singles and doubles method for modeling electronic spectra (excitation energies and transition dipole moments) of the BH, H2O, H2CO, and furan molecules. We also analyze the effect of orbital optimization within pCCD on the resulting LR-pCCD+S transition dipole moments and oscillator strengths and perform a statistical error analysis. We show that the LR-pCCD+S method can correctly reproduce the transition dipole moments features, thus representing a reliable and cost-effective alternative to standard, more expensive electronic structure methods for modeling electronic spectra of simple molecules. Specifically, the proposed models require only mean-field-like computational cost, while excited-state properties may approach the CCSD level of accuracy. Moreover, we demonstrate the capability of our model to simulate electronic transitions with non-negligible contributions of double excitations and the electronic spectra of polyenes of various chain lengths, for which standard electronic structure methods perform purely.

Finite-field Cholesky decomposed coupled-cluster techniques (ff-CD-CC): theory and application to pressure broadening of Mg by a He atmosphere and a strong magnetic field

For the interpretation of spectra of magnetic stellar objects such as magnetic white dwarfs (WDs), highly accurate quantum chemical predictions for atoms and molecules in finite magnetic field are required. Especially the accurate description of electronically excited states and their properties requires established methods such as those from coupled-cluster (CC) theory. However, respective calculations are computationally challenging even for medium-sized systems. Cholesky decomposition (CD) techniques may be used to alleviate memory bottlenecks. In finite magnetic field computations, the latter are increased due to the reduction of permutational symmetry within the electron-repulsion-integrals (ERIs) as well as the need for complex-valued data types. CD enables a memory-efficient, approximate description of the ERIs with rigorous error control and thus the treatment of larger systems at the CC level becomes feasible. In order to treat molecules in a finite magnetic field, we present in this work the working equations of the left and right-hand side equations for finite field (ff)-EOM-CD-CCSD for various EOM flavours as well as for the approximate ff-EOM-CD-CC2 method. The methods are applied to the study of the modification of the spectral lines of a magnesium transition by a helium atmosphere that can be found on magnetic WD stars.

Stochastic Resolution of Identity to CC2 for Large Systems: Excited State Properties

We apply a stochastic resolution of identity approximation (sRI) to the CC2 method for the excitation energy calculations. A set of stochastic orbitals are employed to decouple the crucial 4-index electron repulsion integrals and optimize the contraction steps in CC2 response theory. The CC2 response for excitations builds upon sRI-CC2 ground-state calculations, which scales as O(N3), where N is a measure for the system size. Overall, the current algorithm for excited states also allows a sharp scaling reduction from original O(N5) to O(N3). We test the sRI-CC2 for different molecular systems and basis sets, and we show that our sRI-CC2 method can accurately reproduce the results of the deterministic CC2 approach. Our sRI-CC2 exhibits an experimental scaling of O(N2.59) for a series of olefin chains, allowing us to calculate systems with nearly thousands of electrons.

Coupled-cluster treatment of complex open-shell systems: the case of single-molecule magnets

We investigate the reliability of two cost-effective coupled-cluster methods for computing spin-state energetics and spin-related properties of a set of open-shell transition-metal complexes. Specifically, we employ the second-order approximate coupled-cluster singles and doubles (CC2) method and projection-based embedding that combines equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) with density functional theory (DFT). The performance of CC2 and EOM-CCSD-in-DFT is assessed against EOM-CCSD. The chosen test set includes two hexaaqua transition-metal complexes containing Fe(II) and Fe(III), and a large Co(II)-based single-molecule magnet with a non-aufbau ground state. We find that CC2 describes the excited states more accurately, reproducing EOM-CCSD excitation energies within 0.05 eV. However, EOM-CCSD-in-DFT excels in describing transition orbital angular momenta and spin-orbit couplings. Moreover, for the Co(II) molecular magnet, using EOM-CCSD-in-DFT eigenstates and spin-orbit couplings, we compute spin-reversal energy barriers, as well as temperature-dependent and field-dependent magnetizations and magnetic susceptibilities that closely match experimental values within spectroscopic accuracy. These results underscore the efficiency of CC2 in computing state energies of multi-configurational, open-shell systems and highlight the utility of the more cost-efficient EOM-CCSD-in-DFT for computing spin-orbit couplings and magnetic properties of complex and large molecular magnets.

Cholesky Decomposition-Based Implementation of Relativistic Two-Component Coupled-Cluster Methods for Medium-Sized Molecules

A Cholesky decomposition (CD)-based implementation of relativistic two-component coupled-cluster (CC) and equation-of-motion CC (EOM-CC) methods using an exact two-component Hamiltonian augmented with atomic-mean-field spin-orbit integrals (the X2CAMF scheme) is reported. The present CD-based implementation of X2CAMF-CC and EOM-CC methods employs atomic-orbital-based algorithms to avoid the construction of two-electron integrals and intermediates involving three and four virtual indices. Our CD-based implementation extends the applicability of X2CAMF-CC and EOM-CC methods to medium-sized molecules with the possibility to correlate around 1000 spinors. Benchmark calculations for uranium-containing small molecules were performed to assess the dependence of the CC results on the Cholesky threshold. A Cholesky threshold of 10-4 is shown to be sufficient to maintain chemical accuracy. Example calculations to illustrate the capability of the CD-based relativistic CC methods are reported for the bond-dissociation energy of the uranium hexafluoride molecule, UF6, with up to quadruple-ζ basis sets, and the lowest excitation energy in the solvated uranyl ion [UO22+(H2O)12].

On the performance of second-order approximate coupled-cluster singles and doubles methods for non-valence anions

We investigate the capability of several variants of the second-order approximate coupled-cluster singles and doubles (CC2) method to describe dipole-bound, quadrupole-bound, and correlation-bound molecular anions. The binding energy of anions formed by electron attachment to closed-shell molecules is computed using the electron attachment variant of CC2 (EA-CC2), whereas anions with a closed-shell ground state are treated with the standard CC2 method that preserves the number of particles. We find that EA-CC2 captures the binding energies of dipole-bound radical anions quite well, whereas results for other types of non-valence anions are less reliable. We also test the performance of semi-empirical spin-scaling factors for all types of non-valence anions and observe that the spin-scaled CC2 variants generally do not provide more accurate binding energies for dipole-bound anions, while the binding energies of quadrupole-bound and correlation-bound anions are improved. As exemplary applications of EA-CC2, we investigate the dipole-bound anions of the steroids cortisol, progesterone, and testosterone. In addition, we characterize electron attachment to sym-tetracyanonaphthalene, a molecule that supports five anionic states, two of which can be interpreted as hitherto unobserved π-type quadrupole-bound states.

Near-Exact CCSDT Energetics from Rank-Reduced Formalism Supplemented by Non-iterative Corrections

We introduce a non-iterative energy correction, added on top of the rank-reduced coupled-cluster method with single, double, and triple substitutions, that accounts for excitations excluded from the parent triple excitation subspace. The formula for the correction is derived by employing the coupled-cluster Lagrangian formalism, with an additional assumption that the parent excitation subspace is closed under the action of the Fock operator. Owing to the rank-reduced form of the triple excitation amplitudes tensor, the computational cost of evaluating the correction scales as N⁷, where N is the system size. The accuracy and computational efficiency of the proposed method is assessed for both total and relative correlation energies. We show that the non-iterative correction can fulfill two separate roles. If the accuracy level of a fraction of kJ/mol is sufficient for a given system, the correction significantly reduces the dimension of the parent triple excitation subspace needed in the iterative part of the calculations. Simultaneously, it enables reproducing the exact CCSDT results to an accuracy level below 0.1 kJ/mol, with a larger, yet still reasonable, dimension of the parent excitation subspace. This typically can be achieved at a computational cost only several times larger than required for the CCSD(T) method. The proposed method retains the black-box features of the single-reference coupled-cluster theory; the dimension of the parent excitation subspace remains the only additional parameter that has to be specified.

From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape

Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.

Electron confinement meet electron delocalization: non-additivity and finite-size effects in the polarizabilities and dispersion coefficients of the fullerenes

In this work, we used finite-field derivative techniques and density functional theory (DFT) to compute the static isotropic polarizability series (αl with l = 1, 2, 3) for the C60-C84 fullerenes and quantitatively assess the intrinsic non-additivity in these fundamental response properties. By comparing against classical models of the fullerenes as conducting spherical shells (or solid spheres) of uniform electron density, a detailed critical analysis of the derived effective scaling laws (α1 ∼ N1.2, α2 ∼ N2.0, α3 ∼ N2.7) demonstrates that the electronic structure of finite-sized fullerenes-a unique dichotomy of electron confinement and delocalization effects due to their quasi-spherical cage-like structures and encapsulated void spaces-simultaneously limits and enhances their quantum mechanical response to electric field perturbations. Corresponding frequency-dependent polarizabilities were obtained by inputting the αl series into the hollow sphere model (within the modified single frequency approximation), and used to compute the molecular dispersion coefficients (Cn with n = 6, 8, 9, 10) needed to describe the non-trivial van der Waals (vdW) interactions in fullerene-based systems. Using first-order perturbation theory in conjunction with >140 000 DFT calculations, we also computed the non-negligible zero-point vibrational contributions to α1 in C60 and C70, thereby enabling a more accurate and direct comparison between theory and experiment for these quintessential nanostructures.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

A straightforward a posteriori method for reduction of density-fitting error in coupled-cluster calculations

We present a simple method for a posteriori removal of a significant fraction of the density-fitting error from the calculated total coupled-cluster energies. The method treats the difference between the exact and density-fitted integrals as a perturbation, and simplified response-like equations allow us to calculate improved amplitudes and the corresponding energy correction. The proposed method is tested at the coupled-cluster singles and doubles level of theory for a diverse set of moderately-sized molecules. On average, error reductions by a factor of approximately 10 and 20 are observed in double-zeta and triple-zeta basis sets, respectively. Similar reductions are observed in calculations of interaction energies of several model complexes. The computational cost of the procedure is small in comparison with the preceding coupled-cluster iterations. The applicability of this method is not limited to the density-fitting approximation; in principle, it can be used in conjunction with an arbitrary decomposition scheme of the electron repulsion integrals.

Improving the Efficiency of the Multireference Driven Similarity Renormalization Group via Sequential Transformation, Density Fitting, and the Noninteracting Virtual Orbital Approximation

This study examines several techniques to improve the efficiency of the linearized multireference driven similarity renormalization group truncated to one- and two-body operators [MR-LDSRG(2)]. We propose a sequential MR-LDSRG(2) [sq-MR-LDSRG(2)] scheme, in which one-body substitutions are folded exactly into the Hamiltonian. This new approach is combined with density fitting (DF) to reduce the storage cost of two-electron integrals. To further avoid storage of large four-index intermediates, we propose a noninteracting virtual orbital (NIVO) approximation of the Baker-Campbell-Hausdorff series that neglects commutators terms with three and four virtual indices. The NIVO approximation reduces the computational prefactor of the MR-LDSRG(2), bringing it closer to that of coupled cluster with singles and doubles (CCSD). We test the effect of the DF and NIVO approximations on the MR-LDSRG(2) and sq-MR-LDSRG(2) methods by computing properties of eight diatomic molecules. The diatomic constants obtained by DF-sq-MR-LDSRG(2)+NIVO are found to be as accurate as those from the original MR-LDSRG(2) and coupled cluster theory with singles, doubles, and perturbative triples. Finally, we demonstrate that the DF-sq-MR-LDSRG(2)+NIVO scheme can be applied to chemical systems with more than 550 basis functions by computing the automerization energy of cyclobutadiene with a quintuple-ζ basis set. The predicted automerization energy is found to be similar to the value computed with Mukherjee's state-specific multireference coupled cluster theory with singles and doubles.

Reduced-cost linear-response CC2 method based on natural orbitals and natural auxiliary functions

A reduced-cost density fitting (DF) linear-response second-order coupled-cluster (CC2) method has been developed for the evaluation of excitation energies. The method is based on the simultaneous truncation of the molecular orbital (MO) basis and the auxiliary basis set used for the DF approximation. For the reduction of the size of the MO basis, state-specific natural orbitals (NOs) are constructed for each excited state using the average of the second-order Møller-Plesset (MP2) and the corresponding configuration interaction singles with perturbative doubles [CIS(D)] density matrices. After removing the NOs of low occupation number, natural auxiliary functions (NAFs) are constructed [M. Kállay, J. Chem. Phys. 141, 244113 (2014)], and the NAF basis is also truncated. Our results show that, for a triple-zeta basis set, about 60% of the virtual MOs can be dropped, while the size of the fitting basis can be reduced by a factor of five. This results in a dramatic reduction of the computational costs of the solution of the CC2 equations, which are in our approach about as expensive as the evaluation of the MP2 and CIS(D) density matrices. All in all, an average speedup of more than an order of magnitude can be achieved at the expense of a mean absolute error of 0.02 eV in the calculated excitation energies compared to the canonical CC2 results. Our benchmark calculations demonstrate that the new approach enables the efficient computation of CC2 excitation energies for excited states of all types of medium-sized molecules composed of up to 100 atoms with triple-zeta quality basis sets.

Optical Rotation Calculations for a Set of Pyrrole Compounds

Optical rotation of 14 molecules containing the pyrrole group is calculated by employing both time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the second-order approximate coupled-cluster singles and doubles (CC2) method. All optical rotations have been provided using the aug-cc-pVDZ basis set at λ = 589 nm. The two methods predict similar results for both sign and magnitude for the optical rotation of all molecules. The obtained signs are consistent with experiments as well, although several conformers for four molecules needed to be studied to reproduce the experimental sign. We have also calculated excitation energies and rotatory strengths for the six lowest lying electronic transitions for several conformers of the two smallest molecules and found that each rotatory strength has various contributions for each conformer which can cause different optical rotations for different conformers of a molecule. Our results illustrate that both methods are able to reproduce the experimental optical rotations, and that the CAM-B3LYP functional, the least computationally expensive method used here, is an applicable and reliable method to predict the optical rotation for these molecules in line with previous studies.

A Computational Study of the Interaction and Polarization Effects of Complexes Involving Molecular Graphene and C60 or a Nucleobases

A systematic analysis of the molecular structure, energetics, electronic (hyper)polarizabilities and their interaction-induced counterparts of C60 with a series of molecular graphene (MG) models, CmHn, where m = 24, 84, 114, 222, 366, 546 and n = 12, 24, 30, 42, 54, 66, was performed. All the reported data were computed by employing density functional theory and a series of basis sets. The main goal of the study is to investigate how alteration of the size of the MG model affects the strength of the interaction, charge rearrangement, and polarization and interaction-induced polarization of the complex, C60-MG. A Hirshfeld-based scheme has been employed in order to provide information on the intrinsic polarizability density representations of the reported complexes. It was found that the interaction energy increases approaching a limit of -26.98 kcal/mol for m = 366 and 546; the polarizability and second hyperpolarizability increase with increasing the size of MG. An opposite trend was observed for the dipole moment. Interestingly, the variation of the first hyperpolarizability is relatively small with m. Since polarizability is a key factor for the stability of molecular graphene with nucleobases (NB), a study of the magnitude of the interaction-induced polarizability of C84H24-NB complexes is also reported, aiming to reveal changes of its magnitude with the type of NB. The binding strength of C84H24-NB complexes is also computed and found to be in agreement with available theoretical and experimental data. The interaction involved in C60 B12N12H24-NB complexes has also been considered, featuring the effect of contamination on the binding strength between MG and NBs.

Cholesky Decomposition-Based Multiconfiguration Second-Order Perturbation Theory (CD-CASPT2): Application to the Spin-State Energetics of Co(III)(diiminato)(NPh)

The electronic structure and low-lying electronic states of a Co(III)(diiminato)(NPh) complex have been studied using multiconfigurational wave function theory (CASSCF/CASPT2). The results have been compared to those obtained with density functional theory. The best agreement with ab initio results is obtained with a modified B3LYP functional containing a reduced amount (15%) of Hartree-Fock exchange. A relativistic basis set with 869 functions has been employed in the most extensive ab initio calculations, where a Cholesky decomposition technique was used to overcome problems arising from the large size of the two-electron integral matrix. It is shown that this approximation reproduces results obtained with the full integral set to a high accuracy, thus opening the possibility to use this approach to perform multiconfigurational wave-function-based quantum chemistry on much larger systems relative to what has been possible until now.

A systematic way for the cost reduction of density fitting methods

We present a simple approach for the reduction of the size of auxiliary basis sets used in methods exploiting the density fitting (resolution of identity) approximation for electron repulsion integrals. Starting out of the singular value decomposition of three-center two-electron integrals, new auxiliary functions are constructed as linear combinations of the original fitting functions. The new functions, which we term natural auxiliary functions (NAFs), are analogous to the natural orbitals widely used for the cost reduction of correlation methods. The use of the NAF basis enables the systematic truncation of the fitting basis, and thereby potentially the reduction of the computational expenses of the methods, though the scaling with the system size is not altered. The performance of the new approach has been tested for several quantum chemical methods. It is demonstrated that the most pronounced gain in computational efficiency can be expected for iterative models which scale quadratically with the size of the fitting basis set, such as the direct random phase approximation. The approach also has the promise of accelerating local correlation methods, for which the processing of three-center Coulomb integrals is a bottleneck.

Communication: resolving the three-body contribution to the lattice energy of crystalline benzene: benchmark results from coupled-cluster theory

Coupled-cluster theory including single, double, and perturbative triple excitations [CCSD(T)] has been applied to trimers that appear in crystalline benzene in order to resolve discrepancies in the literature about the magnitude of non-additive three-body contributions to the lattice energy. The present results indicate a non-additive three-body contribution of 0.89 kcal mol(-1), or 7.2% of the revised lattice energy of -12.3 kcal mol(-1). For the trimers for which we were able to compute CCSD(T) energies, we obtain a sizeable difference of 0.63 kcal mol(-1) between the CCSD(T) and MP2 three-body contributions to the lattice energy, confirming that three-body dispersion dominates over three-body induction. Taking this difference as an estimate of three-body dispersion for the closer trimers, and adding an Axilrod-Teller-Muto estimate of 0.13 kcal mol(-1) for long-range contributions yields an overall value of 0.76 kcal mol(-1) for three-body dispersion, a significantly smaller value than in several recent studies.

Energy interactions in amyloid-like fibrils from NNQQNY

We use large-scale MP2 calculations to analyze the interactions appearing in amyloid fibers, which are difficult to determine experimentally. To this end, dimers and trimers of the hexapeptide NNQQNY from the yeast prion-like protein Sup35 were considered as model systems. We studied the energy interactions present in the three levels of organization in which the formation of amyloid fibrils is structured. The structural changes in the hydrogen bonds were studied too. It was found that the most energetic process is the formation of the β-sheet, which is equally due to both hydrogen bonds and van der Waals interactions. The aromatic rings help stabilize these aggregates through stacking of the aromatic rings of tyrosine, the stability produced by the aromatics residues increasing with their aromaticity. The formation of the basic unit of the assembled proto-fiber, the steric zipper, is less energetic and is associated to both dispersion forces and hydrogen bonds. The interactions between pair of β-sheets across the peptide-to-peptide contact through the tyrosine rings are cooperative and due to dispersion effects. Moreover, the strength of this interaction can rationalize the variation of mobility of the aromatic ring in the tyrosine units found in solid NMR experiments.

Determination of the absolute configuration of perylene quinone-derived mycotoxins by measurement and calculation of electronic circular dichroism spectra and specific rotations

Altertoxins I-III, alterlosins I and II, alteichin (alterperylenol), stemphyltoxins I-IV, stemphyperylenol, stemphytriol, 7-epi-8-hydroxyaltertoxin I, and 6-epi-stemphytriol are mycotoxins derived from perylene quinone, for which the absolute configuration was not known. Electronic circular dichroism (ECD) spectra were calculated for these compounds and compared with measured spectra of altertoxins I-III, alteichin, and stemphyltoxin III and with reported Cotton effects. Specific rotations were calculated and compared with reported specific rotations. The absolute configuration of all the toxins, except for stemphyltoxin IV, could thus be determined. The validity of the assignment was high whenever reported ECD data were available for comparison, and the validity was lower when the assignment was based only on the comparison of calculated and reported specific rotations. ECD spectra are intrinsically different for toxins with a biphenyl substructure and for toxins derived from dihydroanthracene.

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, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Calculation of two-photon absorption strengths with the approximate coupled cluster singles and doubles model CC2 using the resolution-of-identity approximation

An implementation of two-photon absorption matrix elements using the approximate second-order coupled-cluster singles and doubles model CC2 is presented. In this implementation we use the resolution-of-the-identity approximation for the two-electron repulsion integrals to reduce the computational cost. To avoid storage of large arrays we introduce in addition a numerical Laplace transformation of orbital energy denominators for the response of the doubles amplitudes. The error due to the numerical Laplace transformation is found to be negligible. Using this new implementation, we performed a series of benchmark calculations on substituted benzene and azobenzene derivatives to get reference values for TD-DFT results. We show that results obtained with the Coulomb-attenuated B3LYP functional are in reasonable agreement with the coupled-cluster results, whereas other density functionals which do not have a long-range correction give considerably less accurate results. Applications to the AF240 dye molecule and a weakly bound molecular tweezer complex demonstrate that this new RI-CC2 implementation allows for the first time to compute two-photon absorption cross sections with a correlated wave function method for molecules with more than 70 atoms and to apply this method for benchmarking TD-DFT calculations on molecules which are of particular relevance for experimental studies of two-photon absorption.