Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

Miscibility-Function Relations in Organic Solar Cells: Significance of Optimal Miscibility in Relation to Percolation

Control over the morphology of the active layer of bulk heterojunction (BHJ) organic solar cells is paramount to achieve high-efficiency devices. However, no method currently available can predict morphologies for a novel donor–acceptor blend. An approach which allows reaching relevant length scales, retaining chemical specificity, and mimicking experimental fabrication conditions, and which is suited for high-throughput schemes has been proven challenging to find. Here, we propose a method to generate atom-resolved morphologies of BHJs which conforms to these requirements. Coarse-grain (CG) molecular dynamics simulations are employed to simulate the large-scale morphological organization during solution-processing. The use of CG models which retain chemical specificity translates into a direct path to the rational design of donor and acceptor compounds which differ only slightly in chemical nature. Finally, the direct retrieval of fully atomistic detail is possible through backmapping, opening the way fo...

https://pubs.acs.org/doi/pdf/10.1021/jacs.6b11717

Bulk Heterojunction Morphologies with Atomistic Resolution from Coarse-Grain Solvent Evaporation Simulations

Accurate single-crystal X-ray diffraction data offer a unique opportunity to compare and contrast the atomistic details of bulk heterojunction photovoltaic small-molecule acceptor structure and packing, as well as provide an essential starting point for computational electronic structure and charge transport analysis. Herein, we report diffraction-derived crystal structures and computational analyses on the n-type semiconductors which enable some of the highest efficiency organic solar cells produced to date, 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-indanone))-5,5,11,11-tetrakis(4-hexylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (ITIC) and seven derivatives (including three new crystal structures: 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-indanone))-5,5,11,11-tetrakis(4-propylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (ITIC-C3), 3,9-bis(2-methylene-(3-(1,1-dicyanomethylene)-indanone))-5,5,11,11-tetrakis(3-hexylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (m-ITIC-C6), and 3,9-bis(2-methylene-((3-(1,1-dicyanomethylene)-6,7-difluoro)-indanone))-5,5,11,11-tetrakis(4-butylphenyl)-dithieno[2,3-d:2',3'-d']-s-indaceno[1,2-b:5,6-b']dithiophene (ITIC-C4-4F). IDTT acceptors typically pack in a face-to-face fashion with π-π distances ranging from 3.28-3.95 A. Additionally, edge-to-face packing is observed with S⋯π interactions as short as 3.21-3.24 A. Moreover, ITIC end group identities and side chain substituents influence the nature and strength of noncovalent interactions (e. g. H-bonding, π-π) and thus correlate with the observed packing motif, electronic structure, and charge transport properties of the crystals. Density functional theory (DFT) calculations reveal relatively large nearest-neighbor intermolecular π-π electronic couplings (5.85-56.8 meV) and correlate the nature of the band structure with the dispersion interactions in the single crystals and core-end group polarization effects. Overall, this combined experimental and theoretical work reveals key insights into crystal engineering strategies for indacenodithienothiophene (IDTT) acceptors, as well as general design rules for high-efficiency post-fullerene small molecule acceptors.

Building Blocks for High-Efficiency Organic Photovoltaics: Interplay of Molecular, Crystal, and Electronic Properties in Post-Fullerene ITIC Ensembles.

Substantial enhancements in the efficiencies of bulk-heterojunction (BHJ) organic solar cells (OSCs) have come from largely trial-and-error-based optimizations of the morphology of the active layers. Further improvements, however, require a detailed understanding of the relationships among chemical structure, morphology, electronic properties, and device performance. On the experimental side, characterization of the local (i.e., nanoscale) morphology remains challenging, which has called for the development of robust computational methodologies that can reliably address those aspects. In this review, we describe how a methodology that combines all-atom molecular dynamics (AA-MD) simulations with density functional theory (DFT) calculations allows the establishment of chemical structure–local morphology–electronic properties relationships. We also provide a brief overview of coarse-graining methods in an effort to bridge local to global (i.e., mesoscale to microscale) morphology. Finally, we give a few examples of machine learning (ML) applications that can assist in the discovery of these relationships.

Organic Photovoltaics: Relating Chemical Structure, Local Morphology, and Electronic Properties

Understanding how additives interact and segregate within bulk heterojunction (BHJ) thin films is critical for exercising control over structure at multiple length scales and delivering improvements in photovoltaic performance. The morphological evolution of poly(3-hexylthiophene) (P3HT) and phenyl-C61-butyric acid methyl ester (PCBM) blends that are commensurate with the size of a BHJ thin film is examined using petascale coarse-grained molecular dynamics simulations. Comparisons between two-component and three-component systems containing short P3HT chains as additives undergoing thermal annealing demonstrate that the short chains alter the morphology in apparently useful ways: they efficiently migrate to the P3HT/PCBM interface, increasing the P3HT domain size and interfacial area. Simulation results agree with depth profiles determined from neutron reflectometry measurements that reveal PCBM enrichment near substrate and air interfaces but a decrease in that PCBM enrichment when a small amount of shor...

Jan-Michael Y. Carrillo

Papers

Petascale Simulations of the Morphology and the Molecular Interface of Bulk Heterojunctions