2.2. Phase Behavior
All title compounds are crystalline at 298.15 K, and their crystal structures identified by XRPD are provided in
Table 2. Subsequently, the phase behavior was investigated in the temperature range from 183 K to the thermal decomposition temperature using heat-flux DSC to confirm/exclude the presence of phase transitions. Records from the first DSC runs are shown in
Figure 2. Temperatures and enthalpies of fusion are compared to TGA results and literature values in
Table 1. The crystal-crystal phase transitions observed in this work are listed and compared with literature values in
Table 3.
For NAVA and NAIA, the onset of decomposition occurs prior to their melting, and the ability of the partially decomposed sample to recrystallize is very poor. Fusion temperatures and enthalpies could be collected from runs with a fresh sample, and their uncertainty was properly increased. However, even for the remaining compounds that decompose far above their melting temperatures according to TGA, slow decomposition was observed when determining fusion temperatures repeatedly. Even here, fusion temperatures and enthalpies were again taken from first runs, where possible. For polymorphs obtained by recrystallization from melt, the uncertainties were increased to account for the lower purity of the sample.
Taking into account the uncertainty of the melting properties due to decomposition, agreement between DSC, TGA, and literature fusion temperatures and enthalpies is generally good for NAGA (phase α), NAAA, and NAVA. Another phase of NAGA (phase β) that was found to preferably crystallize from melt (see
Figure S2 in the Supplementary Materials) was probably not observed before. A positive correlation between the fusion enthalpies and temperatures for a single compound from various sources can be noticed, which is probably a result of using materials with different purities. Repeated melting experiments performed in this work confirmed that the decrease in fusion enthalpy and temperature due to decomposition exhibit a linear correlation in a certain range of impurities.
A small endothermal peak before the melting of NAAA was observed at 424 K (see
Figure S3 in the Supplementary Materials). The reversibility of this enantiotropic crII-crI phase transition is puzzling. At cooling, the reverse transition occurs at around 408 K, but during repeated cycles, areas of both (heating and cooling) peaks decrease more rapidly than would be expected due to the slow decomposition at the given temperatures. After the first melting, the peak was never observed again, although it seemed that it was always the same phase (crI) that melted.
For NAVA, a small endothermic peak was observed at 477.6 ± 2.0 K, but only in one of the two runs with a fresh sample (
Figure S4 in the Supplementary Materials). The initial crystalline structure of the sample was not identified with CSD ref code JEXNAB [
11], but rather with a material containing this structure together with another unresolved polymorph. Therefore, it is possible that the commercial material is a mixture of different crystals, and the occurrence of the transition depends on the polymorph ratio in the used sample.
For NAIA, two endothermic phase transitions at 429.3 ± 1.0 K and 502.0 ± 1.0 K were observed (
Figure S5 in the Supplementary Materials). The first one (crIII-crII) was found to be fairly reversible. The reversibility of the latter one (crII-crI) was not tested since it occurred in the region where the decomposition of the sample proceeds. Neither of these transitions was reported in the previous study [
10]. At the same time, the fusion enthalpy and temperature reported in [
10] differ considerably (+3.1 K and −9%, respectively) from our values. Since the deviations in the fusion enthalpy and temperature do not follow the same direction, the observed difference is likely not due to sample purity, but the previously reported values may belong to another (fourth) polymorph.
The phase transition reported for NALA at 382 K [
7,
10] was not observed in this work using DSC or TGA. Instead, melting of two different polymorphs that can occur separately or concomitantly was observed (see
Figure S6 in the Supplementary Materials). Based on the fusion temperatures and enthalpies of these two polymorphs, they should be related monotropically, and the reported phase transition [
7,
10] may not correspond to the transformation of one of them to the other. Although the XRPD diffractogram of the sample studied did not fit the known structure of JAHZUN10 [
12] and might correspond to a mixture of polymorphs or to an incompletely or imperfectly crystallized sample, the melting peak observed in the first run was smooth and sharp, which is a sign of a crystalline phase with high purity. The fusion enthalpies of NALA reported previously [
7,
9,
10] are lower than our results, which might be explained by our observation that melting of the concomitant mixture exhibited a lower enthalpy than both of the pure phases.
2.3. Heat Capacities
Experimental heat capacities obtained in this work with SETARAM μDSC IIIa, PerkinElmer DSC 8500, and Quantum Design PPMS are listed in
Tables S1–S10 in the Supplementary Materials, including correction scaling factors applied for PerkinElmer DSC 8500 and Quantum Design PPMS results. For NAGA, NAAA, and NALA, liquid phase heat capacity data of reasonable quality were obtained despite the slow decomposition in the liquid phase. For NAIA, heat capacities of both crIII (below 429 K) and crII (above 429 K) were obtained. Available literature data on crystal heat capacities are summarized in
Table 4. To describe the temperature dependence of selected heat capacity data from
Table 4, the following equation proposed by Archer [
14] was used, the parameters of which are given in
Table 5:
where
= 1 K,
= 1 J∙K
−1∙mol
−1, and
Only a single parameter,
di, per each temperature interval is to be optimized, while the values of the other three are imposed by a constraint on the continuity and smoothness of the resulting temperature dependence. Parameter
b can be estimated from the slope of
at temperatures greater than 70 K prior to the optimization procedure [
14]; the universal value
b = 0.13 K
−1 was used in this work.
The experimental heat capacities for all
N-acetyl amino acid amides studied are compared with the smoothed values obtained using Equations (1) and (2) in
Figure 3. The deviations of the selected experimental data (marked bold in
Table 4) from the smoothed values do not exceed 2% (with the exception of the lowest temperatures, where all experimental methods naturally have higher relative uncertainty).
The literature data by Abate et al. [
7] obtained by heat-flux DSC (Mettler DSC 20) are about 3% higher compared to our data, with the exception of NALA, where the agreement is mostly within 1%. The comparison is complicated by the fact that all the compounds seem to occur in several crystalline forms, and in the previous work, the respective crystalline forms were not reported. The same authors resolved the crystalline structures of all five
N-acetyl amino acid amides [
11,
12,
13], but the studied single crystals were prepared specifically for the X-ray experiments and do not need to correspond to the polymorphs studied calorimetrically. Generally, the uncertainty of heat capacity determination by heat-flux DSC is several percent at best, and thus the agreement should be considered reasonable.
The thermodynamic functions of all studied compounds at
T = 298.15 K obtained using Equations (1) and (2) are tabulated in
Table 6; for thermodynamic functions in a wide temperature range (from 0 K to up to 470 K, depending on the thermodynamic stability of the crystalline phase), see
Figure 4 and
Tables S11–S15 in the Supplementary Materials.
The differences in molar isobaric heat capacity Δ
between aliphatic amino acids presented by Pokorný et al. [
1] (L-alanine, L-valine, L-isoleucine, and L-leucine) and their
N-acetyl amides were calculated, and they are shown in
Figure 5a. Although the heat capacities of glycine are also available [
15], glycine was not included because of the existence of several polymorphs and its rather outlying character compared to the other proteinogenic amino acids. An average value of the difference was then fitted to the Archer equation (Equations (1) and (2)), the parameters of which are presented in
Table 7. Data influenced by the lambda transition of L-leucine [
1] were also not considered in the average. For NALA and NAVA, the Δ
trend above 280 and 370 K, respectively, was approximated by a linear function based on a 100 K interval to avoid discontinuity of the curve and a bias towards compounds measured up to higher temperatures.
Deviations of isobaric heat capacities of
N-acetyl amides calculated using the correlation of the average difference between amino acids and their
N-acetyl amides and the experimental heat capacities of the respective amino acids are shown in
Figure 5b. The heat capacity estimated in this way falls within 10 J mol
−1 K
−1 (and within 5% of the experimental values above 50 K) for all tested compounds. The parameters in
Table 7 could therefore serve as a reasonable estimation for the
N-acetyl amides of the remaining amino acids above 50 K.
The predictive ability of this correlation was tested on
N-acetyl-L-proline amide (NAPA) and L-proline. The heat capacities used for these compounds were taken from Abate et al. [
7] and our previous publication [
2], respectively. As shown in
Figure 5b, the predicted heat capacities of NAPA fall within 5 J K
−1 mol
−1 (3% in relative scale) of the estimation, which is well in agreement with the expected accuracy of the estimation and also of the experimental heat capacity data [
7].