Exact solution of the problem of diffuse reflection and transmission by the method of Laplace transform and linear singular operators

The basic integro-differential equation is subjected to a one-sided finite Laplace transform to obtain linear integral equations of angular distribution of bounding faces. These linear integral equations have been transformed into linear singular integral equations which have been solved exactly to get the emergent distributions from the bounding faces by the theory of linear singular operators. Some solutions of linear singular integral equations have also been derived for future use in radiative transfer problems.     

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.     

The transfer of polarized radiation in spectral lines: Formalism and solutions in simple cases

A general treatment of the transfer of polarized radiation in spectral lines assuming a Rayleigh phase function and a general law of frequency redistribution is derived. It is shown how nine families of coupled integral equations for the moments of the radiation field arise which are necessary to fully describe the state of polarization of the emergent radiation from a plane-parallel, semi-infinite atmosphere. The special case of angle independent redistribution functions is derived from the general formalism, and it is shown how the nine families of integral equations reduce to the six linearly independent integral equations derived by Collins (1972). To serve as a test of the formulation, solutions for isothermal atmospheres are given.     

Solution of Riemann–Hilbert problems for determination of new decoupled expressions of Chandrasekhar’s <Emphasis Type="Italic">X</Emphasis>- and <Emphasis Type="Italic">Y</Emphasis>-functions for slab geometry in radiative transfer

In radiative transfer, the intensities of radiation from the bounding faces of a scattering atmosphere of finite optical thickness can be expressed in terms of Chandrasekhar’s X- and Y-functions. The nonlinear nonhomogeneous coupled integral equations which the X- and Y-functions satisfy in the real plane are meromorphically extended to the complex plane to frame linear nonhomogeneous coupled singular integral equations. These singular integral equations are then transformed into nonhomogeneous Riemann–Hilbert problems using Plemelj’s formulae. Solutions of those Riemann–Hilbert problems are obtained using the theory of linear singular integral equations. New forms of linear nonhomogeneous decoupled expressions are derived for X- and Y-functions in the complex plane and real plane. Solutions of these two expressions are obtained in terms of one known N-function and two new unknown functions N ₁- and N ₂- in the complex plane for both nonconservative and conservative cases. The N ₁- and N ₂-functions are expressed in terms of the known N-function using the theory of contour integration. The unknown constants are derived from the solutions of Fredholm integral equations of the second kind uniquely using the new linear decoupled constraints. The expressions for the H-function for a semi-infinite atmosphere are obtained as a limiting case.     

Linear fREDHOLM integral equations for radiative transfer problems in finite plane parallel media. I. Imbedding in an infinite medium

Linear FREDHOLM integral equations are derived for the STOKES vector of the radiation emerging from a scattering plane parallel medium of finite optical thickness. The integral equations are obtained by means of imbedding the slab in an infinite medium. They are formulated in terms of GREEN 's function matrices and renormalized for the asymptotic eigenmode. Explicitly, linear integral equations are given for the reflection and transmission matrices. The reciprocity principle is employed to obtain integral equations also for the mean intensity of the inner radiation field in the case of the slab albedo problem.     

Exact solution of a basic equation in finite atmosphere by the method of Laplace transform and linear singular operators

A finite atmosphere having distribution of intensity at both surfaces with definite form of scattering function and source function is considered here. The basic integro-differential equation for the intensity distribution at any optical depth is subjected to the finite Laplace transform to have linear integral equations for the surface quantities under interest. These linear integral equations are transformed into linear singular integral equations by use of the Plemelj's formulae. The solution of these linear singular integral equations are obtained in terms of theX-Y equations of Chandrasekhar by use of the theory of linear singular operators which is applied in Das (1978a).     

Exact solution of a simple time-dependent integro-differential equation by the method of Laplace transform and the theory of linear singular operators

The simplest form of the equation of transfer for a time dependent radiation field in finite atmosphere is considered. This equation of transfer is an integro-differential equation, the solution of this equation is based on the theory of separation of variables, the Laplace transform and the theory of linear singular operators. The emergent intensities from the bounding faces of the finite atmosphere are determined in terms ofX-Y equations of Chandrasekhar.     

Solution of non-coherent spectral line transfer problem for finite atmosphere by theF n -method

By use of the orthogonality and normalization integrals developed by McCormick and Siewert (1970) a set of singular integral equations suitable forF _n -method is derived for non-coherent spectral line formation problem in finite media.F _n -equations for exit distributions are used to develop some algebraic equations with suitable recursion relations.     

Exact solution of the problem of diffuse reflection for a combination of Rayleigh and isotropic scattering

We consider the basic vector equation of transfer for radiation in a semi-infinite atmosphere for diffuse reflection which scatters radiation in accordance with the phase matrix obtained from a combination of Rayleight and isotropic scattering. This equation will give an integral equation for emergent intensity while subjected to the Laplace transform. The integral equation will give rise to the emergent intensity matrix on application of the Wiener-Hopf technique. This is an exact method.     

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.     

The solution of resonance radiative transfer equation for diffusion by the method of laplace transform and linear singular operators

The equation for radiative transfer in the case of resonance radiation for isotropic scattering has been solved by the method of the Laplace transformation and linear singular operators. The solution for emergent intensities have come out in terms ofX- andY-functions.     

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.     

Solution of a Riemann–Hilbert problem for a new expression of Chandrasekhar’s H-function in radiative transfer

The linear singular integral equation derived from the nonlinear integral equation of Chandrasekhar’s H-function in radiative transfer is considered here to develop a new form of H-function as a solution of a Riemann–Hilbert problem using Plemelj and Cauchy integral formulae for complex domain. This new form of H-function is a simple integral of known functions. Forms of H-function both for conservative and nonconservative cases are obtained. Their numerical evaluations are made by Simpson’s one-third rule to arrive at an accuracy to ninth places of decimals.     

Exact solution of a key equation in a finite planetary atmosphere by the method of laplace transform and linear singular operators

Considering the ground reflection according to Lambert's law, we establish a fundamental equation in finite planetary atmospheres. An exact form of the solution of this equation is obtained for the emergent quantities from the bounding faces in terms ofX-Y equations by the method of Laplace transform, in combination with the theory of linear singular operators.     

A method of computing the emergent radiation by the atmosphere in the region ranging from ultraviolet to infrared

A method of computing the diffuse reflection and transmission radiation by an inhomogeneous, plane-parallel planetary atmosphere with internal emission source is discussed by use of the adding method. If the atmosphere is simulated by a number of homogeneous sub-layers, the radiation diffusely reflected or transmitted by the atmosphere can be expressed in terms of the reflection and transmission matrices of the radiation of sub-layers. The diffusely transmitted radiation due to the internal emission source can be also easily computed in the same manner. These equations for the emergent radiation are in a quite general form and are applicable to radiative transfer in the atmosphere in the region from ultraviolet to infrared radiation. With this method, the tiresome treatment due to the polarity effect of radiation is overcome.     

Exact solution of the equation of radiative transfer for semi-infinite plane-parallel isotropic scattering atmosphere with scattering albedo ω0 ≤ 1 by a method based on Laplace transform and Wiener-Hopf technique

By performing the one-sided Laplace transform on the scalar integro-differential equation for a semi-infinite plane-parallel isotropic scattering atmosphere with a scattering albedo ₀ 1, an integral equation for the emergent intensity has been derived. Application of the Wiener-Hopf technique to this integral equation will give the emergent intensity. The intensity at any optical depth for a positive scattering angle is also derived by inversion. The intensity at any optical depth for a negative scattering angle is also derived in terms of Cauchy's principal value using Plemelj's formulae.     

Atmospheres and spectra of strongly magnetized neutron stars

We construct atmosphere models for strongly magnetized neutron stars with surface fields     and effective temperatures     . The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars, including radio pulsars, soft gamma-ray repeaters, and anomalous X-ray pulsars. In our models, the atmosphere is composed of pure hydrogen or helium and is assumed to be fully ionized. The radiative opacities include free–free absorption and scattering by both electrons and ions computed for the two photon polarization modes in the magnetized electron–ion plasma. Since the radiation emerges from deep layers in the atmosphere with     , plasma effects can significantly modify the photon opacities by changing the properties of the polarization modes. In the case where the magnetic field and the surface normal are parallel, we solve the full, angle-dependent, coupled radiative transfer equations for both polarization modes. We also construct atmosphere models for general field orientations based on the diffusion approximation of the transport equations and compare the results with models based on full radiative transport. In general, the emergent thermal radiation exhibits significant deviation from blackbody, with harder spectra at high energies. The spectra also show a broad feature     around the ion cyclotron resonance     , where Z and A are the atomic charge and atomic mass of the ion, respectively; this feature is particularly pronounced when     . Detection of the resonance feature would provide a direct measurement of the surface magnetic fields on magnetars.     

A method of evaluation of the atmospheric effect on the emergent radiation over the infinite-cliff

To evaluate the effect of the cliff on the radiation field, the upwelling radiation at the top of the atmosphere is computed over the cliff using the reflection and transmission functions derived from the doubling-adding method. The model is defined by the plane-parallel homogeneous atmosphere, which is composed of aerosol and molecules, and is bounded by the top level surface, cliff and low level surface. These surfaces may be assumed to be the Lambertian.In the computational procedure, the equation for the emergent radiation is expanded into a series of radiative interaction modes among atmosphere, surfaces and the cliff. In respective modes, probabilities of respective interactions are firstly evaluated. With the aid of these probabilities, the emergent radiation is calculated using the doubling-adding method for the model atmosphere bounded by the surfaces and cliff, where the above radiative interactions are considered upto twice as large to obtain the enough accuracy of simulation. The multiple scattering is considered.     

The modified Lagrangian equations of general-relativistic hydrodynamics for a nonideal fluid

The general-relativistic equations of hydrodynamics for a nonideal fluid are derived in the modified Lagrangian form. Together with the zeroth and first moment equations of radiative transfer derived by Morita and Kaneko (1986), the equations provide a complete set of the modified Lagrangian equations of radiation hydrodynamics. The equations of hydrodynamics are specialized for a thermally conducting, Newtonian viscous fluid in the modified Lagrangian form, which are the generalization of the specialrelativistic equations of hydrodynamics derived by Greenberg (1975).     

An exact solution of transfer equations for interlocked multiplets

The transfer equations for non-coherent scattering arising from interlocking of principal lines without redistribution is exactly solved in a very simple way by Laplace tranform and Wiener-Hopf technique which are easily applied by the use of the new representation ofH-functions obtained recently by the author (1977). The emergent intensity in therth line is expressed in terms of anH-function and a Cauchy type integral admitting of closed form evaluation.