M. Luz Sánchez

Solvent Effects on the Radiative and Nonradiative Decay of a Model of the Rhodopsin Chromophore

The radiative and nonradiative decay of a model with five double bonds of the 11-cis-retinal protonated Schiff base was studied both in vacuum and in methanol solution using an extended version of the averaged solvent electrostatic potential from molecular dynamics data (ASEP/MD) method that allows the location of crossing points between free energy surfaces both in equilibrium and in frozen solvent conditions. The multireference quantum method CASSCF was used for the description of the states of interest, while the solvent structure was obtained from molecular dynamics simulations. Electron dynamic correlation corrections to the energy were included at CASPT2 level. Unlike in gas phase, where only two states seem to be implicated, in methanol solution, three states are necessary to describe the photoisomerization process. At the Franck-Condon point the S1 and S2 states are almost degenerate; consequently, the S1 surface has a region with an ionic character ((1)Bu-like) and another one with a covalent character ((2)Ag-like). Emission from the ionic minima is responsible for the low-frequency part of the fluorescence band, while emission from the covalent minima originates the high-frequency part. The ionic minimum is separated from the conical intersection yielding the all-trans isomer by an energy barrier that was estimated in 0.7 kcal/mol. The geometry of the optimized conical intersection was found at a torsion angle of the central double bond close to 90° both in vacuum and in methanol solution. This large torsion in addition to the accompanying charge displacements forces a strong solvent reorganization during the de-excitation process which slows down the photoisomerization kinetics in methanol with respect to the gas phase. Solvent fluctuations modulate the minima depth and the barrier height and could explain the multiexponential relaxation time observed in the experiments.

Solvent effects on de-excitation channels in the p-coumaric acid methyl ester anion, an analogue of the photoactive yellow protein (PYP) chromophore

In an attempt to shed light on the environmental effects on the deactivation channels of the PYP chromophore, radiative and non-radiative deactivation mechanisms of the anionic p-coumaric acid methyl ester (pCE-) in the gas phase and water solution are compared at the CASPT2//CASSCF/cc-pVDZ level and, when necessary, at the CASPT2//CASPT2/cc-pVDZ level. We find that the solvent produces dramatic modifications on the free energy profile of the S1 state. Two twisted structures that are minima in the gas phase could not be localized in aqueous solution. Furthermore, the relative stability of minima and conical intersections (CIs) is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond and the non-radiative, non-reactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with lower energy, while in solution the CIβ structure, characterized by a large charge separation, is strongly stabilized by interactions with water molecules and becomes the structure with the lowest energy on S1. These facts explain the low fluorescence signal of pCE- in aqueous solution and the presence of partial trans-cis photoisomerization in this system.

How Methylation Modifies the Photophysics of the Native All-Trans Retinal Protonated Schiff Base: A CASPT2/MD Study in Gas Phase and in Methanol

A comparison between the free-energy surfaces of the all- trans-retinal protonated Schiff base (RPSB) and its 10-methylated derivative in gas phase and methanol solution is performed at CASSCF//CASSCF and CASPT2//CASSCF levels. Solvent effects were included using the average solvent electrostatic potential from molecular dynamics method. This is a QM/MM (quantum mechanics/molecular mechanics) method that makes use of the mean field approximation. It is found that the methyl group bonded to C10 produces noticeable changes in the solution free-energy profile of the S1 excited state, mainly in the relative stability of the minimum energy conical intersections (MECIs) with respect to the Franck-Condon (FC) point. The conical intersections yielding the 9- cis and 11- cis isomers are stabilized while that yielding the 13- cis isomer is destabilized; in fact, it becomes inaccessible by excitation to S1. Furthermore, the planar S1 minimum is not present in the methylated compound. The solvent notably stabilizes the S2 excited state at the FC geometry. Therefore, if the S2 state has an effect on the photoisomerization dynamics, it must be because it permits the RPSB population to branch around the FC point. All these changes combine to speed up the photoisomerization in the 10-methylated compound with respect to the native compound.

QM/MM Study of Substituent and Solvent Effects on the Excited State Dynamics of the Photoactive Yellow Protein Chromophore

Substituent and solvent effects on the excited state dynamics of the Photoactive Yellow Protein chromophore are studied using the average solvent electrostatic potential from molecular dynamics (ASEP/MD) method. Four molecular models were considered: the ester and thioester derivatives of the p-coumaric acid anion and their methylated derivatives. We found that the solvent produces dramatic modifications on the free energy profile of the S1 state: 1) Two twisted structures that are minima in the gas phase could not be located in aqueous solution. 2) Conical intersections (CIs) associated with the rotation of the single bond adjacent to the phenyl group are found for the four derivatives in water solution but only for thio derivatives in the gas phase. 3) The relative stability of minima and CIs is reverted with respect to the gas phase values, affecting the prevalent de-excitation paths. As a consequence of these changes, three competitive de-excitation channels are open in aqueous solution: the fluorescence emission from a planar minimum on S1, the trans-cis photoisomerization through a CI that involves the rotation of the vinyl double bond, and the nonradiative, nonreactive, de-excitation through the CI associated with the rotation of the single bond adjacent to the phenyl group. In the gas phase, the minima are the structures with the lower energy, while in solution these are the conical intersections. In solution, the de-excitation prevalent path seems to be the photoisomerization for oxo compounds, while thio compounds return to the initial trans ground state without emission.

Accelerating QM/MM Calculations by Using the Mean Field Approximation

It is well known that solvents can modify the frequency and intensity of the solute spectral bands, the thermodynamics and kinetics of chemical reactions, the strength of molecular interactions or the fate of solute excited states. The theoretical study of solvent effects is quite complicated since the presence of the solvent introduces additional difficulties with respect to the study of analogous problems in gas phase. The mean field approximation (MFA) is used for many of the most employed solvent effect theories as it permits to reduce the computational cost associated to the study of processes in solution. In this chapter we revise the performance of ASEP/MD, a quantum mechanics/molecular mechanics method developed in our laboratory that makes use of this approximation. It permits to combine state of the art calculations of the solute electron distribution with a detailed, microscopic, description of the solvent. As examples of application of the method we study solvent effects on the absorption spectra of some molecules involved in photoisomerization processes of biological systems.

A new QM/MM method oriented to the study of ionic liquids

The interest on room temperature ionic liquids has grown in the last decades because of their use as all‐purpose solvent and their low environmental impact. In the present work, a new theoretical procedure is developed to study pure ionic liquids within the framework of the quantum mechanics/molecular mechanics method. Each type of ion (cation or anion) is considered as an independent entity quantum mechanically described that follows a differentiated path in the liquid. The method permits, through an iterative procedure, the full coupling between the polarized charge distribution of the ions and the liquid structure around them. The procedure has been tested with 1‐ethyl‐3‐methylimidazolium tetrafluoroborate. It was found that, similar to non‐polar liquids and as a consequence of the low value of the reaction field, the cation and anion charge distributions are hardly polarized by the rest of molecules in the liquid. Their structure is characterized by an alternance between anion and cation shells as evidenced by the coincidence of the first maximum of the anion–anion and cation–cation radial distribution functions with the first minimum of the anion‐cation. Some degree of stacking between the cations is also found.