A Theoretical Study of the Low-Lying States of the Anionic and Protonated Ionic Forms of Urocanic Acid

Active Learning Assisted MCCI to Target Spin States

Strongly correlated systems and their accurate solutions have been challenging to quantum chemistry. Several methods have been developed over the years for the accurate understanding of such systems, and selected configuration interaction and Monte Carlo configuration interaction (MCCI) form important classes of systems in this category. However, MCCI is plagued by slow convergence. This is further exacerbated by the fact that most of the current MCCI implementations do not target specific spin states. In our work, we use active learning assisted MCCI to speed up the convergence manyfold and also develop a method for spin targeting. This method has been tested with several model Hamiltonian systems akin to molecular systems and has shown improved convergence and accuracy.

Photoprotection Mechanism of p-Methoxy Methylcinnamate: A CASPT2 Study

p-Methoxy methylcinnamate (p-MMC) shares the same molecular skeleton with octyl methoxycinnamate sunscreen. It is recently found that adding one water to p-MMC can significantly enhance the photoprotection efficiency. However, the physical origin is elusive. Herein we have employed multireference complete active space self-consistent field (CASSCF) and multistate complete active-space second-order perturbation (MS-CASPT2) methods to scrutinize the photophysical and photochemical mechanism of p-MMC and its one-water complex p-MMC-W. Specifically, we optimize the stationary-point structures on the (1)ππ*, (1)nπ*, and S0 potential energy surfaces to locate the (1)ππ*/S0 and (1)ππ*/(1)nπ* conical intersections and to map (1)ππ* and (1)nπ* excited-state relaxation paths. On the basis of the results, we find that, for the trans p-MMC, the major (1)ππ* deactivation path is decaying to the dark (1)nπ* state via the in-plane (1)ππ*/(1)nπ* crossing point, which only need overcome a small barrier of 2.5 kcal/mol; the minor one is decaying to the S0 state via the (1)ππ*/S0 conical intersection induced by out-of-plane photoisomerization. For the cis p-MMC, these two decay paths are comparable (1)ππ* deactivation paths: one is decaying to the dark (1)nπ* state via the (1)ππ*/(1)nπ* crossing point, and the second is decaying to the ground state via the (1)ππ*/S0 conical intersection. One-water hydration stabilizes the (1)ππ* state and meanwhile destabilizes the (1)nπ* state. As a consequence, the (1)ππ* deactivation path to the dark (1)nπ* state is heavily inhibited. The related barriers are increased to 5.8 and 3.3 kcal/mol for the trans and cis p-MMC-W, respectively. In comparison, the barriers associated with the photoisomerization-induced (1)ππ* decay paths are reduced to 2.5 and 1.3 kcal/mol for the trans and cis p-MMC-W. Therefore, the (1)ππ* decay paths to the S0 state are dominant relaxation channels when adding one water molecule. Finally, the present work contributes a lot of knowledge to understanding the photoprotection mechanism of methylcinnamate derivatives.

Structure, Spectroscopy, and Spectral Tuning of the Gas-Phase Retinal Chromophore: The β-Ionone “Handle” and Alkyl Group Effect

The low-lying singlet states (i.e. S0, S1, and S2) of the chromophore of rhodopsin, the protonated Schiff base of 11-cis-retinal (PSB11), and of its all-trans photoproduct have been studied in isolated conditions by using ab initio multiconfigurational second-order perturbation theory. The computed spectroscopic features include the vertical excitation, the band origin, and the fluorescence maximum of both isomers. On the basis of the S0-->S1 vertical excitation, the gas-phase absorption maximum of PSB11 is predicted to be 545 nm (2.28 eV). Thus, the predicted absorption maximum appears to be closer to that of the rhodopsin pigment (2.48 eV) and considerably red-shifted with respect to that measured in solution (2.82 eV in methanol). In addition, the absorption maxima associated with the blue, green, and red cone visual pigments are tentatively rationalized in terms of the spectral changes computed for PSB11 structures featuring differently twisted beta-ionone rings. More specifically, a blue-shifted absorption maximum is explained in terms of a large twisting of the beta-ionone ring (with respect to the main conjugated chain) in the visual S-cone (blue) pigment chromophore. In contrast, the chromophore of the visual L-cone (red) pigment is expected to have a nearly coplanar beta-ionone ring yielding a six double bond fully conjugated framework. Finally, the M-cone (green) chromophore is expected to feature a twisting angle between 10 and 60 degrees. The spectroscopic effects of the alkyl substituents on the PSB11 spectroscopic properties have also been investigated. It is found that they have a not negligible stabilizing effect on the S1-S0 energy gap (and, thus, cause a red shift of the absorption maximum) only when the double bond of the beta-ionone ring conjugates significantly with the rest of the conjugated chain.

Ground and excited state CASPT2 geometry optimizations of small organic molecules

A method for computing second-order multiconfigurational perturbation theory (CASPT2) energy gradients numerically has been implemented and applied to a range of elementary organic chromophores, including 1,3 butadiene, acrolein, and two protonated Schiff bases. Geometries of ground and excited states-as well as conical intersections-are compared with the corresponding CASSCF structures, illustrating the effect of including the correction for dynamical electron correlation. It is shown that the differences between the two methods are not readily categorized, but that, while individual changes in bond lengths can be quite large ( approximately 0.01-0.02 A), the natures and CASPT2 energetics of the structures remain similar. Exceptions to this tend to be systems that have a strong ionic character and that are not well described at the CASSCF level. Basis set effects (double- vs. triple-zeta) were examined for a limited number of examples, and found to be quite dramatic at both levels of theory.

Reactive Oxygen Species Formation by UV-A Irradiation of Urocanic Acid and the Role of Trace Metals in This Chemistry¶

We have extended our study of the decomposition of urocanic acid (UCA) with ultraviolet A radiation (UV-A) by the self-sensitized generation of singlet oxygen (see Photochem. Photobiol. 75, 565 [2002]). The chemistry has been found to be partially dependent on the presence of trace metal, most likely iron. Rigorous removal of metal impurities from the reaction mixture, using Chelex, retarded (but did not eliminate) the UV-A-initiated UCA degradation. The addition of small amounts of ferric chloride to the Chelex-treated solutions restored reactivity. Chelex treatment had a modest effect on the previously reported ability of UCA photoproducts to photonick supercoiled plasmid DNA. Also, photoinactivation of Sindbis virus on irradiation with the UCA photoproducts is now reported. Inactivation of the virus by a photoproduct mixture derived from a UCA solution that had been pretreated with Chelex was less rapid and gave better behaved time-course plots than was observed for photoproducts from non-Chelex treated solutions. These results are particularly noteworthy in light of the ubiquitous presence of both UCA and iron in the skin.

A TD-DFT study of the photochemistry of urocanic acid in biologically relevant ionic, rotameric, and protomeric forms

The photochemistry of Urocanic acid, a chromophore present in human skin and linked to photoimmunosuppression and skin cancer, is investigated theoretically by means of time-dependent density functional theory. Extensive calculations are carried out for different ionic, rotameric, and protomeric forms of both the trans and cis form. Inclusion of solvation effects, here accounted for by means of a continuum solvent model, are found to be crucial for the correct description of the biologically relevant zwitterionic forms of the molecule. For the trans zwitterionic form, it is found that the planar form usually assumed in the literature is not stable, and that a realistic charge separation cannot be achieved in the gas phase. Calculated vertical excitation energies are in excellent agreement with available experimental data, with a weakly absorbing n --> pi transition around 4.0 eV, and strongly absorbing pi --> pi transitions at 4.5-4.9 eV. The debated intramolecular hydrogen bond is predicted to have a modest impact on the vertical spectra in solution, but improves agreement with experiment when included. In general, we also predict that different rotameric forms have very similar absorption spectra. In addition, we find a candidate absorbing state to link trans-urocanic acid to singlet oxygen production and subsequent photoaging of the skin.