Dust Continuum Emission Research Articles

Multi-wavelength high-resolution polarimetric imaging of second-generation disc around post-AGB binary IRAS 08544-4431 with SPHERE.

Abstract Exploring the formation and evolution of second-generation circumbinary discs around evolved binary stars, such as post-Asymptotic Giant Branch (post-AGB) and post-Red Giant Branch (post-RGB) binaries, provides valuable insights into the complex binary interaction process that concludes the red-giant phase of evolution in these systems. Additionally, it offers a novel opportunity to investigate the formation of second-generation planets within dusty discs surrounding evolved stars. We present a pilot multi-wavelength polarimetric imaging study of the post-AGB binary system IRAS 08544-4431 using the European Southern Observatory-Very Large Telescope/SPHERE instrument. This study is focused on optical V − and I′ −band ZIMPOL data to complement near-infrared H −band IRDIS data presented previously. The study aims to investigate the dust scattering properties and surface morphology of the post-AGB circumbinary disc as a function of wavelength. We successfully resolved the extended disc structure of IRAS 08544-4431, revealing a complex disc morphology, high polarimetric disc brightness (up to ∼1.5%), and significant forward scattering at optical wavelengths. Additionally, we found that the disc shows a grey polarimetric colour in both optical and near-infrared. The findings highlight similarities between post-AGB circumbinary discs and protoplanetary discs, suggesting submicron-size porous aggregates as the dominant surface dust composition, and indicating potential warping within the disc. However, further expansion of the multi-wavelength analysis to a larger sample of post-AGB binary systems, as well as high-resolution observations of dust continuum and gas emission, is necessary to fully explore the underlying structure of post-AGB circumbinary discs and associated physical mechanisms.

REBELS-25: Discovery of a dynamically cold disc galaxy at z = 7.31

Abstract We present high resolution (∼0.14” = 710 pc) ALMA [CII] 158μm and dust continuum follow-up observations of REBELS-25, a [CII]-luminous (L[CII] = (1.7 ± 0.2) × 109L⊙) galaxy at redshift z = 7.3065 ± 0.0001. These high resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{{\rm M}_{\odot }}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{{\rm M}_{\odot }} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionisation. By modelling the kinematics with 3DBAROLO, we find it has a low velocity dispersion ($\bar{\sigma } = 33^{+9}_{-7}$ km s−1) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\bar{\sigma } = 11 ^{+6}_{-5}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [CII] distribution is well fit by a near-exponential disc model, with a Sérsic index, n, of 1.3 ± 0.2, and we see tentative evidence of more complex non-axisymmetric structures suggestive of a bar in the [CII] and dust continuum emission. By comparing to other high spatial resolution cold gas kinematic studies, we find that dynamically cold discs seem to be more common in the high redshift Universe than expected based on prevailing galaxy formation theories, which typically predict more turbulent and dispersion-dominated galaxies in the early Universe as an outcome of merger activity, gas accretion and more intense feedback. This higher degree of rotational support seems instead to be consistent with recent cosmological simulations that have highlighted the contrast between cold and warm ionised gas tracers, particularly for massive galaxies. We therefore show that dynamically settled disc galaxies can form as early as 700 Myr after the Big Bang.

The Spatial Correlation between CN-line- and Dust-continuum-emitting Regions in High-mass Star-forming Clouds

Measuring the strength of a three-dimensional (3D) magnetic field vector is challenging as it is not easy to recognize whether its line-of-sight (LOS) and plane-of-sky (POS) components are obtained from the same region. CN (N = 1–0) emission has been used to get the LOS component of a magnetic field (B LOS) from its Zeeman splitting lines, while dust continuum emission has been used to get the POS component of a magnetic field (B POS). We use the CN (N = 1–0) data observed with the Taeduk Radio Astronomy Observatory 14 m telescope and the dust continuum data from the Herschel archive toward six high-mass star-forming regions in order to test whether CN line and dust continuum emission can trace a similar region and thus can be used for inferring 3D magnetic field strength. Our comparison between CN and H2 column densities for all targets indicates that CN line emission tends to be strong toward bright continuum regions. The positions of peak CN column densities are particularly well correlated with those of peak H2 column densities, at least over the H2 column density of 8.0 × 1022 cm−2 within one or two telescope beam sizes in all targets, implying that CN-line- and dust-continuum-emitting regions are likely spatially coincident. This enabled us to make the reliable measurement of the 3D magnetic field strengths of five targets by taking a vector sum of their B LOS and B POS, helping to decide the magnetical criticality of the targets as supercritical or transcritical.

Early Planet Formation in Embedded Disks. XI. A High-resolution View Toward the BHR 71 Class 0 Protostellar Wide Binary

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary Class 0 protostellar system BHR 71 IRS1 and IRS2 as part of the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program. We describe the 12CO (J = 2–1), 13CO (J = 2–1), C18O (J = 2–1), H2CO (J = 32,1–22,0), and SiO (J = 5–4) molecular lines along with the 1.3 mm continuum at high spatial resolution (∼0.″08 or ∼5 au). Dust continuum emission is detected toward BHR 71 IRS1 and IRS2, with a central compact component and extended continuum emission. The compact components are smooth and show no sign of substructures such as spirals, rings, or gaps. However, there is a brightness asymmetry along the minor axis of the presumed disk in IRS1, possibly indicative of an inclined geometrically and optically thick disk-like component. Using a position–velocity diagram analysis of the C18O line, clear Keplerian motions were not detected toward either source. If Keplerian rotationally supported disks are present, they are likely deeply embedded in their envelope. However, we can set upper limits of the central protostellar mass of 0.46 M ⊙ and 0.26 M ⊙ for BHR 71 IRS1 and BHR 71 IRS2, respectively. Outflows traced by 12CO and SiO are detected in both sources. The outflows can be divided into two components, a wide-angle outflow and a jet. In IRS1, the jet exhibits a double helical structure, reflecting the removal of angular momentum from the system. In IRS2, the jet is very collimated and shows a chain of knots, suggesting episodic accretion events.

Hints of Planet Formation Signatures in a Large-cavity Disk Studied in the AGE-PRO ALMA Large Program

Detecting planet signatures in protoplanetary disks is fundamental to understanding how and where planets form. In this work, we report dust and gas observational hints of planet formation in the disk around 2MASS J16120668-301027, as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Large Program “AGE-PRO: ALMA survey of Gas Evolution in Protoplanetary disks.” The disk was imaged with the ALMA at Band 6 (1.3 mm) in dust continuum emission and four molecular lines: 12CO(J = 2–1), 13CO(J = 2–1), C18O(J = 2–1), and H2CO(J = 3(3,0)–2(2,0)). Resolved observations of the dust continuum emission (angular resolution of ∼150 mas, 20 au) show a ring-like structure with a peak at 0.″57 (75 au), a deep gap with a minimum at 0.″24 (31 au), an inner disk, a bridge connecting the inner disk and the outer ring, along with a spiral arm structure, and a tentative detection (to 3σ) of a compact emission at the center of the disk gap, with an estimated dust mass of ∼2.7−12.9 Lunar masses. We also detected a kinematic kink (not coincident with any dust substructure) through several 12CO channel maps (angular resolution ∼200 mas, 30 au), located at a radius of ∼0.″875 (115.6 au). After modeling the 12CO velocity rotation around the protostar, we identified a purple tentative rotating-like structure at the kink location with a geometry similar to that of the disk. We discuss potential explanations for the dust and gas substructures observed in the disk and their potential connection to signatures of planet formation.

A3COSMOS: Dust mass function and dust mass density at 0.5 < z < 6

Context. Although dust in galaxies represents only a few percent of the total baryonic mass, it plays a crucial role in the physical processes occurring in galaxies. Studying the dust content of galaxies, particularly at high z, is therefore crucial for understanding the link between dust production, obscured star formation, and the build-up of galaxy stellar mass. Aims. We study the dust properties (mass and temperature) of the largest Atacama Large Millimeter/submillimeter Array (ALMA)-selected sample of star-forming galaxies available from the archive (A3COSMOS), and we derive the dust mass function and dust mass density of galaxies from z = 0.5 − 6. Methods. We fit the spectral energy distribution (SED) with the CIGALE code to constrain the dust mass and temperature of the A3COSMOS galaxy sample based on the UV-to-near-infrared photometric coverage of each galaxy combined with the ALMA (and Herschel when available) coverage of the Rayleigh-Jeans tail of their dust-continuum emission. We then computed and fit the dust mass function by combining the A3COSMOS and the most recent Herschel samples in order to obtain the best estimate of the integrated dust mass density up to z ∼ 6. Results. The dust masses in galaxies in A3COSMOS lie between ∼108 and ∼109.5 M⊙. From the SED fitting, we were also able to derive a dust temperature. The distribution of the dust temperature peaks at ∼30 − 35 K. The dust mass function at z = 0.5 − 6 evolves with an increase in M* and a decrease in the number density (Φ*), and it agrees well with literature estimates. The dust mass density decreases smoothly in its evolution from z ∼ 0.5 to z ∼ 6, which is steeper than what is found by models at z ≳ 2.

JWST Observations of Young protoStars (JOYS). HH211: Textbook case of a protostellar jet and outflow

Due to the high visual extinction and lack of sensitive mid-infrared (MIR) telescopes, the origin and properties of outflows and jets from embedded Class\,0 protostars are still poorly constrained. We aim to characterise the physical, kinematic, and dynamical properties of the HH\,211 jet and outflow, one of the youngest protostellar flows. We used the James Webb Space Telescope (JWST) and its Mid-InfraRed Instrument (MIRI) in the 5--28 mu m range to study the embedded HH\,211 flow. We mapped a 0farcm 95times 0farcm 22 region, covering the full extent of the blueshifted lobe, the central protostellar region, and a small portion of the redshifted lobe. We extracted spectra along the jet and outflow and constructed line and excitation maps of both atomic and molecular lines. Additional JWST NIRCam H$_2$ narrow-band images (at 2.122 and 3.235\,mu m) provide a visual-extinction map of the whole flow, and are used to deredden our data. The jet-driving source is not detected even at the longest MIR wavelengths. The overall morphology of the flow consists of a highly collimated jet, which is mostly molecular (H$_2$, HD) with an inner atomic ( structure. The jet shocks the ambient medium, producing several large bow shocks (BSs) that are rich in forbidden atomic ( and molecular lines (H$_2$, HD, CO, OH, H$_2$O, CO$_2$, HCO$^+$), and is driving an H$_2$ molecular outflow that is mostly traced by low-$J$, $v=0$ transitions. Moreover, H$_2$ 0-0\,S(1) uncollimated emission is also detected down to 2arcsec --3arcsec (sim 650--1000\,au) from the source, tracing a cold ($T$=200--400\,K), less dense, and poorly collimated molecular wind. Two H$_2$ components (warm, $T$=300--1000\,K, and hot, $T$=1000--3500\,K) are detected along the jet and outflow. The atomic jet ( at 26\,mu m) is detected down to sim 130\,au from the source, whereas the lack of H$_2$ emission (at 17\,mu m) close to the source is likely due to the large visual extinction ( Dust-continuum emission is detected at the terminal BSs and in the blue- and redshifted jet, and is likely attributable to dust lifted from the disc. The jet shows an onion-like structure, with layers of different size, velocity, temperature, and chemical composition. Moreover, moving from the inner jet to the outer BSs, different physical, kinematic, and excitation conditions for both molecular and atomic gas are observed. The mass-flux rate and momentum of the jet, as well as the momentum flux of the warm H$_2$ component, are up to one order of magnitude higher than those inferred from the atomic jet component. Our findings indicate that the warm H$_2$ red component is the main driver of the outflow, that is to say it is the most significant dynamical component of the jet, in contrast to jets from more evolved YSOs, where the atomic component is dominant.

The ALMA-CRISTAL survey

We present the morphological parameters and global properties of dust-obscured star formation in typical star-forming galaxies at z = 4–6. Among 26 galaxies composed of 20 galaxies observed by the Cycle-8 ALMA Large Program, CRISTAL, and 6 galaxies from archival data, we individually detect rest-frame 158 μm dust continuum emission from 19 galaxies, 9 of which are reported for the first time. The derived far-infrared luminosities are in the range log10LIR [L⊙] = 10.9 − 12.4, an order of magnitude lower than previously detected massive dusty star-forming galaxies (DSFGs). We find the average relationship between the fraction of dust-obscured star formation (fobs) and the stellar mass to be consistent with previous results at z = 4–6 in a mass range of log10M* [M⊙]∼9.5 − 11.0 and to show potential evolution from z = 6 − 9. The individual fobs exhibits significant diversity, and we find a potential correlation with the spatial offset between the dust and UV continuum, suggesting that inhomogeneous dust reddening may cause the source-to-source scatter in fobs. The effective radii of the dust emission are on average ∼1.5 kpc and are about two times more extended than those seen in rest-frame UV. The infrared surface densities of these galaxies (ΣIR ∼ 2.0 × 1010 L⊙ kpc−2) are one order of magnitude lower than those of DSFGs that host compact central starbursts. On the basis of the comparable contribution of dust-obscured and dust-unobscured star formation along with their similar spatial extent, we suggest that typical star-forming galaxies at z = 4 − 6 form stars throughout the entirety of their disks.

Constraints on the Physical Origin of Large Cavities in Transition Disks from Multiwavelength Dust Continuum Emission

The physical origin of the large cavities observed in transition disks is to date still unclear. Different physical mechanisms (e.g., a companion, dead zones, enhanced grain growth) produce disk cavities of different depth, and the expected spatial distribution of gas and solids in each mechanism is not the same. In this work, we analyze the multiwavelength interferometric visibilities of dust continuum observations obtained with Atacama Large Millimeter/submillimeter Array and Very Large Array for six transition disks: CQTau, UXTau A, LkCa15, RXJ1615, SR24S, and DMTau, and calculate brightness radial profiles, where diverse emission morphology is revealed at different wavelengths. The multiwavelength data are used to model the spectral energy distribution and compute constraints on the radial profile of the dust surface density, maximum grain size, and dust temperature in each disk. They are compared with the observational signatures expected from various physical mechanisms responsible for disk cavities. The observational signatures suggest that the cavities observed in the disks around UXTau A, LkCa15, and RXJ1615 could potentially originate from a dust trap created by a companion. Conversely, in the disks around CQTau, SR24S, DMTau, the origin of the cavity remains unclear, although it is compatible with a pressure bump and grain growth within the cavity.

Rossby wave instability and substructure formation in 3D non-ideal MHD wind-launching discs

ABSTRACT Rings and gaps are routinely observed in the dust continuum emission of protoplanetary discs (PPDs). How they form and evolve remains debated. Previous studies have demonstrated the possibility of spontaneous gas rings and gaps formation in wind-launching discs. Here, we show that such gas substructures are unstable to the Rossby wave instability (RWI) through numerical simulations. Specifically, shorter wavelength azimuthal modes develop earlier, and longer wavelength ones dominate later, forming elongated (arc-like) anticyclonic vortices in the rings and (strongly magnetized) cyclonic vortices in the gaps that persist until the end of the simulation. Highly elongated vortices with aspect ratios of 10 or more are found to decay with time in our non-ideal magnetohydrodynamic (MHD) simulation, in contrast with the hydro case. This difference could be caused by magnetically induced motions, particularly strong meridional circulations with large values of the azimuthal component of the vorticity, which may be incompatible with the columnar structure preferred by vortices. The cyclonic and anticyclonic RWI vortices saturate at moderate levels, modifying but not destroying the rings and gaps in the radial gas distribution of the disc. In particular, they do not shut-off the poloidal magnetic flux accumulation in low-density regions and the characteristic meridional flow patterns that are crucial to the ring and gap formation in wind-launching discs. Nevertheless, the RWI and their associated vortices open up the possibility of producing non-axisymmetric dust features observed in a small fraction of PPDs through non-ideal MHD, although detailed dust treatment is needed to explore this possibility.

Orbital dynamics in the GG Tau A system: Investigating its enigmatic disc

Context. GG Tau is one of the most studied young multiple stellar systems: GG Tau A is a hierarchical triple surrounded by a massive disc and its companion, GG Tau B, is also a binary. Despite numerous observational attempts, a comprehensive understanding of the geometry of the GG Tau A system is still elusive. Given the significant role of dynamical interactions in shaping the evolution of these systems, it is relevant to characterise the stellar orbits and the discs’ properties. Aims. To determine the best orbital configuration of the GG Tau A system and its circumtriple disc, we provide new astrometric measures of the system and we run a set of hydrodynamical simulations with two representative orbits to test how they impact a disc composed of dust and gas. Methods. We tested the dynamical evolution of the two scenarios on short and long timescales. We obtained synthetic flux emission from our simulations at different timescales and we compared them with multi-wavelength observations of 1300 µm ALMA dust continuum emission and 1.67 µm SPHERE dust scattering to infer the most likely orbital arrangement. Results. We extend the analysis of the binary orbital parameters using six new epochs from archival data, showing that the current measurements alone (and future observations coming in the next 5–10 yr) are not capable of fully breaking the degeneracy between families of coplanar and misaligned orbits, but finding that a modest misalignment is probable. We find that the timescale for the onset of the disc eccentricity growth, τecc, is a fundamental timescale for the morphology of the system. Results from the numerical simulations obtained using the representative coplanar and misaligned (∆θ = 30°) orbits show that the best match between the position of the stars, the cavity size, and the dust ring size of GG Tau A is obtained with the misaligned configuration on timescales shorter than τecc. The results exhibit an almost circular cavity and dust ring, favouring slightly misaligned (∆θ ~ 10–30°) low-eccentricity (e ~ 0.2–0.4) orbits. However, for both scenarios, the cavity size and its eccentricity quickly grow for timescales longer than τecc and the models do not reproduce the observed morphology anymore. This implies that either the age of the system is shorter than τecc or that the disc eccentricity growth is not triggered or dissipated in the system. This finding raises questions about the future evolution of the GG Tau A system and, more generally, the time evolution of eccentric binaries and their circumbinary discs.

Constraints on PDS 70 b and c from the dust continuum emission of the circumplanetary discs considering in situ dust evolution

Context. The young T Tauri star PDS 70 has two gas accreting planets sharing one large gap in a pre-transitional disc. Dust continuum emission from PDS 70 c has been detected by Atacama Large Millimeter/submillimeter Array (ALMA) Band 7, considered as the evidence of a circumplanetary disc. However, there has been no detection of the dust emission from the CPD of PDS 70 b. Aims. We constrain the planet mass and the gas accretion rate of the planets by introducing a model of dust evolution in the CPDs and reproducing the detection and non-detection of the dust emission. Methods. We first develop a 1D steady gas disc model of the CPDs reflecting the planet properties. We then calculate the radial distribution of the dust profiles considering the dust evolution in the gas disc and calculate the total flux density of dust thermal emission from the CPDs. Results. We find positive correlations between the flux density of dust emission and three planet properties, the planet mass, gas accretion rate, and their product called ‘MMdot’. We then find that the MMdot of PDS 70 c is ≥4 × 10−7 MJ2 yr−1, corresponding to the planet mass of ≥5 MJ and the gas accretion rate of ≥2 × 10−8 MJ yr−1. This is the first case to succeed in obtaining constraints on planet properties from the flux density of dust continuum emission from a CPD. We also find some loose constraints on the properties of PDS 70 b from the non-detection of its dust emission. Conclusions. We propose possible scenarios for PDS 70 b and c explaining the non-detection respectively detection of the dust emission from their CPDs. The first explanation is that planet c has larger planet mass, larger gas accretion rate, or both than planet b. The other possibility is that the CPD of planet c has a larger amount of dust supply, weaker turbulence, or both than that of planet b. If the dust supply to planet c is larger than b due to its closeness to the outer dust ring, it is also quantitatively consistent with that planet c has weaker Hα line emission than planet b considering the dust extinction effect.

Bayesian analysis of the molecular emission and dust continuum of protoplanetary disks

Context. The MIRI instrument on board the James Webb Space Telescope probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail and reveal a rich chemistry with strong diversity between objects. This complicates interpretations that are mainly based on manual continuum subtraction and 0D slab models. Aims. We investigate the physical conditions under which the gas emits in protoplanetary disks. Based on MIRI spectra, we apply a full Bayesian analysis that provides the posterior distributions of dust and molecular properties, such as column densities and emission temperatures. Methods. To do so, we introduced the Dust Continuum Kit with Line emission from Gas (DuCKLinG), a Python-based model simultaneously describing the molecular line emission and the dust continuum of protoplanetary disks without large computational cost. The model describes the dust continuum emission by dust models with precomputed dust opacities. The molecular emission is based on LTE slab models but from extended radial ranges with gradients in column densities and emission temperatures. We compare the model to observations using Bayesian analysis with linear regression techniques to reduce the dimension of the parameter space. We benchmarked this model to a complex thermo-chemical ProDiMo model of AATau and fit the MIRI spectrum of GW Lup. The latter allowed for a comparison to the previous results obtained with single slab models and hand-fitted continuum. Results. We successfully decrease the computational time of the fitting method by a factor of 80 by eliminating linear parameters, such as the emission areas, from the Bayesian run. This approach does not significantly change the retrieved molecular parameters, and only the calculated errors on the optically thin dust masses slightly decrease. For an AA Tau ProDiMo mock observation, we find that the retrieved molecular conditions from DuCKLinG (column densities from 3 × 1018 cm−2 to 4 × 1020 cm−2, radial range from 0.2 au to 1.2 au, and temperature range from about 200 K to 400 K) fall within the true values from ProDiMo (column densities between 4 × 1017 cm-2 to 5 × 1020 cm−2, radial extent 0.1 au to 6.6 au, and temperature range from about 120 to 1000 K). The smaller DuCKLinG ranges can be explained by the relative flux contributions of the different parts of ProDiMo. The parameter posterior of GW Lup reinforces previously found results. The previously determined column densities fall within the retrieved ranges in this study for all examined molecules (CO2, H2O, HCN, and C2H2). Similar overlap is found for the temperatures with only the temperature range of HCN (from 570−60+60 to 750−70+90 K) not including the previously found value (875 K). This discrepancy may be due to the simultaneous fitting of all molecules compared to the step-by-step fitting of the previous study. There is statistically significant evidence for radial temperature and column density gradients for H2O and CO2 compared to the constant temperature and column density assumed in the 0D slab models. Additionally, HCN and C2H2 emit from a small region with near constant conditions. Due to the small selected wavelength range 13.6–16.3 µm, the dust properties are not well constrained for GW Lup. DuCKL inG can become an important tool to analyse the molecular emission and dust mineralogy of large samples based on JWST /MIRI spectra in an automated way.

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

The interplay between T\,Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. In the last years, major progress has been made using instrumentation that probes the dust structure in the mid-plane and at the surface of protoplanetary disks. Observations show a great variety of disk shapes and substructures that are crucial for understanding planet formation. We studied a seemingly old T\,Tauri star, PDS\,111, and its disk. We combined complementary observations of the stellar atmosphere, the circumstellar hot gas, the surface of the disk, and the mid-plane structure. 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 to update its stellar parameters, 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 Halpha variability and with a low mass-accretion rate of $1-5 odot $. We measured an age of the system of 15.9$^ $\,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_ or $ Jup $. The 12CO observations do not show a cavity and the 12CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than often suggested in literature, PDS\,111 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, that "refuse to evolve".

JWST and ALMA Discern the Assembly of Structural and Obscured Components in a High-redshift Starburst Galaxy

We present observations and analysis of the starburst PACS-819 at z = 1.45 (M * = 1010.7 M ⊙), using high-resolution (0.″1; 0.8 kpc) Atacama Large Millimeter/submillimeter Array (ALMA) and multiwavelength JWST images from the COSMOS-Web program. Dissimilar to Hubble Space Telescope (HST) ACS images in the rest-frame UV, the redder NIRCam and MIRI images reveal a smooth central mass concentration and spiral-like features, atypical for such an intense starburst. Through dynamical modeling of the CO (J = 5–4) emission with ALMA, PACS-819 is rotation dominated and thus consistent with having a disk-like nature. However, kinematic anomalies in CO and asymmetric features in the bluer JWST bands (e.g., F150W) support a more disturbed nature likely due to interactions. The JWST imaging further enables us to map the distribution of stellar mass and dust attenuation, thus clarifying the relationships between different structural components not discernible in the previous HST images. The CO (J = 5–4) and far-infrared dust continuum emission are cospatial with a heavily obscured starbursting core (<1 kpc) that is partially surrounded by much less obscured star-forming structures including a prominent arc, possibly a tidally distorted dwarf galaxy, and a massive clump (detected in CO), likely a recently accreted low-mass satellite. With spatially resolved maps, we find a high molecular gas fraction in the central area reaching ∼3 (M gas/M *) and short depletion times (M gas/SFR ∼ 120 Myr, where SFR is star formation rate) across the entire system. These observations provide insights into the complex nature of starbursts in the distant Universe and underscore the wealth of complementary information from high-resolution observations with both ALMA and JWST.

The ALPINE-ALMA [CII] survey: Dust emission effective radius up to 3 kpc in the early Universe

Aims. Measurements of the size of dust continuum emission are an important tool for constraining the spatial extent of star formation, and hence the buildup of stellar mass. Compact dust emission has generally been observed at cosmic noon (z ∼ 2 − 3). However, at earlier epochs, toward the end of the reionization (z ∼ 4 − 6), only the sizes of a handful of infrared (IR) bright galaxies have been measured. In this work, we derive the dust emission sizes of main-sequence (MS) galaxies at z ∼ 5 from the ALPINE survey. Methods. We measured the dust effective radius, re, FIR, in the uv-plane in Band 7 of ALMA for seven ALPINE galaxies with resolved emission and we compared it with rest-frame ultraviolet (UV) and [CII]158 μm measurements. We studied the re, FIR − LIR scaling relation by considering our dust size measurements and all the data in the literature at z ∼ 4 − 6. Finally, we compared our size measurements with predictions from simulations. Results. The dust emission in the selected ALPINE galaxies is rather extended (re, FIR ∼ 1.5 − 3 kpc), similar to [CII]158 μm but a factor of ∼2 larger than the rest-frame UV emission. Putting together all the measurements at z ∼ 5, spanning two decades in luminosity from LIR ∼ 1011 L⊙ to LIR ∼ 1013 L⊙, the data highlight a steeply increasing trend of the re, FIR − LIR relation at LIR &lt; 1012 L⊙, followed by a downturn and a decreasing trend at brighter luminosities. Finally, simulations that extend up to the stellar masses of the ALPINE galaxies considered in the present work predict a subset of galaxies (∼25% at 1010 M⊙ &lt; M⋆ &lt; 1011 M⊙) with sizes as large as those measured.

Self-similarity of the magnetic field at different scales: The case of G31.41+0.31

Context. Dust polarization observations of the massive protocluster G31.41+0.31 carried out at ~1″ (~3750 au) resolution with the SMA at 870 µm have revealed one of the clearest examples to date of an hourglass-shaped magnetic field morphology in the high-mass regime. Additionally, ~O.″24 (~900 au) resolution observations with ALMA at 1.3 mm have confirmed these results. The next step is to investigate whether the magnetic field maintains its hourglass-shaped morphology down to circumstellar scales. Aims. To study the magnetic field morphology toward the four (proto)stars A, B, C, and D contained in G31.41+0.31 and examine whether the self-similarity observed at core scales (1″ and 0.″ 24 resolution) still holds at circumstellar scales, we carried out ALMA observations of the polarized dust continuum emission at 1.3 mm and 3.1 mm at an angular resolution of ~0.″068 (~250 au), sufficient to resolve the envelope emission of the embedded protostars. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6) and 97.5 GHz (Band 3) with a synthesized beam of 0.″072 × 0.″064. We carried out polarization observations at two different wavelengths to confirm that the polarization traces magnetically aligned dust grains and is not due to dust self-scattering. Results. The polarized emission and the direction of the magnetic field obtained at the two wavelengths are basically the same, except for an area between the embedded sources C and B. In such an area, the emission at 1.3 mm could be optically thick and affected by dichroic extinction. In the rest of the core, the similarity of the emission at the two wavelengths suggests that the polarized emission is due to magnetically aligned grains. The polarized emission has been successfully modeled with a poloidal field with a small toroidal component on the order of 10% of the poloidal component, with a position angle ϕ = −63°, an inclination i = 50°, and a mass-to-flux ratio λ = 2.66. The magnetic field axis is oriented perpendicular to the NE-SW velocity gradient detected in the core. The strength of the plane-of-the-sky component of the mean magnetic field, estimated using both the Davis-Chandrasekhar-Fermi and the polarization-intensity gradient methods, is in the range ~10−80 mG, for a density range 1.4 × 107−5 × 108 cm−3. The mass-to-flux ratio is in the range λ~1.9−3.0, which suggests that the core is “supercritical”. The polarization-intensity gradient method indicates that the magnetic field cannot prevent gravitational collapse inside the massive core. The collapse in the external part of the core is (slightly) sub-Alfvénic and becomes super-Alfvénic close to the center. Conclusions. Dust polarization measurements from large core scales to small circumstellar scales, in the hot molecular core G31.41+0.31 have confirmed the presence of a strong magnetic field with an hourglass-shaped morphology. This result suggests that the magnetic field could have a relevant role in regulating the star-forming process of massive stars at all scales, although it cannot prevent the collapse. However, it cannot be ruled out that the large opacity of the central region of the core may hinder the study of the magnetic field at circumstellar scales. Therefore, high-angular resolution observations at longer wavelengths, tracing optically thinner emission, are needed to confirm this self-similarity.

Angular momentum transport via gravitational instability in the Elias 2–27 disc

Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2−27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used archival ALMA observations of 13CO line emission to measure the efficiency of angular momentum transport in the Elias 2−27 system through the kinematic signatures generated by gravitational instability, known as “GI wiggles”. Assuming the angular momentum is transported by the observed spiral structure and leveraging previously-derived dynamical disc mass measurements, the amount of angular momentum transport we found corresponds to an α-viscosity of α = 0.038 ± 0.018. This value implies an accretion rate onto the central star of log10 Ṁ⋆ = −6.99 ± 0.17 M⊙ yr−1, which reproduces the one observed value of log10 Ṁ⋆,obs = −7.2 ± 0.5 M⊙ yr−1 very well. The excellent agreement we have found serves as further proof that gravitational instability is the main driver of angular momentum transport acting in this system.

PRODIGE - Planet-forming disks in Taurus with NOEMA. II. Modeling the CO (2-1) isotopologue emission of the Class II T Tauri disks in Taurus

To understand how planets form in protoplanetary disks, it is necessary to characterize their gas and dust distribution and masses. This requires a combination of high-resolution dust continuum and molecular line interferometric observations, coupled with advanced theoretical models of protoplanetary disk physics, chemical composition, and radiative transfer. We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at $0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution \, PIs: P. Caselli Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use these data to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C$ $O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. The best-fit disk gas masses range between 0.005 and $0.04 M_ sun $. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements via, for example, modeling of the wavelength dependence of the dust continuum emission radii, and HD and CO isotopologue emission.

The Molecular Clouds of M31

Deep interferometric observations of CO and dust continuum emission are obtained with the Submillimeter Array at 230 GHz to investigate the physical nature of the giant molecular cloud (GMC) population in the Andromeda galaxy (M31). We use J = 2 − 1 12CO and 13CO emission to derive the masses, sizes, and velocity dispersions of 162 spatially resolved GMCs. We perform a detailed study of a subset of 117 GMCs that exhibit simple, single-component line profile shapes. Examining the Larson scaling relations for these GMCs, we find (1) a highly correlated mass–size relation in both 12CO and 13CO emission; (2) a weakly correlated 12CO line width–size (LWS) relation along with a weaker, almost nonexistent, 13CO LWS relation, suggesting a possible dependence of the LWS relation on spatial scale; and (3) that only 43% of these GMCs are gravitationally bound. We identify two classes of GMCs based on the strength and extent of their 13CO emission. Examination of the Larson relations finds that both classes are individually characterized by strong 12CO mass–size relations and much weaker 12CO and 13CO LWS relations. The majority (73%) of strong 13CO-emitting GMCs are found to be gravitationally bound. However, only 25% of the weak 13CO-emitting GMCs are bound. The resulting breakdown in the Larson relations in the weak 13CO-emitting population decouples the mass–size and LWS relations, demonstrating that independent physical causes are required to understand the origin of each. Finally, in nearly every aspect, the physical properties of the M31 GMCs are found to be very similar to those of local Milky Way clouds.