Nicolas Ferré

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.

MOLCAS 8: New Capabilities for Multiconfigurational Quantum Chemical Calculations across the Periodic Table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas–Kroll–Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC‐PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large‐scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Tracking the excited-state time evolution of the visual pigment with multiconfigurational quantum chemistry

The primary event that initiates vision is the photoinduced isomerization of retinal in the visual pigment rhodopsin (Rh). Here, we use a scaled quantum mechanics/molecular mechanics potential that reproduces the isomerization path determined with multiconfigurational perturbation theory to follow the excited-state evolution of bovine Rh. The analysis of a 140-fs trajectory provides a description of the electronic and geometrical changes that prepare the system for decay to the ground state. The data uncover a complex change of the retinal backbone that, at ≈60-fs delay, initiates a space saving “asynchronous bicycle-pedal or crankshaft” motion, leading to a conical intersection on a 110-fs time scale. It is shown that the twisted structure achieved at decay features a momentum that provides a natural route toward the photoRh structure recently resolved by using femtosecond-stimulated Raman spectroscopy.

Structure, initial excited-state relaxation, and energy storage of rhodopsin resolved at the multiconfigurational perturbation theory level

We demonstrate that a “brute force” quantum chemical calculation based on an ab initio multiconfigurational second order perturbation theory approach implemented in a quantum mechanics/molecular mechanics strategy can be applied to the investigation of the excited state of the visual pigment rhodopsin (Rh) with a computational error <5 kcal·mol-1. As a consequence, the simulation of the absorption and fluorescence of Rh and its retinal chromophore in solution allows for a nearly quantitative analysis of the factors determining the properties of the protein environment. More specifically, we demonstrate that the Rh environment is more similar to the “gas phase” than to the solution environment and that the so-called “opsin shift” originates from the inability of the solvent to effectively “shield” the chromophore from its counterion. The same strategy is used to investigate three transient structures involved in the photoisomerization of Rh under the assumption that the protein cavity does not change shape during the reaction. Accordingly, the analysis of the initially relaxed excited-state structure, the conical intersection driving the excited-state decay, and the primary isolable bathorhodopsin intermediate supports a mechanism where the photoisomerization coordinate involves a “motion” reminiscent of the so-called bicycle-pedal reaction coordinate. Most importantly, it is shown that the mechanism of the ∼30 kcal·mol-1 photon energy storage observed for Rh is not consistent with a model based exclusively on the change of the electrostatic interaction of the chromophore with the protein/counterion environment.

The ultrafast photoisomerizations of rhodopsin and bathorhodopsin are modulated by bond length alternation and HOOP driven electronic effects.

Rhodopsin (Rh) and bathorhodopsin (bathoRh) quantum-mechanics/molecular-mechanics models based on ab initio multiconfigurational wave functions are employed to look at the light induced π-bond breaking and reconstitution occurring during the Rh → bathoRh and bathoRh → Rh isomerizations. More specifically, semiclassical trajectory computations are used to compare the excited (S(1)) and ground (S(0)) state dynamics characterizing the opposite steps of the Rh/bathoRh photochromic cycle during the first 200 fs following photoexcitation. We show that the information contained in these data provide an unprecedented insight into the sub-picosecond π-bond reconstitution process which is at the basis of the reactivity of the protein embedded 11-cis and all-trans retinal chromophores. More specifically, the data point to the phase and amplitude of the skeletal bond length alternation stretching mode as the key factor switching the chromophore to a bonding state. It is also confirmed/found that the phase and amplitude of the hydrogen-out-of-plane mode controls the stereochemical outcome of the forward and reverse photoisomerizations.

Specific force field parameters determination for the hybrid ab initio QM/MM LSCF method.

The pure quantum mechanics method, called Local Self‐Consistent Field (LSCF), that allows to optimize a wave function within the constraint that some predefined spinorbitals are kept frozen, is discussed. These spinorbitals can be of any shape, and their occupation numbers can be 0 or 1. Any post‐Hartree–Fock method, based on the restricted or unrestricted Hartree–Fock Slater determinant, and Kohn–Sham‐based DFT method are available. The LSCF method is easily applied to hybrid quantum mechanics/molecular mechanics (QM/MM) procedure where the quantum and the classical parts are covalently bonded. The complete methodology of our hybrid QM/MM scheme is detailed for studies of macromolecular systems. Not only the energy but also the gradients are derived; thus, the full geometry optimization of the whole system is feasible. We show that only specific force field parameters are needed for a correct description of the molecule, they are given for some general chemical bonds. A careful analysis of the errors induced by the use of molecular mechanics in hybrid computation show that a general procedure can be derived to obtain accurate results at low computation effort. The methodology is applied to the structure determination of the crambin protein and to Menshutkin reactions between primary amines and chloromethane.

The color of rhodopsins at the ab initio multiconfigurational perturbation theory resolution

We demonstrate that “brute force” quantum-mechanics/molecular-mechanics computations based on ab initio (i.e., first principles) multiconfigurational perturbation theory can reproduce the absorption maxima of a set of modified bovine rhodopsins with an accuracy allowing for the analysis of the factors determining their colors. In particular, we show that the theory accounts for the changes in excitation energy even when the proteins display the same charge distribution. Three color-tuning mechanisms, leading to changes of close magnitude, are demonstrated to operate in these conditions. The first is based on the change of the conformation of the conjugated backbone of the retinal chromophore. The second operates through the control of the distance between the positive charge residing on the chromophore and the carboxylate counterion. Finally, the third mechanism operates through the changes in orientation of the chromophore relative to the protein. These results offer perspectives for the unbiased computational design of mutants or chemically modified proteins with wanted optical properties.

Color-Tuning Mechanism of Firefly Investigated by Multi-Configurational Perturbation Method

This is the first report on a multiconfigurational reference second-order perturbation theory-molecular mechanics study of the color modulation of the observed bioluminescence of the oxyluciferin-luciferase complex of the Japanese genji-botaru firefly using structures according to recent X-ray data. Our theoretical results do not support the experimentally deduced conclusion that the color modulation of the emitted light primarily depends on the size of the compact luciferase protein cavity embedding the excited oxyluciferin molecule. Rather, we find, in agreement with recent experimental observations, that the wavelength of the emitted light depends on the polarity of the microenvironment at the phenol/phenolate terminal of the benzothiazole fragment in oxyluciferin.

The Chemistry of Bioluminescence: An Analysis of Chemical Functionalities

Firefly luciferase is one of the most studied bioluminescent systems, both theoretically and experimentally. Herein we review the current understanding of the bioluminescent process from a chemical functionality perspective based on those investigations. Three key components are emphasized: the chemiluminophore, the electron-donating fragment, and how these are affected by the substrate-enzyme interaction. The understanding is based on details of how the peroxide -O-O- bond supports the production of electronically excited products and how the charge-transfer (CT) mechanism, with the aid of an electron-donating group, lowers the activation barrier to support a reaction occurs in living organisms. For the substrate-enzyme complex it is demonstrated that the enzyme can affect the hydrogen-bonding around the CT-controlling group, resulting in a mechanism for color modulation. Finally, we analyse other luciferin-luciferase systems and compare them to the key chemical functionalities of the fragments of the luciferin-luciferase complex with respect to similarities and differences.

An artificial molecular switch that mimics the visual pigment and completes its photocycle in picoseconds

Single molecules that act as light-energy transducers (e.g., converting the energy of a photon into atomic-level mechanical motion) are examples of minimal molecular devices. Here, we focus on a molecular switch designed by merging a conformationally locked diarylidene skeleton with a retinal-like Schiff base and capable of mimicking, in solution, different aspects of the transduction of the visual pigment Rhodopsin. Complementary ab initio multiconfigurational quantum chemistry-based computations and time-resolved spectroscopy are used to follow the light-induced isomerization of the switch in methanol. The results show that, similar to rhodopsin, the isomerization occurs on a 0.3-ps time scale and is followed by <10-ps cooling and solvation. The entire (2-photon-powered) switch cycle was traced by following the evolution of its infrared spectrum. These measurements indicate that a full cycle can be completed within 20 ps.