Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.

The aim of our study is to explore the possible existence of Earth-mass planets in the habitable zone of 55~Cancri, an effort pursued based on detailed orbital stability simulations. This star is known to possess (at least) five planets with masses ranging between super-Earth and Jupiter-type. Additionally, according to observational constraints, there is a space without planets between $\sim$0.8~au and $\sim$5.7~au, noting that the inner part of this gap largely coincides with 55~Cnc's habitable zone --- a sincere motivation for the search of potentially habitable planets. It has previously been argued that terrestrial habitable planets are able to exist in the 55~Cnc system, including a planet at $\sim$1.5~au. We explore this possibility through employing sets of orbital integrations and assuming an integration time of 50~Myr. We found that the possibility of Earth-mass planets in the system's habitable zone strongly depends on the adopted system parameters, notably the eccentricity of 55~Cnc-f, which is controversial as both a high value ($e \sim 0.32$) and a low value ($e \sim 0.08$) have previously been deduced. In case that the low value is adopted (together with other updates for the system parameters), the more plausible and most recent value, Earth-mass planets would be able to exist in the gap between 1.0~au and 2.0~au, thus implying the possibility of habitable system planets. Thus, 55~Cnc should be considered a favorable target for future habitable planet search missions.

The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are available to the community through $\texttt{Zenodo}$ and $\texttt{GitHub}$, respectively. We employ a grid interpolation method based on $\sim$150 million full N-body simulations of initially circular, coplanar systems and compare to the 9 known Kepler CBP systems. Using a formalism from planet packing studies, we find that 55% of the Kepler CBP systems allow for an additional equal-mass planet to potentially exist on an interior orbit relative to the observed planet. Therefore, we do $\textit{not}$ find strong evidence for a pile-up in the Kepler CBP systems and more detections are needed to adequately characterize the formation mechanisms for the CBP population. Observations from the Transiting Exoplanet Survey Satellite are expected to substantially increase the number of detections using the unique geometry of CBP systems, where multiple transits can occur during a single conjunction.

Stability of planetary orbits around the GJ 832 star system, which contains inner (GJ 832c) and outer (GJ 832b) planets, is investigated numerically and a detailed phase-space analysis is performed. A special emphasis is given to the existence of stable orbits for a planet less than 15M$_\oplus$ which is injected between the inner and outer planets. Thus, numerical simulations are performed for three and four bodies in elliptical orbits (or circular for special cases) by using a large number of initial conditions that cover the selected phase-spaces of the planet's orbital parameters. The results presented in the phase-space maps for GJ 832c indicate the least deviation of eccentricity from its nominal value, which is then used to determine its inclination regime relative to the star-outer planet plane. Also, the injected planet displays stable orbital configurations for at least one billion years. Then, the radial velocity curves based on the signature from the Keplerian motion are generated for the injected planets with masses 1M$_\oplus$ to 15M$_\oplus$ in order to estimate their semimajor axes and mass-limit. The synthetic RV signal suggests that an additional planet of mass $\le$ 15M$_\oplus$ with dynamically stable configuration may be residing between 0.25 - 2.0 AU from the star. We have provided an estimated number of RV observations for the additional planet that is required for further observational verification.

The idea of single exomoon detection due to the radio emissions caused by its interaction with the host exoplanet is extended to multiple-exomoon systems. The characteristic radio emissions are made possible in part by plasma from the exomoons own ionosphere (Noyola et. al. 2014). In this work, it is demonstrated that neighboring exomoons and the exoplanetary magnetosphere could also provide enough plasma to generate a detectable signal. In particular, the plasma-torus-sharing phenomenon is found to be particularly well suited to facilitate the radio detection of plasma-deficient exomoons. The efficiency of this process is evaluated, and the predicted power and frequency of the resulting radio signals are presented.

In the Jupiter-Io system, the moon's motion produces currents along the field lines that connect it to Jupiter's polar regions. The currents generate, and modulate radio emissions along their paths via the electron-cyclotron maser instability. Based on this process, we suggest that such modulation of planetary radio emissions may reveal the presence of exomoons around giant planets in exoplanetary systems. A model explaining the modulation mechanism in the Jupiter-Io system is extrapolated, and used to define criteria for exomoon detectability. A cautiously optimistic scenario of possible detection of such exomoons around Epsilon Eridani b, and Gliese 876 b is provided.

Depending on the planetary orbit around the host star(s), a planet could orbit either one or both stars in a binary system as S-type or P-type, respectively. We have analysed the dynamics of the S-type planetary system in HD 196885 AB with an emphasis on a planet with a higher orbital inclination relative to the binary plane. The mean exponential growth factor of nearby orbits (MEGNO) maps are used as an indicator to determine regions of periodicity and chaos for the various choices of the planet's semimajor axis, eccentricity and inclination with respect to the previously determined observational uncertainties. We have quantitatively mapped out the chaotic and quasi-periodic regions of the system's phase space which indicate a likely regime of the planet's inclination. In addition, we inspect the resonant angle to determine whether alternation between libration and circulation occurs as a consequence of Kozai oscillations, a probable mechanism that can drive the planetary orbit to a very large inclination. Also, we demonstrate the possible higher mass limit of the planet and improve upon the current dynamical model based on our analysis.

The orbits of two individual planets in two known binary star systems, \gamma Cephei and HD 196885 are numerically integrated using various numerical techniques to assess the chaotic or quasi-periodic nature of the dynamical system considered. The Hill stability (HS) function which measures the orbital perturbation of a planet around the primary star due to the secondary star is calculated for each system. The maximum Lyapunov exponent (MLE) time series are generated to measure the divergence/convergence rate of stable manifolds, which are used to differentiate between chaotic and non-chaotic orbits. Then, we calculate dynamical Mean Exponential Growth factor of Nearby Orbits (MEGNO) maps from solving the variational equations along with the equations of motion. These maps allow us to accurately differentiate between stable and unstable dynamical systems. The results obtained from the analysis of HS, MLE, and MEGNO maps are analysed for their dynamical variations and resemblance. The HS test for the planets shows stability and quasi-periodicity for at least ten million years. The MLE and the MEGNO maps have also indicated the local quasi-periodicity and global stability in a relatively short integration period. The orbital stability of the systems is tested using each indicator for various values of planet inclinations (i_{pl} \le 25^\circ) and binary eccentricities. The reliability of HS criterion is also discussed based on its stability results compared with the MLE and MEGNO maps.