Search

The Mid-infrared ELT Imager and Spectrograph (METIS) is one of the first-generation scientific instruments for the ELT, built under the supervision of ESO by a consortium of research institutes across and beyond Europe. Designed to cover the 3 to 13 $\mu$m wavelength range, METIS had its final design reviewed in Fall 2022, and has then entered in earnest its manufacture, assembly, integration, and test (MAIT) phase. Here, we present the final design of the METIS high-contrast imaging (HCI) modes. We detail the implementation of the two main coronagraphic solutions selected for METIS, namely the vortex coronagraph and the apodizing phase plate, including their combination with the high-resolution integral field spectrograph of METIS, and briefly describe their respective backup plans (Lyot coronagraph and shaped pupil plate). We then describe the status of the MAIT phase for HCI modes, including a review of the final design of individual components such as the vortex phase masks, the grayscale ring apodizer, and the apodizing phase plates, as well as a description of their on-going performance tests and of our plans for system-level integration and tests. Using end-to-end simulations, we predict the performance that will be reached on sky by the METIS HCI modes in presence of various environmental and instrumental disturbances, including non-common path aberrations and water vapor seeing, and discuss our strategy to mitigate these various effects. We finally illustrate with mock observations and data processing that METIS should be capable of directly imaging temperate rocky planets around the nearest stars., Comment: 16 pages, paper presented at SPIE Astronomical Telescopes + Instrumentation 2024

New longwave HgCdTe detectors are critical to upcoming plans for ground-based infrared astronomy. These detectors, with fast-readouts and deep well-depths, will be key components of extremely large telescope instruments and therefore must be well understood prior to deployment. We analyze one such HgCdTe detector, a Teledyne Imaging Sensors GeoSnap, at the University of Michigan. We find that the properties of the GeoSnap are consistent with expectations from analysis of past devices. The GeoSnap has a well-depth of 2.75 million electrons per pixel, a read noise of 360 e-/pix, and a dark current of 330,000 e-/s/pix at 45 K. The device experiences 1/f noise which can be mitigated relative to half-well shot noise with modest frequency image differencing. The GeoSnap's quantum efficiency is calculated to be 79.7 +- 8.3 % at 10.6 microns. Although the GeoSnap's bad pixel fraction, on the order of 3%, is consistent with other GeoSnap devices, close to a third of the bad pixels in this detector are clustered in a series of 31 "leopard" spots spread across the detector plane. We report these properties and identify additional analyses that will be performed on future GeoSnap detectors., Comment: 16 pages, 9 figures, Submitted to SPIE Astronomical Telescopes + Instrumentation 2024

We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights into the diversity of the observed exoplanet populations. XUE targets 15 disks in three areas of NGC 6357, which hosts numerous massive OB stars, among which some of the most massive stars in our Galaxy. Thanks to JWST we can, for the first time, study the effect of external irradiation on the inner ($< 10$ au), terrestrial-planet forming regions of proto-planetary disks. In this study, we report on the detection of abundant water, CO, CO$_2$, HCN and C$_2$H$_2$ in the inner few au of XUE 1, a highly irradiated disk in NGC 6357. In addition, small, partially crystalline silicate dust is present at the disk surface. The derived column densities, the oxygen-dominated gas-phase chemistry, and the presence of silicate dust are surprisingly similar to those found in inner disks located in nearby, relatively isolated low-mass star-forming regions. Our findings imply that the inner regions of highly irradiated disks can retain similar physical and chemical conditions as disks in low-mass star-forming regions, thus broadening the range of environments with similar conditions for inner disk rocky planet formation to the most extreme star-forming regions in our Galaxy., Comment: Accepted for publication in ApJ Letters. 20 pages, 7 figures

Hot Jupiters, particularly those with temperature higher than 2000 K are the best sample of planets that allow in-depth characterization of their atmospheres. We present here a thermal emission study of the ultra hot Jupiter WASP-103 b observed in two secondary eclipses with CFHT/WIRCam in J and Ks bands. By means of high precision differential photometry, we determine eclipse depths in J and Ks to an accuracy of 220 and 270 ppm, which are combined with the published HST/WFC3 and Spitzer data to retrieve a joint constraints on the properties of WASP-103 b dayside atmosphere. We find that the atmosphere is best fit with a thermal inversion layer included. The equilibrium chemistry retrieval indicates an enhanced C/O (1.35$^{+0.14}_{-0.17}$) and a super metallicity with [Fe/H]$=2.19^{+0.51}_{-0.63}$ composition. Given the near-solar metallicity of WASP-103 of [Fe/H]=0.06, this planet seems to be $\sim$100 more abundant than its host star. The free chemistry retrieval analysis yields a large abundance of FeH, H$^{-}$, CO$_2$ and CH$_4$. Additional data of better accuracy from future observations of JWST should provide better constraint of the atmospheric properties of WASP-103b.

We use the long-slit spectroscopy mode of the VISIR-NEAR experiment to perform diffraction-limited observations of eight nearby Herbig Ae protoplanetary disks. We extract spectra for various locations along the slit with a spectral resolution of R = 300 and perform a compositional fit at each spatial location using spectral templates of silicates and the four PAH bands. This yields the intensity vs. location profiles of each species. Results. We could obtain spatially-resolved intensity profiles of the PAH emission features in the N-band for five objects (AB Aurigae, HD 97048, HD 100546, HD 163296, and HD 169142). We observe two kinds of PAH emission geometry in our sample: centrally-peaked (HD 97048) and ring-like (AB Aurigae, HD 100546, HD 163296, and potentially HD 169142). Comparing the spatial PAH emission profiles with near-infrared scattered light images, we find a strong correlation in the disk sub-structure but a difference in radial intensity decay rate. The PAH emission shows a less steep decline with distance from the star. Finally, we find a correlation between the presence of (sub-) micron-sized silicate grains leading to the depletion of PAH emission within the inner regions of the disks. In this work, we find the following: (1) PAH emission traces the extent of Herbig Ae disks to a considerable radial distance. (2) The correlation between silicate emission within the inner regions of disks and the depletion of PAH emission can result from dust-mixing and PAH coagulation mechanisms and competition over UV photons. (3) For all objects in our sample, PAHs undergo stochastic heating across the entire spatial extent of the disk and are not saturated. (4) The difference in radial intensity decay rates between the PAHs and scattered-light profiles may be attributed to shadowing and dust-settling effects, which affect the scattering grains more than the PAHs.

Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations and calibrations of the telescope, thirteen specific working groups were created to coordinate the effort between ESO, the instrument consortia, and the wider community. We describe here the goals of these working groups as well as their achievements so far., Comment: 8 pages, 4 figures

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

The high-speed variability of the local water vapor content in the Earth atmosphere is a significant contributor to ground-based wavefront quality throughout the infrared domain. Unlike dry air, water vapor is highly chromatic, especially in the mid-infrared. This means that adaptive optics correction in the visible or near-infrared domain does not necessarily ensure a high wavefront quality at longer wavelengths. Here, we use literature measurements of water vapor seeing, and more recent infrared interferometric data from the Very Large Telescope Interferometer (VLTI), to evaluate the wavefront quality that will be delivered to the METIS mid-infrared camera and spectrograph for the Extremely Large Telescope (ELT), operating from 3 to 13 {\mu}m, after single-conjugate adaptive optics correction in the near-infrared. We discuss how the additional wavefront error due to water vapor seeing is expected to dominate the wavefront quality budget at N band (8-13 {\mu}m), and therefore to drive the performance of mid-infrared high-contrast imaging modes at ELT scale. Then we present how the METIS team is planning to mitigate the effect of water vapor seeing using focal-plane wavefront sensing techniques, and show with end-to-end simulations by how much the high-contrast imaging performance can be improved., Comment: 13 pages, paper presented at SPIE Astronomical Telescopes + Instrumentation 2022

AGB stars are one of the main sources of dust production in the Galaxy. However, it is not clear what this process looks like and where the dust is condensing in the circumstellar environment. By characterizing the location of the dust and the molecules in the close environment of an AGB star, we aim to achieve a better understanding the history of the dust formation process. We observed the carbon star R Scl with the VLTI-MATISSE instrument in L- and N-bands. The high angular resolution of the VLTI observations, combined with a large uv-plane coverage allowed us to use image reconstruction methods. To constrain the dust and molecules' location, we used two different methods: MIRA image reconstruction and the 1D code RHAPSODY. We found evidence of C2H2 and HCN molecules between 1 and 3.4 Rstar which is much closer to the star than the location of the dust (between 3.8 and 17.0 Rstar). We also estimated a mass-loss rate of 1.2+-0.4x10-6 Msun per yr. In the meantime, we confirmed the previously published characteristics of a thin dust shell, composed of amorphous carbon (amC) and silicon carbide (SiC). However, no clear SiC feature has been detected in the MATISSE visibilities. This might be caused by molecular absorption that can affect the shape of the SiC band at 11.3 micron. The appearance of the molecular shells is in good agreement with predictions from dynamical atmosphere models. For the first time, we co-located dust and molecules in the environment of an AGB star. We confirm that the molecules are located closer to the star than the dust. The MIRA images unveil the presence of a clumpy environment in the fuzzy emission region beyond 4.0 Rstar. Furthermore, with the available dynamic range and angular resolution, we did not detect the presence of a binary companion. Additional observations combining MATISSE and SAM-VISIR instrument should enable this detection in future studies., Comment: 19 pages, published in A&A

We present a new optical transmission spectrum of the hot Jupiter HAT-P-1b based on two transits observed with the Double Spectrograph (DBSP) on the Palomar 200-inch (P200) telescope. The DBSP transmission spectrum, covering a wavelength range from 3250 to 10007 \r{A}, is consistent with that observed with the Hubble Space Telescope (HST), but the former has a finer spectral resolution. The DBSP spectrum alone reveals the presence of a pressure broadened line wing for Na, the line cores for both Na and K, and tentative evidence for H$_2$O. We obtain consistent results from the spectral retrieval analyses performed on the DBSP-only dataset and the DBSP, HST, and Spitzer combined dataset. Our retrievals suggest a mostly clear atmosphere for HAT-P-1b, with a cloud coverage of $22^{+5}_{-3}$% that is dominated by enhanced haze. We derive subsolar abundances for Na, K, and C, and subsolar-to-solar for O. Future observations with James Webb Space Telescope and ground-based high-resolution spectrographs should be able to not only confirm the presence of these species but also stringently constrain the formation and migration pathways for HAT-P-1b., Comment: 24 pages, 14 figures. Accepted for publication in AJ

Active galactic nuclei play a key role in the evolution of galaxies, but their inner workings and physical connection to the host are poorly understood due to a lack of angular resolution. Infrared interferometry makes it possible to resolve the circumnuclear dust in the nearby Seyfert 2 galaxy, Circinus. Previous observations have revealed complex structures and polar dust emission but interpretation was limited to simple models. MATISSE makes it possible to image these structures for the first time. We observed the Circinus Galaxy with VLTI/MATISSE, producing 150 correlated flux spectra and 100 closure phase spectra. We reconstructed images in the N-band at ~10 mas resolution. We fit blackbody functions with dust extinction to several aperture-extracted fluxes from the images to produce a temperature distribution of central dusty structures. We find significant substructure in the circumnuclear dust: central unresolved flux of ~0.5 Jy, a thin disk 1.9 pc in diameter oriented along ~45 deg,and a ~4x1.5 pc polar emission extending orthogonal to the disk. The polar emission exhibits patchiness, which we attribute to clumpy dust. Flux enhancements to the east and west of the disk are seen for the first time. We distinguish the temperature profiles of the disk and of the polar emission: the disk shows a steep temperature gradient indicative of denser material; the polar profile is flatter, indicating clumpiness and/or lower dust density. The unresolved flux is fitted with a high temperature, ~370 K. The polar dust remains warm (~200 K) out to 1.5 pc from the disk. The recovered morphology and temperature distribution resembles modeling of accretion disks with radiation-driven winds at large scales, but we placed new constraints on the subparsec dust. The subparsec features imaged here place new constraints on the physical modeling of circumnuclear dust in active galaxies., Comment: Accepted for publication in A&A

Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and near-infrared (SPHERE; 1.5-2.5$\mu$m) wavelengths which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep $K_s$ (2.16$\mu$m) and L' (3.7$\mu$m) band images of AB Aurigae obtained with LMIRCam on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate ($\rho = 0.681", PA = 7.6^{\circ}$) had been previously claimed. The nature of a second innermost planet candidate ($\rho = 0.16'', PA = 203.9^{\circ}$) can not be investigated by the new data. We are able to derive 5$\sigma$ detection limits in both magnitude and mass for the system, going from 14 \Mjup at 0.3'' (49 au) down to 3-4 \Mjup at 0.6'' (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (< 0.5'') resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.5-0.7") showing a clear south-east/north-west asymmetry. This structure, observed for the first time at L'-band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk., Comment: 13 pages, 5 figures, accepted on ApJ

In the widely accepted 'Unified Model' solution of the classification puzzle of Active Galactic Nuclei, the orientation of a dusty accretion torus around the central black hole dominates their appearance. In 'type-1' systems, the bright nucleus is visible at the centre of a face-on torus. In 'type-2' systems the thick, nearly edge-on torus hides the central engine. Later studies suggested evolutionary effects and added dusty clumps and polar winds but left the basic picture intact. However, recent high-resolution images of the archetypal type-2 galaxy NGC 1068 suggested a more radical revision. They displayed a ring-like emission feature which the authors advocated to be hot dust surrounding the black hole at the radius where the radiation from the central engine evaporates the dust. That ring is too thin and too far tilted from edge-on to hide the central engine, and ad hoc foreground extinction is needed to explain the type-2 classification. These images quickly generated reinterpretations of the type 1-2 dichotomy. Here we present new multi-band mid-infrared images of NGC1068 that detail the dust temperature distribution and reaffirm the original model. Combined with radio data, our maps locate the central engine below the previously reported ring and obscured by a thick, nearly edge-on disk, as predicted by the Unified Model. We also identify emission from polar flows and absorbing dust that is mineralogically distinct from that towards the Milky Way centre., Comment: In press at Nature

The Mid-infrared ELT Imager and Spectrograph (METIS) will provide the Extremely Large Telescope (ELT) with a unique window to the thermal- and mid-infrared (3 - 13 microns). Its single-conjugate adaptive optics (SCAO) system will enable high contrast imaging and integral field unit (IFU) spectroscopy (R~100,000) at the diffraction limit of the ELT. This article describes the science drivers, conceptual design, observing modes, and expected performance of METIS., Comment: 5 pages, 3 figures, published in the ESO Messenger

Ground-based thermal-infrared observations have a unique scientific potential, but are also extremely challenging due to the need to accurately subtract the high thermal background. Since the established techniques of chopping and nodding need to be modified for observations with the future mid-infrared ELT imager and spectrograph (METIS), we investigate the sources of thermal background subtraction residuals. Our aim is to either remove or at least minimise the need for nodding in order to increase the observing efficiency for METIS. To this end we need to improve our knowledge about the origin of chop residuals and devise observing methods to remove them most efficiently, i.e. with the slowest possible nodding frequency. Thanks to dedicated observations with VLT/VISIR and GranTeCan/CanariCam, we have successfully traced the origin of three kinds of chopping residuals to (1) the entrance window, (2) the spiders and (3) other warm emitters in the pupil, in particular the VLT M3 mirror cell in its parking position. We conclude that, in order to keep chopping residuals stable over a long time (and therefore allow for slower nodding cycles), the pupil illumination needs to be kept constant, i.e. (imaging) observations should be performed in pupil-stabilised, rather than field-stabilised mode, with image de-rotation in the post-processing pipeline. This is now foreseen as the default observing concept for all METIS imaging modes., Comment: 17 pages, 9 figures, submitted to the SPIE proceedings as paper no. 11447-403

METIS is one of the three first-light instruments planned for the ELT, mainly dedicated to high contrast imaging in the mid-infrared. On the SPHERE high-contrast instrument currently installed at the VLT, we observe that one of the main contrast limitations is the wind driven halo, due to the limited AO running speed with respect to the atmospheric turbulence temporal evolution. From this observation, we extrapolate this signature to the ELT/METIS instrument, which is equipped with a single conjugated adaptive optics system and with several coronagraphic devices. By making use of an analytic AO simulator, we compare the amount of wind driven halo observed with SPHERE and with METIS, under the same turbulence conditions., Comment: AO4ELT6 Conference Proceedings, 12 pages, 8 figures

EX Lup is the prototype of a class of pre-main sequence eruptive stars defined by their repetitive outbursts lasting several months. In 2008 January-September EX Lup underwent its historically largest outburst, brightening by about 4 magnitudes in visual light. In previous studies we discovered on-going silicate crystal formation in the inner disk during the outburst, but also noticed that the measured crystallinity fraction started decreasing after the source returned to the quiescent phase. Here we present new observations of the 10 $\mu$m silicate feature, obtained with the MIDI and VISIR instruments at Paranal Observatory. The observations demonstrate that within five years practically all crystalline forsterite disappeared from the surface of the inner disk. We reconstruct this process by presenting a series of parametric axisymmetric radiative transfer models of an expanding dust cloud that transports the crystals from the terrestrial zone to outer disk regions where comets are supposed to form. Possibly the early Sun also experienced similar flare-ups, and the forming planetesimals might have incorporated crystalline silicate material produced by such outbursts. Finally, we discuss how far the location of the dust cloud could be constrained by future JWST observations., Comment: 12 pages, 4 figures, accepted for publication in ApJ

H$\alpha$ is a powerful tracer of accretion and chromospheric activity, which has been detected in the case of young brown dwarfs and even recently in planetary mass companions (e.g. PDS70 b and c). H$\alpha$ detections and characterisation of brown dwarf and planet companions can further our knowledge of their formation and evolution, and expanding such a sample is therefore our primary goal. We used the Zurich IMaging POLarimeter (ZIMPOL) of the SPHERE instrument at the Very Large Telescope (VLT) to observe the known $38-72$M$_{\mathrm{J}}$ companion orbiting PZ Tel, obtaining simultaneous angular differential imaging observations in both continuum and narrow H$\alpha$ band. We detect H$\alpha$ emission from the companion, making this only the second H$\alpha$ detection of a companion using the SPHERE instrument. We used our newly added astrometric measurements to update the orbital analysis of PZ Tel B, and we used our photometric measurements to evaluate the H$\alpha$ line flux. Given the estimated bolometric luminosity, we obtained an H$\alpha$ activity (log$\mathrm{(L_{H\alpha}/L_{bol}})$) between $-4.16$ and $-4.31$. The H$\alpha$ activity of PZ Tel B is consistent with known average activity levels for M dwarf of the same spectral type. Given the absence of a known gaseous disk and the relatively old age of the system (24 Myr), we conclude that the H$\alpha$ emission around PZ Tel B is likely due to chromospheric activity., Comment: 9 pages, 4 figures, to be published in Astronomy & Astrophysics

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., Comment: White paper for ESA Voyage 2050; 24 pages

Context. The interaction between low-mass companions and the debris discs they reside in is still not fully understood. A debris disc can evolve due to self-stirring, a process in which planetesimals can excite their neighbours to the point of destructive collisions. In addition, the presence of a companion could further stir the disc (companion-stirring). Additional information is necessary to understand this fundamental step in the formation and evolution of a planetary system, and at the moment of writing only a handful of systems are known where a companion and a debris disc have both been detected and studied at the same time. Aims. Our primary goal is to augment the sample of these systems and to understand the relative importance between self-stirring and companion-stirring. Methods. In the course of the VLT/NaCo Imaging Survey for Planets around Young stars (ISPY), we observed HD 193571, an A0 debris disc hosting star at a distance of 68 pc with an age between 60 and 170 Myr. We obtained two sets of observations in L' band and a third epoch in H band using the GPI instrument at Gemini-South. Results. A companion was detected in all three epochs at a projected separation of 11 au (0.17 arcsec), and co-motion was confirmed through proper motion analysis. Given the inferred disc size of 120 au, the companion appears to reside within the gap between the host star and the disc. Comparison between the L' and H band magnitude and evolutionary tracks suggests a mass of 0.31 - 0.39 solar masses . Conclusions. We discovered a previously unknown M-dwarf companion around HD 193571, making it the third low-mass stellar object discovered within a debris disc. A comparison to self- and companion-stirring models suggests that the companion is likely responsible for the stirring of the disc., Comment: 10 pages, 8 figures, accepted for publication in Astronomy&Astrophysics

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction., Comment: Accepted for publication in Experimental Astronomy as part of the topical collection: Future of Optical-infrared Interferometry in Europe

We refine the gap size measurements of the disk surrounding the Herbig Ae star HD 100546 in the N band. Our new mid-infrared interferometric (MIDI) data have been taken with the UT baselines and span the full range of orientations. The correlated fluxes show a wavy pattern in which the minima separation links to a geometrical structure in the disk. We fit each correlated flux measurement with a spline function, deriving the corresponding spatial scale, while assuming that the pattern arises interferometrically due to the bright emission from the inner disk and the opposing sides of the wall of the outer disk. We then fit an ellipse to the derived separations at their corresponding position angles, thereby using the observations to constrain the disk inclination to i =47 +/- 1 degree and the disk position angle to PA =135.0 +/- 2.5 degree East of North, both of which are consistent with the estimated values in previous studies. We also derive the radius of the ellipse to 15.7 +/- 0.8 au. To confirm that the minima separations translate to a geometrical structure in the disk, we model the disk of HD 100546 using a semi-analytical approach taking into account the temperature and optical depth gradients. Using this model, we simultaneously reproduce the level and the minima of the correlated fluxes and constrain the gap size of the disk for each observation. The values obtained for the projected gap size in different orientations are consistent with the separation found by the geometrical model., Comment: Accepted for publication in Astrophysical Journal. 10 pages, 4 figures, 3 tables

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission., Comment: PASP in press

Spatially resolving the immediate surroundings of young stars is a key challenge for the planet formation community. SPHERE on the VLT represents an important step forward by increasing the opportunities offered by optical or near-infrared imaging instruments to image protoplanetary discs. The Guaranteed Time Observation Disc team has concentrated much of its efforts on polarimetric differential imaging, a technique that enables the efficient removal of stellar light and thus facilitates the detection of light scattered by the disc within a few au from the central star. These images reveal intriguing complex disc structures and diverse morphological features that are possibly caused by ongoing planet formation in the disc. An overview of the recent advances enabled by SPHERE is presented., Comment: 6 pages, 5 figures, published by ESO Messenger

We conducted a detailed radiative transfer modeling of the dust emission from the circumstellar disk around HL Tau. The goal of our study is to derive the surface density profile of the inner disk and its structure. In addition to the Atacama Large Millimeter/submillimeter Array images at Band 3 (2.9mm), Band 6 (1.3mm), and Band 7 (0.87mm), the most recent Karl G. Jansky Very Large Array (VLA) observations at 7mm were included in the analysis. A simulated annealing algorithm was invoked to search for the optimum model. The radiative transfer analysis demonstrates that most radial components (i.e., >6AU) of the disk become optically thin at a wavelength of 7mm, which allows us to constrain, for the first time, the dust density distribution in the inner region of the disk. We found that a homogeneous grain size distribution is not sufficient to explain the observed images at different wavelengths simultaneously, while models with a shallower grain size distribution in the inner disk work well. We found clear evidence that larger grains are trapped in the first bright ring. Our results imply that dust evolution has already taken place in the disk at a relatively young (i.e., ~1Myr) age. We compared the midplane temperature distribution, optical depth, and properties of various dust rings with those reported previously. Using the Toomre parameter, we briefly discussed the gravitational instability as a potential mechanism for the origin of the dust clump detected in the first bright ring via the VLA observations., Comment: Accepted for publication in A&A (10 pages)

High-resolution imaging reveals a large morphological variety of protoplanetary disks. To date, no constraints on their global evolution have been found from this census. An evolutionary classification of disks was proposed based on their IR spectral energy distribution, with the Group I sources showing a prominent cold component ascribed to an earlier stage of evolution than Group II. Disk evolution can be constrained from the comparison of disks with different properties. A first attempt of disk taxonomy is now possible thanks to the increasing number of high-resolution images of Herbig Ae/Be stars becoming available. Near-IR images of six Group II disks in scattered light were obtained with VLT/NACO in Polarimetric Differential Imaging, which is the most efficient technique to image the light scattered by the disk material close to the stars. We compare the stellar/disk properties of this sample with those of well-studied Group I sources available from the literature. Three Group II disks are detected. The brightness distribution in the disk of HD163296 indicates the presence of a persistent ring-like structure with a possible connection with the CO snowline. A rather compact (less than 100 AU) disk is detected around HD142666 and AK Sco. A taxonomic analysis of 17 Herbig Ae/Be sources reveals that the difference between Group I and Group II is due to the presence or absence of a large disk cavity (larger than 5 AU). There is no evidence supporting the evolution from Group I to Group II. Group II are not evolved version of the Group I. Within the Group II disks, very different geometries (both self-shadowed and compact) exist. HD163296 could be the primordial version of a typical Group I. Other Group II, like AK Sco and HD142666, could be smaller counterpart of Group I unable to open cavities as large as those of Group I., Comment: 16 pages, 7 figures, published by A&A

The James Webb Space Telescope will enable astronomers to obtain exoplanet spectra of unprecedented precision. Especially the MIRI instrument may shed light on the nature of the cloud particles obscuring planetary transmission spectra in the optical and near-infrared. We provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups. We select targets which span a wide range in planetary temperature and surface gravity, ranging from super-Earths to giant planets, and have a high expected SNR. For all targets we vary the enrichment, C/O ratio, presence of optical absorbers (TiO/VO) and cloud treatment. We calculate atmospheric structures and emission and transmission spectra for all targets and use a radiometric model to obtain simulated observations. We analyze JWST's ability to distinguish between various scenarios. We find that in very cloudy planets such as GJ 1214b less than 10 transits with NIRSpec may be enough to reveal molecular features. Further, the presence of small silicate grains in atmospheres of hot Jupiters may be detectable with a single JWST MIRI transit. For a more detailed characterization of such particles less than 10 transits are necessary. Finally, we find that some of the hottest hot Jupiters are well fitted by models which neglect the redistribution of the insolation and harbor inversions, and that 1-4 eclipse measurements with NIRSpec are needed to distinguish between the inversion models. Wet thus demonstrate the capabilities of JWST for solving some of the most intriguing puzzles in current exoplanet atmospheric research. Further, by publishing all models calculated for this study we enable the community to carry out similar or retrieval analyses for all planets included in our target list., Comment: 24 pages, 7 figures, accepted for publication in A&A

We present scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 0.63, 0.79, 1.24 and 1.62 micron. We also present H2/H3-band ADI observations. Three distinct radial depressions in the polarized intensity distribution are seen, around 85, 21, and 6~au. The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average towards the South and darker towards the North-West. The ADI observations yielded no signifiant detection of point sources in the disk. Our observations have a linear spatial resolution of 1 to 2au, similar to that of recent ALMA dust continuum observations. The sub-micron sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas. We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling. This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations. We find that the gas surface density in the "gaps" is reduced by 50% to 80% relative to an unperturbed model. Should embedded planets be responsible for carving the gaps then their masses are at most a few 10 Mearth. The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets. If forming planetary bodies have undergone collapse and are in the "detachted phase" then they may be directly observable with future facilities such as METIS at the E-ELT., Comment: 24 pages, 20 figures, accepted for publication in the ApJ

The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with JWST and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) Envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O<1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed., Comment: 32 pages, 18 figures. Accepted for publication in ApJ

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings., Comment: Accepted by The Astrophysical Journal Letters

We present new mid-infrared interferometric observations of the massive young stellar object IRAS 13481-6124, using VLTI/MIDI for spectrally-resolved, long-baseline measurements (projected baselines up to $\sim120$ m) and GSO/T-ReCS for aperture-masking interferometry in five narrow-band filters (projected baselines of $\sim1.8-6.4$ m) in the wavelength range of $7.5-13$ $\mu$m. We combine these measurements with previously-published interferometric observations in the $K$ and $N$ bands in order to assemble the largest collection of infrared interferometric observations for a massive YSO to date. Using a combination of geometric and radiative-transfer models, we confirm the detection at mid-infrared wavelengths of the disk previously inferred from near-infrared observations. We show that the outflow cavity is also detected at both near- and mid-infrared wavelengths, and in fact dominates the mid-infrared emission in terms of total flux. For the disk, we derive the inner radius ($\sim1.8$ mas or $\sim6.5$ AU at 3.6 kpc), temperature at the inner rim ($\sim1760$ K), inclination ($\sim48$ deg) and position angle ($\sim107$ deg). We determine that the mass of the disk cannot be constrained without high-resolution observations in the (sub-)millimeter regime or observations of the disk kinematics, and could be anywhere from $\sim10^{-3}$ to $20$ M$_\odot$. Finally, we discuss the prospects of interpreting the spectral energy distributions of deeply-embedded massive YSOs, and warn against attempting to infer disk properties from the SED., Comment: Accepted to Astronomy and Astrophysics

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star., Comment: 6 pages, 4 figures, accepted for publication in ApJL

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen., Comment: 30 pages, 20 figures. Accepted for publication in ApJ

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates., Comment: Accepted for publication in ApJ; 24 pages, 16 figures

We aim to construct a spectral energy distribution (SED) for the emission from the dayside atmosphere of the hot Jupiter WASP-46b and to investigate its energy budget. We observed a secondary eclipse of WASP-46b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2m telescope. Eclipse depths of the acquired light curves were derived to infer the brightness temperatures at multibands that cover the SED peak. We report the first detection of the thermal emission from the dayside of WASP-46b in the K band at 4.2-sigma level and tentative detections in the H (2.5-sigma) and J (2.3-sigma) bands, with flux ratios of 0.253 +0.063/-0.060%, 0.194 +/- 0.078%, and 0.129 +/- 0.055%, respectively. The derived brightness temperatures (2306 +177/-187K, 2462 +245/-302K, and 2453 +198/-258K, respectively) are consistent with an isothermal temperature profile of 2386K, which is significantly higher than the dayside-averaged equilibrium temperature, indicative of very poor heat redistribution efficiency. We also investigate the tentative detections in the g'r'i' bands and the 3-sigma upper limit in the z' band, which might indicate the existence of reflective clouds if these tentative detections do not arise from systematics., Comment: 8 pages, 7 figures, accepted for publication in A&A

(Abridged) WASP-5b is a highly irradiated dense hot Jupiter orbiting a G4V star every 1.6 days. We observed two secondary eclipses of WASP-5b in the J, H and K bands simultaneously. Thermal emission of WASP-5b is detected in the J and K bands. The retrieved planet-to-star flux ratios in the J and K bands are 0.168 +0.050/-0.052% and 0.269+/-0.062%, corresponding to brightness temperatures of 2996 +212/-261K and 2890 +246/-269K, respectively. No thermal emission is detected in the H band, with a 3-sigma upper limit of 0.166%, corresponding to a maximum temperature of 2779K. On the whole, our J, H, K results can be explained by a roughly isothermal temperature profile of ~2700K in the deep layers of the planetary dayside atmosphere that are probed at these wavelengths. Together with Spitzer observations, which probe higher layers that are found to be at ~1900K, a temperature inversion is ruled out in the range of pressures probed by the combined data set. While an oxygen-rich model is unable to explain all the data, a carbon-rich model provides a reasonable fit but violates energy balance., Comment: 13 pages, 9 figures, accepted for publication in A&A

(Abridged) We observed one transit and one occultation of the hot Jupiter WASP-43b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2-meter telescope. From the transit event, we have independently derived WASP-43's system parameters with high precision, and improved the period to be 0.81347437(13) days. No significant variation in transit depths is detected, with the largest deviations coming from the i', H, and K bands. Given the observational uncertainties, the broad-band transmission spectrum can be explained by either a flat featureless straight line that indicates thick clouds, synthetic spectra with absorption signatures of atomic Na/K or molecular TiO/VO that indicate cloud-free atmosphere, or a Rayleigh scattering profile that indicates high-altitude hazes. From the occultation event, we have detected planetary dayside thermal emission in the K-band with a flux ratio of 0.197 +/- 0.042%, which confirms previous detections obtained in the 2.09 micron narrow band and Ks-band. The K-band brightness temperature 1878 +108/-116 K favors an atmosphere with poor day- to night-side heat redistribution. We also have a marginal detection in the i'-band (0.037 +0.023/-0.021%), which is either a false positive, a signature of non-blackbody radiation at this wavelength, or an indication of reflective hazes at high altitude., Comment: 14 pages, 9 figures, accepted for publication in A&A

We aim to characterize the distribution and composition of circumstellar material around young massive stars, and to investigate exactly which physical structures in these objects are probed by long-baseline mid-infrared interferometric observations. We used the two-telescope interferometric instrument MIDI of the Very Large Telescope Interferometer of the European Southern Observatory to observe a sample of 24 intermediate- and high-mass young stellar objects in the N band (8-13 micron). We had successful fringe detections for 20 objects, and present spectrally-resolved correlated fluxes and visibility levels for projected baselines of up to 128 m. We fit the visibilities with geometric models to derive the sizes of the emitting regions, as well as the orientation and elongation of the circumstellar material. Fourteen objects in the sample show the 10 micron silicate feature in absorption in the total and correlated flux spectra. For 13 of these objects, we were able to fit the correlated flux spectra with a simple absorption model, allowing us to constrain the composition and absorptive properties of the circumstellar material. Nearly all of the massive young stellar objects observed show significant deviations from spherical symmetry at mid-infrared wavelengths. In general, the mid-infrared emission can trace both disks and outflows, and in many cases it may be difficult to disentangle these components on the basis of interferometric data alone, because of the sparse spatial frequency coverage normally provided by current long-baseline interferometers. For the majority of the objects in this sample, the absorption occurs on spatial scales larger than those probed by MIDI. Finally, the physical extent of the mid-infrared emission around these sources is correlated with the total luminosity, albeit with significant scatter., Comment: 36 pages, 22 figures. Accepted to Astronomy and Astrophysics

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least \alpha{} >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and \alpha{} = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least \alpha{} = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap., Comment: accepted by A&A. 10 pages, 10 figures, 2 tables

We report new detections of thermal emission from the transiting hot Jupiter WASP-43b in the H and Ks-bands as observed at secondary eclipses. The observations were made with the WIRCam instrument on the CFHT. We obtained a secondary eclipse depth of 0.103$_{-0.017}^{+0.017}%$ and 0.194$_{-0.029}^{+0.029}%$ in the H and Ks-bands, respectively. The Ks band depth is consistent with previous measurement in the narrow band centered at 2.09um by Gillon et al. (2012). Our eclipse depths in both bands are consistent with a blackbody spectrum with a temperature of ~1850 K, slightly higher than the dayside equilibrium temperature without day-night energy redistribution. Based on theoretical models of the dayside atmosphere of WASP-43b, our data constrain the day-night energy redistribution in the planet to be $\lesssim 15-25$%, depending on the metal content in the atmosphere. Combined with energy balance arguments our data suggest that a strong temperature inversion is unlikely in the dayside atmosphere of WASP-43b. However, a weak inversion cannot be strictly ruled out at the current time. Future observations are required to place detailed constraints on the chemical composition of the atmosphere., Comment: 7 pages, 6 figures, to be published in ApJ

We investigate the young stellar objects (YSOs) in the Lynds~1641 (L1641) cloud using multi-wavelength data including Spitzer, WISE, 2MASS, and XMM covering ~1390 YSOs across a range of evolutionary stages. In addition, we targeted a sub-sample of YSOs for optical spectroscopy. We use this data, along with archival photometric data, to derive spectral types, extinction values, masses, ages, as well as accretion rates. We obtain a disk fraction of ~50% in L1641. The disk frequency is almost constant as a function of stellar mass with a slight peak at log(M_*/M_sun)\approx-0.25. The analysis of multi-epoch spectroscopic data indicates that the accretion variability of YSOs cannot explain the two orders of magnitude of scatter for YSOs with similar masses. Forty-six new transition disk (TD) objects are confirmed in this work, and we find that the fraction of accreting TDs is lower than for optically thick disks (40--45% vs. 77--79% respectively). We confirm our previous result that the accreting TDs have a similar median accretion rate to normal optically thick disks. We confirm that two star formation modes (isolated vs. clustered) exist in L1641. We find that the diskless YSOs are statistically older than the YSOs with optically-thick disks and the transition disk objects have a median age which is intermediate between the two populations. We tentatively study the star formation history in L1641 based on the age distribution and find that star formation started to be active 2--3 Myr ago., Comment: 112 pages, 46 figures, 6 tables, accepted by ApJs

Deeply embedded and at distances of several kiloparsecs, massive young stellar objects (MYSOs) present numerous challenges for observation and study. In this work, we present spatially-resolved observations of one MYSO, AFGL 4176, together with survey and literature data, ranging from interferometric observations with VLTI/MIDI in the mid-infrared, to single-dish Herschel measurements in the far-infrared, and sub-millimeter data from APEX. We consider this spatially-resolved, multi-wavelength data set in terms of both radiative transfer and geometric models. We find that the observations are well described by one-dimensional models overall, but there are also substantial deviations from spherical symmetry at scales of tens to hundreds of astronomical units, which are revealed by the mid-infrared interferometric measurements. We use a multiple-component, geometric modeling approach to explain the mid-infrared emission on scales of tens to hundreds of astronomical units, and find the MIDI measurements are well described by a model consisting of a one-dimensional Gaussian halo and an inclined (\theta=60 deg) circumstellar disk extending out to several hundred astronomical units along a position angle of 160 deg. Finally, we compare our results both with previous models of this source, and with those of other MYSOs, and discuss the present situation with mid-infrared interferometric observations of massive stars., Comment: 15 pages, 14 figures. Accepted to Astronomy and Astrophysics

Information about the inner structure of young stellar objects is crucial for understanding how the central forming stars gain their mass. However, especially for observations of (usually further away) high-mass young stellar objects, conventional imaging has limitations in spatial resolution. By means of mid-infrared interferometry, we can peek deeply into the strongly extinguished central 100 AU of such objects. Here, we report on new data we obtained within our programme using MIDI at the VLTI. Shown are preliminary results on the known outflow sources AFGL 2136 IRS 1 and Mon R2 IRS 3A. In particular, we describe how quantities like differential phases give additional geometrical structure information. We demonstrate how the combined interpretation of single-telescope and interferometric data at different wavelength regimes can lead to a more complete picture of the nature of such MYSOs., Comment: 6 pages, 8 figure files; to appear in the ASP conference proceedings of the Magdalena Ridge Observatory Interferometry Workshop "Resolving the Future of Astronomy with Long-Baseline Interferometry" (eds. M. Creech-Eakman et al.), held in Socorro, NM on March 28-31, 2011

We present spatially-resolved mid-infrared interferometric observations of the massive young stellar object AFGL 4176, together with literature and survey data. We model these observations using a simple, spherically-symmetric radiative transfer model, and find that the observational data are consistent with a highly luminous star surrounded by a thick envelope., Comment: To appear in the ASP conference proceedings of the Magdalena Ridge Observatory Interferometry Workshop "Resolving the Future of Astronomy with Long-Baseline Interferometry" (eds. M. Creech-Eakman et al.), held in Socorro, NM on March 28-31, 2011

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra., Comment: Accepted for publication in ApJ, 37 pages, 11 figures, 2 tables

(Note: this is a shortened version of the original A&A-style structured abstract). The physical nature of the strong photometric variability of T Tau Sa, the more massive member of the Southern "infrared companion" to T Tau, has long been debated. Intrinsic luminosity variations due to variable accretion were originally proposed but later challenged in favor of apparent fluctuations due to time-variable foreground extinction. In this paper we use the timescale of the variability as a diagnostic for the underlying physical mechanism. Because the IR emission emerging from Sa is dominantly thermal emission from circumstellar dust at <=1500K, we can derive a minimum size of the region responsible for the time-variable emission. In the context of the variable foreground extinction scenario, this region must be (un-) covered within the variability timescale, which implies a minimum velocity for the obscuring foreground material. If this velocity supercedes the local Kepler velocity we can reject foreground extinction as a valid variability mechanism. The variable accretion scenario allows for shorter variability timescales since the variations in luminosity occur on much smaller scales, essentially at the surface of the star, and the disk surface can react almost instantly on the changing irradiation with a higher or lower dust temperature and according brightness. We have detected substantial variations at long wavelengths in T Tau S: +26% within four days at 12.8 micron. We show that this short-term variability cannot be due to variable extinction and instead must be due to variable accretion. Using a radiative transfer model of the Sa disk we show that variable accretion can in principle also account for the much larger (several magnitude) variations observed on timescales of several years. For the long-term variability, however, also variable foreground extinction is a viable mechanism., Comment: 15 pages, 8 figures, Accepted for publication in Astronomy and Astrophysics

This contribution reviews results from interferometric observations of circumstellar disks around young stars of $\lesssim$3 M$_{\odot}$, performed with the MIDI instrument operating in the 10 $\mu$m spectral region. Two main topics, the disk structure on $\sim$1-10 AU scales and the dust properties in the same region, are illustrated with several examples of MIDI studies, covering various evolutionary stages. The spatially resolved observations largely confirm SED-only based hypotheses on disk structure, yet also reveal degeneracies that may occur in such SED modeling. The properties of the dust on the disk surface show a strong radial dependence: the dust close to the central star has generally larger grain sizes and in particular a much higher crystallinity than the dust in more remote disk regions., Comment: 13 pages, 7 figures, invited review, conference "The Universe under the Microscope - Astrophysics at High Angular Resolution", Bad Honnef, April 2008