Computational Techniques for Predicting Mechanical Properties of Organic Crystals: A Systematic Evaluation

Abstract

Understanding of the structure – mechanical properties relationship in organic crystals can potentially facilitate the design of crystals with desired mechanical properties through crystal engineering. To understand and predict crystal mechanical properties, including tableting behavior, a number of computational methods have been developed to analyze crystal structure. These include visualization, attachment energy calculations, topological analysis, energy framework, and elasticity tensor calculation. However, different methods often lead to conflicting predictions. There is a need for a computational tool kit for predicting crystal mechanical properties from crystal structures. Using α-oxalic acid anhydrous (OAA) and dihydrate (OAD) as a model system, we have systematically compared the predictive accuracy of the experimentally determined mechanical properties using powder compaction and nanoindentation of several methods. We have found that crystal plasticity can be accurately predicted based on energy framework combined with topological analysis and DFT calculated elasticity tensor. Although very useful in characterizing crystal packing features, structure visualization, topology analysis, and attachment energy calculations alone are insufficient for accurately identifying the slip planes and predicting mechanical properties and tableting behavior of organic crystals.