Possible solution to the riddle of HD 82943 multiplanet system: the three-planet resonance 1:2:5?

We carry out a new analysis of the published radial velocity data for the planet-hosting star HD82943. We include the recent Keck/HIRES measurements as well as the aged but much more numerous CORALIE data. We find that the CORALIE radial velocity measurements are polluted by a systematic annual variation which affected the robustness of many previous results. We show that after purging this variation, the residuals still contain a clear signature of an additional $\sim 1100$ days periodicity. The latter variation leaves significant hints in all three independent radial velocity subsets that we analysed: the CORALIE data, the Keck data acquired prior to a hardware upgrade, and the Keck data taken after the upgrade. We mainly treat this variation as a signature of a third planet in the system, although we cannot rule out other interpretations, such as long-term stellar activity. We find it easy to naturally obtain a stable three-planet radial-velocity fit close to the three-planet mean-motion resonance 1:2:5, with the two main planets (those in the 1:2 resonance) in an aligned apsidal corotation. The dynamical status of the third planet is still uncertain: it may reside in as well as slightly out of the 5:2 resonance. We obtain the value of $\sim 1075$ days for its orbital period and of $\sim 0.3 M_{\rm Jup}$ for its minimum mass, while the eccentric parameters are uncertain.Comment: 18 pages, 5 tables, 18 figures; accepted for publication in MNRA

Secular Dynamics of S-type Planetary Orbits in Binary Star Systems: Applicability Domains of First- and Second-Order Theories

We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary. We show that each model has both advantages and limitations. While the first-order analytical model is algebraically simple and easy to implement, it is only applicable in regions of the parameter space where the perturbations are sufficiently small. The second-order model, although more complex, has a larger range of validity and must be taken into account for dynamical studies of some real exoplanetary systems such as $\gamma$-Cephei and HD 41004A. However, in some extreme cases, neither of these analytical models yields quantitatively correct results, requiring either higher-order theories or direct numerical simulations. Finally, we determine the limits of applicability of each analytical model in the parameter space of the system, giving an important visual aid to decode which secular theory should be adopted for any given planetary system in a close binary.Comment: 32 pages, 8 figures, accepted for publication in Celestial Mechanics and Dynamical Astrophysic

Long-term and large-scale hydrodynamical simulations of migrating planets

We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.Comment: Accepted for publication in ApJ. For a movie describing the method, see https://youtu.be/66o0Z2lX8N

Chaotic diffusion in the Gliese-876 planetary system

Chaotic diffusion is supposed to be responsible for orbital instabilities in planetary systems after the dissipation of the protoplanetary disc, and a natural consequence of irregular motion. In this paper, we show that resonant multiplanetary systems, despite being highly chaotic, not necessarily exhibit significant diffusion in phase space, and may still survive virtually unchanged over time-scales comparable to their age. Using the GJ-876 system as an example, we analyse the chaotic diffusion of the outermost (and less massive) planet. We construct a set of stability maps in the surrounding regions of the Laplace resonance. We numerically integrate ensembles of close initial conditions, compute Poincaŕe maps and estimate the chaotic diffusion present in this system. Our results show that, the Laplace resonance contains two different regions: an inner domain characterized by low chaoticity and slow diffusion, and an outer one displaying larger values of dynamical indicators. In the outer resonant domain, the stochastic borders of the Laplace resonance seem to prevent the complete destruction of the system. We characterize the diffusion for small ensembles along the parameters of the outermost planet. Finally, we perform a stability analysis of the inherent chaotic, albeit stable Laplace resonance, by linking the behaviour of the resonant variables of the configurations to the different sub-structures inside the three-body resonance.Instituto de Astrofísica de La PlataFacultad de Ciencias Astronómicas y Geofísica

Chaotic diffusion in the Gliese-876 planetary system

Chaotic diffusion is supposed to be responsible for orbital instabilities in planetary systems after the dissipation of the protoplanetary disc, and a natural consequence of irregular motion. In this paper, we show that resonant multiplanetary systems, despite being highly chaotic, not necessarily exhibit significant diffusion in phase space, and may still survive virtually unchanged over time-scales comparable to their age. Using the GJ-876 system as an example, we analyse the chaotic diffusion of the outermost (and less massive) planet. We construct a set of stability maps in the surrounding regions of the Laplace resonance. We numerically integrate ensembles of close initial conditions, compute Poincaŕe maps and estimate the chaotic diffusion present in this system. Our results show that, the Laplace resonance contains two different regions: an inner domain characterized by low chaoticity and slow diffusion, and an outer one displaying larger values of dynamical indicators. In the outer resonant domain, the stochastic borders of the Laplace resonance seem to prevent the complete destruction of the system. We characterize the diffusion for small ensembles along the parameters of the outermost planet. Finally, we perform a stability analysis of the inherent chaotic, albeit stable Laplace resonance, by linking the behaviour of the resonant variables of the configurations to the different sub-structures inside the three-body resonance.Instituto de Astrofísica de La PlataFacultad de Ciencias Astronómicas y Geofísica

Tidal evolution of a close-in planet with a more massive outer companion

We investigate the motion of a two-planet coplanar system under the combined effects of mutual interaction and tidal dissipation. The secular behavior of the system is analyzed using two different approaches, restricting to the case of a more massive outer planet. First, we solve the exact equations of motion through the numerical simulation of the system evolution. We also compute the stationary solutions of the mean equations of motion based on a Hamiltonian formalism. An application to the real system CoRoT-7 is investigated.Facultad de Ciencias Astronómicas y Geofísica

Tidal evolution of a close-in planet with a more massive outer companion

We investigate the motion of a two-planet coplanar system under the combined effects of mutual interaction and tidal dissipation. The secular behavior of the system is analyzed using two different approaches, restricting to the case of a more massive outer planet. First, we solve the exact equations of motion through the numerical simulation of the system evolution. We also compute the stationary solutions of the mean equations of motion based on a Hamiltonian formalism. An application to the real system CoRoT-7 is investigated.Facultad de Ciencias Astronómicas y Geofísica

On the chaotic diffusion in multidimensional Hamiltonian systems

We present numerical evidence that diffusion in the herein studied multidimensional near-integrable Hamiltonian systems departs from a normal process, at least for realistic timescales. Therefore, the derivation of a diffusion coefficient from a linear fit on the variance evolution of the unperturbed integrals fails. We review some topics on diffusion in the Arnold Hamiltonian and yield numerical and theoretical arguments to show that in the examples we considered, a standard coefficient would not provide a good estimation of the speed of diffusion. However, numerical experiments concerning diffusion would provide reliable information about the stability of the motion within chaotic regions of the phase space. In this direction, we present an extension of previous results concerning the dynamical structure of the Laplace resonance in Gliese-876 planetary system considering variations of the orbital parameters accordingly to the error introduced by the radial velocity determination. We found that a slight variation of the eccentricity of planet c would destabilize the inner region of the resonance that, though chaotic, shows stable when adopting the best fit values for the parameters.Facultad de Ciencias Astronómicas y Geofísica

Tidal evolution of a close-in planet with a more massive outer companion

We investigate the motion of a two-planet coplanar system under the combined effects of mutual interaction and tidal dissipation. The secular behavior of the system is analyzed using two different approaches, restricting to the case of a more massive outer planet. First, we solve the exact equations of motion through the numerical simulation of the system evolution. We also compute the stationary solutions of the mean equations of motion based on a Hamiltonian formalism. An application to the real system CoRoT-7 is investigated.Facultad de Ciencias Astronómicas y Geofísica