arxiv: v1 [cond-mat.mtrl-sci] 7 Mar 2013

Transcription

1 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine Vo. 00, No. 00, 00 Month 202, 2 RESEARCH ARTICLE arxiv: v [cond-mat.mtr-sci] 7 Mar 203 Studies of homogeneous precipitation in very diute iron-coper aoys using kinetic Monte Caro simuations and statistica theory of nuceation V.G. Vaks a,b, F. Soisson c, and I.A. Zhuravev a a Nationa Research Center Kurchatov Institute, 2382 Moscow, Russia; b Moscow Institute of Physics and Technoogy (State University), 7303 Moscow, Russia; c CEA, DEN, Service de Recherches de Métaurgie Physique, 99 Gif-sur-Yvette, France (Received Juy, 207) Abstract Kinetics of homogeneous nuceation and growth of copper precipitates under eectron irradiation of Fe x Cu x aoys at concentrations x from 0.06 to at.% and temperatures T from 290 to 450 C is studied using the kinetic Monte Caro (KMC) simuations and the statistica theory of nuceation (STN). The conventiona assumption about the simiarity of mechanisms of precipitation under eectron irradiation and during therma aging is adopted. The earier-deveoped ab initio mode of interactions in Fe-Cu aoys is used for both the KMC simuations and the STN cacuations. Vaues of the nuceation barrier F c and the prefactor J 0 in the Zedovich-Vomer reation for the nuceation rate J are cacuated for a number of concentrations and temperatures. For the diute aoys with x 0.2%, the STN and the KMC resuts for the nuceation barrier F c do virtuay coincide with each other, which seems to confirm a high reiabiity of the STN for this probem. The STN cacuations are aso used to estimate the temperature dependencies of concentrations which correspond to the homogeneous or the heterogeneous precipitation imit, x hom (T) and x het (T), and both dependencies are found to be rather sharp. Keywords: precipitation, diute Fe-Cu aoys, kinetic Monte Caro simuations, statistica theory of nuceation. Introduction Studies of copper precipitation in irradiated Fe-Cu-based aoys at temperatures T C and ow copper concentrations x between x= % and x= % (in atomic percents, here and beow) attract great attention, first of a in connection with the probem of hardening and embrittement of the nucear reactor pressure vesses (RPVs), see, e. g. []-[9]. Numerous experiments combined with some simuations showed that for the usua conditions of eectron irradiation, that is, when the radiation-induced custers of vacancies and sef-interstitia atoms are not formed and the heterogeneous precipitation on them is basicay absent, SANS and ATP measurements at x 0.% typicay do not revea any copper-rich precipitates [, 2, 6], even though the aoys remain to be strongy supersaturated. The absence of a significant homogeneous copper precipitation in these eectronirradiated Fe-Cu-based aoys is commony expained by the very ow voume fraction of probaby very sma precipitates,... and... one can infer that the imit of observabe copper precipitation at 290 C is ocated between x=0.09 and 0.08% Corresponding author. Emai: vaks@mbsab.kiae.ru ISSN: print/issn onine c 202 Tayor & Francis DOI: 0.080/ xx.xxxxxx

2 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 2 Tayor & Francis and I.T. Consutant [2], or between x=0. and 0.2% [3]. In particuar, custer dynamics simuations of homogeneous precipitation at x=0.088%, T = 290 C described in [6] predict a veryowdensity(5 0 6 m 3 )ofverybig(r 30nm)copperprecipitates. Under neutron or ion irradiations when copper precipitation can occur heterogeneousy at custers of vacancies or sef-interstitia atoms, such precipitation at T = 290 C was observed for x 0.08 [6, 8], but not for x 0.06% [5, 7]. Therefore, theoretica studies of copper precipitation in very diute Fe-Cu-based aoys near the concentration or temperature imits of such precipitation, x hom (T) for the homogeneous nuceation andx het (T) for the heterogeneous nuceation, seem to be interesting from both fundamenta and appied standpoints, even for binary Fe-Cu aoys which can be considered as mode aoys for the RPV stees [, 2, 4, 6]. Such theoretica studies can be particuary interesting if they can sufficienty reiaby predict the temperature dependencies of these precipitation imits at eevated T > 290 C for which experimenta estimates are either uncertain or absent. A consistent ab initio mode for studies of precipitation kinetics in Fe-Cu aoys has been recenty deveoped by Soisson and Fu [0], and detaied kinetic Monte Caro (KMC) simuations of precipitation in a number of Fe-Cu aoys at different T and x based on this mode yied a good agreement with avaiabe experimenta data [0, ]. Therefore, this mode and KMC methods deveoped in [0] seem to be prospective to study the above-discussed probem of copper precipitation in diute Fe-Cu aoys. We wi aso use the statistica theory of nuceation and growth of isoated precipitates (to be abbreviated as STN) deveoped by Dobretsov and Vaks [2, 3]. This theory provides microscopic expressions for both the nuceation barrier (activation barrier for the formation of a critica embryo) F c and the prefactor J 0 in the phenomenoogica Zedovich and Vomer formua for the nuceation rate J (the number of supercritica embryos formed in unit voume per unit time) [4] vaid for the initia, steady-state nuceation stage of precipitation: J = J 0 exp( βf c ) () where β = /T is the reciproca temperature. The STN corresponds to a number of refinements of earier modes of nuceation suggested by Cahn-Hiiard [5] and Langer [6]. It is based on the generaized Gibbs distribution approach [7] and enabes one to quantitativey cacuate vaues of F c and J 0 for the microscopic aoy mode chosen, particuary in the case of high nuceation barriers, βf c, which just corresponds to the diute aoys under consideration. In description of growth of supercritica embryos we aso use the phenomenoogica equation of such growth suggested in the cassica theory of nuceation [4] and access vaidity of this equation by comparison with our KMC simuations. In a our simuations and cacuations we consider ony homogeneous precipitation and accept the conventiona assumption (discussed, in particuar, by Mathon et a. [2]) that for the moderate irradiation intensity (i.e. when the ratio of baistic jump frequency to thermay activated jump frequency is sma), the main effect of irradiation is to increase the point defect concentrations, and hence the Fe and Cu diffusion coefficients. It means that the precipitation kinetics under eectron irradiation is the same as that during therma aging at the same temperature, except for an acceeration factor A irr defined by the reation: where D th Cu and Dirr Cu A irr = D irr Cu/D th Cu (2) are the copper diffusivities during therma aging and under

3 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 3 irradiation, respectivey. As evoution of microstructure under precipitation is determined by the diffusion of copper, Eq. (2) impies an anaogous reation between the evoution times t irr and t th under irradiation and during therma aging: t irr = t th /A irr. (3) Then the precipitation kinetics under eectron irradiation can be described by the KMC codes deveoped by Soisson and Fu [0] for therma aging with repacing the therma aging timet th by the scaed time t irr accordingto Eq. (3). Thisapproach is accepted throughout this work. In more detai, scaing reation (2) and estimates of the acceeration factor A irr are discussed in Sec. 4. At the same time, the main resuts and concusions of this work do not depend on A irr vaues as this scaing factor is canceed in many important kinetic characteristics. In Sec. 2 we present some resuts of cacuations of structure of critica embryos and parameters F c and J 0 in () for diute Fe-Cu aoys obtained using the STN [2, 3]. In Sec. 3 we describe our KMC simuations of nuceation and growth of copper precipitates in diute Fe-Cu aoys. In Sec. 4 we present an estimate of the acceeration factor A irr based on the rate theory modes [22] and find this estimate to reasonaby agree with those obtained using the avaiabe experimenta data [2] and our KMC simuations. In Sec. 5 we use the phenomenoogica equation for growth of a supercritica embryo [4] and our KMC simuations to estimate some important kinetic parameter which enters into the STN expression for the prefactor J 0 in Eq. (). In Sec. 6 we show that the nuceation barrier F c vaues cacuated using the STN agree we with those estimated in our KMC simuations, particuary at arge βf c which correspond to the concentrations and temperatures near the above-mentioned precipitation imits. In Sec. 7 we use the STN cacuations combined with some pausibe physica assumptions to estimate the temperature dependence of both the homogeneous and the heterogeneous precipitation imit, x hom (T) and x het (T), for practicay interesting temperatures T between 290 and 450 C. Our main concusions are summarized in Sec Properties of critica embryos in diute Fe-Cu aoys cacuated using the statistica theory of nuceation In this section we use the STN described in Refs. [2, 3] (to be referred to as I and II) to cacuate some thermodynamic, structura and kinetic characteristics of critica embryos. For these cacuations we use the ab initio mode of Fe-Cu aoys deveoped by Soisson and Fu and described in detai in Ref. [0]. Here we ony note that this mode uses the foowing vaues of the binding energy between two copper atoms and between a copper atom and a vacancy, E bn neighbors (in ev): E b CuCu = T, E b vcu = 0.26, CuCu and Ebn vcu E b2 CuCu = T,, for the n-th E b2 vcu = (4) The high vaues of ECuCu bn in Eqs. (4) correspond to strong thermodynamic driving forces for precipitation, whie a strong attraction between vacancy and copper atoms eads to the strong trapping of vacancies by copper precipitates discussed in detai in [0]. First we consider the structura and thermodynamic characteristics of embryos which are described in the STN in terms of the occupation number n i for each

4 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 4 Tayor & Francis and I.T. Consutant attice site i. The operator n i is unity when a copper atom is at site i and zero otherwise, whie probabiities of various distributions {n i } are described by the distribution function P{n i } given by Eq. (I-3): P{n i } = exp[β(ω+µ i n i H)]. (5) Here H is the configurationa Hamitonian supposed to be pairwise, µ is the chemica potentia, and the grand canonica potentia Ω is determined by the normaization condition: H = i>j v ij n i n j, Ω = T ntrexp[β( i λ i n i H)] (6) where v ij are the configurationa interactions, whie symbo Tr means summation over a configurations {n i }. The free energy F of an aoy with embryo is determined by reations F = Ω+µ i c i, µ = F/ c i = constant (7) where c i = n i = Tr(n i P) is the mean occupation number, or the oca concentration. The free energy F = F{c i } in (7) was cacuated using severa statistica approximations described in detai in I and II: mean-fied, mean-fied with fuctuations, pair-custer, and pair-custer-with fuctuations. The atter was found to be more consistent and accurate than others, thus beow we use ony the paircuster-with fuctuations approximation. In tabes and 2 and Fig. we present some characteristics of critica embryos. To iustrate both concentration and temperature dependencies of these characteristics, we consider two series of Fe x Cu x aoys: those with the same temperature T = 290 C but different concentrations x varying from 0.06 to % (Tabe ), and those with the same concentration x = 0.2% but different temperatures T varying from 290 to 390 C (Tabe 2). In these tabes, N c is the tota number of copper atoms within the embryo, and R c is its effective radius defined as that of the sphere having the same voume as N c copper atoms in the BCC attice of α-iron with the attice constant a 0 =0.288 nm: R c = a 0 (3N c /8π) /3 = 0.42N /3 c nm. (8) Note that this critica radius (as we as the precipitate radius R in Eq. (27) beow) characterizes the tota number of atoms but not the geometrica size of precipitate, in particuar, not its mean squared radius r 2 i /2 defined by Eq. (II- 0). Therefore, for sma and oose precipitates, such as those shown by ines B and C in Fig., these radii can notaby differ from geometrica ones. However, quantities R c and R in (8) and (27) are convenient to describe the precipitation kinetics. Quantities Ω 0 and Ω in tabes and 2 are the zero-order and the first-order terms in fuctuative contributions to the tota nuceation barrier F c = Ω 0 + Ω (defined as the difference between grand canonica potentias of two aoy states, that with the embryo and that without the embryo). These two terms correspond to the iterative treatment of some exact reation of thermodynamics of nonuniform systems given by Eq. (I-9) whichreates thefreeenergy F = F{c i }to thecorreator K ij of fuctuations of site occupations n i : K ij = (n i c i )(n j c j ). As discussed

5 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 5 in II and beow, the accuracy of such iterative treatment of fuctuations for diute aoys considered is usuay rather high. Fig. iustrates variations of structure of critica embryos with concentration or temperature. The resuts presented in tabes, 2 and Fig. ceary iustrate very sharp variations of both structure and thermodynamics of critica embryos with concentration x and temperature T. In particuar, at T = 290 C, the decrease of concentration from x = to x = 0.06% eads to the increase of the embryo size N c and the reduced nuceation barrier βf c by about three times: from N c 8 to N c 25, and from βf c 2 to βf c 35. According to the Zedovich-Vomer reation (), the atter impies decreasing the nuceation rate by about ten orders of magnitude (variation with x of the prefactor J 0 in () wi be shown to be negigibe). Tabe 2 iustrates a simiar sharp increase of N c and βf c under eevating temperature between 290 and 390 C. Fig. shows that this sharp rise of sizes and nuceation barriers is accompanied by notabe changes in the structure of the embryo: its boundary at high βf c 30 is much ess diffuse than that at moderate βf c 20. Discussing the concentration profies shown in Fig. we note that these profies (as we as those shown in Fig. in [] and Fig. 2 in [8]) usuay correspond to rather diffuse interfaces. At the same time, KMC simuations show that ony very sma custers can have a diffuse interface, whie for R 0.3 nm (that is, N 0 copper atoms) one observes amost pure copper custers with a sharp interface, see, e.g., Fig. 8 in [0]. However, one shoud take into account that the profies shown in Fig. (and other simiar figures) correspond to statistica averaging over a orientations of a custer which is typicay rather anisotropic, and this averaging eads to diffuseinterfaces; it is iustrated, in particuar, byfig. 8in [0]. Let usaso note that in phenomenoogica treatments of critica embryos based on the Cahn- Hiiard continuum approach [9, 20], the embryo interfaces are usuay notaby more diffuse than those obtained in the microscopic STN used. In particuar, rather sharp interfaces shown by curves A and D in Fig. can hardy be obtained in the phenomenoogica treatments. Let us now discuss the prefactor J 0 in Eq. (). In the STN, it is given by Eq. (II-3): J 0 = (β γ 0 /2π) /2 ND R (u)d aa. (9) Here the first three factors (discussed in detai in I) have the thermodynamic origin. They describe dependencies of the free energy and the distribution function of the embryo on the critica variabe a characterizing its size, and on three variabes u characterizing tits position R. The ast factor D aa is the generaized diffusivity in the c i -space that corresponds to an increase of a, i. e. to growth of the embryo. This factor is defined by Eq. (I-55) which expresses D aa as a certain inear combination of generaized mobiities M ij which describe the tempora evoution of mean occupations c i via Eq. (I-50): dc i /dt = j [M ij δ ij M ik ]β F/ c j. (0) The M ij vaues can be cacuated using some microscopic modes, for exampe, those empoyed in II. In particuar, for the idea soution (or in the diute aoy imit which corresponds to sma c i ), the mobiity M ij = Mij 0 is given by Eq. (II-8) with u ik = 0: k M 0 ij = δ ij,nnγ 0 nn [c i( c i )c j ( c j )] /2. ()

6 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 6 Tayor & Francis and I.T. Consutant Here δ ij,nn is the Kroneker symbo equa to unity when sites i and j are the nearest neighbors and zero otherwise, whie quantity γ 0 nn has the meaning of the mean rate of exchanges between neighboring copper and iron atoms. It is reated to the copper diffusivity D 0 in a pure iron by the foowing reation (see, e. g., Eqs. (72) and (74) in [23]): γ 0 nn = D 0 /a 2 0 (2) where a 0 is the BCC iron attice constant used in Eq. (8). At the same time, mobiities M ij can be considered as phenomenoogica parameters in the genera Onsager-type equations (0), and for weaky nonuniform aoys under consideration they can be supposed to be nonzero ony for neighboring sites i and j, just as in Eq. (). For our probem, the phenomenoogica mode of M ij can be considered as adequate if it propery describes the diffusion of soute (copper) atoms in processes of nuceation and growth of embryos. In Sec. 5 we show that growth of a supercritica embryo in diute Fe-Cu aoys can be described we using some generaization of the phenomenoogica diffusion equation suggested in the cassica theory of nuceation [4] given by Eq. (26) beow. This generaized diffusion equation differs from those for idea or diute soutions mainy by the vaue of the effective diffusivity D eff which takes into account the acceeration of precipitation kinetics due to the high mobiity of sma copper custers noted in Refs. [0, 2] and discussed beow. The D eff vaues for diute Fe-Cu aoys considered are estimated in Sec. 5 using kinetic Monte Caro simuations. Therefore, we suggest that the adequate phenomenoogica description of precipitation in Fe-Cu aoys under consideration can be obtained if we assume for mobiities M ij in (0) the same simpest form as that in Eq. (), M ij = δ ij,nn γ nn [c i ( c i )c j ( c j )] /2, (3) but assume the effective exchange rate γ nn in (3) to be proportiona not to the diffusivity D 0 in a pure iron as in Eq. (2), but to the phenomenoogica effective diffusivity D eff mentioned above: γ nn = D eff /a 2 0. (4) Vaidity of this assumption wi be checked by comparison with the KMC simuation resuts presented in Secs. 5 and 6. In Tabe 3 we present vaues of various factors in the expression (9) for the prefactor J 0 in Eq. () for some aoys Fe-Cu. Quantity γ 0 in this tabe is the derivative of the free energy F with respect to the effective embryo size a: γ 0 = 2 F/ a 2. As the critica embryo corresponds to the sadde-point of F in the c i -space with respect to this size (I, II), γ 0 vaues in tabe 3 are negative. N is the normaizing constant in the embryo size distribution function; D aa is the generaized diffusivity discussed above, and D R (u) is some geometrica factor which for a arge embryo is proportiona to its surface. For quantities D aa, D R (u), and the tota prefactor J 0 we present their reduced, dimensioness vaues: D aa /γ nn, v a D R (u), and J 0 = v a J 0 /γ nn. Note that the embryo size N c decreases with x and increases with T, as tabes and 2 show. Tabe 3 shows that both concentration and temperature dependencies of quantities N and D aa /γ nn are rather weak. On the contrary, the geometrica factor D R (u) notaby increases with the embryo size N c, whie β γ 0 somewhat decreases with N c. The tota reduced prefactor J 0 = v a J 0 /γ nn notaby increases with N c, by

7 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 7 about 20 or 5 times for the tota x or T intervas shown in tabe 3. However, in the tota nuceation rate J given by Eq. (), these variations of J 0 are negigibe as compared to theabove-mentioned huge changes of the activation factor exp( βf c ) for these intervas of x or T. Now we note that a very sharp decrease of this activation factor with x between x 0.2 and x 0.08% seen in tabe evidenty correates with the above-mentioned experimenta estimates of the homogeneous precipitation imit at T = 290 C given in Refs. [2, 3, 6]. It seems natura to suggest that the position x hom of this imit is mainy determined by the activation factor exp( βf c ), or by the reduced nuceation barrier βf c which at x = x hom reaches a certain high vaue (βf c ). For exampe, if we accept for the experimenta x hom =x at 290 C the estimate from Ref. [2], x exp %, we have: (βf c ) 28, whie taking this estimate from Ref. [3], x exp 0.5%, we obtain: (βf c ) 23. In more detai, estimates of the precipitation imits x are discussed beow in Sec Kinetic Monte Caro simuations of nuceation and growth of precipitates in diute Fe-Cu aoys As mentioned in Sec., our KMC simuations use the scaing assumptions (2) and (3) which enabe us to reate the homogeneous precipitation in irradiated aoys to that during therma aging for which the KMC codes deveoped by Soisson and Fu [0] can be empoyed. In these simuations we consider the foowing Fe x Cu x aoys: (A) those at T = 290 C and x equa to 0.088, 0.6, 0.44, 0.77, 0.2, 0.3 and %, and (B) those at x = 0.2% and temperatures T equa to 290, 305, 320, 335, 350, 365, 380 and 390 C. We use simuation voume V s = L 3 with periodic boundary conditions and, to check the statistica reiabiity of resuts, we usuay empoy severa different vaues of the simuation size L: L = 28a 0, L 2 = 60a 0, L 3 = 92a 0. (5) Before to describe the resuts, we mention a characteristic feature of such simuations for some ow concentrations x or high temperatures T: at such x and T, we observe ony one precipitate within simuation voume, that is, the tota number N p of precipitates within simuation box is unity (whie at further owering x or eevating T we observe no precipitates for the simuation time). In our simuations it was the case for the series (A) aoys with T = 290 C at x 0.44%, and for the series (B) aoys with x = 0.2%, at T 335 C.Ina suchcases, thereation N p = preserves under increasing simuation voume (repacing L by L 2 or L 3 ), though the incubation time t inc (that which precedes to the formation of the embryo [24]) somewhat decreases when V s increases. The presence of changes of physica characteristics (such as the precipitate density d p = N p /V s or the incubation time t inc ) under variations of simuation voume evidenty indicates on the statistica unreiabiity of simuation resuts for these characteristics. At the same time, the growth of a precipitate after its formation seems to be described by such simuations quite propery. Therefore, even though the simuations with N p = can not be used to study the evoution of precipitate density d p (t), they wi be widey used for studies of precipitate growth described in Sec. 5. When the number of precipitates within simuation box significanty exceeds unity, the simuations become statisticay reiabe and can be used to study the precipitate density d p (t). In figures 2 and 3 we show the dependencies d p (t) for a our simuations for which the fina number of precipitates within simuation box is

8 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 8 Tayor & Francis and I.T. Consutant not toosma: N p 3. Thesimuation time on abscissaaxis is given inthe therma aging vaues t th which in our mode are reated to the observed time t irr by Eq. (3). The acceeration factor A irr in this reation and comparison with the avaiabe data [2] about precipitation in an Fe-0.3%Cu aoy under eectron irradiation are discussed beow in Sec. 4. In the ower part of frames 2(e) and 2(f) we aso show tempora dependencies of the mean precipitate radius R m for two simuations with x = 0.3%. For other our simuations with N p 0 these dependencies are simiar and describe usuay a smooth increase of R m by about 20-30%. The dependencies d p (t th ) presented in Figs. 2 and 3 are approximatey inear. Therefore, if we define the nuceation rate J = J KMC for these simuations as the ratio J KMC = d p (τ)/τ (6) where τ = (t th t inc th ) is the evoution time counted off the incubation time tinc th, this nuceation rate for each simuation remains approximatey constant, particuary for the simuations with not sma N p where fuctuations are not too strong. Hence these our simuations seem to describe mainy the first stage of precipitation, the steady-state nuceation, for which the constant vaue of the nuceation rate J is characteristic. The dashed ine in each frame of Figs. 2 and 3 shows our estimate of the J KMC vaue in Eq. (6) assuming it to be constant. For this estimate we tried to use the maximum broad interva of time t th for which the tempora dependence d p (t th ) is cose toinear. Inspiteoftheevidentscatter ofkmcresutspresentedinfigs 2and 3, particuary for simuations with not arge N p, the inear tempora dependence d p (τ th ) can usuay be foowed rather ceary. Therefore, our estimates of J KMC vaues in Eq. (6) seem to be sufficienty definite. At the same time, the Zedovich-Vomer reation () is vaid just for the steadystate nuceation stage. Therefore, we can use our estimates of J KMC to estimate the parameters F c and J 0 in () putting: J 0 exp( βf c ) = J KMC. (7) As discussed in Sec. 2, variations of the activation factor exp( βf c ) with concentration or temperature affect the tota nuceation rate J much stronger than those of the prefactor J 0. Therefore, for the given J KMC, the activation factor exp( βf c ) can be determined from Eq. (7) rather accuratey even if the prefactor J 0 is estimated not too precisey, in particuar, when J 0 is estimated from the STN cacuations described in Sec. 2 which incude mode assumptions (3) and (4). Substitutingfor theprefactor J 0 its expressionviaquantities J 0, v a andγ nn given in the ast ine of tabe 3, and using Eq. (4) which reates γ nn to the effective diffusivity D eff mentioned in Sec. 2, we can re-write reation (7) as the equation for the reduced nuceation barrier βf c : ) (βf c ) KMC = n (2D eff J0 /a 5 0J KMC where index KMC in the eft-hand side indicates that this expression for βf c is based on the KMC simuations. Vaues of the reduced prefactor J 0 can be taken from our STN cacuations iustrated by the ast ine of tabe 3. Hence to find the reduced nuceation barrier βf c from Eq. (8), we need the effective diffusivity D eff. Estimates of D eff based on our KMC simuations of growth of precipitates are described beow in Sec. 5. (8)

9 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 9 4. Estimate of acceeration precipitation factor A irr for eectron irradiation of diute iron-copper aoys The concentrations of point defects (vacancies and interstitias) under permanent irradiation may frequenty exceed their equiibrium vaues by severa orders of magnitudes [22]. In such conditions diffusive phase transformations, such as precipitation or ordering, are strongy acceerated. For the case of eectron irradiation of Fe-Cu aoys this has been observed, in particuar, by Mathon et a. [2] and by Radiguet et a. [6]. If other irradiation effects (such as baistic mixing or radiation induced segregation) can be negected, it is natura to expect that the precipitation kinetics under irradiation is the same as that during therma aging at the same temperature, except for the strong acceeration of this kinetics described by the factor A irr in Eqs. (2) and (3). In this section we estimate this acceeration factor for diute Fe-Cu aoys under typica conditions of eectron irradiation [2] and compare this estimate with both the experimenta and KMC simuation resuts. For diute aoys under consideration, the copper diffusivity under therma aging in Eq. (2) is given by the conventiona expression [25]: D th Cu = α v c eq v D v. (9) Here c eq v is the equiibrium vacancy concentration which is expressed via the enthapy and entropy of vacancy formation, Hv for and Sv for, as: c eq v = exp(sv for βhv for ), and D v is the vacancy diffusion coefficient in a pure iron which is reated to the vacancy migration enthapy and entropy, Hv mig and Sv mig, as foows: D v = a 2 0 ν 0exp(S mig v βh mig v ) (20) where a 0 is the BCC attice constant and ν 0 is the attempt frequency supposed to have the order of the Debye frequency. The coefficient α v in (9) describes vacancy-soute correations, and it is commony written as [25]: ( ) α v = f 2 exp βg bin Cu v where f 2 is the impurity correation factor, and G bin Cu v is the copper-vacancy binding energy. These parameters may be obtained from experimenta measurements, or from ab initio cacuations [0]. Simiary, under irradiation one may write: (2) D irr Cu = α vc irr v D v +α i c irr i D i, (22) where c irr v and c irr i are the vacancy and the interstitia concentrations under irradiation. It seems natura to assume that the copper-vacancy correation effects and binding are not modified by irradiation. Detaied information on the interstitia diffusion (needed for cacuations of α i and D i ) is usuay more difficut to obtain. However, when the point defect concentrations reach their steady-state vaue, one can show that c irr v D v = c irr i D i [22]. If α v and α i are of the same order of magnitude, one finay gets DCu irr 2α vc irr v D v, which corresponds to the foowing acceeration factor in (2) and (3): A irr = 2(c irr v /ceq v ). (23) To compare our KMC simuations during therma aging and the experimenta

10 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 0 Tayor & Francis and I.T. Consutant kinetics of precipitation under irradiation, the key point is therefore to get a reiabe estimation of the vacancy point defect concentrations under irradiation. According to the rate theory modes [22], if point defects created by irradiation disappear by mutua recombination or by annihiation at disocations (which are assumed to be the dominant point defect sinks), the evoution of their concentrations in a pure meta is described by the foowing equations: dc v dt = K Rc ic v ρ d D v (c v c eq v ) dc i dt = K Rc ic v ρ d D i (c i c eq i ) K Rc ic v ρ d D i c i. (24) Here K is the rate of formation of Frenke pairs under eectron irradiation (in dpa/s),ρ d isthedisocation density, andristhevacancy-interstitia recombination rate: R = 4πr vi (D v +D i )/v a (25) where r vi is the recombination radius, and v a = a 3 0 /2 is the atomic voume. For pure iron (for which attice constant is a 0 =0.288 nm), most of the point defect properties have been estimated by ab initio cacuations. In particuar, for the vacancy formation parameters we have: Hv for = 2.8 ev [0] and Sv for = 4.08 [26], which gives: c eq v at 290 C. The migration enthapies are: Hv mig = 0.68 ev [0] and H mig i = 0.34 ev [27]. The Debye frequency is ν 0 = 0 3 s, and with these parameters a vacancy migration entropy Sv mig 2.2 is required to get the experimenta pre-exponentia factor for sef-diffusion of iron (see Refs. [0, 28] for detais). Lacking both experimenta and ab initio estimations, we wi assume the same vaue for the migration entropy of the interstitias. The irradiation conditions of ref. [2] are: K = dpa/s at 290 C, with disocation densities between ρ d = 0 8 cm 2 and ρ d = 0 cm 2. The evoution of point defect concentrations, obtained by numerica integration of Eq. (24) with these two disocations densities, is given on Fig. 4. In both cases, the steady-state vauesarereachedveryrapidy(after essthan0 s),withvacancyconcentrations c irr v and Theacceeration factor is thereforea irr for the highest disocation density and A irr for the owest. In Ref. [2], a precipitate density of m 3 is observed under eectron irradiation of Fe-0.3%Cu aoy after approximatey s. Our KMC simuations of precipitation during therma aging shown in Fig. 2(f) predict a simiar density after s. This corresponds to an acceeration by a factor A irr under irradiation. Considering the strong approximations used in the rate theory mode, and the uncertainties on the point defect properties, the agreement between our estimate of A irr and the combination of experimenta and KMC resuts seems to be reasonabe. A more detaied comparison of our simuations with experiments by Mathon et a. [2] is hindered by a ow resoution in measurements of precipitate sizes R by SANS method used in [2]: R 0.5 nm. It can expain quantitative disagreements between the resuts reported in [2] and those shown in frames 2(e) and 2(f). The maximum vaue of the precipitate density d p (t) reported in [2], d max /m 3, is by about.5 times ower than thefina vaues d p (t) shown in Figs. 2(e) and2(f), even thoughthemaximumvaued max is probabynot reached yet in these simuations. Hence, our simuated d p exceed those reported in [2] by at east severa times. At the same time, vaues of the mean precipitate radius for the nuceation

11 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine stage (which corresponds to d p (t) < d max ) reported in [2], Rm,M n..65 nm, arebyabout2-3timeshigherthantheanaogousvauesrm n nmobserved in both experiments [29, 30] and simuations [0, ] for other Fe-Cu aoys, whie our R m in frames 2(e) and 2(f) are cose to these usua Rm. n Therefore, it seems probabe that because of the above-mentioned ow resoution R 0.5 nm, vaues of the precipitate density d p reported in [2] are underestimated, whie those of the mean precipitate size R m are overestimated. More accurate measurements of d p and R m in diute Fe-Cu aoys under eectron irradiation are evidenty needed for a quantitative comparison with our simuations. 5. Estimates of effective diffusivity of copper atoms for growth of a precipitate using kinetic Monte Caro simuations In the phenomenoogica theory of growth of a supercritica embryo deveoped for a spherica embryo of a radius R with a sharp edge, the tempora dependence R(t) (for the rigid attice aoy mode used) is described by the foowing equation [4]: dr/dt = D eff (R R c )(c c b )/R 2. (26) Here D eff is theeffective diffusivity of a soute atom (in our case, of a copperatom) in an aoy; R c is the critica radius; c = x is the copper concentration far from the embryo; and c b is the binoda concentration (soubiity imit) for temperature under consideration. In more reaistic descriptions, in particuar, in the STN or KMC simuations described in sections 2 and 3, both critica and supercritica embryos have not a sharp edge but a diffuse surface iustrated by Fig.. More important, in the phenomenoogica derivation [4] of Eq. (26), the diffusivity D eff was supposed to be constant independent of the precipitate size N and its mobiity m(n). At the same time, the recent KMC studies [0, 2] have shown that the sma precipitates under consideration which contain N 00 copper atoms are highy mobie, being much more mobie than individua copper atoms. It eads to a great acceeration of precipitation kinetics which quaitativey corresponds to an increase of the effective diffusivity D eff in Eq. (26) with respect to the diffusivity of an individua copper atom. In spite of a these simpifications, the phenomenoogica equation (26) is commony beieved to reaisticay describe growth of supercritica embryos [4]. Accepting this point of view, we first present the expicit soution of Eq. (26) for R(t). Then we use our KMC simuations to access vaidity of Eq. (26) and to estimate the parameter D eff in this equation for diute Fe-Cu aoys under consideration. First we note that empoying the simpest sharp-edge mode of an embryo in usua derivations of Eq. (26) seems to be unessentia, and this equation (at arge N c considered) appears to be vaid for any reaistic description of growth of an embryo, incuding that used in KMC simuations. The radius R(t) in Eq. (26) mainy characterizes the tota number N of soute atoms within the embryo, and for any precipitate this radius can be expressed via N anaogousy to Eq. (8) for the critica embryo: R = a 0 (3N/8π) /3 = 0.42N /3 nm. (27) The right-hand-side of Eq. (26) contains two basic factors (in addition to the diffusivity D eff and a geometrica factor /R 2 ) which naturay describe the driving force forgrowth of theembryo: thefactor (c c b ) that characterizes supersaturation

12 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 2 Tayor & Francis and I.T. Consutant of a metastabe aoy, and the factor (R R c ) that describes vanishing of this force at R = R c. Therefore, the phenomenoogica equation (26) with the generaization (27) seems to be rather pausibe, and comparison to our KMC simuations given beow seems to confirm its vaidity, at east for the R and t intervas studied. First we present the expicit soution R(t) of Eq. (26). As discussed in detai in [4], the embryo can be considered as supercritica ony when its size R exceeds some vaue (R c + ) where determines the scae of critica fuctuations of sizes near R c. Therefore, the time t f when the supercritica embryo has been eventuay formed is defined by the reation: R(t f ) = (R c + ). (28) For the sharp-edge mode of the embryo used in [4], the fuctuation width is expressed via the interfacia energy σ as: = (T/8πσ) /2. To generaize this estimate to the case of rea embryos with diffuse interfaces, we can express σ via the critica radius R c and the nuceation barrier F c. It yieds: = R c α, α = (2T/3F c ) /2. (29) For arge embryos under consideration with βf c, vaues of α are sma; in particuar, for the Fe-Cu aoys isted in tabes and 2, we have: α Integrating Eq. (26) from the initia time t f defined by Eq. (28) to the arbitrary time t when R(t) = R, we obtain: (R R c )(R+3R c + )/2+R 2 c n[(r R c )/ ] = D eff (c c b )τ (30) where τ = (t t f ) is the tota time of growth of a supercritica embryo. If we describe this growth by a reduced variabe y = (R R c )/R c, Eq. (30) takes a universa form which contains ony the reduced dimensioness time ξ: y(2+α+y/2)+n(+y/α) = ξ, ξ = D eff (c c b )τ/r 2 c. (3) At arge R R c, the dependence R(τ) given by Eq. (30) or (3) takes the form R(τ) = [2D eff (c c b )τ] /2 (32) which describes the diffusion-controed growth of a spherica embryo, see, e. g., [3]. Now we compare the description of growth of an embryo by the phenomenoogica equation (30) or (3) with that given by our KMC simuations. First, we note that for our mode of irradiated aoys based on Eqs. (2) and (3), the reduced time ξ in Eq. (3) does not incude the irradiation acceeration factor A irr, just as the righthand side of Eq. (8), and this reduced time is expicity expressed via the therma aging time τ = τ th = (t th t f th ) used in our KMC simuations. Therefore, growth of a supercritica embryo can be described, from one side, by the function R(ξ) determined by Eq. (30) or (3). From the other side, such growth can be foowed in our KMC simuations for any precipitate chosen, which yieds the dependence R = R KMC (t th ) with R defined by Eq. (27). Equating these two quantities, R(ξ) = R KMC (t th ), (33)

13 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 3 with ξ given by Eq. (3) for τ = τ th, we obtain the equaity which contains ony one unknown parameter D eff at a times τ th > 0. Onecanexpect thattheproperchoice ofthissingeparameter canprovideagood accuracy forobeyingequaity (33)at aτ th consideredonyifthephenomenoogica equation (26) hods true for these τ th. Therefore, the check of vaidity of reation (33) enabes us, first, to access reiabiity of the phenomenoogica equations (26) and (30) for description of growth of supercritica embryos and, second, to estimate the effective diffusivity D eff in these equations for diute Fe-Cu aoys. It is convenient to characterize the effective diffusivity D eff (x,t) by its ratio A c to the appropriate diute aoy diffusivity D(0,T) = D 0 (T) writing D eff as D eff (x,t) = A c D 0 (T). (34) The parameter A c = A c (x,t) in (34) characterizes the effective acceeration of growth of precipitates being mainy due to the above-mentioned high mobiity of sma copper custers. To differ it from the irradiation acceeration factor A irr discussed in Sec. 4, A c wi be caed the acceeration diffusion parameter. For the diute aoy diffusivity D 0 (T) of a copper atom in a pure iron at T = C we wi use the conventiona Arrhenius-type expression D 0 (T) = Aexp( βq) (35) with the A and Q vaues suggested by Soisson and Fu [0] basing on combinations of their ab initio cacuations and empirica estimates: A = m 2 /s, Q = 2.67 ev. (36) In Figs. 5 and 6 we show the phenomenoogica functions R(ξ) defined by Eqs. (30) and (3) together with the dependencies R KMC (t th ) observed in our KMC simuations.eachcurver KMC (t th )inthesefiguresstartsfromthetimeofformation of the precipitate chosen observed in the KMC simuation, whie the curve R(ξ) = R[ξ(t th )] starts from the time t f th of formation of the supercritica embryo defined by Eq. (28). Each frame in these figures corresponds to some of our simuations; for other simuations (not shown in Figs. 5 and 6), the resuts are simiar. The parameter A c and the time t f th used to draw the curve R(ξ) in each frame have been estimated from the best fit of this R(ξ) to the R KMC (t th ) curve at this frame. In Fig. 7 we show the concentration and temperature dependencies of the acceeration diffusion parameter A c obtained in these KMC estimates. Let us discuss the resuts presented in Figs First, we note that some of irreguarities in dependencies R KMC (t th ) seen in Figs. 5 and 6 can be reated not to the rea statistica fuctuations of sizes but to the methodica errors due to too ong tempora intervas between sequent savings of KMC data. This seems to be the case, in particuar, for initia stages of simuations shown in frames 5(a), 6(c) and 6(d), and aso for some other frames, e. g., 6(a). At the same time, for the ater stages of growth whenthe precipitate becomes supercritica, R (R c + ), thesemethodica distortions seem to be ess significant. Second, Figs. 5 and 6 ceary iustrate the strong fuctuations of sizes in the critica region R (R c + ) mentioned in the discussion of Eqs. (28) and (29). Such fuctuations are pronounced, in particuar, in frames 5(b), 5(c), 6(c), 6(e), and 6(f). Third and most important, Figs. 5 and 6 show that the description of growth of precipitates by the phenomenoogica equations (26) and (30) seems to agree we with the KMC simuations, at east up to R.5R c, that is, under increase of the number of copper atoms in the growing

14 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 4 Tayor & Francis and I.T. Consutant embryo by about four times. Figs. 5 and 6 aso ceary iustrate the effects of custer mobiity and their direct coaguation mentioned above. In particuar, sudden increases ( jumps ) of R(t) seen in Fig. 5(d) at t th = s, as we as in Fig. 5(f) at t th = and t th = s, occur simutaneousy with the disappearance (evidenty, due to coaguation) of one of supercritica custers in Figs. 2(b) and 2(f). Simiary, jumps of R(t) in Fig. 6(a) at t th = s and in Fig. 6(b) at t th = s occur simutaneousy with the disappearance of one supercritica custer in Figs. 3(b) and 3(d), respectivey. The smaer jumps in R(t) seen in Figs. 5(d), 5(e), 6(a) and 6(b) can correspond to coaguation of subcritica custers which are not registered in Figs. 2 and 3. At the same time, for simuations with N p = shown in Figs. 5(a)-5(c) and 6(d)-6(f), for which ony one supercritica custer is present in the simuation box, the arge and distinct jumps in R(t) are not ceary seen. Inspection of atomic distributions for these simuations shows that growth of this supercritica custer is mainy reaized via its fast diffusion among amost immobie individua copper atoms which are sometimes swept and absorbed by this mobie custer. Some subcritica custers containing severa copper atoms are aso observed in these simuations, and they seem to diffuse notaby sower than the supercritica custer, in a quaitative agreement with Fig. 9 in [0]. Theresutingacceeration diffusionparameter A c for diutefe-cu aoys is shown in Fig. 7. It is rather high: A c 200, having the same scae as the parameters of acceeration of precipitation in Fe-.34%Cu aoys at T = 500 C due to the high mobiity of copper custers studied in Refs. [0] and [2]. Fig. 7 aso shows that the concentration and temperature dependencies of this acceeration are rather sharp. Fig. 7(a) iustrates thenotabe increaseof A c withincrease ofcopperconcentration x; it can be expained by a probabe enhancement of density of mobie copper custers (both supercritica and subcritica) under increase of the copper content in an aoy. Fig. 7(b) showsasignificant decreaseof A c undereevating temperaturet. This aso seems natura, as eevating T shoud ead to the weakening of the coppervacancy binding proportiona to the Mayer functions [exp(βevcu bn from (4) [8]. Hence the strong trapping of vacancies by copper custers (which is the physica origin of their high mobiity [0]) shoud weaken. ) ] with Ebn vcu 6. Comparison of nuceation barriers found using the STN cacuations and the KMC simuations Estimates of the nuceation rate J KMC described in Sec. 3 and iustrated by Figs. 2 and 3, combined with the estimates of the effective diffusivity D eff discussed in Sec. 5 and iustrated by Figs. 5-7, enabe us to find the reduced nuceation barriers βf c in Eq. (8) basing on the kinetic Monte Caro simuations. This equation (8) aso shows that for the arge βf c considered which vary with x and T very sharpy, possibe errors of a reative order of unity in our estimates of parameters J 0, J KMC and D eff make no significant effect on the βf c vaues obtained. In tabe 4 we present the (βf c ) KMC vaues estimated for a simuations shown in Figs. 2 and 3, together with the anaogous (βf c ) STN vaues cacuated in Sec. 2 with the use of the statistica theory of nuceation [2, 3]. In comparison of these KMC and STN resuts we shoud remember that the STN assumes the size and the nuceation barrier of the embryo to be arge: N c, βf c. Therefore, the accuracy of this theory and its agreement with the KMC simuations shoud improve when the size of the embryo increases. The resuts presented in tabes, 2 and 4 agree with these considerations. We see that at ow concentrations x 0.2% when both N c and βf c are sufficienty arge: N c 2, βf c 20, the agreement

15 Juy, 207 3:29 Phiosophica Magazine VSZ-fina Phiosophica Magazine 5 between KMC and STN resuts is virtuay perfect: for a nine simuations with x 0.2% shown in tabe 4, differences between the (βf c ) KMC and the (βf c ) STN vaues have the order of a percent. At the same time, at higher x equa to 0.3 or %, the N c and βf c vaues notaby decrease, and differences between (βf c ) KMC and (βf c ) STN increase, which may refect the owering of accuracy of the STN. Therefore, the resuts presented in tabe 4 enabe us to make an important concusion that at high vaues βf c 20, the STN-based cacuations of nuceation barriers are highy reiabe. At the same time, as mentioned in Sec. 2, the homogeneous precipitation imit in diute Fe-Cu aoys corresponds just to the vaues βf c > 20. Therefore, the nuceation barriers near the homogeneous precipitation imits in these aoys can be reiaby cacuated using the STN. 7. Temperature dependencies of precipitation imits in diute iron-copper aoys Cacuations of nuceation barriers described in Secs. 2 and 6 can be used for tentative estimates of the homogeneous precipitation imit x hom (T) at different temperatures T. As mentioned in Sec. 2, it seems natura to suggest that the position of this imit is mainy determined by the vaue of the activation factor exp( βf c ), that is, by the reduced nuceation barrier βf c which at x = x hom takes a certain high vaue (βf c ). It seems aso natura to assume that these precipitation imits at different temperatures T correspond to the simiar vaues of the activation factor exp( βf c ). It impies the foowing reation: (βf c ) βf c [ x hom ] (T),T C hom (37) where C hom is a constant independent of temperature, and thus it can beestimated using the experimenta vaue x hom (290 C) x exp. For exampe, using estimates of x exp mentioned in Sec. 2, x exp % from [2] or x exp 0.5% from [3], we obtain: C hom 28 or C hom 23. Then equation (37) with the function F c (x,t) cacuated using the STN enabe us to find the x hom (T) vaue. A high accuracy of the STN cacuations of nuceation barriers combined with the above-mentioned physica considerations aow us to expect that such semi-empirica estimates of the homogeneous precipitation imit for temperatures T of practica interest can be sufficienty reiabe. One can aso try to extend this approach to the case of the heterogeneous precipitation of copper on custers of irradiation-induced point defects, vacancies and sef-interstitia atoms (that is, to the case of neutron or ion irradiation []-[9]), to estimate the appropriate heterogeneous precipitation imit x het, which is more interesting for appications. Even though kinetic paths of the heterogeneous and the homogeneous precipitation differ from each other(see, e. g., [32]), the main physica origin for suppressing precipitation at ow soute concentrations x seems to be the same for both processes. In both cases, formation of a precipitate needs overcoming the energetic barrier due to the surface oss in the free energy F which is not compensated by the voume gain in F unti the precipitate voume becomes sufficienty arge. Differences between these two processes are usuay reated mainy to some geometrica factors which depend on the structure of the heterogeneity [32]. Therefore, one may expect that the heterogeneous precipitation imit x het can aso be estimated from a phenomenoogica reation simiar to (37), but the imiting vaue F c in this reation shoud be higher as the nuceation barrier for the heterogeneous nuceation, Fc het, is ower than that for the homogeneous nuceation,

16 Juy, 207 3:29 Phiosophica Magazine VSZ-fina 6 Tayor & Francis and I.T. Consutant Fc hom.however, ifweassumethatthedifferencebetweenfc het andfc hom has mainy a geometrica origin (which is the case for the simpest modes of heterogeneous nuceation [32]), then the ratio Fc het /Fc hom can be supposed to weaky vary with temperature. Then the temperature dependence of the heterogeneous precipitation imit x het (T) can be estimated from a semi-empirica reation simiar to (37): [ ] βf c x het (T),T C het (38) where F c (x,t) is again the STN cacuated nuceation barrier for the homogeneous nuceation, whie the constant C het is estimated using the experimenta vaue x het (290 C) x het,exp. For numerica estimates we use the vaue xhet,exp 0.06% given in Refs. [5] and [7], which yieds: C het 35. In figure 8 we show positions of the homogeneous and the heterogenous precipitation imits in the (x,t) pane, T hom (x) and T het (x), cacuated for diute Fe-Cu aoys using Eqs. (37) and (38). For comparison, in Fig. 8 we aso show the binoda curve T b (x) (the copper soubiity imit) estimated in [0] from experimenta data and ab initio cacuations. The resuts in Fig. 8 are presented for the temperature interva between 290 and 450 C used for anneaing of RPVs [3, 7]. The vaue C hom = 27.6 used in Fig. 8 corresponds to x hom (290 C) = 0.088% for which Radiguet et a. [6] did not observe any homogeneous precipitation, whie the constant C het = 35 used in Fig. 8 corresponds to the estimate x het = 0.06% given by Mier et a. [5, 7]. Our cacuations aso show that a sight variation of the constant C hom or C het eads to an amost rigid shift of the curve T hom (x) or T het (x) in Fig. 7. For exampe, when this variation of C hom eads to the shift δx hom =0.05% to the right at 290 C, the anaogous shift at 450 C is amost the same: δx hom (450 C) 0.06%. ThemostinterestinggenerafeatureseeninFig.8seemstobearathersharptemperature dependence of the precipitation imits, particuary for the heterogeneous precipitation. For exampe, the x het vaue at T = 450 C exceeds that at T = 290 C by amost five times, whie the x hom vaue increases for this temperature interva by about 3.5 times. These quaitative concusions can be usefu, in particuar, for the interpretation of microstructura observations reated to the anneaing of RPVs [7]. Note aso that the variations of both precipitation imits, x het (T) and x hom (T), with temperature for the temperature interva considered are much stronger than those for the soubiity imit x b (T). To concude this section, we again note that our equations (37) and (38) for temperature dependencies x het (T) and x hom (T) are of course tentative and have not been proved formay. At the same time, physica considerations about the dominant roe of the activation factor exp( βf c ) in suppressing precipitation at ow x used in our derivation seem to be rather pausibe, particuary for the homogeneous precipitation. In the derivation of Eq. (38) for the heterogeneous precipitation, we aso used an assumption about a weak temperature dependence of the ratio Fc het /Fc hom for the temperature interva considered which may seem to be ess evident. Therefore, the accuracy of Eq. (38) for x het (T) can, generay, be ower than the accuracy of Eq. (37) for x hom (T). However, one may expect that the main quaitative features of both dependencies, T hom (x) and T het (x), shown in Fig. 8 are correcty described by the simpe mode used.