arxiv: v1 [cond-mat.mtrl-sci] 31 Oct 2012

Transcription

1 A Mult-Component Phase Feld Crystal Model for Structural Transformatons n Metal Alloys arxv: v1 [cond-mat.mtrl-sc] 31 Oct 01 Nana Ofor-Opoku, 1, Vahd Fallah, 1, Mchael Greenwood, 3, 1 Shahrzad Esmael, and Nkolas Provatas 4, 1 1 Department of Materals Scence and Engneerng, McMaster Unversty, 180 Man Street West, Hamlton, Canada L8S-4L7 Department of Mechancal and Mechatroncs Engneerng, Unversty of Waterloo, 00 Unversty Avenue West, Waterloo, Canada NL-3G1 3 CanmetMaterals, NRCan, 183 Longwood Road South, Hamlton, Canada L8P-0A5 4 Department of Physcs and Centre for the Physcs of Materals, McGll Unversty, 3600 Rue Unversty, Montreal, Canada H3A-T8 We present a new phase feld crystal model for structural transformatons n mult-component alloys. The formalsm bulds upon the two-pont correlaton kernel developed n Greenwood et al. for descrbng structural transformatons n pure materals [Phys. Rev. Lett. 105, (010)]. We ntroduce an effectve twopont correlaton functon for mult-component alloys that uses the local speces concentratons to nterpolate between dfferent crystal structures. A smplfed verson of the model s derved for the partcular case of threecomponent (ternary) alloys, and ts equlbrum propertes are demonstrated. Dynamcal equatons of moton for the densty and multple speces concentraton felds are derved, and the robustness of the model s llustrated wth examples of complex mcrostructure evoluton n dendrtc soldfcaton and sold-state precptaton. PACS numbers: K-,61.50.Ah,81.10.Aj, x I. INTRODUCTION Engneerng alloys requre the addtons of multple components to acheve desred propertes. Ths, however, makes the nvestgaton of ther mcrostructure evoluton and defect nteractons dffcult. The propertes, and therefore the resultant behavour, of alloys can drectly be correlated to the chemcal make up, mcrostructure and the phase selecton processes these alloys undergo upon soldfcaton and subsequent downstream processng, such as thermal treatments. In the case of bnary alloys, models of soldfcaton processes such as nucleaton, free growth and coarsenng knetcs, segregaton and second phase formaton have been relatvely well developed. However, for mult-component alloys, the complex nteractons nvolved between the dfferent chemcal speces, dslocatons and other defects make such phenomena far more dffcult to study, even wth advances n characterzaton technques such as conventonal and hgh resoluton transmsson electron mcroscopy. Advances n modellng have sgnfcantly mproved our understandng of the fundamental nature of mcrostructure and phase selecton processes. Notable contrbutons have been made usng the phase feld methodology (PFM), whch has been successful at examnng mesoscale mcrostructure evoluton over dffusve tme scales. The greatest success of the PFM has come n the area of soldfcaton [1 7]. The phase feld concept has gone far beyond ts orgns. It s now capable of descrbng, through the ntroducton of varous auxlary felds, a wealth of phenomena such as multple crystal orentatons [8 10], multple components and phases [11 13], defect-solute nteractons [14], elastcty [15, 16] and plastcty [17]. There has recently emerged an atomc-scale modellng formalsm, called the phase feld crystal model (PFC) [18, 19]. Electronc address: oforon@mcmaster.ca Ths method, operates on atomstc length scales and dffusve tme scales and self-consstently ncorporates elastcty, multple crystal orentatons, gran boundares, dslocatons, and the evoluton of mcrostructure on dffusve tme scales. For both pure materals and bnary alloys, Elder and coworkers [19] and Jn and Khachaturyan [0] have shown that PFC models can be formally derved from classcal densty functonal theory (CDFT), where the order parameter can be related to the atomc probablty densty [1]. As such, many basc mcrostructure phenomena can be seen as arsng selfconsstently from a smple fundamental theory descrbed by a small set of physcally motvated parameters. Wth the ablty of the PFC densty feld to also assume dsordered states, t s also possble to examne amorphous or glassy states [, 3]. Phase feld crystal models are also exceedngly smple to work wth numercally. The use of coarse granng approaches has further shown that PFC-type models can be used as generators of tradtonal phase feld models, as well as so-called ampltude models, essentally phase feld models wth complex order parameters. These make t possble to smulate dfferent crystal orentatons and defect structures on mesoscopc length and tme scales [4 7], and also explot the scalng afforded by adaptve mesh refnement [8]. A weakness of the early PFC models was ther nablty to systematcally descrbe and control complex crystal structures and coexstence between them. Greenwood et al. [9, 30] addressed ths shortcomng by ntroducng a class of multpeaked, two-pont drect correlaton functons that contaned some of the salent features of CDFT, but retaned the smplfcatons that gave the orgnal PFC formalsm ts numercal effcency. Ths so-called XPFC formalsm was later extended to bnary alloys, and appled to phenomena such as eutectc soldfcaton and elastc ansotropy [31], solute drag [3], quas-crystal formaton [33], solute clusterng and precptaton mechansms n smplfed Al-Cu alloys [34] and 3D stackng fault structures n FCC crystals [35]. In ths paper, we generalze the XPFC formalsm of Greenwood et. al to the case of N-component alloys. The approach

2 begns wth the truncated CDFT energy functonal of an N- component system. At the core of our excess free energy are the partcle nteractons of Ref. [9, 30], adapted for dfferent structural phases n alloys by makng the nteracton kernel a functon of the local speces concentratons. We compute the equlbrum propertes of our model for the case of a ternary alloy and compare the resultng model phase dagram to an expermental ternary system. The dynamcs of the model are then demonstrated n the context of dendrtc soldfcaton and sold-state precptaton The remander of ths paper s organzed as follows. We begn wth dervng the full N-component XPFC energy functonal n Secton II from a smplfed, truncated classcal densty functon theory of freezng smlar to that of Ramakrshan and Yussouff [36]. We then derve a second, smplfed verson of the model that s the N-component analogue of the prevous bnary XPFC model n lterature. Secton III calculates the equlbrum propertes of the model for a partcular case of a ternary system va sothermal sectons of the phase dagram. Secton IV presents some numercal examples of mcrostructure evoluton by smulatng dendrtc soldfcaton and sold-state precptaton. We end wth a summary and conclusons. II. XPFC ENERGY FUNCTIONALS FOR N-COMPONENT ALLOYS A general free energy functonal for an N-component alloy s derved startng from the classcal densty functonal theory of freezng energy formalsm of Ramakrshan and Yussouff [36], where each alloy component s wrtten n terms of a densty feld ρ. The model s re-wrtten n terms of total densty and concentraton varables to make contact wth standard models used n the descrpton of alloys. The model s then collapsed to a smplfed form of the free energy, smlar to the smplfed form for the bnary XPFC model of Greenwood et al. [31]. Fnally, equatons of moton for the total densty and each concentraton feld are presented for both versons of the model free energy. A. Dervng an XPFC Energy Functonal for N-component Systems The free energy functonal of an N-component mxture can be descrbed by two contrbutons; a local free energy for each of the N densty felds and an excess free energy due to speces nteractons. The local free energy s treated as an deal energy whch drves the densty felds to become unform. The excess contrbuton drves the densty felds to become perodc by creatng mnma n the free energy for these states. We can wrte the free energy functonal of the mxture as F k B T = dr { Fd k B T + F } ex, (1) k B T where F d denotes the deal energy and F ex s the excess energy whch accounts for nteractons between atoms through correlatve nteractons. Ths latter term, gves rse to structural symmetry, elastcty and nteractons between topologcal defects. The constant k B s the Boltzmann constant and T the temperature. The dfferental dr dxdydz. The deal energy, F d, gves the entropc contrbuton for an N-component system. For small densty changes from a reference densty of each component, t s defned as F N d k B T = ρ ln ( ρ ρ o ) δρ, () where N denotes the number of components, whch are denoted as A, B, C,...,etc., ρ s the densty of component, and ρ o s the reference densty of component n the lqud phase at co-exstence. Followng prevous PFC models [19], we defne a total mass densty ρ = N ρ and the total reference mass densty as ρ o = N ρ o. Followng Refs. [19, 7, 31], we defne concentratons as c = ρ /ρ and the correspondng reference compostons by c o = ρo /ρo. Furthermore, for convenence we defne a dmensonless mass densty of the form n = ρ/ρ o 1. Wth these defntons and the conservaton condton c 1, Eq. () smplfes to the dmensonless form N F d k B T ρ o =(n+1) ln (n + 1) n+(n+1) c ln c c o. (3) The excess energy takes nto account nter-partcle nteractons truncated at two-partcle,.e., A-A,B-B,..., N-N, A- B,...,A-N, nteractons. Ths can be defned as, F ex k B T = 1 dr N j δρ (r) j (r, r ) δρ j (r ), (4) where C j represent all combnatons of two-partcle correlatons, n ths work assumed sotropc (.e., C j (r, r ) = C j ( r r )), between the feld descrbng speces and j, respectvely, where, j = A, B, C,..., N. We wrte Eq. (4) n terms of the reduced densty n and compostons c. As n Refs. [19, 31], we consder only the lowest order contrbutons of the compostons c, whch vary on length scales much larger than the densty n, whch are perodc on the scale of the lattce constant. Ths allows us to smplfy ntegrals arsng from Eq. (4), whch couple c (r ) together wth n(r ) [47]. For example, ( r r )n(r )c (r ) c (r) ( r r )n(r ). To smplfy notaton, the notaton n(r ) n and c (r ) c s used hereafter. Wth these smplfcatons and notatons, the excess energy of Eq.(4) can be wrtten n terms of the dmen-

3 3 sonless varables n and {c } as F ex k B T ρ o = s the Fourer transform of Cj ( r r ), and sats- where Ĉj fes dr [n c c j + c c j c o c j ] dr [n c + c c o ] n c j (5) dr [ c o c o j n c o j c c o j c ] Ĉj ( k =0), Ĉ j ( k =0) = ( r r ), (6) and where we have ntroduced the notaton C j ρ o C j ( r r ), whch s the drect two-pont correlaton functon. Collectng terms from Eqs. (3) and (5) gves the complete N-component free energy functonal, wrtten n dmensonless form, F k B T ρ o = dr (n + 1) ln (n + 1) n+ F mx ({c }) (n+1) 1 dr [n c c j + c c j c o c j ] n 1 1 dr [n c + c c o ] c j (7) dr [ c o c o j n c c o j c o j c ] Ĉj ( k =0), where F mx ({c }) denotes the deal entropy of mxng, F mx ({c }) = c ln c c o. (8) Equaton (7) s the full N-component PFC model n CDFT form. When a form for C j s specfed, t can be used drectly. However, ths form s not convenent to make contact wth other theores and models n the lterature. It wll be transformed nto a smpler form n the next secton. B. Smplfed N-Component XPFC Free Energy It s nstructve to reduce the model of Eq. (7) to a mnmal form that retans the salent features of the orgnal model but can also make contact wth prevous PFC and phase feld models. To do so, certan smplfcatons must be made. Frst, an expanson of the deal free energy term s taken to fourth order n the lmt of small n,.e., around the reference ρ o. The logarthms n the entropy of mxng (Eq. (8)) are left unexpanded for convenence. Secondly, the terms wth correlaton kernels can be smplfed by retanng the long wavelength behavour of all compostons c, where they vary much more slowly than n. Followng the procedures outlned n Refs. [6, 7, 37, 38], t can be shown that upon coarse granng, all terms contanng lnear powers of n or n n Eq. (7) vansh. Also, terms contanng only concentraton felds and a correlaton functon gve rse to local products of c c j (whch arse from the k = 0 part of C j, and look analogous to the last term n Eq. (7)) and products between ther correspondng gradents. The reader s referred to Appendx A for detals of the coarse granng procedure appled to terms of Eq. (7). After some tedous but straghtforward algebra, the above approxmatons lead to the followng smplfed N-component XPFC free energy functonal, { n n3 F = dr η 6 +χn4 1 +ω F mx({c })(n + 1) 1 n dr C eff ( r r ) n + 1 κ j c c j }, (9) where C eff ( r r ) = =1 c c j C j ( r r ). (10) The parameters η, χ and ω are constants, the sgnfcance of whch s dscussed further below. The κ j are gradent energy coeffcents assocated wth compostonal nterfaces nvolvng c and c j. For notatonal convenence, F s used to denote F/k B T ρ o. The parameters η and χ correspondng to Eq. (7) are formally equal to one, but hereafter wll be treated as free parameters that can be used to correct the densty dependence of the deal free energy away from the reference densty ρ o,.e., to match the bulk free energy to materals propertes. Also, t was shown n Ref. [38] that the k = 0 mode of hgher-order correlaton terms n a CDFT expanson wll contrbute local polynomal terms n c and n, analogous to the Ĉj ( k = 0) terms of Eq. (7). These terms can be combned wth an expanson of the F mx term n Eq. (7) to produce a messy polynomal expanson of the local free energy n powers of the elements of {c } and n. To keep the form of the free energy compact, we have found that t s smpler to ntroduce a parameter, ω, whch modfes the mxng free energy from ts deal form, away from the reference compostons c o. Fnally, n the present work, the gradent energy coeffcent tensor wll be assumed to be dagonal for smplcty,.e., κ j = 0 for j and κ > 0 for all. The correlaton functon n Eq. (10) s too basc to capture the propertes of very complex alloys although t can capture some propertes of smple alloys. Guded by the form of the frst term on the second lne of Eq. (7), t can be seen that hgher-order correlaton functons wll contrbute terms of the form c c j c k C jk 3, c c j c k c l C jkl 4, etc. To emulate such

4 4 hgher-order non-local contrbutons effectvely, we ntroduce an effectve correlaton functon of the form C eff ( r r ) = =1 X ({c j }) C ( r r ). (11) The X are as yet undetermned polynomal functons of the elements of {c j }. The role of the X s to determne the resultant local crystallne structure by nterpolatng between the kernels Ĉ (defned below), whch defne the base equlbrum crystal structures of each pure component. The nterpolaton s done through approprately constructed polynomal expansons of the elements of {c j }. The order of X depends on the number of components n the system and can be made as hgh as requred to smoothly nterpolate from one correlaton kernel to another. We have found that Eq. (11), through approprate choces of X, combned wth other model parameters, s robust enough to model a wde varety of alloy systems. The model n Eq. (9) captures the usual features of other PFC models, whle allowng for a very easy control of a wde range of crystal structures n dfferent phases. It s motvated from consderatons of classcal densty functonal theory but smplfed enough to make numercally tractable smulatons possble, as wll be shown below. Fnally, we note that the form of the expanson n Eq. (11) s dmensonally motvated from hgher-order terms n CDFT but s flexble enough to model expermentally relevant mult-component alloys quanttatvely usng, for example, thermodynamc databases. C. Dynamcs Equatons of moton for the densty n and each of the concentraton felds c follow conserved dsspatve dynamcs. Namely the dmensonless densty n obeys ( n t = M n δf ) + ζ n (1) δ n ( { }) = M n n η n + χn3 3 +ω F mx({c }) C eff n + ζ n, whle the dynamcs of each composton feld, c, evolve accordng to ( c t = M c δf ) + ζ c (13) δc ( { }) = + ζ c, M c ω(n + 1) δ F mx 1 δc nδc eff n κ c δc where the followng shorthand notatons have been made, C eff n dr C eff ( r r )n(r ) n δc eff n n(r) dr δc eff ( r r )n(r ). (14) δc δc The coeffcents M n and M c denote the moblty of the densty and each concentraton, respectvely, and strctly speakng can be functons of the felds. The nose terms ζ n and ζ c model coarse graned thermal fluctuatons on densty and concentratons c, respectvely. They formally satsfy ζ q (r, t)ζ q (r, t ) = A χ a (r r )δ(t t ), where q denotes the densty or one of the concentraton felds, A M q k B T and χ a (r r ) s the nverse Fourer transform of a Gaussan functon, whch, followng Tegze and co-workers Ref. [39], can be generalzed to have a hgh frequency cut off for frequences above π/a, where a s the lattce constant. The precse form of A, whch sets the scale of the thermal fluctuatons s not properly understood n the context of PFC modellng but s the object of several nvestgatons. In the applcatons llustrated n ths paper, the nose s left out of smulatons. III. TERNARY SYSTEMS In ths secton, we reduce the smplfed free energy functonal of secton (II B) to the case of three-components, or ternary alloys. We frst descrbe the ternary free energy functonal, followed by a dscusson of the effectve correlaton functon chosen for ternary systems. Wth the free energy and effectve correlaton n hand, we demonstrate the equlbrum propertes of our model by calculatng the ternary phase dagrams for a generc A-B-C system and a smplfed Al-Cu- Mg system. F ter = dr A. Free Energy Functonal Specalzng Eq. (9) for 3-components, denoted here as A, B and C, reduces t to { n where 1 n η n3 6 + χn4 1 + ω F ter mx (n + 1) (15) dr C ter eff ( r r ) n + κ A c A + κ B c B Fmx ter = c A ln c A c o + c B ln c B A c o B + (1 c A c B ) ln (1 c A c B ) 1 c o A, (16) co B },

5 5 and C ter eff ( r r ) = X A (c A, c B )C AA ( r r ) (17) + X B (c A, c B ) C BB ( r r ) + X C (c A, c B ) C CC ( r r ). In arrvng at Eq. (15), the condtons c C = 1 c A c B and c o C = 1 co A co B have been used, and the cross gradent concentraton terms n A and B have been neglected. The effectve ternary correlaton kernel, Ceff ter, s defned by X such that X A + X B + X C 1 at all compostons, analogous to the case for the XPFC bnary model [31]. Ther partcular form, s chosen here to model the generc propertes of eutectc systems. However, by careful alteraton of other parameters, other alloy systems can be modelled, e.g. somorphous and pertectc systems [31]. Here the X used are, X A (c A, c B ) = 3c A + c A c B c 3 A c Ac B c A c B X B (c A, c B ) = c A c B + 3c B c Ac B c A c B c 3 B X C (c A, c B ) = 1 3c A + c 3 A 3c B + c 3 B 4c A c B + 4c Ac B + 4c A c B. (18) B. Correlaton Functons C The XPFC formalsm s best suted for numercal smulaton n Fourer space. The pure component correlaton functons C ( r r ) are thus constructed drectly n Fourer space, where they are denoted Ĉ (k). Each component,, contrbutes a correlaton functon that supports the desred equlbrum crystal structure for a pure component. A Fourer space peak of Ĉ (k) [30], for a gven mode, j, s denoted by Ĉ j = e σ σ Mj (k k j ) e α j. (19) The total correlaton functon for component, Ĉ, s defned by the envelope of all peaks Ĉ j. The frst exponental n Eq. (19) sets the temperature scale va a Debye-Waller prefactor that employs an effectve temperature parameter, σ. We also defne an effectve transton temperature, σ Mj, whch subsumes the effect of planar and atomc denstes assocated wth the famly of planes correspondng to mode j [31]. The second exponental sets the poston of the recprocal space peak at k j, whch defnes the nverse of the nterplanar spacng for the j th famly of planes n the equlbrum unt cell structure of component. Each peak s represented by a Gaussan functon, wth α j beng the wdth of the peak, j. The {α j } have been shown n Ref. [30] to set the elastc and surface energes, as well as ther ansotropc propertes. It s noted that the k = 0 mode of all correlaton functons s essentally zero. In prncple, as dscussed above, the k = 0 mode of these correlaton functons can have ther effects mplctly reflected through local coeffcents n the free energy. In the case of a pure materal, a nonzero peak heght at k = 0 n the correlaton functon merely shfts the local free energy at denstes away from the reference densty, however the stablty of equlbrum structures s typcally unchanged [30]. The stuaton s smlar for alloys, where the k = 0 mode wll have a neglgble contrbuton for phases that reman relatvely close to the reference densty. Devatons of phases away from the reference densty wll be manfested n the average densty dmenson of the phase dagram. Here, t s assumed that the average densty n o = 0 to smplfy the demonstraton of the model. Of course, the more complex stuatons where both the concentraton and average densty need to be modelled can be treated by addng sutable k = 0 contrbutons, or by choosng the approprate coeffcents n the bulk free energy. Thus, wthout loss of generalty, we wll assume no addtonal constant to the correlaton functon Ĉj j here. C. Ternary Dynamcs For the case of 3-component alloys, the dynamcal equatons of motons n Eqs.(1)-(13) reduce to { } n t =M n n η n + χn3 ter + ω Fmx Ceff ter n, (0) 3 c A t =M A c B t =M B { { ω(n+1) ω(n+1) δ F ter mx δc A δ F ter mx δc B } 1 nδcter eff n κ A c A, δc A } 1 nδcter eff n κ B c B, δc B where M n, M ca and M cb are dmensonless moblty coeffcents for densty and compostons felds. They are set to 1 here, snce t s the ntent of ths paper to ntroduce the model and ts physcal features. D. Equlbrum Propertes Ternary equlbrum s defned by co-exstence of bulk phases, e.g. sold α -sold β, lqud-sold α -sold β, etc. The governng propertes, e.g. parttonng, of such an equlbrum state can be determned from standard thermodynamc mnmzaton methods. In general, for 3-component alloys, free energy mnmzaton s defned by a common plane tangent to the free energy wells of any two or three coexstng phases. Ths constructon s a geometrcal representaton of the statement that the chemcal potentals and grand potentals of any two phases are equal wth respect to each component. Here, we construct sothermal ternary phase dagrams by examnng all combnatons of phase coexstence (e.g., sold α -lqud, sold α -sold β, etc.). Procedures for calculatng phase dagrams for PFC models are well-documented [18, 19, 9, 31] and the approach used here wll only be summarzed. For sold phases, the densty feld, whch vares on atomc length scales, s approxmated usng a mult-mode approx-

6 6 maton gven by N n (r ) = j=1 A j N j l=1 ( ) π exp k l,j r, (1) a where a s the lattce spacng of the sold phase and N denotes the number of mode famles (famles of planes) n the unt cell of phase, A j s the ampltude assocated wth the j th famly of planes. Each mode contans N j recprocal lattce peaks, enumerated by the ndex l. Strctly speakng, there s a dstnct ampltude, A l,j, for each recprocal lattce peak. However, for the purposes of smplfyng the constructon of phase dagrams (.e., workng wth the fewest number of varables to mnmze), they are assumed constant leadng to A j. Each ndex l n the famly j has a correspondng recprocal lattce vector k l,j, normalzed by the lattce spacng. Substtutng Eq. (1) nto Eq. (15), and ntegratng over one unt cell, the free energy can be calculated for each phase as a functon of c A, c B and the ampltudes A j. Snce ampltudes are non-conserved felds, the resultng free energy s then mnmzed wth respect to each A j [31]. The result s substtuted back nto the free energy. After ths procedure, we are left wth a free energy landscape Fsol ter (c A, c B ), where Fsol ter represents an ampltude-mnmzed sold free energy. In keepng wth the dscusson of the prevous sectons, we assume that the average densty of all phases s close to the reference densty,.e., n o = 0. For the lqud phase, the free energy Flq ter (c A, c B ) s trvally computed by settng all A j = 0 [48]. Wth the free energy landscapes of lqud and solds, the phase boundary lnes between a combnaton of phases at a gven temperature parameter, σ, are computed by solvng the followng set of equatons smultaneously, µ I c A = µ J c A where the last of these mples, µ I c B = µ J c B Ω I = Ω J, () f I µ I c A c I A µ I c B c I B = f J µ J c A c J A µ J c B c J B. (3) The superscrpts I and J denote any two phases n equlbrum (e.g. lqud-sold α ), respectvely. The expressons µ I c A = f I / c A and µ I c B = f I / c B are the chemcal potentals of phase I wth respect to the concentratons c A and c B, respectvely, wth analogous defntons for µ J c A and µ J c B. The expressons Ω I and Ω J are the grand potentals of phases I and J, respectvely. See Appendx B for further detals on calculatng phase dagrams. The set of condtons n Eq. (), along wth Eq. (B3) defnng the average concentraton, can be solved to fnd the four equlbrum concentratons (two per phase) defnng coexstence on a gven te lne. 1. Generc Ternary Eutectc Alloy A frst example of the equlbrum propertes of the ternary XPFC model are demonstrated for a system where all three components (A, B, and C) are structurally smlar, dfferng only n ther equlbrum lattce spacngs. Here twodmensonal (D) square symmetry s assumed as the equlbrum structure for each pure component, whch n ths context mples that all Ĉ have the same number of peaks, wth the correspondng ratos of ther postons n recprocal space beng the same. However, each structure s dfferentated by the absolute postons (k j ) of each peak. Though t has not been done n ths ntal work, by adjustng the wdths (α j ) of each peak, each element can also be dfferentated by dfferent elastc and surface energes. The full lst of parameters used to construct the phase dagrams n ths subsecton are lsted n the capton of Fg. 1. Allowng all three components to have square structural symmetry, at suffcently low temperature we can construct a bulk sold free energy landscape descrbng multple sold phases, descrbed by an effectve lattce parameter (a ter ) that s a weghted average of the ndvdual lattce parameters of all three components, usng the nterpolaton functons of Eq.(18), namely, a ter = a A X A + a B X B + a C X C. Ths leads to the sold-lqud free energy landscape n Fg. 1(a) for σ = 0.17, where the values of all other parameters are specfed n the fgure capton. The correspondng sothermal phase dagram s llustrated n Fg. 1(c), whch s constructed form the coexstence lnes calculated between the lqud phase and the dfferent sold-soluton phases, usng the set of condtons n Eq. (). Fgure 1(b) shows an sothermal cut at a hgher temperature,.e., σ = 0.18, depctng an ncreased regon where the bulk lqud s stable compared to the sold phases. At suffcently low temperature, the free energy admts eutectc coexstence of three phases. We construct an sothermal cut rght above the eutectc temperature,.e., at σ = 0.164, shown n Fg. 1(d). The correspondng concentratons c A and c B n Fg. 1, are gven as fractons, where unty represents pure A or B, respectvely, along each axs of the phase dagram.. Smplfed Al-Cu-Mg Type Alloy The parameters of the ternary XPFC model can be chosen to produce sectons of expermental phase dagrams qualtatvely, as n the work of Fallah et al. [40], where the present ternary model s used to model precptaton n a D representaton of the Al-Cu-Mg system. Here, we demonstrate how the equlbrum propertes of a porton of the Al-rch (smplfed) part of the Al-Cu-Mg phase dagram can be descrbed quanttatvely by the ternary XPFC model. An expermental phase dagram at 400 C s shown n Fg. (b), taken from Ref. [41]. Consder the part of the phase dagram for (Al) β-θ outlned by the red dashed lne and crcled sold phases n the expermental phase dagram shown n Fg. (b), and gnorng the (Al)+S and (Al)+β+T phase regons. In the dlute-mg regon, a eutectc transton occurs between the Al-rch, (Al)-

7 7 FIG. 1: (Colour onlne) Ternary eutectc system: (a) Sold and lqud energy landscapes of a square-square-square (A-B-C) system at temperature parameter σ = Correspondng phase dagrams at temperatures (b) σ = 0.18, (c) σ = 0.17 and (d) σ = The parameters for deal free energy and entropy of mxng were η = 1.4, χ = 1, ω = 0.005, whle reference concentratons were coa = and cob = Wdths of the correlatons peaks are taken α11 = 0.8 and α10 = α11 for all phases (requred for sotropc elastc constants n a sold phase wth square symmetry [31]). The peak postons for the gven structures are k11a = (81/38)π and k10a = k11a for α, k11b = (54/9)π and k10b = k11b for β and k11c = π and k10c = k11c for γ. The effectve transton temperatures are set to σm j = 0.55 for all famly of planes n all phases. The concentratons on the sothermal phase dagrams are read n a Cartesan coordnate system. fcc phase, and an ntermedate phase θ whch has a tetragonal crystal structure. The eutectc system of (Al)-θ has a small solublty for Mg, however past the maxmum solublty lmt, there exsts other ntermedate phases termnatng at the cubc β-phase. The equlbrum lattce constants (and thus the postons of the recprocal space peaks) of θ and β phases are determned by nterpolatng between those of Al wth 3.5 at.%cu, and Al wth 38.5 at.%mg, respectvely. For smplcty, we assume a square structural symmetry for all three equlbrum phases, and lke the precedng secton, the effectve lattce constant s nterpolated by weghtng by local solute compostons, ccu and cm g. The parameters (η, χ, ω), along wth the peak wdths αj are chosen to gve a satsfactory mappng of the solublty lmts of the (Al)-phase for Cu and Mg to those n the expermental phase dagram of Fg. (b), for a range of temperature parameters (σ). The full lst of parameters used to construct the phase dagrams n ths subsecton are lsted n the capton of Fg.. Fgure (a), shows the free energy landscape for the sold at σ = Fgure (c) shows the correspondng sothermal phase dagram at σ = 0.04, where the nset shows a zoomed n mage of the Al-rch corner. Comparng the nset wth the expermental phase dagram, reasonable agreement s evdent between the calculated and the expermental phase dagram sectons. Fgure (d) shows the sothermal phase dagram for σ = At ths hgher temperature (stll below the eutectc), there s an ncrease n the solublty lmts of the phase boundares. Secton IV B wll use ths phase dagram to demonstrate sold-state precptaton. IV. APPLICATIONS The bnary XPFC approach was prevously demonstrated as a tool wth whch to model the role of defects and elastcty n structural phase transformatons that operate over dffusve tme scales. Further to these capabltes, the ablty to have mult-component nteractons between solute atoms and defects now makes t possble to examne much more complex nteractons of the above atomc-scale effects wth dfferent solutes, and ther dffuson. Ths capablty opens a myrad of possbltes for applcatons for mcrostructure engneerng n materals. Ths secton showcases some applcatons of the XPFC mult-component model presented n ths work.

8 8 FIG. : (Colour onlne) Al-Cu-Mg Phase dagram: (a) Sold and lqud energy landscapes of a square-square-square ((Al) β θ) system at temperature σ = 0.04, (b) The Al-rch sde of an sothermal cut (at 400 C) from the expermental phase dagram of the Al-Cu-Mg system taken from Ref. [41]. Dashed crcles mark the regons of the Al-rch (Al), Cu-rch (θ) and Mg-rch (β) regons consdered for reconstructon by the model phase dagram. Reconstructed phase dagrams at temperatures (c) σ = 0.04 and (d) σ = The parameters for deal free energy and entropy of mxng were η = 1.4, χ = 1, ω = 0.005, cocu = and com g = Wdths of the correlatons peaks are α11 = 0.8 and α10 = α11 for all phases. The peak postons are k11(al) = π, k10(al) = k11(al), k11θ = (.08)π, k10θ = k11θ, k11β = (1.8765)π and k10β = k11β. For all famly of planes, σm j = 0.55, n all phases. The maxma n concentratons ccu and cm g are rescaled from unty, 1, to correspond to the Cu and Mg-content n the θ-phase and β-phase gven by the expermental phase dagram,.e., 3.5 and 38.5 at.%, respectvely. The concentratons on the sothermal phase dagrams are read n a Cartesan coordnate system. In partcular, usng the phase dagrams from the prevous secton, we demonstrate dendrtc soldfcaton and precptaton n the presence of ternary components. These phenomena are paradgms of mcrostructure evoluton of relevance to materals engneerng applcatons and are strongly nfluenced by dffuson of mpurtes, elastc stran, crystal ansotropy and defect structures. A. Dendrtc Soldfcaton Dendrtc soldfcaton arses when a supercooled lqud s quenched nto the sold-lqud coexstence part of the phase dagram. Fgure (3) shows snapshots n tme of a dendrtc crystal n a ternary alloy. The smulaton was done usng the phase dagram n Fg. 1(b). Smulatons were conducted n a D doman of sze 768a 768a, where a s the lattce spacng. A unform grd spacng and dscrete tme of x = 0.15 (whch makes for a doman of sze grd ponts) and t = 3 were used and equatons of moton, Eq. (0), were solved sem-mplctly n Fourer space. The ntal condtons conssted of a small crcular seed of dameter, d = 8a of γ-phase, seeded n lqud at a temperature of σ = The ntal concentraton of solute components A and B was unform n both phases and set to the values c A = 0.1 and c B = 0.1. Several tme slces of the smulaton doman, show-

9 9 (a) (b) (c) FIG. 3: (Colour onlne) Early-tme dendrtc soldfcaton n a ternary alloy, smulated usng the phase dagram of Fg. 1(b). The quench temperature s σ = 0.18 and the ntal solute compostons are unform and set to the alloy averages, c A = 0.1 and c B = 0.1. Each Column of mages represents a dfferent tme durng the smulaton. The tmes shown are: (a) 1000 (b) 3000 and (c)7000 teratons. From bottom to top, each row dsplays the progresson of n, cb and ca, respectvely, wth ca plotted n the colour range from whte (lowest concentraton) to red (hghest concentraton) and cb s plotted n the colour range from whte (lowest concentraton) to blue (hghest concentraton). FIG. 4: (Colour onlne) Dendrtc soldfcaton at tme, t = 10000, dsplayng the values of all three felds. The top left and bottom rght quadrants show the densty feld, n. The nset on the bottom rght, s a zoomed-n mage of the rectangular area marked n black, revealng the structure of the underlyng square lattce and densty wave structure through the nterface. Top rght quadrant shows ca, wth colour range whte (low concentraton) to red (hgh concentraton), whle bottom left shows cb, wth colour range whte (low concentraton) to red (hgh concentraton). Both show the hgh solute content at the nterface. the dynamcal equatons. ng the felds (n, ca and cb ) at early tnes, are shown n Fg. 3. As tme progresses durng the smulaton, Fg. 3(a)-(c), dendrtc growth s evdent. The crystal develops a characterstc 4-fold symmetry of the underlyng square crystal structure, produced wth the correlaton functon for the gven pure component of the γ-phase. The top two rows show the tme evoluton of the concentraton felds (from left to rght), ca and cb, respectvely, ndcatng the nterface boundary layer for each component. Both solutes, A and B, reach ther maxmum solute content at the nterface of the growng dendrte, n agreement wth the solute rejecton mechansm of crystal growth. The bottom row shows the evoluton of the densty. There s also evdence of the assocated densty jump at the nterface between sold and lqud phases as depcted by the lght halo lke regon around the nterface. Fgure 4 shows a composte vew of the dendrte at later tme, hghlghtng n each quadrant one of the three felds. Ths smulaton depcts multple dffusng speces, densty changes and surface tenson ansotropy. In a larger numercal doman (where multple dendrtes can be grown), gran boundares would also naturally emerge. It s noteworthy that these physcal ngredents arse self-consstently and are very straghtforward to smulate numercally. We also note that sde-branchng of the growng dendrte s not observed n Fg. 3 due to the sze of the smulaton doman and the excluson of thermal nose n B. Solute Clusterng and Precptaton Many propertes of engneerng alloys are typcally attaned through downstream processng followng soldfcaton. These downstream processes typcally nvolve ether thermo-mechancal manpulatons or heat treatment of the ascast mcrostructure. One of the most mportant ams s to nduce certan phase transformatons n the as-cast prmary sold matrx to help strengthen alloys, a process known as precptaton hardenng. In ths subsecton we demonstrate ths process usng the ternary XPFC model developed n ths work. In partcular, we llustrate the ntal stages of a heat treatment process leadng to solute clusterng, the precursor stage of precptaton n Al-Cu-Mg alloys. The detals of ths process have been reported elsewhere [40]. Solute clusterng/early-stage precptaton smulatons were preformed usng the equlbrum propertes calculated for the (Al) β-θ system n Fgs. (c) and (d). Smulatons were performed on a D rectangular mesh wth grd spacng x = 0.15 and tme step t = 10. Dynamcal equatons were solved sem-mplctly n Fourer space. Intal condtons conssted of dstorted sngle-phase structures, through the ntroducton of a unform dstrbuton of dslocatons, and a un-

10 10 FIG. 5: (Colour onlne) Tme evoluton of clusters n solutonzed/quenched (a)-(c) Al-1.1Cu and (d)-(f) Al-1.1Cu-0.Mg alloys at σ = The nsets n (a) and (d) show the ntal dstorted/damaged sngle-phase structures, wth dslocatons clearly marked, for each set of smulatons. form composton everywhere of ccu = 1.1 and cm g = 0. at.%. All smulatons were ntally solutonzed for some tme at σ = 0.155, followng whch they were quenched/aged at a temperature σ = Durng ageng, small clusters ntally appear wth hgher Mg and/or Cu-content than that of the matrx. As tme progresses, some of these clusters decrease n sze and Cu-Mg-content, or vansh entrely. A few, however, stablze, as shown by the typcal stablzed clusters a and b n Fg. 5(a)-(c) and (d)-(f) for Al-1.1Cu and Al-1.1Cu0.Mg alloys, respectvely. In contrast, for ether alloy, when we ncrease the ageng temperature wthn the sngle-phase (Al) regon, e.g., σ = 0.145, no clusterng s observed and the ntal dstortons are removed from the matrx. Experments n quenched/aged Al-Cu and Al-Cu-Mg alloys [4 44] have found ncreasng evdence that the nteracton of ternary mpurtes and quenched-n defects such as dslocatons [45] dynamcally reduce the local nucleaton barrer for precptaton at locatons n the matrx. We have also found that the addton of Mg nto an Al-1.1Cu alloy promotes clusterng and refnement of the fnal mcrostructure, as seen n the smulaton data of Fg. (5). The clusterng phenomenon observed n these smulatons can be attrbuted to the propensty for solute segregaton to defects and surroundng areas to releve stresses nduced by the presence of sad defects, n ths case dslocatons. As more solute aggregates to dslocatons, the sze of the cluster ncreases but the structural nature of the cluster also begns to approach that of the next nearest stable sold phase. As ths process contnues and the ever growng cluster attracts more solute, t creates addtonal stresses n the surroundng matrx. Ths n turn draws nearby dslocatons to the cluster n attempts to releve these addtonal stresses caused by solute accumulaton. An extensve nvestgaton of solute clusterng mechansms, n presence of quenchedn bulk crystal defects, has been done through a quanttatve analyss of the system energetcs n bnary alloys n [34] and recently n ternary alloys, usng the present model [40]. V. SUMMARY Ths paper reported a new phase feld crystal model for structural phase transformatons (XPFC) n mult-component alloys. The detals of the model dervaton were dscussed. A smplfed verson of the model was specalzed for ternary alloys and ts equlbrum propertes were shown. The dynamcs of the model were demonstrated on two phenomena of relevance to mcrostructure evoluton n materals scence. Ths s the frst mult-component PFC model, and as such s able to capture the complex knetcs of soldfcaton and elastc and plastc effects on sold state processes, such as clusterng and precptate growth. Ths model has been used n a separate work [40] to support recent experments on the elusve

11 11 mechansms of the early stages of clusterng and precptaton. The phase feld crystal methodology was ntroduced to create a brdge between the atomc and tradtonal phase feld regmes. As a relatvely novel method, many works n ths area of materals scence are workng to valdate the physcs of PFC models. As the frst phase feld crystal model for N-component alloys, ths work has demonstrated some mportant thermodynamc and knetc propertes of the model. Moreover, asde from the model s quanttatve and self-consstent nature, t s partcularly smple to operate numercally. It s expected that ths model can thus be used to elucdate the role of multple solutes n phenomena governed by atomc-scale elastcty and defects operatng on dffusonal tme scales. Acknowledgments We thank the Natural Scence and Engneerng Research Councl of Canada (NSERC) and Ontaro Mnstry of Research and In- novaton (Early Researcher Award Program) for fnancal support, and the CLUMEQ Supercomputng faclty of Compute Canada. Appendx A: Long Wavelength Lmt In secton II B, we reduced the free energy n Eq. (7) nto a smplfed mult-component XPFC energy functonal. In the process of dong ths, we smplfed terms by consderng the long wavelength lmt where the concentraton vares much more slowly than the densty feld. Ths Appendx detals the steps of how some terms of Eq. (7) can be smplfed to derve the smplfed free energy functonal n Eq. (9). 1. Terms Couplng Product of c and c j wth C j We begn frst wth terms nvolvng a couplng of two concentraton felds wth a correlaton functon. As a concrete example, consder the term G = 1 dr c (r) ( r r )c j (r ), (A1) n Eq. (7), where we have used the more explct notaton for clarty. (The other terms follow analogously.) To proceed, we rewrte the correlaton functon n a Fourer seres of the form, C j ( r r ) = Substtutng Eq. (A) nto Eq. (A1) yelds, G = 1 dr c (r) dk Ĉj ( k )ek r e k r. (A) dk Ĉj ( k )ĉ j(k)e k r, (A3) where we defne ĉ j (k) dr c j (r )e k r. (A4) Consderng the long wavelength lmt, we take a Taylor seres expanson of the correlaton functon n powers of k around k = 0. Ths results n, G = 1 dr c (r) dk ( 1) l (k ) l l Ĉ j l! (k ) l l=0 (A5) ĉ j (k)e k r. k=0 We note that to nvoke the long wavelength lmt, we could have also Taylor expanded the concentraton, c j (r ), at r = r as s done n Refs. [5 7] or employed the mult-scale expanson used n Refs. [37, 38]. All these methods, though dfferent and requre dfferent mathematcal treatments, are found to be equvalent. Retanng, to lowest order, terms up to order l = 1, we have G = dr c (r) dr c (r) dk Ĉj ( k ) k=0 ĉ j (k)e k r dk k Ĉj (k ) ĉ j (k)e k r. (A6) k=0 Usng the defnton of the nverse Fourer transform, we recast Eq. (A6) as G = γ j dr c (r) c j (r) j κ j dr c (r) ( ) c j (r), where we have used the followng defntons, and (A7) γ j Ĉj ( k ) k=0 (A8) κ j Ĉj (k ). k=0 (A9) It s thus clear that the frst term n Eq. (A7) wll contrbute terms that renormalze the coeffcent of the c terms n the entropy of mxng, f Eq. (9) were expanded about c = c o. In ths work, the γ j terms are neglected, and ther role s subsumed n an effectve manner, for convenence, nto the prefactor ω n Eq. (9). The second term n Eq. (A7) can be recast nto gradent energy terms analogous to those used n Cahn- Hllard or Gnzburg-Landau theores. To do so, we perform

12 1 ntegraton by parts, yeldng, G = γ j dr c (r) c j (r) κ j dr c (r) c j (r). (A10) Lastly, for clarty of exposton, we separate the gradent terms n Eq. (A10), whch yelds G = γ j dr c (r) c j (r) κ dr c (r) κ j dr c (r) c j (r). (A11) The second term of Eq. (A11) gves gradents terms of the Cahn-Hllard form, whle the thrd lne yelds cross terms. In ths work, for smplcty, we set κ j = 0 for j. We note that such cross terms can become mportant when studyng certan phenomena and/or when hgher-order alloyng nteractons are consdered.. Correlaton Kernels Contanng Lnear Terms n n To demonstrate the long wavelength lmt of terms lnear n densty n Eq. (7), we consder, as an example, the term H= 1 dr c o c j (r) ( r r )n(r ). (A1) Substtutng the Fourer seres expanson of the correlaton functon, Taylor expandng the correlaton as n Eq. (A5) (retanng the lowest order term), and takng the nverse Fourer transform yelds, H = γ j dr c o c j (r)n(r) κ j dr c o c j (r) ( ) n(r), (A13) where γ j and κ j are defned by Eqs. (A8) and (A9), respectvely. The densty, n(r), n Eq. (A13) s rapdly varyng. Its leadng order representaton s defned by a sngle-mode approxmaton of the form n(r) = m A m (r) e qm r, (A14) where q m are the recprocal lattce vectors and A m (r) are slowly varyng ampltudes correspondng to each recprocal lattce vector, m. Substtutng Eq. (A14) nto Eq. (A13) gves H = 1 1 γ dr c o c j (r)a m (r) e qm r m m κ dr c o c j (r) ( A m (r) e qm r). Expandng the Laplacan n Eq. (A15) gves, H = 1 1 γ dr c o c j (r)a m (r) e qm r m m κ dr c o c j (r) e qm r L m A m (r), (A15) (A16) where L m + q m q m s a covarant operator that assures rotatonal nvarance of the free energy n the long wavelength lmt. It s noted that each term n Eq. (A16) only contans one rapdly oscllatng varable,.e., e qm r. If we apply the so-called quck and drty [46] analogue of the volume averagng method employed n Refs. [6, 7] (whch amounts to decouplng slowly varyng felds nsde ntegrals from rapdly varyng phase factors, thus makng the ntegrals effectvely vansh when ntegrated over one unt), we obtan H 0. It s straghtforward to show that all other terms n Eq. (7) that are lnear n n, such as, H = 1 dr n(r)c (r) ( r r )c j (r ), (A17) smlarly vansh upon coarse granng. It should also be evdent from the above consderatons that f Eq. (A1) contaned an n(r) n(r ) combnaton, then Eq. (A16) would contan terms wth phase factors of dfferent combnatons of sums of two recprocal lattce vectors. Some of these two-vector combnatons would add up to zero causng ther correspondng terms to survve upon ntegraton. 3. Volume Averagng Equaton (A16) can more formally be analyzed usng a volume averagng convoluton operator [6], defned by f(r) V 1 πv dr f(r )χ V (r r ), (A18) where f(r ) s the functon beng course graned and V s the coarse granng volume. The functon χ V n the ntegrand of Eq. (A18) s a smoothng functon that s normalzed to unty,

13 13.e., dr χ V (r r ) 1. A convenent form of χ V s gven by (A19) χ V (r r ) = 1 πv e (r r ) V. (A0) In the long wavelength lmt, L c L a where L V 1/d, n d-dmensons, whle L c s the length scale of varaton of the concentraton feld. Ths condton mples that the functon χ V (r) vares on dmensons much larger than the lattce constant a = π/ q m but much less then the length scale of varaton of the concentraton, c (r). Equaton (A18) defnes a nonnvertble lmtng procedure that can be used to average a functon over some volume. It s nstructve to apply the volume averagng procedure to the frst term n Eq. (A16). For convenence we defne φ(r) c o c j(r)a m (r). It s noted that φ(r) vares on scales much larger than the lattce constant snce t s comprsed of slowly varyng functons. Usng the defnton of φ(r), the frst ntegral of Eq. (A16) can be wrtten as H V = 1 = 1 ( γ dr m γ m ( dr ) dr χ V (r r ) φ(r )e qm r dr χ V (r r )φ(r ) e qm r ), (A1) Snce φ(r ) vares more slowly than the scale of varaton of χ V, t s reasonable to expand t n a Taylor seres about r = r. Substtutng φ(r ) = φ(r) φ(r) (r r ) nto the above expresson leads to ( dr φ(r) dr χ V (r r ) e qm r H V == 1 γ m φ(r) dr χ V (r r ) (r r ) e qm r + ). (A) The nonnvertble procedure was ntroduced n the second lne of Eq. (A1). In the long wavelength lmt, when q m L, both ntegrals n Eq. (A) vansh as ( q m L ) 1, makng H V smlarly vansh. Appendx B: Phase Dagram Calculaton Equaton () provdes a system of equatons that are exact when one needs to determne the equlbrum propertes of a gven system. For a bnary system, where the number of equatons n Eq. () s reduced by one, at a specfed temperature and pressure they are suffcent to specfy exactly the unque phase concentratons correspondng to the te lne between two phases. However for mult-component systems, for the present case of a ternary alloy (represented by solute compostons A and B), the set of condtons n Eq. () are under-determned and cannot unquely defne all phase concentratons. Ths s because, for a ternary system at a specfed temperature and pressure, there s not generally a sngle te lne whch specfes phase boundares between coexstng phases but, rather multple te lnes defnng the boundary between any two phases. The under-determned set of condtons n Eq. () contan varables c I A, cj A, ci B and cj B n phases I and J respectvely. To close ths system, an addtonal condton s necessary to provde a fourth equaton relatng the concentratons. A convenent fourth condton s the lever rule, whch relates weght fractons of phases to the average concentraton. For clarty, we specfy t here for ternary sold (α) and lqud (L) phases, and c A = c L Ax L + c α Ax α c B = c L Bx L + c α Bx α, (B1) (B) where c A and c B are the average alloy compostons for components A and B respectvely and x L and x α represent the equlbrum volume fractons of lqud and α respectvely, and satsfy x L + x α 1. Combnng ths last relaton between the volume fractons and Eqs. (B1) and (B) gves, the last equlbrum condton, c A c α A c L A cα A = c B c α B c L B. (B3) cα B Equaton () together wth Eq. (B3) comprses a complete set of equatons whch can admt unque te lne solutons,.e., solutons for c L A, cα A, cl B and cα B n the sold-lqud example just consdered. Wth the free energy functons generally beng hghly nonlnear, t s not possble to fnd analytcal solutons to Eqs. () and (B3), and they must be solved numercally. One approach s to specfy the temperature and then raster through the phase space of average concentratons c A and c B, where the rasterng s done by takng dscrete steps n steps of c A and c B, respectvely (for practcal purposes ts convenent to set c A = c B = c). For each par of c A and c B, Eqs. () and (B3) can be solved numercally. The solutons yeld c L A, cα A, cl B and cα B. A unque soluton for each par of c A and c B defnes one te lne. The collecton of all such te lnes maps out the coexstence phase boundares between any two phases, n the case consdered here, L and α. Where no solutons are admtted correspond to sngle phase regons where no te lnes exst. It s expected that the smoothness of the phase boundares, when plotted for graphcal purposes, wll depend on the step sze, c, chosen to dscretze the average concentraton values. The above mentoned recpe can stll requre ntensve computaton, requrng a soluton of four equatons n four unknowns for M combnatons of average concentraton pars (where M s the dscretzed number of average concentraton