Source author record

Michiel Min

Michiel Min appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2023arXiv

Hydrocarbon chemistry in inner regions of planet forming disks

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 ....

preprint2022arXiv

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

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.

preprint2022arXiv

A retrieval challenge exercise for the Ariel mission

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.

preprint2022arXiv

Convolutional neural networks as an alternative to Bayesian retrievals

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Investigating the detectability of hydrocarbons in exoplanet atmospheres with JWST

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

Dust evolution in pre-stellar cores

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.

preprint2016arXiv

Determining the mid-plane conditions of circumstellar discs using gas and dust modelling: a study of HD 163296

The mass of gas in protoplanetary discs is a quantity of great interest for assessing their planet formation potential. Disc gas masses are, however, traditionally inferred from measured dust masses by applying an assumed standard gas-to-dust ratio of $g/d=100$. Furthermore, measuring gas masses based on CO observations has been hindered by the effects of CO freeze-out. Here we present a novel approach to study the mid-plane gas by combining C$^{18}$O line modelling, CO snowline observations and the spectral energy distribution (SED) and selectively study the inner tens of au where freeze-out is not relevant. We apply the modelling technique to the disc around the Herbig Ae star HD 163296 with particular focus on the regions within the CO snowline radius, measured to be at 90 au in this disc. Our models yield the mass of C$^{18}$O in this inner disc region of $M_{\text{C}^{18}\text{O}}(<90\,\text{au})\sim 2\times10^{-8}$ M$_\odot$. We find that most of our models yield a notably low $g/d<20$, especially in the disc mid-plane ($g/d<1$). Our only models with a more interstellar medium (ISM)-like $g/d$ require C$^{18}$O to be underabundant with respect to the ISM abundances and a significant depletion of sub-micron grains, which is not supported by scattered light observations. Our technique can be applied to a range of discs and opens up a possibility of measuring gas and dust masses in discs within the CO snowline location without making assumptions about the gas-to-dust ratio.

preprint2015arXiv

On the importance of scattering at 8 microns: Brighter than you think

Context. Extinction and emission of dust models need for observational constraints to be validated. The coreshine phenomenon has already shown the importance of scattering in the 3 to 5 micron range and its ability to validate dust properties for dense cores. Aims. We want to investigate whether scattering can also play a role at longer wavelengths and to place even tighter constraints on the dust properties. Methods. We analyze the inversion of the Spitzer 8 micron map of the dense molecular cloud L183, to examine the importance of scattering as a potential contributor to the line-of-sight extinction. Results. The column density deduced from the inversion of the 8 micron map, when we neglect scattering, disagrees with all the other column density measurements of the same region. Modeling confirms that scattering at 8 microns is not negligible with an intensity of several hundred kJy per sr. This demonstrates the need of efficiently scattering dust grains at MIR wavelengths up to 8 microns. Coagulated aggregates are good candidates and might also explain the discrepancy at high extinction between E(J-K) et tau(9.7) toward dense molecular clouds. Further investigation requires considering efficiently scattering dust grains including ices as realistic dust models.

preprint2015arXiv

The snow line in viscous disks around low-mass stars: implications for water delivery to terrestrial planets in the habitable zone

The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.

preprint2014arXiv

SEEDS Adaptive Optics Imaging of the Asymmetric Transition Disk Oph IRS 48 in Scattered Light

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$μ$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

preprint2012arXiv

SPICES: Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System.

preprint2012arXiv

Why circumstellar disks are so faint in scattered light: The case of HD 100546

Scattered light images of circumstellar disks play an important role in characterizing the planet forming environments around young stars. The characteristic size of the scattering dust grains can be estimated from the observed brightness asymmetry between the front and back side of the disk, for example using standard Mie theory. However such models often overpredict their brightness by one or two orders of magnitude, and have difficulty explaining very red disk colors. We aim to develop a dust model that explains simultaneously the observed disk surface brightness, colors and asymmetry in scattered light, focusing on constraining grain sizes. We use the 2D radiative transfer code MCMax with anisotropic scattering to explore the effects of grain size on synthetic scattered light images of circumstellar disks. We compare the results with surface brightness profiles of the protoplanetary disk HD 100546 in scattered light at wavelengths from 0.4 to 2.2 micron. We find that extreme forward scattering by micron sized particles lowers the /effective/ dust albedo and creates a faint and red disk that /appears/ only slightly forward scattering. For the outer (>~100 AU) disk of HD 100546 we derive a minimum grain size of 2.5 micron, likely present in the form of aggregates. Intermediate sized grains are too bright, whereas smaller grains are faint and scatter more isotropically, but also produce disk colors that are too blue. Observed surface brightness asymmetries alone are not sufficient to constrain the grain size in circumstellar disks. Additional information, such as the brightness and colors of the disk are needed to provide additional constraints.

preprint2011arXiv

The effects of disk and dust structure on observed polarimetric images of protoplanetary disks

Imaging polarimetry is a powerful tool for imaging faint circumstellar material. For a correct analysis of observations we need to fully understand the effects of dust particle parameters, as well as the effects of the telescope, atmospheric seeing, and assumptions about the data reduction and processing of the observed signal. Here we study the major effects of dust particle structure, size-dependent grain settling, and instrumental properties. We performed radiative transfer modeling using different dust particle models and disk structures. To study the influence of seeing and telescope diffraction we ran the models through an instrument simulator for the ExPo dual-beam imaging polarimeter mounted at the 4.2m William Herschel Telescope (WHT). Particle shape and size have a strong influence on the brightness and detectability of the disks. In the simulated observations, the central resolution element also contains contributions from the inner regions of the protoplanetary disk besides the unpolarized central star. This causes the central resolution element to be polarized, making simple corrections for instrumental polarization difficult. This effect strongly depends on the spatial resolution, so adaptive optics systems are needed for proper polarization calibration. We find that the commonly employed homogeneous sphere model gives results that differ significantly from more realistic models. For a proper analysis of the wealth of data available now or in the near future, one must properly take the effects of particle types and disk structure into account. The observed signal depends strongly on the properties of these more realistic models, thus providing a potentially powerful diagnostic. We conclude that it is important to correctly understand telescope depolarization and calibration effects for a correct interpretation of the degree of polarization.

preprint2010arXiv

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.