From its discovery, the WASP-18 system with its massive transiting planet on a tight orbit was identified as a unique laboratory for studies on tidal planet-star interactions. In an analysis of Doppler data, which include five new measurements obtained with the HIRES/Keck-I instrument between 2012 and 2018, we show that the radial velocity signal of the photosphere following the planetary tidal potential can be distilled for the host star. Its amplitude is in agreement with both theoretical predictions of the equilibrium tide approximation and an ellipsoidal modulation observed in an orbital phase curve. Assuming a circular orbit, we refine system parameters using photometric time series from TESS. With a new ground-based photometric observation, we extend the span of transit timing observations to 28 years in order to probe the rate of the orbital period shortening. Since we found no departure from a constant-period model, we conclude that the modified tidal quality parameter of the host star must be greater than 3.9x10^6 with 95% confidence. This result is in line with conclusions drawn from studies of the population of hot Jupiters, predicting that the efficiency of tidal dissipation is 1 or 2 orders of magnitude weaker. As the WASP-18 system is one of the prime candidates for detection of orbital decay, further timing observations are expected to push the boundaries of our knowledge on stellar interiors.

We validate the planetary origin of the KOI-1599 transit time variations (TTVs) with statistical and dynamical tests. We re-analysed KEPLER Q1-Q17 light-curves of the star, and we independently derived the TTVs. They appear as strongly anti-correlated, suggestive of two mutually interacting planets. We found similar radii of the candidates, $1.9 \pm 0.2\, \mbox{R}_{\oplus}$ for the inner KOI-1599.02, and $1.9 \pm 0.3 \, \mbox{R}_{\oplus}$ for the outer KOI-1599.01. The standard MCMC TTV analysis constrains the planet masses safely below the dynamical instability limit of $\simeq 3\, \mbox{M}_{\mbox{Jup}} $. The best-fitting MCMC model yields $(9.0\pm 0.3)\, \mbox{R}_{\oplus}$, and $(4.6\pm0.3),\ \mbox{R}_{\oplus}$, for the inner and the outer planet, respectively. The planets are trapped in 3:2 mean motion resonance (MMR) with anti-aligned apsides ($\Delta \varpi = 180^{\circ}$) at low-eccentric ($e\simeq 0.01)$ orbits. However, we found that the TTV mass determination depends on eccentricity priors with the dispersion in the (0.01,0.05) range. They permit a second family of TTV models with smaller masses of $\simeq 7\,\mbox{R}_{\oplus}$, and $\simeq 3.6\,\mbox{R}_{\oplus}$, respectively, exhibiting two modes of $\Delta \varpi = 0^{\circ},180^{\circ}$ librations. The 3:2 MMR is dynamically robust and persists for both modes. In order to resolve the mass duality, we re-analysed the TTV data with a quasi-analytic model of resonant TTV signals. This model favours the smaller masses. We also reproduced this model with simulating the migration capture of the system into the 3:2 MMR.

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.

We study the orbital architecture, physical characteristics of planets, formation and long-term evolution of the Kepler-30 planetary system, detected and announced in 2012 by the KEPLER team. We show that the Kepler-30 system belongs to a particular class of very compact and quasi-resonant, yet long-term stable planetary systems. We re-analyse the light curves of the host star spanning Q1-Q17 quarters of the KEPLER mission. A huge variability of the Transit Timing Variations (TTV) exceeding 2 days is induced by a massive Jovian planet located between two Neptune-like companions. The innermost pair is near to the 2:1 mean motion resonance (MMR), and the outermost pair is close to higher order MMRs, such as 17:7 and 7:3. Our re-analysis of photometric data allows us to constrain, better than before, the orbital elements, planets' radii and masses, which are $9.2 \pm 0.1$, $536 \pm 5$, and $23.7 \pm 1.3$ Earth masses for Kepler-30b, Kepler-30c and Kepler-30d, respectively. The masses of the inner planets are determined within $\sim 1\%$ uncertainty. We infer the internal structures of the Kepler-30 planets and their bulk densities in a wide range from $ (0.19 \pm 0.01)$ g$\cdot$cm$^{-3}$ for Kepler-30d, $(0.96\pm0.15)$ g$\cdot$cm$^{-3}$ for Kepler-30b, to $(1.71 \pm 0.13)$ g$\cdot$cm$^{-3}$ for the Jovian planet Kepler-30c. We attempt to explain the origin of this unique planetary system and a deviation of the orbits from exact MMRs through the planetary migration scenario. We anticipate that the Jupiter-like planet plays an important role in determining the present dynamical state of this system.

We study formation of 9:7 mean motion resonance (MMR) as a result of convergent migration of two planets embedded in a disc. Depending on the migration parameters, initial orbits and planets' masses the system may pass through the resonance or enter it (permanently or temporarily). We illustrate that a stable equilibrium of the averaged system (a periodic orbit of the N-body model) is surrounded by the region of permanent resonance capture, whose size depends on the migration parameters and the planets mass ratio. A system located inside this region tends towards the equilibrium (the capture is permanent), while a system located outside the region evolves away from the equilibrium and leaves the resonance. We verify recent results of Delisle et al. and Xu & Lai where they show that for $m_1 \lesssim m_2$ ($m_1, m_2$ are the inner and outer planets' masses) the equilibrium is unstable when the migration is added, so the capture cannot be permanent. We show that for particular migration parameters the situation may be reversed (the equilibria are unstable for $m_1 \gtrsim m_2$). We illustrate that 9:7 MMR consists of two modes separated with a separatrix. The inner one is centred at the equilibrium and the outer one has no equilibrium in its centre. A system located outside the region of stable capture evolves from the inner into the outer mode. The evolution occurs along families of periodic orbits of the averaged system, that play a crucial role in the dynamics after the resonance capture.

Surveys have revealed many multi-planet systems containing super-Earths and Neptunes in orbits of a few days to a few months. There is debate whether in situ assembly or inward migration is the dominant mechanism of the formation of such planetary systems. Simulations suggest that migration creates tightly packed systems with planets whose orbital periods may be expressed as ratios of small integers (resonances), often in a many-planet series (chain). In the hundreds of multi-planet systems of sub-Neptunes, more planet pairs are observed near resonances than would generally be expected, but no individual system has hitherto been identified that must have been formed by migration. Proximity to resonance enables the detection of planets perturbing each other. Here we report transit timing variations of the four planets in the Kepler-223 system, model these variations as resonant-angle librations, and compute the long-term stability of the resonant chain. The architecture of Kepler-223 is too finely tuned to have been formed by scattering, and our numerical simulations demonstrate that its properties are natural outcomes of the migration hypothesis. Similar systems could be destabilized by any of several mechanisms, contributing to the observed orbital-period distribution, where many planets are not in resonances. Planetesimal interactions in particular are thought to be responsible for establishing the current orbits of the four giant planets in the Solar System by disrupting a theoretical initial resonant chain similar to that observed in Kepler-223.

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.

We study the migration of three-planet systems in an irradiated 1+1D $\alpha$-disc with photoevaporation. We performed $2700$ simulations with various planets' masses and initial orbits. We found that most of the systems which ended up as compact configurations form chains of mean motion resonances (MMRs) of the first and higher orders. Most of the systems involved in chains of MMRs are periodic configurations. The period ratios of such system, though, are not necessarily close to exact commensurability. If a given system resides in a divergent migration zone in the disc, the period ratios increase and evolve along resonant divergent migration paths at ($P_2/P_1, P_3/P_2)$-diagram, where $P_1, P_2, P_3$ are the orbital periods of the first, second and third planet, respectively. The observed systems, though, do not lie on those paths. We show that an agreement between the synthetic and the observed systems distributions could be achieved if the orbital circularization was slower than it results from models of the planet-disc interactions. Therefore, we conclude that most of those systems unlikely formed as a result of divergent migration out of nominal chains of MMRs.

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 dynamics of a system of two super-Earths embedded in a protoplanetary disc. We build a simple model of an irradiated viscous disc and use analytical prescriptions for the planet-disc interactions which lead to migration. We show that depending on the disc parameters, planets' masses and their positions in the disc, the migration of each planet can be inward or outward and the migration of a two-planet system can be convergent (which may lead to formation of a first order mean motion resonance, MMR) or divergent (a system moves away from MMR). We performed $3500$ simulations of the migration of two-planet systems with various masses and initial orbits. Almost all of them end up as resonant configurations, although the period ratios may be very distant from the nominal values of a given MMR. We found that almost all the systems resulting from the migration are periodic configurations.

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.

We introduce the Mechanic, a new open-source code framework. It is designed to reduce the development effort of scientific applications by providing unified API (Application Programming Interface) for configuration, data storage and task management. The communication layer is based on the well-established Message Passing Interface (MPI) standard, which is widely used on variety of parallel computers and CPU-clusters. The data storage is performed within the Hierarchical Data Format (HDF5). The design of the code follows em core-module approach which allows to reduce the user's codebase and makes it portable for single- and multi-CPU environments. The framework may be used in a local user's environment, without administrative access to the cluster, under the PBS or Slurm job schedulers. It may become a helper tool for a wide range of astronomical applications, particularly focused on processing large data sets, such as dynamical studies of long-term orbital evolution of planetary systems with Monte Carlo methods, dynamical maps or evolutionary algorithms. It has been already applied in numerical experiments conducted for Kepler-11 (Migaszewski et al., 2012), and nuOctantis planetary systems (Go\'zdziewski et al., 2013). In this paper we describe the basics of the framework, including code listings for the implementation of a sample user's module. The code is illustrated on a model Hamiltonian introduced by (Froeschle et al., 2000) presenting the Arnold diffusion. The Arnold Web is shown with the help of the MEGNO (Mean Exponential Growth of Nearby Orbits) fast indicator (Go\'zdziewski et al., 2008a) applied onto symplectic SABA integrators family (Laskar and Robutel, 2001).

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.

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.

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.

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.

We study the long-term dynamics of a planetary system composed of a star and a planet. Both bodies are considered as extended, non-spherical, rotating objects. There are no assumptions made on the relative angles between the orbital angular momentum and the spin vectors of the bodies. Thus, we analyze full, spatial model of the planetary system. Both objects are assumed to be deformed due to their own rotations, as well as due to the mutual tidal interactions. The general relativity corrections are considered in terms of the post-Newtonian approximation. Besides the conservative contributions to the perturbing forces, there are also taken into account non-conservative effects, i.e., the dissipation of the mechanical energy. This dissipation is a result of the tidal perturbation on the velocity field in the internal zones with non-zero turbulent viscosity (convective zones). Our main goal is to derive the equations of the orbital motion as well as the equations governing time-evolution of the spin vectors (angular velocities). We derive the Lagrangian equations of the second kind for systems which do not conserve the mechanical energy. Next, the equations of motion are averaged out over all fast angles with respect to time-scales characteristic for conservative perturbations. The final equations of motion are then used to study the dynamics of the non-conservative model over time scales of the order of the age of the star. We analyze the final state of the system as a function of the initial conditions. Equilibria states of the averaged system are finally discussed.

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.

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 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 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 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 present a method of modeling the radial velocity (RV) measurements which can be useful in searching for planets hosted by chromospherically active stars. We assume that the observed RV signal is induced by the reflex motion of a star as well as by distortions of spectral line profiles, measured by the Bisector Velocity Span (BVS). The RVs are fitted with a common planetary model including RV correction term depending linearly on the BVS, which accounts for the stellar activity. The coefficient of correlation is an additional free parameter of the RV model. That approach differs from correcting the RVs before or after fitting the "pure" planetary model. We test the method on simulated data derived for single-planet systems. The results are compared with the outcomes of algorithms found in the literature.

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.

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 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 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.