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.     

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.     

About the Hill stability of extrasolar planets in stellar binary systems

We use a three dimensional generalization of Szebehely’s invariant relation obtained by us (Makó and Szenkovits, Celest. Mech. Dyn. Astron. 90, 51, 2004) in the elliptic restricted three-body problem, to establish more accurate criterion of the Hill stability. By using this criterion, the Hill stability of four extrasolar planets (γ Cephei Ab, Gliese 86 Ab, HD 41004 Ab and HD 41004 Bb) is investigated.     

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.     

A theory of libration

It is shown here that many problems of libration in celestial mechanics can be reduced to a perturbation of anintermediary defined by the Hamiltonian $$F = B\left( y \right) + 2\mu ^2 A\left( y \right)f\left( x \right).$$ This generalization of the Ideal Resonance Problem, with a periodic functionf(x) replacing sin² x, is solved here toO(μ ²) by an algorithm that is essentially the same as the one used in the original formulation. The solution is of the formx=x(u), u=u(t), y=y(x), with the functionx(u) commonly involving the inversion of a hyperelliptic integralu(x), evaluated by quadrature. Libration may be simple or multiple, depending on the nature of the functionf(x) and on the initial conditions. Double libration is illustrated here by the horseshoe-shaped orbits enclosing two libration centers.     

Fast Procedure Solving Universal Kepler's Equation

We developed a procedure to solve a modification of the standard form of the universal Kepler’s equation, which is expressed as a nondimensional equation with respect to a nondimensional variable. After reducing the domain of the variable and the argument by using the symmetry and the periodicity of the equation, the method first separates the case where the solution is so small that it is given an inverted series. Second, it separates the cases where the elliptic, parabolic, or hyperbolic standard forms of Kepler’s equation are suitable. Here the separation is done by judging whether detouring these nonuniversal equations will cause a 1-bit loss of information to their nonuniversal solutions or not. Then the nonuniversal equations are solved by the author’s procedures to solve the elliptic Kepler’s equation (Fukushima, 1997a), Barker’s equation (Fukushima, 1998), and the hyperbolic Kepler’s equation (Fukushima, 1997b), respectively. And their nonuniversal solutions are transformed back to the solution of the universal equation. For the rest of the case, we obtain an approximate solution by solving roughly the approximated cubic equation as we did in solving Barker’s equation. Then the correction to the approximate solution is obtained by Halley’s method precisely. There the special function appeared in the universal equation is rewritten into a combination of similar special functions of small arguments, so that they are efficiently evaluated by their Taylor series. Numerical measurements showed that, in the case of Intel Pentium II processor, the new method is 10–25 times as fast as Shepperd’s method (Shepperd, 1985) and 7–13 times as fast as the standard Newton method. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Dynamic discretization method for solving Kepler’s equation

Kepler’s equation needs to be solved many times for a variety of problems in Celestial Mechanics. Therefore, computing the solution to Kepler’s equation in an efficient manner is of great importance to that community. There are some historical and many modern methods that address this problem. Of the methods known to the authors, Fukushima’s discretization technique performs the best. By taking more of a system approach and combining the use of discretization with the standard computer science technique known as dynamic programming, we were able to achieve even better performance than Fukushima. We begin by defining Kepler’s equation for the elliptical case and describe existing solution methods. We then present our dynamic discretization method and show the results of a comparative analysis. This analysis will demonstrate that, for the conditions of our tests, dynamic discretization performs the best.     

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.     

CCD star images: On the determination of Moffat’s PSF shape parameters

Among the variety of empirical models of optical Point Spread Function used in the astronomical environment, only the Moffat’s (1969) one is able to describe by means of two parameters (in the circular case) both the inner and the outer star image regions. In view of this very important feature, the problem of the simultaneous estimates of Moffat’s PSF shape parameters, off-centring, and the background level in CCD star images has been investigated. The problem does not seem to be rigorously resolvable, but an approximate way to calculate all the parameters except off-centring is shown. It must be stressed that, the Moffat’s PSF model being a softened power law belonging to the family of modified King and Hubble models, the present discussion can be of aid in many other research fields. Also, the integral equation enabling us to convolve a spherical source with Moffat’s PSF is given and applied for comparison to Multi-Gaussian convolution.     

Two exact solutions to the general relativistic Binet’s equation

The general relativistic Binet’s orbit equation is exactly solved using the exp-function method, two exact solutions are obtained which are useful for theoretical analysis.     

Further VLA observations of hydrogen deficient stars

VLA observations at 6 cm have been obtained for three hydrogen-deficient objects υ Sgr, V 348 Sgr, and A bell 58. A bell 58 was also observed at 2 cm. Only upper limits to the flux density could be set for these sources. A new radio source at 6 cm was found in the field of υ Sgr. The upper limit for 6 cm flux density of V348 Sgr sets an upper limit to its reddening asE(B–V) ≤ 0.65. The hydrogen deficient planetary nebula A 58 shows much lower radio flux than expected from the infrared-radio flux density relationship of planetary nebulae. National Radio Astronomy Observatory’s Very Large Array is operated by Associated Universities Inc. under contract with National Science Foundation, USA.     

A new integration theory for conservative two degree-of-freedom systems

A two degree-of-freedom, conservative system is reduced to a single degree-of-freedom, kinematic system with Hamiltonian integral under the change of independent variable: $$dt = \zeta dt (\zeta = \upsilon _x - \upsilon _y )$$ where ζ is the curl (or vorticity) of the velocity field with cartesian inertial componentsu(x, y, t) andv(x, y, t). In the autonomous case whenu _t=v _t=0, orbits are globally represented by the level curves of an autonomous Hamiltonian functionH(x,y) satisfying a second-order quasilinear partial differential equation (Szebehely's Equation): $$2(H + U)\left( {H_{xx} H_y^2 - 2H_{xy} H_x H_y + H_{yy} H_x^2 } \right) + (H_x U_x + H_y U_y )\left( {H_x^2 + H_y^2 } \right) = 0$$ whereU(x, y) is the autonomous potential function. An inversion of dependent and independent variables reduces this equation to a second-order, ordinary differential equation for a function specifying the orbital curve. The true time variable is recovered by evaluating a quadrature. Fundamental differences exist between this approach and Hamilton-Jacobi theory.     

A new class of LRS Bianchi type-II dark energy models with variable EoS parameter

A new class of dark energy models in a Locally Rotationally Symmetric Bianchi type-II (LRS B-II) space-time with variable equation of state (EoS) parameter and constant deceleration parameter have been investigated in the present paper. The Einstein’s field equations have been solved by applying a variation law for generalized Hubble’s parameter given by Berman: Nuovo Cimento 74:182 (1983) which generates two types of solutions for the average scale factor, one is of power-law type and other is of the exponential-law form. Using these two forms, Einstein’s field equations are solved separately that correspond to expanding singular and non-singular models of the universe respectively. The dark energy EoS parameter ω is found to be time dependent and its existing range for both models is in good agreement with the three recent observations of (i) SNe Ia data (Knop et al.: Astrophys. J. 598:102 (2003)), (ii) SNe Ia data collaborated with CMBR anisotropy and galaxy clustering statistics (Tegmark et al.: Astrophys. J. 606:702 (2004)) and latest (iii) a combination of cosmological datasets coming from CMB anisotropies, luminosity distances of high redshift type Ia supernovae and galaxy clustering (Hinshaw et al.: Astrophys. J. Suppl. 180:225 (2009); Komatsu et al. Astrophys. J. Suppl. 180:330 (2009)). The cosmological constant Λ is found to be a positive decreasing function of time and it approaches a small positive value at late time (i.e. the present epoch) which is corroborated by results from recent supernovae Ia observations. The physical and geometric behaviour of the universe have also been discussed in detail.     

Green’s function for the linear Kompaneets equation

Green’s function for the linear Kompaneets equation is calculated; it is expressed in terms of a Whittaker function W_2,iμ(Z) or a MacDonald function K_iμ(z) with a purely imaginary index. A method is proposed for calculating these functions. Langer’s asymptotic solution for large μ is refined in Cherry’s second approximation. With a series expansion for small values of the argument and the asymptotic form for large values, this approximation enables one to calculate Green’s function to five significant figures. Solutions of the Kompaneets equation will be used to estimate the accuracy of numerical methods and to calculate the evolution of the spectrum of a photon gas during Compton scattering, as well as the average frequencies and the dispersion of photon frequencies for different initial spectra. Translated from Astrofizika, Vol. 40, No. 1, pp. 97–116, January–March, 1997.     

Orbit of an Astrometric Binary System

We present a new method to solve the problem of initial orbit determination of any binary system. This method is mainly based on the material available for an observer, for example relative positions at a given time of the couple in the “plane of sky”, namely the tangent plane to the celestial sphere at the position of the primary component. The problem of orbit determination is solved by splitting in successive stages in order to decorrelate the parameters of each other as much as possible. On one hand, the geometric problem is solved using the first Kepler’s law from a single observing run and, on the other hand, dynamical parameters are then inferred from the fit of the Kepler’s equation. At last, the final stage consists in determining the main physical parameters involved in the secular evolution of the system, that is the spin axis and the J₂ parameter of the primary if we assume that it is a quasi-spherical body. As a matter of fact there is no need to make too restrictive initial assumptions (such as circular orbit or zero eccentricity) and initial guesses of parameters required by a non-linear least-squares Levenberg–Marquardt algorithm are finally obtained after each stage. Such a protocol is very useful to study systems like binary asteroids for which all of the parameters should be considered a priori as unknowns. As an example of application, we used our method to estimate the set of the Pluto–Charon system parameters from observations collected in the literature since 1980.     

Ca II K Emission Line Asymmetries Among Red Giants

In the spectra of red giants the chromospheric emission feature found in the core of the Ca II K line often exhibits an asymmetric profile. This asymmetry can be documented by a parameter V/R which is classified as > 1, 1, or < 1 if the violet wing of the emission profile is of greater, equal, or lower intensity than the redward wing. A literature search has been conducted to compile a V/R dataset which builds on the large survey of bright field giants made by Wilson (1976). Among stars of luminosity classes II–III–IV the majority of those with V/R > 1 are found to be bluer than B-V =1.3, while those with V/R < 1 are mostly redder than this colour. Stars with nearly symmetric profiles, V/R≈ 1, are found throughout the colour range 0.8 < B-V < 1.5. There is no sharp transition line separating stars of V/R > 1 and < 1 in the colour-magnitude diagram, but rather a ‘transition zone’ centered at B-V ≈ 1.3. The center of this zone coincides closely with a ‘coronal dividing line’ identified by Haish, Schmitt and Rosso (1991) as the red envelope in the H–R diagram of giants detected in soft x-ray emission by ROSAT. It is suggested that both the transition to a Ca II K emission asymmetry of V/R < 1 and the drop in soft x-ray activity across the coronal dividing line are related to changes in the dynamical state of the chromospheres of red giants. By contrast, the onset of photometric variability due to pulsation occurs among stars of early-M spectral type, that are redward of the mid-point of the Ca II V/R ‘transition zone’, suggesting that the chromospheric motions which produce an asymmetry of V/R < 1 are established prior to the onset of pulsation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mutual phenomena of Jovian satellites

Results of the observations of mutual eclipses of Galilean satellites observed from the Vainu Bappu Observatory during 1985 are presented. Theoretical models assuming a uniform disc, Lambert’s law and Lommel-Seeliger’s law describing the scattering characteristics of the surface of the eclipsed satellite were used to fit the observations. Light curves of the 1E2 event on 1985 September 24 and the 3E1 event on 1985 October 24 observed from VBO and published light curves of the 1E2 event on 1985 September 14, the 3E1 event on 1985 September 26 and the 2E1 event on 1985 October 28 (Arlotet al 1989) were fitted with theoretical light curves using Marquardt’s algorithm. The best fitting was obtained using Lommel-Seeliger’s law to describe the scattering over the surface of Io and Europa. During the fitting, a parameter^δxshift which shifts the theoretical light curve along the direction of relative motion of the eclipsed satellite with respect to the shadow centre, on the sky plane (as seen from the Sun) was determined along with the impact parameter. In absence of other sources like prominent surface features or non perfect sky conditions which could lead to asymmetric light curves,^δxshift would be a measure of the phase correction (Aksnes, Franklin & Magnusson 1986) with an accuracy as that of the midtime. Heliocentric Δα cos (δ) and gDδ at mid times derived from fitted impact parameters are reported     

Brans-Dicke Quantum Cosmology

By solving a Wheeler-De Witt ‘extended’ equation in the Brans-Dicke theory, we have found that the probability distribution predicts: i) An initial value for the Brans-Dicke scalar field φ ∼ ρ^1/2_VAC in the beginning of the inflation, where ρ_VAC is the vacuum density energy (this gives a planck mass ∼ ρ^1/4_VAC) ii) Large values for the Brans-Dicke parameter w. On the other hand it is shown that by taking into account the dynamical behaviour of φ and the matter scalar field σ we can formulate a ‘creation boundary condition’ where in the ‘beginning’ of the Universe (R =0, ‘nothing’ for some authors) we have a dynamical σ already ‘created’. This could be the energetic mechanism which makes Universe tunnels the potential barrier to evolve classically after. Besides we have found the possibility of a cosmological uncertainty principle. This revised version was published online in July 2006 with corrections to the Cover Date.     

Out-of-plane librations of spinning tethered satellite systems

We analyze the out-of-plane librations of a tethered satellite system that is nominally rotating in the orbit plane. To isolate the librational dynamics, the system is modeled as two point masses connected by a rigid rod with the system mass center constrained to an unperturbed circular orbit. For small out-of-plane librations, the in-plane motion is unaffected by the out-of-plane librations and a solution for the in-plane motion is determined in terms of Jacobi elliptic functions. This solution is used in the linearized equation for the out-of-plane librations, resulting in a Hill’s equation. Floquet theory is used to analyze the Hill’s equation, and we show that the out-of-plane librations are unstable for certain ranges of in-plane spin rate. For relatively high in-plane spin rates, the out-of-plane librations are stable, and the Hill’s equation can be approximated by a Mathieu’s equation. Approximate solutions to the Mathieu’s equation are determined, and we analyze the dominant characteristics of the out-of-plane librations for high in-plane spin rates. The results obtained from the analysis of the linearized equations of motion are compared to numerical simulations of the nonlinear equations of motion, as well as numerical simulations of a more realistic system model that accounts for tether flexibility. The instabilities discovered from the linear analysis are present in both the nonlinear system and the more realistic system model. The approximate solutions for the out-of-plane librations compare well to the nonlinear system for relatively high in-plane rotation rates, and also capture the significant qualitative behavior of the flexible system.