The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations

Probing the Photochemistry of Rhodopsin Through Population Dynamics Simulations

The primary event in vision is induced by the ultrafast photoisomerization of rhodopsin, the dim-light visual pigment of vertebrates. While spectroscopic and theoretical studies have identified certain vibrationally coherent atomic motions to promote the rhodopsin photoisomerization, how exactly and to what degree such coherence is biologically related with its isomerizing efficiency (i.e. the photoisomerization quantum yield) remains unknown. In fact, in the past, the computational cost limited the simulation of the rhodopsin photoisomerization dynamics, which could be carried out only for a single molecule or a small set of molecules, therefore lacking the necessary statistical description of a molecular population motion. In this Dissertation I apply a hybrid quantum mechanics/molecular mechanics (QM/MM) models of bovine rhodopsin, the verterbrate visual pigment, to tackle the basic issues mentioned above. Accordingly, my work has been developing along three different lines comprising the development, testing and application of new tools for population dynamics simulation: (I) Development of a suitable protocol to investigate the excited state population dynamics of rhodopsins at room temperature. (II) A correlation between the phase of a hydrogen-out-of-plane (HOOP) motion at the decay point and the outcome of the rhodopsin photoisomerization. (III) A population “splitting” mechanism adopted by the protein to maximize its quantum yield and, therefore, light sensitivity. In conclusion, my Dissertation reports, for the first time, a connection between the initial coherent motion of a population of rhodopsin molecules and the quantum efficiency of their isomerization. The photoisomerization efficiency is ultimately determined by the way in which the degree of coherence of the excited state population motion is modulated by the protein sequence and conformation

Probing the Photochemistry of Rhodopsin Through Population Dynamics Simulations

The primary event in vision is induced by the ultrafast photoisomerization of rhodopsin, the dim-light visual pigment of vertebrates. While spectroscopic and theoretical studies have identified certain vibrationally coherent atomic motions to promote the rhodopsin photoisomerization, how exactly and to what degree such coherence is biologically related with its isomerizing efficiency (i.e. the photoisomerization quantum yield) remains unknown. In fact, in the past, the computational cost limited the simulation of the rhodopsin photoisomerization dynamics, which could be carried out only for a single molecule or a small set of molecules, therefore lacking the necessary statistical description of a molecular population motion. In this Dissertation I apply a hybrid quantum mechanics/molecular mechanics (QM/MM) models of bovine rhodopsin, the verterbrate visual pigment, to tackle the basic issues mentioned above. Accordingly, my work has been developing along three different lines comprising the development, testing and application of new tools for population dynamics simulation: (I) Development of a suitable protocol to investigate the excited state population dynamics of rhodopsins at room temperature. (II) A correlation between the phase of a hydrogen-out-of-plane (HOOP) motion at the decay point and the outcome of the rhodopsin photoisomerization. (III) A population “splitting” mechanism adopted by the protein to maximize its quantum yield and, therefore, light sensitivity. In conclusion, my Dissertation reports, for the first time, a connection between the initial coherent motion of a population of rhodopsin molecules and the quantum efficiency of their isomerization. The photoisomerization efficiency is ultimately determined by the way in which the degree of coherence of the excited state population motion is modulated by the protein sequence and conformation

Electronic State Mixing Controls the Photoreactivity of a Rhodopsin with all- trans Chromophore Analogues

Rhodopsins hosting synthetic retinal protonated Schiff base analogues are important for developing tools for optogenetics and high-resolution imaging. The ideal spectroscopic properties of such analogues include long-wavelength absorption/emission and fast/hindered photoisomerization. While the former may be achieved, for instance, by elongating the chromophore π-system, the latter requires a detailed understanding of the substituent effects (i.e., steric or electronic) on the chromophore light-induced dynamics. In the present letter we compare the results of quantum mechanics/molecular mechanics excited-state trajectories of native and analogue-hosting microbial rhodopsins from the eubacterium Anabaena. The results uncover a relationship between the nature of the substituent on the analogue (i.e., electron-donating (a Me group) or electron-withdrawing (a CF3 group)) and rhodopsin excited-state lifetime. Most importantly, we show that electron-donating or -withdrawing substituents cause a decrease or an increase in the electronic mixing of the first two excited states which, in turn, controls the photoisomerization speed

Comparative Quantum-Classical Dynamics of Natural and Synthetic Molecular Rotors Show How Vibrational Synchronization Modulates the Photoisomerization Quantum Efficiency

We use quantum-classical trajectories to investigate the origin of the different photoisomerization quantum efficiency observed in the dim-light visual pigment Rhodopsin and in the light-driven biomimetic molecular rotor para-methoxy N-methyl indanylidene-pyrrolinium (MeO-NAIP) in methanol. The results reveal that effective light-energy conversion requires, in general, an auxiliary molecular vibration (called promoter) that does not correspond to the rotary motion but that synchronizes with it at specific times. They also reveal that Nature has designed Rhodopsin to exploit two mechanisms working in a vibrationally coherent regime. The first uses a wag promoter to ensure that ca. 75% of the absorbed photons lead to unidirectional rotations. The second mechanism ensures that the same process is fast enough to avoid directional randomization. It is found that MeO-NAIP in methanol is incapable of exploiting the above mechanisms resulting into a 50% quantum efficiency loss. However, when the solvent is removed, MeO-NAIP rotation is predicted to synchronize with a ring-inversion promoter leading to a 30% increase of quantum efficiency and, therefore, biomimetic behavior

Localization and Frequency Identification of Large-Range Wide-Band Electromagnetic Interference Sources in Electromagnetic Imaging System

The identification and localization of large-range, wide-band electromagnetic interference (EMI) sources have always been both costly and time-consuming. The measurements at different times and places are often required before a typical system can locate a target. In this paper, we proposed a 2D electromagnetic imaging system to localize interference sources and identify the EMI frequency in real time. In this system, an offset paraboloid with a diameter of three meters is designed for large-range EMI imaging, while a multi-channel digital signal acquisition system is developed for wide-band EMI localization. The located interference source is segmented by the maximum entropy method based on particle swarm optimization, and the modified generalized regression neural network (MGRNN) is applied to identify the EMI frequency effectively by excluding misleading effects of outliers. The experiment which has been completed on our dataset indicates that our approach not only increases accuracy by 5% compared with the standard generalized regression neural network approaches for identification, but also exerts a large-range wide-band localization of the EMI source detection method

A Space-Variant Deblur Method for Focal-Plane Microwave Imaging

In the research of passive millimetre wave (PMMW) imaging, the focal plane array (FPA) can realize fast, wide-range imaging and detection. However, it has suffered from a limited aperture and off-axis aberration. Thus, the result of FPA is usually blurred by space-variant point spread function (SVPSF) and is hard to restore. In this paper, a polar-coordinate point spread function (PCPSF) model is presented to describe the circle symmetric characteristic of space-variant blur, and a log-polar-coordinate transformation (LPCT) method is propagated as the pre-processing step before the Lucy⁻Richardson algorithm to eliminate the space variance of blur. Compared with the traditional image deblur method, LPCT solves the problem by analyzing the physical model instead of the approximating it, which has proved to be a feasible way to deblur the FPA imaging system

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G∗ level of theory