Explicitly time-dependent coupled cluster singles doubles calculations of laser-driven many-electron dynamics

We report explicitly time-dependent coupled cluster singles doubles (TD-CCSD) calculations, which simulate the laser-driven correlated many-electron dynamics in molecular systems. Small molecules, i.e., HF, H(2)O, NH(3), and CH(4), are treated mostly with polarized valence double zeta basis sets. We determine the coupled cluster ground states by imaginary time propagation for these molecules. Excited state energies are obtained from the Fourier transform of the time-dependent dipole moment after an ultrashort, broadband laser excitation. The time-dependent expectation values are calculated from the complex cluster amplitudes using the corresponding configuration interaction singles doubles wave functions. Also resonant laser excitations of these excited states are simulated, in order to explore the limits for the numerical stability of our current TD-CCSD implementation, which uses time-independent molecular orbitals to form excited configurations.     

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).     

An overview of coupled cluster theory and its applications in physics

Summary What has since become known as the normal coupled cluster method (NCCM) was invented about thirty years ago to calculate ground-state energies of closed-shell atomic nuclei. Coupled cluster (CC) techniques have since been developed to calculate excited states, energies of open-shell systems, density matrices and hence other properties, sum rules, and the sub-sum-rules that follow from imbedding linear response theory within the NCCM. Further extensions deal both with systems at nonzero temperature and with general dynamical behaviour. More recently, a new version of CC theory, the so-called extended coupled cluster method (ECCM) has been introduced. It has the potential to describe such global phenomena as phase transitions, spontaneous symmetry breaking, states of topological excitation, and nonequilibrium behaviour. CC techniques are now widely recognized as providing one of the most universally applicable, most powerful, and most accurate of all microscopicab initio methods in quantum many-body theory. The number of successful applications within physics is now impressively large. In most such cases the numerical results are either the best or among the best available. A typical case is the electron gas, where the CC results for the correlation energy agree over the entire metallic density range to within less than 1 millihartree (or <1%) with the essentially exact Green's function Monte Carlo results. The role of CC theory within modern quantum many-body theory is first surveyed, by a comparison with other techniques. Its full range of applications in physics is then reviewed. These include problems in nuclear physics, both for finite nuclei and infinite nuclear matter; the electron gas; various integrable and nonintegrable models; various relativistic quantum field theories; and quantum spin chain and lattice models. Particular applications of the ECCM include the quantum hydrodynamics of a zero-temperature, strongly-interacting condensed Bose fluid; a charged impurity in a polarizable medium (e.g., positron annihilation in metals); and various anharmonic oscillator and spin systems.     

Theoretical study on spin alignments in ferromagnetic heterospin chains with competing exchange interactions: a generalized ferrimagnetic system containing organic biradicals in the singlet ground state

Spin alignments in heterospin chains are examined from numerical calculations of model spin Hamiltonians. The Hamiltonians of the heterospin chains mimic an open-shell molecular assemblage composed of an organic biradical in a singlet (S = 0) ground state and a doublet (S = 1/2) monoradical, which are coupled by intermolecular ferromagnetic exchange interactions. It is found from numerical calculations of the spin Hamiltonians that the spin value S2 of the ground-state singlet biradical embedded in the exchange-coupled assemblage deviates from zero and contributes to the bulk magnetization. The alternating chain is found to have two kinds of ground spin states, a high- and a low-spin state. All the spins are parallel to each other in the high-spin state, which is characterized by the spin correlation function of (S(i).S(j)) = 0.25. On the other hand, the spin alignment in the low-spin state is found to be dependent on the topology of the intermolecular exchange interactions. The energy preference of the two states depends on the relative amplitude of the exchange interactions in the chain. The intermolecular ferromagnetic couplings are competing in the ground-state singlet biradical with the intramolecular antiferromagnetic interaction. The appearance of the two kinds of ground states is attributed to a quantum spin frustration effect inherent in the triangular motif of the competing interactions. Magnetic properties of a zigzag chain complex composed of a nitronyl nitroxide biradical with a singlet ground state and Cu(hfac)2 are examined on the basis of the theoretical calculations. The vanishing magnetic moments, or the product of susceptibility and temperature chiT, at low temperatures observed for the complex are consistent with those of the low-spin state predicted in the theoretical calculations.     

Two-dimensional H2O-Cl2 and H2O-Br2 potential surfaces: an ab initio study of ground and valence excited electronic states

All electron ab initio calculations for the interaction of H2O with Cl2 and Br2 are reported for the ground state and the lowest triplet and singlet Pi excited states as a function of both the X-X and O-X bond lengths (X = Cl or Br). For the ground state and lowest triplet state, the calculations are performed with the coupled cluster singles, doubles, and perturbative triple excitation level of correlation using an augmented triple-zeta basis set. For the 1Pi state the multireference average quadratic coupled cluster technique was employed. For several points on the potential, the calculations were repeated with the augmented quadruple-zeta basis set. The ground-state well depths were found to be 917 and 1,183 cm-1 for Cl2 and Br2, respectively, with the triple-zeta basis set, and they increased to 982 and 1,273 cm-1 for the quadruple-zeta basis set. At the geometry of the ground-state minimum, the lowest energy state corresponding to the unperturbed 1Pi states of the halogens increases in energy by 637 and 733 cm-1, respectively, relative to the ground-state dissociation limit of the H2O-X2 complex. Adding the attractive ground-state interaction energy to that of the repulsive excited state predicts a blue-shift, relative to that of the free halogen molecules, of approximately 1,600 cm-1 for H2O-Cl2 and approximately 2,000 cm-1 for H2O-Br2. These vertical blue-shifts for the dimers are greater than the shift of the band maximum upon solvation of either halogen in liquid water.     

The fluorobenzene-argon S(1) excited-state intermolecular potential energy surface

We evaluate the first excited-state (S1) intermolecular potential energy surface for the fluorobenzene-Ar van der Waals complex using the coupled cluster method and the augmented correlation-consistent polarized valence double-zeta basis set extended with a set of 3s3p2d1f1g midbond functions. To calculate the S(1) interaction energies, we use ground-state interaction energies evaluated with the same basis set and the coupled cluster singles and doubles (CCSD) including connected triple excitations [CCSD(T)] model and interaction and excitation energies evaluated at the CCSD level. The surface minima are characterized by the Ar atom located above and below the fluorobenzene ring at a distance of 3.5060 A with respect to the fluorobenzene center of mass and at an angle of 5.89 degrees with respect to the axis perpendicular to the fluorobenzene plane. The corresponding interaction energy is -425.226 cm(-1). The surface is used in the evaluation of the intermolecular level structure of the complex, and the results are compared to the experimental data available and to those found in previous theoretical papers on ground-state potentials for similar complexes.     

Quantum-chemical investigation of the structures and electronic spectra of the nucleic acid bases at the coupled cluster CC2 level

We study the ground-state structures and singlet- and triplet-excited states of the nucleic acid bases by applying the coupled cluster model CC2 in combination with a resolution-of-the-identity approximation for electron interaction integrals. Both basis set effects and the influence of dynamic electron correlation on the molecular structures are elucidated; the latter by comparing CC2 with Hartree-Fock and M?ller-Plesset perturbation theory to second order. Furthermore, we investigate basis set and electron correlation effects on the vertical excitation energies and compare our highest-level results with experiment and other theoretical approaches. It is shown that small basis sets are insufficient for obtaining accurate results for excited states of these molecules and that the CC2 approach to dynamic electron correlation is a reliable and efficient tool for electronic structure calculations on medium-sized molecules.     

The coupled-cluster method with a multiconfiguration reference state

The coupled-cluster approach to obtaining the bond-state wave functions of many-electron systems is extended, with a set of physically reasonable approximations, to admit a multiconfiguration reference state. This extension permits electronic structure calculations to be performed on correlated closed- or open-shell systems with potentially uniform precision for all molecular geometries. Explicit coupled cluster working equations are derived using a multiconfiguration reference state for the case in which the so-called cluster operator is approximated by its one- and two-particle components. The evaluation of the requisite matrix elements is facilitated by use of the unitary group generators which have recently received wide attention and use in the quantum chemistry community.     

The symmetrized basis function method for SU(2) and its application to the total spin states of many-electron systems

In this paper, the symmetrized basis function method is extended and used for the classification of the total spin states of many-electron systems. The reduction matrix, which is expressed as a series of products and direct products of matrices is derived. It is very advantageous that the method can be completed by computer. As an example, the reduction matrix for a 5-electron system calculated by IBM PC computer and the classification of the total spin states is presented.     

基于CASSCF参考函数的块相关耦合簇方法对烷烃中单键解离势能面研究

用以完全活化空间自洽场(CASSCF)波函数为参考波函数的块相关耦合簇(BCCC)方法(简称CAS-BCCC)研究了烷烃(甲烷和乙烷)中的单键解离过程的势能面(PES). 与其它理论方法比较的结果表明, 该方法可以对所研究的整个解离势能面给出定量准确的描述.     

Towards multistate multimode landscapes in singlet fission of pentacene: the dual role of charge-transfer states

Singlet fission duplicates triplet excitons for improving light harvesting efficiency. The presence of the interaction between electronic and nuclear degrees of freedom complicates the interpretation of correlated triplet pairs. We report a quantum chemistry study on the significance and subtleties of multistate and multimode pathways in forming triplet pair states of the pentacene dimer through a six-state vibronic-coupling Hamiltonian derived from many-electron adiabatic wavefunctions of an ab initio density matrix renormalization group. The resulting spin values of the singlet manifolds on each pentacene center are computed, and the varying spin nature can be distinguished clearly with respect to dimer stacking and vibronic progression. Our monomer spin assignments reveal the coexistence of both lower-lying weak and higher-lying strong charge transfer states which interact vibronically with the triplet pair state, providing important implications for its generation and separation occurring in vibronic regions. This work conveys the importance of the many-electron process requiring close low-lying singlet manifolds to determine the subtle fission details, and represents an important step for understanding vibronically resolved spin states and conversions underlying efficient singlet fission.

Singlet fission in pentacene necessitates the vibronic progression of weak and strong charge-transfer states with correlated triplet pairs.     

Density functional theory: Approximating quasiparticle density functional

A constructive approach for deriving the approximating quasiparticle energy density functional is proposed. As a matter of fact, the proposed approach is the direct development of the Kohn–Sham quasiparticle concept and the Levy–Valone approach. The approach presented takes into account a pseudopotential character of the exchange-correlation part of the density functional and results in a system of functional equations to obtain ground-state energies of many-electron systems.     

Effective local potentials for excited states

The constrained variational Hartree-Fock method for excited states of the same symmetry as the ground state [Chem. Phys. Lett. 287, 189 (1998)] is combined with the effective local potential (ELP) method [J. Chem. Phys. 125, 081104 (2006)] to generate Kohn-Sham-type exact-exchange potentials for singly excited states of many-electron systems. Illustrative examples include the three lowest (2)S states of the Li and Na atoms and the three lowest (3)S states of He and Be. For the systems studied, excited-state ELPs differ from the corresponding ground-state potentials in two respects: They are less negative and have small additional "bumps" in the outer electron region. The technique is general and can be used to approximate excited-state exchange-correlation potentials for other orbital-dependent functionals.     

Electron correlation studies: Rumer basis approach

A technique for the configuration interaction (CI) study of many-electron systems is developed based on Rumer spin-coupling scheme for the antisymmetrized configuration state functions (CSF). Incorporating a new graphical approach, the primitive configurations have been generated in blocks of definite ionocities to permit ready association of possible spin functions with each of the primitive configurations. Simple as well as extended Hubard model Hamiltonians have been studied to test the efficiency of the method. Procedures have been incorporated to calculate various correlation functions using the spin-adapted CSFs without invoking explicit expansions in terms of slater determinants. © 1996 John Wiley & Sons, Inc.     

On the triplet ground state of tetrahedral X4 clusters (X = Li, Na, K, Cu)

The lowest electronic state of distorted tetrahedral X(4) clusters (with X = Li, Na, K, Cu) is studied at coupled-cluster level using high-quality atomic basis sets. The ground state is found to have a triplet spin symmetry for this kind of geometry and for all the considered atomic species. The equilibrium geometries correspond to Jahn-Teller-distorted oblate tetrahedra having D(2d) symmetry, and tetrahedric structures are local minima on the potential-energy surfaces for the triplet states. Their energies lie between 0.2 eV (for the K(4) cluster) and 0.9 eV (for Cu(4)) above the absolute minimum of the corresponding systems, which is a spin singlet having a rhombus geometry.     

Electronic structures, vibrational and thermochemical properties of neutral and charged niobium clusters Nb(n), n = 7-12

Geometric and electronic structures, vibrational properties, and relative stabilities of niobium clusters Nb(n), n = 7-12, are studied using both DFT (BPW91 and M06 functionals) and CCSD(T) calculations with the cc-pVnZ-PP basis set. In each cluster, various lower-lying states are very close in energy in such a way that the ground state cannot be unambiguously established by DFT computations. Nb clusters tend to prefer the lowest possible spin state as the ground state, except for Nb(12) ((3)A(g)). The optimal structure of the cluster at a certain size does not simply grow from that of the smaller one by adding an atom randomly. Instead, the Nb clusters prefer a close-packed growth behavior. Nb(10) has a spherically aromatic character, high chemical hardness and large HOMO-LUMO gap. Electron affinities, ionization energies, binding energy per atom, and the stepwise dissociation energies are evaluated. Energetic properties exhibit odd-even oscillations. Comparison with experimental values shows that both BPW91 and M06 functionals are reliable in predicting the EA and IE values, but the BPW91 is deficient in predicting the binding and dissociation energies. We re-examine in particular the experimental far IR spectra previously recorded using the IR-MPD and free electron laser spectrometric techniques and propose novel assignments for Nb(7) and Nb(9) systems. The IR spectra of the anions are also predicted.     

Treating molecules in arbitrary spin states using the parametric two-electron reduced-density-matrix method

Minimizing the electronic energy with respect to a parameterized two-electron reduced density matrix (2-RDM) is known as a parametric variational 2-RDM method. The parametric 2-RDM method with the M 2-RDM parametrization [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)] is extended to treat molecules in arbitrary spin states. Like its singlet counterpart, the M parametric 2-RDM method for arbitrary spin states is derived using approximate N-representability conditions, which allow it to capture more correlation energy than coupled cluster with single and double excitations at a lower computational cost. We present energies, optimized bond lengths, potential energy curves, and occupation numbers for a set of molecules in a variety of spin states using the M and K parametric 2-RDM methods as well as several wavefunction methods. We show that the M parametric 2-RDM method can describe bond breaking of open-shell molecules like triplet B(2) and singlet and triplet OH(+) even in the presence of strong correlation. Finally, the computed 2-RDMs are shown to be nearly N-representable at both equilibrium and non-equilibrium geometries.