Interplay of Stereoelectronic and Enviromental Effects in Tuning the Structural and Magnetic Properties of a Prototypical Spin Probe: Further Insights from a First Principle Dynamical Approach

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.

Electronic and Vibrational Manifold of Tetracyanoethylene-Chloronaphthalene Charge Transfer Complex in Solution: Insights from TD-DFT and Ab Initio Molecular Dynamics

The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A 2018, 122 (14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+ polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.

Development and Validation of the B3LYP/N07D Computational Model for Structural Parameter and Magnetic Tensors of Large Free Radicals

Extensive calculations on a large set of free radicals containing atoms of the second and third row show that the B3LYP/N07D computational model provides remarkably accurate structural parameters and magnetic tensors at reasonable computational costs. The key of this success is the optimization of core-valence s functions for hyperfine coupling constants, while retaining (and even improving) the good performances of the parent 6-31+G(d,p) basis set for valence properties through reoptimization of polarization and diffuse p functions.

Toward an Understanding of the Ambiguous Electron Paramagnetic Resonance Spectra of the Iminoxy Radical from o-Fluorobenzaldehyde Oxime: Density Functional Theory and ab Initio Studies

Iminoxy radicals (R1R2C═N—O•) possess an inherent ability to exist as E and Z isomers. Although isotropic hyperfine couplings for the species with R1 = H allow one to distinguish between E and Z, unequivocal assignment of the parameters observed in the EPR spectra of the radicals without the hydrogen atom at the azomethine carbon to the right isomer is not a simple task. The iminoxyl derived from o-fluoroacetophenone oxime (R1 = CH3 and R2 = o-FC6H5) appears to be a case in point. Moreover, for its two isomers the rotation of the o-FC6H5 group brings into existence the syn and anti conformers, depending on the mutual orientation of the F atom and C═N—O• group, making a description of hyperfine couplings to structure even more challenging. To accomplish this, a vast array of theoretical methods (DFT, OO-SCS-MP2, QCISD) was used to calculate the isotropic hyperfine couplings. The comparison between experimental and theoretical values revealed that the E isomer is the dominant radical form, for which a fast interconversion between anti and syn conformers is expected. In addition, the origin of the significant AF increase with solvent polarity was analyzed.

Electronic excitation energies in dimers between radical ions presenting long, multicenter bonding

The formation of long, multicenter dimers between radical ions is usually monitored through UV-vis spectroscopy given the characteristic low-energy absorption band that they exhibit, not observed for the parent monomers. In this work, the performance of CASPT2, RASPT2, and TD-DFT methods for obtaining excitation energies of the long, multicenter bonded π-[TCNE]2(2-) and π-[TTF]2(2+) dimers has been addressed (TCNE = tetracyanoethylene; TTF = tetrathiafulvalene). The impact of the active space on the vertical electronic transitions computed at the RASPT2 and CASPT2 levels has been tested against experimentally observed absorption bands. Analogous tests have been carried out for a wide variety of density functionals within the TD-DFT formalism. Our calculations show that whereas CASPT2 predicts very accurate excitation energies for the π-[TCNE]2(2-), the mean absolute error for π-[TTF]2(2+) is higher for CASPT2 than for TD-DFT calculations, whenever pure density functionals or low % HF exchange hybrid functionals are used. Hybrid functionals with high % HF exchange (and thus RSH functionals) conduct to large errors on the excitation energies in both dimers. Furthermore, vertical electronic transitions are also obtained for 100 configurations extracted from a 45 ps molecular dynamics (CPMD) simulation aimed at providing an accurate description of the thermal fluctuation effects of a π-[TCNE]2(2-) dimer in dichloromethane. These thermal effects explain the shape of the experimental UV-vis spectrum, where the lowest HOMO → LUMO absorption presents a broad band and the following HOMO-1 → LUMO absorption exhibits a narrower band.

Theoretical Characterisation of Phosphinyl Radicals and Their Magnetic Properties: g Matrix

The g matrices (g tensors) of various phosphinyl radicals (R2 P(.) ) were calculated using the DFT and multireference configuration interaction (MRCI) methods. The g matrices were distinctly dependent on the molecular structure of the radical. To thoroughly examine this dependence, the contributions from individual atoms and excited states were calculated. The former revealed the gain from the phosphorus atom to be preeminent unless PO or PS bonds are present in the radical molecule. The contributions owing to excited states arising from electronic transitions between doubly occupied molecular orbitals and the SOMO were clearly positive, as in the case of semiquinone and niroxide radicals. The transitions from the phosphorus lone pair were of paramount importance. Surprisingly, unlike for semiquinones and nitroxides, a significant negative contribution was observed from excitations from the SOMO to unoccupied molecular orbitals. For radicals with PO bonds, this contribution to the g2 component was dominant.

Dynamical effects on the magnetic properties of dithiazolyl bistable materials

The magnetic properties of molecule-based magnets are commonly rationalized by considering only a single nuclear configuration of the system under study (usually an X-ray crystal structure). Here, by means of a computational study, we compare the results obtained using such a static approach with those obtained by explicitly accounting for thermal fluctuations, and uncover the serious limitations of the static perspective when dealing with magnetic crystals whose radicals undergo wide-amplitude motions. As a proof of concept, these limitations are illustrated for the magnetically bistable 1,3,5-trithia-2,4,6-triazapentalenyl (TTTA) material. For its high-temperature phase at 300 K, we show that nuclear dynamics induce large fluctuations in the magnetic exchange interactions (JAB) between spins (up to 1000% of the average value). These deviations result in a ∼20% difference between the 300 K magnetic susceptibility computed by explicitly considering the nuclear dynamics and that computed using the X-ray structure, the former being in better agreement with the experimental data. The unveiled strong coupling between JAB interactions and intermolecular vibrations reveals that considering JAB as a constant value at a given temperature (as always done in molecular magnetism) leads to a flawed description of the magnetism of TTTA. Instead, the physically relevant concept in this case is the statistical distribution of JAB values. The discovery that a single X-ray structure is not adequate enough to interpret the magnetic properties of TTTA is also expected to be decisive in other organic magnets with dominant exchange interactions propagating through labile π-π networks.

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.

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.

Theoretical and experimental studies of the spin trapping of inorganic radicals by 5,5-dimethyl-1-pyrroline N-oxide (DMPO). 3. Sulfur dioxide, sulfite, and sulfate radical anions

Radical forms of sulfur dioxide (SO(2)), sulfite (SO(3)(2-)), sulfate (SO(4)(2-)), and their conjugate acids are known to be generated in vivo through various chemical and biochemical pathways. Oxides of sulfur are environmentally pervasive compounds and are associated with a number of health problems. There is growing evidence that their toxicity may be mediated by their radical forms. Electron paramagnetic resonance (EPR) spin trapping using the commonly used spin trap, 5,5-dimethyl-1-pyrroline N-oxide (DMPO), has been employed in the detection of SO(3)(•-) and SO(4)(•-). The thermochemistries of SO(2)(•-), SO(3)(•-), SO(4)(•-), and their respective conjugate acids addition to DMPO were predicted using density functional theory (DFT) at the PCM/B3LYP/6-31+G**//B3LYP/6-31G* level. No spin adduct was observed for SO(2)(•-) by EPR, but an S-centered adduct was observed for SO(3)(•-)and an O-centered adduct for SO(4)(•-). Determination of adducts as S- or O-centered was made via comparison based on qualitative trends of experimental hfcc's with theoretical values. The thermodynamics of the nonradical addition of SO(3)(2-) and HSO(3)(-) to DMPO followed by conversion to the corresponding radical adduct via the Forrester-Hepburn mechanism was also calculated. Adduct acidities and decomposition pathways were investigated as well, including an EPR experiment using H(2)(17)O to determine the site of hydrolysis of O-centered adducts. The mode of radical addition to DMPO is predicted to be governed by several factors, including spin population density, and geometries stabilized by hydrogen bonds. The thermodynamic data supports evidence for the radical addition pathway over the nucleophilic addition mechanism.

Effects of restricted rotations and dynamic averaging on the calculated isotropic hyperfine coupling constants of the bis-dimethyl and bis-Di(trifluoromethyl) nitroxide radicals

The rotational effects of the CH(3) and CF(3) groups on the electronic structure and nuclear hyperfine coupling constants (HFCCs) in dimethylnitroxide (DMNO*) and ditrifluoro-methynitroxide (TFMNO*) are investigated using the UB1LYP hybrid density functional method. The CH(3) and CF(3) HFCCs of both radicals are found to obey the McConnell relation during rotation. The two CH(3) groups of the DMNO* do not gear with each other, but the rotation of the first CH(3) group induces only a small rocking effect ( approximately 7 degrees ) in the second group. However, in TFMNO*, the fluorine atoms from different CF(3) groups are close enough so that the steric repulsion between them causes them to act as two interlocked gears, where one drives the other. Therefore, both CF(3) groups undergo full rotation. To the best of our knowledge, this is only the second example of "gearing" to be studied. Stabilization due to hyperconjugation is also a major factor that affects the magnitudes of the HFCCs of the CF(3) groups during rotational averaging. Stable configurations at specific CF(3) group orientations have a large overlap with the NO pi-electron cloud because the lobes of the hybridized p(sigma)(F(2)), p(sigma)(F(3)), p(sigma)(F(5)), and p(sigma)(F(6)) orbitals along the F-C bonds have cylindrical symmetry and are of the correct phases for hyperconjugation to occur. The calculated TFMNO* C(1)-N and C(2)-N bond orders range from 0.91 to 0.95 as the CF(3) groups are rotated. Therefore, the C-N bonds are essentially single bonds. This, in conjunction with the low rotational energy barrier of approximately 50 cm(-1), explains why the EPR intensities of the (19)F hyperfine splittings, in the range of 163-297 K, are characteristic of six equivalent rapidly rotating fluorine atoms. The TFMNO* out-of-plane NO vibrations induce additional s character at the (14)N nucleus. This increases the magnitude of the (14)N HFCC and decreases the (19)F HFCCs. As the temperature increases and because of mixing of the first excited out-of-plane vibrational state, the NO vibrational amplitudes also increase. This leads to an increased (14)N HFCC and decreased (19)F HFCCs, which is in agreement with experiment.

On the stability of X2NO radicals (X = F, Cl, Br, I)

The stability of X(2)NO radicals has been studied by investigating the three possible dissociative channels, namely X+XNO, X(2)+NO, and NX(2)+O. While all the radicals have been found stable with respect to the latter, the second pathway shows that Br(2)NO and I(2)NO are unstable with respect to dissociation. The first dissociative channel has been thoroughly investigated with the aim of understanding whether F(2)NO and Cl(2)NO are stable and how much. This implied the molecular structures and energies of X(2)NO and XNO, with X = F, Cl, to be computed at high level of theory. The coupled cluster ansatz in conjunction with hierarchical series of basis sets has been employed, thus accounting for extrapolation to the complete basis-set limit. Core correlation as well as higher excitations in the electronic-correlation treatment have also been taken into account. It is particularly noteworthy that explicit inclusion of quadruple excitations allowed to obtain for the first time equilibrium geometries of FNO and ClNO in full quantitative agreement with their experimental counterparts.

Alignment of the dye's molecular levels with the TiO(2) band edges in dye-sensitized solar cells: a DFT-TDDFT study

We present a theoretical study of the lineup of the LUMO of Ru(II)-polypyridyl (N3 and N719) molecular dyes with the conduction band edge of a TiO(2) anatase nanoparticle. We use density functional theory (DFT) and the Car-Parrinello scheme for efficient optimization of the dye-nanoparticle systems, followed by hybrid B3LYP functional calculations of the electronic structure and time-dependent DFT (TDDFT) determination of the lowest vertical excitation energies. The electronic structure and TDDFT calculations are performed in water solution, using a continuum model. Various approximate procedures to compute the excited state oxidation potential of dye sensitizers are discussed. Our calculations show that the level alignment for the interacting nanoparticle-sensitizer system is very similar, within about 0.1 eV, to that for the separated TiO(2) and dye. The excellent agreement of our results with available experimental data indicates that the approach of this work could be used as an efficient predictive tool to help the optimization of dye-sensitized solar cells.

Quantum mechanical computations and spectroscopy: from small rigid molecules in the gas phase to large flexible molecules in solution

Interpretation of structural properties and dynamic behavior of molecules in solution is of fundamental importance to understand their stability, chemical reactivity, and catalytic action. While information can be gained, in principle, by a variety of spectroscopic techniques, the interpretation of the rich indirect information that can be inferred from the analysis of experimental spectra is seldom straightforward because of the subtle interplay of several different effects, whose specific role is not easy to separate and evaluate. In such a complex scenario, theoretical studies can be very helpful at two different levels: (i) supporting and complementing experimental results to determine the structure of the target molecule starting from its spectral properties; (ii) dissecting and evaluating the role of different effects in determining the observed spectroscopic properties. This is the reason why computational spectroscopy is rapidly evolving from a highly specialized research field into a versatile and widespread tool for the assignment of experimental spectra and their interpretation in terms of chemical physical effects. In such a situation, it becomes important that both computationally and experimentally oriented chemists are aware that new methodological advances and integrated computational strategies are available, providing reliable estimates of fundamental spectral parameters not only for relatively small molecules in the gas phase but also for large and flexible molecules in condensed phases. In this Account, we review the most significant methodological contributions from our research group in this field, and by exploiting some recent results of their application to the computation of IR, UV-vis, NMR, and EPR spectral parameters, we discuss the microscopic mechanisms underlying solvent and vibrational effects on the spectral parameters. After reporting some recent achievements for the study of excited states by first principle quantum mechanical approaches, we focus on the treatment of environmental effects by means of mixed discrete-continuum solvent models and on effective methods for computing vibronic contributions to the spectra. We then discuss some new developments, mainly based on time-dependent approaches, allowing us to go beyond the determination of spectroscopic parameters toward the simulation of line widths and shapes. Although further developments are surely needed to improve the accuracy and effectiveness of several items in the proposed approach, we try to show that the first important steps toward a direct comparison between the results obtained in vitro and those obtained in silico have been made, making easier fruitful crossovers among experiments, computations and theoretical models, which would be decisive for a deeper understanding of the spectral behavior associated with complex systems and processes.

Toward effective and reliable fluorescence energies in solution by a new state specific polarizable continuum model time dependent density functional theory approach

A state specific (SS) model for the inclusion of solvent effects in time dependent density functional theory (TD-DFT) computations of emission energies has been developed and coded in the framework of the so called polarizable continuum model (PCM). The new model allows for a rigorous and effective treatment of dynamical solvent effects in the computation of fluorescence and phosphorescence spectra in solution, and it can be used for studying different relaxation time regimes. SS and conventional linear response (LR) models have been compared by computing the emission energies for different benchmark systems (formaldehyde in water and three coumarin derivatives in ethanol). Special attention is given to the influence of dynamical solvation effects on LR geometry optimizations in solution. The results on formaldehyde point out the complementarity of LR and SS approaches and the advantages of the latter model especially for polar solvents and/or weak transitions. The computed emission energies for coumarin derivatives are very close to their experimental counterparts, pointing out the importance of a proper treatment of nonequilibrium solvent effects on both the excited and the ground state energies. The availability of SS-PCM/TD-DFT models for the study of absorption and emission processes allows for a consistent treatment of a number of different spectroscopic properties in solution.

Halogen bonds between 2,2,6,6-tetramethylpiperidine-N-oxyl radical and CxHyFzI species: DFT calculations of physicochemical properties and comparison with hydrogen bonded adducts

Structural, thermodynamic, and magnetic properties of adducts between the 2,2,6,6-tetramethylpiperidine-N-oxyl radical and representative hydrogen and halogen bond donors in solution have been investigated by an integrated computational tool including hybrid density functionals and discrete-continuum solvent models. From a quantitative point of view, the computed values show a fair agreement with experiment when environmental effects are taken into the proper account. From a more general point of view, our analysis points out a number of analogies, but also some difference, between hydrogen and halogen bond, which have been interpreted in terms of the various effects tuning thermodynamic and spectroscopic parameters.

A quantum mechanical/molecular dynamics/mean field study of acrolein in aqueous solution: analysis of H bonding and bulk effects on spectroscopic properties

A novel molecular dynamics methodology recently proposed by our group [Rega et al., Chem. Phys. Lett. 422, 367 (2006)], which is based on an integrated hybrid potential rooted in high level quantum mechanical methods using localized basis functions and nonperiodic boundary conditions, has been applied to study acrolein in aqueous solution. The solute structural rearrangement and its hydrogen-bonding pattern due to the interactions with water have been analyzed in some detail. Moreover, the solvent effects on the UV n-->pi* vertical transition and on the NMR 13C and 17O shielding constants of acrolein have been investigated theoretically by performing a posteriori quantum mechanical calculations on a statistically significant number of snapshots extracted from both gas-phase and aqueous solution simulations. Results show that such effective computational strategy can be successfully used to improve our understanding, at atomic level, of important spectroscopic observables.

Interplay of Intrinsic, Environmental, and Dynamic Effects in Tuning the EPR Parameters of Nitroxides: Further Insights from an Integrated Computational Approach

The role of stereoelectronic, environmental, and short-time dynamic effects in tuning the hyperfine and gyromagnetic tensors of a prototypical nitroxide spin probe has been investigated by an integrated computational approach based on extended Lagrangian molecular dynamics and discrete-continuum solvent models. Trajectories were generated in two protic solvents as well as in the gas phase for reference; structural analysis of the dynamics, and comparison with optimized solute-solvent clusters, allowed for the identification of the prevailing solute-solvent hydrogen-bonding patterns and helped to define the strategy for the computation of magnetic parameters. This was performed in a separate step, on a large number of frames, by a high-level DFT approach coupling the PBE0 hybrid functional with a tailored basis set and with proper account of specific and bulk solvent effects. Remarkable changes in solvation networks are found on going from aqueous to methanol solution, thus providing a rationalization of indirect experimentally available evidence. The computed magnetic parameters are in satisfactory agreement with the available measured values and allow for an unbiased evaluation of the role of different effects in tuning the overall EPR observables. Apart from their intrinsic interest, our results pave the route toward the development of tunable detection protocols based on specific spectroscopic signatures.

Accurate calculation of the phenyl radical's magnetic inequivalency, relative orientations of its spin hamiltonian tensors, and its electronic spectrum

The phenyl radical's electronic structure, magnetic inequivalency, spin Hamiltonian tensor components, and the relative orientation of their principal axes are computed by Neese's coupled-perturbed Kohn-Sham hybrid density functional (CPKS-HDF) technique in a moderate amount of time without resorting to expensive post-Hartree-Fock techniques. The g tensor component values are in excellent agreement with those determined experimentally and differ by less than 370 ppm. The computed hydrogen nuclear hyperfine tensors, A(1H), are also found to be in very good agreement with their experimental counterparts. The correlation of the radical's electronic structure with its g and A numerical values corroborates that it has a 2A(1) ground state. In accordance with our previous studies on the equivalency of planar radicals that possess C(2v) symmetry, the in-plane g and A(1H) principal axes should not be parallel to one another. Consequently, the spatially equivalent ortho (1H(2), 1H(6)) and meta (1H(3), 1H(5)) proton pairs should be magnetically inequivalent. This was confirmed in both the present computations and the simulation of the EPR solid-state spectrum. To the best of our knowledge, this is the first aromatic in-plane sigma-type radical whose magnetic inequivalency is studied both computationally and experimentally. To properly interpret the radical's electronic excitation spectra, the spectroscopy-oriented dedicated difference configuration interaction (SORCI) procedure was employed. Aside from a slight overestimation, the method seems to be capable of reproducing the C(6)H(5)* electronic vertical excitation energies in the range of 0-50,000 cm(-1). These vertical excitations, in conjunction with the corresponding orbit and spin orbit matrix elements, were also used to compute the g tensor components, employing the sum-over-states technique. Due to the limited number of computed roots and excited states, the results were marginally inferior to those obtained using the CPKS-HDF method.

Hyperfine Coupling Tensors of the Benzosemiquinone Radical Anion from Car–Parrinello Molecular Dynamics

Based on Car-Parrinello ab initio molecular dynamics simulations of the benzosemiquinone radical anion in both aqueous solution and the gas phase, density functional calculations provide the currently most refined EPR hyperfine coupling (HFC) tensors of semiquinone nuclei and solvent protons. For snapshots taken at regular intervals from the molecular dynamics trajectories, cluster models with different criteria for inclusion of water molecules and an additional continuum solvent model are used to analyse the HFCs. These models provide a detailed picture of the effects of dynamics and of different intermolecular interactions on the spin-density distribution and HFC tensors. Comparison with static calculations allows an assessment of the importance of dynamical effects, and of error compensation in static DFT calculations. Solvent proton HFCs depend characteristically on the position relative to the semiquinone radical anion. A point-dipolar model works well for in-plane hydrogen-bonded protons but deviates from the quantum chemical values for out-of-plane hydrogen bonding.