Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

On the Statistical Regime, Coherence versus Incoherence and Ergodicity of Quantum Vibrational Trajectories in Soft Condensed Molecular Systems

A theoretical-computational procedure, recently proposed for modelling Vibrational Energy Relaxation (VER) processes of a molecule (Quantum Center, QC) embedded in a complex atomic-molecular system, is extended and applied for analyzing in detail the features of the QC density matrix (DM) temporal evolution. The results, obtained using aqueous azide ion as a case study, show the total lack of coherence in the DM, when the system is prepared to be initially in a pure vibrational eigenstate. This finding is fully in line with the statistical interpretation of the process typically adopted also in the experimental studies where the relaxation processes are all described within the typical schemes of chemical kinetics. Consistently, when the initial vibrational state corresponds to an eigenstate mixture, although initially coherent, the DM relaxes to a fully incoherent condition with a mean lifetime related to the one of the diagonal elements relaxation. These specific DM features turn out to be essentially governed by the thermal equilibrium condition of the atomic-molecular classical coordinates which drive the ensemble of the quantum-trajectories toward the observed statistical regime. Finally, from the analysis of a single long timescale quantum vibrational trajectory it also clearly emerges its ergodic behaviour.

Theoretical-Computational Modeling of CD Spectra of Aqueous Monosaccharides by Means of Molecular Dynamics Simulations and Perturbed Matrix Method

The electronic circular dichroism (ECD) spectra of aqueous d-glucose and d-galactose were modeled using a theoretical-computational approach combining molecular dynamics (MD) simulations and perturbed matrix method (PMM) calculations, hereafter termed MD-PMM. The experimental spectra were reproduced with a satisfactory accuracy, confirming the good performances of MD-PMM in modeling different spectral features in complex atomic-molecular systems, as already reported in previous studies. The underlying strategy of the method was to perform a preliminary long timescale MD simulation of the chromophore followed by the extraction of the relevant conformations through essential dynamics analysis. On this (limited) number of relevant conformations, the ECD spectrum was calculated via the PMM approach. This study showed that MD-PMM was able to reproduce the essential features of the ECD spectrum (i.e., the position, the intensity, and the shape of the bands) of d-glucose and d-galactose while avoiding the rather computationally expensive aspects, which were demonstrated to be important for the final outcome, such as (i) the use of a large number of chromophore conformations; (ii) the inclusion of quantum vibronic coupling; and (iii) the inclusion of explicit solvent molecules interacting with the chromophore atoms within the chromophore itself (e.g., via hydrogen bonds).