Prediction of the Mechanical Deformation Properties of Organic Crystals Based upon their Crystallographic Structures: Case Studies of Pentaerythritol and Pentaerythritol Tetranitrate.

Abstract

Development of a quantitative model and associated workflow for predicting the mechanical deformation properties (plastic deformation or cleavage fracture) of organic single crystals from their crystallographic structures using molecular and crystallographic modelling. Intermolecular synthons, hydrogen bonding, crystal morphology and surface chemistry are modelled using empirical force fields with the data integrated into the analysis of lattice deformation as computed using a statistical approach. The approach developed comprises three main components. Firstly, the identification of the likely direction of deformation based on lattice unit cell geometry; secondly, the identification of likely lattice planes for deformation through the calculation of the strength and stereochemistry of interplanar intermolecular interactions, surface plane rugosity and surface energy; thirdly, identification of potential crystal planes for cleavage fracture by assessing intermolecular bonding anisotropy. Pentaerythritol is predicted to fracture by brittle cleavage on the {001} lattice planes by strong in-plane hydrogen-bond interactions in the <110>, whereas pentaerythritol tetranitrate is predicted to deform by plastic deformation through the slip system {110} < 001>, with both predictions being in excellent agreement with known experimental data. A crystallographic framework and associated workflow for predicting the mechanical deformation of molecular crystals is developed through quantitative assessment of lattice energetics, crystal surface chemistry and crystal defects. The potential for the de novo prediction of the mechanical deformation of pharmaceutical materials using this approach is highlighted for its potential importance in the design of formulated drug products process as needed for manufacture by direct compression.