Anticorrelation between multiplicity and orbital properties in exoplanetary systems as a possible record of their dynamical histories
^{1} Space Science Data Center (SSDC) – ASI, via del Politecnico snc, 00133 Rome, Italy
email: angelo.zinzi@ssdc.asi.it
^{2} INAF–OAR, via Frascati n. 33, 00078 Monte Porzio Catone (RM), Italy
^{3} INAF–IAPS, via del Fosso del Cavaliere n. 100, 00133 Rome, Italy
email: diego.turrini@iaps.inaf.it
^{4} Departamento de Fisica, Universidad de Atacama, 485 Copayapu, Copiapó, Chile
Received: 19 July 2017
Accepted: 26 August 2017
Context. Previous works focused on exoplanets discovered with the radial velocity (RV) method reported an anticorrelation between the orbital eccentricities of the exoplanets and the multiplicity M (i.e., the number of planets) of the systems they inhabit.
Aims. We further investigate this reported anticorrelation here using a dataset comprising exoplanets discovered with both the RV and transit methods, searching for hints of its causes by exploring the connection between the number of planets and the dynamical state of the exosystems.
Methods. In order to examine the correlation between multiplicity and orbital eccentricity, for every multiplicity case considered (2 ≤ M ≤ 6), we computed the weighted average eccentricities instead of the median eccentricities used previously. The average eccentricities were calculated using the inverse of the uncertainty on the eccentricity values as weights. The analysis of the dynamic state of the exosystems was performed by computing their angular momentum deficit (AMD), which is a diagnostic parameter successfully used in the study of the solar system and recently applied to exosystems as well.
Results. Our results confirm the reported multiplicityeccentricity anticorrelation and show that the use of the uncertainties on the orbital eccentricities in the analysis allows for a better agreement between the data and the fits. Specifically, our best fit reproduces well the behaviour of the average eccentricities for all systems with M> 1, including the additional cases of TRAPPIST1 (M = 7) and of the solar system (M = 8). The AMD analysis, while not conclusive due to the limited number of exosystems that could be analysed, also suggests the existence of an anticorrelation between the multiplicity and the AMD of exosystems. This second anticorrelation, if confirmed by future studies, raises the possibility that the population of lowmultiplicity exosystems is contaminated by former highmultiplicity systems that became dynamically unstable and lost some of their planets.
Key words: methods: statistical / techniques: radial velocities / planets and satellites: dynamical evolution and stability
© ESO, 2017
1. Introduction
The currently known population of exoplanets shows characteristics that are extremely different from those of their counterparts orbiting the Sun and no exact analogue of the solar system has yet been found (see e.g. Hatzes 2016, for a recent review). Albeit it is likely that our picture of exoplanetary systems could be incomplete because of the limitations of our present observational capabilities, this nevertheless reveals the large variety of exosystems across the Galaxy and the diversity of the outcomes of the processes shaping the formation and evolution of planetary systems (see Hatzes 2016, and references therein).
Exploiting the possibility to perform statistical studies for the first time thanks to the large number of discovered exoplanets, recent works (Juric & Tremaine 2008; Limbach & Turner 2015) reported the existence of an anticorrelation between the multiplicity of exosystems (i.e. the number of planets they possess) and the orbital eccentricities of their planets. In particular, Limbach & Turner (2015) analysed a sample of 403 exoplanets showing that, once divided in bins according to the multiplicity of their host systems, the medians of their eccentricities can be fitted by a power law. These authors also verified that the same power law can be used to fit the median of the eccentricities of solar system planets.
In building their sample, Limbach & Turner (2015) selected exoplanets exclusively discovered by means of the radial velocity (RV) method, requiring eccentricity e measurements that are different from zero for all of the planets. This selection criterion ended up in a dataset comprising 276, 81, 25, 12 and 9 exoplanets for multiplicities M equal to 1, 2, 3, 4 and 5, 6, respectively, the last bin comprising two multiplicity values to increase the statistics. Limbach & Turner performed the analysis without considering the uncertainty on eccentricity measurements and calculated a power law fit for the data with M> 1; they obtained the relation e(M) = 0.584M^{1.20}, which they found to reproduce well the median eccentricities for M> 2.
In the present work we further investigated the findings of Limbach & Turner (2015) by taking into account the uncertainties on the eccentricity measurements in our analysis. To gain a deeper insight on the causes of the observed anticorrelation, we also performed a preliminary analysis of the dynamical state of the exosystems (instead of the individual exoplanets) in our sample by evaluating their angular momentum deficits (Laskar 2000).
This work is organized as follows: Sect. 2 is dedicated to the description and reduction of the dataset used; in Sect. 3 these data are discussed; and finally, in Sect. 4 conclusions are given.
2. Dataset and methods
To compose our sample we selected planets around stars with effective temperatures between 7920 K and 2600 K (i.e. as expected for type F, G, K, and M stars) in systems with at least two planets. We used the NASA Exoplanet Archive database (Akeson et al. 2013), which we queried by means of the Exoplanetary Analysis and 3D visualization Tool (ExoplAn3T). This web tool, currently under development at the Space Science Data Center (SSDC), allows users to select between the majority of the parameters provided by the NASA Exoplanet Archive application programming interfaces. However, since the query result also includes the required characteristics for all the objects (i.e. star and planets) that constitute the systems, the output of the tool can be used to perform analyses of the systems as a whole.
In our query we made no distinction on the discovery technique of the exoplanets, with the requirement of knowing the estimated uncertainties on their orbital eccentricities. The resulting dataset consisted of 258 exoplanets. Of these, 227 were detected by means of the RV method and 56 were detected with the transit method; the total number is larger than our sample since some exoplanets were observed with multiple techniques. While our sample is smaller than that used by Limbach & Turner (2015), our selection criterion prevents us from introducing further observational biases in the analysis, allowing us to weight the data according to their reliability and information content.
These 258 planets are distributed among the different multiplicities as described in Table 1 and, even though more than half of the planets belongs to systems with M = 2, all the considered multiplicity cases possess more than 10 planets, thus providing a statically relevant sample. Our dataset is plotted in Fig. 1, where we also show the planets of TRAPPIST1 (Gillon et al. 2017) and of the solar system^{1}. This choice provides us with systems with 7 and 8 planets (M = 7 and M = 8) respectively, which we used to verify how our results extend to higher multiplicities. These two planetary systems, however, have not been included into our analysis and are shown only for comparison purposes.
Number of planets in our dataset divided by the multiplicity (M) of their host system.
Fig. 1 Multiplicity M vs. eccentricity e for all the planets in our sample. The cases for M = 7 and M = 8 are the TRAPPIST1 planets and the planets in the solar system. The individual exoplanets in our sample are represented by black crosses. The blue squares represent the mean for every bin for planets observed with the RV method without considering errors. The red triangles are the means associated with planets detected with the RV method taking into account the uncertainties on the eccentricities measurements. Black circles indicate the means for the entire dataset accounting for measurement errors. 

Open with DEXTER 
2.1. Multiplicity and eccentricity
In analysing the dynamical properties of the exoplanets in our sample, instead of the median values of the eccentricity e in the different multiplicity bins considered by Limbach & Turner (2015), we focused on the relation between M and the weighted average eccentricities.
The weights we used in computing the average values of each multiplicity bin were the inverse of the relative errors of the eccentricities, so that (1)where is the weighted average for a specific value of the multiplicity M (from 2 to 6), e(M)_{i} are the individual eccentricity values of the exoplanets in the multiplicity bin, and σ_{ei} are the associated relative errors, which in the case of different values for positive and negative errors were computed as (2)
2.2. Multiplicity and dynamical state
Our preliminary investigation of the dynamical state of the exoplanetary systems took advantage of the information supplied by the angular momentum deficit AMD (Laskar 2000), which measures the difference between the angular momentum of a planetary system and the angular momentum the system would have if its planets were all on circular and planar orbits.
As demonstrated by Laskar & Petit (2017), the AMD remains stable after totally inelastic collisions, but increases when a body is ejected from the system. Therefore a system that experienced several ejections (hence decreasing in multiplicity M) should have a higher value of AMD with respect to a system that remained dynamically stable over all its history.
Following Laskar (2000), the AMD is defined as (3)where G is the gravitational constant, m_{k} the planetary masses (in Jovian masses), m_{0} is the mass of the host star (in solar masses), and a_{k}, e_{k}, i_{k} are the planetary semimajor axes (expressed in AU), eccentricities and inclinations, respectively. For every system, we considered as the reference inclination plane that of the orbit of the most massive planet and computed the relative inclinations of the other planets with respect to that orbit.
When computing the AMD for the different exosystems in our sample we considered two datasets (see Table 2). The first dataset (hereafter, the “large selection”) comprises 10 exosystems (31 exoplanets) with known values of the mass, semimajor axis, eccentricity, and inclination for all planets.
The second one (hereafter the “strict selection”) comprises 9 exosystems (28 exoplanets) selected with the additional requirement of also having known uncertainties on the four planetary parameters required for the computation of the AMD. To increase the statistics, in this case, we also considered those exosystems with M> 2 for which inclination values were not available for at most one planet.
Number of systems at different multiplicities considered in the analysis of their dynamical state through the AMD.
3. Results and discussion
3.1. Multiplicity and eccentricity
Figure 2 shows the behaviour, computed for three different subsets of our sample of exoplanets, and the associated power law fits. Specifically, we computed the values of as (1) simple averages considering only the 227 exoplanets discovered with the RV method, i.e. similar to the original analysis by Limbach & Turner (2015); (2) weighted averages as per Eq. (1) again considering only the 227 exoplanets discovered with the RV method; and (3) weighted averages as per Eq. (1) considering our entire dataset of 258 exoplanets.
In Table 3 we show the parameters of the fits and the associated Pearson correlation coefficients R^{2}. As shown by Table 3, there is a significant increase in the value of R^{2} when measurement errors are used as weights in computing the averages. Another noticeable effect of the use of weighted averages is the increase in the value of the mean eccentricity for the case M = 2. This results in the estimated power law to reproduce well all cases for M> 1 instead of just M> 2 as reported by Limbach & Turner (2015).
Finally, looking at Fig. 2 one can see that the behaviour of the planets belonging to TRAPPIST1 (M = 7) and the solar system (M = 8) is well reproduced by our fits, where that computed with the entire dataset most closely reproduces the average values associated with these two systems. Using the rootmeansquared deviation (RMSD) as an indicator of the goodness of the fits of the weighted averages for M ≥ 7, the case considering only RV data has RMSD = 0.020, while that considering the entire dataset has rms = 0.013.
Fig. 2 Power law fits. Blue squares indicate RV planets without considering errors, red triangles indicate RV planets with weights as errors, and black circles indicate all planets in the selected dataset with errors used as weights. 

Open with DEXTER 
Fig. 3 Angular momentum deficit of the systems considered. Crosses indicate values for the individual systems, blue squares indicate average values for the large selection sample, red circles indicate (with error bars) averages values for strict selection sample. The TRAPPIST1 value is shown as a black triangle. Error bars are smaller than the symbol sizes. 

Open with DEXTER 
3.2. Multiplicity and dynamical state
Figure 3 shows the average values of the AMD of the exosystems in our large and strict selections for every multiplicity case.
Here we omitted the comparison with the solar system as it possesses masses and distances from the host star that are extremely different from those of the more compact exosystems in our sample. To this regard, Laskar (1997) also pointed out that the AMD of the solar system shows a bimodal behaviour with the inner system (i.e. the terrestrial planets) showing characteristics of an unstable system and the outer system (i.e the giant planets) more resembling a stable system.
The analysis of the AMD shows an inverse trend with multiplicity for both the large and the strict selections. We did not attempt to fit the data as the number of systems analysed is statistically limited. Nevertheless, the fact that lower multiplicity systems display a higher value of the AMD, if confirmed when a larger sample of exosystems becomes available, could indicate that highmultiplicity systems are a common outcome of the planetary formation process, but that a significant portion of these systems form or evolves into unstable orbital configurations.
The loss of planets due to phases of dynamical instability by such systems (e.g. by planetplanet scattering events; see Weidenschilling & Marzari 1996; Rasio & Ford 1996) would cause the increase of their AMD values and, consequently, of the eccentricities of their surviving planets. The resulting systems would then contaminate the population of exosystems with lower native multiplicity values, probably explaining the behaviours observed in Figs. 2 and 3.
4. Conclusions
The results of this study confirm the reported anticorrelation between the orbital eccentricity of exoplanets and the multiplicity of the host system they inhabit. The inclusion of the uncertainties on the measured planetary eccentricities in our analysis not only provides a more robust confirmation, but also allows for fitting with a single power law the average eccentricities for multiplicities comprised between M = 2 and M = 6. Furthermore, both the characteristics of TRAPPIST1 (M = 7) and of the solar system (M = 8) are represented well by the same power law (4)In the attempt to provide a physical explanation of the reported behaviour, we made a preliminary analysis of the correlation between the angular momentum deficit AMD of the systems
included in our sample and their multiplicity values. While the data are not sufficient to provide a conclusive answer, we observed an anticorrelation between the multiplicity of the exosystems and their AMD as well. Should this AMDmultiplicity anticorrelation of exosystems be verified once larger and more precise datasets become available, it would provide us with a deeper insight both on the planetary formation process and on the place of the solar system in an exoplanetary context.
Therefore this could be an interesting field of research when the next generation of space telescopes dedicated to exoplanetary systems is ready, starting with the NASA Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015), which is planned to be operative in the next year.
Specifically, our results raise the possibility that highmultiplicity systems such as the solar system are a relatively common outcome of the planetary formation process. Our results also suggest that, unlike the solar system, a possibly large portion of these systems form or evolve into dynamically unstable orbital configurations. The ejection of one or more planets from these unstable systems would naturally explain our results by decreasing the multiplicities of the systems increasing their AMD and eccentricities and/or inclinations of their surviving planets at the same time. As a consequence, the population of lowmultiplicity exosystems we observe is likely to be contaminated by former highmultiplicity systems that became dynamically unstable and lost some of their planets during their histories.
References
 Akeson, R. L., Chen, X., Ciardi, D., et al. 2013, PASP, 125, 989 [NASA ADS] [CrossRef] [Google Scholar]
 Gillon, M., Triaud, A. H. M. J., Demory, B.O., et al. 2017, Nature, 542, 456 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Hatzes, A. P. 2016, Space Sci. Rev., 205, 267 [NASA ADS] [CrossRef] [Google Scholar]
 Juric, M., & Tremaine, S. 2008, in Extreme Solar Systems, ASP Conf. Ser., 398, 295 [NASA ADS] [Google Scholar]
 Laskar, J. 1997, A&A, 317, L75 [NASA ADS] [Google Scholar]
 Laskar, J. 2000, Phys. Rev. Lett., 84, 3240 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Laskar, J., & Petit, A. 2017, A&A, 605, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Limbach, M. A., & Turner, E. L. 2015, Proc. Nat. Acad. Sci., 112, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Rasio, F. A., & Ford, E. B. 1996, Science, 274, 954 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes, Instruments, and Systems, 1, 014003 [NASA ADS] [CrossRef] [Google Scholar]
 Weidenschilling, S. J., & Marzari, F., 1996, Nature 384, 619 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
All Tables
Number of planets in our dataset divided by the multiplicity (M) of their host system.
Number of systems at different multiplicities considered in the analysis of their dynamical state through the AMD.
All Figures
Fig. 1 Multiplicity M vs. eccentricity e for all the planets in our sample. The cases for M = 7 and M = 8 are the TRAPPIST1 planets and the planets in the solar system. The individual exoplanets in our sample are represented by black crosses. The blue squares represent the mean for every bin for planets observed with the RV method without considering errors. The red triangles are the means associated with planets detected with the RV method taking into account the uncertainties on the eccentricities measurements. Black circles indicate the means for the entire dataset accounting for measurement errors. 

Open with DEXTER  
In the text 
Fig. 2 Power law fits. Blue squares indicate RV planets without considering errors, red triangles indicate RV planets with weights as errors, and black circles indicate all planets in the selected dataset with errors used as weights. 

Open with DEXTER  
In the text 
Fig. 3 Angular momentum deficit of the systems considered. Crosses indicate values for the individual systems, blue squares indicate average values for the large selection sample, red circles indicate (with error bars) averages values for strict selection sample. The TRAPPIST1 value is shown as a black triangle. Error bars are smaller than the symbol sizes. 

Open with DEXTER  
In the text 