We characterize the radii and masses of the star and planets in the Kepler-59 system, as well as their orbital parameters. The star parameters are determined through a standard spectroscopic analysis, resulting in a mass of $1.359\pm 0.155\,M_\odot$ and a radius of $1.367\pm 0.078\,R_\odot$. The planetary radii obtained are $1.5\pm 0.1\,R_\oplus$ for the inner and $2.2\pm 0.1\,R_\oplus$ for the outer planet. The orbital parameters and the planetary masses are determined by the inversion of Transit Timing Variations (TTV) signals. For this, we consider two different data sets, one provided by Holczer et al. 2016, with TTVs only for the planet Kepler-59c, and the other provided by Rowe et al. 2015, with TTVs signals for both planets. The inversion method is carried out by applying an algorithm of Bayesian inference (MultiNest) combined with an efficient N-body integrator (Swift). For each of the data sets, two possible solutions are found, both having the same probability according to their corresponding Bayesian evidences. All four solutions appear to be indistinguishable within their 2-$\sigma$ uncertainties. Nevertheless, statistical analyses show that the solutions from Rowe et al. 2015 data better characterize the data. The first and second solutions identify masses of $5_{-2}^{+4}~M_{\mathrm{\oplus}}$ and $4.6_{-2.0}^{+3.6}~M_{\mathrm{\oplus}}$, and $3.0^{+0.8}_{-0.8}~M_{\mathrm{\oplus}}$ and $2.6^{+1.9}_{-0.8}~M_{\mathrm{\oplus}}$ for the inner and outer planet, respectively. This points to a system with an inner super-Earth and an outer mini-Neptune. Dynamical studies show the planets have almost co-planar orbits with small eccentricities ($e<0.1$), close but not into the 3:2 mean motion resonance. Stability analysis indicates that this configuration is stable over million years of evolution.

We perform dynamical fits to the Transit Timing Variations (TTVs) of Kepler-419b. The TTVs from 17 Kepler quarters are obtained from Holczer et al (2016, ApJS 225,9). The dynamical fits are performed using the MultiNest Bayesian inference tool, coupled to an efficient symplectic N-body integrator. We find that the existing TTV data alone are able to uniquely constrain the planetary masses of Kepler-419b and c. Our estimates are in a good agreement with previous mass determinations that combined different techniques and observations, such as TTVs, radial velocity measurements, and the photoeccentric effect. As expected, however, our mass estimates have larger uncertainty. We study the global stability of the system within the parameters uncertainties to discard possible unstable solutions. We conclude that our method applied to the Holczer et al. data can provide reliable determinations of the planetary parameters. In the fore-coming work, we will use Holczer et al. (2016b) data to determine planetary parameters for a large set of Kepler systems.

The stability of satellites in the solar system is affected by the so-called evection resonance. The moons of Saturn, in particular, exhibit a complex dynamical architecture in which co-orbital configurations occur, especially close to the planet where this resonance is present. We address the dynamics of the evection resonance, with particular focus on the Saturn system, and compare the known behavior of the resonance for a single moon to that of a pair of moons in co-orbital trojan configuration. We developed an analytic expansion of the averaged Hamiltonian of a trojan pair of bodies, including the perturbation from a distant massive body. {The analysis of the corresponding equilibrium points was restricted to the asymmetric apsidal corotation solution of the co-orbital dynamics.} We also performed numerical N-body simulations to construct dynamical maps of the stability of the evection resonance in the Saturn system, and to study the effects of this resonance under the migration of trojan moons caused by tidal dissipation. The structure of the phase space of the evection resonance for trojan satellites is similar to that of a single satellite, differing in that the libration centers are displaced from their standard positions by an angle that depends on the periastron difference $\varpi_2-\varpi_1$ and on the mass ratio $m_2/m_1$ of the trojan pair. The interaction with the inner evection resonance may have been relevant during the early evolution of the Saturn moons Tethys, Dione, and Rhea. In particular, Rhea may have had trojan companions in the past that were lost when it crossed the evection resonance, while Tethys and Dione may either have retained their trojans or have never crossed the evection. This may help to constrain the dynamical processes that led to the migration of these satellites and to the evection itself.

We use 16 quarters of the \textit{Kepler} mission data to analyze the transit timing variations (TTVs) of the extrasolar planet Kepler-46b (KOI-872). Our dynamical fits confirm that the TTVs of this planet (period $P=33.648^{+0.004}_{-0.005}$ days) are produced by a non-transiting planet Kepler-46c ($P=57.325^{+0.116}_{-0.098}$ days). The Bayesian inference tool \texttt{MultiNest} is used to infer the dynamical parameters of Kepler-46b and Kepler-46c. We find that the two planets have nearly coplanar and circular orbits, with eccentricities $\simeq 0.03$ somewhat higher than previously estimated. The masses of the two planets are found to be $M_{b}=0.885^{+0.374}_{-0.343}$ and $M_{c}=0.362^{+0.016}_{-0.016}$ Jupiter masses, with $M_{b}$ being determined here from TTVs for the first time. Due to the precession of its orbital plane, Kepler-46c should start transiting its host star in a few decades from now.