Investigation of topology of intermolecular interactions in the benzene–acetylene co-crystal by different theoretical methods

New Polymorphic Modifications of 6-Methyluracil: An Experimental and Quantum Chemical Study

Polymorphism of 6-methyluracil, which affects the regulation of lipid peroxidation and wound healing, has been studied by experimental and quantum chemical methods. Two known polymorphic modifications and two new crystalline forms were crystallized and characterized by single crystal and powder X-ray diffraction (XRD) methods as well as by the differential scanning calorimetry (DSC) method and infrared (IR) spectroscopy. The calculations of pairwise interaction energies between molecules and lattice energies in periodic boundary conditions have shown that the polymorphic form 6MU_I used in the pharmaceutical industry and two new forms 6MU_III and 6MU_IV, which can be formed due to temperature violations, may be considered as metastable. The centrosymmetric dimer bound by two N-H···O hydrogen bonds was recognized as a dimeric building unit in all of the polymorphic forms of 6-methyluracil. Four polymorphic forms have a layered structure from the viewpoint of interaction energies between dimeric building units. The layers parallel to the (100) crystallographic plane were recognized as a basic structural motif in the 6MU_I, 6MU_III, and 6MU_IV crystals. In the 6MU_II structure, a basic structural motif is a layer parallel to the (001) crystallographic plane. The ratio between the interaction energies within the basic structural motif and between neighboring layers correlates with the relative stability of the studied polymorphic forms. The most stable polymorphic form 6MU_II has the most anisotropic "energetic" structure, while the interaction energies in the least stable form 6MU_IV are very close in various directions. The modeling of shear deformations of layers in the metastable polymorphic structures has not revealed any possibility of these crystals to be deformed under external mechanical stress or pressure influence. These results allow the use of metastable polymorphic forms of 6-methyluracil in the pharmaceutical industry without any limitations.

Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods

Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol^-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.

Quantum Chemical Study on Mefenamic Acid Polymorphic Forms

Three polymorphic structures of mefenamic acid, which is a very popular drug, have been studied using quantum chemical methods. It has been shown that the centrosymmetric dimer formed due to two O-H···O hydrogen bonds is a complex building unit in all of the polymorphic structures under study. On the basis of an analysis of the pairwise interaction energies between molecules, the polymorphic forms I and II are classified as columnar-layered while the polymorphic form III has a columnar structure. The stabilities of the three polymorphic forms of mefenamic acid under ambient conditions (I > II > III) correlate with the degree of anisotropy of the interaction energies between columns (primary basic structural motifs) formed due to stacking interactions. The shear deformation modeling of strongly bound layers in all of the polymorphic structures has not revealed any possibility for deformation of the crystal structure. The construction of the shift energy profiles and calculation of the energy barriers for the displacement along the (100) crystallographic plane in the [100], [010], and [011] crystallographic directions make it possible to explain the experimental data obtained for commercially available polymorphic structure I in a diamond anvil cell. The absence of any local minimum near the starting point on the shift energy profile and the extremely high energy barrier can be considered as criteria for the impossibility of a crystal structure deformation under pressure.

Is it possible to predict the stability of a crystal structure under the influence of pressure? Quantum chemical study of ibuprofen crystals

Quantum chemical modeling was used to analyze the crystalline structure of ibuprofen under atmospheric pressure to determine the structural features, providing its stability under pressure.

Usage of Quantum Chemical Methods to Understand the Formation of Concomitant Polymorphs of Acetyl 2-(N-(2-Fluorophenyl)imino)coumarin-3-carboxamide

Crystallization of concomitant polymorphs is a very intriguing process that is difficult to be studied experimentally. A comprehensive study of two polymorphic modifications of acetyl 2-(N-(2-fluorophenyl)imino)coumarin-3-carboxamide using quantum chemical methods has revealed molecular and crystal structure dependence on crystallization conditions. Fast crystallization associated with a kinetically controlled process results in the formation of a columnar structure with a nonequilibrium molecular conformation and more isotropic topology of interaction energies between molecules. Slow crystallization may be considered as a thermodynamically controlled process and leads to the formation of a layered crystal structure with the conformation of the molecule corresponding to local minima and anisotropic topology of interaction energies. Fast crystallization results in the formation of a lot of weak intermolecular interactions, while slow crystallization leads to the formation of small amounts of stronger interactions.

Conformational polymorphs of 3-cyclopropyl-5-(3-methyl-[1,2,4]triazolo[4,3-a]pyridin-7-yl)-1,2,4-oxadiazole

The dipharmacophore compound 3-cyclopropyl-5-(3-methyl-[1,2,4]triazolo[4,3-a]pyridin-7-yl)-1,2,4-oxadiazole, C₁₂H₁₁N₅O, was studied on the assumption of its potential biological activity. Two polymorphic forms differ in both their molecular and crystal structures. The monoclinic polymorphic form was crystallized from more volatile solvents and contains a conformer with a higher relative energy. The basic molecule forms an abundance of interactions with relatively close energies. The orthorhombic polymorph was crystallized very slowly from isoamyl alcohol and contains a conformer with a much lower energy. The basic molecule forms two strong interactions and a large number of weak interactions. Stacking interactions of the `head-to-head' type in the monoclinic structure and of the `head-to-tail' type in the orthorhombic structure proved to be the strongest and form stacked columns in the two polymorphs. The main structural motif of the monoclinic structure is a double column where two stacked columns interact through weak C-H...N hydrogen bonds and dispersive interactions. In the orthorhombic structure, a single stacked column is the main structural motif. Periodic calculations confirmed that the orthorhombic structure obtained by slow evaporation has a lower lattice energy (0.97 kcal mol^-1) compared to the monoclinic structure.

Revisiting 2-chloro-4-nitroaniline: analysis of intricate supramolecular ordering of a triclinic polymorph featuring a high Z value and strong second harmonic generation

A new polymorph of 2-chloro-4-nitroaniline is more stable than its 55 year-old antecedent.

CrystaLattE: Automated computation of lattice energies of organic crystals exploiting the many-body expansion to achieve dual-level parallelism

We present an algorithm to compute the lattice energies of molecular crystals based on the many-body cluster expansion. The required computations on dimers, trimers, etc., within the crystal are independent of each other, leading to a naturally parallel approach. The algorithm exploits the long-range three-dimensional periodic order of crystals to automatically detect and avoid redundant or unnecessary computations. For this purpose, Coulomb-matrix descriptors from machine learning applications are found to be efficient in determining whether two N-mers are identical. The algorithm is implemented as an open-source Python program, CrystaLattE, that uses some of the features of the Quantum Chemistry Common Driver and Databases library. CrystaLattE is initially interfaced with the quantum chemistry package Psi4. With CrystaLattE, we have applied the fast, dispersion-corrected Hartree-Fock method HF-3c to the lattice energy of crystalline benzene. Including all 73 symmetry-unique dimers and 7130 symmetry-unique trimers that can be formed from molecules within a 15 Å cutoff from a central reference monomer, HF-3c plus an Axilrod-Teller-Muto estimate of three-body dispersion exhibits an error of only -1.0 kJ mol^-1 vs the estimated 0 K experimental lattice energy of -55.3 ± 2.2 kJ mol^-1. The convergence of the HF-3c two- and three-body contributions to the lattice energy as a function of intermonomer distance is examined.

Influence of ortho-substituent on the molecular and crystal structures of 2-(N-arylimino)coumarin-3-carboxamide: isotypic and polymorphic structures

During a comprehensive study of a series of 2-(N-arylimino)coumarin-3-carboxamides with the aryl group substituted in the ortho-position by either a halogen atom, a methyl group or a methoxy group, the existence of three groups of isotypic crystal structures has been revealed. The similarity of crystal structures belonging to the same groups was confirmed by the analysis based on the comparison of pairwise interactions energies obtained from quantum chemical calculations. Group I includes unsubstituted, methyl-substituted and polymorphic modification 1 of fluoro-substituted 2-(N-arylimino)coumarin-3-carboxamide. Structures of polymorphic modification 2 of fluoro-substituted derivative, chloro-substituted and polymorphic modification 1 of bromo-substituted 2-(N-arylimino)coumarin-3-carboxamide may represent group II. Group III contains structures of polymorphic modification 2 of bromo-substituted derivative, iodine- and methoxy-substituted 2-(N-arylimino)coumarin-3-carboxamides. Structures of the same type group have extremely close parameters of the unit cell as well as those of molecular and crystal structures. But they are not identical. Polymorphic modifications of fluoro- and bromo-substituted 2-(N-arylimino)coumarin-3-carboxamides belong to different crystal types mainly due to different arrangement of basic structural motifs separated out using quantum chemical calculations.

The formation of the salt and neutral molecule cocrystal from equimolar solution of heliamine and bicyclo[2.2.1]hept-5-ene-endo-2,3-dicarboxylic acid

The possible interaction of 6,7-dimethoxy-1,2,3,4-tetrahydroisoquinoline (heliamine) with bicyclo[2.2.1]hept-5-ene-endo-2,3-dicarboxylic acid anhydride has been studied. Instead of the reaction with heliamine, the acid anhydride was hydrolyzed into the appropriate dicarboxylic acid. An equimolar mixture of unreacted heliamine and in-situ-generated dicarboxylic acid crystallized in space group P2₁/c. The comprehensive study of the obtained crystals shows that the peculiarities of the crystallization process lead to the formation of the salt-cocrystal structure where the dianion interacts simultaneously with two cations forming a chain as the primary structural motif. The neutral molecules of dicarboxylic acid link the dianions of the neighbouring chains, forming a layer as the secondary structural motif. As a result, the stronger hydrogen bonds formed by the neutral molecules play a secondary role in the crystal structure formation.

Intermolecular interactions in crystals of benzene and its mono- and dinitro derivatives: study from the energetic viewpoint

The weak intermolecular interactions and their role in mono- and dinitrobenzene crystal structure formation have been studied using quantum-chemical calculations.

Hydrogen bonding vs. stacking interaction in the crystals of the simplest coumarin derivatives: a study from the energetic viewpoint

Hydrogen bonded and stacked dimers have very close interaction energies in the crystals of the simplest coumarin derivatives according to the data of the crystal structure analysis based on the comparison of pairwise interaction energies.

Competition between intermolecular hydrogen bonding and stacking in the crystals of 4-Hydroxy-N-(pyridin-2-yl)-2,2-dioxo-1H-2λ6,1-benzothiazine-3-carboxamides

The traditional crystal packing study based on the analysis of geometrical characteristics of intermolecular interactions is found to be not informative enough. The application of quantum-chemical calculations for the evaluation of pairwise interaction energies between molecules allows to get much more information about supramolecular architecture. The staking interactions between π-systems of neighboring molecules form the building unit of the crystal packing in the absence of any strong interactions as well as in the presence of the N–H…O classical hydrogen bond. Unexpectedly the hydrogen bonds play the secondary role in the crystal packing formation as compared to stacking. Analysis of the total interaction energy of the basic molecule with all the molecules of its first coordination sphere and the pairwise interaction energies allows to evaluate the extent of crystal isotropy from point of view of interaction energies.

Acceptor properties of amino groups in aminobenzene crystals: study from the energetic viewpoint

The role of the N–H⋯N hydrogen bonds in the organization of the crystals of the aniline and diaminobenzenes has been studied.

Influence of substituents on the acceptor properties of the amino groups in the diaminobenzene analogues

The influence of substituents on the geometric parameters and acceptor properties of the amino group in the diaminobenzene analogues has been studied using quantum-chemical calculations and the study of crystal packing from the energetic viewpoint.