The photon transport equation for turbid biological media with spatially varying isotropic refractive index
|
|
- Rosemary Houston
- 5 years ago
- Views:
Transcription
1 The photon transport equation for turbid biological media with spatially varying isotropic refractive index Malin Premaratne 1, Erosha Premaratne 2 and Arthur James Lowery 1 1 Advanced Computing and Simulation Laboratory AXL) Department of Electrical and Computer Systems Engineering Monash University, Clayton, Victoria 3800, Australia. malin@eng.monash.edu.au and Arthur.Lowery@eng.monash.edu.au 2 Endocrine Unit, Department of Medicine, The Austin Health The University of Melbourne, Victoria 3010, Australia. ah-endo@unimelb.edu.au Abstract: Using the principle of energy conservation and laws of geometrical optics, we derive the photon transport equation for turbid biological media with spatially varying isotropic refractive index. We show that when the refractive index is constant, our result reduces to the standard radiative transfer equation and when the medium is lossless and free of scattering to the well known geometrical optics equations in refractive media Optical Society of America OCIS codes: ) Medical Optics and Biotechnology; ) Light Propagation in Tissue; ) Turbid Media; ) Photon Migration References and links 1. V. Tuchin, Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis, SPIE Press, Bellingham, 2000). 2. J. C. Hebden, A. Gibson, R. M Yusof, N. Everdell, E. M. C. Hillman, D. T. Deply, S. R. Arridge, T. Austin, J. H. Meek and J. S. Wyatt, Three-dimensional optical tomography of the premature infant brain, Phys. Med. Biol. 47, ). 3. X. J. Wang, T. E. Milner, J. F. De Boer, Y. Zhang, D. H. Pashley and J. S. Nelson Characterization of dentin and enamel by use of optical coherence tomography, Appl. Opt. 38, ). 4. S. A. Boppart, M. E. Brezinski, B. E. Bouma, G. J. Tearney and J. G. Fujimoto, Investigation of developing embryonic morphology using optical coherence tomography, Devel. Biol. 177, ). 5. E. L. Hull, M. N. Ediger, A. H. T. Unione, E. K. Deemer, M. L. Stroman, and J. W. Baynes, Noninvasive, optical detection of diabetes: model studies with porcine skin, Opt. Express 12, ), 6. M. H. Niemz, Laser-Tissue Interactions: Fundamentals and Applications, Third, Revised Edition, Springer, New Jersey, 2004). 7. B. R. Masters ed.) Selected Papers on Optical Low-Coherence Reflectometry & Tomography, SPIE Milestone Series MS 165 SPIE Optical Engineering Press, Bellingham, 2001). 8. A. J. Welch and M. J. C. Van-Gemert, Optical-Thermal Response of Laser-Irradiated Tissue Lasers, Photonics and Electro-Optics), Plenum Publishing Corporation, New York, 1995). 9. L-H. Wang, S. L. Jacques and L-Q Zheng, MCML - Monte Carlo modeling of photon transport in multi-layered tissues, Comp. Meth. and Prog. in Biomed. 47, ). 10. S. Chandrasekhar, Radiative Transfer, Dover Publications, New York, 1960). 11. L-H. Wang and S. L. Jacques, Hybrid model of Monte Carlo simulation and diffusion theory for light reflectance by turbid media, J. Opt. Soc. Am. A 10, ). 12. E. D. Aydin, C. R. E. De Oliveira and A. J. H. Goddard, A comparison between transport and diffusion calculations using a finite element-spherical harmonics radiation transport method, Am. Assoc. Phys. Med. 29, ). C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 389
2 13. A. J. Welch, G. Yoon and M. J. C. Van Gemert, Practical models for light distribution in laser-irradiated tissue, Lasers in Surgery and Medicine 6, ). 14. M. Born and E. Wolf, Principles of Optics, 7th Expanded) Edition, Cambridge University Press, New York, 2002). 15. M. Kline and I. W. Kay, Electromagnetic Theory and Geometrical Optics, Wiley Interscience Publishers, New Jersey, 1965). 16. H. A. Ferwerda, The radiative transfer equation for scattering media with a spatially varying refractive index, J. Opt. A: Pure Appl. Opt. 1, L1 L2 1999). 17. T. Khan and H. Jiang, A new diffusion approximation to the radiative transfer equation for scattering media with spatially varying refractive indices, J. Opt. A: Pure Appl. Opt. 5, ). 18. A. Ishimaru, Wave Propagation and Scattering in Random Media, Academic, New York, 1978). 19. M. Kline, A note on the expansion coefficient of geometrical optics, Comm. Pure and Appl. Maths. 14, ). 20. S. I. Grossman, Calculus, Fifth Edition, Harcourt Brace College Publishers, Philadelphia, 1991). 21. G. Yankovsky, Higher Algebra, Mir Publishers, Moscow, 1980). 22. D. J. Griffiths, Introduction to Electrodynamics, Third Edition, Prentice Hall, New Jersey, 1999). 23. W. L. Burke, Applied Differential Geometry, Cambridge University Press, Cambridge, 1985). 1. Introduction Light can penetrate to several centimeters in biological tissue owing to the absence of chromophores that could absorb radiation in the near infrared region NIR) [1]. Recently, information extracted from such photon transport studies have been used increasingly in medical diagnosis [1]. Notably, photon transport in biological media has been used for interesting applications such as imaging of premature infant brain [2], characterization of dentin and enamel [3], investigation of developing embryonic morphology [4] and noninvasive detection of glucose using spectroscopic methods [5]. However, due to strong scattering relative to absorption, imaging and non-invasive characterization in NIR becomes a very difficult task [1, 6]. Therefore, having a clear understanding and a reliable model of photon transport in turbid biological media plays a significant role in improving techniques used for medical diagnosis [1, 7, 8]. Researchers have used variety of numerical strategies for approximately solving standard photon transport equation in biological media such as Monte Carlo simulation techniques [8, 9], discrete ordinate methods [10], diffusion theory [8, 11], finite element methods [12] and various flux models [8] Kubelka-Munk theory [6], four flux model [13], seven flux model [13]). In Monte Carlo techniques, photons are injected from a source either individually or as a bunch or a packet) and traced through the scattering/absorptive medium following photon ray paths [9]. Each photon path incurs a random set of scattering and absorption events. A large number of photons is used to predict an average path. Even though Monte Carlo techniques can provide high accuracy result, due to high computation cost, they are used mainly for validating alternative, faster techniques. In the discrete ordinate method, the transport equation is discretized in different directions or ordinates) and solved using accurate quadrature techniques [10]. However, similar to Monte Carlo techniques, prohibitive computational cost has reduced its usage in practice. The diffusion equation provides an approximation to the transport equation when the number of scattering events is high. Given that biological media is strongly scattering relative to absorption, diffusion equation has the ability to give very accurate results in practice [8, 11]. Apart from requiring a low computational effort, the diffusion equation can also be solved analytically [8]. Recently, finite element methods have also attracted much attention because of their ability to account for complex geometries and inhomogeneous media [12]. It is interesting to observe that up to now, these techniques have been mainly used for solving homogeneous media with constant refractive indices. However, most of the these techniques have the ability to account for refractive index discontinuities using Snell s Law and Fresnel equations for plane boundaries [8, 9]. At a refractive index discontinuity, the incident photon beam or packet is resolved into reflected and refracted wave using Fresnel equations [14] and C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 390
3 traced forward as two separate photon beams or packets with modified intensity or photon number. However, this technique fails to satisfactorily represent the photon transport in turbid media where refractive index continuously changes along the spatial dimensions [14]. This is mainly because when the refractive index is continuously changing along a certain direction, light rays propagating in that direction do not reflect any energy in opposite direction but tend to follow a curved path depending on the gradient of the refractive index along the path [14, 15]. Recently, Ferwerda [16] suggested a modification to standard radiative transport equation to account for media with isotropic, spatially varying continuous refractive index profiles. Khan and Jian [17] identified some redundant terms in Ferwerda s derivation [16] and hence gave a modified, but equivalent form of equation. They also derived the analogous diffusion approximation for Ferwerda s [16] formulation. In this paper, we show that Ferwerda s [16] and hence Khan and Jian s [17] ) formulations do not satisfy energy conservation even in lossless, nonscattering media and fail to account for wavefront curvature. Using the principle of energy conservation [10, 14] and laws of geometrical optics [14, 15], we derive the photon transport equation for turbid biological media with spatially varying isotropic refractive index. We show that when the refractive index is constant, our result reduces to the standard radiative transfer equation [18, 16] and when the medium is lossless and free of scattering to the well known geometrical optics equations in refractive media [19, 15, 14]. This paper is organized as follows. In the Section 2, we track a photon packet propagating through isotropic, continuously spatially varying refractive index profile in scattering/absorptive turbid biological media. In Section 2.1, we show that photon propagation can be viewed as a coordinate transformation or mapping through the turbid medium and hence the Jacobian [20] of the transformation plays a vital part in the characterization of the photon flow. In Section 2.2, using the definition of the Jacobian, we derive its differential properties when the coordinate transformation is governed by the laws of geometrical optics. Using the properties of the Jacobian for ray congruences [14], we derive two different but equivalent formulations of photon transport equation in turbid biological media. The first format given in Section 2.2 is based on the arc length parameter of the photon trajectory or light path) and also demonstrates the conservative features of the result. The second formulation given in Section 2.3 is the standard form used in literature [18, 16, 17, 10] where Cartesian position coordinates and spherical propagation direction coordinates are used. In Section 3, we look at special cases where results are known from independent analysis. As expected, we show in Section 3.1 that the current formulation reduces to the well-known standard photon transport equation [18, 10] when the refractive index is homogeneous. In Section 3.2, we investigate the behavior of the photon transport equation for a lossless and non-scattering medium and show the results completely agrees with what is expected from geometrical optics analysis [19, 15, 14]. In Section 4, we conclude the paper after summarizing results and discussing further research directions. 2. The photon transport equation with spatially varying isotropic refractive index To derive the photon transport equation TPE) with spatially varying refractive index, we track the evolution of an infinitesimally small convex elementary scattering volume along a infinitesimally small ray tube surrounding a central ray in phase space. Figure 1 shows the evolution of initial photon packet with phase space volume, δv S0 at arc length s = s 0. The arc length s is measure from an arbitrary but fixed reference point relative to the laboratory frame) of the central light ray. To proceed further we introduce the rectangular cartesian coordinate system shown in Fig. 2. Positions along the right ray are given by the vector, r = xˆx + yŷ + zẑ where ˆx, ŷ and ẑ are unit vectors along x, y and z directions, respectively. We also assume that tangent vector at position r C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 391
4 { Infinitesimally Small Ray Tube s=0 s=s 0 s s+ s { δv S+ S Central Light Ray { δv S0 δv S Fig. 1. Transport of an infinitesimally small photon packet with phase space volume, V S0, along an infinitesimally small ray tube surrounding a central light ray. ^z θ Ω^ z r s Central Light Ray 0 y ^ y s=s 0 x ϕ ^x Fig. 2. Details of the stationary and moving coordinate systems C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 392
5 is given by ˆΩ Ω x ˆx + Ω y ŷ + Ω z ẑ with polar angle, θ with the ẑ direction and azimuthal angle, ϕ with the ˆx direction as shown in Fig. 2. The unit tangent vector, ˆΩ, can be then written as: ˆΩ = dr ds = 1 µ 2 cosϕ)ˆx + 1 µ 2 sinϕ)ŷ + µẑ 1) where ϕ [0, 2π), θ [0, pi] and µ = cosθ) giving the range µ [ 1, 1]. Suppose L r, ˆΩ,t ) L x,y,z, µ,ϕ,t) is the radiance of the photon flow [10, 18], which describes the amount of energy in the photon beam at position r, flowing per second through a unit area perpendicular to the tangential unit vector ˆΩ. Therefore, the integral L x,y,z, µ,ϕ,t)dxdydzdµdϕ 2) V is proportional to the change in photon number inside the phase space volume V per unit time. Applying this result to the case given in Fig. 1, we track an infinitesimally small photon packet at s = s 0 with phase space volume, δv S0 along the trajectory of a light ray. We basically follow the approach taken by Ferwerda [16] in his derivation of photon transport equation with spatially varying, isotropic refractive index media. However, our approach is systematic and rigorous, leading to a different result. We observe that the photon number in this initial volume can grow due to incoming photons from distributed sources along the light path and photons scattering to the volume from outside, whereas photon number in this initial volume can reduce due to attenuation and scattering of photons to outside. Application of energy conservation principle [10, 18, 16] i.e. photon number conservation in this case due to tracking of a single wavelength) to this moving volume leads, for a small s length along the trajectory, the following result: L x,y,z, µ,ϕ,t + t)dxdydzdµdϕ L x,y,z, µ,ϕ,t)dxdydzdµdϕ + s = s µ a + µ s )L x,y,z, µ,ϕ,t)dxdydzdµdϕ 1 2π + s µ s f µ,ϕ, `µ, `ϕ)l x,y,z, `µ, `ϕ,t)d `µd `ϕ dxdydzdµdϕ 3) s ε x,y,z, µ,ϕ,t)dxdydzdµdϕ + O s 2 ) where µ a µ a x,y,z) is the absorption coefficient per unit length and µ s µ s x,y,z) is the scattering coefficient per unit length. ε εx,y,z, µ,ϕ,t) represents the source distribution per unit volume [16, 18]. The probability of scattering a photon moving in direction `µ, `ϕ) to µ,ϕ) is give by the normalized phase function, f µ,ϕ, `µ, `ϕ) [18, 16]. Noting that this infinitesimal photon packet move along ray congruences [14], the photon packet volume sizes, + s and,atsand s + s, respectively, can be seen as a mapping of initial volume 0 along the ray congruences [14, 15]. In multidimensional space, such mapping happens through the Jacobian, J Jx,y,z, µ,ϕ) [20]. Using this relation, we could get the following expression for the left hand side of Eq.3): L x,y,z, µ,ϕ,t + t)dxdydzdµdϕ L x,y,z, µ,ϕ,t)dxdydzdµdϕ + s C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 393
6 = 0 [L J) s+ s,t+ t L J) s,t ]dx 0 dy 0 dz 0 dµ 0 dϕ 0 4) where x 0,y 0,z 0, µ 0,ϕ 0 ) is the coordinate variables x,y,z, µ,ϕ) at s = s 0. However, for a medium with spatially varying refractive index, n nx,y,z), we have the relation, ds/dt = c/n, where c is the speed of light in vacuum [14]. Using, this, we could get the following relation for the integrand of right hand side of Eq.4): L J) s+ s,t+ t L J) s,t = J n L s + J L s + L J c t s + O s2 ) 5) Now to evaluate Eq.4), we need to evaluate differential of Jacobian with respect to arc length at a certain time instant i.e. J ) Evaluation of term J Suppose x represents the determinant of matrix x. Then we could write the Jacobian J as [20]: x x x x x x 0 y 0 z 0 µ 0 0 y y y y J = x 0 y 0 z 0 0 y z z z µ 0 z z x 0 y 0 z 0 µ 0 0 µ µ µ µ µ x 0 y 0 z 0 µ 0 0 x 0 y 0 z 0 µ 0 0 where x 0,y 0,z 0, µ 0,ϕ 0 ) is the coordinate variables x,y,z, µ,ϕ) at s = s 0. To evaluate this we introduce the permutation function, σ, of the integer set {1,2,,n} [21]. Transposition i, j) of σ is an operation which interchanges i and j but leaves other numbers fixed. It can be shown that every permutation mapping of σ is a combination of transpositions [21]. It can be shown easily that if one representation of a permutation by transposition involves an even number of transpositions, then all representations involving transpositions result in even number of transpositions [21]. An analogous statement applies to the combination of odd numbers of transpositions. Using this property, we denote the sign of a permutation σ, denoted sgnσ) is +1 ifσ is a combination of even number of transpositions and is 1 ifσ is a combination of an odd number of transpositions. Using these definitions, we could write J in Eq.6) as a sum of permutations of σ = {1,2,3,4,5} [21]: J = sgnσ)j 1,σ1) J 2,σ2) J 3,σ3) J 4,σ4) J 5,σ5) 7) σ If ϑ 0 is a one-to-one, onto map from index set {1,2,3,4,5} to variable set x 0,y 0,z 0, µ 0,ϕ 0 ) and ϑ is a similar map to variable set x,y,z, µ,ϕ) then we could write Eq.7) in the following expanded differential form: ϑ j) J = sgnσ) σ j=1 ϑ 0 σ j)) Differentiating Eq.8) respect to s when t is held constant, we get: J 5 ) Ω v ϑk) 5 = sgnσ) v {x 0,y 0,z 0 } σ k=1 ϑk) v j=1 ϑ 0 j) v + v {µ 0,ϕ 0 } 5 sgnσ) σ k=1 5 ϑk) v ϑk) v ϑ j) ϑ 0 σ j)) ) ) 5 j=1 ϑ 0 j) v ϑ j) ϑ 0 σ j)) C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 394 6) 8) 9)
7 Careful inspection of each of the terms under summation index v in Eq.9) shows that only a single term out of five expressions contributes to the final sum due to repeats of identical rows. Thus, we have the following result after some tedious but straightforward calculations: J = Ωx x + Ω y y + Ω z z + µ µ ) + )) J 10) To write this equation in compact form we could use the definitions of divergence operator in a rectangular coordinate system. For the cartesian system ˆx, ŷ, ẑ) the divergence operator of the vector function V = V x ˆx +V y ŷ +V z ẑ has the form [22]: r V = V x x + V y y + V z z 11) Using this definition, we could write Eq.10) as: J = J r ˆΩ + µ µ ) + )) J 12) However, ˆΩ denotes a unit vector in the light propagation direction of the ray s trajectory. It is a well known result in geometric optics that the divergence of ˆΩ is given as the sum of principal radii of curvatures R 1 s) and R 2 s) of the geometrical wave-fronts [19, 14, 15]: r ˆΩ = 1 R 1 s) + 1 R 2 s) 13) It is also well known from geometric optics that the light ray trajectory is given by [14, 15]: n ˆΩ ) = r n 14) where n nx, y, z) is the spatially varying refractive index of the propagation medium. After carrying out the differentiation with respect to s, Eq.14) can be rearranged to get ˆΩ = rn n 1 dn n ds ˆΩ 15) Because ˆΩ is a unit vector i.e. ˆΩ ˆΩ = 1), it is straightforward to show that ˆΩ/ is perpendicular to the vector ˆΩ [20, 22]. Therefore we introduce following unit vectors, ˆµ and ˆϕ in the plane perpendicular to the unit vector ˆΩ: ˆµ = µ cosϕ)ˆx µ sinϕ)ŷ + 1 µ 2 ẑ 16) ˆϕ = sinϕ)ˆx + cosϕ)ŷ 17) Differentiating Eq.1) with respect to s, it possible to derive the following result: ˆΩ = 1 µ 1 µ 2 ˆµ + 1 µ 2 ˆϕ 18) Substituting Eq.18) to Eq.15) and projecting results along ˆµ and ˆϕ, we derive the following expressions for the rate of change of µ and ϕ along the light trajectory: µ = 1 µ 2 rn ˆµ n 19) C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 395
8 = rn ˆϕ n 20) 1 µ 2 The unit vectors ˆµ and ˆϕ in Eq.16) and Eq.17), respectively have the following rate changes along µ and ϕ directions. ˆµ µ = ˆΩ 21) 1 µ 2 ˆϕ = µ ˆµ 1 µ 2 ˆΩ 22) Differentiating Eq.19) and Eq.20) with respect to µ and ϕ, respectively, we get the following results: ) µ µ = µ 1 µ 2 ) µ = 1 µ 2 r n ˆµ rn ˆΩ n n r n ˆµ rn ˆΩ n n Substitution of these results and Eq.13)to Eq.12), we finally get the following compact equation for the differential of Jacobian J along the ray trajectory: 23) 24) J 1 = R 1 s) + 1 ) R 2 s) 2 rn ˆΩ J 25) n Noting that dn ds = rn ˆΩ, Eq.26) can be written in the following equivalent form: J 1 = R 1 s) + 1 R 2 s) 2 ) dn J 26) n ds 2.2. Photon transport equation in light path trajectory arc length format We use the results in previous section to find the photon transport equation along a light path using arc length variable, s. Substitution of Eq.26) and Eq.5) in Eq.4), we get the following relation, integrated over the initial phase space volume, 0 : L x,y,z, µ,ϕ,t + t)dxdydzdµdϕ L x,y,z, µ,ϕ,t)dxdydzdµdϕ + s [ n L 1 = s + c t R 1 s) + 1 ) L + n 2 )] L R 2 s) n 2 Jdx 0 dy 0 dz 0 dµ 0 dϕ 0 + O s 2 ) 27) 0 Recognizing that jacobian, J maps the infinitesimal volume 0 to, we could rewrite Eq.27) as: L x,y,z, µ,ϕ,t + t)dxdydzdµdϕ L x,y,z, µ,ϕ,t)dxdydzdµdϕ + s [ n L 1 = s + c t R 1 s) + 1 ) L + n 2 )] L R 2 s) n 2 dxdydzdµdϕ + O s 2 ) 28) C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 396
9 Using this result in Eq.3) and taking the limit as s 0, we get the integral equation: [ n L 1 + c t R 1 s) + 1 ) L + n 2 )] L R 2 s) n 2 dxdydzdµdϕ µ a + µ s )L + ε)dxdydzdµdϕ 1 2π f µ,ϕ, `µ, `ϕ)l x,y,z, `µ, `ϕ,t)d `µd `ϕ dxdydzdµdϕ = 0 29) µ s 1 0 This result is valid for any arbitrary infinitesimal phase space volume,. Therefore, the integrand of the above equation must be identically equal to zero [22]. nr) L r, ˆΩ,t ) 1 + c t R 1 s) + 1 ) L r, ˆΩ,t ) + n 2 r) R 2 s) n 2 r) =µ a r)+µ s r))l r, ˆΩ,t ) + ε r, ˆΩ,t ) + µ s r) f ˆΩ,Ὼ)L r,ὼ,t ) dὼ 30) where we used the following equivalent relations in function arguments: x,y,z) r; µ,ϕ) ˆΩ and `µ, `ϕ) Ὼ. Eq.30) is the main result of this paper. It is important to note that both r and ˆΩ are functions of s with functional relation given by Eq.14). Therefore, once the spatially dependent, isotropic refractive index relation, nr) is given, photon trajectory can be constructed by solving Eq.14). This result can then be used to calculate the eikonal, S using the relation [14, 15]: r S = nr) ˆΩ 31) Once the eikonal, S is known, principal radii of curvatures, R 1 s) and R 2 s) can be calculated using standard techniques in differential geometry [23] Photon transport equation in r, ˆΩ ) format It is also useful to express the photon transport equation given in Eq.30) using r, ˆΩ ) variables because this is the standard form that researchers have reported results in the past [16, 17, 18]. Moreover, this format is much suited for numerical calculations [10]. Noting that L r, ˆΩ ) /n 2 r) is a function of both r and ˆΩ, we could write: 4π n 2 r) n 2 r) = n 2 r) r L r, ˆΩ,t ) ) r n 2 + n 2 r) ˆΩ L r, r) ˆΩ,t ) ) ˆΩ n 2 r) 32) where we used the definition for gradient operator in the moving orthogonal coordinate system ˆΩ, ˆµ, ˆϕ ) [22]: ˆΩ n 2 = ˆµ 1 µ r) 2 µ n 2 + r) ˆϕ 1 µ 2 n 2 r) 33) C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 397
10 Using the Eq.15) for ˆΩ/ and noting the identities ˆµ ˆϕ = 0, Ω ˆ ˆΩ = 0 and / = ˆΩ r, we could write Eq.32) in following equivalent form: n 2 r) n 2 r) = ˆΩ r L r, ˆΩ,t ) + 1 nr) rnr) ˆΩ L r, ˆΩ,t ) 2 nr) ˆΩ r nr) 34) Substitution of this result to Eq.30) gives the photon transport equation in r, ˆΩ ) format as: nr) L r, ˆΩ,t ) + ˆΩ r L r, ˆΩ,t ) + 1 c t nr) rnr) ˆΩ L r, ˆΩ,t ) 1 + R 1 s) + 1 ) L r, ˆΩ,t ) 2 R 2 s) nr) ˆΩ r nr) =µ a r)+µ s r))l r, ˆΩ,t ) + ε r, ˆΩ ) + µ s r) f ˆΩ,Ὼ)L r,ὼ,t ) dὼ 35) It is instructive compare this result with Ferwerda [16]and so Khang and Jiang [17]). It is easy show that their [16, 17] calculation of r Ω is incorrect because it depends only on differentials of the refractive index and hence does not take wavefront curvature into account. It is a well known fact that even in constant refractive index media, depending on the source configuration, wave fronts can have finite radius of curvature [14, 15]. However, Ferwerda s r Ω term reduces to zero in such cases. Also, his result does not have the term 2/nr)) ˆΩ r nr). This term is essential for making the Eq.35) conserve energy. 3. Special cases To confirm the validity of proposed photon transport equation given in Eq.30) or Eq.35), we investigate its behavior for some special cases where expected outcome is known from independent analysis Photon transport equation in constant refractive index media It is instructive to see the behavior of Eq.35) in a constant refractive index media i.e. nr)=n 0 ) where wavefront curvatures are assumed to be plane. Noting that for such medium, r nr)=0 and both principal radii of curvatures, R 1 s) and R 2 s) become infinite, we get: n 0 L r, ˆΩ,t ) + ˆΩ r L r, ˆΩ,t ) c t = µ a r)+µ s r))l r, ˆΩ,t ) + ε r, ˆΩ,t ) + µ s r) f ˆΩ,Ὼ)L r,ὼ,t ) dὼ 36) This is the well know photon transport equation for homogeneous turbid media with constant refractive index) [10, 18, 16] Lossless non-scattering medium - steady state case In a lossless, non-scattering medium, the photon transport equation in Eq.30) simplifies to: 1 R 1 s) + 1 ) L r, ˆΩ ) + n 2 r) d L r, ˆΩ ) ) R 2 s) ds n 2 = 0 37) r) 4π 4π C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 398
11 If the radiance of the photon beam has a value of L = L 0 at s = s 0 corresponding to r = r 0, then Eq.37) gives: L r, ˆΩ ) = n2 r) n 2 r 0 ) L 0 exp s s 0 1 R 1 s) + 1 ) ) ds R 2 s) This result agrees with the geometric optics result given by Born and Wolf [14] and also Kline [19, 15] for wavefronts with finite radius of curvature. Especially when the wavefronts are plane, both principal radii of curvatures, R 1 s) and R 2 s) become infinite [14, 15]. Therefore, the integral of Eq.38) becomes zero, giving the following result: 38) L r, ˆΩ ) = n2 r) n 2 r 0 ) L 0 39) This result agrees with Ishimaru s result [18] on plane wave photon transport in turbid media. It is interesting to note that none of the results given in this section can be derived from the photon transport equation reported by Ferwerda [16] or it subsequent modified form by Khan and Jiang [17]. 4. Conclusions In this paper, using the principle of energy conservation and laws of geometrical optics, we derived the photon transport equation for turbid biological media with spatially varying isotropic refractive index. We showed that photon propagation through the medium can be viewed as a coordinate transformation or mapping through the turbid medium and hence the Jacobian [20] of the transformation plays vital part in the characterization of the photon flow. Using the properties of the Jacobian for ray congruences [14], we presented two equivalent representation of the transport equation. The first format is based on the arc length parameter of the photon trajectory or the light path) and also demonstrates the conservative features of the result. The second formulation is the standard form used in literature [18, 16, 17, 10] where cartesian position coordinates and spherical propagation direction coordinates are used. Recently, Ferwerda [16] and Khan and Jian [17] made attempts to derive such an equation. However we showed that Ferwerda s [16] and hence Khan and Jian s [17] ) formulations do not satisfy energy conservation and fail to account for wavefront curvature. We showed that when the refractive index is constant, our result reduces to the standard radiative transfer equation [18, 16] and when the medium is lossless and free of scattering to the well known geometrical optics equations in refractive media [19, 15, 14]. C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 399
Multiple-source optical diffusion approximation for a multilayer scattering medium
Multiple-source optical diffusion approximation for a multilayer scattering medium Joseph. Hollmann 1 and ihong V. Wang 1,2, * 1 Optical Imaging aboratory, Department of Biomedical Engineering, Texas A&M
More informationAmbiguity of optical coherence tomography measurements due to rough surface scattering
Ambiguity of optical coherence tomography measurements due to rough surface scattering Y. Ashtamker, 1 V Freilikher, 1,* and J C Dainty 2 1 Department of Physics, Bar-Ilan University, Ramat Gan 52900,
More informationImproved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid medium
246 J. Opt. Soc. Am. A/Vol. 14, No. 1/January 1997 A. Kienle and M. S. Patterson Improved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid
More informationRadiative transfer equations with varying refractive index: a mathematical perspective
Radiative transfer equations with varying refractive index: a mathematical perspective Guillaume Bal November 23, 2005 Abstract This paper reviews established mathematical techniques to model the energy
More informationVector diffraction theory of refraction of light by a spherical surface
S. Guha and G. D. Gillen Vol. 4, No. 1/January 007/J. Opt. Soc. Am. B 1 Vector diffraction theory of refraction of light by a spherical surface Shekhar Guha and Glen D. Gillen* Materials and Manufacturing
More information2 The Radiative Transfer Equation
9 The Radiative Transfer Equation. Radiative transfer without absorption and scattering Free space or homogeneous space I (r,,) I (r,,) r -r d da da Figure.: Following a pencil of radiation in free space
More informationModeling microlenses by use of vectorial field rays and diffraction integrals
Modeling microlenses by use of vectorial field rays and diffraction integrals Miguel A. Alvarez-Cabanillas, Fang Xu, and Yeshaiahu Fainman A nonparaxial vector-field method is used to describe the behavior
More informationAPPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION AND REFLECTION BY A CONDUCTING SHEET
In: International Journal of Theoretical Physics, Group Theory... ISSN: 1525-4674 Volume 14, Issue 3 pp. 1 12 2011 Nova Science Publishers, Inc. APPLICATION OF THE MAGNETIC FIELD INTEGRAL EQUATION TO DIFFRACTION
More informationMCRT L10: Scattering and clarification of astronomy/medical terminology
MCRT L10: Scattering and clarification of astronomy/medical terminology What does the scattering? Shape of scattering Sampling from scattering phase functions Co-ordinate frames Refractive index changes
More informationDiscrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media
Discrete-ordinates solution of short-pulsed laser transport in two-dimensional turbid media Zhixiong Guo and Sunil Kumar The discrete-ordinates method is formulated to solve transient radiative transfer
More informationROINN NA FISICE Department of Physics
ROINN NA FISICE Department of 1.1 Astrophysics Telescopes Profs Gabuzda & Callanan 1.2 Astrophysics Faraday Rotation Prof. Gabuzda 1.3 Laser Spectroscopy Cavity Enhanced Absorption Spectroscopy Prof. Ruth
More informationGreen's functions for the two-dimensional radiative transfer equation in bounded media
Home Search Collections Journals About Contact us My IOPscience Green's functions for the two-dimensional radiative transfer equation in bounded media This article has been downloaded from IOPscience Please
More informationEquivalent isotropic scattering formulation for transient
Equivalent isotropic scattering formulation for transient short-pulse radiative transfer in anisotropic scattering planar media Zhixiong Guo and Sunil Kumar An isotropic scaling formulation is evaluated
More informationAngular distribution of diffuse reflectance in biological tissue
Angular distribution of diffuse reflectance in biological tissue Jinjun Xia and Gang Yao* Department of Biological Engineering, University of Missouri Columbia, Columbia, Missouri 65211, USA *Corresponding
More informationPhoton migration in turbid media using a cumulant approximation to radiative transfer
PHYSICAL REVIEW E, VOLUME 65, 066609 Photon migration in turbid media using a cumulant approximation to radiative transfer M. Xu,* W. Cai, M. Lax, and R. R. Alfano Institute for Ultrafast Spectroscopy
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationFundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation
Fundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation Klaus Jockers November 11, 2014 Max-Planck-Institut für Sonnensystemforschung
More informationarxiv: v1 [physics.class-ph] 8 Apr 2019
Representation Independent Boundary Conditions for a Piecewise-Homogeneous Linear Magneto-dielectric Medium arxiv:1904.04679v1 [physics.class-ph] 8 Apr 019 Michael E. Crenshaw 1 Charles M. Bowden Research
More informationLichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele
Lichtausbreitung in streuenden Medien: Prinzip und Anwendungsbeispiele Alwin Kienle 06.12.2013 Institut für Lasertechnologien in der Medizin und Meßtechnik an der Universität Ulm Overview 1) Theory of
More informationTheoretical Analysis of Thermal Damage in Biological Tissues Caused by Laser Irradiation
Copyright c 007 Tech Science Press MCB, vol.4, no., pp.7-39, 007 Theoretical Analysis of Thermal Damage in Biological Tissues Caused by Laser Irradiation Jianhua Zhou,J.K.Chen and Yuwen Zhang Abstract:
More informationAnalysis of second-harmonic generation microscopy under refractive index mismatch
Vol 16 No 11, November 27 c 27 Chin. Phys. Soc. 19-1963/27/16(11/3285-5 Chinese Physics and IOP Publishing Ltd Analysis of second-harmonic generation microscopy under refractive index mismatch Wang Xiang-Hui(
More informationULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY
Heat Transfer Research 46(9), 861 879 (2015) ULTRAFAST LASER PULSE TRAIN RADIATION TRANSFER IN A SCATTERING-ABSORBING 3D MEDIUM WITH AN INHOMOGENEITY Masato Akamatsu 1,* & Zhixiong Guo 2 1 Graduate School
More informationTowards quantitative tissue absorption imaging by combining photoacoustics and acousto optics
Towards quantitative tissue absorption imaging by combining photoacoustics and acousto optics K. Daoudi, and W. Steenbergen * Biomedical Photonic Imaging group, MIRA Institute for Biomedical Technology
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More information6. LIGHT SCATTERING 6.1 The first Born approximation
6. LIGHT SCATTERING 6.1 The first Born approximation In many situations, light interacts with inhomogeneous systems, in which case the generic light-matter interaction process is referred to as scattering
More informationDetermining how uncertainties in optical properties affect light dose calculations
Determining how uncertainties in optical properties affect light dose calculations Julia Sandell 1, Jarod Finlay, Timothy Zhu 1 Department of Physics, University of Pennsylvania Radiation Oncology, Hospital
More informationLaser Beam Interactions with Solids In absorbing materials photons deposit energy hc λ. h λ. p =
Laser Beam Interactions with Solids In absorbing materials photons deposit energy E = hv = hc λ where h = Plank's constant = 6.63 x 10-34 J s c = speed of light Also photons also transfer momentum p p
More informationScattering of light from quasi-homogeneous sources by quasi-homogeneous media
Visser et al. Vol. 23, No. 7/July 2006/J. Opt. Soc. Am. A 1631 Scattering of light from quasi-homogeneous sources by quasi-homogeneous media Taco D. Visser* Department of Physics and Astronomy, University
More informationPolarized light propagation and scattering in random media
Polarized light propagation and scattering in random media Arnold D. Kim a, Sermsak Jaruwatanadilok b, Akira Ishimaru b, and Yasuo Kuga b a Department of Mathematics, Stanford University, Stanford, CA
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationSTOCHASTIC & DETERMINISTIC SOLVERS
STOCHASTIC & DETERMINISTIC SOLVERS Outline Spatial Scales of Optical Technologies & Mathematical Models Prototype RTE Problems 1-D Transport in Slab Geometry: Exact Solution Stochastic Models: Monte Carlo
More informationNanoscale shift of the intensity distribution of dipole radiation
Shu et al. Vol. 26, No. 2/ February 2009/ J. Opt. Soc. Am. A 395 Nanoscale shift of the intensity distribution of dipole radiation Jie Shu, Xin Li, and Henk F. Arnoldus* Department of Physics and Astronomy,
More informationModeling Focused Beam Propagation in scattering media. Janaka Ranasinghesagara, Ph.D.
Modeling Focused Beam Propagation in scattering media Janaka Ranasinghesagara, Ph.D. Teaching Objectives The need for computational models of focused beam propagation in scattering media Introduction to
More informationAdaptability analysis of radiative transport diffusion approximation in planar-graded-index media
Research Article Adaptability analysis of radiative transport diffusion approximation in planar-graded-index media Advances in Mechanical Engineering 2018, ol. 10(11) 1 6 Ó The Author(s) 2018 DOI: 10.1177/1687814018809613
More informationEITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity
EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity Daniel Sjöberg Department of Electrical and Information Technology Spring 2018 Outline 1 Basic reflection physics 2 Radar cross section definition
More informationLecture 3: Asymptotic Methods for the Reduced Wave Equation
Lecture 3: Asymptotic Methods for the Reduced Wave Equation Joseph B. Keller 1 The Reduced Wave Equation Let v (t,) satisfy the wave equation v 1 2 v c 2 = 0, (1) (X) t2 where c(x) is the propagation speed
More informationSEAFLOOR MAPPING MODELLING UNDERWATER PROPAGATION RAY ACOUSTICS
3 Underwater propagation 3. Ray acoustics 3.. Relevant mathematics We first consider a plane wave as depicted in figure. As shown in the figure wave fronts are planes. The arrow perpendicular to the wave
More informationAnalysis of Scattering of Radiation in a Plane-Parallel Atmosphere. Stephanie M. Carney ES 299r May 23, 2007
Analysis of Scattering of Radiation in a Plane-Parallel Atmosphere Stephanie M. Carney ES 299r May 23, 27 TABLE OF CONTENTS. INTRODUCTION... 2. DEFINITION OF PHYSICAL QUANTITIES... 3. DERIVATION OF EQUATION
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More informationElectrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic
Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an
More informationAT622 Section 2 Elementary Concepts of Radiometry
AT6 Section Elementary Concepts of Radiometry The object of this section is to introduce the student to two radiometric concepts intensity (radiance) and flux (irradiance). These concepts are largely geometrical
More informationReview Article Mathematical Methods in Biomedical Optics
ISRN Biomedical Engineering Volume 2013, Article ID 464293, 8 pages http://dx.doi.org/10.1155/2013/464293 Review Article Mathematical Methods in Biomedical Optics Macaveiu Gabriela Academic Center of Optical
More informationUniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation
Uniform Plane Waves Page 1 Uniform Plane Waves 1 The Helmholtz Wave Equation Let s rewrite Maxwell s equations in terms of E and H exclusively. Let s assume the medium is lossless (σ = 0). Let s also assume
More informationImplicitly causality enforced solution of multidimensional transient photon transport equation
Implicitly causality enforced solution of multidimensional transient photon transport equation Chintha C Handapangoda and Malin Premaratne Advanced Computing and Simulation Laboratory (AχL), Department
More informationRadiative transport in the delta-p 1 approximation: accuracy of fluence rate and optical penetration depth predictions in turbid semi-infinite media
Journal of Biomedical Optics 9(3), 632 647 (May/June 24) Radiative transport in the delta-p 1 approximation: accuracy of fluence rate and optical penetration depth predictions in turbid semi-infinite media
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationAkira Ishimaru, Sermsak Jaruwatanadilok and Yasuo Kuga
INSTITUTE OF PHYSICS PUBLISHING Waves Random Media 4 (4) 499 5 WAVES IN RANDOMMEDIA PII: S959-774(4)789- Multiple scattering effects on the radar cross section (RCS) of objects in a random medium including
More informationOptical Properties of Tissues after Laser Treatments in the Wavelength Range of nm
Optical Properties of Tissues after Laser Treatments in the Wavelength Range of 350-1000 nm Katsunori Ishii, Akinori Kimura, Kunio Awazu Sustainable Energy and Environmental Engineering, Graduate School
More informationUltrafast-laser-radiation transfer in heterogeneous tissues with the discrete-ordinates method
Ultrafast-laser-radiation transfer in heterogeneous tissues with the discrete-ordinates method Zhixiong Guo and Kyunghan Kim Here light propagation and radiation transfer of ultrafast laser pulses in heterogeneous
More informationLecture 19 Optical MEMS (1)
EEL6935 Advanced MEMS (Spring 5) Instructor: Dr. Huikai Xie Lecture 19 Optical MEMS (1) Agenda: Optics Review EEL6935 Advanced MEMS 5 H. Xie 3/8/5 1 Optics Review Nature of Light Reflection and Refraction
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationMAGNETIC FIELDS OF STEADY ELECTRIC CURRENTS, BIOT SAVART LAPLACE LAW, AND ITS APPLICATIONS
MAGNETIC FIELDS OF STEADY ELECTRIC CURRENTS, BIOT SAVART LAPLACE LAW, AND ITS APPLICATIONS In these notes I explain the magnetic fields of steady electric currents. The fields of time-dependent currents
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationLight propagation in structural anisotropic media in the steady-state and time domains
Light propagation in structural anisotropic media in the steady-state and time domains Alwin Kienle, Florian Foschum, Ansgar Hohmann Institut für Lasertechnologien in der Medizin und Meßtechnik, Helmholtzstraße
More informationLecture notes 5: Diffraction
Lecture notes 5: Diffraction Let us now consider how light reacts to being confined to a given aperture. The resolution of an aperture is restricted due to the wave nature of light: as light passes through
More informationDetection and characterization of optical inhomogeneities with diffuse photon density waves: a signal-to-noise analysis
Detection and characterization of optical inhomogeneities with diffuse photon density waves: a signal-to-noise analysis D. A. Boas, M. A. O Leary, B. Chance, and A. G. Yodh Diffusing photons provide information
More informationPC4262 Remote Sensing Scattering and Absorption
PC46 Remote Sensing Scattering and Absorption Dr. S. C. Liew, Jan 003 crslsc@nus.edu.sg Scattering by a single particle I(θ, φ) dφ dω F γ A parallel beam of light with a flux density F along the incident
More informationp(θ,φ,θ,φ) = we have: Thus:
1. Scattering RT Calculations We come spinning out of nothingness, scattering stars like dust. - Jalal ad-din Rumi (Persian Poet, 1207-1273) We ve considered solutions to the radiative transfer equation
More informationModeling Focused Beam Propagation in a Scattering Medium. Janaka Ranasinghesagara
Modeling Focused Beam Propagation in a Scattering Medium Janaka Ranasinghesagara Lecture Outline Introduction Maxwell s equations and wave equation Plane wave and focused beam propagation in free space
More informationCHAPTER 3 CYLINDRICAL WAVE PROPAGATION
77 CHAPTER 3 CYLINDRICAL WAVE PROPAGATION 3.1 INTRODUCTION The phase and amplitude of light propagating from cylindrical surface varies in space (with time) in an entirely different fashion compared to
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationUnit-1 Electrostatics-1
1. Describe about Co-ordinate Systems. Co-ordinate Systems Unit-1 Electrostatics-1 In order to describe the spatial variations of the quantities, we require using appropriate coordinate system. A point
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationSimulation of diffuse reflectance for characterisation of particle suspensions
Thomson, Kelly and Stoliarskaia, Daria and Tiernan-Vandermotten, Sarra and Lue, Leo and Chen, Yi-Chieh (2017) Simulation of diffuse reflectance for characterisation of particle suspensions. In: Optical
More informationScattering of Electromagnetic Radiation. References:
Scattering of Electromagnetic Radiation References: Plasma Diagnostics: Chapter by Kunze Methods of experimental physics, 9a, chapter by Alan Desilva and George Goldenbaum, Edited by Loveberg and Griem.
More informationAdvanced Microscopy. Wintersemester 2012/13 Alexander Heisterkamp Tissue Optics, Optical Clearing, Optical Tomography
Advanced Microscopy Wintersemester 2012/13 Alexander Heisterkamp Tissue Optics, Optical Clearing, Optical Tomography 21.01.2013 [Welch] Radiometric variables Goal: calculation of the distribution of the
More informationLecture Notes on Wave Optics (03/05/14) 2.71/2.710 Introduction to Optics Nick Fang
Outline: A. Electromagnetism B. Frequency Domain (Fourier transform) C. EM waves in Cartesian coordinates D. Energy Flow and Poynting Vector E. Connection to geometrical optics F. Eikonal Equations: Path
More informationInverse Transport Problems and Applications. II. Optical Tomography and Clear Layers. Guillaume Bal
Inverse Transport Problems and Applications II. Optical Tomography and Clear Layers Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University http://www.columbia.edu/ gb23 gb23@columbia.edu
More informationMie theory for light scattering by a spherical particle in an absorbing medium
Mie theory for light scattering by a spherical particle in an absorbing medium Qiang Fu and Wenbo Sun Analytic equations are developed for the single-scattering properties of a spherical particle embedded
More informationPhotonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique
Photonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique Zhigang Chen and Allen Taflove Department of Electrical and Computer
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationWave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces
Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph
More informationTHE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION
Progress In Electromagnetics Research Letters, Vol. 5, 137 149, 2008 THE SCATTERING FROM AN ELLIPTIC CYLINDER IRRADIATED BY AN ELECTROMAGNETIC WAVE WITH ARBITRARY DIRECTION AND POLARIZATION Y.-L. Li, M.-J.
More information1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationNON LINEAR ANOMALOUS SKIN EFFECT IN METALS
www.arpapress.com/volumes/vol7issue3/ijrras_7_3_14.pdf NON LINEAR ANOMALOUS SKIN EFFECT IN METALS Arthur Ekpekpo Department of Physics, Delta State University, Abraka, Nigeria E-mail: arthurekpekpo@yahoo.com
More informationMeasurement of the optical absorption coefficient of a liquid by use of a time-resolved photoacoustic technique
Measurement of the optical absorption coefficient of a liquid by use of a time-resolved photoacoustic technique Yaochun Shen, Zuhong Lu, Stephen Spiers, Hugh A. MacKenzie, Helen S. Ashton, John Hannigan,
More informationSpectral Degree of Coherence of a Random Three- Dimensional Electromagnetic Field
University of Miami Scholarly Repository Physics Articles and Papers Physics 1-1-004 Spectral Degree of Coherence of a Random Three- Dimensional Electromagnetic Field Olga Korotkova University of Miami,
More informationMultiple scattering of light by water cloud droplets with external and internal mixing of black carbon aerosols
Chin. Phys. B Vol. 21, No. 5 (212) 5424 Multiple scattering of light by water cloud droplets with external and internal mixing of black carbon aerosols Wang Hai-Hua( 王海华 ) and Sun Xian-Ming( 孙贤明 ) School
More informationMicrowave-induced thermoacoustic tomography using multi-sector scanning
Microwave-induced thermoacoustic tomography using multi-sector scanning Minghua Xu, Geng Ku, and Lihong V. Wang a) Optical Imaging Laboratory, Biomedical Engineering Program, Texas A&M University, 3120
More information7 Curvilinear coordinates
7 Curvilinear coordinates Read: Boas sec. 5.4, 0.8, 0.9. 7. Review of spherical and cylindrical coords. First I ll review spherical and cylindrical coordinate systems so you can have them in mind when
More informationLecture 21 Reminder/Introduction to Wave Optics
Lecture 1 Reminder/Introduction to Wave Optics Program: 1. Maxwell s Equations.. Magnetic induction and electric displacement. 3. Origins of the electric permittivity and magnetic permeability. 4. Wave
More informationModelling the Directionality of Light Scattered in Translucent Materials
Modelling the Directionality of Light Scattered in Translucent Materials Jeppe Revall Frisvad Joint work with Toshiya Hachisuka, Aarhus University Thomas Kjeldsen, The Alexandra Institute March 2014 Materials
More information(a) Show that the amplitudes of the reflected and transmitted waves, corrrect to first order
Problem 1. A conducting slab A plane polarized electromagnetic wave E = E I e ikz ωt is incident normally on a flat uniform sheet of an excellent conductor (σ ω) having thickness D. Assume that in space
More informationReflection/Refraction
Reflection/Refraction Page Reflection/Refraction Boundary Conditions Interfaces between different media imposed special boundary conditions on Maxwell s equations. It is important to understand what restrictions
More informationSimplified Collector Performance Model
Simplified Collector Performance Model Prediction of the thermal output of various solar collectors: The quantity of thermal energy produced by any solar collector can be described by the energy balance
More informationUniversity of Cyprus. Reflectance and Diffuse Spectroscopy
University of Cyprus Biomedical Imaging and Applied Optics Reflectance and Diffuse Spectroscopy Spectroscopy What is it? from the Greek: spectro = color + scope = look at or observe = measuring/recording
More informationSolutions: Homework 7
Solutions: Homework 7 Ex. 7.1: Frustrated Total Internal Reflection a) Consider light propagating from a prism, with refraction index n, into air, with refraction index 1. We fix the angle of incidence
More informationWave-Front Curvature And Resolution Of The Abraham-Minkowski Controversy
Wave-Front Curvature And Resolution Of The Abraham-Minkowski Controversy Devashish Vikas Gupta Abstract: This paper discusses about the momenta of photons in dielectric media. It is one of the unsolved
More informationElectromagnetic Waves
Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic
More informationMotion control of the wedge prisms in Risley-prism-based beam steering system for precise target tracking
Motion control of the wedge prisms in Risley-prism-based beam steering system for precise target tracking Yuan Zhou,,2, * Yafei Lu, 2 Mo Hei, 2 Guangcan Liu, and Dapeng Fan 2 Department of Electronic and
More informationOffset Spheroidal Mirrors for Gaussian Beam Optics in ZEMAX
Offset Spheroidal Mirrors for Gaussian Beam Optics in ZEMAX Antony A. Stark and Urs Graf Smithsonian Astrophysical Observatory, University of Cologne aas@cfa.harvard.edu 1 October 2013 This memorandum
More information- 1 - θ 1. n 1. θ 2. mirror. object. image
TEST 5 (PHY 50) 1. a) How will the ray indicated in the figure on the following page be reflected by the mirror? (Be accurate!) b) Explain the symbols in the thin lens equation. c) Recall the laws governing
More informationSolar radiation / radiative transfer
Solar radiation / radiative transfer The sun as a source of energy The sun is the main source of energy for the climate system, exceeding the next importat source (geothermal energy) by 4 orders of magnitude!
More informationOptimized diffusion approximation
356 Vol. 35, No. 2 / February 208 / Journal of the Optical Society of America A Research Article Optimized diffusion approximation UGO TRICOLI, CALLUM M. MACDONALD, ANABELA DA SILVA, AND VADIM A. MARKEL*
More informationExperiment 8. Fresnel Coefficients. 8.1 Introduction. References
Experiment 8 Fresnel Coefficients References Optics by Eugene Hecht, Chapter 4 Introduction to Modern Optics by Grant Fowles, Chapter 2 Principles of Optics by Max Born and Emil Wolf, Chapter 1 Optical
More informationChapter 1 Introduction
Introduction This monograph is concerned with two aspects of optical engineering: (1) the analysis of radiation diffraction valid for all wavelengths of the incident radiation, and (2) the analysis of
More informationMonte Carlo Methods:
Short Course on Computational Monte Carlo Methods: Fundamentals Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine Shuang Zhao 1 Teaching Objective Introducing
More information