The relative energies of the energetically low-lying tautomers of pyridone, cytosine, uracil, thymine, guanine, and iso-cytosine are studied by a variety of different quantum chemical methods. In particular, we employ density functional theory (DFT) using the six functionals HCTH407, PBE, BP86, B-LYP, B3-LYP, and BH-LYP, and the ab initio methods Hartree-Fock (HF), standard second-order Møller-Plesset perturbation theory (MP2), an improved version of it (SCS-MP2), and quadratic configuration interaction including single and double excitations (QCISD) and perturbative triple corrections [QCISD(T)]. A detailed basis set study is performed for the formamide/formamidic acid tautomeric pair. In general, large AO basis sets of at least valence triple-zeta quality including f-functions (TZV) are employed, which are found to be necessary for an accurate energetic description of the various structures. The performance of the more approximate methods is evaluated with QCISD(T)/TZV(2df,2dp) data taken as reference. In general it is found that DFT is not an appropriate method for the problem. For the tautomers of pyridone and cytosine, most density functionals, including the popular B3-LYP hybrid, predict a wrong energetic order, and only for guanine, the correct sequence of tautomers is obtained with all functionals. Out of the density functionals tested, BH-LYP, which includes a rather large fraction of HF exchange, performs best. A consistent description of the nonaromatic versus aromatic tautomers seems to be a general problem especially for pure, nonhybrid functionals. Tentatively, this could be assigned to the exchange potentials used while the functional itself, including the correlation part, seems to be appropriate. Out of the ab initio methods tested, the new SCS-MP2 approach seems to perform best because it effectively reduces some outliers obtained with standard MP2. It outperforms the much more costly QCISD method and seems to be a very good compromise between computational effort and accuracy.

The absorption and fluorescence properties in a class of oligothiophene push-pull biomarkers are investigated with a long-range-corrected (LC) density functional method. Using linear response timedependent density functional theory (TDDFT), we calculate excitation energies, fluorescence energies, oscillator strengths, and excited-state dipole moments. To benchmark and assess the quality of the LC-TDDFT formalism, an extensive comparison is made between LC-BLYP excitation energies and approximate coupled cluster singles and doubles (CC2) calculations. When using a properly-optimized value of the range parameter, , we find that the LC technique provides an accurate description of charge-transfer excitations as a function of biomarker size and chemical functionalization. In contrast, we find that re-optimizing the fraction of Hartree Fock exchange in conventional hybrid functionals still yields an inconsistent description of excitation energies and oscillator strengths for the two lowest excited states in our series of biomarkers. The results of the present study emphasize the importance of a distance-dependent contribution of exchange in TDDFT for investigating excited-state properties.2

Full geometry optimizations at the dispersion-corrected DFT-BLYP level of theory were carried out for dimers and trimers of pyridine. The DFT-D interaction energies were checked against results from single-point SCS-MP2/aug-cc-pVTZ calculations. Three stacked structures and a planar H-bonded dimer were found to be very close in energy (interaction energies in the range from -3.4 to -4.0 kcal mol(-1)). Two T-shaped geometries are higher lying, by about 1 kcal mol(-1), which is explained by the more favorable electrostatic interactions in the stacked and H-bonded arrangements. The DFT-D approach has proved to be a reliable and efficient tool to explore the conformational space of aromatic van der Waals complexes and furthermore provides interaction energies with errors of less than 10-20 % of DeltaE. Comparisons with previous results obtained by using only partially optimized model geometries strongly indicate that unconstrained optimizations are mandatory in such weakly bonded low-symmetry systems.

The magnetic configuration of a nanostructure can be altered by an external magnetic field, by spin-transfer torque or by its magnetoelastic response. Here, we explore an alternative route, namely the possibility of switching the sign of the exchange coupling between two magnetic centres by means of an electric potential. This general effect, which we name electrostatic spin crossover, occurs in insulating molecules with super-exchange magnetic interaction and inversion symmetry breaking. As an example we present the case of a family of di-cobaltocene-based molecules. The critical fields for switching, calculated from first principles, are of the order of 1 V nm(-1) and can be achieved in two-terminal devices. More crucially, such critical fields can be engineered with an appropriate choice of substituents to add to the basic di-cobaltocene unit. This suggests that an easy chemical strategy for achieving the synthesis of suitable molecules is possible.

Full geometry optimizations at the dispersion corrected DFT-BLYP/TZV2P level of theory have been performed for dimers of azulene that may serve as a model system for the van der Waals complexes of polar pi systems. The structures and binding energies for 11 dimers are investigated in detail. The DFT-D interaction energies have been successfully checked against results from the accurate SCS-MP2/aug-cc-pVTZ approach. Out of the nine investigated stacked complexes, eight have binding energies larger than 7.4 kcal/mol (SCS-MP2) that exceed the value of 7.1 kcal/mol for the best naphthalene dimer. T-shaped arrangements (CH...pi) are significantly less stable. Two out of the three best structures have an antiparallel alignment of the monomer dipole moments in the complex, although the best ones with a parallel orientation are only about 0.5 kcal/mol less strongly bound which points to a minor importance of dipole-dipole interactions to binding. Quite surprisingly, the energetically lowest structure (DeltaE = -9.2 kcal/mol) corresponds to a situation where the two seven-membered rings are located almost on top of each other (7-7) and the long molecular axes are rotated against each other by 130 degrees. The 7-7 structural motif is found also in other energetically low-lying structures, and the expected 5-7 (two-side) arrangement is less strongly bound by about 2 kcal/mol. This can be explained by the electrostatic potential of azulene that only partially reflects the charge separation according to the common 4n + 2 pi electron rule. General rules for predicting stable van der Waals complexes of polar pi systems are discussed.