• 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.
  • 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.
  • 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.
  • 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 have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system be the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by integrating the entire 5 bodies (binary + planets) for 10 Myr. We identified several stable regions between the two known planets as well as a region beyond the orbit of Kepler-47c where the orbit of the third planet could be stable. To constrain the orbit of this planet, we used the measured duration of the unexplained transit event ($\sim 4.15$ hours) and compared that with the transit duration of the third planet in an ensemble of stable orbits. To remove the degeneracy among the orbits with similar transit durations, we considered the planet to be in a circular orbit and calculated its period analytically. The latter places an upper limit of 424 days on the orbital period of the third planet. Our analysis suggests that if the unexplained transit event detected during the discovery of the Kepler-47 circumbinary system is due to a planetary object, this planet will be in a low eccentricity orbit with a semimajor axis smaller than 1.24 AU. Further constraining of the mass and orbital elements of this planet requires a re-analysis of the entire currently available data including those obtained post-announcement of the discovery of this system. We present details of our methodology and discuss the implication of the results.
  • 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).
  • 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.
  • 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 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.
  • 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.