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.

The influence of the calculation method in mimicking experimental (13)C NMR chemical shifts of 15 low-polarity natural products singularly containing 10-20 carbon atoms was investigated by employing different quantum chemistry approaches and basis sets, both in the preliminary geometry optimizations and in the following single-point (13)C GIAO calculations of the NMR chemical shifts. The geometries of the involved species were optimized at the PM3, HF, B3LYP and mPW1PW91 levels whereas the (13)C NMR parameters were determined at the HF, B3LYP and mPW1PW91 levels. Different combinations of basis sets were also tested. The consistency and efficiency of the considered combinations of geometry optimizations and GIAO (13)C NMR calculations were thoroughly checked by the analysis of statistical parameters concerning computed and experimental (13)C NMR chemical shift values.

The structures, relative stabilities as well as the rotational and vibrational spectra of the three lowenergy conformers of Pyruvic acid (PA) have been characterized using a state-of-the-art quantummechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semi-experimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg. for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol −1 . Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Ttand Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm −1 are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc-and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/ density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.

No abstract

The nitrogen isotropic hyperfine coupling constant (hcc) and the g tensor of a prototypical spin probe (di-tert-butyl nitroxide, DTBN) in aqueous solution have been investigated by means of an integrated computational approach including Car-Parrinello molecular dynamics and quantum mechanical calculations involving a discrete-continuum embedding. The quantitative agreement between computed and experimental parameters fully validates our integrated approach. Decoupling of the structural, dynamical, and environmental contributions acting onto the spectral observables allows an unbiased judgment of the role played by different effects in determining the overall experimental observables and highlights the importance of finite-temperature vibrational averaging. Together with their intrinsic interest, our results pave the route toward more reliable interpretations of EPR parameters of complex systems of biological and technological relevance.

The popular AMBER force-field has been extended to provide an accurate description of large and flexible nitroxide free-radicals in condensed phases. New atom types have been included, and relevant parameters have been fitted based on geometries, vibrational frequencies and potential energy surfaces computed at the DFT level for several different classes of nitroxides, both in vacuo and in different solvents. The resulting computational tool is capable of providing reliable structures, vibrational frequencies, relative energies and spectroscopic observables for large and flexible nitroxide systems, including those typically used as spin labels. The modified force field has been employed in the context of an integrated approach, based on classical molecular dynamics and discrete-continuum solvent models, for the investigation of environmental and short-time dynamic effects on the hyperfine and gyromagnetic tensors of PROXYL, TEMPO and INDCO spin probes. The computed magnetic parameters are in very good agreement with the available experimental values, and the procedure allows for an unbiased evaluation of the role of different effects in tuning the overall EPR observables.

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.

Calculations on a large set of free radicals containing atoms of the second and third row show that the computational model defined by the new N07D basis set and hybrid density functionals (B3LYP and PBE0) provides remarkably accurate g-tensor values at reasonable computational costs. Since in previous works it has been shown that the same computational model delivers reliable results also for structural parameters and hyperfine couplings, the route seems paved toward full a priori computation of EPR spectra of large free radicals both in vacuo and in condensed phases.