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.     

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.     

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.     

Time-dependentX- andY-functions by integral equation method

The time-dependent equation of radiative transfer for isotropic scattering has been solved by integral equation technique in terms ofX- andY-functions appropriate for the problem. It is seen thatX- andY-functions are reducible to the corresponding function for steady-state problems by simply changing the Laplace transform parameters-i.e., byS0.     

Time-dependentH-,X- andY-functions for an anisotropically-scattering atmosphere

We have considered a homogeneous atmosphere scattering anisotropically with Dirac -function type time-dependent incidence. We used the method of integral operator developed by Ambartsumian and the theory ofN-solutions developed by Busbridge to find the correspondingH-function (in semi-infinite atmosphere) andX- andY-functions (in finite atmosphere).     

Free molecular flow of rarefied gas over an oscillating plate under a periodic external force

The free molecular flow over an infinite oscillating plane wall under external periodic force is considered. The Boltzmann equation is solved by using moments method with two-stream distribution functions. The boundary condition is obtained by assuming that the reflection of the particles from the solid surface takes place with complete energy accommodation. An analytical form for the velocity (X) and shear stress (Y) at any point is obtained. The results show that the amplitude of both the velocity change (X ₁) and the shear stress change (Y ₁) due to the periodic external force at the boundary (y=0) is an increasing function of time (t).     

5D space-time-mass gravity theory with off-diagonal components in the metric tensor

The bi-variational technique is used to calculated Chandrasekhar'sX- andY-functions and their first two moments for an isotropic homogeneous finite slab. Numerical results obtained are compared with published results.     

Maximum-entropy for chandrasekhar'sX- andY-functions

The maximum-entropy approach is used to calculate the Chandrasekhar'sX- andY-function and their moments. Numerical results are obtained and compared.     

An exact linearization and decoupling of the integral equations satisfied by time-dependent X- and Y-functions

We discuss a simple method of linearization and decoupling of the integral equations satisfied by time-dependentX - andY -functions which play an important rôle in the study of non-stationary radiative transfer problems.     

Time-dependentX- andY-functions for anisotropically scattering inhomogeneous atmosphere by principle of invariance

The nonlinear integral equations forX- andY-functions have been developed for an inhomogeneous atmosphere scattering anisotropically using the principle of invariance. The anisotropy is represented by means of a phase function expressed in terms of finite-order Legendre polynomials.     

Time-dependentX- andY-functions for an isotropically scattering neutron transport problem in a finite slab

Expressions for time-dependentX- andY-functions for a one-speed neutron transport problem in a finite slab have been derived using a technique combining invariant imbedding method and eigenfunction expansion method. The atmosphere has been considered to scatter isotropically.     

External radiation field in Rayleigh-Cabannes atmospheres with constant and linear sources

The external field of radiation in Rayleigh-Cabannes atmospheres with constant and linear sources is found using the resolvent matrix approach. If the internal sources are constant the external field may be described by theX-, Y-, andH-matrices. For the case with linear sources we need the derivatives of these matrices with respect to angular variable. The respective scheme for their determination is given.A set of integro-differential equations for theX- andY-matrices is derived and solved numerically. Some relations between the moments of theH-matrix are given and a sample of results for external fields are provided.     

Average photon path-length in isotropically scattering finite atmospheres

The determination of the average path-length of photons in a finite isotropically scattering plane-parallel homogeneous atmosphere is discussed. To solve this problem we have used the kernel approximation method which easily allows us to find the derivatives of the intensity with respect to optical depth, optical thickness and albedo of single scattering.In order to check the results we have used another approach by exploiting the set of integrodifferential equations of Chandrasekhar for theX- andY-functions. This approach allows us to find the average path length only at the boundaries of the atmosphere but on the other hand it gives also the dispersion of the path-length distribution function, thus generating the input parameters for determining the approximate path-length distribution function. It occurred that the set so obtained is stable and the results are highly accurate.As a by-product we obtain the first two derivatives of theX- andY-functions with respect to the albedo of single scattering and optical thickness, and the mixed derivative.     

Derivation of frequency-dependentX- andY-functions for a non-coherent scattering problem

The equation of transfer for the case of non-coherent scattering (Hummer, 1968; Ivanov, 1973; McCormick and Siewert, 1970) has been considered. The correspondingX- andY-functions have been derived by a combination of eigenfunction method developed by Case, and from the principle of invariance as developed by Chandrasekhar (1960).     

Characteristics of plasma temperature variation during two-current-loop collisions

Characteristics of plasma temperature variations during two-current-loop collisions are described. It is revealed that plasma temperature has an oscillatory feature with damping amplitude and growing quasi-period in the case of anI-type collision. In the case of aY-type collision, if the initial current becomes strong enough, there also occur pulsations of the temperature. However, the temperature profile of anX-type collision is characterized by a single pulse only.     

Evolutionary considerations involving the internal density concentration parameter of binary stars

The evolutionary changes that occur in the internal density concentration parameterk ₂ (called the apsidal constant for brevity) for a star of given mass and initial composition are examined in detail. The purpose is to ascertain whether or not such an approach leads to a reduction in the differences now noted between the theoretically derived values ofk ₂ and the observed values derived from the secular advance of the periastron in close eclipsing binary systems.A series of stellar models of mass 2.0, 5.0, 10.0 and 20.0M were employed, with an initial compositional mixture ofX=0.739,Y=0.24 andZ=0.021. These models cover an evolutionary range from a point in time where the star has just completed the Hayashi phase of its pre-Main Sequence contraction through its entire Main Sequence phase to the point where hydrogen depletion in the core is complete.For each model, a value ofk ₂ is determined by numerically integrating Radau's equation and using the values of , the ratio of the star's density at pointa to its mean density, as taken from the models. The result is the time history ofk ₂ for each stellar mass over the evolutionary range of interest. The results are then summarized in the (logk ₂, logT _e) plane which, for the first time, quantitatively indicates the variation ink ₂ as a function of the evolutionary state of the star.A comparison between these theoretically derived values ofk ₂ and a selected set of observationally determined ones immediately indicates that the secular variation ink ₂ plays an extremely important part in any comparison between theory and observation. For most of the cases studied, the difference between the theoretically and observationally determined values ofk ₂ can be reconciled in terms of the evolutionary history of the binary system.While tentatively providing a satisfactory explanation for the previously noted differences in the determination ofk ₂, there now exists the problem of accurately pinpointing the evolutionary state of the observed binary system.     

Pulsations of stellar models in H and He burning phases

A study of pulsational properties with evolution has been done for a 15.6M star withX _e =0.90 andY _e =0.08. Pulsational properties in the hydrogen-burning stages have been compared with those in helium-burning stages. A comparison with observed characteristics of Cepheids, classical Cepheids and supergiant variables has been made during the course of its evolution. In addition, models of 5,9 and 15M withX _e =0.708,Y _e =0.272 have also been studied for pulsational properties during the helium burning stage. It is also seen that pulsational instability is sensitive to changes in initial chemical composition and opacity parameters,n ands. A low helium abundance could be a reason for the stability of the models, even when lying in the instability strip of the H-R diagram.     

Bivariate least squares linear regression: Towards a unified analytic formalism. I. Functional models

Concerning bivariate least squares linear regression, the classical approach pursued for functional models in earlier attempts (York, 1966, York, 1969) is reviewed using a new formalism in terms of deviation (matrix) traces which, for unweighted data, reduce to usual quantities leaving aside an unessential (but dimensional) multiplicative factor. Within the framework of classical error models, the dependent variable relates to the independent variable according to the usual additive model. The classes of linear models considered are regression lines in the general case of correlated errors in X and in Y for weighted data, and in the opposite limiting situations of (i) uncorrelated errors in X and in Y, and (ii) completely correlated errors in X and in Y. The special case of (C) generalized orthogonal regression is considered in detail together with well known subcases, namely: (Y) errors in X negligible (ideally null) with respect to errors in Y; (X) errors in Y negligible (ideally null) with respect to errors in X; (O) genuine orthogonal regression; (R) reduced major-axis regression. In the limit of unweighted data, the results determined for functional models are compared with their counterparts related to extreme structural models i.e. the instrumental scatter is negligible (ideally null) with respect to the intrinsic scatter (Isobe et al., 1990, Feigelson and Babu, 1992). While regression line slope and intercept estimators for functional and structural models necessarily coincide, the contrary holds for related variance estimators even if the residuals obey a Gaussian distribution, with the exception of Y models. An example of astronomical application is considered, concerning the [O/H]–[Fe/H] empirical relations deduced from five samples related to different stars and/or different methods of oxygen abundance determination. For selected samples and assigned methods, different regression models yield consistent results within the errors (?σ) for both heteroscedastic and homoscedastic data. Conversely, samples related to different methods produce discrepant results, due to the presence of (still undetected) systematic errors, which implies no definitive statement can be made at present. A comparison is also made between different expressions of regression line slope and intercept variance estimators, where fractional discrepancies are found to be not exceeding a few percent, which grows up to about 20% in the presence of large dispersion data. An extension of the formalism to structural models is left to a forthcoming paper.     

Rayleigh scattering: Volterra equation for the matrix source function

Multiple Rayleigh scattering is examined in a semi-infinite atmosphere with uniformly distributed primary sources of partially polarized radiation. The resulting linear polarization is described by a 2×2 matrix transfer equation. A matrix generalization of Rybicki's two point Q-integral is obtained for this case. It is shown that the Volterra equation for the matrix source function for this problem is a particular case of our Q integral. Applying the Laplace transform to it yields the matrix form of the Ambartsumyan-Chandrasekhar H-equation. The Volterra equation for Sobolev's matrix resolvent function is another simple consequence of this equation. Translated from Astrofizika, Vol. 52, No. 2, pp. 301–310 (May 2009).