Search

Dimethyl disulfide (DMDS), one of the smallest organic molecules with an S–S bond, serves as a model system for understanding photofragmentation in polypeptides and proteins. Prior studies of DMDS photodissociation excited at ∼266 nm and ∼248 nm have elucidated the mechanisms of S–S and C–S bond cleavage, which involve the lowest excited electronic states S1 and S2. Far less is known about the dissociation mechanisms and electronic structure of relevant excited states of DMDS excited at ∼200 nm. Herein we present calculations of the electronic structure and properties of electronic states S1–S6 accessed when DMDS is excited at ∼200 nm. Our analysis includes a comparison of theoretical and experimental UV spectra, as well as theoretically predicted one-dimensional cuts through the singlet and triplet potential energy surfaces along the S–S and C–S bond dissociation coordinates. Finally, we present calculations of spin–orbit coupling constants at the Franck–Condon geometry to assess the likelihood of ultrafast intersystem crossing. We show that choosing an accurate yet computationally efficient electronic structure method for calculating the S0–S6 potential energy surfaces along relevant dissociation coordinates is challenging due to excited states with doubly excited character and/or mixed Rydberg-valence character. Our findings demonstrate that the extended multi-state complete active space second-order perturbation theory (XMS-CASPT2) balances this computational efficiency and accuracy, as it captures both the Rydberg character of states in the Franck–Condon region and multiconfigurational character toward the bond-dissociation limits. We compare the performance of XMS-CASPT2 to a new variant of equation of motion coupled cluster theory with single, double, and perturbative triple corrections, EOM-CCSD(T)(a)*, finding that EOM-CCSD(T)(a)* significantly improves the treatment of doubly excited states compared to EOM-CCSD, but struggles to quantitatively capture asymptotic energies along bond dissociation coordinates for these states. [ABSTRACT FROM AUTHOR]

The Massively Parallel Quantum Chemistry (MPQC) program is a 30-year-old project that enables facile development of electronic structure methods for molecules for efficient deployment to massively parallel computing architectures. Here, we describe the historical evolution of MPQC's design into its latest (fourth) version, the capabilities and modular architecture of today's MPQC, and how MPQC facilitates rapid composition of new methods as well as its state-of-the-art performance on a variety of commodity and high-end distributed-memory computer platforms. [ABSTRACT FROM AUTHOR]

The Distinguishable Cluster (DC) approximation to the coupled cluster (CC) doubles, proposed by Kats and Manby [J. Chem. Phys. 139, 021102 (2013)], can semiquantitatively describe multiple bond dissociation (which is traditionally considered a paradigm of strongly correlated electronic structure methods based on the multideterminant approach) without leaving the single-reference coupled cluster framework. DC is just one of many internally corrected (CC) methods that improve on the parent CC method by approximation. To build on the success of the DC methods, it is important to probe whether they can be systematically improved. To answer this question, we considered a set of methods in the distinguishable cluster family, culminating in the DC singles, doubles, and triples (DCSDT), a DC modification of the parent CC singles, doubles, and triples that maintains exactness for 3-electron systems. Inclusion of the complete treatment of triples leads to systematic improvement over the DCSD for equilibrium correlation energy estimates. However, this improvement is not matched by enhanced accuracy for multiple bond breaking processes. [ABSTRACT FROM AUTHOR]

Simplifications or modifications of coupled cluster methods such as the CCSD (coupled cluster singles and doubles) model often perform better than the original method in providing the total energy, equilibrium geometries, and harmonic vibration frequencies for the ground state. Three such methods that have been recently proposed include 2CC, parameterized CCSD generalizations, and distinguishable cluster singles and doubles (DCSD) approach. In this paper, we lay the theoretical foundation needed to treat excited states via the equation of motion (EOM) approach using these ground state CC methods. As these ground state approximations to CCSD share its property of being exact for two-electron systems, so will their excited state extensions. These methods are tested for two complementary benchmark sets of excited states for a wide range of organic molecules with focus on singlet and triplet excited states of both valence and Rydberg nature. We also test these methods for doubly excited states, taking CH+ as an example to test their performance at equilibrium and stretched bond geometries. Finally, we assess if any of these methods perform consistently better than EOM CCSD. [ABSTRACT FROM AUTHOR]

Obtaining the correct potential energy curves for the dissociation of multiple bonds is a challenging problem for ab initio methods which are affected by the choice of a spin-restricted reference function. Coupled cluster (CC) methods such as CCSD (coupled cluster singles and doubles model) and CCSD(T) (CCSD + perturbative triples) correctly predict the geometry and properties at equilibrium but the process of bond dissociation, particularly when more than one bond is simultaneously broken, is much more complicated. New modifications of CC theory suggest that the deleterious role of the reference function can be diminished, provided a particular subset of terms is retained in the CC equations. The Distinguishable Cluster (DC) approach of Kats and Manby [J. Chem. Phys. 139, 021102 (2013)], seemingly overcomes the deficiencies for some bond-dissociation problems and might be of use in quasi-degenerate situations in general. DC along with other approximate coupled cluster methods such as ACCD (approximate coupled cluster doubles), ACP-D45, ACP-D14, 2CC, and pCCSD(α, β) (all defined in text) falls under a category of methods that are basically obtained by the deletion of some quadratic terms in the double excitation amplitude equation for CCD/CCSD (coupled cluster doubles model/coupled cluster singles and doubles model). Here these approximate methods, particularly those based on the DC approach, are studied in detail for the nitrogen molecule bond-breaking. The N2 problem is further addressed with conventional single reference methods but based on spatial symmetry-broken restricted Hartree-Fock (HF) solutions to assess the use of these references for correlated calculations in the situation where CC methods using fully symmetry adapted SCF solutions fail. The distinguishable cluster method is generalized: 1) to different orbitals for different spins (unrestricted HF based DCD and DCSD), 2) by adding triples correction perturbatively (DCSD(T)) and iteratively (DCSDT-n), and 3) via an excited state approximation through the equation of motion (EOM) approach (EOM-DCD, EOM-DCSD). The EOM-CC method is used to identify lower-energy CC solutions to overcome singularities in the CC potential energy curves. It is also shown that UHF based CC and DC methods behave very similarly in bond-breaking of N2, and that using spatially broken but spin preserving SCF references makes the CCSD solutions better than those for DCSD. [ABSTRACT FROM AUTHOR]

Modified coupled-cluster (CC) methods such as linearized coupled-cluster doubles (LinCCD), approximate coupled pair (ACP D14), 2CC (from nCC family), parameterized CCSD (pCCSD) and distinguishable cluster (DCSD) can have their advantages over general CC methods. Though these methods include connected clusters of single and double excitations at most, distinguishable cluster, parameterized CC and approximate coupled pair methods, in particular, have been shown to produce quantitatively correct results in benchmark studies. To put these methods on a stronger foothold, it is essential to understand the rationale for their success: mimicking the effect of connected triple excitations. We exploit the relation between CC and many body perturbation theory (MBPT) in general, and between CCSD and MBPT(4)/MP4 in particular, to take a step towards bringing clarity to this persisting conundrum. Our aim here is to look for numerical signs of 'addition by subtraction' or 'inclusion by deletion' effect that is likely behind the success of these modified CCD or CCSD methods. We achieve this by revisiting well-studied examples of single and multiple bond dissociation and comparing the performance of these modified CCSD methods with higher-level CC methods. Though our results are qualitative in nature, we hope this would lead to more rigorous analysis in future studies. [ABSTRACT FROM AUTHOR]