A unified presentation of the Voigt functions

This paper aims at presenting a unified study of the Voigt functionsK(x,y) andL(x,y) which play a rather important role in several diverse fields of physics such as astrophysical spectroscopy and the theory of neutron reactions. Explicit expressions for these functions are given in terms of relatively more familiar special functions of one and two variables; indeed, each of these representations will naturally lead to various other needed properties of the Voigt functions.     

Astrophysical spectroscopy and neutron reactions: Integral transforms and Voigt functions

The Voigt functionsK(x, y) andL(x, y) which play an essential role in astrophysical spectroscopy and neutron physics are investigated and generalized from the viewpoint of integral operators. Unified representations and series expansions involving classical functions of mathematical physics and multivariable hypergeometric functions are established. From the delicate asymptotic analysis of Laplace and Hankel integral transforms we extract complete and rigorous asymptotic expansions of the generalized Voigt functions for large values of the variablesx andy which are of great value in the theory of spectral line profiles.     

Some unified presentations of the Voigt functions

The principal object of this note is to provide a natural further step toward the unified presentations of the Voigh functionsK(x, y) andL(x, y) which play a rather important role in such diverse fields of physics as astrophysical spectroscopy and the theory of neutron reactions. Explicit representations for these functions, given in terms of some relatively more familiar special functions of one and two variables, are potentially useful in finding many other needed (numerical or analytical) properties of the Voigt functions. Several erroneous recent contributions to the theory of Voigt functions, including (for example) the main result of A. Siddiqui (1990), are also corrected here.     

Inverse problem with two-parametric families of planar orbits

As a possible extension of recent work we study the following version of the inverse problem in dynamics: Given a two-parametric familyf(x, y, b)=c of plane curves, find an autonomous dynamical system for which these curves are orbits.We derive a new linear partial differential equation of the first order for the force componentsX(x, y) andY(x, y) corresponding to the given family. With the aid of this equation we find that, depending on the given functionf, the problem may or may not have a solution. Based on given criteria, we present a full classification of the various cases which may arise.     

Spectral line profiles and neutron cross sections: New results concerning the analysis of Voigt functions

The Voigt functions, so important in spectroscopy and neutron physics, are represented as generalized hypergeometric functions (G-functions) of two real variables. A system of partial differential equations for the Voigt functions is derived. By applying Hölder's inequality to an integral representation of the Voigt functions apparently not known in the literature until now, lower and upper bounds are obtained. Moreover, from this representation an asymptotic expansion of Voigt functions for large values of one variable is extracted.     

Solution of the three-dimensional inverse problem

The three dimensional inverse problem for a material point of unit mass, moving in an autonomous conservative field, is solved. Given a two-parametric family of space curvesf(x, y, z)=c ₁,g(x, y, z)=c ₂, it is shown that, in general, no potentialU=U(x, y, z) exists which can give rise to this family. However, if the given functionsf(x, y, z) andg(x, y, z) satisfy certain conditions, the corresponding potentialU(x, y, z), as well as the total energyE=E(f, g) are determined uniquely, apart from a multiplicative and an additive constant.     

The gravitational field of a disk

This note gives the gravitational potential of the disk {(x, y, z):x ² +y ² p ² , z=0} and the gravitational field at the point (x, y, z). Formulas for a ring can be obtained as the difference of our results for two different values ofp. Results are obtained in terms of elliptic integrals and we indicate how these functions can be computed efficiently. Formulas necessary for the computation of partial derivatives are also given.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract NAS7-100, sponsored by the National Aeronautics and Space Administration.     

Generalization of Szebehely's equation

Szebehely's partial differential equation for the force functionU=U(x,y) which gives rise to a given family of planar orbitsf(x,y)=Constant is generalized to account for velocity-dependent potentials V^*=V^*(x,y, ). The new partial differential equation is quasi-linear and of the first order. An example is given and a comparison is made of the two equations.     

The angle-dependent redistribution functionsR III andR IV: Line absorption probability distribution functions and reemission coefficients

The line absorption probability distribution functions and the reemission coefficients are derived for the non-coherent scattering functionsR _III andR _IV. The appropriate line profile function forR _III is shown to be a simple Voigt function, while forR _IV, the line absorption probability distribution function is more complex involving a linear combination of two Voigt functions and another more complex probability distribution. The structure of the reemission coefficients forR _III andR _IV is then discussed.     

An equation for the characteristic curve of a family of symmetric periodic orbits

Szebehely's equation for the inverse problem of Dynamics is used to obtain the equation of the characteristic curve of a familyf(x,y)=c of planar periodic orbits (crossing perpendicularly thex-axis) created by a certain potentialV(x,y). Analytic expressions for the characteristic curves are found both in sideral and synodic systems. Examples are offered for both cases. It is shown also that from a given characteristic curve, associated with a given potential, one can obtain an analytic expression for the slope of the orbit at any point.     

Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

Compatibility conditions for a non-quadratic integral of motion

Conditions are found which are satisfied by the coefficients of the expression being a second integral of the motion of an autonomous dynamical system with two degrees of freedom. The coefficientsA, B. , ,E are differentiable functions of the cartesian position coordinatesx, y. The velocity components are denoted by . It is shown that must be constant andB must be of the formB =f(x+y) +g(x-y) wheref, g are arbitrary.Given andB one can always find the remaining coefficientsA, E and also the corresponding potential and second integral. Depending on the specifica case at hand a certain number of arbitrary constants (or arbitrary functions) enter into the potential and the second integral. To each potential (which may be of the separable or nonseparable type in the coordinatesx andy)there corresponds one integral of the above form.     

Intrinsic stability of periodic orbits

Families of orbits of a conservative, two degree-of-freedom system are represented by an unsteady velocity field with componentsu(x, y, t) andv(x, y, t). Intrinsic stability properties depend on velocity field divergence and curl, whose dynamical evolution is determined by a matrix Riccati equation. Near equilibrium, divergence-free or irrotational fields are dynamically compatible with the conservative force field. It is shown that a necessary condition for stable periodic orbits is satisfied when the orbitaveraged divergence is zero, which results in bounded normal variations. A sufficient condition for stability is derived from the requirement that tangential variations do not exhibit secular growth.In a steady, divergence-free field, velocity component functionsu(x, y) andv(x, y) may be continuedanalytically from any initial condition, except when velocity is parallel to U or at equilibria. In an unsteady field, the orbit-averaged divergence is zero when the vorticity function is periodic. When such a field exists, initial conditions for stable periodic orbits (i.e., characteristic loci) may be determinedanalytically.     

A systematic method for the analysis of high-resolution fraunhofer line profiles

A fraunhofer line profile depends on various parameters, partly related to the photospheric structure (T, P _g, P _e, v _conv, v _turb), partly to the atom or ion involved (such as oscillator strength, energy levels), partly also resulting from the interaction of the relevant kind of particles with the photosphere, and the photospheric radiation field. In this paper we shall mainly pay attention to the determination of: the macroturbulent (convective) velocities, v _conv (); the damping constant (); the abundance, A _el; the distribution function (v _conv, ) of the convective velocities at each depth ; the source function, S (); the microturbulent velocities, v _turb ().The particular difficulty with these unknowns is that they are, as a rule, coupled in the resulting line profiles, that is: the shapes and intensities in these profiles are determined by the combined influence of these unknowns (together with the other above-given parameters).In this paper we describe a method to determine these six unknowns empirically by separating them, in analysing accurate high-resolution observations of line profiles of a multiplet. The unknown functions and quantities are consecutively determined in the above given succession. For each determination another, appropriate part of the line profile is used. In some cases the influence of the mutual coupling of the various parameters cannot be completely eliminated, and an iterative method has to be used.The method is summarized in Table II and section 2, and is further explained in sections 3 to 8. It is applied to an infrared Ci multiplet. The main results are the following:     

The solution of radiation transfer problem for a slab with space-dependent albedo and anisotropic scattering

The transport of thermal radiation has been considered within a finite slab which absorb and scatter anisotropically. The problem involves the space-dependent single-scattering albedow(x). Two approximations are taken forw(x). In the first it is represented in exponential form asw(x)=w ₀ exp(–x/s), wherew ₀ ands are given constants andx is the optical variable. The second approximation assumes the formw(x) = _r=0 ^R d _r ^* p _r (x/a), whered _r ^* are known expansion coefficients anda is the half optical thickness of the slab. Analytic expressions for the forward, backward radiation intensities and fluxes are given in each approximation. The solution of the linear transport equation is performed on the basis of integral Fourier transforms.     

Compatible potentials and orbits in synodic systems

For a given family of orbits f(x,y) = c ^* which can be traced by a material point of unit in an inertial frame it is known that all potentials V(x,y) giving rise to this family satisfy a homogeneous, linear in V(x,y), second order partial differential equation (Bozis,1984). The present paper offers an analogous equation in a synodic system Oxy, rotating with angular velocity . The new equation, which relates the synodic potential function (x,y), = –V(x, y) + % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSqaaSqaai% aaigdaaeaacaaIYaaaaaaa!3780!\[\tfrac{1}{2}\]²(x ² + y ²) to the given family f(x,y) = c ^*, is again of the second order in (x,y) but nonlinear.As an application, some simple compatible pairs of functions (x,y) and f(x, y) are found, for appropriate values of , by adequately determining coefficients both in and f.     

On the adelphic potentials compatible with a set of planar orbits

Given a planar potentialB=B(x, y), compatible with a monoparametric family of planar orbitsf(x, y)=c, we face the problem of producing potentialsA=A(x, y), adelphic toB(x, y), i.e. nontrivial potentials which have in common withB(x, y) the given set of orbits. We establish a linear, second order partial differential equation for a functionP(x, y) and we prove that, to any definite positive solution of this equation, there corresponds a potentialA(x, y) adelphic toB(x, y).     

The anomalous extinction law

In this paper we present the results of the calculations of anomalous extinction laws characterized by different values of the ratio of the total to selective extinctionR=A _v/E(B-V). As far as we are aware, no such theoretical extinction laws have been published before. The calculation is based on the Mie theory of light scattering by small particles assuming a two-grain dust model of Mathiset al. (1977) and by adopting theri particle size (a) distribution, given byn(a)=a ^–3.5, in which the minimum and maximum size is taken to be 0.005 and 0.22 m, respectively. Furthermore, we use the dielectric functions for the two types of grains, graphite and astronomical silicate, as derived by Draine and Lee (1984), published in tabular form by Draine (1985). Following Mathis and Wallenhorst (1981), we then obtained the anomalous extinction laws by changing the upper size cutoff of the particles.A comparison shows that the calculated anomalous extinction laws agree quite well with the corresponding laws derived observationally using the colour-difference method. Finally, a trial and error method for the determination of the anomalous extinction law based on a comparison of the observed extinction-free spectral energy distribution with the corresponding theoretical Kurucz (1979) model, is explained.     

Linear analysis of the light curves of eclipsing variables

In this paper, general model analysis of the form c _{i i} () for the fractional loss of lightf() exhibited by a close binary system will be considered in the sense of the least-squares criterion. A general recursive method will be established for constructing the normal equations for the most useful functions _i (), and an economical storage of the method on a digital computer, with its computational steps, will also be given. Moreover, a full recursive computational algorithm for the least-squares approximation will also be established. By means of this algorithm, all the solution vectors, the variance for different orders of fit and the corresponding variance-covariance matrix could be computed once and for all and, moreover, recursively. The economical storage of the algorithm and its computational steps will be given. Finally, some of the practical difficulties encountered in the application of the least-squares criterion will be analysed, and some techniques for detecting and controlling these difficulties are also given. Numerical examples on the Algol system using a Fourier cosine series of the form c _i cos [(j - 1) /] will be given for illustration.     

The influence of a resisting medium on 3-dimensional periodic orbits of the restricted 3-body problem

Summary The evolution of the extreme values of the functionsx(t),y(t) andz(t) which are the coordinates of the third bodyP in the barycentric rotating frame of reference, when friction is present, is discussed. These values, which are constant for periodic orbits, change due to the presence of the resisting medium. It is shown that either the orbit tends to become circular and coplanar with the two primaries or to collide with one of the primaries.