Morphological Stability during Electrodeposition

Size: px
Start display at page:

Download "Morphological Stability during Electrodeposition"

Transcription

1 C708 Journal of The Electrochemical Society, C708-C /2003/15010/C708/9/$7.00 The Electrochemical Society, Inc. Morphological Stability uring Electroeposition II. Aitive Effects Mikko Haataja, a,z Davi J. Srolovitz, a an Anrew B. Bocarsly b, * a Princeton Materials Institute an Department of Mechanical an Aerospace Engineering an b Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA Common experience shows that electroeposite ED metallic films exhibit rough surfaces unless the electrochemical bath contains small quantities of molecular aitives. We evelop a formalism for escribing the effects of aitives on surface morphology evolution, which buils on that in a companion paper for the aitive-free case. We emonstrate that the aitives suppress the morphological instability that leas to roughening by preferentially accumulating near surface protrusions an blocking growth sites. Our chemically base moel shows that aitives which reaily asorb onto the surface an have a strong tenency to complex with the metal cations reuce the riving force for the instability an thus enhance leveling. Furthermore, polar aitives provie an aitional stabilizing effect, in accor with experimental observations. It is also shown that linearly stable growth can be achieve over a wie range of eposition flux at sufficiently large aitive bulk concentrations. We preict the ED conitions necessary for growing flat films an emonstrate that these are in goo agreement with experimental observations The Electrochemical Society. DOI: / All rights reserve. Manuscript submitte January 28, 2003; revise manuscript receive April 7, Available electronically August 18, Extensive experience shows 1,2 that the properties of electroeposite ED metallic films can be effectively controlle by introucing small quantities of aitives into the electrochemical bath. For example, morphological instabilities uring eposition commonly prouce poor quality films with rough surfaces 3-7 in the absence of aitives. Aition of even small quantities of aitives can rastically ecrease the amplitue of the surface roughness as well as increase the characteristic wavelength of the roughness. 7 In aition to controlling the large-scale roughness leveling, aitives have been use to control crystallographic texture an grain size of ED films. The goal of the present paper is to provie a formalism which can be use to evaluate propose mechanisms for surface leveling by the introuction of aitives into the electrochemical baths. We apply this formalism to explicitly evaluate the effects of complex formation between cations an aitives, aitive ipole moments, surface segregation, an aitive incorporation into the growing electroeposit. Typically, aitives consist of organic compouns or small inorganic species that can act as ligans with respect to the metal being eposite. The ionic nature of the aitive along with the exact moe of ligation is extremely variable. Commonly employe aitives can be organic e.g., thiourea, CH 4 N 2 S in nickel plating baths or inorganic e.g., potassium nitrate, KNO 3 in silver plating baths, ionic e.g., soium sulfocyanate, NaSCN in silver electroeposition, or nonionic, polar e.g., naphthoquinone, C 10 H 6 O 2 in lea electrorefining or nonpolar e.g., anthraquinone, C 14 H 9 NO 2 in lea electrorefining. 2 This suggests that the leveling effect of the aitives epens on only a few key physical characteristics. While there are many commercial levelers, a complete unerstaning of their effects an the physical an chemical mechanisms that allow these molecules to effectively reuce the roughness of films remains elusive. This is in contrast with the emerging unerstaning of the role of aitives in vapor phase eposition; e.g., the surfactant effect in semiconuctor growth see Ref. 8 an references cite therein. In Part I 9 we presente a formalism that provies quantitative preictions for the evolution of the surface morphology in the absence of aitives. In this paper we exten the formalism of Part I to the case where ED occurs in the presence of aitives. We present a physically an chemically motivate moel for ED with aitives, then investigate the influence of the aitives on the morphological evolution an stability of surfaces uring ED. A brief account of this work has appeare in Ref. 10. * Electrochemical Society Active Member. z mhaataja@princeton.eu Morphology evolution in the presence of aitives has been previously aresse through both experiments 2,7 an theory Qualitatively, there is some experimental evience suggesting that polar aitives or aitives which ten to form complexes with the metal cations lea to more effective leveling. 1,2 An experimental stuy of the growth of copper films from a copper sulfate bath prouce increasingly rough surfaces uring growth in the absence of aitives, while aing thiourea le to smoother films. 7 Increases in thiourea concentration le to smoother films an longer wavelength roughness. A range of thiourea concentration exists for which the growth front remains essentially flat. We argue that this is attributable to the complete suppression of growth instability for aitive concentrations larger than a critical value. Several moels have been propose to explain the leveling effects of aitives. They all assume that the aitives asorb onto the growing surface an interfere with the eposition process through a specific mechanism, such as blocking the surface growth sites; 1,2,11 the leveling effect results from preferential growth at surface epressions instea of protrusions. Maore et al. 11 introuce a continuum moel which inclue the blocking effect of charge aitives as well as aitive reuction at the electroe, couple to iffusive transport of the aitives to the surface. They foun that the leveling effect of aitives is controlle by the ratio of the aitive iffusivity to the rate of aitive consumption at the cathoe. Georgiaou et al. 12 stuie the evolution of film morphology uring ED of copper in microtrenches by employing a moel where the aitive suppresses metal eposition an is consume by an electrochemical reaction. By numerically solving the full nonlinear moving bounary problem they emonstrate that the aitives promote voi-free filling of trenches. Similarly, Cao et al. 13 numerically solve a moel in which the asorption of aitives locally inhibits eposition an showe that this leas to effective leveling an voi-free filling of trenches. Recently, Pricer et al. 14 introuce a couple Monte-Carlo/continuum approach for stuying the effect of a blocking aitive in copper ED in rectangular trenches. Finally, Josell et al. 15,16 have recently introuce an accelerator-aitive moel which quantitatively explains superconformal ED in submicrometer features. While the aforementione theoretical approaches have successfully explaine some of the experimental results, they have neglecte several common experimental observations an have base their moels on an oversimplifie picture of transport in the electrolyte. For example, effective aitives are known to ligan to the electroactive metal ions, forming complexes in solution. The aitives alone or in a complex can prouce neutral species which interact with the electric fiel through strong molecular ipole mo-

2 Journal of The Electrochemical Society, C708-C C709 ments. We consier several types of aitive/metal ion interactions an explicitly account for metal cation an aitive rift in the electric fiel within electrolytes that are not fully supporte. Incluing the effects of spectator ions allows us to interpolate between the limiting cases of fully supporte purely iffusive transport an unsupporte electrolytes. The resultant linear stability theory allows us to examine the interplay between aitive an metal ion transport, rift, an iffusive transport an make preictions about the resulting morphology in terms of the external control parameters e.g., aitive bulk concentration, spectator ion concentration, eposition flux an the chemical properties of the metal cations an aitives e.g., charge, ipole moment, heat of segregation, an complex formation. In the next section we very briefly review the relevant aspects of the formalism an the results on growth instabilities in the absence of aitives from the companion paper, Part I. 9 We then escribe several physical an chemical moels for the mechanisms by which the aitives influence the eposition process. Then we perform a self-consistent, uniform steay-state analysis of ED in the presence of aitives incluing the concentration an electrostatic fiels. Morphology evolution is analyze using perturbation theory, where we linearize aroun the uniform steay-state result. The resultant preictions for the wavelength an growth rate of the roughness as a function of chemical an eposition parameters are collecte in the form of a stability map, which shows uner which conitions the surface instability can be entirely suppresse. Finally, we compare our preictions with experimental observations. Moel for Electroeposition in the Absence of Aitives Formalism. We briefly review the relevant aspects of the morphology evolution formalism put forth in Part I. The number ensity for each ionic species obeys a continuity equation C i t j i 0 where i C, A, XC, XA enotes the metal cations, anions, spectator cations, an spectator anions, respectively. The ion fluxes j i (r) follow from 1 j i r D i C i q ied i C i 2 where D i enotes the iffusivity an is the electrostatic potential. is etermine from the Poisson equation 17 2 e ˆ i q i C i r 3 where ˆ enotes the constant local permittivity of the ionic solution. We consier a rectangular system of imensions L x L z W L, where W enotes the linear imension with of the planar surface an L enotes the thickness of the mass-transfer layer over which the concentrations vary. Beyon L the concentrations of all species are very nearly constant. Perioic bounary conitions are employe in the x irection, an the bounary conitions in the z irection become C C (x, L) C A (x, L) C 0, C XC (x, L) C XA (x, L) C 1, an (x, L) 0. At the growing surface (x, 0) V ext, an the ion fluxes i (A,XC,XA) vanish as these species o not coeposit: j i j i nˆ 0, where nˆ is the surface normal pointing into the bath. The magnitue of the local metal cation flux onto the growing surface j C taken to be positive for a growing film is given by the Butler-Volmer B-V equation j C j C nˆ j 0 C C C 0 e 1 F/RT e 2 F/RT 4 where F an R are the Faraay s an gas constants, an ( 1, 2 ) enote the so-calle symmetry factors relate to the potential barrier for metal cation reuction. 18 The overpotential is given by 19 V eq V ext RT F ln C C C C eq F qc The equilibrium potential of the metal-solution interface is V eq, an varying the external electroe potential V ext away from V eq leas to either eposition ( 0) or issolution ( 0). The thir term in Eq. 5 accounts for the metal cation concentration epenence of V eq, an the last term escribes the effects of surface curvature an surface tension, i.e., the Gibbs-Thomson effect. The quantity q C ej 0 enotes the exchange current ensity, an enotes the atomic volume of the metal in the eposit. Finally, the local growth velocity in the normal irection nˆ, v n, follows from the massbalance relation v n j C nˆ. Morphology evolution. In Part I we introuce a physically motivate continuum moel for the morphological stability of surfaces uring ED in the absence of aitives. The moel explicitly accounts for the electric fiel in the electrolyte, the metal cations an anions, an any aitional ions from a supporting electrolyte. By explicitly solving the steay-state equations for a planar growth front, we emonstrate that the moel naturally gives rise to the Gouy-Chapman G-C bounary layer an that increasing the concentration of the spectator ions leas to a rapi ecrease in the magnitue of the electric fiel in the bulk. In this limit the metal cation transport is ominate by iffusion in the concentration graient, as expecte. A first unerstaning of the morphology evolution was obtaine by carrying out a perturbation analysis in the surface profile. In agreement with previous stability analysis, the surface was shown always to be linearly unstable against perturbations with a sufficiently large wavelength for any finite external eposition flux. This result hols for both unsupporte an fully supporte baths an can be unerstoo as follows. Consier a small perturbation to the planar surface. Because C C increases away from the growing surface, the metal cation concentration C C is larger smaller near protrusions epressions than near the flat surface. In the absence of surface tension an increase in C C leas to a larger metal cation flux via the B-V equation, Eq. 4. The larger metal cation flux leas to an increase in v n an hence, to faster slower growth at a protrusion epression. This positive feeback leas to unstable growth of the perturbations. Capillarity, however, effectively suppresses perturbations on small scales, an these two competing effects eposition an capillarity set the lateral scale of the surface roughness observe in experiment. In particular, increasing the eposition flux makes the surface rough on smaller scales an leas to faster roughening. Interestingly, for a fixe eposition flux j the system without the supporting electrolyte was foun to be more stable than one with a supporting electrolyte, in agreement with previous stuies. 22 This was shown to be a consequence of the increasing local metal cation concentration graient as more spectator ions are ae into the bath. Moels for Aitive Effects in Electroeposition We now turn to the escription of aitives within the framework introuce in Part I. Qualitatively, there are three main questions that must be consiere in assessing the role of aitives in ED: i how o they interact with the ions in the bulk, (ii) how o they get to the surface, an (iii) once on the surface, how o they moify the eposition of the metal cations. It has been suggeste that a metal cation an an aitive molecule can form a complex in the bulk solution; 1,2 in particular, it has been note that increasing the tenency of metal ions to form complexes with the aitives correlates with their susceptibility to the effect of aitives. 1 The transport of the aitives to the surface is usually assume to be 5

3 C710 Journal of The Electrochemical Society, C708-C iffusion-limite, 2,7,14 although it is conceivable that some of the aitives are brought onto the surface through the motion of the metal cations as a part of a complex. Several moels have been propose to explain the leveling effects of aitives, the most popular of which are i blocking the surface, (ii) changes in the bounary layer potential, an (iii) ion briging. 1,2 The first mechanism, also known as the irt mechanism, is attribute to the tenency of the aitive molecules to asorb onto the surface an interfere with either the attachment of metal cations or their iffusion along the surface to preferential growth sites. This type of aitive causes a ecrease in the rate of the electroe reaction at fixe potential. The leveling effect of such aitives can be unerstoo by consiering the aitive concentration aroun a protrusion: if the protrusion collects more aitives than o epressions e.g., ue to an aitive concentration graient, it tens to grow more slowly than the epressions, implying leveling. 2,7,14 The secon mechanism is base on the iea that the reuction reaction oes not necessarily have to occur at the surface but can take place anywhere within the G-C bounary layer. 2 Therefore, any foreign ions or neutral molecules which asorb onto the surface moify the potential within the G-C bounary layer an thus influence the eposition rate. Finally, the aitives may increase the rate of reuction; 2 this coul apply to aitives which bon to the surface an meiate the charge transfer between the metal cation an the eposit ion briging. It is noteworthy that the effect of a particular aitive species can be a combination of all these mechanisms. Furthermore, there is some evience that the leveling is brought about through the interaction between reaction proucts of the aitive at the surface an the growing eposit. 2 Finally, in some cases the growth of a level an bright eposit requires the simultaneous use of multiple types of aitives. Such synergistic effects are not consiere here. In light of this iscussion it is clear that unraveling the unerlying mechanisms for the leveling effect of a given aitive molecular species is a aunting task. Therefore, the goal of the present work on aitive effects is by necessity somewhat more limite, namely, to construct a minimal moel for aitive effects, to quantify the role of ifferent interactions in the morphology evolution of the surface, an finally, to relate this to leveling. In particular, we consier aitives which i form complexes with the metal cations, (ii) interact with the electric fiel via a ipole moment, an (iii) interfere with the eposition process by blocking surface growth sites. Although some of the microscopic etails are necessarily ignore in this approach, many of the common properties of aitives are accounte for. This approach provies a means to assess the importance of ifferent types of interactions an therefore provies a guie for the evelopment of aitives. Coupling between the aitives an the eposition process. Complex formation. We assume that complex formation between the aitives an the metal cations is local an pairwise. That is, we moel complex formation between cations an aitive in a small volume of space as proportional to the prouct of the concentration of each in this volume. To obtain the contribution of complex formation to the total Helmholtz free energy, we a the contributions of each local volume element as F C C C rc I rr This type of expression is commonly use in statistical mechanics treatments of solution theory. 26 In this equation is a parameter which characterizes the strength of the local attractive interaction between the aitives an the metal cations. A simple relation between the heat of complex formation Q I per aitive molecule an the parameter can be establishe by assuming that one aitive molecule interacts with a single metal cation: Q I /C 0. In this moel any inhomogeneities in the metal cation concentration are accompanie by a variation in the aitive concentration C I (r). 6 Equation 6 accurately escribes the interaction between the aitives an the metal cations through complex formation, given that the complexes are relatively weakly boun, which implies that Q I. Polar aitives. There is experimental evience 2 that polarity of the aitives enhances the stability of the growing surface. It has been observe in lea electrorefining that naphthoquinone, which has a molecular ipole moment, is a very goo leveler, while a chemically similar complex molecule anthraquinone with very small molecular ipole moment has a very small leveling effect. Hence, it has been hypothesize 2 that the leveling capacity of naphthoquinone can be attribute to the finite molecular ipole moment which allows the aitives to interact with the inhomogeneous electric fiel at the metal/solution interface. The energy of a ipole in an electric fiel is simply the ot prouct of the ipole moment vector an the local electric fiel. 17 The total interaction energy of the ipole moment on the aitives an the electric fiel F E is proportional to the local concentration of ipoles aitives an the local electric fiel summe over all small volume elements F E C I p r C I pr rc I pr r 7 where p is the molecular ipole vector of the aitives, an the last two equalities assume that the molecules quickly orient with their ipoles along the fiel. Note that this coupling enhances the concentration of aitives in regions where the magnitue of the electric fiel is large. Transport equations an fluxes. Base on the aitional interaction terms introuce we rewrite the metal cation flux as j C r D C C C qed CC C D CC C C I 8 The last term is ue to the coupling between the aitives an the metal cations, as escribe by Eq. 6. Because the anions are assume not to form complexes with the aitives their flux is still written as Finally, the aitive flux reas j A r D A C A qed AC A j I r D I C I D IpC I D IC I C C 10 The aitive flux has contributions from iffusion, graients in the electric fiel if the molecules are polar (p 0), asorption to the surface, an complex formation with the metal cations. Diffusion tens to smooth out the inhomogeneities in the concentration profile, whereas the aitives ten to accumulate in regions of large electric fiel graients because the aitives are polar, large metal cation concentration, or both. Aitives an surface site blocking. We incorporate the effect of aitives in blocking potential surface growth sites by aopting a moifie B-V equation for the local metal cation flux j C which satisfies j C D C C C D CqeC C nˆ j 0 1 I D CC C I C FRT e 1 F/RT e 2 C 0 C C 9 11

4 Journal of The Electrochemical Society, C708-C C711 where I enotes the aitive surface coverage. We escribe the surface coverage using a Langmuir isotherm I KC I 1 KC I 12 where the asorption equilibrium constant K exp(q/) an Q enotes the heat of asorption. The tenency of the aitives to locally block some of the available growth sites on the surface thus leaing to slower eposition rates for high local aitive coverage is phenomenologically incorporate by the multiplicative term (1 I ) in the local flux. In particular, we assume that a sufficiently high aitive concentration can lea to a total saturation of the growth sites on the surface, an hence to a suppression of eposition. This moel qualitatively captures the salient features of the interaction between aitives an the metal cations at the growth surface. We also consier the possibility that the aitives on the surface may be incorporate into the growing film by coeposition along with the metal cations. This coeposition is moele by imposing a imensionless aitive flux at the surface which is proportional to the local metal cation flux: j I C I j, where the proportionality constant escribes the extent of coeposition an C I enotes the imensionless aitive bulk concentration. While this amittely constitutes an oversimplification of the true surface physics, it is consistent with experimental observations an allows us to assess the importance of coeposition on the growth of ED films, of which little is generally known. Note that 1 implies that metal cations an aitives are incorporate into the growing eposit in proportions given by their bulk concentrations. It is noteworthy that a finite ipole moment effectively increases the equilibrium constant K as the polar aitives are attracte to the surface by the large electric fiel strengths. Unfortunately, it is generally not easy to estimate the contribution of this interaction to K quantitatively because it epens on the etaile atomic structure of the surface an the structure of the aitive in question. However, as a simple estimate we consier a typical electric fiel strength of 0.3 V/nm at the interface an a typical ipole moment of 5e Å for a polar molecular species, the prouct of which is the ipole moment contribution to Q. Using these values we fin that this contribution is ev, which is larger than by a factor of orer 3 or 4 at room temperature. This translates into a change in the value of K by a factor of nearly 25, as compare with its value if no ipole moment were present. A ipole moment on the aitives may also play a minor role on the metal cation-aitive complex formation. However, in orer to assess the significance of the aforementione interactions on the growth front stability, we treat K, p, an as inepenent parameters in the remainer of the paper. Steay-state properties. We again examine the steay-state properties of the moel as a prerequisite to the morphology analysis. In steay state the metal cation flux is constant across the length of the cell, while the anion an spectator ion fluxes vanish, respectively. Therefore, upon measuring all lengths in units of the masstransfer layer thickness L an thus all wavenumbers in units of 1/L), time in units of iffusion time L 2 /D C, ionic an aitive concentrations in units of the bulk metal cation concentration C 0, electric potential in units of the equilibrium voltage V eq, an fluxes in units of the iffusion-limite flux D C C 0 /L as employe in Part I, we may write the new imensionless steay-state equations as z C C C C z C C z C I j z C A C A 0 z z C XC C XC 0 z z C XA C XA 0 z z C I C I z z C I z C C D C I C I j 2 2 z 2 C C C XC C A C XA 13 where D C I D C /D I, an where we have use two important imensionless parameters C 0 /() Q I /() ln K c an V eq p/(l) in the aitive flux. escribes the interaction strength for complex formation an escribes the interaction strength between the aitives an the inhomogeneous electric fiels through the finite molecular ipole moment. Furthermore, K c enotes the equilibrium constant for complex formation between the aitives an the metal cations. In the subsequent calculations we employ L m, appropriate for free convection electroplating conitions. Written in a imensionless form, it can be inferre from the metal cation steay-state equation in Eq. 13 that the term C C C I /z, arising from complex formation, is typically very small in magnitue (10 6 ) compare to the iffusion an rift terms 1, because the bulk aitive concentration C I 10 6 C 0. Therefore, we may rop the last term in the metal cation flux to a very goo approximation. Note, however, that the complexing term cannot be neglecte in the aitive flux an is thus retaine in the following analysis. Because the aitives now affect the metal cation concentration only through the moifie B-V bounary conition, the steay-state solutions for C C an are exactly as given in Part I in the absence of aitives C C 1 C 1 exp C 1 exp z 2 where z, j 1 * ln 1 j 2 2C 1 1 tanh1 4 e 21 C 1 z/ tanh * 4 1 ln 1 j 2 2C 1 jz 2 2C 1 15 an * V ext /V eq. This solution is accurate to lowest orer in everywhere an to lowest orer in j/(1 C 1 ) within the bounary layer 0 z an to all orers in j in the bulk. Similarly, the steay-state solutions for C A, C XA, C XC, an C I are simply C A e, C XC C 1 e, an C XA C 1 e, an C I C I e C C /z 1 D I C j D I C jz 16 appropriate for 1. We note that the aitive concentration at the surface is given by cf. Eq. 16 C I 0 C I e C C 0 0/z 1 D I C j 17 It is noteworthy that the large electric fiels in the bounary layer effectively attract polar aitive particles. In orer to get a better unerstaning of the effects of the various parameters in the aitive concentration profile C I, we plot some representative profiles in steay state in the absence of spectator ions in Fig. 1. In particular, Fig. 1a shows the equilibrium ( j 0) concentration profiles for several values of the interaction parameter that escribes the strength of complex formation. Notice how an

5 C712 Journal of The Electrochemical Society, C708-C Figure 1. Aitive concentration profiles in steay state in the unsupporte electrolyte i.e., no spectator ions for a, top left j 0, 0, 0, an several values of. b, top right Concentration profiles at finite eposition rate j 0.2, 0, 0, an several values of. c, bottom The effect of the ipole moment on the concentration profiles for j 0.2, 0, an 0. increase in the metal cation concentration in the vicinity of the bounary layer leas to a corresponing increase in C I ue to this attractive interaction. Aitionally, the aitive concentration attains a constant value away from the bounary layer. Upon applying an external eposition flux, the aitive concentration profile is influence by the metal cation concentration graient in the bulk, as shown in Fig. 1b for j 0.2, 0, 0, an several values of. In particular, an aitive concentration graient in the bulk evelops in steay state ue to the interaction between the metal cations an aitives. Finally, in Fig. 1c we plot the aitive concentration profiles in the absence of complex formation for j 0.2, 0, 0, an several values of the parameter which escribes the coupling between the electric fiels an the polar aitives. Note the accumulation of aitives within the bounary layer as a consequence of their interaction with the electric fiel. Linear stability analysis: Aitive effects. Next we apply the linear stability analysis to etermine how the surface evolves uring ED in an electrolyte containing aitives. The analysis follows the proceure escribe in Part I in the absence of aitives. In the fully supporte electrolyte case the analysis is carrie out analytically an is vali for all values of the eposition flux j, while in the partially supporte electrolyte case the analysis is carrie out as a perturbation expansion in the flux j, as outline in Appenix B of Part I in the absence of aitives. Physically, the aitives occupy surface sites where metal cations can attach, thus slowing growth. Therefore, if protrusions collect more aitives than epressions they grow more slowly, hence favoring the growth of smooth surfaces. 7,14,18 Increase aitive concentrations occur where the aitive flux to the surface is greatest. This flux is controlle by i aitive-metal cation complexing represente by in which the metal cations rag the aitives to the growth front uring eposition, (ii) interactions between the polar aitives an the inhomogeneous electric fiel associate with the surface perturbation, an by (iii) aitive coeposition represente by which gives rise to an aitive concentration graient. The linear stability analysis shows that the perturbation growth rate in the presence of aitives (k) in a fully supporte electrolyte is given by

6 Journal of The Electrochemical Society, C708-C C713 Figure 2. Linear ispersion relation (k)/ vs. k for a fully supporte electrolyte with a, top left j 0.5, 4.0, 0.0, 1 0.5, an ˆ an several values of KC I an b, top right j 0.5, KC I 0.25, 0.0, 1 0.5, an ˆ an several values of the parameter escribing complex formation. Linear ispersion relation (k)/ vs. k for an unsupporte electrolyte for c, bottom left j 0.5, KC I 0.5, 0.0, 1 0.5, an ˆ an several values of the parameter an, bottom right j 0.5, KC I 0.25, 0.0, 1 0.5, ˆ 0.001, an 1.0, an several values of the parameter which escribes the ipole moment of the aitives. k jk j eff1 1 1 j eff j1 1 1 j 2 ˆ 1 k j1 D C I j k 18 where the effective current j eff (1 1 )/(1 j eff ) j(1 1 )/(1 j) (j D I C j j 2 D I C ), an the aitive effects are conveniently expresse by as KC I e 1j 1 KC I e 1j 1 D I C j 19 Qualitatively, Eq. 18 is very similar to the corresponing result in the absence of aitives, Eq. 17 in Part I. In particular, the growth rate (k) is positive for small k an becomes negative for k k 0, implying stability on scales s 2/k 0, where s enotes the stabilization length. Equation 18 has a very appealing physical interpretation. The presence of aitives gives rise to an effective flux j eff, the magnitue of which epens on the aitive bulk concentration C I, the strength of complex formation, the egree of coeposition, an the susceptibility to segregate onto the surface of the growing film, escribe by the equilibrium constant K. Because j eff is always less than j, the metal cation flux riven instability is reuce by the presence of the aitives. This follows from the key observation in Part I, i.e., that the morphological instability in the absence of aitives is ue to the finite eposition flux j, an ecreasing j makes the surface smoother. In Fig. 2a we plot (k) vs. k, several aitive concentrations C I, an a fixe metal cation eposition flux j. The other parameters were set to j 0.6, 4.0, 0.0, 1 0.5, an ˆ Note that increasing C I leas to reuce k max an instability growth rates, (k). Thus, increasing the aitive concentration increases the range of surface stability implying leveling for the same external eposition flux j. Closer examination of Eq. 18 also reveals that increasing the egree of aitive coeposition (

7 C714 Journal of The Electrochemical Society, C708-C ), the equilibrium aitive surface concentration represente by K, an the egree of coupling between aitives an metal cations all enhance stability, consistent with experimental observations. 2 This is illustrate in Fig. 2b, where we show (k) vs. k for several values of the interaction parameter with j 0.5, KC I 0.25, 0.0, 1 0.5, an ˆ We note that the aitives can increase the surface stability even without aitive incorporation into the growing film as clearly seen from this figure. Let us now explicitly compare our preictions for surface stabilization with existing experimental ata. In particular, in Ref. 7 the growth of copper eposits from a copper sulfate bath in the presence of thiourea was stuie. For example, employing the current ensity J 200 A/m 2 an cation bulk concentration C M from Ref. 7, an employing the reasonable values D C 10 9 m 2 /s, 1.6 J/m 2, K 10 6 I, 1 1, an 0.5, an experimentally etermine value , c appropriate for thiourea-cupric complexes, we fin that an aitive bulk concentration of mm leas to the stabilization length 0.1 mm s 0.2 mm, in goo quantitative agreement with experimental result of s 0.5 mm in Ref. 7. It is noteworthy that the parameter which escribes the strength of the interaction between the polar aitives an the electric fiel oes not explicitly appear in Eq. 18. This is a irect consequence of the fact that for a sufficiently large spectator ion concentration, the electric fiel is almost completely screene in the bulk, an therefore it is practically unaffecte by the surface perturbation. This implies that the electric fiel aroun a protrusion is the same as for a flat part of the surface or a epression, an therefore the interaction between the ipole moment on the aitives an the perturbe electric fiel oes not moify the aitive efficacy in stabilizing the surface in the fully supporte electrolyte. Increasing the aitive ipole moment effectively increases the equilibrium constant K as the aitives are attracte to the surface by the large electric fiel strengths, an this helps stabilize the surface, as iscusse previously. In the partially supporte electrolyte case, the polar nature of an aitive can play a further stabilizing role, as iscusse shortly. We have also performe a stability analysis for ED in the presence of aitives for the case in which the electrolyte is unsupporte or only partially supporte using the same approach as outline in the previous section an in Appenix B. Figure 2c shows the perturbation growth rate (k)/ as a function of wavenumber for the case of the completely unsupporte electrolyte (C 1 0), for j 0.5, KC I 0.5, 0.0, 1 0.5, ˆ 0.001, an several values of. As in the fully supporte electrolyte case, increasing is beneficial for stability. Increasing the bulk aitive concentration an the surface coverage of the aitives through increasing the equilibrium constant K also increases the surface stability not shown, as in the fully supporte electrolyte case. A novel feature arises in the unsupporte or partially supporte electrolyte case when the ipole moment of the aitives is finite, as anticipate previously. During eposition, a finite electric fiel is set up in the bulk; the presence of this fiel is ue to the absence of spectator ions an it persists for small spectator ion concentrations, as shown in Part I. Consier next the nonzero electric fiel associate with a surface protrusion. Since the growing surface is an equipotential one, it follows that the equipotential lines are crowe together ahea of the protrusion, an the corresponing electric fiel becomes large this is shown explicitly in Part I, Appenix B; the opposite is true for a epression. Because the polar aitives are attracte to regions of large electric fiels, a protrusion collects more aitives, which tens to stabilize the surface. These physical arguments are confirme by calculating the linear stability curves for an c Doona an Stanbury 27 measure the equilibrium constant K c for thiourea-cupric ion complex formation, from which we fin ln K c. Figure 3. A stability map showing fixe sample imensions an experimental conitions for which the ED surface is smooth stable or rough unstable for 4.0, W 1.0, 0.0, ˆ 0.001, an unsupporte electrolyte with polar aitives. The results are shown in Fig. 2, where we plot (k)/ for j 0.5, KC I 0.25, 0.0, 1 0.5, ˆ 0.001, an 1.0, an several values of the parameter which is relate to the ipole moment of the aitives. It can be seen that increasing is beneficial for leveling. However, the main aitional stabilizing effect of ipolar aitives comes from a larger effective equilibrium constant K, as 10 5 for L 10 4 m, appropriate for typical electroplating conitions. Aitive-assiste growth of linearly stable surfaces. The results presente can be use to choose eposition conitions corresponing to the maximum film growth rate for which the surface remains smooth over the requisite length scale, W. To this en, we require that all perturbation moes with wavelengths smaller than W ecay; i.e., (k) 0 for k 2/W. Equation 18 implies j eff (1 1 )/(1 j eff ) ˆ 1 (2/W) 2 for the case of a fully supporte electrolyte. In orer to guie experiment, we express the stability conitions in terms of the experimentally accessible parameters without changing aitives, metal cation flux j which is proportional to the growth rate, an aitive concentration C I. In Fig. 3 we construct a stability map for realistic growth conitions. Such a map conveniently isplays the full range of stable an unstable growth conitions. While introucing aitives can ramatically increase the eposition currents that can be use, the effect saturates at large aitive concentrations. There exists a critical current j* 1 (1 1 )/ above which the aitives are incapable of leveling the surface. This can be unerstoo by incorporating the blocking effect into the secon term of an effective overpotential eff RT/Fln(1 I ), obtaine from the moifie B-V equation. In orer to prouce large eposition fluxes, the metal cation concentration graient ahea of the surface is large an the metal cation concentration is low. To maintain these high fluxes the overpotential must be large, such that RT/Fln(1 I ). This implies that the aitives represente by I ) are incapable of stabilizing the smooth surface. Finally, let us relate the stability map to experiments. The surface is linearly stable at sufficiently small currents j even in the absence of aitives. However, for sufficiently large j, instabilities with wavenumbers greater than 2/W appear, implying roughening. Increasing C I at fixe j j* leas first to slower roughening an then to a transition from a rough to a smooth surface. This behavior

8 Journal of The Electrochemical Society, C708-C C715 is consistent with experimental observations of Schilari et al. 7 which show that the ominant wavelength of the surface instability increases with increasing aitive concentration up to a critical concentration beyon which the surface of the growing film remains level. Similarly, we preict that increasing the eposition flux j for a fixe aitive bulk concentration C I leas to a transition from planar to rough morphology. Conclusions We have introuce a chemically motivate continuum moel for the morphological stability of surfaces uring ED an the effects of aitives on that stability. The ED moel explicitly accounts for the electric fiel in the electrolyte, the metal cations an anions, the aitives, an the spectator ions from the supporting electrolyte. The moel for aitives employe here accounts for three main chemical/physical phenomena: i the aitives form complexes with the metal cations in the electrolyte, (ii) the aitive molecules have a finite ipole moment an hence interact with inhomogeneities in the electric fiels in the electrolyte, an (iii) the aitives occupy sites on the surface, thereby blocking eposition of metal cations. The stability of the ED surface in the presence of aitives was analyze using perturbation theory. In the fully supporte electrolyte case, we foun that aitives enhance leveling by making the surface stable on longer length scales an ecreasing the rate at which roughness grows; our quantitative preictions for the stabilization length are in very goo agreement with existing experimental observations. The aitives stabilize the surface by ecreasing the effective riving force for the instability. The effective riving force epens on the bulk concentration of the aitives, the tenency of the aitives to segregate onto the surface, the strength of complex formation, an the local aitive flux; increasing the magnitue of any of these terms increases the surface stability. The tenency of the aitives to segregate onto the surface is enhance in the case of polar aitives, which are attracte to the surface ue to the presence of strong electric fiel graients. This occurs for any bath composition. Aitionally, the ipole moment of the aitive may help further stabilize the surface in the case of an unsupporte or partially supporte electrolyte, as the polar aitives also interact with the electric fiels associate with surface protrusions. Uner typical eposition conitions, however, we expect this effect to be small. Leveling is promote by increasing the concentration of the aitives, the tenency for the aitives to thermoynamically segregate to the surface, the tenency for aitive-metal cation complexing, incorporation of the aitives into the growing electroeposit, an the tenency of the aitives on the surface to block metal cation attachment. At a minimum, aitives must be able to segregate to the surface in equilibrium an be able to block surface sites once they get there. In aition, at least one of the following must also be present: aitive-metal cation complexing, aitive incorporation in the growing eposit or, equivalently, aitive reuction at the surface, or a nonzero ipole moment if the electrolyte is not fully supporte. These ifferent effects act in concert in a well-esigne aitive. Much of these results can be summarize in stability maps that show uner which experimental conitions a surface is stable or unstable for a given feature size below which the surface must be flat. We have constructe such a map in metal cation flux j i.e., growth rate an aitive concentration C I space. This map clearly shows the existence of a sharp bounary between the stable an unstable growth regimes as well as ientifies conitions uner which increasing the aitive concentration is unable to suppress this natural instability. The experiments of Schilari et al. 7 show that incresing the aitive concentration leas to first ecrease in roughness rather than an increase in roughness above some critical concentration. Those authors attribute this to a phase transition within an aitive layer. Acknowlegments We thank Dr. Corbett Battaile an Dr. John Hamilton for fruitful iscussions, an Sania National Laboratories for financial support. Princeton University assiste in meeting the publication costs of this article. List of Symbols C i concentration of species i C 0 metal cation bulk concentration C 1 spectator cation bulk concentration, normalize by C 0 C I aitive bulk concentration, normalize by C 0 D i iffusivity of species i F Faraay constant h position of the surface in the lab frame j imensionless steay-state metal cation flux j i flux of the species i j 0 exchange flux ensity k initial perturbation wavenumber, in units of 1/L K equilibrium aitive surface asorption constant k B Boltzmann constant K c equilibrium constant for aitive-metal cation complex formation k max maximally unstable moe k 0 neutral moe L thickness of the mass-transfer layer s stabilization length nˆ surface normal, pointing into the solution p strength of the ipole moment per aitive q i e charge of the species i R gas constant T temperature V eq equilibrium potential of the metal-solution interface V ext electric potential of the metal-solution interface uring growth v n normal velocity of the growth front z istance from the surface Greek imensionless equilibrium potential of the metal-solution interface surface tension of the metal-solution interface ˆ imensionless surface tension of the metal-solution interface imensionless aitive ipole moment C i perturbation of the species i j perturbation of the local metal cation flux ˆ initial surface perturbation amplitue imensionless thickness of the G-C bounary layer ˆ permittivity of the solution overpotential I aitive surface coverage local curvature of the metal-solution interface GC thickness of the G-C bounary layer imensionless strength of complex formation electric potential * potential of the metal-solution interface, normalize by the equilibrium potential parameter escribing the extent of aitive coeposition parameter which escribes the strength of complexing atomic volume of the metal in the eposit imensionless atomic volume of the metal in the eposit (k) growth rate of perturbation with wavenumber k References 1. J. W. Dini, Electroeposition: The Materials Science of Coatings an Substrates, Noyes Publications, Park Rige, NJ J. Oniciu an L. Muresan, J. Appl. Electrochem., 21, ; T. C. Franklin, Surf. Coat. Technol., 30, F. Argoul, A. Arneo, G. Grasseau, an H. L. Swinney, Phys. Rev. Lett., 61, G. L. M. K. S. Kahana, X.-Q. Zou, R. Farrell, an P.-O. Wong, Phys. Rev. Lett., 68, D. Barkey, F. Oberholtzer, an Q. Wu, Phys. Rev. Lett., 75, J. M. Pastor an M. A. Rubio, Phys. Rev. Lett., 76, P. L. Schilari, O. Azzoni, an R. C. Salvarezza, Phys. Rev. B, 62, D. Kanel an E. Kaxiras, Soli State Phys., 54, M. Haataja, D. J. Srolovitz, an A. B. Bocarsly, J. Electrochem. Soc., 150, C M. Haataja an D. J. Srolovitz, Phys. Rev. Lett., 89, C. Maore, M. Matlosz, an D. Lanolt, J. Electrochem. Soc., 143, M. Georgiaou, D. Veyret, R. L. Sani, an R. C. Alkire, J. Electrochem. Soc., 148, C Y. Cao, P. Taephaisitphongse, R. Chalupa, an A. C. West, J. Electrochem. Soc., 148, C

9 C716 Journal of The Electrochemical Society, C708-C T. J. Pricer, M. J. Kushner, an R. C. Alkire, J. Electrochem. Soc., 149, C D. Josell, D. Wheeler, W. H. Huber, an T. P. Moffat, Phys. Rev. Lett., 87, D. Josell, B. Baker, C. Witt, D. Wheeler, an T. P. Moffat, J. Electrochem. Soc., 149, C J. D. Jackson, Classical Electroynamics, John Wiley & Sons, New York M. Paunovic an M. Schlesinger, Funamentals of Electrochemical Deposition, John Wiley & Sons, New York M. D. Pritzker an T. Z. Fahiy, Electrochim. Acta, 37, R. Aogaki, K. Kitazawa, Y. Kose, an K. Fueki, Electrochim. Acta, 25, D. P. Barkey, R. H. Muller, an C. W. Tobias, J. Electrochem. Soc., 136, C.-P. Chen an J. Jorne, J. Electrochem. Soc., 138, L.-G. Sunstrom an F. H. Bark, Electrochim. Acta, 40, J. Elezgaray, C. Léger, an F. Argoul, J. Electrochem. Soc., 145, R. Cuerno an M. Castro, Phys. Rev. Lett., 87, C. H. P. Lupis, Chemical Thermoynamics of Materials, North-Hollan, New York C. J. Doona an D. M. Stanbury, Inorg. Chem., 35,

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

water adding dye partial mixing homogenization time

water adding dye partial mixing homogenization time iffusion iffusion is a process of mass transport that involves the movement of one atomic species into another. It occurs by ranom atomic jumps from one position to another an takes place in the gaseous,

More information

Critical Size and Particle Growth

Critical Size and Particle Growth Critical Size an article Growth rof. Sotiris E. ratsinis article Technology Laboratory Department of Mechanical an rocess Engineering, ETH Zürich, Switzerlan www.ptl.ethz.ch 1 Nucleation-Conensation A

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

The Coupling between the Hydration and Double Layer Interactions

The Coupling between the Hydration and Double Layer Interactions 7584 Langmuir 00, 18, 7584-7593 The Coupling between the Hyration an ouble Layer Interactions Eli Ruckenstein* an Marian Manciu epartment of Chemical Engineering, State University of New York at Buffalo,

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

Unified kinetic model of dopant segregation during vapor-phase growth

Unified kinetic model of dopant segregation during vapor-phase growth PHYSICAL REVIEW B 72, 195419 2005 Unifie kinetic moel of opant segregation uring vapor-phase growth Craig B. Arnol 1, * an Michael J. Aziz 2 1 Department of Mechanical an Aerospace Engineering an Princeton

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016 Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]

More information

Modeling and analysis of hydrogen permeation in mixed proton electronic conductors

Modeling and analysis of hydrogen permeation in mixed proton electronic conductors Chemical Engineering Science 58 (2003 1977 1988 www.elsevier.com/locate/ces Moeling an analysis of hyrogen permeation in mixe proton electronic conuctors Lin Li a;b;1, Enrique Iglesia a;b; a Department

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

A Model of Electron-Positron Pair Formation

A Model of Electron-Positron Pair Formation Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating droplets tracking by holographic high speed video in turbulent flow Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

Efficient Macro-Micro Scale Coupled Modeling of Batteries

Efficient Macro-Micro Scale Coupled Modeling of Batteries A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue

PCCP PAPER. 1 Introduction. A. Nenning,* A. K. Opitz, T. M. Huber and J. Fleig. View Article Online View Journal View Issue PAPER View Article Online View Journal View Issue Cite this: Phys. Chem. Chem. Phys., 2014, 16, 22321 Receive 4th June 2014, Accepte 3r September 2014 DOI: 10.1039/c4cp02467b www.rsc.org/pccp 1 Introuction

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

C. Montella * Abstract. Keywords: Potential step; chronoamperometry; chronocoulometry; absorption; intercalation; diffusion coefficient.

C. Montella * Abstract. Keywords: Potential step; chronoamperometry; chronocoulometry; absorption; intercalation; diffusion coefficient. J. Electroanal. Chem., sous presse Discussion of the potential step metho for the etermination of the iffusion coefficients of guest species in host materials. I. Influence of charge transfer kinetics

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field

Quantum optics of a Bose-Einstein condensate coupled to a quantized light field PHYSICAL REVIEW A VOLUME 60, NUMBER 2 AUGUST 1999 Quantum optics of a Bose-Einstein conensate couple to a quantize light fiel M. G. Moore, O. Zobay, an P. Meystre Optical Sciences Center an Department

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

Applications of First Order Equations

Applications of First Order Equations Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation

More information

Interaction force in a vertical dust chain inside a glass box

Interaction force in a vertical dust chain inside a glass box Interaction force in a vertical ust chain insie a glass box Jie Kong, Ke Qiao, Lorin S. Matthews an Truell W. Hye Center for Astrophysics, Space Physics, an Engineering Research (CASPER) Baylor University

More information

arxiv:physics/ v4 [physics.class-ph] 9 Jul 1999

arxiv:physics/ v4 [physics.class-ph] 9 Jul 1999 AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly

More information

Project 3 Convection

Project 3 Convection 12.307 Project 3 Convection Loovica Illari an John Marshall March 2010 1 Backgroun The Earth is bathe in raiation from the Sun whose intensity peaks in the visible. In orer to maintain energy balance the

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION Long-istance coherent coupling in a quantum ot array Floris R. Braakman 1, Pierre Barthelemy 1, Christian Reichl, Werner Wegscheier, Lieven M.K. Vanersypen 1 1 Kavli Institute of Nanoscience, TU Delft,

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

arxiv:cond-mat/ v1 12 Oct 2005

arxiv:cond-mat/ v1 12 Oct 2005 Creation an estruction of a spin gap in weakly couple quarter-fille laers B. Eegger,, H.G. Evertz, an R.M. Noack Institut für Theoretische Physik, Technische Universität Graz, A-8 Graz, Austria Institut

More information

(3-3) = (Gauss s law) (3-6)

(3-3) = (Gauss s law) (3-6) tatic Electric Fiels Electrostatics is the stuy of the effects of electric charges at rest, an the static electric fiels, which are cause by stationary electric charges. In the euctive approach, few funamental

More information

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0

Extinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0 Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

LeChatelier Dynamics

LeChatelier Dynamics LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a

More information

Lecture contents. Metal-semiconductor contact

Lecture contents. Metal-semiconductor contact 1 Lecture contents Metal-semiconuctor contact Electrostatics: Full epletion approimation Electrostatics: Eact electrostatic solution Current Methos for barrier measurement Junctions: general approaches,

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty Plasma Void Dynamics in Unmoving and Moving Flows 7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological

More information

Light-Soaking Effects on the Open-Circuit Voltage of a-si:h Solar Cells

Light-Soaking Effects on the Open-Circuit Voltage of a-si:h Solar Cells Syracuse University SURFACE Physics College of Arts an Sciences 2005 Light-Soaking Effects on the Open-Circuit Voltage of a-si:h Solar Cells Jianjun Liang Syracuse University Eric A. Schiff Syracuse University

More information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012 arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity

More information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

ABCD42BEF F2 F8 5 4D658 CC89

ABCD42BEF F2 F8 5 4D658 CC89 ABCD BEF F F D CC Vetri Velan GSI, Physics 7B Miterm 2: Problem Solution. Outsie sphere, E looks like a point charge. E = The total charge on the sphere is Q sphere = ρ 4 3 πr3 Thus, outsie the sphere,

More information

(NaCl) x (KCl) y-x (KBr) 1-y single crystals: study of the ac conductivity activation energy.

(NaCl) x (KCl) y-x (KBr) 1-y single crystals: study of the ac conductivity activation energy. (NaCl) x (KCl) y-x (KBr) 1-y single crystals: stuy of the ac conuctivity ivation energy. Vassiliki Katsika-Tsigourakou * Department of Soli State Physics, Faculty of Physics, University of Athens, Panepistimiopolis,

More information

v r 1 E β ; v r v r 2 , t t 2 , t t 1 , t 1 1 v 2 v (3) 2 ; v χ αβγδ r 3 dt 3 , t t 3 ) βγδ [ R 3 ] exp +i ω 3 [ ] τ 1 exp i k v [ ] χ αβγ , τ 1 dτ 3

v r 1 E β ; v r v r 2 , t t 2 , t t 1 , t 1 1 v 2 v (3) 2 ; v χ αβγδ r 3 dt 3 , t t 3 ) βγδ [ R 3 ] exp +i ω 3 [ ] τ 1 exp i k v [ ] χ αβγ , τ 1 dτ 3 ON CLASSICAL ELECTROMAGNETIC FIELDS PAGE 58 VII. NONLINEAR OPTICS -- CLASSICAL PICTURE: AN EXTENDED PHENOMENOLOGICAL MODEL OF POLARIZATION : As an introuction to the subject of nonlinear optical phenomena,

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an

More information

Numerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number

Numerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number Journal of Applie Flui Mechanics, Vol. 9, No. 2, pp. 643-651, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Investigation of Non-Stationary Parameters on Effective

More information

Collective optical effect in complement left-handed material

Collective optical effect in complement left-handed material Collective optical effect in complement left-hane material S.-C. Wu, C.-F. Chen, W. C. Chao, W.-J. Huang an H. L. Chen, National Nano Device Laboratories, 1001-1 Ta-Hsueh Roa, Hsinchu, Taiwan R.O.C A.-C.

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses.

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses. Physics 195a Course Notes The K 0 : An Interesting Example of a Two-State System 021029 F. Porter 1 Introuction An example of a two-state system is consiere. involve: 1. Treatment of a ecaying particle.

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Thermodynamic basis for a variational model for crystal growth

Thermodynamic basis for a variational model for crystal growth PHYSICL REIEW E OLUME 60, NUMBER 1 JULY 1999 Thermoynamic basis for a variational moel for crystal growth Bayar K. Johnson* an Robert F. Sekerka Carnegie Mellon University, Pittsburgh, Pennsylvania 1513

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Thermal runaway during blocking

Thermal runaway during blocking Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application

More information

Magnetic helicity evolution in a periodic domain with imposed field

Magnetic helicity evolution in a periodic domain with imposed field PHYSICAL REVIEW E 69, 056407 (2004) Magnetic helicity evolution in a perioic omain with impose fiel Axel Branenburg* Norita, Blegamsvej 17, DK-2100 Copenhagen Ø, Denmark William H. Matthaeus University

More information

Formulation of statistical mechanics for chaotic systems

Formulation of statistical mechanics for chaotic systems PRAMANA c Inian Acaemy of Sciences Vol. 72, No. 2 journal of February 29 physics pp. 315 323 Formulation of statistical mechanics for chaotic systems VISHNU M BANNUR 1, an RAMESH BABU THAYYULLATHIL 2 1

More information

Chapter 11: Feedback and PID Control Theory

Chapter 11: Feedback and PID Control Theory Chapter 11: Feeback an D Control Theory Chapter 11: Feeback an D Control Theory. ntrouction Feeback is a mechanism for regulating a physical system so that it maintains a certain state. Feeback works by

More information

Magnetic field generated by current filaments

Magnetic field generated by current filaments Journal of Phsics: Conference Series OPEN ACCESS Magnetic fiel generate b current filaments To cite this article: Y Kimura 2014 J. Phs.: Conf. Ser. 544 012004 View the article online for upates an enhancements.

More information

Study of the composition dependence of the ionic conductivity of LiH x D 1-x alloys

Study of the composition dependence of the ionic conductivity of LiH x D 1-x alloys Stuy of the composition epenence of the ionic conuctivity of LiH x D 1-x alloys Vassiliki Katsika-Tsigourakou Section of Soli State Physics, Department of Physics, National an Kapoistrian University of

More information

Advanced Partial Differential Equations with Applications

Advanced Partial Differential Equations with Applications MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.

More information