Issue |
A&A
Volume 625, May 2019
|
|
---|---|---|
Article Number | A71 | |
Number of page(s) | 16 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201935356 | |
Published online | 14 May 2019 |
The CORALIE survey for southern extrasolar planets
XVIII. Three new massive planets and two low-mass brown dwarfs at greater than 5 AU separation★,★★
1
Département d’astronomie, Université de Genève,
51 ch. des Maillettes,
1290 Versoix, Switzerland
e-mail: emily.rickman@unige.ch
2
Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP,
Rua das Estrelas,
4150-762 Porto, Portugal
3
Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto,
Rua do Campo Alegre,
4169-007 Porto, Portugal
4
Institute of Physics, Laboratory of Astrophysics, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny,
1290, Versoix, Switzerland
5
European Southern Observatory, Alonso de Córdova 3107,
Vitacura, Santiago, Chile
6
Department of Astronomy and Astrophysics, University of Chicago,
5640 South Ellis Avenue,
Chicago,
IL 60637, USA
7
Department of Physics, University of Warwick,
Coventry CV4 7AL, UK
8
Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326,
13388 Marseille, France
9
Carnegie Institution for Science, Las Campanas Observatory,
Casilla 601, Colina El Pino s/n, La Serena, Chile
10
Physikalisches Institut, Universität Bern,
Gesellschaftsstr. 6,
3012 Bern, Switzerland
11
Instituto de Astrofísica de Canarias, 38205 La Laguna,
Tenerife,
Spain
12
School of Physics & Astronomy, University of Birmingham, Edgbaston,
Birmingham,
B15 2TT, UK
13
Leiden Observatory, Leiden University,
Postbus 9513,
2300 RA Leiden, Netherlands
Received:
25
February
2019
Accepted:
2
April
2019
Context. Since 1998, a planet-search around main sequence stars within 50 pc in the southern hemisphere has been underway with the CORALIE spectrograph at La Silla Observatory.
Aims. With an observing time span of more than 20 yr, the CORALIE survey is able to detect long-term trends in data with masses and separations large enough to select ideal targets for direct imaging. Detecting these giant companion candidates will allow us to start bridging the gap between radial-velocity-detected exoplanets and directly imaged planets and brown dwarfs.
Methods. Long-term precise Doppler measurements with the CORALIE spectrograph reveal radial-velocity signatures of massive planetary companions and brown dwarfs on long-period orbits.
Results. In this paper, we report the discovery of new companions orbiting HD 181234, HD 13724, HD 25015, HD 92987 and HD 50499. We also report updated orbital parameters for HD 50499b, HD 92788b and HD 98649b. In addition, we confirm the recent detection of HD 92788c. The newly reported companions span a period range of 15.6–40.4 yr and a mass domain of 2.93–26.77 MJup, the latter of which straddles the nominal boundary between planets and brown dwarfs.
Conclusions. We report the detection of five new companions and updated parameters of four known extrasolar planets. We identify at least some of these companions to be promising candidates for imaging and further characterisation.
Key words: techniques: radial velocities / planets and satellites: detection / binaries: visual / planetary systems
The radial velocity measurements and additional data products discussed in this paper are available on the DACE web platform at https://dace.unige.ch/radialVelocities. See the appendix for a direct link to the individual target data products. A copy of the data is also available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/625/A71
© ESO 2019
1 Introduction
Little is known about massive giant planets and brown dwarfs at orbital separations between 5 and 50 AU due to their low occurrence rate (Bowler 2016) and to the lower sensitivity of the different observing methods in this separation range. Indeed, radial velocity (RV) and transit techniques are extremely efficient at detecting planets around older stars at short separations (Fischer et al. 2014). On the other hand, direct imaging is most efficient at detecting younger planets at separations larger than several times the diffraction limit of the telescope (typically 5–10 λ∕D). This translates into several tens of astronomical units for the closest young stellar associations (e.g. β Pic and 51 Eri as part of the β Pic moving group (Zuckerman et al. 2001; Feigelson et al. 2006) and HR 8799 as part of the Columba association (Zuckerman et al. 2011). Nevertheless, the population of massive giant exoplanets at intermediate orbital separations between 5 and 50 AU is an important puzzle piece needed for constraining the uncertainties that exist in planet formation and evolution models.
The historical CORALIE planet-search survey has been ongoing for more than 20 yr in the southern hemisphere and monitors a volume-limited sample of 1647 main sequence (MS) stars from F8 down to K0 located within 50 pc of the Sun (Udry et al. 2000). With an individual measurement precision ranging between 3.5 and 6 ms −1, CORALIE has permitted (or has contributed to) the detection of more than 140 extra-solar planet candidates (Pepe et al. 2002; Udry et al. 2002; Tamuz et al. 2008; Ségransan et al. 2010; Marmier et al. 2013). Such a long and continuous monitoring of nearby MS stars is unique among all planet search surveys; it allows us to detect massive giant planets at separations larger 5 AU as well as to identify small RV drifts hinting at the presence of low-mass companions at even wider separations.
These are indeed very useful targets for direct imaging, as such old and very-low-mass companions are rare and very difficult to search for blindly. Cheetham et al. (2018) has shown with the discovery of the ultra-cool brown dwarf companion orbiting the planet host star HD 4113 A that long-term RV surveys are an extremely useful tool to select targets to image. Not only does it allow us to start filling in a largely unexplored parameter space, but through combining RV and direct imaging we can now expect to measure the masses of these companions using Kepler’s laws. By constraining the mass, we are able to place additional constraints on the evolution of the companion, both in terms of temperature and atmospheric composition (Crepp et al. 2018; Peretti et al. 2019).
In this paper we report the discovery of four new giant planets and brown dwarfs orbiting HD 181234, HD 13724, HD 25015, and HD 92987, together with the updated CORALIE orbital elements for an already known exoplanet around HD 98649 (Marmier et al. 2013). We also report updated orbital parameters for HD 50499b (Vogt et al. 2005), as well as the detection of HD 50499c, which has previously been noted by Vogt et al. (2005), Butler et al. (2017), and Barbato et al. (2018). We also report the updated orbital parameters of HD 92788b detected by Fischer et al. (2001) and confirm the recent detection of HD 92788c (Wittenmyer et al. 2019).
The paper is organised as follows. The properties of the host stars are summarised in Sect. 2. In Sect. 3 we present our RV data and the inferred orbital solution of the newly detected companions. In Sect. 4 we present the CORALIE updated parameters of already known exoplanets with new detections in two of these systems. The results are discussed in Sect. 5 with some concluding remarks.
2 Stellar characteristics
Spectral types, V band magnitude and colour indices are taken from the HIPPARCOS catalogue (Perryman et al. 1997) while astrometric parallaxes (π) and luminosities are taken from the second Gaia date release (Gaia Collaboration 2018). Effective temperatures, gravities, and metallicities are derived using the same spectroscopic methods as applied in Santos et al. (2013), whilst the v sin (i) is computed using the calibration of CORALIE’s cross correlation function (CCF; Santos et al. 2001; Marmier 2014).
The mean chromospheric activity index – – of each star is computed by co-adding the corresponding CORALIE spectra to improve the signal-to-noise ratio (S/N) which allows us to measure the Ca II re-emission at λ = 3933.66 Å. We derived an estimate of the rotational period of the star from the mean
activity index using the calibration of Mamajek & Hillenbrand (2008).
Stellar radii and their uncertainties are derived from the Gaia luminosities and the effective temperatures obtained from the spectroscopic analysis. A systematic error of 50 K was quadratically added to the effective temperature error bars and was propagated in the radius uncertainties.
The mass and the age of the stars, as well as their uncertainties, are derived using the Geneva stellar-evolution models (Ekström et al. 2012; Georgy et al. 2013). The interpolation in the model grid was made through a Bayesian formalism using observational Gaussian priors on Teff, MV, log g, and [Fe/H] (Marmier 2014).
The observed and inferred stellar parameters for newly detected host stars to planetary companions are summarised in Table 1 and host stars to brown dwarf companions in Table 2.
3 Radial velocities and orbital solutions
The CORALIE observations span over more than 20 yr, from June 1998 to December 2018. During that time, CORALIE went through two major upgrades in June 2007 (Ségransan et al. 2010) and in November 2014 to increase overall efficiency and accuracy of the instrument. These changes introduced small offsets in the measured RVs that depend on several parameters such as the spectral type of the star and its systemic velocity. For this reason, we decided to consider CORALIE as three different instruments, corresponding to the different upgrades: the original CORALIE as CORALIE-98 (C98), the first upgrade as CORALIE-07 (C07), and the latest upgrade as CORALIE-14 (C14).
In addition to the RV time series, the CORALIE automated pipeline also provides several useful indicators that help pinpoint the origin ofobserved periodic signals. These are the CCF full width at half maximum (FWHM), the bisector, and the Hα chromospheric activity indicator.
We are also using published RVs taken with other spectrographs, namely HARPS (Mayor et al. 2003), HIRES (Vogt et al. 1994), and HAMILTON (Vogt 1987). The data products presented in this paper are available at the Data and Analysis Center for Exoplanets (DACE)1.
We perform an initial modelling of the RV time series using the online DACE platform. Keplerian model initial conditions are computed using the formalism described in Delisle et al. (2016). The stellar activity detrending and the modelling of the instrumental noise and the stellar jitter follow the formalism described in Díaz et al. (2016) and Delisle et al. (2018). Analytical false-alarm probabilities (FAPs) are computed on the periodogram of the residuals following Baluev (2008) and numerical FAP values which are used in this paper are computed by permutation of the calendar. Periodic signals with a FAP lower than 0.1% are considered significant and are added to the model. For each periodogram shown, the three lines represent the 10, 1, and 0.1% FAP in ascendingorder.
Once the RV time series is fully modelled using DACE online tools, we run a Markov chain Monte Carlo (MCMC) analysis of each system using the algorithm described in Díaz et al. (2014, 2016) and Delisle et al. (2018) to obtain the posterior distributions of the model parameters. Each MCMC simulation is run with 10 000 000 iterations drawing the proposal solution obtained using DACE. The parameter confidence intervals are computed for a 68.27% confidence level.
Gaussian priors are set for the instrument offsets and stellar mass with uniform priors for the orbital elements. In the cases where the minimum RV () or the maximum RV (
) is well sampled, we perform the fit using either of these instead of fitting the phase.
In this section we present the orbital solutions for newly reported giant planets and brown dwarfs from the CORALIE survey. A summary of the orbital solutions can be found in Table 3 for the newly detected companions and the fully probed MCMC parameter space is shown in the appendix.
Observed and inferred stellar parameters for host stars to the planet candidates: HD 181234, HD 25015, HD 50499, HD 92788, and HD 98649.
Observed and inferred stellar parameters for host stars to the brown dwarf candidates: HD 13724 and HD 92987.
3.1 HD 181234 (LTT 5654, HIP 95015)
HD 181234 was observed with CORALIE at La Silla Observatory since May 2000. Fifteen measurements were taken with CORALIE-98, 21 additional RV measurements were obtained with CORALIE-07, and 59 additional RV measurements were obtained with CORALIE-14. HD 181234 has also been observed with Keck/HIRES (Butler et al. 2017) with 20 RV measurements from June 1999 to August 2014.
The best-fit Keplerian, as shown in Fig. 1, shows that we are looking at a highly eccentric system with an eccentricity of 0.73. It has an orbital period of 20.4 yr with a minimum mass of 8.4 MJup. The orbital solutions are summarised in Table 3. Figure 1 shows the CORALIE RVs and the corresponding best-fit Keplerian model along with the RV residuals and a periodogram of the residuals. The results from the fully probed parameter space from the MCMC are shown in the appendix.
3.2 HD 92987 (HIP 52472)
HD 92987 was observed with CORALIE at La Silla Observatory since January 1999. Fifty-three measurements were taken with CORALIE-98, 18 additional RV measurements were obtained with CORALIE-07, and 22 additional RV measurements were obtained with CORALIE-14.
HD 92987b is one of the brown dwarf candidates with a minimum mass of 16.88 MJup and a semi-major axis of 9.62 AU, making it a promising candidate for direct imaging. The orbital solutions for HD 92987 are summarised in Table 3. Figure 2 shows the CORALIE RVs and the corresponding best-fit Keplerian model along with the RV residuals and a periodogram of the residuals. The results from the fully probed parameter space from the MCMC are shown in the appendix.
3.3 HD 25015 (HIP 18527)
HD 25015 was observed with CORALIE at La Silla Observatory since May 2001. Twenty-two measurements were taken with CORALIE-98, 32 additional RV measurements were obtained with CORALIE-07, followed by 56 additional RV measurements obtained with CORALIE-14.
We note here that all of the stars in our sample are very quiet with the exception of HD 25015. When we plot the Hα versus time, we see that some periodicity exists in the periodogram at a FAP of greater than 10%, as seen in Fig 3 at approximately 370.26 days. We do not see a periodicity in the FWHM or bisector, therefore we cannot exclude that there is telluric line contamination in the Hα index with CORALIE. We detrend the RV from the stellar activity using the Hα indicators. A scale factor is adjusted to the smoothed Hα (using a Gaussian filter at 0.1 yr), as well as an additional jitter term proportional to the activity trend.
The orbital solutions for HD 25015 are summarised in Table 3. Figure 3 shows the CORALIE RVs and the corresponding best-fit Keplerian model along with the RV residuals, the periodogram of the residuals and a periodogram for Hα before detrending. The results from the fully probed parameter space from the MCMC are shown in the appendix.
Best-fitted solutions for the substellar companions orbiting HD 13724, HD 181234, HD 25015, and HD 92987.
![]() |
Fig. 1 Top: HD 181234 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple) and HIRES data (Butler et al. 2017) in orange. The bestsingle-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 181234. Bottom: periodogram of the residuals for HD 181234. The three black lines represent the 10, 1, and 0.1% FAPs in ascending order. |
![]() |
Fig. 2 Top: HD 92987 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The fitted single-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 92987. Bottom: periodogram of the residuals for HD 92987. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
![]() |
Fig. 3 Top: HD 25015 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The fitted single-planet Keplerian model is represented as a black curve. Second figure: RV residuals of HD 25015. Third figure: periodogram of the residuals for HD 25015 after the signal has been removed showing no significant signals. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. Bottom: periodogram of Hα before detrending. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order, showing a significant peak above the 10% FAP at 370.26 days. |
3.4 HD 13724 (HIP 10278)
HD 13724 has been observed with CORALIE at La Silla Observatory since August 1999. It is a relatively young star at 0.76± 0.71 Gyr old.
During the past 19.3 yr, 167 Doppler measurements were taken on this target with 70 RV measurements taken with CORALIE-98, 19 with CORALIE-07, 48 with CORALIE-14 and 30 with HARPS. HD 13724b is a brown dwarf companion with a minimum mass of 26.77 MJup and a semi-major axis of 12.4 AU, making it a promising candidate for direct imaging.
The orbital solutions for HD 13724 are summarised in Table 3. Figure 4 shows the CORALIE RVs and the corresponding best-fit Keplerian model along with the RV residuals and a periodogram of the residuals. The results from the fully probed parameter space from the MCMC are shown in the appendix.
![]() |
Fig. 4 Top: HD 13724 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple) and HARPS (red). The best single-planet Keplerian model is representedas a black curve. Middle: RV residuals of HD 13724. Bottom: periodogram of the residuals for HD 13724 after the signal has been removed showing no significant signals. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
4 Updated parameters for known exoplanets
We present updated orbital parameters for a known exoplanet around HD 98649 (Marmier et al. 2013). We also report updated orbital parameters for HD 92788b (Fischer et al. 2001) and confirm the detection of HD 92788c (Wittenmyer et al. 2019). Moreover, we report the discovery of a new planet around HD 50499 and report updated orbital parameters for HD 50499b (Vogt et al. 2005).
The RV data is fitted in the same way as described in Sect. 3. The stellar parameters for these systems are summarised in Table 1. For HD 50499 we also report the discovery of a new exoplanet.
The orbital parameters for HD 98649, HD 50499, and HD 92788 are summarised in Table 4. The probed physical parameters using the MCMC for each target are shown in the appendix.
4.1 A companion of 7 MJup on an eccentric orbit of 15 yr around HD 98649 (LTT 4199, HIP 55409)
We report updated parameters for a known exoplanet detected by Marmier et al. (2013). The star HD 98649 is of G3/G5V type at 42.19 ± 0.09 pc from the Sun. The star properties are summarised in Table 4.
HD 98649 has been observed with CORALIE at La Silla Observatory since February 2003. Fourteen measurements were taken with CORALIE-98, 42 additional RV measurements were obtained with CORALIE-07, followed by 12 additional RV measurements obtained with CORALIE-14.
HD 98649b is an extremely eccentric planet with an eccentricity of . The orbital parameters agree well with Marmier et al. (2013) who reported a period of
yr and a mass of Msini = 6.8 ± 0.5 MJup.
The orbital solutions for HD 98649 are summarised in Table 4. Figure 5 shows the CORALIE RVs and the corresponding best-fit Keplerian model along with the RV residuals and a periodogram of the residuals. The results from the fully probed parameter space from the MCMC are shown in the appendix.
Best-fitted solution for the substellar companions orbiting HD 50499, HD 92788, and HD 98649.
![]() |
Fig. 5 Top: HD 98649 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The best single-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 98649. Bottom: periodogram of the residuals for HD 98649 showing no significant signals remaining. The three black lines represent the 10, 1 and 0.1% FAP in ascending order. |
4.2 HD 50499 (HIP 32970)
HD 50499 is a G1V star at 46.34 ± 0.06 pc from the Sun. The star properties are summarised in Table 1.
We report updated orbital parameters for a known exoplanet around HD 50499b previously detected by Vogt et al. (2005). We also report the discovery of a new exoplanet in this system.
HD 50499 was observed with CORALIE at La Silla Observatory since January 1999. Forty-four measurements were taken with CORALIE-98, 39 additional RV measurements with CORALIE-07, followed by 40 additional RV measurements obtained with CORALIE-14. There are also 5 measurements taken with HARPS, as well as an additional 86 RV points publicly available from HIRES.
The outer signal was previously noted by Vogt et al. (2005). It was also reported by Butler et al. (2017) who noted that the outer trend is parabolic with additional data points from HIRES. In addition, Barbato et al. (2018) fits the RV data for a single Keplerian orbit plus a quadratic term obtaining a best-fit curve where lower limits for the outer companion are derived. Here we provide additional data points from CORALIE, which allows us to further constrain this outer planet.
We derive the orbital solution for HD 50499 using two Keplerians as seen in Fig. 6. The orbital solutions for HD 50499 are summarised in Table 4. Figure 6 shows the CORALIE, HARPS, and HIRES RVs and the corresponding best-fit Keplerian models along with the periodogram of the residuals. Figure 7 shows the phase-foldedRV diagrams for HD 50499b and HD 50499c. The results from the fully probed parameter space from the MCMC are shown in the appendix.
These parameters agree well with Vogt et al. (2005) who reports that HD 50499b has a period of 6.80± 0.30 yr and a mass of 1.71 ± 0.2 MJup.
Because not all of the outer planet period is covered, the uncertainties on the period of the second Keplerian remain relatively large with an unconstrained ω.
The parameters we obtain for HD 50499c agree with Barbato et al. (2018) who report that the quadratic trend corresponds to a planet with an orbital period of P ≥ 22.61 yr and a minimum mass of Msini ≥ 0.942 MJup.
![]() |
Fig. 6 Top: HD 50499 RV curves. Blue: CORALIE-98 data; green: CORALIE-07; purple: CORALIE-14 data; red: HARPS data; and orange: HIRES data (Butler et al. 2017). The Keplerian models are represented by black curves. Middle: RV residuals of HD 50499. Bottom: periodogram of the residuals of HD 50499 after the two planetary signals were removed, indicating that there are no more significant signals remaining in the data. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
![]() |
Fig. 7 Phase-folded RV curves for HD 50499. Blue: CORALIE-98 data; green: CORALIE-07; purple: CORALIE-14 data; red: HARPS data; orange: HIRES data. The Keplerian models are represented by black curves. Top: phase-folded curve for HD 50499b. Bottom: phase-folded curve for HD 50499c. |
![]() |
Fig. 8 Top: HD 92788 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple), HARPS (red), HIRES (orange; Butler et al. 2017) and HAMILTON (pink; Butler et al. 2006). The best single-planet Keplerian model is represented by black curve. Middle: RV residuals of HD 92788. Bottom: periodogram of the residuals for HD 92788 after thesignal was removed showing no significant signals. The three black lines represent the 10, 1, and 0.1% FAPs in ascending order. |
4.3 HD 92788 (HIP 52409)
We report updated orbital parameters for HD 92788b (Fischer et al. 2001) and confirm the detection of HD 92788c (Wittenmyer et al. 2019). HD 92788 is a G6V star at 28.83 ± 0.05 pc from the Sun. The star properties are summarised in Table 4.
HD 92788 was observed with CORALIE at La Silla Observatory since March 1999, 59 RV measurements were obtained with CORALIE-98, an additional 10 RV measurements were obtained with CORALIE-07, and an additional 11 RV measurements were obtained with CORALIE-14. There are also 61 measurements taken with HARPS, 42 RV points publicly available from HIRE, and an additional 31 RV points from HAMILTON.
The orbital solutions for HD 92788 are summarised in Table 4. Figure 8 shows the CORALIE, HARPS, HIRES, and HAMILTON RVs and the corresponding best-fit Keplerian models along with a time series of the residuals and a periodogram of the residuals. Figure 9 shows the phase-folded RV diagram for HD 92788b and the time series for HD 92788c. The results from the fully probed parameter space from the MCMC are shown in the appendix.
The orbital parameters for HD 92788b agree well with Fischer et al. (2001) who reported a period of P = 0.894 ± 0.009 yr and a minimum mass of Msini = 3.34 MJup and also agree well with the recently reported parameters by Wittenmyer et al. (2019) who report HD 92788b to have an orbital period of 0.892 ± 0.00008 yr and a minimum mass of Msini = 3.78 ± 0.18MJup.
Wittenmyer et al. (2013) tested the system for a potential additional planet with a period of 162 days which they found by limiting the eccentricity of the Keplerian model of HD 92788b. This is because fitting Keplerians can be biased towards fitting an increased eccentricity, especially when the semi-amplitude is small or the system has not been sampled well (Shen & Turner 2008). Although Wittenmyer et al. (2013) suggested that there is a possible additional planet in the system, when we fitted the Keplerian the signal was not significant above the noise (lower than a 10% FAP) to claim an additional planet at this period, as seen in Fig. 8 and therefore we do not detect this signal.
We do however detect the signal of another planet in the system, HD 92788c, as was recently discovered by Wittenmyer et al. (2019). We confirm this planet with orbital parameters that agree with Wittenmyer et al. (2019), who report a period of 26.99 ± 2.54 yr and a minimum mass of Msini of 3.64 ± 0.69 MJup.
![]() |
Fig. 9 Blue: CORALIE-98 data; green: CORALIE-07 data; purple: CORALIE-14 data; red: HARPS data; orange: HIRES data; pink: HAMILTON data. The Keplerian models are represented by black curves. Top: phase-folded curve for HD 92788b.Bottom: time series for HD 92788c. |
5 Discussion and conclusion
Here we report the discovery of five new giant planets and brown dwarf candidates discovered with the CORALIE spectrograph mounted on 1.2 m Euler Swiss telescope at La Silla Observatory as well as updated orbital parameters for four previously detected planets. CORALIE time series combined with the published data sets independently confirms the existence of these four already published companions. In addition, we do not find any significant evidence for the exoplanet HD 92788c suggested by Wittenmyer et al. (2013). The newly reported companions span a period range of 15.6–40.4 yr and a mass domain of 2.93–26.77 MJup.
Most of the parent stars in this paper have a metallicity excess, as seen in Fig. 10. Despite the small size of our sample, our results seem to agree with previous observations that giant planets appear to occur significantly around stars that are more metal-rich, which has been noted before by Santos et al. (2004), Fischer & Valenti (2005), Mayor et al. (2011) and Boisse et al. (2012).
The focus of this paper is on long-period exoplanets, where all of the newly reported planets have periods over 15 yr. This contributes to the relatively small number of previously known planets with periods in this range, where according to the NASA Exoplanet2 Archive3 there are only 26 known exoplanets with a period greater than 15 yr.
As seen in Fig. 11, the planets and brown dwarfs presented in this paper are bridging the gap between the RV detected exoplanets and the directly imaged exoplanets. As we achieve deeper detection limits and smaller inner working angles in imaging with new instrumentation and telescopes, and span longer base lines with RV techniques, this gap in separation will decrease.
Combining RV and direct imaging data has previously been done, for example by Kane et al. (2014), Rodigas et al. (2016a,b) and Crepp et al. (2012). More recently, these two techniques have been combined for two candidates from the CORALIE RV survey with the detection of HD 4747 B by Peretti et al. (2019) and with the discovery of an ultra-cool brown dwarf companion to HD 4113 A by Cheetham et al. (2018).
With future direct imaging detections of these candidates presented here, we aim to follow the same method as Cheetham et al. (2018) and Peretti et al. (2019), and perform an atmospheric retrieval analysis of carbon and oxygen abundances, as demonstrated by Lavie et al. (2017). The CORALIE survey plays an important step in identifying promising targets for such observations, making CORALIE a unique instrument in being able to carry out such a long continuous survey at high precision. Furthermore, the stars in the CORALIE sample are older than the commonly and directly imaged targets, andso the substellar companions probed in this paper represent a new and complementary parameter space.
A great many attempts have been made to detect these long-period companions through imaging as part of a 15-yr effort using VLT/NACO with little success. Now with the capabilities of VLT/SPHERE with a contrast limit of ~ 10−3–10−4 at a separation of 0.1 arcsec (Beuzit et al. 2019) we are able to start detecting these massive planets and brown dwarfs.
Some of these targets may be challenging for SPHERE, but where this is the case, they should be within the capabilities of ELT/METIS, where METIS should achieve a 10−5 contrast at 0.1 arcseconds (Carlomagno et al. 2016). In addition, combining astrometry that will be available from Gaia with RV data will allow us to further constrain the range of possible masses of these massive companions.
![]() |
Fig. 10 Companion mass in the limits 1–50 MJup as a function of the host star metallicity. The new companions presented in this paper are shown by the orange stars. Previously detected planets and brown dwarfs are shown in blue2. The black dashed line shows the metallicity of the Sun. Most of the stars with detected companions in this paper have a significant metallicity excess. |
![]() |
Fig. 11 Mass of detected exoplanets and brown dwarfs as a function of separation. The new companions presented in this paper are shown by the orange stars. Previously detected imaged planets and brown dwarfs are shown in blue and the RV-detected planets are shown in red2. The black dashed lines show the limits for companions that are in a potential detectable parameter space with imaging, i.e., giant planets more massive than 2 MJup and planets at a separation larger than 5 AU. |
Acknowledgements
This work has been carried out within the framework of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation. The authors acknowledge the financial support of the SNSF. This publication makes use of the Data & Analysis Center for Exoplanets (DACE), which is a facility based at the University of Geneva (CH) dedicated to extrasolar planets data visualisation, exchange and analysis. DACE is a platform of the Swiss National Centre of Competence in Research (NCCR) PlanetS, federating the Swiss expertise in Exoplanet research. The DACE platform is available at https://dace.unige.ch. L.A.dS. and J.V.S. acknowledge the support from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (project FOUR ACES, grant agreement No 724427). N.C.S. was supported by FCT – Fundação para a Ciência e a Tecnologia – through national funds and by FEDER through COMPETE2020 – Programa Operacional Competitividade e Internacionalização by these grants: UID/FIS/04434/2013 & POCI-01-0145-FEDER-007672; PTDC/FIS-AST/28953/2017 & POCI-01-0145-FEDER-028953 and PTDC/FIS-AST/32113/2017 & POCI-01-0145-FEDER-032113. This work has made use of data from the European Space Agency (ESA) mission Gaia (http://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, http://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.
Appendix A Direct access to the RVs and other data products
The RV measurements and additional data products discussed in this paper are available in electronic form on the DACE web platform for each individual target with each link. A copy of the data is also available at the CDS.
HD 13724: https://dace.unige.ch/radialVelocities/?pattern=HD13724
HD 181234: https://dace.unige.ch/radialVelocities/?pattern=HD181234
HD 25015: https://dace.unige.ch/radialVelocities/?pattern=HD25015
HD 50499: https://dace.unige.ch/radialVelocities/?pattern=HD50499
HD 92788: https://dace.unige.ch/radialVelocities/?pattern=HD92788
HD 92987: https://dace.unige.ch/radialVelocities/?pattern=HD92987
HD 98649: https://dace.unige.ch/radialVelocities/?pattern=HD98649
Appendix B MCMC tables
We probed the model parameter space with the MCMC sampler described in (Díaz et al. 2014, 2016). For each MCMC simulation, we performed 10 000 000 iterations with initial conditions drawn from the solution obtained using DACE. The corresponding parameters and confidence intervals for each star are listed in Tables A.1–B.6 for HD 181234, HD 92987, HD 13724, HD 25015, HD 98649, HD 50499 and HD 92788 respectively.
Parameters probed by the MCMC used to fit the RV measurements of HD 181234.
Parameters probed by the MCMC used to fit the RV measurements of HD 92987.
Parameters probed by the MCMC used to fit the RV measurements of HD 13724.
Parameters probed by the MCMC used to fit the RV measurements of HD 25015.
Parameters probed by the MCMC used to fit the RV measurements of HD 98649.
Parameters probed by the MCMC used to fit the RV measurements of HD 50499.
Parameters probed by the MCMC used to fit the RV measurements of HD 92788.
References
- Baluev, R. V. 2008, MNRAS, 385, 1279 [NASA ADS] [CrossRef] [Google Scholar]
- Barbato, D., Sozzetti, A., Desidera, S., et al. 2018, A&A, 615, A175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Beuzit, J.-L., Vigan, A., Mouillet, D., et al. 2019, A&A, submitted, [arXiv:1902.04080] [Google Scholar]
- Boisse, I., Pepe, F., Perrier, C., et al. 2012, A&A, 545, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bond, J. C., Tinney, C. G., Butler, R. P., et al. 2006, MNRAS, 370, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bowler, B. P. 2016, PASP, 128, 102001 [NASA ADS] [CrossRef] [Google Scholar]
- Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505 [NASA ADS] [CrossRef] [Google Scholar]
- Butler, R. P., Vogt, S. S., Laughlin, G., et al. 2017, AJ, 153, 208 [NASA ADS] [CrossRef] [Google Scholar]
- Carlomagno, B., Absil, O., Kenworthy, M., et al. 2016, in Adaptive Optics Systems V, Proc. SPIE, 9909, 990973 [CrossRef] [Google Scholar]
- Cheetham, A., Ségransan, D., Peretti, S., et al. 2018, A&A, 614, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Crepp, J. R., Johnson, J. A., Howard, A. W., et al. 2012, ApJ, 761, 39 [NASA ADS] [CrossRef] [Google Scholar]
- Crepp, J. R., Principe, D. A., Wolff, S., et al. 2018, ApJ, 853, 192 [NASA ADS] [CrossRef] [Google Scholar]
- Delisle, J. B., Ségransan, D., Buchschacher, N., & Alesina, F. 2016, A&A, 590, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Delisle, J.-B., Ségransan, D., Dumusque, X., et al. 2018, A&A, 614, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Díaz, R. F., Almenara, J. M., Santerne, A., et al. 2014, MNRAS, 441, 983 [NASA ADS] [CrossRef] [Google Scholar]
- Díaz, R. F., Ségransan, D., Udry, S., et al. 2016, A&A, 585, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Feigelson, E. D., Lawson, W. A., Stark, M., Townsley, L., & Garmire, G. P. 2006, AJ, 131, 1730 [NASA ADS] [CrossRef] [Google Scholar]
- Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 [NASA ADS] [CrossRef] [Google Scholar]
- Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2001, ApJ, 551, 1107 [NASA ADS] [CrossRef] [Google Scholar]
- Fischer, D. A., Howard, A. W., Laughlin, G. P., et al. 2014, Protostars and Planets VI (Tucson: University of Arizona Press), 715 [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Georgy, C., Ekström, S., Eggenberger, P., et al. 2013, A&A, 558, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gomes da Silva, J., Santos, N. C., Boisse, I., Dumusque, X., & Lovis, C. 2014, A&A, 566, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kane, S. R., Howell, S. B., Horch, E. P., et al. 2014, ApJ, 785, 93 [NASA ADS] [CrossRef] [Google Scholar]
- Lavie, B., Mendonça, J. M., Mordasini, C., et al. 2017, AJ, 154, 91 [NASA ADS] [CrossRef] [Google Scholar]
- Mamajek, E. E., & Hillenbrand, L. A. 2008, ApJ, 687, 1264 [NASA ADS] [CrossRef] [Google Scholar]
- Marmier, M. 2014, Ph.D. Thesis, Geneva Observatory, University of Geneva, Switzerland [Google Scholar]
- Marmier, M., Ségransan, D., Udry, S., et al. 2013, A&A, 551, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
- Mayor, M., Marmier, M., Lovis, C., et al. 2011, ArXiv e-prints [arXiv:1109.2497] [Google Scholar]
- Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Peretti, S., Ségransan, D., Lavie, B., et al. 2019, A&A, accepted [arXiv:1805.05645] [Google Scholar]
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
- Porto de Mello, G. F., da Silva, R., da Silva, L., & de Nader, R. V. 2014, A&A, 563, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rodigas, T. J., Bergeron, P., Simon, A., et al. 2016a, ApJ, 831, 177 [NASA ADS] [CrossRef] [Google Scholar]
- Rodigas, T. J., Arriagada, P., Faherty, J., et al. 2016b, ApJ, 818, 106 [NASA ADS] [CrossRef] [Google Scholar]
- Santos, N. C., Israelian, G., & Mayor, M. 2001, A&A, 373, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Sousa, S. G., Mortier, A., et al. 2013, A&A, 556, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ségransan, D., Udry, S., Mayor, M., et al. 2010, A&A, 511, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Shen, Y., & Turner, E. L. 2008, ApJ, 685, 553 [NASA ADS] [CrossRef] [Google Scholar]
- Sousa, S. G., Santos, N. C., Mayor, M., et al. 2008, A&A, 487, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Tamuz, O., Ségransan, D., Udry, S., et al. 2008, A&A, 480, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Udry, S., Mayor, M., Naef, D., et al. 2000, A&A, 356, 590 [NASA ADS] [Google Scholar]
- Udry, S., Mayor, M., Naef, D., et al. 2002, A&A, 390, 267 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vogt, S. S. 1987, PASP, 99, 1214 [NASA ADS] [CrossRef] [Google Scholar]
- Vogt, S. S., Allen, S. L., Bigelow, B. C., et al. 1994, in Instrumentation in Astronomy VIII, eds. D. L. Crawford, & E. R. Craine, Proc. SPIE, 2198, 362 [Google Scholar]
- Vogt, S. S., Butler, R. P., Marcy, G. W., et al. 2005, ApJ, 632, 638 [NASA ADS] [CrossRef] [Google Scholar]
- Wittenmyer, R. A., Wang, S., Horner, J., et al. 2013, ApJS, 208, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Wittenmyer, R. A., Clark, J. T., Zhao, J., et al. 2019, MNRAS, 484, 5859 [NASA ADS] [CrossRef] [Google Scholar]
- Zuckerman, B., Song, I., Bessell, M. S., & Webb, R. A. 2001, ApJ, 562, L87 [NASA ADS] [CrossRef] [Google Scholar]
- Zuckerman, B., Rhee, J. H., Song, I., & Bessell, M. S. 2011, ApJ, 732, 61 [NASA ADS] [CrossRef] [Google Scholar]
The data are available at the Data and Analysis Center for Exoplanets (DACE) which can be accessed at https://dace.unige.ch
The NASA Exoplanet Archive can be accessed at https://exoplanetarchive.ipac.caltech.edu/index.html
All Tables
Observed and inferred stellar parameters for host stars to the planet candidates: HD 181234, HD 25015, HD 50499, HD 92788, and HD 98649.
Observed and inferred stellar parameters for host stars to the brown dwarf candidates: HD 13724 and HD 92987.
Best-fitted solutions for the substellar companions orbiting HD 13724, HD 181234, HD 25015, and HD 92987.
Best-fitted solution for the substellar companions orbiting HD 50499, HD 92788, and HD 98649.
All Figures
![]() |
Fig. 1 Top: HD 181234 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple) and HIRES data (Butler et al. 2017) in orange. The bestsingle-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 181234. Bottom: periodogram of the residuals for HD 181234. The three black lines represent the 10, 1, and 0.1% FAPs in ascending order. |
In the text |
![]() |
Fig. 2 Top: HD 92987 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The fitted single-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 92987. Bottom: periodogram of the residuals for HD 92987. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
In the text |
![]() |
Fig. 3 Top: HD 25015 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The fitted single-planet Keplerian model is represented as a black curve. Second figure: RV residuals of HD 25015. Third figure: periodogram of the residuals for HD 25015 after the signal has been removed showing no significant signals. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. Bottom: periodogram of Hα before detrending. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order, showing a significant peak above the 10% FAP at 370.26 days. |
In the text |
![]() |
Fig. 4 Top: HD 13724 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple) and HARPS (red). The best single-planet Keplerian model is representedas a black curve. Middle: RV residuals of HD 13724. Bottom: periodogram of the residuals for HD 13724 after the signal has been removed showing no significant signals. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
In the text |
![]() |
Fig. 5 Top: HD 98649 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green) and CORALIE-14 (purple). The best single-planet Keplerian model is represented as a black curve. Middle: RV residuals of HD 98649. Bottom: periodogram of the residuals for HD 98649 showing no significant signals remaining. The three black lines represent the 10, 1 and 0.1% FAP in ascending order. |
In the text |
![]() |
Fig. 6 Top: HD 50499 RV curves. Blue: CORALIE-98 data; green: CORALIE-07; purple: CORALIE-14 data; red: HARPS data; and orange: HIRES data (Butler et al. 2017). The Keplerian models are represented by black curves. Middle: RV residuals of HD 50499. Bottom: periodogram of the residuals of HD 50499 after the two planetary signals were removed, indicating that there are no more significant signals remaining in the data. The three black lines represent the 10, 1 and 0.1% FAPs in ascending order. |
In the text |
![]() |
Fig. 7 Phase-folded RV curves for HD 50499. Blue: CORALIE-98 data; green: CORALIE-07; purple: CORALIE-14 data; red: HARPS data; orange: HIRES data. The Keplerian models are represented by black curves. Top: phase-folded curve for HD 50499b. Bottom: phase-folded curve for HD 50499c. |
In the text |
![]() |
Fig. 8 Top: HD 92788 RV measurements as a function of Julian Date obtained with CORALIE-98 (blue), CORALIE-07 (green), CORALIE-14 (purple), HARPS (red), HIRES (orange; Butler et al. 2017) and HAMILTON (pink; Butler et al. 2006). The best single-planet Keplerian model is represented by black curve. Middle: RV residuals of HD 92788. Bottom: periodogram of the residuals for HD 92788 after thesignal was removed showing no significant signals. The three black lines represent the 10, 1, and 0.1% FAPs in ascending order. |
In the text |
![]() |
Fig. 9 Blue: CORALIE-98 data; green: CORALIE-07 data; purple: CORALIE-14 data; red: HARPS data; orange: HIRES data; pink: HAMILTON data. The Keplerian models are represented by black curves. Top: phase-folded curve for HD 92788b.Bottom: time series for HD 92788c. |
In the text |
![]() |
Fig. 10 Companion mass in the limits 1–50 MJup as a function of the host star metallicity. The new companions presented in this paper are shown by the orange stars. Previously detected planets and brown dwarfs are shown in blue2. The black dashed line shows the metallicity of the Sun. Most of the stars with detected companions in this paper have a significant metallicity excess. |
In the text |
![]() |
Fig. 11 Mass of detected exoplanets and brown dwarfs as a function of separation. The new companions presented in this paper are shown by the orange stars. Previously detected imaged planets and brown dwarfs are shown in blue and the RV-detected planets are shown in red2. The black dashed lines show the limits for companions that are in a potential detectable parameter space with imaging, i.e., giant planets more massive than 2 MJup and planets at a separation larger than 5 AU. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.