We study proximity of the Kepler-25 planetary system to a periodic configuration, which is known to be the final state of a system that undergoes smooth migration resulting from the planet-disc interaction. We show that the system is close to the periodic configuration of 2:1 mean motion resonance (MMR) what indicates that its past migration was neither disturbed significantly by turbulence in the disc nor the orbits were perturbed by planetesimals that left after the disc dispersal. We show that, because of the TTV model degeneracy, a periodic configuration is difficult to be found when the standard modelling of the transit timing variations (TTVs) is used. The TTV signal of a periodic configuration (with anti-aligned apsidal lines) may be misinterpreted as an aligned non-resonant system. We demonstrate that the standard MCMC modelling of the Kepler-25 TTVs is very sensitive to an a~priori information on the eccentricities (prior probability distributions). Wide priors (of the order of the ones typically used in the literature) result in favouring the aligned non-resonant configurations with small planets' masses and moderate eccentricities, while for the narrower priors the most likely are the anti-aligned resonant systems with larger masses and very low eccentricities.

The HR8799 planetary system with four ~10 mJup planets in wide orbits up to 70 au, and periods up to 500 yr has been detected with the direct imaging. Its intriguing orbital architecture is not fully resolved due to time-limited astrometry covering ~20 years. Earlier, we constructed a heuristic model of the system based on rapid, convergent migration of the planets. We developed a better structured and CPU-efficient variant of this model. We re-analyzed the self-consistent, homogeneous astrometric dataset in Konopacky (2016). The best-fitting configuration agrees with our earlier findings. The planets are likely involved in dynamically robust Laplace resonance chain. Hypothetical planets with masses below the current detection limit of 0.1-3 Jupiter masses, within the observed inner, or beyond the outer orbit, respectively, do not influence the long term stability of the system. We predict positions of such non-detected objects. The long-term stable orbital model of the observed planets helps to simulate the dynamical structure of debris disks in the system. A CPU-efficient fast indicator technique makes it possible to reveal their complex, resonant shape in 10^6 particles scale. We examine the inner edge of the outer disk detected between 90-145 au. We also reconstruct the outer disk assuming that it has been influenced by convergent migration of the planets. A complex shape of the disk strongly depends on various dynamical factors, like orbits and masses of non-detected planets. It may be highly non-circular and its models are yet non-unique, regarding both observational constraints, as well as its origin.

The cyclic behaviour of (O-C) residuals of eclipse timings in the sdB+M eclipsing binary NSVS 14256825 was previously attributed to one or two Jovian-type circumbinary planets. We report 83 new eclipse timings that not only fill in the gaps in those already published but also extend the time span of the (O-C) diagram by three years. Based on the archival and our new data spanning over more than 17 years we re-examined the up to date system (O-C). The data revealed systematic, quasi-sinusoidal variation deviating from an older linear ephemeris by about 100 s. It also exhibits a maximum in the (O-C) near JD 2,456,400 that was previously unknown. We consider two most credible explanations of the (O-C) variability: the light propagation time due to the presence of an invisible companion in a distant circumbinary orbit, and magnetic cycles reshaping one of the binary components, known as the Applegate or Lanza-Rodono effect. We found that the latter mechanism is unlikely due to the insufficient energy budget of the M-dwarf secondary. In the framework of the third-body hypothesis, we obtained meaningful constraints on the Keplerian parameters of a putative companion and its mass. Our best-fitting model indicates that the observed quasi-periodic (O-C) variability can be explained by the presence of a brown dwarf with the minimal mass of 15 Jupiter masses rather than a planet, orbiting the binary in a moderately elliptical orbit (~ 0.175) with the period of ~ 10 years. Our analysis rules out two planets model proposed earlier.

We present multi-epoch astrometric observations of the M9 ultra-cool dwarf TVLM513-46546 that is placed at the brown dwarf boundary. The new observations have been performed with the European VLBI Network (EVN) at 6cm band. The target has been detected at 7~epochs spanning three years, with measured quiescent emission flux in the range 180-300 muJy. We identified four short-duration flaring events (0.5-2 mJy) with very high circular polarization (~75%-100%). Properties of the observed radio flares support the physical model of the source that is characterized by the electron cyclotron maser instability responsible for outbursts of radio emission. Combined with Very Long Baseline Array (VLBA) earlier data, our detections make it possible to refine the absolute parallax pi=93.27(+0.18,-0.17) mas. Our measurements rule out TVLM513-46546 companions more massive than Jupiter in orbits with periods longer than ~1 yr.

We analyse the Transit Timing Variation (TTV) measurements of a~system of two super-Earths detected as Kepler-29, in order to constrain the planets' masses and orbital parameters. A dynamical analysis of the best-fitting configurations constrains the masses to be $\sim 6$ and $\sim 5$ Earth masses for the inner and the outer planets, respectively. The analysis also reveals that the system is likely locked in the 9:7~mean motion resonance. However, a variety of orbital architectures regarding eccentricities and the relative orientation of orbits is permitted by the observations as well as by stability constraints. We attempt to find configurations preferred by the planet formation scenarios as an additional, physical constraint. We show that configurations with low eccentricities and anti-aligned apsidal lines of the orbits are a natural and most likely outcome of the convergent migration. However, we show that librations of the critical angles are not necessary for the Kepler-29 system to be dynamically resonant, and such configurations may be formed on the way of migration as well. We argue, on the other hand, that aligned configurations with $e \gtrsim 0.03$ may be not consistent with the migration scenario.

The main aim of this work is to estimate possible radio GHz emission of extrasolar planets and brown dwarfs and to check if such radiation can be detected by Very Large Baseline Interferometers (VLBI). In the estimation we assume that the emission may originate in processes similar to those observed in the Jupiter system. The frequency of the radio emission that is produced in this system depends mostly on the magnetic field strength. Jupiter's magnetic field ($\sim 9$ G on average) allows for radiation from kHz frequencies up to 40 MHz. This is well below the frequency range of VLBI. However, it was demonstrated that the magnetic field strength in massive and young object may be up to two orders of magnitude higher than for Jupiter, which is especially relevant for planets around short-lived A type stars. This should extend the range of the emission up to GHz frequencies. We calculated expected flux densities of radio emission for a variety of hypothetical young planetary systems. We analysed two different emission scenarios, and found that the radiation induced by moons (process similar to Jupiter-Io interactions) appears to be less efficient than the emission generated by a stellar wind on a planetary magnetosphere. We also estimated hypothetical emission of planets and brown dwarfs located around relatively young and massive main sequence A-type stars. Our results show that the emission produced by stellar winds could be detected by currently operating VLBI networks.

We investigate the dynamical stability of the Kepler-60 planetary system with three super-Earths. We first determine their orbital elements and masses by Transit Timing Variation (TTV) data spanning quarters Q1-Q16 of the KEPLER mission. The system is dynamically active but the TTV data constrain masses to ~4 Earth masses and orbits in safely wide stable zones. The observations prefer two types of solutions. The true three-body Laplace MMR exhibits the critical angle librating around 45 degrees and aligned apsides of the inner and outer pair of planets. In the Laplace MMR formed through a chain of two-planet 5:4 and 4:3 MMRs, all critical angles librate with small amplitudes of ~30 degrees and apsidal lines in planet's pairs are anti-aligned. The system is simultaneously locked in a three-body MMR with librations amplitude of ~10 degrees. The true Laplace MMR can evolve towards a chain of two-body MMRs in the presence of planetary migration. Therefore the three-body MMR formed in this way seems to be more likely state of the system. However, the true three-body MMR cannot be disregarded a priori and it remains a puzzling configuration that may challenge the planet formation theory.

We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the world. They reveal a systematic deviation of ~60 - 120 seconds from the older ephemeris. We re-analyse the whole set of the timing data available. Our results provide an erratum to the previous HU Aqr planetary models, indicating that the hypothesis for a third and fourth body in this system is uncertain. The dynamical stability criterion and a particular geometry of orbits rule out coplanar 2-planet configurations. A putative HU Aqr planetary system may be more complex, e.g., highly non-coplanar. Indeed, we found examples of 3-planet configurations with the middle planet in a retrograde orbit, which are stable for at least 1Gyr, and consistent with the observations. The (O-C) may be also driven by oscillations of the gravitational quadrupole moment of the secondary, as predicted by the Lanza et al. modification of the Applegate mechanism. Further systematic, long-term monitoring of HU Aqr is required to interpret the (O-C) residuals.

HR 8799 is a nearby star hosting at least four ~10 Jovian mass planets in wide orbits up to ~70au, detected through the direct, high-contrast infrared imaging. Large companions and debris disks reported interior to ~10au, and exterior to ~100au indicate massive protoplanetary disc in the past. The dynamical state of the HR 8799 system is not yet fully resolved, due to limited astrometric data covering tiny orbital arcs. We construct a new, orbital model of the HR 8799 system, assuming rapid migration of the planets after their formation in wider orbits. We found that the HR 8799 planets are likely involved in double Laplace resonance, 1e:2d:4c:8b MMR. Quasi-circular planetary orbits are coplanar with the stellar equator and inclined by ~25 degrees to the sky plane. This best-fit orbital configuration matches astrometry, debris disk models, and mass estimates from cooling models. The multiple MMR is stable for the age of the star ~160Myr, for at least 1 Gyr unless significant perturbations to the N-body dynamics are present. We predict four configurations with the fifth hypothetical innermost planet HR 8799f in ~9.7au, or ~7.5au orbit, extending the MMR chain to triple Laplace resonance 1f:2e:4d:8c:16b MMR or to the 1f:3e:6d:12c:24b MMR, respectively. Our findings may establish strong boundary conditions for the system formation and its early history.

In this work we carry out an analysis of the observed times of primary and secondary eclipses of the post-common envelope binary NSVS14256825. Recently, \cite{Almeida2013} proposed that two circumbinary companions orbit this short-period eclipsing binary, in order to explain observed variations in the timing of mutual eclipses between the two binary components. Using a standard weighted least-squares minimisation technique, we have extensively explored the topology of $\chi^2$ parameter space of a single planet model. We find the data set to be insufficient to reliably constrain a one-companion model. Various models, each with similar statistical significance, result in substantially different orbital architectures for the additional companion. No evidence is seen for a second companion of planetary nature. We suspect insufficient coverage (baseline) of timing data causing the best-fit parameters to be unconstrained.

Red Dwarfs (RDs) are the most common, low-mass stars (<=0.5 M_sun) in the Solar neighbourhood, and probably in the Universe as well. Most (likely all) young RDs are magnetically active, and therefore it is impossible to measure their radial velocities (RVs) with the accuracy required by contemporary planetary searches (5-50 m/s). Some of RDs are known as a source of variable radio emission at centimetre wavelengths. This radiation is relatively weak (~0.2 to 1 mJy). However, recent development of the EVN and e-VLBI systems in terms of sensitivity gives a new possibility to investigate such stars. Here, we present first results from our ongoing RISARD project (Radio Interferometric Survey of Active Red Dwarfs). The main goal of this project is to detect an exoplanet by direct, precise measurements of a RD position and possible changes to this position caused by the planet.

We report a linear ordering of orbits in a sample of multiple extrasolar planetary systems with super-Earth planets. We selected 20 cases, mostly discovered by the Kepler mission, hosting at least four planets within \sim 0.5 au. The semi-major axis a_n of an n-th planet in each system of this sample obeys a(n) = a_1 + (n-1) \Delta a, where a_1 is the semi-major axis of the innermost orbit and \Delta a is a spacing between subsequent planets, which are specific for a particular system. For instance, the Kepler-33 system hosting five super-Earth planets exhibits the relative deviations between the observed and linearly predicted semi-major axes of only a few percent. At least half of systems in the sample fulfill the linear law with a similar accuracy. We explain the linear distribution of semi-major axes as a natural implication of multiple chains of mean motion resonances between subsequent planets, which emerge due to planet--disk interactions and convergent migration at early stages of their evolution.

The transiting planet WASP-12 b was identified as a potential target for transit timing studies because a departure from a linear ephemeris was reported in the literature. Such deviations could be caused by an additional planet in the system. We attempt to confirm the existence of claimed variations in transit timing and interpret its origin. We organised a multi-site campaign to observe transits by WASP-12 b in three observing seasons, using 0.5-2.6-metre telescopes. We obtained 61 transit light curves, many of them with sub-millimagnitude precision. The simultaneous analysis of the best-quality datasets allowed us to obtain refined system parameters, which agree with values reported in previous studies. The residuals versus a linear ephemeris reveal a possible periodic signal that may be approximated by a sinusoid with an amplitude of 0.00068+/-0.00013 d and period of 500+/-20 orbital periods of WASP-12 b. The joint analysis of timing data and published radial velocity measurements results in a two-planet model which better explains observations than single-planet scenarios. We hypothesize that WASP-12 b might be not the only planet in the system and there might be the additional 0.1 M_Jup body on a 3.6-d eccentric orbit. A dynamical analysis indicates that the proposed two-planet system is stable over long timescales.

We investigate the orbital stability of a putative Jovian planet in a compact binary \nu Octantis reported by Ramm et al. We re-analyzed published radial velocity data in terms of self-consistent Newtonian model and we found stable best-fit solutions that obey observational constraints. They correspond to retrograde orbits, in accord with an earlier hypothesis of Eberle & Cuntz, with apsidal lines anti-aligned with the apses of the binary. The best-fit solutions are confined to tiny stable regions of the phase space. These regions have a structure of the Arnold web formed by overlapping low-order mean motion resonances and their sub-resonances. The presence of a real planet is still questionable, because its formation would be hindered by strong dynamical perturbations. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework MECHANIC developed to run massive numerical experiments on CPU clusters.

The Kepler-11 star hosts at least six transiting super-Earth planets detected through the precise photometric observations of the Kepler mission (Lissauer et al.). In this paper, we re-analyze the available Kepler data, using the direct N-body approach rather than an indirect TTV method in the discovery paper. The orbital modeling in the realm of the direct approach relies on the whole data set, not only on the mid-transits times. Most of the results in the original paper are confirmed and extended. We constrained the mass of the outermost planet g to less than 30 Earth masses. The mutual inclinations between orbits b and c as well as between orbits d and e are determined with a good precision, in the range of [1,5] degrees. Having several solutions to four qualitative orbital models of the Kepler-11 system, we analyze its global dynamics with the help of dynamical maps. They reveal a sophisticated structure of the phase space, with narrow regions of regular motion. The dynamics are governed by a dense net of three- and four-body mean motion resonances, forming the Arnold web. Overlapping of these resonances is a main source of instability. We found that the Kepler-11 system may be long-term stable only in particular multiple resonant configurations with small relative inclinations. The mass-radius data derived for all companions reveal a clear anti-correlation between the mean density of the planets with their distance from the star. This may reflect the formation and early evolution history of the system.

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 develop the Mechanic package, which is a new numerical framework for dynamical astronomy. The aim of our software is to help in massive numerical simulations by efficient task management and unified data storage. The code is built on top of the Message Passing Interface (MPI) and Hierarchical Data Format (HDF5) standards and uses the Task Farm approach to manage numerical tasks. It relies on the core-module approach. The numerical problem implemented in the user-supplied module is separated from the host code (core). The core is designed to handle basic setup, data storage and communication between nodes in a computing pool. It has been tested on large CPU-clusters, as well as desktop computers. The Mechanic may be used in computing dynamical maps, data optimization or numerical integration. The code and sample modules are freely available at http://git.astri.umk.pl/projects/mechanic.

In a recent study, Lee et al. presented new photometric follow-up timing observations of the semi-detached binary system SZ Herculis and proposed the existence of two hierarchical cirumbinary companions. Based on the light-travel time effect, the two low-mass M-dwarf companions are found to orbit the binary pair on moderate to high eccentric orbits. The derived periods of these two companions are close to a 2:1 mean-motion orbital resonance. We have studied the stability of the system using the osculating orbital elements as presented by Lee et al. Results indicate an orbit-crossing architecture exhibiting short-term dynamical instabilities leading to the escape of one of the proposed companions. We have examined the system's underlying model parameter-space by following a Monte Carlo approach and found an improved fit to the timing data. A study of the stability of our best-fitting orbits also indicates that the proposed system is generally unstable. If the observed anomalous timing variations of the binary period is due to additional circumbinary companions, then the resulting system should exhibit a long-term stable orbital configuration much different from the orbits suggested by Lee et al. We, therefore, suggest that based on Newtonian-dynamical considerations, the proposed quadruple system cannot exist. To uncover the true nature of the observed period variations of this system, we recommend future photometric follow-up observations that could further constrain eclipse-timing variations and/or refine light-travel time models.

In this work we propose a new orbital architecture for the two proposed circumbinary planets around the polar eclipsing binary HU Aquarii. We base the new two-planet, light-travel time model on the result of a Monte Carlo simulation driving a least-squares Levenberg-Marquardt minimisation algorithm on the observed eclipse egress times. Our best-fitting model with $\chi_{r}^2=1.43$ resulted in high final eccentricities for the two companions leading to an unstable orbital configuration. From a large ensemble of initial guesses we examined the distribution of final eccentricities and semi-major axes for different $\chi_{r}^2$ parameter intervals and encountered qualitatively a second population of best-fitting parameters. The main characteristic of this population is described by low-eccentric orbits favouring long-term orbital stability of the system. We present our best-fitting model candidate for the proposed two-planet system and demonstrate orbital stability over one million years using numerical integrations.

We investigate the non-resonant, 3-D (spatial) model of the hierarchical system composed of point-mass stellar (or sub-stellar) binary and a low-mass companion (a circumbinary planet or a brown dwarf). We take into account the leading relativistic corrections to the Newtonian gravity. The secular model of the system relies on the expansion of the perturbing Hamiltonian in terms of the ratio of semi-major axes $\alpha$, averaged over the mean anomalies. We found that the low-mass object in a distant orbit may excite large eccentricity of the inner binary when the mutual inclination of the orbits is larger than about of 60 deg. This is related to strong instability caused by a phenomenon which acts similarly to the Lidov-Kozai resonance (LKR). The secular system may be strongly chaotic and its dynamics unpredictable over the long-term time scale. Our study shows that in the Jupiter-- or brown dwarf-- mass regime of the low-massive companion, the restricted model does not properly describe the long-term dynamics in the vicinity of the LKR. The relativistic correction is significant for the parametric structure of a few families of stationary solutions in this problem, in particular, for the direct orbits configurations (with the mutual inclination less than 90 degrees). We found that the dynamics of hierarchical systems with small $\alpha \sim 0.01$ may be qualitatively different in the realm of the Newtonian (classic) and relativistic models. This holds true even for relatively large masses of the secondaries.

More than 450 extrasolar planets are known to date. To detect these intriguing objects, many photometric and radial velocity (RV) surveys are in progress. We developed the Keplerian FITting (KFIT) code, to model published and available RV data. KFIT is based on a hybrid, quasi-global optimization technique relying on the Genetic Algorithms and simplex algorithm. Here, we re-analyse the RV data of evolved K3III star HD 240210. We confirm three equally good solutions which might be interpreted as signals of 2-planet systems. Remarkably, one of these best-fits describes the long-term stable two-planet system, involved in the 2:1~mean motion resonance (MMR). It may be the first instance of this strong MMR in a multi-planet system hosted by evolved star, as the 2:1 MMR configurations are already found around a few solar dwarfs.

The detailed study of the exoplanetary systems HD189733 and HD209458 has given rise to a wealth of exciting information on the physics of exoplanetary atmospheres. To further our understanding of the make-up and processes within these atmospheres we require a larger sample of bright transiting planets. We have began a project to detect more bright transiting planets in the southern hemisphere by utilising precision radial-velocity measurements. We have observed a constrained sample of bright, inactive and metal-rich stars using the HARPS instrument and here we present the current status of this project, along with our first discoveries which include a brown dwarf/extreme-Jovian exoplanet found in the brown dwarf desert region around the star HD191760 and improved orbits for three other exoplanetary systems HD48265, HD143361 and HD154672. Finally, we briefly discuss the future of this project and the current prospects we have for discovering more bright transiting planets.

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 report the discovery of a brown dwarf on an eccentric orbit and with a semimajor axis that places it in the brown dwarf desert region around the star HD191760. The star has a spectral type of G3IV/V and a metallicity ([Fe/H]) of 0.29 dex. HD191760 adds to the small number of metal-rich stars with brown dwarf companions. The brown dwarf (HD191760b) is found to have an orbital period of 505.57+/-0.40 days and semimajor axis of 1.35+/-0.01 AU, placing it firmly in the brown dwarf desert. The eccentricity of HD191760b is found to be 0.63+/-0.01, meaning it reaches as close as 0.5 AU from the host star. Dynamical simulations indicate that no inner planets could reside at separations beyond ~0.17 AU due to the disastrous gravity imposed by HD191760b. In addition to these first results we also refine the orbits found for the exoplanets around the stars HD48265, HD143361 and HD154672. All 1-planet solutions are in agreement with those previously published by the Magellan Planet Search.

The recent discovery of a three-planet extrasolar system of HR 8799 by Marois et al. is a breakthrough in the field of the direct imaging. This great achievement raises questions on the formation and dynamical stability of the HR 8799 system, because Keplerian fits to astrometric data are strongly unstable during ~0.2Myr. We search for stable, self-consistent N-body orbits with the so called GAMP method that incorporates stability constraints into the optimization algorithm. Our searches reveal only small regions of stable motions in the phase space of three-planet, coplanar configurations. Most likely, if the planetary masses are in 10-Jupiter-mass range, they may be stable only if the planets are involved in two- or three-body mean motion resonances (MMRs). We found that 80% systems found by GAMP that survived 30Myr backwards integrations, eventually become unstable after 100Myr. It could mean that the HR 8799 system undergo a phase of planet-planet scattering. We test a hypothesis that the less certain detection of the innermost object is due to a blending effect. In such a case, two-planet best-fit systems are mostly stable, on quasi-circular orbits and close to the 5:2 MMR, resembling the Jupiter-Saturn pair.

We investigate the secular dynamics of a planetary system composed of the parent star and two massive planets in mutually inclined orbits. The dynamics are investigated in wide ranges of semi-major axes ratios (0.1-0.667), and planetary masses ratios (0.25-2) as well as in the whole permitted ranges of the energy and total angular momentum. The secular model is constructed by semi-analytic averaging of the three-body system. We focus on equilibria of the secular Hamiltonian (periodic solutions of the full system), and we analyze their stability. We attempt to classify families of these solutions in terms of the angular momentum integral. We identified new equilibria, yet unknown in the literature. Our results are general and may be applied to a wide class of three-body systems, including configurations with a star and brown dwarfs and sub-stellar objects. We also describe some technical aspects of the semi-numerical averaging. The HD 12661 planetary system is investigated as an example configuration.

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analog of the Solar system. Our work is inspired by work of Malhotra and Minton (2008). Using the linear Laplace-Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the sub-system of Jupiters may be found in dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semi-major axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.

We consider the secular dynamics of a binary and a planet in terms of non-restricted, hierarchical three-body problem, including the general relativity corrections to the Newtonian gravity. We determine regions in the parameter space where the relativistic corrections may be important for the long-term dynamics. We try to constrain the inclinations of putative Jovian planets in recently announced binary systems of HD 4113 and HD 156846.

We construct a secular theory of a coplanar system of N-planets not involved in strong mean motion resonances, and which are far from collision zones. Besides the point-to-point Newtonian mutual interactions, we consider the general relativity corrections to the gravitational potential of the star and the innermost planet, and also a modification of this potential by the quadrupole moment and tidal distortion of the star. We focus on hierarchical planetary systems. A survey regarding model parameters (the masses, semi-major axes, spin rate of the star) reveals a rich and non-trivial dynamics of the secular system. Our study is focused on its equilibria. Such solutions predicted by the classic secular theory, which correspond to aligned (mode I) or anti-aligned (mode II) apsides, may be strongly affected by the gravitational corrections. The so called true secular resonance, which is a new feature of the classic two-planet problem discovered by Michtchenko & Malhotra (2004), may appear in other, different regions of the phase space of the generalized model. We found bifurcations of mode II which emerge new, yet unknown in the literature, secularly unstable equilibria and a complex structure of the phase space. These equilibria may imply secularly unstable orbital configurations even for nitially moderate eccentricities. The point mass gravity corrections can affect the long term-stability in the secular time scale, which may directly depend on the age of the host star through its spin rate. We also analyze the secular dynamics of the upsilon Andromede system in the realm of the generalized model. Also in this case of the three-planet system, new secular equilibria may appear.

We present the secular theory of coplanar $N$-planet system, in the absence of mean motion resonances between the planets. This theory relies on the averaging of a perturbation to the two-body problem over the mean longitudes. We expand the perturbing Hamiltonian in Taylor series with respect to the ratios of semi-major axes which are considered as small parameters, without direct restrictions on the eccentricities. Next, we average out the resulting series term by term. This is possible thanks to a particular but in fact quite elementary choice of the integration variables. It makes it possible to avoid Fourier expansions of the perturbing Hamiltonian. We derive high order expansions of the averaged secular Hamiltonian (here, up to the order of 24) with respect to the semi-major axes ratio. The resulting secular theory is a generalization of the octupole theory. The analytical results are compared with the results of numerical (i.e., practically exact) averaging. We estimate the convergence radius of the derived expansions, and we propose a further improvement of the algorithm. As a particular application of the method, we consider the secular dynamics of three-planet coplanar system. We focus on stationary solutions in the HD 37124 planetary system.

Most of the first-discovered extrasolar multi-planet systems were found to lie close to dynamically unstable configurations. However a few observed multi-planet systems (e.g. HD 74156) did not show this trait. Those systems could share this property if they contain an additional planet in between those that are known. Previous investigations identified the properties of hypothetical planets that would place these systems near instability. The hypothetical planet in HD 74156 was expected to have a mass about equal to that of Saturn, a semi-major axis between 0.9 and 1.4 AU, and an eccentricity less than 0.2. HD 74156 d, a planet with a mass of 1.3 Saturn masses at 1.04 AU with an eccentricity of 0.25, was recently announced. We have reanalyzed all published data on this system in order to place tighter constraints on the properties of the new planet. We find two possible orbits for this planet, one close to that already identified and another (with a slightly better fit to the data) at ~0.89 AU. We also review the current status of other planet predictions, discuss the observed single planet systems, and suggest other systems which may contain planets in between those that are already known. The confirmation of the existence of HD 74156 d suggests that planet formation is an efficient process, and planetary systems should typically contain many planets.

We present an analysis of high precision radial velocity (RV) observations of stars hosting multi-planet systems with Jovian companions. We use dynamical stability constraints and quasi-global methods of optimization. As an illustration, we present new results derived for the RV data of the Sun-like dwarfs HD 155358 and $\tau^1$ Gruis.

Precision radial velocity (RV) measurements of the Sun-like dwarf 14 Herculis published by Naef et. al (2004), Butler et. al (2006) and Wittenmyer et al (2007) reveal a Jovian planet in a 1760 day orbit and a trend indicating the second distant object. On the grounds of dynamical considerations, we test a hypothesis that the trend can be explained by the presence of an additional giant planet. We derive dynamical limits to th orbital parameters of the putative outer Jovian companion in an orbit within ~13 AU. In this case, the mutual interactions between the Jovian planets are important for the long-term stability of the system. The best self-consistent and stable Newtonian fit to an edge-on configuration of Jovian planets has the outer planet in 9 AU orbit with a moderate eccentricity ~0.2 and confined to a zone spanned by the low-order mean motion resonances 5:1 and 6:1. This solution lies in a shallow minimum of \Chi and persists over a wide range of the system inclination. Other stable configurations within 1\sigma confidence interval of the best fit are possible for the semi-major axis of the outer planet in the range of (6,13) AU and the eccentricity in the range of (0,0.3). The orbital inclination cannot yet be determined but when it decreases, both planetary masses approach ~10 Jupiter masses and for ~30 deg the hierarchy of the masses is reversed.

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best fit solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best fit solutions are found in dynamically active region of the phase space. The long term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myrs of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behavior.

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.

Two nearby stars, HD 128311 and HD 82943, are believed to host pairs of Jupiter-like planets involved in a strong first order 2:1 mean motion resonance (MMR). In this work we reanalyze available radial velocity (RV) measurements and demonstrate that it is also possible to explain the observed RV variations of the parents stars as being induced by a pair of Trojan planets (i.e., in a 1:1 MMR). We show that these Trojan configurations reside in extended zones of stability in which such systems may easily survive in spite of large masses of the planets, large eccentricities and nonzero mutual inclinations of their orbits. We also show that HD 82943 could harbor a previously unknown third planet about of 0.5 Jupiter masses in 2.1 AU orbit.

The kinematic structure of a sample of planetary nebulae, consisting of 23 [WR] central stars, 21 weak emission line stars (wels) and 57 non-emission line central stars, is studied. The [WR] stars are shown to be surrounded by turbulent nebulae, a characteristic shared by some wels but almost completely absent from the non-emission line stars. The fraction of objects showing turbulence for non-emission-line stars, wels and [WR] stars is 7%, 24% and 91%, respectively. The [WR] stars show a distinct IRAS 12-micron excess, indicative of small dust grains, which is not found for wels. The [WR]-star nebulae are on average more centrally condensed than those of other stars. On the age-temperature diagram, the wels are located on tracks of both high and low stellar mass, while [WR] stars trace a narrow range of intermediate masses. Emission-line stars are not found on the cooling track. One group of wels may form a sequence wels--[WO] stars with increasing temperature. For the other groups both the wels and the [WR] stars appear to represent several, independent evolutionary tracks. We find a discontinuity in the [WR] stellar temperature distribution and suggest different evolutionary sequences above and below the temperature gap. One group of cool [WR] stars has no counterpart among any other group of PNe and may represent binary evolution. A prime factor distinguishing wels and [WR] stars appears to be stellar luminosity. We find no evidence for an increase of nebular expansion velocity with time.

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 re-analyze the precision radial velocity (RV) data of HD188015, HD114729, HD190360, HD147513 and HD208487. All these stars are supposed to host Jovian companions in long-period orbits. We test a hypothesis that the residuals of the 1-planet model of the RV or an irregular scatter of the measurements about the synthetic RV curve may be explained by the existence of additional planets in short-period orbits. We perform a global search for the best fits in the orbital parameters space with genetic algorithms and simplex method. This makes it possible to verify and extend the results obtained with an application of commonly used FFT-based periodogram analysis for identifying the leading periods. Our analysis confirms the presence of a periodic component in the RV data of HD190360 which may correspond to a hot-Neptune planet. We found four new cases when the 2-planet model yields significantly better fits to the RV data than the best 1-planet solutions. If the periodic variability of the residuals of single-planet fits has indeed a planetary origin then hot-Neptune planets may exist in these extrasolar systems. We estimate their orbital periods in the range of 7-20d and minimal masses about of 20 masses od the Earth.

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.

We study the long-term dynamics of the PSR 1257+12 planetary system. Using the recently determined accurate initial condition by Konacki & Wolszczan (2003) who derived the orbital inclinations and the absolute masses of the planets B and C, we investigate the system stability by long-term, 1Gyr direct integrations. No secular changes of the semi-major axes eccentricities and inclinations appear during such an interval. This stable behavior is confirmed with the fast indicator MEGNO. The analysis of the orbital stability in the neighborhood of the nominal initial condition reveals that the PSR 1257+12 system is localized in a wide stable region of the phase space but close to a few weak 2 and 3-body mean motion resonances. The long term stability is additionally confirmed by a negligible exchange of the Angular Momentum Deficit between the innermost planet A and the pair of the outer planets B and C. An important feature of the system that helps sustain the stability is the secular apsidal resonance (SAR) between the planets B and C with the center of libration about 180deg. We find useful limits on the elements of the innermost planet A which are otherwise unconstrained by the observations. We also perform a preliminary study of the short-term dynamics of massless particles in the system.

The new 2-planetary system around HD169830 has been announced during the XIX-th IAP Colloquium "Extrasolar Planets: Today & Tomorrow" (Paris, June 30 - July 4, 2003) by the Geneva Extrasolar Planet Search team. We study the orbital dynamics of this system in the framework of the $N$-body problem. The analysis of its orbital stability is performed using the long-term integrations and the fast indicators, the Mean Exponential Growth factor of Nearby Orbits and the Frequency Map Analysis. The HD169830 appears to be located in a wide stable region of the phase space. The ratio of the mean motions of the planets HD169830b and c is between low-order mean motion resonances, 9:1 and 10:1. The long-term integration of the coplanar configurations, conducted over 1Gyr, reveals that the eccentricities of the companions vary with a large amplitude about 0.4-0.5 but there is no sign of instability. The orbital parameters of the planets resemble those of another 2-planetary system, around HD12661. Both of them can be classified as hierarchical planetary systems. We investigate whether these two exosystems are dynamically similar. Such similarities may be important for finding out if the formation and subsequent orbital evolution of exoplanetary systems obey common rules.

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)