
We study proximity of the Kepler25 planetary system to a periodic
configuration, which is known to be the final state of a system that undergoes
smooth migration resulting from the planetdisc 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 antialigned apsidal lines) may be
misinterpreted as an aligned nonresonant system. We demonstrate that the
standard MCMC modelling of the Kepler25 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 nonresonant configurations with small planets' masses
and moderate eccentricities, while for the narrower priors the most likely are
the antialigned resonant systems with larger masses and very low
eccentricities.

We study the orbital architecture, physical characteristics of planets,
formation and longterm evolution of the Kepler30 planetary system, detected
and announced in 2012 by the KEPLER team. We show that the Kepler30 system
belongs to a particular class of very compact and quasiresonant, yet longterm
stable planetary systems. We reanalyse the light curves of the host star
spanning Q1Q17 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 Neptunelike 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 reanalysis 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 Kepler30b, Kepler30c and Kepler30d, respectively. The
masses of the inner planets are determined within $\sim 1\%$ uncertainty. We
infer the internal structures of the Kepler30 planets and their bulk densities
in a wide range from $ (0.19 \pm 0.01)$ g$\cdot$cm$^{3}$ for Kepler30d,
$(0.96\pm0.15)$ g$\cdot$cm$^{3}$ for Kepler30b, to $(1.71 \pm 0.13)$
g$\cdot$cm$^{3}$ for the Jovian planet Kepler30c. 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
Jupiterlike planet plays an important role in determining the present
dynamical state of this system.

The cyclic behaviour of (OC) residuals of eclipse timings in the sdB+M
eclipsing binary NSVS 14256825 was previously attributed to one or two
Joviantype 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 (OC) diagram by three years. Based on the archival and our new data
spanning over more than 17 years we reexamined the up to date system (OC).
The data revealed systematic, quasisinusoidal variation deviating from an
older linear ephemeris by about 100 s. It also exhibits a maximum in the (OC)
near JD 2,456,400 that was previously unknown. We consider two most credible
explanations of the (OC) 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 LanzaRodono effect. We found that the latter mechanism is unlikely due to
the insufficient energy budget of the Mdwarf secondary. In the framework of
the thirdbody hypothesis, we obtained meaningful constraints on the Keplerian
parameters of a putative companion and its mass. Our bestfitting model
indicates that the observed quasiperiodic (OC) 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 multiepoch astrometric observations of the M9 ultracool dwarf
TVLM51346546 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 180300 muJy. We identified four
shortduration flaring events (0.52 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
TVLM51346546 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
superEarths detected as Kepler29, in order to constrain the planets' masses
and orbital parameters. A dynamical analysis of the bestfitting 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 antialigned 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 Kepler29 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 Kepler60 planetary system with
three superEarths. We first determine their orbital elements and masses by
Transit Timing Variation (TTV) data spanning quarters Q1Q16 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 threebody 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 twoplanet 5:4
and 4:3 MMRs, all critical angles librate with small amplitudes of ~30 degrees
and apsidal lines in planet's pairs are antialigned. The system is
simultaneously locked in a threebody MMR with librations amplitude of ~10
degrees. The true Laplace MMR can evolve towards a chain of twobody MMRs in
the presence of planetary migration. Therefore the threebody MMR formed in
this way seems to be more likely state of the system. However, the true
threebody 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
Kepler47 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 quasiperiodic orbits, and assessed the longterm stability of the
threeplanet 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 Kepler47c 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 Kepler47
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
reanalysis of the entire currently available data including those obtained
postannouncement 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, highcontrast 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. Quasicircular planetary
orbits are coplanar with the stellar equator and inclined by ~25 degrees to the
sky plane. This bestfit 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 Nbody 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 opensource 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 wellestablished
Message Passing Interface (MPI) standard, which is widely used on variety of
parallel computers and CPUclusters. The data storage is performed within the
Hierarchical Data Format (HDF5). The design of the code follows em coremodule
approach which allows to reduce the user's codebase and makes it portable for
single and multiCPU 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 longterm orbital evolution of planetary systems
with Monte Carlo methods, dynamical maps or evolutionary algorithms. It has
been already applied in numerical experiments conducted for Kepler11
(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 postcommon envelope binary NSVS14256825. Recently,
\cite{Almeida2013} proposed that two circumbinary companions orbit this
shortperiod eclipsing binary, in order to explain observed variations in the
timing of mutual eclipses between the two binary components. Using a standard
weighted leastsquares 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 onecompanion 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 bestfit parameters to be
unconstrained.

Red Dwarfs (RDs) are the most common, lowmass 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 (550 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 eVLBI 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 superEarth planets. We selected 20 cases, mostly
discovered by the Kepler mission, hosting at least four planets within \sim 0.5
au. The semimajor axis a_n of an nth planet in each system of this sample
obeys a(n) = a_1 + (n1) \Delta a, where a_1 is the semimajor axis of the
innermost orbit and \Delta a is a spacing between subsequent planets, which are
specific for a particular system. For instance, the Kepler33 system hosting
five superEarth planets exhibits the relative deviations between the observed
and linearly predicted semimajor 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 semimajor axes as a natural implication of
multiple chains of mean motion resonances between subsequent planets, which
emerge due to planetdisk 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 reanalyzed published radial
velocity data in terms of selfconsistent Newtonian model and we found stable
bestfit solutions that obey observational constraints. They correspond to
retrograde orbits, in accord with an earlier hypothesis of Eberle & Cuntz, with
apsidal lines antialigned with the apses of the binary. The bestfit solutions
are confined to tiny stable regions of the phase space. These regions have a
structure of the Arnold web formed by overlapping loworder mean motion
resonances and their subresonances. 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 Kepler11 star hosts at least six transiting superEarth planets detected
through the precise photometric observations of the Kepler mission (Lissauer et
al.). In this paper, we reanalyze the available Kepler data, using the direct
Nbody 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 midtransits 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 Kepler11 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 fourbody mean motion resonances, forming the Arnold
web. Overlapping of these resonances is a main source of instability. We found
that the Kepler11 system may be longterm stable only in particular multiple
resonant configurations with small relative inclinations. The massradius data
derived for all companions reveal a clear anticorrelation 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 (OC) 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
reanalyse the whole currently available data set. We add almost 60 new, yet
unpublished, mostly precision light curves obtained using the time
highresolution photopolarimeter OPTIMA, as well as photometric observations
performed at the MONET/N, PIRATE and TCS telescopes. We determine new
midegress 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 nonhomogeneous with respect to spectral windows (ultraviolet, Xray,
visual, polarimetric mode) and the reported midegress 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 bestfit 2planet configurations,
including selfconsistent, Newtonian Nbody solutions may be able to explain
the data. However, using high resolution, precision OPTIMA light curves, we
find that the (OC) 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 Jupitermasses. 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 Kepler16, Kepler34 and Kepler35.

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 coremodule approach. The numerical problem implemented
in the usersupplied 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 CPUclusters, 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 twoplanet, lighttravel time model on the result of a Monte Carlo
simulation driving a leastsquares LevenbergMarquardt minimisation algorithm
on the observed eclipse egress times. Our bestfitting 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 semimajor
axes for different $\chi_{r}^2$ parameter intervals and encountered
qualitatively a second population of bestfitting parameters. The main
characteristic of this population is described by loweccentric orbits
favouring longterm orbital stability of the system. We present our
bestfitting model candidate for the proposed twoplanet system and demonstrate
orbital stability over one million years using numerical integrations.

We investigate the nonresonant, 3D (spatial) model of the hierarchical
system composed of pointmass stellar (or substellar) binary and a lowmass
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 semimajor axes $\alpha$, averaged over the mean anomalies. We
found that the lowmass 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 LidovKozai resonance (LKR). The secular system may
be strongly chaotic and its dynamics unpredictable over the longterm time
scale. Our study shows that in the Jupiter or brown dwarf mass regime of
the lowmassive companion, the restricted model does not properly describe the
longterm 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 nonlinear, 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 coplanar", and "conic" precessions,
respectively. Moreover, we also show that in the case of conic precession,
although the relative orbit is circular, the pointmass 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 metalrich 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 1planet
solutions are in agreement with those previously published by the Magellan
Planet Search.

The recent discovery of a threeplanet 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, selfconsistent Nbody 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 threeplanet, coplanar configurations. Most
likely, if the planetary masses are in 10Jupitermass range, they may be
stable only if the planets are involved in two or threebody 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 planetplanet scattering. We test a
hypothesis that the less certain detection of the innermost object is due to a
blending effect. In such a case, twoplanet bestfit systems are mostly stable,
on quasicircular orbits and close to the 5:2 MMR, resembling the
JupiterSaturn 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 semimajor axes ratios (0.10.667), and
planetary masses ratios (0.252) as well as in the whole permitted ranges of
the energy and total angular momentum. The secular model is constructed by
semianalytic averaging of the threebody 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
threebody systems, including configurations with a star and brown dwarfs and
substellar objects. We also describe some technical aspects of the
seminumerical averaging. The HD 12661 planetary system is investigated as an
example configuration.

We investigate the dynamics of putative Earthmass planets in the habitable
zone (HZ) of the extrasolar planetary system OGLE2006BLG109L, a close analog
of the Solar system. Our work is inspired by work of Malhotra and Minton
(2008). Using the linear LaplaceLagrange 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 subsystem
of Jupiters may be found in dynamically active region of the phase space
spanned by loworder meanmotion resonances. To generalize this secular model,
we construct a semianalytical 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 highorder analytic secular theory developed
in our previous works. We found regions in the parameter space of the problem
in which stable, quasicircular 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 semimajor axes
consistent with the observations, a terrestrial planet may survive in low
eccentric orbits. We also study the effect of postNewtonian gravity correction
on the innermost secular resonance.

We consider the secular dynamics of a binary and a planet in terms of
nonrestricted, hierarchical threebody 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
longterm 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 Nplanets not involved
in strong mean motion resonances, and which are far from collision zones.
Besides the pointtopoint 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, semimajor
axes, spin rate of the star) reveals a rich and nontrivial 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 antialigned (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 twoplanet 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 termstability 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 threeplanet
system, new secular equilibria may appear.