Can we Predict Plastic and Elastic Properties of Monocrystalline Materials? Now we Can.


A crystal is a solid material whose constituent atoms, molecules, or ions are arranged in an orderly repeating pattern extending in all three spatial dimensions. The ability to comprehend and control the mechanical properties of crystalline solids holds great promise for new advancements in crystal engineering and material design. Among other classes, the organic materials offer especially wide range of stress response. For instance, the “organosuperplasticity” effect, defined as superplastic deformation without heating, was found recently in an organic crystal. Another advantage of organic materials is the opportunity for systematic improvement. Recent publications have further demonstrated how density functional theory (DFT) in conjunction with global minimization methods for the lattice energy may serve as powerful tools in crystal structure predictions (CSP). In particular, the latter include evolutionary algorithms that help to avoid computationally expensive exhaustive enumeration of all the local minima and direct the search toward lower energy structures by retaining and combining energetically favorable structural features. However, while CSP methods have been widespread in recent years, developments of computational tools for predicting mechanical properties have received less attention.

Many organic compounds form two or more solid phases that differ only in the packing of molecules in the unit cell. These phases are known as polymorphs and can have dissimilar physical properties. Therefore, there is the need to develop a methodology that can discern various polymorphs. In this regard, University of Central Florida and South Ural State University researchers led by Professors Artem Masunov and Ekaterina Bartashevich developed a novel methodology for computational predictions of the mechanical properties for single crystals. The proposed methodology was based on constrained optimization using dispersion-corrected density functional theory level, and can be dubbed the virtual tensile test. Their work is currently published in the research journal, Crystal Growth & Design.

This approach was tested on 4-bromophenyl 4-bromobenzoate (BPBB), which has been shown to possess different mechanical properties in its three known polymorphs: form I, form II, and form III. Each one of these polymorphic crystal structures was stretched stepwise along each crystallographic axis, while the remaining lattice parameters and atomic coordinates were relaxed. The geometrical properties of halogen bonds and the other noncovalent interactions were monitored at each step to understand the nature of mechanical response. Lastly, the unit cell volumes and lattice energies were plotted as functions of the stretching parameter, and these curves were analyzed in terms of mechanical properties of the brittle, plastic, and elastic polymorphs.

The authors reported that when elastic form I was stretched along the c axis, the local cavities between molecules in the crystal were formed, which explained the increase in cell volume. At the 40% stretch, the cavity volume was seen to reach its maximum. With further increase in strain, the cavities “collapsed” and Br atoms formed new short contacts, resulting in decrease of lattice energy into a new local minimum. For brittle form II, even under small deformations along all cell axes, the lattice energy spiked. The researchers also reported that when the plastic form III was stretched along the c axis, the layers of molecules slipped along each other.

In summary, the study reported on the development of a novel methodology for computational predictions of the mechanical properties for single crystals. The approach described here can be called the virtual tensile test. Remarkably, it seems to be quite universal and applicable to a wide range of crystals. In a statement to Advances in Engineering, Professor Artem Masunov explained: “Combining crystal structure prediction techniques with our virtual tensile test holds a great promise for crystal engineering. We envision no limits on the type and composition of a crystal, since only the stepwise elongation of one lattice parameter followed by its freezing in the relaxation of the altered structure is required in our approach.”

Can we Predict Plastic and Elastic Properties of Monocrystalline Materials? Now we Can. - Advances in Engineering

About the author

Artëm E. Masunov studied chemistry at Moscow State University, where he received his B.S./M.S. degree in 1988. He earned his Ph.D. in 2000 from the The Graduate School, City University of New York. After two years as Postdoctoral Research Associate at City College of New York and three years as Technical Staff Member at Los Alamos National Laboratory, he joined the faculty at University of Central Florida NanoScience Technology Center in 2005, where he presently is an Associate Professor. Dr. Masunov’s group uses Density Functional Theory and Molecular Dynamics approaches for computational materials design and crystal engineering, predictions of nonlinear optical properties, photophysics and photochemistry, and studies in chemical kinetics of combustion.

One of the main research topics of Dr. Masunov’s group is the aggregation and self-assembly phenomena. Aggregation determines organization of many organic materials, and Molecular Dynamics simulations of self-assembly is a promising tool for computational materials engineering. Dr. Masunov discovered counterintuitive blue shift in two-photon absorption spectra of J-aggregates, and effect of polymorphism on exciton localization in conducting polymers. The nature of the conducting polymer was found to affect the competition between the polymer aggregation and polymer wrapping in of carbon nanotube in the composites prepared for organic photovoltaic. Masunov’s group also investigated the role of aggregation inhibitors and mutations on formation and disassembly of amyloid aggregates responsible for Altzheimer’s and Parkinson’s diseases, determined the size of the critical nucleolus, and established the role of dynamic cooperative effects during the amyloid growth. More recently, Dr. Masunov turned his attention to prediction of structure and physical properties of organic crystals.

Dr. Masunov is a 2010 Wiley Young Investigator, and 2008 HP Outstanding Junior Faculty in Computational Chemistry section of the American Chemical Society. His group has attracted funding from the NSF, DOE, NASA, and DOD. During his career, Masunov supervised and graduated seven Ph.D. and five M.S. students. To date, he published 145 papers in peer-reviewed journals. He also serves on the editorial board for a Hindawi’s Journal of Spectroscopy.


Artem E. Masunov, Meryl Wiratmo, Alexander A. Dyakov, Yury V. Matveychuk, Ekaterina V. Bartashevich. Virtual Tensile Test for Brittle, Plastic, and Elastic Polymorphs of 4‑Bromophenyl 4‑Bromobenzoate. Crystal Growth & Design 2020: volume 20, page 6093−6100.

Go To Crystal Growth & Design

Check Also

Method to penetrate through the overdense plasma sheath: Mode conversion in nonuniform plasmas with two-temperature electron population - Advances in Engineering

Method to penetrate through the overdense plasma sheath