Thickness and refractive index measurements using multiple beam interference fringes (FECO)

Size: px
Start display at page:

Download "Thickness and refractive index measurements using multiple beam interference fringes (FECO)"

Transcription

1 Journal of Colloid and Interface Science Note Thickness and refractive index measurements using multile beam interference fringes FECO Rafael Tadmor, 1 Nianhuan Chen, and Jacob N. Israelachvili Materials Research Laboratory and Deartment of Chemical Engineering, University of California, Santa Barbara, CA 93106, USA Received 28 January 2002; acceted 11 Aril 2003 Abstract We reort on the use of otical interferometry emloying fringes of equal chromatic order FECO in a surface force aaratus SFA to determine film thicknesses and refractive indices of confined media for a wide range of searations. In articular, we show how to calculate the surface searation film thickness based on two fringes whose contact osition was not measured. We discuss the measurement accuracy, and though the theoretical accuracy is 1 Å for all searations, we show that in ractice, for large searations, it is very hard to get to this accuracy Elsevier Inc. All rights reserved. Keywords: Large searations; Thick films; Three-layer interferometer; Disersion; Disersive hase change In a SFA surface force aaratus [1] exeriment, there is usually a three- or five-layer interferometer [2]. For simlicity we assume here a three-layer interferometer, although the analysis alies to any number of layers. When the transarent substrate layers tyically mica sheets are brought into contact, the ositions λ 0 of the th-order fringes is usually measured in the visible range of wavelengths. After the substrate surfaces are searated by a distance D, a three-layer interferometer is formed. The fringes shift to longer wavelengths and their new ositions λ D are given see Ref. [3] by tan 2πµ med D/λ D 2 µ sin4πµy/λ D = 1 µ µ 2 cos4πµy/λ D, which, assuming no disersion we discuss this oint later in the text, reduces to the more useful form tan 2πµ med D/λ D 2 µ sin π 1 λ0 /λd = 1 + µ 2 cos 1 λ 0 /λ0 1 π 1 λ0 /λd 1 λ 0 /λ0 1 ± µ 2 1, 1a * Corresonding author. address: jacob@engineering.ucsb.edu J.N. Israelachvili. 1 Present address: Det. of Chemical Engineering, Lamar University, Beaumont, TX where Y is the otical thickness of each substrate we discuss later the difference between the otical and hysical thicknesses and µ = µ/µ med,whereµµ γ or µ β and µ med are, resectively, the refractive indices of the mica and the intervening medium at λ D.InEq.1a,+ and refer to odd and even order fringes, resectively. Other ways of calculating thicknesses are available see, for examle, Ref. [4]; however, the simlicity of Eq. 1a made it the most common method. Indeed, due to the increased use of Eq. 1a in SFA measurements, eole have addressed some secial cases which were not addressed by the original ublication [3], such as reflecting metallic media [5], nonsymmetric three-layer interferometry [6], and most recently the slitting of the fringes as a result of a liquid crystalline medium [7]. Currently Eq. 1a is limited due to the fact that the fringes at contact need to be of the same order as the fringes of measurements. Aart from imosing inconvenience, this also limits the range of measurable searations, since at far searations, the contact fringes are no longer in the visible range. Different aroaches to this roblem were addressed by Horn and Smith [6], by Heuberger et al. [4], and by Farrell et al. [8]. In this study, however, we resent a very simle account for this roblem. An even simler account than the one resented here is given in the thesis of H. Christenson [9]. This solution, however, is extremely sensitive to the osition of the fringes, sub-ångstrom errors in the wavelength measurement can roduce a micrometer-scale /$ see front matter 2003 Elsevier Inc. All rights reserved. doi: /s

2 R. Tadmor et al. / Journal of Colloid and Interface Science error in the searation. This aer describes an easy and accurate way of calculating the searation from any order of fringes at any surface searation. From Eq. 1a it is clear that if the refractive indices of the mica substrate and the medium are known then the searation between the surfaces D can be obtained from λ D, λ0, and λ 0 1. In the case where the medium refractive index µ is not known [10,11], there is a need to solve two equations with two unknowns, D and µ med. The two equations are Eq. 1a for an odd fringe and Eq. 1a for an even fringe. In order to solve these simultaneous equations, we need to know the location of two adjacent fringes, λ D and λd 1,at finite D and the contact ositions of three adjacent fringes, λ 0, λ0 1,andλ0 2,atD = 0. The limitation of the above conventional method is that λ D, λd 1, λ0, λ0 1,andλ0 2 need to be in the visible for them to be measured. This limits the range of measurable searations D obtained using this aroach to D<1µm for tyical substrate thicknesses thicker substrates would increase the range but decrease the accuracy of the measurements. In this aer we describe a rocedure that allows accurate and unambiguous measurements of distance and refractive index to be made using the contact osition of only two adjacent fringes, λ 0 and λ0 1, and more imortantly using two adjacent fringes of any order, λ D m and λd m 1, and not necessarily λ D and λd General rocedure for measuring D and µ med of films of arbitrary thickness The contact ositions, λ 0, corresond to a single-layer interferometer. A one-layer interferometer made of a substrate with a refractive index µ and hysical thickness Y, obeys the relation 2µY = λ, = 1, 2, 3,...,, 2 where λ is the wavelength that corresonds to constructive interference and is a natural number the fringe order. In ractice there is never a one-layer interferometer, since the reflecting layer results in a hase change of the light at the substrate reflector interface say a mica silver interface. Such a hase change for the mica silver interface reflection has been studied [12], and it has been shown that it can be viewed as an aarent small change in the otical thickness of the substrate to Y instead of Y [13]: 2µY = λ, = 1, 2, 3,...,. 2a It is intuitive that the solution to the roosed roblem may be rovided using this equation. Secifically, using Eq. 2a to substitute the thickness of the layer Y for the contact osition of the fringe in Eq. 1a, we should be able to relate one contact fringe osition to any other contact fringe osition. Hence, the roblem becomes one of how to exress a singlelayer interferometer in terms of a relation between different fringe orders rather than the resentation of Eq. 2a. Using Eq. 2a and assuming constant Y [13], we write for three arbitrary fringes of order, 1, and m 2µ Y = λ 0, 3a 3b 2µ 1 Y = 1λ 0 1, 2µ m Y = mλ 0 m, 3c where we write µ to emhasize the disersive nature of the substrate, i.e., that each refractive index µ corresonds to a wavelength λ 0. Equations 3a and 3b have just two unknowns, Y and, and solving for them we get 1 = 1 µ 1λ 0 µ λ 0. 1 From Eqs. 3a and 3c and using Eq. 4 we obtain λ 0 m = µ m λ 0. 5 µ 1 m 1 µ 1λ 0 µ λ 0 1 Equation 5 relates the osition of the mth fringe to that of the th fringe as a function of the number of fringes m between them, as well as five other arameters that need to be known, including µ m. Note that m is an integer, which may be ositive, negative, or zero. Now, we may use any fringe in order to use Eq. 1a for the calculation of searations and refractive indices: we simly need to count the number of fringes from m to, ut this number as m into Eq. 5, and calculate a new set of contact ositions λ 0 m, λ 0 m 1,andλ0 m 2. As we show later in the worked Examle 1, the value of µ m can be obtained by substitution of an aroximated λ m given by Eq. 6 into Eq. 12, which we introduce below. Exerimentally, if a measurement is done at a large searation on a fringe λ D m whose contact osition λ 0 m is not visible, one way of finding m after comleting a measurement using the mth fringe is to bring the surfaces to substrate substrate contact while counting the assing fringes until one finds the -fringe. This number is the integer m. Equation 5 requires knowledge of µ m. To obtain this, note that Eq. 5 has a very weak deendence on the refractive index of the substrate. Neglecting the disersion, for examle, results in λ 0 m = λ m1 λ 0 /λ0 1, which is an excellent aroximation for many cases. Equation 6 can also be used for a first guess for the wavelength at which one needs to look for µ m in Eq. 5. In worked Examle 1, we show how this can be done. Similarly, Eq. 5 or its aroximated version, Eq. 6, which is written for any λ 0 m, can also give λ0 2, thus reducing the minimal number of required contact fringes to 2. In fact, there are cases where only two fringes are in the visible range for very thin mica, for which Eq. 5 may be very useful. In worked Examle 2 we show that for a three-layer interferometer of known 4

3 550 R. Tadmor et al. / Journal of Colloid and Interface Science medium refractive index and substrate disersion [14], one can theoretically calculate the ositions of the contacting fringes even without counting the integer m. However, in ractice this can be used only for small m; therefore, for large searations, m usually needs to be counted from the assing fringes as described above. The simlified result of Eq. 6 can also be obtained by utting D = 0 in Eq. 1a. This gives 1 λ 0 /λd 7 1 λ 0 = z, /λ0 1 where z is an integer. In Eq. 7, λ D is not the contact osition λ 0,butany contact osition λ0 m. This is because any such fringe would corresond to tan2πµ med D/λ D = 0in Eq. 1a. Thus, solving for Eq. 7, we again get Eq. 6. It is clear from the above treatment that Eq. 1a does not take into account the disersion of the substrate as was noted already in the original aer [3]. The disersion of mica has been studied [15,16], and below we derive the equation corresonding to Eq. 1 for disersive substrates and media. Using the more fundamental form of Eq. 1, rather than 1a, we note that a one-layer interferometer imlies that 4πµY 8 λ D = π λ 0 /λd, where µ is the refractive index of the mica at λ D as noted in Eq. 1. Substituting Eqs. 4 and 8 into Eq. 1, we obtain the disersive version of 1a as tan 2πµ med D/λ D = 2 µ sin π 1 λ 0 /λd 1 µ 1λ 0 1 µ λ 0 1 [ 1 + µ 2 cos π 1 λ 0 /λd 1 µ 1λ 0 1 µ λ 0 1 and then use Eq. 5 to relace the unknowns λ 0 m and λ0 m 1 with the measured λ 0 and λ0 1 and the calculated λ0 2 in the case where the refractive index of the medium is not known. In certain situations, the integer number m may not be available or measurable. Since we have two different equations 10 for odd and even fringes, we can use the measured λ D m and λd m 1, and simly guess a value for m. Then, using Eq. 5 with λ 0 and λ0 1, we calculate the contact ositions which corresond to fringes m, m 1, and m 2. Aarently, only if the guess is correct will one obtain the same searation D using both forms of Eq. 10 for even and for odd fringes. In ractice, however, this method which anyway works only for a three-layer interferometer can be done only for very small m values, while for large m other factors such as the disersion of the medium µ med and the disersive hase change at the reflector substrate interface should also be known very accurately. This is discussed further after Examle 2. The examles below are based on measured values. 2. Worked examles 2.1. Examle 1 Estimating the value of λ 0 2 from λ0 and λ0 1. In an exeriment using brownish mica as substrate, the contact wavelengths of three adjacent fringes of unknown order, 1, and 2 are measured Fig. 1 and found to be at λ 0 = ± 0.11 Å, λ0 1 = ± 0.16 Å, and λ 0 2 = ± 0.21 Å. If we want to estimate the osition of λ 0 2 from λ0 and λ0 1, then utting m = 2 into Eq. 6 we obtain λ 0 2 = / = Å, neglecting disersion. 11 ± µ 2 1] 1. 9 Equation 9 is the analogue of Eq. 1a with a correction for the disersion in the substrate refractive index, and as in Eq. 1a, + and refer to odd and even order fringes, resectively. Finally, we may write the analogue to Eq. 9 for any order of fringe, tan 2πµ med D/λ D m = 2 µ sin π 1 λ 0 m /λd m 1 µ m 1λ 0 1 m µ m λ 0 [ m µ 2 cos π 1 λ 0 m /λd m 1 µ m 1λ 0 1 m µ m λ 0 m 1 ± µ 2 1] 1, 10 Fig. 1. Schematic reresentation of the locations and shaes of the FECO in Examles 1 and 2. The flat arts in λ 0, λ0 1,andλ0 2, result from the elastic distortion of the glue which suorts the mica [3]. For such a flat contact of mica in air, odd...,, 2,... fringes always have a similar shae, which is different from that of the even..., 1, 3,... fringes. Secifically, the lace at which the fringes sto being flat and start bending has a bigger discontinuity in its derivative for odd fringes than for even fringes. Also note that since m>, λ D m, λd m 1,andλD m 2 are closer together than λ 0, λ0 1,andλ0 2.

4 R. Tadmor et al. / Journal of Colloid and Interface Science If we include disersion in the refractive index, then we need to know the relation for the refractive index of the fringe β or γ that is being measured as a function of λ. For brownish mica, this is [17] µ β = /λ 2 for β fringes, 12 where λ is the wavelength in Å. We first calculate the refractive indices at the three contact wavelengths, and obtain µ β λ 0 = , µ βλ 0 1 = , µ βλ 0 2 = Now, using Eq. 5, we obtain λ 0 2 = [ ] = Å, 13 which is in excellent agreement with the measured value of ± 0.21 Å, indeed, within the exerimental error. Actually, one could use the measured value of λ 0 2 to determine one of the constants in Eq. 12 or, if we also measured λ 0 +1 or λ0 3, to get both constants. Note that it is not necessary to know the value of, nor does it matter if is odd or even the equations are identical. Only for the searation measurements of the three-layer interferometer is it imortant to know whether is odd or even it haened to be odd for this secific examle. We also note that this method can be used to calculate the contact wavelength of any fringe order, not just λ Examle 2 Calculating distances D from λ 0 and λ0 1, and the measured ositions λ D m and λd m 1 of any two adjacent fringes of unknown order m and m 1. A measured contact osition of a fringe of unknown order is at λ 0 = Å and that of order 1isat λ 0 1 = Å. The substrate is brownish mica, whose refractive index is given by Eq. 12. The surfaces are well searated and we want to calculate the distance D between the two surfaces. We know that the medium between the two mica substrates has a refractive index of [14], and we erform simultaneous measurements of λ D m and λd m 1 and obtain λ D m = Å and λd m 1 = Å see Fig. 1 for the qualitative relative ositions of λ 0, λ0 1, λ D m, λd m 1 in this examle. For a quick estimate of D, one can use Eq. 6 to calculate the contact ositions of various fringes + 1,+ 2,... see Examle 1 and Eq. 1a to calculate the distances of our three-layer interferometer until + m is found. A more accurate aroach is to use Eqs. 5 and 6 for the calculation of the ositions + 1, + 2,..., as described in Examle 1, and Eq. 9 for the distances. In this examle we show a case where Eq. 5 is used with Eq. 1a [14]. Obviously m ; hence our first guess is m =. Using this guess, we ut the values for λ 0, λ0 1,andλD m which we guess to be λd in the odd version of Eq. 1a and obtain D = 1363 Å. We then erform the same calculation using λ D m 1 which we guess to be λ D 1 and the values for λ0 1 and λ0 2 using the even version of Eq. 1a, and obtain D = 1310 Å. Since , our guess was wrong. Our second guess would be m = + 1wherem is now even, and we use the even version of Eq. 1a with the calculated value of λ 0 +1, measured λ0 and λd m which we guess to be λd +1. We obtain D = 3288 Å. The odd version of Eq. 1a with λ 0, λ0 1,andλD m 1 which we guess to be λd,nowresults in D = 3280 Å. Since 3288 = 3280, we conclude that m = + 1. The difference of 8 Å is in art because the refractive index of the medium was assumed to be nondisersive, and in art because the mica silver hase change was also assumed to be nondisersive. We should also note that for m = + 2, we obtain D = 5217 Å using the odd version of Eq. 1a with λ 0 +2, λ0 +1,andλD m guessedtobe λ D +2,andD = 5248 Å using the even version of Eq. 1a with λ 0 +1, λ0,andλd m 1 guessed to be λd +1,andwe get Dλ D m,λ0 +2,λ0 +1 <DλD m 1,λ0 +1,λ0. Indeed, for all guesses above + 1 the second estimate is higher than the first, while for all guesses below + 1 the reverse is true, and this observation is general. In other words, the searation calculated using a higher order contact fringe coule is bigger than that using a lower order contact fringe coule for m guesses which are smaller than the real m, whereas the searation calculated using a higher order contact fringe coule is smaller than that using a lower order contact fringe coule for m guesses which are bigger than the real m. Formulating this would look like this: D λ D m,λ0 +i +i 1,λ0 >D λ D m 1,λ 0 +i 1,λ0 +i 2 for i<m, D λ D m,λ0 +i +i 1,λ0 = D λ D m 1,λ 0 +i 1,λ0 +i 2 for i = m, D λ D m,λ0 +i +i 1,λ0 >D λ D m 1,λ 0 +i 1,λ0 +i 2 for i>m. The above examle uses only m = + 1 i = 1 and would work for most cases only for very small values of m. Thus, in general it is advised to count the number of assing fringes in order to determine m. Theoretically, it should be ossible to include the disersive refractive index of the medium and the disersive hase change at the reflector substrate interface. However, as we discuss below, incororating a disersive hase change, while theoretically desirable, is difficult in ractice. Since the above did not include the disersive hase change, the similarity obtained for the two calculations of D does not suggest that this is the error in the searation, which could be larger. As we shall see, incororation of the hase change gives a formally accurate solution; however, there are good reasons to use only the substrate disersive refractive index and neglect the disersive hase change.

5 552 R. Tadmor et al. / Journal of Colloid and Interface Science Comaring disersion and hase change If we consider a normal SFA exeriment with mica as a substrate and a silver layer of 550 Å as a reflector, then out of the two corrections to the ideal osition of the fringes substrate refractive index disersion and substrate reflector disersive hase change, the substrate disersion is usually the bigger contributor, as demonstrated in Fig. 2. In addition to having a smaller effect, the hase change adds comlexity and requires the use of henomenological relations. The above discussion does not neglect hase change comletely, but rather uses a constant value for the hase change at D = 0 as exlained by Bailey and co-workers [8], since the wavelength ositions measured in contact corresond to the reflector hase changing system that roduced them. Thus, for many systems, this constant value correction suffices. However, there are situations in which one would like to incororate the disersive hase change. 4. Incororating the hase change at the substrate reflector mica silver interface For better accuracy there is a need to incororate the hase change at the substrate reflector interface. We use the idea suggested by Bailey and co-workers [8], and we follow their methodology to incororate the hase change contribution exlicitly. The hase change has the effect of adding or subtracting an aarent thickness y to the real thickness, and hence mica of real thickness Y will, due to the hase change, aear to be of thickness Y = Y + y. Note that y is a function of the wavelength. We start again with the single layer interferometer, and following Bailey and co-workers [8], we can rewrite Eq. 3 as 2µ Y + y = λ 0, 14a 14b 2µ 1 Y + y 1 = 1λ 0 1, 2µ m Y + y m = mλ 0 m. 14c Solving Eq. 14 for, weget 1 = λ0 1 λ0 µ 1 µ λ 0 1 µ 1y y 1, 15 where y i is the value of y for the wavelength λ 0 i that is, yλ 0 i and it can easily be seen that for a constant nondisersive y, Eq. 15 reduces to Eq. 4. Generally yλ i and articularly yλ 0 i is a function of the disersive hase change φλ i, λ i yλ i = φλ i 16 4πµλ i, where φλ i was shown by Bailey and co-workers [8] to be given by the exression kλ i φλ i = arctan nλ i µλ i kλ i arctan, 17 nλ i + µλ i where kλ i and nλ i are the imaginary and real arts of the comlex refractive index of the reflector usually silver in our case. Tabulated values for kλ i and nλ i can be found in ordinary handbooks see, for examle, Ref. [18]. In Fig. 2, we used the heuristic relation from Ref. [19], which was obtained from measurements with a silver layer Å thick deosited by evaoration: Fig. 2. The difference [Å] between the calculated and measured exerimental contact wavelengths for three calculations: square use Eq. 6 neglecting both disersion of mica and of mica silver interface; triangle use Eq. 5 neglecting only disersion of mica silver interface; circle use Eq. 20 considering disersion both of mica and of mica silver interface. The filled squares on the vertical line across the 0 of the ordinate corresond to the measured exerimental value. All the calculations were done using the two rightmost measured exerimental data oints. kλ i = λ i λ 2 i, nλ i = In Ref. [19], silver films of thickness Å were measured, and no deendence on film thickness was observed for films thicker than 250 Å. Note that Bailey and co-workers also do not consider the hase change at the mica silver interface to be a strong function of the silver thickness [8]; rather, the silver roughness, which deends on the way of silver deosition, the temerature and ressure, and the duration of the deosition, is the dominant factor in the hase change. Now, the value of y y 1 can easily be calculated from y y 1 = y λ 0 y λ 0 1 = φ λ 0 λ πµ 0 φ λ 0 λ πµ 0, 1 where φλ i is given by Eq. 17. This in turn can be substituted into Eq. 15 to calculate 1/.

6 R. Tadmor et al. / Journal of Colloid and Interface Science We can now use Eq. 14 to calculate the osition of any contact fringe λ 0 m based on two measured fringes λ0 and λ 0 1, λ 0 m = µ m λ 0 µ 1 m λ 0 1 µ 1 y y 1 λ 0 1 λ0 µ 1/µ 2µ m y y m λ 0, 1 µ 1 y y 1 λ 0 1 λ0 µ m 1/µ 20 and one can see that for a nondisersive y, Eq. 20 reduces to Eq. 5. Figure 2 shows the calculated ositions of λ 0 m using Eqs. 20, 5, and 6. Equation 6 is an aroximate solution which neglects both mica disersion and mica silver hase disersion; Eq. 5 accounts for the substrate mica disersion but neglects the hase change at the substrate reflector mica silver interface; and Eq. 20 accounts for both substrate and substrate reflector hase disersions. We can see that as the equation used is more exact the calculated wavelength osition becomes indeed increasingly close to the measured one. Aarently, however, Eq. 20 fails to calculate the exact osition of wavelength and does not rovide an absolutely accurate solution. This is not surrising for two reasons. One, as noted earlier, is that the relations in Eq. 18 are deendent on the exact exerimental conditions that formed the silver layer, which can be different from exeriment to exeriment. Another reason is that Eq. 20 has two oosite large corrective terms, which eventually result in a small correction. Hence a small inaccuracy in the measurement of the wavelength or in Eq. 18 could result in a rather big error if Eq. 20 were used. Unlike Eq. 5, where the corrective term merely multilies an existing term by a small number. Thus, relying on literature values for the reflector hase change is roblematic since they may use different conditions than the exerimental ones, which cannot necessarily be accounted for e.g., silver roughness at the interface. We therefore suggest measuring the values of kλ i and nλ i using the known m values in Eq. 20 for every exeriment in which it may be necessary and recalculating the constants using Eq. 18. The only roblem with this aroach is that unlike the disersive refractive index, which uses the rather well-established Cauchy relation, the equations for the disersive hase change are urely heuristic emirical relations, and hence the extraolation and erhas even the interolation of the exerimental values is not well-founded and could be erroneous. To comlete the icture, we write the equation for the three-layer interferometer that accounts for the hase change. Here again we follow Bailey and co-workers [8], who noted that in Eq. 1 the term 4πµY/λ D can be relaced by [ 4πµ D Y/λ D = πλ0 21 λ D µ D 4µ D y µ 0 λ 0 y D ], where µ D is µλ D, µ 0 is µλ 0, yd = yλd, and y = yλ 0. Substituting from Eq. 15 into Eq. 21 and substituting that into Eq. 1 would give the correct exression for the three-layer interferometer, which accounts for both the disersion and the hase change. Acknowledgments We thank Qi Lin, Peter Ercius, and Linh Tran for useful discussions and comutational assistance. NASA Grant NAG is gratefully acknowledged. References [1] J.N. Israelachvili, Intermolecular and Surface Forces, Academic Press, London, [2] S. Tolansky, Interferometry, Longman, London, [3] J.N. Israelachvili, J. Colloid Interface Sci [4] M. Heuberger, G. Luengo, J. Israelachvili, Langmuir [5] M.T. Clarkson, J. Phys. D Al. Phys [6] R.G. Horn, D.T. Smith, Al. Ot [7] V. Kitaev, E. Kumacheva, J. Phys. Chem. B [8] B. Farrell, A.I. Bailey, D. Chaman, Al. Ot [9] H.K. Christenson, Ph.D. thesis, Australian National University, [10] J. Janik, R. Tadmor, J. Klein, Langmuir [11] J. Janik, R. Tadmor, J. Klein, Langmuir [12] E. Eisner, Research [13] The otical thickness Y which takes into account the hase change of light at the reflector substrate interface, as oosed to the hysical thickness Y, is a function of the wavelength and the thickness of the reflector silver layer. Throughout the first art of the aer the aroximation of a constant Y is made. Although later in the aer the disersive Y is accounted for, note that in ractice, the use of a constant Y is more useful for most cases. [14] In rincile, one should use the disersive form of the refractive index for the medium as well; however, from Eq. 10 one sees that the refractive index of the medium is only relevant for the measured fringes which, by definition, are all in the visible range. Thus, the error in the disersion of the medium is limited to the narrow visible range. On the other hand, the refractive index of the substrate should be known at the measured visible fringe as well as at the contact ositions which could be well outside the visible range. Thus, the disersion of the substrate does need to be known. [15] A.I. Bailey, S.M. Kay, Br. J. Al. Phys [16] S.Y. El-Zaiat, Ot. Laser Technol [17] J.N. Israelachvili, G.E. Adams, J. Chem. Soc. Faraday Trans [18] D.R. Lide Ed., CRC Handbook of Chemistry and Physics, CRC Press, London, [19] P.B. Johnson, R.W. Christy, Phys. Rev. B

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system

More information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: v1 [physics.data-an] 26 Oct 2012 Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch

More information

Effective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes

Effective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes hys. stat. sol. b 36, 65-633 003 Effective conductivity in a lattice model for binary disordered media with comlex distributions of grain sizes R. PIASECKI Institute of Chemistry, University of Oole, Oleska

More information

An Improved Calibration Method for a Chopped Pyrgeometer

An Improved Calibration Method for a Chopped Pyrgeometer 96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum

More information

State Estimation with ARMarkov Models

State Estimation with ARMarkov Models Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,

More information

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.

More information

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands

More information

MATH 2710: NOTES FOR ANALYSIS

MATH 2710: NOTES FOR ANALYSIS MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite

More information

Statics and dynamics: some elementary concepts

Statics and dynamics: some elementary concepts 1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and

More information

The Quark-Parton Model

The Quark-Parton Model The Quark-Parton Model Before uarks and gluons were generally acceted Feynman roosed that the roton was made u of oint-like constituents artons Both Bjorken Scaling and the Callan-Gross relationshi can

More information

Linear diophantine equations for discrete tomography

Linear diophantine equations for discrete tomography Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,

More information

A Model for Randomly Correlated Deposition

A Model for Randomly Correlated Deposition A Model for Randomly Correlated Deosition B. Karadjov and A. Proykova Faculty of Physics, University of Sofia, 5 J. Bourchier Blvd. Sofia-116, Bulgaria ana@hys.uni-sofia.bg Abstract: A simle, discrete,

More information

Optical Fibres - Dispersion Part 1

Optical Fibres - Dispersion Part 1 ECE 455 Lecture 05 1 Otical Fibres - Disersion Part 1 Stavros Iezekiel Deartment of Electrical and Comuter Engineering University of Cyrus HMY 445 Lecture 05 Fall Semester 016 ECE 455 Lecture 05 Otical

More information

Convergence performance of the coupled-wave and the differential methods for thin gratings

Convergence performance of the coupled-wave and the differential methods for thin gratings Convergence erformance of the couled-wave and the differential methods for thin gratings Philie Lalanne To cite this version: Philie Lalanne. Convergence erformance of the couled-wave and the differential

More information

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College

More information

PRELIMINARY SELECTION OF BASIC PARAMETERS OF DIFFERENT TORSIONAL VIBRATION DAMPERS INTENDED FOR USE IN MEDIUM-SPEED DIESEL ENGINES

PRELIMINARY SELECTION OF BASIC PARAMETERS OF DIFFERENT TORSIONAL VIBRATION DAMPERS INTENDED FOR USE IN MEDIUM-SPEED DIESEL ENGINES Ivan Filiović Dževad Bibić Aleksandar Milašinović Almir Blažević Adnan Pecar ISSN 333-4 PRELIMINARY SELECTION OF BASIC PARAMETERS OF DIFFERENT TORSIONAL VIBRATION DAMPERS INTENDED FOR USE IN MEDIUM-SPEED

More information

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites

More information

The directivity of the forced radiation of sound from panels and openings including the shadow zone

The directivity of the forced radiation of sound from panels and openings including the shadow zone The directivity of the forced radiation of sound from anels and oenings including the shadow zone J. Davy RMIT University, Alied Physics, GPO Box 476V, 3001 Melbourne, Victoria, Australia john.davy@rmit.edu.au

More information

5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:

5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell: 5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum

More information

A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAPE

A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAPE '. ' 220 HILlS TECHNICAL REVIEW VOLUME 20 A M,ETHOD OF MEASURING THE RESISTIVITY AND HALL' COEFFICIENT ON LAMELLAE OF ARBITRARY SHAE 621.317.331:538.632.083 Resistivity and Hall-coefficient measurements

More information

arxiv:cond-mat/ v2 25 Sep 2002

arxiv:cond-mat/ v2 25 Sep 2002 Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,

More information

Chapter 7 Rational and Irrational Numbers

Chapter 7 Rational and Irrational Numbers Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers

More information

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional

More information

All-fiber Optical Parametric Oscillator

All-fiber Optical Parametric Oscillator All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical

More information

Feedback-error control

Feedback-error control Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller

More information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment

More information

A. G. Falkowski Chrysler Corporation, Detroit, Michigan 48227

A. G. Falkowski Chrysler Corporation, Detroit, Michigan 48227 Three-ass mufflers with uniform erforations A. Selamet V. Easwaran The Ohio State University, Deartment of Mechanical Engineering Center for Automotive Research, 06 West 18th Avenue, Columbus, Ohio 4310

More information

A Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split

A Bound on the Error of Cross Validation Using the Approximation and Estimation Rates, with Consequences for the Training-Test Split A Bound on the Error of Cross Validation Using the Aroximation and Estimation Rates, with Consequences for the Training-Test Slit Michael Kearns AT&T Bell Laboratories Murray Hill, NJ 7974 mkearns@research.att.com

More information

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment

More information

3.4 Design Methods for Fractional Delay Allpass Filters

3.4 Design Methods for Fractional Delay Allpass Filters Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or

More information

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures Chater 2 Introductory Concets of Wave Proagation Analysis in Structures Wave roagation is a transient dynamic henomenon resulting from short duration loading. Such transient loadings have high frequency

More information

A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands

A Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands A Secial Case Solution to the Persective -Point Problem William J. Wolfe California State University Channel Islands william.wolfe@csuci.edu Abstract In this aer we address a secial case of the ersective

More information

On split sample and randomized confidence intervals for binomial proportions

On split sample and randomized confidence intervals for binomial proportions On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have

More information

arxiv: v1 [nucl-ex] 28 Sep 2009

arxiv: v1 [nucl-ex] 28 Sep 2009 Raidity losses in heavy-ion collisions from AGS to RHIC energies arxiv:99.546v1 [nucl-ex] 28 Se 29 1. Introduction F. C. Zhou 1,2, Z. B. Yin 1,2 and D. C. Zhou 1,2 1 Institute of Particle Physics, Huazhong

More information

Intrinsic Approximation on Cantor-like Sets, a Problem of Mahler

Intrinsic Approximation on Cantor-like Sets, a Problem of Mahler Intrinsic Aroximation on Cantor-like Sets, a Problem of Mahler Ryan Broderick, Lior Fishman, Asaf Reich and Barak Weiss July 200 Abstract In 984, Kurt Mahler osed the following fundamental question: How

More information

Optical self-energy of superconducting Pb in the terahertz region

Optical self-energy of superconducting Pb in the terahertz region Otical self-energy of suerconducting Pb in the terahertz region T. Mori, 1 E. J. Nicol, 2, * S. Shiizuka, 1 K. Kuniyasu, 3 T. Nojima, 3 N. Toyota, 1 and J. P. Carbotte 4 1 Physics Deartment, Graduate School

More information

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372

More information

Highly improved convergence of the coupled-wave method for TM polarization

Highly improved convergence of the coupled-wave method for TM polarization . Lalanne and G. M. Morris Vol. 13, No. 4/Aril 1996/J. Ot. Soc. Am. A 779 Highly imroved convergence of the couled-wave method for TM olarization hilie Lalanne Institut d Otique Théorique et Aliquée, Centre

More information

Pulse Propagation in Optical Fibers using the Moment Method

Pulse Propagation in Optical Fibers using the Moment Method Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment

More information

Positive decomposition of transfer functions with multiple poles

Positive decomposition of transfer functions with multiple poles Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.

More information

16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE

16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE 16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE H. Yamasaki, M. Abe and Y. Okuno Graduate School at Nagatsuta, Tokyo Institute of Technology 459, Nagatsuta, Midori-ku, Yokohama,

More information

CMSC 425: Lecture 4 Geometry and Geometric Programming

CMSC 425: Lecture 4 Geometry and Geometric Programming CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas

More information

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS 4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis

More information

Probability Estimates for Multi-class Classification by Pairwise Coupling

Probability Estimates for Multi-class Classification by Pairwise Coupling Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics

More information

Dimensional perturbation theory for Regge poles

Dimensional perturbation theory for Regge poles Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West

More information

How to Estimate Expected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty

How to Estimate Expected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty How to Estimate Exected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty Christian Servin Information Technology Deartment El Paso Community College El Paso, TX 7995, USA cservin@gmail.com

More information

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam

More information

Introduction to Landau s Fermi Liquid Theory

Introduction to Landau s Fermi Liquid Theory Introduction to Landau s Fermi Liquid Theory Erkki Thuneberg Deartment of hysical sciences University of Oulu 29 1. Introduction The rincial roblem of hysics is to determine how bodies behave when they

More information

Estimation of Separable Representations in Psychophysical Experiments

Estimation of Separable Representations in Psychophysical Experiments Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)

More information

Asymptotically Optimal Simulation Allocation under Dependent Sampling

Asymptotically Optimal Simulation Allocation under Dependent Sampling Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee

More information

Multiscale Surfaces and Amontons Law of Friction

Multiscale Surfaces and Amontons Law of Friction Tribol Lett (23) 49:539 543 DOI.7/s249-2-94-6 ORIGINAL PAPER Multiscale Surfaces and Amontons Law of Friction J. R. Barber Received: 4 December 22 / Acceted: 24 December 22 / Published online: 26 January

More information

Buckling of Spherical Shells Revisited. John W. Hutchinson. School of Engineering and Applied Sciences, Harvard University.

Buckling of Spherical Shells Revisited. John W. Hutchinson. School of Engineering and Applied Sciences, Harvard University. Buckling of Sherical Shells Revisited John W. Hutchinson School of Engineering and Alied Sciences, Harvard University Abstract A study is resented of the ost-buckling behavior and imerfection-sensitivity

More information

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete

More information

On the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined

On the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined Downloaded 0/9/3 to 99.86.4.8. Redistribution subject to SEG license or coyright; see Terms of Use at htt://library.seg.org/ On the luid Deendence of Rock Comressibility: Biot-Gassmann Refined Leon Thomsen,

More information

Lower bound solutions for bearing capacity of jointed rock

Lower bound solutions for bearing capacity of jointed rock Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying

More information

Notes on Instrumental Variables Methods

Notes on Instrumental Variables Methods Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of

More information

Notes on pressure coordinates Robert Lindsay Korty October 1, 2002

Notes on pressure coordinates Robert Lindsay Korty October 1, 2002 Notes on ressure coordinates Robert Lindsay Korty October 1, 2002 Obviously, it makes no difference whether the quasi-geostrohic equations are hrased in height coordinates (where x, y,, t are the indeendent

More information

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu

More information

Comparative study on different walking load models

Comparative study on different walking load models Comarative study on different walking load models *Jining Wang 1) and Jun Chen ) 1), ) Deartment of Structural Engineering, Tongji University, Shanghai, China 1) 1510157@tongji.edu.cn ABSTRACT Since the

More information

ARTICLE IN PRESS. Physica B

ARTICLE IN PRESS. Physica B Physica B 45 (21) 316 321 Contents lists available at ScienceDirect Physica B journal homeage: www.elsevier.com/locate/hysb On the retrieval of article size from the effective otical roerties of colloids

More information

Pb nanoprecipitates in Al: Magic-shape effects due to elastic strain

Pb nanoprecipitates in Al: Magic-shape effects due to elastic strain Downloaded from orbit.dtu.dk on: Nov 04, 018 nanoreciitates in Al: Magic-shae effects due to elastic strain Hamilton, J.C.; Leoard, F.; Johnson, Erik; Dahmen, U. Published in: Physical Review Letters Link

More information

Section 4: Electromagnetic Waves 2

Section 4: Electromagnetic Waves 2 Frequency deendence and dielectric constant Section 4: Electromagnetic Waves We now consider frequency deendence of electromagnetic waves roagating in a dielectric medium. As efore we suose that the medium

More information

Positivity, local smoothing and Harnack inequalities for very fast diffusion equations

Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Positivity, local smoothing and Harnack inequalities for very fast diffusion equations Dedicated to Luis Caffarelli for his ucoming 60 th birthday Matteo Bonforte a, b and Juan Luis Vázquez a, c Abstract

More information

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,

More information

Equivalence of Wilson actions

Equivalence of Wilson actions Prog. Theor. Ex. Phys. 05, 03B0 7 ages DOI: 0.093/te/tv30 Equivalence of Wilson actions Physics Deartment, Kobe University, Kobe 657-850, Jaan E-mail: hsonoda@kobe-u.ac.j Received June 6, 05; Revised August

More information

Distribution of populations in excited states of electrodeless discharge lamp of Rb atoms

Distribution of populations in excited states of electrodeless discharge lamp of Rb atoms Article Atomic & Molecular Physics June 2013 Vol.58 No.16: 18761881 doi: 10.1007/s11434-013-5789-z Distribution of oulations in excited states of electrodeless discharge lam of Rb atoms TAO ZhiMing 1,2,

More information

An Analysis of Reliable Classifiers through ROC Isometrics

An Analysis of Reliable Classifiers through ROC Isometrics An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit

More information

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009 University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental

More information

Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model

Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Q. Sun, Y.J. Liu, R.H, Jin School of Aeronautics, Northwestern Polytechnical University, Xi an 710072,

More information

The thermal wind 1. v g

The thermal wind 1. v g The thermal win The thermal win Introuction The geostrohic win is etermine by the graient of the isobars (on a horizontal surface) or isohyses (on a ressure surface). On a ressure surface the graient of

More information

Time Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique

Time Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique 2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Time Domain Calculation of Vortex Induced Vibration of Long-San Bridges by Using a Reduced-order Modeling

More information

VISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES

VISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES VISCOELASTIC PROPERTIES OF INHOMOGENEOUS NANOCOMPOSITES V. V. Novikov ), K.W. Wojciechowski ) ) Odessa National Polytechnical University, Shevchenko Prosekt, 6544 Odessa, Ukraine; e-mail: novikov@te.net.ua

More information

Dufour-Driven Thermosolutal Convection of the Veronis Type

Dufour-Driven Thermosolutal Convection of the Veronis Type JOURNAL OF MATEMATICAL ANALYI AND APPLICATION 8, 56958 998 ARTICLE AY975753 Dufour-Driven Thermosolutal Convection of the Veronis Tye ari Mohan Deartment of Mathematics, imachal Pradesh Uniersity, ummer

More information

ON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE

ON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE

More information

Developing A Deterioration Probabilistic Model for Rail Wear

Developing A Deterioration Probabilistic Model for Rail Wear International Journal of Traffic and Transortation Engineering 2012, 1(2): 13-18 DOI: 10.5923/j.ijtte.20120102.02 Develoing A Deterioration Probabilistic Model for Rail Wear Jabbar-Ali Zakeri *, Shahrbanoo

More information

An Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices

An Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem

More information

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is

More information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO) Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment

More information

arxiv: v1 [quant-ph] 22 Apr 2017

arxiv: v1 [quant-ph] 22 Apr 2017 Quaternionic Quantum Particles SERGIO GIARDINO Institute of Science and Technology, Federal University of São Paulo (Unifes) Avenida Cesare G. M. Lattes 101, 147-014 São José dos Camos, SP, Brazil arxiv:1704.06848v1

More information

Multi-Operation Multi-Machine Scheduling

Multi-Operation Multi-Machine Scheduling Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job

More information

Casimir Force Between the Two Moving Conductive Plates.

Casimir Force Between the Two Moving Conductive Plates. Casimir Force Between the Two Moving Conductive Plates. Jaroslav Hynecek 1 Isetex, Inc., 95 Pama Drive, Allen, TX 751 ABSTRACT This article resents the derivation of the Casimir force for the two moving

More information

Kinetics of Protein Adsorption and Desorption on Surfaces with Grafted Polymers

Kinetics of Protein Adsorption and Desorption on Surfaces with Grafted Polymers 1516 Biohysical Journal Volume 89 Setember 2005 1516 1533 Kinetics of Protein Adsortion and Desortion on Surfaces with Grafted Polymers Fang Fang,* Javier Satulovsky, y and Igal Szleifer* *Deartment of

More information

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino

More information

Outline. EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Simple Error Detection Coding

Outline. EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Simple Error Detection Coding Outline EECS150 - Digital Design Lecture 26 Error Correction Codes, Linear Feedback Shift Registers (LFSRs) Error detection using arity Hamming code for error detection/correction Linear Feedback Shift

More information

A Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars

A Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars Aalborg Universitet A Simle And Efficient FEM-Imlementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars Published in: Proceedings of the 9th Nordic Seminar on Comutational

More information

For q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i

For q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:

More information

Author(s) Tohei, T; Kuwabara, A; Yamamoto, T; Citation PHYSICAL REVIEW LETTERS (2005), 94(

Author(s) Tohei, T; Kuwabara, A; Yamamoto, T; Citation PHYSICAL REVIEW LETTERS (2005), 94( Title General rule for dislacive comounds revisited by first hase rinci t Author(s) Tohei, T; Kuwabara, A; Yamamoto, T; Citation PHYSICAL REVIEW LETTERS (2005), 94( Issue Date 2005-01-28 URL htt://hdl.handle.net/2433/39924

More information

Contribution of the cosmological constant to the relativistic bending of light revisited

Contribution of the cosmological constant to the relativistic bending of light revisited PHYSICAL REVIEW D 76, 043006 (007) Contribution of the cosmological constant to the relativistic bending of light revisited Wolfgang Rindler and Mustaha Ishak* Deartment of Physics, The University of Texas

More information

Characteristics of Beam-Based Flexure Modules

Characteristics of Beam-Based Flexure Modules Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced

More information

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management

More information

HENSEL S LEMMA KEITH CONRAD

HENSEL S LEMMA KEITH CONRAD HENSEL S LEMMA KEITH CONRAD 1. Introduction In the -adic integers, congruences are aroximations: for a and b in Z, a b mod n is the same as a b 1/ n. Turning information modulo one ower of into similar

More information

Elliptic Curves and Cryptography

Elliptic Curves and Cryptography Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a

More information

Residual Stress and Deformation Modelling for Metal Additive Manufacturing Processes

Residual Stress and Deformation Modelling for Metal Additive Manufacturing Processes Proceedings of the World Congress on Mechanical, Chemical, and Material Engineering (MCM 215) Barcelona, Sain July 2-21, 215 Paer No. 245 Residual Stress and Deformation Modelling for Metal Additive Manufacturing

More information

A Closed-Form Solution to the Minimum V 2

A Closed-Form Solution to the Minimum V 2 Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month

More information

FUGACITY. It is simply a measure of molar Gibbs energy of a real gas.

FUGACITY. It is simply a measure of molar Gibbs energy of a real gas. FUGACITY It is simly a measure of molar Gibbs energy of a real gas. Modifying the simle equation for the chemical otential of an ideal gas by introducing the concet of a fugacity (f). The fugacity is an

More information

arxiv: v3 [physics.data-an] 23 May 2011

arxiv: v3 [physics.data-an] 23 May 2011 Date: October, 8 arxiv:.7v [hysics.data-an] May -values for Model Evaluation F. Beaujean, A. Caldwell, D. Kollár, K. Kröninger Max-Planck-Institut für Physik, München, Germany CERN, Geneva, Switzerland

More information

Node-voltage method using virtual current sources technique for special cases

Node-voltage method using virtual current sources technique for special cases Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological

More information

A compression line for soils with evolving particle and pore size distributions due to particle crushing

A compression line for soils with evolving particle and pore size distributions due to particle crushing Russell, A. R. (2011) Géotechnique Letters 1, 5 9, htt://dx.doi.org/10.1680/geolett.10.00003 A comression line for soils with evolving article and ore size distributions due to article crushing A. R. RUSSELL*

More information

Mathematics as the Language of Physics.

Mathematics as the Language of Physics. Mathematics as the Language of Physics. J. Dunning-Davies Deartment of Physics University of Hull Hull HU6 7RX England. Email: j.dunning-davies@hull.ac.uk Abstract. Courses in mathematical methods for

More information

Modeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation

Modeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation 6.3 Modeling and Estimation of Full-Chi Leaage Current Considering Within-Die Correlation Khaled R. eloue, Navid Azizi, Farid N. Najm Deartment of ECE, University of Toronto,Toronto, Ontario, Canada {haled,nazizi,najm}@eecg.utoronto.ca

More information