Showing 1–50 of 65 results for author: Marleau, G

Accreting companion occurrence rates using a new method to compute emission-line survey sensitivity: Application to the H-alpha Giant Accreting Protoplanet Survey (GAPlanetS)

Authors: Cailin Plunkett, Katherine B. Follette, Gabriel-Dominique Marleau, Eric Nielsen

Abstract: A key scientific goal of exoplanet surveys is to characterize the underlying population of planets in the local galaxy. In particular, the properties of accreting protoplanets can inform the rates and physical processes of planet formation. We develop a novel method to compute sensitivity to protoplanets in emission-line direct-imaging surveys, enabling estimates of protoplanet population properti… ▽ More A key scientific goal of exoplanet surveys is to characterize the underlying population of planets in the local galaxy. In particular, the properties of accreting protoplanets can inform the rates and physical processes of planet formation. We develop a novel method to compute sensitivity to protoplanets in emission-line direct-imaging surveys, enabling estimates of protoplanet population properties under various planetary accretion and formation theories. In this work, we specialize to the case of H-alpha and investigate three formation models governing the planetary-mass-to-mass-accretion-rate power law, and two accretion models that describe the scaling between total accretion luminosity and observable H-alpha line luminosity. We apply our method to the results of the Magellan Giant Accreting Protoplanet Survey (GAPlanetS) to place the first constraints on accreting companion occurrence rates in systems with transitional circumstellar disks. We compute the posterior probability for transitional disk systems to host an accreting companion (-8< log MMdot MJ^2/yr < -2) within 2 arcseconds (~200 au). Across accretion models, we find consistent accreting companion rates, with median and one-sigma credible intervals of 0.15 (+0.18, -0.10) and 0.19 (+0.23, -0.12). Our technique enables studying protoplanet populations under flexible assumptions about planet formation. This formalism provides the statistical underpinning necessary for protoplanet surveys to discriminate among formation and accretion theories for planets and brown dwarfs. △ Less

Submitted 27 August, 2024; v1 submitted 2 August, 2024; originally announced August 2024.

Comments: Submitted to ApJ. Comments welcome. 18 pages, 7 figures. Updated bibliography

Analyzing JWST/NIRSpec Hydrogen Line Detections at TWA 27B: Constraining Accretion Properties and Geometry

Authors: Yuhiko Aoyama, Gabriel-Dominique Marleau, Jun Hashimoto

Abstract: Hydrogen lines from forming planets are crucial for understanding planet formation. However, the number of planetary hydrogen line detections is still limited. Recent JWST/NIRSpec observations have detected Paschen and Brackett hydrogen lines at TWA 27 B (2M1207b). TWA 27 B is classified as a planetary-mass companison (PMC) rather than a planet due to its large mass ratio to the central object (… ▽ More Hydrogen lines from forming planets are crucial for understanding planet formation. However, the number of planetary hydrogen line detections is still limited. Recent JWST/NIRSpec observations have detected Paschen and Brackett hydrogen lines at TWA 27 B (2M1207b). TWA 27 B is classified as a planetary-mass companison (PMC) rather than a planet due to its large mass ratio to the central object ($\approx 5 M_\mathrm{J}$ compared to $25 M_\mathrm{J}$). Nevertheless, TWA 27 B's hydrogen line emission is expected to be same as for planets, given its small mass. We aim to constrain the accretion properties and accretion geometry of TWA 27 B, contributing to our understanding of hydrogen line emission mechanism common to both PMCs and planets. We conduct spectral fitting of four bright hydrogen lines (Pa-$α$, Pa-$β$, Pa-$γ$, Pa-$δ$) with an accretion-shock emission model tailored for forming planets. We estimate the mass accretion rate at $\dot{M} \approx 3 \times 10^{-9}\, M_\mathrm{J}\,\mathrm{yr}^{-1}$ with our fiducial parameters, though this is subject to an uncertainty of up to factor of ten. Our analysis also indicates a dense accretion flow, $n\gtrsim 10^{13}\,\mathrm{cm^{-3}}$ just before the shock, implying a small accretion-shock filling factor $f_\mathrm{f}$ on the planetary surface ($f_\mathrm{f} \lesssim 5\times10^{-4}$). This finding suggests that magnetospheric accretion is occurring at TWA 27 B. Additionally, we carry out a comparative analysis of hydrogen-line emission color to identify the emission mechanism, but the associated uncertainties proved too large for definitive conclusions. This underscores the need for further high-precision observational studies to elucidate these emission mechanisms fully. △ Less

Submitted 24 July, 2024; v1 submitted 22 July, 2024; originally announced July 2024.

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

High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light

Authors: N. Pourré, T. O. Winterhalder, J. -B. Le Bouquin, S. Lacour, A. Bidot, M. Nowak, A. -L. Maire, D. Mouillet, C. Babusiaux, J. Woillez, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube , et al. (151 additional authors not shown)

Abstract: Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working… ▽ More Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 16 pages, 14 figures. Submitted to A&A

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Authors: Asmita Bhandare, Benoît Commerçon, Guillaume Laibe, Mario Flock, Rolf Kuiper, Thomas Henning, Andrea Mignone, Gabriel-Dominique Marleau

Abstract: Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. W… ▽ More Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. We explore if thermal and dynamical conditions developed during disc formation can generate gas flows that efficiently mix and transport well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a 2D axisymmetric geometry and follow the azimuthal component of velocity. Dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10 and 100 micron-sized neutral spherical dust. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 Msun and explore effects of initial rotation rates and cloud core sizes. Our study underlines mechanisms for early transport of dust from inner hot disc regions via the occurrence of meridional flows and outflow. The vortical flow fosters dynamical mixing and retention of dust while thermal pressure driven outflow replenishes dust in the outer disc. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse. [Abridged] △ Less

Submitted 12 May, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: Language edited version

Journal ref: A&A 687, A158 (2024)

Four-of-a-kind? Comprehensive atmospheric characterisation of the HR 8799 planets with VLTI/GRAVITY

Authors: E. Nasedkin, P. Mollière, S. Lacour, M. Nowak, L. Kreidberg, T. Stolker, J. J. Wang, W. O. Balmer, J. Kammerer, J. Shangguan, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli , et al. (73 additional authors not shown)

Abstract: With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels… ▽ More With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5σ$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings. △ Less

Submitted 17 July, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: 45 pages, 25 figures, accepted by Astronomy & Astrophysics

Journal ref: A&A 687, A298 (2024)

Combining Gaia and GRAVITY: Characterising five new Directly Detected Substellar Companions

Authors: T. O. Winterhalder, S. Lacour, A. Mérand, A. -L. Maire, J. Kammerer, T. Stolker, N. Pourré, C. Babusiaux, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (74 additional authors not shown)

Abstract: Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observati… ▽ More Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+. △ Less

Submitted 24 June, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 688, A44 (2024)

MINDS: JWST/NIRCam imaging of the protoplanetary disk PDS 70

Authors: V. Christiaens, M. Samland, Th. Henning, B. Portilla-Revelo, G. Perotti, E. Matthews, O. Absil, L. Decin, I. Kamp, A. Boccaletti, B. Tabone, G. -D. Marleau, E. F. van Dishoeck, M. Güdel, P. -O. Lagage, D. Barrado, A. Caratti o Garatti, A. M. Glauser, G. Olofsson, T. P. Ray, S. Scheithauer, B. Vandenbussche, L. B. F. M. Waters, A. M. Arabhavi, S. L. Grant , et al. (6 additional authors not shown)

Abstract: Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $μ$m and 4.83 $μ$m, assess the presence of Pa-$α$ emission due to accretion onto the protoplan… ▽ More Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths difficult to access from the ground. Aims. We aim to image the circumstellar environment of PDS 70 at 1.87 $μ$m and 4.83 $μ$m, assess the presence of Pa-$α$ emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtain non-coronagraphic JWST/NIRCam images of PDS 70 within the MINDS (MIRI mid-INfrared Disk Survey) program. We leverage the Vortex Image Processing (VIP) package for data reduction, and develop dedicated routines for optimal stellar PSF subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk is used to disentangle the contributions from the disk and the protoplanets. Results. We re-detect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream feeding its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal consistent with a previously proposed protoplanet candidate enshrouded in dust, near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 $μ$m flux of planet b is consistent with atmospheric model predictions, but not that of planet c. We discuss potential origins for this discrepancy, including significant Pa-$α$ line emission. The 4.83 $μ$m fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 6+11 pages, 3+10 figures (text+appendix). Accepted for publication in A&A

Journal ref: A&A 685, L1 (2024)

A catalogue of dual-field interferometric binary calibrators

Authors: M. Nowak, S. Lacour, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, B. Charnay, G. Chauvin, A. Chavez, E. Choquet, V. Christiaens, Y. Clénet, V. Coudé du Foresto, A. Cridland , et al. (75 additional authors not shown)

Abstract: Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to with… ▽ More Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 6 figures. Submitted to A&A

Journal ref: A&A 687, A248 (2024)

Revisiting the Helium and Hydrogen Accretion Indicators at TWA 27B: Weak Mass Flow at Near-Freefall Velocity

Authors: Gabriel-Dominique Marleau, Yuhiko Aoyama, Jun Hashimoto, Yifan Zhou

Abstract: TWA 27B (2M1207b) is the first directly-imaged planetary-mass (MP ~ 5 MJ) companion (Chauvin et al. 2004) and was observed at 0.9--5.3 micron with JWST/NIRSpec (Luhman et al. 2023). To understand the accretion properties of TWA 27B, we search for continuum-subtracted near-infrared helium and hydrogen emission lines and measure their widths and luminosities. We detect the He I triplet at 4.3 sigma… ▽ More TWA 27B (2M1207b) is the first directly-imaged planetary-mass (MP ~ 5 MJ) companion (Chauvin et al. 2004) and was observed at 0.9--5.3 micron with JWST/NIRSpec (Luhman et al. 2023). To understand the accretion properties of TWA 27B, we search for continuum-subtracted near-infrared helium and hydrogen emission lines and measure their widths and luminosities. We detect the He I triplet at 4.3 sigma and all Paschen-series lines covered by NIRSpec (Pa alpha, Pa beta, Pa gamma, Pa delta) at 4--5 sigma. The three brightest Brackett-series lines (Br alpha, Br beta, Br gamma) as well as Pf gamma and Pf delta are tentative detections at 2--3 sigma. We provide upper limits on the other hydrogen lines, including on H alpha through Hubble Space Telescope archival data. Three lines can be reliably deconvolved to reveal an intrinsic width Delta v = 67+-9 km/s, which is 60% of the surface freefall velocity. The line luminosities seem significantly too high to be due to chromospheric activity. Converting line luminosities to an accretion rate yields Mdot ~ 5e-9 MJ/yr when using scaling relationships for planetary masses, and Mdot ~ 0.1e-9 MJ/yr with extrapolated stellar scalings. Several of these lines represent first detections at an accretor of such low mass. The weak accretion rate implies that formation is likely over. This analysis shows that JWST can be used to measure low line-emitting mass accretion rates onto planetary-mass objects, motivates deeper searches for the mass reservoir feeding TWA 27B, and hints that other young directly-imaged objects might -- hitherto unbeknownst -- also be accreting. △ Less

Submitted 19 March, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: In press at ApJ. Main text: 10 pages, 5 figures, 1 table, 1 equation. Appendix: 12 figures (line profiles). Comments welcome. v2: added one reference. v3: added RV velocity (probable systemic velocity) of primary to Fig. 3b. v4: no change, only added DOI of open-access, published version (2024 March 18)

VLTI/GRAVITY Provides Evidence the Young, Substellar Companion HD 136164 Ab formed like a "Failed Star"

Authors: William O. Balmer, L. Pueyo, S. Lacour, J. J. Wang, T. Stolker, J. Kammerer, N. Pourré, M. Nowak, E. Rickman, S. Blunt, A. Sivaramakrishnan, D. Sing, K. Wagner, G. -D. Marleau, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, M. S. Bordoni , et al. (71 additional authors not shown)

Abstract: Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or c… ▽ More Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity

Authors: S. Blunt, W. O. Balmer, J. J. Wang, S. Lacour, S. Petrus, G. Bourdarot, J. Kammerer, N. Pourré, E. Rickman, J. Shangguan, T. Winterhalder, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (73 additional authors not shown)

Abstract: Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging… ▽ More Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426. △ Less

Submitted 6 October, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 28 pages, 14 figures. Revised and resubmitted to AJ

The Comprehensive Archive of Substellar and Planetary Accretion Rates

Authors: S. K. Betti, K. B. Follette, K. Ward-Duong, A. E. Peck, Y. Aoyama, J. Bary, B. Dacus, S. Edwards, G. -D. Marleau, K. Mohamed, J. Palmo, C. Plunkett, C. Robinson, H. Wang

Abstract: Accretion rates ($\dot{M}$) of young stars show a strong correlation with object mass ($M$); however, extension of the $\dot{M}-M$ relation into the substellar regime is less certain. Here, we present the Comprehensive Archive of Substellar and Planetary Accretion Rates (CASPAR), the largest to-date compilation of substellar accretion diagnostics. CASPAR includes: 658 stars, 130 brown dwarfs, and… ▽ More Accretion rates ($\dot{M}$) of young stars show a strong correlation with object mass ($M$); however, extension of the $\dot{M}-M$ relation into the substellar regime is less certain. Here, we present the Comprehensive Archive of Substellar and Planetary Accretion Rates (CASPAR), the largest to-date compilation of substellar accretion diagnostics. CASPAR includes: 658 stars, 130 brown dwarfs, and 10 bound planetary mass companions. In this work, we investigate the contribution of methodological systematics to scatter in the $\dot{M}-M$ relation, and compare brown dwarfs to stars. In our analysis, we rederive all quantities using self-consistent models, distances, and empirical line flux to accretion luminosity scaling relations to reduce methodological systematics. This treatment decreases the original $1σ$ scatter in the $\log \dot{M}-\log M$ relation by $\sim17$%, suggesting that it makes only a small contribution to the dispersion. CASPAR rederived values are best fit by $\dot{M}\propto M^{2.02\pm0.06}$ from 10~$M_\mathrm{J}$ to 2~$M_\odot$, confirming previous results. However, we argue that the brown dwarf and stellar populations are better described separately and by accounting for both mass and age. Therefore, we derive separate age-dependent $\dot{M}-M$ relations for these regions, and find a steepening in the brown dwarf $\dot{M}-M$ slope with age. Within this mass regime, the scatter decreases from 1.36 dex to 0.94 dex, a change of $\sim$44%. This result highlights the significant role that evolution plays in the overall spread of accretion rates, and suggests that brown dwarfs evolve faster than stars, potentially as a result of different accretion mechanisms. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: Accepted for publication in AJ. 38 pages, 21 figures, 11 tables. The database can be found online at Zenodo: https://doi.org/10.5281/zenodo.8393054

VLTI/GRAVITY Observations and Characterization of the Brown Dwarf Companion HD 72946 B

Authors: W. O. Balmer, L. Pueyo, T. Stolker, H. Reggiani, S. Lacour, A. -L. Maire, P. Mollière, M. Nowak, D. Sing, N. Pourré, S. Blunt, J. J. Wang, E. Rickman, Th. Henning, K. Ward-Duong, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet , et al. (74 additional authors not shown)

Abstract: Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. W… ▽ More Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,σ$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis. △ Less

Submitted 15 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

An imaged 15Mjup companion within a hierarchical quadruple system

Authors: A. Chomez, V. Squicciarini, A. -M. Lagrange, P. Delorme, G. Viswanath, M. Janson, O. Flasseur, G. Chauvin, M. Langlois, P. Rubini, S. Bergeon, D. Albert, M. Bonnefoy, S. Desidera, N. Engler, R. Gratton, T. Henning, E. E. Mamajek, G. -D. Marleau, M. R. Meyer, S. Reffert, S. C. Ringqvist, M. Samland

Abstract: Since 2019, the direct imaging B-star Exoplanet Abundance Study (BEAST) at SPHERE@VLT has been scanning the surroundings of young B-type stars in order to ascertain the ultimate frontiers of giant planet formation. Recently, the $17^{+3}_{-4}$ Myr HIP 81208 was found to host a close-in (~50 au) brown dwarf and a wider (~230 au) late M star around the central 2.6Msun primary. Alongside the continua… ▽ More Since 2019, the direct imaging B-star Exoplanet Abundance Study (BEAST) at SPHERE@VLT has been scanning the surroundings of young B-type stars in order to ascertain the ultimate frontiers of giant planet formation. Recently, the $17^{+3}_{-4}$ Myr HIP 81208 was found to host a close-in (~50 au) brown dwarf and a wider (~230 au) late M star around the central 2.6Msun primary. Alongside the continuation of the survey, we are undertaking a complete reanalysis of archival data aimed at improving detection performances so as to uncover additional low-mass companions. We present here a new reduction of the observations of HIP 81208 using PACO ASDI, a recent and powerful algorithm dedicated to processing high-contrast imaging datasets, as well as more classical algorithms and a dedicated PSF-subtraction approach. The combination of different techniques allowed for a reliable extraction of astrometric and photometric parameters. A previously undetected source was recovered at a short separation from the C component of the system. Proper motion analysis provided robust evidence for the gravitational bond of the object to HIP 81208 C. Orbiting C at a distance of ~20 au, this 15Mjup brown dwarf becomes the fourth object of the hierarchical HIP 81208 system. Among the several BEAST stars which are being found to host substellar companions, HIP 81208 stands out as a particularly striking system. As the first stellar binary system with substellar companions around each component ever found by direct imaging, it yields exquisite opportunities for thorough formation and dynamical follow-up studies. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 12 pages, 9 figures, 5 tables. Accepted for publication as a Letter in Astronomy and Astrophysics, section 1. Letters to the Editor

Journal ref: A&A 676, L10 (2023)

Emission line variability of young 10-30 Mjup companions : I. The case of GQ Lup b and GSC 06214-00210 b

Authors: Dorian Demars, Mickael Bonnefoy, Catherine Dougados, Yuhiko Aoyama, Thanawuth Thanathibodee, Gabriel-Dominique Marleau, Pascal Tremblin, Philippe Delorme, Paulina Palma-Bifani, Simon Petrus, Brendan P. Bowler, Gael Chauvin, Anne-Marie Lagrange

Abstract: Emission lines indicative of active accretion have been seen on a handful of low-mass companions (M < 30 MJup) to stars. Line variability is ubiquitous on stellar accretors but has never been characterized in detail on low-mass companions and can give insights on the accretion mechanism at play. We investigate the emission line variability of two low-mass companions (M<30 MJup) to stars to underst… ▽ More Emission lines indicative of active accretion have been seen on a handful of low-mass companions (M < 30 MJup) to stars. Line variability is ubiquitous on stellar accretors but has never been characterized in detail on low-mass companions and can give insights on the accretion mechanism at play. We investigate the emission line variability of two low-mass companions (M<30 MJup) to stars to understand their accretion mechanisms. Using J-band observations, we analyze the short to long-term variability of the HI Paschen β emission line (1.282 μm) for GQ Lup b and GSC 06214-00210 b. Archival spectroscopic observations are also examined to extend the time span. We compare their line profiles and intensities to more massive accretors and magnetospheric accretion and shock models. Both objects have HI Paschen β flux variability that is moderate at short timescales (< 50 %) and increases at longer timescales (~1000 % on decade timescales), resembling classical T Tauri stars. GQ Lup b's line profiles are compatible with magnetospheric accretion. GSC 06214-00210 b's profiles are reproduced by both magnetospheric accretion and shock models, except for the brightest epoch for which the shock model is highly favored. Both companions have C/O values broadly consistent with solar values. While magnetospheric accretion is favored for GQ Lup b, higher resolution (R > 10000) observations are required to disentangle the two (non-exclusive) line formation mechanisms. The similarity in variability behavior may support similar accretion mechanisms between these low-mass companions and classical T Tauri stars. The significant variability observed at months and longer timescales could explain the low yield of Hα imaging campaigns. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 21 pages, 10 figures

Journal ref: A&A 676, A123 (2023)

The Planetary Accretion Shock. III. Smoothing-free 2.5D simulations and calculation of H alpha emission

Authors: Gabriel-Dominique Marleau, Rolf Kuiper, William Béthune, Christoph Mordasini

Abstract: Surveys have looked for H alpha emission from accreting gas giants but found very few objects. Analyses of the detections and non-detections have assumed that the entire gas flow feeding the planet is in radial free-fall. However, hydrodynamical simulations suggest that this is far from reality. We calculate the H alpha emission from multidimensional accretion onto a gas giant, following the gas f… ▽ More Surveys have looked for H alpha emission from accreting gas giants but found very few objects. Analyses of the detections and non-detections have assumed that the entire gas flow feeding the planet is in radial free-fall. However, hydrodynamical simulations suggest that this is far from reality. We calculate the H alpha emission from multidimensional accretion onto a gas giant, following the gas flow from Hill-sphere scales down to the circumplanetary disc (CPD) and the planetary surface. We perform azimuthally-symmetric radiation-hydrodynamics simulations around the planet and use modern tabulated gas and dust opacities. Crucially, contrasting with most previous simulations, we do not smooth the gravitational potential and do follow the flow down to the planetary surface, where grid cells are 0.01 Jupiter radii small radially. We find that only roughly one percent of the net gas inflow into the Hill sphere reaches directly the planet. As expected for ballistic infall trajectories, most of the gas falls at too large a distance on the CPD to generate H alpha. Including radiation transport removes the high-velocity sub-surface flow previously seen in hydrodynamics-only simulations, so that only the free planet surface and the inner regions of the CPD emit substantially H alpha. Unless magnetospheric accretion, which we neglect here, additionally produces H alpha, the corresponding H alpha production efficiency is much smaller than usually assumed, which needs to be taken into account when analysing (non-)detection statistics. △ Less

Submitted 19 July, 2023; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: Main text: 16 pages, 7 figures; appendices: 4 pages, 4 figures; references: 1.7 pages. Printed at ApJ on 2023 July 20th. Version 4: This information and correcting R_Hill -> R_Jup (obvious) in the axis labels of Fig. 9

Journal ref: The Astrophysical Journal, Volume 952, Issue 1, article id. 89, 18 pp. (2023)

Planetary line-to-accretion luminosity scaling relations: Extrapolating to higher-order hydrogen lines

Authors: Gabriel-Dominique Marleau, Yuhiko Aoyama

Abstract: Aoyama et al. (2021, ApJL) provided scaling relations between hydrogen-line luminosities and the accretion luminosity for planetary-mass objects. These fits should be an improvement over blind extrapolations of stellar relations. The fits go up only to the n = 8 electron energy level, but higher-n Balmer lines have been observed in the near-UV at Delorme 1 (AB)b with UVES (Ringqvist et al. 2023).… ▽ More Aoyama et al. (2021, ApJL) provided scaling relations between hydrogen-line luminosities and the accretion luminosity for planetary-mass objects. These fits should be an improvement over blind extrapolations of stellar relations. The fits go up only to the n = 8 electron energy level, but higher-n Balmer lines have been observed in the near-UV at Delorme 1 (AB)b with UVES (Ringqvist et al. 2023). We extend the scaling relations to higher-n levels for the Balmer and other series by fitting the fit coefficients (a, b) themselves and extrapolating them. Within the assumption of an accretion shock as the source of line emission, these fits should be robust for accreting planetary-mass objects. △ Less

Submitted 2 March, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 3 pages, 1 four-panel figure. Published on December 9th, 2022. v1: updated links to published papers, added thanks. Fits of scaling relations available in Tomas Stolker's "species" toolkit (https://species.readthedocs.io). v2: corrected the accidental shuffle of the author names (!)

Journal ref: RNAAS (2022) 6, 262

Forbidden emission lines in protostellar outflows and jets with MUSE

Authors: Lizxandra Flores-Rivera, Mario Flock, Nicolás Kurtovic, Bernd Husemann, Andrea Banzatti, Simon C. Ringqvist, Sebastian Kamann, André Müller, Christian Fendt, Rebeca Garcia Lopez, Gabriel-Dominique Marleau, Thomas Henning, Carlos Carrasco-Gonzalez, Roy van Boekel, Miriam Keppler, Ralf Launhardt, Yuhiko Aoyama

Abstract: Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral fiel… ▽ More Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral field spectrograph of the Multi Unit Spectroscopic Explorer (MUSE), at the Very Large Telescope (VLT). With a resolution of 0.025 X 0.025 arcsec$^{2}$, we aim to derive the position angle of the outflow/jet (PA$_{outflow/jet}$) that is connected with the inner disk. The forbidden emission lines analyzed here have their origin at the inner parts of the protoplanetary disk. From the maximum intensity emission along the outflow/jet in DL Tau, CI Tau, DS Tau, IP Tau, and IM Lup, we were able to reliably measure the PA$_{outflow/jet}$ for most of the identified lines. We found that our estimates agree with PA$_{dust}$ for most of the disks. These estimates depend on the signal-to-noise level and the collimation of the outflow (jet). The outflows/jets in CIDA 9, GO Tau, and GW Lup are too compact for a PA$_{outflow/jet}$ to be estimated. Based on our kinematics analysis, we confirm that DL Tau and CI Tau host a strong outflow/jet with line-of-sight velocities much greater than 100 km s$^{-1}$, whereas DS Tau, IP Tau, and IM Lup velocities are lower and their structures encompass low-velocity components to be more associated with winds. Our estimates for the mass-loss rate, $\dot{M}_{loss}$, range between (1.1-6.5)X10$^{-7}$-10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for the disk-outflow/jet systems analyzed here. The outflow/jet systems analyzed here are aligned within around 1 degree between the inner and outer disk. Further observations are needed to confirm a potential misalignment in IM Lup. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted on December 19 2022

Resolved near-UV hydrogen emission lines at 40-Myr super-Jovian protoplanet Delorme 1 (AB)b: Indications of magnetospheric accretion

Authors: Simon C. Ringqvist, Gayathri Viswanath, Yuhiko Aoyama, Markus Janson, Gabriel-Dominique Marleau, Alexis Brandeker

Abstract: We have followed up on our observations of the ~ 40-Myr, and still accreting, PMC Delorme 1 (AB)b. We used high-resolution spectroscopy to characterise the accretion process further by accessing the wealth of emission lines in the near-UV. With VLT/UVES, we obtained R ~ 50000 spectroscopy at 330--452 nm. After separating the emission of the companion from that of the M5 low-mass binary, we perform… ▽ More We have followed up on our observations of the ~ 40-Myr, and still accreting, PMC Delorme 1 (AB)b. We used high-resolution spectroscopy to characterise the accretion process further by accessing the wealth of emission lines in the near-UV. With VLT/UVES, we obtained R ~ 50000 spectroscopy at 330--452 nm. After separating the emission of the companion from that of the M5 low-mass binary, we performed a detailed emission-line analysis, which included planetary accretion shock modelling. We reaffirm ongoing accretion in Delorme 1 (AB)b and report the first detections in a (super-Jovian) protoplanet of resolved hydrogen line emission in the near-UV (H-gamma, H-delta, H-epsilon, H8 and H9). We tentatively detect H11, H12, He I and Ca II H/K. The analysis strongly favours a planetary accretion shock with a line-luminosity-based accretion rate dMp/dt = 2e-8 MJ/yr. The lines are asymmetric and well described by sums of narrow and broad components with different velocity shifts. Overall line shapes are best explained by a pre-shock velocity v0 = 170+-30 km/s, implying a planetary mass Mp = 13+-5 MJ, and number densities n0 ~ 1e13/cc or n0 ~ 1e11/cc. The higher density implies a small line-emitting area of ~ 1% relative to the planetary surface. This favours magnetospheric accretion, a case potentially strengthened by the presence of blueshifted emission in the asymmetrical profiles.High-resolution spectroscopy offers the opportunity to resolve line profiles, crucial for studying the accretion process in depth. The super-Jovian protoplanet Delorme 1 (AB)b is still accreting at ~ 40 Myr. Thus, Delorme 1 belongs to the growing family of Peter Pan disc systems with protoplanetary and/or circumplanetary disc(s) far beyond the typically assumed disc lifetimes. Further observations of this benchmark companion, and its presumed disc(s), will help answer key questions about the accretion geometry in PMCs. △ Less

Submitted 22 March, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: Published in A&A 669, L12, 11 pages, abbreviated abstract

Journal ref: A&A 669, L12 (2023)

Radial Velocity Survey for Planets around Young stars (RVSPY) A transiting warm super-Jovian planet around HD 114082, a young star with a debris disk

Authors: O. Zakhozhay, R. Launhardt, T. Trifonov, M. Kürster, S. Reffert, Th. Henning, R. Brahm, J. Vinés, G. -D. Marleau, J. Patel

Abstract: Aiming to detect planetary companions to young stars with debris disks via the radial velocity method, we observed HD114082 during April 2018 - August 2022 as one of the targets of our RVSPY program (Radial Velocity Survey for Planets around Young stars). We used the FEROS spectrograph, mounted to the MPG/ESO 2.2 m telescope in Chile, to obtain high signal-to-noise spectra and time series of preci… ▽ More Aiming to detect planetary companions to young stars with debris disks via the radial velocity method, we observed HD114082 during April 2018 - August 2022 as one of the targets of our RVSPY program (Radial Velocity Survey for Planets around Young stars). We used the FEROS spectrograph, mounted to the MPG/ESO 2.2 m telescope in Chile, to obtain high signal-to-noise spectra and time series of precise radial velocities (RVs). Additionally, we analyzed archival HARPS spectra and TESS photometric data. We used the CERES, CERES++ and SERVAL pipelines to derive RVs and activity indicators and ExoStriker for the independent and combined analysis of the RVs and TESS photometry. We report the discovery of a warm super-Jovian companion around HD114082 based on a 109.8$\pm$0.4 day signal in the combined RV data from FEROS and HARPS, and on one transit event in the TESS photometry. The best-fit model indicates a 8.0$\pm$1.0 Mjup companion with a radius of 1.00$\pm$0.03 Rjup in an orbit with a semi-major axis of 0.51$\pm$0.01 au and an eccentricity of 0.4$\pm$0.04. The companions orbit is in agreement with the known near edge-on debris disk located at about 28 au. HD114082b is possibly the youngest (15$\pm$6 Myr), and one of only three younger than 100 Myr giant planetary companions for which both their mass and radius have been determined observationally. It is probably the first properly model-constraining giant planet that allows distinguishing between hot and cold-start models. It is significantly more compatible with the cold-start model. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 10 pages, 9 figures, 5 tables; Accepted for publication in A&A Letters

Journal ref: A&A 667, L14 (2022)

Near-infrared Accretion Signatures from the Circumbinary Planetary Mass Companion Delorme 1 (AB)b

Authors: S. K. Betti, K. B. Follette, K. Ward-Duong, Y. Aoyama, G. -D. Marleau, J. Bary, C. Robinson, M. Janson, W. Balmer, G. Chauvin, P. Palma-Bifani

Abstract: Accretion signatures from bound brown dwarf and protoplanetary companions provide evidence for ongoing planet formation, and accreting substellar objects have enabled new avenues to study the astrophysical mechanisms controlling formation and accretion processes. Delorme 1 (AB)b, a ~30-45 Myr circumbinary planetary mass companion, was recently discovered to exhibit strong H$α$ emission. This sugge… ▽ More Accretion signatures from bound brown dwarf and protoplanetary companions provide evidence for ongoing planet formation, and accreting substellar objects have enabled new avenues to study the astrophysical mechanisms controlling formation and accretion processes. Delorme 1 (AB)b, a ~30-45 Myr circumbinary planetary mass companion, was recently discovered to exhibit strong H$α$ emission. This suggests ongoing accretion from a circumplanetary disk, somewhat surprising given canonical gas disk dispersal timescales of 5-10 Myr. Here, we present the first NIR detection of accretion from the companion in Pa$β$, Pa$γ$, and Br$γ$ emission lines from SOAR/TripleSpec 4.1, confirming and further informing its accreting nature. The companion shows strong line emission, with $L_{line} \approx 1-6 \times 10^{-8}~L_\odot$ across lines and epochs, while the binary host system shows no NIR hydrogen line emission ($L_{line} <0.32-11\times10^{-7}\ L_\odot$). Observed NIR hydrogen line ratios are more consistent with a planetary accretion shock than with local line excitation models commonly used to interpret stellar magnetospheric accretion. Using planetary accretion shock models, we derive mass accretion rate estimates of $\dot{M}_{\mathrm{pla}}\sim3$-$4\times 10^{-8}\ M_\mathrm{J}$ yr$^{-1}$, somewhat higher than expected under the standard star formation paradigm. Delorme 1 (AB)b's high accretion rate is perhaps more consistent with formation via disk fragmentation. Delorme 1 (AB)b is the first protoplanet candidate with clear (S/N$\sim$5) NIR hydrogen line emission. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJL

Direct discovery of the inner exoplanet in the HD206893 system. Evidence for deuterium burning in a planetary-mass companion

Authors: S. Hinkley, S. Lacour, G. -D. Marleau, A. M. Lagrange, J. J. Wang, J. Kammerer, A. Cumming, M. Nowak, L. Rodet, T. Stolker, W. -O. Balmer, S. Ray, M. Bonnefoy, P. Mollière, C. Lazzoni, G. Kennedy, C. Mordasini, R. Abuter, S. Aigrain, A. Amorim, R. Asensio-Torres, C. Babusiaux, M. Benisty, J. -P. Berger, H. Beust , et al. (89 additional authors not shown)

Abstract: Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoi… ▽ More Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $μ$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au. △ Less

Submitted 3 April, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to A&A

X-SHYNE: X-shooter spectra of young exoplanet analogs. I. A medium-resolution 0.65-2.5$\mathrm{μm}$ one-shot spectrum of VHS\,1256-1257 b

Authors: Simon Petrus, Gaël Chauvin, Mickaël Bonnefoy, Pascal Tremblin, Benjamin Charnay, Philippe Delorme, Gabriel-Dominique Marleau, Amelia Bayo, Elena Manjavacas, Anne-Marie Lagrange, Paul Mollière, Paulina Palma-Bifani, Beth Biller James-S. Jenkins

Abstract: We present simultaneous 0.65-2.5 microns medium resolution (3300 < R < 8100) VLT/X-Shooter spectra of the young low-mass (19+/-5MJup) L-T transition object VHS 1256-1257 b, a known spectroscopic analogue of HR8799d. The companion is a prime target for the JWST Early Release Science (ERS) and one of the highest-amplitude variable brown-dwarf known to date. We compare the spectrum to the custom grid… ▽ More We present simultaneous 0.65-2.5 microns medium resolution (3300 < R < 8100) VLT/X-Shooter spectra of the young low-mass (19+/-5MJup) L-T transition object VHS 1256-1257 b, a known spectroscopic analogue of HR8799d. The companion is a prime target for the JWST Early Release Science (ERS) and one of the highest-amplitude variable brown-dwarf known to date. We compare the spectrum to the custom grids of cloudless ATMO models exploring different atmospheric composition with the Bayesian inference tool ForMoSA. We also re-analyze low-resolution HST/WFC3 1.10-1.67 microns spectra at minimum and maximum variability to contextualize the X-Shooter data interpretation. The models reproduce the slope and most molecular absorption from 1.10 to 2.48 microns self-consistently but fail to provide a radius consistent with evolutionary model predictions. They do not reproduce consistently the optical spectrum and the depth of the K I doublets in the J-band. We derive Teff = 1380+/-54 K, log(g) = 3.97+/-0.48 dex, [M/H] = 0.21+/-0.29, and C/O > 0.63. Our inversion of the HST/WFC3 spectra suggests a relative change of 27+6-5 K of the disk-integrated Teff correlated with the near-infrared brightness. Our data anchor the characterization of that object in the near-infrared and could be used jointly to the ERS mid-infrared data to provide the most detailed characterization of an ultracool dwarf to date. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 9 pages, 4 figures, 4 tables

Journal ref: A&A 670, L9 (2023)

A scaled-up planetary system around a supernova progenitor

Authors: V. Squicciarini, R. Gratton, M. Janson, E. E. Mamajek, G. Chauvin, P. Delorme, M. Langlois, A. Vigan, S. C. Ringqvist, G. Meeus, S. Reffert, M. Kenworthy, M. R. Meyer, M. Bonnefoy, M. Bonavita, D. Mesa, M. Samland, S. Desidera, V. D'Orazi, N. Engler, E. Alecian, A. Miglio, T. Henning, S. P. Quanz, L. Mayer , et al. (2 additional authors not shown)

Abstract: Virtually all known exoplanets reside around stars with $M<2.3~M_\odot$; to clarify if the dearth of planets around more massive stars is real, we launched the direct-imaging B-star Exoplanet Abundance STudy (BEAST) survey targeting B stars ($M>2.4~M_\odot$) in the young (5-20 Myr) Scorpius-Centaurus association (Sco-Cen). Here we present the case of a massive ($M \sim 9~M_\odot$) BEAST target,… ▽ More Virtually all known exoplanets reside around stars with $M<2.3~M_\odot$; to clarify if the dearth of planets around more massive stars is real, we launched the direct-imaging B-star Exoplanet Abundance STudy (BEAST) survey targeting B stars ($M>2.4~M_\odot$) in the young (5-20 Myr) Scorpius-Centaurus association (Sco-Cen). Here we present the case of a massive ($M \sim 9~M_\odot$) BEAST target, $μ^2$ Sco. Based on kinematic information, we found that $μ^2$ Sco is a member of a small group which we label Eastern Lower Scorpius, refining in turn the precision on stellar parameters. Around this star we identified a robustly detected substellar companion ($14.4\pm 0.8 M_J$) at a projected separation of $290\pm 10$ au, and a probable second object ($18.5\pm 1.5 M_J$) at $21\pm 1$ au. The planet-to-star mass ratios of these objects are similar to that of Jupiter to the Sun, and their irradiation is similar to those of Jupiter and Mercury, respectively. The two companions of $μ^2$ Sco are naturally added to the giant planet b Cen b recently discovered by BEAST; although slightly more massive than the deuterium burning limit, their properties resemble those of giant planets around less massive stars and they are better reproduced by a formation under a planet-like, rather than a star-like scenario. Irrespective of the (needed) confirmation of the inner companion, $μ^2$ Sco is the first star that would end its life as a supernova that hosts such a system. The tentative high frequency of BEAST discoveries shows that giant planets or small-mass brown dwarfs can form around B stars. When putting this finding in the context of core accretion and gravitational instability, we conclude that the current modeling of both mechanisms is not able to produce this kind of companion. BEAST will pave the way for the first time to an extension of these models to intermediate and massive stars. (abridged) △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 26 pages, 19 figures. Accepted for publication in A&A

Journal ref: A&A 664, A9 (2022)

A wide-orbit giant planet in the high-mass b Centauri binary system

Authors: Markus Janson, Raffaele Gratton, Laetitia Rodet, Mickael Bonnefoy, Philippe Delorme, Eric E. Mamajek, Sabine Reffert, Lukas Stock, Gabriel-Dominique Marleau, Maud Langlois, Gael Chauvin, Silvano Desidera, Simon Ringqvist, Lucio Mayer, Gayathri Viswanath, Vito Squicciarini, Michael R. Meyer, Matthias Samland, Simon Petrus, Ravit Helled, Matthew A. Kenworthy, Sascha P. Quanz, Beth Biller, Thomas Henning, Dino Mesa , et al. (2 additional authors not shown)

Abstract: Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until… ▽ More Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until a turnover point at 1.9 solar masses, above which the frequency rapidly decreases. This could potentially imply that planet formation is impeded around more massive stars, and that giant planets around stars exceeding 3 solar masses may be rare or non-existent. However, the methods used to detect planets in small orbits are insensitive to planets in wide orbits. Here we demonstrate the existence of a planet at 560 times the Sun-Earth distance from the 6-10 solar mass binary b Centauri through direct imaging. The planet-to-star mass ratio of 0.10-0.17% is similar to the Jupiter-Sun ratio, but the separation of the detected planet is ~100 times wider than that of Jupiter. Our results show that planets can reside in much more massive stellar systems than what would be expected from extrapolation of previous results. The planet is unlikely to have formed in-situ through the conventional core accretion mechanism, but might have formed elsewhere and arrived to its present location through dynamical interactions, or might have formed via gravitational instability. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: Manuscript version. Published in Nature 9 December 2021

Accreting protoplanets: Spectral signatures and magnitude of gas and dust extinction at H alpha

Authors: G. -D. Marleau, Y. Aoyama, R. Kuiper, K. Follette, N. J. Turner, G. Cugno, C. F. Manara, S. Y. Haffert, D. Kitzmann, S. C. Ringqvist, K. R. Wagner, R. van Boekel, S. Sallum, M. Janson. T. O. B. Schmidt, L. Venuti, Ch. Lovis, C. Mordasini

Abstract: Accreting planets have been seen at Ha (H alpha), but targeted searches have not been fruitful. For planets, accretion tracers should come from the shock itself, exposing them to extinction by the accreting material. High-resolution (R>5e4) spectrographs at Ha should soon allow studying how the incoming material shapes the line profile. We calculate how much the gas and dust accreting onto a plane… ▽ More Accreting planets have been seen at Ha (H alpha), but targeted searches have not been fruitful. For planets, accretion tracers should come from the shock itself, exposing them to extinction by the accreting material. High-resolution (R>5e4) spectrographs at Ha should soon allow studying how the incoming material shapes the line profile. We calculate how much the gas and dust accreting onto a planet reduce the Ha flux from the shock at the planetary surface and how they affect line shapes. We also study the absorption-modified relationship between Ha luminosity and Mdot. We compute the high-resolution radiative transfer of the Ha line using a 1D velocity-density-temperature structure for the inflowing matter in three representative accretion geometries: spherical symmetry, polar inflow, and magnetospheric accretion. For each, we explore wide ranges of Mdot and planet mass M. We use detailed gas opacities and estimate dust opacities. At Mdot<3e-6 MJ/yr, gas extinction is negligible for spherical or polar inflow and at most A_Ha<0.5 mag for magnetospheric accretion. Up to Mdot~3e-4 MJ/yr, the gas has A_Ha<4 mag. This decreases with M. We estimate realistic dust opacities at Ha as ~0.01-10 cm^2/g, i.e., 10-1e4 times lower than in the ISM. Extinction flattens the L_Ha-Mdot relationship, which becomes non-monotonic with a maximum L_Ha~1e-4 LSun near Mdot~1e-4 MJ/yr for M~10 MJ. In magnetospheric accretion, the gas can introduce features in line profiles, but the velocity gradient smears them out in other geometries. For most of parameter space, extinction by the accreting matter should be negligible, simplifying interpretation of observations, especially for planets in gaps. At high Mdot, strong absorption reduces the Ha flux, and some measurements can be interpreted as two Mdot values. Line profiles at R~1e5 can provide complex constraints on the accretion flow's thermal-dynamical structure. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 25 pages, 11 figures (main text) + 3 pages, 4 figures (three appendices). In press at A&A. Comments welcome

Journal ref: A&A 657, A38 (2022)

Characterizing the protolunar disk of the accreting companion GQ Lupi B

Authors: Tomas Stolker, Sebastiaan Y. Haffert, Aurora Y. Kesseli, Rob G. van Holstein, Yuhiko Aoyama, Jarle Brinchmann, Gabriele Cugno, Julien H. Girard, Gabriel-Dominique Marleau, Gabriele Cugno, Michael R. Meyer, Julien Milli, Sascha P. Quanz, Ignas A. G. Snellen, Kamen O. Todorov

Abstract: GQ Lup B is a young and accreting, substellar companion that appears to drive a spiral arm in the circumstellar disk of its host star. We report high-contrast imaging observations of GQ Lup B with VLT/NACO at 4-5 $μ$m and medium-resolution integral field spectroscopy with VLT/MUSE. The optical spectrum is consistent with an M9 spectral type, shows characteristics of a low-gravity atmosphere, and e… ▽ More GQ Lup B is a young and accreting, substellar companion that appears to drive a spiral arm in the circumstellar disk of its host star. We report high-contrast imaging observations of GQ Lup B with VLT/NACO at 4-5 $μ$m and medium-resolution integral field spectroscopy with VLT/MUSE. The optical spectrum is consistent with an M9 spectral type, shows characteristics of a low-gravity atmosphere, and exhibits strong H$α$ emission. The $H-M'$ color is $\gtrsim$1 mag redder than field dwarfs with similar spectral types and a detailed analysis of the spectral energy distribution (SED) from optical to mid-infrared wavelengths reveals excess emission in the $L'$, NB4.05, and $M'$ bands. The excess flux is well described by a blackbody component with $T_\mathrm{disk} \approx 460$ K and $R_\mathrm{disk} \approx 65\,R_\mathrm{J}$ and is expected to trace continuum emission from small grains in a protolunar disk. We derive an extinction of $A_V \approx 2.3$ mag from the broadband SED with a suspected origin in the vicinity of the companion. We also combine 15 yr of astrometric measurements and constrain the mutual inclination with the circumstellar disk to $84 \pm 9$ deg, indicating a tumultuous dynamical evolution or a stellar-like formation pathway. From the measured H$α$ flux and the estimated companion mass, $M_\mathrm{p} \approx 30\,M_\mathrm{J}$, we derive an accretion rate of $\dot{M} \approx 10^{-6.5}\,M_\mathrm{J}\,\mathrm{yr}^{-1}$. We speculate that the disk is in a transitional stage in which the assembly of satellites from a pebble reservoir has opened a central cavity while GQ Lup B is in the final stages of its formation. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: 28 pages, 11 figures, accepted for publication in AJ

Comparison of planetary Hα-emission models: A new correlation with accretion luminosity

Authors: Yuhiko Aoyama, Gabriel-Dominique Marleau, Masahiro Ikoma, Christoph Mordasini

Abstract: Accreting planets have been detected through their hydrogen-line emission, specifically H$α$. To interpret this, stellar-regime empirical correlations between the H$α$ luminosity $L_\mathrm{Hα}$ and the accretion luminosity $L_\mathrm{acc}$ or accretion rate $\dot{M}$ have been extrapolated to planetary masses, however without validation. We present a theoretical $L_\mathrm{acc}$--$L_\mathrm{Hα}$… ▽ More Accreting planets have been detected through their hydrogen-line emission, specifically H$α$. To interpret this, stellar-regime empirical correlations between the H$α$ luminosity $L_\mathrm{Hα}$ and the accretion luminosity $L_\mathrm{acc}$ or accretion rate $\dot{M}$ have been extrapolated to planetary masses, however without validation. We present a theoretical $L_\mathrm{acc}$--$L_\mathrm{Hα}$ relationship applicable to a shock at the surface of a planet. We consider wide ranges of accretion rates and masses and use detailed spectrally-resolved, non-equilibrium models of the postshock cooling. The new relationship gives a markedly higher $L_\mathrm{acc}$ for a given $L_\mathrm{Hα}$ than fits to young stellar objects, because Ly-$α$, which is not observable, carries a large fraction of $L_\mathrm{acc}$. Specifically, an $L_\mathrm{Hα}$ measurement needs ten to 100 times higher $L_\mathrm{acc}$ and $\dot{M}$ than previously predicted, which may explain the rarity of planetary H$α$ detections. We also compare the $\dot{M}$--$L_\mathrm{Hα}$ relationships coming from the planet-surface shock or implied by accretion-funnel emission. Both can contribute simultaneously to an observed H$α$ signal but at low (high) $\dot{M}$ the planetary-surface shock (heated funnel) dominates. Only the shock produces Gaussian line wings. Finally, we discuss accretion contexts in which different emission scenarios may apply, putting recent literature models in perspective, and also present $L_\mathrm{acc}$--$L_\mathrm{line}$ relationships for several other hydrogen lines. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 13 pages, accepted to ApJL

Journal ref: The Astrophysical Journal Letters, Volume 917, Issue 2, article id. L30 (2021)

Molecular mapping of the PDS70 system: No molecular absorption signatures from the forming planet PDS70 b

Authors: G. Cugno, P. Patapis, T. Stolker, S. P. Quanz, A. Boehle, H. J. Hoeijmakers, G. -D. Marleau, P. Mollière, E. Nasedkin, I. A. G. Snellen

Abstract: Aims. We aim to detect molecules in the atmosphere of the young forming companion PDS70 b by searching for atmospheric absorption features typical of substellar objects. Methods. We obtained medium-resolution (R$\approx$5075) spectra of the PDS70 planetary system with the SINFONI integral field spectrograph at the Very Large Telescope. We applied molecular mapping, based on cross-correlation with… ▽ More Aims. We aim to detect molecules in the atmosphere of the young forming companion PDS70 b by searching for atmospheric absorption features typical of substellar objects. Methods. We obtained medium-resolution (R$\approx$5075) spectra of the PDS70 planetary system with the SINFONI integral field spectrograph at the Very Large Telescope. We applied molecular mapping, based on cross-correlation with synthetic spectra, to identify signatures of molecular species in the atmosphere of the planet. Results. Although the planet emission is clearly detected when resampling the data to lower resolution, no molecular species could be identified with the cross-correlation technique. We estimated upper limits on the abundances of H$_2$O, CO and CH$_4$ ($\log(X_\mathrm{mol}) < -4.0$, $-4.1$ and $-4.9$, respectively) assuming a clear atmosphere, and we explored the impact of clouds, which increase the upper limits by a factor up to 0.7 dex. Assuming that the observations directly probe the planet's atmosphere, we found a lack of molecular species compared to other directly imaged companions or field objects. Under the assumption that the planet atmosphere presents similar characteristics to other directly imaged planets, we conclude that a dusty environment surrounds the planet, effectively obscuring any feature generated in its atmosphere. We quantify the extinction necessary to impede the detection ($A_V\approx16-17$ mag), pointing to the possibility of higher optical thickness than previously estimated from other studies. Finally, the non-detection of molecular species conflicts with atmospheric models previously proposed to describe the forming planet. Conclusions. To unveil how giant planets form, a comprehensive approach that includes constraints from multiple techniques needs to be undertaken. Molecular mapping emerges as an alternative to more classical techniques like SED fitting. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A (16 pages, 10 figures)

Journal ref: A&A 653, A12 (2021)

The First Dynamical Mass Measurement in the HR 8799 System

Authors: G. Mirek Brandt, Timothy D. Brandt, Trent J. Dupuy, Daniel Michalik, Gabriel-Dominique Marleau

Abstract: HR 8799 hosts four directly imaged giant planets, but none has a mass measured from first principles. We present the first dynamical mass measurement in this planetary system, finding that the innermost planet HR~8799~e has a mass of $9.6^{+1.9}_{-1.8} \, M_{\rm Jup}$. This mass results from combining the well-characterized orbits of all four planets with a new astrometric acceleration detection (… ▽ More HR 8799 hosts four directly imaged giant planets, but none has a mass measured from first principles. We present the first dynamical mass measurement in this planetary system, finding that the innermost planet HR~8799~e has a mass of $9.6^{+1.9}_{-1.8} \, M_{\rm Jup}$. This mass results from combining the well-characterized orbits of all four planets with a new astrometric acceleration detection (5$σ$) from the Gaia EDR3 version of the Hipparcos-Gaia Catalog of Accelerations. We find with 95\% confidence that HR~8799~e is below $13\, M_{\rm Jup}$, the deuterium-fusing mass limit. We derive a hot-start cooling age of $42^{+24}_{-16}$\,Myr for HR~8799~e that agrees well with its hypothesized membership in the Columba association but is also consistent with an alternative suggested membership in the $β$~Pictoris moving group. We exclude the presence of any additional $\gtrsim$5-$M_{\rm Jup}$ planets interior to HR~8799~e with semi-major axes between $\approx$3-16\,au. We provide proper motion anomalies and a matrix equation to solve for the mass of any of the planets of HR~8799 using only mass ratios between the planets. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ Letters

Perturbers: SPHERE detection limits to planetary-mass companions in protoplanetary disks

Authors: R. Asensio-Torres, Th. Henning, F. Cantalloube, P. Pinilla, D. Mesa, A. Garufi, S. Jorquera, R. Gratton, G. Chauvin, J. Szulagyi, R. van Boekel, R. Dong, G. -D. Marleau, M. Benisty, M. Villenave, C. Bergez-Casalou, C. Desgrange, M. Janson, M. Keppler, M. Langlois, F. Menard, E. Rickman, T. Stolker, M. Feldt, T. Fusco , et al. (3 additional authors not shown)

Abstract: The detection of a wide range of substructures such as rings, cavities and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in… ▽ More The detection of a wide range of substructures such as rings, cavities and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in the disk (perturbers), but so far the only bona-fide detection has been the two giant planets around PDS 70. Here, we collect a sample of 15 protoplanetary disks showing substructures in SPHERE scattered light images and present a homogeneous derivation of planet detection limits in these systems. We also estimate the mass of these perturbers through a Hill radius prescription and a comparison to ALMA data. Assuming that one single planet carves each substructure in scattered light, we find that more massive perturbers are needed to create gaps within cavities than rings, and that we might be close to a detection in the cavities of RX J1604, RX J1615, Sz Cha, HD 135344B and HD 34282. We reach typical mass limits in these cavities of 3-10 Mjup. For planets in the gaps between rings, we find that the detection limits of SPHERE are about an order of magnitude away in mass, and that the gaps of PDS 66 and HD 97048 seem to be the most promising structures for planet searches. The proposed presence of massive planets causing spiral features in HD 135344B and HD 36112 are also within SPHERE's reach assuming hot-start models.These results suggest that current detection limits are able to detect hot-start planets in cavities, under the assumption that they are formed by a single perturber located at the centre of the cavity. More realistic planet mass constraints would help to clarify whether this is actually the case, which might point to perturbers not being the only way of creating substructures. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 26 pages, 33 figures, accepted for publication in A&A

Journal ref: A&A 652, A101 (2021)

A faint companion around CrA-9: protoplanet or obscured binary?

Authors: V. Christiaens, M. -G. Ubeira-Gabellini, H. Cánovas, P. Delorme, B. Pairet, O. Absil, S. Casassus, J. H. Girard, A. Zurlo, Y. Aoyama, G-D. Marleau, L. Spina, N. van der Marel, L. Cieza, G. Lodato, S. Pérez, C. Pinte, D. J. Price, M. Reggiani

Abstract: Understanding how giant planets form requires observational input from directly imaged protoplanets. We used VLT/NACO and VLT/SPHERE to search for companions in the transition disc of 2MASS J19005804-3645048 (hereafter CrA-9), an accreting M0.75 dwarf with an estimated age of 1-2 Myr. We found a faint point source at $\sim$0.7'' separation from CrA-9 ($\sim$108 au projected separation). Our 3-epoc… ▽ More Understanding how giant planets form requires observational input from directly imaged protoplanets. We used VLT/NACO and VLT/SPHERE to search for companions in the transition disc of 2MASS J19005804-3645048 (hereafter CrA-9), an accreting M0.75 dwarf with an estimated age of 1-2 Myr. We found a faint point source at $\sim$0.7'' separation from CrA-9 ($\sim$108 au projected separation). Our 3-epoch astrometry rejects a fixed background star with a $5σ$ significance. The near-IR absolute magnitudes of the object point towards a planetary-mass companion. However, our analysis of the 1.0-3.8$μ$m spectrum extracted for the companion suggests it is a young M5.5 dwarf, based on both the 1.13-$μ$m Na index and comparison with templates of the Montreal Spectral Library. The observed spectrum is best reproduced with high effective temperature ($3057^{+119}_{-36}$K) BT-DUSTY and BT-SETTL models, but the corresponding photometric radius required to match the measured flux is only $0.60^{+0.01}_{-0.04}$ Jovian radius. We discuss possible explanations to reconcile our measurements, including an M-dwarf companion obscured by an edge-on circum-secondary disc or the shock-heated part of the photosphere of an accreting protoplanet. Follow-up observations covering a larger wavelength range and/or at finer spectral resolution are required to discriminate these two scenarios. △ Less

Submitted 20 February, 2021; originally announced February 2021.

Comments: 24 pages, 14 figures, 4 tables, to be published in MNRAS

BEAST begins: Sample characteristics and survey performance of the B-star Exoplanet Abundance Study

Authors: Markus Janson, Vito Squicciarini, Philippe Delorme, Raffaele Gratton, Mickael Bonnefoy, Sabine Reffert, Eric E. Mamajek, Simon C. Eriksson, Arthur Vigan, Maud Langlois, Natalia Engler, Gael Chauvin, Silvano Desidera, Lucio Mayer, Gabriel-Dominique Marleau, Alexander J. Bohn, Matthias Samland, Michael Meyer, Valentina d'Orazi, Thomas Henning, Sascha Quanz, Matthew Kenworthy, Joseph C. Carson

Abstract: While the occurrence rate of wide giant planets appears to increase with stellar mass at least up through the A-type regime, B-type stars have not been systematically studied in large-scale surveys so far. It therefore remains unclear up to what stellar mass this occurrence trend continues. The B-star Exoplanet Abundance Study (BEAST) is a direct imaging survey with the extreme adaptive optics ins… ▽ More While the occurrence rate of wide giant planets appears to increase with stellar mass at least up through the A-type regime, B-type stars have not been systematically studied in large-scale surveys so far. It therefore remains unclear up to what stellar mass this occurrence trend continues. The B-star Exoplanet Abundance Study (BEAST) is a direct imaging survey with the extreme adaptive optics instrument SPHERE, targeting 85 B-type stars in the young Scorpius-Centaurus (Sco-Cen) region with the aim to detect giant planets at wide separations and constrain their occurrence rate and physical properties. The statistical outcome of the survey will help determine if and where an upper stellar mass limit for planet formation occurs. In this work, we describe the selection and characterization of the BEAST target sample. Particular emphasis is placed on the age of each system, which is a central parameter in interpreting direct imaging observations. We implement a novel scheme for age dating based on kinematic sub-structures within Sco-Cen, which complements and expands upon previous age determinations in the literature. We also present initial results from the first epoch observations, including the detections of ten stellar companions, of which six were previously unknown. All planetary candidates in the survey will need follow up in second epoch observations, which are part of the allocated observational programme and will be executed in the near future. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 22 pages, 15 figures, accepted for publication in A&A

A medium-resolution spectrum of the exoplanet HIP 65426 b

Authors: Simon Petrus, Mickaël Bonnefoy, Gaël Chauvin, Benjamin Charnay, Gabriel-Dominique Marleau, Raffaele Gratton, Anne-Marie Lagrange, Julien Rameau, Chistoph Mordasini, Mathias Nowak, Philippe Delorme, Anthony Boccaletti, Alexis Carlotti, Mathis Houllé, Arthur Vigan, France Allard, Silvano Desidera, Valentina D'Orazi, Herman Jens Hoeijmakers, Aurélien Wyttenbach, Baptiste Lavie

Abstract: Medium-resolution integral-field spectrographs (IFS) coupled with adaptive-optics such as Keck/OSIRIS, VLT/MUSE, or SINFONI are appearing as a new avenue for enhancing the detection and characterization capabilities of young, gas giant exoplanets at large heliocentric distances (>5 au). We analyzed K-band VLT/SINFONI medium-resolution (R_lambda~5577) observations of the young giant exoplanet HIP 6… ▽ More Medium-resolution integral-field spectrographs (IFS) coupled with adaptive-optics such as Keck/OSIRIS, VLT/MUSE, or SINFONI are appearing as a new avenue for enhancing the detection and characterization capabilities of young, gas giant exoplanets at large heliocentric distances (>5 au). We analyzed K-band VLT/SINFONI medium-resolution (R_lambda~5577) observations of the young giant exoplanet HIP 65426 b. Our dedicated IFS data analysis toolkit (TExTRIS) optimized the cube building, star registration, and allowed for the extraction of the planet spectrum. A Bayesian inference with the nested sampling algorithm coupled with the self-consistent forward atmospheric models BT-SETTL15 and Exo-REM using the ForMoSA tool yields Teff=1560 +/- 100K, log(g)<4.40dex, [M/H]=0.05 +/- 0.22dex, and an upper limit on the C/O ratio (<0.55). The object is also re-detected with the so-called "molecular mapping" technique. The technique yields consistent atmospheric parameters, but the loss of the planet pseudo-continuum in the process degrades or modifies the constraints on these parameters. The solar to sub-solar C/O ratio suggests an enrichment by solids at formation if the planet was formed beyond the water snowline (>20 au) by core-accretion. However, a formation by gravitational instability can not be ruled out. The metallicity is compatible with the bulk enrichment of massive Jovian planets from the Bern planet population models. Finally, we measure a radial velocity of 26 +/- 15km/s compatible with our revised measurement on the star. This is the fourth imaged exoplanet for which a radial velocity can be evaluated, illustrating the potential of such observations for assessing the coevolution of imaged systems belonging to star forming regions, such as HIP 65426. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 21 pages, 16 figures, 3 Appendix

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

Spectral appearance of the planetary-surface accretion shock: Global spectra and hydrogen-line profiles and fluxes

Authors: Yuhiko Aoyama, Gabriel-Dominique Marleau, Christoph Mordasini, Masahiro Ikoma

Abstract: Hydrogen-line emission from an accretion shock has recently been observed at planetary-mass objects. Our previous work predicted the shock spectrum and luminosity for a shock on the circumplanetary disc. We extend this to the planet-surface shock. We calculate the global spectral energy distribution (SED) of accreting planets by combining our model emission spectra with photospheric SEDs, and pred… ▽ More Hydrogen-line emission from an accretion shock has recently been observed at planetary-mass objects. Our previous work predicted the shock spectrum and luminosity for a shock on the circumplanetary disc. We extend this to the planet-surface shock. We calculate the global spectral energy distribution (SED) of accreting planets by combining our model emission spectra with photospheric SEDs, and predict the line-integrated flux for several hydrogen lines, especially H alpha, but also H beta, Pa alpha, Pa beta, Pa gamma, Br alpha, and Br gamma. We apply our non-equilibrium emission model to the surface accretion shock for a wide range of accretion rates Mdot and masses M_p . Fits to formation calculations provide radii and effective temperatures. Extinction by the surrounding material is neglected, which is arguably often relevant. We find that the line luminosity increases monotonically with Mdot and M_p , depending mostly on Mdot and weakly on M_p for the relevant range of parameters. The Lyman, Balmer, and Paschen continua can exceed the photosphere. The H beta line is fainter by 0 to 1 dex than H alpha, whereas other lines are weaker (by 1 to 3 dex). Shocks on the planet or the CPD surface are distinguishable at very high spectral resolution, but the planet surface shock likely dominates if both are present. Applied to recent non-detections of H alpha, our models imply looser constraints on the Mdot of putative planets than when extrapolating fits from the stellar regime. These hydrogen-line luminosity predictions are useful for interpreting (non-)detections of accreting planets. △ Less

Submitted 21 April, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 20 pages, 8 figures (main text) + 14 pages, 7 figures (Appendix + references), re-submitted to ApJ. V.3: Model details moved to Appendix, Sect. 6.5 simplified, Appendix A on Storey & Hummer (1995) corrected, small changes. Results unchanged. Fits to radii and Teff's of accreting planets at https://github.com/gabrielastro/St-Moritz (in gnuplot/IDL/python)

The influence of infall on the properties of protoplanetary discs

Authors: O. Schib, C. Mordasini, N. Wenger, G. -D. Marleau, R. Helled

Abstract: We perform a population synthesis of protoplanetary discs including infall with a total of $50\,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a differen… ▽ More We perform a population synthesis of protoplanetary discs including infall with a total of $50\,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of $\sim 0.3$ and $0.1\,\mathrm{M_\odot}$ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems with final stellar mass $\sim 1\,\mathrm{M_\odot}$, we find disc masses of $\sim 0.7\,\mathrm{M_\odot}$ for the `hydro' case and $\sim 0.2\,\mathrm{M_\odot}$ for the `MHD' case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, $\approx 5-7\,\mathrm{Myr}$ on average, despite our choice of a high value of $10^{-2}$ for the background viscosity $α$-parameter. In addition, fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation is suppressed entirely. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: 26 pages, 14 figures. Accepted to A&A

Journal ref: A&A 645, A43 (2021)

MIRACLES: atmospheric characterization of directly imaged planets and substellar companions at 4-5 $μ$m. II. Constraints on the mass and radius of the enshrouded planet PDS 70 b

Authors: Tomas Stolker, Gabriel-Dominique Marleau, Gabriele Cugno, Paul Mollière, Sascha P. Quanz, Kamen O. Todorov, Jonas Kühn

Abstract: The circumstellar disk of PDS 70 hosts two forming planets, which are actively accreting gas from their environment. In this work, we report the first detection of PDS 70 b in the Br$α$ and $M'$ filters with VLT/NACO, a tentative detection of PDS 70 c in Br$α$, and a reanalysis of archival NACO $L'$ and SPHERE $H23$ and $K12$ imaging data. The near side of the disk is also resolved with the Br$α$… ▽ More The circumstellar disk of PDS 70 hosts two forming planets, which are actively accreting gas from their environment. In this work, we report the first detection of PDS 70 b in the Br$α$ and $M'$ filters with VLT/NACO, a tentative detection of PDS 70 c in Br$α$, and a reanalysis of archival NACO $L'$ and SPHERE $H23$ and $K12$ imaging data. The near side of the disk is also resolved with the Br$α$ and $M'$ filters, indicating that scattered light is non-negligible at these wavelengths. The spectral energy distribution of PDS 70 b is well described by blackbody emission, for which we constrain the photospheric temperature and photospheric radius to $T_\mathrm{eff}=1193 \pm 20$ K and $R=3.0 \pm 0.2$ $R_\mathrm{J}$. The relatively low bolometric luminosity, $\log(L/L_\odot) = -3.79 \pm 0.02$, in combination with the large radius, is not compatible with standard structure models of fully convective objects. With predictions from such models, and adopting a recent estimate of the accretion rate, we derive a planetary mass and radius in the range of $M_\mathrm{p}\approx 0.5-1.5$ $M_\mathrm{J}$ and $R_\mathrm{p}\approx 1-2.5$ $R_\mathrm{J}$, independently of the age and post-formation entropy of the planet. The blackbody emission, large photospheric radius, and the discrepancy between the photospheric and planetary radius suggests that infrared observations probe an extended, dusty environment around the planet, which obscures the view on its molecular composition. Finally, we derive a rough upper limit on the temperature and radius of potential excess emission from a circumplanetary disk, $T_\mathrm{eff}\lesssim256$ K and $R\lesssim245$ $R_\mathrm{J}$, but we do find weak evidence that the current data favors a model with a single blackbody component. △ Less

Submitted 19 November, 2020; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: 19 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 644, A13 (2020)

Evidence of Three Mechanisms Explaining the Radius Anomaly of Hot Jupiters

Authors: Paula Sarkis, Christoph Mordasini, Thomas Henning, Gabriel D. Marleau, Paul Mollière

Abstract: The radii of hot Jupiters are still not fully understood and all of the proposed explanations are based on the idea that these close-in giant planets possess hot interiors. We approach the radius anomaly problem by adopting a statistical approach. We infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is t… ▽ More The radii of hot Jupiters are still not fully understood and all of the proposed explanations are based on the idea that these close-in giant planets possess hot interiors. We approach the radius anomaly problem by adopting a statistical approach. We infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is the dominant one. We develop a flexible and robust hierarchical Bayesian model that couples the interior structure of exoplanets to their observed properties. We apply the model to 314 hot Jupiters and infer the internal luminosity distribution for each planet and study at the population level ({\it i}) the mass-luminosity-radius distribution and as a function of equilibrium temperature the distributions of the ({\it ii}) heating efficiency, ({\it iii}) internal temperature, and the ({\it iv}) pressure of the radiative-convective-boundary (RCB). We find that hot Jupiters tend to have high internal luminosity leading to hot interiors. This has important consequences on the cooling rate and we find that the RCB is located at low pressures. Assuming that the ultimate source of the extra heating is the irradiation from the host star, we illustrate that the heating efficiency follows a Gaussian distribution, in agreement with previous results. We discuss our findings in the context of the proposed heating mechanisms and illustrate that ohmic dissipation, advection of potential temperature, and thermal tides are in agreement with certain trends inferred from our analysis and thus all three models can explain aspects of the observations. We provide new insights on the interior structure of hot Jupiters and show that with our current knowledge it is still challenging to firmly identify the universal mechanism driving the inflated radii. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 27 pages and 12 figures. Accepted in A&A. Source code can be found at https://github.com/psarkis/bloatedHJs and data at https://www.space.unibe.ch/research/research_groups/planets_in_time/numerical_data/index_eng.html

Journal ref: A&A 645, A79 (2021)

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)

The Dynamic Proto-atmospheres around Low-Mass Planets with Eccentric Orbits

Authors: Chuhong Mai, Steven J. Desch, Rolf Kuiper, Gabriel-Dominique Marleau, Cornelis Dullemond

Abstract: Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We inves… ▽ More Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We investigate the proto-atmospheres on low-mass planets with eccentric orbits with radiation-hydrodynamics simulations. A 2D radiative model of the proto-atmosphere is established with tabulated opacities for the gas and dust. The solutions reveal large-scale gas recycling inside a bow shock structure. The atmospheres on eccentric planets are typically three to four orders of magnitude less massive than those of planets with circular orbits. Overall, however, a supersonic environment is favorable for planets to keep an early stable atmosphere, rather than harmful, due to the steady gas supply through the recycling flow. We also quantitatively explore how such atmospheres are affected by the relative velocity of the planet to the gas, the planet mass, and the background gas density. Our time-dependent simulations track the orbital evolution of the proto-atmosphere with the planet--disk parameters changing throughout the orbit. Atmospheric properties show oscillatory patterns as the planet travels on an eccentric orbit, with a lag in phase. To sum up, low-mass eccentric planets can retain small proto-atmospheres despite the stripping effects of bow shocks. The atmospheres are always connected to and interacting with the disk gas. These findings provide important insights into the impacts of migration and scattering on planetary proto-atmospheres. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJ (23 pages, 15 figures, 2 tables)

Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

Authors: P. Mollière, T. Stolker, S. Lacour, G. P. P. L. Otten, J. Shangguan, B. Charnay, T. Molyarova, M. Nowak, Th. Henning, G. -D. Marleau, D. A. Semenov, E. van Dishoeck, F. Eisenhauer, P. Garcia, R. Garcia Lopez, J. H. Girard, A. Z. Greenbaum, S. Hinkley, P. Kervella, L. Kreidberg, A. -L. Maire, E. Nasedkin, L. Pueyo, I. A. G. Snellen, A. Vigan , et al. (3 additional authors not shown)

Abstract: Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our r… ▽ More Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets. △ Less

Submitted 18 June, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A, 28 pages, 13 figures, updated author list

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

Strong H$α$ emission and signs of accretion in a circumbinary planetary mass companion from MUSE

Authors: Simon C. Eriksson, Rubén Asensio Torres, Markus Janson, Yuhiko Aoyama, Gabriel-Dominique Marleau, Mickael Bonnefoy, Simon Petrus

Abstract: Context. Intrinsic H$α$ emission can be advantageously used to detect substellar companions because it improves contrasts in direct imaging. Characterising this emission from accreting exoplanets allows for the testing of planet formation theories. Aims. We characterise the young circumbinary planetary mass companion 2MASS J01033563-5515561 (AB)b (Delorme 1 (AB)b) through medium-resolution spect… ▽ More Context. Intrinsic H$α$ emission can be advantageously used to detect substellar companions because it improves contrasts in direct imaging. Characterising this emission from accreting exoplanets allows for the testing of planet formation theories. Aims. We characterise the young circumbinary planetary mass companion 2MASS J01033563-5515561 (AB)b (Delorme 1 (AB)b) through medium-resolution spectroscopy. Methods. We used the new narrow-field mode (NFM) for the MUSE integral-field spectrograph, located on the ESO Very Large Telescope, during science verification time to obtain optical spectra of Delorme 1 (AB)b. Results. We report the discovery of very strong H$α$ and H$β$ emission, accompanied by He I emission. This is consistent with an active accretion scenario. We provide accretion rate estimates obtained from several independent methods and find a likely mass of $12-15$ M$_{\rm Jup}$ for Delorme 1 (AB)b. This is also consistent with previous estimates. Conclusions. Signs of active accretion in the Delorme 1 system might indicate a younger age than the $\sim 30-40$ Myr expected from a likely membership in Tucana-Horologium (THA). Previous works have also shown the central binary to be overluminous, which gives further indication of a younger age. However, recent discoveries of active discs in relatively old ($\sim 40$ Myr), very low-mass systems suggests that ongoing accretion in Delorme 1 (AB)b might not require in and of itself that the system is younger than the age implied by its THA membership. △ Less

Submitted 22 August, 2020; v1 submitted 24 May, 2020; originally announced May 2020.

Comments: 9 pages. Supplementary material and figures in appendix

Journal ref: A&A 638, L6 (2020)

VLT/SPHERE survey for exoplanets around young, early-type stars including systems with multi-belt architectures

Authors: M. Lombart, G. Chauvin, P. Rojo, E. Lagadec, P. Delorme, H. Beust, M. Bonnefoy, R. Galicher, R. Gratton, D. Mesa, M. Bonavita, F. Allard, A. Bayo, A. Boccaletti, S. Desidera, J. Girard, J. S. Jenkins, H. Klahr, G. Laibe, A. -M. Lagrange, C. Lazzoni, G. -D. Marleau, D. Minniti, C. Mordasini

Abstract: Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substa… ▽ More Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers. △ Less

Submitted 5 June, 2020; v1 submitted 18 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Astronomy & Astrophysics. 12 pages, 9 figures

Journal ref: A&A 639, A54 (2020)

Birth of convective low-mass to high-mass second Larson cores

Authors: Asmita Bhandare, Rolf Kuiper, Thomas Henning, Christian Fendt, Mario Flock, Gabriel-Dominique Marleau

Abstract: Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the depen… ▽ More Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 22 pages, 20 figures, and 2 tables. Accepted for publication in A&A. Collapse movie: https://keeper.mpdl.mpg.de/f/f04abdeabdf3472fb56d/

Journal ref: A&A 638, A86 (2020)

SPHERE+: Imaging young Jupiters down to the snowline

Authors: A. Boccaletti, G. Chauvin, D. Mouillet, O. Absil, F. Allard, S. Antoniucci, J. -C. Augereau, P. Barge, A. Baruffolo, J. -L. Baudino, P. Baudoz, M. Beaulieu, M. Benisty, J. -L. Beuzit, A. Bianco, B. Biller, B. Bonavita, M. Bonnefoy, S. Bos, J. -C. Bouret, W. Brandner, N. Buchschache, B. Carry, F. Cantalloube, E. Cascone , et al. (108 additional authors not shown)

Abstract: SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with S… ▽ More SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei. △ Less

Submitted 13 March, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: White paper submitted to ESO on Feb. 20th, 2020

A new take on the low-mass brown dwarf companions on wide-orbits in Upper-Scorpius

Authors: Simon Petrus, Mickaël Bonnefoy, Gaël Chauvin, Carine Babusiaux, Philippe Delorme, Anne-Marie Lagrange, Nathan Florent, Amelia Bayo, Markus Janson, Beth Biller, Elena Manjavacas, Gabriel-Dominique Marleau, Taisiya G. Kopytova

Abstract: The Upper-Scorpius association (5-11 Myr) contains a unique population of low-mass (M<30 MJup) brown-dwarfs. The detailed relative characterization of their physical properties (mass, radius, temperature, composition, ongoing accretion) offers the opportunity to potentially explore their origin and their mechanisms of formation. We aim at characterizing the chemical and physical properties of thre… ▽ More The Upper-Scorpius association (5-11 Myr) contains a unique population of low-mass (M<30 MJup) brown-dwarfs. The detailed relative characterization of their physical properties (mass, radius, temperature, composition, ongoing accretion) offers the opportunity to potentially explore their origin and their mechanisms of formation. We aim at characterizing the chemical and physical properties of three young, late-M brown-dwarfs claimed to be companions of the Upper-Scorpius stars USco161031.9-16191305, HIP77900, and HIP78530 using medium resolution spectroscopy at UV (R~3300), optical (R~5400), and near-infrared (R~4300) wavelengths. The spectra of six free-floating analogues from the same association are analyzed for comparison and to explore the potential physical differences between these substellar objects. We also aim at analyzing emission lines at UV and optical wavelengths to investigate the presence of ongoing accretion processes. The X-Shooter spectrograph at VLT was used to obtain the spectra of our targets over the 0.3-2.5 $μ$m range simultaneously. Performing a forward modelling of the observed spectra with the ForMoSA code (Nested Sampling Bayesian inference method in using BT-SETTL15 models), we infer the Teff, log(g), and radius of our objects. Our results are compatible with physical parameters found for free floating analogues of the association and with evolutionary-model predictions. However the final accuracy on the Teff estimates is strongly limited by non-reproducibilities of the BT-SETTL15 models for late-M brown-dwarfs. We identified emission lines in the spectrum of several objects attributed to chromospheric activity except for USco1608-2315 for which they are indicative of active accretion. We confirm the x4 over-luminosity of USco161031.9-16191305B. It could be explained by the object activity and if the companion is an unresolved multiple system. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 21 pages, 15 figures, 10 tables

Report number: AA/2019/35732

Journal ref: A&A 633, A124 (2020)

The Planetary Accretion Shock. II. Grid of Post-Shock Entropies and Radiative Shock Efficiencies for Non-Equilibrium Radiation Transport

Authors: Gabriel-Dominique Marleau, Christoph Mordasini, Rolf Kuiper

Abstract: In the core-accretion formation scenario of gas giants, most of the gas accreting onto a planet is processed through an accretion shock. In this series of papers we study this shock since it is key in setting the forming planet's structure and thus its post-formation luminosity, with dramatic observational consequences. We perform one-dimensional grey radiation-hydrodynamical simulations with non-… ▽ More In the core-accretion formation scenario of gas giants, most of the gas accreting onto a planet is processed through an accretion shock. In this series of papers we study this shock since it is key in setting the forming planet's structure and thus its post-formation luminosity, with dramatic observational consequences. We perform one-dimensional grey radiation-hydrodynamical simulations with non-equilibrium (two-temperature) radiation transport and up-to-date opacities. We survey the parameter space of accretion rate, planet mass, and planet radius and obtain post-shock temperatures, pressures, and entropies, as well as global radiation efficiencies. We find that usually, the shock temperature T_shock is given by the "free-streaming" limit. At low temperatures the dust opacity can make the shock hotter but not significantly. We corroborate this with an original semi-analytical derivation of T_shock . We also estimate the change in luminosity between the shock and the nebula. Neither T_shock nor the luminosity profile depend directly on the optical depth between the shock and the nebula. Rather, T_shock depends on the immediate pre-shock opacity, and the luminosity change on the equation of state (EOS). We find quite high immediate post-shock entropies (S ~ 13-20 kB/mH), which makes it seem unlikely that the shock can cool the planet. The global radiation efficiencies are high (eta^phys > 97%) but the remainder of the total incoming energy, which is brought into the planet, exceeds the internal luminosity of classical cold starts by orders of magnitude. Overall, these findings suggest that warm or hot starts are more plausible. △ Less

Submitted 13 June, 2019; originally announced June 2019.

Comments: 27 pages, 11 figures. In press at ApJ

The B-Star Exoplanet Abundance Study: a co-moving 16-25 Mjup companion to the young binary system HIP 79098

Authors: Markus Janson, Ruben Asensio-Torres, Damien Andre, Mickael Bonnefoy, Philippe Delorme, Sabine Reffert, Silvano Desidera, Maud Langlois, Gael Chauvin, Raffaele Gratton, Alexander J. Bohn, Simon C. Eriksson, Gabriel-Dominique Marleau, Eric E. Mamajek, Arthur Vigan, Joseph C. Carson

Abstract: Wide low-mass substellar companions are known to be very rare among low-mass stars, but appear to become increasingly common with increasing stellar mass. However, B-type stars, which are the most massive stars within ~150 pc of the Sun, have not yet been examined to the same extent as AFGKM-type stars in that regard. In order to address this issue, we launched the ongoing B-star Exoplanet Abundan… ▽ More Wide low-mass substellar companions are known to be very rare among low-mass stars, but appear to become increasingly common with increasing stellar mass. However, B-type stars, which are the most massive stars within ~150 pc of the Sun, have not yet been examined to the same extent as AFGKM-type stars in that regard. In order to address this issue, we launched the ongoing B-star Exoplanet Abundance Study (BEAST) to examine the frequency and properties of planets, brown dwarfs, and disks around B-type stars in the Scorpius-Centaurus (Sco-Cen) association; we also analyzed archival data of B-type stars in Sco-Cen. During this process, we identified a candidate substellar companion to the B9-type spectroscopic binary HIP 79098 AB, which we refer to as HIP 79098 (AB)b. The candidate had been previously reported in the literature, but was classified as a background contaminant on the basis of its peculiar colors. Here we demonstrate that the colors of HIP 79098 (AB)b are consistent with several recently discovered young and low-mass brown dwarfs, including other companions to stars in Sco-Cen. Furthermore, we show unambiguous common proper motion over a 15-year baseline, robustly identifying HIP 79098 (AB)b as a bona fide substellar circumbinary companion at a 345+/-6 AU projected separation to the B9-type stellar pair. With a model-dependent mass of 16-25 Mjup yielding a mass ratio of <1%, HIP 79098 (AB)b joins a growing number of substellar companions with planet-like mass ratios around massive stars. Our observations underline the importance of common proper motion analysis in the identification of physical companionship, and imply that additional companions could potentially remain hidden in the archives of purely photometric surveys. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: 9 pages, 6 figures, accepted for publication in A&A

Self-luminous and irradiated exoplanetary atmospheres explored with HELIOS

Authors: Matej Malik, Daniel Kitzmann, João M. Mendonça, Simon L. Grimm, Gabriel-Dominique Marleau, Esther F. Linder, Shang-Min Tsai, Kevin Heng

Abstract: We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sou… ▽ More We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80°, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 33 pages, 18 figures, accepted for publication in AJ

Exploring the formation by core accretion and the luminosity evolution of directly imaged planets: The case of HIP 65426 b

Authors: Gabriel-Dominique Marleau, Gavin A. L. Coleman, Adrien Leleu, Christoph Mordasini

Abstract: A low-mass companion to the two-solar mass star HIP65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. We want to derive statistical constraints on the mass and initial entropy of HIP65426b and to explore possible formation pathways of directl… ▽ More A low-mass companion to the two-solar mass star HIP65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. We want to derive statistical constraints on the mass and initial entropy of HIP65426b and to explore possible formation pathways of directly imaged objects within the core-accretion paradigm, focusing on HIP65426b. Constraints on the planet's mass and post-formation entropy are derived from its age and luminosity combined with cooling models. For the first time, the results of population synthesis are also used to inform the results. Then, a formation model that includes N-body dynamics with several embryos per disc is used to study possible formation histories and the properties of possible additional companions. Finally, the outcomes of two- and three-planet scattering in the post-disc phase are analysed, taking tides into account. The mass of HIP65426b is found to be Mp = 9.9 +1.1 -1.8 MJ using the hot population and Mp = 10.9 +1.4 -2.0 MJ with the cold-nominal population. Core formation at small separations from the star followed by outward scattering and runaway accretion at a few hundred AU succeeds in reproducing the mass and separation of HIP65426b. Alternatively, systems having two or more giant planets close enough to be on an unstable orbit at disc dispersal are likely to end up with one planet on a wide HIP65426b-like orbit with a relatively high eccentricity (>~ 0.5). If this scattering scenario explains its formation, HIP65426b is predicted to have a high eccentricity and to be accompanied by one or several roughly Jovian-mass planets at smaller semi-major axes, which also could have a high eccentricity. This could be tested by further direct-imaging as well as radial-velocity observations. △ Less

Submitted 20 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Comments: 17 pages, 11 figures. A&A in press. Bern EXoplanet cooling curves (BEX) available upon request. v2: Language and other minor changes; Fig. 4 now has labels summarising a possible formation pathway discussed in the text

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