π-Stacking in Charged Thiophene Oligomers

Crystal-structure of active layers of small molecule organic photovoltaics before and after solvent vapor annealing

It is demonstrated by a detailed structural analysis that the crystallinity and the efficiency of small molecule based organic photovoltaics can be tuned by solvent vapor annealing (SVA). Blends made of the small molecule donor 2,2′-[(3,3′″,3″″,4′-tetraoctyl[2,2′:5′,2″:5″,2′″:5′″,2″″-quinquethiophene]-5,5″″-diyl)bis[(Z)-methylidyne(3-ethyl-4-oxo-5,2-thiazolidinediylidene)]]bis-propanedinitrile (DRCN5T) and the acceptor [6,6]-phenyl C71 butyric acid methyl ester (PC71BM) were annealed using solvent vapors with either a high solubility for the donor (tetrahydrofuran), the acceptor (carbon disulfide) or both (chloroform). The samples were analyzed by grazing-incidence wide-angle X-ray scattering (GIWAXS), electron diffraction, X-ray pole figures, and time-of-flight secondary ion mass spectrometry (ToF-SIMS). A phase separation of DRCN5T and PC71BM is induced by SVA leading to a crystallization of DRCN5T and the formation of a DRCN5T enriched layer. The DRCN5T crystallites possess the two dimensional oblique crystal system with the lattice parameters a = 19.2 Å, c = 27.1 Å, and β = 111.1° for the chloroform case. No major differences in the crystal structure for the other solvent vapors were observed. However, the solvent choice strongly influences the size of the DRCN5T enriched layer. Missing periodicity in the [010]-direction leads to the extinction of all Bragg reflections with k ≠ 0. The annealed samples are randomly orientated with respect to the normal of the substrate (fiber texture).

Self-assembled structures of bent-shaped π-conjugated compounds: effect of siloxane groups for nano-segregation

The tuning of nanostructures is successfully achieved by introduction of siloxane unit to bithiophene-modified bent-shaped skeleton.

Revealing the Sequence-Structure-Electronic Property Relation of Self-Assembling π-Conjugated Oligopeptides by Molecular and Quantum Mechanical Modeling

Self-assembled nanoaggregates of π-conjugated synthetic peptides present a biocompatible and highly tunable alternative to silicon-based optical and electronic materials. Understanding the relationship between structural morphology and electronic properties of these assemblies is critical for understanding and controlling their mechanical, optical, and electronic responses. In this work, we combine all-atom classical molecular simulations with quantum mechanical electronic structure calculations to ascertain the sequence-structure-electronic property relationship within a family of Asp-X-X-quaterthiophene-X-X-Asp (DXX-OT4-XXD) oligopeptides in which X is one of the five amino acids {Ala, Phe, Gly, Ile, Val} ({A, F, G, I, V}). Molecular dynamics simulations reveal that smaller amino acid substituents (A, G) favor linear stacking within a peptide dimer, whereas larger groups (F, I, V) induce larger twist angles between the peptides. Density functional theory calculations on the dimer show the absorption spectrum to be dominated by transitions between carbon and sulfur p orbitals. Although the absorption spectrum is largely insensitive to the relative twist angle, the highest occupied molecular orbital strongly localizes onto one molecule within the dimer at large twist angles, impeding the efficiency of transport between molecules. Our results provide a fundamental understanding of the relation between peptide orientation and electronic structure and offer design precepts for rational engineering of these systems.

Electrochromism of alkylene-linked discrete chromophore polymers with broad radical cation light absorption

This paper details the theoretical calculations, synthesis, and electrochromism of polymeric materials consisting of conjugated chromophores separated by nonconjugated linkers.

Distribution of dopant ions around poly(3,4-ethylenedioxythiophene) chains: a theoretical study

The effect of counterions and multiple polymer chains on the properties and structure of poly(3,4-ethylenedioxythiophene) (PEDOT) doped with ClO₄⁻ has been examined using density functional theory (DFT) calculations with periodic boundary conditions (PBCs).

Steric control in the π-dimerization of oligothiophene radical cations annelated with bicyclo[2.2.2]octene units

A Two series of oligothiophenes 2(nT) (n=4,5), annelated with bicyclo[2.2.2]octene (BCO) units at both ends, and quaterthiophenes 3 a-c, annelated with various numbers of BCO units at different positions, were newly synthesized to investigate the driving forces of π-dimerization and the structure-property relationships of the π-dimers of oligothiophene radical cations. Their radical-cation salts were prepared through chemical one-electron oxidation by using nitrosonium hexafluoroantimonate. From variable-temperature electron spin resonance and electronic absorption measurements, the π-dimerization capability was found to vary among the members of the 2(nT)(+)(·)SbF6(-) series and 3(+)(·)SbF6(-) series of compounds. To examine these results, density functional theory (DFT) calculations at the M06-2X/6-31G(d) level were conducted for the π-dimers. This level of theory was found to successfully reproduce the previously reported X-ray structure of (2(3T))2(2+) having a bent π-dimer structure with cis-cis conformations. The absorption bands obtained by time-dependent DFT calculations for the π-dimers were in reasonable agreement with the experimental spectra. The attractive and repulsive forces for the π-dimerization were divided into four factors: 1) SOMO-SOMO interactions, 2) van der Waals forces, 3) solvation, and 4) Coulomb repulsion, and the effects of each factor on the structural differences and chain-length dependence are discussed in detail.

Stacking interaction and its role in kynurenic acid binding to glutamate ionotropic receptors

Stacking interaction is known to play an important role in protein folding, enzyme-substrate and ligand-receptor complex formation. It has been shown to make a contribution into the aromatic antagonists binding with glutamate ionotropic receptors (iGluRs), in particular, the complex of NMDA receptor NR1 subunit with the kynurenic acid (KYNA) derivatives. The specificity of KYNA binding to the glutamate receptors subtypes might partially result from the differences in stacking interaction. We have calculated the optimal geometry and binding energy of KYNA dimers with the four types of aromatic amino acid residues in Rattus and Drosophila ionotropic iGluR subunits. All ab initio quantum chemical calculations were performed taking into account electron correlations at MP2 and MP4 perturbation theory levels. We have also investigated the potential energy surfaces (PES) of stacking and hydrogen bonds (HBs) within the receptor binding site and calculated the free energy of the ligand-receptor complex formation. The energy of stacking interaction depends both on the size of aromatic moieties and the electrostatic effects. The distribution of charges was shown to determine the geometry of polar aromatic ring dimers. Presumably, stacking interaction is important at the first stage of ligand binding when HBs are weak. The freedom of ligand movements and rotation within receptor site provides the precise tuning of the HBs pattern, while the incorrect stacking binding prohibits the ligand-receptor complex formation.

A unified electrostatic and cavitation model for first-principles molecular dynamics in solution

The electrostatic continuum solvent model developed by [Fattebert and Gygi J. Comput. Chem. 23, 662 (2002); Int. J. Quantum Chem. 93, 139 (2003)] is combined with a first-principles formulation of the cavitation energy based on a natural quantum-mechanical definition for the surface of a solute. Despite its simplicity, the cavitation contribution calculated by this approach is found to be in remarkable agreement with that obtained by more complex algorithms relying on a large set of parameters. Our model allows for very efficient Car-Parrinello simulations of finite or extended systems in solution and demonstrates a level of accuracy as good as that of established quantum-chemistry continuum solvent methods. We apply this approach to the study of tetracyanoethylene dimers in dichloromethane, providing valuable structural and dynamical insights on the dimerization phenomenon.

Theoretical Analysis of Intermolecular Covalent π−π Bonding and Magnetic Properties of Phenalenyl and spiro-Biphenalenyl Radical π-Dimers

Singlet-triplet splittings DeltaEST and intermolecular covalent pi-pi bonding characteristics of the prototypical phenalenyl pi-dimer and eight spiro-biphenalenyl radical pi-dimer structures are analyzed with the aid of restricted and unrestricted density functional theory calculations and paramagnetic susceptibility data fitted using the Bleaney-Bowers dimer model and the Curie-Weiss model. Single determinant approximations for DeltaEST as a function of transfer integrals and on-site Coulomb repulsion energy are presented for the two-electron two-site pi-dimers of phenalenyls and the two-electron four-site pi-dimers of spiro-biphenalenyl radicals. Within the range of intermolecular separation of 3.12<D<3.51 A, for the shorter separations, restricted theory works quite well and indicates the presence of a relatively strong intermolecular covalent pi-pi bonding interaction. For the longer separations, the singlet-triplet splittings are small; electron correlation plays a significant role, and only the unrestricted theory provides results that are in qualitative agreement with experiments. The bonding interactions in the pi-dimers are gradually weakened with increasing D, showing a transition from low D values with significant intermolecular pi-pi bonding and electron delocalization to high D values with localized spins and a biradicaloid character.

Solvent Effects on Electronic Structures and Chain Conformations of α-Oligothiophenes in Polar and Apolar Solutions

Solvent effects on electronic structures and chain conformations of alpha-oligothiophenes nTs (n = 1 to 10) are investigated in solvents of n-hexane, 1,4-dioxane, carbon tetrachloride, chloroform, and water by using density functional theory (DFT) and molecular dynamics (MD) simulations. Both implicit and explicit solvent models are employed. The polarized continuum model (PCM) calculations and MD simulations demonstrate the weak solvent effects on the electronic structures of alpha-oligothiophenes. The lowest dipole-allowed vertical excitation energies of nTs, obtained from time-dependent DFT/PCM calculations at the B3LYP/6-31G(d) level, exhibit a red shift as the solvent polarity increases, in agreement with experiments. The studied solvents have little impact on the state order of the low-lying excited states provided that the nTs are kept in C2h or C2v symmetry. The MD simulations demonstrate that the chain conformations are distorted to some extent in polar and nonpolar solvents. A qualitative picture of the distribution of solvent molecules around the solvated nTs is drawn by means of radial and spatial distribution functions. The S...H-O and pi...H-O solute-solvent interactions are insignificant in aqueous solution.

Intermolecular Covalent π−π Bonding Interaction Indicated by Bond Distances, Energy Bands, and Magnetism in Biphenalenyl Biradicaloid Molecular Crystal

Density-functional theory (DFT) calculations were performed for energy band structure and geometry optimizations on the stepped pi-chain, the isolated molecule and (di)cations of the chain, and various related molecules of a neutral biphenalenyl biradicaloid (BPBR) organic semiconductor 2. The dependence of the geometries on crystal packing provides indirect evidence for the intermolecular covalent pi-pi bonding interaction through space between neighboring pi-stacked phenalenyl units along the chain. The two phenalenyl electrons on each molecule, occupying the singly occupied molecular orbitals (SOMOs), are participating in the intermolecular covalent pi-pi bonding making them partially localized on the phenalenyl units and less available for intramolecular delocalization. The band structure shows a relatively large bandwidth and small band gap indicative of good pi-pi overlap and delocalization between neighboring pi-stacked phenalenyl units. A new interpretation is presented for the magnetism of the stepped pi-chain of 2 using an alternating Heisenberg chain model, which is consistent with DFT total energy calculations for 2 and prevails against the previous interpretation using a Bleaney-Bowers dimer model. The obtained transfer integrals and the magnetic exchange parameters fit well into the framework of a Hubbard model. All presented analyses on molecular geometries, energy bands, and magnetism provide a coherent picture for 2 pointing toward an alternating chain with significant intermolecular through-space covalent pi-pi bonding interactions in the molecular crystal. Surprisingly, both the intermolecular transfer integrals and exchange parameters are larger than the intramolecular through-bond values indicating the effectiveness of the intermolecular overlap of the phenalenyl SOMO electrons.

Photoluminescence and Conductivity of Self-Assembled π–π Stacks of Perylene Bisimide Dyes

The self-assembly of a new, highly fluorescent perylene bisimide dye 2 into pi stacks, both in solution and condensed phase, has been studied in detail by NMR spectroscopy, vapor pressure osmometry (VPO), UV/Vis and fluorescence spectroscopy, differential scanning calorimetry (DSC), optical polarizing microscopy (OPM) and X-ray diffraction. The NMR and VPO measurements revealed the formation of extended pi-pi stacks of the dye molecules in solution. The aggregate size determined from VPO and DOSY NMR measurements agree well with that obtained from the concentration and temperature-dependent UV/Vis spectral data by employing the isodesmic model (equal K model). In the condensed state, dye 2 possesses a hexagonal columnar liquid crystalline (LC) phase as confirmed by X-ray diffraction analysis. The columnar stacking of this dye has been further explored by atomic force microscopy (AFM). Well-resolved columnar nanostructures of the compound are observed on graphite surface. A color-tunable luminescence from green to red has been observed upon aggregation which is accompanied by an increase of the fluorescence lifetime and depolarization. The observed absorption properties can be explained in terms of molecular exciton theory. The charge transport properties of dye 2 have been investigated by pulse radiolysis-time resolved microwave conductivity measurements and a 1D charge carrier mobility up to 0.42 cm(2) V(-1) s(-1) is obtained. Considering the promising self-assembly, semiconducting, and luminescence properties of this dye, it might serve as a useful functional material for nano(opto)electronics.

Effect of Counterions on the Interactions of Charged Oligothiophenes

The functionality of conjugated polymer systems often relies on oxidations or reductions, in most cases mediated by the presence of counterions. The effect that the common counterion hexafluorophosphate (PF6-) has on the intermolecular interactions between charged oligothiophenes is investigated here using ab initio quantum chemistry methods. Counterions are explicitly included in the simulations of oxidized oligothiophenes and in the dimerization process. Our calculations provide quantitative and qualitative insight into the intermolecular interactions in oligothiophene-counterion systems and show that the intermolecular pi-stacking of oligothiophenes is not adversely affected by the presence of counterions and that in fact oligothiophene dimerization is further stabilized by their presence.

Solvent-Induced Effects: Self-Association of Positively Charged π Systems

Antimalarial cationic drugs, such as chloroquine (CQ) and ferroquine (FQ), form stable dimer structures not only in the solid state but also in solution. The short distances (3.3-3.5 A) observed between the positively charged quinolinium rings suggest that this self-association process is driven by pi/pi stacking interactions. Nevertheless, the strength of these dispersive forces is likely not sufficient to overcome the strong repulsive +/+ electrostatic effects. The question of the exact role of the environment, particularly the solvent, is clearly raised here. Characterization of these unconventional stabilizing nonbonding interactions which we have named +-pi/+-pi is therefore of great importance. In the present work, we describe theoretical calculations and NMR experiments undertaken to probe the nature and the strength of +-pi/+-pi interactions occurring upon self-association of FQ and CQ molecules in water.

Theoretical Study of Long Oligothiophene Dications: Bipolaron vs Polaron Pair vs Triplet State

A series of oligothiophene dications (from the sexithiophene dication to the 50-mer oligothiophene dication, nT2+, n = 6-50) were studied. Density functional theory (DFT) at the B3LYP/6-31G(d) level and, in some cases, also at BLYP/6-31Gd, was applied to study the singlet and triplet states of the whole series. We found that the singlet state is the ground state for all oligothiophene dications up to the 20-mer, and that the singlet and triplet states are degenerate for longer oligomers. Thus, the triplet state is never a pure ground state for these dications. We found that, for short oligothiophenes dication (e.g., 6T2+), the bipolaron state is the more important state, with only a small contribution made by the polaron pair state. For medium size oligothiophene dications (e.g., 14T2+), both the bipolaron state and the polaron-pair state contribute to the electronic structure. Finally, in long oligothiophene dications, such as 30T2+ and 50T2+, the contribution from the polaron pair state becomes dominant, and these molecules can be considered as consisting of two independent cation radicals or a polaron pair. Results from isodesmic reactions show that the stability of oligothiophene cation radicals over dications is inversely proportional to chain length. Small oligothiophene dications (n = 6-12) were studied at the CASSCF(m,m)/6-31G(d) (m = 4, 6, and 10) level. The major conclusions of this paper regarding the relative energy of the singlet state versus the triplet state and regarding the relative stability of the bipolaron versus the polaron pair were also supported by CASSCF calculations.

Stacking of oligo- and polythiophene cations in solution: Surface tension and dielectric saturation

The stacking of positively charged (or doped) terthiophene oligomers and quaterthiophene polymers in solution is investigated applying a recently developed unified electrostatic and cavitation model for first-principles calculations in a continuum solvent. The thermodynamic and structural patterns of the dimerization are explored in different solvents, and the distinctive roles of polarity and surface tension are characterized and analyzed. Interestingly, we discover a saturation in the stabilization effect of the dielectric screening that takes place at rather small values of epsilon(0). Moreover, we address the interactions in trimers of terthiophene cations, with the aim of generalizing the results obtained for the dimers to the case of higher-order stacks and nanoaggregates.

Intermolecular Charge Transfer between Heterocyclic Oligomers. Effects of Heteroatom and Molecular Packing on Hopping Transport in Organic Semiconductors

For electron or hole transfer between neighboring conducting polymer strands or oligomers, the intrinsic charge-transfer rate is dictated by the charge-resonance integral and by the reorganization energy due to geometric relaxation. To explain conduction anisotropy and other solid-state effects, a multivariate, systematic analysis of bandwidth as a function of intermolecular orientations is undertaken for a series of oligoheterocycles, using first-principles methods. While cofacial oligomers show the greatest bandwidths at a given intermolecular C-C contact distance, for a fixed center-to-center intermolecular distance, tilted pi-stacking increases pi-overlap (particularly for LUMO orbitals) and decreases electrostatic repulsion, yielding optimum tilt angles for packing of approximately 40-60 degrees at small intermolecular separations. The calculations also reveal that bandwidths and intrinsic mobilities of holes and electrons in conjugated oligoheterocycles can be quite comparable.

π-Stacking in Thiophene Oligomers as the Driving Force for Electroactive Materials and Devices

The pi-stacking between aromatic oligomers has been extensively studied for many years, although the notion of exploiting this phenomenon as the driving force for molecular actuation has only recently emerged. In this work we examine with MP2 and Car-Parrinello ab initio calculations the actuation properties of a novel class of thiophene-based materials introduced by Swager et al. (Adv. Mater. 2002, 14, 368; J. Am. Chem. Soc. 2003, 125, 1142). The chemical ingredients of the assembly, calix[4]arenes and oligothiophenes, are screened separately to characterize the actuation mechanisms and design optimal architectures. In particular, ab initio methods are used to study pi-stacking in mixed-valence oligothiophene dimers, revealing strong interactions that can be turned on and off as a function of the electrochemical potential. We show how these interactions could be harnessed to achieve molecular actuation and investigate the response of an active unit in real time with first-principles molecular dynamics simulations.