This paper bridges the gap between high-level ab initio computations of gas-phase models of 1 : 1 arene–arene complexes and calculations of the two-component (binary) organic crystals using atom–atom potentials. The studied crystals consist of electron-rich and electron-deficient compounds, which form infinite stacks (columns) of heterodimers. The sublimation enthalpy of crystals has been evaluated by DFT periodic calculations, while intermolecular interactions have been characterized by Bader analysis of the periodic electronic density. The consideration of aromatic compounds without a dipole moment makes it possible to reveal the contribution of quadrupole–quadrupole interactions to the π-stacking energy. These interactions are significant for heterodimers formed by arenes with more than 2 rings, with absolute values of the traceless quadrupole moment (Q_zz) larger than 10 D Å. The further aggregation of neighboring stacks is due to the C–H⋯F interactions in arene/perfluoroarene crystals. In crystals consisting of arene and an electron-deficient compound such as pyromellitic dianhydride, aggregation occurs due to the C–H⋯O interactions. The C–H⋯F and C–H⋯O inter-stacking interactions make the main contribution to the sublimation enthalpy, which exceeds 150 kJ mol⁻¹ for the two-component crystals formed by arenes with more than 2 rings.

The influence of pressure on the thermodynamic, structural, and elastic properties of the 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) molecular crystal at T = 0 is systematically studied. Calculations are carried out using density functional theory methods in a plane wave basis set with dispersion corrections for the exchange-correlation part of total energy, and ultrasoft pseudopotentials. The equilibrium unit cell parameters, the cold compression curve in the pressure range of 0–50 GPa and the sound speeds are computed. The effect of finite pressure on the molecular structure of TATB is elucidated from the analysis of relative changes in the intra- and intermolecular geometrical parameters. For the first time, the full set of elastic constants of this crystal at zero and non-zero pressures is determined from ab initio calculations. The resulted structural, elastic, and acoustic properties of TATB are shown to be in a good agreement with available experimental and theoretical data.

Synthesis of multicomponent solid forms is an important method of modifying and fine-tuning the most critical physicochemical properties of drug compounds. The design of new multicomponent pharmaceutical materials requires reliable information about the supramolecular arrangement of molecules and detailed description of the intermolecular interactions in the crystal structure. It implies the use of a combination of different experimental and theoretical investigation methods. Organic salts present new challenges for those who develop theoretical approaches describing the structure, spectral properties, and lattice energy E_latt. These crystals consist of closed-shell organic ions interacting through relatively strong hydrogen bonds, which leads to E_latt > 200 kJ/mol. Some technical problems that a user of periodic (solid-state) density functional theory (DFT) programs encounters when calculating the properties of these crystals still remain unsolved, for example, the influence of cell parameter optimization on the E_latt value, wave numbers, relative intensity of Raman-active vibrations in the low-frequency region, etc. In this work, various properties of a new two-component carbendazim maleate crystal were experimentally investigated, and the applicability of different DFT functionals and empirical Grimme corrections to the description of the obtained structural and spectroscopic properties was tested. Based on this, practical recommendations were developed for further theoretical studies of multicomponent organic pharmaceutical crystals.