-
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.
-
The ESO SupJup Survey I: Chemical and isotopic characterisation of the late L-dwarf DENIS J0255-4700 with CRIRES$^+$
Authors:
S. de Regt,
S. Gandhi,
I. A. G. Snellen,
Y. Zhang,
C. Ginski,
D. González Picos,
A. Y. Kesseli,
R. Landman,
P. Mollière,
E. Nasedkin,
A. Sánchez-López,
T. Stolker
Abstract:
It has been proposed that the distinct formation and evolution of exoplanets and brown dwarfs may affect the chemical and isotopic content of their atmospheres. Recent work has indeed shown differences in the $^{12}$C/$^{13}$C isotope ratio, provisionally attributed to the top-down formation of brown dwarfs and the core accretion pathway of super-Jupiters. The ESO SupJup Survey aims to disentangle…
▽ More
It has been proposed that the distinct formation and evolution of exoplanets and brown dwarfs may affect the chemical and isotopic content of their atmospheres. Recent work has indeed shown differences in the $^{12}$C/$^{13}$C isotope ratio, provisionally attributed to the top-down formation of brown dwarfs and the core accretion pathway of super-Jupiters. The ESO SupJup Survey aims to disentangle the formation pathways of isolated brown dwarfs and planetary-mass companions using chemical and isotopic tracers. The survey uses high-resolution spectroscopy with the recently upgraded VLT/CRIRES$^+$ spectrograph, covering a total of 49 targets. Here, we present the first results: an atmospheric characterisation of DENIS J0255-4700, an isolated brown dwarf near the L-T transition. We analyse its K-band spectrum using a retrieval framework where the radiative transfer code petitRADTRANS is coupled to PyMultiNest. Gaussian Processes are employed to model inter-pixel correlations and we adopt an updated parameterisation of the PT-profile. Abundances of CO, H$_2$O, CH$_4$, and NH$_3$ are retrieved for this fast-rotating L-dwarf. The ExoMol H$_2$O line list provides a significantly better fit than that of HITEMP. A free-chemistry retrieval is strongly favoured over equilibrium chemistry, caused by an under-abundance of CH$_4$. The free-chemistry retrieval constrains a super-solar C/O-ratio of $\sim0.68$ and a solar metallicity. We find tentative evidence ($\sim3σ$) for the presence of $^{13}$CO, with a constraint on the isotope ratio of $\mathrm{^{12}C/^{13}C}=184^{+61}_{-40}$ and a lower limit of $\gtrsim97$, suggesting a depletion of $^{13}$C compared to the interstellar medium ($\sim68$). High-resolution, high signal-to-noise K-band spectra provide an excellent means to constrain the chemistry and isotopic content of sub-stellar objects, as is the main objective of the ESO SupJup Survey.
△ Less
Submitted 17 May, 2024;
originally announced May 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 4 April, 2024;
originally announced April 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.
-
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.
-
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.
-
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.
-
Atmospheric retrievals with petitRADTRANS
Authors:
Evert Nasedkin,
Paul Mollière,
Doriann Blain
Abstract:
petitRADTRANS (pRT) is a fast radiative transfer code used for computing emission and transmission spectra of exoplanet atmospheres, combining a FORTRAN back end with a Python based user interface. It is widely used in the exoplanet community with 222 references in the literature to date, and has been benchmarked against numerous similar tools. The spectra calculated with pRT can be used as a forw…
▽ More
petitRADTRANS (pRT) is a fast radiative transfer code used for computing emission and transmission spectra of exoplanet atmospheres, combining a FORTRAN back end with a Python based user interface. It is widely used in the exoplanet community with 222 references in the literature to date, and has been benchmarked against numerous similar tools. The spectra calculated with pRT can be used as a forward model for fitting spectroscopic data using Monte Carlo techniques, commonly referred to as an atmospheric retrieval. The new retrieval module combines fast forward modelling with nested sampling codes, allowing for atmospheric retrievals on a large range of different types of exoplanet data. Thus it is now possible to use pRT to easily and quickly infer the atmospheric properties of exoplanets in both transmission and thermal emission.
△ Less
Submitted 3 April, 2024; v1 submitted 13 September, 2023;
originally announced September 2023.
-
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.
-
Impacts of high-contrast image processing on atmospheric retrievals
Authors:
Evert Nasedkin,
Paul Mollière,
Jason Wang,
Faustine Cantalloube,
Laura Kreidberg,
Laurent Pueyo,
Tomas Stolker,
Arthur Vigan
Abstract:
Many post-processing algorithms have been developed in order to better separate the signal of a companion from the bright light of the host star, but the effect of such algorithms on the shape of exoplanet spectra extracted from integral field spectrograph data is poorly understood. The resulting spectra are affected by noise that is correlated in wavelength space due to both optical and data proc…
▽ More
Many post-processing algorithms have been developed in order to better separate the signal of a companion from the bright light of the host star, but the effect of such algorithms on the shape of exoplanet spectra extracted from integral field spectrograph data is poorly understood. The resulting spectra are affected by noise that is correlated in wavelength space due to both optical and data processing effects. Within the framework of Bayesian atmospheric retrievals, we aim to understand how these correlations and other systematic effects impact the inferred physical parameters. We consider three algorithms (KLIP, PynPoint and ANDROMEDA), optimizing the choice of algorithmic parameters using a series of injection tests into archival SPHERE and GPI data of the HR 8799 system. The wavelength-dependent covariance matrix is calculated to provide a measure of instrumental and algorithmic systematics. We perform atmospheric retrievals using petitRADTRANS on optimally extracted spectra to measure how these data processing systematics influence the retrieved parameter distributions. The choice of data processing algorithm and parameters significantly impact the accuracy of retrieval results, with the mean posterior parameter bias ranging from 1 to 3 $σ$ from the true input parameters. Including the full covariance matrix in the likelihood improves the accuracy of inferred parameters, and cannot be accounted for using ad hoc scaling parameters in the retrieval framework. Using the Bayesian information criterion and other statistical measures as a heuristic goodness-of-fit metrics, the retrievals including the full covariance matrix are favoured when compared to using only the diagonal elements.
△ Less
Submitted 2 August, 2023;
originally announced August 2023.
-
Neural posterior estimation for exoplanetary atmospheric retrieval
Authors:
Malavika Vasist,
François Rozet,
Olivier Absil,
Paul Mollière,
Evert Nasedkin,
Gilles Louppe
Abstract:
Retrieving the physical parameters from spectroscopic observations of exoplanets is key to understanding their atmospheric properties. Exoplanetary atmospheric retrievals are usually based on approximate Bayesian inference and rely on sampling-based approaches to compute parameter posterior distributions. Accurate or repeated retrievals, however, can result in very long computation times due to th…
▽ More
Retrieving the physical parameters from spectroscopic observations of exoplanets is key to understanding their atmospheric properties. Exoplanetary atmospheric retrievals are usually based on approximate Bayesian inference and rely on sampling-based approaches to compute parameter posterior distributions. Accurate or repeated retrievals, however, can result in very long computation times due to the sequential nature of sampling-based algorithms. We aim to amortize exoplanetary atmospheric retrieval using neural posterior estimation (NPE), a simulation-based inference algorithm based on variational inference and normalizing flows. In this way, we aim (i) to strongly reduce inference time, (ii) to scale inference to complex simulation models with many nuisance parameters or intractable likelihood functions, and (iii) to enable the statistical validation of the inference results. We evaluate NPE on a radiative transfer model for exoplanet spectra petitRADTRANS, including the effects of scattering and clouds. We train a neural autoregressive flow to quickly estimate posteriors and compare against retrievals computed with MultiNest. NPE produces accurate posterior approximations while reducing inference time down to a few seconds. We demonstrate the computational faithfulness of our posterior approximations using inference diagnostics including posterior predictive checks and coverage, taking advantage of the quasi-instantaneous inference time of NPE. Our analysis confirms the reliability of the approximate posteriors produced by NPE. The accuracy and reliability of the inference results produced by NPE establishes it as a promising approach for atmospheric retrievals. Amortization of the posterior inference makes repeated inference on several observations computationally inexpensive since it does not require on-the-fly simulations, making the retrieval efficient, scalable, and testable.
△ Less
Submitted 10 February, 2023; v1 submitted 16 January, 2023;
originally announced January 2023.
-
Photochemically-produced SO$_2$ in the atmosphere of WASP-39b
Authors:
Shang-Min Tsai,
Elspeth K. H. Lee,
Diana Powell,
Peter Gao,
Xi Zhang,
Julianne Moses,
Eric Hébrard,
Olivia Venot,
Vivien Parmentier,
Sean Jordan,
Renyu Hu,
Munazza K. Alam,
Lili Alderson,
Natalie M. Batalha,
Jacob L. Bean,
Björn Benneke,
Carver J. Bierson,
Ryan P. Brady,
Ludmila Carone,
Aarynn L. Carter,
Katy L. Chubb,
Julie Inglis,
Jérémy Leconte,
Mercedes Lopez-Morales,
Yamila Miguel
, et al. (60 additional authors not shown)
Abstract:
Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WA…
▽ More
Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $μ$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$σ$) and G395H (4.5$σ$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations.
△ Less
Submitted 24 March, 2023; v1 submitted 18 November, 2022;
originally announced November 2022.
-
Exoplanet Imaging Data Challenge, phase II: Characterization of exoplanet signals in high-contrast images
Authors:
F. Cantalloube,
V. Christiaens,
C. Cantero,
E. Nasedkin,
A. Cioppa,
O. Absil,
J. M. Bonse,
P. Delorme,
C. Gomez-Gonzalez,
S. Juillard,
J. Mazoyer,
M. Samland Ruffio J. -B. i,
Van Droogenbroeck M. c
Abstract:
Today, there exists a wide variety of algorithms dedicated to high-contrast imaging, especially for the detection and characterisation of exoplanet signals. These algorithms are tailored to address the very high contrast between the exoplanet signal(s), which can be more than two orders of magnitude fainter than the bright starlight residuals in coronagraphic images. The starlight residuals are in…
▽ More
Today, there exists a wide variety of algorithms dedicated to high-contrast imaging, especially for the detection and characterisation of exoplanet signals. These algorithms are tailored to address the very high contrast between the exoplanet signal(s), which can be more than two orders of magnitude fainter than the bright starlight residuals in coronagraphic images. The starlight residuals are inhomogeneously distributed and follow various timescales that depend on the observing conditions and on the target star brightness. Disentangling the exoplanet signals within the starlight residuals is therefore challenging, and new post-processing algorithms are striving to achieve more accurate astrophysical results. The Exoplanet Imaging Data Challenge is a community-wide effort to develop, compare and evaluate algorithms using a set of benchmark high-contrast imaging datasets. After a first phase ran in 2020 and focused on the detection capabilities of existing algorithms, the focus of this ongoing second phase is to compare the characterisation capabilities of state-of-the-art techniques. The characterisation of planetary companions is two-fold: the astrometry (estimated position with respect to the host star) and spectrophotometry (estimated contrast with respect to the host star, as a function of wavelength). The goal of this second phase is to offer a platform for the community to benchmark techniques in a fair, homogeneous and robust way, and to foster collaborations.
△ Less
Submitted 16 September, 2022;
originally announced September 2022.
-
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.
-
Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions
Authors:
Paul Mollière,
Tamara Molyarova,
Bertram Bitsch,
Thomas Henning,
Aaron Schneider,
Laura Kreidberg,
Christian Eistrup,
Remo Burn,
Evert Nasedkin,
Dmitry Semenov,
Christoph Mordasini,
Martin Schlecker,
Kamber R. Schwarz,
Sylvestre Lacour,
Mathias Nowak,
Matthäus Schulik
Abstract:
Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric comp…
▽ More
Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations.
△ Less
Submitted 28 April, 2022;
originally announced April 2022.
-
Direct emission spectroscopy of exoplanets with the medium resolution imaging spectrometer on board JWST MIRI: I. Molecular mapping and sensitivity to instrumental effects
Authors:
P. Patapis,
E. Nasedkin,
G. Cugno,
A. M. Glauser,
I. Argyriou,
N. P. Whiteford,
P. Mollière,
A. Glasse,
S. P. Quanz
Abstract:
The Medium Resolution Spectrometer on board JWST/MIRI will give access to mid-IR spectra while retaining spatial information. With the unparalleled sensitivity of JWST and the MIRI detectors, the MRS has the potential to revolutionise our understanding of giant exoplanet atmospheres. Molecular mapping is a promising detection and characterisation technique used to study the spectra of directly ima…
▽ More
The Medium Resolution Spectrometer on board JWST/MIRI will give access to mid-IR spectra while retaining spatial information. With the unparalleled sensitivity of JWST and the MIRI detectors, the MRS has the potential to revolutionise our understanding of giant exoplanet atmospheres. Molecular mapping is a promising detection and characterisation technique used to study the spectra of directly imaged exoplanets. We aim to examine the feasibility and application of this technique to MRS observations. We used the instrument simulator MIRISIM to create mock observations of resolved star and exoplanet systems. As an input for the simulator, we used stellar and planet parameters from literature, with the planet spectrum being modelled with the radiative transfer code petitRADTRANS. After processing the raw data with the JWST pipeline, we high pass filter the data to account for the stellar point spread function, and used a forward modelling approach to detect the companions and constrain the chemical composition of their atmospheres through their molecular signatures. We identified limiting factors in spectroscopic characterisation of directly imaged exoplanets with the MRS and simulated observations of two representative systems, HR8799 and GJ504. In both systems, we could detect the presence of multiple molecules that were present in the input model of their atmospheres. We used two different approaches with single molecule forward models, used in literature, that are sensitive to detecting mainly H$_2$O, CO, CH$_4$, and NH$_3$, and a log-likelihood ratio test that uses full atmosphere forward models and is sensitive to a larger number of less dominant molecular species. We show that the MIRI MRS can be used to characterise widely separated giant exoplanets in the mid-IR using molecular mapping.
△ Less
Submitted 29 October, 2021;
originally announced October 2021.
-
The mass of Beta Pictoris c from Beta Pictoris b orbital motion
Authors:
S. Lacour,
J. J. Wang,
L. Rodet,
M. Nowak,
J. Shangguan,
H. Beust,
A. -M. Lagrange,
R. Abuter,
A. Amorim,
R. Asensio-Torres,
M. Benisty,
J. -P. Berger,
S. Blunt,
A. Boccaletti,
A. Bohn,
M. -L. Bolzer,
M. Bonnefoy,
H. Bonnet,
G. Bourdarot,
W. Brandner,
F. Cantalloube,
P. Caselli,
B. Charnay,
G. Chauvin,
E. Choquet
, et al. (74 additional authors not shown)
Abstract:
We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was…
▽ More
We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $β$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $β$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $β$ Pictoris c are now constrained as precisely as those of $β$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
GRAVITY K-band spectroscopy of HD 206893 B: brown dwarf or exoplanet
Authors:
J. Kammerer,
S. Lacour,
T. Stolker,
P. Mollière,
D. K. Sing,
E. Nasedkin,
P. Kervella,
J. J. Wang,
K. Ward-Duong,
M. Nowak,
R. Abuter,
A. Amorim,
R. Asensio-Torres,
M. Bauböck,
M. Benisty,
J. -P. Berger,
H. Beust,
S. Blunt,
A. Boccaletti,
A. Bohn,
M. -L. Bolzer,
M. Bonnefoy,
H. Bonnet,
W. Brandner,
F. Cantalloube
, et al. (72 additional authors not shown)
Abstract:
We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectra…
▽ More
We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~μ\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~μ\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$.
△ Less
Submitted 15 June, 2021;
originally announced June 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.
-
The ExoGRAVITY project: using single mode interferometry to characterize exoplanets
Authors:
S. Lacour,
J. J. Wang,
M. Nowak,
L. Pueyo,
F. Eisenhauer,
A. -M. Lagrange,
P. Mollière,
R. Abuter,
A. Amorim,
R. Asensio-Torres,
M. Bauböck,
M. Benisty,
J. P. Berger,
H. Beust,
S. Blunt,
A. Boccaletti,
A. Bohn,
M. Bonnefoy,
H. Bonnet,
W. Brandner,
F. Cantalloube,
P. Caselli,
B. Charnay,
G. Chauvin,
E. Choquet
, et al. (67 additional authors not shown)
Abstract:
Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orb…
▽ More
Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.
△ Less
Submitted 19 January, 2021; v1 submitted 18 January, 2021;
originally announced January 2021.
-
Exoplanet Imaging Data Challenge: benchmarking the various image processing methods for exoplanet detection
Authors:
F. Cantalloube,
C. Gomez-Gonzalez,
O. Absil,
C. Cantero,
R. Bacher,
M. J. Bonse,
M. Bottom,
C. -H. Dahlqvist,
C. Desgrange,
O. Flasseur,
T. Fuhrmann,
Th. Henning,
R. Jensen-Clem,
M. Kenworthy,
D. Mawet,
D. Mesa,
T. Meshkat,
D. Mouillet,
A. Mueller,
E. Nasedkin,
B. Pairet,
S. Pierard,
J. -B. Ruffio,
M. Samland,
J. Stone
, et al. (1 additional authors not shown)
Abstract:
The Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is host…
▽ More
The Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is hosted on the CodaLab competition platform, where participants can upload their results. The specifications of the data challenge are published on our website. The first phase, launched on the 1st of September 2019 and closed on the 1st of October 2020, consisted in detecting point sources in two types of common data-set in the field of high-contrast imaging: data taken in pupil-tracking mode at one wavelength (subchallenge 1, also referred to as ADI) and multispectral data taken in pupil-tracking mode (subchallenge 2, also referred to as ADI mSDI). In this paper, we describe the approach, organisational lessons-learnt and current limitations of the data challenge, as well as preliminary results of the participants submissions for this first phase. In the future, we plan to provide permanent access to the standard library of data sets and metrics, in order to guide the validation and support the publications of innovative image processing algorithms dedicated to high-contrast imaging of planetary systems.
△ Less
Submitted 13 January, 2021;
originally announced January 2021.
-
Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY
Authors:
J. J. Wang,
A. Vigan,
S. Lacour,
M. Nowak,
T. Stolker,
R. J. De Rosa,
S. Ginzburg,
P. Gao,
R. Abuter,
A. Amorim,
R. Asensio-Torres,
M. Baubck,
M. Benisty,
J. P. Berger,
H. Beust,
J. -L. Beuzit,
S. Blunt,
A. Boccaletti,
A. Bohn,
M. Bonnefoy,
H. Bonnet,
W. Brandner,
F. Cantalloube,
P. Caselli,
B. Charnay
, et al. (79 additional authors not shown)
Abstract:
We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbi…
▽ More
We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.
△ Less
Submitted 3 February, 2021; v1 submitted 11 January, 2021;
originally announced January 2021.
-
Direct confirmation of the radial-velocity planet $β$ Pic c
Authors:
M. Nowak,
S. Lacour,
A. -M. Lagrange,
P. Rubini,
J. Wang,
T. Stolker,
A. Amorim,
R. Asensio-Torres,
M. Bauböck,
M. Benisty,
J. P. Berger,
H. Beust,
S. Blunt,
A. Boccaletti,
M. Bonnefoy,
H. Bonnet,
W. Brandner,
F. Cantalloube,
B. Charnay,
E. Choquet,
V. Christiaens,
Y. Clénet,
V. Coudé du Foresto,
A. Cridland,
P. T. de Zeeuw
, et al. (68 additional authors not shown)
Abstract:
Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. %, and thus on long orbital period. This dichotomy makes it difficult to combine the two techniques on a single ta…
▽ More
Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. %, and thus on long orbital period. This dichotomy makes it difficult to combine the two techniques on a single target at once. Simultaneous measurements made by direct and indirect techniques offer the possibility of determining the mass and luminosity of planets and a method of testing formation models. Here, we aim to show how long-baseline interferometric observations guided by radial-velocity can be used in such a way. We observed the recently-discovered giant planet $β$ Pictoris c with GRAVITY, mounted on the Very Large Telescope Interferometer (VLTI). This study constitutes the first direct confirmation of a planet discovered through radial velocity. We find that the planet has a temperature of $T = 1250\pm50$\,K and a dynamical mass of $M = 8.2\pm0.8\,M_{\rm Jup}$. At $18.5\pm2.5$\,Myr, this puts $β$ Pic c close to a 'hot start' track, which is usually associated with formation via disk instability. Conversely, the planet orbits at a distance of 2.7\,au, which is too close for disk instability to occur. The low apparent magnitude ($M_{\rm K} = 14.3 \pm 0.1$) favours a core accretion scenario. We suggest that this apparent contradiction is a sign of hot core accretion, for example, due to the mass of the planetary core or the existence of a high-temperature accretion shock during formation.
△ Less
Submitted 9 October, 2020;
originally announced October 2020.
-
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.
-
A holistic perspective on the dynamics of G035.39-00.33: the interplay between gas and magnetic fields
Authors:
Tie Liu,
Pak Shing Li,
Mika Juvela,
Kee-Tae Kim,
Neal J. Evans II,
James Di Francesco,
Sheng-Yuan Liu,
Jinghua Yuan,
Ken'ichi Tatematsu,
Qizhou Zhang,
Derek Ward-Thompson,
Gary Fuller,
Paul F. Goldsmith,
P. M. Koch,
Patricio Sanhueza,
I. Ristorcelli,
Sung-ju Kang,
Huei-Ru Chen,
N. Hirano,
Yuefang Wu,
Vlas Sokolov,
Chang Won Lee,
Glenn J. White,
Ke Wang,
David Eden
, et al. (32 additional authors not shown)
Abstract:
Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at…
▽ More
Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $μ$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.
△ Less
Submitted 18 June, 2018; v1 submitted 26 March, 2018;
originally announced March 2018.
