An expansion of the Voigt function

The Voigt function and its derivatives are represented by means of series in Hermite polynomials. The equations obtained can be used both for numerical calculations of these functions and for analytical research.Translated from Astrofizika, Vol. 39, No. 4, pp. 619–625, October–December, 1996.     

Quantitative predictions with detuned normal forms

The phase-space structure of two families of galactic potentials is approximated with a resonant detuned normal form. The normal form series is obtained by a Lie transform of the series expansion around the minimum of the original Hamiltonian. Attention is focused on the quantitative predictive ability of the normal form. We find analytical expressions for bifurcations of periodic orbits and compare them with other analytical approaches and with numerical results. The predictions are quite reliable even outside the convergence radius of the perturbation and we analyze this result using resummation techniques of asymptotic series.     

YORP torque as the function of shape harmonics

The second-order analytical approximation of the mean Yarkovsky–O'Keefe–Radzievskii–Paddack (YORP) torque components is given as an explicit function of the shape spherical harmonics coefficients for a sufficiently regular minor body. The results are based upon a new expression for the insolation function, significantly simpler than in previous works. Linearized plane-parallel model of the temperature distribution derived from the insolation function allows us to take into account a non-zero conductivity. Final expressions for the three average components of the YORP torque related with rotation period, obliquity and precession are given in a form of the Legendre series of the cosine of obliquity. The series have good numerical properties and can be easily truncated according to the degree of the Legendre polynomials or associated functions, with first two terms playing the principal role.     

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.     

GRB 060418 and 060607A: the medium surrounding the progenitor and the weak reverse shock emission

We provide a remedy for a recently published formulation of the Voigt function by reformulating the function into a single proper integral with a damped sine integrand. The present formulation clears up concerns highlighted about the original formulation. The reduction of the Voigt function to a single proper integral enables the use of algorithms available in the literature and included in many software packages to integrate the function and to evaluate the line profile with relative simplicity and superior accuracy. Evidence of the usefulness and superior accuracy of the new formulation is provided.     

The Expansion of the Hamiltonian of the Planetary Problem into the Poisson Series in All Keplerian Elements (Theory)

The study of the evolution of planetary systems, primarily of the Solar System, is one of the basic problems of celestial mechanics. The stability of motion of giant planets on cosmogonic time scales was established by numerical and analytical methods, but the question about the evolution of orbits of terrestrial planets and arbitrary solar-type planetary systems remained open. This work initiates a series of papers allowing one to advance in solving the problem of the evolution of the solar-type planetary systems on cosmogonic time scales by using powerful analytical tools. In the first paper of this series, we choose the optimum reference system and obtain the Poisson series expansion of the Hamiltonian of the problem in all Keplerian elements. We propose to use the integral representation of the corresponding coefficients or the Poisson processor means instead of conventionally addressing any possible special functions. This approach extremely simplifies the algorithm. The next paper of this series deals with the calculation of the expansion coefficients.     

Voigt profile fitting to quasar absorption lines: an analytic approximation to the Voigt–Hjerting function

The Voigt–Hjerting function is fundamental in order to correctly model the profiles of absorption lines imprinted on the spectra of bright background sources by intervening absorbing systems. In this work, we present a simple analytic approximation to this function in the context of absorption-line profiles of intergalactic H  i absorbers. Using basic calculus tools, we derive an analytic expression for the Voigt–Hjerting function that contains only fourth-order polynomial and Gaussian functions. In connection with the absorption coefficient of intergalactic neutral hydrogen, this approximation is suitable for modelling Voigt profiles with an accuracy of 10⁻⁴ or better for an arbitrary wavelength baseline, for column densities up to   N _H I= 10 ²² cm⁻²  , and for damping parameters   a ≲ 10⁻⁴  , that is, the entire range of parameters characteristic to all Lyman transitions arising in a variety of H  i absorbing systems such as Lyman α (Lyα) forest clouds, Lyman limit systems and damped Lyα systems. We hence present an approximation to the Voigt–Hjerting function that is both accurate and flexible to implement in various types of programming languages and machines, and with which Voigt profiles can be calculated in a reliable and very simple manner.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Analytical evaluation of the relativistic thermodynamic functions using binomial expansion theorem and incomplete Gamma functions

A simple evaluation algorithm for calculation of the fermions and bosons relativistic thermodynamic functions with high accuracy is presented. This method is based on the use of the binomial expansion theorems and incomplete Gamma functions. The importance of this result is that it enables a series analytical formula to determine the thermodynamic functions of fermions and bosons in cases of degeneracy and relativity. The formulas obtained are numerically stable for ψ < 0. Examples of comparative calculations are given. Numerical results verify the effectiveness of the present approach and they are in good agreement with literature.     

Line formation in moving media: Asymptotic expansions of some special functions

Line formation in the spectrum of a moving medium with a spherical geometry is considered. In the Sobolev approximation there are some special functions that determine the source function and the force of radiation pressure in the line. The most important case is that of a small dimensionless velocity gradient (i.e., a large dimensionless Sobolev length τ) and a small ratio β of the opacity in the continuum to the opacity in the line. Until now there has been no detailed analytical information about the asymptotic behavior of these functions. For the case of a Doppler profile of the absorption coefficient, we clarify the nontrivial structure of their total asymptotic expansions for τ » 1, β « 1, and arbitrary Βτ. We give an algorithm for obtaining all the coefficients of these expansions and give explicit expressions for the first few coefficients. We also compare the asymptotic expansions with the numerical calculations of these functions available in the literature. We also briefly consider the case of a power-law decrease in the absorption coefficient in the line wing (and, in more detail, the case of Lorentz wings of the Voigt profile).     

Radius of the Roche equipotential surfaces

In literature, there is no exact analytical solution available for determining the radius of Roche equipotential surfaces of distorted close binary systems in synchronous rotation. However, Kopal (Roche Model and Its Application to Close Binary Systems, Advances in Astronomy and Astrophysics, Academic Press, New York 1972) and Morris (Publ. Astron. Soc. Pac. 106:154, 1994) have provided the approximate analytical solutions in the form of infinite mathematical series. These series expressions have been commonly used by various authors to determine the radius of the Roche equipotential surfaces, and hence the equilibrium structures of rotating stars and stars in the binary systems. However, numerical results obtained from these approximating series expressions are not very accurate. In the present paper, we have expanded these series expressions to higher orders so as to improve their accuracy. The objective of this paper is to check, whether, there is any effect on the accuracy of these series expressions when the terms of higher orders are considered. Our results show that in most of the cases these expanded series give better results than the earlier series. We have further used these expanded series to find numerically the volume radius of the Roche equipotential surfaces. The obtained results are in good agreement with the results available in literature. We have also presented simple and accurate approximating formulas to calculate the radius of the primary component in a close binary system. These formulas give very accurate results in a specified range of mass ratio.     

Gravitational potential energy of interpenetrating spherical galaxies in Hernquist’s model

Hernquist’s (1990) mass model for spherical galaxies and bulges described by the deVaucouleur’s profile gives analytical expressions for the density profile and the potential. These have been used to derive a simple and exact analytical expression for the gravitational potential energy of a pair of interpene-trating spherical galaxies represented by this model. The results are compared with those for polytropic and Plummer models of galaxis.     

Integration of the elliptic restricted three-body problem with Lie series

The method of Lie series is used to construct a solution for the elliptic restricted three body problem. In a synodic pulsating coordinate system, the Lie operator for the motion of the third infinitesimal body is derived as function of coordinates, velocities and true anomaly of the primaries. The terms of the Lie series for the solution are then calculated with recurrence formulae which enable a rapid successive calculation of any desired number of terms. This procedure gives a very useful analytical form for the series and allows a quick calculation of the orbit.The project is supported by the Austrian Fonds zur Förderung der wissénschaftlichen Forschung under Project No. 4471.     

Gamma-ray bursts: optical afterglows in the deep Newtonian phase

Gamma-ray burst remnants become trans-relativistic typically in days to tens of days, and they enter the deep Newtonian phase in tens of days to months, during which the majority of shock-accelerated electrons will no longer be highly relativistic. However, a small portion of electrons are still accelerated to ultra-relativistic speeds and are capable of emitting synchrotron radiation. The distribution function for electrons is re-derived here so that synchrotron emission from these relativistic electrons can be calculated. Based on the revised model, optical afterglows from both isotropic fireballs and highly collimated jets are studied numerically, and compared to analytical results. In the beamed cases, it is found that, in addition to the steepening due to the edge effect and the lateral expansion effect, the light curves are universally characterized by a flattening during the deep Newtonian phase.     

An analytical expression for apsidal superhump precession and comparisons with numerical simulations and dwarf nova observations

We compare analytical expressions of precession rates from apsidal (positive) superhumps in close binary systems with numerical disc simulation results and observed values. In the analytical expressions, we include both the dynamical effects on the precession of the disc and effects caused by pressure forces that have been theorized to provide a retrograde effect (i.e. slowing) on the prograde disc precession. We establish new limits on density wave pitch angle to a normalized disc sound speed 60≥Ω_orb  d  tan  i / c >2.214 . Using average values for the density wave pitch angle i and speed of sound c , we find good correlation between numerical simulations and the analytical expression for the apsidal superhump period excess, which includes both the prograde and retrograde effects, for mass ratios of 0.025≤ q ≤0.33 . We also show good correlations with the four known eclipsing systems, OY Car, Z Cha, HT Cas, and WZ Sge. Our analytical expression for apsidal superhump period excess as a function of orbital period is consistent with the trend found in observed systems.     

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.     

任意倾角的共轨运动研究

共轨运动天体与摄动天体的半长径相同,处于1:1平运动共振中.太阳系内多个行星的特洛伊天体即为处于蝌蚪形轨道的共轨运动天体,其中一些高轨道倾角特洛伊天体的轨道运动与来源仍未被完全理解.利用一个新发展的适用于处理1:1平运动共振的摄动函数展开方式,对三维空间中的共轨运动进行考察,计算不同初始轨道根数情况下共轨轨道的共振中心、共振宽度,分析轨道类型与初始轨道根数的关系.并将分析方法所得结果与数值方法的结果相互比较验证,得到了广阔初始轨道根数空间内共轨运动的全局图景.     

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.     

On the atmosphere with exponential turbulence

We present an exact calculation for the photon-average bending angle and intensity produced during occultation of a point source by an isothermal atmosphere with an exponential height dependence of turbulent power. The calculation is performed via an expansion in powers of the gradients of average refractivity and turbulent power. Conditions for the validity of the expansion are discussed. We compare results with other expressions, for the same distribution of refractivity, recently presented by V. R. Eshleman and B. S. Haugstad. Discrepancies are noted. Certain properties of photon averages, which may be of importance for the reduction of occultation data, are discussed.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, the literal analytical expressions for the Fourier expansion of the Earth's spherical harmonic potential will be explored in terms of the sectorial variables _j ⁽ⁱ⁾ introduced in Paper IV (Sharaf, 1982) to regularize the highly-oscillating perturbation force of some orbital systems.Now at the Department of Astronomy, King Abdulaziz University Jeddah, Saudi Arabia.