Javier Segarra-Martí

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Ab initio determination of the ionization potentials of water clusters (H2O)n (n = 2-6).

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic ionization potentials of several water clusters: dimer, trimer, tetramer, pentamer, hexamer book, hexamer ring, hexamer cage, and hexamer prism. The present results establish reference values at a level not reported before for these systems, calibrating different computational strategies and helping to discard less reliable theoretical and experimental data. The systematic study with the increasing size of the water cluster allows obtaining some clues on the structure and reductive properties of liquid water.

Modeling the high-energy electronic state manifold of adenine: Calibration for nonlinear electronic spectroscopy

Pump-probe electronic spectroscopy using femtosecond laser pulses has evolved into a standard tool for tracking ultrafast excited state dynamics. Its two-dimensional (2D) counterpart is becoming an increasingly available and promising technique for resolving many of the limitations of pump-probe caused by spectral congestion. The ability to simulate pump-probe and 2D spectra from ab initio computations would allow one to link mechanistic observables like molecular motions and the making/breaking of chemical bonds to experimental observables like excited state lifetimes and quantum yields. From a theoretical standpoint, the characterization of the electronic transitions in the visible (Vis)/ultraviolet (UV), which are excited via the interaction of a molecular system with the incoming pump/probe pulses, translates into the determination of a computationally challenging number of excited states (going over 100) even for small/medium sized systems. A protocol is therefore required to evaluate the fluctuations of spectral properties like transition energies and dipole moments as a function of the computational parameters and to estimate the effect of these fluctuations on the transient spectral appearance. In the present contribution such a protocol is presented within the framework of complete and restricted active space self-consistent field theory and its second-order perturbation theory extensions. The electronic excited states of adenine have been carefully characterized through a previously presented computational recipe [Nenov et al., Comput. Theor. Chem. 1040-1041, 295-303 (2014)]. A wise reduction of the level of theory has then been performed in order to obtain a computationally less demanding approach that is still able to reproduce the characteristic features of the reference data. Foreseeing the potentiality of 2D electronic spectroscopy to track polynucleotide ground and excited state dynamics, and in particular its expected ability to provide conformational dependent fingerprints in dimeric systems, the performances of the selected reduced level of calculations have been tested in the construction of 2D electronic spectra for the in vacuo adenine monomer and the unstacked adenine homodimer, thereby exciting the Lb/La transitions with the pump pulse pair and probing in the Vis to near ultraviolet spectral window.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems.

A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.

On the photophysics and photochemistry of the water dimer.

The photochemistry of the water dimer irradiated by UV light is studied by means of the complete active space perturbation theory//complete active space self-consistent field (CASPT2//CASSCF) method and accurate computational approaches like as minimum energy paths. Both electronic structure computations and ab initio molecular dynamics simulations are carried out. The results obtained show small shifts relative to a single water molecule on the vertical excitation energies of the dimer due to the hydrogen bond placed between the water donor (W(D)) and the water acceptor (W(A)). A red-shift and a blue-shift are predicted for the W(D) and W(A), respectively, supporting previous theoretical and experimental results. The photoinduced chemistry of the water dimer is described as a process occurring between two single water molecules in which the effect of the hydrogen bond plays a minor role. Thus, the photoinduced decay routes correspond to two photodissociation processes, one for each water molecule. The proposed mechanism for the decay channels of the lowest-lying excited states of the system is established as the photochemical production of a hydrogen-bonded H(2)O...HO species plus a hydrogen H atom.

Ultraviolet vision: photophysical properties of the unprotonated retinyl Schiff base in the Siberian hamster cone pigment

Abstract The Siberian hamster ultraviolet (SHUV) visual pigment has an unprotonated Schiff-base (SB) retinyl chromophore in the dark state, which becomes protonated after photoexcitation during the early stages of the photobleaching cycle. While the photochemical relaxation processes of the SHUV remain poorly understood, they are expected to show significant differences when compared to those of the protonated SB (PSB) chromophore in visual rhodopsin. Here, we report a study of the photophysical properties of the SHUV unprotonated SB (SHUV-USB), based on multiconfigurational and multireference perturbative methods within a hybrid quantum mechanics/molecular mechanics scheme. Comparisons of multireference and time-dependent density functional theory results indicate that both methodologies predict an ionic excited state (S1), similar to the PSB of rhodopsin, although its minimum has even bond-lengths in the central region of the retinyl polyene chain. The analysis of excited-state manifolds at the Franck–Condon region and S1 minimum configuration indicates that the skeletal relaxation initiated in the S1 surface is likely to involve S1/S2 surface crossing. These results provide valuable insights for future studies of the SHUV-USB photoisomerization mechanism.

On the N1-H and N3-H Bond Dissociation in Uracil by Low Energy Electrons: A CASSCF/CASPT2 Study.

The dissociative electron-attachment (DEA) phenomena at the N1-H and N3-H bonds observed experimentally at low energies (<3 eV) in uracil are studied with the CASSCF/CASPT2 methodology. Two valence-bound π(-) and two dissociative σ(-) states of the uracil anionic species, together with the ground state of the neutral molecule, are proven to contribute to the shapes appearing in the experimental DEA cross sections. Conical intersections (CI) between the π(-) and σ(-) are established as the structures which activate the DEA processes. The N1-H and N3-H DEA mechanisms in uracil are described, and experimental observations are interpreted on the basis of two factors: (1) the relative energy of the (U-H)(-) + H fragments obtained after DEA with respect to the ground-state equilibrium structure (S0) of the neutral molecule (threshold for DEA) and (2) the relative energy of the CIs also with respect to S0 (band maxima). The π1(-) state is found to be mainly responsible for the N1-H bond breaking, whereas the π2(-) state is proved to be involved in the cleavage of the N3-H bond.

Multiple Decay Mechanisms and 2D-UV Spectroscopic Fingerprints of Singlet Excited Solvated Adenine-Uracil Monophosphate.

Abstract The decay channels of singlet excited adenine uracil monophosphate (ApU) in water are studied with CASPT2//CASSCF:MM potential energy calculations and simulation of the 2D‐UV spectroscopic fingerprints with the aim of elucidating the role of the different electronic states of the stacked conformer in the excited state dynamics. The adenine 1La state can decay without a barrier to a conical intersection with the ground state. In contrast, the adenine 1Lb and uracil S(U) states have minima that are separated from the intersections by sizeable barriers. Depending on the backbone conformation, the CT state can undergo inter‐base hydrogen transfer and decay to the ground state through a conical intersection, or it can yield a long‐lived minimum stabilized by a hydrogen bond between the two ribose rings. This suggests that the 1Lb, S(U) and CT states of the stacked conformer may all contribute to the experimental lifetimes of 18 and 240 ps. We have also simulated the time evolution of the 2D‐UV spectra and provide the specific fingerprint of each species in a recommended probe window between 25 000 and 38 000 cm−1 in which decongested, clearly distinguishable spectra can be obtained. This is expected to allow the mechanistic scenarios to be discerned in the near future with the help of the corresponding experiments. Our results reveal the complexity of the photophysics of the relatively small ApU system, and the potential of 2D‐UV spectroscopy to disentangle the photophysics of multichromophoric systems.

Deciphering the photochemical mechanisms describing the UV-induced processes occurring in solvated guanine monophosphate.

The photophysics and photochemistry of water-solvated guanine monophosphate (GMP) are here characterized by means of a multireference quantum-chemical/molecular mechanics theoretical approach (CASPT2//CASSCF/AMBER) in order to elucidate the main photo-processes occurring upon UV-light irradiation. The effect of the solvent and of the phosphate group on the energetics and structural features of this system are evaluated for the first time employing high-level ab initio methods and thoroughly compared to those in vacuo previously reported in the literature and to the experimental evidence to assess to which extent they influence the photoinduced mechanisms. Solvated electronic excitation energies of solvated GMP at the Franck-Condon (FC) region show a red shift for the ππ* La and Lb states, whereas the energy of the oxygen lone-pair nπ* state is blue-shifted. The main photoinduced decay route is promoted through a ring-puckering motion along the bright lowest-lying La state toward a conical intersection (CI) with the ground state, involving a very shallow stationary point along the minimum energy pathway in contrast to the barrierless profile found in gas-phase, the point being placed at the end of the minimum energy path (MEP) thus endorsing its ultrafast deactivation in accordance with time-resolved transient and photoelectron spectroscopy experiments. The role of the nπ* state in the solvated system is severely diminished as the crossings with the initially populated La state and also with the Lb state are placed too high energetically to partake prominently in the deactivation photo-process. The proposed mechanism present in solvated and in vacuo DNA/RNA chromophores validates the intrinsic photostability mechanism through CI-mediated non-radiative processes accompanying the bright excited-state population toward the ground state and subsequent relaxation back to the FC region.

Resolving Ultrafast Photoinduced Deactivations in Water-Solvated Pyrimidine Nucleosides

For the first time, ultrafast deactivations of photoexcited water-solvated pyrimidine nucleosides are mapped employing hybrid QM(CASPT2)/MM(AMBER) optimizations that account for explicit solvation, sugar effects, and dynamically correlated potential energy surfaces. Low-energy S1/S0 ring-puckering and ring-opening conical intersections (CIs) are suggested to drive the ballistic coherent subpicosecond (<200 fs) decays observed in each pyrimidine, the energetics controlling this processes correlating with the lifetimes observed. A second bright 1π2π* state, promoting excited-state population branching and leading toward a third CI with the ground state, is proposed to be involved in the slower ultrafast decay component observed in Thd/Cyd. The transient spectroscopic signals of the competitive deactivation channels are computed for the first time. A general unified scheme for ultrafast deactivations, spanning the sub- to few-picosecond time domain, is eventually delivered, with computed data that matches the experiments and elucidates the intrinsic photoprotection mechanism in solvated pyrimidine nucleosides.