Validation of the B3LYP/N07D and PBE0/N07D Computational Models for the Calculation of Electronic g-Tensors

Benchmarking Basis Sets for Density Functional Theory Thermochemistry Calculations: Why Unpolarized Basis Sets and the Polarized 6-311G Family Should Be Avoided

Basis sets are a crucial but often largely overlooked choice in setting up quantum chemistry calculations. The choice of the basis set can be critical in determining the accuracy and calculation time of your quantum chemistry calculations. Clear recommendations based on thorough benchmarking are essential but not readily available currently. This study investigates the relative quality of basis sets for general properties by benchmarking basis set performance for a diverse set of 139 reactions (from the diet-150-GMTKN55 data set). In our analysis, we find the distributions of errors are often significantly non-Gaussian, meaning that the joint consideration of median errors, mean absolute errors, and outlier statistics is helpful to provide a holistic understanding of basis set performance. Our direct comparison of performance between most modern basis sets provides quantitative evidence for basis set recommendations that broadly align with the established understanding of basis set experts and is evident in the design of modern basis sets. For example, while zeta is a good measure of quality, it is not the only determining factor for an accurate calculation with unpolarized double- and triple-ζ basis sets (like 6-31G and 6-311G) having very poor performance. Appropriate use of polarization functions (e.g., 6-31G*) is essential to obtain the accuracy offered by double- or triple-ζ basis sets. In our study, the best performances for double- and triple-ζ basis sets are 6-31++G** and pcseg-2, respectively. However, the performances of singly polarized double-ζ and doubly polarized triple-ζ basis sets are quite similar with one key exception: the polarized 6-311G basis set family has poor parametrization, which means its performance is more like a double-ζ than a triple-ζ basis set. All versions of the 6-311G basis set family should be avoided entirely for valence chemistry calculations moving forward.

Probing Liquid-Ordered and Disordered Phases in Lipid Model Membranes: A Combined Theoretical and Spectroscopic Study of a Fluorescent Molecular Rotor

An integrated theoretical/experimental strategy has been applied to the study of environmental effects on the spectroscopic parameters of 4-(diphenylamino)phtalonitrile (DPAP), a fluorescent molecular rotor. The computational part starts from the development of an effective force field for the first excited electronic state of DPAP and proceeds through molecular dynamics simulations in solvents of different polarities toward the evaluation of Stokes shifts by quantum mechanics/molecular mechanics (QM/MM) approaches. The trends of the computed results closely parallel the available experimental results thus giving confidence to the interpretation of new experimental studies of the photophysics of DPAP in lipid bilayers. In this context, results show unambiguously that both flexible dihedral angles and global rotations are significantly retarded in a cholesterol/DPPC lipid matrix with respect to the DOPC matrix, thus confirming the sensitivity of DPAP to probe different environments and, therefore, its applicability as a probe for detecting different structures and levels of plasma membrane organization.

A radical approach to radicals

Nitroxide radicals are characterized by a long-lived spin-unpaired electronic ground state and are strongly sensitive to their chemical surroundings. Combined with electron paramagnetic resonance spectroscopy, these electronic features have led to the widespread application of nitroxide derivatives as spin labels for use in studying protein structure and dynamics. Site-directed spin labelling requires the incorporation of nitroxides into the protein structure, leading to a new protein-ligand molecular model. However, in protein crystallographic refinement nitroxides are highly unusual molecules with an atypical chemical composition. Because macromolecular crystallography is almost entirely agnostic to chemical radicals, their structural information is generally less accurate or even erroneous. In this work, proteins that contain an example of a radical compound (Chemical Component Dictionary ID MTN) from the nitroxide family were re-refined by defining its ideal structural parameters based on quantum-chemical calculations. The refinement results show that this procedure improves the MTN ligand geometries, while at the same time retaining higher agreement with experimental data.

Controlling pairing of π-conjugated electrons in 2D covalent organic radical frameworks via in-plane strain

Controlling the electronic states of molecules is a fundamental challenge for future sub-nanoscale device technologies. π-conjugated bi-radicals are very attractive systems in this respect as they possess two energetically close, but optically and magnetically distinct, electronic states: the open-shell antiferromagnetic/paramagnetic and the closed-shell quinoidal diamagnetic states. While it has been shown that it is possible to statically induce one electronic ground state or the other by chemical design, the external dynamical control of these states in a rapid and reproducible manner still awaits experimental realization. Here, via quantum chemical calculations, we demonstrate that in-plane uniaxial strain of 2D covalently linked arrays of radical units leads to smooth and reversible conformational changes at the molecular scale that, in turn, induce robust transitions between the two kinds of electronic distributions. Our results pave a general route towards the external control, and thus technological exploitation, of molecular-scale electronic states in organic 2D materials.

Controlling the electronic states of molecules is a fundamental challenge for future sub-nanoscale device technologies but the external dynamical control of these states still awaits experimental realization. Here, via quantum chemical calculations, the authors demonstrate that in-plane uniaxial strain of 2D covalently linked arrays of radical units induces controlled pairing of π-conjugated electrons in a reversible way.

High-resolution EPR distance measurements on RNA and DNA with the non-covalent Ǵ spin label

Pulsed electron paramagnetic resonance (EPR) experiments, among them most prominently pulsed electron-electron double resonance experiments (PELDOR/DEER), resolve the conformational dynamics of nucleic acids with high resolution. The wide application of these powerful experiments is limited by the synthetic complexity of some of the best-performing spin labels. The recently developed $\bf\acute{G}$ (G-spin) label, an isoindoline-nitroxide derivative of guanine, can be incorporated non-covalently into DNA and RNA duplexes via Watson-Crick base pairing in an abasic site. We used PELDOR and molecular dynamics (MD) simulations to characterize $\bf\acute{G}$, obtaining excellent agreement between experiments and time traces calculated from MD simulations of RNA and DNA double helices with explicitly modeled $\bf\acute{G}$ bound in two abasic sites. The MD simulations reveal stable hydrogen bonds between the spin labels and the paired cytosines. The abasic sites do not significantly perturb the helical structure. $\bf\acute{G}$ remains rigidly bound to helical RNA and DNA. The distance distributions between the two bound $\bf\acute{G}$ labels are not substantially broadened by spin-label motions in the abasic site and agree well between experiment and MD. $\bf\acute{G}$ and similar non-covalently attached spin labels promise high-quality distance and orientation information, also of complexes of nucleic acids and proteins.

Electronic g-tensors of nanodiamonds: Dependence on the size, shape, and surface functionalization

The electronic g-tensor dependence on the size, shape, and surface functionalization of nanodiamonds (NDs) is theoretically investigated by selecting dangling bonds and single substitutional nitrogen atoms as a main source of the unpaired electrons. The performed g-tensor calculations reveal that aforementioned paramagnetic impurities introduced into octahedrally shaped ND of C₈₄H₆₄ size behave in a very similar manner as those embedded into a smaller octahedral model of C₃₅H₃₆ size. Since cubic and tetrahedral NDs-C₅₄H₄₈ and C₅₁H₅₂-demonstrate a wider range of g-shift values than octahedral systems, the g-tensor dependence on different shapes can be considered as more pronounced. However, a different surface functionalization scheme, namely, fluorination, results in a much larger variation of the g-shifts, pointing to a significant impact the F atoms have on the local environment of the unpaired electrons in C₃₅F₃₆. A partial surface functionalization of C₃₅H₃₆ with benzoic acid and aniline groups indicates that, in some special cases, these linkers might induce a noticeable spin density redistribution which in turn substantially modifies the g-shift values of the system. Additional infrared (IR) spectra calculations show that some of paramagnetic defects in C₃₅H₃₆ and C₃₅F₃₆ possess clearly expressed signatures which could be useful while analyzing the experimental IR spectra of NDs.

On the magnetic properties of nanodiamonds: Electronic g-tensor calculations

The electronic g-tensor calculations are carried out for various paramagnetic defects introduced into hydrogenated diamond nanocrystal C₃₅H₃₆, showing that such a system can be successfully used to model magnetic properties of nanodiamonds (NDs) with paramagnetic centers containing no vacancies. In addition, it is revealed that, depending on the geometric positions in ND, paramagnetic centers of the same type produce noticeable variations of the g-tensor values. A side-by-side comparison of the performance of effective nuclear charge and spin-orbit mean field (SOMF) approaches indicates that the latter is more sensitive to the quality of basis sets, especially concerning diffuse functions, the inclusion of which is found to be nonbeneficial. What is more, the SOMF method also exhibits a much more pronounced gauge-origin dependence. Compared to electronic charge centroid, spin centers (SCs) demonstrate a superior suitability as gauge origins, providing a better agreement with diamagnetic and paramagnetic contributions of g-tensor obtained employing gauge-including atomic orbitals (GIAOs). Therefore, SCs can be recommended for the g-tensor calculations of NDs whenever GIAOs are not available.

Computational study of the DPAP molecular rotor in various environments: from force field development to molecular dynamics simulations and spectroscopic calculations

Fluorescent molecular rotors (FMRs) belong to an important class of environment-sensitive dyes capable of acting as nanoprobes in the measurement of viscosity and polarity of their micro-environment. FMRs have found widespread applications in various research fields, ranging from analytical to biochemical sciences, for example in intracellular imaging studies or in volatile organic compound detection. Here, a computational investigation of a recently proposed FMR, namely 4-(diphenylamino)phthalonitrile (DPAP), in various chemical environments is presented. A purposely developed molecular mechanics force field is proposed and then applied to simulate the rotor in a high- and low-polar solvent (i.e., acetonitrile, tetrahydrofuran, o-xylene and cyclohexane), a polymer matrix and a lipid membrane. Subtle effects of the molecular interactions with the embedding medium, the structural fluctuations of the rotor and its rotational dynamics are analyzed in some detail. The results correlate with a previous work, thus supporting the reliability of the model, and provide further insights into the environment-specific properties of the dye. In particular, it is shown how molecular diffusion and rotational correlation times of the FMR are affected by the surrounding medium and how the molecular orientation of the dye becomes anisotropic once immersed in the lipid bilayer. Moreover, a qualitative correlation between the FMR rotational dynamics and the fluorescence lifetime is detected, a result in line with the observed viscosity dependence of its emission. Finally, optical absorption spectra are computed and successfully compared with their experimental counterparts.

All-atom molecular dynamics simulations of spin labelled double and single-strand DNA for EPR studies

Prediction of motional EPR spectra of spin labelled DNA structures from fully atomistic MD simulations.

Simulation of Vacuum UV Absorption and Electronic Circular Dichroism Spectra of Methyl Oxirane: The Role of Vibrational Effects

Vibrationally resolved one-photon absorption and electronic circular dichroism spectra of (R)-methyl oxirane were calculated with different electronic and vibronic models selecting, through an analysis of the convergence of the results, the best compromise between reliability and computational cost. Linear-response TD-DFT/CAM-B3LYP/SNST electronic computations in conjunction with the simple vertical gradient vibronic model were chosen and employed for systematic comparison with the available experimental data. Remarkable agreement between simulated and experimental spectra was achieved for both one-photon absorption and circular dichroism concerning peak positions, relative intensities, and general spectral shapes considering the computational efficiency of the chosen theoretical approach. The significant improvement of the results with respect to smearing of vertical electronic transitions by phenomenological Gaussian functions and the possible inclusion of solvent effects by polarizable continuum models at a negligible additional cost paves the route toward the simulation and analysis of spectral shapes of complex molecular systems in their natural environment.

Tuning of dye optical properties by environmental effects: a QM/MM and experimental study

The present work is aimed at a deeper investigation of two recently synthesized heteroaromatic fluorophores by means of a computational multilayer approach, integrating quantum mechanics (QM) and molecular mechanics (MM).

The present work is aimed at a deeper investigation of two recently synthesized heteroaromatic fluorophores by means of a computational multilayer approach, integrating quantum mechanics (QM) and molecular mechanics (MM). In particular, dispersion of the title dyes in a polymer matrix is studied in connection with potential applications as photoactive species in luminescent solar concentrators (LSCs). Molecular dynamics simulations, based on accurate QM-derived force fields, reveal increased stiffness of these organic dyes when going from CHCl₃ solution to the polymer matrix. QM/MM computations of UV spectra for snapshots extracted from MD simulations show that this different flexibility permits explaining the different spectral shapes obtained experimentally for the two different environments. Moreover, the general spectroscopic trends are reproduced well by static computations employing a polarizable continuum description of environmental effects.

Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases

The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.

The Virtual Multifrequency Spectrometer: a new paradigm for spectroscopy

On going developments of hardware and software are changing computational spectroscopy from a strongly specialized research area to a general tool in the inventory of most researchers. Increased interactions between experimentally-oriented users and theoretically-oriented developers of new methods and models would result in more robust, flexible and reliable tools and studies for the systems of increasing complexity, which are of current scientific and technological interest. This is the philosophy behind this review, which presents the development of a so-called virtual multi-frequency spectrometer (VMS) including state-of-the-art approaches in a user-friendly frame. The current status of the VMS tool will be illustrated by a number of case studies with special reference to infrared and UV-vis regions of the electro-magnetic spectrum including also chiral spectroscopies. Only the basic theoretical background will be provided avoiding explicit equations as far as possible, and pointing out the most recent advancements beyond the standard rigid-rotor harmonic-oscillator model coupled to vertical electronic excitation energies.

Anharmonic simulations of the vibrational spectrum of sulfated compounds: application to the glycosaminoglycan fragment glucosamine 6-sulfate

Anharmonic behavior of sulfated glucosamine resolved by hybrid GVPT2 approach.

A general time-dependent route to resonance-Raman spectroscopy including Franck-Condon, Herzberg-Teller and Duschinsky effects

We present a new formulation of the time-dependent theory of Resonance-Raman spectroscopy (TD-RR). Particular attention has been devoted to the generality of the framework and to the possibility of including different effects (Duschinsky mixing, Herzberg-Teller contributions). Furthermore, the effects of different harmonic models for the intermediate electronic state are also investigated. Thanks to the implementation of the TD-RR procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects have been included in an effective way. The reliability and stability of our TD-RR implementation are validated against our previously proposed and well-tested time-independent procedure. Practical applications are illustrated with some closed- and open-shell medium-size molecules (anthracene, phenoxyl radical, benzyl radical) and the simulated spectra are compared to the experimental results. More complex and larger systems, not limited to organic compounds, can be also studied, as shown for the case of Tris(bipyridine)ruthenium(II) chloride.

New developments of a multifrequency virtual spectrometer: stereo-electronic, dynamical, and environmental effects on chiroptical spectra

Computational spectroscopy has recently evolved from a field reserved for specialists toward a general tool allowing interpretations and analyses of experimental results. However, the current practice of providing tables of transitions for rigid geometries, possibly tuned by phenomenological broadening, is by far too naive. In order to improve this situation in the last few years we have been developing a general, robust, and user-friendly virtual spectrometer (VS) able to complement experimental studies for complex systems in condensed phases. The VS is based on flexible graphical pre- and postprocessing tools interfaced with general number-crunching software. This last tool is rooted in several electronic structure methodologies (DFT, TD-DFT, post-Hartree-Fock), powerful discrete/continuum models for describing environmental effects, and general vibrational and vibronic models. These last topics are the main focus of this work, which sketches our latest developments related to effective inclusion of anharmonic contributions, together with time-independent and/or time-dependent descriptions of vibronic transitions including Franck-Condon, Herzberg-Teller, and Duschinsky effects. Some test cases are described in some detail with the aim of showing the role of different effects in ruling vibrational (VCD) and electronic (ECD, CPL) chiral spectroscopies.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

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.

First principles calculation of inhomogeneous broadening in solid-state cw-EPR spectroscopy

We present a scheme for the first-principles calculation of EPR lineshapes for continuous-wave-EPR spectroscopy (cw-EPR) of spin centers in complex chemical environments. We specifically focus on poorly characterized systems, e.g. powders and frozen glasses with variable microsolvation structures. Our approach is based on ab initio molecular dynamics simulations and ab initio calculations of the ensemble of g- and A-tensors along the trajectory. The method incorporates temperature effects as well as the full anharmonicity of the intra- and intermolecular degrees of freedom of the system. We apply this scheme to compute the lineshape of a prototypical spin probe, the nitrosodisulfonate dianionic radical (Fremy's salt), dissolved in a 50 : 50 mixture of water and methanol. We are able to determine the specific effect of variations of local solvent composition and microsolvation structure on the cw-EPR lineshape. Our molecular dynamics reveal a highly anisotropic solvation structure with distinct spatial preferences for water and methanol around Fremy's salt that can be traced back to a combination of steric and polar influences. The overall solvation structure and conformational preferences of Fremy's salt as found in our MD simulations agree very well with the results obtained from EPR and orientation-selective ENDOR spectroscopy performed on the frozen glass. The simulated EPR lineshapes show good agreement with the experimental spectra. When combined with our MD results, they characterize the lineshape dependence on local morphological fluctuations.

General Time Dependent Approach to Vibronic Spectroscopy Including Franck-Condon, Herzberg-Teller, and Duschinsky Effects

An effective time-dependent (TD) approach to compute vibrationally resolved optical spectra from first principles is presented for the computation of one-photon electronic spectra induced by either electric or magnetic transition dipoles or by their mutual interaction, namely absorption, emission, and circular dichroism. Particular care has been devoted to generality, modularity, and numerical stability including all the contributions that play a role at the harmonic level of approximation, namely Franck-Condon, Herzberg-Teller, and Dushinsky (i.e., mode mixing) effects. The implementation shares the same general framework of our previous time-independent (TI) model, thus allowing an effective integration between both approaches with the consequent enhancement of their respective strengths (e.g., spectrum completeness and straightforward account of temperature effects for the TD route versus band resolution and assignment for the TI route) using a single set of starting data. Implementation of both models in the same general computer program allows comprehensive studies using several levels of electronic structure description together with effective account of environmental effects by atomistic and/or continuum models of different sophistication. A few medium-size molecules (furan, phenyl radical, anthracene, dimethyloxirane, coumarin 339) have been studied in order to fully validate the approach.

Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds

The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.

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.

π-Stacking effects on the EPR parameters of a prototypical DNA spin label

The character and value of spin labels for probing environments like double-stranded DNA depend on the degree of change in the spin Hamiltonian parameters of the spin label induced by the environment. Herein we provide a systematic theoretical investigation of this issue, based on a density functional theory method applied to a spin labeled DNA model system, focusing on the dependence of the EPR properties of the spin label on the π stacking and hydrogen bonding that occur upon incorporating the spin label into the selected base pair inside DNA. It is found that the EPR spin Hamiltonian parameters of the spin label are only negligibly affected by its incorporation into DNA, when compared to its free form. This result gives a theoretical ground for the common empirical assumption regarding the behaviour of spin Hamiltonian parameters made in EPR based measurements of the distance between spin labels incorporated into DNA.

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.

Toward Ab Initio Anharmonic Vibrational Circular Dichroism Spectra in the Condensed Phase

The first implementation and calculation of anharmonic VCD rotational strengths for solvated systems is reported. Our approach, rooted in the polarizable continuum model (PCM) and in the second-order vibrational perturbation theory (VPT2), permits not only correction for anharmonicity in the signals associated with fundamental transitions but also calculation of rotational strengths of overtones and combination bands. This allows for a more physically consistent comparison between experiment and calculations together with the analysis of spectral regions dominated by anharmonic effects. The developed model is applied to a few test cases, and the computational outcomes are directly compared with experimental data.

Molecular structure refinement by direct fitting of atomic coordinates to experimental ESR spectra

An attempt is made to bypass spectral analysis and fit internal coordinates of radicals directly to experimental liquid- and solid-state electron spin resonance (ESR) spectra. We take advantage of the recently introduced large-scale spin dynamics simulation algorithms and of the fact that the accuracy of quantum mechanical calculations of ESR parameters has improved to the point of quantitative correctness. Partial solutions are offered to the local minimum problem in spectral fitting and to the problem of spin interaction parameters (hyperfine couplings, chemical shifts, etc.) being very sensitive to vibrational excursions from the equilibrium geometry.

Theoretical study of the conformational and optical properties of a fluorescent dye. A step toward modeling sensors grafted on polymer structures

The conformational and photophysical properties in toluene solution of a naphthoic acid derivative, namely 4-naphthoyloxy-1-methoxy-2,2,6,6-tetramethylpiperidine, which can simulate a pendant active group grafted on apolar polymer structures, were studied by means of calculations based on the density functional theory (DFT) and its time-dependent extension (TD-DFT) coupled to a polarizable continuum model (PCM) of the surrounding medium. The conformational landscape was exhaustively explored and the conformers responsible for the absorption and emission spectra were identified through well established computational procedures based on a tuned combination of functional (PBE0) and basis set (N07D) which reproduced satisfactorily the experimental absorption and emission wavelengths.

Reliable structural, thermodynamic, and spectroscopic properties of organic molecules adsorbed on silicon surfaces from computational modeling: the case of glycine@Si(100)

Chemisorption of glycine on Si(100) has been studied by an integrated computational strategy based on perturbative anharmonic computations employing geometries and harmonic force fields evaluated by hybrid density functionals coupled to purposely tailored basis sets. It is shown that such a strategy allows the prediction of spectroscopic properties of isolated and chemisorbed molecules with comparable accuracy, paving the route toward a detailed analysis of surface-induced changes of glycine vibrational spectra.

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.

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.