Massive Giant Planets Research Articles

Revisiting planetary systems in the Okayama Planet Search Program: A new long-period planet, RV astrometry joint analysis, and a multiplicity-metallicity trend around evolved stars

Abstract In this study, we revisit 32 planetary systems around evolved stars observed within the framework of the Okayama Planet Search Program (OPSP) and its collaborative framework of the East Asian Planet Search Network to search for additional companions and investigate the properties of stars and giant planets in multiple-planet systems. With our latest radial velocities obtained from Okayama Astrophysical Observatory, we confirm an additional giant planet in the wide orbit of the 75 Cet system ($P_{\rm {c}} = 2051.62_{-40.47}^{+45.98}\ \rm {d}$, $M_{\rm {c}}\sin i=0.912_{-0.090}^{+0.088}\,\,M_{\rm {J}}$, and $a_{\rm {c}}=3.929_{-0.058}^{+0.052}\ \rm {au}$), along with five stars exhibiting long-term radial velocity accelerations, which indicates massive companions in the wide orbits. We also find that the radial velocity variations of several planet-harboring stars may indicate additional planet candidates, stellar activities, or other understudied sources. These stars include ϵ Tau, 11 Com, 24 Boo, 41 Lyn, 14 And, HD 32518, and ω Ser. We further constrain the orbital configuration of the HD 5608, HD 14067, HD 120084, and HD 175679 systems by combining radial velocities with astrometry, as their host central stars exhibit significant astrometric accelerations. For other systems, we simply refine their orbital parameters. Moreover, our study indicates that the OPSP planet-harboring stars are more metal poor compared to the currently known planet-harboring stars, and this is likely due to the B − V color upper limit at 1.0 for star selection in the beginning of the survey. Finally, by investigating the less massive giant planets (<5 MJ) around currently known planet-harboring evolved stars, we find that metallicity positively correlates with the multiplicity and total planet mass of the system, which can be evidence for the core-accretion planet formation model.

An inflationary disk phase to explain extended protoplanetary dust disks

Context. Understanding planetesimal formation is an essential first step towards understanding planet formation. The distribution of these first solid bodies drives the locations where planetary embryos, which eventually form fully-fledged planets, grow. Aims. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at least the following three criteria: (1) It must produce an extended gas and dust disk (e.g. 45 au for the dust); (2) within the disk, the local dust-to-gas ratio in at least two distinct locations must sufficiently increase to explain the early formation of the parent bodies of non-carbonaceous and carbonaceous iron meteorites; and (3) dust particles, which have condensed at high temperatures (i.e. calcium–aluminium-rich inclusions), must be transported to the outer disk. Though current protoplanetary disk models are able to satisfy one or two of these criteria, none have been successful in recreating all three. We aim to find scenarios that satisfy all three. Methods. In this study we used a 1D disk model that tracks the evolution of the gas and dust disks. Planetesimals are formed within the disk at locations where the streaming instability can be triggered. We explored a large parameter space to study the effect of the disk viscosity, the timescale of infall of material into the disk, the distance within which material is deposited into the disk, and the fragmentation threshold of dust particles. Results. We find that scenarios with a large initial disk viscosity (α > 0.05), a relatively short infall timescale (Tinfall < 100–200kyr), and a small centrifugal radius (RC ~ 0.4 au; i.e. the distance within which material falls into the disk) result in disks that satisfy all three criteria needed to represent the protoplanetary disk of the Solar System. The large initial viscosity and short infall timescale result in a rapid initial expansion of the disk, which we dub the ‘inflationary phase’ of the disk. Furthermore, a temperature-dependent fragmentation threshold, which accounts for cold icy particles breaking more easily, results in larger and more massive disks. This, in turn, results in more ‘icy’ than ‘rocky’ planetesimals. Such scenarios are also better in line with our Solar System, which has small terrestrial planets and massive giant planet cores. Finally, we find that scenarios with large RC cannot transport calcium–aluminium-rich inclusions to the outer disk and do not produce planetesimals at two locations within the disk.

Preparation for an unsupervised massive analysis of SPHERE high-contrast data with PACO

Context. Despite tremendous progress in the detection and characterization of extrasolar planetary systems in the last 25 yr, we have not pinpointed any Solar System analogues. In particular, Jupiter-like planets (either mature or old) are barely detectable beyond 5 au with indirect techniques and they are still out of the reach of direct imaging techniques. Aims. Our study is aimed at a search for exoplanets throughout the whole ESO/VLT-SPHERE archive with an improved and unsupervised data analysis algorithm that could allow us to detect massive giant planets at 5 au. To prepare, test, and optimize our approach, we gathered a sample of 24 solar-type stars observed with SPHERE using angular and spectral differential imaging modes. Methods. We used PACO, a recently developed new-generation algorithm that has been shown to outperform classical methods. We also improved the SPHERE pre-reduction pipeline and optimized the outputs of PACO to enhance the detection performance. We developed custom-built spectral prior libraries to optimize the detection capability of the ASDI mode for both IRDIS and IFS. Results. Compared to previous works conducted with more classical algorithms, the contrast limits we derived with PACO are more reliable and significantly improved, especially at short angular separations, where a gain by a factor ten has been obtained between 0.2 and 0.5 arcsec. Under good observing conditions, planets down to 5 MJup, orbiting at 5 au could be detected around stars within 60 parsec. We identified two exoplanet candidates that will require a follow-up to test for a common proper motion. Conclusions. In this work, we use a small sample to demonstrate the benefits of PACO in terms of achievable contrast and of control of the confidence levels. In addition, we have developed custom tools to take full advantage of this algorithm and to quantity the total error budget on the estimated astrometry and photometry. This work paves the way towards an end-to-end, homogeneous, and unsupervised massive re-reduction of archival direct imaging surveys in the quest for new exo-Jupiters.

Discovery of a massive giant planet with extreme density around the sub-giant star TOI-4603

We present the discovery of a transiting massive giant planet around TOI-4603, a sub-giant F-type star from NASA’s Transiting Exoplanet Survey Satellite (TESS). The newly discovered planet has a radius of 1.042−0.035+0.038 RJand an orbital period of 7.24599−0.00021+0.00022days. Using radial velocity measurements with the PARAS and TRES spectrographs, we determined the planet’s mass to be 12.89−0.57+0.58 MJ, resulting in a bulk density of 14.1−1.6+1.7g cm−3. This makes it one of the few known massive giant planets with an extreme density. It lies in the transition mass region of massive giant planets and low-mass brown dwarfs, an important addition to the population of fewer than five known objects in this mass range. The eccentricity of 0.325 ± 0.020 and an orbital separation of 0.0888 ± 0.0010 AU from its host star suggest that the planet is likely undergoing high-eccentricity tidal migration. We find a fraction of heavy elements of 0.13−0.06+0.05and metal enrichment of the planet (ZP/Zstar) of 4.2−2.0+1.6. Detection of such systems will enable us to gain valuable insights into the governing mechanisms of massive planets and improve our understanding of their dominant formation and migration mechanisms.

TOI-2525 b and c: A Pair of Massive Warm Giant Planets with Strong Transit Timing Variations Revealed by TESS* ∗ Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under MPG program 0104.A-9007. This paper includes data gathered with the 6.5 m Magellan Telescopes located at Las Campanas Observatory, Chile.

The K-type star TOI-2525 has an estimated mass of M = M ⊙ and radius of R = R ⊙ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semiamplitude of ∼6 hr. We performed TTV dynamical and photodynamical light-curve analysis of the TESS data combined with radial velocity measurements from FEROS and PFS, and we confirmed the planetary nature of these companions. The TOI-2525 system consists of a transiting pair of planets comparable to Neptune and Jupiter with estimated dynamical masses of m b = and m c = M Jup, radii of r b = and r c = R Jup, and orbital periods of P b = and P c = days for the inner and outer planet, respectively. The period ratio is close to the 2:1 period commensurability, but the dynamical simulations of the system suggest that it is outside the mean-motion resonance (MMR) dynamical configuration. Object TOI-2525 b is among the lowest-density Neptune-mass planets known to date, with an estimated median density of ρ b = g cm−3. The TOI-2525 system is very similar to the other K dwarf systems discovered by TESS, TOI-2202 and TOI-216, which are composed of almost identical K dwarf primaries and two warm giant planets near the 2:1 MMR.

Gamma-ray flux limits from brown dwarfs: Implications for dark matter annihilating into long-lived mediators

Brown dwarfs (BDs) are celestial objects representing the link between the least massive main-sequence stars and giant gas planets. In the first part of this article, we perform a model-independent search of a gamma-ray signal from the direction of nine nearby BDs in 13 years of \Fermi-LAT data. We find no significant excess of gamma rays, and we, therefore, set 95\% confidence level upper limits on the gamma-ray flux with a binned-likelihood approach. In the second part of the paper, we interpret these bounds within an exotic mechanism proposed for gamma-ray production in BDs: If the dark matter (DM) of the universe is constituted of particles with non-negligible couplings to the standard model, BDs may efficiently accumulate them through scatterings. DM particles eventually thermalize, and can annihilate into light, long-lived, mediators which later decay into photons outside the BD. Within this framework, we set a stacked upper limit on the DM-nucleon elastic scattering cross section at the level $\sim 10^{-38}$ cm$^{2}$ for DM masses below 10 GeV. Our limits are comparable to similar bounds from the capture of DM particles in celestial objects, but have the advantage of covering a larger portion of the parameter space in mediator decay length and DM mass. They also depend only on the local DM abundance, as opposed to the inner Galaxy profile, and are thus more robust.

A mixed stirring mechanism for debris discs with giant and dwarf planetary perturbations

ABSTRACT Debris discs consist of belts of bodies ranging in size from dust grains to planetesimals; these belts are visible markers of planetary systems around other stars that can reveal the influence of extrasolar planets through their shape and structure. Two key stirring mechanisms – self-stirring by planetesimals and secular perturbation by an external giant planet – have been identified to explain the dynamics of planetesimal belts; their relative importance has been studied independently, but are yet to be considered in combination. In this work, we perform a suite of 286 N-body simulations exploring the evolution of debris discs over 1 Gyr, combining the gravitational perturbations of both dwarf planets embedded in the discs, and an interior giant planet. Our systems were somewhat modelled after the architecture of the outer Solar system: a Solar mass star, a single massive giant planet at 30 au (MGP = 10 to 316 M⊕), and a debris disc formed by 100 massive dwarf planets and 1000 massless particles (MDD = 3.16 to 31.6 M⊕). We present the evolution of both the disc and the giant planet after 1 Gyr. The time evolution of the average eccentricity and inclination of the disc is strongly dependent on the giant planet mass as well as on the remaining disc mass. We also found that efficient stirring is achieved even with small disc masses. In general, we find that a mixed mechanism is more efficient in the stirring of cold debris discs than either mechanism acting in isolation.

The interplay between forming planets and photoevaporating discs I: forbidden line diagnostics

ABSTRACT Disc winds and planet formation are considered to be two of the most important mechanisms that drive the evolution and dispersal of protoplanetary discs and in turn define the environment in which planets form and evolve. While both have been studied extensively in the past, we combine them into one model by performing three-dimensional radiation-hydrodynamic simulations of giant planet hosting discs that are undergoing X-ray photoevaporation, with the goal to analyse the interactions between both mechanisms. In order to study the effect on observational diagnostics, we produce synthetic observations of commonly used wind-tracing forbidden emission lines with detailed radiative transfer and photoionization calculations. We find that a sufficiently massive giant planet carves a gap in the gas disc that is deep enough to affect the structure and kinematics of the pressure-driven photoevaporative wind significantly. This effect can be strong enough to be visible in the synthetic high-resolution observations of some of our wind diagnostic lines, such as the [O i] 6300 Å or [S ii] 6730 Å lines. When the disc is observed at inclinations around 40° and higher, the spectral line profiles may exhibit a peak in the redshifted part of the spectrum, which cannot easily be explained by simple wind models alone. Moreover, massive planets can induce asymmetric substructures within the disc and the photoevaporative wind, giving rise to temporal variations of the line profiles that can be strong enough to be observable on time-scales of less than a quarter of the planet’s orbital period.

A quarter century of spectroscopic monitoring of the nearby M dwarf Gl 514

Context. Statistical analyses based on Kepler data show that most of the early-type M dwarfs host multi-planet systems consisting of Earth- to sub-Neptune-sized planets with orbital periods of up to ~250 days, and that at least one such planet is likely located within the habitable zone. M dwarfs are therefore primary targets to search for potentially habitable planets in the solar neighbourhood. Aims. We investigated the presence of planetary companions around the nearby (7.6 pc) and bright (V = 9 mag) early-type M dwarf Gl 514, analysing 540 radial velocities collected over nearly 25 yr with the HIRES, HARPS, and CARMENES spectrographs. Methods. The data are affected by time-correlated signals at the level of 2–3 m s−1 due to stellar activity, which we filtered out, testing three different models based on Gaussian process regression. As a sanity cross-check, we repeated the analyses using HARPS radial velocities extracted with three different algorithms. We used HIRES radial velocities and Hipparcos-Gaia astrometry to put constraints on the presence of long-period companions, and we analysed TESS photometric data. Results. We find strong evidence that Gl 514 hosts a super-Earth on a likely eccentric orbit, residing in the conservative habitable zone for nearly 34% of its orbital period. The planet Gl 514b has minimum mass mb sin ib = 5.2 ± 0.9 M⊕, orbital period Pb = 140.43 ± 0.41 days, and eccentricity eb = 0.45−0.14+0.15. No evidence for transits is found in the TESS light curve. There is no evidence for a longer period companion in the radial velocities and, based on astrometry, we can rule out a ~0.2 MJup planet at a distance of ~3–10 astronomical units, and massive giant planets and brown dwarfs out to several tens of astronomical units. We discuss the possible presence of a second low-mass companion at a shorter distance from the host than Gl 514 b. Conclusions. Gl 514 b represents an interesting science case for studying the habitability of planets on eccentric orbits. We advocate for additional spectroscopic follow-up to get more accurate and precise planetary parameters. Further follow-up is also needed to investigate the presence of additional planetary signals of less than 1 m s−1.

The California Legacy Survey. III. On the Shoulders of (Some) Giants: The Relationship between Inner Small Planets and Outer Massive Planets

We use a high-precision radial velocity survey of FGKM stars to study the conditional occurrence of two classes of planets: close-in small planets (0.023–1 au, 2–30 M ⊕) and distant giant planets (0.23–10 au, 30–6000 M ⊕). We find that of systems with a close-in, small planet also host an outer giant, compared to for stars irrespective of small planet presence. This implies that small planet hosts may be enhanced in outer giant occurrences compared to all stars with 1.7σ significance. Conversely, we estimate that of cold giant hosts also host an inner small planet, compared to of stars irrespective of cold giant presence. We also find that more massive and close-in giant planets are not associated with small inner planets. Specifically, our sample indicates that small planets are less likely to have outer giant companions more massive than approximately 120 M ⊕ and within 0.3–3 au, than to have less massive or more distant giant companions, with ∼2.2σ confidence. This implies that massive gas giants within 0.3–3 au may suppress inner small planet formation. Additionally, we compare the host-star metallicity distributions for systems with only small planets and those with both small planets and cold giants. In agreement with previous studies, we find that stars in our survey that only host small planets have a metallicity distribution that is consistent with the broader solar-metallicity-median sample, while stars that host both small planets and gas giants are distinctly metal rich with ∼2.3σ confidence.

TOI-530b: a giant planet transiting an M-dwarf detected by TESS

ABSTRACT We report the discovery of TOI-530b, a transiting Saturn-like planet around an M0.5V dwarf, delivered by the Transiting Exoplanet Survey Satellite (TESS). The host star is located at a distance of 147.7 ± 0.6 pc with a radius of R* = 0.54 ± 0.03 R⊙ and a mass of M* = 0.53 ± 0.02 M⊙. We verify the planetary nature of the transit signals by combining ground-based multiwavelength photometry, high-resolution spectroscopy from SPIRou as well as high-angular-resolution imaging. With V = 15.4 mag, TOI-530b is orbiting one of the faintest stars accessible by ground-based spectroscopy. Our model reveals that TOI-530b has a radius of 0.83 ± 0.05 RJ and a mass of 0.37 ± 0.08 MJ on a 6.39-d orbit. TOI-530b is the sixth transiting giant planet hosted by an M-type star, which is predicted to be infrequent according to core accretion theory, making it a valuable object to further study the formation and migration history of similar planets. Furthermore, we identify a potential dearth of hot massive giant planets around M-dwarfs with separation distance smaller than 0.1 au and planet-to-star mass ratio between 2 × 10−3 and 10−2. We also find a possible correlation between hot giant planet formation and the metallicity of its parent M-dwarf. We discuss the potential formation channel of such systems.

Two Massive Jupiters in Eccentric Orbits from the TESS Full-frame Images

Abstract We report the discovery of two short-period massive giant planets from NASA’s Transiting Exoplanet Survey Satellite (TESS). Both systems, TOI-558 (TIC 207110080) and TOI-559 (TIC 209459275), were identified from the 30 minute cadence full-frame images and confirmed using ground-based photometric and spectroscopic follow-up observations from TESS’s follow-up observing program working group. We find that TOI-558 b, which transits an F-dwarf (M * = 1.349 − 0.065 + 0.064 M ⊙, R * = 1.496 − 0.040 + 0.042 R ⊙, T eff = 6466 − 93 + 95 K, age 1.79 − 0.73 + 0.91 Gyr) with an orbital period of 14.574 days, has a mass of 3.61 ± 0.15 M J, a radius of 1.086 − 0.038 + 0.041 R J, and an eccentric (e = 0.300 − 0.020 + 0.022 ) orbit. TOI-559 b transits a G dwarf (M * = 1.026 ± 0.057 M ⊙, R * = 1.233 − 0.026 + 0.028 R ⊙, T eff = 5925 − 76 + 85 K, age 6.8 − 2.0 + 2.5 Gyr) in an eccentric (e = 0.151 ± 0.011) 6.984 days orbit with a mass of 6.01 − 0.23 + 0.24 M J and a radius of 1.091 − 0.025 + 0.028 R J. Our spectroscopic follow up also reveals a long-term radial velocity trend for TOI-559, indicating a long-period companion. The statistically significant orbital eccentricity measured for each system suggests that these planets migrated to their current location through dynamical interactions. Interestingly, both planets are also massive (>3 M J), adding to the population of massive giant planets identified by TESS. Prompted by these new detections of high-mass planets, we analyzed the known mass distribution of hot and warm Jupiters but find no significant evidence for multiple populations. TESS should provide a near magnitude-limited sample of transiting hot Jupiters, allowing for future detailed population studies.

KMT-2019-BLG-0371 and the Limits of Bayesian Analysis

Abstract We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of q ∼ 0.08. Due to the short event duration (t E ∼ 6.5 days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass M h = 0.09 − 0.05 + 0.14 M ⊙ and a massive giant planet with a mass M p = 7.70 − 3.90 + 11.34 M Jup . However, the interpretation of the secondary as a planet (i.e., as having M p < 13M Jup) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of θ E and are relatively insensitive to other factors such as the direction of the event (ℓ, b), the lens–source relative proper motion μ rel, or the specific Galactic model prior.

NGTS-13b: a hot 4.8 Jupiter-mass planet transiting a subgiant star

We report the discovery of the massive hot Jupiter NGTS-13b by the Next Generation Transit Survey (NGTS). The V = 12.7 host star is likely in the subgiant evolutionary phase with logg* = 4.04 ± 0.05, Teff = 5819 ± 73 K, M* = 1.30−0.18+0.11 M⊙, and R* = 1.79 ± 0.06 R⊙. The NGTS detected a transiting planet with a period of P = 4.12 days around the star, which was later validated with the Transiting Exoplanet Survey Satellite (TESS; TIC 454069765). We confirm the planet using radial velocities from the CORALIE spectrograph. Using NGTS and TESS full-frame image photometry combined with CORALIE radial velocities, we determine NGTS-13b to have a radius of RP = 1.142 ± 0.046 RJup, a mass of MP = 4.84 ± 0.44 MJup, and an eccentricity of e = 0.086 ± 0.034. Previous studies have suggested that ~4 MJup may be the border separating two formation scenarios (e.g., core accretion and disk instability) and that massive giant planets share similar formation mechanisms as lower-mass brown dwarfs. NGTS-13b is just above 4 MJup, making it an important addition to the statistical sample needed to understand the differences between various classes of substellar companions. The high metallicity of NGTS-13, [Fe/H] = 0.25 ± 0.17, does not support previous suggestions that massive giants are found preferentially around lower metallicity host stars, but NGTS-13b does support findings that more massive and evolved hosts may have a higher occurrence of close-in massive planets than lower-mass unevolved stars.

Formation of intermediate-mass planets via magnetically controlled disk fragmentation

Intermediate mass planets, from Super-Earth to Neptune-sized bodies, are the most common type of planets in the galaxy. The prevailing theory of planet formation, core-accretion, predicts significantly fewer intermediate-mass giant planets than observed. The competing mechanism for planet formation, disk instability, can produce massive gas giant planets on wide-orbits, such as HR8799, by direct fragmentation of the protoplanetary disk. Previously, fragmentation in magnetized protoplanetary disks has only been considered when the magneto-rotational instability is the driving mechanism for magnetic field growth. Yet, this instability is naturally superseded by the spiral-driven dynamo when more realistic, non-ideal MHD conditions are considered. Here we report on MHD simulations of disk fragmentation in the presence of a spiral-driven dynamo. Fragmentation leads to the formation of long-lived bound protoplanets with masses that are at least one order of magnitude smaller than in conventional disk instability models. These light clumps survive shear and do not grow further due to the shielding effect of the magnetic field, whereby magnetic pressure stifles local inflow of matter. The outcome is a population of gaseous-rich planets with intermediate masses, while gas giants are found to be rarer, in qualitative agreement with the observed mass distribution of exoplanets.

The eccentricity distribution of giant planets and their relation to super-Earths in the pebble accretion scenario

Observations of the population of cold Jupiter planets (r >1 AU) show that nearly all of these planets orbit their host star on eccentric orbits. For planets up to a few Jupiter masses, eccentric orbits are thought to be the outcome of planet–planet scattering events taking place after gas dispersal. We simulated the growth of planets via pebble and gas accretion as well as the migration of multiple planetary embryos in their gas disc. We then followed the long-term dynamical evolution of our formed planetary system up to 100 Myr after gas disc dispersal. We investigated the importance of the initial number of protoplanetary embryos and different damping rates of eccentricity and inclination during the gas phase for the final configuration of our planetary systems. We constrained our model by comparing the final dynamical structure of our simulated planetary systems to that of observed exoplanet systems. Our results show that the initial number of planetary embryos has only a minor impact on the final orbital eccentricity distribution of the giant planets, as long as the damping of eccentricity and inclination is efficient. If the damping is inefficient (slow), systems with a larger initial number of embryos harbour larger average eccentricities. In addition, for slow damping rates, we observe that scattering events are already common during the gas disc phase and that the giant planets that formed in these simulations match the observed giant planet eccentricity distribution best. These simulations also show that massive giant planets (above Jupiter mass) on eccentric orbits are less likely to host inner super-Earths as they get lost during the scattering phase, while systems with less massive giant planets on nearly circular orbits should harbour systems of inner super-Earths. Finally, our simulations predict that giant planets are not single, on average, but they live in multi-planet systems.

Unveiling the β Pictoris system, coupling high contrast imaging, interferometric, and radial velocity data

Context. The nearby and young β Pictoris system hosts a well resolved disk, a directly imaged massive giant planet orbiting at ≃9 au, as well as an inner planet orbiting at ≃2.7 au, which was recently detected through radial velocity (RV). As such, it offers several unique opportunities for detailed studies of planetary system formation and early evolution. Aims. We aim to further constrain the orbital and physical properties of β Pictoris b and c using a combination of high contrast imaging, long base-line interferometry, and RV data. We also predict the closest approaches or the transit times of both planets, and we constrain the presence of additional planets in the system. Methods. We obtained six additional epochs of SPHERE data, six additional epochs of GRAVITY data, and five additional epochs of RV data. We combined these various types of data in a single Markov-chain Monte Carlo analysis to constrain the orbital parameters and masses of the two planets simultaneously. The analysis takes into account the gravitational influence of both planets on the star and hence their relative astrometry. Secondly, we used the RV and high contrast imaging data to derive the probabilities of presence of additional planets throughout the disk, and we tested the impact of absolute astrometry. Results. The orbital properties of both planets are constrained with a semi-major axis of 9.8 ± 0.4 au and 2.7 ± 0.02 au for b and c, respectively, and eccentricities of 0.09 ± 0.1 and 0.27 ± 0.07, assuming the HIPPARCOS distance. We note that despite these low fitting error bars, the eccentricity of β Pictoris c might still be over-estimated. If no prior is provided on the mass of β Pictoris b, we obtain a very low value that is inconsistent with what is derived from brightness-mass models. When we set an evolutionary model motivated prior to the mass of β Pictoris b, we find a solution in the 10–11 MJup range. Conversely, β Pictoris c’s mass is well constrained, at 7.8 ± 0.4 MJup, assuming both planets are on coplanar orbits. These values depend on the assumptions on the distance of the β Pictoris system. The absolute astrometry HIPPARCOS-Gaia data are consistent with the solutions presented here at the 2σ level, but these solutions are fully driven by the relative astrometry plus RV data. Finally, we derive unprecedented limits on the presence of additional planets in the disk. We can now exclude the presence of planets that are more massive than about 2.5 MJup closer than 3 au, and more massive than 3.5 MJup between 3 and 7.5 au. Beyond 7.5 au, we exclude the presence of planets that are more massive than 1–2 MJup. Conclusions. Combining relative astrometry and RVs allows one to precisely constrain the orbital parameters of both planets and to give lower limits to potential additional planets throughout the disk. The mass of β Pictoris c is also well constrained, while additional RV data with appropriate observing strategies are required to properly constrain the mass of β Pictoris b.

TW Hya: an old protoplanetary disc revived by its planet

Abstract Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of a ∼40 M⊕ Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way – by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

Super-Earth masses sculpted by pebble isolation around stars of different masses

We developed a pebble-driven core accretion model to study the formation and evolution of planets around stars in the stellar mass range of 0.08 M⊙–1 M⊙. By Monte Carlo sampling of the initial conditions, the growth and migration of a large number of individual protoplanetary embryos were simulated in a population synthesis manner. We tested two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplanetary disks. Two types of disks with different turbulent viscous parametersαtof 10−3and 10−4are also investigated to shed light on the role of outward migration of protoplanets. The forming planets are compared with the observed exoplanets in terms of mass, semimajor axis, metallicity, and water content. We find that gas giant planets are likely to form when the characteristic disk sizes are larger, the disk accretion rates are higher, the disks are more metal rich, and/or their stellar hosts are more massive. Our model shows that first, the characteristic mass of super-Earth is set by the pebble isolation mass. Super-Earth masses increase linearly with the mass of its stellar host, which corresponds to one Earth mass around a late M-dwarf star and 20 Earth masses around a solar-mass star. Second, the low-mass planets, up to 20 M⊕, can form around stars with a wide range of metallicities, while massive gas giant planets are preferred to grow around metal rich stars. Third, super-Earth planets that are mainly composed of silicates, with relatively low water fractions, can form from protoplanetary embryos at the water ice line in weakly turbulent disks where outward migration is suppressed. However, if the embryos are formed over a wide range of radial distances, the super-Earths would end up having a distinctive, bimodal composition in water mass. Altogether, our model succeeds in quantitatively reproducing several important observed properties of exoplanets and correlations with their stellar hosts.

An Eccentric Massive Jupiter Orbiting a Subgiant on a 9.5-day Period Discovered in the Transiting Exoplanet Survey Satellite Full Frame Images

Abstract We report the discovery of TOI-172 b from the Transiting Exoplanet Survey Satellite (TESS) mission, a massive hot Jupiter transiting a slightly evolved G star with a 9.48-day orbital period. This is the first planet to be confirmed from analysis of only the TESS full frame images, because the host star was not chosen as a two-minute cadence target. From a global analysis of the TESS photometry and follow-up observations carried out by the TESS Follow-up Observing Program Working Group, TOI-172 (TIC 29857954) is a slightly evolved star with an effective temperature of T eff = 5645 ± 50 K, a mass of M ⋆ = M ⊙, radius of R ⋆ = R ⊙, a surface gravity of log g ⋆ = , and an age of . Its planetary companion (TOI-172 b) has a radius of R P = R J, a mass of M P = M J, and is on an eccentric orbit ( ). TOI-172 b is one of the few known massive giant planets on a highly eccentric short-period orbit. Future study of the atmosphere of this planet and its system architecture offer opportunities to understand the formation and evolution of similar systems.