arxiv: v1 [astro-ph.sr] 17 Sep 2014

Transcription

1 Depolarization of multiple scattered light in atmospheres due to anisotropy of small grains and molecules. II. The problems with sources arxiv:49.547v [astro-ph.sr] 7 Sep 4 N. A. Silant ev, G. A. Alekseeva, V. V. Novikov Central Astronomical Observatory at Pulkovo of Russian Academy of Sciences, 964, Saint-Petersburg, Pulkovskoe shosse 65, Russia September, 4 Abstract In the previous paper we considered two classic problems - the diffuse reflection of the light beam from semi-infinite atmosphere, and the Milne problem. For both problems we used the technique of invariance principle. In this paper we consider the solution of the problem when in semi-infinite atmosphere the sources of unpolarized radiation are exist. Here we used the technique of the Green matrices. We consider only continuum radiation. Keywords: Radiative transfer, scattering, polarization Introduction Real atmospheres of stars and planets, and optically thick accretion disks consist of anisotropic particles. Usually these particles dust grains and molecules are small compared to the wavelengths of radiation. To consider the process of multiple light scattering in such atmospheres we have to remember the features of single scattering. The incident electromagnetic wave Eω induces dipole moment in the scattering particle, p i ω = β ij ωe j ω. This dipole moment gives rise to scattered radiation. The polarizability tensor β ij ω characterizes the structure of the dust grains and molecules. Anisotropic particle with axial symmetry is characterized by two polarizabilities - along the symmetry axis β ω, and in transverse direction - β ω. For isotropic particles electron, spherical dust grains polarizabilities β = β. It is known see, for example, Dolginov et al. 995, Gnedin and Silant ev 997 that for β = β the polarization of scattered radiation is greater than that for anisotropic particles with β β. The anisotropic scattering particles due to chaotic thermal motions are freely oriented in an atmosphere. The radiative transfer equation for this case firstly was derived by Chandrasekhar 96. The average over all orientations of particles results to the additional term, depending on anisotropy of particles. This term depends integrally on the intensity Iτ, of radiation, and gives rise to additional isotropic unpolarized part of scattered radiation. For this reason, the depolarizing effect of anisotropy of particles mostly reveals by consideration of axially symmetric problems. One of such problem is the multiple scattering of light in an semi-infinite atmosphere with sources of unpolarized radiation. In such problem the radiation is described by two intensities - I l τ, and I r τ,. Here τ is the optical depth below the surface of semi-infinite plane-parallel atmosphere, = cos ϑ with ϑ being the angle between the outer normal N to the surface and the direction of light propagation n. The intensity I l describes the light linearly polarized in the plane nn, and I r has polarization perpendicular to this plane. The total intensity I = I l + I r, and the Stokes parameter Q = I l I r. The Stokes parameter U. The degree of linear polarization is equal to p = I l I r /I l +I r. Circularly polarized light is described by the separate equation. We do not consider this equation. A detailed consideration of problems without parameter b is presented in many papers see, for example, Chandrasekhar 96; Horak & Chandrasekhar 96; Lenoble 97; Abhyankar & Fymat 97. nsilant@bk.ru

2 The consideration of diffuse reflection of the light beam from semi-infinite atmosphere and the Milne problem, which take into account parameter b, is presented in paper Silant ev et al. 4. These problems were considered using the principle of invariance. The problems with sources of nonpolarized radiation is impossible consider in such a manner. We consider these problems using the Green matrix approach. Below we present of theorethical part of this consideration. Radiative transfer equation for intensities I l τ, and I r τ, Introducing the column vector I with the components I l,i r we obtain the following matrix transfer equation for multiple scattering of continuum radiation on small anisotropic particles see Chandrasekhar 96; Dolginov et al. 995: diτ, = Iτ, b ˆM Kτ b K τ q Qτ, where q = N a σ a /N s σ s +N a σ a is the probability oflight absorption, N a and N s are the number densities of absorbing the grains and scattering particles, respectively; = N s σ s + N a σ a dz determines the dimensionless optical depth, σ s and σ a are the cross sections of scattering and absorption. The scattering cross section of small particles dust grains, molecules is σ s = π/ω/c 4 b +b, whereω = πν is cyclic frequencyoflight, c is the speed oflight. For freelychaoticallyorientedparticles σ s is independent of the polarization of incident electromagnetic wave Eω. The values b and b are related to polarizability tensor β ij ω of a particle as a whole. Induced dipole moment of a particle, as whole, is equal to p i ω = β ij ωe j ω. Anisotropic particle with axial symmetry is characterized by two polarizabilities - along the symmetry axis β ω, and in transverse direction β ω. For such particles b = 9 β +β + 45 β β, b = 5 β β. In transfer equation we use the dimensionless parameters b = b b, b =, b +b =. b +b b +b For needle like particles β β parameters b =.4, b =., and for plate like particles β β we have b =.7, b =.. Parameter b describes the depolarization of radiation, scattered by freely oriented anisotropic particles. The matrix ˆM has the form see Lenoble 97; Abhyankar and Fymat 97: ˆM =,. 4, Note that the superscript T will be used for matrix transpose. The vector Kτ and scalar K τ are determined by the expressions: Kτ = q K τ = q 4 It is easy seen from these expressions that ˆM T Iτ,, 5 I l τ,+i r τ,. 6 K τ = K lτ+ K r τ. 7

3 Usual source of radiation is the thermal radiation Qτ = B ω Tτ, where B ω T is the Planck function with the temperature, depending on the optical depth τ. This dependence appears as a result of consideration of various models for atmospheres. The main features of our investigation are consideration of problems with depolarization parameter b, and taking into account the true absorbtion parameter q i.e., the existence of absorbing grains in an atmosphere. The integral equation for Kτ The radiative transfer equation can be written in the form: diτ, = Iτ, Sτ,, 9 where the vector Sτ, is described by the expression: Sτ, =  Kτ+ q Qτ. Here the matrix  appears from the relation:  Kτ b ˆM Kτ+b K lτ+ K r τ. The explicit form of matrix  is the following:  b = +b, b +b b +b, b. The vector Iτ,, according to Eq.9, can be written as: Iτ, = Iτ, = τ τ exp τ τ Sτ,,, exp τ τ Sτ,,. Substitution these expressions into Eq.5, and taking into account that matrix ˆM T and vector S are even functions of, we derive the expression: Kτ = q ˆM T Sτ, exp τ τ. 4 Using the formulas and, we obtain the integral equation for Kτ: Kτ = gτ+ The matrix kernel of this integral equation is the following: ˆL τ τ = q ˆL τ τ Kτ. 5 exp τ τ â. 6

4 Here ˆM T  = /4â. The explicit form of matrix â is the following: â = b + 4 +b +, b +b b + b, b +b. 7 The free term gτ has the form: gτ = q q Qτ exp τ τ +. 4 Solution of integral equation 5 using the Green matrix According to the standard theory of integral equations see, for example, Tricomi 957, the solution of equation 5 can be presented in the form: Kτ = gτ+ where the Green matrix ˆRτ,τ obeys the integral equation: It follows from this equation that ˆRτ,τ = ˆL τ τ + ˆRτ, = ˆLτ+ ˆRτ,τ gτ, 9 ˆL τ τ ˆRτ,τ. ˆL τ τ ˆRτ,. This equation coincides formally with Eq.5, and, consequently, ˆRτ, can be written in the form: ˆRτ, = ˆLτ+ ˆRτ,τ ˆLτ q [ exp τ + ˆRτ, ] â, where ˆRa denotes the Laplace transform: ˆRa = As is known, the Green matrix ˆRτ,τ also obeys the another equation: ˆRτ,τ = ˆL τ τ + From this equation we obtain ˆR,τ = ˆLτ+ exp aτˆrτ. ˆRτ,τ ˆL τ τ. 4 ˆR,τ ˆL τ τ. 5 4

5 Analogously to equation, the ˆR,τ can be written in the form: ˆR,τ = ˆLτ+ q Differentiation of Eq. or 4 gives rise to relation: ˆRτ,τ τ ˆLτ ˆRτ,τ [ â exp τ + ˆR ],τ. 6 + ˆRτ,τ τ = ˆRτ,ˆR,τ. 7 This relation means that the total Green matrix ˆRτ,τ can be obtained, if we know two auxiliary matrices, depending on one variable - τ or τ. The Laplace transform over τ and over τ of the matrix ˆRτ,τ we denote as ˆRa,b = exp aτexp bτ ˆRτ,τ. Taking this double Laplace transform in equation 7, we obtain the relation: ˆRa,b = a+b [ ˆRa,+ ˆR,b+ ˆRa, ˆR,b]. 9 Taking the Laplace transform over variable τ in equations and 6, and using the relation 9, we obtain the system equations for ˆRa, and ˆR,a. It is useful to introduce the new matrices Ĥ b = Î + ˆR,b, For these new matrices we obtain the following system of equations: Ĥ b = Î + q Ĥ b = Î + q Ĥ b = Î + ˆRb,. â b+ĥ Ĥ b, b+ĥ b Ĥ â. For the sake of convenience, we introduce the notation /b = z. Then equations acquire more simple form: Ĥ z = Î + q z â +Ĥ z Ĥ z, Ĥ z = Î + q z +Ĥ z zĥ â. It is seen, that to know Ĥ z and Ĥ z we are to know Ĥ and Ĥ. Taking z =, we derive from the system the separate system of equations for Ĥ and Ĥ : Ĥ = Î + q â + Ĥ Ĥ, 5

6 Ĥ = Î + q + ĤĤ â. Thus, the system plays the main role in the Green matrix approach. Note that for the case b = the matrix â = â T. In this case we have Ĥ = Ĥ T. The system is nonlinear system of matrix equations. It is known see Sobolev 969 that scalar H-function also obeys the linear equation. Let us derive the linear equations for our matrix functions Ĥ and Ĥ. 5 Linear equations for matrix functions Ĥ and Ĥ Let us multiply equation for H by the value â /z from the right side of equation and take the integration over in the interval,. For brevity, we introduced λ = q. As a result, we obtain the expression: Ĥ â z = The same operation with z + gives: Ĥ â z + = â z +λ z â z + +λ z â + z Ĥ Ĥ â. 4 â + z +Ĥ Ĥ â. 5 Now let us multiply equation for Ĥ by the value â /z from the left side of equation and take the usual integration over in the interval,. As a result, we obtain the expression: â Ĥ z = The same operation with z + gives: â Ĥ z + = â z +λ z â z + +λ z The summation of Eqs.4 and 7 gives rise: â + z Ĥ Ĥ â. 6 â + z +Ĥ Ĥ â. 7 Ĥ â z + â Ĥ z + = â z + â z + + λ z λ z â + z Ĥ Ĥ â + â + z +Ĥ Ĥ â. In last double integral we take the substitutions and. After that the sum of two last integrals in Eq. acquires the form: λ z â z + z Ĥ Ĥ â. 9 6

7 From first equation in system it follows that â Ĥ z + = Î Ĥ z. 4 This integration exists in Eq.9 and in left part of Eq.. Using the relation 4 in Eqs. and 9, we obtain the following linear equation: Ĥ z Î â Ĥ z = Î â. 4 z The summation of equations 5 and 6 and using the relation Ĥ â z + = Î Ĥ z, 4 whichfollowsfromthe secondequationin system, weobtainthe followinglinearequationforĥ z: Î â Ĥ z z = Î â Ĥ. 4 z Equations 4 and 4 are the equations with singular kernels. The solution such scalar equations is presented in the book Gakhov 977. Some times the expression Ĥ z is used. To derive this expression, we represent formula 7 in the form: 44 â = â +â +â 4 4, 45 where the matrices â,â and â 4 are independent of : b +b â =, b, b, b +b b, b â = b, b +b, 4b b b, b â 4 = b, b. 46 Using these relations, we obtain the following expression: Ĥ z = Ĥ z Î + â âz z âz [ ] Ĥ +zĥ â +Ĥ +zĥ +z Ĥ â 4 âz. 47 Here Ĥ n denote the moments of the matrix Ĥ : Ĥ n = Remember that inverse X matrix has very simple form: ˆB = B, B B, B 4 n Ĥ. 4 = B B 4 B B B4, B B, B. 49 7

8 6 Formulas for outgoing radiation According to Eq. 9 and, the vector I,, describing the outgoing radiation, has the form: I, = exp τ Sτ, exp τ  Kτ+ q Qτ, 5 i.e. this expression proportional to the Laplace transform over variable τ/. The vector Kτ is presented in Eq.9, where the vector gτ is given in Eq.. The Laplace transform of vector Kτ can be written in the form: K = exp τ Kτ exp τ Î + ˆR,τgτ. 5 Thus, expression 5 acquires the form: I, = g  + ˆR,τgτ + q Q. 5 Note that the value Qτ characterises the distribution of sources of nonpolarized radiation in an atmosphere. For source function Q h τ = Q h exp hτ the expression for I, acquires comparatively simple form. In this case the expression depends on Ĥ and Ĥ, not on matrix ˆRτ,τ. The sources of types Q n τ = Q n τ n are related with exponential source by the simple formula: It means that the sources of type gτ = q qq h Q n τ = n Q n d n dh n exp hτ h=. 5 Qτ = Q h exp hτ+q +Q τ +Q τ can be considered on the base of exponential source. Consider this case in detail. Taking Q h τ = Q h exp hτ in Eq. we obtain exp τ exp hτ + h + exp hτ. h+ 55 Substituting this expression into Eq. 5 and taking into account formula 9 for double Laplace transform of matrix ˆRτ,τ, we derive the following expression: q I, = Q h h+ + q qâ Ĥ Ĥ h h h+ + Ĥ h ] + h +h+ [ Ĥ h + For homogeneous source Q we are take h =. In this case we obtain from Eq. 56: q [ I, = Q + q qâ Ĥ Ĥ + + ] + Ĥ }. 56 }. 57

9 Taking z = in first equation of system, we obtain: Ĥ = [ Î q ] â Ĥ +â Ĥ +â 4 Ĥ 4. 5 The matrices â, â and â 4 are presented in Eq. 46. Definition of moments Ĥ n 4. The outgoing radiation for the source Q τ can be derived according to Eq.5: I, = Q q + q qâ Ĥ Ĥ + + Analogously can be derived formulas for outgoing radiation for sources of type τ,τ etc. 7 Conclusion is given in Eq.. 59 In this paper we consider the role of anisotropy of molecules and dust grains in the process of depolarization of multiple scattered light in semi-infinite plane-parallel atmospheres with the sources of unpolarized radiation. The consideration uses the matrix Green function Ĥτ,τ approach. This matrix is connected with two auxiliar matrices, depending on one variable τ or τ, which obey the system of nonlinear equations similar to known system of equations for scalar H- functions of Chandrasekhar. It is shown that for sources of type exp hτ and sources of type τ,τ,τ etc. the formulas for outgoing radiation depend on matrix H-functions and have rather simple form. This research was supported by the Program of Prezidium of RAS No, the Program of the Department of Physical Sciences of RAS No7, the Federal Target Program Science and Scientific-Pedagogical Personnel of Innovative Russia XXXVII in turn - the action.., and by the Grant from President of the Russian Federation The Basic Scientific Schools NSh References [] Abhyankar, K. D., Fymat, A. L. Astrophys. J. Suppl. Ser., 5 97 [] Chandrasekhar, S.: Radiative transfer. Dover, New York 96 [] Dolginov, A. Z., Gnedin, Yu. N., Silant ev, N. A. : Propagation and Polarization of Radiation in Cosmic Media. Gordon & Breach Publ., Amsterdam 995 [4] Gakhov, F. D.: Boundary problems, Moscow, Nauka, 977 in russian [5] Gnedin, Yu. N., Silant ev, N. A. Astronomy and Astrophys. Review, 997 [6] Horak, H. G., Chandrasekhar, S. Astropys. J. 4, [7] Lenoble, J., J. Quant. Spectrosc. Radiat. Transf., 5 97 [] Silant ev, N. A., Alekseeva, G. A., Novikov, V. V. Astropys. Space Sci. 5, [9] Sobolev, V. V. : Course in theorethical astrophysics. NASA Thechnical Translation F-5, Washington 969 [] Tricomi, F. G.: Integral equations, Intersci. Publishers, inc., New York 957 9