Examination of stabilty of molecular agglomerates in molecular crystals

We report an examination algorithm of stability of molecular aggregates based on the estimation of rigidity of intermolecular contacts in a crystal structure. The algorithm includes the intermolecular interaction energy calculation (in the atom-atom potential approximation) of a pair of molecules selected in the crystal structure. Further, the energy is minimized using a least-squares technique by varying the position and orientation of one of the molecules. The contact rigidity is quantitatively assessed by the minimized rms difference between the positions of the atoms in the original and optimized structures (Zorkii’s criterion). Every rigid contact revealed in the structure determines finite or infinite stable agglomerates. The paper presents the results of testing the computer program based on this algorithm with a number of real crystal structures previously determined by single crystal X-ray diffraction, and also the examples of the most common stable molecular agglomerates found with the aid of the program.