Inherent structures of crystalline pentacene

Charge mobility calculation of organic semiconductors without use of experimental single-crystal data

Prediction of material properties of newly designed molecules is a long-term goal in organic electronics. In general, it is a difficult problem, because the material properties are dominated by the unknown packing structure. We present a practical method to obtain charge transport properties of organic single crystals, without use of experimental single-crystal data. As a demonstration, we employ the promising molecule C10-DNBDT. We succeeded in quantitative evaluation of charge mobility of the single crystal using our quantum wave-packet dynamical simulation method. Here, the single-crystal data is computationally obtained by searching possible packing structures from structural formula of the molecule. We increase accuracy in identifying the actual crystal structure from suggested ones by using not only crystal energy but also similarity between calculated and experimental powder X-ray diffraction patterns. The proposed methodology can be a theoretical design technique for efficiently developing new high-performance organic semiconductors, since it can estimate the charge transport properties at early stage in the process of material development.

Computational modelling of solvent effects in a prolific solvatomorphic porous organic cage

Crystal structure prediction methods can enable the in silico design of functional molecular crystals, but solvent effects can have a major influence on relative lattice energies, sometimes thwarting predictions. This is particularly true for porous solids, where solvent included in the pores can have an important energetic contribution. We present a Monte Carlo solvent insertion procedure for predicting the solvent filling of porous structures from crystal structure prediction landscapes, tested using a highly solvatomorphic porous organic cage molecule, CC1. Using this method, we can understand why the predicted global energy minimum structure for CC1 is never observed from solvent crystallisation. We also explain the formation of three different solvatomorphs of CC1 from three structurally-similar chlorinated solvents. Calculated solvent stabilisation energies are found to correlate with experimental results from thermogravimetric analysis, suggesting a future computational framework for a priori materials design that factors in solvation effects.

Energy-Structure-Function Maps: Cartography for Materials Discovery

Some of the most successful approaches to structural design in materials chemistry have exploited strong directional bonds, whose geometric reliability lends predictability to solid-state assembly. For example, metal-organic frameworks are an important design platform in materials chemistry. By contrast, the structure of molecular crystals is defined by a balance of weaker intermolecular forces, and small changes to the molecular building blocks can lead to large changes in crystal packing. Hence, empirical rules are inherently less reliable for engineering the structures of molecular solids. Energy-structure-function (ESF) maps are a new approach for the discovery of functional organic crystals. These maps fuse crystal-structure prediction with the computation of physical properties to allow researchers to choose the most promising molecule for a given application, prior to its synthesis. ESF maps were used recently to discover a highly porous molecular crystal that has a high methane deliverable capacity and the lowest density molecular crystal reported to date (r = 0.41 g cm^-3 , SA_BET = 3425 m² g^-1 ). Progress in this field is reviewed, with emphasis on the future opportunities and challenges for a design strategy based on computed ESF maps.

Crystal structure prediction of rigid molecules

A non-polarizable force field based on atomic multipoles fit to reproduce experimental crystal properties andab initiogas-phase dimers is described. The Ewald method is used to calculate both long-range electrostatic and 1/r6dispersion energies of crystals. The dispersion energy of a crystal calculated by a cutoff method is shown to converge slowly to the exact Ewald result. A method for constraining space-group symmetry during unit-cell optimization is derived. Results for locally optimizing 4427 unit cells including volume, cell parameters, unit-cell r.m.s.d. and CPU timings are given for both flexible and rigid molecule optimization. An algorithm for randomly generating rigid molecule crystals is described. Using the correct experimentally determined space group, the average and maximum number of random crystals needed to find the correct experimental structure is given for 2440 rigid single component crystals. The force field energy rank of the correct experimental structure is presented for the same set of 2440 rigid single component crystals assuming the correct space group. A complete crystal prediction is performed for two rigid molecules by searching over the 32 most probable space groups.

Convergence Properties of Crystal Structure Prediction by Quasi-Random Sampling

Generating sets of trial structures that sample the configurational space of crystal packing possibilities is an essential step in the process of ab initio crystal structure prediction (CSP). One effective methodology for performing such a search relies on low-discrepancy, quasi-random sampling, and our implementation of such a search for molecular crystals is described in this paper. Herein we restrict ourselves to rigid organic molecules and, by considering their geometric properties, build trial crystal packings as starting points for local lattice energy minimization. We also describe a method to match instances of the same structure, which we use to measure the convergence of our packing search toward completeness. The use of these tools is demonstrated for a set of molecules with diverse molecular characteristics and as representative of areas of application where CSP has been applied. An important finding is that the lowest energy crystal structures are typically located early and frequently during a quasi-random search of phase space. It is usually the complete sampling of higher energy structures that requires extended sampling. We show how the procedure can first be refined, through targetting the volume of the generated crystal structures, and then extended across a range of space groups to make a full CSP search and locate experimentally observed and lists of hypothetical polymorphs. As the described method has also been created to lie at the base of more involved approaches to CSP, which are being developed within the Global Lattice Energy Explorer (Glee) software, a few of these extensions are briefly discussed.

Epitaxial growth of π-stacked perfluoropentacene on graphene-coated quartz

Chemical-vapor-deposited large-area graphene is employed as the coating of transparent substrates for the growth of the prototypical organic n-type semiconductor perfluoropentacene (PFP). The graphene coating is found to cause face-on growth of PFP in a yet unknown substrate-mediated polymorph, which is solved by combining grazing-incidence X-ray diffraction with theoretical structure modeling. In contrast to the otherwise common herringbone arrangement of PFP in single crystals and "standing" films, we report a π-stacked arrangement of coplanar molecules in "flat-lying" films, which exhibit an exceedingly low π-stacking distance of only 3.07 Å, giving rise to significant electronic band dispersion along the π-stacking direction, as evidenced by ultraviolet photoelectron spectroscopy. Our study underlines the high potential of graphene for use as a transparent electrode in (opto-)electronic applications, where optimized vertical transport through flat-lying conjugated organic molecules is desired.

Are crystal polymorphs predictable? The case of sexithiophene

Using sexithiophene as a benchmark compound, we present a very effective strategy for searching the potential energy minima of a crystalline material, described in terms of rigid molecules with Coulombic and atom-atom interactions. The strategy involves uniform sampling of the many-body energy hypersurface, mechanical identification of all constraints deriving from the crystallographic symmetry, and a "sight-resight" method, originally introduced in wildlife ecology, for assessing the completeness of the search. Thousands of distinct potential energy minima, with a surprising variety of structural arrangements, are identified for sexithiophene. Despite the large number of competing minima, the system presents a small number of deep minima, with very different structures and not particularly congested in energy or density. The two deepest minima correspond to the structures of the two known experimental polymorphs, which are satisfactorily described.

Inherent structures of crystalline tetracene

We have systematically sampled the potential energy surface of crystalline tetracene to identify its local minima. These minima represent all possible stable configurations and constitute the "inherent structures" of the system. The crystal is described in terms of rigid molecules with Coulombic and atom-atom interactions. Hundreds of distinct minima are identified, mostly belonging to the space groups P (triclinic) and P2(1)/c (monoclinic), with a variety of structural arrangements. The deepest minimum corresponds to the high temperature-low pressure polymorph. This is the only polymorph with a completely described X-ray structure, which is satisfactorily described by the calculations. The next deep minimum is likely to correspond to the low temperature-high pressure polymorph, which has been experimentally identified but not yet fully described.

Microscopic Description of Elementary Growth Processes and Classification of Structural Defects in Pentacene Thin Films

Elementary growth processes such as kink initiation, adding a molecule to a kink, and adding a molecule between two neighboring kinks and between two grains are theoretically studied in pentacene films by adding one molecule at a time to a predefined aggregate. For each molecule, the potential energy surface is calculated using the MM3 molecular mechanics force field, which allowed one to identify useful parameters like the energy barrier for diffusion and the energy to create kinks, as well as defect configurations. Depending on the properties of the potential energy surface and the resulting growth-condition-dependent probabilities of initiating defect configurations in the film, three types of pentacene defects are identified: a thermally activated defect, an intrinsic defect, and a kinetic defect. Upon film growth, most defects relax into the ideal crystal configuration. Bulk defects that resist relaxation have densities lower than 10(16) defects/cm3 at typical growth conditions. Grain boundary defects, on the other hand, are very stable. Moreover, interstitial molecules at grain boundaries are identified as a source of compressive stress.

Dynamics of bimodal growth in pentacene thin films

Previous studies have established that pentacene films deposited on silicon oxide consist of a substrate-induced "thin-film" phase, with the bulk phase of pentacene detected in thicker films only. We show that the bulk phase nucleates as early as the first monolayer, and continues to nucleate as film growth progresses, shadowing the growth of the thin-film phase. Moreover, we find that the transition between the "thin-film" and the bulk phase is not a continuous one, as observed in heteroepitaxial systems, but rather the two phases nucleate and grow independently.

Band Structure of the Four Pentacene Polymorphs and Effect on the Hole Mobility at Low Temperature

The band structure of the four known polymorphs of pentacene is computed from first principles using the accurate molecular orbitals of the isolated molecule as the basis for the calculation of the crystalline orbitals. The computed bands are remarkably different for each polymorph, but their diversity can be easily rationalized using a simple analytical model that employs only three parameters. The effect of the electronic structure on the hole mobility was evaluated using a simple model based on the constant relaxation time approximation. It is found that the mobility tensor is highly anisotropic for three of the four considered polymorphs. The practical implication of this prediction on the technology of thin-film organic transistors is discussed.

Inter-residue and solvent-residue interactions in proteins: a statistical study on experimental structures

A large set of protein structures resolved by X-ray or NMR techniques has been extracted from the Protein Data Bank and analyzed using statistical methods. In particular, we investigate the interactions between side chains and the interactions between solvent and side chains, pointing out on the possibility of including the solvent as part of a knowledge-based potential. The solvent-residue contacts are accounted for on the basis of the Voronoi's polyhedron analysis. Our investigation confirms the importance of hydrophobic residues in determining the protein stability. We observe that in general hydrophobic-hydrophobic interactions and, more specifically, aromatic-aromatic contacts tend to be increasingly distally separated in the primary sequence of proteins, thus connecting distinct secondary structure elements. A simple relation expressing the dependence of the protein free energy by the number of residues is proposed. Such a relation includes both the residue-residue and the solvent-residue contributions. The former is dominant for large size proteins, whereas for small sizes (number of residues less than 100) the two terms are comparable. Gapless threading experiments show that the solvent-residue knowledge-based potential yields a significant contribution with respect to discriminating the native structure of proteins. Such contribution is important especially for proteins of small size and is similar to that given by the most favorable residue-residue knowledge-based potential referring to hydrophobic-hydrophobic interactions such as isoleucine-leucine. In general, the inclusion of the solvent-residue interaction produces a relevant increase of the free energy gap between the native structures and decoys.