Radiative transfer in spectral lines by noncoherent scattering

A functional analytic method of solution of the operator equations is used to obtain a solution of the radiative transfer equation in spectral lines. A problem of scattering in the spectral line with frequency redistribution in an isotropic scattering medium is considered.     

The highly excited hydrogen b n -factors calculated

Radio recombination line intensities are functions of highly excited atom nonequilibrium populations. These populations are determined by the spontaneous, collision and induced radiation transition rates that arise from the balance equations. The population balance equations are solved in a diffusion approximation which is proved for the departure factors from LTE. Diffusion equation departure coefficients are far from LTE due to the spontaneous transition flow, found as a second-order differential equation general solution. This solution is chosen in an explicit form from a physical situation with the temperature and density determined. The background radiation induced transition rates strongly influence the population flow to the equilibrium distribution. Radio recombination line amplification coefficients are calculated for the different medium temperatures and electron densities; these coefficients permit interpretation of the line intensities in the actual experiments.     

The three-dimensional motion of Trojan asteroids

The problem is considered within the framework of the elliptic restricted three-body problem. The asymptotic solution is derived by a three-variable expansion procedure. The variables of the expansion represent three time-scales of the asteroids: the revolution around the Sun, the libration around the triangular Lagrangian pointsL ₄,L ₅, and the motion of the perihelion. The solution is obtained completely in the first order and partly in the second order. The results are given in explicit form for the coordinates as functions of the true anomaly of Jupiter. As an example for the perturbations of the orbital elements the main perturbations of the eccentricity, the perihelion longitude and the longitude of the ascending node are given. Conditions for the libration of the perihelion are also discussed.     

Volterra equation for the matrix source function at resonance scattering

A matrix transfer equation for multiple resonance scattering of radiation in a spectral line in a semiinfinite atmosphere with a uniform distribution of primary radiation sources is examined. A nonlinear matrix integral is obtained for this equation as a generalization of the Rybicki two-point Q-integral. One special case of the matrix [^(Q)] {\mathbf{\hat{Q}}} -integral is the Volterra equation for the matrix source function of the problem discussed here. The Volterra equation is solved numerically for a Doppler profile of the absorption coefficient. Several polarization characteristics of the emerging radiation are obtained.     

On nonstationary radiative transfer in a spectral line in stellar atmospheres

We consider nonstationary radiative transfer in a line in stellar atmospheres simulated as a stationary semi-infinite plane-parallel medium. We assume complete frequency redistribution in the elementary act of scattering. We assume that the time a photon spends in the medium is determined only by the mean time spent in the absorbed state. We obtain an explicit expression for the resolvent of the nonstationary integral equation of transfer, which is a bilinear expansion with respect to the eigenfunctions found in [12] for the corresponding stationary transfer equation.Translated fromAstrofizika, Vol. 37, No. 3, 1994.I am grateful to the American Astronomical Society for financial support of this work.     

An algorithm for solving the pulsar equation

We present an algorithm of finding numerical solutions of pulsar equation. The problem of finding the solutions was reduced to finding expansion coefficients of the source term of the equation in a base of orthogonal functions defined on the unit interval by minimizing a multi-variable mismatch function defined on the light cylinder. We applied the algorithm to Scharlemann and Wagoner boundary conditions by which a smooth solution is reconstructed that by construction passes successfully the Gruzinov’s test of the source function exponent.      

Calibration of Vector Magnetograms with the Chromospheric Mg b2 Line

Stokes polarization profiles of the Mg?b₂ 5172.68 Å spectral line on two simple sunspots are obtained with the Multi-Channel Solar Telescope (MCST) at the Huairou Solar Observing Station (HSOS). This is done by means of scanning this line over the wavelength interval from 200 mÅ redward of the line center to 200 mÅ blueward, in steps of 10 mÅ. A generalized analytic solution to the transfer equation for polarized radiation is presented. With a nonlinear least-square fitting technique, the linear calibration coefficients for the low-chromospheric longitudinal magnetic field is obtained in the weak-field case. We also discuss the problems in calibrating the transverse field with this line. It is shown that the weak-field approximation is not applicable to the chromospheric Mg?b₂ line for the transverse component of the magnetic field.     

Laser radiation in active amplifying media treated as a transport problem: Transfer equation derived and exactly solved

The flow of laser radiation in a plane-parallel cylindrical slab of active amplifying medium with axial symmetry is treated as a problem in radiative transfer. The appropriate one-dimensional transfer equation describing the transfer of laser radiation has been derived by an appeal to Einstein'sA, B coefficients (describing the processes of stimulated line absorption, spontaneous line emission, and stimulated line emission sustained by population inversion in the medium) and considering the rate equations to completely establish the rational of the transfer equation obtained. The equation is then exactly solved and the angular distribution of the emergent laser beam intensity is obtained; its numerically computed values are given in tables and plotted in graphs showing the nature of peaks of the emerging laser beam intensity about the axis of the laser cylinder.     

The elimination of short and intermediate period terms from the problem of a high altitude earth satellite

The Delaunay-Similar elements of Scheifele are applied to the problem of an Earth satellite that is perturbed by the Sun, Moon andJ ₂. All three effects are assumed to be the same order of magnitude. Since the external body terms depend explicitly on time, the time element appears as an additional angle variable. Also, the eccentric anomaly is used as a noncanonical auxiliary variable. A solution to the first Von Zeipel equation allows simultaneous elimination of short and intermediate period terms. The canonical transformation to mean elements is defined by a generating function that is a series involving Bessel coefficients.     

Radiative transfer in a cylinder. II. Special problems. Asymptotics

The method of computing the radiation field in an infinite circular cylinder proposed in Part I is now applied to the case of isotropic scattering with sources on the boundary and axis of the cylinder, as well as for a uniform distribution of sources inside the cylinder. For the simplest kernel we obtain exact solutions of the basic integral equation in explicit form. For scattering in a spectral line with complete frequency redistribution and a power absorption profile we develop an asymptotic theory for the case when the optical radius of the cylinder is large. We solve the asymptotic equations for the basic characteristics of the scattering in closed form for conservative scattering and find its asymptotics. We obtain estimates of the mean number of scatterings with a layered source, and also the mean and variance of the number of scatterings with a uniform source distribution.Translated fromAstrofizika, Vol. 37, No. 4, 1994.This work was carried out with the financial support of the Russian Basic Research Fund (grant 93-02-2957).     

A new analytic approach to the Sitnikov problem

A new analytic approach to the solution of the Sitnikov Problem is introduced. It is valid for bounded small amplitude solutions (z _max = 0.20) (in dimensionless variables) and eccentricities of the primary bodies in the interval (–0.4 < e < 0.4). First solutions are searched for the limiting case of very small amplitudes for which it is possible to linearize the problem. The solution for this linear equation with a time dependent periodic coefficient is written up to the third order in the primaries eccentricity. After that the lowest order nonlinear amplitude contribution (being of order z ³) is dealt with as perturbation to the linear solution. We first introduce a transformation which reduces the linear part to a harmonic oscillator type equation. Then two near integrals for the nonlinear problem are derived in action angle notation and an analytic expression for the solution z(t) is derived from them. The so found analytic solution is compared to results obtained from numeric integration of the exact equation of motion and is found to be in very good agreement. CERN SL/AP     

Non-Abelian brane cosmology

We discuss isotropic and homogeneous D-brane-world cosmology with non-Abelian Born-Infeld (NBI) matter on the brane. In the usual Friedmann-Robertson-Walker (FRW) model the scale non-invariant NBI matter gives rise to an equation of state which asymptotes to the string gas equation p=-ε/3 and ensures a start-up of the cosmological expansion with zero acceleration. We show that the same state equation in the brane-world setup leads to the Tolman type evolution as if the conformal symmetry was effectively restored. This is not precisely so in the NBI model with symmetrized trace, but the leading term in the expansion law is still the same. A cosmological sphaleron solution on the D-brane is presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Unified and Very Fast Way for Computing the Global Potential and Linear Force-Free Fields

We present a fast solver for computing potential and linear force-free fields (LFFF) above the full solar disk with a synoptic magnetic map as input. The global potential field and the LFFF are dealt with in a unified way by solving a three-dimensional Helmholtz equation in a spherical shell and a two-dimensional Poisson equation on the solar surface. The solver is based on a combination of the spectral method and the finite-difference scheme. In the longitudinal direction the equation is transformed into the Fourier spectral space, and the resulting two-dimensional equations in the r?C?? plane for the Fourier coefficients are solved by finite differencing. The solver shows an extremely fast computing speed, e.g., the computation for a magnetogram with a resolution of 180(??)×360(?) is completed in less than 2 s. Even on a high-resolution 600×1200 grid, the solution can be obtained within only about one minute on a single CPU. The solver can potentially be applied directly to the original resolution of observed magnetograms from SDO/HMI for routinely analyzing daily full-disk data.     

A nonlinear stability problem in the three dimensional restricted three body problem

The question of whether or not there is a transfer of energy between the in-plane motion and out-of-plane motion in the neighborhood ofL ₄ in the restricted problem of three bodies is investigated in this paper. The in-plane motion is assumed to be finite and the out-of-plane motion to be infinitesimal. The equation governing the out-of-plane motion becomes one with time varying coefficients. The stability of this equation is then investigated using Lie Series.Presented as a paper AAS No. 70-313, at the AAS/AIAA Astrodynamics Specialists Conference 1971 at Fort Lauderdale Fla., U.S.A.     

Empirical NLTE analyses of solar spectral lines

A method is suggested for empirical NLTE analyses of solar spectral lines. Special depthdependent departure coefficients (x) are introduced and the formulas are given for further application. From test calculations it is shown that the separation of the departure coefficients for the upper and the lower level from each other and from the uncertainties in several input parameters (oscillator strength times abundance, turbulent velocities and damping constant) is greatly facilitated when spectral lines are analysed on the disk, around the limb, as well as in the flash spectrum. Therefore it is necesssary that all accessible line data are used, from accurate line profiles to equivalent widths or integrated intensities. To reduce the number of independent variables the analysis should include many multiplets between a few spectroscopic terms.We discuss the advantages of translating the departure coefficients into a set of temperatures. Some published LTE and NLTE analyses are compared in terms of these temperatures. The results do not show inexplicable contradictions. For several atomic spectra the departures from LTE seem large enough to be established quantitatively by the detailed method suggested in this paper.     

An extended canonical perturbation method

In this investigation, a procedure is described for extending the application of canonical perturbation theories, which have been applied previously to the study of conservative systems only, to the study of non-conservative dynamical systems. The extension is obtained by imbedding then-dimensional non-conservative motion in a 2n-dimensional space can always be specified in canonical form, and, consequently, the motion can be studied by direct application of any canonical perturbation method. The disadvantage of determining a solution to the 2n-dimensional problem instead of the originaln-dimensional problem is minimized if the canonical transformation theory is used to develop the perturbation solution. As examples to illustrate the application of the method, Duffing's equation, the equation for a linear oscillator with cubic damping and the van der Pol equation are solved using the Lie-Hori perturbation algorithm.This research was supported by the Office of Naval Research under Contract N00014-67-a-0126-0013.     

On a class of variational equations transformable to the Gauss hypergeometric equation

A new class of linear ordinary differential equations with periodic coefficients is found which can be transformed to the Gauss hypergeometric equation, and therefore the monodromy matrices are computable explicitly. These equations appear as the variational equations around a straight-line solution in Hamiltonian systems of the form H = T(p) + V(q), where T(p) and V(q) are homogeneous functions of p and q, respectively.     

An approximate integral and solution for eccentric planetary type orbits in the restricted problem of three bodies

The circular restricted problem of three bodies is investigated analytically with respect to the problem of deriving a second integral of motion besides the well known Jacobian Integral. The second integral is searched for as a correction the angular momentum integral valid in the two body case. A partial differential equation equivalent to the problem is derived and solved approximately by an asymptotic Fourier method assuming either sufficiently small values for the dimensionless mass parameter or sufficiently large distances from the barycentre. The solution of the partial equation then leads to a function of the coordinates, velocities and time being nearly constant, which means that its variation with time is about 40–300 times less than that of the pure angular momentum. By averaging over the remaining fluctuating part of the quasi-integral we are able to integrate the first order equations using a renormalization transformation. This leads to an explicit expression for the approximate solution of the circular problem which describes the motion of the third body orbiting both primaries with nonvanishing initial eccentricity (eccentric planetary type orbits). One of the main results is an explicit formula for the frequency of the perihelion motion of the third body which depends on the mass parameter, the initial distance of the third body from the barycentre and the initial eccentricity. Finally we study orbits of the P-Type, being defined as solutions of the restricted problem with circular initial conditions (vanishing initial eccentricity).     

On the solving of the Hill's differential equation through the method of operational calculus

Using a method previously applied to the treatment of the Mathieu differential equation, we solve the Hill's differential equation of lunar theory through the way of operational calculus, which avoids the cumbersome infinite determinants of the classical procedure. The one-sided Laplace transformation changes it into a difference equation with an infinite number of terms and variable coefficients. When its first member is divided by a suitable factor, this difference equation is the image of an integral equation of the Volterra type which is equivalent to the initial Hill's differential equation. Solution of this Volterra integral equation is unique and it is the general solution of the Hill's differential equation. This solution is a series in the powers of a small dimensionless parameter ² which appears as a factor in the second member of the Hill's differential equation. We reduce it to the sum of its terms of degree 12 with respect to which is the precision usually required in a lunar theory and we write down effectively the coefficients of the terms in ², ( ²)² and the coefficient of the term in ( ²)³ which depends upon the initial valuey(0) of the Hill's differential equation.     

K/S two-point-boundary-value problems

A method for developing the missing general K/S (Kustaanheimo/Stiefel) boundary conditions is presented, with use of the formalism of optimal control theory. As an illustrative example, the method is applied to the K/S Lambert problem to derive the missing terminal condition. The necessary equations are developed for a solution to this problem with the generalized eccentric anomaly,E, as the independent variable. This formulation, requiring the solution of only one nonlinear, well-behaved equation in one unknown,E, results in considerable simplification of the problem.Presented also at the AAS/AIAA Astrodynamics Specialist Conference, Nassau, Bahamas, July 1975 (AAS Paper No. 75-032).