Investigating the (Poly)Radicaloid Nature of Real-World Organic Compounds with DFT-Based Methods.

Recent advances in the synthesis of stable organic (open-shell) polyradicaloids have opened their application as active compounds for emerging technologies. These systems typically exhibit small energy differences between states with different spin multiplicities, which are intrinsically difficult to calculate by theoretical methods. We thus apply here some DFT-based variants (FT-DFT, SF-DFT, and SF-TDDFT) on a test set of large and real-world molecules, as test systems for which such energy differences are experimentally available, also comparing systematically with RAS-SF results to infer if shortcomings of previous DFT applications are corrected. Additionally, we explore the spin-spin contribution to the ZFS tensor, of high interest for EPR spectroscopy, and derive the spatial extent of the corresponding (photoexcited) triplet state.

Quadratic integrand double-hybrid made spin-component-scaled.

We propose two analytical expressions aiming to rationalize the spin-component-scaled (SCS) and spin-opposite-scaled (SOS) schemes for double-hybrid exchange-correlation density-functionals. Their performances are extensively tested within the framework of the nonempirical quadratic integrand double-hybrid (QIDH) model on energetic properties included into the very large GMTKN30 benchmark database, and on structural properties of semirigid medium-sized organic compounds. The SOS variant is revealed as a less computationally demanding alternative to reach the accuracy of the original QIDH model without losing any theoretical background.

Tuning the Electron Acceptor in Phthalocyanine-Based Electron Donor-Acceptor Conjugates.

Zinc phthalocyanines (ZnPc) have been attached to the peri-position of a perylenemonoimide (PMI) and a perylenemonoanhydride (PMA), affording electron donor-acceptor conjugates 1 and 2, respectively. In addition, a perylene-monoimide-monoanhydride (PMIMA) has been connected to a ZnPc through its imido position to yield the ZnPc-PMIMA conjugate 10. The three conjugates have been studied for photoinduced electron transfer. For ZnPc-PMIMA 10, electron transfer occurs upon both ZnPc and PMIMA excitation, giving rise to a long-lived (340 ps) charge-separated state. For ZnPc-PMI 1 and ZnPc-PMA 2, stabilization of the radical ion pair states by using polar media is necessary. In THF, photoexcitation of either ZnPc or PMI/PMA produces charge-separated states with lifetimes of 375 and 163 ps, respectively.

Single Lévy states-disorder induced energy funnels in molecular aggregates.

Using fluorescence super-resolution microscopy we studied simultaneous spectral, spatial localization, and blinking behavior of individual 1D J-aggregates. Excitons migrating 100 nm are funneled to a trap appearing as an additional red-shifted blinking fluorescence band. We propose that the trap is a Frenkel exciton state formed much below the main exciton band edge due to an environmentally induced heavy-tailed Lévy disorder. This points to disorder engineering as a new avenue in controlling light-harvesting in molecular ensembles.

Assembling a phthalocyanine and perylenediimide donor-acceptor hybrid through a platinum(II) diacetylide linker.

In a novel electron-donor-acceptor conjugate, phthalocyanine (Pc) and perylenediimide (PDI) are connected through a trans-platinum(II) diacetylide linker to yield Pc-Pt-PDI 1. In the ground state, the presence of Pt(II) disrupts the electronic communication between the two electroactive components, as revealed by UV/Vis spectroscopy and electrochemical studies. The photophysical behavior of 1 is compared with that of the corresponding Pc-PDI electron-donor-acceptor conjugate 2 in terms of charge separation and charge recombination. The insertion of Pt(II) between Pc and PDI impacts the results in a longer-lived Pc(.) (+) /PDI(.) (-) radical ion-pair state. In addition, the intermediately formed Pc triplet excited state is formed with higher quantum yields in 1 than in 2.

Assembling phthalocyanine dimers through a platinum(II) acetylide linker.

Phthalocyanine (Pc) dimers connected through trans-platinum(II) diacetylide linkers have been prepared by reaction of the corresponding ethynylphthalocyanines with trans-bis(triethylphosphine)platinum(II) chloride. Special emphasis was placed on the analysis of the ground- and excited-state features of these compounds in relation to butadiyne-bridged Pc dimers and the corresponding monomers. Both Zn(II)-containing Pc dimers exhibit long-lived triplet excited states. The insertion of σ-bonded trans-platinum(II) diacetylide spacers decoupled the two Pc groups and led to an appreciable acceleration (by a factor of up to 10) of the radiative and nonradiative decay rate of the singlet and triplet excited states.

Synthesis, characterization, and photoinduced energy and electron transfer in a supramolecular tetrakis (ruthenium(II) phthalocyanine) perylenediimide pentad.

Metal coordination was probed as a versatile approach for designing a novel electron donor/acceptor hybrid [PDIpy(4){Ru(CO)Pc}(4)] (1), in which four pyridines placed at the bay region of a perylenediimides (PDIpy(4)) coordinate with four ruthenium phthalocyanine units [Ru(CO)Pc]. This structural motif was expected to promote strong electronic coupling between the electron donors and the electron acceptor, a hypothesis that was confirmed in a full-fledged physicochemical investigation focusing on the ground and excited state reactivities. As far as the ground state is concerned, absorption and electrochemical assays indeed reveal a notable redistribution of electron density, that is, from the electron-donating [Ru(CO)Pc] to the electron-accepting PDIpy(4). The most important thing to note in this context is that both the [Ru(CO)Pc] oxidation and the PDIpy(4) reduction are rendered more difficult in 1 than in the individual building blocks. Likewise, in the excited state, strong electronic communication is the inception for a rapid charge-transfer process in photoexcited 1. Regardless of exciting [Ru(CO)Pc] or PDIpy(4), spectral characteristics of the [RuPc] radical cation (broad absorptive features from 425 to 600 nm with a maximum at 575 nm, as well as a band centered at 725 nm) and of the PDI radical anion (780 nm maximum) emerge. The correspondingly formed radical ion pair state lasts for up to several hundred picoseconds in toluene, for example. On the other hand, employing more polar solvents, such as dichloromethane, destabilizes the radical ion pair state.

Conductance enhancement in nanographene-gold junctions by molecular pi-stacking.

First-principles calculations on gold-pentacene-gold and several gold-circumacene-gold nanojunctions indicate that their low-bias conductance is due to the onset of a HOMO-derived resonance, thus being quite sensitive to the detailed interaction between the molecule and the gold leads. It is also found that such interaction is dominated by the electrophilic binding of Au to the (circum)acene, in agreement with previous theoretical and experimental results on pentacene. Therefore, the alignment of the HOMO resonance with the Fermi level, and thus the conductance, increases as the ionization potential and the HOMO-LUMO gap of the molecular arrangement diminish. It is shown here that both quantities are inversely proportional to the molecular length and the number of molecules present on a pi stack. It is also found that the conductance depends dramatically on the amount of pi overlap between the molecules in the stack, as well as on the particular disposition of the metallic tips with respect to the molecule. The conclusions reached point toward pi-stacked arrangements of large circumacenes as potential candidates to build useful nanodevices for molecular electronics made out of nanographene-based materials.

Photoinduced electron and energy transfer in aryldihydropyridines.

Dimethyl 2,6-dimethyl-1,4-dihydropyridine-3,5-dicarboxylates (Hantzsch DHPs) fluoresce weakly in fluid solution. However, these compounds exhibit an efficient fluorescence both in a viscous medium (glycerin) at room temperature and in a glassy matrix at 77 K (but no phosphorescence, since ISC is negligible). DHPs bearing an aryl group in position 4 have been synthesized. These contain two different pi systems separated by an sp(3) carbon (DHP-Ar dyads). The occurrence of energy and electron transfer processes between the chromophores is investigated through luminescence measurements. In particular, when (3)Ar emits at a slow rate (e.g., Ar = phenanthryl) or not at all (Ar = nitrophenyl) the intradyad forward/backward electron transfer sequence offers a path for arriving at the DHP-localized triplet and the corresponding phosphorescence is observed. When (3)Ar emits at a faster rate (Ar = acylphenyl), the phosphorescence from either of the two localized triplets, (3)Ar or (3)DHP, can be observed depending on lambda(exc). When the aryl group has a triplet energy lower than that of (3)DHP, this functions as emitting (4-cyano-1-naphthyl) or nonemitting (MeO(2)CCH horizontal lineCHC(6)H(4)) energy sink. The results document the possibility of building tailor-made Hantzsch aryldihydropyridines as versatile photoactivated dyads.

Using circumacenes to improve organic electronics and molecular electronics: design clues.

Theoretical modeling is used here to ascertain the potential use of circumacenes to improve the transport parameters of pi-conjugated materials acting as: (i) the layered molecular constituent for organic electronic devices; and (ii) the molecular component of gold-molecule-gold nanobridges for molecular electronic device use. It is concluded that, to a first approximation, the molecular length or, alternatively, the HOMO-LUMO gap (HOMO: highest occupied molecular orbital; LUMO: lowest unoccupied molecular orbital) can be used to relate the two transport regimes usually found in these two fields, thus serving as a key design parameter for guaranteeing good performance of circumanthracene for both regimes. It is also clearly established that going beyond this simple relationship requires knowledge of the detailed molecule-contact geometry of the molecular nanobridge, and how its tremendous impact on the binding strength and the conductance prevents blind extrapolation of results obtained for molecular nanobridges built by means of different experimental set-ups.

A tightly coupled bis(zinc(II) phthalocyanine)-perylenediimide ensemble to yield long-lived radical ion pair states.

A conjugated donor-acceptor array composed of two phthalocyanines connected to the bay region of a perylenediimide is assembled by using palladium chemistry. Excitation of the phthalocyanine produces a nanosecond lived charge-separated state.

Uncovering transport properties of 4,4'-bipyridine/gold molecular nanobridges.

Here, the fascinating connection between the chemical and the transport properties of recently fabricated 4,4'-bipyridine/gold nanobridges is addressed. By means of first-principles ab initio calculations, the remarkable reproducibility of the 4,4'-bipyridine conductance properties is explained as the combined result of (i) the bonding of the molecule to the metallic leads through hybridization between the 4,4'-bipyridine highest occupied molecular orbitals and lowest unoccupied molecular orbitals (LUMOs) with s and d orbitals at low-coordination gold atoms, (ii) the limited number of molecule-lead arrangements due to gold-hydrogen steric repulsions, and (iii) the electron transmission through a LUMO-derived resonance, whose positioning with respect to the Fermi level determines which of the above arrangements yields nonnegligible conductance. Structural and electronic interpretations to the stepped dependence reported for the electronic transport of 4,4'-bipyridine as a function of the distance between the gold tips are also given.

Self-assembly of octachloroperylene diimide into 1D rods and 2D plates by manipulating the growth kinetics for waveguide applications.

One-dimensional (1D) rods and 2D hexagonal plates of octachloroperylene diimide (Cl8-PTCDI) have been selectively prepared by controlling the growth kinetic processes. Both ensemble and single-particle spectroscopy clarify that 1D rods and 2D plates have shape dependent optical waveguiding properties.

Bay-substituted perylene bisimide dye with an undistorted planar scaffold and outstanding solid state fluorescence properties.

Single crystal X-ray analysis reveals a flat perylene π-scaffold for a perylene bisimide bearing bulky 2,6-diphenylphenoxy substituents at the 1,7 bay positions. The flat structure provides sharp vibronic progressions in the absorption and fluorescence spectra, the sterical shielding outstanding fluorescence quantum yields in the solid state.

Synthesis and photophysical properties of a hydrogen-bonded phthalocyanine-perylenediimide assembly.

A supramolecular phthalocyanine-perylenediimide donor-acceptor array has been assembled by using a melamine/perylenediimide motif. Photoexcitation of the perylenediimide component affords transduction of singlet excited state energy to the energetically lower lying phthalocyanine.

Merging multiconfigurational wavefunctions and correlation functionals to predict magnetic coupling constants.

We study the performance of different approaches that combine multiconfigurational wavefunctions with correlation functionals for the calculation of magnetic coupling constants of several materials and molecules. The systems under study include four antiferromagnetic materials: NiO, KNiF(3), K(2)NiF(4) and La(2)CuO(4); two biradicals: alpha-4-Dehydrotoluene and 1,1',5,5'-Tetramethyl-6,6'-dioxo-3,3'-biverdazyl; two molecular complexes: [Cu(2)Cl(6)](-2) and Copper(II) acetate monohidrate; and the prototypical H-He-H system. On average, the best results are obtained with a recently proposed method [Phys. Rev. A 75, 012503 (2007)] that estimates the correlation energy of density functionals from a pair of alternative spin densities built from the natural orbitals and occupation numbers of the multiconfigurational wavefunction.

Combining two-body density correlation functionals with multiconfigurational wave functions using natural orbitals and occupation numbers.

We propose a procedure that combines multiconfigurational (MC) wave functions with two-body density correlation functionals by transforming the latter into functionals of the MC natural orbitals and occupation numbers. The method is tested with the spectroscopic constants of a set of 11 diatomics, the diradical-involved automerization barrier of cyclobutadiene, the energy difference between triplet and open-shell singlet states in He and the methylene molecule, and the magnetic coupling constants of several systems, such as NiO, KNiF(3), K(2)NiF(4), La(2)CuO(4), alpha-4-dehydrotoluene, 1,1('),5,5(')-tetramethyl-6,6(')-dioxo-3,3(')-biverdazyl, [Cu(2)Cl(6)](-2), copper(II) acetate monohidrate and H-He-H. The procedure is applied to the Colle-Salvetti [Theor. Chim. Acta 37, 329 (1975); 53, 55 (1979)], functional and to a size-consistent functional depending on the on-top pair density (F1-5-N(eff)). On average, the best results are provided by the transformed F1-5-N(eff) [J. Chem. Phys. 114, 2022 (2001)] functional.

Torsional potential of 4,4'-bipyridine: ab initio analysis of dispersion and vibrational effects.

Ab initio calculations using restricted Hartree-Fock, second-order Møller-Plesset perturbation theory (MP2), density-functional theory (DFT), and coupled-cluster methods have been done to obtain the torsional potential-energy profile of the aza-aromatic molecule 4,4'-bipyridine. The torsional potential is evaluated adiabatically by fixing the normalized sum of the dihedral angles through the C-C inter-ring bond at several values along the torsional path and relaxing the remaining degrees of freedom. Previous discrepancies between MP2 and DFT internal rotation barrier heights are removed, and seen to be mostly due to the underestimation of the dispersion energy in the coplanar conformer by MP2 when using relatively small basis sets. The calculations indicate that the barrier height between the twisted global minimum and the 0 degrees conformer is around 1.5-1.8 kcal mol-1 while that corresponding to the 90 degrees one is about 2.0-2.2 kcal mol-1. This same relative energy ordering of the coplanar and perpendicular conformers was experimentally derived from nuclear magnetic resonance (NMR) measurements of 1H dipolar couplings on 4,4'-bipyridine solutions in a nematic liquid crystal, although the barrier heights are much lower than those estimated from NMR experiments in the gas phase. The DFT infrared spectrum and zero-point vibrational energy corrections to the torsional energy profile have also been calculated, the latter having a small influence on the torsional potential-energy profiles.

Correlation energy functionals dependent on an effective number of electrons: charged species and equilibrium geometries.

Recently proposed spin-dependent and spin-independent correlation energy functionals [Perez-Jimenez et al., J. Chem. Phys. 116, 10571 (2002)] based on an effective number of electrons N are extended to deal with charged systems. By introducing the concept of an effective atomic number Z analogous to N, the spin-dependent functional in combination with Becke's exchange [Becke, Phys. Rev. A 38, 3098 (1988)] yields a mean absolute error (MAE) of 5.4 kcal/mol for the 88 ionization potentials and 58 electron affinities included in the extended G2 set, and a MAE of 4.1 kcal/mol for the 312 data comprising the above plus the 148 enthalpies of formation of the extended G2 set and the 18 total energies of the neutral atoms H through Ar. Geometry optimizations performed on the 53 molecules of the G2-1 test set with the above combination of exchange and correlation functionals yield MAEs of 0.017 A and 1.5 degrees for the 68 bond lengths and 29 angles analyzed as compared with the experimental estimates.

Density functional theory with alternative spin densities: application to magnetic systems with localized spins.

A new method to improve the excess spin density obtained from unrestricted Hartree-Fock wave functions in terms of natural orbitals is proposed. Using this modified excess spin density to evaluate the correlation energy by means of density functionals leads to large improvements in the computed magnetic coupling constants of several materials without need to modify the exchange contribution. This is important because it reconciles the density functional theory description with the one provided by multi-determinant wave functions. Using the present approach, the leading contribution to the magnetic coupling constant arises from electron correlation effects. The performance of the new method is illustrated on various materials including high-critical-temperature superconductors parent compounds.