The accuracy of the numerical findings of Lukes-Gerakopoulos PRD 86, 044013 (2012), regarding the existence of additional first integrals in the Zipoy-Voorhees space-time, was recently questioned by Maciejewski et al. PRD 88, 064003 (2013). In this comment, it is shown that the discrepancy between the results of Lukes-Gerakopoulos (2012) and Maciejewski et al. (2013) is not due to issues related to numerical accuracy, as claimed in Maciejewski et al (2013), but due to different choice of coordinates used in Maciejewski et al.(2013).

In this work, we investigate the eclipse timing of the polar binary HU Aquarii that has been observed for almost two decades. Recently, Qian et al. attributed large (O-C) deviations between the eclipse ephemeris and observations to a compact system of two massive jovian companions. We improve the Keplerian, kinematic model of the Light Travel Time (LTT) effect and re-analyse the whole currently available data set. We add almost 60 new, yet unpublished, mostly precision light curves obtained using the time high-resolution photo-polarimeter OPTIMA, as well as photometric observations performed at the MONET/N, PIRATE and TCS telescopes. We determine new mid--egress times with a mean uncertainty at the level of 1 second or better. We claim that because the observations that currently exist in the literature are non-homogeneous with respect to spectral windows (ultraviolet, X-ray, visual, polarimetric mode) and the reported mid--egress measurements errors, they may introduce systematics that affect orbital fits. Indeed, we find that the published data, when taken literally, cannot be explained by any unique solution. Many qualitatively different and best-fit 2-planet configurations, including self-consistent, Newtonian N-body solutions may be able to explain the data. However, using high resolution, precision OPTIMA light curves, we find that the (O-C) deviations are best explained by the presence of a single circumbinary companion orbiting at a distance of ~4.5 AU with a small eccentricity and having ~7 Jupiter-masses. This object could be the next circumbinary planet detected from the ground, similar to the announced companions around close binaries HW Vir, NN Ser, UZ For, DP Leo or SZ Her, and planets of this type around Kepler-16, Kepler-34 and Kepler-35.

We consider the unrestricted problem of two mutually attracting rigid bodies, an uniform sphere (or a point mass) and an axially symmetric body. We present a global, geometric approach for finding all relative equilibria (stationary solutions) in this model, which was already studied by Kinoshita (1970). We extend and generalize his results, showing that the equilibria solutions may be found by solving at most two non-linear, algebraic equations, assuming that the potential function of the symmetric rigid body is known explicitly. We demonstrate that there are three classes of the relative equilibria, which we call "cylindrical", "inclined co-planar", and "conic" precessions, respectively. Moreover, we also show that in the case of conic precession, although the relative orbit is circular, the point-mass and the mass center of the body move in different parallel planes. This solution has been yet not known in the literature.

We re-analyze the precision radial velocity (RV) observations of HD160691 (mu Ara) by the Anglo-Australian Planet Search Team. The star is supposed to host two Jovian companions (HD160691b, HD160691c) in long-period orbits about 630 days and about 2500 days, respectively) and a hot-Neptune (HD160691d) in about 9 days orbit. We perform a global search for the best fits in the orbital parameters space with a hybrid code employing the genetic algorithm and simplex method. The stability of Keplerian fits is verified with the N-body model of the RV signal that takes into account the dynamical constraints (so called GAMP method). Our analysis reveals a signature of the fourth, yet unconfirmed, Jupiter-like planet HD160691e in about 307 days orbit. Overall, the global architecture of four-planet configuration recalls the Solar system. All companions of mu Ara move in quasi-circular orbits. The orbits of two inner Jovian planets are close to the 2:1 mean motion resonance. The alternative three-planet system involves two Jovian planets in eccentric orbits (e about 0.3), close to the 4:1 MMR, but it yields a significantly worse fit to the data. We also verify a hypothesis of the 1:1 MMR in the subsystem of two inner Jovian planets in the four-planet model.

We perform a dynamical analysis of the recently published radial velocity (RV) measurements of a few solar type stars which host multiple Jupiter-like planets. In particular, we re-analyze the data for HD 202206, 14 Her, HD 37124 and HD 108874. We derive dynamically stable configurations which reproduce the observed RV signals using our method called GAMP (an acronym of the Genetic Algorithm with MEGNO Penalty). The GAMP relies on the N-body dynamics and makes use of genetic algorithms merged with a stability criterion. For this purpose, we use the maximal Lyapunov exponent computed with the dynamical fast indicator MEGNO. Through a dynamical analysis of the phase-space in a neighborhood of the obtained best-fit solutions, we derive meaningful limits on the parameters of the planets. We demonstrate that GAMP is especially well suited for the analysis of the RV data which only partially cover the longest orbital period and/or correspond to multi-planet configurations involved in low-order mean motion resonances (MMRs). In particular, our analysis reveals a presence of a second Jupiter-like planet in the 14 Her system (14 Her c) involved in a 3:1 or 6:1 MMR with the known companion b. We also show that the dynamics of the HD 202206 system may be qualitatively different when coplanar and mutually-inclined orbits of the companions are considered. We demonstrate that the two outer planets in the HD 37124 system may reside in a close neighborhood of the 5:2 MMR. Finally, we found a clear indication that the HD 108874 system may be very close to, or locked in an exact 4:1 MMR.

We perform a dynamical analysis of the recently updated set of the radial velocity (RV) measurements of the HD160691 (mu Arae). The kinematic, 2-Keplerian model of the measurements leads to the best-fit solution in which the eccentricity of the outer planet is about 0.7 and its semi-major axis is about 4AU. The parameters of the inner planet are well determined. The eccentricity is about 0.3 and the semi-major axis is about 1.65AU. The best 2-Keplerian model leads to a catastrophically unstable configuration, self disrupting in less than 20,000yr. To derive dynamically stable configurations which are simultaneously consistent with the RV data, we use the so called GAMP (Genetic Algorithm with MEGNO Penalty). Using this method, we derive meaningful limits on the parameters of the outer planet which provide a stable behavior of the system. The best-fit solutions are located in a shallow valley of Chi^2, in the (a_c,e_c)-plane, extending over 2AU (for the formal 1sigma confidence interval of the best fit). We find two equally good best-fit solutions leading to the qualitatively different orbital configurations. One of them corresponds to the center of the 5:1 mean motion resonance (MMR) and the second one describes a configuration between the 6:1 and 17:2MMRs. The mu Arae system can be found in the zone confined to other low-order MMRs of the type p:1 with p>5. Our results support the classification of the mu Arae as a hierarchical planetary system, dynamically similar to other known multi-planet systems around HD12661 and HD169830. The results of the GAMP on the extended data set are fully consistent with our previous conclusions concerning a much shorter observational window.

It is proposed to use the Lie group theory of symmetries of differential equations to investigate the system of equations describing a static star in a radiative and convective equilibrium. It is shown that the action of an admissible group induces a certain algebraic structure in the set of all solutions, which can be used to find a family of new solutions. We have demonstrated that, in the most general case, the equations admit an infinite parameter group of quasi-homologous transformations. We have found invariants of the symmetries group which correspond to the fundamental relations describing a physical characteristic of the stars such as the Hertzsprung-Russell diagram or the mass-luminosity relation. In this way we can suggest that group invariants have not only purely mathematical sense, but their forms are closely associated with the basic empirical relations.

The set of radial velocity measurements of the HD 160691 has been recently published by Jones et.al 2002 (MNRAS). It reveals a linear trend that indicates a presence of the second planet in this system. The preliminary double-Keplerian orbital fit to the observations, announced by the discovery team, describes a highly unstable, self-disrupting configuration. Because the observational window of the HD 160691 system is narrow, the orbital parameters of the hypothetical second companion are unconstrained. In this paper we try to find out whether a second giant planet can exist up to the distance of Jupiter and search for the dynamical constraints on its orbital parameters. Our analysis employs a combination of fitting algorithms and simultaneous examination of the dynamical stability of the obtained orbital fits. It reveals that if the semi-major axis of the second planet is smaller than $\simeq 5.2$ AU, the observations are consistent with quasi-periodic, regular motions of the system confined to the islands of various low-order mean motion resonances, e.g., 3:1, 7:2, 4:1, 5:1, or to their vicinity. In such cases the second planet has smaller eccentricity $\simeq 0.2-0.5$ than estimated in the previous works. We show that the currently available Doppler data rather preclude the 2:1 mean motion resonance expected by some authors to be present in the HD 160691 planetary system. We also demonstrate that the MEGNO-penalty method, developed in this paper, which is a combination of the genetic minimization algorithm and the MEGNO analysis, can be efficiently used for predicting stable planetary configurations when only a limited number of observations is given or the data do not provide tight constraints on the orbital elements.

In this work we perform a global analysis of the radial velocity curve of the HD 12661 system. Orbital fits that are obtained by the genetic and gradient algorithms of minimization reveal the proximity of the system to the 6:1 mean motion resonance. The orbits are locked in the secular resonance with apsidal axes librating about 180 deg with the full amplitude $\simeq (90,180) deg$. Our solution incorporates the mutual interaction between the companions. The stability analysis with the MEGNO fast indicator shows that the system is located in an extended stable zone of quasi-periodic motions. These results are different from those obtained on the basis of the orbital fit published by Fischer et al. (2003)

We present a new analysis of the dynamics of the planetary system around the pulsar B1257+12. A semi-analytical theory of perturbation between terrestrial-mass planets B and C is developed and applied to improve multi-orbit timing formula for this object. We use numerical simulations of the pulse arrival times for PSR B1257+12 to demonstrate that our new timing model can serve as a toll to determine the masses of the two planets.

We present theoretical analysis of the astrometric searches for extrasolar planets with the Space Interferometry Mission (SIM). Particularly, we derive a model for the future measurements with SIM and discuss the problem of reliable estimation of orbital elements of planets. For this purpose we propose a new method of data analysis and present a numerical test of its application on simulated SIM astrometric measurements of $\upsilon$ Andromedae planetary system. We demonstrate that our approach allows successfull determination of its orbital elements.

Since it has become possible to discovery planets orbiting nearby solar-type stars through very precise Doppler-shift measurements, the role of methods used to analyze such observations has grown significantly. The widely employed model-dependent approach based on the least-squares fit of the Keplerian motion to the radial-velocity variations can be, as we show, unsatisfactory. Thus, in this paper, we propose a new method that may be easily and successfully applied to the Doppler-shift measurements. This method allows us to analyze the data without assuming any specific model and yet to extract all significant features of the observations. This very simple idea, based on the subsequent subtraction of all harmonic components from the data, can be easily implemented. We show that our method can be used to analyze real 16 Cygni B Doppler-shift observations with a surprising but correct result which is substantially different from that based on the least-squares fit of a Keplerian orbit. Namely, using frequency analysis we show that with the current accuracy of this star's observations it is not possible to determine the value of the orbital eccentricity which is claimed to be as high as 0.6.

In this paper we present a new method that can be used for analysis of time of arrival of a pulsar pulses (TOAs). It is designated especially to detect quasi-periodic variations of TOAs. We apply our method to timing observations of PSR B1257+12 and demonstrate that using it it is possible to detect not only first harmonics of a periodic variations, but also the presence of a resonance effect. The resonance effect detected, independently of its physical origin, can appear only when there is a non-linear interaction between two periodic modes. The explanation of TOAs variations as an effect of the existence of planets is, till now, the only known and well justified. In this context, the existence of the resonance frequency in TOAs is the most significant signature of the gravitational interaction of planets.