Quantum-classical calculation of the absorption and emission spectral shapes of oligothiophenes at low and room temperature by first-principle calculations.

We report a thorough computational characterization of the low- and room-temperature absorption and emission spectra of a series of oligothiophenes that contain between three and seven thiophene units. Our computational approach is based on time-dependent (TD) density functional calculations with the CAM-B3LYP functional. The effect of vibrations is included without resorting to any empirical parameters either at a fully quantum level or with a hybrid quantum-classical protocol. This latter approach is introduced to describe the relevant broadening effects in absorption at room temperature and is based on the partition of the vibrational modes into two sets: the inter-ring torsions treated at the anharmonic level in a classical way and the remaining modes described at the quantum level. The contribution of the quantum modes to the spectrum is computed by using a harmonic approximation, which accounts for Duschinsky mixing and changes in the vibrational frequencies associated with the electronic transition; a path-integral TD approach is adopted to account for the effect of temperature. The spectra simulated at low temperatures are in very good agreement with their experimental counterparts, which indicates that our calculations can quantitatively reproduce the effect of chain lengthening on the position and the shape of the spectra. Good agreement is also obtained at room temperature, for which we show that the classical description of the broadening, owing to the inter-ring torsions, reproduces the loss of the vibronic structure observed in the experiment and introduces only a slight overestimation of the spectral width.

Insights for an Accurate Comparison of Computational Data to Experimental Absorption and Emission Spectra: Beyond the Vertical Transition Approximation.

In this work we carefully investigate the relationship between computed data and experimental electronic spectra. To that end, we compare both vertical transition energies, EV, and characteristic frequencies of the spectrum like the maximum, ν(max), and the center of gravity, M(1), taking advantage of an analytical expression of M(1) in terms of the parameters of the initial- and final-state potential energy surfaces. After pointing out that, for an accurate comparison, experimental spectra should be preliminarily mapped from wavelength to frequency domain and transformed to normalized lineshapes, we simulate the absorption and emission spectra of several prototypical chromophores, obtaining lineshapes in very good agreement with experimental data. Our results indicate that the customary comparison of experimental ν(max) and computational EV, without taking into account vibrational effects, is not an adequate measure of the performance of an electronic method. In fact, it introduces systematic errors that, in the investigated systems, are on the order of 0.1-0.3 eV, i.e., values comparable to the expected accuracy of the most accurate computational methods. On the contrary, a comparison of experimental and computed M(1) and/or 0-0 transition frequencies provides more robust results. Some rules of thumbs are proposed to help rationalize which kind of correction one should expect when comparing EV, M(1), and ν(max).

Extension of the AMBER force-field for the study of large nitroxides in condensed phases: an ab initio parameterization.

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.

Interplay of stereo-electronic, environmental, and dynamical effects in determining the EPR parameters of aromatic spin-probes: INDCO as a test case.

An integrated computational strategy for the evaluation of reliable structures and magnetic properties of spin probes and spin labels has been extended to aromatic species. From an electronic point of view, delocalization of the unpaired electron density over aromatic moieties reduces significantly the computed nitrogen isotropic hyperfine coupling constant (A(N)) with respect to values characteristic of aliphatic nitroxides. Solvent effects in not too high polarity media are quite small, but not negligible. At this stage computed A(N) are lower than their experimental counterparts by more than 1 G. Inclusion of vibrational averaging effects by molecular dynamics simulations with a new reliable force field restores full agreement with experiment pointing out the limits of static approaches irrespective of the sophistication of the electronic quantum mechanical method. The generality and computational effectiveness of the proposed integrated approach paves the route toward a reliable analysis of the interplay of stereo-electronic, environmental, and dynamical effects in tuning the properties of large flexible magnetic systems of biological and technological interest.

Insight into the mechanism of action of plakortins, simple 1,2-dioxane antimalarials.

A multidisciplinary approach, based on molecular dynamics/mechanics, ab initio calculations, dynamic docking studies, and chemical reactions, has been employed to gain insight into the mechanism of the antimalarial action of plakortin and dihydroplakortin, simple 1,2-dioxanes isolated from the sponge Plakortis simplex. Our results show that these molecules, after interaction of the endoperoxide bond with Fe(ii), likely coming from the heme molecule, give rise to the formation of an oxygen radical, followed by rearrangement to give a carbon radical centered on the "western" alkyl side-chain. The carbon radicals generated on the side-chain, amenable for intermolecular reactions, should represent the toxic intermediates responsible for subsequent reactions leading to plasmodium death. The minimal structural requirements necessary for the activity of this class of antimalarial agents have been identified and discussed throughout the paper.

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.

Artarborol, a nor-caryophyllane sesquiterpene alcohol from Artemisia arborescens. stereostructure assignment through concurrence of NMR data and computational analysis.

A nor-caryophyllane derivative, artarborol, has been isolated from wormwood (Artemisia arborescens) and its stereostructure established by using a combination of chemical derivatization, NMR data, molecular modeling, and quantum-mechanical calculations. In particular, comparison of experimental 13C NMR data with a Boltzmann-weighed average of 13C NMR chemical shifts, calculated by ab initio DFT method, supported the stereochemical assignment.