22 November 1999
Physics Letters A 263 Ž1999. 70–77 www.elsevier.nlrlocaterphysleta
A new model of periodic precipitation incorporating nucleation, growth and ripening Maurice Chacron, Ivan L’Heureux
)
Ottawa-Carleton Institute for Physics, Department of Physics, UniÕersity of Ottawa, Ottawa, Ontario, Canada K1N 6N5 Received 27 May 1997; received in revised form 28 September 1999; accepted 28 September 1999 Communicated by C.R. Doering
Abstract We present and solve a one-dimensional model of periodic precipitation which includes nucleation, growth and ripening processes. This model thus generalizes two important models: the prenucleation-based model of Dee and the postnucleation competitive growth model ŽCGM. of Feeney et al. By tuning a simple phenomenological parameter, our model smoothly bridges the gap between a nucleation-growth dominated regime and one where ripening is active. q 1999 Published by Elsevier Science B.V. All rights reserved. PACS: 36.40.Sx; 47.54.q r; 81.10.Aj Keywords: Periodic precipitation; Liesegang bands; Competitive growth model; Prenucleation model; Postnucleation model
1. Introduction Spatially periodic precipitation phenomena have generated much interest w1–3x ever since their first investigation by Liesegang at the end of last century. In a typical Liesegang experiment, a salt solution diffuses through an aqueous gel and reacts with another salt. According to the system geometry, the product then precipitates by forming a series of concentric rings or bands parallel to the diffusion front. Liesegang band formation have also been applied to explain various spatial patterns observed in geological systems w4–7x. The presence of the gel is not necessary as Liesegang bands have been observed in water-filled capillaries or in the gas phase ) Corresponding author. Tel.: q1-613-562-5800r6770; fax: q1-613-562-5190; e-mail:
[email protected]
w8x. In many cases, the bands obey the Jablczinsky’s spacing law w9x: X nq 1 s Pn X n
Ž 1.
and the time scaling law w10x: Tn s Q n X n2 .
Ž 2.
Here, X n , Tn denote the position and time of formation of the nth band, and Pn , Q n tend to constants as n increases. Models for periodic precipitation can be classified in two broad categories. Ža. In prenucleation models w11–15x, a feedback mechanism exists between the nucleation kinetics and diffusional transport. Nucleation of the precipitate particles takes place when a saturation threshold is reached. Since the process is diffusion-limited,
0375-9601r99r$ - see front matter q 1999 Published by Elsevier Science B.V. All rights reserved. PII: S 0 3 7 5 - 9 6 0 1 Ž 9 9 . 0 0 7 0 9 - 4
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
precipitation results in local depletion of reactant concentration and the nucleation process is stopped. As the reaction front moves further away, a new saturation threshold is eventually reached, nucleation of precipitate is triggered and a new band is formed. The cycle repeats itself in this fashion. The models in Refs. w12–15x are based on the classical atomistic theory of nucleation. Other simpler models w3,16–19x use diffusional transport coupled with a phenomenological macroscopic expression for the saturation threshold. Cellular automata reduction of such models are also available w20x. Although the models of Refs. w14,15x incorporate growth kinetics, they are considered w21x pre-nucleation models, since the nucleation phase is essential in determining the basic features of the resulting precipitation pattern. Prenucleation models generate bands which are consistent with the scaling laws but do not explain other often observed features, such as pattern formation in absence of an initial concentration gradient, the time-evolution of the bands and the existence of revert spacing or other irregular patterns. Žb. In postnucleation models w4,22,23x, a feedback between growth, diffusion and surface tension effects generates the pattern through a variant of the Lifshitz–Slyozov instability w24x. Here, the bands evolve through ripening after the nucleation phase is over. Homogeneous concentration profiles are unstable in such systems. Consider for example a small localized increase in particle size. This generates a local depletion of reactant. Diffusive transport then results in a net mass current towards the perturbation which increases its size. As the solution component precipitates, the particle in the zone adjacent to the perturbation will start to dissolve, thus generating an enhancement in the local concentration. This causes a net mass flux further away from the initial perturbation and growth of particle size there. The model has thus the potential to generate a succession of bands characterized by an oscillation in the particle size. Except for some general comments in Ref. w19x, a model that generalizes the two previous classes of models and reduces to them in the appropriate limits is not available. It is the purpose of this paper to present and solve such a model, in which nucleation, growth and ripening processes are present. We will specifically generalize Dee’s prenucleation model w14x and the postnucleation competitive growth model
71
ŽCGM. of Ref. w23x. We find that bridging the gap between a nucleation-growth dominated regime and one where ripening is important can be achieved by tuning a phenomenological parameter D, which characterizes the thickness of the Gibbs surface of a precipitate particle Žwhich in turn is of the order of the critical nucleation radius.. In contrast to what is assumed in many applications of the CGM, our model indicates that the number density of precipitate particles is not uniform.
2. The model We consider a one-dimensional finite reactor of length l along the spatial dimension x. Initially, the aqueous solution A Že.g., HCl. occupies the space x - 0 with a uniform concentration A 0 , whereas the solution B Že.g., AgNO 3 . is localized in the reactor 0 - x - l with a uniform concentration B0 . At x s 0, A and B react to form the product C in solution Že.g., AgCl., which can then precipitate. If A 0 4 B0 , the reaction front propagates in the x ) 0 direction. Let a, b and c be the concentration Žmoles per unit volume. of the reactants A and B and of the product C, respectively. Let D 1 , D 2 and D 3 denote the diffusion constant of A, B and C, respectively and k the reaction rate coefficient for the reaction A q B ™ C, considered irreversible. Then, the following diffusion-reaction equations describe the concentrations dynamics:
Ea Et Eb Et Ec
s D1 s D2 s D3
E 2a E x2 E 2b E x2 E 2c
y kab, y kab,
q kab y u. Ž 3. Et E x2 Here t is time and u is the precipitation rate, i.e., the number of moles of C precipitated per unit volume per unit time. In the spirit of the classical nucleation mechanism w25x, we assume that the precipitate particles are spherical and immobile. Let J Ž x,t . be the nucleation rate at position x and time t, r Ž x,t,t X . be the radius of a precipitate particle at time t, given that it was nucleated at time tX , and let Õ denote the molar
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
72
volume of precipitate. Then, the molar concentration of precipitate is given by w14x: f s Ž 4 pr3Õ .
t
X
H0 J Ž x ,t . r
3
Ž x ,t ,tX . d tX .
Ž 4.
The precipitation rate is thus: Ef us Et s Ž 4 pr3Õ . J Ž x ,t . rn3 Ž x ,t . q Ž 4 prÕ . Er t = J Ž x ,tX . r 2 Ž x ,t ,tX . Ž x ,t ,tX . d tX . Ž 5. E t 0 The first term on the right hand side of Eq. Ž5. represents the contribution due to the instantaneous nucleation process, r Ž x,t,t . being the radius of the critical particle rnŽ x,t .. The second term represents the growth of particles which nucleated in the past. We introduce the supersaturation s through: s ' Ž N0 c y C0 . rC0 Ž 6. where C0 is the concentration Žin molecules per unit volume. of the product C in equilibrium with precipitate particles of large radius and N0 is Avogadro’s number. Using classical nucleation theory, Dee w14x obtained J s Jc F Ž s . ,
H
rn s wg Ž s .
Ž 7.
where d is a typical molecular size Žtaken as twice the diameter of the C molecule.. Note that our definitnion of Jc is different from the one reported by Dee. We now generalize this nucleation mechanism by including growth and ripening dynamics in the model. This is accomplished through a nonlinear coupling in the particle growth between the particle radius and the supersaturation, as in post-nucleation models. Although other growth laws could also be considered w15x, we take the growth as interface-controlled. One has:
Er Et
sG
N0 c y Ceq Ž r .
where G is a kinetic coefficient and the Gibbs– Thomson relation: Ceq Ž r . s C0 exp wc Ž r . rr
F Ž s. '
½
Ž 1 q s . exp Ž y b g Ž s .
2
c Ž r . s exp y
0, s - 0 g Ž s . ' 1rln Ž 1 q s . . Ž 8. A similar expression for the precipitation rate is used in Ref. w15x. In these relations, the capillary length w is w22x 2 Õs ws Ž 9. pN0 k B T where s is the surface tension Žassumed independent of the particle radius here., p is a stoichiometric coefficient Žequal to 2 in our case., T is the temperature and k B is Boltzmann’s constant. The constants Jc and b are Jc s 4 p D 3 w 2 C02rd ,
bs
ž
4 ps w 2 3k B T
1r2
/
drr
H0
. , sG0
Ž 10 .
Ž 12 .
gives the molecular concentration of the product C in equilibrium with a precipitate particle of radius r. The dimensionless function c Ž r . expresses the radius dependence of the surface tension, which has an important effect on the growth of particles of small radius, but is neglected in the nucleation term ŽEq. Ž8... Koenig w26x has derived a thermodynamicallybased expression for this function:
where 2
Ž 11 .
C0
dx
2 Ž 1 q x q x 2r3 . 1 q 2 x Ž 1 q x q x 2r3 .
Ž 13 . where d f O Ž w . is a parameter characterizing the thickness of the Gibbs surface. For computational purposes, it is convenient to approximate c by a rational function having the correct asymptotic behavior for large and small r. One finds to a good approximation Žwith a maximum relative error of 0.8%.:
c Ž r. (
r2qd r r 2 q 3 d r q d 2rq
Ž 14 .
with q s 0.304359. Our expression for the equilibrium concentration is new and different from the one used in the CGM. Following Le Van and Ross w15x, it is convenient to transform the integro-differential Eq. Ž5. to a set of coupled differential equations by defining the
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
average surface of the precipitate particles per unit volume S :
S'4p
t
X
H0 J Ž x ,t ,t . r
2
Ž x ,t ,tX . d tX ,
Ž 15 .
the average radius of the particles per unit volume r :
r'
t
X
X
H0 J Ž x ,t ,t . r Ž x ,t ,t . d t
X
Ž 16 .
and the average particle number density N: N'
t
X
X
H0 J Ž x ,t ,t . d t .
Ž 17 .
We nondimensionalize the variables according to: w 1r2 t s t ) , x s Ž D 1wrG . x ) , a s A 0 a ) , G ÕC0 ) ÕC0 b s A0 b ) , S s S , rs r), N0 w 4 p N0 w 2 Ns
ÕC0 4 p N0 w
3
N ) , r s wr ) ,
Ž 18 .
and we introduce the six dimensionless parameters: K ' kA 0 wrG, K X ' KA 0 N0rC0 , D ' D 2rD1 , DX ' D 3rD1 ,
a ' 4 p Jc w 4 N0rGÕC0 , D ' drw
Ž 19 .
together with the dimensionless system size L s Ž GrD1w .1r2 l . In terms of dimensionless variables, we thus obtain the following generalized set of seven coupled nonlinear differential equations Ždropping the ).:
E 2a
Ea Et Eb Et
s
E x2
sD
y Kab
E 2b E x2
y Kab
E 2a s DX 2 q K X ab y Ž 1 q s y e c r r . S Et Ex a y g 3 Ž s. F Ž s. 3 ES s a g 2 Ž s . F Ž s . q 2Ž 1 q s y ec r r . r Et Er s a g Ž s . F Ž s . q Ž 1 q s y ec r r . N Et
Ž 20a. Ž 20b.
EN Et Er Et
Ž 20d. Ž 20e.
Ž 20f.
s 1 q s y ec r r .
Ž 20g.
fs
4pN 3Õ
r 3 Ž x ,t .
Ž 21 .
where N is the particle number density, assumed constant and uniform. In most applications of the CGM, N acts as a simple external parameter. It is easily seen that, if the nucleation rate decays on a time short compared to the growth and ripening time scales of interest, and if we neglect the spatial dependence of J, then we have approximately from Eqs. Ž4. and Ž21.: `
N(
X
H0 J Ž t . d t
X
Ž 22 .
in agreement with Eq. Ž17.. With N constant in time, the coefficient a F Ž s . must be set equal to zero for consistency with our model Eq. Ž20f.. Also, one finds by integration of Eqs. Ž20d. and Ž20e. that S s Nr 2 . Scaling the variables a, b, r, t, x and N as in Eq. Ž18., the supersaturation evolution Eq. Ž20c. becomes:
Et
Ž 20c.
sa F Ž s.
c is as in Eq. Ž14. with D substituting d . Dee’s prenucleation model is obtained as a special case by formally setting c s 0. Our model also reduces to the CGM in the appropriate limit. Indeed, in the CGM, the nucleation phase is assumed terminated and only production of C by reaction, growth of the precipitate particles and ripening occur. In the CGM, the reaction-diffusion rate laws are still the same as in Eq. Ž3., except that the concentration of precipitate ŽEq. Ž4.. is replaced by
Es
Es
73
s DX
E 2a E x2
q K X ab y Ž 1 q s y e c r r . Nr 2 .
Ž 23 .
This is equivalent to using Eq. Ž21. in Eq. Ž5., as in the CGM. Thus, the dynamics of the post-nucleation CGM constitutes a special case of our generalized model ŽEqs. Ž20a., Ž20b., Ž20c., Ž20d., Ž20e. and Ž20f.. under the condition that N is constant and uniform. The numerical solution of our generalized model indeed shows that N becomes constant after a short time but is far from being spatially uniform.
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
74
3. Stability of a uniform sol In order to see if the model allows for coarsening, it is useful to consider the dynamics of a spatially uniform sol in the absence of a source term Žchemical reaction.. Strictly, the system Ž20. with K X s 0 has no uniform steady state since Eqs. Ž20f. and Ž20g. have no common roots. However, for F Ž s . small, N is almost constant during the times of interest. To find this uniform ‘pseudo-steady state’, we integrate Eq. Ž20c. and, with the use of Eq. Ž20g., we get: sqN
r3
sH Ž 24 . 3 where H is a constant of the motion. Combining
Fig. 2. Time evolution of an uniform initial condition in the absence of reaction exhibiting coarsening for D s 0.16. The initial condition was chosen as r s 0.25, sse c r r y1, N s1.0, S s Nr 2 , r s Nr. A small perturbation of amplitude 10y3 was applied on r at t s 0 and at x s 0.5. This perturbation resulted in a local increase of r Ž0.5,t . and a smooth decrease Ževentually to zero. in the neighborhood of the initial perturbation.
with Eq. Ž20g., we see that the pseudo-steady state is given by the roots of e
crr
qN
r3 3
s H q 1.
Ž 25 .
Depending on the value of H, there can be two, one or no steady-states. The linear stability of the steady state Ž ss , rs ., if it exists, can be investigated by linearizing the system about Ž ss , rs .. We assume perturbations of the form, s s ss q ´ 1exp Ž i kx q v t . , r s rs q ´ 2 exp Ž i kx q v t .
Ž 26 .
where ´ 1,2 are small amplitudes, k is a wave-vector and v is an eigenvalue. The solution for v is given by
v "s y " Fig. 1. Schematic representation of the dispersion curves v Ž k . of Eq. Ž27. describing the linear stability of the uniform pseudo-steady state for the system Ž20.: Ža. Eq Nrs2 - 0; and Žb. Eq Nrs2 ) 0.
1 2 1
Ž E q DX k 2 q Nrs2 .
(Ž E q D k q Nr . y 4 D Ek 2 X
2
2 2 s
X
2
Ž 27 .
where E ' Žde c rrd r .< r s . The dispersion relations v " Ž k . are schematized in Fig. 1Ža. and Fig. 1Žb. for
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
the cases where E q Nrs2 - 0 or ) 0, respectively. In the first case, the steady-state is unstable to homogeneous Ž k s 0. perturbations, whereas it is stable to homogeneous perturbations in the second case. However, in both situations, the steady state is unstable to spatially-dependent perturbations Ž k / 0.. Moreover, the maximum perturbation growth rate occurs for large k. No finite length-scale is selected by the instability, at least at this linear level of description. Thus, the system has a tendency to coarsen. These findings are analogous to the ones for the CGM. Fig. 2 illustrates the numerical solution of the full model ŽEqs. Ž20a., Ž20b., Ž20c., Ž20d., Ž20e. and Ž20f.. for an initial condition corresponding to the pseudo-steady state solution with rs s 0.25 superposed to a small local perturbation at x s 0.5. This situation corresponds to the case of Fig. 1Ža.. The value of the radius tends to increase at the position of the perturbation and to decrease eventually to zero in the neighborhood of the initial perturbation.
75
7.48 = 10 15 cmy3 ; k s 10 5 moly1 sy1 ; D1 s D 2 s D 3 s 10y5 cm2 sy1 ; G s 10y7 cmrs; s s 170 erg cmy2 ; p s 2; w s 10y7 cm; T s 300 K; d s 8.7 = 10y8 cm; and Õ s 25.29 cm3 mol.y1 . These values are typical of the system AgNO 3 q HCl ™ AgClq inactive products and, except for B0 , were taken from Ref. w14x. Fig. 3 illustrates the scaled precipitate molar concentration f N0rC0 ŽEq. Ž3.. and the scaled radius r
4. Results and discussion The initial conditions are: aŽ x s 0,0. s 1, aŽ x / 0,0. s 0, bŽ x,0. s B0rA 0 , sŽ x,0. s y1, S Ž x,0. s r Ž x,0. s N Ž x,0. s r Ž x,0. s 0, i.e., no precipitate exists at t s 0. The system dimensionless size L s 100 is chosen large enough with the following boundary conditions: aŽ0,t . s 1, aŽ L,t . s 0, bŽ L,t . s B0rA 0 and a zero flux condition at the origin for the components B and C: Eb Ž0,t .rE x s EsŽ0,t .rE x s 0. According to w27,28x, the parameters A 0 y B0 and C0 are important parameters that characterizes the type of precipitation pattern obtained. A complete analysis of our model in terms of these parameters is outside the scope of this work. Instead, we will fix the value of these parameters and vary the phenomenological parameter D. We have solved the set of parabolic partial differential Eqs. Ž20a., Ž20b., Ž20c., Ž20d., Ž20e. and Ž20f. by a non-iterative Crank–Nicholson method, with forward projection on the nonlinear terms w29x. The scheme was stable and convergent. We used dimensionless discrete time and space steps of 0.001 and 0.1, respectively. The parameter values are given by: A 0 s 10y5 mol cmy3 ; B0 s 10y6 mol cmy3 ; C0 s
Fig. 3. Ža. Scaled molar concentration of precipitate fN0 r C0 as a function of scaled position for D s 0.16 and for two different scaled times: t s 300 Žcontinuous line. and t s900 Ždotted line.; and Žb. scaled particle radius r for the same parameter values. One dimensionless unit in x and t corresponds to 3.16=10y3 cm and 1 s, respectively. The scale of f is in 1.24=10y8 mol cmy3 , whereas one dimensionless unit of r corresponds to 10y7 cm.
76
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
as a function of the scaled distance for two different times. The value D s 0.16 was chosen, but the time evolution turns out to be very nearly independent of D, as long as D G 0.156. Thus, the solution exhibits the formation of precipitate bands. The density function is similar to the one obtained by Dee w14x and it is easy to verify that our precipitation pattern is consistent with the spacing and time laws of Eqs. Ž1. and Ž2.. Note that the position for which the radius is a relative minimum corresponds to a maximum in the concentration, in accordance with the fact that many nuclei of small radius are nucleated. In fact, in the limit of large D, the ripening factor c Ž r .rr is negligible for all values of r and Dee’s prenucleation model is recovered.
Fig. 4. Same as Fig. 3 with D s 0.136 for t s 300 Žcontinuous line. and t s900 Ždashed line..
Fig. 5. Scaled particle number density for a scaled time t s 300 Žcontinuous line. and t s900 Ždashed line., and for D s 0.136. One dimensionless unit in N corresponds to 2.5=10 13 cmy3 . The density is very nearly constant in time. The times taken to reach the first four maxima Žat x s 0.0, 2.0, 4.7 and 8.2. are 0.71, 2.65, 6.70 and 14.18, respectively.
Fig. 4 illustrates the same quantities, but using a smaller value of DŽs 0.136.. The solution is quantitatively different. For the later stages of the process and for a position larger than some threshold, the bands continue to evolve after their formation: they create a doublet Žsatellites.. Thus, each band is eventually characterized by two narrow sub-bands of high-density aggregate, separated by a wider zone of low-density precipitate. During this phase of the process, a minimum in the radius tends to occur at the same position as a minimum in the precipitation, indicating that ripening is taking place. Another feature of the solution is to exhibit interbands regions where no precipitate exist. These characteristics are typically found in postnucleation models of periodic precipitation and are also observed in nature w4x. As the value of D is further decreased, the solution is similar to the case of Fig. 4 except that the doublet formation is more pronounced: the contrast between the high-density edges of the doublet and its lowdensity middle region is enhanced. A CGM solution with a typical, constant and uniform value of the particle number density N, also generates doublets in the molar concentration of precipitate. However, in contrast to our model where both ripening and nucleation take place, the solution
M. Chacron, I. L’Heureuxr Physics Letters A 263 (1999) 70–77
obtained from the CGM exhibits much higher and narrower peaks in the molar concentration of precipitate than those in Fig. 4. Moreover, our results ŽFig. 5. show that, whereas N is indeed nearly constant in time, it is far from uniform. A banding pattern is also reflected in the spatial dependence of N.
77
Acknowledgements This research was supported by grants from the Natural Sciences and Engineering Research Council of Canada. I.L. acknowledges fruitful discussions with A.D. Fowler. References
5. Conclusion In this paper, we have presented a periodic precipitation model which smoothly bridges the gap between a regime where nucleation and growth are dominant and one where ripening is important by varying a simple phenomenological parameter D. This parameter is a by-product of our expression for the radius-dependence of the surface tension and characterizes the thickness of the Gibbs surface. For large values of D, the system exhibits the formation of bands of precipitate mostly via nucleation and diffusional transport as in Dee’s prenucleation model. For small values of D, the ripening mechanism becomes important and further evolution of the bands is characterized by the formation of doublets, as in the CGM. One basic assumption underlying the applications of the CGM is that the precipitate number density is constant, uniform and a priori known. By direct computation, we have found that, although the number density is very nearly constant, it is far from uniform. This feature and the fact that nucleation is not neglected in our model explains why our solutions are quantitatively different from those obtained in the CGM. As in many periodic precipitation models, ours neglects fluctuations and thus constitutes a mean-field approach to the problem. A more complete analysis would require the calculation of the probability density of finding a crystal with radius between r and r q dr at position x and time t. This could be achieved by using our Gibbs–Thomson relation in existing molecular theories of coarsening in periodic precipitation patterns Žsee for example w30,31x. or by incorporating the precipitation reactions to existing theories of phase separation combining nucleation, growth and ripening w32x.
w1x R.E. Liesegang, Photogr. Archiv. 21 Ž1896. 221. w2x K.H. Stern, Chem. Rev. 54 Ž1954. 79. w3x H.K. Henisch, Periodic Precipitation, Pergamon Press, Oxford, 1991. w4x A.E. Boudreau, Min. Pet. 54 Ž1995. 55. w5x A.E. Boudreau, S. Afr. J. Geol. 97 Ž1994. 473. w6x P. Ortoleva, Geochemical Self-Organization, Oxford University Press, New York, 1994. w7x I. L’Heureux, A.D. Fowler, in: B. Jamtveit, P. Meakin ŽEds.., Growth, Dissolution and Pattern Formation in Geosystems, Kluwer Academic Publishers, Amsterdam, 1999. w8x A. Buki, M. Zrınyi, Physica A 220 ¨ E. Karpati-Smidroczki, ´ ´ ´ ´ Ž1995. 357. w9x K. Jablczynski, Bull. Soc. Chim. France 33 Ž1923. 1592. w10x H.W. Morce, G.W. Pierce, Physik. Chem. 45 Ž1903. 589. w11x W. Ostwald, Lehrbuch der Allgemeinem Chemie, Engelman, Leipzig, 1897. w12x S. Prager, J. Chem. Phys. 25 Ž1956. 279. w13x D.A. Smith, J. Chem. Phys. 81 Ž1984. 3102. w14x G.T. Dee, Phys. Rev. Lett. 57 Ž1986. 275. w15x M.E. Le Van, J. Ross, J. Phys. Chem. 91 Ž1987. 6300. w16x Y. Brechet, J.S. Kirkaldy, J. Chem. Phys. 90 Ž1989. 1499. w17x V. Talanquer, J. Chem. Educ. 71 Ž1994. 58. w18x A. Buki, M. Zrınyi, J. Chem. Phys. ¨ E. Karpati-Smidroczki, ´ ´ ´ ´ 103 Ž1995. 10387. w19x D.S. Chernavskii, A.A. Polezhaev, S.C. Muller, Physica D ¨ 54 Ž1991. 160. w20x B. Chopard, P. Luthi, Phys. Rev. Lett. 72 Ž1994. 1384. w21x R. Lovett, P. Ortoleva, J. Ross, J. Chem. Phys. 69 Ž1978. 947. w22x A.A. Polezhaev, S.C. Muller, Chaos 4 Ž1994. 631. ¨ w23x R. Feeney, S.L. Scmidt, P. Strickholm, J. Chadam, P. Ortoleva, J. Chem. Phys. 78 Ž1983. 1293. w24x I.M. Lifshitz, V.V. Slyozov, J. Phys. Chem. Solids 19 Ž1961. 35. w25x I.V. Markov, Crystal Growth for Beginners, World Scientific, Singapore, 1996. w26x F.O. Koenig, J. Chem. Phys. 18 Ž1950. 449. w27x S.C. Muller, S. Kai, J. Ross, J. Phys. Chem. 86 Ž1982. 4078. ¨ w28x S. Kai, S.C. Muller, J. Ross, J. Phys. Chem. 87 Ž1983. 806. ¨ w29x D.U. von Rosenberg, Methods for the Numerical Solution of Partial Differential Equations, Elsevier, New York, 1969. w30x G. Venzl, Phys. Rev. A 31 Ž1985. 3431. w31x G. Venzl, J. Chem. Phys. 85 Ž1986. 1996. w32x C. Sagui, M. Grant, Phys. Rev. E 59 Ž1999. 4175.