Effects of scattering and dust grain size on the temperature structure of protoplanetary discs: a three-layer approach

The temperature in the optically thick interior of protoplanetary discs is essential for the interpretation of millimetre observations of the discs, for the vertical structure of the discs, for models of the disc evolution and the planet formation, and for the chemistry in the discs. Since large icy grains have a large albedo even in the infrared, the effect of scattering of the diffuse radiation in the discs on the interior temperature should be examined. We have performed a series of numerical radiation transfer simulations, including isotropic scattering by grains with various typical sizes for the diffuse radiation as well as for the incident stellar radiation. We also have developed an analytic model including isotropic scattering to understand the physics concealed in the numerical results. With the analytic model, we have shown that the standard two-layer approach is valid only for grey opacity (i.e. grain size ≳10 μm) even without scattering. A three-layer interpretation is required for grain size ≲10 μm. When the grain size is 0.1–10 μm, the numerical simulations show that the isotropic scattering reduces the temperature of the disc interior. This reduction is nicely explained by the analytic three-layer model as a result of the energy loss by scatterings of the incident stellar radiation and of the warm diffuse radiation in the disc atmosphere. For grain size ≳10 μm (i.e. grey scattering), the numerical simulations show that the isotropic scattering does not affect the interior temperature. This is nicely explained by the analytic two-layer model; the energy loss by scattering in the disc atmosphere is exactly offset by the 'green-house effect' due to the scattering of the cold diffuse radiation in the interior.     

Polarization of resonance lines in the case of a Voigt absorption profile

Multiple resonance scattering of radiation in a spectral line is considered in the case of a Voigt absorption profile. The scattering is assumed to take place in a nonmagnetic semi-infinite atmosphere with uniformly distributed sources of unpolarized radiation. Polarization characteristics have been obtained for the emergent radiation by numerically solving the Ambartsumian-Chandrasekhar matrix integral equation.     

Depolarization of multiple scattered light in atmospheres due to anisotropy of small grains and molecules

Freely oriented small anisotropic grains and molecules depolarize radiation both in single scattering and in the process of multiple scattering. Especially large depolarization occurs for resonant scattering corresponding to the electron transitions between the energy levels with very different quantum numbers. The existence of light absorption also changes essentially the angular distribution and polarization of radiation, outgoing from an atmosphere. In the present paper we consider these effects in detail both for continuum radiation and for resonant lines. Because the term describing the depolarization deals with isotropic radiation, we consider the axially symmetric part of radiation. We derived the formulas for observed intensity and polarization using the invariance-principles both for continuum and resonant scattering. We confine ourselves to two problems—the diffuse reflection of the light beam from semi-infinite atmosphere, and the Milne problem.     

Polarized line formation by resonance scattering: Lorentz profile

Multiple resonance scattering of spectral line radiation is examined in atmospheres with uniformly distributed sources of unpolarized radiation. It is assumed that the profile of the absorption coefficient is lorentzian and that scattering involves complete frequency redistribution. The polarization characteristics of the emerging radiation are determined by iterative solution of a nonlinear Ambartsumyan-Chandrasekhar matrix integral equation. In particular, it is found that for pure scattering the maximum polarization at the limb of the disk is 5.421%. The polarization characteristics of the emerging radiation are compared for three different absorption profiles: Lorentz, Doppler, and rectangular (monochromatic radiation). __________ Translated from Astrofizika, Vol. 50, No. 2, pp. 199–217 (May 2007).     

Radiative transfer in a semi-infinite plane parallel atmosphere with an inhomogeneous source function

We investigate the multiple scattering of radiation in semi-infinite homogeneous atmosphere when the sources of the radiation are distributed inhomogeneous, for example, are created by restricted beams penetrating into the medium. The case of isotropic scattering is considered. It is shown that the density of radiation and the intensity of outgoing radiation for any forms of the sources can be represented as some integrals with the real and imaginary parts of the universal H-function, which satisfies the nonlinear integral equation. We calculated the intensity of radiation emerging from the surface after multiple scattering for the case when a beam with a finite radius incident perpendicular on the medium surface. The results allowed us to estimate quantitatively when the intensity of outgoing radiation in the center of a beam coincides with that for the classical case of unbounded flux (the case considered by Chandrasekhar et al.). We compared our exact solutions with those in the diffusion approximation. For conservative medium the difference is ?20–30%, depending on the particular forms of the radiation sources. For absorbing medium the difference is much larger. Our exact semi-analytical solution can be generalized for the cases of multiple anisotropic scattering of the polarized beams. The presented simple theory can be used at the consideration of close binary systems, flare stars etc.     

Stokes parameters for Thomson scattering in a strong magnetic field with radiation damping

The effect of a strong magnetic field on neutron stars or white dwarfs is calculated for Thomson scattering in a strong magnetic field with radiation damping. The Stokes's parameters for the scattered radiation are computed explicitly in terms of the state of polarization of the incident wave, the electrocyclotron frequency, the angle of incidence, and the angle of scattering. The criterion for the magnetic field to substantially affect the Stokes's parameters is that the photon frequency be less than the electrocyclotron frequency. The effects of classical radiation damping are explicitly taken into account by using the Abraham-Lorentz equation. The corresponding quantum mechanical treatment for Compton scattering in a relativistic electron gas is now under investigation.     

One more explanation of superluminal motion

The occurrence of superluminal motion in extragalactic radio sources is believed to be quite common. Among others, the geometrical scattering of radio radiation can also cause superluminal expansion and or motion and halo formation, In this paper, the effectiveness of the stimulated Raman scattering in producing these features is investigated. The scattering medium is a plasma whose position, density and temperature decide the rate and angle of scattering. When the radiation from a stationary and constant source gets scattered from a stationary plasma, a halo is formed around the source. However, the scattering of a rotating radiation beam does produce superluminal motion of the virtual source. It is found that the plasma should have the characteristics of the emission-line regions and the intercloud medium in order to Raman scatter the radiation. Since the scattering is polarization dependent, it is possible to estimate the rotation of the electric vector along the direction of the apparent motion of a radio source.     

Average photon path-length of radiation emerging from finite atmospheres

The determination of the average path-length of photons emerging from a finite planeparallel atmosphere with molecular scattering is discussed. We examine the effects of polarisation on the average path-length of the emergent radiation by comparing the results with those obtained for the atmosphere where the scattering obeys the scalar Rayleigh function. Only the axial radiation field is considered for both cases.To solve this problem we have used the integro-differential equations of Chandrasekhar for the diffuse scattering and transmission functions (or matrices). By differentiation of these equations with respect to the albedo of single scattering we obtain new equations the solution of which gives us the derivatives of the intensities of the emergent radiation at the boundaries.As in the case of scalar transfer the principles of invariance by Chandrasekhar may be used to find an adding scheme to obtain both the scattering and transmission matrices and their derivatives with respect to the albedo of single scattering. These derivatives are crucial in determining the average path length.The numerical experiments have shown that the impact of the polarisation on the average pathlength of the emergent radiation is the largest in the atmospheres with optical thickness less than, or equal to, three, reaching 6.9% in the reflected radiation.     

Effect of continuum absorption on the polarization of resonance lines

Multiple resonance scattering of radiation in a spectrum line with absorption in the continuum is examined. It is assumed that the scattering atmosphere is semi-infinite and that there is no magnetic field or continuum emission at the frequencies of the spectrum line. The polarization characteristics of the emerging radiation are determined for unpolarized primary radiation sources distributed uniformly within an atmosphere in the case of a Voigt absorption profile. The calculations employ an iterative solution of the Ambartsumyan-Chandrasekhar matrix integral equation.     

Asymmetry parameter and efficiency for radiation pressure and scattering of electromagnetic radiation by a very large dielectric sphere

The asymmetry parameter of the scattering diagram and the efficiency for radiation pressure for scattering by a very large dielectric sphere conforming to geometrical optics and diffraction, have been investigated beyond the limit of existing literature. A new feature of shallow broad extremum (resonance) has been found in the asymmetry parameter as well as in the radiation pressure both as function of the index of refraction.     

A new representation of Güttler's theory of electromagnetic scattering by a composite sphere

The multipole expansion coefficients for electromagnetic scattering by a stratified core-mantle sphere, suitable for fast and accurate evaluation of the radiation scattering parameters, have been recast.     

On the theory of X-ray pulsar radiation

The origin of hard X-ray spectrum (continuum and cyclotron lines) of pulsars in binary systems is discussed. A model of the polar region of a neutron star consisting of a hot spot in a dense plasma atmosphere with a quasi-homogeneous magnetic field and an extended accreting column in an inhomogeneous dipolar field is investigated. In the hot spot bremsstrahlung and Thomson scattering form continuum radiation, while bremsstrahlung and cyclotron scattering produce the absorption cyclotron lines. By the observed continuum intensity one can estimate the maximum distances to pulsars. Cyclotron scattering in gyro-resonant layers localized in the accreting column leads to a general attenuation of the radiation of a hot spot, but is unable to ensure the formation of cyclotron lines. For strong accretion the hot spot radiation becomes insignificant, the lines disappear and the pulsating component of an X-ray pulsar is produced by the accreting column bremsstrahlung transformed by Thomson scattering.     

Scattering of waves in a magnetoactive plasma

A general theory of scattering of waves in a magnetoactive plasma by particles of arbitrary energy is presented. The cross-section for the scattering of magnetoionic waves by thermal particles is derived and discussed. Conditions under which the effect of the spiralling motion of the scattering electron can be neglected in treating inverse Compton radiation are found.     

Radiative transfer in an infinite cylinder

The problem of radiation transfer in a cylinder with diffuse reflectivity and containing an energy source is connected with the source-free radiation transfer problem with isotropic boundary condition. Equation for the radiation heat flux is obtained for a polynomial source. In the special case of isotropic scattering, the radiation heat flux is given in terms of the albedo of the second problem. An expression is also given for the net radiation heat flux.     

A new method for exact solution of transfer equations in finite media

In this paper we develop a new exact method combined with finite Laplace transform and theory of linear singular operators to obtain a solution of transport equation in finite plane-parallel steady-state scattering atmosphere both for angular distribution of radiation from the bounding faces of the atmosphere and for intensity of radiation at any depth of the atmosphere. The emergent intensity of radiation from the bounding faces are determined from simultaneous linear integral equations of the emergent intensity of radiation in terms ofX andY equations of Chandrasekhar. The intensity of radiation at any optical depth for a positive and negative direction parameter is derived by inversion of the Laplace transform in terms of intergrals of the emergent intensity of radiation. A new expression of theX andY equation is also derived for easy numerical computation. This is a new and exact method applicable to all problems in finite plane parallel steady scattering atmosphere.     

Polarization in spectral lines

The general formalism, presented in a previous paper of this series (Landi Degl'Innocenti, 1983a), is particularized to deduce the radiative transfer equations for polarized radiation and the statistical equilibrium equations for a multi-level atom in the zero-magnetic field, collisionless regime. The formulae are developed both in the standard representation and in the representation of the statistical tensors. For resonance scattering in a two-level atom, in the limiting case of complete depolarization of the ground level, we recover the classical results for Rayleigh scattering and we derive the expression of the phase matrix in terms of ordinary rotation matrices. The law of scattering is then generalized to take properly into account the influence of the anisotropy of the radiation field on the atomic polarization of the ground level (depopulation pumping).     

Stokes parameters for Thomson scattering in a strong magnetic field

The effect of a strong magnetic field on neutron stars or white dwarfs is calculated for Thomson scattering. The Stokes parameters for the scattered radiation are computed explicitly in terms of the state of polarization of the incident wave, the electron-cyclotron frequency, the angle of incidence, and the angle of scattering. The effect of the magnetic field on these parameters is considered in detail for incoming radiation propagating along or across the external magnetic field. The criterion for the magnetic field to substantially affect the Stokes parameters is that the photon frequency be less than the electron-cyclotron frequency.     

Order-of-scattering of partially polarized radiation in inhomogeneous,anisotropically scattering atmospheres

A complete set of transfer equations required for the order-of-scattering analysis of partially polarized radiation in inhomogeneous, anisotropically scattering atmospheres is provided. The equations have been derived for both a local study using the radiative transfer equation and its associated auxiliary equation for the source-matrix, and a global study in terms of the scattering and transmission matrices; they account for the polarity of the scattering medium. Their derivations for the finite order scattering and the finitely cumulative scattering, in particular, have yielded important new equations expressing the invariance principles and the integro-differential recurrences for the scattering and transmission matrices. These novel expressions contain as a special case Bellmanet al's (1972) equations for the simpler case of isotropic scattering of unpolarized light in homogeneous atmospheres.     

强磁场中的逆Compton散射与γ射线爆能谱

本文认为强磁场中的逆Compton散射可能是γ射线爆的主要辐射机制.其能谱是由源区质子产生的低频光子经强磁场中非热电子的Compton散射形成的.我们利用非相对论情形(B/B_(cr)≤1,hv_i/m_ec~2≤1)下强磁场中的Compton散射微分截面,导出了上述Compton散射的辐射谱公式,由此很好地拟合了典型γ射线爆GB811016的观测能谱.     

Elemental abundance analyses with DAO spectrograms — XVIII. The double-lined spectroscopic binary 46 Draconis

Equipartition magnetic fields can dramatically affect the polarization of radiation emerging from accretion disc atmospheres in active galactic nuclei. We extend our previous work on this subject by exploring the interaction between Faraday rotation and absorption opacity in local, plane-parallel atmospheres with parameters appropriate for accretion discs. Faraday rotation in pure scattering atmospheres acts to depolarize the radiation field by rotating the polarization planes of photons after last scattering. Absorption opacity in an unmagnetized atmosphere can increase or decrease the polarization compared to the pure scattering case, depending on the thermal source function gradient. Combining both Faraday rotation and absorption opacity, we find the following results. If absorption opacity is much larger than scattering opacity throughout the atmosphere, then Faraday rotation generally has only a small effect on the emerging polarization because of the small electron column density along a photon mean free path. However, if the absorption opacity is not too large and it acts alone to increase the polarization, then the effects of Faraday rotation can be enhanced over those in a pure scattering atmosphere. Finally, while Faraday rotation often depolarizes the radiation field, it can in some cases increase the polarization when the thermal source function does not rise too steeply with optical depth. We confirm the correctness of the analytic calculation by Silant'ev of the high magnetic field limit of the pure scattering atmosphere, which we incorrectly disputed in our previous paper.