The formation of hybrid light–molecule states (polaritons) offers a new strategy to manipulate the photochemistry of molecules. To fully exploit its potential, one needs to build a toolbox of polaritonic phenomenologies that supplement those of standard photochemistry. By means of a state-of-the-art computational photochemistry approach extended to the strong-coupling regime, here we disclose various mechanisms peculiar of polaritonic chemistry: coherent population oscillations between polaritons, quenching by trapping in dead-end polaritonic states and the alteration of the photochemical reaction pathway and quantum yields. We focus on azobenzene photoisomerization, that encompasses the essential features of complex photochemical reactions such as the presence of conical intersections and reaction coordinates involving multiple internal modes. In the strong coupling regime, a polaritonic conical intersection arises and we characterize its role in the photochemical process. Our chemically detailed simulations provide a framework to rationalize how the strong coupling impacts the photochemistry of realistic molecules.

Diradical molecules are essential species involved in many organic and inorganic chemical reactions. The computational study of their electronic structure is often challenging, because a reliable description of the correlation, and in particular of the static one, requires multireference techniques. The Jastrow correlated antisymmetrized geminal power (JAGP) is a compact and efficient wave function ansatz, based on the valence-bond representation, which can be used within quantum Monte Carlo (QMC) approaches. The AGP part can be rewritten in terms of molecular orbitals, obtaining a multideterminant expansion with zero-seniority number. In the present work we demonstrate the capability of the JAGP ansatz to correctly describe the electronic structure of two diradical prototypes: the orthogonally twisted ethylene, C2H4, and the methylene, CH2, representing respectively a homosymmetric and heterosymmetric system. In the orthogonally twisted ethylene, we find a degeneracy of π and π* molecular orbitals, as correctly predicted by multireference procedures, and our best estimates of the twisting barrier, using respectively the variational Monte Carlo (VMC) and the lattice regularized diffusion Monte Carlo (LRDMC) methods, are 71.9(1) and 70.2(2) kcal/mol, in very good agreement with the high-level MR-CISD+Q value, 69.2 kcal/mol. In the methylene we estimate an adiabatic triplet-singlet (X̃(3)B1-ã(1)A1) energy gap of 8.32(7) and 8.64(6) kcal/mol, using respectively VMC and LRDMC, consistently with the experimental-derived finding for Te, 9.363 kcal/mol. On the other hand, we show that the simple ansatz of a Jastrow correlated single determinant (JSD) wave function is unable to provide an accurate description of the electronic structure in these diradical molecules, both at variational level (VMC torsional barrier of C2H4 of 99.3(2) kcal/mol, triplet-singlet energy gap of CH2 of 13.45(10) kcal/mol) and, more remarkably, in the fixed-nodes projection schemes (LRDMC torsional barrier of 97.5(2) kcal/mol, triplet-singlet energy gap of 13.36(8) kcal/mol) showing that a poor description of the static correlation yields an inaccurate nodal surface. The suitability of JAGP to correctly describe diradicals with a computational cost comparable with that of a JSD calculation, in combination with a favorable scalability of QMC algorithms with the system size, opens new perspectives in the ab initio study of large diradical systems, like the transition states in cycloaddition reactions and the thermal isomerization of biological chromophores.

We explore the computation of high-harmonic generation spectra by means of Gaussian basis sets in approaches propagating the time-dependent Schrödinger equation. We investigate the efficiency of Gaussian functions specifically designed for the description of the continuum proposed by Kaufmann et al. [J. Phys. B 22, 2223(1989]. We assess the range of applicability of this approach by studying the hydrogen atom, i.e. the simplest atom for which "exact" calculations on a grid can be performed. We notably study the effect of increasing the basis set cardinal number, the number of diffuse basis functions, and the number of Gaussian pseudo-continuum basis functions for various laser parameters. Our results show that the latter significantly improve the description of the low-lying continuum states, and provide a satisfactory agreement with grid calculations for laser wavelengths λ0 = 800 and 1064 nm. The Kaufmann continuum functions therefore appear as a promising way of constructing Gaussian basis sets for studying molecular electron dynamics in strong laser fields using time-dependent quantum-chemistry approaches.

The penta-2,4-dieniminium cation (PSB3) displays similar ground state and first excited state potential energy features as those of the retinal protonated Schiff base (RPSB) chromophore in rhodopsin. Recently, PSB3 has been used to benchmark several electronic structure methods, including highly correlated multireference wave function approaches, highlighting the necessity to accurately describe the electronic correlation in order to obtain reliable properties even along the ground state (thermal) isomerization paths. In this work, we apply two quantum Monte Carlo approaches, the variational Monte Carlo and the lattice regularized diffusion Monte Carlo, to study the energetics and electronic properties of PSB3 along representative minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo is used in combination with the Jastrow antisymmetrized geminal power ansatz, which guarantees an accurate and balanced description of the static electronic correlation thanks to the multiconfigurational nature of the antisymmetrized geminal power term, and of the dynamical correlation, due to the presence of the Jastrow factor explicitly depending on electron–electron distances. Along the two ground state isomerization minimum energy paths of PSB3, CASSCF calculations yield wave functions having either charge transfer or diradical character in proximity of the two transition state configurations. Here, we observe that at the quantum Monte Carlo level of theory, only the transition state with charge transfer character can be located. The conical intersection, which becomes highly sloped, is observed only if the path connecting the two original CASSCF transition states is extended beyond the diradical one, namely by increasing the bond-length-alternation (BLA). These findings are in good agreement with the results obtained by MRCISD+Q calculations, and they demonstrate the importance of having an accurate description of the static and dynamical correlation when studying isomerization and transition states of conjugated systems.

The accurate determination of the geometrical details of the dark state of 11- retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the -ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [. ,, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S and S states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.

In this letter, we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP, and CAM-B3LYP) and shorter than that calculated at the Hartree-Fock (HF) level. TDDFT and EOM-CCSD calculations on a reduced peridinin model confirm the HOMO-LUMO major contribution of the B u + -like (S 2 ) bright excited state. Many Body Green's Function Theory (MBGFT) calculations of the vertical excitation energy of the B u + -like state for the VMC structure (VMC/ MBGFT) provide an excitation energy of 2.62 eV, in agreement with experimental results in nhexane (2.72 eV). The dependence of the excitation energy on the bond length alternation in the MBGFT and TDDFT calculations with different functionals is discussed.Carotenoids are among the most abundant chemical species in biological systems. Similarly to chlorophylls, they play a fundamental role in light harvesting and energy transfer mechanism in photosynthetic organisms.1 In addition, carotenoids get the necessary photoprotective activity of quenching triplet chlorophylls and singlet oxygen molecules.2 Although carotenoids are characterized by a large variety in terms of conjugation length of the polyenic chain and of substituents, they share some properties like the state ordering and the character of the singlet and triplet excitations. In the present letter, we consider the carotenoid peridinin (PID, Chart 1), present in the Peridinin-Chlorophyll-a Protein (PCP), a water-soluble complex deriving from marine dinoflagellate Amphidinium carterae containing the highest peridinin to chlorophyll-a ratio in nature, namely 4:1 for each domain.3 After a single-photon singlet excitation S 0 → S 2 (B u + -like, one-photon allowed), the chromophore couples with an adjacent chlorophyll-a according to a resonant energy transfer mechanism. Two energy transfer channels have been experimentally detected:2-7 the first, * leonardo.guidoni@univaq.it. Author ContributionsThe manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. NotesThe authors declare no competing financial interest. The accurate determination of the relaxed ground state geometry of carotenoids still represents a challenge for quantum chemistry calculations, due to the difficulties in the correct description of electronic correlation of the conjugated polyenic chain. Nevertheless, structural effects play an important role on the spectral tuning of carotenoids in the gas phase, in solution, or in their protein environment, and a reference high-level structure would be desirable. A key structural parameter in the spectral tuning of linear chromophores is the average bond length alternation (BLA), defined as the difference between the average of single bond and double bond distances, excluding (in the case of PID) the terminal double bond ...

Variational MonteCarlo and Diffusion MonteCarlo calculations have been carried out for cations like Li + , Na + and K + as dopants of small helium clusters over a range of cluster sizes up to about 12 solvent atoms.The interaction has been modelled through a sum-of-potential picture that disregards higher order effects beyond atom-atom and atom-ion contributions. The latter were obtained from highly correlated ab-initio calculations over a broad range of interatomic distances.This study focuses on two of the most striking features of the microsolvation in a quantum solvent of a cationic dopant: electrostriction and snowball effects. They are here discussed in detail and in relation with the nanoscopic properties of the interaction forces at play within a fully quantum picture of the clusters features.

In this work, we study the electronic and geometrical properties of the ground state of the Retinal Minimal Model C(5)H(6)NH(2)(+) using the variational Monte Carlo (VMC) method by means of the Jastrow antisymmetrized geminal power (JAGP) wavefunction. A full optimization of all wavefunction parameters, including coefficients, and exponents of the atomic basis, has been achieved, giving converged geometries with a compact and correlated wavefunction. The relaxed geometries of the cis and trans isomers present a pronounced bond length alternation pattern characterized by a C=C central double bond slightly shorter than that reported by the CASPT2 structures. The comparison between different basis sets indicates converged values of geometrical parameters, energy differences, and dipole moments even when the smallest wavefunction is used. The compactness of the wavefunction as well as the scalability of VMC optimization algorithms on massively parallel computers opens the way to perform full structural optimizations of conjugated biomolecules of hundreds of electrons by correlated methods like Quantum Monte Carlo.