arXiv:0704.1381v1 [physics.chem-ph] 11 Apr 2007 ∗ xlnto fteeet ete h ieeasdnor elapsed predictable. time being the temperature neither freezing effective effect, the statistical the a of toward explanation by points supported investigations, apparently As scenario, experimental a earlier. temperature such freezing local consequence, spontaneously higher a thus hot a water, to Initially cold effect than supercool water. to Mpemba both non-preheated seems the supercooling in water of where and occurrence [3], preheated the in to Ref. related been in has opinion, our in support to seem references not [3]). (see effect cases) in do the these of calculations appearance in the case, quantitatively even any observed in been and, has the (whilst effect experiments Mpemba the during taken are precautions should water in [5]. [2] change heated some when that literature, occur the out in pointing proposed mainly been have observations, such an original pure make not (the to tried cream). for liquids he ice even when different of place occurred or curve observation takes Mpemba cooling solutions effect the with The follow water, initially one. closely the last then of the temperature and the one, to down water cold cool hot to initially which first to has according the view simple adopting his- naive, The when short [3]). counter-intuitive a in sound references (for the observations effect see Mpemba survey scientific the and as torical to which referred phenomenon now a is one, cold initially than quickly more rtto fte 2,[] 3,[] hs eut consist results These [5]. that, [3], fact [4], the inter- in [2], physical them accepted of generally the pretation a and of [1] lacking to results due apparent times experimental recent counter-intuitive in of some interest much freezing the attracted of has that water as such phenomenon well-known A pmaeetadpaetastosi h daai cooling adiabatic the in transitions phase and effect Mpemba orsodn uhr [email protected] author, Corresponding ee epee ofc h rbe ysatn from starting by problem the face to prefer we Here, discussion, the in introduced been has light novel Some some if applied be cannot explanations such However, explain to aimed phenomena, physical possible Several n 1 and ieetodrn fcutr fmlclsi ae,i give is water, in experi transi molecules the phase of by clusters such of Supported and ordering supercooling different earlier. of detected analysis been theoretical not has 0.21, oei lydb uecoiga ela ypaetransition phase by the as in well that as showing supercooling out, by carried played is is one) role cold before water hot nacrt xeietlivsiaino h pmaeffect Mpemba the on investigation experimental accurate An aytimes many . 3 ± opes nvriai iMneS neo i iti,I- Cinthia, Via Angelo, S. Monte di Universitario Complesso 0 . iatmnod cez iih,Uiestad aoi“Fe Napoli Universit`a di Fisiche, Scienze di Dipartimento 6 o ,rsetvl.Tels rniin curn ihannn non a with occurring transition, last The respectively. C, ntal o ae freezes water hot initially , siuoNzoaed iiaNcer,Szoed Napoli di Sezione Nucleare, Fisica di Nazionale Istituto .Esposito S. ∗ .D ii n .Somma L. and Risi, De R. , hti nw bu the about known is what the h aeo ae ht ttepressure just the example, at for that, is, water it of occurs; case transition fulfilled, the phase not more a are or that conditions so two stability into when them. happens breaks between This equilibrium it which mutual in but in parts state homogeneous, homogeneous a not in exist is the may it and system physical volume a the energy), example (for quantities modynamic temperature eoe agr n h eprtr ftesse rises system the speed of 0 average temperature the to the that and so dis- larger, water, Maxwell becomes liquid given the the for are to tribution molecules attaining these bond those form, of from to energy removed begins and potential ice cluster When into that energies formation. join kinetic to their able put are only molecules furthermore, moving and, slow nucleus initial the veloci- arrange lower water, can statistically ties pure with liquid in the fact, in In molecules only freezing. effective its before quired rcse r ohfvrda eprtrslwrthan lower temperatures growth at crystal favored 0 and non-ordered both Nucleation the are liquid. with the processes of the nucleus means in the by molecules ice, of loca- of characteristic forces structure the interaction crystalline in the of molecules starting tions further such hap- attract Second, may a state, to form randomly. liquid to the ice liquid only in order the that pen this, well-defined of order and a molecules crystal in in some minimum fact, all arrange of In should first water form, start. con- to sufficient to begins a freezing not is for this lower, dition or coexist) phases solid em ftoejiigadtoelaigtefre crys- formed in the leaving equal surface). those is tal and joining molecules those of of exchange terms continuing the where t oeue erae u,ee ftetmeauegoes 0 temperature the to if even down but, decreases molecules its tt fpeaiu qiiru wtri nametastable a in is a (water in equilibrium exists precarious water of state supercooled although process, freezing o ,s htsproln flqi ae sgnrlyre- generally is water liquid of supercooling that so C, ngnrli skonta,frgvnvle ftether- the of values given for that, known is it general In hssproln is, supercooling Thus hnlqi ae scoe,teaeaevlcte of velocities average the cooled, is water liquid When cooling o ovosy h eprtr ssta h value the at set is temperature the (obviously, C n. in,wihaeitrrtdi em of terms in interpreted are which tions, aigpaea 6 at place taking s daai oln fwtrarelevant a water of cooling adiabatic etlrslsahee,athorough a achieved, results mental o one. tefie eprtr hr iudand liquid where temperature fixed (the C T ta s h reigo initially of freezing the is, (that f0 of 02 aoi Italy Napoli, 80126 eioI”and II” derico o ,bcmsunstable. becomes C, giil rbblt of probability egligible fwtrbfr freezing before water of efacto de ± freezing 1 o ,3 C, e nrdeti the in ingredient key a , . 5 rcs,rte than rather process, ± 0 p . 5 f1amadat and atm 1 of o C nucleus has 2 state). Minor perturbations such as impurities or other V = 20 cm3 V = 50 cm3 V = 65 cm3 V = 80 cm3 can trigger the sudden appearance of the stable crys- PSC 1 0.28 0 0.46 o o o o talline phase for the whole liquid mass, again with the Tc = −8 C Tc = −14 C Tc = −22 C Tc = −26 C release of the entire crystallization heat (melting heat) PSC 0.75 0.50 0 0.11 which increases the temperature of the freezing liquid to the normal 0oC one. TABLE I: Probabilities for the occurrence of supercooling for V In general, when a system is in a metastable state, different volumes of the sample and different temperatures of the cryostat Tc. sooner or later it will pass to another stable state. In water, density and entropy fluctuations favor the forma- tion of crystallization nuclei but, if the liquid constitutes a stable state, such nuclei are always unstable and will liquid instead of pure water, or by fluctuations of the ex- disappear with time being. However, since the fluctua- ternal pressure or temperature, caused in the commercial tions become more pronounced the lower the tempera- refrigerator itself. This explains why no appreciable su- ture, if water is supercooled, for sufficiently large nuclei percooling is observed in normal situations. Obviously, they will result to be stable and grew with time, becom- the most direct way to induce freezing in supercooled ing freezing centers. The starting of the phase transition water is to introduce an external body in it, in order to is thus determined by the probability of appearance of directly lower the surface tension. those nuclei, and the reported Mpemba effect could be We have thus performed an accurate experimental in- simply a manifestation of this process. vestigation, accounting for a total of about one hundred We have calculated just this probability P as function runs, aimed to clarify the phenomenology of the Mpemba of the absolute temperature T of the metastable phase effect and its interpretation. In the first part of our ex- (the one at which the nucleus is in equilibrium with the periments we have tested all the above qualitative pre- liquid), obtaining the following result [11]: dictions about supercooling, by studying the cooling and freezing of tens of cm3 of normal water in a commer- 2 α − T∗ cial refrigerator, in daily operation conditions. The key P = exp β 2 . (1) T∗ (T − T∗) point, in fact, is not to obtain the most favorable physical conditions, employing sophisticated setups, but rather Here T∗ is the equilibrium temperature of the liquid-solid to reproduce the Mpemba conditions, that is adiabatic phase, α is a dimensionless normalization factor and β is cooling (with commercial refrigerators) of not extremely a constant whose expression is given by small quantities of water. We have used an Onofri re- 16πτ 3v2 frigerator for the cooling of double distilled water and a β = 2 , (2) NiCr-Ni thermocouple as a temperature sensor (Leybold 3Q kT∗ 666193), interfaced with a Cassy Lab software for data where τ is the surface tension, v the molecular volume acquisition. of the crystallization nucleus, Q the molecular heat of For fixed temperatures of the cryostat we have indeed the transition from the metastable phase to the nucleus observed supercooling in our samples, with the freezing phase, and k is the Boltzmann constant. Just to give occurring just along the lines predicted above. In par- an idea of the macroscopic value of the constant β, let 3 2 3 ticular, during the supercooling phase we have induced us note that τ v = Wsurf is the cube of the work done a number of small perturbations in our samples, namely, ∼ by the surface forces, and by assuming that Q kT∗ we variations of external pressure or temperature, mechan- may write: ical perturbations or introduction of an external macro- 3 scopic body (a glass thermometer held at the same tem- 16π Wsurf β ≃ , (3) perature of the sample). In all these cases we have regis- 3 Q tered the sudden interruption of the supercooling phase
surf and a practically instantaneous increase of the temper- that is the constant β is ruled by the ratio W /Q. o The probability P has a minimum at the liquid-solid ature to the value of 0 C, denoting the starting of the equilibrium temperature T∗ and increases for decreasing freezing process. Conversely, if no perturbation is in- temperature, as expected. From the formulae above it duced (or takes places) the water reached an equilibrium − o is clear that the probability for nucleation, and thus the with the cryostat at temperatures up to about 30 C onset of the freezing process as well, is enhanced if the (lasting also for several thousands of seconds). work done by the surface forces (or the surface tension We have then verified that when the freezing process itself) is lowered in some way. In normal daily conditions started from the supercooling phase, the Mpemba effect when a commercial refrigerator is employed, this is eas- took place with a probability in agreement with that re- ily induced in two simple ways: either by the presence ported in Ref. [3]. of impurities, when solutions (such as an ice cream so- In about half (with a total probability of 0.47) of the lution, as in the Mpemba case) are used as the freezing runs performed we have detected a supercooling phase. 3
o Tc = −8 ± 2 C V = 20 cm3 V = 50 cm3 V = 65 cm3 V = 80 cm3 ∆t1 (s) 7 ± 1 ∆t2 (s) 11 ± 6 220 ± 100 500 ± 170 630 ± 160 ∆t3 (s) 12 ± 6 70 ± 30 o Tc = −14 ± 2 C V = 20 cm3 V = 50 cm3 V = 65 cm3 V = 80 cm3 ∆t1 (s) 37 ± 1 ∆t2 (s) 8 ± 3 130 ± 80 480 ± 160 500 ± 60 ∆t3 (s) 7 ± 4 o 3 o Tc = −22 ± 1 C V Tc − FIG. 1: Left: Cooling curves for = 20 cm and = 8 C. V 3 V 3 V 3 V 3 Right: The fitted time duration of the phase transition at = 20 cm = 50 cm = 65 cm = 80 cm o t ± ± 3.5 C as function of the volume V of the samples, for different ∆ 1 (s) 63 1 7 1 t ± ± temperatures Tc of the cryostat. ∆ 2 (s) 170 100 130 70 ∆t3 (s) o Tc = −26 ± 1 C In Table I we report the observed probability PSC for V = 20 cm3 V = 50 cm3 V = 65 cm3 V = 80 cm3 the occurrence of supercooling for different volumes V of ∆t1 (s) 3.5 ± 0.7 the water sample and for different temperatures Tc of the ∆t2 (s) 3 ± 1 320 ± 70 210 ± 170 cryostat. We find the data to be fitted by a straight line, ∆t3 (s) 200 ± 70 1 ± 1 denoting (in the range considered) a linear decreasing TABLE II: Time duration of the phase transitions at 6oC of PSC for decreasing temperatures of the cryostat and o o for increasing volumes of the samples, this probability (∆t1), 3.5 C (∆t2) and 1.3 C (∆t3) for different volumes V o of the sample and different temperatures of the cryostat Tc. reaching the maximum PSC = 1 for Tc =0 C (and V = 0). An interesting feature of what we have observed is the for the temperature as function of time t is obtained: sensible appearance of iced water in our samples. In fact, when supercooling did not occur, the ice started to form −t/τ T = T − (T − T0)e , (5) around the walls of the beaker, while the inner parts were c c still in a liquid form, as usually expected. Instead the im- where τ = C/δ is a time constant measuring the cool- mediate freezing of supercooled water involved the whole ing rate of the sample. However, although the overall sample, this showing a very peculiar symmetric form. dependence of T on time is that expressed by Eq. (5), We have used cylindrical beakers with the temperature our experimental data clearly reveal the presence of three sensor in their periphery, near the walls; the observed transition points before freezing (or supercooling), where structure was a pure radial (planar) one, with no liquid τ changes its value. This transitions occur at tempera- water and radial filaments of ice from the center of the o o o tures T1 =6±1 C, T2 =3.5±0.5 C and T3 =1.3±0.6 C beakers to the walls (in one case we have been also able with a probability of P1 =0.11, P2 =0.84 and P3 =0.21, to take a low resolution picture of this, before its destruc- respectively. The time duration ∆t of each phase tran- tion outside the refrigerator). sition, during which the temperature keeps practically However, although supercooling plays a relevant role in constant [12], depends on the volume of the sample and the manifestation of the Mpemba effect, the things are on the temperature of the cryostat. The data we have made more complicated by the occurrence of other statis- collected are summarized in Table II. For the phase tical effects before the temperature of the water reaches transition at T2 these data show a linear increase of the value of 0oC. This comes out when an accurate mea- ∆t2 with Tc and a quadratic one with V ; in Fig. 1 we surement of the cooling curves is performed (some exam- give the fitting curves corresponding to best fit function ples of what we have obtained during the second part of 2 ∆t2 = (a+bTc)V . Instead, for the other two phase tran- our experiments are reported in Fig. 1). sitions no sufficient data are available in order to draw According to a simple naive model, the heat exchange any definite conclusion on the dependence on V and Tc, from the water sample (at initial temperature T0) to the though ∆t1 and ∆t3 appear to be shorter than ∆t2. cryostat (at fixed temperature Tc) is described by the The occurrence of these phase transitions is likely re- equation lated to the formation of more or less ordered structures in water, resulting from the competition between long- C dT = δ (Tc − T0) dt, (4) range density ordering and local bond ordering maxi- where C and δ are the thermal capacity and the heat mizing the number of local bonds [8]. The anomalous conductivity of the water, respectively. Thus by solving density maximum at about 4oC (which we observe here o the differential equation in (4), the following expression at T2 = 3.5 ± 0.5 C) is, for example, explained just in 4
o o o o Tc = −8 C Tc = −14 C Tc = −22 C Tc = −26 C tially hot water before cold one. τ1 (s) 600 ± 110 680 ± 100 1000 ± 110 950 ± 190 From the data we have collected we have verified that, τ ± ± ± ± 2 (s) 1080 260 1060 170 530 90 570 3 for given V and Tc, in many cases no inversion between τ3 (s) 1590 ± 930 1520 ± 730 270 ± 50 220 ± 80 the cooling curves happens before the freezing point, ir- τ4 (s) 620 ± 480 500 ± 180 150 ± 30 640 ± 490 respective of the change in the value of τ or the time duration of the phase transitions. Nevertheless we have TABLE III: Time constants τ1 (T