Visualization of freezing process in situ upon cooling and warming of aqueous solutions.

The freezing of aqueous solutions and reciprocal distribution of ice and a freeze-concentrated solution (FCS) are poorly understood in spite of their importance in fields ranging from biotechnology and life sciences to geophysics and climate change. Using an optical cryo-microscope and differential scanning calorimetry, we demonstrate that upon cooling of citric acid and sucrose solutions a fast freezing process results in a continuous ice framework (IF) and two freeze-concentrated solution regions of different concentrations, FCS1 and FCS2. The FCS1 is maximally freeze-concentrated and interweaves with IF. The less concentrated FCS2 envelops the entire IF/FCS1. We find that upon further cooling, the FCS1 transforms to glass, whereas the slow freezing of FCS2 continues until it is terminated by a FCS2-glass transition. We observe the resumed slow freezing of FCS2 upon subsequent warming. The net thermal effect of the resumed freezing and a reverse glass-FCS1 transition produces the Ttr2-transition which before has only been observed upon warming of frozen hydrocarbon solutions and which nature has remained misunderstood for decades.

Liquid water, arguably the most important solvent on Earth, rarely occurs in pure state but rather as a component of aqueous solutions. In contrast, the solid form of water, ice, is highly intolerant to impurities1. Hence, upon freezing aqueous solutions separate into pure ice and a FCS which vitrifies23456789101112 or freezes1314 upon further cooling. This phase separation and FCS distribution within the ice play an important role in various natural, industrial and biotechnological processes. For example, FCS veins/pockets within the ice affect microbial activity in ice sheets15, hydromechanics of freezing soils16, rheology and transport properties of glaciers1718 and sea ice1920. In the atmosphere, FCS around cloud ice particles affects physical and chemical properties of cirrus clouds2122 and the rate of stratospheric ozone destruction2223 and, consequently, impacts the climate. When living matter freezes, growing extra- and intracellular ice disrupts cell membranes and this together with other freeze-induced stresses (the formation of FCS, cellular dehydration, etc.) is fatal to cells2425. Freeze-induced separation is crucial in freeze-desalination of sea water26, freeze-purification of waste water27, food industry89101128293031 and biotechnology, particularly, in freeze-drying (lyophilization) which is used to extend the stability and shelf life of foods89101128293031 and labile drugs, especially pharmaceutical proteins23456732, because degradation reactions are decelerated in lyophilized products2345678910112829303132.

Lyophilization is a time- and energy-intensive process which besides freezing consists of primary drying, sometimes preceded by annealing33, and secondary drying2345673334353637 performed upon subsequent warming. The duration of drying is largely determined by the freezing step3334353637. Vitrified or crystallized FCS creates a solid matrix suitable for drying. The morphology of the ice/FCS-matrix controls product resistance to vapour flow of sublimated ice during primary drying, desorption of residual water from a resulting porous cake during secondary drying, and the quality attributes of final lyophilized products such as product porous structure, physical state, residual moisture, reconstitution time, etc.2345673334353637. Freezing methods impose constraints on ice/FCS-matrix morphology. Methods which involve small formulation supercoiling and small cooling rate, produce fewer and larger ice crystals which makes primary drying faster and leaves larger pores in a cake after ice sublimation2345673334353637. However, the genuine ice/FCS morphology formed during freezing is not known. Currently, it is believed that freezing produces ice crystals embedded and dispersed in a matrix of glassy and/or crystallized FCS233334353637. However, such seeming picture of ice/FCS morphology cannot account for the appearance of two transitions, T and T, calorimetrically observed upon warming of frozen carbohydrate solutions91011122930313839. The cold transition, T, is usually related to a glass transition of FCS. Hitherto the nature of the warm transition, T, and the question of whether T or T should be related to the glass transition of maximally FCS, T', has remained a subject of debate for decades91011122930313839. The knowledge of T' is important for the determination of collapse temperature, T, at which lyophilized products start losing their amorphous structure4041. The primary drying is performed at a product temperature, T, which is slightly below T ≈ T + 2 K340.

Visualization of the freezing process in situ would reveal the genuine morphology of ice/FCS. Unfortunately, using an optical cryo-microscope (OC-M), the freezing of bulk solutions is seen as an abrupt black flash because of light scattering from numerous rapidly formed ice crystals42. Finding methods for observing the freezing process in situ is challenging but would crucially improve our knowledge of the freezing phenomenon and understanding of the variety of natural and biotechnological processes. In this work, we observe the freezing process of ‘2-dimensional' samples (5–10 μm films) of citric acid (CA) and sucrose solutions in situ with an OC-M. The term “2-dimensional” is used in the following solely as a shorthand to discriminate the samples used in OC-M from those in DSC. It does not imply a different dimensionality of the physics, such as different nucleation behavior or different dimensionality of the growing ice-network, since a thickness of 5–10 μm is still large in comparison to the size of a critical nucleus and large in comparison to the characteristic dimension of the ice structures observed in OC-M. The choice of solutes was motivated by the fact that CA is widely used in food industry, pharmaceutics4344, tissue engineering45, and sucrose, being a natural lyoprotectant, is important in life sciences12, food industry91011293031, biotechnology3839 etc. We also investigate ‘3-dimensional' (bulk) samples of the same solutions with differential scanning calorimetry (DSC). The obtained OC-M and DSC results are mutually complementary and give a clear picture of the freezing process and formed ice/FCS morphology.

Results

Figure 1A displays the thermograms of ‘3-dimensional' 10, 30, and 55 wt% CA solutions. Exothermic, T, and endothermic, T, peaks are produced by the enthalpy of fusion emitted during the freezing to pure ice and absorbed during equilibrium ice melting46, respectively. Regions without transitions are seen as a straight baseline. The different shape of T-peaks shows that freezing is a fast process in diluted solutions and is hindered by increasing viscosity in a concentrated 55 wt% CA solution. The long low-temperature tail of T-peaks indicates that ice starts to melt gradually from an ice/FCS interface where FCS concentration is largest. In OC-M observations of the freezing process of ‘3-dimensional' solutions, we always observe an abrupt dark flash produced by freezing. In Figure 1B, an OC-M image of a frozen ‘3-dimensional' solution shows that ice morphology and FCS are not distinguishable.

Figure 1

DSC thermograms and OC-M images of CA/H2O.

(A). Upper blue and lower red lines are cooling and warming thermograms, respectively. Skewed lines truncate freezing peaks, Tf, to fit the figure. Concentration (wt%), heat flow (mW), and direction of temperature change (3 K/min) are indicated. (B). Image of frozen ‘3-dimensional' 20 wt% CA at ~211 K. Si marks a silicon substrate. (C). Magnified thermograms from panel (A). The T-transition is a net thermal effect produced by the resumed slow freezing of FCS2 and reverse glass-FCS1 transition, T (see text). Open arrow marks the temperature at which the resumed slow freezing ceases (see text). The meaning of other symbols is given in the text. (D, E). Images of frozen ‘2-dimensional' 10 wt% CA and 52 wt% CA at ~210 and 200 K, respectively. Bright spots are the parts of ice in contact with a cover glass. Arrows mark the channels of FCS1, ice, and a borderline of FCS2 (see text and movies S1, S2). (F). Image of frozen ‘2-dimensional' 55 wt% CA at ~221 K shows that freezing begins from multiple ice nucleating events (movie S3).

Magnified thermograms in Figure 1C for CA/H2O (and in Figure 2A for sucrose/H2O) are more informative than the thermograms in Figure 1A. In addition to the fast freezing process, T-peak, which produces the majority of ice, the magnified cooling thermograms reveal a slow freezing process13, which manifests itself through an inclined exotherm on the cold side of T-peak. We observed the fast and slow freezing processes also in OC-M measurements of ‘2-dimensional' solutions of all concentrations, including the solutions whose thermograms are presented in Figures 1 and 2 (movies S1 and S4). The cooling thermograms also reveal two liquid-glass transitions, T and T, which are recognized by the appearance of two steps, ΔC and ΔC, produced by the heat capacity change46. The ΔC and ΔC steps are only very subtle in the thermograms of 10 wt% CA (and 10 wt% sucrose in Figure 2A) because of the small amount of FCS formed. The warming thermograms reveal a reverse glass-liquid transition, T, and the T-transition91011122930313839. The existence of two liquid-glass transitions upon cooling of CA/H2O and sucrose/H2O and the T-transition during the warming of frozen CA/H2O, to our best knowledge, has not been reported before.

Figure 2

DSC thermograms and OC-M images of sucrose/H2O.

(A). The thermograms are obtained from ‘3-dimensional' 10 and 40 wt% sucrose. All symbols have the same meaning as in Figure 1C. (B, C) Images of ‘2-dimensional' 10 wt% and 5 wt% sucrose taken at ~253 K demonstrate that freezing is triggered heterogeneously from a single ice nucleating event. Spherulitic IF is seen as bright tortuous ice needles/plates interweaved with the dark spots/channels of FCS1 (see also Figure 1D). The images of 10 wt% (B), 5 wt% (C), and 40 wt% sucrose (D) show how IF morphology changes with increasing concentration.

Discussion

The existence of T and T upon cooling requires the existence of two reverse glass-liquid transitions upon warming and, consequently, the formation of two FCS regions of different concentrations during freezing. The fact that T-transition is on the inclined thermogram (Figures 1C and 2A) suggests that the slow freezing and T-transition occur simultaneously, which also requires the existence of two FCS regions of different concentrations. In OC-M images of frozen ‘2-dimensional' CA/H2O and sucrose/H2O, the first region, FCS1, is seen as dark tortuous channels/pockets in between bright tortuous ice twigs (Figures 1D and 2B) or ice needles/plates (Figure 2C). OC-M images also demonstrate that supercooled diluted solutions freeze heterogeneously from a single ice nucleating event triggered by a foreign particle (Figure 2B) or substrate (Figure 2C). After nucleation, ice propagates rapidly as radial ~2–4 μm-thick tortuous twigs which form a continuous ice framework (IF) immersed into FCS1. We also observe that, as concentration increases, IF becomes a dendritic multi-branching pattern (Figures 1E and 2D). This dendritic morphology arises from growth instabilities brought about by insufficiently fast latent heat conduction and solute exclusion from ice during fast freezing. Concentrated solutions freeze from multiple ice nucleating events (Figure 1F and movie S3). Thus, our OC-M observations demonstrate that freezing supercooled solutions produce a continuous IF immersed into FCS1 and not isolated ice crystals, as has previously been believed.

The second, less concentrated (see below) region, FCS2, is formed in front of the advancing IF/FCS1 front and envelops the entire IF/FCS1 (movie S2, Figures 1D, 1E and 2D). Due to the limited rate of low-temperature diffusion of H2O to ice, a concentration gradient is established between FCS1 and FCS2. However, the volume of the transition region is much smaller than that of FCS1 and FCS2 and, consequently, only T and T are visible in the thermograms.

In OC-M measurements, we observe that, as temperature decreases, the slow freezing of FCS2 slows down due to increasing viscosity and ultimately ceases at ~208 K in CA/H2O and ~230 K in sucrose/H2O (movies S1 and S4). The fact that these temperatures coincide with the onset of liquid-glass transition, T, (Figures 1C and 2A), indicates that the FCS2 is associated with the T-transition and, consequently, it is less concentrated than FCS1, which itself vitrifies at T. Upon subsequent warming, the slow freezing resumes also at ~208 K and ~230 K in CA/H2O and sucrose/H2O, respectively (movies S1 and S4). In warming thermograms, these temperatures are the end of reverse glass-FCS2 transition, T (Figures 1C and 2A), where the viscosity of FCS2 has decreased sufficiently for resumed slow ice growth.

In Figure 3, we present the images which captured the onset and end of the resumed slow freezing of FCS2 upon warming. They show that the resumed freezing continues to ~230 K in CA/H2O and ~245 K in sucrose/H2O (movies S1 and S4) i.e., it completely covers the temperature region of the T-transition (Figures 1C and 2A). From this fact and from what was stated above, we conclude that the T-transition is a net thermal effect produced by the resumed slow freezing of FCS2, which is responsible for the exothermic feature of the T-transition, and reverse glass-FCS1 transition, T, which produces the ΔC-step. This solves the long-standing problem of the T-transition and accounts for the appearance of “non-reversing” (crystallization) and “reversing” (glass transition) events in modulated DSC scans of the T-transition93038. We determine the onset temperature of T-transition at ~217 for CA/H2O and ~239 K for sucrose/H2O.

Figure 3

OC-M images of frozen CA/H2O and sucrose/H2O.

(A). Pairs of images (1, 2) taken from ‘2-dimensional' 10 wt% CA and (3, 4) from 20 wt% CA show that upon warming of frozen solutions the resumed slow freezing of FCS2 continues to ~230 K. (B). Pairs of images (5, 6) taken from 40 wt% sucrose and (7, 8) from 45 wt% sucrose show that the resumed slow freezing continues to ~245 K. Squares show locations where the resumed ice growth is visible best.

Both Tw and Tw are characteristic and reproducible temperatures which are independent of the initial solution concentration, as has been seen before for carbohydrate solutions10. Since T increases with concentration and Tw > Tw, we relate the Tw to the glass transition of maximally FCS, T, and the concentration of FCS1 to the maximal freeze-concentration, C'38. This solves another long-standing problem, namely, the problem of T and C'. We calculate C' ≈ 81 wt% and C ≈ 75 wt% for CA/H2O and C' ≈ 85 wt% and C ≈ 81 wt% for sucrose/H2O using the Gordon-Taylor approach4748. In the calculations, we use T = 284 K for pure CA43, T = 335 K for pure sucrose48, and the Gordon-Taylor parameter of kGT ≈ 5.43 for sucrose/H2O48 and our calculated kGT ≈ 3.46 for CA/H2O. Our value of C' ≈ 85 wt% for sucrose/H2O is larger than the literature data of C' ≈ 82 wt%38.

The fact that the resumed slow freezing continues to ~230 and 245 K (Figures 1C and 2A) implies that upon warming of frozen CA/H2O and sucrose/H2O, the fraction of FCS2 remains in liquid phase above T ≈ T + 2 K. This suggests that if the freezing behaviour of pharmaceutical formulations is similar to that described above, then the remaining liquid FCS2 can form a ‘skin' on top of formulations and resist the vapour flow of sublimated ice during primary drying at T < Tc. Further, in the case of CA/H2O, the resumed freezing of FCS2 and ice melting at ice/FCS1 interface occur simultaneously between ~220 and 230 K, because this temperature range is well on the ice melting endotherm (Figure 1C). The resumed freezing increases FCS2 concentration, whereas ice melting at the ice/FCS1 interface, which starts at the end of the Tw-transition at ~220 K, decreases FCS1 concentration. Above ~230 K, when the two concentrations become equal, only ice melting continues. In sucrose/H2O, the annihilation of the concentration gradient between FCS1 and FCS2 takes place between ~242 and 245 K (Figure 2A). Simultaneous freezing and ice melting are best seen upon warming of ‘2-dimensional' 62 wt% CA previously cooled to 173 K (Figure 4 and movie S5).

Figure 4

OC-M images and DSC thermograms of 62 wt% CA.

(A). Images (1–4) are taken upon warming of ‘2-dimensional' 62 wt% CA previously cooled to 173 K. We observed no freezing upon cooling. Upon warming, freezing starts at ~220 K from multiple ice nucleating events (movie S5) and continues to ~240 K (image 4). Ice melting in FCS1, beginning at ~233 K (image 1), is seen as increasing brightness of IF/FSC1 region. Between ~233 and 240 K freezing and ice melting occur simultaneously. (B). The cooling thermogram of ‘3-dimensional' 62 wt% CA shows no indication of freezing, as in ‘2-dimensional' 62 wt% CA. A liquid-glass transition occurs at ~189 K. Upon warming, freezing starts also at ~220 K. Arrows show the temperature region in which freezing and melting occur simultaneously.

A natural question may arise concerning the extent to which conclusions about the ice/FCS morphology of the bulk ‘3-dimensional' solutions can be drawn from the OC-M data of ~5–10 μm-thick solutions, which we call ‘2-dimensional' solutions only in order to distinguish them from large drops. The necessity to use micrometer-scaled thick solutions arises because one can only focus on approximately one micron-thick layer in the optical microscopy technique. We emphasize that our ‘2-dimensional' solutions are very different from the thin films of just a few molecular layers thickness. Whereas the physics and chemistry in such thin films is dominated by surface processes, our ‘2-dimensional' solutions are large (1 cm in diameter, see Methods section) and thick enough to behave as bulk solutions and, consequently, produce the ice/FCS morphology similar to that in bulk samples. Besides the parallels in the DSC and OC-M data discussed above, this is, e.g., also confirmed by the similarity of our pictures in Figures 1D and 2B and pictures obtained with a cryo-scanning electron microscope (C-SEM) in Figure 6 in Ref. 9. Our pictures show that below Tg', FCS1 is amorphous (glassy) and represents a porous matrix (cake) with the pores filled with ice. Similarly, pictures in Figure 6 in Ref. 9 show an amorphous porous cake which was obtained after ice sublimation at 238 K < Tg' ≈ 239 K from frozen ‘3-dimensional' 40 wt% sucrose. Thus our ‘2-dimensional' and ‘3-dimensional' solutions freeze similarly and produce a continuous IF (not isolated ice crystals as previously believed) immersed into FCS1 + FCS2.

In conclusion, this study introduces the ‘2-dimensional-solution' strategy as an approach for the visualization of freezing process in situ and, consequently, the determination of ice/FCS morphology of frozen biopharmaceutical formulations. Together with DSC measurements of ‘3-dimensional' bulk solutions this strategy solves the long-standing problems of the T-transition of frozen hydrocarbon solutions and allows the practically precise determination of the critical formulation parameters T, T and T the knowledge of which is crucial for the optimization of lyophilisation process435. Our findings suggest that a continuous IF may also be formed upon freezing of biological cells and organs that may give a new impetus to investigation of the resistance of living matter to freezing and its survival at low temperatures.

Methods

We prepared 10–62 wt% citric acid, (C6H8O7), and 10–45 wt% sucrose (C12H22O11) solutions by dissolution of 99.5% CA (Merck) and 99.5% sucrose (Sigma) in ultrapure water. The freezing behaviour of approximately semi-spherical drops (5.5–6.5 mg), which in the text are referred to as ‘3-dimensional' solutions, were studied with a Mettler DSC822 calorimeter. The drops were cold-sealed in Al crucibles of 40 μl and studied at a cooling/warming rate of 3K/min between 320 and 163 K. Calorimeter calibration and details about DSC measurements are described elsewhere131423. We also employed ~5–10 μm-thick films of CA/H2O and sucrose/H2O (referred to as ‘2-dimensional' solutions) for in situ observation of freezing/melting processes with an optical cryo-microscope (Olympus BX51) equipped with a Linkam cold stage and Linksys32 temperature control and video capture software. The solution films were formed between a Si-wafer and a cover-glass of 1 cm in diameter. Cryo-microscope measurements were performed at cooling/warming rate of 3 and 5 K/min between 320 and 163 K. The temperatures of heterogeneous freezing and melting of ‘3-dimensional' drops and ‘2-dimensional' solutions of the same concentration were quite similar, as expected, because ~5–10 μm-thick solutions are large/thick enough to behave as bulk solutions. More than 300 measurements performed with DSC and OC-M showed very good reproducibility of results.

Author Contributions

A.B. designed the research, performed DSC and OC-M measurements and calculations, collected and analyzed data, and wrote the manuscript. N.B. performed some of DSC measurements, discussed results, and contributed to writing the manuscript. M.J.M., H.T., E.B. and T.L. discussed results and commented on the manuscript.