---
author:
- 'J. Lillo-Box, D. Barrado, A. C. M. Correia'
bibliography:
- 'biblio2.bib'
date: 'Accepted on March 11th, 2016'
subtitle: 'Lack of hot-Jupiters and prevalence of multi-planetary systems'
title: 'Close-in planets around giant stars'
---

[Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that $70.0\pm6.6$% of the planets in this regime are inner components of multiplanetary systems. This value is 4.2$\sigma$ higher than for main-sequence hosts, which we find to be $42.4\pm0.1$%. [The properties of the known planets seem to indicate that the closest-in planets ($a<0.06$ au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits $0.06<a<0.5$ au. We discuss the implications of this emerging observational trend in the context of formation and evolution of hot Jupiters.]{}]{}

Introduction
============

The large crop of extrasolar planets discovered so far offers a new possibility of studying the properties of these systems from a statistical point of view [e.g., @batalha14]. This bounty of planetary systems and their widely diverse properties allow us to start answering the question of how these bodies form and evolve from an observational point of view. However, while hundreds of planets have been found around solar-mass stars, there is still a scarcity of planets orbiting more massive hosts. Some of these hosts have already left the main sequence, now being in the giant or subgiant phase. In consequence, looking for planets around these more evolved stars (with much sharper absorption lines) can help to better constrain the demography of planets around early-type stars and so probing the efficiency of planet formation mechanisms for the different mass ranges of the host star. The discovery of planets around K and G giants and subgiants is therefore crucial for planet formation theories.

From the observational point of view, there is a dearth of planets with short periods around stars ascending the red giant branch (RGB, @johnson07). The reason could be twofold. On one hand, it could be explained by a scarcity of close-in planets around early-type main-sequence stars ($M_{\star} > 1.2 M_{\odot}$). This scarcity has been hypothesized to be related to different migration mechanisms for planets around stars of different masses (Udry et al., 2003), owing to the shorter dissipation timescales of the protoplanetary disks of these stars, which prevents the formed planets from migrating to close-in orbits [see, e.g., @burkert07; @currie09]. On the other hand, the paucity of planets around these evolved stars has been considered by theoretical studies as evidence of the planet engulfment or disruption even in the first stages of the evolution off the main sequence of their parents. [@villaver09] calculated how tidal interactions in the subgiant and giant stages can lead to the final engulfment of the close-in planets and how this process is more efficient for more massive planets. Likewise, [@villaver14] have analyzed the effects of the evolution of the planet’s orbital eccentricity, mass-loss rate, and planetary mass on the survivability of planets orbiting massive stars ($M_{\star} > 1.5 M_{\odot}$). They conclude that the planet mass is a key parameter for the engulfment during the subgiant phase, with the more massive planets more likely falling into the stellar envelope during this phase (for the same initial orbital parameters). Also, planets located at 2-3 $R_{\star}$, when the star begins to leave the main sequence, may suffer orbital decay from the influence of stellar tides. Any planet closer than this orbital distance may be engulfed by the star in the subgiant phase. However, these results are based on a limited sample of confirmed exoplanets around RGB stars. Therefore the detection of close-in planets around post main-sequence stars is crucial for constraining theoretical models of planet engulfment. Several long-term projects have therefore focused on finding these planets (e.g., TAPAS: @niedzielski15; EXPRESS: @jones11).

In this paper, we analyze the growing sample of close-in planets around giant stars [from an observational point of view]{}. [We analyze their properties and multiplicity and compare them to planets around main-sequence stars, trying to determine the evolution of these properties. ]{}[We discuss the possible explanations for the observed trend in the context of the formation and evolution of hot-Jupiter planets and the consequences of their inward migration.]{}

Properties of the known close-in planets around giant stars {#sec:properties}
===========================================================

From the sample of confirmed or validated extrasolar planets, we have selected those with close-in orbits around stars in the subgiant phase or ascending the RGB. We consider planets in this regime as those revolving in orbits closer than $a<0.5$ au around host stars with $\log{g}<3.8$. The limit in semi-major axis was chosen as the apparent limit closer than which no planets have been found until recent years. The limit in surface gravity was chosen to select both giant and subgiant stars. In Fig. \[fig:logg\_sma\], we show the location of all known planets in the $[a,\log{g}]$ parameter space. In total, 13 systems have been identified so far in the mentioned regime (including uncertainties): Kepler-56, Kepler-91, Kepler-108, Kepler-278, Kepler-368, Kepler-391, Kepler-432, 8Umi, 70Vir, HD11964, HD38529 , HD102956, and HIP67851. In Table \[tab:properties\] we present the main properties of the planets and host stars in these systems. We include all systems that lie inside the above-mentioned boundaries within their $3\sigma$ uncertainties. In total, the sample is composed of four su-giant hosts and nine stars already ascending the RGB.

![Semi-major axis of known planets against their host surface gravity. We include single planets (red plus symbols) and multiplanetary systems (circles) with the inner component being the dark blue and the outer planet(s) in light blue. [For reference, we also plot the stellar radius of the host star at every evolutionary stage (i.e., for the different surface gravities) for two masses that comprise the range of masses of the hosts analyzed in this paper, $1~M_{\odot}$(dotted) and $1.6~M_{\odot}$(dash-dotted).]{} []{data-label="fig:logg_sma"}](sma_logg2.pdf){width="50.00000%"}

Results {#sec:results}
=======

We can see from Fig. \[fig:logg\_sma\] that most of the planets with $\log{g}<3.8$ and $a<0.5$ au are multiplanetary systems with at least the inner component in this regime. In four cases, no additional planets are reported by the discovery papers. These are Kepler-91, HD102956, 8Umi, and 70Vir. In this section we calculate the actual rate (and uncertainty) of multiplanetary systems in this regime and analyze the RV observations of the single systems to put constraints on possibly undetected outer planets.

We retrieved the host masses and radii from the Exoplanet Catalogue[^1] and the NASA Exoplanet Archive[^2]. The surface gravity and its uncertainty is then computed from those values. The semi-major axis was also retrieved from the same source. We divide the $[\log{g},a]$ parameter space into three regimes (see Fig. \[fig:logg\_sma\]): close-in planets around main-sequence stars (MS$_{\rm close}$, $\log{g}>3.8$ and $a<0.5$ au), close-in planets around giant stars (GS$_{\rm close}$, $\log{(g)}<3.8$ and $a<0.5$ au), and far away planets around giant stars (GS$_{\rm far}$, $\log{g}<3.8$ and $a>0.5$ au). We bootstrap the values of all planets by considering Gaussian distribution of probability with standard deviation equal to the uncertainty in the parameters. In each step, we count the number of planets in each particular regime that are inner components in multiplanetary systems with respect to the total number of systems in that regime ($\zeta$). We ran $10^3$ steps and obtained the distribution of $\zeta$ for each regime (see Fig. \[fig:distributions\]).

The median and standard deviation of the distribution of the bootstrapped populations for each regime are kept. These values represent the ratio of inner planets in multi-planetary systems in each of the three regimes. We find that $\zeta_{\rm MS,close}=42.4\pm0.1$% of the planetary systems hosted by main sequence stars with at least one component closer than 0.5 au are multiplanetary systems. Regarding systems in orbits farther than 0.5 au, we find that $\zeta_{\rm GS,far}=7.1\pm0.1$% of the systems in this regime are multiple systems with their inner component farther than 0.5 au. Finally, we find that in $\zeta_{\rm GS,close}=70.0\pm6.6$% of the systems with the closer planet more interior than 0.5 au, this planet is the inner component of a multi-planetary system. According to this, the ratio of close-in planets that are inner components in multiplanetary systems is significantly greater ($>4.2\sigma$) in the case of giant stars than it is during the main-sequence stage. The implications of this significant difference are discussed in Sect. \[sec:discussion\].

We must note that we are aware that these numbers may suffer from some observational biases (specially $\zeta_{\rm GS,far}$) since, [for instance, detecting outer components to planets with $a>0.5$ au around giant stars by]{} both radial velocity and transit methods is much more difficult or, at the least, less likely. Also, detecting planets with the transit method is more difficult in the case of giants owing to the contrast radius ratio. These biases have not been taken into account in this calculation. However, it is important to point out that both biases would favor an increase of $\zeta_{\rm GS,close}$ over $\zeta_{\rm MS,close}$.

![image](all_percentages.pdf){width="100.00000%"}

![image](99_min_Mp_hd102956.pdf){width="33.00000%"} ![image](99_min_Mp_8umi.pdf){width="33.00000%"} ![image](99_min_Mp_70vir.pdf){width="33.00000%"}

Among the four planetary systems in close-in orbits around giant stars found to host only one planet, we have investigated the mass-period parameter space yet unexplored by the discovery RV data to put constraints on additional undetected planet. We have used the RV data published in the discovery papers of HD102956 (22 epochs in a 1164-days timespan; @johnson10), 8 Ursa Minor (21 measurements in a 1888-day timespan; @lee15), and 70 Vir (@naef03) to check for the limitations of these observations for detecting additional planets in the system. We used these data and the derived parameters in the discovery works to remove the contribution of the known planet and look for a second planet in the residuals.

This was done by following the prescriptions in [@lagrange09], by simulating a planet signal in a circular orbit with the corresponding $[m_p,P_{\rm orb}]$ and observed at the same dates, and by introducing the typical noise for the particular instrument used in the discovery of the inner planet. The results of the simulation are provided in Fig. \[fig:minimumMp\], where we set the properties of the planets that could have been missed by the discovery papers. Regarding Kepler-91 [@lillo-box14; @lillo-box14c], subsequent follow-up observations after its confirmation pointed out the possibility that additional bodies exist by detecting a possible radial velocity drift [@sato15] or by detecting additional dips in the *Kepler* phase-folded light curve.

In these four cases, the possibility of an additional outer planet is not discarded among the presented limits, but we count them as single for the purposes of this work.

Discussion and conclusions {#sec:discussion}
==========================

The results presented in the previous section indicate that [(according to observations)]{} the large majority of planetary systems in close-in orbits around giant stars are multiple. We have found that the ratio of these systems against the total number of systems is significantly higher in giant and subgiant stars than in main-sequence hosts. Even though we are in the low number statistics, there seems to be an increasing sample of multiplanetary systems around evolved stars with their inner component in an orbit closer than 0.5 au. We thus may wonder about the reasons for this higher rate against single planets [and about how this is related to the evolution of planetary systems.]{}

![image](MS_Radius_single_vs_multi.pdf){width="45.00000%"} ![image](MS_Mass_single_vs_multi.pdf){width="45.00000%"}

It has been theorized that giant planets formed in the outer regions of the disk are capable of destroying or swallowing smaller planets in their inward migration to the inner parts of the system [e.g., @mustill15]. This could explain the high concentration of single planets in orbits closer than 0.06 au around main-sequence stars (see Fig. \[fig:logg\_sma\]). [Indeed, 68.1% of the planets found in this regime are found to be single planets, while inner components of multiplanetary systems represent the 31.0% of the systems. By contrast, planets at 0.06-0.5 au are mostly multiple systems with at least the inner component in this range. According to the sample of known planets used in this paper, we find that 64.5% of the systems in this regime are inner planets in multiplanetary systems, while just 35.5% are single-planet systems.\
By focusing on the properties of the extremely close-in planets ($a<0.06$ au) around main-sequence stars, we see a clear segregation, with the single planets being Jupiter-like and the inner components of multiplanetary systems being small (rocky) planets. This segregation is clear in Fig. \[fig:closein\], where we show the mass and radius distribution of extremely close-in planets around main-sequence stars for the sample of single planetary systems and inner planets in multiplanetary systems. From this reasoning, it is clear that during the main-sequence stage, planets in these extremely close-in orbits are mainly single Jupiter-like planets (i.e., hot Jupiters). This observational result is in good agreement with the idea that the migration of massive planets from the outer regions of the disk tends to destroy any already formed rocky planets in the inner regions.\
It is important to point out that this destruction of the smaller planets may not always take place as shown in the case of the WASP-47 system, where a Neptune-sized outer planet and a super-Earth inner companion to a hot-Jupiter planet were found by [@becker15]. This discovery highlights that several effects may be taking place during this process.]{}

[The previously summarized conclusions are relevant for the results of this work since we have not found extremely close-in planets around giant or subgiant stars. However, we find planets in the 0.06-0.5 au regime around these evolved stars. And, more importantly, we find that they are mostly inner components in multiplanetary systems in a fraction of $70.0\pm6.6$ %, which is compatible with the 64.5% previously mentioned for the main-sequence systems in the same regime. The observational conclusion that can be extracted from these numbers is that once the star evolves off the main-sequence, the closest planets (usually massive and single) are swallowed by their host, leaving only the population of multiplanetary systems in orbits beyond 0.06 au. Numerical simulations by [@frewen16] have shown that this could indeed explain the deficiency of warm Jupiters in stars larger than two solar radii. They hypothesize that these gaseous planets with periods of 10-100 days may migrate owing to Kozai-Lidov oscillations with a subsequent decay via tidal interactions with the star. In the current sample, the single planets Kepler-91b and HD102956b are indeed single hot Jupiters in a very close-in orbit ($a<0.1$ au), so they seem to have (by now) survived this phase.]{}

[According to [@villaver14], more massive planets are more easily and more rapidly engulfed during the evolution of their host star off the main sequence. It is thus clear from observations that hot Jupiters should practically disappear at these stages. Low-mass planets, however, could survive at least during the first stages of the evolution of the star. But, their detection is much more challenging given the small radius ratio compared to the giant planets (i.e., it makes it difficult to detect their transits) and the low-mass ratios combined with stellar pulsations during this evolutionary stage (what difficult the detection by the radial velocity technique). However, small planets should be there, and future instrumentation (e.g., ESPRESSO/VLT) will be capable of detecting them, filling the new gap of extremely close-in planets around evolved stars.]{}

[There are other possibilities for explaining the prevalence of close-in multiplanetary systems around giant stars. For instance, gravitational interactions between the planets can play an important role during the subgiant phase, although they tend to destabilize the system [@voyatzis13; @veras13; @mustill14]. However, mean motion resonances (MMR) can help the inner orbit to migrate outward and thus escape being engulfed by the star, as in the formation of the solar system [e.g., @levison03; @dangelo12]. Among multiplanetary systems with the inner component being closer than 0.5 au, we checked for possible MMR. Only Kepler-56 has a near MMR of 2:1, as pointed out in the discovery paper [@steffen13]. In two other cases, Kepler-108 and Kepler-432, the planets are close to MMR 4:1 and 8:1. The remaining systems do not present any resonant or near-resonant periodicities. Although in a small sample, it is therefore more difficult to explain the major presence of multiplanetary systems than single planets around subgiant and giant stars only in terms of gravitational interactions between the planets in the system.]{}

Finally, we must note that the results published here are purely observational. They can feed theoretical studies of the evolution of planetary systems across the different stellar stages to explain the planet-star interactions once the star leaves the main sequence. However, these studies are beyond the scope of this observational work, which aims to highlight the need for detecting more close-in planets around giant stars and to provide some hints for theoretical analysis.

This research was partially funded by Spanish grant AYA2012-38897-C02-01. This work was co-funded under the Marie Curie Actions of the European Commission (FP7-COFUND). A.C. acknowledges support from CIDMA strategic project UID/MAT/04106/2013.

[^1]: [exoplanet.eu](exoplanet.eu)

[^2]: [exoplanetarchive.ipac.caltech.edu](exoplanetarchive.ipac.caltech.edu)
