We present 2D hydrodynamical simulations of hot Jupiters orbiting near the inner edge of protoplanetary discs. We systemically explore how the accretion rate at the inner disc edge is regulated by a giant planet of different mass, orbital separation and eccentricity. We find that a massive (with planet-to-star mass ratio $\gtrsim 0.003$) eccentric ($e_p\gtrsim 0.1$) planet drives a pulsed accretion at the inner edge of the disc, modulated at one or two times the planet's orbital frequency. The amplitude of accretion variability generally increases with the planet mass and eccentricity, although some non-monotonic dependences are also possible. Applying our simulation results to the T Tauri system CI Tau, where a young hot Jupiter candidate has been detected, we show that the observed luminosity variability in this system can be explained by pulsed accretion driven by an eccentric giant planet.

We develop a simplified model for studying the long-term evolution of giant planets in protoplanetary discs. The model accounts for the eccentricity evolution of the planets and the dynamics of eccentric discs under the influences of secular planet-disc interactions and internal disc pressure, self-gravity and viscosity. Adopting the ansatz that the disc precesses coherently with aligned apsides, the eccentricity evolution equations of the planet-disc system reduce to a set of linearized ODEs, which allows for fast computation of the evolution of planet-disc eccentricities over long timescales. Applying our model to "giant planet + external disc" systems, we are able to reproduce and explain the secular behaviours found in previously published hydrodynamical simulations. We re-examine the possibility of eccentricity excitation (due to secular resonance) of multiple planets embedded in a dispersing disc, and find that taking into account the dynamics of eccentric discs can significantly affect the evolution of the planets' eccentricities.

The population of giant planets on short-period orbits can potentially be explained by some flavours of high-eccentricity migration. In this paper we investigate one such mechanism involving "secular chaos", in which secular interactions between at least three giant planets push the inner planet to a highly eccentric orbit, followed by tidal circularization and orbital decay. In addition to the equilibrium tidal friction, we incorporate dissipation due to dynamical tides that are excited inside the giant planet. Using the method of Gaussian rings to account for planet-planet interactions, we explore the conditions for extreme eccentricity excitation via secular chaos and the properties of hot Jupiters formed in this migration channel. Our calculations show that once the inner planet reaches a sufficiently large eccentricity, dynamical tides quickly dissipate the orbital energy, producing an eccentric warm Jupiter, which then decays in semi-major axis through equilibrium tides to become a hot Jupiter. Dynamical tides help the planet avoid tidal disruption, increasing the chance of forming a hot Jupiter, although not all planets survive the process. We find that the final orbital periods generally lie in the range of 2-3 days, somewhat shorter than those of the observed hot Jupiter population. We couple the planet migration to the stellar spin evolution to predict the final spin-orbit misalignments. The distribution of the misalignment angles we obtain shows a lack of retrograde orbits compared to observations. Our results suggest that high-eccentricity migration via secular chaos can only account for a fraction of the observed hot Jupiter population.

During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($t\sim 3\times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{\rm d}/M_{\rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{\rm d}/M_{\rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.

We explore the evolution of the eccentricity of an accretion disc perturbed by an embedded planet whose mass is sufficient to open a large gap in the disc. Various methods for representing the orbit-averaged motion of an eccentric disc are discussed. We characterize the linear instability which leads to the growth of eccentricity by means of hydrodynamical simulations. We numerically recover the known result that eccentricity growth in the disc is possible when the planet-to-star mass ratio exceeds 0.003. For mass ratios larger than this threshold, the precession rates and growth rates derived from simulations, as well as the shape of the eccentric mode, compare well with the predictions of a linear theory of eccentric discs. We study mechanisms by which the eccentricity growth eventually saturates into a non-linear regime.

Following the recent discovery of the first radial velocity planet in a star still possessing a protoplanetary disc (CI Tau), we examine the origin of the planet's eccentricity (e $\sim 0.3$). We show through long timescale ($10^5$ orbits) simulations that the planetary eccentricity can be pumped by the disc, even when its local surface density is well below the threshold previously derived from short timescale integrations. We show that the disc may be able to excite the planet's orbital eccentricity in $<$ a Myr for the system parameters of CI Tau. We also perform two planet scattering experiments and show that alternatively the observed planet may plausibly have acquired its eccentricity through dynamical scattering of a migrating lower mass planet, which has either been ejected from the system or swallowed by the central star. In the latter case the present location and eccentricity of the observed planet can be recovered if it was previously stalled within the disc's magnetospheric cavity.

We formulate a set of linear equations that describe the behaviour of small eccentricities in a protoplanetary system consisting of a gaseous disc and a planet. Eccentricity propagates through the disc by means of pressure and self-gravity, and is exchanged with the planet via secular interactions. Excitation and damping of eccentricity can occur through Lindblad and corotation resonances, as well as viscosity. We compute normal modes of the coupled disc-planet system in the case of short-period giant planets orbiting inside an inner cavity, possibly carved by the stellar magnetosphere. Three-dimensional effects allow for a mode to be trapped in the inner parts of the disc. This mode can easily grow within the disc's lifetime. An eccentric mode dominated by the planet can also grow, although less rapidly. We compute the structure and growth rates of these modes and their dependence on the assumed properties of the disc.

Winds from short-period Earth and Neptune mass exoplanets, driven by high energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super Earth mass planet that loses a large fraction of its mass in a wind could drift a few percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

We determine, analytically and numerically, the conditions needed for a system of two migrating planets trapped in a 2:1 mean motion resonance to enter an inclination-type resonance. We provide an expression for the asymptotic equilibrium value that the eccentricity $e_{\rm i}$ of the inner planet reaches under the combined effects of migration and eccentricity damping. We also show that, for a ratio $q$ of inner to outer masses below unity, $e_{\rm i}$ has to pass through a value $e_{\rm i,res}$ of order 0.3 for the system to enter an inclination-type resonance. Numerically, we confirm that such a resonance may also be excited at another, larger, value $e_{\rm i, res} \simeq 0.6$, as found by previous authors. A necessary condition for onset of an inclination-type resonance is that the asymptotic equilibrium value of $e_{\rm i}$ is larger than $e_{\rm i,res}$. We find that, for $q \le 1$, the system cannot enter an inclination-type resonance if the ratio of eccentricity to semimajor axis damping timescales $t_e/t_a$ is smaller than 0.2. This result still holds if only the eccentricity of the outer planet is damped and $q \lesssim 1$. As the disc/planet interaction is characterized by $t_e/t_a \sim 10^{-2}$, we conclude that excitation of inclination through the type of resonance described here is very unlikely to happen in a system of two planets migrating in a disc.

Observations of exoplanets over the last two decades have revealed a new class of Jupiter-size planets with orbital periods of a few days, the so-called "hot Jupiters". Recent measurements using the Rossiter-McLaughlin effect have shown that many (~ 50%) of these planets are misaligned; furthermore, some (~ 15%) are even retrograde with respect to the stellar spin axis. Motivated by these observations, we explore the possibility of forming retrograde orbits in hierarchical triple configurations consisting of a star-planet inner pair with another giant planet, or brown dwarf, in a much wider orbit. Recently Naoz et al. (2011) showed that in such a system, the inner planet's orbit can flip back and forth from prograde to retrograde, and can also reach extremely high eccentricities. Here we map a significant part of the parameter space of dynamical outcomes for these systems. We derive strong constraints on the orbital configurations for the outer perturber that could lead to the formation of hot Jupiters with misaligned or retrograde orbits. We focus only on the secular evolution, neglecting other dynamical effects such as mean-motion resonances, as well as all dissipative forces. For example, with an inner Jupiter-like planet initially on a nearly circular orbit at 5 AU, we show that a misaligned hot Jupiter is likely to be formed in the presence of a more massive planetary companion (> 2 MJ) within 140 AU of the inner system, with mutual inclination 50 degrees and eccentricity above 0.25. This is in striking contrast to the test-particle approximation, where an almost perpendicular configuration can still cause large eccentricity excitations, but flips of an inner Jupiter-like planet are much less likely to occur. The constraints we derive can be used to guide future observations, and, in particular, searches for more distant companions in systems containing a hot Jupiter.

The secular approximation for the evolution of hierarchical triple configurations has proven to be very useful in many astrophysical contexts, from planetary to triple-star systems. In this approximation the orbits may change shape and orientation, on time scales longer than the orbital time scales, but the semi major axes are constant. For example, for highly inclined triple systems, the Kozai-Lidov mechanism can produce large-amplitude oscillations of the eccentricities and inclinations. Here we revisit the secular dynamics of hierarchical triple systems. We derive the secular evolution equations to octupole order in Hamiltonian perturbation theory. Our derivation corrects an error in some previous treatments of the problem that implicitly assumed a conservation of the z-component of the angular momentum of the inner orbit (i.e., parallel to the total angular momentum of the system). Already to quadrupole order, our results show new behaviors including the possibility for a system to oscillate from prograde to retrograde orbits. At the octupole order, for an eccentric outer orbit, the inner orbit can reach extremely high eccentricities and undergo chaotic flips in its orientation. We discuss applications to a variety of astrophysical systems, from stellar triples to merging compact binaries and planetary systems. Our results agree with those of previous studies done to quadrupole order only in the limit in which one of the inner two bodies is a massless test particle and the outer orbit is circular;our results agree with previous studies at octupole order for the eccentricity evolution, but not for the inclination evolution.

We study the dynamics of a planet on an orbit inclined with respect to a disc. If the initial inclination of the orbit is larger than some critical value, the gravitational force exerted by the disc on the planet leads to a Kozai cycle in which the eccentricity of the orbit is pumped up to large values and oscillates with time in antiphase with the inclination. On the other hand, both the inclination and the eccentricity are damped by the frictional force that the planet is subject to when it crosses the disc. We show that, by maintaining either the inclination or the eccentricity at large values, the Kozai effect provides a way of delaying alignment with the disc and circularization of the orbit. We find the critical value to be characteristically as small as about 20 degrees. Typically, Neptune or lower mass planets would remain on inclined and eccentric orbits over the disc lifetime, whereas orbits of Jupiter or higher mass planets would align and circularize. This could play a significant role in planet formation scenarios.

About 25 per cent of `hot Jupiters' (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbit's angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet-star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.