Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Q_x and Q_y states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Q_x-Q_y crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Role of specific solute–solvent interactions on the photophysical properties of distyryl substituted BODIPY derivatives

The role played by specific solute–solvent interactions on the spectroscopic properties of experimentally available BODIPY derivatives has been investigated.

Fatty Acid Binding to Human Serum Albumin in Blood Serum Characterized by EPR Spectroscopy

One of the functions of Human Serum Albumin (HSA) is binding and transport of fatty acids. This ability could be altered by the presence of several blood components such as toxins or peptides - which in turn alters the functionality of the protein. We aim at characterizing HSA and its fatty acid binding in native serum environment. Native ligand binding and deviations from normal function can be monitored by electron paramagnetic resonance (EPR) spectroscopy using spin labeled fatty acids (FAs). Blood serum from healthy individuals is used to examine healthy HSA in its natural physiological conditions at different loading ratios of protein to FAs. Among the EPR spectroscopic parameters (like hyperfine coupling, line shape, rotational correlation time and population of different binding sites) the rotational correlation time is found to differ significantly between binding sites of the protein, especially at loading ratios of four FAs per HSA. Although differences are observed between individual samples, a general trend regarding the dynamics of healthy HSA at different loading ratios could be obtained and compared to a reference of purified commercially available HSA in buffer.

Performance of DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals

In the present study, we tested the widely used density functionals BP86, PBE, OLYP, TPSS, M06-L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, TPSS0, M06, M06-2X, CAM-B3LYP, ωB97x, and B2PLYP with the cc-pCVQZ basis set in calculations on a set of 23 nitroxide radicals with well-resolved ¹⁴N anisotropic hyperfine coupling (HFC) constants. The results were compared with those obtained using the B3LYP/N07D and PBE/N07D methods. The convergence of the HFC values to the complete basis set limit is briefly discussed. The best results were obtained using the M06/COSMO method, with a mean absolute deviation (MAD) of 0.4 G for the dipole-dipole contribution and MAD = 0.6 G for the contact coupling contribution (as compared to 1.1 G and 1.0 G, respectively, for the B3LYP/N07D/COSMO method and 1.7 G and 0.5 G, respectively, for the B3LYP/N07D method). The majority of the functionals yielded satisfactory results for the dipole-dipole contribution, but only the M06 functional yielded similar errors for both the dipole-dipole and isotropic contributions. The RIJCOSX and RI approximations introduced errors equal to or smaller than 0.01 G.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

Density functional theory studies of MTSL nitroxide side chain conformations attached to an activation loop

A quantum mechanical (QM) method rooted on density functional theory (DFT) has been employed to determine conformations of the methane-thiosulfonate spin label (MTSL) attached to a fragment extracted from the activation loop of Aurora-A kinase. The features of the calculated energy surface revealed low energy barriers between isoenergetic minima, and the system could be described in a population of 76 rotamers that can be also considered for other systems since it was found that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{3} $$\end{document}χ3, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{4} $$\end{document}χ4 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{5} $$\end{document}χ5 do not depend on the previous two dihedral angles. Conformational states obtained were seen to be comparable to those obtained in the α-helix systems studied previously, indicating that the protein backbone does not affect the torsional profiles significantly and suggesting the possibility to use determined conformations for other protein systems for further modelling studies.

Modeling EPR parameters of nitrogen containing conjugated radical cations

DFT investigation on conjugated radical cations containing¹⁴N nucleus to obtain accurate isotropic hyperfine coupling constants.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Extension of the AMBER Force Field for Nitroxide Radicals and Combined QM/MM/PCM Approach to the Accurate Determination of EPR Parameters of DMPO-H in Solution

A computational strategy that combines both time-dependent and time-independent approaches is exploited to accurately model molecular dynamics and solvent effects on the isotropic hyperfine coupling constants of the DMPO-H nitroxide. Our recent general force field for nitroxides derived from AMBER ff99SB is further extended to systems involving hydrogen atoms in β-positions with respect to NO-moiety. The resulting force-field has been employed in a series of classical molecular dynamics simulations, comparing the computed EPR parameters from selected molecular configurations to the corresponding experimental data in different solvents. The effect of vibrational averaging on the spectroscopic parameters is also taken into account, by second-order vibrational perturbation theory involving semidiagonal third energy derivatives together first and second property derivatives.

The world as viewed by and with unpaired electrons

Recent advances in electron paramagnetic resonance (EPR) include capabilities for applications to areas as diverse as archeology, beer shelf life, biological structure, dosimetry, in vivo imaging, molecular magnets, and quantum computing. Enabling technologies include multifrequency continuous wave, pulsed, and rapid scan EPR. Interpretation is enhanced by increasingly powerful computational models.

The electronic structure of the lutein triplet state in plant light-harvesting complex II

Carotenoid molecules are essential for the life of photosynthetic organisms in that they protect the cell from the photo-oxidative damage induced by light-stress conditions. One of the photo-protective mechanisms involves triplet-triplet energy transfer from the chlorophyll molecules to the carotenoids: a process that is strongly dependent on the electronic properties of the triplet states involved. Here, we obtain a clear description of the triplet state of lutein in LHCII from higher plants for the first time by density functional theory (DFT) calculations. DFT predictions have been validated by comparison with hyperfine couplings obtained with pulsed-ENDOR spectroscopy. Knowledge of the spin density distribution, the frontier orbitals and orbital excitations forms a basis for discussing the requirements for an efficient triplet-triplet energy transfer. The results obtained for the lutein in LHCII are compared with those of the highly-substituted carotenoid peridinin in PCP from Amphidinium carterae [Di Valentin et al., Biochim. Biophys. Acta, 2008, 1777, 295-307]. The presence of substituents in the peridinin molecule does not alter significantly the triplet state electronic structure compared to lutein. Despite the unusual spectroscopic behaviour of the peridinin excited singlet state, lutein and peridinin have similar triplet state properties. In both molecules the unpaired spins are delocalized uniformly over the whole π-conjugated system in an alternating even-odd pattern.

Absorption and emission spectra of fluorescent silica nanoparticles from TD-DFT/MM/PCM calculations

A multi-scale computational protocol, which combines Quantum Mechanics and Molecular Mechanics (QM/MM) calculations with the polarisable continuum model (PCM), has been used to study the tetramethylrhodamine isothiocyanate (TRITC) fluorophore, embedded in three different environments, namely in water, on an amorphous silica surface and covalently encapsulated in a silica nanoparticle (C dot). Absorption and emission spectra have been simulated by using TD-B3LYP/PCM calculations, performed on the TRITC ground and excited state geometries, optimized at the QM/MM level. The results are in good agreement with experimental data confirming the caging effect played by the silica shell on the mobility of the TRITC molecule when covalently encapsulated in silica nanoparticles. This could result in a decrease of the nonradiative decay rate and thus an increase of the quantum yield of the molecule.

DFT Calculations of Isotropic Hyperfine Coupling Constants of Nitrogen Aromatic Radicals: The Challenge of Nitroxide Radicals

The performance of DFT methodology to predict with accuracy the isotropic hyperfine coupling constants (hfccs) of aromatic radicals containing (14)N nucleus is investigated by an extensive study in which 165 hfccs, belonging to 38 radical species, are obtained from calculations with B3LYP and PBE0 functionals combined with 6-31G*, N07D, TZVP, and EPR-III basis sets, and are compared to the reported experimental data. The results indicate that the selection of the basis set is of fundamental importance in the calculation of (14)N hfccs, whereas there is not so great an influence on the accurate computation of that parameter for (1)H nuclei. The values of the calculated (14)N coupling constants of aromatic nitroxide radicals using DFT methodology are noticeably lower than the experimental ones. A very simple relation to predict these hfccs with high accuracy is proposed on the basis of the present results, as an interesting alternative to the highly computationally demanding integrated approaches so far used.

An integrated computational protocol for the accurate prediction of EPR and PNMR parameters of aminoxyl radicals in solution

Magnetic spectroscopic techniques such as electron paramagnetic resonance (EPR) and paramagnetic NMR (PNMR) are valuable tools for understanding the structure and dynamics of complex systems such as, for example, biomolecules or nanomaterials labeled with suitable free radicals. Unfortunately, such spectra do not give direct access to the radical structure because of the subtle interplay between several different effects not easily separable and evaluable by experimentalists alone. In this respect, computational spectroscopy is becoming an essential and versatile tool for the assignment and interpretation of experimental spectra. In this article, the new integrated computational approaches developed in the recent years in our research group are reviewed. Such approaches have been applied to two widely used spin probes showing that proper account of stereo-electronic, environmental and dynamical effects leads to magnetic properties in remarkable agreement with experimental results.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.