Quantum monte carlo study of the retinal minimal model C5H6NH2+

Toward quantum Monte Carlo forces on heavier ions: Scaling properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Z_eff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z_eff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Z_eff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z_eff ^6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z_eff ^-2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z_eff.

Variational Multistate Density Functional Theory for a Balanced Treatment of Static and Dynamic Correlations

Here, an approach to variational multistate density functional theory (vMSDFT) is explored. In this approach, the Kohn-Sham orbitals as well as configuration coefficients were simultaneously optimized, thus yielding a full variational minimum. Furthermore, this work also proposes two important improvements on the MSDFT framework. First, a "point-to-point correction" is used to correct the static correlation present in the DFT framework. Therefore, double counting of static correlation in vMSDFT is mitigated. Second, a general form to construct the transition density functional in the vMSDFT framework is proposed, which allows for the properties of vMSDFT wave functions to be standardized to the complete active space self-consistent field properties. The utility of vMSDFT is illustrated on molecular systems of interest including bond breaking, diradicals, excited states, and conical intersections. The numerical results suggest that the accuracy of vMSDFT is in close agreement with the high-level multireference methods.

Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2

This study provides gas-phase S₁ excited-state geometries along with emission and adiabatic energies for methylated/demethylated and ring-locked analogues of protonated Schiff base retinal models comprising system of five conjugated double bonds (PSB5), using second order multiconfiguration perturbation theory (CASPT2). CASPT2 results serve as reference data to assess the performance of CC2 (second-order approximate coupled cluster singles and doubles) and a commonly used CASSCF/CASPT2 protocol, that is, complete active space self-consistent field (CASSCF) geometry optimization followed by CASPT2 energy calculation. We find that the CASSCF methodology fails to locate planar S₁ minimum energy structures for four out of five investigated planar models in contrast to CC2 and CASPT2 methods. However, for those which were found: one planar and two twisted minima, there is an excellent agreement between CASSCF and CASPT2 results in terms of geometrical parameters, one-electron properties, as well as emission and adiabatic energies. CC2 performs well for in-plane S₁ minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S₂ electronic state, with strong multireference and significant double excitation character, is very poor, exhibiting errors in transition energies exceeding 1 eV. They may be substantially diminished by recalculating transition energies with CASPT2 method. Our work shows that CASSCF/CASPT2 and CC2 shortcomings may influence gas-phase retinal analogues' excited state description in a dramatic way. © 2017 Wiley Periodicals, Inc.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model

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.

Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics

Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.

Ab Initio Geometry and Bright Excitation of Carotenoids: Quantum Monte Carlo and Many Body Green's Function Theory Calculations on Peridinin

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 Bu(+)-like (S2) bright excited state. Many Body Green's Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu(+)-like state for the VMC structure (VMC/MBGFT) provide an excitation energy of 2.62 eV, in agreement with experimental results in n-hexane (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.

Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule

Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.

Protein Field Effect on the Dark State of 11-cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics

The accurate determination of the geometrical details of the dark state of 11-cis 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. [J. Am. Chem. Soc. 2004, 126, 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 S1 and S2 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.

Reaction pathways by quantum Monte Carlo: insight on the torsion barrier of 1,3-butadiene, and the conrotatory ring opening of cyclobutene

Quantum Monte Carlo (QMC) methods are used to investigate the intramolecular reaction pathways of 1,3-butadiene. The ground state geometries of the three conformers s-trans, s-cis, and gauche, as well as the cyclobutene structure are fully optimised at the variational Monte Carlo (VMC) level, obtaining an excellent agreement with the experimental results and other quantum chemistry high level calculations. Transition state geometries are also estimated at the VMC level for the s-trans to gauche torsion barrier of 1,3-butadiene and for the conrotatory ring opening of cyclobutene to the gauche-1,3-butadiene conformer. The energies of the conformers and the reaction barriers are calculated at both variational and diffusional Monte Carlo levels providing a precise picture of the potential energy surface of 1,3-butadiene and supporting one of the two model profiles recently obtained by Raman spectroscopy [Boopalachandran et al., J. Phys. Chem. A 115, 8920 (2011)]. Considering the good scaling of QMC techniques with the system's size, our results also demonstrate how variational Monte Carlo calculations can be applied in the future to properly investigate the reaction pathways of large and correlated molecular systems.

Optimized Structure and Vibrational Properties by Error Affected Potential Energy Surfaces

The precise theoretical determination of the geometrical parameters of molecules at the minima of their potential energy surface and of the corresponding vibrational properties are of fundamental importance for the interpretation of vibrational spectroscopy experiments. Quantum Monte Carlo techniques are correlated electronic structure methods promising for large molecules, which are intrinsically affected by stochastic errors on both energy and force calculations, making the mentioned calculations more challenging with respect to other more traditional quantum chemistry tools. To circumvent this drawback in the present work, we formulate the general problem of evaluating the molecular equilibrium structures, the harmonic frequencies, and the anharmonic coefficients of an error affected potential energy surface. The proposed approach, based on a multidimensional fitting procedure, is illustrated together with a critical evaluation of systematic and statistical errors. We observe that the use of forces instead of energies in the fitting procedure reduces the statistical uncertainty of the vibrational parameters by 1 order of magnitude. Preliminary results based on variational Monte Carlo calculations on the water molecule demonstrate the possibility to evaluate geometrical parameters and harmonic and anharmonic coefficients at this level of theory with an affordable computational cost and a small stochastic uncertainty (<0.07% for geometries and <0.7% for vibrational properties).