Showing 1–50 of 147 results for author: Rice, K

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Gliese 12 b, A Temperate Earth-sized Planet at 12 Parsecs Discovered with TESS and CHEOPS

Authors: Shishir Dholakia, Larissa Palethorpe, Alexander Venner, Annelies Mortier, Thomas G. Wilson, Chelsea X. Huang, Ken Rice, Vincent Van Eylen, Emma Nabbie, Ryan Cloutier, Walter Boschin, David Ciardi, Laetitia Delrez, Georgina Dransfield, Elsa Ducrot, Zahra Essack, Mark E. Everett, Michaël Gillon, Matthew J. Hooton, Michelle Kunimoto, David W. Latham, Mercedes López-Morales, Bin Li, Fan Li, Scott McDermott , et al. (11 additional authors not shown)

Abstract: We report on the discovery of Gliese 12 b, the nearest transiting temperate, Earth-sized planet found to date. Gliese 12 is a bright ($V=12.6$ mag, $K=7.8$ mag) metal-poor M4V star only $12.162\pm0.005$ pc away from the Solar System with one of the lowest stellar activity levels known for an M-dwarf. A planet candidate was detected by TESS based on only 3 transits in sectors 42, 43, and 57, with a… ▽ More We report on the discovery of Gliese 12 b, the nearest transiting temperate, Earth-sized planet found to date. Gliese 12 is a bright ($V=12.6$ mag, $K=7.8$ mag) metal-poor M4V star only $12.162\pm0.005$ pc away from the Solar System with one of the lowest stellar activity levels known for an M-dwarf. A planet candidate was detected by TESS based on only 3 transits in sectors 42, 43, and 57, with an ambiguity in the orbital period due to observational gaps. We performed follow-up transit observations with CHEOPS and ground-based photometry with MINERVA-Australis, SPECULOOS, and Purple Mountain Observatory, as well as further TESS observations in sector 70. We statistically validate Gliese 12 b as a planet with an orbital period of $12.76144\pm0.00006$ days and a radius of $1.0\pm{0.1}$ R$_\oplus$, resulting in an equilibrium temperature of $\sim$315K. Gliese 12 b has excellent future prospects for precise mass measurement, which may inform how planetary internal structure is affected by the stellar compositional environment. Gliese 12 b also represents one of the best targets to study whether Earth-like planets orbiting cool stars can retain their atmospheres, a crucial step to advance our understanding of habitability on Earth and across the Galaxy. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 19 pages, 7 figures, Accepted for publication in MNRAS, Authors Shishir Dholakia and Larissa Palethorpe contributed equally

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

Authors: Alison K. Young, Maggie Celeste, Richard A. Booth, Ken Rice, Adam Koval, Ethan Carter, Dimitris Stamatellos

Abstract: The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate… ▽ More The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths respectively. The second, the "Modified Lombardi" method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the Modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted to MNRAS. 10 pages

Confronting compositional confusion through the characterisation of the sub-Neptune orbiting HD 77946

Authors: L. Palethorpe, A. Anna John, A. Mortier, J. Davoult, T. G. Wilson, K. Rice, A. C. Cameron, Y. Alibert, L. A. Buchhave, L. Malavolta, J. Cadman, M. López-Morales, X. Dumusque, A. M. Silva, S. N. Quinn, V. Van Eylen, S. Vissapragada, L. Affer, D. Charbonneau, R. Cosentino, A. Ghedina, R. D. Haywood, D. W. Latham, F. Lienhard, A. F. Martínez Fiorenzano , et al. (7 additional authors not shown)

Abstract: We report on the detailed characterization of the HD 77946 planetary system. HD 77946 is an F5 ($M_*$ = 1.17 M$_{\odot}$, $R_*$ = 1.31 R$_{\odot}$) star, which hosts a transiting planet recently discovered by NASA's Transiting Exoplanet Survey Satellite (TESS), classified as TOI-1778 b. Using TESS photometry, high-resolution spectroscopic data from HARPS-N, and photometry from CHEOPS, we measure t… ▽ More We report on the detailed characterization of the HD 77946 planetary system. HD 77946 is an F5 ($M_*$ = 1.17 M$_{\odot}$, $R_*$ = 1.31 R$_{\odot}$) star, which hosts a transiting planet recently discovered by NASA's Transiting Exoplanet Survey Satellite (TESS), classified as TOI-1778 b. Using TESS photometry, high-resolution spectroscopic data from HARPS-N, and photometry from CHEOPS, we measure the radius and mass from the transit and RV observations, and find that the planet, HD 77946 b, orbits with period $P_{\rm b}$ = $6.527282_{-0.000020}^{+0.000015}$ d, has a mass of $M_{\rm b} = 8.38\pm{1.32}$M$_\oplus$, and a radius of $R_{\rm b} = 2.705_{-0.081}^{+0.086}$R$_\oplus$. From the combination of mass and radius measurements, and the stellar chemical composition, the planet properties suggest that HD 77946 b is a sub-Neptune with a $\sim$1\% H/He atmosphere. However, a degeneracy still exists between water-world and silicate/iron-hydrogen models, and even though interior structure modelling of this planet favours a sub-Neptune with a H/He layer that makes up a significant fraction of its radius, a water-world composition cannot be ruled out, as with $T_{\rm eq} = 1248^{+40}_{-38}~$K, water may be in a supercritical state. The characterisation of HD 77946 b, adding to the small sample of well-characterised sub-Neptunes, is an important step forwards on our journey to understanding planetary formation and evolution pathways. Furthermore, HD 77946 b has one of the highest transmission spectroscopic metrics for small planets orbiting hot stars, thus transmission spectroscopy of this key planet could prove vital for constraining the compositional confusion that currently surrounds small exoplanets. △ Less

Submitted 1 May, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 19 pages, 13 figures, 5 tables. Published in MNRAS

Journal ref: MNRAS, 529, 3323-3341 (2024)

Populating The Milky Way: Characterising Planet Demographics by Combining Galaxy Formation Simulations and Planet Population Synthesis Models

Authors: C. Boettner, P. Dayal, M. Trebitsch, N. Libeskind, K. Rice, C. Cockell, B. I. Tieleman

Abstract: Stellar populations and their distribution differ widely across the Galaxy, which is likely to affect planet demographics. Our local neighbourhood is dominated by young, metal-rich stars in the galactic thin disc, while the stellar halo and galactic bulge host a large fraction of older, metal-poor stars. We study the impact of these variations on planet populations in different regions of the Gala… ▽ More Stellar populations and their distribution differ widely across the Galaxy, which is likely to affect planet demographics. Our local neighbourhood is dominated by young, metal-rich stars in the galactic thin disc, while the stellar halo and galactic bulge host a large fraction of older, metal-poor stars. We study the impact of these variations on planet populations in different regions of the Galaxy by combining a high-resolution galaxy formation simulation with state-of-the-art planet population synthesis models. We construct a population model to estimate occurrence rates of different planet types, based on the New Generation Planet Population Synthesis by Emsenhuber et al., 2021. We apply this to a simulated Milky Way--Analogue in the HESTIA galaxy formation simulation. We study the planet occurrence rate in the metal-rich regions of the inner Galaxy, i.e. in the galactic bulge and thin disc, and contrast them to the frequencies in the more distant, metal-poor region like the thick disc and stellar halo. We find that the planet demographics in the metal-poor regions of the Milky Way-Analogue, differ strongly from the planet populations in the more distant, metal-poor regions. The occurrence rate of giant planets ($>300 M_\bigoplus$) is 10 to 20 times larger in the thin disc compared to the thick disc, driven by the low amounts of solid material available for planet formation around metal-poor stars. Similarly, low-mass Earth-like planets around Sun-like stars are most abundant in the thick disc, being 1.5 times more frequent than in the thin disc. Moreover, low-mass planets are expected to be abundant throughout the galaxy, from the central regions to the outer halo, due to their formation processes being less dependent on stellar metallicity. The planet populations differ more strongly around Sun-like stars compared to dwarfs with masses 0.3 - 0.5 $M_\odot$. △ Less

Submitted 2 April, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 19 pages, 13 figures, accepted by A&A. Code available at https://github.com/ChrisBoettner/skaro

Journal ref: A&A 686, A167 (2024)

The Role of Drag and Gravity on Dust Concentration in a Gravitationally Unstable Disc

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru, James Wurster, Hossam Aly, Richard Alexander, Ken Rice, Richard A. Booth

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If w… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If we include both, grains that would otherwise be partially decoupled will become well coupled and trace the spirals. For the dust grains most influenced by drag (with Stokes numbers near unity), the dust disc quickly becomes gravitationally unstable and rapidly forms clumps with masses between 0.15 - 6 Earth masses. A large fraction of clumps are below the threshold where runaway gas accretion can occur. However, if dust self-gravity is neglected, the dust is unable to form clumps, despite still becoming trapped in the gas spirals. When large dust grains are unable to feel either gas gravity or drag, the dust is unable to trace the gas spirals. Hence, full physics is needed to properly simulate dust in gravitationally unstable discs. Dust trapping of large grains in spiral arms of discs stable to gas fragmentation could explain planet formation in very young discs by a population of planetesimals formed due to the combined roles of drag and gravity in the earliest stages of a disc's evolution. Furthermore, it highlights that gravitationally unstable discs are not just important for forming gas giants quickly, it can also rapidly form Earth mass bodies. △ Less

Submitted 24 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. Updated references

The magnetically quiet solar surface dominates HARPS-N solar RVs during low activity

Authors: Ben S. Lakeland, Tim Naylor, Raphaëlle Haywood, Nadège Meunier, Federica Rescigno, Shweta Dalal, Annelies Mortier, Samantha J. Thompson, Andrew Collier Cameron, Xavier Dumusque, Mercedes López-Morales, Francesco Pepe, Ken Rice, Alessandro Sozzetti, Stéphane Udry, Eric Ford, Adriano Ghedina, Marcello Lodi

Abstract: Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities obser… ▽ More Using images from the Helioseismic and Magnetic Imager aboard the \textit{Solar Dynamics Observatory} (SDO/HMI), we extract the radial-velocity (RV) signal arising from the suppression of convective blue-shift and from bright faculae and dark sunspots transiting the rotating solar disc. We remove these rotationally modulated magnetic-activity contributions from simultaneous radial velocities observed by the HARPS-N solar feed to produce a radial-velocity time series arising from the magnetically quiet solar surface (the 'inactive-region radial velocities'). We find that the level of variability in the inactive-region radial velocities remains constant over the almost 7 year baseline and shows no correlation with well-known activity indicators. With an RMS of roughly 1 m/s, the inactive-region radial-velocity time series dominates the total RV variability budget during the decline of solar cycle 24. Finally, we compare the variability amplitude and timescale of the inactive-region radial velocities with simulations of supergranulation. We find consistency between the inactive-region radial-velocity and simulated time series, indicating that supergranulation is a significant contribution to the overall solar radial velocity variability, and may be the main source of variability towards solar minimum. This work highlights supergranulation as a key barrier to detecting Earth twins. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 12 pages, 11 figures, accepted to MNRAS

Binary planet formation through tides

Authors: C. Lazzoni, K. W. Rice, A. Zurlo, S. Hinkley, S. Desidera

Abstract: The search for satellites around exoplanets represents one of the greatest challenges in advancing the characterization of planetary systems. Currently, we can only detect massive satellites, which resemble additional planetary companions rather than rocky moons. It is not yet well understood whether such substellar pairs, known as binary planets, are common or how they form. In this study, we inv… ▽ More The search for satellites around exoplanets represents one of the greatest challenges in advancing the characterization of planetary systems. Currently, we can only detect massive satellites, which resemble additional planetary companions rather than rocky moons. It is not yet well understood whether such substellar pairs, known as binary planets, are common or how they form. In this study, we investigated the formation scenario for binary planets resulting from tidal dissipation during close encounters in the gravitational instability scenario. We conducted seven sets of simulations, varying the number of initial planets injected into the system from two to five, as well as the amount of energy lost due to tides. Our results demonstrate that this formation mechanism is quite efficient in producing binary planets, with an average occurrence rate for the simulated systems of 14.3%. Additionally, we present the distribution of relevant physical parameters (semi-major axis, eccentricity, mass ratios, and formation time) for planet-planet pairs. We also provide comprehensive statistics for single planets and planet-planet pairs. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 11 pages, 12 figures

HD152843 b & c: the masses and orbital periods of a sub-Neptune and a super-puff Neptune

Authors: B. A. Nicholson, S. Aigrain, N. L. Eisner, M. Cretignier, O. Barragán, L. Kaye, J. Taylor, J. Owen, A. Mortier, L. Affer, W. Boschin, A. Collier Cameron, M. Damasso, L. Di Fabrizio, V. DiTomasso, X. Dumusque, A. Gehdina, A. Harutyunyan, D. W. Latham, M. Lopez-Morales, V. Lorenzi, A. F. Martínez Fiorenzano, E. Molinari, M. Pedani, M. Pinamonti , et al. (2 additional authors not shown)

Abstract: We present the characterisation of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 M$_\oplus$ but differing radii of $3.05 \pm 0.11$ R$_\oplus$ and $5.94^{+0.18}_{-0.16}$ R$_\oplus$ , re… ▽ More We present the characterisation of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 M$_\oplus$ but differing radii of $3.05 \pm 0.11$ R$_\oplus$ and $5.94^{+0.18}_{-0.16}$ R$_\oplus$ , respectively, and orbital periods of $11.62071^{+9.6e-05}_{-0.000106}$ days and $19.502104^{+7.4e-05}_{-8.5e-05}$ days. This indicates that HD 152843 c is in the lowest fifth percentile in density of the known exoplanet population, and has the longest orbital period among these low density planets. Further, HD 152843 c's radius places it in the Saturn valley, the observed lack of planets larger than Neptune, but smaller than Saturn. The orbital periods of these planets indicate they are near a 5:3 mean motion resonance, indicating the possibility of transit timing variations, and hints at the possibility of interaction with a third planet at some point in the evolution of this system. Further, the brightness of the host star and the low density of HD 152843 c make it a key target for atmospheric characterisation. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Submitted for review to the Monthly Notices of the Royal Astronomical Society. 13 pages, 16 figures

A hot mini-Neptune and a temperate, highly eccentric sub-Saturn around the bright K-dwarf TOI-2134

Authors: F. Rescigno, G. Hébrard, A. Vanderburg, A. W. Mann, A. Mortier, S. Morrell, L. A. Buchhave, K. A. Collins, C. R. Mann, C. Hellier, R. D. Haywood, R. West, M. Stalport, N. Heidari, D. Anderson, C. X. Huang, M. López-Morales, P. Cortés-Zuleta, H. M. Lewis, X. Dumusque, I. Boisse, P. Rowden, A. Collier Cameron, M. Deleuil, M. Vezie , et al. (42 additional authors not shown)

Abstract: We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.… ▽ More We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.27$\pm$0.42 R$_{e}$ for the outer one. We acquired 111 radial-velocity spectra with HARPS-N and 108 radial-velocity spectra with SOPHIE. After careful periodogram analysis, we derive masses for both planets via Gaussian Process regression: 9.13$^{+0.78}_{-0.76}$ M$_{e}$ for TOI-2134b and 41.86$^{+7.69}_{-7.83}$ M$_{e}$ for TOI-2134c. We analysed the photometric and radial-velocity data first separately, then jointly. The inner planet is a mini-Neptune with density consistent with either a water-world or a rocky core planet with a low-mass H/He envelope. The outer planet has a bulk density similar to Saturn's. The outer planet is derived to have a significant eccentricity of 0.67$^{+0.05}_{-0.06}$ from a combination of photometry and RVs. We compute the irradiation of TOI-2134c as 1.45$\pm$0.10 times the bolometric flux received by Earth, positioning it for part of its orbit in the habitable sone of its system. We recommend further RV observations to fully constrain the orbit of TOI-2134c. With an expected Rossiter-McLaughlin (RM) effect amplitude of 7.2$\pm$1.3 m/s, we recommend TOI-2134c for follow-up RM analysis to study the spin-orbit architecture of the system. We calculate the Transmission Spectroscopy Metric, and both planets are suitable for bright-mode NIRCam atmospheric characterisation. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 24 pages, 7 tables, 14 figures

Masses, Revised Radii, and a Third Planet Candidate in the "Inverted" Planetary System Around TOI-1266

Authors: Ryan Cloutier, Michael Greklek-McKeon, Serena Wurmser, Collin Cherubim, Erik Gillis, Andrew Vanderburg, Sam Hadden, Charles Cadieux, Étienne Artigau, Shreyas Vissapragada, Annelies Mortier, Mercedes López-Morales, David W. Latham, Heather Knutson, Raphaëlle D. Haywood, Enric Pallé, René Doyon, Neil Cook, Gloria Andreuzzi, Massimo Cecconi, Rosario Cosentino, Adriano Ghedina, Avet Harutyunyan, Matteo Pinamonti, Manu Stalport , et al. (18 additional authors not shown)

Abstract: Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a larg… ▽ More Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a large sub-Neptune ($P_b=10.9$ days, $R_{p,b}=2.62\pm 0.11\, \mathrm{R}_{\oplus}$) orbiting interior to that of the system's smaller planet ($P_c=18.8$ days, $R_{p,c}=2.13\pm 0.12\, \mathrm{R}_{\oplus}$). Here we present revised planetary radii based on new TESS and diffuser-assisted ground-based transit observations, and characterize both planetary masses using a set of 145 radial velocity measurements from HARPS-N ($M_{p,b}=4.23\pm 0.69\, \mathrm{M}_{\oplus}, M_{p,c}=2.88\pm 0.80\, \mathrm{M}_{\oplus}$). Our analysis also reveals a third planet candidate ($P_d=32.3$ days, $M_{p,d}\sin{i} = 4.59^{+0.96}_{-0.94}\, \mathrm{M}_{\oplus}$), which if real, would form a chain of near 5:3 period ratios, although the system is likely not in a mean motion resonance. Our results indicate that TOI-1266 b and c are among the lowest density sub-Neptunes around M dwarfs and likely exhibit distinct bulk compositions of a gas-enveloped terrestrial ($X_{\mathrm{env},b}=5.5\pm 0.7$%) and a water-rich world (WMF$_c=59\pm 14$%), which is supported by hydrodynamic escape models. If distinct bulk compositions are confirmed through atmospheric characterization, the system's unique architecture would represent an interesting test case of inside-out sub-Neptune formation at pebble traps. △ Less

Submitted 6 November, 2023; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 21 pages. Our spectroscopic time series are included in the arXiv source files as table6.csv

Eruptivity Criteria for Solar Coronal Flux Ropes in Magnetohydrodynamic and Magnetofrictional Models

Authors: Oliver E. K. Rice, Anthony R. Yeates

Abstract: We investigate which scalar quantity or quantities can best predict the loss of equilibrium and subsequent eruption of magnetic flux ropes in the solar corona. Our models are initialized with a potential magnetic arcade, which is then evolved by means of two effects on the lower boundary: firstly a gradual shearing of the arcade, modelling differential rotation on the solar surface, and secondly s… ▽ More We investigate which scalar quantity or quantities can best predict the loss of equilibrium and subsequent eruption of magnetic flux ropes in the solar corona. Our models are initialized with a potential magnetic arcade, which is then evolved by means of two effects on the lower boundary: firstly a gradual shearing of the arcade, modelling differential rotation on the solar surface, and secondly supergranular diffusion. These result in flux cancellation at the polarity inversion line and the formation of a twisted flux rope. We use three model setups: full magnetohydrodynamics (MHD) in cartesian coordinates, and the magnetofrictional model in both cartesian and polar coordinates. The flux ropes are translationally-invariant, allowing for very fast computational times and thus a comprehensive parameter study, comprising hundreds of simulations and thousands of eruptions. Similar flux rope behavior is observed using either magnetofriction or MHD, and there are several scalar criteria that could be used as proxies for eruptivity. The most consistent predictor of eruptions in either model is the squared current in the axial direction of the rope, normalised by the relative helicity, although a variation on the previously proposed 'eruptivity index' is also found to perform well in both the magnetofrictional and MHD simulations. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 20 pages, 11 figures, 1 table. Accepted to the Astrophysical Journal on 13th August 2023

A review of planetary systems around HD 99492, HD 147379 and HD 190007 with HARPS-N

Authors: M. Stalport, M. Cretignier, S. Udry, A. Anna John, T. G. Wilson, J. -B. Delisle, A. S. Bonomo, L. A. Buchhave, D. Charbonneau, S. Dalal, M. Damasso, L. Di Fabrizio, X. Dumusque, A. Fiorenzano, A. Harutyunyan, R. D. Haywood, D. W. Latham, M. López-Morales, V. Lorenzi, C. Lovis, L. Malavolta, E. Molinari, A. Mortier, M. Pedani, F. Pepe , et al. (4 additional authors not shown)

Abstract: The Rocky Planet Search (RPS) program is dedicated to a blind radial velocity (RV) search of planets around bright stars in the Northern hemisphere, using the high-resolution echelle spectrograph HARPS-N installed on the Telescopio Nazionale Galileo (TNG). The goal of this work is to revise and update the properties of three planetary systems by analysing the HARPS-N data with state-of-the-art s… ▽ More The Rocky Planet Search (RPS) program is dedicated to a blind radial velocity (RV) search of planets around bright stars in the Northern hemisphere, using the high-resolution echelle spectrograph HARPS-N installed on the Telescopio Nazionale Galileo (TNG). The goal of this work is to revise and update the properties of three planetary systems by analysing the HARPS-N data with state-of-the-art stellar activity mitigation tools. The stars considered are HD 99492 (83Leo B), HD 147379 (Gl617 A) and HD 190007. We employ a systematic process of data modelling, that we selected from the comparison of different approaches. We use YARARA to remove instrumental systematics from the RV, and then use SPLEAF to further mitigate the stellar noise with a multidimensional correlated noise model. We also search for transit features in the Transiting Exoplanets Survey Satellite (TESS) data of these stars. We report on the discovery of a new planet around HD 99492, namely HD 99492 c, with an orbital period of 95.2 days and a minimum mass of msin i = 17.9 M_Earth, and refine the parameters of HD 99492 b. We also update and refine the Keplerian solutions for the planets around HD 147379 and HD 190007, but do not detect additional planetary signals. We discard the transiting geometry for the planets, but stress that TESS did not exhaustively cover all the orbital phases. The addition of the HARPS-N data, and the use of advanced data analysis tools, has allowed us to present a more precise view of these three planetary systems. It demonstrates once again the importance of long observational efforts such as the RPS program. Added to the RV exoplanet sample, these planets populate two apparently distinct populations revealed by a bimodality in the planets minimum mass distribution. The separation is located between 30 and 50 M_Earth. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 27 pages, 26 figures (13 in Appendix); Accepted for publication in A&A

Journal ref: A&A 678, A90 (2023)

On the conditions for warping and breaking protoplanetary discs

Authors: Alison K. Young, Struan Stevenson, C. J. Nixon, Ken Rice

Abstract: Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc mo… ▽ More Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc morphology in the wave-like regime of warp propagation that is appropriate to protoplanetary discs. The simulations confirm that disc tearing is possible in protoplanetary discs as long as the mass ratio, $μ$, and disc-binary inclination angle, $i$, are not too small. For the simulations presented here this corresponds to $μ> 0.1$ and $i \gtrsim 40^\circ$. For highly eccentric binaries, tearing can occur for discs with smaller misalignment. Existing theoretical predictions provide an estimate of the radial extent of the disc in which we can expect breaking to occur. However, there does not seem to be a simple relationship between the disc properties and the radius within the circumbinary disc at which the breaks appear, and furthermore the radius at which the disc breaks can change as a function of time in each case. We discuss the implications of our results for interpreting observations and suggest some considerations for modelling misaligned discs in the future. △ Less

Submitted 11 August, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: Accepted to MNRAS. 16 pages

TOI-4010: A System of Three Large Short-Period Planets With a Massive Long-Period Companion

Authors: Michelle Kunimoto, Andrew Vanderburg, Chelsea X. Huang, M. Ryleigh Davis, Laura Affer, Andrew Collier Cameron, David Charbonneau, Rosario Cosentino, Mario Damasso, Xavier Dumusque, A. F. Martnez Fiorenzano, Adriano Ghedina, R. D. Haywood, Florian Lienhard, Mercedes López-Morales, Michel Mayor, Francesco Pepe, Matteo Pinamonti, Ennio Poretti, Jesús Maldonado, Ken Rice, Alessandro Sozzetti, Thomas G. Wilson, Stéphane Udry, Jay Baptista , et al. (31 additional authors not shown)

Abstract: We report the confirmation of three exoplanets transiting TOI-4010 (TIC-352682207), a metal-rich K dwarf observed by TESS in Sectors 24, 25, 52, and 58. We confirm these planets with HARPS-N radial velocity observations and measure their masses with 8 - 12% precision. TOI-4010 b is a sub-Neptune ($P = 1.3$ days, $R_{p} = 3.02_{-0.08}^{+0.08}~R_{\oplus}$, $M_{p} = 11.00_{-1.27}^{+1.29}~M_{\oplus}$)… ▽ More We report the confirmation of three exoplanets transiting TOI-4010 (TIC-352682207), a metal-rich K dwarf observed by TESS in Sectors 24, 25, 52, and 58. We confirm these planets with HARPS-N radial velocity observations and measure their masses with 8 - 12% precision. TOI-4010 b is a sub-Neptune ($P = 1.3$ days, $R_{p} = 3.02_{-0.08}^{+0.08}~R_{\oplus}$, $M_{p} = 11.00_{-1.27}^{+1.29}~M_{\oplus}$) in the hot Neptune desert, and is one of the few such planets with known companions. Meanwhile, TOI-4010 c ($P = 5.4$ days, $R_{p} = 5.93_{-0.12}^{+0.11}~R_{\oplus}$, $M_{p} = 20.31_{-2.11}^{+2.13}~M_{\oplus}$) and TOI-4010 d ($P = 14.7$ days, $R_{p} = 6.18_{-0.14}^{+0.15}~R_{\oplus}$, $M_{p} = 38.15_{-3.22}^{+3.27}~M_{\oplus}$) are similarly-sized sub-Saturns on short-period orbits. Radial velocity observations also reveal a super-Jupiter-mass companion called TOI-4010 e in a long-period, eccentric orbit ($P \sim 762$ days and $e \sim 0.26$ based on available observations). TOI-4010 is one of the few systems with multiple short-period sub-Saturns to be discovered so far. △ Less

Submitted 19 June, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 26 pages, 16 figures, published in AJ; (v3) added missing citation

Journal ref: AJ, 166, 7 (2023)

Cold Jupiters and improved masses in 38 Kepler and K2 small planet systems from 3661 HARPS-N radial velocities. No excess of cold Jupiters in small planet systems

Authors: A. S. Bonomo, X. Dumusque, A. Massa, A. Mortier, R. Bongiolatti, L. Malavolta, A. Sozzetti, L. A. Buchhave, M. Damasso, R. D. Haywood, A. Morbidelli, D. W. Latham, E. Molinari, F. Pepe, E. Poretti, S. Udry, L. Affer, W. Boschin, D. Charbonneau, R. Cosentino, M. Cretignier, A. Ghedina, E. Lega, M. López-Morales, M. Margini , et al. (9 additional authors not shown)

Abstract: The exoplanet population characterized by relatively short orbital periods ($P<100$ d) around solar-type stars is dominated by super-Earths and sub-Neptunes. However, these planets are missing in our Solar System and the reason behind this absence is still unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inw… ▽ More The exoplanet population characterized by relatively short orbital periods ($P<100$ d) around solar-type stars is dominated by super-Earths and sub-Neptunes. However, these planets are missing in our Solar System and the reason behind this absence is still unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inward migration of sub-Neptunes from beyond the water iceline; alternatively, Jupiter may have reduced considerably the inward flux of material (pebbles) required to form super-Earths inside that iceline. Both scenarios predict an anti-correlation between the presence of small planets (SPs) and that of cold Jupiters (CJs) in exoplanetary systems. To test that prediction, we homogeneously analyzed the radial-velocity (RV) measurements of 38 Kepler and K2 transiting SP systems gathered over almost 10 years with the HARPS-N spectrograph, as well as publicly available RVs collected with other facilities. We detected five CJs in three systems, two in Kepler-68, two in Kepler-454, and a very eccentric one in K2-312. We derived an occurrence rate of $9.3^{+7.7}_{-2.9}\%$ for CJs with $0.3-13~M_{Jup}$ and 1-10 AU, which is lower but still compatible at $1.3σ$ with that measured from RV surveys for solar-type stars, regardless of the presence or absence of SPs. The sample is not large enough to draw a firm conclusion about the predicted anti-correlation between SPs and CJs; nevertheless, we found no evidence of previous claims of an excess of CJs in SP systems. As an important by-product of our analyses, we homogeneously determined the masses of 64 Kepler and K2 small planets, reaching a precision better than 5, 7.5 and 10$σ$ for 25, 13 and 8 planets, respectively. Finally, we release the 3661 HARPS-N radial velocities used in this work to the scientific community. [Abridged] △ Less

Submitted 6 September, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 21 pages, 10 figures, 10 tables, published in Astronomy and Astrophysics. The updated version of the article takes into account the A&A language editing and guidelines. Tables 1, A.1 and full Table 2 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/677/A33

Journal ref: A&A 677, A33 (2023)

Retrieval study of cool, directly imaged exoplanet 51 Eri b

Authors: Niall Whiteford, Alistair Glasse, Katy L. Chubb, Daniel Kitzmann, Shrishmoy Ray, Mark W. Phillips, Beth A. Biller, Paul I. Palmer, Ken Rice, Ingo P. Waldmann, Quentin Changeat, Nour Skaf, Jason Wang, Billy Edwards, Ahmed Al-Refaie

Abstract: Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophot… ▽ More Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$σ$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures. △ Less

Submitted 15 February, 2023; originally announced February 2023.

TOI-1695 b: A Water World Orbiting an Early M Dwarf in the Planet Radius Valley

Authors: Collin Cherubim, Ryan Cloutier, David Charbonneau, Bill Wohler, Chris Stockdale, Keivan G. Stassun, Richard P. Schwarz, Boris Safonov, Annelies Mortier, David W. Latham, Keith Horne, Raphaëlle D. Haywood, Erica Gonzales, Maria V. Goliguzova, Karen A. Collins, David R. Ciardi, Allyson Bieryla, Alexander A. Belinski, Christopher A. Watson, Rolands Vanderspek, Stéphane Udry, Alessandro Sozzetti, Damien Ségransan, Dimitar Sasselov, George R. Ricker , et al. (16 additional authors not shown)

Abstract: Characterizing the bulk compositions of transiting exoplanets within the M dwarf radius valley offers a unique means to establish whether the radius valley emerges from an atmospheric mass loss process or is imprinted by planet formation itself. We present the confirmation of such a planet orbiting an early M dwarf (… ▽ More Characterizing the bulk compositions of transiting exoplanets within the M dwarf radius valley offers a unique means to establish whether the radius valley emerges from an atmospheric mass loss process or is imprinted by planet formation itself. We present the confirmation of such a planet orbiting an early M dwarf ($T_{\rm mag} = 11.0294 \pm 0.0074, M_s = 0.513 \pm 0.012\ M_\odot, R_s = 0.515 \pm 0.015\ R_\odot, T_{\rm eff} =3690\pm 50 K$): TOI-1695 b ($P = 3.13$ days, $R_p = 1.90^{+0.16}_{-0.14}\ R_\oplus$). TOI-1695 b's radius and orbital period situate the planet between model predictions from thermally-driven mass loss versus gas depleted formation, offering an important test case for radius valley emergence models around early M dwarfs. We confirm the planetary nature of TOI-1695 b based on five sectors of TESS data and a suite of follow-up observations including 49 precise radial velocity measurements taken with the HARPS-N spectrograph. We measure a planetary mass of $6.36 \pm 1.00\ M_\oplus$, which reveals that TOI-1695 b is inconsistent with a purely terrestrial composition of iron and magnesium silicate, and instead is likely a water-rich planet. Our finding that TOI-1695 b is not terrestrial is inconsistent with the planetary system being sculpted by thermally driven mass loss. We present a statistical analysis of seven well-characterized planets within the M dwarf radius valley demonstrating that a thermally-driven mass loss scenario is unlikely to explain this population. △ Less

Submitted 13 February, 2023; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 22 pages, 11 figures. Accepted in AJ

Kepler-102: Masses and Compositions for a Super-Earth and Sub-Neptune Orbiting an Active Star

Authors: Casey Brinkman, James Cadman, Lauren Weiss, Eric Gaidos, Ken Rice, Daniel Huber, Zachary R. Claytor, Aldo S. Bonomo, Lars A. Buchhave, Andrew Collier Cameron, Rosario Cosentino, Xavier Dumusque, Aldo F Martinez Fiorenzano, Adriano Ghedina, Avet Harutyunyan, Andrew Howard, Howard Isaacson, David W. Latham, Mercedes Lopez-Morales, Luca Malavolta, Giuseppina Micela, Emilio Molinari, Francesco Pepe, David F Philips, Ennio Poretti , et al. (2 additional authors not shown)

Abstract: Radial velocity (RV) measurements of transiting multiplanet systems allow us to understand the densities and compositions of planets unlike those in the Solar System. Kepler-102, which consists of 5 tightly packed transiting planets, is a particularly interesting system since it includes a super-Earth (Kepler-102d) and a sub-Neptune-sized planet (Kepler-102e) for which masses can be measured using… ▽ More Radial velocity (RV) measurements of transiting multiplanet systems allow us to understand the densities and compositions of planets unlike those in the Solar System. Kepler-102, which consists of 5 tightly packed transiting planets, is a particularly interesting system since it includes a super-Earth (Kepler-102d) and a sub-Neptune-sized planet (Kepler-102e) for which masses can be measured using radial velocities. Previous work found a high density for Kepler-102d, suggesting a composition similar to that of Mercury, while Kepler-102e was found to have a density typical of sub-Neptune size planets; however, Kepler-102 is an active star, which can interfere with RV mass measurements. To better measure the mass of these two planets, we obtained 111 new RVs using Keck/HIRES and TNG/HARPS-N and modeled Kepler-102's activity using quasi-periodic Gaussian Process Regression. For Kepler-102d, we report a mass upper limit of M$_{d} < $5.3 M$_{\oplus}$ [95\% confidence], a best-fit mass of M$_{d}$=2.5 $\pm$ 1.4 M$_{\oplus}$, and a density of $ρ_{d}$=5.6 $\pm$ 3.2 g/cm$^{3}$ which is consistent with a rocky composition similar in density to the Earth. For Kepler-102e we report a mass of M$_{e}$=4.7 $\pm$ 1.7 M$_{\oplus}$ and a density of $ρ_{e}$=1.8 $\pm$ 0.7 g/cm$^{3}$. These measurements suggest that Kepler-102e has a rocky core with a thick gaseous envelope comprising 2-4% of the planet mass and 16-50% of its radius. Our study is yet another demonstration that accounting for stellar activity in stars with clear rotation signals can yield more accurate planet masses, enabling a more realistic interpretation of planet interiors. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Accepted to AJ 11/08/2022

Independent validation of the temperate Super-Earth HD79211 b using HARPS-N

Authors: Victoria DiTomasso, Chantanelle Nava, Mercedes López-Morales, Allyson Bieryla, Ryan Cloutier, Luca Malavolta, Annelies Mortier, Lars A. Buchhave, Keivan G. Stassun, Alessandro Sozzetti, Aldo Stefano Bonomo, David Charbonneau, Andrew Collier Cameron, Rosario Cosentino, Mario Damasso, Xavier Dumusque, A. F. Martínez Fiorenzano, Adriano Ghedina, Avet Harutyunyan, R. D. Haywood, David Latham, Emilio Molinari, Francesco A. Pepe, Matteo Pinamonti, Ennio Poretti , et al. (6 additional authors not shown)

Abstract: We present high-precision radial velocities (RVs) from the HARPS-N spectrograph for HD79210 and HD79211, two M0V members of a gravitationally-bound binary system. We detect a planet candidate with a period of $24.421^{+0.016}_{-0.017}$ days around HD79211 in these HARPS-N RVs, validating the planet candidate originally identified in CARMENES RV data alone. Using HARPS-N, CARMENES and HIRES RVs spa… ▽ More We present high-precision radial velocities (RVs) from the HARPS-N spectrograph for HD79210 and HD79211, two M0V members of a gravitationally-bound binary system. We detect a planet candidate with a period of $24.421^{+0.016}_{-0.017}$ days around HD79211 in these HARPS-N RVs, validating the planet candidate originally identified in CARMENES RV data alone. Using HARPS-N, CARMENES and HIRES RVs spanning a total of 25 years, we further refine the planet candidate parameters to $P=24.422\pm0.014$ days, $K=3.19\pm0.27$ m/s, $M$ sin $i = 10.6 \pm 1.2 M_\oplus$, and $a = 0.142 \pm0.005$ au. We do not find any additional planet candidate signals in the data of HD79211 nor do we find any planet candidate signals in HD79210. This system adds to the number of exoplanets detected in binaries with M dwarf members, and serves as a case study for planet formation in stellar binaries. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: Accepted for publication in AJ, 29 pages, 17 figures

Eruptivity Criteria for Two-dimensional Magnetic Flux Ropes in the Solar Corona

Authors: Oliver E. K. Rice, Anthony R. Yeates

Abstract: We apply the magneto-frictional approach to investigate which quantity or quantities can best predict the loss of equilibrium of a translationally-invariant magnetic flux rope. The flux rope is produced self-consistently by flux cancellation combined with gradual footpoint shearing of a coronal arcade which is open at the outer boundary. This models the magnetic field in decaying active regions on… ▽ More We apply the magneto-frictional approach to investigate which quantity or quantities can best predict the loss of equilibrium of a translationally-invariant magnetic flux rope. The flux rope is produced self-consistently by flux cancellation combined with gradual footpoint shearing of a coronal arcade which is open at the outer boundary. This models the magnetic field in decaying active regions on the Sun. Such a model permits two types of eruption: episodic small events caused by shearing and relaxation of the overlying arcade, and major eruptions of the main low-lying coronal flux rope. Through a parameter study, we find that the major eruptions are best predicted not by individual quantities but by thresholds in the ratios of squared rope current to either magnetic energy or relative magnetic helicity. We show how to appropriately define the latter quantity for translationally-invariant magnetic fields, along with a related eruptivity index that has recently been introduced for three-dimensional magnetic fields. In contrast to previous configurations studied, we find that the eruptivity index has only a weak predictive skill, and in fact is lower prior to eruption, rather than higher. This is because the overlying background magnetic field has the same direction as the arcade itself. Thus we propose that there are a whole class of solar eruptions that cannot be predicted by a high eruptivity index. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: 25 pages, 11 figures. To be published in 'Frontiers in Astronomy and Space Sciences', accepted for publication 8th March 2022

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N and TESS

Authors: G. Lacedelli, T. G. Wilson, L. Malavolta, M. J. Hooton, A. Collier Cameron, Y. Alibert, A. Mortier, A. Bonfanti, R. D. Haywood, S. Hoyer, G. Piotto, A. Bekkelien, A. M. Vanderburg, W. Benz, X. Dumusque, A. Deline, M. López-Morales, L. Borsato, K. Rice, L. Fossati, D. W. Latham, A. Brandeker, E. Poretti, S. G. Sousa, A. Sozzetti , et al. (93 additional authors not shown)

Abstract: We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$,… ▽ More We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($ρ_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 22 pages, 19 figures. Accepted for publication in MNRAS

Binary companions triggering fragmentation in self-gravitating discs

Authors: James Cadman, Cassandra Hall, Clémence Fontanive, Ken Rice

Abstract: Observations of systems hosting close in ($<1$ AU) giant planets and brown dwarfs ($M\gtrsim7$ M$_{\rm Jup}$) find an excess of binary star companions, indicating that stellar multiplicity may play an important role in their formation. There is now increasing evidence that some of these objects may have formed via fragmentation in gravitationally unstable discs. We present a suite of 3D smoothed p… ▽ More Observations of systems hosting close in ($<1$ AU) giant planets and brown dwarfs ($M\gtrsim7$ M$_{\rm Jup}$) find an excess of binary star companions, indicating that stellar multiplicity may play an important role in their formation. There is now increasing evidence that some of these objects may have formed via fragmentation in gravitationally unstable discs. We present a suite of 3D smoothed particle hydrodynamics (SPH) simulations of binary star systems with circumprimary self-gravitating discs, which include a realistic approximation to radiation transport, and extensively explore the companion's orbital parameter space for configurations which may trigger fragmentation. We identify a "sweet spot" where intermediate separation binary companions ($100$ AU $\lesssim a\lesssim400$ AU) can cause a marginally stable disc to fragment. The exact range of ideal binary separations is a function of the companion's eccentricity, inclination and mass. Heating is balanced by efficient cooling, and fragmentation occurs inside a spiral mode driven by the companion. Short separation, disc penetrating binary encounters ($a\lesssim100$ AU) are prohibitive to fragmentation, as mass stripping and disc heating quench any instability. This is also true of binary companions with high orbital eccentricities ($e\gtrsim0.75$). Wide separation companions ($a\gtrsim500$ AU) have little effect on the disc properties for the setup parameters considered here. The sweet spot found is consistent with the range of binary separations which display an excess of close in giant planets and brown dwarfs. Hence we suggest that fragmentation triggered by a binary companion may contribute to the formation of these substellar objects. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Comments: 15 pages, 8 figures, accepted for publication in MNRAS

K2-79b and K2-222b: Mass measurements of two small exoplanets with periods beyond 10 days that overlap with periodic magnetic activity signals

Authors: Chantanelle Nava, Mercedes López-Morales, Annelies Mortier, Li Zeng, Helen A. C. Giles, Allyson Bieryla, Andrew Vanderburg, Lars A. Buchhave, Ennio Poretti, Steven H. Saar, Xavier Dumusque, David W. Latham, David Charbonneau, Mario Damasso, Aldo S. Bonomo, Christophe Lovis, Andrew Collier Cameron, Jason D. Eastman, Alessandro Sozzetti, Rosario Cosentino, Marco Pedani, Francesco Pepe, Emilio Molinari, Dimitar Sasselov, Michel Mayor , et al. (6 additional authors not shown)

Abstract: We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars. Their respective 10.99d and 15.39d orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity… ▽ More We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars. Their respective 10.99d and 15.39d orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity and activity indicator observations are available over multiple observing seasons and the orbital period of the exoplanet is known. We perform correlation and periodogram analyses on sub-sets composed of each target's two observing seasons, in addition to the full data sets. For both targets, these analyses reveal an optimal season with little to no interference at the orbital period of the known exoplanet. We make a confident mass detection of each exoplanet by confirming agreement between fits to the full radial velocity set and the optimal season. For K2-79b, we measure a mass of 11.8 $\pm$ 3.6 $M_{Earth}$ and a radius of 4.09 $\pm$ 0.17 $R_{Earth}$. For K2-222b, we measure a mass of 8.0 $\pm$ 1.8 $M_{Earth}$ and a radius of 2.35 $\pm$ 0.08 $R_{Earth}$. According to model predictions, K2-79b is a highly irradiated Uranus-analog and K2-222b hosts significant amounts of water ice. We also present an RV solution for a candidate second companion orbiting K2-222 at 147.5d. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: 33 pages, 17 figures, 13 tables

Global Coronal Equilibria with Solar Wind Outflow

Authors: Oliver E. K. Rice, Anthony R. Yeates

Abstract: Given a known radial magnetic field distribution on the Sun's photospheric surface, there exist well-established methods for computing a potential magnetic field in the corona above. Such potential fields are routinely used as input to solar wind models, and to initialize magneto-frictional or full magnetohydrodynamic simulations of the coronal and heliospheric magnetic fields. We describe an impr… ▽ More Given a known radial magnetic field distribution on the Sun's photospheric surface, there exist well-established methods for computing a potential magnetic field in the corona above. Such potential fields are routinely used as input to solar wind models, and to initialize magneto-frictional or full magnetohydrodynamic simulations of the coronal and heliospheric magnetic fields. We describe an improved magnetic field model which calculates a magneto-frictional equilibrium with an imposed solar wind profile (which can be Parker's solar wind solution, or any reasonable equivalent). These `outflow fields' appear to approximate the real coronal magnetic field more closely than a potential field, take a similar time to compute, and avoid the need to impose an artificial source surface. Thus they provide a practical alternative to the potential field model for initializing time-evolving simulations or modeling the heliospheric magnetic field. We give an open-source Python implementation in spherical coordinates and apply the model to data from Solar Cycle 24. The outflow tends to increase the open magnetic flux compared to the potential field model, reducing the well known discrepancy with in situ observations. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 16 Pages, 6 Figures. To be published in the Astrophysical Journal (accepted 1/10/21)

A HARPS-N mass for the elusive Kepler-37d: a case study in disentangling stellar activity and planetary signals

Authors: V. M. Rajpaul, L. A. Buchhave, G. Lacedelli, K. Rice, A. Mortier, L. Malavolta, S. Aigrain, L. Borsato, A. W. Mayo, D. Charbonneau, M. Damasso, X. Dumusque, A. Ghedina, D. W. Latham, M. López-Morales, A. Magazzù, G. Micela, E. Molinari, F. Pepe, G. Piotto, E. Poretti, S. Rowther, A. Sozzetti, S. Udry, C. A. Watson

Abstract: To date, only 18 exoplanets with radial velocity (RV) semi-amplitudes $<2$ m/s have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 - know… ▽ More To date, only 18 exoplanets with radial velocity (RV) semi-amplitudes $<2$ m/s have had their masses directly constrained. The biggest obstacle to RV detection of such exoplanets is variability intrinsic to stars themselves, e.g. nuisance signals arising from surface magnetic activity such as rotating spots and plages, which can drown out or even mimic planetary RV signals. We use Kepler-37 - known to host three transiting planets, one of which, Kepler-37d, should be on the cusp of RV detectability with modern spectrographs - as a case study in disentangling planetary and stellar activity signals. We show how two different statistical techniques - one seeking to identify activity signals in stellar spectra, and another to model activity signals in extracted RVs and activity indicators - can enable detection of the hitherto elusive Kepler-37d. Moreover, we show that these two approaches can be complementary, and in combination, facilitate a definitive detection and precise characterisation of Kepler-37d. Its RV semi-amplitude of $1.22\pm0.31$ m/s (mass $5.4\pm1.4$ $M_\oplus$) is formally consistent with TOI-178b's $1.05^{+0.25}_{-0.30}$ m/s, the latter being the smallest detected RV signal of any transiting planet to date, though dynamical simulations suggest Kepler-37d's mass may be on the lower end of our $1σ$ credible interval. Its consequent density is consistent with either a water-world or that of a gaseous envelope ($\sim0.4\%$ by mass) surrounding a rocky core. Based on RV modelling and a re-analysis of Kepler-37 TTVs, we also argue that the putative (non-transiting) planet Kepler-37e should probably be stripped of its 'confirmed' status. △ Less

Submitted 30 July, 2021; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS. 22 pages, 10 figures, 13 tables

AB Aurigae: Possible evidence of planet formation through the gravitational instability

Authors: James Cadman, Ken Rice, Cassandra Hall

Abstract: Recent observations of the protoplanetary disc surrounding AB Aurigae have revealed the possible presence of two giant planets in the process of forming. The young measured age of $1-4$Myr for this system allows us to place strict time constraints on the formation histories of the observed planets. Hence we may be able to make a crucial distinction between formation through core accretion (CA) or… ▽ More Recent observations of the protoplanetary disc surrounding AB Aurigae have revealed the possible presence of two giant planets in the process of forming. The young measured age of $1-4$Myr for this system allows us to place strict time constraints on the formation histories of the observed planets. Hence we may be able to make a crucial distinction between formation through core accretion (CA) or the gravitational instability (GI), as CA formation timescales are typically Myrs whilst formation through GI will occur within the first $\approx10^4-10^5$yrs of disc evolution. We focus our analysis on the $4-13$M$_{\rm Jup}$ planet observed at $R\approx30$AU. We find CA formation timescales for such a massive planet typically exceed the system's age. The planet's high mass and wide orbit may instead be indicative of formation through GI. We use smoothed particle hydrodynamic simulations to determine the system's critical disc mass for fragmentation, finding $M_{\rm d,crit}=0.3$M$_{\odot}$. Viscous evolution models of the disc's mass history indicate that it was likely massive enough to exceed $M_{\rm d,crit}$ in the recent past, thus it is possible that a young AB Aurigae disc may have fragmented to form multiple giant gaseous protoplanets. Calculations of the Jeans mass in an AB Aurigae-like disc find that fragments may initially form with masses $1.6-13.3$M$_{\rm Jup}$, consistent with the planets which have been observed. We therefore propose that the inferred planets in the disc surrounding AB Aurigae may be evidence of planet formation through GI. △ Less

Submitted 27 March, 2021; originally announced March 2021.

Comments: 12 pages, 5 figures

TOI-1634 b: an Ultra-Short Period Keystone Planet Sitting Inside the M Dwarf Radius Valley

Authors: R. Cloutier, D. Charbonneau, K. G. Stassun, F. Murgas, A. Mortier, R. Massey, J. J. Lissauer, D. W. Latham, J. Irwin, R. D. Haywood, P. Guerra, E. Girardin, S. A. Giacalone, P. Bosch-Cabot, A. Bieryla, J. Winn, C. A. Watson, R. Vanderspek, S. Udry, M. Tamura, A. Sozzetti, A. Shporer, D. Ségransan, S. Seager, A. B. Savel , et al. (41 additional authors not shown)

Abstract: Studies of close-in planets orbiting M dwarfs have suggested that the M dwarf radius valley may be well-explained by distinct formation timescales between enveloped terrestrials, and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by… ▽ More Studies of close-in planets orbiting M dwarfs have suggested that the M dwarf radius valley may be well-explained by distinct formation timescales between enveloped terrestrials, and rocky planets that form at late times in a gas-depleted environment. This scenario is at odds with the picture that close-in rocky planets form with a primordial gaseous envelope that is subsequently stripped away by some thermally-driven mass loss process. These two physical scenarios make unique predictions of the rocky/enveloped transition's dependence on orbital separation such that studying the compositions of planets within the M dwarf radius valley may be able to establish the dominant physics. Here, we present the discovery of one such keystone planet: the ultra-short period planet TOI-1634 b ($P=0.989$ days, $F=121 F_{\oplus}$, $r_p = 1.790^{+0.080}_{-0.081} R_{\oplus}$) orbiting a nearby M2 dwarf ($K_s=8.7$, $R_s=0.45 R_{\odot}$, $M_s=0.50 M_{\odot}$) and whose size and orbital period sit within the M dwarf radius valley. We confirm the TESS-discovered planet candidate using extensive ground-based follow-up campaigns, including a set of 32 precise radial velocity measurements from HARPS-N. We measure a planetary mass of $4.91^{+0.68}_{-0.70} M_{\oplus}$, which makes TOI-1634 b inconsistent with an Earth-like composition at $5.9σ$ and thus requires either an extended gaseous envelope, a large volatile-rich layer, or a rocky portion that is not dominated by iron and silicates to explain its mass and radius. The discovery that the bulk composition of TOI-1634 b is inconsistent with that of the Earth favors the gas-depleted formation mechanism to explain the emergence of the radius valley around M dwarfs with $M_s\lesssim 0.5 M_{\odot}$. △ Less

Submitted 18 May, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 27 pages, 13 figures, accepted to AAS journals. Our time series are included as a csv file in the arXiv source files

Identifying Exoplanets with Deep Learning. IV. Removing Stellar Activity Signals from Radial Velocity Measurements Using Neural Networks

Authors: Zoe L. de Beurs, Andrew Vanderburg, Christopher J. Shallue, Xavier Dumusque, Andrew Collier Cameron, Christopher Leet, Lars A. Buchhave, Rosario Cosentino, Adriano Ghedina, Raphaëlle D. Haywood, Nicholas Langellier, David W. Latham, Mercedes López-Morales, Michel Mayor, Giusi Micela, Timothy W. Milbourne, Annelies Mortier, Emilio Molinari, Francesco Pepe, David F. Phillips, Matteo Pinamonti, Giampaolo Piotto, Ken Rice, Dimitar Sasselov, Alessandro Sozzetti , et al. (2 additional authors not shown)

Abstract: Exoplanet detection with precise radial velocity (RV) observations is currently limited by spurious RV signals introduced by stellar activity. We show that machine learning techniques such as linear regression and neural networks can effectively remove the activity signals (due to starspots/faculae) from RV observations. Previous efforts focused on carefully filtering out activity signals in time… ▽ More Exoplanet detection with precise radial velocity (RV) observations is currently limited by spurious RV signals introduced by stellar activity. We show that machine learning techniques such as linear regression and neural networks can effectively remove the activity signals (due to starspots/faculae) from RV observations. Previous efforts focused on carefully filtering out activity signals in time using modeling techniques like Gaussian Process regression (e.g. Haywood et al. 2014). Instead, we systematically remove activity signals using only changes to the average shape of spectral lines, and no information about when the observations were collected. We trained our machine learning models on both simulated data (generated with the SOAP 2.0 software; Dumusque et al. 2014) and observations of the Sun from the HARPS-N Solar Telescope (Dumusque et al. 2015; Phillips et al. 2016; Collier Cameron et al. 2019). We find that these techniques can predict and remove stellar activity from both simulated data (improving RV scatter from 82 cm/s to 3 cm/s) and from more than 600 real observations taken nearly daily over three years with the HARPS-N Solar Telescope (improving the RV scatter from 1.753 m/s to 1.039 m/s, a factor of ~ 1.7 improvement). In the future, these or similar techniques could remove activity signals from observations of stars outside our solar system and eventually help detect habitable-zone Earth-mass exoplanets around Sun-like stars. △ Less

Submitted 13 June, 2022; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: 28 pages, 14 figures, Accepted for publication in the Astronomical Journal

K2-111: an old system with two planets in near-resonance

Authors: A. Mortier, M. R. Zapatero Osorio, L. Malavolta, Y. Alibert, K. Rice, J. Lillo-Box, A. Vanderburg, M. Oshagh, L. Buchhave, V. Adibekyan, E. Delgado Mena, M. Lopez-Morales, D. Charbonneau, S. G. Sousa, C. Lovis, L. Affer, C. Allende Prieto, S. C. C. Barros, S. Benatti, A. S. Bonomo, W. Boschin, F. Bouchy, A. Cabral, A. Collier Cameron, R. Cosentino , et al. (42 additional authors not shown)

Abstract: This paper reports on the detailed characterisation of the K2-111 planetary system with K2, WASP, and ASAS-SN photometry as well as high-resolution spectroscopic data from HARPS-N and ESPRESSO. The host, K2-111, is confirmed to be a mildly evolved ($\log g=4.17$), iron-poor ([Fe/H]$=-0.46$), but alpha-enhanced ([$α$/Fe]$=0.27$), chromospherically quiet, very old thick disc G2 star. A global fit, p… ▽ More This paper reports on the detailed characterisation of the K2-111 planetary system with K2, WASP, and ASAS-SN photometry as well as high-resolution spectroscopic data from HARPS-N and ESPRESSO. The host, K2-111, is confirmed to be a mildly evolved ($\log g=4.17$), iron-poor ([Fe/H]$=-0.46$), but alpha-enhanced ([$α$/Fe]$=0.27$), chromospherically quiet, very old thick disc G2 star. A global fit, performed by using PyORBIT shows that the transiting planet, K2-111b, orbits with a period $P_b=5.3518\pm0.0004$ d, and has a planet radius of $1.82^{+0.11}_{-0.09}$ R$_\oplus$ and a mass of $5.29^{+0.76}_{-0.77}$ M$_\oplus$, resulting in a bulk density slightly lower than that of the Earth. The stellar chemical composition and the planet properties are consistent with K2-111b being a terrestrial planet with an iron core mass fraction lower than the Earth. We announce the existence of a second signal in the radial velocity data that we attribute to a non-transiting planet, K2-111c, with an orbital period of $15.6785\pm 0.0064$ days, orbiting in near-3:1 mean-motion resonance with the transiting planet, and a minimum planet mass of $11.3\pm1.1$ M$_\oplus$. Both planet signals are independently detected in the HARPS-N and ESPRESSO data when fitted separately. There are potentially more planets in this resonant system, but more well-sampled data are required to confirm their presence and physical parameters. △ Less

Submitted 21 December, 2020; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS on 28 Sept 2020. Paper is 18 pages with an additional 12 pages of supplementary material. Data is available at https://vizier.u-strasbg.fr/viz-bin/VizieR?-source=J/MNRAS/499/5004

Journal ref: Monthly Notices of the Royal Astronomical Society, 2020, Volume 499, Issue 4, pp.5004-5021

An unusually low density ultra-short period super-Earth and three mini-Neptunes around the old star TOI-561

Authors: G. Lacedelli, L. Malavolta, L. Borsato, G. Piotto, D. Nardiello, A. Mortier, M. Stalport, A. Collier Cameron, E. Poretti, L. A. Buchhave, M. López-Morales, V. Nascimbeni, T. G. Wilson, S. Udry, D. W. Latham, A. S. Bonomo, M. Damasso, X. Dumusque, J. M. Jenkins, C. Lovis, K. Rice, D. Sasselov, J. N. Winn, G. Andreuzzi, R. Cosentino , et al. (16 additional authors not shown)

Abstract: Based on HARPS-N radial velocities (RVs) and TESS photometry, we present a full characterisation of the planetary system orbiting the late G dwarf TOI-561. After the identification of three transiting candidates by TESS, we discovered two additional external planets from RV analysis. RVs cannot confirm the outer TESS transiting candidate, which would also make the system dynamically unstable. We d… ▽ More Based on HARPS-N radial velocities (RVs) and TESS photometry, we present a full characterisation of the planetary system orbiting the late G dwarf TOI-561. After the identification of three transiting candidates by TESS, we discovered two additional external planets from RV analysis. RVs cannot confirm the outer TESS transiting candidate, which would also make the system dynamically unstable. We demonstrate that the two transits initially associated with this candidate are instead due to single transits of the two planets discovered using RVs. The four planets orbiting TOI-561 include an ultra-short period (USP) super-Earth (TOI-561 b) with period $P_{\rm b} = 0.45$ d, mass $M_{\rm b} =1.59 \pm 0.36$ M$_\oplus$ and radius $R_{\rm b}=1.42 \pm 0.07$ R$_\oplus$, and three mini-Neptunes: TOI-561 c, with $P_{\rm c} = 10.78$ d, $M_{\rm c} = 5.40 \pm 0.98$ M$_\oplus$, $R_{\rm c}= 2.88 \pm 0.09$ R$_\oplus$; TOI-561 d, with $P_{\rm d} = 25.6$ d, $M_{\rm d} = 11.9 \pm 1.3$ M$_\oplus$, $R_{\rm d} = 2.53 \pm 0.13$ R$_\oplus$; and TOI-561 e, with $P_{\rm e} = 77.2$ d, $M_{\rm e} = 16.0 \pm 2.3$ M$_\oplus$, $R_{\rm e} = 2.67 \pm 0.11$ R$_\oplus$. Having a density of $3.0 \pm 0.8$ g cm$^{-3}$, TOI-561 b is the lowest density USP planet known to date. Our N-body simulations confirm the stability of the system and predict a strong, anti-correlated, long-term transit time variation signal between planets d and e. The unusual density of the inner super-Earth and the dynamical interactions between the outer planets make TOI-561 an interesting follow-up target. △ Less

Submitted 27 November, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 20 pages, 12 figures. Accepted for publication in MNRAS

Three Years of HARPS-N High-Resolution Spectroscopy and Precise Radial Velocity Data for the Sun

Authors: X. Dumusque, M. Cretignier, D. Sosnowska, N. Buchschacher, C. Lovis, D. F. Phillips, F. Pepe, F. Alesina, L. A. Buchhave, J. Burnier, M. Cecconi, H. M. Cegla, R. Cloutier, A. Collier Cameron, R. Cosentino, A. Ghedina, M. Gonzalez, R. D. Haywood, D. W. Latham, M. Lodi, M. Lopez-Morales, J. Maldonado, L. Malavolta, G. Micela, E. Molinari , et al. (13 additional authors not shown)

Abstract: The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such high-cadence, long-baseline data set is crucial for understanding spurious radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub-\ms\,systematics that needed to be corrected for. The goal of this manuscript is to i) present… ▽ More The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such high-cadence, long-baseline data set is crucial for understanding spurious radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub-\ms\,systematics that needed to be corrected for. The goal of this manuscript is to i) present a new data reduction software for HARPS-N, ii) demonstrate the improvement brought by this new software on the first three years of the HARPS-N solar data set, and iii) release all the obtained solar products, from extracted spectra to precise radial velocities. To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimized recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS. The most significant improvement brought by the new DRS is a strong decrease in the day-to-day radial-velocity scatter, from 1.27 to 1.07\ms; this is thanks to a more robust method to derive wavelength solutions, but also to the use of calibrations closer in time. The newly derived solar radial-velocities are also better correlated with the chromospheric activity level of the Sun on the long-term, with a Pearson correlation coefficient of 0.93 compared to 0.77 before, which is expected from our understanding of stellar signals. Finally, we also discuss how HARPS-N spectral ghosts contaminate the measurement of the calcium activity index, and present an efficient technique to derive an index free of instrumental systematics. This paper presents a new data reduction software for HARPS-N, and demonstrates its improvements [...] △ Less

Submitted 6 May, 2021; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 20 pages, 15 figures, 2 tables, version accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 648, A103 (2021)

The observational impact of dust trapping in self-gravitating discs

Authors: James Cadman, Cassandra Hall, Ken Rice, Tim J. Harries, Pamela D. Klaassen

Abstract: We present a 3D semi-analytic model of self-gravitating discs, and include a prescription for dust trapping in the disc spiral arms. Using Monte-Carlo radiative transfer we produce synthetic ALMA observations of these discs. In doing so we demonstrate that our model is capable of producing observational predictions, and able to model real image data of potentially self-gravitating discs. For a dis… ▽ More We present a 3D semi-analytic model of self-gravitating discs, and include a prescription for dust trapping in the disc spiral arms. Using Monte-Carlo radiative transfer we produce synthetic ALMA observations of these discs. In doing so we demonstrate that our model is capable of producing observational predictions, and able to model real image data of potentially self-gravitating discs. For a disc to generate spiral structure that would be observable with ALMA requires that the disc's dust mass budget is dominated by millimetre and centimetre-sized grains. Discs in which grains have grown to the grain fragmentation threshold may satisfy this criterion, thus we predict that signatures of gravitational instability may be detectable in discs of lower mass than has previously been suggested. For example, we find that discs with disc-to-star mass ratios as low as $0.10$ are capable of driving observable spiral arms. Substructure becomes challenging to detect in discs where no grain growth has occurred or in which grain growth has proceeded well beyond the grain fragmentation threshold. We demonstrate how we can use our model to retrieve information about dust trapping and grain growth through multi-wavelength observations of discs, and using estimates of the opacity spectral index. Applying our disc model to the Elias 27, WaOph 6 and IM Lup systems we find gravitational instability to be a plausible explanation for the observed substructure in all 3 discs, if sufficient grain growth has indeed occurred. △ Less

Submitted 22 August, 2020; originally announced August 2020.

Comments: 19 pages, 21 figures, accepted for publication in MNRAS

Detection Limits of Low-mass, Long-period Exoplanets Using Gaussian Processes Applied to HARPS-N Solar RVs

Authors: N. Langellier, T. W. Milbourne, D. F. Phillips, R. D. Haywood, S. H. Saar, A. Mortier, L. Malavolta, S. Thompson, A. Collier Cameron, X. Dumusque, H. M. Cegla, D. W. Latham, J. Maldonado, C. A. Watson, N. Buchschacher, M. Cecconi, D. Charbonneau, R. Cosentino, A. Ghedina, M. Gonzalez, C-H. Li, M. Lodi, M. López-Morales, G. Micela, E. Molinari , et al. (7 additional authors not shown)

Abstract: Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression… ▽ More Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s$^{-1}$. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass ($\sim$ 0.5 m s$^{-1}$) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus ($\sim$ 0.1 m s$^{-1}$) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets. △ Less

Submitted 2 June, 2021; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 12 pages, 4 figures, 1 table, accepted in AJ

Journal ref: AJ 161 287 (2021)

The search for disks or planetary objects around directly imaged companions: A candidate around DH Tau B

Authors: C. Lazzoni, A. Zurlo, S. Desidera, D. Mesa, C. Fontanive, M. Bonavita, S. Ertel, K. Rice, A. Vigan, A. Boccaletti, M. Bonnefoy, G. Chauvin, P. Delorme, R. Gratton, M. Houllé, A. L. Maire, M. Meyer, E. Rickman, E. A. Spalding, R. Asensio-Torres, M. Langlois, A. Müller, J-L. Baudino, J. -L. Beuzit, B. Biller , et al. (23 additional authors not shown)

Abstract: In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. In this paper, we focus our attention on substellar companions detected with the direct imaging technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery w… ▽ More In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. In this paper, we focus our attention on substellar companions detected with the direct imaging technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery would shed light on many unresolved questions, particularly with regard to their possible formation mechanisms. To reveal bound features of directly imaged companions we need to suppress the contribution from the source itself. Therefore, we developed a method based on the negative fake companion (NEGFC) technique that first estimates the position in the field of view (FoV) and the flux of the imaged companion, then subtracts a rescaled model point spread function (PSF) from the imaged companion. Next it performs techniques, such as angular differential imaging (ADI), to further remove quasi-static patterns of the star. We applied the method to the sample of substellar objects observed with SPHERE during the SHINE GTO survey. Among the 27 planets and brown dwarfs we analyzed, we detected a possible point source close to DH Tau B. This candidate companion was detected in four different SPHERE observations, with an estimated mass of $\sim 1$ M\textsubscript{Jup}, and a mass ratio with respect to the brown dwarf of $1/10$. This binary system, if confirmed, would be the first of its kind, opening up interesting questions for the formation mechanism, evolution, and frequency of such pairs. In order to address the latter, the residuals and contrasts reached for 25 companions in the sample of substellar objects observed with SPHERE were derived. If the DH Tau Bb companion is real, the binary fraction obtained is $\sim 7\%$, which is in good agreement with the results obtained for field brown dwarfs. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Comments: 19 pages, 17 figures

Journal ref: A&A 641, A131 (2020)

The SPHERE infrared survey for exoplanets (SHINE). III. The demographics of young giant exoplanets below 300 au with SPHERE

Authors: A. Vigan, C. Fontanive, M. Meyer, B. Biller, M. Bonavita, M. Feldt, S. Desidera, G. -D. Marleau, A. Emsenhuber, R. Galicher, K. Rice, D. Forgan, C. Mordasini, R. Gratton, H. Le Coroller, A. -L. Maire, F. Cantalloube, G. Chauvin, A. Cheetham, J. Hagelberg, A. -M. Lagrange, M. Langlois, M. Bonnefoy, J. -L. Beuzit, A. Boccaletti , et al. (86 additional authors not shown)

Abstract: The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses betwee… ▽ More The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 24 pages, 14 figures, 3 tables. Accepted for publication in A&A

Journal ref: A&A 651, A72 (2021)

TOI-1235 b: a keystone super-Earth for testing radius valley emergence models around early M dwarfs

Authors: Ryan Cloutier, Joseph E. Rodriguez, Jonathan Irwin, David Charbonneau, Keivan G. Stassun, Annelies Mortier, David W. Latham, Howard Isaacson, Andrew W. Howard, Stéphane Udry, Thomas G. Wilson, Christopher A. Watson, Matteo Pinamonti, Florian Lienhard, Paolo Giacobbe, Pere Guerra, Karen A. Collins, Allyson Beiryla, Gilbert A. Esquerdo, Elisabeth Matthews, Rachel A. Matson, Steve B. Howell, Elise Furlan, Ian J. M. Crossfield, Jennifer G. Winters , et al. (63 additional authors not shown)

Abstract: Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we presen… ▽ More Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we present the confirmation of TOI-1235 b ($P=3.44$ days, $r_p=1.738^{+0.087}_{-0.076}$ R$_{\oplus}$), a planet whose size and period are intermediate between the competing model predictions thus making the system an important test case for emergence models of the rocky/non-rocky transition around early M dwarfs ($R_s=0.630\pm 0.015$ R$_{\odot}$, $M_s=0.640\pm 0.016$ M$_{\odot}$). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial-velocities from HARPS-N and HIRES. We measure a planet mass of $6.91^{+0.75}_{-0.85}$ M$_{\oplus}$, which implies an iron core mass fraction of $20^{+15}_{-12}$% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like and we constrain a H/He envelope mass fraction to be $<0.5$% at 90% confidence. Our results are consistent with model predictions from thermally-driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remain efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically-determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin ($P=21.8^{+0.9}_{-0.8}$ days, $m_p\sin{i}=13.0^{+3.8}_{-5.3}$ M$_{\oplus}$) that cannot be firmly ruled out by our data. △ Less

Submitted 21 May, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Comments: Accepted to The Astronomical Journal. 8 figures & 5 tables. Table 2 is provided in the arXiv source code

A pair of TESS planets spanning the radius valley around the nearby mid-M dwarf LTT 3780

Authors: Ryan Cloutier, Jason D. Eastman, Joseph E. Rodriguez, Nicola Astudillo-Defru, Xavier Bonfils, Annelies Mortier, Christopher A. Watson, Manu Stalport, Matteo Pinamonti, Florian Lienhard, Avet Harutyunyan, Mario Damasso, David W. Latham, Karen A. Collins, Robert Massey, Jonathan Irwin, Jennifer G. Winters, David Charbonneau, Carl Ziegler, Elisabeth Matthews, Ian J. M. Crossfield, Laura Kreidberg, Samuel N. Quinn, George Ricker, Roland Vanderspek , et al. (62 additional authors not shown)

Abstract: We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbita… ▽ More We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley. △ Less

Submitted 12 May, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: Accepted to AJ. 8 figures, 6 tables. CSV file of the RV measurements (i.e. Table 2) are included in the source code

Massive discs around low-mass stars

Authors: Thomas J. Haworth, James Cadman, Farzana Meru, Cassandra Hall, Emma Albertini, Duncan Forgan, Ken Rice, James E. Owen

Abstract: We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. T… ▽ More We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation. △ Less

Submitted 27 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, Accepted for publication MNRAS

Fragmentation favoured in discs around higher mass stars

Authors: James Cadman, Ken Rice, Cassandra Hall, Thomas J. Haworth, Beth Biller

Abstract: We investigate how a protoplanetary disc's susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increas… ▽ More We investigate how a protoplanetary disc's susceptibility to gravitational instabilities and fragmentation depends on the mass of its host star. We use 1D disc models in conjunction with 3D SPH simulations to determine the critical disc-to-star mass ratios at which discs become unstable against fragmentation, finding that discs become increasingly prone to the effects of self-gravity as we increase the host star mass. The actual limit for stability is sensitive to the disc temperature, so if the disc is optically thin stellar irradiation can dramatically stabilise discs against gravitational instability. However, even when this is the case we find that discs around $2$M$_{\odot}$ stars are prone to fragmentation, which will act to produce wide-orbit giant planets and brown dwarfs. The consequences of this work are two-fold: that low mass stars could in principle support high disc-to-star mass ratios, and that higher mass stars have discs that are more prone to fragmentation, which is qualitatively consistent with observations that favour high-mass wide-orbit planets around higher mass stars. We also find that the initial masses of these planets depends on the temperature in the disc at large radii, which itself depends on the level of stellar irradiation. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS

An ultra-short period rocky super-Earth orbiting the G2-star HD 80653

Authors: G. Frustagli, E. Poretti, T. Milbourne, L. Malavolta, A. Mortier, Vikash Singh, A. S. Bonomo, L. A. Buchhave, L. Zeng, A. Vanderburg, S. Udry, G. Andreuzzi, A. Collier-Cameron, R. Cosentino, M. Damasso, A. Ghedina, A. Harutyunyan, R. D. Haywood, D. W. Latham, M. López-Morales, V. Lorenzi, A. F. Martinez Fiorenzano, M. Mayor, G. Micela, E. Molinari , et al. (4 additional authors not shown)

Abstract: Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise mas… ▽ More Ultra-short period (USP) planets are a class of exoplanets with periods shorter than one day. The origin of this sub-population of planets is still unclear, with different formation scenarios highly dependent on the composition of the USP planets. A better understanding of this class of exoplanets will, therefore, require an increase in the sample of such planets that have accurate and precise masses and radii, which also includes estimates of the level of irradiation and information about possible companions. Here we report a detailed characterization of a USP planet around the solar-type star HD 80653 $\equiv$ EP 251279430 using the K2 light curve and 108 precise radial velocities obtained with the HARPS-N spectrograph, installed on the Telescopio Nazionale Galileo. From the K2 C16 data, we found one super-Earth planet ($R_{b}=1.613\pm0.071 R_{\oplus}$) transiting the star on a short-period orbit ($P_{\rm b}=0.719573\pm0.000021$ d). From our radial velocity measurements, we constrained the mass of HD 80653 b to $M_{b}=5.60\pm0.43 M_{\oplus}$. We also detected a clear long-term trend in the radial velocity data. We derived the fundamental stellar parameters and determined a radius of $R_{\star}=1.22\pm0.01 R_{\odot}$ and mass of $M_{\star}=1.18\pm0.04 M_{\odot}$, suggesting that HD 80653, has an age of $2.7\pm1.2$ Gyr. The bulk density ($ρ_{b} = 7.4 \pm 1.1$ g cm$^{-3}$) of the planet is consistent with an Earth-like composition of rock and iron with no thick atmosphere. Our analysis of the K2 photometry also suggests hints of a shallow secondary eclipse with a depth of 8.1$\pm$3.7 ppm. Flux variations along the orbital phase are consistent with zero. The most important contribution might come from the day-side thermal emission from the surface of the planet at $T\sim3480$ K. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 15 pages, 12 figures; accepted by A&A

An extremely low-density and temperate giant exoplanet

Authors: A. Santerne, L. Malavolta, M. R. Kosiarek, F. Dai, C. D. Dressing, X. Dumusque, N. C. Hara, T. A. Lopez, A. Mortier, A. Vanderburg, V. Adibekyan, D. J. Armstrong, D. Barrado, S. C. C. Barros, D. Bayliss, D. Berardo, I. Boisse, A. S. Bonomo, F. Bouchy, D. J. A. Brown, L. A. Buchhave, R. P. Butler, A. Collier Cameron, R. Cosentino, J. D. Crane , et al. (46 additional authors not shown)

Abstract: Transiting extrasolar planets are key objects in the study of the formation, migration, and evolution of planetary systems. In particular, the exploration of the atmospheres of giant planets, through transmission spectroscopy or direct imaging, has revealed a large diversity in their chemical composition and physical properties. Studying these giant planets allows one to test the global climate mo… ▽ More Transiting extrasolar planets are key objects in the study of the formation, migration, and evolution of planetary systems. In particular, the exploration of the atmospheres of giant planets, through transmission spectroscopy or direct imaging, has revealed a large diversity in their chemical composition and physical properties. Studying these giant planets allows one to test the global climate models that are used for the Earth and other solar system planets. However, these studies are mostly limited either to highly-irradiated transiting giant planets or directly-imaged giant planets at large separations. Here we report the physical characterisation of the planets in a bright multi-planetary system (HIP41378) in which the outer planet, HIP41378 f is a Saturn-sized planet (9.2 $\pm$ 0.1 R$_\oplus$) with an anomalously low density of 0.09 $\pm$ 0.02 g cm$^{-3}$ that is not yet understood. Its equilibrium temperature is about 300 K. Therefore, it represents a planet with a mild temperature, in between the hot Jupiters and the colder giant planets of the Solar System. It opens a new window for atmospheric characterisation of giant exoplanets with a moderate irradiation, with the next-generation space telescopes such as JWST and ARIEL as well as the extremely-large ground-based telescopes. HIP41378 f is thus an important laboratory to understand the effect of the irradiation on the physical properties and chemical composition of the atmosphere of planets. △ Less

Submitted 17 November, 2019; originally announced November 2019.

Comments: Preprint submitted to Nature Astronomy. The results have not been peer-reviewed yet. Supplementary informations available as ancillary file

Using HARPS-N to characterise the long-period planets in the PH-2 and Kepler-103 systems

Authors: Sophie C. Dubber, Annelies Mortier, Ken Rice, Chantanelle Nava, Luca Malavolta, Helen Giles, Adrien Coffinet, David Charbonneau, Andrew Vanderburg, Aldo S. Bonomo, Walter Boschin, Lars A. Buchhave, Andrew Collier Cameron, Rosario Cosentino, Xavier Dumusque, Adriano Ghedina, Avet Harutyunyan, Raphaelle D. Haywood, David Latham, Mercedes Lopez-Morales, Giusi Micela, Emilio Molinari, Francesco A. Pepe, David Phillips, Giampaolo Piotto , et al. (4 additional authors not shown)

Abstract: We present confirmation of the planetary nature of PH-2b, as well as the first mass estimates for the two planets in the Kepler-103 system. PH-2b and Kepler-103c are both long-period and transiting, a sparsely-populated category of exoplanet. We use {\it Kepler} light-curve data to estimate a radius, and then use HARPS-N radial velocities to determine the semi-amplitude of the stellar reflex motio… ▽ More We present confirmation of the planetary nature of PH-2b, as well as the first mass estimates for the two planets in the Kepler-103 system. PH-2b and Kepler-103c are both long-period and transiting, a sparsely-populated category of exoplanet. We use {\it Kepler} light-curve data to estimate a radius, and then use HARPS-N radial velocities to determine the semi-amplitude of the stellar reflex motion and, hence, the planet mass. For PH-2b we recover a 3.5-$σ$ mass estimate of $M_p = 109^{+30}_{-32}$ M$_\oplus$ and a radius of $R_p = 9.49\pm0.16$ R$_\oplus$. This means that PH-2b has a Saturn-like bulk density and is the only planet of this type with an orbital period $P > 200$ days that orbits a single star. We find that Kepler-103b has a mass of $M_{\text{p,b}} = 11.7^{+4.31}_{-4.72}$ M$_{\oplus}$ and Kepler-103c has a mass of $M_{\text{p,c}} = 58.5^{+11.2}_{-11.4}$ M$_{\oplus}$. These are 2.5$σ$ and 5$σ$ results, respectively. With radii of $R_{\text{p,b}} = 3.49^{+0.06}_{-0.05}$ R$_\oplus$, and $R_{\text{p,c}} = 5.45^{+0.18}_{-0.17}$ R$_\oplus$, these results suggest that Kepler-103b has a Neptune-like density, while Kepler-103c is one of the highest density planets with a period $P > 100$ days. By providing high-precision estimates for the masses of the long-period, intermediate-mass planets PH-2b and Kepler-103c, we increase the sample of long-period planets with known masses and radii, which will improve our understanding of the mass-radius relation across the full range of exoplanet masses and radii. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

An 11 Earth-Mass, Long-Period Sub-Neptune Orbiting a Sun-like Star

Authors: Andrew W. Mayo, Vinesh M. Rajpaul, Lars A. Buchhave, Courtney D. Dressing, Annelies Mortier, Li Zeng, Charles D. Fortenbach, Suzanne Aigrain, Aldo S. Bonomo, Andrew Collier Cameron, David Charbonneau, Adrien Coffinet, Rosario Cosentino, Mario Damasso, Xavier Dumusque, A. F. Martinez Fiorenzano, Raphaëlle D. Haywood, David W. Latham, Mercedes López-Morales, Luca Malavolta, Giusi Micela, Emilio Molinari, Logan Pearce, Francesco Pepe, David Phillips , et al. (5 additional authors not shown)

Abstract: Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) w… ▽ More Although several thousands of exoplanets have now been detected and characterized, observational biases have led to a paucity of long-period, low-mass exoplanets with measured masses and a corresponding lag in our understanding of such planets. In this paper we report the mass estimation and characterization of the long-period exoplanet Kepler-538b. This planet orbits a Sun-like star (V = 11.27) with M_* = 0.892 +/- (0.051, 0.035) M_sun and R_* = 0.8717 +/- (0.0064, 0.0061) R_sun. Kepler-538b is a 2.215 +/- (0.040, 0.034) R_earth sub-Neptune with a period of P = 81.73778 +/- 0.00013 d. It is the only known planet in the system. We collected radial velocity (RV) observations with HIRES on Keck I and HARPS-N on the TNG. We characterized stellar activity by a Gaussian process with a quasi-periodic kernel applied to our RV and cross correlation function full width at half maximum (FWHM) observations. By simultaneously modeling Kepler photometry, RV, and FWHM observations, we found a semi-amplitude of K = 1.68 +/- (0.39, 0.38) m s^-1 and a planet mass of M_p = 10.6 +/- (2.5, 2.4) M_earth. Kepler-538b is the smallest planet beyond P = 50 d with an RV mass measurement. The planet likely consists of a significant fraction of ices (dominated by water ice), in addition to rocks/metals, and a small amount of gas. Sophisticated modeling techniques such as those used in this paper, combined with future spectrographs with ultra high-precision and stability will be vital for yielding more mass measurements in this poorly understood exoplanet regime. This in turn will improve our understanding of the relationship between planet composition and insolation flux and how the rocky to gaseous transition depends on planetary equilibrium temperature. △ Less

Submitted 11 December, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: Published in AJ, 17 pages, 8 figures, 3 tables

Journal ref: 2019 AJ 158 4

Temporal evolution and correlations of optical activity indicators measured in Sun-as-a-star observations

Authors: J. Maldonado, D. F. Phillips, X. Dumusque, A. Collier Cameron, R. D. Haywood, A. F. Lanza, G. Micela, A. Mortier, S. H. Saar, A. Sozzetti, K. Rice, T. Milbourne, M. Cecconi, H. M. Cegla, R. Cosentino, J. Costes, A. Ghedina, M. Gonzalez, J. Guerra, N. Hernández, C. -H. Li, M. Lodi, L. Malavolta, E. Molinari, F. Pepe , et al. (7 additional authors not shown)

Abstract: (Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$δ$. The der… ▽ More (Abridged) We perform a detailed study of the main optical activity indicators (Ca II H & K, Balmer lines, Na I D$_{\rm 1}$ D$_{\rm 2}$, and He I D$_{\rm 3}$) measured for the Sun using the data provided by the HARPS-N solar-telescope feed at the Telescopio Nazionale Galileo. The value of the solar rotation period is found in all the activity indicators, with the only exception being H$δ$. The derived values vary from 26.29 days (H$γ$ line) to 31.23 days (He I). From an analysis of sliding periodograms we find that in most of the activity indicators the spectral power is split into several "bands" of periods around 26 and 30 days, that might be explained by the migration of active regions between the equator and a latitude of $\sim$ 30$^{\circ}$, spot evolution or a combination of both effects. In agreement with previous works a typical lifetime of active regions of $\sim$ ten rotation periods is inferred from the pooled variance diagrams. We find that H$α$, H$β$, H$γ$, H$ε$, and He I show a significant correlation with the S index. Significant correlations between the contrast, bisector span, and the heliocentric radial velocity with the activity indexes are also found. We show that the full width at half maximum, the bisector, and the disc-integrated magnetic field correlate with the radial velocity variations. The correlation of the S index and H$α$ changes with time, increasing with larger sun spot numbers and solar irradiance. A similar tendency with the S index - radial velocity correlation is also present in the data. Our results are consistent with a scenario in which higher activity favours the correlation between the S index and the H$α$ activity indicators and between the S index and radial velocity variations. △ Less

Submitted 7 June, 2019; originally announced June 2019.

Comments: Accepted by Astronomy & Astrophysics

Journal ref: A&A 627, A118 (2019)

Three years of Sun-as-a-star radial-velocity observations on the approach to solar minimum

Authors: A. Collier Cameron, A. Mortier, D. Phillips, X. Dumusque, R. D. Haywood, N. Langellier, C. A. Watson, H. M. Cegla, J. Costes, D. Charbonneau, A. Coffinet, D. W. Latham, M. Lopez-Morales, L. Malavolta, J. Maldonado, G. Micela, T. Milbourne, E. Molinari, S. H. Saar, S. Thompson, N. Buchschacher, M. Cecconi, R. Cosentino, A. Ghedina, A. Glenday , et al. (11 additional authors not shown)

Abstract: The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-veloci… ▽ More The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity. △ Less

Submitted 27 April, 2019; originally announced April 2019.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

A high binary fraction for the most massive close-in giant planets and brown dwarf desert members

Authors: C. Fontanive, K. Rice, M. Bonavita, E. Lopez, K. Muzic, B. Biller

Abstract: Stellar multiplicity is believed to influence planetary formation and evolution, although the precise nature and extent of this role remain ambiguous. We present a study aimed at testing the role of stellar multiplicity in the formation and/or evolution of the most massive, close-in planetary and substellar companions. Using direct imaging observations, as well as the Gaia DR2 catalogue, we search… ▽ More Stellar multiplicity is believed to influence planetary formation and evolution, although the precise nature and extent of this role remain ambiguous. We present a study aimed at testing the role of stellar multiplicity in the formation and/or evolution of the most massive, close-in planetary and substellar companions. Using direct imaging observations, as well as the Gaia DR2 catalogue, we searched for wide binary companions to 38 stars hosting massive giant planets or brown dwarfs (M > 7 MJup) on orbits shorter than ~1 AU. We report the discovery of a new component in the WASP-14 system, and present an independent confirmation of a comoving companion to WASP-18. From a robust Bayesian statistical analysis, we derived a binary fraction of 79.0+13.2-14.7% between 20-10,000 AU for our sample, twice as high as for field stars with a 3-σ significance. This binary frequency was found to be larger than for lower-mass planets on similar orbits, and we observed a marginally higher binary rate for inner companions with periods shorter than 10 days. These results demonstrate that stellar companions greatly influence the formation and/or evolution of these systems, suggesting that the role played by binary companions becomes more important for higher-mass planets, and that this trend may be enhanced for systems with tighter orbits. Our analysis also revealed a peak in binary separation at 250 AU, highlighting a shortfall of close binaries among our sample. This indicates that the mechanisms affecting planet and brown dwarf formation or evolution in binaries must operate from wide separations, although we found that the Kozai-Lidov mechanism is unlikely to be the dominant underlying process. We conclude that binarity plays a crucial role in the existence of very massive short-period giant planets and brown dwarf desert inhabitants, which are almost exclusively observed in multiple systems. △ Less

Submitted 28 March, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS. 30 pages, 20 figures. Updated to include proof corrections

HARPS-N Solar Radial-Velocity Variations Are Dominated By Large, Bright Magnetic Regions

Authors: T. W. Milbourne, R. D. Haywood, D. F. Phillips, S. H. Saar, H. M. Cegla, A. C. Cameron, J. Costes, X. Dumusque, N. Langellier, D. W. Latham, J. Maldonado, L. Malavolta, A. Mortier, M. L. Palumbo Iii, S. Thompson, C. A. Watson, F. Bouchy, N. Buchschacher, M. Cecconi, D. Charbonneau, R. Cosentino, A. Ghedina, A. G. Glenday, M. Gonzalez, C-H. Li , et al. (15 additional authors not shown)

Abstract: State of the art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3-year period during the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (H… ▽ More State of the art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3-year period during the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and $\log{R'_{\rm HK}}$. The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N $\log{R'_{\rm HK}}$ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than $60 \ \rm Mm^2$ do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 16 pages, 7 figures, accepted by ApJ

So close, so different: characterization of the K2-36 planetary system with HARPS-N

Authors: M. Damasso, L. Zeng, L. Malavolta, A. Mayo, A. Sozzetti, A. Mortier, L. A. Buchhave, A. Vanderburg, M. Lopez-Morales, A. S. Bonomo, A. C. Cameron, A. Coffinet, P. Figueira, D. W. Latham, M. Mayor, E. Molinari, F. Pepe, D. F. Phillips, E. Poretti, K. Rice, S. Udry, C. A. Watson

Abstract: K2-36 is a K dwarf orbited by two small ($R_{\rm b}=1.43\pm0.08$ $R_\oplus$ and $R_{\rm c}=3.2\pm0.3$ $R_\oplus$), close-in ($a_{\rm b}$=0.022 AU and $a_{\rm c}$=0.054 AU) transiting planets discovered by Kepler/K2. They are representatives of two families of small planets ($R_{\rm p}$<4 $R_\oplus$) recently emerged from the analysis of Kepler data, with likely a different structure, composition a… ▽ More K2-36 is a K dwarf orbited by two small ($R_{\rm b}=1.43\pm0.08$ $R_\oplus$ and $R_{\rm c}=3.2\pm0.3$ $R_\oplus$), close-in ($a_{\rm b}$=0.022 AU and $a_{\rm c}$=0.054 AU) transiting planets discovered by Kepler/K2. They are representatives of two families of small planets ($R_{\rm p}$<4 $R_\oplus$) recently emerged from the analysis of Kepler data, with likely a different structure, composition and evolutionary pathways. We revise the fundamental stellar parameters and the sizes of the planets, and provide the first measurement of their masses and bulk densities, which we use to infer their structure and composition. We observed K2-36 with the HARPS-N spectrograph over $\sim$3.5 years, collecting 81 useful radial velocity measurements. The star is active, with evidence for increasing levels of magnetic activity during the observing time span. The radial velocity scatter is $\sim$17 \ms due to the stellar activity contribution, which is much larger that the semi-amplitudes of the planetary signals. We tested different methods for mitigating the stellar activity contribution to the radial velocity time variations and measuring the planet masses with good precision. We found that K2-36 is likely a $\sim$1 Gyr old system, and by treating the stellar activity through a Gaussian process regression, we measured the planet masses $m_{\rm b}$=3.9$\pm$1.1 $M_\oplus$ and $m_{\rm c}$=7.8$\pm$2.3 $M_\oplus$. The derived planet bulk densities $ρ_{\rm b}$=7.2$^{+2.5}_{-2.1}$ $g/cm^{3}$ and $ρ_{\rm c}$=1.3$^{+0.7}_{-0.5}$ $g/cm^{3}$ point out that K2-36\,b has a rocky, Earth-like composition, and K2-36\,c is a low-density sub-Neptune. Composed of two planets with similar orbital separations but different densities, K2-36 represents an optimal laboratory for testing the role of the atmospheric escape in driving the evolution of close-in, low-mass planets after $\sim$1 Gyr from their formation. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: Accepted for publication on Astronomy $\&$ Astrophysics

Journal ref: A&A 624, A38 (2019)

A giant impact as the likely origin of different twins in the Kepler-107 exoplanet system

Authors: Aldo S. Bonomo, Li Zeng, Mario Damasso, Zoë M. Leinhardt, Anders B. Justesen, Eric Lopez, Mikkel N. Lund, Luca Malavolta, Victor Silva Aguirre, Lars A. Buchhave, Enrico Corsaro, Thomas Denman, Mercedes Lopez-Morales, Sean M. Mills, Annelies Mortier, Ken Rice, Alessandro Sozzetti, Andrew Vanderburg, Laura Affer, Torben Arentoft, Mansour Benbakoura, François Bouchy, Jørgen Christensen-Dalsgaard, Andrew Collier Cameron, Rosario Cosentino , et al. (27 additional authors not shown)

Abstract: Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the pl… ▽ More Measures of exoplanet bulk densities indicate that small exoplanets with radius less than 3 Earth radii ($R_\oplus$) range from low-density sub-Neptunes containing volatile elements to higher density rocky planets with Earth-like or iron-rich (Mercury-like) compositions. Such astonishing diversity in observed small exoplanet compositions may be the product of different initial conditions of the planet-formation process and/or different evolutionary paths that altered the planetary properties after formation. Planet evolution may be especially affected by either photoevaporative mass loss induced by high stellar X-ray and extreme ultraviolet (XUV) flux or giant impacts. Although there is some evidence for the former, there are no unambiguous findings so far about the occurrence of giant impacts in an exoplanet system. Here, we characterize the two innermost planets of the compact and near-resonant system Kepler-107. We show that they have nearly identical radii (about $1.5-1.6~R_\oplus$), but the outer planet Kepler-107c is more than twice as dense (about $12.6~\rm g\,cm^{-3}$) as the innermost Kepler-107b (about $5.3~\rm g\,cm^{-3}$). In consequence, Kepler-107c must have a larger iron core fraction than Kepler-107b. This imbalance cannot be explained by the stellar XUV irradiation, which would conversely make the more-irradiated and less-massive planet Kepler-107b denser than Kepler-107c. Instead, the dissimilar densities are consistent with a giant impact event on Kepler-107c that would have stripped off part of its silicate mantle. This hypothesis is supported by theoretical predictions from collisional mantle stripping, which match the mass and radius of Kepler-107c. △ Less

Submitted 4 February, 2019; originally announced February 2019.

Comments: Published in Nature Astronomy on 4 February 2019, 35 pages including Supplementary Information material

Journal ref: Nature Astronomy 02/2019