Empirical potential for molecular simulation of graphene nanoplatelets.

A new empirical potential for layered graphitic materials is reported. Interatomic interactions within a single graphene sheet are modeled using a Stillinger-Weber potential. Interatomic interactions between atoms in different sheets of graphene in the nanoplatelet are modeled using a Lennard-Jones interaction potential. The potential is validated by comparing molecular dynamics simulations of tensile deformation with the reported elastic constants for graphite. The graphite is found to fracture into graphene nanoplatelets when subjected to ∼15% tensile strain normal to the basal surface of the graphene stack, with an ultimate stress of 2.0 GPa and toughness of 0.33 GPa. This force field is useful to model molecular interactions in an important class of composite systems comprising 2D materials like graphene and multi-layer graphene nanoplatelets.

Aggregation and Solubility of a Model Conjugated Donor-Acceptor Polymer.

In conjugated polymers, solution-phase structure and aggregation exert a strong influence on device morphology and performance, making understanding solubility crucial for rational design. Using atomistic molecular dynamics (MD) and free-energy sampling algorithms, we examine the aggregation and solubility of the polymer PTB7, studying how side-chain structure can be modified to control aggregation. We demonstrate that free-energy sampling can be used to effectively screen polymer solubility in a variety of solvents but that solubility parameters derived from MD are not predictive. We then study the aggregation of variants of PTB7 including those with linear (octyl), branched (2-ethylhexyl), and cleaved (methyl) side chains, in a selection of explicit solvents and additives. Energetic analysis demonstrates that while side chains do disrupt polymer backbone stacking, solvent exclusion is a critical factor controlling polymer solubility.

Heterogeneous Nucleation of an n-Alkane on Tetrahedrally Coordinated Crystals.

Heterogeneous nucleation refers to the propensity for phase transformations to initiate preferentially on foreign surfaces, such as vessel walls, dust particles, or formulation additives. In crystallization, the form of the initial nucleus has ramifications for the crystallographic form, morphology, and properties of the resulting solid. Nevertheless, the discovery and design of nucleating agents remains a matter of trial and error because of the very small spatiotemporal scales over which the critical nucleus is formed and the extreme difficulty of examining such events empirically. Using molecular dynamics simulations, we demonstrate a method for the rapid screening of entire families of materials for activity as nucleating agents and for characterizing their mechanism of action. The method is applied to the crystallization of n-pentacontane, a model surrogate for polyethylene, on the family of tetrahedrally coordinated crystals, including diamond and silicon. A systematic variation of parameters in the interaction potential permits a comprehensive, physically based screening of nucleating agents in this class of materials, including both real and hypothetical candidates. The induction time for heterogeneous nucleation is shown to depend strongly on crystallographic registry between the nucleating agent and the critical nucleus, indicative of an epitaxial mechanism in this class of materials. Importantly, the severity of this registry requirement weakens with decreasing rigidity of the substrate and increasing strength of attraction to the surface of the nucleating agent. Employing this method, a high-throughput computational screening of nucleating agents becomes possible, facilitating the discovery of novel nucleating agents within a broad "materials genome" of possible additives.