Conical Intersection And Potential Energy Surface Features Of A Model Retinal Chromophore: Comparison Of Eom-cc And Multireference Methods

This work investigates the performance of equation-of-motion coupled-cluster (EOM-CC) methods for describing the changes in the potential energy surfaces of the penta-2,4-dieniminium cation, a reduced model of the retinal chromophore of visual pigments, due to dynamical electron correlation effects. The ground-state wave function of this model includes charge-transfer and diradical configurations whose weights vary along different displacements and are rapidly changing at the conical intersection between the ground and the first excited states, making the shape of the potential energy surface sensitive to a balanced description of nondynamical and dynamical correlation. Recently, variational (MRCISD) and perturbative (MRPT2) approaches for including dynamical correlation in CASSCF-based calculations were tested along three representative ground state paths. Here, we use the same three paths to compare the performance of single-reference EOM-CC methods against MRCISD and MRCISD+Q We find that the spin-flip variant of EOM-CCSD with perturbative inclusion of triple excitations (dT or IT) produces potential energy profiles of the two lowest electronic states in quantitative agreement with MRCISD+Q (our highest-quality reference method). The nonparallelity errors and differences in vertical energy differences of the two surfaces along these scans are less than 1.4 kcal/mol (EOM-SF-CCSD(dT) versus MRCISD+Q). For comparison, the largest error of MRCISD versus MRCISD+Q is 1.7 kcal/mol. Our results show that the EOM-CC methods provide an alternative to multireference approaches and may be used to study photochemical systems like the one used in this work

To Be or Not To Be a Molecular Ion: The Role of the Solvent in Photoionization of Arginine.

Application of photoionization mass spectroscopy, a technique capable of assessing protonation states in complex molecules in the gas phase, is challenging for arginine due to its fragility. We report photoionization efficiencies in the valence region of aqueous aerosol particles produced from arginine solutions under various pH and vaporization conditions. By using ab initio calculations, we investigate the stability of different conformers. Our results show that neutral arginine fragments upon ionization in the gas phase but solvation stabilizes the molecular ion, resulting in different photoionization dynamics. We also report the valence-band photoelectron spectra of the aerosol solutions obtained at different pH values

Two-photon absorption spectroscopy of stilbene and phenanthrene: Excited-state analysis and comparison with ethylene and toluene

The following article appeared in The Journal of Chemical Physics 146, 174102 (2017); doi: 10.1063/1.498204 and may be found at http://doi.org/10.1063/1.4982045Two-photon absorption (2PA) spectra of several prototypical molecules (ethylene, toluene, trans- and cis-stilbene, and phenanthrene) are computed using the equation-of-motion coupled-cluster method with single and double substitutions. The states giving rise to the largest 2PA cross sections are analyzed in terms of their orbital character and symmetry-based selection rules. The brightest 2PAtransitions correspond to Rydberg-like states from fully symmetric irreducible representations. Symmetry selection rules dictate that totally symmetric transitions typically have the largest 2PA cross sections for an orientationally averaged sample when there is no resonance enhancement via one-photon accessible intermediate states. Transition dipole arguments suggest that the strongest transitions also involve the most delocalized orbitals, including Rydberg states, for which the relative transition intensities can be rationalized in terms of atomic selection rules. Analysis of the 2PA transitions provides a foundation for predicting relative 2PA cross sections of conjugated molecules based on simple symmetry and molecular orbital arguments

Using the charge-stabilization technique in the double ionization potential equation-of-motion calculations with dianion references

The charge-stabilization method is applied to double ionization potential equation-of-motion (EOM-DIP) calculations to stabilize unstable dianion reference functions. The auto-ionizing character of the dianionic reference states spoils the numeric performance of EOM-DIP limiting applications of this method. We demonstrate that reliable excitation energies can be computed by EOM-DIP using a stabilized resonance wave function instead of the lowest energy solution corresponding to the neutral + free electron(s) state of the system. The details of charge-stabilization procedure are discussed and illustrated by examples. The choice of optimal stabilizing Coulomb potential, which is strong enough to stabilize the dianion reference, yet, minimally perturbs the target states of the neutral, is the crux of the approach. Two algorithms of choosing optimal parameters of the stabilization potential are presented. One is based on the orbital energies, and another -on the basis set dependence of the total Hartree-Fock energy of the reference. Our benchmark calculations of the singlet-triplet energy gaps in several diradicals show a remarkable improvement of the EOM-DIP accuracy in problematic cases. Overall, the excitation energies in diradicals computed using the stabilized EOM-DIP are within 0.2 eV from the reference EOM spin-flip values

Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework:Theory, implementation, and validation against multi-reference methods

We report an implementation of non-adiabatic coupling (NAC) forces within the equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) framework via the summed-state approach. Using illustrative examples, we compare NAC forces computed with EOM-CCSD and multi-reference (MR) wave functions (for selected cases, we also consider configuration interaction singles). In addition to the magnitude of the NAC vectors, we analyze their direction, which is important for the calculations of the rate of non-adiabatic transitions. Our benchmark set comprises three doublet radical-cations (hexatriene, cyclohexadiene, and uracil), neutral uracil, and sodium-doped ammonia clusters. When the characters of the states agree among different methods, we observe good agreement between the respective NAC vectors, both in the Franck-Condon region and away. In the cases of large discrepancies between the methods, the disagreement can be attributed to the difference in the states' character, which, in some cases, is very sensitive to electron correlation, both within single-reference and multi-reference frameworks. The numeric results confirm that the accuracy of NAC vectors depends critically on the quality of the underlying wave functions. Within their domain of applicability, EOM-CC methods provide a viable alternative to MR approaches

Analytic evaluation of non-adiabatic couplings within the complex absorbing potential equation-of-motion coupled-cluster method

We present the theory for the evaluation of non-adiabatic couplings (NACs) involving resonance states within the complex absorbing potential equation-of-motion coupled-cluster (CAP-EOM-CC) framework implemented within the singles and doubles approximation. Resonance states are embedded in the continuum and undergo rapid decay through autodetachment. In addition, nuclear motions can facilitate transitions between different resonances and between resonances and bound states. These non-adiabatic transitions affect the chemical fate of resonances and have distinct spectroscopic signatures. The NAC vector is a central quantity needed to model such effects. In the CAP-EOM-CC framework, resonance states are treated on the same footing as bound states. Using the example of fumaronitrile, which supports a bound radical anion and several anionic resonances, we analyze the non-adiabatic coupling between bound states and pseudocontinuum states, between bound states and resonances and between two resonances. We find that the NAC between a bound state and a resonance is nearly independent of the CAP strength and thus straightforward to evaluate whereas the NAC between two resonance states or between a bound state and a pseudocontinuum state is more difficult to evaluate

Two-photon absorption spectroscopy of trans-stilbene, cis-stilbene, and phenanthrene: Theory and experiment

The following article appeared in The Journal of Chemical Physics 146, 144305 (2017); doi: 10.1063/1.4979651 and may be found at http://dx.doi.org/10.1063/1.4979651Two-photon absorption (2PA) spectroscopy provides complementary, and sometimes more detailed, information about the electronic structure of a molecule relative to one-photon absorption (1PA) spectroscopy. However, our understanding of the 2PA processes is rather limited due to technical difficulties in measuring experimental 2PA spectra and theoretical challenges in computing higher-order molecular properties. This paper examines the 2PA spectroscopy of trans-stilbene, cis-stilbene, and phenanthrene by a combined experimental and theoretical approach. The broadband 2PA spectra of all three compounds are measured under identical conditions in order to facilitate a direct comparison of the absolute 2PA cross sections in the range 3.5-6.0 eV. For comparison, the theoretical 2PA cross sections are computed using the equation-of-motion coupled-cluster method with single and double substitutions. Simulated 2PA spectra based on the calculations reproduce the main features of the experimental spectra in solution, although the quantitative comparison is complicated by a number of uncertainties, including limitations of the theoretical model, vibronic structure, broadening of the experimental spectra, and solvent effects. The systematic comparison of experimental and theoretical spectra for this series of structurally similar compounds provides valuable insight into the nature of 2PA transitions in conjugated molecules. Notably, the orbital character and symmetry-based selection rules provide a foundation for interpreting the features of the experimental 2PA spectra in unprecedented detail

Mapping The Excited State Potential Energy Surface Of A Retinal Chromophore Model With Multireference And Equation-of-motion Coupled-cluster Methods

The photoisomerization of the retinal chromophore of visual pigments proceeds along a complex reaction coordinate on a multidimensional surface that comprises a hydrogen-out-of-plane (HOOP) coordinate, a bond length alternation (BLA) coordinate, a single bond torsion and, finally, the reactive double bond torsion. These degrees of freedom are coupled with changes in the electronic structure of the chromophore and, therefore, the computational investigation of the photochemistry of such systems requires the use of a methodology capable of describing electronic structure changes along all those coordinates. Here, we employ the penta-2,4-dieniminium (PSB3) cation as a minimal model of the retinal chromophore of visual pigments and compare its excited state isomerization paths at the CASSCF and CASPT2 levels of theory. These paths connect the cis isomer and the trans isomer of PSB3 with two structurally and energetically distinct conical intersections (CIs) that belong to the same intersection space. MRCISD+Q energy profiles along these paths provide benchmark values against which other ab initio methods are validated. Accordingly, we compare the energy profiles of MRPT2 methods (CASPT2, QD-NEVPT2, and XMCQDPT2) and EOM-SF-CC methods (EOM-SF-CCSD and EOM-SF-CCSD(dT)) to the MRCISD+Q reference profiles. We find that the paths produced with CASSCF and CASPT2 are topologically and energetically different, partially due to the existence of a locally excited region on the CASPT2 excited state near the Franck-Condon point that is absent in CASSCF and that involves a single bond, rather than double bond, torsion. We also find that MRPT2 methods as well as EOM-SF-CCSD(dT) are capable of quantitatively describing the processes involved in the photoisomerization of systems like PSB3

The effect of pi-stacking, h-bonding, and electrostatic interactions on the ionization energies of nucleic acid bases: adenine-adenine, thymine-thymine and adenine-thymine dimers

A combined theoretical and experimental study of the ionized dimers of thymine and adenine, TT, AA, and AT, is presented. Adiabatic and vertical ionization energies(IEs) for monomers and dimers as well as thresholds for the appearance of the protonated species are reported and analyzed. Non-covalent interactions stronglyaffect the observed IEs. The magnitude and the nature of the effect is different for different isomers of the dimers. The computations reveal that for TT, the largestchanges in vertical IEs (0.4 eV) occur in asymmetric h-bonded and symmetric pi- stacked isomers, whereas in the lowest-energy symmetric h-bonded dimer the shiftin IEs is much smaller (0.1 eV). The origin of the shift and the character of the ionized states is different in asymmetric h-bonded and symmetric stacked isomers. Inthe former, the initial hole is localized on one of the fragments, and the shift is due to the electrostatic stabilization of the positive charge of the ionized fragment by thedipole moment of the neutral fragment. In the latter, the hole is delocalized, and the change in IE is proportional to the overlap of the fragments' MOs. The shifts in AAare much smaller due to a less effcient overlap and a smaller dipole moment. The ionization of the h-bonded dimers results in barrierless (or nearly barrierless) protontransfer, whereas the pi-stacked dimers relax to structures with the hole stabilized by the delocalization or electrostatic interactions