Luis Serrano-Andrés

Software news and update MOLCAS 7 : The Next Generation

Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two‐electron integrals and in the generation of so‐called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self‐consistent field, density functional theory, 2nd order perturbation theory, complete‐active space self‐consistent field multiconfigurational reference 2nd order perturbation theory, and coupled‐cluster methods. The report further elaborates on the implementation of a restricted‐active space self‐consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas‐Kroll‐Hess transformation for one‐component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so‐called picture‐change‐free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.

The multi-state CASPT2 method

Abstract An extension of the multiconfigurational second-order perturbation approach CASPT2 is suggested, where several electronic states are coupled at second order via an effective-Hamiltonian approach. The method has been implemented into the MOLCAS-4 program system, where it will replace the single-state CASPT2 program. The accuracy of the method is illustrated through calculations of the ionic-neutral avoided crossing in the potential curves for LiF and of the valence-Rydberg mixing in the V-state of the ethylene molecule.

Does density functional theory contribute to the understanding of excited states of unsaturated organic compounds

A comparative study has been performed on the electronic spectra of a number of unsaturated organic molecules, using on the one hand density functional linear response theory and on the other multiconfigurational second-order perturbation theory, in order to establish the accuracy that the density functional based methods can give for excitation energies and energy surfaces for excited states. The following molecules are included in the study: tetrazine; the five-membered ring systems cyclopentadiene, furan, pyrrole, and thiophene; acetone; and a dipeptide. The results show that DFT valence excited states have errors that vary between 0 and 1 eV, while Rydberg states are accurate to about 0.2eV in most cases. The use of an asymptotically corrected exchange-correlation potential was essential for the latter result. However, transitions which involve a considerable charge transfer have much larger errors. The results show in some cases a surprisingly strong interaction between valence and Rydberg excited st...

Towards an accurate molecular orbital theory for excited states: Ethene, butadiene, and hexatriene

A newly proposed quantum chemical approach for ab initio calculations of electronic spectra of molecular systems is applied to the molecules ethene, trans-1,3-butadiene, and trans-trans-1,3,5-hexat ...

A Three-State Model for the Photophysics of Guanine

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated 1(pi pi* L(a)) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L(a))CI. On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, 1(pi pi* L(b)) and, in particular, 1(n(O) pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.

Computation of conical intersections by using perturbation techniques

Multiconfigurational second-order perturbation theory, both in its single-state multiconfigurational second-order perturbation theory (CASPT2) and multistate (MS-CASPT2) formulations, is used to search for minima on the crossing seams between different potential energy hypersurfaces of electronic states in several molecular systems. The performance of the procedures is tested and discussed, focusing on the problem of the nonorthogonality of the single-state perturbative solutions. In different cases the obtained structures and energy differences are compared with available complete active space self-consistent field and multireference configuration interaction solutions. Calculations on different state crossings in LiF, formaldehyde, the ethene dimer, and the penta-2,4-dieniminium cation illustrate the discussions. Practical procedures to validate the CASPT2 solutions in polyatomic systems are explored, while it is shown that the application of the MS-CASPT2 procedure is not straightforward and requires a careful analysis of the stability of the results with the quality of the reference wave functions, that is, the size of the active space.

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.

Electrostatic Control of the Photoisomerization Efficiency and Optical Properties in Visual Pigments: On the Role of Counterion Quenching

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophores optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Interpretation of the electronic absorption spectrum of free base porphin by using multiconfigurational second-order perturbation theory

Abstract Multiconfigurational second-order perturbation (CASPT2) calculations have been performed on the low-lying optically allowed valence excited states of the free base porphin molecule in order to assign the four lowest bands of the spectrum. The low-lying triplet states have also been characterized. A basis set of the atomic natural orbital type of split-valence plus polarization quality for first-row atoms has been employed. Polarization functions are important for an accurate description of the transitions. These CASPT2 results provide a consistent picture of the experimental spectrum. Each band of the spectrum up to 4.5 eV is composed of a pair of states, which become degenerate in the spectra of metal–porphyrins. Our interpretation deviates from assignments based on other types of ab initio calculations.

The Role of Adenine Excimers in the Photophysics of Oligonucleotides

Energies and structures of different arrangements of the stacked adenine homodimer have been computed at the ab initio CASPT2 level of theory in isolation and in an aqueous environment. Adenine dimers are shown to form excimer singlet states with different degrees of stacking and interaction. A model for a 2-fold decay dynamics of adenine oligomers can be supported in which, after initial excitation in the middle UV range, unstacked or slightly stacked pairs of nucleobases will relax by an ultrafast internal conversion to the ground state, localizing the excitation in the monomer and through the corresponding conical intersection with the ground state. On the other hand, long-lifetime intrastrand stacked excimer singlet states will be formed in different conformations, including neutral and charge transfer dimers, which originate the red-shifted emission observed in the oligonucleotide chains and that will evolve toward the same monomer decay channel after surmounting an energy barrier. By computing the transient absorption spectra for the different structures considered and their relative stability in vacuo and in water, it is concluded that in the adenine homodimers the maximum-overlap face-to-face orientations are the most stable excimer conformations observed in experiment.