Magneto-convection in sunspot penumbrae

Finite amplitude convection in the presence of a horizontal magnetic field has been investigated in a region where thermal diffusivity (κ) is less than magnetic diffusivity (η) and whenκ/η > 1,Q ≤Q _c, where $$Q_c = \frac{{(1 + \sigma _1 )(\pi ^2 + q_c^2 )^2 }}{{q_c^2 (\sigma _2 - \sigma _1 )}}$$ ,Q is the Chandrasekhar number,σ ₁ the Prandtl number,σ ₂ the magnetic Prandtl number, andq _c the critical wave number at the onset of stationary convection. We have derived a nonlinear time-dependent Landau—Ginzburg equation near the onset of supercritical stationary convection and a nonlinear, second-order equation at the Takens—Bogdanov bifurcation. We have obtained steady-state solutions of these equations, which describe the nonlinear behaviour near the onset of stationary convection.     

The seasonal variation of the newtonian constant of gravitation and its relationship with the “classical tests” of general relativity

The ratio between the Earth's perihelion advance (Δθ) _E and the solar gravitational red shift (GRS) (Δø _s e)a ₀/c ² has been rewritten using the assumption that the Newtonian constant of gravitationG varies seasonally and is given by the relationship, first found by Gasanalizade (1992b) for an aphelion-perihelion difference of (ΔG)_a?p . It is concluded that $$\begin{gathered} (\Delta \theta )_E = \frac{{3\pi }}{e}\frac{{(\Delta \phi _{sE} )_{A_0 } }}{{c^2 }}\frac{{(\Delta G)_{a - p} }}{{G_0 }} = 0.038388 \sec {\text{onds}} {\text{of}} {\text{arc}} {\text{per}} {\text{revolution,}} \hfill \\ \frac{{(\Delta G)_{a - p} }}{{G_0 }} = \frac{e}{{3\pi }}\frac{{(\Delta \theta )_E }}{{(\Delta \phi _{sE} )_{A_0 } /c^2 }} = 1.56116 \times 10^{ - 4} . \hfill \\ \end{gathered} $$ The results obtained here can be readily understood by using the Parametrized Post-Newtonian (PPN) formalism, which predicts an anisotropy in the “locally measured” value ofG, and without conflicting with the general relativity.     

Black Hole Dynamic Potentials

In the following paper, certain black hole dynamic potentials have been developed definitively on the lines of classical thermodynamics. These potentials have been refined in view of the small differences in the equations of the laws of black hole dynamics as given by Bekenstein and those of thermodynamics. Nine fundamental black hole dynamical relations have been developed akin to the four fundamental thermodynamic relations of Maxwell. The specific heats C _??,?? and C _J,Q have been defined. For a black hole, these quantities are negative. The ??dA equation has been obtained as an application of these fundamental relations. Time reversible processes observing constancy of surface gravity are considered and an equation connecting the internal energy of the black hole E, the additional available energy defined as the first free energy function K, and the surface gravity ??, has been obtained. Finally as a further application of the fundamental relations, it has been proved for a homogeneous gravitational field in black hole space times or a de Sitter black hole that $C_{\Omega ,\Phi } -C_{J,Q} =\kappa \left[ {\left( {\frac{\partial J}{\partial \kappa }} \right)_{\Omega ,\Phi } \left( {\frac{\partial \Omega }{\partial \kappa }} \right)_{J,Q} +\left( {\frac{\partial Q}{\partial \kappa }} \right)_{\Omega ,\Phi } \left( {\frac{\partial \Phi }{\partial \kappa }} \right)_{J,Q} } \right]$ . This is dubbed as the homogeneous fluid approximation in context of the black holes.     

A few observations of and remarks on V1500 Cygni

The development of the post-nova light curve of V1500 Cyg inUBV andHβ, for 15 nights in September and October 1975 are presented. We confirm previous reports that superimposed on the steady decline of the light curve are small amplitude cyclic variations. The times of maxima and minima are determined. These together with other published values yield the following ephemerides from JD 2 442 661 to JD 2 442 674: $$\begin{gathered} {\text{From}} 17 {\text{points:}} {\text{JD}}_{ \odot \min } = 2 442 661.4881 + 0_{^. }^{\text{d}} 140 91{\text{n}} \hfill \\ \pm 0.0027 \pm 0.000 05 \hfill \\ {\text{From}} 15 {\text{points:}} {\text{JD}}_{ \odot \max } = 2 442 661.5480 + 0_{^. }^{\text{d}} 140 89{\text{n}} \hfill \\ \pm 0.0046 \pm 0.0001 \hfill \\ \end{gathered} $$ with standard errors of the fits of ±0 _. ^d 0052 for the minima and ±0 _. ^d 0091 for the maxima. Assuming V1500 Cyg is similar to novae in M31, we foundr=750 pc and a pre-nova absolute photographic magnitude greater than 9.68.     

Effect of oblateness on the non-linear stability of L 4 in the restricted three-body problem

Non-linear stability of the libration point L ₄ of the restricted three-body problem is studied when the more massive primary is an oblate spheroid with its equatorial plane coincident with the plane of motion, Moser's conditions are utilised in this study by employing the iterative scheme of Henrard for transforming the Hamiltonian to the Birkhoff's normal form with the help of double D'Alembert's series. It is found that L ₄ is stable for all mass ratios in the range of linear stability except for the three mass ratios: $$\begin{gathered} \mu _{c1} = 0.0242{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.1790{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c2} = 0.0135{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0993{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c3} = 0.0109{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0294{\text{ }}...{\text{ }}A_1 . \hfill \\ \end{gathered} $$      

A note on theH-functions of transfer problems in multiplying media

Some useful results and remodelled representations ofH-functions corresponding to the dispersion function $$T\left( z \right) = 1 - 2z^2 \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x/\left( {z^2 - x^2 } \right)} $$ are derived, suitable to the case of a multiplying medium characterized by $$\gamma _0 = \sum\limits_1^n {\int_0^{\lambda r} {Y_r } \left( x \right){\text{d}}x > \tfrac{1}{2} \Rightarrow \xi = 1 - 2\gamma _0< 0} $$      

A Widebinary Solar Companion as a Possible Origin of Sedna-like Objects

Sedna is the first inner Oort cloud object to be discovered. Its dynamical origin remains unclear, and a possible mechanism is considered here. We investigate the parameter space of a hypothetical solar companion which could adiabatically detach the perihelion of a Neptune-dominated TNO with a Sedna-like semimajor axis. Demanding that the TNO’s maximum value of osculating perihelion exceed Sedna’s observed value of 76 AU, we find that the companion’s mass and orbital parameters (m _c , a _c , q _c , Q _c , i _c ) are restricted to $$m_c>rapprox 5\hskip.25em\hbox{M}_{\rm J}\left(\frac{Q_c}{7850\hbox{ AU}} \frac{q_c}{7850\hbox{ AU}}\right)^{3/2}$$ during the epoch of strongest perturbations. The ecliptic inclination of the companion should be in the range $45{\deg}\lessapprox i_c\lessapprox 135{\deg}$ if the TNO is to retain a small inclination while its perihelion is increased. We also consider the circumstances where the minimum value of osculating perihelion would pass the object to the dynamical dominance of Saturn and Jupiter, if allowed. It has previously been argued that an overpopulated band of outer Oort cloud comets with an anomalous distribution of orbital elements could be produced by a solar companion with present parameter values $$m_c\approx 5\hskip.25em\hbox{M}_{\rm J}\left(\frac{9000\hbox{ AU}}{a_c}\right)^{1/2}.$$ If the same hypothetical object is responsible for both observations, then it is likely recorded in the IRAS and possibly the 2MASS databases.     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

Dimensional approach to mass loss: O- and B-type stars

An attempt has been made to obtain an expression for the rate of stellar mass loss using dimensional analysis. The best expression for O and B stars is of the form: $$\dot M = A'{\text{ }}\left( {\frac{1}{{G^{1/2} c^4 }}} \right){\text{ }}L^{\text{2}} {\text{ (}}R/M)^{{\text{3/2}}} .$$ It is also found thatA′ increases as one goes from B→O stars and from O→O(f)→O(f)), but is not sensitive to luminosity.     

Solution of the equation of transfer for interlocked multiplets by the method of Laplace transform and Winer-Hopf technique with planck function as a nonlinear function of optical depth

The equation of transfer for interlocked multiplets has been solved by Laplace transformation and the Wiener-Hopf technique developed by Dasgupta (1978) considering two nonlinear forms of Planck function: i.e., (a) $$B{\text{ }}_{\text{v}} (T) = B(t) = b_0 + b_1 {\text{ }}e^{ - \alpha t} ,$$ (b) $$B{\text{ }}_{\text{v}} (T) = B(t) = b_0 + b_1 t + b_2 E_2 (t).$$ Solutions obtained by Dasgupta (1978) or by Chandrasekhar (1960) may be obtained from our solutions by dropping the nonlinear terms.     

Nonlinear newtonian theory and the effect of time dilatation

The fact that the energy density ρ_g of a static spherically symmetric gravitational field acts as a source of gravity, gives us a harmonic function \(f\left( \varphi \right) = e^{\varphi /c^2 } \) , which is determined by the nonlinear differential equation $$\nabla ^2 \varphi = 4\pi k\rho _g = - \frac{1}{{c^2 }}\left( {\nabla \varphi } \right)^2 $$ Furthermore, we formulate the infinitesimal time-interval between a couple of events measured by two different inertial observers, one in a position with potential φ-i.e., dt _φ and the other in a position with potential φ=0-i.e., dt ₀, as $${\text{d}}t_\varphi = f{\text{d}}t_0 .$$ When the principle of equivalence is satisfied, we obtain the well-known effect of time dilatation.     

Normality condition in the ideal resonance problem

The Ideal Resonance Problem, defined by the Hamiltonian $$F = B(y) + 2\mu ^2 A(y)\sin ^2 x,\mu \ll 1,$$ has been solved in Garfinkelet al. (1971). As a perturbed simple pendulum, this solution furnishes a convenient and accurate reference orbit for the study of resonance. In order to preserve the penduloid character of the motion, the solution is subject to thenormality condition, which boundsAB" andB' away from zero indeep and inshallow resonance, respectively. For a first-order solution, the paper derives the normality condition in the form $$pi \leqslant max(|\alpha /\alpha _1 |,|\alpha /\alpha _1 |^{2i} ),i = 1,2.$$ Herep _i are known functions of the constant ‘mean element’y', α is the resonance parameter defined by $$\alpha \equiv - {\rm B}'/|4AB\prime \prime |^{1/2} \mu ,$$ and $$\alpha _1 \equiv \mu ^{ - 1/2}$$ defines the conventionaldemarcation point separating the deep and the shallow resonance regions. The results are applied to the problem of the critical inclination of a satellite of an oblate planet. There the normality condition takes the form $$\Lambda _1 (\lambda ) \leqslant e \leqslant \Lambda _2 (\lambda )if|i - tan^{ - 1} 2| \leqslant \lambda e/2(1 + e)$$ withΛ ₁, andΛ ₂ known functions of λ, defined by $$\begin{gathered} \lambda \equiv |\tfrac{1}{5}(J_2 + J_4 /J_2 )|^{1/4} /q, \hfill \\ q \equiv a(1 - e). \hfill \\ \end{gathered}$$      

Changes in binding energy of interacting galaxies

It is shown that the fractional increase in binding energy of a galaxy in a fast collision with another galaxy of the same size can be well represented by the formula $$\xi _2 = 3({G \mathord{\left/ {\vphantom {G {M_2 \bar R}}} \right. \kern-\nulldelimiterspace} {M_2 \bar R}}) ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {V_p }}} \right. \kern-\nulldelimiterspace} {V_p }})^2 e^{ - p/\bar R} = \xi _1 ({{M_1 } \mathord{\left/ {\vphantom {{M_1 } {M_2 }}} \right. \kern-\nulldelimiterspace} {M_2 }})^3 ,$$ whereM ₁,M ₂ are the masses of the perturber and the perturbed galaxy, respectively,V _p is the relative velocity of the perturber at minimum separationp, and \(\bar R\) is the dynamical radius of either galaxy.     

Sur un modele explicitant les differents types morphologiques des galaxies

In the now classical Lindblad-Lin density-wave theory, the linearization of the collisionless Boltzmann equation is made by assuming the potential functionU expressed in the formU=U ₀ + \(\tilde U\) +... WhereU ₀ is the background axisymmetric potential and \(\tilde U<< U_0 \) . Then the corresponding density distribution is \(\rho = \rho _0 + \tilde \rho (\tilde \rho<< \rho _0 )\) and the linearized equation connecting \(\tilde U\) and the component \(\tilde f\) of the distribution function is given by $$\frac{{\partial \tilde f}}{{\partial t}} + \upsilon \frac{{\partial \tilde f}}{{\partial x}} - \frac{{\partial U_0 }}{{\partial x}} \cdot \frac{{\partial \tilde f}}{{\partial \upsilon }} = \frac{{\partial \tilde U}}{{\partial x}}\frac{{\partial f_0 }}{{\partial \upsilon }}.$$ One looks for spiral self-consistent solutions which also satisfy Poisson's equation $$\nabla ^2 \tilde U = 4\pi G\tilde \rho = 4\pi G\int {\tilde f d\upsilon .} $$ Lin and Shu (1964) have shown that such solutions exist in special cases. In the present work, we adopt anopposite proceeding. Poisson's equation contains two unknown quantities \(\tilde U\) and \(\tilde \rho \) . It could be completelysolved if a second independent equation connecting \(\tilde U\) and \(\tilde \rho \) was known. Such an equation is hopelesslyobtained by direct observational means; the only way is to postulate it in a mathematical form. In a previouswork, Louise (1981) has shown that Poisson's equation accounted for distances of planets in the solar system(following to the Titius-Bode's law revised by Balsano and Hughes (1979)) if the following relation wasassumed $$\rho ^2 = k\frac{{\tilde U}}{{r^2 }} (k = cte).$$ We now postulate again this relation in order to solve Poisson's equation. Then, $$\nabla ^2 \tilde U - \frac{{\alpha ^2 }}{{r^2 }}\tilde U = 0, (\alpha ^2 = 4\pi Gk).$$ The solution is found in a classical way to be of the form $$\tilde U = cte J_v (pr)e^{ - pz} e^{jn\theta } $$ wheren = integer,p =cte andJ _v (pr) = Bessel function with indexv (v ² =n ² + α²). By use of the Hankel function instead ofJ _v (pr) for large values ofr, the spiral structure is found to be given by $$\tilde U = cte e^{ - pz} e^{j[\Phi _v (r) + n\theta ]} , \Phi _v (r) = pr - \pi /2(v + \tfrac{1}{2}).$$ For small values ofr, \(\tilde U\) = 0: the center of a galaxy is not affected by the density wave which is onlyresponsible of the spiral structure. For various values ofp,n andv, other forms of galaxies can be taken into account: Ring, barred and spiral-barred shapes etc. In order to generalize previous calculations, we further postulateρ ₀ =kU ₀/r ², leading to Poisson'sequation which accounts for the disc population $$\nabla ^2 U_0 - \frac{{\alpha ^2 }}{{r^2 }}U_0 = 0.$$ AsU ₀ is assumed axisymmetrical, the obvious solution is of the form $$U_0 = \frac{{cte}}{{r^v }}e^{ - pz} , \rho _0 = \frac{{cte}}{{r^{2 + v} }}e^{ - pz} .$$ Finally, Poisson's equation is completely solvable under the assumptionρ =k(U/r ². The general solution,valid for both disc and spiral arm populations, becomes $$U = cte e^{ - pz} \left\{ {r^{ - v} + } \right.\left. {cte e^{j[\Phi _v (r) + n\theta ]} } \right\},$$ The density distribution along the O _z axis is supported by Burstein's (1979) observations.     

On equivalent widths in a penumbral spectrum (a test of the Moe-Maltby model)

From new observational material we made a curve of growth analysis of the penumbra of a large, stable sunspot. The analysis was done relative to the undisturbed photosphere and gave the following results (⊙ denotes photosphere, * denotes penumbra): $$\begin{gathered} (\theta ^ * - \theta ^ \odot )_{exe} = 0.051 \pm 0.007 \hfill \\ {{\xi _t ^ * } \mathord{\left/ {\vphantom {{\xi _t ^ * } {\xi _t }}} \right. \kern-\nulldelimiterspace} {\xi _t }}^ \odot = 1.3 \pm 0.1 \hfill \\ {{P_e ^ * } \mathord{\left/ {\vphantom {{P_e ^ * } {P_e ^ \odot = 0.6 \pm 0.1}}} \right. \kern-\nulldelimiterspace} {P_e ^ \odot = 0.6 \pm 0.1}} \hfill \\ {{P_g ^ * } \mathord{\left/ {\vphantom {{P_g ^ * } {P_g }}} \right. \kern-\nulldelimiterspace} {P_g }}^ \odot = 1.0 \pm 0.2 \hfill \\ \end{gathered} $$ The results of the analysis are in satisfactory agreement with the penumbral model as published by Kjeldseth Moe and Maltby (1969). Additionally we tested this model by computing the equivalent widths of 28 well selected lines and comparing them with our observations.     

Static spherically-symmetric scalar-field theory in general relativity

A spherically-symmetric static scalar field in general relativity is considered. The field equations are defined by $$\begin{gathered} R_{ik} = - \mu \varphi _i \varphi _k ,\varphi _i = \frac{{\partial \varphi }}{{\partial x^i }}, \varphi ^i = g^{ik} \varphi _k , \hfill \\ \hfill \\ \end{gathered} $$ where ?=?(r,t) is a scalar field. In the past, the same problem was considered by Bergmann and Leipnik (1957) and Buchdahl (1959) with the assumption that ?=?(r) be independent oft and recently by Wyman (1981) with the assumption ?=?(r, t). The object of this paper is to give explicit results with a different approach and under a more general condition $$\phi _{;i}^i = ( - g)^{ - 1/2} \frac{\partial }{{\partial x^i }}\left[ {( - g)^{1/2} g^{ik} \frac{\partial }{{\partial x^k }}} \right] = - 4\pi ( -g )^{ - 1/2} \rho $$ where ?=?(r, t) is the mass or the charge density of the sources of the field.     

On pulsar-driven isothermal blast-wave models of supernova remanants

We investigate the ‘equilibrium’ and stability of spherically-symmetric self-similar isothermal blast waves with a continuous post-shock flow velocity expanding into medium whose density varies asr ^?ω ahead of the blast wave, and which are powered by a central source (a pulsar) whose power output varies with time ast ^ω?3. We show that:

1. for ω<0, no physically acceptable self-similar solution exists;
2. for ω>3, no solution exists since the mass swept up by the blast wave is infinite;
3. ? must exceed zero in order that the blast wave expand with time, but ?<2 in order that the central source injects a finite total energy into the blast wave;
4. for 3>ω_min(?)>ω>ω_max(?)>0, where $$\begin{gathered} \omega _{\min } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} + {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} + {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \omega _{\max } (\varphi ){\text{ }} = {\text{ }}2[5{\text{ }} - {\text{ }}\varphi {\text{ }} - {\text{ }}(10{\text{ }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 )^{1/2} ]^2 [2{\text{ }} - {\text{ (10 }} + {\text{ 4}}\varphi {\text{ }} - {\text{ 2}}\varphi ^2 {\text{)}}^{{\text{1/2}}} ]^{ - 2} , \hfill \\ \end{gathered} $$ two critical points exist in the flow velocity versus position plane. The physically acceptable solution must pass through the origin with zero flow speed and through the blast wave. It must also pass throughboth critical points if \(\varphi > \tfrac{5}{3}\) , while if \(\varphi< \tfrac{5}{3}\) it must by-pass both critical points. It is shown that such a solution exists but a proper connection at the lower critical point (for ?>5/3) (through whichall solutions pass with thesame slope) has not been established;
5. for 3>ω>ω_min(?) it is shown that the two critical points of (iv) disappear. However a new pair of critical points form. The physically acceptable solution passing with zero flow velocity through the origin and also passing through the blast wave mustby-pass both of the new critical points. It is shown that the solution does indeed do so;
6. for 3>ω_min(?)>ω_max(?)>ω it is shown that the dependence of the self-similar solution on either ω or ? is non-analytic and therefore, inferences drawn from any solutions obtained in ω>ω_max(?) (where the dependence of the solutionis analytic on ω and ?) are not valid when carried over into the domain 3>ω_min(?)>ω_max(?)>ω;
7. all of the physically acceptable self-similar solutions obtained in 3>ω>0 are unstable to short wavelength, small amplitude but nonself-similar radial velocity perturbations near the origin, with a growth which is a power law in time;
8. the physical self-similar solutions are globally unstable in a fully nonlinear sense to radial time-dependent flow patterns. In the limit of long times, the nonlinear growth is a power law in time for 5<ω+2?, logarithmic in time for 5>ω+2?, and the square of the logarithm in time for 5=ω+2?.

The results of (vii) and (viii) imply that the memory of the system to initial and boundary values does not decay as time progresses and so the system does not tend to a self-similar form. These results strongly suggest that the evolution of supernova remnants is not according to the self-similar form.     

The diffuse cosmic X-ray background spectrum in the 2–10 keV energy range

We report a measurement of the background spectrum based on 10000 counts observed in the energy range 2–10 keV. The rocketborne detector system was optimised for cosmic ray noise rejection. A best fit power law spectrum $$\frac{{dN}}{{dE}} = 16E^{ - 1.8} photons{\text{ }}cm^{ - 2} s^{ - 1} sr^{ - 1} keV^{ - 1} .$$ resulted from the analysis. At 10 keV this result is consistent with recently assessed higher energy data. We show therefore that the change in spectral index between 5 and 50 ke V is approximately ?0.2.     

On the integrability cases of the equation of motion for a satellite in an axially symmetric gravitational field

The projection of an axially symmetric satellite's orbit on a plane perpendicular to the rotation axis (z=const.) is given by the second-order differential equation. $$\frac{{y''}}{{1 + y'^2 }} = \bar \Psi _y - y'\bar \Psi _{x,}$$ where the prime denotes the derivative with respect tox and \(\bar \Psi (x,y)\) is a known function. Two integrability cases have been investigated and it has been shown that for these two cases the integration can be carried out either by quadratures or reduced to a first-order differential equation. Analytical and physical properties are expressed, and it is shown that the equation can be derived from the calssical plane eikonal equation of geometric optics.     

A simple procedure to extend the Gauss method of determining orbital parameters from three to N points

A simple procedure is developed to determine orbital elements of an object orbiting in a central force field which contribute more than three independent celestial positions. By manipulation of formal three point Gauss method of orbit determination, an initial set of heliocentric state vectors r _i and $\dot{\mathbf{r}}_{i}$ is calculated. Then using the fact that the object follows the path that keep the constants of motion unchanged, I derive conserved quantities by applying simple linear regression method on state vectors r _i and $\dot{\mathbf{r}}_{i}$ . The best orbital plane is fixed by applying an iterative procedure which minimize the variation in magnitude of angular momentum of the orbit. Same procedure is used to fix shape and orientation of the orbit in the plane by minimizing variation in total energy and Laplace Runge Lenz vector. The method is tested using simulated data for a hypothetical planet rotating around the sun.