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 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 conducted radio-interferometric observations of the well known binary cataclysmic system AM Herculis. This particular system is formed by the magnetic white dwarf (primary) and the red dwarf (secondary), and is the prototype of so-called polars. Our observations were conducted with the European VLBI Network (EVN) in the e-EVN mode at 5 GHz. We obtained six astrometric measurements spanning one year, which make it possible to update the annual parallax for this system with the best precision to date ($\pi=11.29\pm0.08$ mas), equivalent to the distance of $88.6\pm0.6$ pc. The system was observed mostly in the quiescent phase (visual magnitude $m_v\sim15.3$), when the radio emission was at the level of about $300~\mu$Jy. Our analysis suggests that the radio flux of AM Herculis is modulated with the orbital motion. Such specific properties of the radiation can be explained using the emission mechanism similar to the scenario proposed for V471 Tauri and, in general, for RS CVn type stars. In this scenario the radio emission arises near the surface of the red dwarf, where the global magnetic field strength may reach a few kG. We argue that the quiescent radio emission distinguishes AM Herculis together with AR Ursae Majoris (the second known persistent radio polar) from other polars, as the systems with a magnetized secondary star.

We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the time reversibility and the displacement from the initial condition, occurring when we integrate it forward and backward for the same time interval, is related to the dynamical character of the trajectory. If the motion is chaotic, in the sense of non-zero maximal Characteristic Lyapunov Exponent (mLCE), then REM increases exponentially with time, as exp \lambda t, while when the motion is regular (quasi-periodic) then REM increases as a power law in time, as t^{\alpha}, where {\alpha} and {\lambda} are real coefficients. We compare the REM with a variant of mLCE, the Mean Exponential Growth factor of Nearby Orbits (MEGNO). The test set includes the restricted three body problem and five resonant planetary systems: HD 37124, Kepler-60, Kepler-36, Kepler-29 and Kepler-26. We found a very good agreement between the outcomes of these algorithms. Moreover, the numerical implementation of REM is astonishing simple, and is based on solid theoretical background. The REM requires only a symplectic and time-reversible (symmetric) integrator of the equations of motion. This method is also CPU efficient. It may be particularly useful for the dynamical analysis of multiple planetary systems in the KEPLER sample, characterized by low-eccentricity orbits and relatively weak mutual interactions. As an interesting side-result, we found a possible stable chaos occurrence in the Kepler-29 planetary system.

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.

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 focus on short--period eclipsing binaries that belong to a class of Cataclysmic Variables (CVs). They are known as polars and intermediate polars, closely resembling their prototype AM Herculis. These binaries consist of a red dwarf and a strongly magnetic white dwarf, having orbital periods of only a few hours. Monitoring eclipses of these typically faint sources demands high-time resolution photometry. We describe the very recent results obtained for two CVs, HU Aqr and DQ Her, which were observed with the Optical Pulsar Timing Analyzer (OPTIMA). The new observations of HU Aqr confirm that the O--C (Observed minus Calculated) diagrams exhibit variations known for this binary which can be explained by a single, massive Jupiter--like planet, possibly accompanied by a very distant companion.

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.