An extended multireference study of the electronic states of para-benzyne

A state-averaged, multireference complete active space (CAS) approach was used for the determination of the vertical excitation energies of valence and Rydberg states of para-benzyne. Orbitals were generated with a 10- and 32-state averaged multiconfigurational self-consistent field approach. Electron correlation was included using multireference configuration interaction with singles and doubles, including the Pople correction for size extensivity, multireference averaged quadratic coupled cluster (MR-AQCC), and MR-AQCC based on linear response theory. There is a very high density of electronic states in this diradical system-there are more than 17 states within 7 eV of the ground state including two 3s Rydberg states. All excitations, except 2 (1)A(g), are from the pi system to the sigmasigma(*) system. Of the 32 states characterized, 15 were multiconfigurational, including the ground (1)A(g) state, providing further evidence for the necessity of a multireference approach for p-benzyne. The vertical singlet-triplet splitting was also characterized using a two-state averaged approach. A CAS(2,2) calculation was shown to be inadequate due to interaction with the pi orbitals.     

Linear-response theory for Mukherjee's multireference coupled-cluster method: Excitation energies

The recently presented linear-response function for Mukherjee's multireference coupled-cluster method (Mk-MRCC) [T.-C. Jagau and J. Gauss, J. Chem. Phys. 137, 044115 (2012)] is employed to determine vertical excitation energies within the singles and doubles approximation (Mk-MRCCSD-LR) for ozone as well as for o-benzyne, m-benzyne, and p-benzyne, which display increasing multireference character in their ground states. In order to assess the impact of a multireference ground-state wavefunction on excitation energies, we compare all our results to those obtained at the single-reference coupled-cluster level of theory within the singles and doubles as well as within the singles, doubles, and triples approximation. Special attention is paid to the artificial splitting of certain excited states which arises from the redundancy intrinsic to Mk-MRCC theory and hinders the straightforward application of the Mk-MRCC-LR method.     

Exploiting regularity in systematic sequences of wavefunctions which approach the full CI limit

Summary The ground state total energy and related 1-electron properties are computed for three small molecules (N₂, H₂O, and H₂CN) using several systematic sequences of wavefunctions which approach the full CI. These sequences include multireference CI, averaged coupled pair functional and quasidegenerate variational perturbation theory wavefunctions. It is demonstrated that sufficient regularity exists in the sequence of variationally computed energies to permit extrapolation to the full CI limit using simple analytic expressions. It is furthermore demonstrated that a subset of the original list of configurations employed in the normal singles and doubles CI procedure can be selected using second order perturbation theory without adversely affecting the extrapolation to the full CI limit. This significantly broadens the range of applicability of the method. Along these lines, a scheme is proposed for the extrapolation of the selected CI results to the zero threshold (i.e. unselected) values in cases where the numbers of configurations associated with the latter would render the calculations intractable. Due to the vast reduction in the number of configurations which are handled variationally, the proposed scheme makes it possible to derive estimates of the full CI limit in cases where explicit full CI is either very difficult or currently impossible.Dedicated to Prof. Klaus RuedenbergThe Pacific Northwest Laboratory is operated for the U.S. Department of Energy by Battelle Memorial Institute under contract DE-AC06-76RLO 1830     

Basis set limit electronic excitation energies, ionization potentials, and electron affinities for the 3d transition metal atoms: Coupled cluster and multireference methods

Recently developed correlation consistent basis sets for the first row transition metal elements Sc-Zn have been utilized to determine complete basis set (CBS) scalar relativistic electron affinities, ionization potentials, and 4s(2)3d(n-2)-4s(1)d(n-1) electronic excitation energies with single reference coupled cluster methods [CCSD(T), CCSDT, and CCSDTQ] and multireference configuration interaction with three reference spaces: 3d4s, 3d4s4p, and 3d4s4p3d'. The theoretical values calculated with the highest order coupled cluster techniques at the CBS limit, including extrapolations to full configuration interaction, are well within 1 kcal/mol of the corresponding experimental data. For the early transition metal elements (Sc-Mn) the internally contracted multireference averaged coupled pair functional method yielded excellent agreement with experiment; however, the atomic properties for the late transition metals (Mn-Zn) proved to be much more difficult to describe with this level of theory, even with the largest reference function of the present work.     

Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules

We report an extension of the local correlation concept to electronically excited states via the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) method. We apply the same orbital domain structure used successfully for ground-state CCSD by Werner and co-workers and find that the resulting localized excitation energies are in error generally by less than 0.2 eV relative to their canonical EOM-CCSD counterparts, provided the basis set is flexible and includes Rydberg-like functions. In addition, we account for weak-pair contributions efficiently using a correction to local-EOM-CCSD transition energies based on the perturbative (D) correction used with configuration interaction singles (CIS).     

A simple correction to final state energies of doublet radicals described by equation-of-motion coupled cluster theory in the singles and doubles approximation

A computationally inexpensive energy correction is suggested for radicals described by the equation-of-motion coupled cluster method for ionized states in the singles and doubles approximation (EOMIP-CCSD). The approach is primarily intended for doublet states that are qualitatively described by Koopmans' approximation. Following a strategy similar to those used in multireference coupled cluster theory, the proposed correction accounts for all correlation effects through third order in perturbation theory and also includes selected contributions to higher-order energies. As an initial test of the numerical performance of the method, total energies and energy splittings are calculated for some small prototype radicals.     

Generalization of coupled-cluster response theory to multireference expansion spaces: application of the coupled-cluster singles and doubles effective Hamiltonian

Employing separate cluster ansatz in time-independent and time-dependent wave-operators, coupled-cluster (CC) response theory is generalized to multireference (MR) expansion spaces. For state energies, this corresponds to the MR secular problem with an arbitrary similarity-transformed effective Hamiltonian, H˜=Ω⁻¹ HΩ. The effective Hamiltonian can be generated via size-extensive CC methods. Thus the states in MR linear response theory (MRLRT) maintain the usual CC core-extensive properties. We have used the Gelfand unitary group basis of the spin-adapted configurations to construct the matrix of H˜ in the MR excitation space. As a preliminary application, the CC singles and doubles effective Hamiltonian is applied to excitation and photoionization energies of the CH⁺ and N₂ molecules, and is compared with experimental results and results from other numerical procedures including conventional CC linear response theory (CC-LRT), MR and full configuration interaction (MRCI and FCI) methods. The numerical results indicate that MRLRT reproduces valence and external excited states quantitatively, combining the best features of CC-LRT and MRCI. Received: 2 July 1998 / Accepted: 28 August 1998 / Published online: 11 November 1998     

Multireference perturbation theory with optimized partitioning. II. Applications to molecular systems

The second‐order multireference perturbation theory using an optimized partitioning, denoted as MROPT(2), is applied to calculations of various molecular properties—excitation energies, spectroscopic parameters, and potential energy curves—for five molecules: ethylene, butadiene, benzene, N₂, and O₂. The calculated results are compared with those obtained with second‐ and third‐order multireference perturbation theory using the traditional partitioning techniques. We also give results from computations using the multireference configuration interaction (MRCI) method. The presented results show very close resemblance between the new method and MRCI with renormalized Davidson correction. The accuracy of the new method is good and is comparable to that of second‐order multireference perturbation theory using Møller‐Plesset partitioning. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1390–1400, 2003     

Quantum chemical comparison of vertical,adiabatic, and 0‐0 excitation energies: The PYP and GFP chromophores

A number of benchmark studies investigating the performance of quantum chemical methods for calculating vertical excitation energies are today available in the literature. However, less established is the variation between methods in their estimates of the differences between vertical, adiabatic, and 0‐0 excitation energies. To this end, such excitation energies are here calculated for the bright S₁ states of the anionic chromophores of the photoactive yellow protein (PYP) and the green fluorescent protein (GFP) in the gas phase using configuration interaction singles, complete active space self‐consistent field, coupled‐cluster singles and doubles, and time‐dependent density functional theory methods. Although the estimates of the excitation energies vary by more than 1 eV between the methods, the differences between the different types of excitation energies are found to be relatively method‐insensitive, varying by ～0.1 eV only for these particular chromophores. Specifically, the adiabatic energies are uniformly 0.10–0.17 (PYP) and 0.06–0.17 eV (GFP) lower than the vertical energies, and the 0‐0 energies are similarly 0.09–0.14 (PYP) and 0.07–0.17 eV (GFP) lower than the adiabatic energies. © 2012 Wiley Periodicals, Inc.     

The lowest 2A' excited state of the water-hydroxyl complex

Vertical and adiabatic excitation energies of the lowest (2)A(') excited state in the water-hydroxyl complex have been determined using coupled cluster, multireference configuration interaction, multireference perturbation theory, and density-functional methods. A significant redshift of about 0.4 eV in the vertical excitation energy of the complex compared to that of the hydroxyl radical monomer is found with the coupled cluster calculations validating previous results. Electronic excitation leads to a structure with near-equal sharing of the hydroxyl hydrogen by both oxygen atoms and a concomitantly large redshift of the adiabatic excitation energy of approximately 1 eV relative to the vertical excitation energy. The combination of redshifts ensures that the electronic transition in the complex lies well outside the equivalent excitation in the hydroxyl radical monomer. The complex is approximately five times more strongly bound in the excited state than in the ground state.     

Torsional effects on excitation energies of thiophene derivatives induced by beta-substituents: comparison between time-dependent density functional theory and approximated coupled cluster approaches

The influence of methyl or phenyl substitution in beta-position of dioxygenated terthiophene and diphenylthiophene on the optical properties is investigated by first-principles calculations. We compare the approximated singles and doubles coupled cluster (CC2) approach with time-dependent density functional theory methods. CC2 reproduces experimental excitation energies with an accuracy of 0.1 eV. We find that the different substituents modify the inter-ring torsional angle which in turn strongly influences the excitation energies. The steric contribution to the excitation energies have been separated from the total substituent effects.     

A high‐accuracy theoretical study of the CHnP Systems n = 1–3

We have performed high‐level electronic structure computations on the most important species of the CH_nP systems n = 1–3 to characterize them and provide reliable information about the equilibrium and vibrationally averaged molecular structures, rotational constants, vibrational frequencies (harmonic and anharmonic), formation enthalpies, and vertical excitation energies. Those chemical systems are intermediates for several important reactions and also prototypical phosphorus‐carbon compounds; however, they are often elusive to experimental detection. The present results significantly complement their knowledge and can be used as an assessment of the experimental information when available. The explicitly correlated coupled‐cluster RCCSD(T)‐F12 method has been used for geometry optimizations and vibrational frequency calculations. Vibrational configuration interaction theory has been used to account for anharmonicity effects. Basis‐set limit extrapolations have been carried out to determine accurate thermochemical quantities. Electronic excited states have been calculated with coupled‐cluster approaches and also by means of the multireference configuration interaction method. © 2013 Wiley Periodicals, Inc.     

Electronic spectra of nitroethylene

A systematic study of the electronic excited states of nitroethylene (C₂H₃NO₂) was carried out using the approximate coupled‐cluster singles‐and‐doubles approach with the resolution of the identity (RI‐CC2), the time dependent density functional theory with the CAMB3LYP functional (TDDFT/CAMB3LYP) and the DFT multireference configuration interaction (DFT/MRCI) method. Vertical transition energies and optical oscillator strengths were computed for a maximum of 20 singlet transitions. Semiclassical simulations of the ultraviolet (UV) spectra were performed at the RI‐CC2 and DFT/MRCI levels. The main features in the UV spectrum were assigned to a weak n‐π* transition, and two higher energy π_CC+O‐π* bands. These characteristics are common to molecules containing NO₂ groups. Simulated spectra are in good agreement with the experimental spectrum. The energy of the bands in the DFT/MRCI simulation agrees quite well with the experiment, although it overestimates the band intensities. RI‐CC2 produced intensities comparable to the experiment, but the bands were blue shifted. A strong π_CC+O‐π* band, not previously measured, was found in the 8–9 eV range. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Coupled-cluster methods including noniterative corrections for quadruple excitations

A new method is presented for treating the effects of quadruple excitations in coupled-cluster theory. In the approach, quadruple excitation contributions are computed from a formula based on a non-Hermitian perturbation theory analogous to that used previously to justify the usual noniterative triples correction used in the coupled cluster singles and doubles method with a perturbative treatment of the triple excitations (CCSD(T)). The method discussed in this paper plays a parallel role in improving energies obtained with the full coupled-cluster singles, doubles, and triples method (CCSDT) by adding a perturbative treatment of the quadruple excitations (CCSDT(Q)). The method is tested for an extensive set of examples, and is shown to provide total energies that compare favorably with those obtained with the full singles, doubles, triples, and quadruples (CCSDTQ) method.     

Comprehensive relativistic ab initio and density functional theory studies on PtH, PtF, PtCl, and Pt(NH(3))(2)Cl(2)

Platinum monohydride is taken as an example to compare the performance of various relativistic and correlation approaches, such as all-electron DPT (direct perturbation theory), ECP (effective core potential); RSPT2, RSPT3 (second- and third-order multireference Rayleigh-Schr?dinger perturbation theory), CCSD(T) (coupled-cluster with singles, doubles, and perturbative triples), as well as the four-component relativistic density functional theory. It is shown that first-order DPT performs significantly better than the (first-order) Breit-Pauli Hamiltonian. The performance of different approaches for the excitation energies of the platinum diatomics is discussed critically. The molecular spectroscopic constants for PtF and PtCl are predicted for the first time. The geometric data for several isomers of cis- and trans-Pt(NH(3))(2)Cl(2) are reported. The corresponding energetic data are calculated at relativistic all-electron and ECP-CCSD(T) as well as four-component relativistic density functional levels of theory. Contrary to previous results, it is found that the two C(2v) isomers of cis-Pt(NH(3))(2)Cl(2) are marginally separated in energy, which could be ascribed to Cl-H interactions.     

Hybrid coupled cluster and molecular dynamics approach: application to the excitation spectrum of cytosine in the native DNA environment

Evolution of the excited state energies of cytosine base in the native DNA environment was investigated using a hybrid coupled cluster and classical molecular dynamics approach. The time averaged excitation energies obtained with the variant of the completely renormalized equation-of-motion with singles, doubles, and non-iterative triples approach that includes a bulk of the correlation effects for excited states, are compared with the analogous calculations in the gas phase. Significant blue shifts for the two lowest singlet excitation energies can be observed as a result of the interaction of the quantum system with the surrounding environment.     

Theoretical investigation of the excited states of coumarin dyes for dye-sensitized solar cells

Using time-dependent density functional theory (TD-DFT), configuration interaction single (CIS) method, and approximate coupled cluster singles and doubles (CC2) method, we investigated the absorption spectra of coumarin derivative dyes (C343, NKX-2388, NKX-2311, NKX-2586, and NKX-2677), which have been synthesized for efficient dye-sensitized solar cells. The CC2 calculations are found in good agreement with the experimental results except for the smallest coumarin dye (C343). TD-DFT underestimates the vertical excitation energy of the larger coumarin dyes (NKX-2586 and -2677). Solvents (methanol) are found to induce a red shift of the vertical excitation energies, and their effects on the molecular geometry and the electronic structure are examined in detail. The deprotonated form of coumarin is also investigated, where a blue shift of the vertical excitation energies is observed.     

High-order excitations in state-universal and state-specific multireference coupled cluster theories: model systems

For the first time high-order excitations (n>2) have been studied in three multireference couple cluster (MRCC) theories built on the wave operator formalism: (1) the state-universal (SU) method of Jeziorski and Monkhorst (JM) (2) the state-specific Brillouin-Wigner (BW) coupled cluster method, and (3) the state-specific MRCC approach of Mukherjee (Mk). For the H4, P4, BeH(2), and H8 models, multireference coupled cluster wave functions, with complete excitations ranging from doubles to hextuples, have been computed with a new arbitrary-order string-based code. Comparison is then made to corresponding single-reference coupled cluster and full configuration interaction (FCI) results. For the ground states the BW and Mk methods are found, in general, to provide more accurate results than the SU approach at all levels of truncation of the cluster operator. The inclusion of connected triple excitations reduces the nonparallelism error in singles and doubles MRCC energies by a factor of 2-10. In the BeH(2) and H8 models, the inclusion of all quadruple excitations yields absolute energies within 1 kcal mol(-1) of the FCI limit. While the MRCC methods are very effective in multireference regions of the potential energy surfaces, they are outperformed by single-reference CC when one electronic configuration dominates.     

A new size extensive multireference perturbation theory

A new multireference perturbation series is derived based on the Rayleigh–Schrödinger perturbation theory. It is orbitally invariant. Its computational cost is comparable to the single reference Møller–Plesset perturbation theory. It is demonstrated numerically that the present multireference second‐ and third‐order energies are size extensive by two types of supermolecules composed of H₂ and BH monomers. Spectroscopic constants of as well as the ground state energies of H₂O, NH₂, and CH₂ at three bond lengths have been calculated with the second multireference perturbation theory. The dissociation behaviors of CH₄ and HF have also been investigated. Comparisons with other approximate theoretical models as well as the experimental data have been carried out to show their relative performances. © 2013 Wiley Periodicals, Inc.