Search

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

Context. Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions (~0.1–10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions (~1–10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Aims. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. Methods. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry. Results. Our temperature-gradient model can consistently reproduce our dataset. We confirmed a spatially extended inner 10 au emission in H- and L-bands, with a homogeneously high temperature (~1700 K), which we associate with the presence of stochastically heated nano-grains. On the other hand, the N-band emitting region presents a ring-like geometry that starts at about 10 au with a temperature of 400 K. Moreover, the existing low resolution MATISSE data exclude the presence of aromatic carbon grains (i.e., producing the 3.3 μm feature) in close proximity tothe star (≲1 au). Future medium spectral resolution MATISSE data will confirm their presence at larger distances. Conclusions. Our best-fit model demonstrates the presence of two separated dust populations: nano-grains that dominate the near- to mid-IR emission in the inner 10 au region and larger grains that dominate the emission outward. The presence of such nano-grains in the highly irradiated inner 10 au region of HD 179218 requires a replenishment process. Considering the expected lifetime of carbon nano-grains from The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS model), the estimated disk accretion inflow of HD 179218 could significantly contribute to feed the inner 10 au region in nano-grains.Moreover, we also expect a local regeneration of those nano-grains by the photo-fragmentation of larger aggregates.

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r, accepted for publication in A&A

The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars., Comment: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

These protoplanetary disks in T Tauri stars play a central role in star and planet formation. We spatially resolve at sub-au scales the innermost regions of a sample of T Tauri's disks to better understand their morphology and composition. We extended our homogeneous data set of 27 Herbig stars and collected near-IR K-band observations of 17 T Tauri stars, spanning effective temperatures and luminosities in the ranges of ~4000-6000 K and ~0.4-10 Lsun. We focus on the continuum emission and develop semi-physical geometrical models to fit the interferometric data and search for trends between the properties of the disk and the central star. The best-fit models of the disk's inner rim correspond to wide rings. We extend the Radius-luminosity relation toward the smallest luminosities (0.4-10 Lsun) and find the R~L^(1/2) trend is no longer valid, since the K-band sizes measured with GRAVITY are larger than the predicted sizes from sublimation radius computation. No clear correlation between the K-band half-flux radius and the mass accretion rate is seen. Having magnetic truncation radii in agreement with the K-band GRAVITY sizes would require magnetic fields as strong as a few kG, which should have been detected, suggesting that accretion is not the main process governing the location of the half-flux radius of the inner dusty disk. Our measurements agree with models that take into account the scattered light. The N-to-K band size ratio may be a proxy for disentangling disks with silicate features in emission from disks with weak and/or in absorption silicate features. When comparing inclinations and PA of the inner disks to those of the outer disks (ALMA) in nine objects of our sample, we detect misalignments for four objects., Comment: Accepted for publication in A&A; 20 pages, 15 figures, 9 tables

The formation and evolution of planetary systems impact the primordial accretion disk. HD141569 is the only known pre-main sequence star characterized by a hybrid disk. Observations probed the outer-disk structure showing a complex system of rings and interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. The goal of this work was to explore with new high-resolution interferometric observations the properties of the dust and gas in the internal regions of HD141569. We observed HD141569 on mas scales with GRAVITY/VLTI in the near-infrared at low and high spectral resolution. We interpreted the visibilities and spectral energy distribution with geometrical models and radiative transfer techniques to constrain the dust emission. We analyzed the high spectral resolution quantities to investigate the properties of the Br-Gamma line emitting region. Thanks to the combination of three different epochs, GRAVITY resolves the inner dusty disk in the K band. Data modeling shows that an IR excess of about 6% is spatially resolved and that the origin of this emission is confined in a ring of material located at a radius of 1 au from the star with a width smaller than 0.3 au. The MCMax modeling suggests that this emission could originate from a small amount of QHPs, while large silicate grain models cannot reproduce at the same time the observational constraints on the properties of near-IR and mid-IR fluxes. The differential phases in the Br-Gamma line clearly show an S-shape that can be best reproduced witha gas disk in Keplerian rotation, confined within 0.09 au. This is also hinted at by the double-peaked Br-Gamma emission line shape. The modeling of the continuum and gas emission shows that the inclination and position angle of these two components are consistent with a system showing relatively coplanar rings on all scales., Comment: Accepted for publication in A&A; 25 pages, 15 figures, 5 tables

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim., Comment: accepted for publication in A&A

We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the Spectral Energy Distribution and match the interferometric observations. To match the SED , our model requires a stellar luminosity of 2.5 Lsun, higher than any previously determined values. Such a high value is needed to accommodate the circumstellar extinction caused by the highly inclined disk, which has been neglected in previous studies. While using an effective temperature of 4800 K determined through high-resolution spectroscopy, we derive a stellar radius of 2.29 Rsun. These revised fundamental parameters, when combined with the mass estimates available , lead to an age of 0.5-2.0 Ma for RY Lup, in better agreement with the age of the Lupus association than previous determinations. Our disk model nicely reproduces the interferometric GRAVITY data and is in good agreement with the PIONIER ones. We derive an inner rim location at 0.12~au from the central star. This model corresponds to an inclination of the inner disk of 50deg, which is in mild tension with previous determinations of a more inclined outer disk from SPHERE (70 deg in NIR) and ALMA(67 $\pm$5 deg) images, but consistent with the inclination determination from the ALMA CO spectra (55$\pm$5deg). Increasing the inclination of the inner disk to 70 deg leads to a higher line-of-sight extinction and therefore requires a higher stellar luminosity of 4.65 Lsun to match the observed flux levels. This luminosity would translate to a stellar radius of 3.13~Rsun, leading to an age of 2-3~Ma, and a stellar mass of about 2 Msun, in disagreement with the observed dynamical mass estimate of 1.3-1.5 Msun. Critically, this high-inclination inner disk model also fails to reproduce the visibilities observed with GRAVITY., Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results., Comment: Accepted for publication in A&A letters

Protoplanetary disks drive some of the formation process (e.g., accretion, gas dissipation, formation of structures, etc.) of stars and planets. Understanding such physical processes is one of the main astrophysical questions. HD 163296 is an interesting young stellar object for which infrared and sub-millimeter observations have shown a prominent circumstellar disk with gaps plausibly created by forming planets. This study aims at characterizing the morphology of the inner disk in HD 163296 with multi-epoch near-infrared interferometric observations performed with GRAVITY at the Very Large Telescope Interferometer (VLTI). Our goal is to depict the K-band (lambda_0 ~ 2.2 um) structure of the inner rim with milliarcsecond (sub-au) angular resolution. Our data is complemented with archival PIONIER (H-band; lambda_0 ~ 1.65 um) data of the source. We performed a Gradient Descent parametric model fitting to recover the sub-au morphology of our source. Our analysis shows the existence of an asymmetry in the disk surrounding the central star of HD 163296. We confirm variability of the disk structure in the inner ~2 mas (0.2 au). While variability of the inner disk structure in this source has been suggested by previous interferometric studies, this is the first time that it is confirmed in the H- and K-bands by using a complete analysis of the closure phases and squared visibilities over several epochs. Because of the separation from the star, position changes, and persistence of this asymmetric structure on timescales of several years, we argue that it is a dusty feature (e.g., a vortex or dust clouds), probably, made by a mixing of sillicate and carbon dust and/or refractory grains, inhomogeneously distributed above the mid-plane of the disk., Accepted to be published in Astronomy and Astrophysics; main-body: 11 pages, 3 figures and 3 tables

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation., Comment: Accepted for publication in A&A; 23 pages, 16 figures, 7 tables

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

13 pags., 18 figs., 2 tabs., Aims. This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. Methods. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ice features at 45 and 63 μm change with key disk properties: disk size, grain sizes, disk dust mass, dust settling, and ice thickness. In addition, a series of models is devised to calculate the water ice emission features from warmup, direct deposit, and cooldown scenarios of the water ice in disks. Results. Photodesorption from icy grains in disk surfaces weakens the mid-IR water ice features by factors of 4-5. The far-IR water ice emission features originate from small grains at the surface snow line in disks at distance of 10-100 au. Unless this reservoir is missing in disks (e.g., transitional disks with large cavities), the feature strength does not change. Grains larger than 10 μm do not contribute to the features. Grain settling (using turbulent description) affects the strength of the ice features by at most 15%. The strength of the ice feature scales with the disk dust mass and water ice fraction on the grains, but saturates for dust masses higher than 10 M and for ice mantles that increase the dust mass by more than 50%. The various thermal histories of water ice leave an imprint on the shape of the features (crystalline and/or amorphous) and on the peak strength and position of the 45 μm feature. SOFIA/HIRMES can only detect crystalline ice features that are much stronger than those simulated in our standard T Tauri disk model in deep exposures (1 h). SPICA/SAFARI can detect the typical ice features in our standard T Tauri disk model in short exposures (10 min). Conclusions. The sensitivity of SPICA/SAFARI will allow the detailed study of the 45 and 63 μm water ice feature in unbiased surveys of T Tauri stars in nearby star forming regions and an estimate of the mass of their ice reservoir. The water ice emission features carry an imprint of the thermal history of the ice, and thus can distinguish between various formation and transport scenarios. Amorphous ice at 45 μm that has a much broader and flatter peak could be detected in deep surveys if the underlying continuum can be well characterized and the baseline stability of SAFARI is better than a few percent., IK and MM acknowledge funding from the EU FP7-2011 under Grant Agreement no. 284405. LK is supported by a grant from the Netherlands Research School for Astronomy (NOVA).

To understand the chemical composition of planets, it is important to know the chemical composition of the region where they form in protoplanetary disks. Due to its fundamental role in chemical and biological processes, carbon is a key element to trace. We aim to identify the carriers and processes behind the extended NIR flux observed around several Herbig stars. We compare the extended NIR flux from objects in the PIONIER Herbig Ae/Be survey with their flux in the PAH features. HD 100453 is used as a benchmark case to investigate the influence of quantum heated particles, like PAHs or very small carbonaceous grains, in more detail. We use the Monte Carlo radiative transfer code MCMax to do a parameter study of the QHP size and scale- height and examine the influence of quantum heating on the amount of extended flux in the NIR visibilities. There is a correlation between the PAH feature flux of a disk and the amount of its extended NIR flux. We find that very small carbonaceous grains create the observed extended NIR flux around HD 100453 and still lead to a realistic SED. These results can not be achieved without using quantum heating effects, e.g. only with scattered light and grains in thermal equilibrium. It is possible to explain the extended NIR emission around Herbig stars with the presence of carbonaceous, quantum heated particles. Interferometric observations can be used to constrain the spatial distribution and typical size of carbonaceous material in the terrestrial planet forming region., 10 pages, 6 figures, 2 tables, accepted for publication at A&A

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks., Accepted for publication in A&A (including acknowledgments)