Showing 1–50 of 262 results for author: Dominik, C

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): PDS 111, an old T Tauri star with a young-looking disk

Authors: Annelotte Derkink, Christian Ginski, Paola Pinilla, Nicolas Kurtovic, Lex Kaper, Alex de Koter, Per-Gunnar Valegård, Eric Mamajek, Frank Backs, Myriam Benisty, Til Birnstiel, Gabriele Columba, Carsten Dominik, Antonio Garufi, Michiel Hogerheijde, Rob van Holstein, Jane Huang, François Ménard, Christian Rab, María Claudia Ramírez-Tannus, Álvaro Ribas, Jonathan P. Williams, Alice Zurlo

Abstract: The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary… ▽ More The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$α$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 23 pages, 23 figures, accepted by A&A (abstract shortened)

SPHERE RefPlanets: Search for epsilon Eridani b and warm dust

Authors: C. Tschudi, H. M. Schmid, M. Nowak, H. Le Coroller, S. Hunziker, R. G. van Holstein, C. Perrot, D. Mouillet, J. -C. Augereau, A. Bazzon, J. L. Beuzit, A. Boccaletti, M. J. Bonse, G. Chauvin, S. Desidera, K. Dohlen, C. Dominik, N. Engler, M. Feldt, J. H. Girard, R. Gratton, Th. Henning, M. Kasper, P. Kervella, A. -M. Lagrange , et al. (13 additional authors not shown)

Abstract: We carried out very deep VLT/SPHERE imaging polarimetry of the nearby system Eps Eri based on 38.5 hours of integration time with a 600 - 900 nm broadband filter to search for polarized scattered light from a planet or from circumstellar dust using AO, coronagraphy, high precision differential polarimetry, and angular differential imaging. We have improved several data reduction and post-processin… ▽ More We carried out very deep VLT/SPHERE imaging polarimetry of the nearby system Eps Eri based on 38.5 hours of integration time with a 600 - 900 nm broadband filter to search for polarized scattered light from a planet or from circumstellar dust using AO, coronagraphy, high precision differential polarimetry, and angular differential imaging. We have improved several data reduction and post-processing techniques and also developed new ones to further increase the sensitivity of SPHERE/ZIMPOL. The data provide unprecedented contrast limits, but no significant detection of a point source or an extended signal from circumstellar dust. For each observing epoch, we obtained a point source contrast for the polarized intensity between $2\cdot 10^{-8}$ and $4\cdot 10^{-8}$ at the expected separation of the planet Eps Eri b of 1'' near quadrature phase. The polarimetric contrast limits are about six to 50 times better than the intensity limits because polarimetric imaging is much more efficient in speckle suppression. Combining the entire 14-month data set to the search for a planet moving on a Keplerian orbit with the K-Stacker software further improves the contrast limits by a factor of about two, to about $8 \cdot 10^{-9}$ at 1''. This would allow the detection of a planet with a radius of about 2.5 Jupiter radii. The surface brightness contrast limits achieved for the polarized intensity from an extended scattering region are about 15 mag arcsec$^{-2}$ at 1'', or up to 3 mag arcsec$^{-2}$ deeper than previous limits. For Eps Eri, these limits exclude the presence of a narrow dust ring and they constrain the dust properties. This study shows that the polarimetric contrast limits for reflecting planets with SPHERE/ZIMPOL can be improved to a level $<10^{-8}$ simply by collecting more data over many nights and using the K-Stacker software. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 26 pages, 20 figures, 5 tables. Accepted for publication in A&A

A broken debris cascade as a possible source of hot dust emission in transitioning planet-forming disks

Authors: Niels Swinkels, Carsten Dominik

Abstract: Planet-forming disks turn from gas-rich, massive disks made of dust and gas into planetary systems containing only small amounts dust produced by collisions between smaller planetary objects like planetesimals, asteroids, or comets. Traditionally we talk about protoplanetary (age $\sim$1 Myr), transitional ($\sim$ 5-10 Myr), and debris disks ($\sim$ 10-hundreds of Myr), even though the overlap bet… ▽ More Planet-forming disks turn from gas-rich, massive disks made of dust and gas into planetary systems containing only small amounts dust produced by collisions between smaller planetary objects like planetesimals, asteroids, or comets. Traditionally we talk about protoplanetary (age $\sim$1 Myr), transitional ($\sim$ 5-10 Myr), and debris disks ($\sim$ 10-hundreds of Myr), even though the overlap between these phases may be relevant. We aim to show that in the transition phase of a disk, when gas surface densities are reduced but not yet negligible, a seemingly small amount of collisional activity may lead to the production of dust on a level that is observationally relevant by creating regions of optical depth 1 or above. In particular, we aim to show that the hot dust emission component of transitional disks may in fact be debris dust produced in such collisions. We develop an analytical model to derive the conditions under which observationally relevant amounts of dust can be produced. We focus on the effect of the gas surface density during the transition phase of a disk from fully gas-dominated to the gas-poor debris stage. We show that the decrease of the gas surface density has an important effect. It allows smaller planetesimals to become collisional, starting a cascade. At the same time, small particles are not destroyed by collisions. That interrupted cascade is critical to preserve the produced dust for significant time intervals, allowing a seemingly minor amount of planetesimal collisions to be effective in producing detectable amounts of dust. The warm emission of low amounts of dust in many transitional planet-forming disks might be caused entirely by second-generation dust, exposing planetary material at a time much earlier than the age of what are traditionally considered as debris disks. △ Less

Submitted 7 May, 2024; v1 submitted 27 April, 2024; originally announced April 2024.

Comments: To be published in Astronomy & Astrophysics

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

The SPHERE view of the Taurus star-forming region

Authors: A. Garufi, C. Ginski, R. G. van Holstein, M. Benisty, C. F. Manara, S. Pérez, P. Pinilla, Á. Ribas, P. Weber, J. Williams, L. Cieza, C. Dominik, S. Facchini, J. Huang, A. Zurlo, J. Bae, J. Hagelberg, Th. Henning, M. R. Hogerheijde, M. Janson, F. Ménard, S. Messina, M. R. Meyer, C. Pinte, S. P. Quanz , et al. (9 additional authors not shown)

Abstract: The sample of planet-forming disks observed by high-contrast imaging campaigns over the last decade is mature enough to enable the demographical analysis of individual star-forming regions. We present the full census of Taurus sources with VLT/SPHERE polarimetric images available. The whole sample sums up to 43 targets (of which 31 have not been previously published) corresponding to one-fifth of… ▽ More The sample of planet-forming disks observed by high-contrast imaging campaigns over the last decade is mature enough to enable the demographical analysis of individual star-forming regions. We present the full census of Taurus sources with VLT/SPHERE polarimetric images available. The whole sample sums up to 43 targets (of which 31 have not been previously published) corresponding to one-fifth of the Class II population in Taurus and about half of such objects that are observable. A large fraction of the sample is apparently made up of isolated faint disks (equally divided between small and large self-shadowed disks). Ambient signal is visible in about one-third of the sample. This probes the interaction with the environment and with companions or the outflow activity of the system. The central portion of the Taurus region almost exclusively hosts faint disks, while the periphery also hosts bright disks interacting with their surroundings. The few bright disks are found around apparently older stars. The overall picture is that the Taurus region is in an early evolutionary stage of planet formation. Yet, some objects are discussed individually, as in an intermediate or exceptional stage of the disk evolution. This census provides a first benchmark for the comparison of the disk populations in different star forming regions. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted by A&A

The SPHERE view of the Orion star-forming region

Authors: P. -G. Valegard, C. Ginski, A. Derkink, A. Garufi, C. Dominik, A. Ribas, J. P. Williams, M. Benisty, T. Birnstiel, S. Facchini, G. Columba, M. Hogerheijde, R. G. Van Holstein, J. Huang, M. Kenworthy, C. F. Manara, P. Pinilla, Ch. Rab, R. Sulaiman, A. Zurlo

Abstract: We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analys… ▽ More We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analyse the disks in context of the stellar parameters and the environment of the Orion star-forming region. We use ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore use a combination of new and archival ALMA mm-continuum observations to characterize the dust masses present in the circumstellar disks. Within our sample we detect extended circumstellar disks in 10 of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori and V1012 Ori) and show scattered light asymmetries which may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high resolution imaging observations are also sensitive to close (sub)stellar companions and we detect 9 such objects in our sample of which 5 were previously unknown. We find in particular a possible sub-stellar companion (either a very low mass star or a high mass brown dwarf) 137 au from the star RY Ori. We find a strong anti-correlation between disk detection and multiplicity, with only 2 of our 10 disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimetre flux suggesting that disks that have a high dust content are typically bright in near-infrared scattered light. Conversely we do not find significant correlations between scattered light contrast of the disks and the stellar mass or age. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted by A&A

The SPHERE view of the Chamaeleon I star-forming region

Authors: C. Ginski, A. Garufi, M. Benisty, R. Tazaki, C. Dominik, A. Ribas, N. Engler, T. Birnstiel, G. Chauvin, G. Columba, S. Facchini, A. Goncharov, J. Hagelberg, T. Henning, M. Hogerheijde, R. G. van Holstein, J. Huang, T. Muto, P. Pinilla, K. Kanagawa, S. Kim, N. Kurtovic, M. Langlois, C. Manara, J. Milli , et al. (10 additional authors not shown)

Abstract: We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems… ▽ More We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems. In 13 of the 20 systems included in this study we detected resolved scattered light signals from circumstellar dust. For the CR Cha, CT Cha, CV Cha, SY Cha, SZ Cha, and VZ Cha systems we present the first detailed descriptions of the disks in scattered light. The observations found typically smooth or faint disks, often with little substructure, with the notable exceptions of SZ Cha, which shows an extended multiple-ringed disk, and WW Cha, which shows interaction with the cloud environment. New high S/N K- band observations of the HD 97048 system in our survey reveal a significant brightness asymmetry that may point to disk misalignment and subsequent shadowing of outer disk regions, possibly related to the suggested planet candidate in the disk. We resolve for the first time the stellar binary in the CS Cha system. Multiple wavelength observations of the disk around CS Cha have revealed that the system contains small, compact dust grains that may be strongly settled, consistent with numerical studies of circumbinary disks. We find in our sample that there is a strong anti-correlation between the presence of a (close) stellar companion and the detection of circumstellar material with five of our seven nondetections located in binary systems. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted by A&A

Mid-infrared evidence for iron-rich dust in the multi-ringed inner disk of HD 144432

Authors: J. Varga, L. B. F. M. Waters, M. Hogerheijde, R. van Boekel, A. Matter, B. Lopez, K. Perraut, L. Chen, D. Nadella, S. Wolf, C. Dominik, Á. Kóspál, P. Ábrahám, J. -C. Augereau, P. Boley, G. Bourdarot, A. Caratti o Garatti, F. Cruz-Sáenz de Miera, W. C. Danchi, V. Gámez Rosas, Th. Henning, K. -H. Hofmann, M. Houllé, J. W. Isbell, W. Jaffe , et al. (18 additional authors not shown)

Abstract: Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interfe… ▽ More Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interferometric data taken by the Very Large Telescope Interferometer (VLTI), combining PIONIER, GRAVITY, and MATISSE observations. Methods. We introduced a new physical disk model, TGMdust, to image the interferometric data, and to fit the disk structure and dust composition. We also performed equilibrium condensation calculations with GGchem. Results. Our best-fit model has three disk zones with ring-like structures at 0.15, 1.3, and 4.1 au. Assuming that the dark regions in the disk at ~0.9 au and at ~3 au are gaps opened by planets, we estimate the masses of the putative gap-opening planets to be around a Jupiter mass. We find evidence for an optically thin emission ($τ<0.4$) from the inner two disk zones ($r<4$ au) at $λ>3\ μ$m. Our silicate compositional fits confirm radial mineralogy gradients. To identify the dust component responsible for the infrared continuum emission, we explore two cases for the dust composition, one with a silicate+iron mixture and the other with a silicate+carbon one. We find that the iron-rich model provides a better fit to the spectral energy distribution. Conclusions. We propose that in the warm inner regions ($r<5$ au) of typical planet-forming disks, most if not all carbon is in the gas phase, while iron and iron sulfide grains are major constituents of the solid mixture along with forsterite and enstatite. Our analysis demonstrates the need for detailed studies of the dust in inner disks with new mid-infrared instruments such as MATISSE and JWST/MIRI. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: 29 pages, 24 figures

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

The bouncing barrier revisited: Impact on key planet formation processes and observational signatures

Authors: Carsten Dominik, Cornelis Dullemond

Abstract: Context. A leading paradigm in planet formation is currently the streaming instability and pebble accretion scenario. For this scenario, dust must grow into sizes in a specific regime of Stokes numbers in order to make these processes viable and sufficiently effective. The dust growth models currently in use do not implement some of the growth barriers suggested to be relevant in the literature. A… ▽ More Context. A leading paradigm in planet formation is currently the streaming instability and pebble accretion scenario. For this scenario, dust must grow into sizes in a specific regime of Stokes numbers in order to make these processes viable and sufficiently effective. The dust growth models currently in use do not implement some of the growth barriers suggested to be relevant in the literature. Aims. We want to investigate if the bouncing barrier, when effective, has impact on the time scales and efficiencies of processes like the streaming instability and pebble accretion, as well as on the observational appearance of planet-forming disks. Methods. We implement a formalism for the bouncing barrier into the publicly available dust growth model DustPy and run a series of models to understand the impact. Results. We find that the bouncing barrier has significant effect on the dust evolution in planet-forming disks. It reduces in many cases the size of the typical or largest particles available in the disk, it produces a very narrow, almost mono-disperse size distribution and removes most micrometer-sized grains in the process, with impact on scattered light images. It modifies the settling and therefore the effectiveness of and timescales for the streaming instability and for pebble accretion. An active bouncing barrier may well have observational consequences. It may reduce the strength if signatures of small particles (e.g. the 10 micron silicate feature), and it may create additional shadowed regions visible in scattered light images. Conclusions. Modeling of planet formation leaning heavily on the streaming instability and on pebble accretion should take the bouncing barrier into account. The complete removal of small grains in our model is not consistent with observations. However, this could be resolved by incomplete vertical mixing or some level of erosion in collisions. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 20 pages, 19 figures. Accepted for publication by Astronomy and Astrophysics

The polarisation properties of the HD 181327 debris ring. Evidence for sub-micron particles from scattered light observations

Authors: Julien Milli, Elodie Choquet, Ryo Tazaki, François Ménard, Jean-Charles Augereau, Johan Olofsson, Philippe Thébault, Olivier Poch, Anny-Chantal Levasseur-Regourd, Jérémie Lasue, Jean-Baptiste Renard, Edith Hadamcik, Clément Baruteau, Hans Martin Schmid, Natalia Engler, Rob G. van Holstein, Evgenij Zubko, Anne-Marie Lagrange, Sebastian Marino, Chirstophe Pinte, Carsten Dominik, Anthony Boccaletti, Maud Langlois, Alice Zurlo, Célia Desgrange , et al. (4 additional authors not shown)

Abstract: Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We… ▽ More Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We obtained high-contrast polarimetric images of HD 181327 in the H band with the SPHERE / IRDIS instrument on the Very Large Telescope (ESO). We complemented them with archival data from HST / NICMOS in the F110W filter reprocessed in the context of the Archival Legacy Investigations of Circumstellar Environments (ALICE) project. We developed a combined forward-modelling framework to simultaneously retrieve the scattering phase function in polarisation and intensity. We detected the debris disc around HD 181327 in polarised light and total intensity. We measured the scattering phase function and the degree of linear polarisation of the dust at 1.6 micron in the birth ring. The maximum polarisation is 23.6% +/- 2.6% and occurs between a scattering angle of 70 deg and 82 deg. We show that compact spherical particles made of a highly refractive and relatively absorbing material in a differential power-law size distribution of exponent $-3.5$ can simultaneously reproduce the polarimetric and total intensity scattering properties of the dust. This type of material cannot be obtained with a mixture of silicates, amorphous carbon, water ice, and porosity, and requires a more refracting component such as iron-bearing minerals. We reveal a striking analogy between the near-infrared polarisation of comets and that of HD 181327. The methodology developed here combining VLT/SPHERE and HST/NICMOS may be applicable in the future to combine the polarimetric capabilities of SPHERE with the sensitivity of JWST. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 21 pages, 22 figures, accepted in Astronomy & Astrophysics

The need for spatially resolved observations of PAHs in protoplanetary discs

Authors: K. Lange, C. Dominik, A. G. G. M. Tielens

Abstract: The signatures of polycyclic aromatic hydrocarbons (PAHs) have been observed in protoplanetary discs, and their emission features obtained from spectral energy distributions (SED) have been used in the literature to characterise their size and determine their abundance. Two simple disc models (uniform PAH distribution against a PAH gap in the inner disc) are compared to investigate the difference… ▽ More The signatures of polycyclic aromatic hydrocarbons (PAHs) have been observed in protoplanetary discs, and their emission features obtained from spectral energy distributions (SED) have been used in the literature to characterise their size and determine their abundance. Two simple disc models (uniform PAH distribution against a PAH gap in the inner disc) are compared to investigate the difference of their SED and obtainable information. We used the radiative transfer code RADMC-3D to model the SED of two protoplanetary discs orbiting a typical Herbig star, one of which features a depletion of PAHs in the inner disc. We further created artificial images of the discs at face-on view to extract radial profiles of the PAH emission in the infrared. We find that the extracted PAH features from an SED provide limited information about the PAHs in protoplanetary disc environments, except for the ionisation state. The distribution of PAHs in a protoplanetary disc influences the total observed PAH luminosity in a non-linear fashion and alters the relative strength between the 3.3\,$μ$m and 11.3\,$μ$m features. Furthermore, we produced radial profiles at the 3\,$μ$m, 6\,$μ$m and, 11\,$μ$m PAH emission features and find that they follow a double power-law profile where the slope reflects the radiative environment (single photon regime vs. multi-photon regime) in which the PAHs lie. Using spatially resolved techniques such as IFU or imaging in the era of the James Webb Space Telescope, we find that multi-wavelength radial emission profiles will not only provide information on the spatial distribution of the PAHs, but may also provide information on their size and underlying UV environment, which is crucial for photo-evaporative disc wind models. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 9 pages, 5 figures, accepted for Astronomy and Astrophysics

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): HD 34700 A unveils an inner ring

Authors: G. Columba, E. Rigliaco, R. Gratton, D. Mesa, V. D'Orazi, C. Ginski, N. Engler, J. P. Williams, J. Bae, M. Benisty, T. Birnstiel, P. Delorme, C. Dominik, S. Facchini, F. Menard, P. Pinilla, C. Rab, Á. Ribas, V. Squicciarini, R. G. van Holstein, A. Zurlo

Abstract: Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms. Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures. Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and… ▽ More Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms. Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures. Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and of ALMA observations, we analyse this system at multiple wavelengths. We study the rings and spiral arms morphology and the scattering properties of the dust. We discuss the possible causes of all the observed features. Results. We detect for the first time, in the H$α$ band, a ring extending from $\sim$65 au to ${\sim}$120 au, inside the ring already known from recent studies. These two have different physical and geometrical properties. Based on the scattering properties, the outer ring may consist of grains of typical size $a_{out} > 4 μm$, while the inner ring of smaller grains ($a_{in} <= 0.4 {μm}$). Two extended logarithmic spiral arms stem from opposite sides of the disk. The outer ring appears as a spiral arm itself, with a variable radial distance from the centre and extended substructures. ALMA data confirm the presence of a millimetric dust substructure centred just outside the outer ring, and detect misaligned gas rotation patterns for HD 34700 A and B. Conclusions. The complexity of HD 34700 A, revealed by the variety of observed features, suggests the existence of one or more disk-shaping physical mechanisms. Possible scenarios, compatible with our findings, involve the presence inside the disk of a yet undetected planet of several Jupiter masses and the system interaction with the surroundings by means of gas cloudlet capture or flybys. Further observations with JWST/MIRI or ALMA (gas kinematics) could shed more light on these. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: Accepted for publication on A&A. 14 + 5 pages, 9 + 7 figures (text + appendix)

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

The effect of Jupiter on the CAI storage problem

Authors: Stefan Jongejan, Carsten Dominik, Cornelis Dullemond

Abstract: By studying the distribution of calcium-aluminium-rich inclusions (CAIs) that are embedded within meteorites, we can learn about the dynamical history of the protoplanetary disk from which our Solar System formed. A long-standing problem concerning CAIs is the CAI storage problem. CAIs are thought to have formed at high temperatures near the Sun, but they are primarily found in carbonaceous chondr… ▽ More By studying the distribution of calcium-aluminium-rich inclusions (CAIs) that are embedded within meteorites, we can learn about the dynamical history of the protoplanetary disk from which our Solar System formed. A long-standing problem concerning CAIs is the CAI storage problem. CAIs are thought to have formed at high temperatures near the Sun, but they are primarily found in carbonaceous chondrites, which formed much further out, beyond the orbit of Jupiter. Additionally, radial drift of CAI particles should have removed them from the solar protoplanetary disk several million years before the parent bodies of meteorites in which they are encountered would have accreted. We revisit a previously suggested solution to the CAI storage problem by Desch, Kalyaan, and Alexander which proposed that CAIs were mixed radially outward through the disk and subsequently got trapped in a pressure maximum created by Jupiter's growing core opening a planet gap. Our aim is to investigate whether their solution still works when we take into account the infall phase during which the disk builds up from the collapse of a molecular cloud core. We build a 1D numerical code in Python using the DISKLAB package to simulate the evolution of the solar protoplanetary disk, starting with a collapsing molecular cloud. We find that outward transport of CAIs during the infall phase is very efficient, possibly mixing them all the way into the far outer disk. Subsequent inward radial drift collects CAIs in the pressure maximum beyond Jupiter's orbit while draining the inner disk, roughly reproducing parts of the result by Desch et al. By introducing CAI formation so early, abundances out to 100 AU remain significant, possibly not consistent with some meteoritic data. It is possible to create a disk that does not expand as far out and also does not push CAIs as far out by using a very slowly rotating cloud. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Journal ref: A&A 679, A45 (2023)

The influence of a static planetary atmosphere on spin transfer during pebble accretion

Authors: M. J. Yzer, R. G. Visser, C. Dominik

Abstract: We study the effect an atmosphere has on pebble orbits and spin build-up on a planet's surface during pebble accretion in the extreme case of a static atmosphere. We numerically integrate the equations of motion of pebbles in a planar, global frame with a planet, a central star and gas from a protoplanetary disc. An adiabatic atmosphere is then placed around the planet, and the spin deposited onto… ▽ More We study the effect an atmosphere has on pebble orbits and spin build-up on a planet's surface during pebble accretion in the extreme case of a static atmosphere. We numerically integrate the equations of motion of pebbles in a planar, global frame with a planet, a central star and gas from a protoplanetary disc. An adiabatic atmosphere is then placed around the planet, and the spin deposited onto the planet's surface is measured. These simulations are evaluated for different distances to the star, Stokes numbers, and planet masses. Pebble feedback to the gas is not taken into account. We find that a static atmosphere dampens the spin the planet's surface receives by absorbing part of the angular momentum of the pebbles and circularising their orbits. This could prevent the excessive spin values predicted in some 3D pebble accretion simulations without an atmosphere. For planets larger than 0.5 Earth masses, a stationary atmosphere absorbs all angular momentum, leaving no spin for the surface. Significant quantities of angular momentum are stored in the inner and intermediate atmosphere ($<0.3$ Bondi radii). Depending on the atmospheric and disc model, this spin could be transported either to the disc through atmospheric recycling or to the planet through drag between the surface and the atmosphere. Further research is required to quantify the spin transfer within the atmosphere. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 18 pages, 14 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 678, A37 (2023)

SHAMPOO: A stochastic model for tracking dust particles under the influence of non-local disk processes

Authors: M. Oosterloo, I. Kamp, W. van Westrenen, C. Dominik

Abstract: The abundances of CHNOS are crucial for the composition of planets. At the onset of planet formation, large amounts of these elements are stored in ices on dust grains in planet-forming disks. The evolution of this ice is affected by dynamical transport, collisional processes, and the formation and sublimation of ice. We aim to constrain the disk regions where these processes are fully coupled, an… ▽ More The abundances of CHNOS are crucial for the composition of planets. At the onset of planet formation, large amounts of these elements are stored in ices on dust grains in planet-forming disks. The evolution of this ice is affected by dynamical transport, collisional processes, and the formation and sublimation of ice. We aim to constrain the disk regions where these processes are fully coupled, and develop a flexible modelling approach that is able to predict the effects of these processes acting simultaneously on the CHNOS budgets of the dust in these regions. We compared timescales associated with these disk processes to constrain the disk regions where this approach is necessary, and developed the SHAMPOO code, which tracks the CHNOS abundances in the ice mantle of a single monomer dust particle, embedded in a larger aggregate and undergoing these processes simultaneously. The adsorption and photodesorption of monomer ices depend on the depth of the monomer in the aggregate. We investigated the effect of fragmentation velocity and aggregate filling factor on the amount of ice on monomers residing at r = 10 AU. The locations where disk processes are fully coupled depend on both grain size and ice species. Monomers embedded in aggregates with fragmentation velocities of 1 m/s are able to undergo adsorption and photodesorption more often compared to a fragmentation velocity of 5 m/s or 10 m/s. Aggregates with a filling factor of $10^{-3}$ are able to accumulate ice 22 times faster on average than aggregates with a filling factor of 1. As different grain sizes are coupled through collisions and the grain ice consists of multiple ice species, it is difficult to isolate the locations where disk processes are fully coupled, necessitating the development of the SHAMPOO code. The processing of ice may not be spatially limited to dust aggregate surfaces for either fragile or porous aggregates. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 30 pages, 24 figures, 4 tables, to appear in Astronomy & Astrophysics

Journal ref: A&A 674, A124 (2023)

Turbulent processing of PAHs in protoplanetary discs -- Coagulation and freeze-out leading to depletion of gas-phase PAH

Authors: K. Lange, C. Dominik, A. G. G. M. Tielens

Abstract: Polycyclic aromatic hydrocarbons (PAHs) have been detected in numerous circumstellar discs. We propose the continuous processing of PAHs through clustering, adsorption on dust grains, and their reverse-processes as key mechanisms to reduce the emission-capable PAH abundance in protoplanetary discs. This cycle of processing is driven by vertical turbulence in the disc mixing PAHs between the disc m… ▽ More Polycyclic aromatic hydrocarbons (PAHs) have been detected in numerous circumstellar discs. We propose the continuous processing of PAHs through clustering, adsorption on dust grains, and their reverse-processes as key mechanisms to reduce the emission-capable PAH abundance in protoplanetary discs. This cycle of processing is driven by vertical turbulence in the disc mixing PAHs between the disc midplane and the photosphere. We used a theoretical Monte Carlo model for photodesorption and a coagulation code in the disc midplane to estimate the relevance and timescale of these processes in a Herbig Ae/Be disc environment. By combining these components in a 1D vertical model, we calculated the gas-phase depletion of PAHs that stick as clusters on dust grains. Our results show that the clustering of gas-phase PAHs is very efficient, and that clusters with more than 100 monomers can grow for years before they are able to freeze out in the disc midplane. Once a PAH cluster is frozen on the dust grain surface, the large heat capacity of these clusters prevents them from evaporating off the grains in UV-rich environments such as the photosphere. Therefore, the clustering of PAHs followed by freeze-out can lead to a depletion of gas-phase PAHs in discs. Evaluated over the lifetime of protoplanetary discs, we find a depletion of PAHs by a factor that ranges between 50 and 1000 compared to the standard ISM abundance of PAHs in the inner disc through turbulent processing. Through these processes, we favour PAHs smaller than circumovalene as the major gas-phase emitters of the disc photosphere as larger PAH monomers cannot photodesorb from the grain surface. These gas-phase PAHs co-exist with large PAH clusters sticking on dust grains. We find a close relation between the amount of PAHs frozen out on dust grains and the dust population, as well as the strength of the vertical turbulence. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Accepted for Astronomy and Astrophysics. 18 pages, 14 figures

Fractal aggregates of sub-micron-sized grains in the young planet-forming disk around IM Lup

Authors: Ryo Tazaki, Christian Ginski, Carsten Dominik

Abstract: Despite rapidly growing disk observations, it remains a mystery what primordial dust aggregates look like and what the physical and chemical properties of their constituent grains (monomers) are in young planet-forming disks. Confrontation of models with observations to answer this mystery has been a notorious task because we have to abandon a commonly used assumption, perfectly spherical grains,… ▽ More Despite rapidly growing disk observations, it remains a mystery what primordial dust aggregates look like and what the physical and chemical properties of their constituent grains (monomers) are in young planet-forming disks. Confrontation of models with observations to answer this mystery has been a notorious task because we have to abandon a commonly used assumption, perfectly spherical grains, and take into account particles with complex morphology. In this Letter, we present the first thorough comparison between near-infrared scattered light of the young planet-forming disk around IM Lup and the light-scattering properties of complex-shaped dust particles. The availability of scattering observations at multiple wavelengths and over a significant range of scattering angles allows for the first determination of the monomer size, fractal dimension, and size of dust aggregates in a planet-forming disk. We show that the observations are best explained by fractal aggregates with a fractal dimension of 1.5 and a characteristic radius larger than $\sim2~μ$m. We also determined the radius of the monomer to be $\sim200$ nm, and monomers much smaller than this size can be ruled out on the premise that the fractal dimension is less than 2. Furthermore, dust composition comprising amorphous carbon is found to be favorable to simultaneously account for the faint scattered light and the flared disk morphology. Our results support that planet formation begins with fractal coagulation of sub-micron-sized grains. All the optical properties of complex dust particles computed in this study are publicly available. △ Less

Submitted 9 February, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 19 pages, 11 figures; Accepted for publication in ApJ Letters. The optical properties of dust particles are available at https://github.com/rtazaki1205/AggScatVIR . Updated the caption of fig. 8 to improve clarity

Observed polarized scattered light phase functions of planet-forming disks

Authors: Christian Ginski, Ryo Tazaki, Carsten Dominik, Tomas Stolker

Abstract: Dust particles are the building blocks from which planetary bodies are made. A major goal of the studies of planet-forming disks is to constrain the properties of dust particles and aggregates in order to trace their origin, structure, and the associated growth and mixing processes in the disk. Observations of scattering and/or emission of dust in a location of the disk often lead to degenerate in… ▽ More Dust particles are the building blocks from which planetary bodies are made. A major goal of the studies of planet-forming disks is to constrain the properties of dust particles and aggregates in order to trace their origin, structure, and the associated growth and mixing processes in the disk. Observations of scattering and/or emission of dust in a location of the disk often lead to degenerate information about the kind of particles, such as size, porosity, or fractal dimension of aggregates. Progress can be made by deriving the full (polarizing) scattering phase function of such particles at multiple wavelengths. This has now become possible by careful extraction from scattered light images. Such an extraction requires knowledge about the shape of the scattering surface in the disk and we discuss how to obtain such knowledge as well as the associated uncertainties. We use a sample of disk images from observations with VLT/SPHERE to, for the first time, extract the phase functions of a whole sample of disks with broad phase angle coverage. We find that polarized phase functions come in two categories. Comparing the extracted functions with theoretical predictions from rigorous T-Matrix computations of aggregates, we show that one category can be linked back to fractal, porous aggregates, while the other is consistent with more compact, less porous aggregates. We speculate that the more compact particles become visible in disks where embedded planets trigger enhanced vertical mixing. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ with parallel submission to arxiv for logistical reasons; We will update the arxiv submission with the final version once the manuscript was refereed and accepted by ApJ. Data files for all phase functions will be made available once refereeing is concluded. We welcome comments from the community

Forming equal mass planetary binaries via pebble accretion

Authors: T. J. Konijn, R. G. Visser, C. Dominik, C. W. Ormel

Abstract: Binary solar system objects are common and range from satellite systems with very large mass ratios $M_1/M_2$ to mass ratios very close to unity. A well-known example of a binary is the Pluto-Charon system. With Charon only eight times less massive than Pluto the question arises as for many other systems, why the mass-ratio is still close to unity. There is much evidence that (binary) planet(esima… ▽ More Binary solar system objects are common and range from satellite systems with very large mass ratios $M_1/M_2$ to mass ratios very close to unity. A well-known example of a binary is the Pluto-Charon system. With Charon only eight times less massive than Pluto the question arises as for many other systems, why the mass-ratio is still close to unity. There is much evidence that (binary) planet(esimal) formation happened early, when the protoplanetary gas disk was still around. It is likely that (some of) these binaries grew up together subject to pebble accretion. Here we focus on the question of how the mass arriving in the gravitational influence zone of the binary during pebble accretion, is distributed over the binary components. Does the accretion through time lead to a converging mass ratio, or to a diverging mass ratio? We numerically integrate pebble paths in the same well-known fashion as for a single mass subject to pebble accretion and track what the efficiency of accretion is for the two separate binary components, compared to a single body with the same mass. These numerical simulations are done for a range of binary mass-ratios, mutual separations, Stokes numbers and two orbital distances, 2.5 and 39 au. We find that in the limit where pebbles start to spiral around the primary (this holds for relatively large pebbles), the pebble preferentially collides with the secondary, causing the mass ratio to converge towards unity on Myr timescales. In this regime the total sweep-up efficiency can lower to half that of a pebble-accreting single body because pebbles that are thrown out of the system, after close encounters with the system. The results show that systems such as Pluto-Charon and other larger equal mass binaries could well have co-accreted by means of pebble accretion in the disk phase without producing binaries with highly diverging mass-ratios. △ Less

Submitted 17 February, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 10 pages, 8 figures, submitted to Astronomy and Astrophysics

Journal ref: A&A 670, A116 (2023)

Razor-thin dust layers in protoplanetary disks: Limits on the vertical shear instability

Authors: C. P. Dullemond, A. Ziampras, D. Ostertag, C. Dominik

Abstract: Context: Recent observations with the Atacama Large Millimeter Array (ALMA) have shown that the large dust aggregates observed at millimeter wavelengths settle to the midplane into a remarkably thin layer. Aims: We intend to find out if the geometric thinness of these layers is evidence against the vertical shear instability (VSI) operating in these disks. Methods: We performed hydrodynamic simula… ▽ More Context: Recent observations with the Atacama Large Millimeter Array (ALMA) have shown that the large dust aggregates observed at millimeter wavelengths settle to the midplane into a remarkably thin layer. Aims: We intend to find out if the geometric thinness of these layers is evidence against the vertical shear instability (VSI) operating in these disks. Methods: We performed hydrodynamic simulations of a protoplanetary disk with a locally isothermal equation of state, and let the VSI fully develop. We sprinkled dust particles and followed their motion as they got stirred up by the VSI. We determined for which grain size the layer becomes geometrically thin enough to be consistent with ALMA observations. We then verified if, with these grain sizes, it is still possible to generate a moderately optically thick layer at millimeter wavelengths, as observations appear to indicate. Results: We found that even very large dust aggregates with Stokes numbers close to unity get stirred up to relatively large heights above the midplane by the VSI, which is in conflict with the observed geometric thinness. For grains so large that the Stokes number exceeds unity, the layer can be made to remain thin, but we show that it is hard to make dust layers optically thick at ALMA wavelengths (e.g., tau(1.3mm)>=1) with such large dust aggregates. Conclusions: We conclude that protoplanetary disks with geometrically thin midplane dust layers cannot be VSI unstable, at least not down to the disk midplane. Explanations for the inhibition of the VSI include a reduced dust-to-gas ratio of the small dust grains that are responsible for the radiative cooling of the disk. A reduction of small grains by a factor of between 10 and 100 is sufficient to quench the VSI. Such a reduction is plausible in dust growth models, and still consistent with observations at optical and infrared wavelengths. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Scattered light detection of a possible disk wind in RY Tau

Authors: P. -G. Valegård, C. Ginski, C. Dominik, J. Bae, M. Benisty, T. Birnstiel, S. Facchini, A. Garufi, M. Hogerheijde, R. G. van Holstein, M. Langlois, C. F. Manara, P. Pinilla, Ch. Rab, Á. Ribas, L. B. F. M. Waters, J. Williams

Abstract: Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-i… ▽ More Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-infrared but not seen at optical wavelengths. In the near-infrared, dark wedges that separates the horns from the disk, indicating we may see the scattered light from a disk wind. We use archived ALMA and SPHERE/ZIMPOL I-band observations combined with newly acquired SPEHRE/IRDIS H-band observations and available literature to build a simple geometric model of the RY Tau disk and disk wind. We use Monte Carlo radiative transfer modelling \textit{MCMax3D} to create comparable synthetic observations that test the effect of a dusty wind on the optical effect in the observations. We constrain the grain size and dust mass needed in the disk wind to reproduce the effect from the observations. A model geometrically reminiscent of a dusty disk wind with small micron to sub-micron size grains elevated above the disk can reproduce the optical effect seen in the observations. The mass in the obscuring component of the wind has been constrained to $1\times10^{-9} M_{\odot} \leq M \leq 5\times10^{-8} M_{\odot}$ which corresponds to a lower limit mass loss rate in the wind of about $\sim 1\times10^{-8}M_{\odot}\mathrm{yr}^{-1}$. While an illuminate dust cavity cannot be ruled out without measurements of the gas velocity, we argue that a magnetically launched disk wind is the most likely scenario. △ Less

Submitted 5 October, 2022; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 11 pages, 5 figures

Journal ref: A&A 668, A25 (2022)

Locating dust and molecules in the inner circumstellar environment of R~Sculptoris with MATISSE

Authors: Julien Drevon, Florentin Millour, Pierre Cruzalèbes, Claudia Paladini, Josef Hron, A. Meilland, F. Allouche, K. -H. Hofmann, S. Lagarde, B. Lopez, A. Matter, R. Petrov, S. Robbe-Dubois, D. Schertl, M. Wittkowski, G. Zins, P. Ábrahám, P. Antonelli, U. Beckmann, P. Berio, F. Bettonvil, A. Glindemann, U. Graser, M. Heininger, Thomas Henning , et al. (27 additional authors not shown)

Abstract: 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 carbo… ▽ More 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. △ Less

Submitted 23 August, 2022; originally announced August 2022.

Comments: 19 pages, published in A&A

The disk of FU Orionis viewed with MATISSE/VLTI: first interferometric observations in $L$ and $M$ bands

Authors: F. Lykou, P. Ábrahám, L. Chen, J. Varga, Á. Kóspál, A. Matter, M. Siwak, Zs. M. Szabó, Z. Zhu, H. B. Liu, B. Lopez, F. Allouche, J. -C. Augereau, P. Berio, P. Cruzalèbes, C. Dominik, Th. Henning, K. -H. Hofmann, M. Hogerheijde, W. J. Jaffe, E. Kokoulina, S. Lagarde, A. Meilland, F. Millour, E. Pantin , et al. (8 additional authors not shown)

Abstract: The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data i… ▽ More The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data indicate that the dusty passive disk is silicate-rich. Only the innermost region of said dusty disk is found to emit strongly in the $N$ band, and it is resolved at an angular size of $\sim5$ mas, which translates to a diameter of about 2 au. The observations therefore place stringent constraints for the outer radius of the inner accretion disk. Dust radiative transfer simulations with RADMC-3D provide adequate fits to the spectral energy distribution from the optical to the submillimeter and to the interferometric observables when opting for an accretion rate $\dot{M}\sim 2\times 10^{-5}\, M_\odot$ yr$^{-1}$ and assuming $M_*=0.6\, M_\odot$. Most importantly, the hot inner accretion disk's outer radius can be fixed at 0.3 au. The outer radius of the dusty disk is placed at 100 au, based on constraints from scattered-light images in the literature. The dust mass contained in the disk is $2.4\times10^{-4}\, M_\odot$, and for a typical gas-to-dust ratio of 100, the total mass in the disk is approximately 0.02 $M_\odot$. We did not find any evidence for a nearby companion in the current interferometric data, and we tentatively explored the case of disk misalignment. For the latter, our modeling results suggest that the disk orientation is similar to that found in previous imaging studies by ALMA. Should there be an asymmetry in the very compact, inner accretion disk, this might be resolved at even smaller spatial scales ($\leq1$ mas). △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 17 pages, 14 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 663, A86 (2022)

The dusty heart of Circinus: I. Imaging the circumnuclear dust in N-band

Authors: Jacob W. Isbell, Klaus Meisenheimer, Jörg-Uwe Pott, Marko Stalevski, Konrad R. W. Tristram, Joel Sanchez-Bermudez, Karl-Heinz Hofmann, Violeta Gámez Rosas, Walter Jaffe, Leonard Burtscher, James Leftley, Romain Petrov, Bruno Lopez, Thomas Henning, Gerd Weigelt, Fatme Allouche, Philippe Berio, Felix Bettonvil, Pierre Cruzalebes, Carsten Dominik, Matthias Heininger, Michiel Hogerheijde, Stéphane Lagarde, Michael Lehmitz, Alexis Matter , et al. (6 additional authors not shown)

Abstract: 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 interpreta… ▽ More 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. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 663, A35 (2022)

How large are the monomers of dust aggregates in planet-forming disks?: Insights from quantitative optical and near-infrared polarimetry

Authors: Ryo Tazaki, Carsten Dominik

Abstract: Context: The size of the constituent particles (monomers) of dust aggregates is one of the most uncertain parameters directly affecting collisional growth of aggregates in planet-forming disks. Despite its importance, the monomer size has not yet been meaningfully constrained by disk observations. Aims: We attempt to derive the monomer size from optical and near-infrared (IR) polarimetric observat… ▽ More Context: The size of the constituent particles (monomers) of dust aggregates is one of the most uncertain parameters directly affecting collisional growth of aggregates in planet-forming disks. Despite its importance, the monomer size has not yet been meaningfully constrained by disk observations. Aims: We attempt to derive the monomer size from optical and near-infrared (IR) polarimetric observations of planet-forming disks. Methods: We perform a comprehensive parameter survey on the degree of linear polarization of light scattered by dust aggregates, using an exact numerical method called the $T$-matrix method. We investigate the effect of the monomer size, aggregate size, porosity, and composition on the degree of polarization. The obtained results are then compared with observed polarization fractions of several planet-forming disks at optical and near-IR wavelengths. Results: It is shown that the degree of polarization of aggregates depends sensitively on the monomer size unless the monomer size parameter is smaller than one or two. Comparing the simulation results with the disk observations, we find that the monomer radius is no greater than $0.4~μ$m. The inferred monomer size is therefore similar to subunit sizes of the solar system dust aggregates and the maximum size of interstellar grains. Conclusions: Optical and near-IR quantitative polarimetry will provide observational grounds on the initial conditions for dust coagulation and thereby planetesimal formation in planet-forming disks. △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: 17 pages, 12 figures, 1 table; Accepted for publication in A&A

Optical and Near-infrared View of Planet-forming Disks and Protoplanets

Authors: M. Benisty, C. Dominik, K. Follette, A. Garufi, C. Ginski, J. Hashimoto, M. Keppler, W. Kley, J. Monnier

Abstract: In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spi… ▽ More In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spirals, arcs, and shadows, are seen in all imaged infrared-bright disks to date, and are ubiquitous across spectral types. These substructures are believed to be the result of disk evolution processes, and in particular disk-planet interactions. Since small dust grains that scatter light are tightly bound to the disk's gas, these observations closely trace disk structures predicted by hydrodynamical models and serve as observational tests of the predictions of planet formation theories. We argue that the results of current and next-generation high-contrast imaging surveys will, when combined with complementary data from ALMA, lead to a much deeper understanding of the co-evolution of disks and planets, and the mechanisms by which planets form. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Accepted version, before community feedback. 31 pages (37 with references and appendix), 17 figures. This Chapter is dedicated to our co-author Willy Kley. An obituary written by Richard Nelson can be found at the end of the Chapter

Thermal imaging of dust hiding the black hole in the Active Galaxy NGC 1068

Authors: Violeta Gamez Rosas, Jacob W. Isbell, Walter Jaffe, Romain G. Petrov, James H. Leftley, Karl-Heinz Hofmann, Florentin Millour, Leonard Burtscher, Klaus Meisenheimer, Anthony Meilland, Laurens B. F. M. Waters, Bruno Lopez, Stephane Lagarde, Gerd Weigelt, Philippe Berio, Fatme Allouche, Sylvie Robbe-Dubois, Pierre Cruzalebes, Felix Bettonvil, Thomas Henning, Jean-Charles Augereau, Pierre Antonelli, Udo Beckmann, Roy van Boekel, Philippe Bendjoya , et al. (27 additional authors not shown)

Abstract: 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 evolutio… ▽ More 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. △ Less

Submitted 27 December, 2021; originally announced December 2021.

Comments: In press at Nature

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Authors: A. Chiavassa, K. Kravchenko, M. Montargès, F. Millour, A. Matter, B. Freytag, M. Wittkowski, V. Hocdé, P. Cruzalèbes, F. Allouche, B. Lopez, S. Lagarde, R. G. Petrov, A. Meilland, S. Robbe-Dubois, K. -H. Hofmann, G. Weigelt, P. Berio, P. Bendjoya, F. Bettonvil, A. Domiciano de Souza, M. Heininger, Th. Henning, J. W. Isbell, W. Jaffe , et al. (28 additional authors not shown)

Abstract: Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reco… ▽ More Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: Accepted for publication on Astronomy and Astrophysics

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

VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris

Authors: K. -H. Hofmann, A. Bensberg, D. Schertl, G. Weigelt, S. Wolf, A. Meilland, F. Millour, L. B. F. M. Waters, S. Kraus, K. Ohnaka, B. Lopez, R. G. Petrov, S. Lagarde, Ph. Berio, F. Allouche, S. Robbe-Dubois, W. Jaffe, Th. Henning, C. Paladini, M. Schöller, A. Mérand, A. Glindemann, U. Beckmann, M. Heininger, F. Bettonvil , et al. (36 additional authors not shown)

Abstract: 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 di… ▽ More 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. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: accepted for publication in A&A

An extended scattered light disk around AT Pyx -- Possible planet formation in a cometary globule

Authors: C. Ginski, R. Gratton, A. Bohn, C. Dominik, S. Jorquera, G. Chauvin, J. Milli, M. Rodriguez, M. Benisty, R. Launhardt, A. Mueller, G. Cugno, R. G. van Holstein, A. Boccaletti, G. A. Muro-Arena, S. Desidera, M. Keppler, A. Zurlo, E. Sissa, T. Henning, M. Janson, M. Langlois, M. Bonnefoy, F. Cantalloube, V. D'Orazi , et al. (13 additional authors not shown)

Abstract: To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have obser… ▽ More To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. We used the extreme adaptive optics imager VLT/SPHERE to observe AT Pyx in polarized light as well as total intensity in the J, H and K-band. Additionally we employed VLT/NACO to observe the system in the L-band. We resolve the disk around AT Pyx in scattered light across multiple wavelengths. We find an extended (>126 au) disk, with an intermediate inclination between 35 deg and 42 deg. The disk shows complex sub-structure and we identify 2 and possibly 3 spiral-like features. Depending on the precise geometry of the disk (which we can not unambiguously infer from our data) the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet equal in mass to Jupiter. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are not sensitive enough to detect gas giants. AT Pyx is the first disk in a cometray globule in the Gum Nebula which is spatially resolved. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with higher sensitivity. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 11 pages, 9 figures, accepted for publication in A&A

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

A SPHERE survey of self-shadowed planet-forming disks

Authors: A. Garufi, C. Dominik, C. Ginski, M. Benisty, R. G. van Holstein, Th. Henning, N. Pawellek, C. Pinte, H. Avenhaus, S. Facchini, R. Galicher, R. Gratton, F. Menard, G. Muro-Arena, J. Milli, T. Stolker, A. Vigan, M. Villenave, T. Moulin, A. Origne, F. Rigal, J. -F. Sauvage, L. Weber

Abstract: To date, nearly two hundred planet-forming disks have been imaged with high resolution. Our propensity to study bright and extended objects is however biasing our view of the disk demography. In this work, we contribute to alleviate this bias by analyzing fifteen disks targeted with VLT/SPHERE that look faint in scattered light. Sources were selected based on a low far-IR excess from the spectral… ▽ More To date, nearly two hundred planet-forming disks have been imaged with high resolution. Our propensity to study bright and extended objects is however biasing our view of the disk demography. In this work, we contribute to alleviate this bias by analyzing fifteen disks targeted with VLT/SPHERE that look faint in scattered light. Sources were selected based on a low far-IR excess from the spectral energy distribution. The comparison with the ALMA images available for a few sources shows that the scattered light surveyed by these datasets is only detected from a small portion of the disk extent. The mild anti-correlation between the disk brightness and the near-IR excess demonstrates that these disks are self-shadowed: the inner disk rim intercepts much starlight and leaves the outer disk in penumbra. Based on the uniform distribution of the disk brightness in scattered light across all spectral types, self-shadowing would act similarly for inner rims at a different distance from the star. We discuss how the illumination pattern of the outer disk may evolve with time. Some objects in the sample are proposed to be at an intermediate stage toward bright disks from the literature with either no shadow or with sign of azimuthally confined shadows. △ Less

Submitted 29 November, 2021; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A\&A

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

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

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

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

Submitted 30 October, 2021; originally announced November 2021.

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

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

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

Authors: B. Lopez, S. Lagarde, R. G. Petrov, W. Jaffe, P. Antonelli, F. Allouche, P. Berio, A. Matter, A. Meilland, F. Millour, S. Robbe-Dubois, Th. Henning, G. Weigelt, A. Glindemann, T. Agocs, Ch. Bailet, U. Beckmann, F. Bettonvil, R. van Boekel, P. Bourget, Y. Bresson, P. Bristow, P. Cruzalèbes, E. Eldswijk, Y. Fanteï Caujolle , et al. (128 additional authors not shown)

Abstract: Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have develope… ▽ More Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now. △ Less

Submitted 2 March, 2022; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: 24 pages, 26 figures, submitted to Astronomy & Astrophysics

Journal ref: A&A 659, A192 (2022)

Stability of Polycyclic Aromatic Hydrocarbon Clusters in Protoplanetary Disks

Authors: K. Lange, C. Dominik, A. G. G. M. Tielens

Abstract: The infrared signature of polycyclic aromatic hydrocarbons (PAHs) are present in many protostellar disks and these speciesare thought to play an important role in heating of the gas in the photosphere. We aim to consider PAH cluster formation as one possible cause for non-detections of PAH features in protoplanetary disks. We test the necessary conditions for cluster formation and cluster dissocia… ▽ More The infrared signature of polycyclic aromatic hydrocarbons (PAHs) are present in many protostellar disks and these speciesare thought to play an important role in heating of the gas in the photosphere. We aim to consider PAH cluster formation as one possible cause for non-detections of PAH features in protoplanetary disks. We test the necessary conditions for cluster formation and cluster dissociation by stellar optical and FUV photons in protoplanetarydisks using a Herbig Ae/Be and a T Tauri star disk model. We perform Monte-Carlo (MC) and statistical calculations to determine dissociation rates for coronene, circumcoronene and circumcoronene clusters with sizes between 2 and 200 cluster members. By applying general disk models to our Herbig Ae/Be and T Tauri star model, we estimate the formation rate of PAH dimers and compare these with the dissociation rates. We show that the formation of PAH dimers can take place in the inner 100 AU of protoplanetary disks in sub-photospheric layers. Dimer formation takes seconds to years allowing them to grow beyond dimer size in a short time. We further demonstrate that PAH cluster increase their stability while they grow if they are located beyond a critical distance that depends on stellar properties and PAH species. The comparison with the local vertical mixing time scale allows a determination of the minimum cluster size necessaryfor survival of PAH clusters. Considering the PAH cluster formation sites, cluster survival in the photosphere of the inner disk of Herbig stars isunlikely because of the high UV radiation. For the T Tauri stars, survival of coronene, circumcoronene and circumcircumcoronene clusters is possible and cluster formation should be considered as one possible explanation for low PAH detection rates in T Tauri star disks. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: Accepted for Astronomy and Astrophysics. 12 pages, 8 figures

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

A Circumplanetary Disk Around PDS70c

Authors: Myriam Benisty, Jaehan Bae, Stefano Facchini, Miriam Keppler, Richard Teague, Andrea Isella, Nicolas T. Kurtovic, Laura M. Perez, Anibal Sierra, Sean M. Andrews, John Carpenter, Ian Czekala, Carsten Dominik, Thomas Henning, Francois Menard, Paola Pinilla, Alice Zurlo

Abstract: PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$α$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter… ▽ More PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$α$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 855$μ$m at high angular resolution ($\sim$20mas, 2.3au) that aim to resolve the circumplanetary disks and constrain their dust masses. Our observations confirm the presence of a compact source of emission co-located with PDS70c, spatially separated from the circumstellar disk and less extended than $\sim$1.2au in radius, a value close to the expected truncation radius of the cicumplanetary disk at a third of the Hill radius. The emission around PDS70c has a peak intensity of $\sim$86$\pm$16 $μ\mathrm{Jy}~\mathrm{beam}^{-1}$ which corresponds to a dust mass of $\sim$0.031M$_{\oplus}$ or $\sim$0.007M$_{\oplus}$, assuming that it is only constituted of 1 $μ$m or 1 mm sized grains, respectively. We also detect extended, low surface brightness continuum emission within the cavity near PDS70b. We observe an optically thin inner disk within 18au of the star with an emission that could result from small micron-sized grains transported from the outer disk through the orbits of b and c. In addition, we find that the outer disk resolves into a narrow and bright ring with a faint inner shoulder. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: ApJ Letters, in press; 19 pages, 9 figures

Journal ref: 2021, ApJ, 916L, 2B

The formation of planetary systems with SPICA

Authors: I. Kamp, M. Honda, H. Nomura, M. Audard, D. Fedele, L. B. F. M. Waters, Y. Aikawa, A. Banzatti, J. E. Bowey, M. Bradford, C. Dominik, K. Furuya, E. Habart, D. Ishihara, D. Johnstone, G. Kennedy, M. Kim, Q. Kral, S. P. Lai, B. Larsson, M. McClure, A. Miotello, M. Momose, T. Nakagawa, D. Naylor , et al. (16 additional authors not shown)

Abstract: In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared spac… ▽ More In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated) △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: accepted for publication in PASA

First MATISSE L-band observations of HD 179218. Is the inner 10 au region rich in carbon dust particles?

Authors: E. Kokoulina, A. Matter, B. Lopez, E. Pantin, N. Ysard, G. Weigelt, E. Habart, J. Varga, A. Jones, A. Meilland, E. Dartois, L. Klarmann, J. -C. Augereau, R. van Boekel, M. Hogerheijde, G. Yoffe, L. B. F. M. Waters, C. Dominik, W. Jaffe, F. Millour, Th. Henning, K. -H. Hofmann, D. Schertl, S. Lagarde, R. G. Petrov , et al. (36 additional authors not shown)

Abstract: 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 ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as… ▽ More 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 ($\sim $0.1 to 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 (PAHs)), 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 ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. 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. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) 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. △ Less

Submitted 29 July, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Journal ref: A&A 652, A61 (2021)

What happened before? -- The disks around the precursors of young Herbig Ae/Be stars

Authors: P. G. Valegård, L. B. F. M. Waters, C. Dominik

Abstract: We seek to find the precursors of the Herbig Ae/Be stars in the solar vicinity within 500 pc from the Sun. We do this by creating an optically selected sample of intermediate mass T-Tauri stars (IMTT stars) here defined as stars of masses $1.5 M_{\odot}\leq M_* \leq 5 M_{\odot}$ and spectral type between F and K3, from literature. We use literature optical photometry (0.4-1.25$μ$m) and distances d… ▽ More We seek to find the precursors of the Herbig Ae/Be stars in the solar vicinity within 500 pc from the Sun. We do this by creating an optically selected sample of intermediate mass T-Tauri stars (IMTT stars) here defined as stars of masses $1.5 M_{\odot}\leq M_* \leq 5 M_{\odot}$ and spectral type between F and K3, from literature. We use literature optical photometry (0.4-1.25$μ$m) and distances determined from Gaia DR2 parallax measurements together with Kurucz stellar model spectra to place the stars in a HR-diagram. With Siess evolutionary tracks we identify intermediate mass T-Tauri stars from literature and derive masses and ages. We use Spitzer spectra to classify the disks around the stars into Meeus Group I and Group II disks based on their [F$_{30}$/F$_{13.5}$] spectral index. We also examine the 10$μ$m silicate dust grain emission and identify emission from Polycyclic Aromatic Hydrocarbons (PAH). From this we build a qualitative picture of the disks around the intermediate mass T-Tauri stars and compare this with available spatially resolved images at infrared and at sub-millimeter wavelengths to confirm our classification. We find 49 intermediate mass T-Tauri stars with infrared excess. The identified disks are similar to the older Herbig Ae/Be stars in disk geometries and silicate dust grain population. Spatially resolved images at infra-red and sub-mm wavelengths suggest gaps and spirals are also present around the younger precursors to the Herbig Ae/Be stars. Comparing the timescale of stellar evolution towards the main sequence and current models of protoplanetary disk evolution the similarity between Herbig Ae/Be stars and the intermediate mass T-Tauri stars points towards an evolution of Group I and Group II disks that are disconnected, and that they represent two different evolutionary paths. △ Less

Submitted 30 April, 2021; v1 submitted 29 April, 2021; originally announced April 2021.

Comments: 25 pages, 9 figures, 8 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 652, A133 (2021)

Mid-infrared circumstellar emission of the long-period Cepheid l Carinae resolved with VLTI/MATISSE

Authors: V. Hocdé, N. Nardetto, A. Matter, E. Lagadec, A. Mérand, P. Cruzalèbes, A. Meilland, F. Millour, B. Lopez, P. Berio, G. Weigelt, R. Petrov, J. W. Isbell, W. Jaffe, P. Kervella, A. Glindemann, M. Schöller, F. Allouche, A. Gallenne, A. Domiciano de Souza, G. Niccolini, E. Kokoulina, J. Varga, S. Lagarde, J. -C. Augereau , et al. (129 additional authors not shown)

Abstract: 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. W… ▽ More 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. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

New binaries from the SHINE survey

Authors: M. Bonavita, R. Gratton, S. Desidera, V. Squicciarini, V. D'Orazi, A. Zurlo, B. Biller, G. Chauvin, C. Fontanive, M. Janson, S. Messina, F. Menard, M. Meyer, A. Vigan, H. Avenhaus, R. Asensio Torres, J. -L. Beuzit, A. Boccaletti, M. Bonnefoy, W. Brandner, F. Cantalloube, A. Cheetham, M. Cudel, S. Daemgen, P. Delorme , et al. (45 additional authors not shown)

Abstract: We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets obser… ▽ More We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose. △ Less

Submitted 28 July, 2022; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 59 pages, 11 tables, 18 figures

Journal ref: A&A 663, A144 (2022)

HD142527: Quantitative disk polarimetry with SPHERE

Authors: S. Hunziker, H. M. Schmid, J. Ma, F. Menard, H. Avenhaus, A. Boccaletti, J. L. Beuzit, G. Chauvin, K. Dohlen, C. Dominik, N. Engler, C. Ginski, R. Gratton, T. Henning, M. Langlois, J. Milli, D. Mouillet, C. Tschudi, R. G. van Holstein, A. Vigan

Abstract: We present high-precision photometry and polarimetry for the protoplanetary disk around HD142527, with a focus on determining the light scattering parameters of the dust. We re-reduced polarimetric differential imaging data of HD142527 in the VBB (735 nm) and H-band (1625 nm) from the ZIMPOL and IRDIS subinstruments of SPHERE/VLT. With polarimetry and photometry based on reference star differentia… ▽ More We present high-precision photometry and polarimetry for the protoplanetary disk around HD142527, with a focus on determining the light scattering parameters of the dust. We re-reduced polarimetric differential imaging data of HD142527 in the VBB (735 nm) and H-band (1625 nm) from the ZIMPOL and IRDIS subinstruments of SPHERE/VLT. With polarimetry and photometry based on reference star differential imaging, we were able to measure the linearly polarized intensity and the total intensity of the light scattered by the circumstellar disk with high precision. We used simple Monte Carlo simulations of multiple light scattering by the disk surface to derive constraints for three scattering parameters of the dust: the maximum polarization of $P_{\rm max}$, the asymmetry parameter $g$, and the single-scattering albedo $ω$. We measure a reflected total intensity of $51.4\pm1.5$ mJy and $206\pm12$ mJy and a polarized intensity of $11.3\pm0.3$ mJy and $55.1\pm3.3$ mJy in the VBB and H-band, respectively. We also find in the visual range a degree of polarization that varies between $28\%$ on the far side of the disk and $17\%$ on the near side. The disk shows a red color for the scattered light intensity and the polarized intensity, which are about twice as high in the near-infrared when compared to the visual. We determine with model calculations the scattering properties of the dust particles and find evidence for strong forward scattering ($g\approx 0.5-0.75$), relatively low single-scattering albedo ($ω\approx 0.2-0.5$), and high maximum polarization ($P_{\rm max} \approx 0.5-0.75$) at the surface on the far side of the disk for both observed wavelengths. The optical parameters indicate the presence of large aggregate dust particles, which are necessary to explain the high maximum polarization, the strong forward-scattering nature of the dust, and the observed red disk color. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: 20 pages, 14 figures

Journal ref: A&A 648, A110 (2021)

The SPHERE infrared survey for exoplanets (SHINE)- I Sample definition and target characterization

Authors: S. Desidera, G. Chauvin, M. Bonavita, S. Messina, H. LeCoroller, T. Schmidt, R. Gratton, C. Lazzoni, M. Meyer, J. Schlieder, A. Cheetham, J. Hagelberg, M. Bonnefoy, M. Feldt, A-M. Lagrange, M. Langlois, A. Vigan, T. G. Tan, F. -J. Hambsch, M. Millward, J. Alcala, S. Benatti, W. Brandner, J. Carson, E. Covino , et al. (83 additional authors not shown)

Abstract: Large surveys with new-generation high-contrast imaging instruments are needed to derive the frequency and properties of exoplanet populations with separations from $\sim$5 to 300 AU. A careful assessment of the stellar properties is crucial for a proper understanding of when, where, and how frequently planets form, and how they evolve. The sensitivity of detection limits to stellar age makes this… ▽ More Large surveys with new-generation high-contrast imaging instruments are needed to derive the frequency and properties of exoplanet populations with separations from $\sim$5 to 300 AU. A careful assessment of the stellar properties is crucial for a proper understanding of when, where, and how frequently planets form, and how they evolve. The sensitivity of detection limits to stellar age makes this a key parameter for direct imaging surveys. We describe the SpHere INfrared survey for Exoplanets (SHINE), the largest direct imaging planet-search campaign initiated at the VLT in 2015 in the context of the SPHERE Guaranteed Time Observations of the SPHERE consortium. In this first paper we present the selection and the properties of the complete sample of stars surveyed with SHINE, focusing on the targets observed during the first phase of the survey (from February 2015 to February 2017). This early sample composed of 150 stars is used to perform a preliminary statistical analysis of the SHINE data, deferred to two companion papers presenting the survey performance, main discoveries, and the preliminary statistical constraints set by SHINE. Based on a large database collecting the stellar properties of all young nearby stars in the solar vicinity (including kinematics, membership to moving groups, isochrones, lithium abundance, rotation, and activity), we selected the original sample of 800 stars that were ranked in order of priority according to their sensitivity for planet detection in direct imaging with SPHERE. The properties of the stars that are part of the early statistical sample were revisited, including for instance measurements from the GAIA Data Release 2. △ Less

Submitted 7 March, 2021; originally announced March 2021.

Journal ref: A&A 651, A70 (2021)

The SPHERE infrared survey for exoplanets (SHINE) -- II. Observations, Data reduction and analysis Detection performances and early-results

Authors: M. Langlois, R. Gratton, A. -M. Lagrange, P. Delorme, A. Boccaletti, M. Bonnefoy, A. -L. Maire, D. Mesa, G. Chauvin, S. Desidera, A. Vigan, A. Cheetham, J. Hagelberg, M. Feldt, M. Meyer, P. Rubini, H. Le Coroller, F. Cantalloube, B. Biller, M. Bonavita, T. Bhowmik, W. Brandner, S. Daemgen, V. D'Orazi, O. Flasseur , et al. (96 additional authors not shown)

Abstract: Over the past decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) from their host stars. To understand their formation and evolution mechanisms, we have initiated in 2015 the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars to explore their demographics.} {We aim to… ▽ More Over the past decades, direct imaging has confirmed the existence of substellar companions (exoplanets or brown dwarfs) on wide orbits (>10 au) from their host stars. To understand their formation and evolution mechanisms, we have initiated in 2015 the SPHERE infrared survey for exoplanets (SHINE), a systematic direct imaging survey of young, nearby stars to explore their demographics.} {We aim to detect and characterize the population of giant planets and brown dwarfs beyond the snow line around young, nearby stars. Combined with the survey completeness, our observations offer the opportunity to constrain the statistical properties (occurrence, mass and orbital distributions, dependency on the stellar mass) of these young giant planets.} {In this study, we present the observing and data analysis strategy, the ranking process of the detected candidates, and the survey performances for a subsample of 150 stars, which are representative of the full SHINE sample. The observations were conducted in an homogeneous way from February 2015 to February 2017 with the dedicated ground-based VLT/SPHERE instrument equipped with the IFS integral field spectrograph and the IRDIS dual-band imager covering a spectral range between 0.9 and 2.3 $μ$m. We used coronographic, angular and spectral differential imaging techniques to reach the best detection performances for this study down to the planetary mass regime.} △ Less

Submitted 5 March, 2021; originally announced March 2021.

Journal ref: A&A 651, A71 (2021)

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Late infall causing disk misalignment and dynamic structures in SU Aur

Authors: C. Ginski, S. Facchini, J. Huang, M. Benisty, D. Vaendel, L. Stapper, C. Dominik, J. Bae, F. Menard, G. Muro-Arena, M. Hogerheijde, M. McClure, R. G. van Holstein, T. Birnstiel, Y. Boehler, A. Bohn, M. Flock, E. E. Mamajek, C. F. Manara, P. Pinilla, C. Pinte, A. Ribas

Abstract: Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resol… ▽ More Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 18 pages, 12 figures, published in ApJL on 18-02-2021

Journal ref: ApJL, 908, L25 (2021)

The radial structure of planetary bodies formed by the streaming instability

Authors: Rico G. Visser, Joanna Drążkowska, Carsten Dominik

Abstract: Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global an… ▽ More Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global and lack the resolution to follow the final collapse stage of a pebble cloud within its Hill radius. We aim to track the collapse of a gravitationally bound pebble cloud subject to mutual collisions and gas drag with the representative particle approach. We determine the radial pebble size distribution of the collapsed core and the impact of mutual pebble collisions on the pebble size distribution. We find that virial equilibrium is never reached during the cloud evolution and that, in general, pebbles with given Stokes number (St) collapse towards an optically thick core in a sequence from aerodynamically largest to aerodynamically smallest. We show that at the location for which the core becomes optically thick, the terminal velocity is well below the fragmentation threshold velocity. While collisional processing is negligible during cloud evolution, the collisions that do occur are sticking. These results support the observations that comets and small planetary bodies are composed of primordial pebbles in the milimeter to centimeter size range △ Less

Submitted 4 February, 2021; v1 submitted 22 January, 2021; originally announced January 2021.

Comments: 15 pages, accepted in A&A on January 21, 2021

Journal ref: A&A 647, A126 (2021)

A survey of the linear polarization of directly imaged exoplanets and brown dwarf companions with SPHERE-IRDIS. First polarimetric detections revealing disks around DH Tau B and GSC 6214-210 B

Authors: R. G. van Holstein, T. Stolker, R. Jensen-Clem, C. Ginski, J. Milli, J. de Boer, J. H. Girard, Z. Wahhaj, A. J. Bohn, M. A. Millar-Blanchaer, M. Benisty, M. Bonnefoy, G. Chauvin, C. Dominik, S. Hinkley, C. U. Keller, M. Keppler, M. Langlois, S. Marino, F. Ménard, C. Perrot, T. O. B. Schmidt, A. Vigan, A. Zurlo, F. Snik

Abstract: Young giant planets and brown dwarf companions emit near-infrared radiation that can be linearly polarized up to several percent. This polarization can reveal the presence of a circumsubstellar accretion disk, rotation-induced oblateness of the atmosphere, or an inhomogeneous distribution of atmospheric dust clouds. We measured the near-infrared linear polarization of 20 known directly imaged exop… ▽ More Young giant planets and brown dwarf companions emit near-infrared radiation that can be linearly polarized up to several percent. This polarization can reveal the presence of a circumsubstellar accretion disk, rotation-induced oblateness of the atmosphere, or an inhomogeneous distribution of atmospheric dust clouds. We measured the near-infrared linear polarization of 20 known directly imaged exoplanets and brown dwarf companions with the high-contrast imager SPHERE-IRDIS at the VLT. We reduced the data using the IRDAP pipeline to correct for the instrumental polarization and crosstalk with an absolute polarimetric accuracy <0.1% in the degree of polarization. We report the first detection of polarization originating from substellar companions, with a polarization of several tenths of a percent for DH Tau B and GSC 6214-210 B in H-band. By comparing the measured polarization with that of nearby stars, we find that the polarization is unlikely to be caused by interstellar dust. Because the companions have previously measured hydrogen emission lines and red colors, the polarization most likely originates from circumsubstellar disks. Through radiative transfer modeling, we constrain the position angles of the disks and find that the disks must have high inclinations. The presence of these disks as well as the misalignment of the disk of DH Tau B with the disk around its primary star suggest in situ formation of the companions. For the 18 other companions, we do not detect significant polarization and place subpercent upper limits on their degree of polarization. These non-detections may indicate the absence of circumsubstellar disks, a slow rotation rate of young companions, the upper atmospheres containing primarily submicron-sized dust grains, and/or limited cloud inhomogeneity. Finally, we present images of the circumstellar disks of DH Tau, GQ Lup, PDS 70, Beta Pic, and HD 106906. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in A&A. Shortened abstract. 29 pages, 22 figures

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

The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?

Authors: J. Varga, M. Hogerheijde, R. van Boekel, L. Klarmann, R. Petrov, L. B. F. M. Waters, S. Lagarde, E. Pantin, Ph. Berio, G. Weigelt, S. Robbe-Dubois, B. Lopez, F. Millour, J. -C. Augereau, H. Meheut, A. Meilland, Th. Henning, W. Jaffe, F. Bettonvil, P. Bristow, K. -H. Hofmann, A. Matter, G. Zins, S. Wolf, F. Allouche , et al. (111 additional authors not shown)

Abstract: 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. W… ▽ More 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<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $μ$m-sized grains. Refractory grains or large ($\gtrsim 10\ μ$m-sized) grains could be the origin for this emission. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: accepted for publication in A&A

Journal ref: A&A 647, A56 (2021)

A triple star in disarray -- Multi-epoch observations of T Tauri with VLT-SPHERE and LBT-LUCI

Authors: M. Kasper, K. K. R. Santhakumari, T. M. Herbst, R. van Boekel, F. Menard, R. Gratton, R. G. van Holstein, M. Langlois, C. Ginski, A. Boccaletti, J. de Boer, P. Delorme, S. Desidera, C. Dominik, J. Hagelberg, T. Henning, R. Koehler, D. Mesa, S. Messina, A. Pavlov, C. Petit, E. Rickman, A. Roux, F. Rigal, A. Vigan , et al. (2 additional authors not shown)

Abstract: T Tauri remains an enigmatic triple star for which neither the evolutionary state of the stars themselves, nor the geometry of the complex outflow system is completely understood. Eight-meter class telescopes equipped with state-of-the-art adaptive optics provide the spatial resolution necessary to trace tangential motion of features over a timescale of a few years, and they help to associate them… ▽ More T Tauri remains an enigmatic triple star for which neither the evolutionary state of the stars themselves, nor the geometry of the complex outflow system is completely understood. Eight-meter class telescopes equipped with state-of-the-art adaptive optics provide the spatial resolution necessary to trace tangential motion of features over a timescale of a few years, and they help to associate them with the different outflows. We used J-, H-, and K-band high-contrast coronagraphic imaging with VLT-SPHERE recorded between 2016 and 2018 to map reflection nebulosities and obtain high precision near-infrared (NIR) photometry of the triple star. We also present molecular hydrogen emission maps of the 1-0 S(1) line at 2.122 micron obtained with LBT-LUCI during its commissioning period at the end of 2016. The data reveal a number of new features in the system, some of which are seen in reflected light and some are seen in H2 emission; furthermore, they can all be associated with the main outflows. The tangential motion of the features provides compelling evidence that T Tauri Sb drives the southeast-northwest outflow. T Tauri Sb has recently faded probably because of increased extinction as it passes through the southern circumbinary disk. While T Tauri Sb is approaching periastron, T Tauri Sa instead has brightened and is detected in all our J-band imagery for the first time. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 10 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 644, A114 (2020)

A single-armed spiral in the protoplanetary disk around HD34282 ?

Authors: J. de Boer, C. Ginski, G. Chauvin, F. Menard, M. Benisty, C. Dominik, K. Maaskant, J. H. Girard, G. van der Plas, A. Garufi, C. Perrot, T. Stolker, H. Avenhaus, A. Bohn, A. Delboulbe, M. Jaquet, T. Buey, O. Moller-Nilsson, J. Pragt, T. Fusco

Abstract: During the evolution of protoplanetary disks into planetary systems we expect to detect signatures that trace mechanisms such as planet-disk interaction. Protoplanetary disks display a large variety of structures in recently published high-spatial resolution images. However, the three-dimensional morphology of these disks is often difficult to infer from the two-dimensional projected images we obs… ▽ More During the evolution of protoplanetary disks into planetary systems we expect to detect signatures that trace mechanisms such as planet-disk interaction. Protoplanetary disks display a large variety of structures in recently published high-spatial resolution images. However, the three-dimensional morphology of these disks is often difficult to infer from the two-dimensional projected images we observe. We spatially resolve the disk around HD 34282 using VLT/SPHERE in polarimetric imaging mode. We retrieve a profile for the height of the scattering surface to create a height-corrected deprojection, which simulates a face-on orientation. The disk displays a complex scattering surface. An inner clearing or cavity extending up to r<0.28" (88 au) is surrounded by a bright inclined (i = 56 deg) ring with a position angle of 119 deg. The center of this ring is offset from the star along the minor axis with 0.07", which can be explained with a disk-height of 26 au above the mid-plane. Outside this ring, beyond its south-eastern ansa we detect an azimuthal asymmetry or blob at r ~ 0.4". At larger separation, we detect an outer disk structure that can be fitted with an ellipse, compatible with a circular ring seen at r = 0.62" (190 au) and height of 77 au. After applying a height-corrected deprojection we see a circular ring centered on the star at 88 au, while what seemed to be a separate blob and outer ring, now both could be part of a single-armed spiral. Based on the current data it is not possible to conclude decisively whether $H_{\rm scat} / r$ remains constant or whether the surface is flared with at most $H_{\rm scat} \propto r^{1.35}$ , although we favor the constant ratio based on our deprojections. The height-corrected deprojection allows a more detailed interpretation of the observed structures, after which we discern the detection of a single-armed spiral. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: 10 pages, 8 figures, accepted for publication by A&A

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

Dust dynamics in planet-driven spirals

Authors: Ardjan Sturm, Giovanni Rosotti, Carsten Dominik

Abstract: Context. Protoplanetary disks are known to host spiral features that are observed in scattered light, ALMA continuum and more recently in CO gas emission and gas dynamics. It is however unknown if spirals in gas and dust trace the same morphology. Aims. We study the morphology and amplitude of dusty spirals as function of Stokes number and the underlying mechanisms causing a difference from gas… ▽ More Context. Protoplanetary disks are known to host spiral features that are observed in scattered light, ALMA continuum and more recently in CO gas emission and gas dynamics. It is however unknown if spirals in gas and dust trace the same morphology. Aims. We study the morphology and amplitude of dusty spirals as function of Stokes number and the underlying mechanisms causing a difference from gas spirals. We then construct a model to relate the deviation from Keplerian rotation in the gas to a perturbation in surface density of gas and dust. Methods. We use FARGO-3D with dust implementation to numerically study the spirals, after which the results are interpreted using a semi-analytical model. This model is tested on observational data to predict the perturbation of the spiral in gas dynamics based on the continuum data. Results. We find that the pitch angle of a spiral does not differ significantly between gas and dust. The amplitude of the dust spiral decreases with Stokes number (St) and starts to fade out at a typical St > 0.1 as the dust becomes decoupled from the gas. The semianalytical model provides an accurate and fast representation of the spiral in surface density of the dust from the gas. We find a spiral in the TW Hya velocity residual map, never seen before, which is a feature in the vertical velocity and has a kink at the continuum gap, yielding strong evidence for a planet at 99 au. Conclusions. We built a model that gives an estimate of the underlying dynamics of dust in a spiral, which can serve as a proof of planetary origin of spirals and can be a probe for the Stokes number in the disk. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 11 pages, 8 figures, appendices, accepted to A&A 14 September 2020

Journal ref: A&A 643, A92 (2020)

Interpreting high spatial resolution line observations of planet-forming disks with gaps and rings -- The case of HD 163296

Authors: Ch. Rab, I. Kamp, C. Dominik, C. Ginski, G. A. Muro-Arena, W. -F. Thi, L. B. F. M. Waters, P. Woitke

Abstract: Spatially resolved continuum observations of planet-forming disks show prominent ring and gap structures in their dust distribution. However, the picture from gas observations is much less clear and constraints on the radial gas density structure (i.e. gas gaps) remain rare and uncertain. We want to investigate the importance of thermo-chemical processes for the interpretation of high-spatial-reso… ▽ More Spatially resolved continuum observations of planet-forming disks show prominent ring and gap structures in their dust distribution. However, the picture from gas observations is much less clear and constraints on the radial gas density structure (i.e. gas gaps) remain rare and uncertain. We want to investigate the importance of thermo-chemical processes for the interpretation of high-spatial-resolution gas observations of planet-forming disks and their impact on derived gas properties. We apply the radiation thermo-chemical disk code ProDiMo (PROtoplanetary DIsk MOdel) to model self-consistently the dust and gas disk of HD 163296, using the DSHARP gas and dust observations. With this model we investigate the impact of dust gaps and gas gaps, considering chemistry and heating/cooling processes, on the observables and the derived gas properties. We find distinct peaks in the radial line intensity profiles of the CO line data of HD 163296 at the location of the dust gaps. Our model indicates that those peaks are not only a consequence of a gas temperature increase within the gaps but are mainly caused by the absorption of line emission from the back side of the disk by the dust rings. For two of the three prominent dust gaps in HD 163296, we find that thermo-chemical effects are negligible for deriving density gradients via measurements of the rotation velocity. However, for the gap with the highest dust depletion, the temperature gradient can be dominant and needs to be considered to derive accurate gas density profiles. Self-consistent gas and dust thermo-chemical modelling in combination with high-quality observations of multiple molecules are necessary to accurately derive gas gap depths and shapes. This is crucial to determine the origin of gaps and rings in planet-forming disks and to improve the mass estimates of forming planets if they are the cause of the gap. △ Less

Submitted 3 September, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

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

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