Showing 1–50 of 153 results for author: Min, M

The PLATO Mission

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

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

Submitted 8 June, 2024; originally announced June 2024.

Constraining the formation of WASP-39b using JWST transit spectroscopy

Authors: N. Khorshid, M. Min, J. Polman, L. B. F. M. Waters

Abstract: Understanding the formation history of planets is one of the goals of studying exoplanet atmospheres. The atmospheric composition of planets can provide insights into the formation pathways of planets. Even though the mapping of the atmospheric composition onto a formation pathway is not unambiguous, with the increasing sensitivity of modern instruments, we can derive promising constraints. In thi… ▽ More Understanding the formation history of planets is one of the goals of studying exoplanet atmospheres. The atmospheric composition of planets can provide insights into the formation pathways of planets. Even though the mapping of the atmospheric composition onto a formation pathway is not unambiguous, with the increasing sensitivity of modern instruments, we can derive promising constraints. In this work, we aim to understand the formation pathway of WASP-39b. We discuss whether the detection of SO2 in its atmosphere would impact our understanding of the formation of the planet and whether it enables us to determine the formation pathway of the planet with greater accuracy. We used the JWST transit observation of the planet together with the available HST and Spitzer observations. We used a formation model coupled with a radiative transfer retrieval model to derive the planet's atmospheric characteristics and formation history. Furthermore, we used a photochemical model to derive the impact of photochemistry on the atmosphere of the planet. In this work, we show that the planet is most likely to have initiated beyond the CO2 ice line of its natal disk. Furthermore, the planet is likely to have have accreted some planetesimals during its formation. We show that the sulfur abundance in the atmosphere of the planet is probably lower than $2.27 \times 10^{-4}$. This abundance indicates that the planet is likely to exhibit a lower metallicity than suggested by the retrievals. Furthermore, such an abundance for sulfur is more likely if WASP-39b had been formed beyond the CO ice line of its natal disk. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 10 pages, 7 figures, 5 tables

Journal ref: A&A, volume 685, Article A64, Year 2024

Bayesian Analysis of Molecular Emission and Dust Continuum of Protoplanetary Disks

Authors: T. Kaeufer, M. Min, P. Woitke, I. Kamp, A. M. Arabhavi

Abstract: The Mid-InfraRed Instrument (MIRI) on board the James Webb Space Telescope (JWST) probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail, complicating interpretations which are mainly based on manual continuum subtraction and 0D slab models. We investigate the physical conditions under which the gas emits in pro… ▽ More The Mid-InfraRed Instrument (MIRI) on board the James Webb Space Telescope (JWST) probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail, complicating interpretations which are mainly based on manual continuum subtraction and 0D slab models. We investigate the physical conditions under which the gas emits in protoplanetary disks. Based on MIRI spectra, we apply a full Bayesian analysis that provides the posterior distributions of dust and molecular properties. For doing so, we introduce the Dust Continuum Kit with Line emission from Gas (DuCKLinG), a model describing the molecular line emission and the dust continuum simultaneously without large computational cost. The dust model is based on work by Juhasz et al. (2009, 2010). The molecular emission is based on LTE slab models, but with radial gradients in column densities and temperatures. The model is compared to observations using Bayesian analysis. We benchmark this model to a complex thermo-chemical ProDiMo model and fit the MIRI spectrum of GWLup. We find that the retrieved molecular conditions from DuCKLinG fall within the true values from ProDiMo. The column densities retrieved by Grant et al. (2023) fall within the retrieved ranges in this study for all examined molecules (CO2, H2O, HCN, and C2H2). Similar overlap is found for the temperatures with only the temperature range of HCN not including the previously found value. This discrepancy may be due to the simultaneous fitting of all molecules compared to the step-by-step fitting of the previous study. There is statistically significant evidence for radial temperature and column density gradients for H2O and CO2 compared to the constant temperature and column density assumed in the 0D slab models. Additionally, HCN and C2H2 emit from a small region with near constant conditions. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: accepted by Astronomy & Astrophysics

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

FlopPITy: Enabling self-consistent exoplanet atmospheric retrievals with machine learning

Authors: Francisco Ardévol Martínez, Michiel Min, Daniela Huppenkothen, Inga Kamp, Paul I. Palmer

Abstract: Interpreting the observations of exoplanet atmospheres to constrain physical and chemical properties is typically done using Bayesian retrieval techniques. Because these methods require many model computations, a compromise is made between model complexity and run time. Reaching this compromise leads to the simplification of many physical and chemical processes (e.g. parameterised temperature stru… ▽ More Interpreting the observations of exoplanet atmospheres to constrain physical and chemical properties is typically done using Bayesian retrieval techniques. Because these methods require many model computations, a compromise is made between model complexity and run time. Reaching this compromise leads to the simplification of many physical and chemical processes (e.g. parameterised temperature structure). Here we implement and test sequential neural posterior estimation (SNPE), a machine learning inference algorithm, for exoplanet atmospheric retrievals. The goal is to speed up retrievals so they can be run with more computationally expensive atmospheric models, such as those computing the temperature structure using radiative transfer. We generate 100 synthetic observations using ARCiS (ARtful Modeling Code for exoplanet Science, an atmospheric modelling code with the flexibility to compute models in varying degrees of complexity) and perform retrievals on them to test the faithfulness of the SNPE posteriors. The faithfulness quantifies whether the posteriors contain the ground truth as often as we expect. We also generate a synthetic observation of a cool brown dwarf using the self-consistent capabilities of ARCiS and run a retrieval with self-consistent models to showcase the possibilities that SNPE opens. We find that SNPE provides faithful posteriors and is therefore a reliable tool for exoplanet atmospheric retrievals. We are able to run a self-consistent retrieval of a synthetic brown dwarf spectrum using only 50,000 forward model evaluations. We find that SNPE can speed up retrievals between $\sim2\times$ and $\geq10\times$ depending on the computational load of the forward model, the dimensionality of the observation, and the signal-to-noise ratio of the observation. We make the code publicly available for the community on Github. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted for publication at A&A

Retrieving planet formation parameters of WASP-77Ab using SimAb

Authors: N. Khorshid, M. Min, J. M. Désert

Abstract: The atmospheric compositions of planets offer a unique view into their respective formation processes. State-of-the-art observatories and techniques are finally able to provide high-precision data on atmospheric composition that can be used to constrain planet formation. In this context, we focus on the formation of WASP-77Ab based on previous observations of its atmosphere, which have provided pr… ▽ More The atmospheric compositions of planets offer a unique view into their respective formation processes. State-of-the-art observatories and techniques are finally able to provide high-precision data on atmospheric composition that can be used to constrain planet formation. In this context, we focus on the formation of WASP-77Ab based on previous observations of its atmosphere, which have provided precise C/O and metallicity measurements. We use the SimAb planet formation simulation to model the formation of WASP-77Ab. We assume two compositions for the disk WASP-77Ab was formed within: one of a solar composition and one that represents the composition of WASP-77A. In addition, we considered two different scenarios regarding the migration of the planet and we study the possible planet formation paths that reproduce the composition of WASP-77Ab. This work shows that the planet is expected to have formed in a disk where not many planetesimals could be accreted. Moreover, we demonstrate that the most likely migration scenario is disk-free migration, whereby the planet initiates its Type II migration within the CO ice line and ends it beyond the water ice line. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 10 pages, 9 figures

Journal ref: volum number:675 Year: 2023 Pages:A95

SO$_2$, silicate clouds, but no CH$_4$ detected in a warm Neptune

Authors: Achrène Dyrek, Michiel Min, Leen Decin, Jeroen Bouwman, Nicolas Crouzet, Paul Mollière, Pierre-Olivier Lagage, Thomas Konings, Pascal Tremblin, Manuel Güdel, John Pye, Rens Waters, Thomas Henning, Bart Vandenbussche, Francisco Ardevol Martinez, Ioannis Argyriou, Elsa Ducrot, Linus Heinke, Gwenael Van Looveren, Olivier Absil, David Barrado, Pierre Baudoz, Anthony Boccaletti, Christophe Cossou, Alain Coulais , et al. (22 additional authors not shown)

Abstract: WASP-107b is a warm ($\sim$740 K) transiting planet with a Neptune-like mass of $\sim$30.5 $M_{\oplus}$ and Jupiter-like radius of $\sim$0.94 $R_{\rm J}$, whose extended atmosphere is eroding. Previous observations showed evidence for water vapour and a thick high-altitude condensate layer in WASP-107b's atmosphere. Recently, photochemically produced sulphur dioxide (SO$_2$) was detected in the at… ▽ More WASP-107b is a warm ($\sim$740 K) transiting planet with a Neptune-like mass of $\sim$30.5 $M_{\oplus}$ and Jupiter-like radius of $\sim$0.94 $R_{\rm J}$, whose extended atmosphere is eroding. Previous observations showed evidence for water vapour and a thick high-altitude condensate layer in WASP-107b's atmosphere. Recently, photochemically produced sulphur dioxide (SO$_2$) was detected in the atmosphere of a hot ($\sim$1,200 K) Saturn-mass planet from transmission spectroscopy near 4.05 $μ$m, but for temperatures below $\sim$1,000 K sulphur is predicted to preferably form sulphur allotropes instead of SO$_2$. Here we report the 9$σ$-detection of two fundamental vibration bands of SO$_2$, at 7.35 $μ$m and 8.69 $μ$m, in the transmission spectrum of WASP-107b using the Mid-Infrared Instrument (MIRI) of the JWST. This discovery establishes WASP-107b as the second irradiated exoplanet with confirmed photochemistry, extending the temperature range of exoplanets exhibiting detected photochemistry from $\sim$1,200 K down to $\sim$740 K. Additionally, our spectral analysis reveals the presence of silicate clouds, which are strongly favoured ($\sim$7$σ$) over simpler cloud setups. Furthermore, water is detected ($\sim$12$σ$), but methane is not. These findings provide evidence of disequilibrium chemistry and indicate a dynamically active atmosphere with a super-solar metallicity. △ Less

Submitted 21 November, 2023; originally announced November 2023.

15NH3 in the atmosphere of a cool brown dwarf

Authors: David Barrado, Paul Mollière, Polychronis Patapis, Michiel Min, Pascal Tremblin, Francisco Ardevol Martinez, Niall Whiteford, Malavika Vasist, Ioannis Argyriou, Matthias Samland, Pierre-Olivier Lagage, Leen Decin, Rens Waters, Thomas Henning, María Morales-Calderón, Manuel Guedel, Bart Vandenbussche, Olivier Absil, Pierre Baudoz, Anthony Boccaletti, Jeroen Bouwman, Christophe Cossou, Alain Coulais, Nicolas Crouzet, René Gastaud , et al. (18 additional authors not shown)

Abstract: Brown dwarfs serve as ideal laboratories for studying the atmospheres of giant exoplanets on wide orbits as the governing physical and chemical processes in them are nearly identical. Understanding the formation of gas giant planets is challenging, often involving the endeavour to link atmospheric abundance ratios, such as the carbon-to-oxygen (C/O) ratio, to formation scenarios. However, the comp… ▽ More Brown dwarfs serve as ideal laboratories for studying the atmospheres of giant exoplanets on wide orbits as the governing physical and chemical processes in them are nearly identical. Understanding the formation of gas giant planets is challenging, often involving the endeavour to link atmospheric abundance ratios, such as the carbon-to-oxygen (C/O) ratio, to formation scenarios. However, the complexity of planet formation requires additional tracers, as the unambiguous interpretation of the measured C/O ratio is fraught with complexity. Isotope ratios, such as deuterium-to-hydrogen and 14N/15N, offer a promising avenue to gain further insight into this formation process, mirroring their utility within the solar system. For exoplanets only a handful of constraints on 12C/13C exist, pointing to the accretion of 13C-rich ice from beyond the disks' CO iceline. Here we report on the mid-infrared detection of the 14NH3 and 15NH3 isotopologues in the atmosphere of a cool brown dwarf with an effective temperature of 380 K in a spectrum taken with the Mid-InfraRed Instrument of the James Webb Space Telescope. As expected, our results reveal a 14N/15N value consistent with star-like formation by gravitational collapse, demonstrating that this ratio can be accurately constrained. Since young stars and their planets should be more strongly enriched in the 15N isotope, we expect that 15NH3 will be detectable in a number of cold, wide-separation exoplanets. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Accepted by Nature. 28 pages, 7 figures, uses nature3.cls

The sulphur species in hot rocky exoplanet atmospheres

Authors: L. J. Janssen, P. Woitke, O. Herbort, M. Min, K. L. Chubb, Ch. Helling, L. Carone

Abstract: The first JWST observations of hot Jupiters showed an unexpected detection of SO2 in their hydrogen-rich atmospheres. We investigate how much sulphur can be expected in the atmospheres of rocky exoplanets and which sulphur molecules can be expected to be most abundant and detectable by transmission spectroscopy. We run thermo-chemical equilibrium models at the crust-atmosphere interface, consideri… ▽ More The first JWST observations of hot Jupiters showed an unexpected detection of SO2 in their hydrogen-rich atmospheres. We investigate how much sulphur can be expected in the atmospheres of rocky exoplanets and which sulphur molecules can be expected to be most abundant and detectable by transmission spectroscopy. We run thermo-chemical equilibrium models at the crust-atmosphere interface, considering surface temperatures 500 to 5000 K, surface pressures 1 to 100 bar, and various sets of element abundances based on common rock compositions. Between 1000 K and 2000 K, we find gaseous sulphur concentrations of up to 25 percent above the rock in our models. SO2, SO, H2S and S2 are by far the most abundant sulphur molecules. SO2 shows potentially detectable features in transmission spectra at about 4 micron, between 7 and 8 micron, and beyond 15 micron. In contrast, the sometimes abundant H2S molecule is difficult to detect in these spectra, which are mostly dominated by H2O and CO2. Although the molecule PS only occurs with concentrations below 300 ppm, it can cause a strong absorption feature between 0.3 and 0.65 micron in some of our models for high surface pressures. The detection of sulphur molecules would enable a better characterisation of the planetary surface. △ Less

Submitted 31 October, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Typos have been corrected compared to version 1

MSC Class: 85

Hydrocarbon chemistry in inner regions of planet forming disks

Authors: Jayatee Kanwar, Inga Kamp, Peter Woitke, Christian Rab, Wing-Fai Thi, Michiel Min

Abstract: The analysis of the mid-infrared spectra helps understanding the composition of the gas in the inner, dense and warm terrestrial planet forming region of disks around young stars. ALMA has detected hydrocarbons in the outer regions of the planet forming disk and Spitzer detected \ce{C2H2} in the inner regions. JWST- MIRI provides high spectral resolution observations of \ce{C2H2} and a suite of mo… ▽ More The analysis of the mid-infrared spectra helps understanding the composition of the gas in the inner, dense and warm terrestrial planet forming region of disks around young stars. ALMA has detected hydrocarbons in the outer regions of the planet forming disk and Spitzer detected \ce{C2H2} in the inner regions. JWST- MIRI provides high spectral resolution observations of \ce{C2H2} and a suite of more complex hydrocarbons are now reported. Interpreting the fluxes observed in the spectra is challenging and radiation thermo-chemical codes are needed to properly take into account the disk structure, radiative transfer, chemistry and thermal balance. Various disk physical parameters like the gas-to-dust ratio, dust evolution including radial drift, dust growth and settling can affect the fluxes observed in the mid-IR. Still, thermo-chemical disk models were not always successful in matching all observed molecular emission bands simultaneously. The goal of this project is two-fold. We analyse the warm carbon chemistry in the inner regions of the disk, i.e. within 10 au to find pathways forming \ce{C2H2} potentially missing from the existing chemical networks. Second, we analyse the effect of the new chemistry on the line fluxes of acetylene. We use radiative thermo-chemical disk code {P{\small RO}D{\small I}M{\small O}} to expand the hydrocarbon chemistry that occurs in a typical standard T Tauri disks. We used the UMIST and the KIDA rate databases for collecting reactions for the species. We include a number of three-body and thermal decomposition reactions from STAND2020 network. We included isotopomers for the species that were present in the databases. The chemistry is then analysed in the regions that produce observable features in the mid-infrared spectra. The effect of expanding the hydrocarbon chemistry on the mid-infrared spectra is studied. Acetylene is formed via two .... △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: accepted for publication in A&A

Journal ref: A&A 681, A22 (2024)

Analysing the SEDs of protoplanetary disks with machine learning

Authors: T. Kaeufer, P. Woitke, M. Min, I. Kamp, C. Pinte

Abstract: ABRIDGED. The analysis of spectral energy distributions (SEDs) of protoplanetary disks to determine their physical properties is known to be highly degenerate. Hence, a Bayesian analysis is required to obtain parameter uncertainties and degeneracies. The challenge here is computational speed, as one radiative transfer model requires a couple of minutes to compute. We performed a Bayesian analysis… ▽ More ABRIDGED. The analysis of spectral energy distributions (SEDs) of protoplanetary disks to determine their physical properties is known to be highly degenerate. Hence, a Bayesian analysis is required to obtain parameter uncertainties and degeneracies. The challenge here is computational speed, as one radiative transfer model requires a couple of minutes to compute. We performed a Bayesian analysis for 30 well-known protoplanetary disks to determine their physical disk properties, including uncertainties and degeneracies. To circumvent the computational cost problem, we created neural networks (NNs) to emulate the SED generation process. We created two sets of radiative transfer disk models to train and test two NNs that predict SEDs for continuous and discontinuous disks. A Bayesian analysis was then performed on 30 protoplanetary disks with SED data collected by the DIANA project to determine the posterior distributions of all parameters. We ran this analysis twice, (i) with old distances and additional parameter constraints as used in a previous study, to compare results, and (ii) with updated distances and free choice of parameters to obtain homogeneous and unbiased model parameters. We evaluated the uncertainties in the determination of physical disk parameters from SED analysis, and detected and quantified the strongest degeneracies. The NNs are able to predict SEDs within 1ms with uncertainties of about 5% compared to the true SEDs obtained by the radiative transfer code. We find parameter values and uncertainties that are significantly different from previous values obtained by $χ^2$ fitting. Comparing the global evidence for continuous and discontinuous disks, we find that 26 out of 30 objects are better described by disks that have two distinct radial zones. Also, we created an interactive tool that instantly returns the SED predicted by our NNs for any parameter combination. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 40 pages, 22 figures, the online tool is available at https://tillkaeufer.github.io/sedpredictor

Journal ref: A&A 672, A30 (2023)

What does a typical full-disc around a post-AGB binary look like? -- Radiative transfer models reproducing PIONIER, GRAVITY, and MATISSE data

Authors: A. Corporaal, J. Kluska, H. Van Winckel, D. Kamath, M. Min

Abstract: (abridged) Stable circumbinary discs around evolved post-Asymptotic Giant branch (post-AGB) binary systems show many similarities with protoplanetary discs around young stellar objects. These discs can provide constraints on both binary evolution and the formation of macrostructures within circumstellar discs. Here we focus on one post-AGB binary system: IRAS08544-4431. We aim to refine the physic… ▽ More (abridged) Stable circumbinary discs around evolved post-Asymptotic Giant branch (post-AGB) binary systems show many similarities with protoplanetary discs around young stellar objects. These discs can provide constraints on both binary evolution and the formation of macrostructures within circumstellar discs. Here we focus on one post-AGB binary system: IRAS08544-4431. We aim to refine the physical model of IRAS08544-4431 with a radiative transfer treatment and continue the near-infrared and mid-infrared interferometric analysis covering the H-, K-, L-, and N-bands. We aim to capture the previously detected amount of over-resolved flux and the radial intensity profile at and beyond the inner dust disc rim to put constraints on the physical processes in the inner disc regions. We used a Monte Carlo radiative transfer code to investigate the physical structure of the disc by reproducing both the photometry and the multi-wavelength infrared interferometric data set. We developed a strategy to identify the models which perform best to reproduce our data set. We found a family of models that successfully fit the infrared photometric and interferometric data in all bands. Some over-resolved flux component was recovered in all bands but the optimised models still fall short to explain all the over-resolved flux. This suggests that another dusty structure within the system plays a role. Multi-wavelength infrared interferometric observations of circumstellar discs allow to study the inner disc regions in unprecedented detail. The refined physical models can reproduce most of the investigated features, including the photometric characteristics, the radial extent, and the overall shape of the visibility curves. Our multi-wavelength interferometric observations combined with photometry show that the disc is similar to protoplanetary discs with similar dust masses and efficient dust growth. △ Less

Submitted 12 January, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 18 pages, 13 figures (including apppendix), accepted for publication in A&A

Journal ref: A&A 671, A15 (2023)

Clouds form on the hot Saturn JWST ERO target WASP-96b

Authors: Dominic Samra, Christiane Helling, Katy Chubb, Michiel Min, Ludmila Carone, Aaron Schneider

Abstract: WASP-96b is a hot Saturn exoplanet, with an equilibrium temperature well within the regime of thermodynamically expected extensive cloud formation. Prior observations with Hubble/WFC3, Spitzer/IRAC, and VLT/FORS2 have been combined into a single spectra for which retrievals suggest a cold but cloud-free atmosphere. Recently, the planet was observed with the James Webb Space Telescope (JWST) as par… ▽ More WASP-96b is a hot Saturn exoplanet, with an equilibrium temperature well within the regime of thermodynamically expected extensive cloud formation. Prior observations with Hubble/WFC3, Spitzer/IRAC, and VLT/FORS2 have been combined into a single spectra for which retrievals suggest a cold but cloud-free atmosphere. Recently, the planet was observed with the James Webb Space Telescope (JWST) as part of the Early Release Observations (ERO). 1D profiles are extracted from the 3D GCM expeRT/MITgcm results and used as input for a kinetic, non-equilibrium model to study the formation of mineral cloud particles of mixed composition. The ARCiS retrieval framework is applied to the pre-JWST WASP-96b transit spectra to investigate the apparent contradiction between cloudy models and assumed cloud-free transit spectra. Clouds are predicted to be ubiquitous throughout the atmosphere of WASP-96b. Silicate materials contribute between 40% and 90%, hence, also metal oxides contribute with up to 40% in the low-pressure regimes that effect the spectra. We explore how to match these cloudy models with currently available atmospheric transit spectra. A reduced vertical mixing acts to settle clouds to deeper in the atmosphere, and an increased cloud particles porosity reduces the opacity of clouds in the near-IR and optical region. Both these effects allow for clearer molecular features to be observed, while still allowing clouds to be in the atmosphere. The atmosphere of WASP-96b is unlikely to be cloud free. Also retrievals of HST, Spitzer and VLT spectra show that multiple cloudy solutions reproduce the data. JWST observations will be affected by clouds, where within even the NIRISS wavelength range the cloud top pressure varies by an order of magnitude. The long wavelength end of NIRSpec and short end of MIRI may probe atmospheric asymmetries between the limbs of the terminator on WASP-96b. △ Less

Submitted 1 November, 2022; originally announced November 2022.

Comments: 17 pages, 13 figures, accepted for publication in A&A

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

H2S and SO2 detectability in Hot Jupiters: Sulfur species as indicator of metallicity and C/O ratio

Authors: J. Polman, L. B. F. M. Waters, M. Min, Y. Miguel, N. Khorshid

Abstract: The high cosmic abundance and the intermediate volatility and chemical properties of sulfur allow the use of sulfur-bearing species as a tracer of the chemical processes in the atmospheres of hot Jupiter exoplanets. Nevertheless, despite its properties and relevance as a tracer of the giant planets' formation history, little attention has been paid to this species in the context of hot Jupiter atm… ▽ More The high cosmic abundance and the intermediate volatility and chemical properties of sulfur allow the use of sulfur-bearing species as a tracer of the chemical processes in the atmospheres of hot Jupiter exoplanets. Nevertheless, despite its properties and relevance as a tracer of the giant planets' formation history, little attention has been paid to this species in the context of hot Jupiter atmospheres. Here we provide an overview of the abundances of sulfur-bearing species in hot Jupiter atmospheres under different conditions and explore their observability. We use the photochemical kinetics code VULCAN to model hot Jupiter atmospheric disequilibrium chemistry. Transmission spectra for these atmospheres are created using the modelling framework ARCiS. We vary model parameters such as the diffusion coefficient, and we study the importance of photochemistry on the resulting mixing ratios. Furthermore, we vary the chemical composition of the atmosphere by increasing the metallicity from solar to ~10 times solar. We also explore different C/O ratios. We find that H2S and SO2 are the best candidates for detection between 1 and 10 micron, using a spectral resolution that is representative of the instruments on board the JWST. H2S is easiest to detect at an equilibrium temperature of ~1500 K and C/O ratios between 0.7 and 0.9, with the ideal value increasing slightly for increasing metallicity. SO2 is most likely to be detected at an equilibrium temperature of ~1000 K at low C/O ratios and high metallicities. Nevertheless, among these two molecules, we expect SO2 detection to be more common, as it is detectable in scenarios more favoured by formation models. We conclude that H2S and SO2 will most likely be detected in the coming years with the JWST and that the detection of these species will provide information on atmospheric processes and planet formation scenarios. △ Less

Submitted 16 November, 2022; v1 submitted 31 July, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A

Exoplanet atmosphere retrievals in 3D using phase curve data with ARCiS: application to WASP-43b

Authors: Katy L. Chubb, Michiel Min

Abstract: Aims. To create a retrieval framework which encapsulates the 3D nature of exoplanet atmospheres, and to apply it to observed emission phase curve and transmission spectra of hot Jupiter exoplanet WASP-43b. Methods. We present our 3D framework, which is freely available as a stand-alone module from GitHub. We use the atmospheric modelling and Bayesian retrieval package ARCiS (ARtful modelling Code… ▽ More Aims. To create a retrieval framework which encapsulates the 3D nature of exoplanet atmospheres, and to apply it to observed emission phase curve and transmission spectra of hot Jupiter exoplanet WASP-43b. Methods. We present our 3D framework, which is freely available as a stand-alone module from GitHub. We use the atmospheric modelling and Bayesian retrieval package ARCiS (ARtful modelling Code for exoplanet Science) to perform 8 3D retrievals on simultaneous transmission (HST/WFC3) and phase-dependent emission (HST/WFC3, Spitzer/IRAC) observations of WASP-43b as a case study. We assess how input assumptions affect our retrieval outcomes. In particular we look at constraining equilibrium chemistry vs a free molecular retrieval, the case of no clouds vs parametrised clouds, and using Spitzer phase data that have been reduced from two different literature sources. For the free chemistry retrievals, we retrieve abundances of H2O, CH4, CO, CO2, AlO, and NH3 as a function of phase, with many more species considered for the equilibrium chemistry retrievals. Results. We find consistent super-solar C/O (0.6-0.9) and super-solar metallicities (1.7-2.9 dex) for all retrieval setups that assume equilibrium chemistry. Atmospheric heat distribution, hotspot shift (15.6 vs 4.5 for different Spitzer datasets), and temperature structure are influenced by the Spitzer emission phase data. Comparisons are made with other studies of WASP-43b, including available GCM simulations. Conclusions. The parametrised 3D setup we have developed provides a valuable tool to analyse extensive observational datasets such as spectroscopic phase curves. Near-future observations with missions such as the James Webb Space Telescope (JWST) and Ariel will greatly improve our understanding of the atmospheres of exoplanets such as WASP-43b. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 31 pages; published in A&A

Journal ref: A&A 665, A2 (2022)

Convolutional neural networks as an alternative to Bayesian retrievals

Authors: Francisco Ardevol Martinez, Michiel Min, Inga Kamp, Paul I. Palmer

Abstract: Exoplanet observations are currently analysed with Bayesian retrieval techniques. Due to the computational load of the models used, a compromise is needed between model complexity and computing time. Analysis of data from future facilities, will need more complex models which will increase the computational load of retrievals, prompting the search for a faster approach for interpreting exoplanet o… ▽ More Exoplanet observations are currently analysed with Bayesian retrieval techniques. Due to the computational load of the models used, a compromise is needed between model complexity and computing time. Analysis of data from future facilities, will need more complex models which will increase the computational load of retrievals, prompting the search for a faster approach for interpreting exoplanet observations. Our goal is to compare machine learning retrievals of exoplanet transmission spectra with nested sampling, and understand if machine learning can be as reliable as Bayesian retrievals for a statistically significant sample of spectra while being orders of magnitude faster. We generate grids of synthetic transmission spectra and their corresponding planetary and atmospheric parameters, one using free chemistry models, and the other using equilibrium chemistry models. Each grid is subsequently rebinned to simulate both HST/WFC3 and JWST/NIRSpec observations, yielding four datasets in total. Convolutional neural networks (CNNs) are trained with each of the datasets. We perform retrievals on a 1,000 simulated observations for each combination of model type and instrument with nested sampling and machine learning. We also use both methods to perform retrievals on real WFC3 transmission spectra. Finally, we test how robust machine learning and nested sampling are against incorrect assumptions in our models. CNNs reach a lower coefficient of determination between predicted and true values of the parameters. Nested sampling underestimates the uncertainty in ~8% of retrievals, whereas CNNs estimate them correctly. For real WFC3 observations, nested sampling and machine learning agree within $2σ$ for ~86% of spectra. When doing retrievals with incorrect assumptions, nested sampling underestimates the uncertainty in ~12% to ~41% of cases, whereas this is always below ~10% for the CNN. △ Less

Submitted 3 March, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 662, A108 (2022)

A retrieval challenge exercise for the Ariel mission

Authors: Joanna K. Barstow, Quentin Changeat, Katy L. Chubb, Patricio E. Cubillos, Billy Edwards, Ryan J. MacDonald, Michiel Min, Ingo P. Waldmann

Abstract: The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a c… ▽ More The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a critical preparatory step for the mission. Here, we present the results of a retrieval challenge. We use five different exoplanet retrieval codes to analyse the same synthetic datasets, and test a) the ability of each to recover the correct input solution and b) the consistency of the results. We find that generally there is very good agreement between the five codes, and in the majority of cases the correct solutions are recovered. This demonstrates the reproducibility of retrievals for transit spectra of exoplanets, even when codes are not previously benchmarked against each other. △ Less

Submitted 15 March, 2022; v1 submitted 1 March, 2022; originally announced March 2022.

Comments: 28 pages, 14 figures. Accepted in Experimental Astronomy (2022)

Investigating the detectability of hydrocarbons in exoplanet atmospheres with JWST

Authors: Danny Gasman, Michiel Min, Katy L. Chubb

Abstract: We investigate at what abundances various hydrocarbon molecules (e.g. acetylene (C$_2$H$_2$), ethylene (C$_2$H$_4$), and methane (CH$_4$)) become detectable when observing the atmospheres of various planets using JWST. We focused on atmospheric models based on the parameters of a small sample of planets: HD 189733b, HD 209458b, HD 97658b, and Kepler-30c. We computed model transmission spectra usin… ▽ More We investigate at what abundances various hydrocarbon molecules (e.g. acetylene (C$_2$H$_2$), ethylene (C$_2$H$_4$), and methane (CH$_4$)) become detectable when observing the atmospheres of various planets using JWST. We focused on atmospheric models based on the parameters of a small sample of planets: HD 189733b, HD 209458b, HD 97658b, and Kepler-30c. We computed model transmission spectra using the Bayesian retrieval package ARCiS. We simulated observed spectra using the PandExo package. We subsequently ran retrievals on these spectra to determine whether the molecular abundances can be accurately retrieved from these simulated spectra. We find that generally we can detect and retrieve abundances of the hydrocarbon species as long as they have a VMR above approximately 1x10$^{-7}$-1x10$^{-6}$, at least for the brighter targets. There are variations based on planet type and instrument(s) used, and these limits will likely change depending on the abundance of other strong absorbers. We also find scenarios where the presence of one hydrocarbon is confused with another; this is often improved when two instruments are combined. C$_2$H$_2$, CH$_4$, and C$_2$H$_4$ will all be detectable with JWST, provided they are present in high enough abundances, and that the optimal instruments are chosen. A combination of two instruments, either NIRSpec G395M and MIRI LRS, or NIRCam F322W2 and MIRI LRS, is best for observing these species in bright exoplanet systems, with NIRSpec G395M and MIRI LRS the best option for the HD 189733b-like atmosphere with clouds included. The use of NIRSpec Prism is tentatively found to be best for fainter targets, potentially in combination with the MIRI LRS slit mode, although the target we test is too faint to draw any strong conclusions. Instrument sensitivity, noise, and wavelength range are all thought to play a role in being able to distinguish spectral features. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: accepted A&A

A population of transition disks around evolved stars: Fingerprints of planets? Catalog of disks surrounding Galactic post-AGB binaries

Authors: Jacques Kluska, Hans Van Winckel, Quentin Coppée, Glenn-Michael Oomen, Devika Kamath, Valentin Bujarrabal, Michiel Min

Abstract: Post-asymptotic giant branch (post-AGB) binaries are surrounded by massive disks of gas and dust that are similar to protoplanetary disks around young stars. We assembled a catalog of all known Galactic post-AGB binaries with disks. We explore the correlations between the different observables with the aim of learning more about potential disk-binary interactions. We compiled spectral energy distr… ▽ More Post-asymptotic giant branch (post-AGB) binaries are surrounded by massive disks of gas and dust that are similar to protoplanetary disks around young stars. We assembled a catalog of all known Galactic post-AGB binaries with disks. We explore the correlations between the different observables with the aim of learning more about potential disk-binary interactions. We compiled spectral energy distributions of 85 Galactic post-AGB binary systems. We built a color-color diagram to differentiate between the different disk morphologies traced by the infrared excess. We categorized the different disk types and searched for correlations with other observational characteristics of these systems. Between 8 and 12% of our targets are surrounded by transition disks, that is, disks having no or low near-infrared excess. We find a strong link between these transition disks and the depletion of refractory elements on the surface of the post-AGB star. We interpret this link as evidence of a mechanism that stimulates the dust and gas separation within the disk and that also produces the transition disk structure. We propose that such a mechanism is likely to be due to a giant planet carving a hole in the disk, effectively trapping the dust in the outer disk parts. We propose disk evolutionary scenarios, depending on the presence of such a giant planet in the disk. We advocate that giant planets can successfully explain the link between the transition disks and the observed depletion of refractory materials. If the planetary scenario is confirmed, disks around post-AGB binaries could be a unique laboratory for testing planet-disk interactions and their influence on the late evolution of binary stars. The question of whether such planets are first- or second-generation bodies remains to be considered. We argue that these disks are ideal for studying planet formation in unprecedented parameter space. △ Less

Submitted 31 January, 2022; originally announced January 2022.

Comments: Published in A&A. 27 pages, 25 figures

Journal ref: A&A 658, A36 (2022)

Self-consistent modelling of the dust component in protoplanetary and circumplanetary disks: the case of PDS 70

Authors: B. Portilla-Revelo, I. Kamp, Ch. Rab, E. F. van Dishoeck, M. Keppler, M. Min, G. A. Muro-Arena

Abstract: Direct observations of young stellar objects are important to test established theories of planet formation. PDS 70 is one of the few cases where robust evidence favours the presence of two planetary mass companions inside the gap of the transition disk. Those planets are believed to be going through the last stages of accretion from the protoplanetary disk, a process likely mediated by a circumpl… ▽ More Direct observations of young stellar objects are important to test established theories of planet formation. PDS 70 is one of the few cases where robust evidence favours the presence of two planetary mass companions inside the gap of the transition disk. Those planets are believed to be going through the last stages of accretion from the protoplanetary disk, a process likely mediated by a circumplanetary disk (CPD). We aim to develop a three dimensional radiative transfer model for the dust component of the PDS 70 system which reproduces the system's global features observed at two different wavelengths: 855 $μ\, \mathrm{m}$ with ALMA and 1.25 $μ\, \mathrm{m}$ with VLT/SPHERE. We use this model to investigate the physical properties of the planetary companion PDS 70 c and its potential circumplanetary disk. We select initial values for the physical properties of the planet and CPD through appropriate assumptions about the nature and evolutionary stage of the object. We modify iteratively the properties of the protoplanetary disk until the predictions retrieved from the model are consistent with both data sets. We provide a model that jointly explains the global features of the PDS 70 system seen in submillimeter and polarised-scattered light. Our model suggests that spatial segregation of dust grains is present in the protoplanetary disk. The submillimeter modelling of the PDS 70 c source favours the presence of an optically thick CPD and places an upper limit to its dust mass of 0.7 $M_\oplus$. Furthermore, analysis of the thermal structure of the CPD demonstrates that the planet luminosity is the dominant heating mechanism of dust grains inside 0.6 au from the planet while heating by stellar photons dominates at larger planetocentric distances. △ Less

Submitted 25 April, 2022; v1 submitted 16 November, 2021; originally announced November 2021.

Comments: accepted for publication in A&A

Journal ref: A&A 658, A89 (2022)

SimAb: A simple, fast and flexible model to assess the effects of planet formation on the atmospheric composition of gas giants

Authors: N. Khorshid, M. Min, J. M. Désert, P. Woitke, C. Dominik

Abstract: We present a basic, fast, and flexible planet formation model, called SimAb (Simulating Abundances), to form giant planets and study their primary atmospheric composition soon after their formation. In SimAb we introduce parameters to simplify the assumptions about the complex physics involved in the formation of a planet. This approach allows us to trace and understand the influence of complex ph… ▽ More We present a basic, fast, and flexible planet formation model, called SimAb (Simulating Abundances), to form giant planets and study their primary atmospheric composition soon after their formation. In SimAb we introduce parameters to simplify the assumptions about the complex physics involved in the formation of a planet. This approach allows us to trace and understand the influence of complex physical processes on the formed planets. We focus on the C/O ratio and the metallicity of the planetary atmosphere as an indicator of their compositions. We show that the initial protoplanet core mass does not influence the final composition of the planetary atmosphere in the context of our model. The initial orbital distance affects the C/O ratio due to the different C/O ratios in the gas phase and the solid phase at different orbital distances. Additionally, the initial orbital distance together with the amount of accreted planetesimals cause the planet to have sub-solar or super-solar metallicity. Furthermore, the C/O ratio is affected by the dust grain fraction and the planetesimal fraction. Planets that accrete most of their heavy elements through dust grains will have a C/O ratio close to the solar C/O ratio, while planets that accrete most of their heavy elements from the planetesimals in the disk will end up with a C/O ratio closer to the C/O ratio in the solid phase of the disk. By using the C/O ratio and metallicity together we can put a lower and upper boundary on the initial orbital distance where super-solar metallicity planets are formed. We show that planetesimals are the main source for reaching super-solar metallicity planets. On the other hand, planets that mainly accrete dust grains will show a more solar composition. Super-solar metallicity planets that initiate their formation farther than the CO ice line have a C/O ratio closer to the solar value. △ Less

Submitted 30 October, 2021; originally announced November 2021.

Comments: 14 pages,11 figures, Submitted to A&A

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

Implementation of disequilibrium chemistry to spectral retrieval code ARCiS and application to sixteen exoplanet transmission spectra: Indication of disequilibrium chemistry for HD 209458b and WASP-39b

Authors: Yui Kawashima, Michiel Min

Abstract: The retrieval approach is currently a standard method for deriving atmospheric properties from observed spectra of exoplanets. However, the approach ignores disequilibrium chemistry in most current retrieval codes, which can lead to misinterpretation of the metallicity or elemental abundance ratios of the atmosphere. We have implemented the disequilibrium effect of vertical mixing or quenching for… ▽ More The retrieval approach is currently a standard method for deriving atmospheric properties from observed spectra of exoplanets. However, the approach ignores disequilibrium chemistry in most current retrieval codes, which can lead to misinterpretation of the metallicity or elemental abundance ratios of the atmosphere. We have implemented the disequilibrium effect of vertical mixing or quenching for the major species in hydrogen/helium-dominated atmospheres, namely $\mathrm{CH_4}$, $\mathrm{CO}$, $\mathrm{H_2O}$, $\mathrm{NH_3}$, $\mathrm{N_2}$, and $\mathrm{CO_2}$, for the spectral retrieval code ARCiS with a physical basis using the chemical relaxation method. Then, using ARCiS updated with this module, we have performed retrievals of the observed transmission spectra of 16 exoplanets with sizes ranging from Jupiter to mini-Neptune. As a result, we find indications of disequilibrium chemistry for HD 209458b ($\geq 4.1σ$) and WASP-39b ($\geq 2.7σ$). The retrieved spectrum of HD 209458b exhibits a strong $\mathrm{NH_3}$ absorption feature at 10.5 $μ$m accessible by JWST owing to an enhanced abundance of $\mathrm{NH_3}$ due to the quenching effect. This feature is absent in the spectrum retrieved assuming equilibrium chemistry, which makes HD 209458b an ideal target for studying disequilibrium chemistry in exoplanet atmospheres. Moreover, for HAT-P-11b and GJ 436b, we obtain relatively different results than for the retrieval with the equilibrium assumption, such as a $2.9σ$ difference for the C/O ratio. We have also examined the retrieved eddy diffusion coefficient, but could not identify a trend over the equilibrium temperature, possibly due to the limits of the current observational precision. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 37 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 656, A90 (2021)

Multi-wavelength VLTI study of the puffed-up inner rim of a circumbinary disc

Authors: A. Corporaal, J. Kluska, H. Van Winckel, D. Bollen, D. Kamath, M. Min

Abstract: The presence of stable, compact circumbinary discs of gas and dust around post-asymptotic giant branch (post-AGB) binary systems has been well established. We focus on one such system: IRAS 08544-4431. We present an interferometric multi-wavelength analysis of the circumstellar environment of IRAS 08544-4431. The aim is to constrain different contributions to the total flux in the H, K, L, and N-b… ▽ More The presence of stable, compact circumbinary discs of gas and dust around post-asymptotic giant branch (post-AGB) binary systems has been well established. We focus on one such system: IRAS 08544-4431. We present an interferometric multi-wavelength analysis of the circumstellar environment of IRAS 08544-4431. The aim is to constrain different contributions to the total flux in the H, K, L, and N-bands in the radial direction. The data from VLTI/PIONIER, VLTI/GRAVITY, and VLTI/MATISSE range from the near-infrared, where the post-AGB star dominates, to the mid-infrared, where the disc dominates. We fitted two geometric models to the visibility data to reproduce the circumbinary disc: a ring with a Gaussian width and a flat disc model with a temperature gradient. The flux contributions from the disc, the primary star (modelled as a point-source), and an over-resolved component are recovered along with the radial size of the emission, the temperature of the disc as a function of radius, and the spectral dependencies of the different components. The trends of all visibility data were well reproduced with the geometric models. The near-infrared data were best fitted with a Gaussian ring model while the mid-infrared data favoured a temperature gradient model. This implies that a vertical structure is present at the disc inner rim, which we attribute to a rounded puffed-up inner rim. The N-to-K size ratio is 2.8, referring to a continuous flat source, analogues to young stellar objects. By combining optical interferometric instruments operating at different wavelengths we can resolve the complex structure of circumstellar discs and study the wavelength-dependent opacity profile. A detailed radial, vertical, and azimuthal structural analysis awaits a radiative transfer treatment in 3D to capture all non-radial complexity. △ Less

Submitted 11 June, 2021; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 8 pages, 5 figures, accepted for publication in A&A Letters

Journal ref: A&A 650, L13 (2021)

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

Cloud property trends in hot and ultra-hot giant gas planets (WASP-43b, WASP-103b, WASP-121b, HAT-P-7b, and WASP-18b)

Authors: Ch. Helling, D. Lewis, D. Samra, L. Carone, V. Graham, O. Herbort, K. L. Chubb, M. Min, R. Waters, V. Parmentier, N. Mayne

Abstract: Ultra-hot Jupiters are the hottest exoplanets discovered so far. Observations begin to provide insight into the composition of their extended atmospheres and their chemical day/night asymmetries. Both are strongly affected by cloud formation. We explore trends in cloud properties for a sample of five giant gas planets: WASP-43b, WASP-18b, HAT-P-7b, WASP-103b, and WASP-121b. This provides a referen… ▽ More Ultra-hot Jupiters are the hottest exoplanets discovered so far. Observations begin to provide insight into the composition of their extended atmospheres and their chemical day/night asymmetries. Both are strongly affected by cloud formation. We explore trends in cloud properties for a sample of five giant gas planets: WASP-43b, WASP-18b, HAT-P-7b, WASP-103b, and WASP-121b. This provides a reference frame for cloud properties for the JWST targets WASP-43b and WASP-121b. We further explore chemically inert tracers to observe geometrical asymmetries, and if the location of inner boundary of a 3D GCM matters for the clouds that form. The large day/night temperature differences of ultra-hot Jupiters cause large chemical asymmetries: cloud-free days but cloudy nights, atomic vs. molecular gases and respectively different mean molecular weights, deep thermal ionospheres vs. low-ionised atmospheres, undepleted vs enhanced C/O. WASP-18b, as the heaviest planet in the sample, has the lowest global C/O. The global climate may be considered as similar amongst ultra-hot Jupiters, but different to that of hot gas giants. The local weather, however, is individual for each planet since the local thermodynamic conditions, and hence the local cloud and gas properties, differ. The morning and the evening terminator of ultra-hot Jupiters will carry signatures of their strong chemical asymmetry such that ingress/egress asymmetries can be expected. An increased C/O ratio is a clear sign of cloud formation, making cloud modelling a necessity when utilizing C/O (or other mineral ratios) as tracer for planet formation. The changing geometrical extension of the atmosphere from the day to the nightside may be probed through chemically inert species like helium. Ultra-hot Jupiters are likely to develop deep atmospheric ionospheres which may impact the atmosphere dynamics through MHD processes. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 38 pages, accepted for publication in A&A

Journal ref: A&A 649, A44 (2021)

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

Authors: E. Cannon, M. Montarges, A. de Koter, L. Decin, M. Min, E. Lagadec, P. Kervella, J. O. Sundqvist, H. Sana

Abstract: The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHER… ▽ More The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: 15 pages, accepted for publication in MNRAS

VLTI images of circumbinary disks around evolved stars

Authors: Jacques Kluska, Rik Claes, Akke Corporaal, Hans Van Winckel, Javier Alcolea, Narsireddy Anugu, Jean-Philippe Berger, Dylan Bollen, Valentin Bujarrabal, Robert Izzard, Devika Kamath, Stefan Kraus, Jean-Baptiste Le Bouquin, Michiel Min, John D. Monnier, Hans Olofsson

Abstract: The new generation of VLTI instruments (GRAVITY, MATISSE) aims to produce routinely interferometric images to uncover the morphological complexity of different objects at high angular resolution. Image reconstruction is, however, not a fully automated process. Here we focus on a specific science case, namely the complex circumbinary environments of a subset of evolved binaries, for which interfero… ▽ More The new generation of VLTI instruments (GRAVITY, MATISSE) aims to produce routinely interferometric images to uncover the morphological complexity of different objects at high angular resolution. Image reconstruction is, however, not a fully automated process. Here we focus on a specific science case, namely the complex circumbinary environments of a subset of evolved binaries, for which interferometric imaging provides the spatial resolution required to resolve the immediate circumbinary environment. Indeed, many binaries where the main star is in the post-asymptotic giant branch (post-AGB) phase are surrounded by circumbinary disks. Those disks were first inferred from the infrared excess produced by dust. Snapshot interferometric observations in the infrared confirmed disk-like morphology and revealed high spatial complexity of the emission that the use of geometrical models could not recover without being strongly biased. Arguably, the most convincing proof of the disk-like shape of the circumbinary environment came from the first interferometric image of such a system (IRAS08544-4431) using the PIONIER instrument at the VLTI. This image was obtained using the SPARCO image reconstruction approach that enables to subtract a model of a component of the image and reconstruct an image of its environment only. In the case of IRAS08544-4431, the model involved a binary and the image of the remaining signal revealed several unexpected features. Then, a second image revealed a different but also complex circumstellar morphology around HD101584 that was well studied by ALMA. To exploit the VLTI imaging capability to understand these targets, we started a large program at the VLTI to image post-AGB binary systems using both PIONIER and GRAVITY instruments. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 9 pages, Draft of the SPIE Proceedings Volume 11446, Optical and Infrared Interferometry and Imaging VII; 114460D (2020)

Neural network based image reconstruction with astrophysical priors

Authors: R. Claes, J. Kluska, H. Van Winckel, M. Min

Abstract: With the advent of interferometric instruments with 4 telescopes at the VLTI and 6 telescopes at CHARA, the scientific possibility arose to routinely obtain milli-arcsecond scale images of the observed targets. Such an image reconstruction process is typically performed in a Bayesian framework where the function to minimize is made of two terms: the datalikelihood and the Bayesian prior. This prio… ▽ More With the advent of interferometric instruments with 4 telescopes at the VLTI and 6 telescopes at CHARA, the scientific possibility arose to routinely obtain milli-arcsecond scale images of the observed targets. Such an image reconstruction process is typically performed in a Bayesian framework where the function to minimize is made of two terms: the datalikelihood and the Bayesian prior. This prior should be based on our prior knowledge of the observed source. Up to now,this prior was chosen from a set of generic and arbitrary functions, such as total variation for example. Here, we present an image reconstruction framework using generative adversarial networks where the Bayesian prior is defined using state-of-the-art radiative transfer models of the targeted objects. We validate this new image reconstruction algorithm on synthetic data with added noise. The generated images display a drastic reduction of artefacts and allow a more straight forward astrophysical interpretation. The results can be seen as a first illustration of how neural networks can provide significant improvements to the image reconstruction of a variety of astrophysical sources. △ Less

Submitted 18 December, 2020; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: SPIE Astronomical Telescopes + Instrumentation 2020 conference Paper No. 11446-110

Journal ref: Proc. SPIE 11446, Optical and Infrared Interferometry and Imaging VII, 114461U (13 December 2020)

Evidence for chromium hydride in the atmosphere of Hot Jupiter WASP-31b

Authors: Marrick Braam, Floris F. S. van der Tak, Katy L. Chubb, Michiel Min

Abstract: Context. The characterisation of exoplanet atmospheres has shown a wide diversity in compositions. Hot Jupiters have the appropriate temperatures to host metallic compounds, which should be detectable through transmission spectroscopy. Aims. We aim to detect exotic species in transmission spectra of hot Jupiters, specifically WASP-31b, by testing a variety of chemical species to explain the spectr… ▽ More Context. The characterisation of exoplanet atmospheres has shown a wide diversity in compositions. Hot Jupiters have the appropriate temperatures to host metallic compounds, which should be detectable through transmission spectroscopy. Aims. We aim to detect exotic species in transmission spectra of hot Jupiters, specifically WASP-31b, by testing a variety of chemical species to explain the spectrum. Methods. We conduct a re-analysis of publicly available transmission data of WASP-31b using the Bayesian retrieval framework TauREx II. We retrieve various combinations of the opacities of 25 species to determine the minimum set needed to fit the observed spectrum. Results. We report evidence for the spectroscopic signatures of chromium hydride (CrH), H2O and K in WASP-31b. Compared to a flat model without any signatures, a CrH-only model is preferred with a statistical significance of ~3.9 sigma. A model consisting of both CrH+H2O is found with ~2.6 sigma and ~3 sigma confidence over a CrH-only and H2O-only model respectively. Weak evidence for the addition of K is found at ~2.2 sigma over the H2O+CrH model, although the fidelity of the data point associated with this signature was questioned in earlier studies. Finally, the inclusion of collision-induced absorption and a Rayleigh scattering slope is found with ~3.5 sigma confidence over a flat model. This analysis presents the first evidence for signatures of CrH in a hot Jupiter atmosphere. At a retrieved temperature of ~1481 K, the atmosphere of WASP-31b is hot enough to host gaseous Cr-bearing species. The retrieved abundances agree well with predictions from thermal equilibrium chemistry and with the abundance found in an L-type brown dwarf atmosphere. However, additional retrievals using VLT FORS2 data lead to a non-detection of CrH. Future coverage with JWST has the potential to confirm the detection and discover other CrH features. △ Less

Submitted 22 December, 2020; v1 submitted 20 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A (Nov 24, 2020). 12 pages, 6 figures, 9 tables

Journal ref: A&A 646, A17 (2021)

The ExoMolOP Database: Cross-sections and k-tables for Molecules of Interest in High-Temperature Exoplanet Atmospheres

Authors: Katy L. Chubb, Marco Rocchetto, Sergei N. Yurchenko, Michiel Min, Ingo Waldmann, Joanna K. Barstow, Paul Mollière, Ahmed F. Al-Refaie, Mark W. Phillips, Jonathan Tennyson

Abstract: A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures… ▽ More A publicly available database of opacities for molecules of astrophysical interest, ExoMolOP, has been compiled for over 80 species, based on the latest line list data from the ExoMol, HITEMP and MoLLIST databases. These data are generally suitable for characterising high temperature exoplanet or cool stellar/substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room-temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS and petitRADTRANS, and include both cross-sections (at R~=~$\fracλ{Δλ}$~=~15,000) and k-tables (at R~=~$\fracλ{Δλ}$~=~1000) for the 0.3~-~50$μ$m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or from those for a known molecule with similar molecular properties where no broadening data are available. The data are available from www.exomol.com. △ Less

Submitted 1 September, 2020; originally announced September 2020.

Comments: 25 pages, 9 figures. Accepted for publication in A&A

Journal ref: A&A 646, A21 (2021)

The ARCiS framework for Exoplanet Atmospheres: Modelling Philosophy and Retrieval

Authors: Michiel Min, Chris W. Ormel, Katy Chubb, Christiane Helling, Yui Kawashima

Abstract: Aims: ARCiS, a novel code for the analysis of exoplanet transmission and emission spectra is presented. The aim of the modelling framework is to provide a tool able to link observations to physical models of exoplanet atmospheres. Methods: The modelling philosophy chosen in this paper is to use physical and chemical models to constrain certain parameters while keeping free the parts where our phys… ▽ More Aims: ARCiS, a novel code for the analysis of exoplanet transmission and emission spectra is presented. The aim of the modelling framework is to provide a tool able to link observations to physical models of exoplanet atmospheres. Methods: The modelling philosophy chosen in this paper is to use physical and chemical models to constrain certain parameters while keeping free the parts where our physical understanding is still more limited. This approach, in between full physical modelling and full parameterisation, allows us to use the processes we understand well and parameterise those less understood. A Bayesian retrieval framework is implemented and applied to the transit spectra of a set of 10 hot Jupiters. The code contains chemistry and cloud formation and has the option for self consistent temperature structure computations. Results: The code presented is fast and flexible enough to be used for retrieval and for target list simulations for e.g. JWST or the ESA Ariel missions. We present results for the retrieval of elemental abundance ratios using the physical retrieval framework and compare this to results obtained using a parameterised retrieval setup. Conclusions: We conclude that for most of the targets considered the current dataset is not constraining enough to reliably pin down the elemental abundance ratios. We find no significant correlations between different physical parameters. We confirm that planets in our sample with a strong slope in the optical transmission spectrum are the planets where we find cloud formation to be most active. Finally, we conclude that with ARCiS we have a computationally efficient tool to analyse exoplanet observations in the context of physical and chemical models. △ Less

Submitted 28 July, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 642, A28 (2020)

Mineral cloud and hydrocarbon haze particles in the atmosphere of the hot Jupiter JWST target WASP-43b

Authors: Ch. Helling, Y. Kawashima, V. Graham, D. Samra, K. L. Chubb, M. Min, L. B. F. M. Waters, V. Parmentier

Abstract: Having a short orbital period and being tidally locked makes WASP-43b an ideal candidate for JWST observations. Phase curve observations of an entire orbit will enable the mapping of the atmospheric structure across the planet, with different wavelengths of observation allowing different atmospheric depths to be seen. We provide insight into the details of the clouds that may form on WASP-43b in o… ▽ More Having a short orbital period and being tidally locked makes WASP-43b an ideal candidate for JWST observations. Phase curve observations of an entire orbit will enable the mapping of the atmospheric structure across the planet, with different wavelengths of observation allowing different atmospheric depths to be seen. We provide insight into the details of the clouds that may form on WASP-43b in order to prepare the forthcoming interpretation of the JWST and follow-up data. We utilize 3D GCM results as input for a kinetic, non-equilibrium model for mineral cloud particles, and for a kinetic model to study a photochemicaly-driven hydrocarbon haze component. Mineral condensation seeds form throughout the atmosphere of WASP-43b. This is in stark contrast to the ultra-hot Jupiters, like WASP-18b and HAT-P-7b. The dayside is loaded with few but large mineral cloud particles in addition to hydrocarbon haze particles of comparable abundance. Photochemically driven hydrocarbon haze appears on the dayside, but does not contribute to the cloud formation on the nightside. The geometrical cloud extension differs across the globe due to the changing thermodynamic conditions. Day and night differ by 6000km in pressure scale height. As reported for other planets, the C/O is not constant throughout the atmosphere. The mean molecular weight is approximately constant in a H2-dominated WASP-43b. WASP-43b is expected to be fully covered in clouds which are not homogeneously distributed throughout the atmosphere. The dayside and the terminator clouds will be a combination of mineral particles of locally varying size and composition, and of hydrocarbon hazes. The optical depth of hydrocarbon hazes is considerably lower than that of mineral cloud particles such that a wavelength-dependent radius measurement of WASP-43b would be determined by the mineral cloud particles but not by hazes. △ Less

Submitted 2 June, 2020; v1 submitted 28 May, 2020; originally announced May 2020.

Comments: 22 pages, accepted for publication A&A (correcting auto-correct)

Aluminium oxide in the atmosphere of hot Jupiter WASP-43b

Authors: Katy L. Chubb, Michiel Min, Yui Kawashima, Christiane Helling, Ingo Waldmann

Abstract: We have conducted a re-analysis of publicly available Hubble Space Telescope Wide Field Camera 3 (HST WFC3) transmission data for the hot-Jupiter exoplanet WASP-43b, using the Bayesian retrieval package Tau-REx. We report evidence of AlO in transmission to a high level of statistical significance (> 5-sigma in comparison to a flat model, and 3.4-sigma in comparison to a model with H2O only). We fi… ▽ More We have conducted a re-analysis of publicly available Hubble Space Telescope Wide Field Camera 3 (HST WFC3) transmission data for the hot-Jupiter exoplanet WASP-43b, using the Bayesian retrieval package Tau-REx. We report evidence of AlO in transmission to a high level of statistical significance (> 5-sigma in comparison to a flat model, and 3.4-sigma in comparison to a model with H2O only). We find no evidence of the presence of CO, CO2, or CH4 based on the available HST WFC3 data or on Spitzer IRAC data. We demonstrate that AlO is the molecule that fits the data to the highest level of confidence out of all molecules for which high-temperature opacity data currently exists in the infrared region covered by the HST WFC3 instrument, and that the subsequent inclusion of Spitzer IRAC data points in our retrieval further supports the presence of AlO. H2O is the only other molecule we find to be statistically significant in this region. AlO is not expected from the equilibrium chemistry at the temperatures and pressures of the atmospheric layer that is being probed by the observed data. Its presence therefore implies direct evidence of some disequilibrium processes with links to atmospheric dynamics. Implications for future study using instruments such as the James Webb Space Telescope (JWST) are discussed, along with future opacity needs. Comparisons are made with previous studies into WASP-43b. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: Accepted for publication in A&A. 16 pages, 7 figures

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

Mineral snowflakes on exoplanets and brown dwarfs: Effects of micro-porosity, size distributions, and particle shape

Authors: Dominic Samra, Christiane Helling, Michiel Min

Abstract: Exoplanet atmosphere characterisation has become an important tool in understanding exoplanet formation, evolution. However, clouds remain a key challenge for characterisation: upcoming space telescopes (e.g. JWST, ARIEL) and ground-based high-resolution spectrographs (e.g. CRIRES+) will produce data requiring detailed understanding of cloud formation and cloud effects. We aim to understand how th… ▽ More Exoplanet atmosphere characterisation has become an important tool in understanding exoplanet formation, evolution. However, clouds remain a key challenge for characterisation: upcoming space telescopes (e.g. JWST, ARIEL) and ground-based high-resolution spectrographs (e.g. CRIRES+) will produce data requiring detailed understanding of cloud formation and cloud effects. We aim to understand how the micro-porosity of cloud particles affects the cloud structure, particle size, and material composition. We examine the spectroscopic effects of micro-porous particles, the particle size distribution, and non-spherical cloud particles. We expanded our kinetic non-equilibrium cloud formation model and use a grid of prescribed 1D (T_gas-p_gas) DRIFT-PHOENIX profiles. We applied the effective medium theory and the Mie theory to model the spectroscopic properties of cloud particles with micro-porosity and a derived particle size distribution. We used a statistical distribution of hollow spheres to represent the effects of non-spherical cloud particles. Highly micro-porous cloud particles (90% vacuum) have a larger surface area, enabling efficient bulk growth higher in the atmosphere than for compact particles. Increases in single-scattering albedo and cross-sectional area for these mineral snowflakes cause the cloud deck to become optically thin only at a wavelength of ~100 ${\rm μm}$ instead of at the ~20 ${\rm μm}$ for compact cloud particles. A significant enhancement in albedo is also seen when cloud particles occur with a locally changing Gaussian size distribution. Non-spherical particles increase the opacity of silicate spectral features, which further increases the wavelength at which the clouds become optically thin. JWST MIRI will be sensitive to signatures of micro-porous and non-spherical cloud particles based on the wavelength at which clouds are optically thin. △ Less

Submitted 28 April, 2020; originally announced April 2020.

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

Journal ref: A&A 639, A107 (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

ALMA and NACO observations towards the young exoring transit system J1407 (V1400 Cen)

Authors: M. A. Kenworthy, P. D. Klaassen, M. Min, N. van der Marel, A. J. Bohn, M. Kama, A. Triaud, A. Hales, J. Monkiewicz, E. Scott, E. E. Mamajek

Abstract: Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimete… ▽ More Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the field centred on J1407 in the 340 GHz (Band 7) continuum in order to determine the flux and astrometric location of the ring system relative to the star. We used the VLT/NACO camera to observe the J1407 system in March 2019 and to search for the central planetary mass object at thermal infrared wavelengths. We detect no point source at the expected location of J1407, and derive an upper limit $3σ$ level of $57.6~μ\rm{Jy}$. There is a point source detected at an angular separation consistent with the expected location for a free-floating ring system that occulted J1407 in May 2007, with a flux of $89~μ\rm{Jy}$ consistent with optically thin dust surrounding a massive substellar companion. At 3.8 microns with the NACO camera, we detect the star J1407 but no other additional point sources within 1.3 arcseconds of the star, with a lower bound on the sensitivity of $6M_{Jup}$ at the location of the ALMA source, and down to $4M_{Jup}$ in the sky background limit. The ALMA upper limit at the location of J1407 implies that a hypothesised bound ring system is composed of dust smaller than $1\rm{~mm}$ in size, implying a young ring structure. The detected ALMA source has multiple interpretations, including: (i) it is an unbound substellar object surrounded by warm dust in Sco-Cen with an upper mass limit of $6M_{Jup}$, or (ii) it is a background galaxy. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: Accepted for publication in A&A (6 pages, 5 figures). Reduced data and reduction scripts on GitHub at https://github.com/mkenworthy/j1407_alma_detection

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

Dust evolution in pre-stellar cores

Authors: Charlène Lefèvre, Laurent Pagani, Bilal Ladjelate, Michiel Min, Hiroyuki Hirashita, Robert Zylka

Abstract: Dust grains are the building {blocks} of future planets. They evolve in size, shape and composition during the life cycle of the interstellar medium. We seek to understand the process which leads from diffuse medium grains to dust grains in the vicinity of protostars inside disks. As a first step, we propose to characterize the dust evolution inside pre-stellar cores thanks to multi-wavelength obs… ▽ More Dust grains are the building {blocks} of future planets. They evolve in size, shape and composition during the life cycle of the interstellar medium. We seek to understand the process which leads from diffuse medium grains to dust grains in the vicinity of protostars inside disks. As a first step, we propose to characterize the dust evolution inside pre-stellar cores thanks to multi-wavelength observations. We will present how NIKA2 maps are crucial to better constrain dust properties and {we will} introduce SIGMA: a new flexible dust model in open access. △ Less

Submitted 14 November, 2019; v1 submitted 8 November, 2019; originally announced November 2019.

Comments: 6 pages, 5 figures, to appear in the proceedings of the international conference entitled mm Universe @ NIKA2, Grenoble (France), June 2019, EPJ Web of conferences

VLTI/PIONIER survey of disks around post-AGB binaries. Dust sublimation physics rules!

Authors: J. Kluska, H. Van Winckel, M. Hillen, J. -P. Berger, D. Kamath, J. -B. Le Bouquin, M. Min

Abstract: Post-AGB binaries are surrounded by circumbinary disks of gas and dust that are similar to protoplanetary disks found around young stars. We aim to understand the structure of these disks and identify the physical phenomena at play in their very inner regions. We want to understand the disk-binary interaction and to further investigate the comparison with protoplanetary disks. We have conducted an… ▽ More Post-AGB binaries are surrounded by circumbinary disks of gas and dust that are similar to protoplanetary disks found around young stars. We aim to understand the structure of these disks and identify the physical phenomena at play in their very inner regions. We want to understand the disk-binary interaction and to further investigate the comparison with protoplanetary disks. We have conducted an interferometric snapshot survey of 23 post-AGB binaries in the near-infrared (H-band) using VLTI/PIONIER. We have fitted the multiwavelength visibilities and closure phases with purely geometrical models with an increasing complexity in order to retrieve the sizes, temperatures and flux ratios of the different components All sources are resolved and the different components contributing to the H-band flux are dissected. The environment of these targets is very complex: 13/23 targets need models with thirteen or more parameters to fit the data. We find that the inner disk rims follow and extend the size-luminosity relation established for disks around young stars with an offset toward larger sizes. The measured temperature of the near-infrared circumstellar emission of post-AGB binaries is lower (Tsub~1200K) than for young stars, probably due to a different dust mineralogy and/or gas density in the dust sublimation region. The dusty inner rims of the circumbinary disks around post-AGB binaries are ruled by dust sublimation physics. Additionally, a significant amount of the circumstellar H-band flux is over-resolved (14 targets have more than 10% of their non-stellar flux over-resolved) hinting for more structure from a yet unknown origin (disk structure or outflow). The amount of over-resolved flux is larger than around young stars. Due to the complexity of these targets, interferometric imaging is a necessary tool to reveal the interacting inner regions in a model-independent way. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: Accepted for publication in A&A, 40 pages

Journal ref: A&A 631, A108 (2019)

ESA Voyage 2050 White Paper: Detecting life outside our solar system with a large high-contrast-imaging mission

Authors: Ignas Snellen, Simon Albrecht, Guillem Anglada-Escude, Isabelle Baraffe, Pierre Baudoz, Willy Benz, Jean-Luc Beuzit, Beth Biller, Jayne Birkby, Anthony Boccaletti, Roy van Boekel, Jos de Boer, Matteo Brogi, Lars Buchhave, Ludmila Carone, Mark Claire, Riccardo Claudi, Brice-Olivier Demory, Jean-Michel Desert, Silvano Desidera, Scott Gaudi, Raffaele Gratton, Michael Gillon, John Lee Grenfell, Olivier Guyon , et al. (42 additional authors not shown)

Abstract: In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestr… ▽ More In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: White paper for ESA Voyage 2050; 24 pages

Optical polarized phase function of the HR\,4796A dust ring

Authors: J. Milli, N. Engler, H. M. Schmid, J. Olofsson, F. Menard, Q. Kral, A. Boccaletti, P. Thebault, E. Choquet, D. Mouillet, A. -M. Lagrange, J. C. Augereau, C. Pinte, G. Chauvin, C. Dominik, C. Perrot, A. Zurlo, T. Henning, M. Min, J. L. Beuzit, H. Avenhaus, A. Bazzon, T. Moulin, M. Llored, O. Moeller-Nilsson , et al. (2 additional authors not shown)

Abstract: The scattering properties of the dust originating from debris discs are still poorly known. The analysis of scattered light is however a powerful remote-sensing tool to understand the physical properties of dust particles orbiting other stars. Scattered light is indeed widely used to characterise the properties of cometary dust in the solar system. We aim to measure the morphology and scattering… ▽ More The scattering properties of the dust originating from debris discs are still poorly known. The analysis of scattered light is however a powerful remote-sensing tool to understand the physical properties of dust particles orbiting other stars. Scattered light is indeed widely used to characterise the properties of cometary dust in the solar system. We aim to measure the morphology and scattering properties of the dust from the debris ring around HR4796A in polarised optical light. We obtained high-contrast polarimetric images of HR4796A in the wavelength range 600-900nm with the SPHERE / ZIMPOL instrument on the Very Large Telescope. We measured for the first time the polarised phase function of the dust in a debris system over a wide range of scattering angles in the optical. We confirm that it is incompatible with dust particles being compact spheres under the assumption of the Mie theory, and propose alternative scenarios compatible with the observations, such as particles with irregular surface roughness or aggregate particles. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: 12 pages, 11 figures

Journal ref: A&A 626, A54 (2019)

Consistent dust and gas models for protoplanetary disks IV. A panchromatic view of protoplanetary disks

Authors: O. Dionatos, P. Woitke, M. Guedel, P. Degroote, A. Liebhart, F. Anthonioz, S. Antonellini, C. Baldovin-Saavedra, A. Carmona, C. Dominik, J. Greaves, J. D. Ilee, I. Kamp, F. Menard, M. Min, C. Pinte, C. Rab, L. Rigon, W. F. Thi, L. B. F. M. Waters

Abstract: Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present… ▽ More Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disc mass and shape parameters. dust size and PAH properties are also provided for completeness. Targets were selected based on their properties data availability. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Compiled datasets along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow to backtrack them to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing. △ Less

Submitted 28 February, 2019; originally announced February 2019.

Comments: 32 pages, 11 figures, A&A accepted

Journal ref: A&A 625, A66 (2019)

The ARCiS framework for Exoplanet Atmospheres: The Cloud Transport Model

Authors: Chris W. Ormel, Michiel Min

Abstract: Understanding of clouds is instrumental in interpreting current and future spectroscopic observations of exoplanets. Modelling clouds consistently is complex, since it involves many facets of chemistry, nucleation theory, condensation physics, coagulation, and particle transport. We develop a simple physical model for cloud formation and transport, efficient and versatile enough that it can be use… ▽ More Understanding of clouds is instrumental in interpreting current and future spectroscopic observations of exoplanets. Modelling clouds consistently is complex, since it involves many facets of chemistry, nucleation theory, condensation physics, coagulation, and particle transport. We develop a simple physical model for cloud formation and transport, efficient and versatile enough that it can be used in modular fashion for parameter optimization searches of exoplanet atmosphere spectra. The transport equations are formulated in 1D, accounting for sedimentation and diffusion. The grain size is obtained through a moment method. For simplicity, only one cloud species is considered and the nucleation rate is parametrized. From the resulting physical profiles we simulate transmission spectra covering the visual to mid-IR wavelength range. We apply our models towards KCl clouds in the atmosphere of GJ1214 b and towards MgSiO3 clouds of a canonical hot-Jupiter. We find that larger cloud diffusivity $K_{zz}$ increases the thickness of the cloud, pushing the $τ=1$ surface to a lower pressure layer higher in the atmosphere. A larger nucleation rate also increases the cloud thickness while it suppresses the grain size. Coagulation is most important at high nuclei injection rates ($\dotΣ_n$) and low $K_{zz}$. We find that the investigated combinations of $K_{zz}$ and $\dotΣ_n$ greatly affect the transmission spectra in terms of the slope at near-IR wavelength (a proxy for grain size), the molecular features seen at ~1\micr (which disappear for thick clouds, high in the atmosphere), and the 10\micr silicate feature, which becomes prominent for small grains high in the atmosphere. The result of our hybrid approach -- aimed to provide a good balance between physical consistency and computational efficiency -- is ideal towards interpreting (future) spectroscopic observations of exoplanets. △ Less

Submitted 23 January, 2019; v1 submitted 12 December, 2018; originally announced December 2018.

Comments: language and other tiny corrections

Journal ref: A&A 622, A121 (2019)

Consistent dust and gas models for protoplanetary disks III. Models for selected objects from the FP7 DIANA project

Authors: P. Woitke, I. Kamp, S. Antonellini, F. Anthonioz, C. Baldovin-Saveedra, A. Carmona, O. Dionatos, C. Dominik, J. Greaves, M. Güdel, J. D. Ilee, A. Liebhardt, F. Menard, M. Min, C. Pinte, C. Rab, L. Rigon, W. -F. Thi, N. Thureau, L. B. F. M. Waters

Abstract: The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line ma… ▽ More The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana). △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 70 pages, 35 figures, 22 tables, accepted by PASP

Chromatic transit light curves of disintegrating rocky planets

Authors: A. R. Ridden-Harper, C. U. Keller, M. Min, R. van Lieshout, I. A. G. Snellen

Abstract: Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition a… ▽ More Context. Kepler observations have revealed a class of short period exoplanets, of which Kepler-1520 b is the prototype, which have comet-like dust tails thought to be the result of small, rocky planets losing mass. The shape and chromaticity of the transits constrain the properties of the dust particles originating from the planet's surface, offering a unique opportunity to probe the composition and geophysics of rocky exoplanets. Aims. We aim to approximate the average Kepler long-cadence light curve of Kepler-1520 b and investigate how the optical thickness and transit cross-section of a general dust tail can affect the observed wavelength dependence and depth of transit light curves. Methods. We developed a new 3D model that ejects sublimating particles from the planet surface to build up a dust tail, assuming it to be optically thin, and used 3D radiative transfer computations that fully treat scattering using the distribution of hollow spheres (DHS) method, to generate transit light curves between 0.45 and 2.5 $μ$m. Results. We show that the transit depth is wavelength independent for optically thick tails, potentially explaining why only some observations indicate a wavelength dependence. From the 3D nature of our simulated tails, we show that their transit cross-sections are related to the component of particle ejection velocity perpendicular to the planet's orbital plane and use this to derive a minimum ejection velocity of 1.2 kms$^{-1}$. To fit the average transit depth of Kepler-1520 b of 0.87%, we require a high dust mas-loss rate of 7 $-$ 80 M$_\oplus$ Gyr$^{-1}$ which implies planet lifetimes that may be inconsistent with the observed sample. Therefore, these mass-loss rates should be considered to be upper limits. △ Less

Submitted 20 July, 2018; originally announced July 2018.

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

Journal ref: A&A 618, A97 (2018)

Modelling mid-infrared molecular emission lines from T Tauri stars

Authors: P. Woitke, M. Min, W. -F. Thi, C. Roberts, A. Carmona, I. Kamp, F. Menard, C. Pinte

Abstract: We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission… ▽ More We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission lines of SI, SII, SIII, SiII, FeII, NeII, NeIII, ArII and ArIII. In contrast to previously published works, we do not treat the abundances of the molecules nor the temperature in the disc as free parameters, but use the complex results of detailed 2D ProDiMo disc models concerning gas and dust temperature structure, and molecular concentrations. FLiTs computes the line emission spectra by ray tracing in an efficient, fast and reliable way. The results are broadly consistent with R=600 Spitzer/IRS observational data of T Tauri stars concerning line strengths, colour, and line ratios. In order to achieve that agreement, however, we need to assume either a high gas/dust mass ratio of order 1000, or the presence of illuminated disc walls at distances of a few au. The molecules in these walls cannot be photo-dissociated easily by UV because of the large densities in the walls which favour their re-formation. Most observable molecular emission lines are found to be optically thick, rendering a standard analysis with column densities difficult. We find that the difference between gas and dust temperatures in the disc surface is important for the line formation. We briefly discuss the effects of C/O ratio and choice of chemical rate network on these results. Our analysis offers new ways to infer the chemical and temperature structure of T Tauri discs from future JWST/MIRI observations, and to possibly detect secondary illuminated disc walls based on their specific mid-IR molecular signature. △ Less

Submitted 25 July, 2018; v1 submitted 16 July, 2018; originally announced July 2018.

Comments: accepted by A&A, 13 figures, 4 tables

Journal ref: A&A 618, A57 (2018)

The perturbed sublimation rim of the dust disk around the post-AGB binary IRAS08544-4431

Authors: J. Kluska, M. Hillen, H. Van Winckel, R. Manick, M. Min, S. Regibo, P. Royer

Abstract: Context: Post-Asymptotic Giant Branch (AGB) binaries are surrounded by stable dusty and gaseous disks similar to the ones around young stellar objects. Whereas significant effort is spent on modeling observations of disks around young stellar objects, the disks around post-AGB binaries receive significantly less attention, even though they pose significant constraints on theories of disk physics a… ▽ More Context: Post-Asymptotic Giant Branch (AGB) binaries are surrounded by stable dusty and gaseous disks similar to the ones around young stellar objects. Whereas significant effort is spent on modeling observations of disks around young stellar objects, the disks around post-AGB binaries receive significantly less attention, even though they pose significant constraints on theories of disk physics and binary evolution. Aims: We want to examine the structure of and phenomena at play in circumbinary disks around post-AGB stars. We continue the analysis of our near-infrared interferometric image of the inner rim of the circumbinary disk around IRAS08544-4431. We want to understand the physics governing this inner disk rim. Methods: We use a radiative transfer model of a dusty disk to reproduce simultaneously the photometry as well as the near-infrared interferometric dataset on IRAS08544-4431. The model assumes hydrostatic equilibrium and takes dust settling self-consistently into account. Results: The best-fit radiative transfer model shows excellent agreement with the spectral energy distribution up to mm wavelengths as well as with the PIONIER visibility data. It requires a rounded inner rim structure, starting at a radius of 8.25 au. However, the model does not fully reproduce the detected over-resolved flux nor the azimuthal flux distribution of the inner rim. While the asymmetric inner disk rim structure is likely to be the consequence of disk-binary interactions, the origin of the additional over-resolved flux remains unclear. Conclusions: As in young stellar objects, the disk inner rim of IRAS08544-4431 is ruled by dust sublimation physics. Additional observations are needed to understand the origin of the extended flux and the azimuthal perturbation at the inner rim of the disk. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: Accepted for publication in A&A, 13 figures, 13 pages

Journal ref: A&A 616, A153 (2018)

The contribution of the ARIEL space mission to the study of planetary formation

Authors: D. Turrini, Y. Miguel, T. Zingales, A. Piccialli, R. Helled, A. Vazan, F. Oliva, G. Sindoni, O. Panić, J. Leconte, M. Min, S. Pirani, F. Selsis, V. Coudé du Foresto, A. Mura, P. Wolkenberg

Abstract: The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanet… ▽ More The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanetary population, however, do not provide enough information to break the intrinsic degeneracy of their histories, as different evolutionary tracks can result in the same final configurations. The lessons learned from the Solar System indicate us that the solution to this problem lies in the information contained in the composition of planets. The goal of the Atmospheric Remote-Sensing Infrared Exoplanet Large-survey (ARIEL), one of the three candidates as ESA M4 space mission, is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres, which should show minimal condensation and sequestration of high-Z materials and thus reveal their bulk composition across all main cosmochemical elements. In this work we will review the most outstanding open questions concerning the way planets form and the mechanisms that contribute to create habitable environments that the compositional information gathered by ARIEL will allow to tackle △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: 20 pages, 6 figures, accepted for publication on Experimental Astronomy, to appear in the special issue on the ESA space mission ARIEL

Diagnostic value of far-IR water ice features in T Tauri disks

Authors: Inga Kamp, Anton Scheepstra, Michiel Min, Lucia Klarmann, Pablo Riviere-Marichalar

Abstract: This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ic… ▽ More This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ice features at 45 and 63 micron change with key disk properties: disk size, grain sizes, disk dust mass, dust settling, and ice thickness. In addition, a series of models is devised to calculate the water ice emission features from warmup, direct deposit and cooldown scenarios of the water ice in disks. Photodesorption from icy grains in disk surfaces weakens the mid-IR water ice features by factors 4-5. The far-IR water ice emission features originate from small grains at the surface snow line in disks at distance of 10-100 au. Unless this reservoir is missing in disks (e.g. transitional disks with large cavities), the feature strength is not changing. Grains larger than 10 micron do not contribute to the features. Grain settling (using turbulent description) is affecting the strength of the ice features by at most 15%. The strength of the ice feature scales with the disk dust mass and water ice fraction on the grains, but saturates for dust masses larger than 1.e-4 Msun and for ice mantles that increase the dust mass by more than 50%. The various thermal histories of water ice leave an imprint on the shape of the features (crystalline/amorphous) as well as on the peak strength and position of the 45 micron feature. SOFIA/HIRMES can only detect crystalline ice features much stronger than simulated in our standard T Tauri disk model in deep exposures (1 hr). SPICA/SAFARI can detect the typical ice features in our standard T Tauri disk model in short exposures (10 min). (abbreviated) △ Less

Submitted 15 April, 2018; originally announced April 2018.

Comments: accepted for publication in A&A

Dust modeling of the combined ALMA and SPHERE datasets of HD163296. Is HD163296 really a Meeus group II disk?

Authors: G. A. Muro-Arena, C. Dominik, L. B. F. M. Waters, M. Min, L. Klarmann, C. Ginski, A. Isella, M. Benisty, A. Pohl, A. Garufi, J. Hagelberg, M. Langlois, F. Menard, C. Pinte, E. Sezestre, G. van der Plas, M. Villenave, A. Delboulbé, Y. Magnard, O. Möller-Nilsson, J. Pragt, P. Rabou, R. Roelfsema

Abstract: Context. Multi-wavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aimed to construct a 3-dimensional model of HD 163296 capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We… ▽ More Context. Multi-wavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aimed to construct a 3-dimensional model of HD 163296 capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We want to determine why the SED of HD 163296 is intermediary between the otherwise well-separated group I and group II Herbig stars. Methods. The disk was modelled using the Monte Carlo radiative transfer code MCMax3D. The radial dust surface density profile was modelled after the ALMA observations, while the polarized scattered light observations were used to constrain the inclination of the inner disk component and turbulence and grain growth in the outer disk. Results. While three rings are observed in the disk midplane in millimeter thermal emission at $\sim$80, 124 and 200 AU, only the innermost of these is observed in polarized scattered light, indicating a lack of small dust grains on the surface of the outer disk. We provide two models capable of explaining this difference. The first model uses increased settling in the outer disk as a mechanism to bring the small dust grains on the surface of the disk closer to the midplane, and into the shadow cast by the first ring. The second model uses depletion of the smallest dust grains in the outer disk as a mechanism for decreasing the optical depth at optical and NIR wavelengths. In the region outside the fragmentation-dominated regime, such depletion is expected from state-of-the-art dust evolution models. We studied the effect of creating an artificial inner cavity in our models, and conclude that HD 163296 might be a precursor to typical group I sources. △ Less

Submitted 9 February, 2018; originally announced February 2018.

Comments: 15 pages, 12 figures

Journal ref: A&A 614, A24 (2018)

Laboratory mid-IR spectra of equilibrated and igneous meteorites. Searching for observables of planetesimal debris

Authors: B. L. de Vries, H. Skogby, L. B. F. M. Waters, M. Min

Abstract: Meteorites contain minerals from Solar System asteroids with different properties (like size, presence of water, core formation). We provide new mid-IR transmission spectra of powdered meteorites to obtain templates of how mid-IR spectra of asteroidal debris would look like. This is essential for interpreting mid-IR spectra of past and future space observatories, like the James Webb Space Telescop… ▽ More Meteorites contain minerals from Solar System asteroids with different properties (like size, presence of water, core formation). We provide new mid-IR transmission spectra of powdered meteorites to obtain templates of how mid-IR spectra of asteroidal debris would look like. This is essential for interpreting mid-IR spectra of past and future space observatories, like the James Webb Space Telescope. We show that the transmission spectra of wet and dry chondrites, carbonaceous and ordinary chondrites and achondrite and chondrite meteorites are distinctly different in a way one can distinguish in astronomical mid-IR spectra. The two observables that spectroscopically separate the different meteorites groups (and thus the different types of parent bodies) are the pyroxene-olivine feature strength ratio and the peak shift of the olivine spectral features due to an increase in the iron concentration of the olivine. △ Less

Submitted 30 November, 2017; originally announced November 2017.

Journal ref: Published in Icarus: https://doi.org/10.1016/j.icarus.2017.10.029

X-ray radiative transfer in protoplanetary disks - The role of dust and X-ray background fields

Authors: Ch. Rab, M. Güdel, P. Woitke, I. Kamp, W. -F. Thi, M. Min, G. Aresu, R. Meijerink

Abstract: The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the c… ▽ More The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 16 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 609, A91 (2018)