Search

We present a hierarchy of local coupled cluster (CC) linear response (LR) methods to calculate ionization potentials (IPs), i.e., excited states with one electron annihilated relative to a ground state reference. The time-dependent perturbation operator V(t), as well as the operators related to the first-order (with respect to V(t)) amplitudes and multipliers, thus are not number conserving and have half-integer particle rank m. Apart from calculating IPs of neutral molecules, the method offers also the possibility to study ground and excited states of neutral radicals as ionized states of closed-shell anions. It turns out that for comparable accuracy IPs require a higher-order treatment than excitation energies; an IP-CC LR method corresponding to CC2 LR or the algebraic diagrammatic construction scheme through second order performs rather poorly. We therefore systematically extended the order with respect to the fluctuation potential of the IP-CC2 LR Jacobian up to IP-CCSD LR, keeping the excitation space of the first-order (with respect to V(t)) cluster operator restricted to the m = 1/2 ⊕ 3/2 subspace and the accuracy of the zero-order (ground-state) amplitudes at the level of CC2 or MP2. For the more expensive diagrams beyond the IP-CC2 LR Jacobian, we employ local approximations. The implemented methods are capable of treating large molecular systems with hundred atoms or more. [ABSTRACT FROM AUTHOR]

We introduce orbital specific virtuals (OSVs) to represent the truncated pair-specific virtual space in periodic local Møller-Plesset perturbation theory of second order (LMP2). The OSVs are constructed by diagonalization of the LMP2 amplitude matrices which correspond to diagonal Wannier-function (WF) pairs. Only a subset of these OSVs is adopted for the subsequent OSV-LMP2 calculation, namely, those with largest contribution to the diagonal pair correlation energy and with the accumulated value of these contributions reaching a certain accuracy. The virtual space for a general (non diagonal) pair is spanned by the union of the two OSV sets related to the individual WFs of the pair. In the periodic LMP2 method, the diagonal LMP2 amplitude matrices needed for the construction of the OSVs are calculated in the basis of projected atomic orbitals (PAOs), employing very large PAO domains. It turns out that the OSVs are excellent to describe short range correlation, yet less appropriate for long range van der Waals correlation. In order to compensate for this bias towards short range correlation, we augment the virtual space spanned by the OSVs by the most diffuse PAOs of the corresponding minimal PAO domain. The Fock and overlap matrices in OSV basis are constructed in the reciprocal space. The 4-index electron repulsion integrals are calculated by local density fitting and, for distant pairs, via multipole approximation. New procedures for determining the fit-domains and the distant-pair lists, leading to higher efficiency in the 4-index integral evaluation, have been implemented. Generally, and in contrast to our previous PAO based periodic LMP2 method, the OSV-LMP2 method does not require anymore great care in the specification of the individual domains (to get a balanced description when calculating energy differences) and is in that sense a black box procedure. Discontinuities in potential energy surfaces, which may occur for PAO-based calculations if one is not careful, virtually disappear for OSV-LMP2. Moreover, due to much increased compactness of the pair-specific virtual spaces, the OSV-LMP2 calculations are faster and require much less memory than PAO-LMP2 calculations, despite the noticeable overhead of the initial OSV construction procedure. [ABSTRACT FROM AUTHOR]

We describe theory and implementation of oscillator strengths, orbital-relaxed first-order properties, and nuclear gradients for the local algebraic diagrammatic construction scheme through second order. The formalism is derived via time-dependent linear response theory based on a second-order unitary coupled cluster model. The implementation presented here is a modification of our previously developed algorithms for Laplace transform based local time-dependent coupled cluster linear response (CC2LR); the local approximations thus are state specific and adaptive. The symmetry of the Jacobian leads to considerable simplifications relative to the local CC2LR method; as a result, a gradient evaluation is about four times less expensive. Test calculations show that in geometry optimizations, usually very similar geometries are obtained as with the local CC2LR method (provided that a second-order method is applicable). As an exemplary application, we performed geometry optimizations on the low-lying singlet states of chlorophyllide a. [ABSTRACT FROM AUTHOR]

In order to arrive at linear scaling of the computational cost with molecular size, local coupled cluster methods discriminate pairs of local molecular orbitals according to the spatial separation R of the latter. Only strong pairs are treated at the full coupled cluster level, whereas for weak pairs a lower level of theory (usually Møller-Plesset perturbation theory of second order,MP2) is used. Yet anMP2 treatment of weak pairs is inadequate in certain situations (for example, for describing p-stacking), which calls for an improved but still inexpensive method for dealing with the weak pairs. In a previous contribution, we proposed as a substituent for MP2 the LrCCD3 method, which is based on ring coupled cluster doubles (ring-CCD) and includes all third-order diagrams with energy contributions decaying not quicker than R-6. In the present work, we explore a still more accurate method, which is based on the same principles. It turned out to be essential to abandon the restriction to ring-CCD, i.e., to include further CCD diagrams beyond the ring approximation. The occurring intermediates turn out to be formally very similar to LMP2 density matrices, such that an efficient evaluation of these non-ring CCD diagrams is possible. Furthermore, a computationally cheap a posteriori estimate for the fourth-order singles contribution to the weak pair energy, which also exhibits a decay behavior of R-6, is introduced. The resulting method, denoted as LCCD[S]-R-6, indeed provides a substantial improvement in accuracy over the previous LrCCD3 method at a relatively modest additional computational cost. [ABSTRACT FROM AUTHOR]

A multistate local CC2 response method for the calculation of analytic energy gradients with respect to nuclear displacements is presented for ground and electronically excited states. The gradient enables the search for equilibrium geometries of extended molecular systems. Laplace transform is used to partition the eigenvalue problem in order to obtain an effective singles eigenvalue problem and adaptive, state-specific local approximations. This leads to an approximation in the energy Lagrangian, which however is shown (by comparison with the corresponding gradient method without Laplace transform) to be of no concern for geometry optimizations. The accuracy of the local approximation is tested and the efficiency of the new code is demonstrated by application calculations devoted to a photocatalytic decarboxylation process of present interest. [ABSTRACT FROM AUTHOR]

A multistate local CC2 response method for the calculation of orbital-relaxed first order properties is presented for ground and electronically excited states. It enables the treatment of excited state properties including orbital relaxation for extended molecular systems and is a major step on the way towards analytic gradients with respect to nuclear displacements. The Laplace transform method is employed to partition the eigenvalue problem and the lambda equations, i.e., the doubles parts of these equations are inverted on-the-fly, leaving only the corresponding effective singles equations to be solved iteratively. Furthermore, the state specific local approximations are adaptive. Density-fitting is utilized to decompose the electron-repulsion integrals. The accuracy of the local approximation is tested and the efficiency of the new code is demonstrated on the example of an organic sensitizer for solar-cell applications, which consists of about 100 atoms. [ABSTRACT FROM AUTHOR]

A local method based on orbital specific virtuals (OSVs) for calculating the perturbative triples correction in local coupled cluster calculations is presented. In contrast to the previous approach based on projected atomic orbitals (PAOs), described by Schütz [J. Chem. Phys. 113, 9986 (2000)], the new scheme works without any ad hoc truncations of the virtual space to domains. A single threshold defines the pair and triple specific virtual spaces completely and automatically. It is demonstrated that the computational cost of the method scales linearly with molecular size. Employing the recommended threshold a similar fraction of the correlation energy is recovered as with the original PAO method at a somewhat lower cost. A benchmark for 52 reactions demonstrates that for reaction energies the intrinsic accuracy of the coupled cluster with singles and doubles excitations and a perturbative treatment of triples excitations method can be reached by OSV-local coupled cluster theory with singles and doubles and perturbative triples, provided a MP2 correction is applied that accounts for basis set incompleteness errors as well as for remaining domain errors. As an application example the interaction energies of the guanine-cytosine dimers in the Watson-Crick and stacked arrangements are investigated at the level of local coupled cluster theory with singles and doubles and perturbative triples. Based on these calculations we propose new complete-basis-set-limit estimates for these interaction energies at this level of theory. [ABSTRACT FROM AUTHOR]

We present a density fitted local configuration interaction singles (CIS) method for calculating optical bandgaps in 3D-periodic systems. We employ an Ewald technique to carry out infinite lattice summations for the exciton-exciton interaction, and robust product-density specific local density fitting in direct space for the electron-hole interaction. Moreover, we propose an alternative to the usual cyclic model with Born-von Karman periodic boundary conditions, the so called Wigner-Seitz supercell truncated infinite model, which exhibits much improved convergence of the CIS excitation energy with respect to the size of the supercell. Test calculations on a series of prototypical systems demonstrate that the method at the present stage can be used to calculate the excitonic bandgaps of 3D periodic systems with up to a dozen atoms in the unit cell, ranging from wide-gap insulators to semiconductors. [ABSTRACT FROM AUTHOR]

An efficient method for the calculation of nuclear magnetic resonance (NMR) shielding tensors is presented, which treats electron correlation at the level of second-order Mo\ller-Plesset perturbation theory. It uses spatially localized functions to span occupied and virtual molecular orbital spaces, respectively, which are expanded in a basis of gauge including atomic orbitals (GIAOs or London atomic orbitals). Doubly excited determinants are restricted to local subsets of the virtual space and pair energies with an interorbital distance beyond a certain threshold are omitted. Furthermore, density fitting is employed to factorize the electron repulsion integrals. Ordinary Gaussians are employed as fitting functions. It is shown that the errors in the resulting NMR shielding constant, introduced (i) by the local approximation and (ii) by density fitting, are very small or even negligible. The capabilities of the new program are demonstrated by calculations on some extended molecular systems, such as the cyclobutane pyrimidine dimer photolesion with adjacent nucleobases in the native intrahelical DNA double strand (ATTA sequence). Systems of that size were not accessible to correlated ab initio calculations of NMR spectra before. The presented method thus opens the door to new and interesting applications in this area. [ABSTRACT FROM AUTHOR]

We extend the orbital-specific-virtual tensor factorization, introduced for local Mo\ller-Plesset perturbation theory in Ref. [J. Yang, Y. Kurashige, F. R. Manby and G. K. L. Chan, J. Chem. Phys. 134, 044123 (2011)], to local coupled cluster singles and doubles theory (OSV-LCCSD). The method is implemented by modifying an efficient projected-atomic-orbital local coupled cluster program (PAO-LCCSD) described recently, [H.-J. Werner and M. Schütz, J. Chem. Phys. 135, 144116 (2011)]. By comparison of both methods we find that the compact representation of the amplitudes in the OSV approach affords various advantages, including smaller computational time requirements (for comparable accuracy), as well as a more systematic control of the error through a single energy threshold. Overall, the OSV-LCCSD approach together with an MP2 correction yields small domain errors in practical calculations. The applicability of the OSV-LCCSD is demonstrated for molecules with up to 73 atoms and realistic basis sets (up to 2334 basis functions). [ABSTRACT FROM AUTHOR]

The atomic orbital basis set limit is approached in periodic correlated calculations for solid LiH. The valence correlation energy is evaluated at the level of the local periodic second order Mo\ller-Plesset perturbation theory (MP2), using basis sets of progressively increasing size, and also employing 'bond'-centered basis functions in addition to the standard atom-centered ones. Extended basis sets, which contain linear dependencies, are processed only at the MP2 stage via a dual basis set scheme. The local approximation (domain) error has been consistently eliminated by expanding the orbital excitation domains. As a final result, it is demonstrated that the complete basis set limit can be reached for both HF and local MP2 periodic calculations, and a general scheme is outlined for the definition of high-quality atomic-orbital basis sets for solids. [ABSTRACT FROM AUTHOR]

We present a density fitted local configuration interaction singles (CIS) method for calculating optical band gaps in 1D-periodic systems. The method is based on the Davidson diagonalization procedure, carried out in the reciprocal space. The one-electron part of the matrix-vector products is also evaluated in the reciprocal space, where the diagonality of the Fock matrix can be exploited. The contraction of the CIS vectors with the two electron integrals is performed in the direct space in the basis of localized occupied (Wannier) and virtual (projected atomic) orbitals. The direct space approach allows to utilize the sparsity of the integrals due to the local representation and locality of the exciton. The density fitting approximation employed for the two electron integrals reduces the nominal scaling with unit cell size to O(N4). Test calculations on a series of prototypical systems demonstrate that the method in its present stage can be used to calculate the excitonic band gaps of polymers with up to a few dozens of atoms in the cell. The computational cost depends on the locality of the exciton, but even relatively delocalized excitons occurring in the polybiphenyl in the parallel orientation, can be routinely treated with this method. [ABSTRACT FROM AUTHOR]

A new multistate local CC2 response method for calculating excitation energies and first-order properties of excited triplet states in extended molecular systems is presented. The Laplace transform technique is employed to partition the left/right local CC2 eigenvalue problems as well as the linear equations determining the Lagrange multipliers needed for the properties. The doubles part in the equations can then be inverted on-the-fly and only effective equations for the singles part must be solved iteratively. The local approximation presented here is adaptive and state-specific. The density-fitting method is utilized to approximate the electron-repulsion integrals. The accuracy of the new method is tested by comparison to canonical reference values for a set of 12 test molecules and 62 excited triplet states. As an illustrative application example, the lowest four triplet states of 3-(5-(5-(4-(bis(4-(hexyloxy)phenyl)amino)phenyl)thiophene-2-yl)thiophene-2-yl)-2-cyanoacrylic acid, an organic sensitizer for solar-cell applications, are computed in the present work. No triplet charge-transfer states are detected among these states. This situation contrasts with the singlet states of this molecule, where the lowest singlet state has been recently found to correspond to an excited state with a pronounced charge-transfer character having a large transition strength. [ABSTRACT FROM AUTHOR]

The lowest electronically excited states of the aniline dimer and trimer related to the lowest π*←π transition of the monomer are investigated by applying time-dependent coupled cluster theory, primarily at the level of the (spin-component-scaled) CC2 model. Minimum energy structures in the vicinity of the Franck-Condon points were determined on the individual potential energy surfaces. For the dimer we find an excimer and a head-to-tail configuration (with the monomers substantially displaced relative to the ground state minimum) for the lowest (dark) and second lowest (bright) states, respectively. The excitation is delocalized on both chromophores for both of these states. For the trimer three distinct minima with quite different hydrogen-bonding arrangements are found for the three lowest states. In strong contrast to the dimer the excitation here is clearly localized on the individual aniline chromophores for each of these three states. One of the three geometries is rather similar to the ground state minimum, while the two others are rather different and thus have presumably quite small Franck-Condon factors. It can be expected that only the electronic origin of the first conformer can eventually be detected in the absorption spectrum of the trimer, provided that it is separated by high-enough barriers from other, energetically lower configurations. [ABSTRACT FROM AUTHOR]

The aniline dimer and trimer are investigated in their electronic ground state. The potential energy surface was thoroughly searched for low lying minima by applying global optimizations on a model potential, which is recalibrated on-the-fly by ab initio calculations (spin-component-scaled LMP2) at relevant configurations. The most stable structure of the dimer corresponds to a head-to-tail arrangement with both aniline monomers being nearly equivalent. DFT-SAPT calculations reveal that the interaction energy is dominated by van der Waals dispersion, which is of comparable size as for the benzene dimer, but with a much larger total interaction energy than for the latter. The global minimum of the aniline trimer corresponds to a hydrogen bonding arrangement involving three directional NH–N hydrogen bonds, with the individual monomers being clearly distinguishable. Nonadditive three-body dispersion contributions appear to play a minor role for the trimer. [ABSTRACT FROM AUTHOR]

We have calculated the equilibrium geometry, formation energy, and bulk modulus of two molecular bulk crystals, NH3 and CO2, at the periodic post-Hartree–Fock correlated level. The dependence of the results on the basis set has been analyzed, by employing basis sets up to aug-cc-pVQZ quality. In the calculations, we used the periodic local Mo\ller–Plesset second order perturbation theory (LMP2), implemented in the CRYSCOR program. Multipolar expansion techniques, as well as density fitting, are employed in this code to reduce the number of and to factorize the required electron repulsion integrals; as a consequence of that, the computational cost for the correlation part of the calculations is comparable to that of the Hartree–Fock. Auxiliary calculations performed on molecular dimers are also reported to verify the accuracy of the LMP2 approach and of the basis sets used. Furthermore, the effect of spin-component scaling has been investigated for the two crystals. One intention of the present paper is also to lay out and specify the computational setup, which is generally applicable for accurate CRYSCOR calculations on molecular crystals. [ABSTRACT FROM AUTHOR]

A new Laplace transform based multistate local CC2 response method for calculating excitation energies of extended molecular systems is presented. By virtue of the Laplace transform trick, the eigenvalue problem involving the local CC2 Jacobian is partitioned along the doubles-doubles block (which is diagonal in the parental canonical method) without losing the sparsity in the integral, amplitude, and amplitude response supermatrices. Hence, only an effective eigenvalue problem involving singles vectors has to be solved, while the doubles part can be computed on-the-fly. Within this framework, a multistate treatment of excited states with state specific and adaptive local approximations imposed on the doubles part is straightforwardly possible. Furthermore, in the context of the density fitting approximation of the two-electron integrals, a procedure to specify the local approximation, i.e., the restricted pair lists and domains, on the basis of an analysis of the object to be approximated itself is proposed. Performance and accuracy of the new Laplace transformed density fitted local CC2 (LT-DF-LCC2) response method are tested for set of different test molecules and states. It turns out that LT-DF-LCC2 response is much more robust than the earlier local CC2 response method proposed before, which failed to find some excited states in difficult cases. [ABSTRACT FROM AUTHOR]

The potential energy surfaces of the S0 and S1(π*←π) states of the 2-naphthol(H2O)n, n∈{1,2} clusters were explored at the level of coupled cluster (CC2) response theory. In the electronic ground state two different types of hydrogen-bonding networks coexist for n=2, (i) a cyclic one [similar to those of the water trimer and phenol(H2O)2] where the hydroxy group of the aryl alcohol acts simultaneously as H donor for the first, and as H acceptor for the second water molecule, and (ii) a hydrogen-bonding arrangement where the aromatic π system is taking over the role as H acceptor. In the S1 state, on the other hand, the cyclic conformers are unstable. Consequently, the first group of cyclic ground state conformers gives rise to broad unstructured band shapes in the absorption spectrum, whereas the second group of conformers involving the aromatic π system gives rise to nicely structured band shapes. Based on these results the puzzling absorption spectrum of the n=2 cluster can properly be interpreted. [ABSTRACT FROM AUTHOR]