Two-component molecular crystals: evaluation of the formation thermodynamics based on melting points and sublimation data †

Based on literature analysis, we have built up a database containing the fusion temperatures of two-component molecular co-crystals and individual compounds (1175 co-crystals/salts). In order to estimate the thermodynamics of two-component crystal formation, we have created a database on the basis of values reported in the literature from 1900 till 2016 inclusive. The database includes values of the en-thalpies and Gibbs energies of individual molecular crystals obtained by various methods. The distribution functions of two-component crystals have been analysed using their fusion temperatures, both for the full sample and separately for the salts and the co-crystals. A comparative analysis was conducted to determine the differences in the melting temperatures of monotropic polymorphic forms, as well as a similar analysis of two-component crystals with the same composition and different stoichiometries. Correlation equations have been obtained, connecting the melting points of co-crystals/salts and individual components for 74 active pharmaceutical ingredients and coformers, which enabled us to design co-crystals with predictable melting temperatures. An approach to estimating co-crystal sublimation thermodynamic characteristics has been developed. The thermodynamic functions of the formation process of 281 co-crystals have been obtained and analyzed. The diagram method has been used to analyze the parameters under study. Analysis of experimental data distribution in the diagram sectors has shown that the number of two-component crystals with enthalpy determined processes of co-crystal/salt formation corresponds to 70.9%, whereas the number of those with entropy determined processes corresponds to 29.1%. A general algorithm for estimating the thermodynamics of the formation of two-component crystals is proposed.


Introduction
Co-crystallization has become an important research area in recent years because of the great potential of fine-tuning the physical properties of multi-component crystal components. The term "compound" has long been used to denote cocrystals with a fixed component ratio, while those with variable stoichiometry were named "mixed crystals". 1,2 Despite numerous studies, the nature of co-crystallization remains inadequately understood. Moreover, understanding and controlling the main driving forces of co-crystal formation is a significant fundamental task. In the past decades, considerable efforts have been made to apply co-crystallization techniques to improving factors such as dissolution rate, thermal stabil-ity, and bioavailability in the pharmaceutical industry. [3][4][5][6] Encouraged by the great success of pharmaceutical co-crystallization, the development of new functional materials by cocrystallization techniques in other fields is attracting more and more interest. Nowadays, energetic co-crystallization becomes very promising as an alternative method of tuning explosive materials and properties by rearranging existing molecules, instead of synthesizing new compounds. 7,8 Another application of co-crystals is in organic nonlinear optical materials with high second-order nonlinear susceptibilities. These materials are used in electro-optical devices such as optical waveguides or frequency modulators. 9,10 Obtaining co-crystals is a spontaneous, unpredictable process that determines the existence of different screening algorithms. In turn, the need for screening procedures makes the co-crystal technology itself expensive and time-consuming. All this makes it necessary to develop algorithms for predicting thermodynamically stable co-crystals/salts (from the available set of individual substances), in order to minimize the number of screening operations and to know in advance the thermodynamic stability of the new materials.
There is a number of works analyzing the co-crystal formation enthalpy. One of the approaches relies on crystal structure prediction (CSP) using anisotropic potential 11,12 and quick methods of energy estimation based on molecular electrostatic potential surfaces. 13 Abramov et al. 14 applied COSMO-RS fluid-phase thermodynamics computations describing the miscibility of co-crystal formers in a supercooled liquid (melt) phase to virtual coformer screening. Moreover, Hansen solubility parameters were recently used to describe the miscibility of APIs and coformers to predict co-crystal formation in order to guide co-crystal screening. 15 An essential fault of all these approaches is that they analyze the enthalpic characteristics of co-crystal formation, while the entropic terms are not taken into account. The existence of such approaches can be attributed to the fact that there are almost no experimental data about the Gibbs energies of co-crystal formation.
Development of approaches for estimating the thermodynamic characteristics of multi-component crystal formation requires experimental data. Unfortunately, there is a limited number of experimental works in this field so far, and they can be conditionally divided into several groups. The most informative works are those containing data on all the thermodynamic functions of the co-crystal/salt formation (Gibbs energy, enthalpy and entropy). As a rule, such information is obtained by co-crystal equilibrium solubility experiments conducted at different temperatures. [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31] However, there are works where two independent techniquessolubility calorimetry (to analyze the enthalpic term) and solubility measured by the saturation method (to analyze the Gibbs energy)are used to determine all the thermodynamic parameters. For example, Oliveira et al. 32 studied the thermodynamics of [Carbamazepine + Saccharin] co-crystal formation by this method. The second group of experimental approaches focuses on studying the co-crystal component complex formation at one temperature. [33][34][35] In this case, only the Gibbs energy of the co-crystal/salt formation can be calculated.
Several years ago, we proposed an approach 36 to estimating the thermodynamic characteristics of twocomponent crystal formation based on the melting temperatures of the co-crystal and its individual components, as well as on the analysis of the sublimation thermodynamic characteristics of the individual substances. The main sources of such analysis were the databases of the melting points of the co-crystals/salts and the sublimation thermodynamics of the molecular crystals created by us. Because the number of occurrences in the designated databases has risen significantly since the publication of our article, we have attempted to systematize the obtained material and analyze the processes of formation of two-component crystals from the standpoint of the main contributions (enthalpy or entropic) to the Gibbs energy. We also tried to present and justify a general scheme/algorithm for estimating the thermodynamic stability of any co-crystal of interest.  36 It should be noted that in order to simplify the process of finding regularities, we selected only two-component systems/co-crystals for the study. The relevant literature was selected through the SciFinder 37 system by using two key word forms: "cocrystal" and "co-crystal". After that, we chose the papers where the melting points of two-component crystals were analysed by different methods (DSC, Kofler and others). Crystals and individual compounds which melt with decomposition (or partial decomposition) were excluded from the study. The melting points obtained by DSC experiments corresponded to T onset . The melting points obtained by alternative methods were used as they were.
The database included both co-crystals and salts (proton transfer between the co-crystal components). In contrast to our previous analysis, in this paper we tried to make a clear distinction between salts and co-crystals, using the single crystal X-ray diffraction results. Where this was impossible (there were no single crystal X-ray diffraction data), twocomponent systems were referred to as co-crystals. In addition, for each two-component crystal, we included the following information in the database: 1) data from single crystal X-ray (refcode/cif file); 2) evidence (yes or no) of the existence of a co-crystal/salt using XRPD (comparative analysis of diffractograms for a two-component crystal and crystals of individual compounds); 3) data about the presence of DSC experiment results (yes or no) and evidence of the existence of a co-crystal/salt (comparative analysis of DSC curves for twocomponent crystals and crystals of individual compounds); 4) data about the presence of the results of experiments aimed at determining the crystal/salt melting temperature using Kofler's method; 5) data about the presence of solubility experiments, and 6) data about the presence of the results of biological experiments. Special attention was paid to the presence of polymorphic modifications of two-component crystals: each polymorphic phase was described as an independent system (occurrence). The data on the melting points of the co-crystals/salts obtained from crystalline hydrates/ crystalline solvates should be mentioned in particular. The literature contains data on DSC and TG experiments which describe the dehydration/desolvation processes of co-crystals/ salts with the temperatures substantially different from the melting points of the non-solvated two-component crystals. For such systems, the resulting melting temperatures were analyzed together with a common data set. For the analysis, we selected only those co-crystals/salts in which the melting temperature of the individual components was over 25°C. Based on the literature analysis, we have obtained 1175 cocrystals/salts. For 941 (80.2%) crystals, single crystal X-ray diffraction experiments with the solved crystal structures were carried out. Whereas 233 (19.8%) two-component crystals were characterized using only XRPD. The data of single X-ray diffraction experiments allow us to unambiguously interpret the state of the molecules in the crystal: they either take the form of salts (proton transfer between two organic molecules) or a cocrystal (no proton transfer). It should be noted that 233 two-component crystals (described by XRPD experiments) were referred to as co-crystals as this technique does not allow us to make an unambiguous interpretation of the molecule state in the crystal. With such assumptions, the database under consideration contained 184 salts (15.7%) and 991 cocrystals (84.3%). There were 86 occurrences of polymorphic forms (7.3%) (different polymorphic modifications of the same two-component crystal were considered different occurrences). Of these, 70 were monotropic and 16enantiotropic phases (81%/19%).
All the data can be divided into several groups based on their reliability levels: a) those describing the crystal structure of single crystals by using X-ray analysis and the DSC melting curve (745 crystals (63.4%)) (the highest reliability); b) those describing the crystal structure of single crystals by using X-ray analysis and the melting temperature (Kofler's method) (196 occurrences (16.7%)); c) those proving co-crystal formation by comparing the X-ray powder diffraction patterns of the co-crystal/salt and its individual components, and providing the DSC melting curve of the co-crystal/salt (219 occurrences (18.7%)); d) those proving co-crystal/salt formation by comparing the X-ray powder diffraction patterns of the cocrystal/salt and its individual components, and providing the co-crystal/salt melting temperature without a DSC curve (14 occurrences (1.2%)). When different sources described the same two-component crystal, the choice was made in favour of systems belonging to a higher reliability level.
2.1.2. Database of thermodynamic characteristics of sublimation processes of molecular crystals. In order to analyze the experimental data, we used a database created by us on the basis of the values published in the literature from 1900 till 2016 inclusive. The database contains the following information: the method used to measure the thermodynamic characteristics; the experimental temperature interval; the equation describing the temperature dependence of saturated vapor pressure; the sublimation enthalpies both at the experimental temperature, ΔH T sub (exp), and at 298.15 K, ΔH 298 sub (exp), in kJ mol −1 ; the standard Gibbs energy values, ΔG 298 sub (exp), in kJ mol −1 ; the fusion enthalpies and melting temperatures of molecular crystals and, finally, the refcodes 38 of the compounds described by the X-ray diffraction experiments. The database includes enthalpies and Gibbs energies obtained by various methods and at different experimental temperatures; therefore, we used a special algorithm to reduce the experimental values to comparable conditions. For this purpose, we calculated the standard Gibbs energy values at 298 K using the temperature dependences of vapour pressure. In a similar way, if the sublimation enthalpies were obtained at elevated temperatures, these values were recalculated at 298 K based on the heat capacities of the studied compounds. If the same substance had been earlier described in the literature by several methods, preference was given to methods that allowed us to obtain both the sublimation enthalpy and Gibbs energy (saturated vapor pressure) data at the same time: transpiration, torsion effusion, mass effusion, isoteniscope, mercury manometer methods, etc. Moreover, with other conditions being identical, preference was given to the data received at temperatures closest to the standard condition (298.15 K). The saturated vapor pressure at 298.15 K was calculated from the temperature dependence. If the data were only taken for sublimation enthalpies, preference was given to the values obtained by the calorimetric method. In the final analysis, the database included 2899 sublimation enthalpy values. Among these compounds, only 1780 values of the Gibbs energies were known. So, it was evident that the number of the sublimation enthalpy values was 1.6 times bigger than the Gibbs energy ones. This relationship can be explained by the fact that the earlier works described a lot of calorimetric experiments which did not allow for saturated vapor pressure estimation. The experimental data for the sublimation enthalpies were derived by the following methods: 1136mass effusion -Knudsen effusion; 494transpiration; 306calorimetry; 111torsion; 79manometer (mercury, piston); 181 calculation/estimation from vaporization and fusion enthalpies (DSC); 18head space analysis; 7isoteniscope, and 567others. The melting points were known for 2542 molecular crystals. It should be noted that the main reason for the lack of melting temperatures was the decomposition of the compounds under fusion (porphyrins and other substances with a high molecular weight). However, there were a few newly synthesized compounds for which there was no information about the melting points in the original papers. Because the main idea of the work was to find out the relationship between the sublimation Gibbs energies and melting temperatures, we selected 1515 compounds having sublimation Gibbs energies, enthalpies and melting points at the same time.

Algorithm for co-crystal thermodynamic function estimation
The literature data show that in structurally similar compounds (according to Tanimoto similarity coefficients), there is a linear correlation between ΔG 298 sub and T fus . 39 The crystal structure of a two-component co-crystal (with any stoichiometry) can be considered structurally similar to that of the individual compounds. In other words, the experimental values of ΔG 298 sub vs. T fus of the individual compounds and the cocrystal belong to one (and the same) cluster. Therefore, the coefficients A and B of eqn (1) can be calculated if ΔG 298 sub and T fus of the individual compounds are known.
Knowing coefficients A, B and T fus ĲCC), we can calculate the ΔG 298 sub (CC)value: In order to calculate the Gibbs energy of co-crystal formation with the stoichiometric ratio (API) n ĲCF) m ΔG 298 f (CC), we should use the following equations: where (API) n ĲCF) m : X 1 = n/(n + m); X 2 = m/(n + m). ΔG 298 sub (PM), ΔG 298 sub (API), ΔG 298 sub (CF) are the sublimation Gibbs energies of the physical mixture, API and CF, respectively.
The co-crystal formation enthalpy value, ΔH 298 f (CC), was obtained using the following algorithm. It is well known that there is a linear dependence between ΔG 298 sub and ΔH 298 sub (the so-called compensation effect). 40 Thus, knowing the experimental values ΔG 298 sub and ΔH 298 sub , it is possible to calculate the coefficients of eqn (5): and then calculate the necessary value: ΔH 298 sub (PM) = X 1 ·ΔH 298 sub (API) + X 2 ·ΔH 298 sub (CF) (8)

Comparative analysis of the melting points of co-crystals and individual compounds
It was interesting to analyze the relationship between the melting temperatures of the co-crystal and its individual components. It is worth mentioning that the division of the cocrystal components into API and CF is done only for convenience in order to separate the components from each other. And in some cases, API molecules are not the molecules for pharmaceutical application. Nevertheless, we used this symbol because most of the analyzed literature reports are devoted to pharmaceutical co-crystals. The melting points in the figures considered below are presented in Celsius for clarity of the analysis. Fig. 1 shows the results of this analysis. Fig. 1a summarizes the results for all the co-crystals with all possible stoichiometric compositions. These results show that most of the co-crystals/salts have the melting point between the melting points of the individual components (54.6%). In turn, 30.6% (one-third) of the co-crystals/salts have the melting point below the melting points of the individual components (and these are the potential candidates for pharmaceutical co-crystals/salts as their Gibbs energies of crystal lattices are lower than those of the crystals of the individual compounds). However, in this case it should be kept in mind that a melting temperature decrease will reduce the crystal thermodynamic stability. Therefore, pharmaceutical applications require a compromise between the crystal thermodynamic stability and solubility. Finally, 14.8% of co-crystals/salts have the melting temperatures above those of the individual compounds which are the co-crystal/salt components. It is this fact/circumstance that classifies this category of the systems as unfavorable in terms of improving the solubility of poorly soluble components. Fig. 1(b-d) present similar histograms detailed based on different stoichiometric compositions. Each group has its variations, but the general trend remains.
As noted above, the database includes 184 salts and 774 co-crystals (taken from single crystal X-ray diffraction data). It was interesting to compare the distribution of the melting temperatures of two-component crystals among the dedicated groups (I-III) for the co-crystals and salts. The analysis of the results is shown in Fig. 2.
It is evident that for the salts, the proportion of twocomponent crystals with melting temperatures higher than that of the most high-melting component (Higher-II) increases from 13.3% (co-crystals) to 23.4% (1.8 times). In turn, for the salts, the share of two-component crystals of the group (Low-III) remains practically the same in comparison with the co-crystals: 27.2 and 27.7%, respectively. Finally, for the salts, the share of two-component crystals of the group (Between-I) decreases by 9.6% as compared to the co-crystals (from 59.0 to 49.4%).
for the full set obtained by single crystal diffraction experiments (958) and separately for the salts (189) and co-crystals (774). The results are given in Fig. 3(a-c). It is clearly seen that the most frequently encountered melting point for the two-component crystals from the analyzed database for both salts and co-crystals lies in the range of 140-170°C (about 25-28% of the systems based on the total number of samples).
Thus, there is no shift in the maximum of the melting temperature distribution function for the salts in the highermelting region as compared to co-crystals (as is often asserted in the literature).
30 two-component molecular crystals in the analyzed database had two monotropic polymorphic modifications (i.e. polymorphs did not transform into one another when heated and melted at different temperatures). As mentioned above in section 2.1.1, there are 70 occurrences in the database which are associated with monotropic phases. It should be noted that 10 occurrences are associated with three or more polymorphic modifications, so to simplify the analysis we excluded them from consideration. Therefore, it was useful to compare how strongly the melting temperatures of monotropic forms differed from each other. Fig. 4 shows a histogram illustrating the relation between the distribution of the number of polymorphic modifications (two polymorphic forms for one two-component crystal) and the difference in the melting temperatures of these forms (absolute value). The results show that the maximum of the distribution function is in the temperature range from 0 to 5°C (50% of the polymorph pairs under consideration). However, there are polymorphic modifications (20% of the sample used), in which the difference in the melting temperatures lies between 10 and 20°C. This circumstance essentially distinguishes the polymorphic forms of two-component crystals from those of single-component crystals, where the differences in melting temperatures of monotropic modifications of the latter rarely exceed 10°C. 41 Another characteristic feature of the two-component systems (in contrast to the onecomponent crystals) is that there are monotropic forms in which the melting points differ by more than 40 and 50°C.
Thus, for the two-component crystals, polymorphism gives more possibilities to change the Gibbs energy (the energy of the crystal lattice) in comparison with the individual compounds. This, in turn, is one of the ways to obtain pharmaceutical crystals with improved API solubility.  A similar analysis was carried out for two-component crystals with the same composition and different stoichiometries. 42 co-crystals/salts with stoichiometries of (1 : 1) and (1 : 2)/ (2 : 1) were chosen. Fig. 5 shows the distribution function of two-component crystals based on the difference in melting temperatures between co-crystals/salts with the selected stoichiometric compositions. As in the case of monotropic polymorphic forms, the maximum of the distribution function falls within the temperature range of 0-5°C.
However, unlike the previous case, the maximum value is much lower: 31% compared to 50% for the polymorphic forms. The distribution of the co-crystals in the interval from 5 to 10°C is comparable to polymorphs and co-crystals with different stoichiometries (about 23%). However, for the latter (co-crystals), the number of two-component crystals in the interval from 10 to 30°C considerably exceeds that for the former (polymorphs). The number of co-crystals/salts with melting points differing by more than 40°C for polymorphous modifications is comparable with the case under consideration. Thus, a change in the stoichiometric composition of two-component co-crystals makes it possible to lower the melting temperatures much more efficiently than the production of various monotropic polymorphic modifications. This fact becomes one more way of obtaining pharmaceutical cocrystals/salts with improved solubility of the poorly soluble component.
In conclusion, we analyzed the difference in the melting temperatures of co-crystals/salts, where the racemate or its enantiomer is one of the components, while the second component remains unchanged. Fig. 6 presents the analysis results.
The black dots correspond to the difference in melting temperatures of the co-crystals/salts if one of the components is a racemate or an (S)-enantiomer, the red dots represent a racemate or an (R)-enantiomer, and finally, the blue dots correspond to different enantiomers. Evidently, the analyzed temperature differences can reach significant valuesup to 40°C. Therefore, this variant of co-crystal/salt design can be also considered as a way of obtaining pharmaceutical systems with improved solubility.

Thermodynamics of two-component crystal formation
Earlier, we developed an approach to estimating the thermodynamic parameters (Gibbs energy, enthalpy, entropy term) for the formation of co-crystals. 36 The basic equations for calculations are given above (see section 2.2). Thus, in order to evaluate the thermodynamic characteristics of the formation of a two-component crystal, it is necessary to know the following parameters: the melting points of the individual components (T fus ĲAPI), T fus ĲCF)) and the two-component crystal (T fus ĲCC)), and also the Gibbs energy and the sublimation enthalpy of the individual compounds (ΔG 298 sub (API), ΔG 298 sub (CF), ΔH 298 sub (API), ΔH 298 sub (CF)). To begin with, let us consider the   simplest case, when the melting points of the individual substances and two-component crystals are known. There are 1175 such two-component crystals in the database we created. As the next step, we need to select only the co-crystals for which there are data about the Gibbs energies and sublimation enthalpies of the individual compoundsco-crystal/ salt components. For this, we will use the molecular crystal sublimation thermodynamics database created by us (described in section 2.1.2). Analysis of the intersections of the two databases (two-component crystals and sublimation thermodynamics) has shown that for 281 co-crystals/salts, all the necessary information for evaluating the thermodynamic characteristics of their formation is available. characteristics of their formation is available. In order to simplify the analysis of the thermodynamic functions of the cocrystal formation we used the diagram method applied by us in a previous work. 36 In short, it can be described as follows. Every process can be described by three thermodynamic functions (Gibbs energy, enthalpy and entropy). If we plot the enthalpy term along the OX axis in the coordinate plane, and the entropy term along the OY axis, the obtained point accurately describes the thermodynamic process. In this case the Gibbs energy can be evaluated by means of the point position on the isoenergetic curves (straight lines in this case) that are parallel to the YOX angle bisector. Fig. 7 shows the results of calculating the entropy term versus the enthalpy term in the coordinates.
The blue dots denote two-component crystals with a stoichiometry of (1 : 1), whereas the red ones represent the remaining stoichiometry (mainly (1 : 2) and (2 : 1)). The regions of triangular sectors C, D, G and H in Fig. 7 Fig. 7 Thermodynamic functions of two-component crystal formation in the entropy vs. enthalpy co-ordinates. The isoenergetic curves of the ΔG 298 f function are marked by dotted lines. The blue points correspond to the cocrystals/salts with a stoichiometry of (1 : 1), whereas the red points correspond to the cocrystals/salts with a stoichiometry of (2 : 1).
The isoenergetic curves of the ΔG 298 f function are marked as dotted lines in Fig. 7. The calculation accuracy of ΔG 298 f (CC)values was 6-10%. However, several reservations should be made. Firstly, both the obtained co-crystals and the individual compounds can have different polymorphic modifications. It is difficult to predict which polymorphic modifications of individual compounds should be chosen to make the suggested calculations. Moreover, there are no experimental data about the sublimation thermodynamic parameters of some polymorphic modifications. Therefore, we assume that the criterion of thermodynamically stable co-crystal formation must meet the following requirement: It is the threshold value of 5 kJ mol −1 that allows assuming that it is possible to obtain different polymorphic modifications. It should also be noted that the proposed algorithm for estimating the thermodynamics of the formation of twocomponent crystals does not work if the melting points of individual substances of the co-crystal are close to each other. Based on this limitation, 40 out of 281 two-component crystals were excluded (14.2%). Taking into account the proposed criterion of the thermodynamic stability of co-crystals/salts, 196 crystals out of 241 were predicted correctly (81.3%). If we take into account the crystals excluded because it was impossible to estimate the melting points of the co-crystal components at close melting temperatures, then the prediction accuracy is 69.8%.
Each sector of the diagram corresponds to a definite relationship between the enthalpy and entropy terms of the two-component crystal formation. Therefore, we tried to construct the distribution functions for crystals (thermodynamically stable) by sectors. Fig. 8 illustrates these functions for all the crystals (a) and the crystals with different stoichiometries (b). It is evident that for all the crystals considered, the maximum distribution is in sector D (41.8%). The next (in descending order) are sectors C, A and H with the number of points about 3.6 times smaller than in sector D. If all the analyzed crystals are grouped into enthalpy and entropy determined ones, then the following regularity is obtained: the number of two-component crystals with enthalpy determined processes of co-crystals/salts formation corresponds to 70.9% (139), whereas the entropy determined processes -29.1% (57). Apparently, the obtained regularity partly justifies the criteria of thermodynamic stability of multicomponent crystals used in the literature, analyzing the enthalpy terms only. 14,42 However, such approaches are "not impeccable", and as our analysis shows, the formation of every third multicomponent crystal is determined by entropy factors. Fig. 8b represents a similar analysis for different crystal stoichiometries separately. The observed regularities are similar to those of the general sample with small variations for each group. The number of crystals with enthalpy determined processes of co-crystals/salts formation corresponds to 72.4% (102) for (1 : 1) and 67.2% (37) for (2 : 1)/(1 : 2).
It was interesting to analyze whether there is a regularity in the thermodynamic stability of multicomponent crystals relative to different stoichiometric compositions. For this purpose, 12 co-crystals/salts with the same composition and different stoichiometries: (1 : 1) and (2 : 1)/(1 : 2) were selected for 241 two-component crystals. For the convenience of the analysis, the thermodynamic functions of the phase transition ΔY 298 tr ĲCC(1 : 1) → CC((2 : 1)/(1 : 2))) (where Y corresponds to G, H, S) were calculated. Fig. 9 shows the results of the calculations in the coordinates of entropy versus enthalpy terms of the phase transition being analyzed. It is evident that for 75% of the examined co-crystals/salts (9), the two-component crystals with a stoichiometry of (2 : 1)/(1 : 2) are more thermodynamically stable than those with (1 : 1).
3.2.2. Two-component crystals with known melting points and unknown sublimation thermodynamic characteristics of individual components. As shown above, for 198 twocomponent crystals the databases contain only the melting points, and there are no data about the sublimation thermodynamic functions of the individual substances constituting the co-crystal/salt. For 523 two-component crystals, the databases have records about the melting points and the values of the sublimation Gibbs energy and enthalpy for only one component of the co-crystal/salt. Finally, for 173 twocomponent crystals, they contain the melting points and values of the sublimation enthalpy for one or two components of the co-crystal/salt. For such co-crystals with unknown values of the thermodynamic sublimation functions of the individual compounds, we are developing an approach based on using the database collected by us for the sublimation processes of molecular crystals (see section 2.1.2). The main idea is to select (using the Tanimoto similarity coefficient) from the database the structurally similar compounds with known sublimation thermodynamic functions to the compound with the unknown sublimation characteristics. Within this cluster, a linear correlation model that relates the sublimation Gibbs energies and the melting temperatures of the compounds is constructed. Knowing the melting point of our compound, we can estimate the sublimation Gibbs energy. At the next stage, using the correlation between the sublimation Gibbs energies and the molecular crystal enthalpy in the selected cluster ("compensation effect"), one can estimate the sublimation enthalpy of the unknown compound. The more representative the sample for each cluster is, the more accurately the values of the thermodynamic functions can be estimated. Some ideas of the presented approach are described in the work by Perlovich and Raevsky. 39 3.2.3. Two-component crystals with unknown melting points. The melting temperature of a two-component crystal can only be determined if the crystal itself is obtained. However, it would be very useful to estimate the melting temperature of a co-crystal/salt and its thermodynamic parameters of formation before the screening procedure. In this case, it would be possible to make conclusions about the reasonability of starting the screening procedures. In our previous paper, 36 we presented/described an algorithm for estimating the melting temperature of an unknown co-crystal/salt (T fus ĲCC)) using the database of the melting points of twocomponent crystals. In order to do that, we divided all the co-crystals/salts into groups/clusters. Each of the groups contained two-component co-crystals, in which one component was the same (conventionally marked as API) for all the co-crystals, while the second one was different (conventionally marked as CF). Within each of the clusters, we correlated the melting temperatures of co-crystals with those of coformers using the following equation: By substantially increasing the database from 729 to 1175 two-component crystals, we were able to analyze 112 independent clusters and to obtain for 74 of them the correlation equations for the individual compoundscomponents of the co-crystal/salt. The results are presented in Table 1. Unfortunately, for 38 clusters, no correlations were found. Table 1 clearly shows that the pair correlation coefficients for the obtained correlation equations are relatively low but high enough to predict T fus ĲCC) values with mean square deviations from 3 to 24°(only for 1,2-phenylenediamine, this value reaches 38.9°due to a small number of points in the cluster). The number of points in the cluster varies from 4 to 48. Thus, we can estimate the melting temperature of the unknown cocrystal/salt, one of the components of which is the compound selected for screening, and the second component is the compound for which there is a correlation equation given in Table 1. After that, we can begin to determine the thermodynamic characteristics of the co-crystal formation according to the algorithms described in sections 3.2.1 and 3.2.2.
3.2.4. General algorithm for estimating the thermodynamics of two-component crystal formation. In general, the application of the developed approach might look like this (Fig. 10). Let us assume that before obtaining the co-crystal/ salt experimentally, we need to evaluate the possibility of its formation between the pre-selected API and the coformer. At the initial stage, we must analyze the available correlation equations in order to find out if there is an equation relating the melting temperatures of two-component crystals including the chosen coformer among them. If such an equation exists, then by knowing the melting temperature of the API we can estimate the co-crystal/salt melting temperature. The next step is to find out if the database of the sublimation thermodynamic parameters contains the values of the sublimation thermodynamic functions for both the API and the coformer. If such values are included in the database, we can calculate the thermodynamics of co-crystal formation. If not,  then there are two possible ways to solve the problem. Firstly, it is possible to perform the sublimation experiments for the co-crystal components. It should be noted that these experiments are quite routine. As a rule, it takes one to several weeks (depending on the compound being studied) to obtain a complete thermodynamic picture of the sublimation of a single molecular crystal. However, it should be emphasized that the set of coformers used to design co-crystals is usually the same. Therefore, over time, a library of the sublimation thermodynamic parameters for these compounds will be formed and it will not be necessary to waste time on carrying out the indicated experiments. Secondly, if it is not possible to conduct the sublimation experiments, we can estimate the thermodynamic parameters of sublimation using the sublimation database for molecular crystals and the clusterization approach described briefly in section 3.2.2. After that, one can calculate the necessary thermodynamic functions for the co-crystal/salt formation. If the Gibbs energy of the crystal formation is lower than the value chosen as the criterion (5 kJ mol −1 ), then the probability of obtaining a co-crystal/salt is high, and one can proceed to its production. Otherwise, you need to use a different coformer. It is of interest that if ΔG 298 f (CC) < −5 kJ mol −1 (i.e. greater than the energies of polymorphic transitions), then any method can be safely chosen for screening (dry grinding techniques, with an addition of solvent, slurry, recrystallization). If ΔG 298 f (CC) does not fit into the indicated criterion, it is necessary to carefully choose the method with which to start. In our opinion, for the reaction to run at low driving force values, the crystals of the initial substances must have an excessive number of defects. Such conditions can be achieved by dry grinding techniques and grinding with an addition of solvent.

Conclusions
Based on literature analysis, we have built up a database including the melting temperatures of two-component molecular co-crystals and individual compounds (1175 co-crystals/ salts). In order to estimate the thermodynamics of two-component crystal formation, we created a database based on the values published in the literature from 1900 till 2016 inclusive. The database includes enthalpies and Gibbs energies of individual molecular crystals obtained by various methods.
Having analyzed all the co-crystals/salts, we have found out that the melting points of 54.6% of the co-crystals are in the range between (I), of 14.8% higher (II) than, and of 30.6% lower (I) than those of the individual compounds.
A comparative analysis of the distribution of the melting temperatures of two-component crystals has been carried out for the selected groups (I-III) of the co-crystals and salts. It has been found that for the salts, the share of twocomponent crystals with melting temperatures higher than those of the highest-melting component (Higher-II) increases from 13.3% (co-crystals) to 23.4% (1.8 times). In turn, for the salts, the share of two-component crystals of the group (Low-III) remains practically the same in comparison with the co-