1
|
Lin X, Zhang C, Hu S, Chen R. Heterogeneous ice nucleation of salt solution in porous media. J Chem Phys 2024; 160:094501. [PMID: 38426515 DOI: 10.1063/5.0190862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2023] [Accepted: 02/12/2024] [Indexed: 03/02/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Xin Lin
- College of Civil Engineering, Hunan University, Changsha 410082, China
| | - Chao Zhang
- Ministry of Education Key Laboratory of Building Safety and Energy Efficiency, College of Civil Engineering, Hunan University, Changsha 410082, China
| | - Shaojie Hu
- Ministry of Education Key Laboratory of Building Safety and Energy Efficiency, College of Civil Engineering, Hunan University, Changsha 410082, China
| | - Renpeng Chen
- Research Center for Advanced Underground Space Technologies, Hunan University, Changsha 410082, China; Ministry of Education Key Laboratory of Building Safety and Energy Efficiency, Hunan University, Changsha 410082, China; and College of Civil Engineering, Hunan University, Changsha 410082, China
| |
Collapse
|
2
|
Kumada T, Nakagawa H, Miura D, Sekine Y, Motokawa R, Hiroi K, Inamura Y, Oku T, Ohishi K, Morikawa T, Kawamura Y, Kawai K. Polarized Neutrons Observed Nanometer-Thick Crystalline Ice Plates in Frozen Glucose Solution. J Phys Chem Lett 2023; 14:7638-7643. [PMID: 37605312 PMCID: PMC10476184 DOI: 10.1021/acs.jpclett.3c01448] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Accepted: 07/21/2023] [Indexed: 08/23/2023]
Abstract
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.
Collapse
Affiliation(s)
- Takayuki Kumada
- Materials
Sciences Research Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
- J-PARC
Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Hiroshi Nakagawa
- Materials
Sciences Research Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
- J-PARC
Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Daisuke Miura
- Materials
Sciences Research Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Yurina Sekine
- Materials
Sciences Research Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Ryuhei Motokawa
- Materials
Sciences Research Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Kosuke Hiroi
- J-PARC
Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Yasuhiro Inamura
- J-PARC
Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Takayuki Oku
- J-PARC
Center, Japan Atomic Energy
Agency, Tokai, Ibaraki 319-1195, Japan
| | - Kazuki Ohishi
- Neutron
Science and Technology Center, Comprehensive
Research Organization for Science and Society (CROSS), Tokai, Ibaraki 319-1106, Japan
| | - Toshiaki Morikawa
- Neutron
Science and Technology Center, Comprehensive
Research Organization for Science and Society (CROSS), Tokai, Ibaraki 319-1106, Japan
| | - Yukihiko Kawamura
- Neutron
Science and Technology Center, Comprehensive
Research Organization for Science and Society (CROSS), Tokai, Ibaraki 319-1106, Japan
| | - Kiyoshi Kawai
- Graduate
School of Integrated Sciences for Life, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8528, Japan
| |
Collapse
|
3
|
Kilani M, Ahmed M, Mayyas M, Wang Y, Kalantar-Zadeh K, Mao G. Toward Precision Deposition of Conductive Charge-Transfer Complex Crystals Using Nanoelectrochemistry. SMALL METHODS 2023; 7:e2201198. [PMID: 36856170 DOI: 10.1002/smtd.202201198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Revised: 12/12/2022] [Indexed: 06/18/2023]
Abstract
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.
Collapse
Affiliation(s)
- Mohamed Kilani
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| | - Mostak Ahmed
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| | - Mohannad Mayyas
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| | - Yifang Wang
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| | - Kourosh Kalantar-Zadeh
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| | - Guangzhao Mao
- School of Chemical Engineering, University of New South Wales (UNSW Sydney), Sydney, New South Wales, 2052, Australia
| |
Collapse
|
4
|
P Lamas C, Vega C, G Noya E, Sanz E. The water cavitation line as predicted by the TIP4P/2005 model. J Chem Phys 2023; 158:124504. [PMID: 37003766 DOI: 10.1063/5.0139470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Cintia P Lamas
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eva G Noya
- Instituto de Química Física Rocasolano, CSIC, C/ Serrano 119, 28006 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Química-Física I (Unidad de I+D+i Asociada al CSIC), Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
5
|
Grabowska J, Blazquez S, Sanz E, Noya EG, Zeron IM, Algaba J, Miguez JM, Blas FJ, Vega C. Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations. J Chem Phys 2023; 158:114505. [PMID: 36948790 DOI: 10.1063/5.0132681] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
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 log10 (J/(m3 s)) = -7(5).
Collapse
Affiliation(s)
- J Grabowska
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E G Noya
- Instituto de Química Física Rocasolano, CSIC, C/ Serrano 119, 28006 Madrid, Spain
| | - I M Zeron
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J Algaba
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - J M Miguez
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - F J Blas
- Laboratorio de Simulación Molecular y Química Computacional, CIQSO-Centro de Investigación en Química Sostenible and Departamento de Ciencias Integradas, Universidad de Huelva, 21006 Huelva, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
6
|
Pereyra RG, Sebastianelli P, Ávila EE. Homogeneous nucleation in supercooled liquid water. Determination of ice germ size and activation energy barrier in Molecular Dynamics simulations. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2068801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Rodolfo G. Pereyra
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- IFEG-CONICET, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
| | - Paolo Sebastianelli
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- School of Chemistry, University of New South Wales, Sydney, Australia
| | - Eldo E. Ávila
- Fa.M.A.F., Universidad Nacional de Córdoba, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
- IFEG-CONICET, Av. Medina Allende s/n, Ciudad Universitaria Córdoba, Argentina
| |
Collapse
|
7
|
Nikiforidis VM, Datta S, Borg MK, Pillai R. Impact of surface nanostructure and wettability on interfacial ice physics. J Chem Phys 2021; 155:234307. [PMID: 34937379 DOI: 10.1063/5.0069896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Vasileios-Martin Nikiforidis
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Saikat Datta
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Matthew K Borg
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| | - Rohit Pillai
- School of Engineering, Institute for Multiscale Thermofluids, The University of Edinburgh, Edinburgh EH9 3FB, United Kingdom
| |
Collapse
|
8
|
Investigation of the Freezing Phenomenon in Vials Using an Infrared Camera. Pharmaceutics 2021; 13:pharmaceutics13101664. [PMID: 34683957 PMCID: PMC8539506 DOI: 10.3390/pharmaceutics13101664] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Revised: 09/30/2021] [Accepted: 10/05/2021] [Indexed: 01/19/2023] Open
Abstract
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.
Collapse
|
9
|
Zhang S, Han J, Luo X, Wang Z, Gu X, Li N, de Souza NR, Garcia Sakai V, Chu XQ. Investigations of structural and dynamical mechanisms of ice formation regulated by graphene oxide nanosheets. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2021; 8:054901. [PMID: 34549075 PMCID: PMC8443304 DOI: 10.1063/4.0000111] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2021] [Accepted: 07/20/2021] [Indexed: 05/14/2023]
Abstract
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 H2O and D2O. 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.
Collapse
Affiliation(s)
| | - Jingjing Han
- Graduate School of China Academy of Engineering Physics, Beijing, China
| | - Xiang Luo
- Graduate School of China Academy of Engineering Physics, Beijing, China
| | - Zhixin Wang
- Graduate School of China Academy of Engineering Physics, Beijing, China
| | - Xudong Gu
- Graduate School of China Academy of Engineering Physics, Beijing, China
| | - Na Li
- National Facility for Protein Science in Shanghai, Zhangjiang Laboratory, Shanghai Advanced Research Institute, CAS, No.333, Haike Road, Shanghai 201210, China
| | - Nicolas R. de Souza
- Australian Nuclear Science and Technology Organization, Lucas Heights, NSW 2234, Australia
| | - Victoria Garcia Sakai
- ISIS Facility, Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, United Kingdom
| | - Xiang-Qiang Chu
- Graduate School of China Academy of Engineering Physics, Beijing, China
- Author to whom correspondence should be addressed:
| |
Collapse
|
10
|
Goswami A, Dalal IS, Singh JK. Universal Nucleation Behavior of Sheared Systems. PHYSICAL REVIEW LETTERS 2021; 126:195702. [PMID: 34047572 DOI: 10.1103/physrevlett.126.195702] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Accepted: 04/15/2021] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Amrita Goswami
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Indranil Saha Dalal
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh 208016, India
| |
Collapse
|
11
|
Bianco V, de Hijes PM, Lamas CP, Sanz E, Vega C. Anomalous Behavior in the Nucleation of Ice at Negative Pressures. PHYSICAL REVIEW LETTERS 2021; 126:015704. [PMID: 33480790 DOI: 10.1103/physrevlett.126.015704] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 10/14/2020] [Accepted: 12/08/2020] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Valentino Bianco
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - P Montero de Hijes
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Cintia P Lamas
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Quimica, Universidad Complutense de Madrid, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|
12
|
Goswami A, Singh JK. Homogeneous nucleation of sheared liquids: advances and insights from simulations and theory. Phys Chem Chem Phys 2021; 23:15402-15419. [PMID: 34279013 DOI: 10.1039/d1cp02617h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Amrita Goswami
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India.
| | - Jayant K Singh
- Department of Chemical Engineering, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India.
| |
Collapse
|
13
|
Sanchez-Burgos I, de Hijes PM, Rosales-Pelaez P, Vega C, Sanz E. Equivalence between condensation and boiling in a Lennard-Jones fluid. Phys Rev E 2020; 102:062609. [PMID: 33466022 DOI: 10.1103/physreve.102.062609] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Accepted: 11/30/2020] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- I Sanchez-Burgos
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - P Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - P Rosales-Pelaez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
14
|
Kang T, Hoptowit R, Jun S. Effects of an oscillating magnetic field on ice nucleation in aqueous iron‐oxide nanoparticle dispersions during supercooling and preservation of beef as a food application. J FOOD PROCESS ENG 2020. [DOI: 10.1111/jfpe.13525] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Taiyoung Kang
- Department of Molecular Biosciences and Bioengineering University of Hawaii at Manoa Honolulu Hawaii USA
| | - Raymond Hoptowit
- Department of Molecular Biosciences and Bioengineering University of Hawaii at Manoa Honolulu Hawaii USA
| | - Soojin Jun
- Department of Human Nutrition, Food and Animal Sciences University of Hawaii at Manoa Honolulu Hawaii USA
| |
Collapse
|
15
|
Bianco V, Espinosa JR, Vega C. Antifreeze proteins and homogeneous nucleation: On the physical determinants impeding ice crystal growth. J Chem Phys 2020; 153:091102. [PMID: 32891082 DOI: 10.1063/5.0023211] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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 ΔTmax. Here, we present an extensive computational study for a model protein adsorbed on a TIP4P/Ice crystal, computing ΔTmax 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 ΔTmax. Then, we find that the water-ice-protein contact angle along the line ΔTmax(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 ≤ ΔTmax, 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.
Collapse
Affiliation(s)
- Valentino Bianco
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0H3, United Kingdom
| | - Carlos Vega
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|
16
|
Rosales-Pelaez P, Sanchez-Burgos I, Valeriani C, Vega C, Sanz E. Seeding approach to nucleation in the NVT ensemble: The case of bubble cavitation in overstretched Lennard Jones fluids. Phys Rev E 2020; 101:022611. [PMID: 32168559 DOI: 10.1103/physreve.101.022611] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 01/24/2020] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- P Rosales-Pelaez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - I Sanchez-Burgos
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Valeriani
- Departamento de Estructura de la Materia, Física Termica y Electronica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
17
|
Huang H, Rey-Bedón C, Yarmush ML, Usta OB. Deep-supercooling for extended preservation of adipose-derived stem cells. Cryobiology 2020; 92:67-75. [PMID: 31751557 PMCID: PMC7195234 DOI: 10.1016/j.cryobiol.2019.11.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 11/12/2019] [Accepted: 11/13/2019] [Indexed: 01/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Haishui Huang
- Center for Engineering in Medicine, Massachusetts General Hospital, Harvard Medical School and Shriners Hospitals for Children, Boston, MA, 02114, United States
| | - Camilo Rey-Bedón
- Center for Engineering in Medicine, Massachusetts General Hospital, Harvard Medical School and Shriners Hospitals for Children, Boston, MA, 02114, United States
| | - Martin L Yarmush
- Center for Engineering in Medicine, Massachusetts General Hospital, Harvard Medical School and Shriners Hospitals for Children, Boston, MA, 02114, United States; Department of Biomedical Engineering, Rutgers University, Piscataway, NJ, 08854, United States.
| | - O Berk Usta
- Center for Engineering in Medicine, Massachusetts General Hospital, Harvard Medical School and Shriners Hospitals for Children, Boston, MA, 02114, United States.
| |
Collapse
|
18
|
Cho Y, Lee EJ, Lee J, Lee S, Yun YC, Hong GP. Pressure Induced Structural Changes of Proteins Affecting the Ice Nucleation Temperature of Pork Loins. Food Sci Anim Resour 2020; 39:1008-1014. [PMID: 31950116 PMCID: PMC6949529 DOI: 10.5851/kosfa.2019.e89] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Revised: 11/18/2019] [Accepted: 11/19/2019] [Indexed: 12/05/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Youngjae Cho
- Department of Food Science and Biotechnology of Animal Resources, Konkuk University, Seoul 05029, Korea
| | - Eun-Jung Lee
- Department of Food Science and Biotechnology, Sejong University, Seoul 05006, Korea
| | - Jiseon Lee
- Department of Food Science and Biotechnology of Animal Resources, Konkuk University, Seoul 05029, Korea
| | - SangYoon Lee
- Department of Food Science and Biotechnology of Animal Resources, Konkuk University, Seoul 05029, Korea
| | - Young-Chan Yun
- Department of Food Science and Biotechnology, Sejong University, Seoul 05006, Korea
| | - Geun-Pyo Hong
- Department of Food Science and Biotechnology, Sejong University, Seoul 05006, Korea
| |
Collapse
|
19
|
Probing the critical nucleus size for ice formation with graphene oxide nanosheets. Nature 2019; 576:437-441. [DOI: 10.1038/s41586-019-1827-6] [Citation(s) in RCA: 152] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2018] [Accepted: 09/17/2019] [Indexed: 12/24/2022]
|
20
|
Montero de Hijes P, Espinosa JR, Sanz E, Vega C. Interfacial free energy of a liquid-solid interface: Its change with curvature. J Chem Phys 2019; 151:144501. [PMID: 31615240 DOI: 10.1063/1.5121026] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- P Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge CB3 0H3, United Kingdom
| | - Eduardo Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
21
|
Montero de Hijes P, Espinosa JR, Vega C, Sanz E. Ice growth rate: Temperature dependence and effect of heat dissipation. J Chem Phys 2019; 151:044509. [DOI: 10.1063/1.5103273] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P. Montero de Hijes
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, United Kingdom
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
22
|
Teboul V, Rajonson G. Simulations of supercooled water under passive or active stimuli. J Chem Phys 2019; 150:214505. [DOI: 10.1063/1.5093353] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Victor Teboul
- Laboratoire de Photonique d’Angers EA 4464, Physics Department, Université d’Angers, 2 Bd Lavoisier, 49045 Angers, France
| | - Gabriel Rajonson
- Laboratoire de Photonique d’Angers EA 4464, Physics Department, Université d’Angers, 2 Bd Lavoisier, 49045 Angers, France
| |
Collapse
|
23
|
Teboul V, Kerasidou AP. Specific properties of supercooled water in light of water anomalies. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1505045] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Victor Teboul
- Physics Department, Université d'Angers, Angers, France
| | | |
Collapse
|
24
|
Soria GD, Espinosa JR, Ramirez J, Valeriani C, Vega C, Sanz E. A simulation study of homogeneous ice nucleation in supercooled salty water. J Chem Phys 2018; 148:222811. [DOI: 10.1063/1.5008889] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Guiomar D. Soria
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge R. Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jorge Ramirez
- Departamento de Ingenieria Quimica Industrial y Medio Ambiente, Escuela Tecnica Superior de Ingenieros Industriales, Universidad Politecnica de Madrid, 28006 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
- Departamento de Fisica Aplicada I, Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
25
|
Espinosa JR, Navarro C, Sanz E, Valeriani C, Vega C. On the time required to freeze water. J Chem Phys 2018; 145:211922. [PMID: 28799362 DOI: 10.1063/1.4965427] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- J R Espinosa
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Navarro
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E Sanz
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Valeriani
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain and Departamento de Fisica Aplicada I , Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
26
|
Curland S, Meirzadeh E, Cohen H, Ehre D, Maier J, Lahav M, Lubomirsky I. The Contribution of Pyroelectricity of AgI Crystals to Ice Nucleation. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201802291] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Sofia Curland
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Elena Meirzadeh
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Hagai Cohen
- Department of Chemical Research Support; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - David Ehre
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Joachim Maier
- Max Planck Institute for Solid State Research; 70569 Stuttgart Germany
| | - Meir Lahav
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Igor Lubomirsky
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| |
Collapse
|
27
|
Curland S, Meirzadeh E, Cohen H, Ehre D, Maier J, Lahav M, Lubomirsky I. The Contribution of Pyroelectricity of AgI Crystals to Ice Nucleation. Angew Chem Int Ed Engl 2018; 57:7076-7079. [DOI: 10.1002/anie.201802291] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2018] [Indexed: 11/05/2022]
Affiliation(s)
- Sofia Curland
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Elena Meirzadeh
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Hagai Cohen
- Department of Chemical Research Support; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - David Ehre
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Joachim Maier
- Max Planck Institute for Solid State Research; 70569 Stuttgart Germany
| | - Meir Lahav
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| | - Igor Lubomirsky
- Department of Materials and Interfaces; The Weizmann Institute of Science; 76100 Rehovot Israel
| |
Collapse
|
28
|
Mandal T, Larson RG. Nucleation of urea from aqueous solution: Structure, critical size, and rate. J Chem Phys 2017; 146:134501. [DOI: 10.1063/1.4979141] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
|
29
|
Sosso G, Chen J, Cox SJ, Fitzner M, Pedevilla P, Zen A, Michaelides A. Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem Rev 2016; 116:7078-116. [PMID: 27228560 PMCID: PMC4919765 DOI: 10.1021/acs.chemrev.5b00744] [Citation(s) in RCA: 379] [Impact Index Per Article: 47.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2015] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Gabriele
C. Sosso
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Ji Chen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | | | - Martin Fitzner
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Philipp Pedevilla
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Andrea Zen
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| | - Angelos Michaelides
- Thomas Young Centre, London
Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street WC1E
6BT London, U.K.
| |
Collapse
|
30
|
Espinosa JR, Vega C, Valeriani C, Sanz E. Seeding approach to crystal nucleation. J Chem Phys 2016; 144:034501. [PMID: 26801035 DOI: 10.1063/1.4939641] [Citation(s) in RCA: 122] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jorge R Espinosa
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Quimica Fisica I, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
31
|
Jung J, Jang E, Shoaib MA, Jo K, Kim JS. Droplet formation and growth inside a polymer network: A molecular dynamics simulation study. J Chem Phys 2016; 144:134502. [PMID: 27059575 DOI: 10.1063/1.4944965] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jiyun Jung
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Eunseon Jang
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Mahbubul Alam Shoaib
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Kyubong Jo
- Department of Chemistry and Interdisciplinary Program of Integrated Biotechnology, Sogang University, Seoul 04107, Republic of Korea
| | - Jun Soo Kim
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| |
Collapse
|
32
|
Zielke SA, Bertram AK, Patey GN. Simulations of Ice Nucleation by Model AgI Disks and Plates. J Phys Chem B 2016; 120:2291-9. [DOI: 10.1021/acs.jpcb.5b06605] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Stephen A. Zielke
- Department
of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - Allan K. Bertram
- Department
of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| | - G. N. Patey
- Department
of Chemistry, University of British Columbia, Vancouver, British Columbia V6T 1Z1, Canada
| |
Collapse
|
33
|
Espinosa JR, Vega C, Valeriani C, Sanz E. The crystal-fluid interfacial free energy and nucleation rate of NaCl from different simulation methods. J Chem Phys 2015; 142:194709. [DOI: 10.1063/1.4921185] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Jorge R. Espinosa
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Chantal Valeriani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Eduardo Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
34
|
Kuffel A, Czapiewski D, Zielkiewicz J. Unusual structural properties of water within the hydration shell of hyperactive antifreeze protein. J Chem Phys 2015; 141:055103. [PMID: 25106616 DOI: 10.1063/1.4891810] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Anna Kuffel
- Department of Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Dariusz Czapiewski
- Department of Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| | - Jan Zielkiewicz
- Department of Chemistry, Gdansk University of Technology, Narutowicza 11/12, 80-233 Gdansk, Poland
| |
Collapse
|
35
|
Kuiper MJ, Morton CJ, Abraham SE, Gray-Weale A. The biological function of an insect antifreeze protein simulated by molecular dynamics. eLife 2015; 4. [PMID: 25951514 PMCID: PMC4442126 DOI: 10.7554/elife.05142] [Citation(s) in RCA: 72] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2014] [Accepted: 05/06/2015] [Indexed: 11/16/2022] Open
Abstract
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
Collapse
Affiliation(s)
- Michael J Kuiper
- Victorian Life Sciences Computation Initiative, The University of Melbourne, Carlton, Australia
| | - Craig J Morton
- ACRF Rational Drug Discovery Centre, St Vincent's Institute of Medical Research, Fitzroy, Australia
| | - Sneha E Abraham
- School of Chemistry, The University of Melbourne, Melbourne, Australia
| | - Angus Gray-Weale
- School of Chemistry, The University of Melbourne, Melbourne, Australia
| |
Collapse
|
36
|
Gladich I, Oswald A, Bowens N, Naatz S, Rowe P, Roeselova M, Neshyba S. Mechanism of anisotropic surface self-diffusivity at the prismatic ice–vapor interface. Phys Chem Chem Phys 2015; 17:22947-58. [DOI: 10.1039/c5cp01330e] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Models of mesoscopic ice roughening require improved understanding of surface attachment kinetics.
Collapse
Affiliation(s)
- Ivan Gladich
- Qatar Environment and Energy Research Institute
- Hamad Bin Khalifa University
- Qatar Foundation
- Doha
- Qatar
| | - Amrei Oswald
- Chemistry Department
- University of Puget Sound
- Tacoma
- USA
| | | | - Sam Naatz
- Chemistry Department
- University of Puget Sound
- Tacoma
- USA
| | - Penny Rowe
- Department of Physics
- Universidad de Santiago de Chile
- Santiago
- Chile
- Geography Department
| | - Martina Roeselova
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- 16610 Prague 6
- Czech Republic
| | | |
Collapse
|
37
|
Espinosa JR, Sanz E, Valeriani C, Vega C. Homogeneous ice nucleation evaluated for several water models. J Chem Phys 2014; 141:18C529. [DOI: 10.1063/1.4897524] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J. R. Espinosa
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - E. Sanz
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Valeriani
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
38
|
Clausse D, Wardhono E, Lanoiselle JL. Formation and determination of the amount of ice formed in water dispersed in various materials. Colloids Surf A Physicochem Eng Asp 2014. [DOI: 10.1016/j.colsurfa.2014.06.032] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
39
|
Kilbride P, Morris GJ, Milne S, Fuller B, Skepper J, Selden C. A scale down process for the development of large volume cryopreservation. Cryobiology 2014; 69:367-75. [PMID: 25219980 PMCID: PMC4271741 DOI: 10.1016/j.cryobiol.2014.09.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2014] [Revised: 09/01/2014] [Accepted: 09/02/2014] [Indexed: 12/13/2022]
Abstract
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.
Collapse
Affiliation(s)
- Peter Kilbride
- Institute of Liver and Digestive Health, Royal Free Hospital Campus, UCL, London NW3 2PF, UK.
| | - G John Morris
- Asymptote Ltd., St. John's Innovation Centre, Cowley Road, Cambridge CB4 0WS, UK
| | - Stuart Milne
- Asymptote Ltd., St. John's Innovation Centre, Cowley Road, Cambridge CB4 0WS, UK
| | - Barry Fuller
- Department of Surgery, Royal Free Hospital Campus, UCL, London NW3 2PF, UK
| | - Jeremy Skepper
- Multi Imaging Centre, Anatomy Building, Downing Site, Cambridge University, CB2 3DY, UK
| | - Clare Selden
- Institute of Liver and Digestive Health, Royal Free Hospital Campus, UCL, London NW3 2PF, UK
| |
Collapse
|
40
|
Reinhardt A, Doye JPK. Effects of surface interactions on heterogeneous ice nucleation for a monatomic water model. J Chem Phys 2014; 141:084501. [DOI: 10.1063/1.4892804] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Affiliation(s)
- Aleks Reinhardt
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Jonathan P. K. Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
41
|
Sanz E, Vega C, Espinosa JR, Caballero-Bernal R, Abascal JLF, Valeriani C. Homogeneous Ice Nucleation at Moderate Supercooling from Molecular Simulation. J Am Chem Soc 2013; 135:15008-17. [DOI: 10.1021/ja4028814] [Citation(s) in RCA: 214] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- E. Sanz
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. R. Espinosa
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - R. Caballero-Bernal
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - J. L. F. Abascal
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Valeriani
- Departamento
de Química
Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
42
|
Lee HS. Size of a crystal nucleus in the isothermal crystallization of supercooled liquid. J Chem Phys 2013; 139:104909. [PMID: 24050367 DOI: 10.1063/1.4820560] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Heon Sang Lee
- Department of Chemical Engineering, Dong-A University, Nakdong-daero 550, Sahagu, Busan 604-701, South Korea
| |
Collapse
|
43
|
Reinhardt A, Doye JPK. Note: Homogeneous TIP4P/2005 ice nucleation at low supercooling. J Chem Phys 2013; 139:096102. [DOI: 10.1063/1.4819898] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Malaspina DC, di Lorenzo AJB, Pereyra RG, Szleifer I, Carignano MA. The water supercooled regime as described by four common water models. J Chem Phys 2013; 139:024506. [DOI: 10.1063/1.4812928] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Reinhardt A, Doye JPK, Noya EG, Vega C. Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water. J Chem Phys 2012. [DOI: 10.1063/1.4766362] [Citation(s) in RCA: 81] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|
46
|
Reinhardt A, Doye JPK. Free energy landscapes for homogeneous nucleation of ice for a monatomic water model. J Chem Phys 2012; 136:054501. [DOI: 10.1063/1.3677192] [Citation(s) in RCA: 86] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
|
47
|
Vega C, Abascal JLF. Simulating water with rigid non-polarizable models: a general perspective. Phys Chem Chem Phys 2011; 13:19663-88. [DOI: 10.1039/c1cp22168j] [Citation(s) in RCA: 658] [Impact Index Per Article: 50.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|