New coupled-cluster methods with singles, doubles, and noniterative triples for high accuracy calculations of excited electronic states

Fragment molecular orbital calculations on red fluorescent proteins (DsRed and mFruits)

We have performed a series of fragment molecular orbital (FMO) calculations for a family of red fluorescent proteins, DsRed and mFruits. The electronic transition energies were evaluated by the method of configuration interaction singles with perturbative doubles [CIS(D)] including higher-order corrections. The calculated values were in good agreement with the corresponding experimental peak values of spectra. Additionally, the chromophore environment was systematically analyzed in terms of the interaction energies between the pigment moiety and neighboring residues. It was theoretically revealed that the electrostatic interactions play a dominant role in the DsRed chromophore, whereas the color tunings in mFruits are controlled in a more delicate fashion.

Extension of the renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed hamiltonian to open-shell systems: a benchmark study

The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; Włoch, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; Włoch, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the F2+ ion and singlet-triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree--Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H --> C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to show that the dependence of the ROHF-based CR-CC(2,3) energies on the method of canonicalization of the ROHF orbitals is, for all practical purposes, negligible.

Structure of Low-Lying Excited States of Guanine in DNA and Solution: Combined Molecular Mechanics and High-Level Coupled Cluster Studies

High-level ab-initio equation-of-motion coupled-cluster methods with singles, doubles, and noniterative triples are used, in conjunction with the combined quantum mechanical molecular mechanics approach, to investigate the structure of low-lying excited states of the guanine base in DNA and solvated environments. Our results indicate that while the excitation energy of the first excited state is barely changed compared to its gas-phase counterpart, the excitation energy of the second excited state is blue-shifted by 0.24 eV.

Reaction pathways and excited states in H2O2+OH→HO2+H2O: A new ab initio investigation

The mechanism of the hydrogen abstraction reaction H(2)O(2)+OH-->HO(2)+H(2)O in gas phase was revisited using density functional theory and other highly correlated wave function theories. We located two pathways for the reaction, both going through the same intermediate complex OH-H(2)O(2), but via two distinct transition state structures that differ by the orientation of the hydroxyl hydrogen relative to the incipient hydroperoxy hydrogen. The first two excited states were calculated for selected points on the pathways. An avoided crossing between the two excited states was found on the product side of the barrier to H transfer on the ground state surface, near the transition states. We report on the calculation of the rate of the reaction in the gas phase for temperatures in the range of 250-500 K. The findings suggest that the strong temperature dependence of the rate at high temperatures is due to reaction on the low-lying excited state surface over a barrier that is much larger than on the ground state surface.

Assessment of time-dependent density functional schemes for computing the oscillator strengths of benzene, phenol, aniline, and fluorobenzene

In present study the relevance of using the time-dependent density functional theory (DFT) within the adiabatic approximation for computing oscillator strengths (f) is assessed using different LDA, GGA, and hybrid exchange-correlation (XC) functionals. In particular, we focus on the lowest-energy valence excitations, dominating the UV/visible absorption spectra and originating from benzenelike HOMO(pi)-->LUMO(pi(*)) transitions, of several aromatic molecules: benzene, phenol, aniline, and fluorobenzene. The TDDFT values are compared to both experimental results obtained from gas phase measurements and to results determined using several ab initio schemes: random phase approximation (RPA), configuration interaction single (CIS), and a series of linear response coupled-cluster calculations, CCS, CC2, and CCSD. In particular, the effect of the amount of Hartree-Fock (HF) exchange in the functional is highlighted, whereas a basis set investigation demonstrates the need of including diffuse functions. So, the hybrid XC functionals--and particularly BHandHLYP--provide f values in good agreement with the highly correlated CCSD scheme while these can be strongly underestimated using pure DFT functionals. These results also display systematic behaviors: (i) larger f and squares of the transition dipole moments (mid R:mumid R:(2)) are associated with larger excitation energies (DeltaE); (ii) these relationships present generally a linear character with R>0.9 in least-squares fit procedures; (iii) larger amounts of HF exchange in the XC functional lead to larger f, R:mumid R:(2), as well as DeltaE values; (iv) these increases in f, mid R:mumid R:(2), and DeltaE are related to increased HOMO-LUMO character; and (v) these relationships are, however, not universal since the linear regression parameters (the slopes and the intercepts at the origin) depend on the system under investigation as well as on the nature of the excited state.

On the origin of the ultrafast internal conversion of electronically excited pyrimidine bases

The ultrafast radiationless decay of photoexcited uracil and cytosine has been investigated by ab initio quantum chemical methods based on CIS and CR-EOM-CCSD(T) electronic energy calculations at optimized CIS geometries. The calculated potential energy profiles indicate that the S(1) --> S(0) internal conversion of the pyrimidine bases occurs through a barrierless state switch from the initially excited (1)pipi state to the out-of-plane deformed excited state of biradical character, which intersects the ground state at a lower energy. This three-state nonradiative decay mechanism predicts that replacement of the C5 hydrogen by fluorine introduces an energy barrier for the initial state switch, whereas replacement of the C6 hydrogen by fluorine does not. These predictions are borne out by the very different fluorescence yields of 5-fluorinated bases relative to the corresponding 6-fluorinated bases. It is concluded from these results that the origin of the ultrafast radiationless decay is the same for the two pyrimidine bases.

Ab initio study of a biradical radiationless decay channel of the lowest excited electronic state of cytosine and its derivatives

A theoretical model for the ultrafast S1-->S0 internal conversion of cytosine is presented, in which a state switch from the initially prepared 1pipi* state to the out-of-plane deformed excited state of biradical character controls the rate of the S1(1pipi*) decay. This mechanism successfully accounts for the dramatically longer S1 lifetimes of 5-fluorocytosine and N-acetylcytosine relative to cytosine. The replacement of the C5 hydrogen atom by a methyl group is predicted to lead to a substantial, but not dramatic, increase in the S1 lifetime, also consistent with experiment. It is this ability to correctly predict the substituent effects that distinguishes the present model from the previously proposed mechanisms.

Second- and third-order triples and quadruples corrections to coupled-cluster singles and doubles in the ground and excited states

Second- and third-order perturbation corrections to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) incorporating excited configurations in the space of triples [EOM-CCSD(2)T and (3)T] or in the space of triples and quadruples [EOM-CCSD(2)TQ] have been implemented. Their ground-state counterparts--third-order corrections to coupled-cluster singles and doubles (CCSD) in the space of triples [CCSD(3)T] or in the space of triples and quadruples [CCSD(3)TQ]--have also been implemented and assessed. It has been shown that a straightforward application of the Rayleigh-Schrodinger perturbation theory leads to perturbation corrections to total energies of excited states that lack the correct size dependence. Approximations have been introduced to the perturbation corrections to arrive at EOM-CCSD(2)T, (3)T, and (2)TQ that provide size-intensive excitation energies at a noniterative O(n(7)), O(n(8)), and O(n(9)) cost (n is the number of orbitals) and CCSD(3)T and (3)TQ size-extensive total energies at a noniterative O(n(8)) and O(n(10)) cost. All the implementations are parallel executable, applicable to open and closed shells, and take into account spin and real Abelian point-group symmetries. For excited states, they form a systematically more accurate series, CCSD<CCSD(2)T<CCSD(2)TQ<CCSD(3)T<CCSDT, with the second- and third-order corrections capturing typically approximately 80% and 100% of such effects, when those effects are large (>1 eV) and the ground-state wave function has single-determinant character. In other cases, however, the corrections tend to overestimate the triples and quadruples effects, the origin of which is discussed. For ground states, the third-order corrections lead to a rather small improvement over the highly effective second-order corrections [CCSD(2)T and (2)TQ], which is a manifestation of the staircase convergence of perturbation series.

Calculations of Molecular Properties in Hybrid Coupled-Cluster and Molecular Mechanics Approach

We report benchmark calculations obtained with our new coupled-cluster singles and doubles (CCSD) code for calculating the first- and second-order molecular properties. This code can be easily incorporated into combined [Valiev, M.; Kowalski, K. J. Chem. Phys. 2006, 125, 211101] classical molecular mechanics (MM) and ab initio coupled-cluster (CC) calculations using NWChem, enabling us to study molecular properties in a realistic environment. To test this methodology, we discuss the results of calculations of dipole moments and static polarizabilities for the Cl2O system in the CCl4 solution using the CCSD (CC with singles and doubles) linear response approach. We also discuss the application of the asymptotic extrapolation scheme (AES) [Kowalski, K.; Valiev, M. J. Phys. Chem. A 2006, 110, 13106] in reducing the numerical cost of CCSD calculations.