Vibrational Dynamics of a Chiral Smectic Liquid Crystal Undergoing Vitriﬁcation and Cold Crystallization

: Vibrational dynamics in the glass transition and the cold crystallization process of (S)-4 -(1-methyloctyloxycarbonyl) biphenyl-4-yl 4-[7-(2,2,3,3,4,4,4-heptaﬂuorobutoxy) heptyl-1-oxy]-benzoate (3F7HPhH7) was studied by Fourier transform infrared spectroscopy (FTIR) during cooling / heating experimental runs. The measured spectra processing was supported by quantum chemical density functional theory (DFT) calculations (frequency assignments). The perturbation-correlation moving window two-dimensional analysis (PCMW2D) was performed to examine how the height of individual absorption bands change under with temperature. Two-dimensional correlation analysis (2D-COS) was used to detect freezing-in or activation of the stochastic movements during the vitriﬁcation and the cold crystallization processes. Upon cooling, the vitriﬁcation process involves freezing-in of the stochastic movements of ester groups. Upon heating, as the cold crystallization process begins, the ﬁrst to respond are the vibrations of the C–O–C and C = O groups in the rigid core.


Introduction
Although the history of liquid crystals goes back more than 130 years [1,2], these substances are still readily studied in hot pursuit for new materials that would meet ever growing application needs, at the same time building up specific knowledge in the field of fundamental science of soft condensed matter. Liquid-crystal displays (LCD) or temperature sensors are probably the best-known examples of their industrial application. Liquid crystals in their nematic phase are mainly used in LCD but displays based on smectic phases with ferro-or antiferro-electric properties are also known [3,4]. In addition, if the liquid crystal phase undergoes a glass transition, then the resulting glass with anisotropic properties can be used in electro-optical devices for selective reflection of light, for manufacturing amorphous pharmaceuticals or drug carriers, assuring a better bioavailability of active substances [5,6].
Liquid crystal compounds often exhibit a fairly rich polymorphism. It is not uncommon to observe vitrification and cold crystallization [7][8][9][10]. One of such compounds is (S)-4 -(1-methyloctyloxycarbonyl) biphenyl-4-yl 4-[7-(2,2,3,3,4,4,4-heptafluorobutoxy) heptyl-1-oxy]-benzoate, abbreviated as 3F7HPhH7 ( Figure 1). Differential scanning calorimetry (DSC) thermograms as well as the images of polarized light microscopic (POM) textures of 3F7HPhH7 are shown in Figure 2. During cooling, the compound under study transits at 389 K from the isotropic phase to the SmA phase with the characteristic fan-shaped texture, then (at 388 K) it immediately proceeds to the SmC* phase with the broken fan-shaped texture. At 375 K, the transition to the SmCA* phase with the striated fan-shaped texture was visible. At approximately 259 K, the SmCA* phase vitrified. The characteristic signature of glass transition are sharp discontinuities in the heat flow, seen on the DSC thermogram between 250 and 180 K, as well as the appearance of cracks on POM images [11], just as it was noticed in References [12,13]. During heating, the process of a cold crystallization at 272 K was observed. After the transition between crystal phase and the SmCA* phase (at 331 K), the enantiotropic system previously detected during cooling was observed in the reverse order: SmCA*-380 K-SmC*-389 K-SmA-390 K-Iso. It is worth adding that in the vitrification process of the SmCA* phase, the active stochastic movements are frozen. Upon heating, the glass softens and passes to the supercooled SmCA* phase, which then crystallizes. Thus, during the cycle of the vitrification-glass softening-cold crystallization, we do not go back to the starting point. While the vitrification and softening processes involve the same stochastic movements, the same cannot be said of cold crystallization. More information on polymorphism of 3F7HPhH7 is available in Reference [14].  Differential scanning calorimetry (DSC) thermograms as well as the images of polarized light microscopic (POM) textures of 3F7HPhH7 are shown in Figure 2. During cooling, the compound under study transits at 389 K from the isotropic phase to the SmA phase with the characteristic fan-shaped texture, then (at 388 K) it immediately proceeds to the SmC* phase with the broken fan-shaped texture. At 375 K, the transition to the SmC A * phase with the striated fan-shaped texture was visible. At approximately 259 K, the SmC A * phase vitrified. The characteristic signature of glass transition are sharp discontinuities in the heat flow, seen on the DSC thermogram between 250 and 180 K, as well as the appearance of cracks on POM images [11], just as it was noticed in References [12,13]. During heating, the process of a cold crystallization at 272 K was observed. After the transition between crystal phase and the SmC A * phase (at 331 K), the enantiotropic system previously detected during cooling was observed in the reverse order: SmC A *-380 K-SmC*-389 K-SmA-390 K-Iso. It is worth adding that in the vitrification process of the SmC A * phase, the active stochastic movements are frozen. Upon heating, the glass softens and passes to the supercooled SmC A * phase, which then crystallizes. Thus, during the cycle of the vitrification-glass softening-cold crystallization, we do not go back to the starting point. While the vitrification and softening processes involve the same stochastic movements, the same cannot be said of cold crystallization. More information on polymorphism of 3F7HPhH7 is available in Reference [14].
Fourier transform infrared spectroscopy (FTIR) has been widely used to elucidate changes that on a microscopic, molecular level, accompany the phase transitions [15][16][17][18]. Such changes affect the vibrational dynamics and become detected as variations of the parameters of absorption bands (e.g., wavenumbers, full width at half maxima (FWHM's), areas and heights) observed in infrared spectra consecutively measured as a function of temperature [19][20][21]. As far as 3F7HPhH7 is concerned, direct detailed analysis of such properties is difficult because vibrations within the core part, i.e., biphenyl benzoate, especially in the region of 1300-1000 cm −1 associated with ester groups, produce overlapping bands. The two-dimensional correlation analysis (2D-COS), proposed by Noda, opens up new possibilities of interpreting series of vibrational spectra collected as a function of varying external perturbation, such as temperature, pressure, etc., i.e., hyperspectroscopic data [22]. Morita et al. went further, developing a perturbation-correlation moving window two-dimensional analysis (PCMW2D), which is based on 2D-COS analysis and a moving-window two-dimensional correlation analysis (MW2D) [23]. PCMW2D is an elegant method of detecting changes in band characteristics and providing analysis thereof by means of statistical tools [24]. Both correlation analyses were successfully used to interpret the vibrational dynamics of compounds showing glass transition or cold crystallization processes [25][26][27][28].
The aim of the present paper is to study the vibrational dynamics of 3F7HPhH7 during its vitrification and then cold crystallization, using FTIR, supported by quantum chemical density functional theory (DFT) calculations. In the earlier work [14], we introduced the mesomorphic behavior and the kinetics of the isothermal and non-isothermal cold crystallization of 3F7HPhH7, and this work is a continuation of that research. Hyperspectral FTIR data were examined by means of 2D-COS and PCMW2D analyses. FTIR spectra were measured for different cooling-heating cycles to assess how it affects the occurrence of the cold crystallization process in the compound under study.
SmA-390 K-Iso. It is worth adding that in the vitrification process of the SmCA* phase, the active stochastic movements are frozen. Upon heating, the glass softens and passes to the supercooled SmCA* phase, which then crystallizes. Thus, during the cycle of the vitrification-glass softening-cold crystallization, we do not go back to the starting point. While the vitrification and softening processes involve the same stochastic movements, the same cannot be said of cold crystallization. More information on polymorphism of 3F7HPhH7 is available in Reference [14].

Material
3F7HPhH7 was obtained in the Institute of Chemistry, Military University of Technology, Warsaw, Poland [29]. The purity of the sample was higher than 99.4% (using the ASTM E928 Standard Test Method for Determining Purity by Differential Scanning Calorimetry with the use of DSC TA 2500) and additionally confirmed by high performance liquid chromatography-photo diode array-mass spectrometry (HPLC-PDA-MS) analysis (APCI-ESI-dual source) Shimadzu Prominence HPLC-LCMS 2010EV equipped with a polychromatic UV-VIS detector and MS (ESI/APCI) analyzer using a Kinetex C18 and PFP column.

Fourier Transform Infrared Absorption Spectra (FTIR)
FTIR absorption spectra were recorded using a Bio-Rad Digilab FTS 3000 Excalibur spectrometer in the wavenumber range of 4000-400 cm −1 , with a resolution close to 0.4 cm −1 and taking 64 scans. The sample in the form of a thin film was placed between two zinc selenide ZnSe window discs. The film thickness was approximately 20 µm. During the measurements, the spectrometer was purged with dry nitrogen. Temperature investigations have been performed in three sequences: (i) fast cooling from 403 to 173 K with the rate of 5 K/min, (ii) heating with the rate of 2 K/min (after a previous slow cooling with the rate of 2 K/min) and (iii) heating with the rate of 2 K/min (after a previous fast cooling with the rate of 10 K/min). Additionally, the isothermal crystallization at 261 K for 60 min was examined.

Two-Dimensional Analysis (2D-COS)
In many cases, the experimenter is interested not only in how the system under study evolves when subjected to an external perturbation, e.g., to changing temperature, pressure, etc., but also how its specific subsystems (such as molecular moieties, functional groups, etc.) are coupled to each other, by unveiling whether and how their response to such perturbation (or dynamic variable) is correlated. The idea of how to do this came from two-dimensional spectroscopy, originally developed, e.g., for the nuclear magnetic resonance (NMR), and after Noda and co-workers, is known as two-dimensional correlation analysis [17,30]. When FTIR spectra are being collected one after another, as the dynamic variable (here: temperature) is changed, a two-dimensional matrix is built, known as hyperspectral data, whose elements are spectral intensities y(ν, t), being functions of the spectral variable ν, and dynamic variable t. Let the dynamic spectrum be: where y(ν) is t-averaged, or otherwise determined reference spectrum defined over the dynamic variable range of interest. Then, the generalized 2D correlation map is calculated for any pair of spectral variable values (ν 1 , ν 2 ) via Fourier transform in the domain of the dynamic variable t, and hence it becomes a complex function of (ν 1 , ν 2 ) : The real part, Φ(ν 1 , ν 2 ), is commonly referred to as synchronous, while the imaginary one, Ψ(ν 1 , ν 2 ), as asynchronous.
The synchronous map carries information about simultaneous (hence the name) changes of spectral intensity corresponding to each pair (ν 1 , ν 2 ). Diagonal elements represent the autocorrelation function and are therefore always positive. The maxima on a diagonal are then called autopeaks, and they indicate bands whose intensity changes under the influence of external perturbation. Positive off-diagonal peaks (so called positive cross-peaks) reveal the simultaneous increase or decrease of intensities of two bands under external perturbation, while negative off-diagonal cross-peaks show up whenever the intensity of one band increases while that of the other decreases [31]. The asynchronous map (imaginary part) represents out-of-phase changes of spectral intensities. Therefore, it has no diagonal part. A positive cross-peak in the asynchronous map appears at (ν 1 , ν 2 ) when the change of maximum band intensity at the spectral variable ν 1 occurs before changing the maximum band intensity at the spectral variable ν 2 . Should the opposite happen, a negative cross-peak shows up.
In this work, the 2D-COS analysis of hyperspectral data was carried out by means of the 2D-SHIGE program [32].

Perturbation-Correlation Moving Window Two-Dimensional Analysis (PCMW2D)
PCMW2D is a step forward and a significant improvement over conventional MW2D analysis, which operates on autocorrelations only. Basically, MW2D consists in plotting autocorrelation intensities calculated within a small window in the dynamic variable domain, which is then "slid" along the dynamic variable range. PCMW2D also employs the concept of a finite-width window sliding down the dynamic variable (perturbation) axis, but this time, both synchronous and asynchronous correlation maps are calculated. Graphical representation can then be either in the form of two three-dimensional (3D) surface plots (synchronous and asynchronous PCMW2D correlations) or the corresponding contour maps on a 2D-plane: the perturbation variable axis (e.g., temperature) vs. the spectral variable axis (e.g., wavenumber) [33]. This type of hyperspectral data analysis allows one to perform an in-depth study of the temperature dependence of vibrational dynamics of the compound under interest.
In the present work, PCMW2D maps were plotted by the 2D-SHIGE program [32]. In PCMW2D correlation maps, the red areas represent positive correlation intensity, which means the increasing of the intensity at a given wavenumber as a function of temperature, while the blue ones represent negative correlation intensity (and thus decreasing of the intensity).

Computational Details
The quantum chemical DFT calculations were performed using TURBOMOLE v7.3 at the Academic Computer Centre CYFRONET (Cracow) [34]. The calculations were carried out for an isolated molecule of 3F7HPhH7, with the basis set of def2-TZVPPD for all atoms [35], the B3-LYP functional [36] and using the dispersion corrections of D4 [37].
Def2-TZVPPD, the basis set used, is a triple-zeta, property-optimized augmented one. It provides adequate precision at reasonable computational costs. The choice of B3-LYP functional was justified for a number of reasons. It is a hybrid functional (correlation/exchange) known to be one of the best, and certainly the most frequently used one if these type of calculations are expected to deliver precision results. Altogether, def2-TZVPPD + B3-LYP guarantee results that are not only good, but also comparable with those of other researchers, carried out on similar systems. From the many years of experience of the authors, the use of def2-TZVPPD + B3-LYP should be combined with semi-empirical dispersion corrections, as "pure" DFT, for fundamental theoretical reasons, does not include such interactions. They are meaningful in the case of large, flexible molecules, such as rod-like liquid crystals, where middle-and long-range weak non-covalent intra-molecular interactions influence the conformational state of molecules as well as vibrational frequencies, especially those involving light atoms (hydrogen).

DFT Calculations vs. FTIR Results
Theoretical calculations for an isolated 3F7HPhH7 molecule in energetically most-favorable conformation ( Figure 1) were performed in order to interpret measured infrared absorption spectra, as shown in Figure 3. The description of absorption bands with the assigned types of vibrations is contained in Table 1. The calculated spectrum remains in an acceptable agreement with the measured one (including a scaling factor of 0.98) [38], except the bands associated with stretching vibrations of C-H bonds (the changes of the position of these bands on both spectra are connected with the anharmonic effects occurring in these vibrations). Other slight differences in the position of absorption bands in the region of stretching vibrations may be due to the fact that the spectrum was calculated for a single molecule of a specific conformation, while the measured spectrum was obtained for the real sample in which molecules can take a variety of conformations, affecting the effective frequencies of bands [39]. Since such a complex spectrum consists of many overlapping bands, not all vibrations can be further studied in detail. The most significant bands, taken for the analysis and discussed further on, are those in the region of stretching vibrations of carbonyl groups in the chiral center, ν(C=O) chiral_c. , and in the rigid core, ν(C=O) core (at 1750-1700 cm −1 ). These bands are sensitive to various intra-and inter-molecular interactions, and they can be easily analyzed as they are well-separated from other bands [40,41]. Moreover, also other bands (e.g., in the wavenumber regions of asymmetric ν(C-H) asym and symmetric ν(C-H) sym C-H stretching vibrations, of ρ(C-H) rocking vibrations and of asymmetric ν(C-O-C) asym and symmetric ν(C-O-C) sym C-O-C stretching vibrations) were taken into account.

Fast Cooling
The FTIR spectra changes in the wavenumber regions of ν(C-H) (at 3100-2800 cm −1 ), ν(C=O) (at 1800-1580 cm −1 ), ν(C-O-C) (at 1520-1050 cm −1 ) and ρ(C-H) (at 1050-680 cm −1 ) vibrations during the cooling process at the rate of 5 K/min from 403 to 173 K are shown in Figure 4. These regions proved sensitive to the structural changes occurring after the vitrification process, but insensitive to the transition between smectic phases.
Comparing the shape of bands in the SmC A * phase and in the gSmC A * state it can be noticed that the bands of ν(C-H) asym (at 2936 cm −1 ) and ν(C-H) sym (at 2858 cm −1 ) vibrations split into two partially overlapping ones (Figure 4a). Furthermore, in the glassy state, the band at 1712 cm −1 dominates over that at 1736 cm  To examine how the parameters of individual absorption bands change under the influence of temperature, the PCMW2D analysis was performed. The red and blue areas in Figure 5 represent positive and negative correlation intensities, respectively. Changes in intensity between the SmCA* phase and the gSmCA* state as a result of cooling are observed on PCMW2D maps at temperatures higher than glass transition temperature = 259 , determined from DSC results [14]. These differences probably result from the nature of vitrification, which is a kinetic transition.
The ν(C-H)asym (at 2936 cm -1 ) band shows positive correlation peaks, while the ν(C-H)sym (at 2858 cm -1 ) band exhibits negative ones (Figure 5a). Also, negative correlation peaks have been assigned to βPh(C-H) bands (at 1308, 1120, 1073 cm -1 ). While 3F7HPhH7 is cooled to , the mobility of hydrocarbon chain segments decreases to reach the thermodynamically stable state, which is in 3100 3050 3000 2950 2900 2850 2800  To examine how the parameters of individual absorption bands change under the influence of temperature, the PCMW2D analysis was performed. The red and blue areas in Figure 5 represent positive and negative correlation intensities, respectively. Changes in intensity between the SmC A * phase and the gSmC A * state as a result of cooling are observed on PCMW2D maps at temperatures higher than glass transition temperature T g = 259 K, determined from DSC results [14]. These differences probably result from the nature of vitrification, which is a kinetic transition.

Heating after Slow Cooling
After cooling with the rate of 2 K/min, the sample was heated with the same rate to check how the parameters of previously analyzed bands change during the cold crystallization and other phase transitions of 3F7HPhH7. It turned out that when heating the sample after slow cooling, no changes in the characteristics of the absorption bands on the FTIR spectra could be seen in the range corresponding to the cold crystallization or the transitions between smectic phases. However, a change in the intensity of the analyzed bands was visible during transition from the crystal phase to the SmC A * phase ( Figure A1, Appendix A). This effect was observed on infrared spectra at temperatures higher than the Cr-SmC A * transition temperature T cr = 331 K determined by DSC [14]. The reason is that the DSC technique measures the heat flow during the crystallization of the sample, while FTIR spectroscopy is sensitive to a local molecular environment.

Heating After Fast Cooling
During heating with the rate of 2 K/min (after previous fast cooling with the rate of 10 K/min), cold crystallization from the gSmC A * state occurs before the Cr-SmC A * transition, known from the heating after slow cooling. Both are associated with the changes in the FTIR spectra. Those corresponding to the cold crystallization were observed at temperatures lower than T cc = 272 K determined by DSC [14]. The reason is that the DSC responds to the three-dimensional order of the crystal structure, whereas the FTIR technique is sensitive to the short-range correlations that show up in the degree to which a given vibrational mode is coupled to adjacent vibrations or affected by a local mean field created by neighboring molecules.
while the negative ones for ν(C=O)core (at 1736 cm -1 ), ν(C=O)chiral_c. (at 1712 cm -1 ) and βPh(C-H) (at 1308, 1120, 1073 cm -1 ) bands (Figure 10b,c). On the asynchronous maps, the so-called "butterfly pattern" appeared in the C=O stretching region (Figure 10e). Usually, the appearance of such a pattern in an asynchronous map is attributed to a peak shift combined with the intensity changes [45,46]. Moreover, the ν(C-H)asym (at 2936 cm -1 ) and βPh(C-H) (at 1308, 1120, 1073 cm -1 ) vibrations changed before the ν(C-H)sym (at 2858 cm -1 ), ν(C-O-C)asym (at 1283-1264 cm -1 ) and ν(C-O-C)asym (at 1185-1172 cm -1 ) vibrations (Figure 10d,f), respectively. The generalized 2D-COS correlation analysis is presented in Figure 11, in the region of 1750-1710 cm -1 vs. 3000-2840 cm -1 (Figure 11a,c) and 1750-1710 cm -1 vs. 1400-1000 cm -1 (Figure 11b,d). The analysis of asynchronous correlation maps suggests the following order of band changes: ν(C-O-C)asym (at 1283-1264 cm -1 ) and ν(C-O-C)asym (at 1185-1172 cm -1 ) > ν(C=O)core (at 1736 cm -1 ) > ν(C-H)asym (at 2936 cm -1 ), ν(C-H)sym (at 2858 cm -1 ) and βPh(C-H) (at 1308, 1120, 1073 cm -1 ) > ν(C=O)chiral_c. (at 1712 cm -1 ). As the cold crystallization starts, the first to respond are the vibrations of the C-O-C and C=O groups in the rigid core, and then those of C-O-C groups prior to the C=O group. We were also interested in the range of ν(C-H) sym and ν(C-H) asym vibrations, which are sensitive to the cold crystallization. The ν(C-H) asym band shows a small, high-wavenumber shift, and the intensity decreases during the cold crystallization, while the second ν(C-H) sym band behaves oppositely ( Figure 12). The shifts are probably caused by the splitting of peaks at T cc = 261 K, which corresponds to the development of the crystal phase. Also, other bands related to the C-H vibrations exhibit different characteristics during the cold crystallization, e.g., the increase of intensity of the β(C-H) asym band assigned to the CH 3 asymmetric deformation mode (at approximately 1450 cm −1 ), a change of shape of ρ(CH 3 ) and ω(CH 3 ) bands associated with the CH 3 rocking and C-CH 3 wagging modes (at approximately 1130 and 1040 cm −1 , respectively) and the appearance of a new β(C-H) sym band connected with the CH 3 symmetric deformation mode (at approximately 1380 cm −1 ). The new features of previously mentioned bands should be related to the structural order of the methyl side groups during the cold crystallization.
The generalized 2D-COS correlation analysis is presented in Figure 11, in the region of 1750-1710 cm -1 vs. 3000-2840 cm -1 (Figure 11a,c) and 1750-1710 cm -1 vs. 1400-1000 cm -1 (Figure 11b,d). The analysis of asynchronous correlation maps suggests the following order of band changes: ν(C-O-C)asym (at 1283-1264 cm -1 ) and ν(C-O-C)asym (at 1185-1172 cm -1 ) > ν(C=O)core (at 1736 cm -1 ) > ν(C-H)asym (at 2936 cm -1 ), ν(C-H)sym (at 2858 cm -1 ) and βPh(C-H) (at 1308, 1120, 1073 cm -1 ) > ν(C=O)chiral_c. (at 1712 cm -1 ). As the cold crystallization starts, the first to respond are the vibrations of the C-O-C and C=O groups in the rigid core, and then those of C-O-C groups prior to the C=O group.   Stretching vibrations of the C=O groups prove an especially sensitive probe of changes to the local surrounding of a molecule, including weak non-covalent interactions, hydrogen bonds, etc. This feature was already made use of above in this work. To continue, also in this case, we turn our attention to the changes of wavenumbers, FWHM's, areas and heights of ν(C=O)chiral_c. (Figure 13a) and ν(C=O)core (Figure 13b) bands, as a function of temperature. An example of fitting two Lorentz-shaped absorption bands in the carbonyl groups stretching vibrations region was presented in Reference [14]. During cooling with the rate of 5 K/min, a gradual decrease of wavenumber was Stretching vibrations of the C=O groups prove an especially sensitive probe of changes to the local surrounding of a molecule, including weak non-covalent interactions, hydrogen bonds, etc. This feature was already made use of above in this work. To continue, also in this case, we turn our attention to the changes of wavenumbers, FWHM's, areas and heights of ν(C=O) chiral_c. (Figure 13a) and ν(C=O) core (Figure 13b) bands, as a function of temperature. An example of fitting two Lorentz-shaped absorption bands in the carbonyl groups stretching vibrations region was presented in Reference [14]. During cooling with the rate of 5 K/min, a gradual decrease of wavenumber was seen for both bands in the SmC A * phase. Next, the ν(T) of both bands exhibited inflection points, and so did FWHM of the ν(C=O) chiral_c. band at the temperature of about 320 K, i.e., at the glass transition. At the same temperature, the extremes were observed for the remaining characteristics, i.e., the area under the ν(C=O) chiral_c. band as well as the height of both bands reach their minima, while the area and FWHM of the (C=O) core band are at their maxima. Then, the values of the ν(C=O) chiral_c. band parameters slowly increase during cooling, while for ν(C=O) core , the opposite is observed (except the wavenumber) in the gSmC A * state. During heating with the rate of 2 K/min (after previous fast cooling with the rate of 10 K/min), a rapid decrease of both band parameters was seen at approximately 260 K, i.e., at temperature of the sample crystallization. This effect is associated with the increase of intermolecular interactions that affect the carbonyl groups vibrations. Next, the ν(C=O) chiral_c. band parameters increase upon heating, while the opposite situation is visible for the ν(C=O) core band parameters in the crystal phase. The observed changes could be explained by the sensitivity of C=O groups to the short-range order of molecules, which is first built up during crystallization process. Then, a rapid increase of both band parameters upon heating at approximately 350 K was noticed, which is associated with the Cr-SmC A * phase transition. This growth may be associated with the movements of C=O groups relative to each other, and thus with change of the sensitivity of carbonyl groups to inter-molecular interactions with temperature.
Additionally, we examined the spectral evolution in the carbonyl stretching region during the cold crystallization process at 261 K of 3F7HPhH7, which is shown in Figure 14. The intensity of the ν(C=O) chiral_c. band increases gradually with time at T cc = 261 K, while that of the ν(C=O) core band decreases. We concluded that, during the cold crystallization, the vibrations related with carbonyl groups behave differently. It seems rational to assign one component (ν(C=O) chiral_c. ) as characteristic for the crystal phase, while the other (ν(C=O) core ) for the gSmC A * state [28].
Crystals 2020, 10, x FOR PEER REVIEW 15 of 20 approximately 350 K was noticed, which is associated with the Cr-SmCA* phase transition. This growth may be associated with the movements of C=O groups relative to each other, and thus with change of the sensitivity of carbonyl groups to inter-molecular interactions with temperature.
(a) (b) Figure 13. Temperature dependence of wavenumbers, FWHM's, areas and heights of the bands associated with stretching vibrations of: (a) (C=O)chiral_c. and (b) (C=O)core registered during cooling (blue) with the rate of 5 K/min and heating (red) with the rate of 2 K/min (after previous fast cooling with the rate of 10 K/min).
Additionally, we examined the spectral evolution in the carbonyl stretching region during the cold crystallization process at 261 K of 3F7HPhH7, which is shown in Figure 14. The intensity of the ν(C=O)chiral_c. band increases gradually with time at = 261 , while that of the ν(C=O)core band decreases. We concluded that, during the cold crystallization, the vibrations related with carbonyl groups behave differently. It seems rational to assign one component (ν(C=O)chiral_c.) as characteristic for the crystal phase, while the other (ν(C=O)core) for the gSmCA* state [28].

Conclusions
In this paper, the FTIR spectroscopy followed by the PCMW2D and the 2D-COS analyses were applied to investigate the vitrification and the cold crystallization processes of 3F7HPhH7. DFT quantum chemical calculations (def2-TZVPPD + B3-LYP + dispersion corrections) were carried out to find minimum-energy molecule conformation, and then to assist unambiguous band frequency assignment in the experimental spectra. The infrared spectra were measured in the different heating modes (after slow and fast cooling) to see how the change in rate affects the occurrence of the cold crystallization in the compound under study. During heating after slow cooling, only changes in the absorption band characteristics associated with the Cr-SmCA* transition were observed, while during sample heating after previous fast cooling, additional changes corresponding to the cold crystallization process were additionally visible. The FTIR spectra changes were noticed in the wavenumber regions of C-H stretching vibrations (at 3100-2800 cm -1 ), C=O stretching vibrations (at 1800-1580 cm -1 ), C-O-C stretching vibrations (at 1520-1050 cm -1 ) and C-H rocking vibrations (at 1050-680 cm -1 ), because these regions are sensitive to the structural changes occurring in the vitrification as well as in the cold crystallization process.
The correlation analysis was performed to study the changes of the stochastic movements (characteristic for partially disordered phases) during the vitrification and cold crystallization processes. The glass transition consists in changing the stochastic movements involving: first, the ester groups, then, the C=O groups and finally, the side methyl groups. As the cold crystallization begins, the first to respond are the vibrations of the C-O-C and C=O groups in the rigid core, and next is the structural order of the methyl side groups.

Conclusions
In this paper, the FTIR spectroscopy followed by the PCMW2D and the 2D-COS analyses were applied to investigate the vitrification and the cold crystallization processes of 3F7HPhH7. DFT quantum chemical calculations (def2-TZVPPD + B3-LYP + dispersion corrections) were carried out to find minimum-energy molecule conformation, and then to assist unambiguous band frequency assignment in the experimental spectra. The infrared spectra were measured in the different heating modes (after slow and fast cooling) to see how the change in rate affects the occurrence of the cold crystallization in the compound under study. During heating after slow cooling, only changes in the absorption band characteristics associated with the Cr-SmC A * transition were observed, while during sample heating after previous fast cooling, additional changes corresponding to the cold crystallization process were additionally visible. The FTIR spectra changes were noticed in the wavenumber regions of C-H stretching vibrations (at 3100-2800 cm −1 ), C=O stretching vibrations (at 1800-1580 cm −1 ), C-O-C stretching vibrations (at 1520-1050 cm −1 ) and C-H rocking vibrations (at 1050-680 cm −1 ), because these regions are sensitive to the structural changes occurring in the vitrification as well as in the cold crystallization process.
The correlation analysis was performed to study the changes of the stochastic movements (characteristic for partially disordered phases) during the vitrification and cold crystallization processes. The glass transition consists in changing the stochastic movements involving: first, the ester groups, then, the C=O groups and finally, the side methyl groups. As the cold crystallization begins, the first to respond are the vibrations of the C-O-C and C=O groups in the rigid core, and next is the structural order of the methyl side groups.