Temperature dependence of ice critical nucleus size

Heterogeneous ice nucleation of salt solution in porous media

Water ubiquitously exists with dissolved salt in both natural and engineered porous media, such as soil, rock, concrete, and tissue; therefore, its freezing temperature depression behavior is of particular interest to various scientific communities tackling with mechanics and physics of porous media. To date, it remains elusive which physical mechanism accounts for its freezing temperature depression and how dissolved ions affect it. Herein, a series of pore-scale experiments were designated to investigate the freezing temperature of salt solutions in tubes with varying pore diameters, pore solution volumes, solid-liquid interfacial areas, ion concentrations, and ion types. The results reveal two main findings: (i) the freezing temperature depression of pore solutions is governed by the heterogeneous ice nucleation (HIN) at the water-solid interface, as evidenced by the observation that the freezing temperature decreases with the decreasing solid-liquid interfacial areas, regardless of pore diameter and pore solution volume; (ii) the dissolved salts alter HIN processes via changing the osmotic potential across the ice embryo-liquid water interface, as indicated by the observation that the freezing temperature is mainly determined by the salt concentration irrespective of salt types. Furthermore, the classical nucleation theory model is adapted for the freezing behavior of pore solutions by including an osmotic potential term. The model shows excellent performance in capturing experimental data with various pore solution concentrations, further substantiating the HIN as the physical mechanism governing pore solution freezing.

Polarized Neutrons Observed Nanometer-Thick Crystalline Ice Plates in Frozen Glucose Solution

Spin-contrast-variation (SCV) small-angle neutron scattering (SANS) is a technique to determine the nanostructure of composite materials from the scattering of polarized neutrons that changes with proton polarization of samples. The SCV-SANS enabled us to determine structure of nanoice crystals that were generated in rapidly frozen sugar solutions by separating the overlapped signals from the nanoice crystals and frozen amorphous solutions. In the frozen glucose solution, we found that the nanoice crystals formed a planar structure with a radius larger than several tens of nanometers and a thickness of 2.5 ± 0.5 nm, which was close to the critical nucleation size of ice crystals in supercooled water. This result suggests that the glucose molecules were preferentially bound to a specific face of nanoice crystals and then blocked the crystal growth perpendicular to that face.

Toward Precision Deposition of Conductive Charge-Transfer Complex Crystals Using Nanoelectrochemistry

The lack of understanding for precise synthesis and assembly of nano-entities remains a major challenge for nanofabrication. Electrocrystallization of a charge-transfer complex (CTC), tetrathiafulvalene bromide (TTF)Br, is studied on micro/nanoelectrodes for precision deposition of functional materials. The study reveals new insights into the entire CTC electrocrystallization process from the initial nanocluster nucleation to the final elongated crystals with hollow ends grown from the working electrode to the neighboring receiving electrode. On microelectrodes, the number of nucleation sites is reduced to one by lowering the applied overpotential or precursor concentration. Certain current-time transients exhibit significant induction periods prior to stable nucleus growth. The induction regime contains small fluctuating current spikes consistent with stochastic formation of precritical nanoclusters with lifetimes of 0.1-30 s and sizes of 20-160 nm. Electrochemical analyses further reveal rate, size distribution, and formation/dissipation dynamics of the nanoclusters. Crystal growth of (TTF)Br is further studied on triangular nanoelectrode patterns with thickness of 5-500 nm, which shows a mass-transfer-controlled process applicable for precision deposition of functional (TTF)Br crystals. This study, for the first time, establishes CTC nanoelectrochemistry as a platform technology for precise deposition of conductive crystal assemblies spanning the source and drain electrode for sensing applications.

The water cavitation line as predicted by the TIP4P/2005 model

The formation of vapor bubbles in a metastable liquid, cavitation, is an activated process due to the free energy cost of having both phases at contact. Such an energetic penalty enables the existence of the liquid beyond its thermodynamic borders. Establishing the stability limits of a liquid as ubiquitous as water has important practical implications and has thereby attracted a lot of attention. Different experimental strategies and theoretical analyses have been employed to measure and predict the cavitation line, or the pressure-temperature kinetic stability border of liquid water. Understanding the location of the cavitation line requires knowing the cavitation rate dependence on pressure and temperature. Such dependency is difficult to obtain in experiments, and we use molecular simulations with the TIP4P/2005 model to fill this gap. By deeply overstretching liquid water below the saturation pressure, we are able to observe and quantify spontaneous cavitation. To deal with a lower overstretching regime, we resort to the Seeding technique, which consists of analyzing simulations of a liquid containing a vapor bubble under the theoretical framework of Classical Nucleation Theory. Combining spontaneous cavitation with Seeding, we get a wide overview of the cavitation rate. We study two different temperatures (450 and 550 K) and complement our perspective with the results previously obtained at 296.4 K [Menzl et al., Proc. Natl. Acad. Sci. 113, 13582 (2016)] to establish a broad simulation-experiment comparison. We find a good agreement between simulations and both isobaric heating and isochoric cooling experiments using quartz inclusions. We are, however, unable to reconcile simulations with other experimental techniques. Our results predict a decrease in the solid-liquid interfacial free energy as the liquid becomes increasingly overstretched with a temperature independent Tolman length of 0.1 nm. Therefore, the capillarity approximation underestimates the nucleation rate. Nonetheless, it provides a fair indication of the location of the cavitation line given the steep rate vs pressure dependence. Overall, our work provides a comprehensive view of the water cavitation phenomenon and sets an efficient strategy to investigate it with molecular simulations.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log₁₀ (J/(m³ s)) = -7(5).

Impact of surface nanostructure and wettability on interfacial ice physics

Ice accumulation on solid surfaces is a severe problem for safety and functioning of a large variety of engineering systems, and its control is an enormous challenge that influences the safety and reliability of many technological applications. The use of molecular dynamics (MD) simulations is popular, but as ice nucleation is a rare event when compared to simulation timescales, the simulations need to be accelerated to force ice to form on a surface, which affects the accuracy and/or applicability of the results obtained. Here, we present an alternative seeded MD simulation approach, which reduces the computational cost while still ensuring accurate simulations of ice growth on surfaces. In addition, this approach enables, for the first time, brute-force all-atom water simulations of ice growth on surfaces unfavorable for nucleation within MD timescales. Using this approach, we investigate the effect of surface wettability and structure on ice growth in the crucial surface-ice interfacial region. Our main findings are that the surface structure can induce a flat or buckled overlayer to form within the liquid, and this transition is mediated by surface wettability. The first overlayer and the bulk ice compete to structure the intermediate water layers between them, the relative influence of which is traced using density heat maps and diffusivity measurements. This work provides new understanding on the role of the surface properties on the structure and dynamics of ice growth, and we also present a useful framework for future research on surface icing simulations.

Investigation of the Freezing Phenomenon in Vials Using an Infrared Camera

The freezing phenomenon has a dramatic impact on the quality of freeze-dried products. Several freezing models applied to solutions in vials have been proposed to predict the resulting product morphology and describe heat transfer mechanisms. However, there is a lack of detailed experimental observations of the freezing phenomenon in vials in the literature. Thus, the present work offers new experimental observations of the freezing phenomenon in vials by infrared (IR) thermography. IR imaging allowed each vial’s whole axial temperature profile to be collected during freezing, providing significant insights into the process. Spontaneous nucleation and vacuum-induced surface freezing (VISF), as a controlled nucleation technique, are investigated. Batches having vials in direct contact with the shelf (exchanging heat mainly through conduction) as well as suspended (exchanging heat mainly through natural convection and radiation) were tested. The study used three solutions: sucrose 5%, mannitol 5%, and dextran 10%. SEM images coupled with an automated image segmentation technique were also performed to examine possible correlations between the freezing observations and the resulting pore size distributions. IR thermography was found to be a promising tool for experimentally predicting the resulting product morphology in-line.

Investigations of structural and dynamical mechanisms of ice formation regulated by graphene oxide nanosheets

Recent research indicates that graphene oxide (GO) nanosheets can be used to regulate ice formation by controlling critical ice nucleus growth in water at supercooling temperatures. In addition, the study of ice formation mechanisms regulated by GO nanosheets, a good model system for antifreeze proteins (AFPs), will shed light on how AFPs regulate ice formation in nature. In this work, time-resolved small-angle x-ray scattering (TR-SAXS) and quasi-elastic neutron scattering (QENS) experiments were carried out to investigate the structural and dynamical mechanisms of ice formation regulated by GO nanosheets. Strikingly, a transient intermediate state was observed in TR-SAXS experiments that only exists in the aqueous dispersions with a larger GO size (11 nm). This serves as evidence that the size of GO is critical for regulating ice formation. Elastic neutron scattering results indicate that ice is formed in all samples and thermal hysteresis occurs in GO aqueous dispersions in both H₂O and D₂O. The structural and dynamics information about water molecules in GO, extracted from QENS, reveals different dynamical behaviors of water molecules in GO aqueous dispersions when approaching the ice formation temperature.

Universal Nucleation Behavior of Sheared Systems

Using molecular simulations and a modified classical nucleation theory, we study the nucleation, under flow, of a variety of liquids: different water models, Lennard-Jones, and hard sphere colloids. Our approach enables us to analyze a wide range of shear rates inaccessible to brute-force simulations. Our results reveal that the variation of the nucleation rate with shear is universal. A simplified version of the theory successfully captures the nonmonotonic temperature dependence of the nucleation behavior, which is shown to originate from the violation of the Stokes-Einstein relation.

Anomalous Behavior in the Nucleation of Ice at Negative Pressures

Ice nucleation is a phenomenon that, despite the relevant implications for life, atmospheric sciences, and technological applications, is far from being completely understood, especially under extreme thermodynamic conditions. In this work we present a computational investigation of the homogeneous ice nucleation at negative pressures. By means of the seeding technique we estimate the size of the ice critical nucleus N_{c} for the TIP4P/Ice water model. This is done along the isotherms 230, 240, and 250 K, from positive to negative pressures until reaching the liquid-gas kinetic stability limit (where cavitation cannot be avoided). We find that N_{c} is nonmonotonic upon depressurization, reaching a minimum at negative pressures in the doubly metastable region of water. According to classical nucleation theory we establish the nucleation rate J and the surface tension γ, revealing a retracing behavior of both when the liquid-gas kinetic stability limit is approached. We also predict a reentrant behavior of the homogeneous nucleation line. The reentrance of these properties is related to the reentrance of the coexistence line at negative pressure, revealing new anomalies of water. The results of this work suggest the possibility of having metastable samples of liquid water for long times at negative pressure provided that heterogeneous nucleation is suppressed.

Homogeneous nucleation of sheared liquids: advances and insights from simulations and theory

One of the most ubiquitous and technologically important phenomena in nature is the nucleation of homogeneous flowing systems. The microscopic effects of shear on a nucleating system are still imperfectly understood, although in recent years a consistent picture has emerged. The opposing effects of shear can be split into two major contributions for simple atomic and molecular liquids: increase of the energetic cost of nucleation, and enhancement of the kinetics. In this perspective, we describe the latest computational and theoretical techniques which have been developed over the past two decades. We collate and unify the overarching influences of shear, temperature, and supersaturation on the process of homogeneous nucleation. Experimental techniques and capabilities are discussed, against the backdrop of results from simulations and theory. Although we primarily focus on simple systems, we also touch upon the sheared nucleation of more complex systems, including glasses and polymer melts. We speculate on the promising directions and possible advances that could come to fruition in the future.

Equivalence between condensation and boiling in a Lennard-Jones fluid

Condensation and boiling are phase transitions highly relevant to industry, geology, and atmospheric science. These phase transitions are initiated by the nucleation of a drop in a supersaturated vapor and of a bubble in an overstretched liquid, respectively. The surface tension between both phases, liquid and vapor, is a key parameter in the development of such nucleation stage. Whereas the surface tension can be readily measured for a flat interface, there are technical and conceptual limitations to obtain it for the curved interface of the nucleus. On the technical side, it is quite difficult to observe a critical nucleus in experiments. From a conceptual point of view, the interfacial free energy depends on the choice of the dividing surface, being the surface of tension the one relevant for nucleation. We bypass the technical limitation by performing simulations of a Lennard-Jones fluid where we equilibrate critical nuclei (both drops and bubbles). Regarding the conceptual hurdle, we find the relevant cluster size by searching the radius that correctly predicts nucleation rates and nucleation free energy barriers when combined with Classical Nucleation Theory. With such definition of the cluster size we find the same value of the surface tension for drops and bubbles of a given radius. Thus, condensation and boiling can be viewed as two sides of the same coin. Finally, we combine the data coming from drops and bubbles to obtain, via two different routes, estimates of the Tolman length, a parameter that allows describing the curvature dependence of the surface tension in a theoretical framework.

Antifreeze proteins and homogeneous nucleation: On the physical determinants impeding ice crystal growth

Antifreeze proteins (AFPs) are biopolymers capable of interfering with ice growth. Their antifreeze action is commonly understood considering that the AFPs, by pinning the ice surface, force the crystal-liquid interface to bend forming an ice meniscus, causing an increase in the surface free energy and resulting in a decrease in the freezing point ΔT^max. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔT^max as a function of the average distance d between AFPs, with simulations spanning over 1 µs. First, we show that the lower the d, the larger the ΔT^max. Then, we find that the water-ice-protein contact angle along the line ΔT^max(d) is always larger than 0°, and we provide a theoretical interpretation. We compute the curvature radius of the stable solid-liquid interface at a given supercooling ΔT ≤ ΔT^max, connecting it with the critical ice nucleus at ΔT. Finally, we discuss the antifreeze capability of AFPs in terms of the protein-water and protein-ice interactions. Our findings establish a unified description of the AFPs in the contest of homogeneous ice nucleation, elucidating key aspects of the antifreeze mechanisms and paving the way for the design of novel ice-controlling materials.

Seeding approach to nucleation in the NVT ensemble: The case of bubble cavitation in overstretched Lennard Jones fluids

Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a metastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.

Deep-supercooling for extended preservation of adipose-derived stem cells

Cell preservation is an enabling technology for widespread distribution and applications of mammalian cells. Traditional cryopreservation via slow-freezing or vitrification provides long-term storage but requires cytotoxic cryoprotectants (CPA) and tedious CPA loading/unloading, cooling, and recovering procedures. Hypothermic storage around 0-4 °C is an alternative method but only works for a short period due to its high storage temperatures. Here, we report on the deep-supercooling (DSC) preservation of human adipose-derived stem cells at deep subzero temperatures without freezing for extended storage. Enabled by surface sealing with an immiscible oil phase, cell suspension can be preserved in a liquid state at -13 °C and -16 °C for 7 days with high cell viability, retention of stemness, attachment, and multilineage differentiation capacities. These results demonstrate that DSC is an improved short-term preservation approach to provide off-the-shelf cell sources for booming cell-based medicine and bioengineering.

Pressure Induced Structural Changes of Proteins Affecting the Ice Nucleation Temperature of Pork Loins

This study investigated the effects of pressure-mediated protein changes on the ice nucleation temperature of pork loins. To variate chemical state of meat proteins, pork loin was pressurized at varying pressure levels (100–500 MPa) for 3 min, and moisture content, expressible moisture (EM) and differential scanning calorimetry (DSC) were analyzed. Although, all treatments showed similar moisture content, EM and degree of protein unfolding of pork loin showed different features as of 300 MPa. At moderate pressure treatments (100–200 MPa), all protein fractions were detected in DSC experiments, and pork loin had lower EM than control (p<0.05). Meanwhile, myosin and actin of pork loin treated at greater than 300 MPa were completely unfolded, and the treatments showed high EM compared to control (p<0.05). Unfolding of meat proteins was a factor suppressing ice nucleation, and the ice nucleation temperature tended to decrease with increasing applied pressure level. The ice nucleation characteristics of pressurized pork loin exhibited a potential application in freezing storage of pressurized meat with less tissue damage comparing to freeze fresh meat, and further exploration regarding the quality change after freezing of fresh and pressurized meat was warranted.

Interfacial free energy of a liquid-solid interface: Its change with curvature

We analyze the changes in the interfacial free energy between a spherical solid cluster and a fluid due to the change of the radius of the solid. Interfacial free energies from nucleation studies using the seeding technique for four different systems, being hard spheres, Lennard-Jones, and two models of water (mW and TIP4P/ICE), were plotted as a function of the inverse of the radius of the solid cluster. In all cases, the interfacial free energy was a linear function of the inverse of the radius of the solid cluster and this is consistent with Tolman's equation. This linear behavior is shown not only in isotherms but also along isobars. The effect of curvature on the interfacial free energy is more pronounced in water, followed by hard spheres, and smaller for Lennard-Jones particles. We show that it is possible to estimate nucleation rates of Lennard-Jones particles at different pressures by using information from simple NpT simulations and taking into account the variation of the interfacial free energy with the radius of the solid cluster. Neglecting the effects of the radius on the interfacial free energy (capillarity approximation) leads to incorrect values of the nucleation rate. For the Lennard-Jones system, the homogeneous nucleation curve is not parallel to the melting curve as was found for water in previous work. This is due to the increase in the interfacial free energy along the coexistence curve as the pressure increases. This work presents a simple and relatively straightforward way to approximately estimate nucleation rates.

On the time required to freeze water

By using the seeding technique the nucleation rate for the formation of ice at room pressure will be estimated for the TIP4P/ICE model using longer runs and a smaller grid of temperatures than in the previous work. The growth rate of ice will be determined for TIP4P/ICE and for the mW model of water. Although TIP4P/ICE and mW have a similar melting point and melting enthalpy, they differ significantly in the dynamics of freezing. The nucleation rate of mW is lower than that of TIP4P/ICE due to its higher interfacial free energy. Experimental results for the nucleation rate of ice are between the predictions of these two models when obtained from the seeding technique, although closer to the predictions of TIP4P/ICE. The growth rate of ice for the mW model is four orders of magnitude larger than for TIP4P/ICE. Avrami's expression is used to estimate the crystallization time from the values of the nucleation and growth rates. For mW the minimum in the crystallization time is found at approximately 85 K below the melting point and its value is of about a few ns, in agreement with the results obtained from brute force simulations by Moore and Molinero. For the TIP4P/ICE the minimum is found at about 55 K below the melting point, but its value is about ten microseconds. This value is compatible with the minimum cooling rate required to avoid the formation of ice and obtaining a glass phase. The crossover from the nucleation controlled crystallization to the growth controlled crystallization will be discussed for systems of finite size. This crossover could explain the apparent discrepancy between the values of J obtained by different experimental groups for temperatures below 230 K and should be considered as an alternative hypothesis to the two previously suggested: internal pressure and/or surface freezing effects. A maximum in the compressibility was found for the TIP4P/ICE model in supercooled water. The relaxation time is much smaller than the crystallization time at the temperature at which this maximum occurs, so this maximum is a real thermodynamic feature of the model. At the temperature of minimum crystallization time, the crystallization time is larger than the relaxation time by just two orders of magnitude.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Seeding approach to crystal nucleation

We present a study of homogeneous crystal nucleation from metastable fluids via the seeding technique for four different systems: mW water, Tosi-Fumi NaCl, Lennard-Jones, and Hard Spheres. Combining simulations of spherical crystal seeds embedded in the metastable fluid with classical nucleation theory, we are able to successfully describe the nucleation rate for all systems in a wide range of metastability. The crystal-fluid interfacial free energy extrapolated to coexistence conditions is also in good agreement with direct calculations of such parameter. Our results show that seeding is a powerful technique to investigate crystal nucleation.

Droplet formation and growth inside a polymer network: A molecular dynamics simulation study

We present a molecular dynamics simulation study that focuses on the formation and growth of nanoscale droplets inside polymer networks. Droplet formation and growth are investigated by the liquid-vapor phase separation of a dilute Lennard-Jones (LJ) fluid inside regularly crosslinked, polymer networks with varying mesh sizes. In a polymer network with small mesh sizes, droplet formation can be suppressed, the extent of which is dependent on the attraction strength between the LJ particles. When droplets form in a polymer network with intermediate mesh sizes, subsequent growth is significantly slower when compared with that in bulk without a polymer network. Interestingly, droplet growth beyond the initial nucleation stage occurs by different mechanisms depending on the mesh size: droplets grow mainly by diffusion and coalescence inside polymer networks with large mesh sizes (as observed in bulk), whereas Ostwald ripening becomes a more dominant mechanism for droplet growth for small mesh sizes. The analysis of droplet trajectories clearly reveals the obstruction effect of the polymer network on the movement of growing droplets, which leads to Ostwald ripening of droplets. This study suggests how polymer networks can be used to control the growth of nanoscale droplets.

Unusual structural properties of water within the hydration shell of hyperactive antifreeze protein

Many hypotheses can be encountered explaining the mechanism of action of antifreeze proteins. One widespread theory postulates that the similarity of structural properties of solvation water of antifreeze proteins to ice is crucial to the antifreeze activity of these agents. In order to investigate this problem, the structural properties of solvation water of the hyperactive antifreeze protein from Choristoneura fumiferana were analyzed and compared with the properties of solvation water present at the surface of ice. The most striking observations concerned the temperature dependence of changes in water structure. In the case of solvation water of the ice-binding plane, the difference between the overall structural ordering of solvation water and bulk water diminished with increasing temperature; in the case of solvation water of the rest of the protein, the trend was opposite. In this respect, the solvation water of the ice-binding plane roughly resembled the hydration layer of ice. Simultaneously, the whole solvation shell of the protein displayed some features that are typical for solvation shells of many other proteins and are not encountered in the solvation water of ice. In the first place, this is an increase in density of water around the protein. The opposite is true for the solvation water of ice - it is less dense than bulk water. Therefore, even though the structure of solvation water of ice-binding plane and the structure of solvation water of ice seem to share some similarities, densitywise they differ.

The biological function of an insect antifreeze protein simulated by molecular dynamics

Antifreeze proteins (AFPs) protect certain cold-adapted organisms from freezing to death by selectively adsorbing to internal ice crystals and inhibiting ice propagation. The molecular details of AFP adsorption-inhibition is uncertain but is proposed to involve the Gibbs–Thomson effect. Here we show by using unbiased molecular dynamics simulations a protein structure-function mechanism for the spruce budworm Choristoneura fumiferana AFP, including stereo-specific binding and consequential melting and freezing inhibition. The protein binds indirectly to the prism ice face through a linear array of ordered water molecules that are structurally distinct from the ice. Mutation of the ice binding surface disrupts water-ordering and abolishes activity. The adsorption is virtually irreversible, and we confirm the ice growth inhibition is consistent with the Gibbs–Thomson law.

DOI:http://dx.doi.org/10.7554/eLife.05142.001

Water expands as it freezes. If this happens to the water inside plants and animals, the resulting ice crystals can rupture cells. To prevent this, many plants and animals that live in cold climates have evolved ‘antifreeze proteins’. When a small particle of ice first starts to form, the antifreeze proteins bind to it and prevent the water around it freezing, hence preventing the growth of an ice crystal.

There are many different types of antifreeze protein, and some are more active than others. For example, some insects including the spruce budworm are exposed to extremely cold temperatures—sometimes below −30°C—and these insects have antifreeze proteins that are highly active.

It is not fully understood how different antifreeze proteins interact with ice and prevent the growth of ice crystals. This is largely because, as yet, there are no experimental techniques that make it possible to see how antifreeze proteins and water molecules arrange themselves at the surface of a growing particle of ice. Instead, scientists have developed computer simulations to investigate this process. While many of these studies have provided valuable information, the computational methods used have only recently become powerful enough to analyze how the antifreeze proteins approach the surface of the ice particle.

Kuiper et al. carried out simulations involving a highly active antifreeze protein from the spruce budworm. The results of these simulations revealed that this antifreeze protein does not bind directly to ice; instead, water molecules at the surface of the protein act as a bridge between the protein and the ice. These water molecules are highly ordered and though they have similarities with how water is structured in the ice, they are distinct from the ice lattice itself. Furthermore, this arrangement appears to be important for allowing the spruce budworm antifreeze protein to interact with the ice.

This study provides detailed insights as to how a highly active antifreeze protein helps to prevent ice crystals forming. In the future, the computational simulations used here may be extended to study the dynamics of other antifreeze proteins, and also how crystals of other materials form.

DOI:http://dx.doi.org/10.7554/eLife.05142.002

Mechanism of anisotropic surface self-diffusivity at the prismatic ice–vapor interface

Models of mesoscopic ice roughening require improved understanding of surface attachment kinetics.

A scale down process for the development of large volume cryopreservation

The process of ice formation and propagation during cryopreservation impacts on the post-thaw outcome for a sample. Two processes, either network solidification or progressive solidification, can dominate the water–ice phase transition with network solidification typically present in small sample cryo-straws or cryo-vials. Progressive solidification is more often observed in larger volumes or environmental freezing. These different ice phase progressions could have a significant impact on cryopreservation in scale-up and larger volume cryo-banking protocols necessitating their study when considering cell therapy applications.

This study determines the impact of these different processes on alginate encapsulated liver spheroids (ELS) as a model system during cryopreservation, and develops a method to replicate these differences in an economical manner.

It was found in the current studies that progressive solidification resulted in fewer, but proportionally more viable cells 24 h post-thaw compared with network solidification. The differences between the groups diminished at later time points post-thaw as cells recovered the ability to undertake cell division, with no statistically significant differences seen by either 48 h or 72 h in recovery cultures.

Thus progressive solidification itself should not prove a significant hurdle in the search for successful cryopreservation in large volumes. However, some small but significant differences were noted in total viable cell recoveries and functional assessments between samples cooled with either progressive or network solidification, and these require further investigation.

Size of a crystal nucleus in the isothermal crystallization of supercooled liquid

We present an alternative to classical nucleation theory (CNT). We introduce a size-dependent surface energy into the total Gibbs free-energy of formation of a crystal (ΔG). We consider the free-energy in the core part of the total volume of crystal and the free-energy in the surface-layer part of it, separately, for the evaluation of ΔG. As a result, we present an explicit model to evaluate a characteristic size of an initial nucleus that differs from the critical nucleus of CNT, but whose temperature dependence agrees well with that reported for the temperature dependency initial fold length of isotactic polystyrene and polyethylene in the literature. Our model has fitted the experimental data in the literature with only one adjustable parameter that is defined as nucleation constant. The nucleation constant is the Gibbs free-energy difference between the crystal and supercooled liquid phases for the volume of initial nucleus. We also present an expression to approximate the evolution of free-energy in the surface-layer part of crystal during the crystal growth.