Atacama Large Millimeter/submillimeter Array Research Articles

Complex angular structure of three elliptical galaxies from high-resolution ALMA observations of strong gravitational lenses

The large-scale mass distributions of galaxy-scale strong lenses have long been assumed to be well described by a singular ellipsoidal power-law density profile with external shear. However, the inflexibility of this model could lead to systematic errors in astrophysical parameters inferred with gravitational lensing observables. Here, we present observations with the Atacama Large (sub-)Millimetre Array (ALMA) of three strongly lensed dusty star-forming galaxies at $ mas angular resolution and investigate the sensitivity of these data to angular structure in the lensing galaxies. We jointly infer the lensing mass distribution and the full surface brightness of the lensed sources with multipole expansions of the power-law density profile up to the fourth order using a technique developed for interferometric data. All three datasets strongly favour third and fourth-order multipole amplitudes of $ percent of the convergence. While the infrared stellar isophotes and isodensity shapes agree for one lens system, for the other two the isophotes disagree to varying extents, suggesting contributions to the angular structure from dark matter intrinsic or extrinsic to the lensing galaxy.

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

We use 0.1″ observations from the Atacama Large Millimeter Array (ALMA), Hubble Space Telescope (HST), and JWST to study young massive clusters (YMCs) in their embedded “infant” phase across the central starburst ring in NGC 3351. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (∼105 M ⊙), small radii (≲ 5 pc), large escape velocities (6–10 km s−1), and short freefall times (0.5–1 Myr). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed H ii region–cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline, ranging from ∼1–2 Myr before cluster formation as starless clumps, to ∼4–6 Myr after as exposed H ii region–cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit a nonuniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows.

CON-quest

Context. Nuclear regions of ultraluminous and luminous infrared galaxies (U/LIRGs) are powered by starbursts and/or active galactic nuclei (AGNs). These regions are often obscured by extremely high columns of gas and dust. Molecular lines in the submillimeter windows have the potential to determine the physical conditions of these compact obscured nuclei (CONs). Aims. We aim to reveal the distributions of HCN and HCO+ emission in local U/LIRGs and investigate whether and how they are related to galaxy properties. Methods. Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted sensitive observations of the HCN J = 3−2 and HCO+J = 3−2 lines toward 23 U/LIRGs in the local Universe (z < 0.07) with a spatial resolution of ∼0.3″ (∼50−400 pc). Results. We detected both HCN and HCO+ in 21 galaxies, only HCN in one galaxy, and neither in one galaxy. The global HCN/HCO+ line ratios, averaged over scales of ∼0.5−4 kpc, range from 0.4 to 2.3, with an unweighted mean of 1.1. These line ratios appear to have no systematic trend with bolometric AGN luminosity or star formation rate. The line ratio varies with position and velocity within each galaxy, with an average interquartile range of 0.38 on a spaxel-by-spaxel basis. In eight out of ten galaxies known to have outflows and/or inflows, we found spatially and kinematically symmetric structures of high line ratios. These structures appear as a collimated bicone in two galaxies and as a thin spherical shell in six galaxies. Conclusions. Non-LTE analysis suggests that the high HCN/HCO+ line ratio in outflows is predominantly influenced by the abundance ratio. Chemical model calculations indicate that the enhancement of HCN abundance in outflows is likely due to high-temperature chemistry triggered by shock heating. These results imply that the HCN/HCO+ line ratio can aid in identifying the outflow geometry when the shock velocity of the outflows is sufficiently high to heat the gas.

The ALMA-QUARKS Survey. II. The ACA 1.3 mm Continuum Source Catalog and the Assembly of Dense Gas in Massive Star-Forming Clumps

Leveraging the high resolution, sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for “Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures”, is observing 139 massive star-forming clumps at ALMA Band 6 (λ ∼ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7 m data of the QUARKS survey, describing the ACA observations and data reduction. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16–84 percentile values of 6–160 M ⊙), gravity-dominated (M ∝ R 1.1) fragments within massive clumps, with supersonic turbulence ( M>1 ) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio (L/M), we find that the DGF increases significantly with L/M, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by a dynamic global collapse process.

Chemistry across dust and gas gaps in protoplanetary disks. Modelling the co-spatial molecular rings in the HD 100546 disk

Chemistry across dust and gas gaps in protoplanetary disks. Modelling the co-spatial molecular rings in the HD 100546 disk

EMISSA: Exploring millimetre indicators of solar-stellar activity III. Comparison of Ca II indices and millimetre continua in a 3D model atmosphere

Amongst several spectral lines, some of the strongest chromospheric diagnostics are offered by the Ca II H K lines. These lines can be used to gauge the temperature stratification of the atmosphere since the line core and wings are formed in different regions of the solar atmosphere. Furthermore, the Ca II lines act as tracers for the magnetic structure of the solar atmosphere, as the line cores are formed in the upper chromosphere even though they are formed in non-local thermodynamic equilibrium (NLTE). In contrast, the formation of millimetre (mm) continuum radiation occurs under local thermodynamic equilibrium (LTE) conditions. As a result, the brightness temperatures obtained from observations with the Atacama Large Millimetre/Submillimetre Array (ALMA) offer a complementary perspective on the activity and thermal structure of stellar atmospheres. The overall aim is to establish more robust solar/stellar activity indicators using ALMA observations in comparison with classical diagnostics, such as the s index and infrared triplet (IRT) index. We employed the 1.5D radiative transfer codes RH1.5D and advanced radiative transfer (ART) to compute the synthetic spectra for the Ca II lines and the millimetre (mm) continua, respectively. These calculations were performed using an enhanced network atmosphere model, which incorporates non-equilibrium hydrogen ionisation generated by the state-of-the-art 3D radiation magnetohydrodynamics (rMHD) Bifrost code. To account for the limited spatial resolution of ALMA, we simulated the effect using a Gaussian point spread function (PSF). Additionally, we analysed the correlations and slopes of scatter plots between the Ca II indices and mm continuum for the original and degraded resolutions, focusing on the entire simulation box, quiet Sun regions, and enhanced network patches separately. The activity indices generated from these lines could further be used to compare the spectra of Sun-like stars with the solar spectrum. We present a comparative study between synthetic continuum brightness temperature maps at mm wavelengths (0.3 mm to 8.5 mm) and the Ca II activity indices; namely, the s index and infrared triplet (IRT) index. The Ca II activity indices and mm brightness temperatures are weakly correlated at the high resolution, with the highest correlation observed at a wavelength of 0.3 mm, corresponding to ALMA band 10. As the resolution decreases, the correlation consistently increases. Conversely, the slopes exhibit a decreasing trend with increasing wavelength, while the degradation of resolution does not noticeably affect the calculated slopes. As the spatial resolution decreases, the standard deviations of the Ca II activity indices and brightness temperatures decrease, while the correlations between them increase. However, the slopes do not exhibit significant changes. Consequently, these relationships could be valuable for calibrating the mm continuum maps obtained through ALMA observations.

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

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. 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. 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. 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 $ M_ J yr $, corresponding to the planet mass of $ J $ and the gas accretion rate of $ M_ J yr $. 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. 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 Halpha line emission than planet b considering the dust extinction effect.

Studying [C ii] emission in low-mass galaxies at z ∼ 7

ABSTRACT We report on a $\rm {[CII]}_{158\mu \rm {m}}$ search using the Atacama Large Millimeter/submillimeter Array (ALMA) on three lensed, confirmed Ly α emitting galaxies at z ∼ 7. Our targets are ultraviolet (UV) faint systems with stellar masses on the order of M* ∼ 109 M⊙. We detect a single [C ii] line emission (4σ) from the brightest (L ∼ 2.4 × 1010L⊙) galaxy in our sample, MACS0454-1251. We determine a systemic redshift (z[C ii] = 6.3151 ± 0.0005) for MACS0454-1251 and measure a Ly α velocity offset of $\Delta v \approx 300 \pm 70 \rm {km\, s}^{-1}$. The remaining two galaxies we detect no [C ii] but provide 3σ upper limits on their [C ii] line luminosities which we use to investigate the $L_{\textrm {[CII]}} - \rm {SFR}$ relation. Overall our single [C ii] detection shows agreement with the relation for local dwarf galaxies. Our [C ii] deficient galaxies could potentially be exhibiting low metallicities (Z < Z⊙). Another possible explanation for weaker [C ii] emission could be strong feedback from star formation disrupting molecular clouds. We do not detect continuum emission in any of the sources, placing upper limits on their dust masses. Assuming a single dust temperature of $T_{d}=35 \rm {K}$ dust masses (Mdust) range from <4.8 × 107 M⊙ to 2.3 × 108 M⊙. Collectively, our results suggest faint reionization era sources could be metal poor and/or could have strong feedback suppressing [C ii] emission.

Probing the dynamical and kinematical structures of detached shells around AGB stars

The chemical evolution of asymptotic giant branch (AGB) stars is driven by repeated thermal pulses (TPs). The duration of a TP is only a few hundred years, whereas an inter-pulse period lasts $10^4-10^5$ years. Direct observations of TPs are hence unlikely. However, the detached shells seen in CO line emission that are formed as a result of a TP provide indirect constraints on the changes experienced by the star during the pulse. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. We used observations of the emission towards five carbon-AGB stars with ALMA (Atacama Large Millimeter/submillimeter Array), including previously published observations of the carbon AGB star U Ant. The data have angular resolutions of 0 to 1 and a velocity resolution of 0.3\ This enabled us to quantify spatial and kinematic structures in the shells. Combining the ALMA data with single-dish observations of the to emission towards the sources, we used radiative transfer models to compare the observed structures with previous estimates of the shell masses and temperatures. The observed emission is separated into two distinct components: a more coherent, bright outer shell and a more filamentary, fainter inner shell. The kinematic information shows that the inner sub-shells move at a higher velocity relative to the outer sub-shells. The observed sub-structures reveal a negative velocity gradient outwards across the detached shells, confirming the predictions from hydrodynamical models. However, the models do not predict a double-shell structure, and the CO emission likely only traces the inner and outer edges of the shell, implying a lack of CO in the middle layers of the detached shell. Previous estimates of the masses and temperatures are consistent with originating mainly from the brighter subshell, but the total shell masses are likely lower limits. Also, additional structures in the form of partial shells outside the detached shell around V644 Sco, arcs within the shell of R Scl, and a partially filled shell for DR Ser indicate a more complicated evolution of the shells and mass-loss process throughout the TP cycle than previously assumed. The observed spatial and kinematical splittings of the shells appear consistent with results from the hydrodynamical models, provided the CO emission does not trace the density distribution in the shell but rather traces the edges of the shells. The hydrodynamical models predict very different density profiles depending on the evolution of the shells and the different physical processes involved in the wind-wind interaction (e.g. heating and cooling processes). It is therefore not possible to constrain the total shell mass based on the CO observations alone. Additional features outside and inside the shells complicate the interpretation of the data. Complementary observations of, for example, CI as a dissociation product of CO would be necessary to understand the distribution of CO compared to in addition to new detailed hydrodynamical models of the pre-pulse, pulse, and post-pulse wind. Only a comprehensive combination of observations and models will allow us to constrain the evolution of the shells and the changes in the star during the thermal-pulse cycle.

FAUST. XIV. Probing the Flared Disk in L1527 with Sulfur-bearing Molecules

IRAS04368+2557 in L1527 is a Class 0/I protostar with a clear disk-envelope system revealed by previous Atacama Large Millimeter/submillimeter Array (ALMA) observations. In this paper, we discuss the flared structure of this source with observed sulfur-bearing molecules included in the FAUST ALMA large program. The analyses of molecular distributions and kinematics have shown that CS, SO, and OCS trace different regions of the disk-envelope system. To evaluate the temperature across the disk, we derive rotation temperature with the two observed SO lines. The temperature profile shows a clear, flared “butterfly” structure with the higher temperature being ∼50 K and the central lower temperature region (<30 K) coinciding with the continuum peak, suggesting dynamically originated heating rather than radiation heating from the central protostar. Other physical properties, including column densities, are also estimated and further used to demonstrate the vertical structure of the disk-envelope system. The “warped” disk structure of L1527 is confirmed with our analyses, showing that sulfur-bearing molecules are not only effective material probes but also sufficient for structural studies of protostellar systems.

TEMPLATES: Characterization of a Merger in the Dusty Lensing SPT0418–47 System

We present JWST and Atacama Large Millimeter/submillimeter Array (ALMA) imaging for the lensing system SPT0418−47, which includes a strongly lensed, dusty, star-forming galaxy at redshift z = 4.225 and an associated multiply imaged companion. The JWST NIRCam and MIRI imaging observations presented in this paper were acquired as part of the Early Release Science program Targeting Extremely Magnified Panchromatic Lensed Arcs and Their Extended Star formation (TEMPLATES). This data set provides robust mutiwavelength detections of stellar light in both the main (SPT0418A) and companion (SPT0418B) galaxies, while the ALMA detection of [C ii] emission confirms that SPT0418B lies at the same redshift as SPT0418A. We infer that the projected physical separation of the two galaxies is 4.42 ± 0.05 kpc. We derive total magnifications of μ = 29 ± 1 and μ = 4.1 ± 0.7 for SPT0418A and SPT0418B, respectively. We use both prospector and cigale to derive stellar masses. We find that SPT0418A has a stellar mass of M*=3.4−0.6+1.1×1010M⊙ from prospector or M * = 1.5 ± 0.3 × 1010 M ⊙ from cigale. The stellar mass ratio of SPT0418A and SPT0418B is roughly between 4 and 7 ( 4.2−1.6+1.9 for prospector and 7.5 ± 3.7 for cigale). We see evidence of extended structure associated with SPT0418A that is suggestive of a tidal feature. These features, along with the close projected proximity, imply that the system is interacting. Interestingly, the star formation rates and stellar masses of both galaxies are consistent with the main sequence of star-forming galaxies at this epoch, indicating that this ongoing interaction has not noticeably elevated the star formation levels.

Exploring the Mpc Environment of the Quasar ULAS J1342+0928 at z = 7.54

Theoretical models predict that z ≳ 6 quasars are hosted in the most massive halos of the underlying dark matter distribution and thus would be immersed in protoclusters of galaxies. However, observations report inconclusive results. We investigate the 1.1 proper-Mpc2 environment of the z = 7.54 luminous quasar ULAS J1342+0928. We search for Lyman-break galaxy (LBG) candidates using deep imaging from the Hubble Space Telescope (HST) in the Advanced Camera for Surveys (ACS)/F814W, Wide Field Camera 3 (WFC3)/F105W/F125W bands, and Spitzer/Infrared Array Camera at 3.6 and 4.5 μm. We report a zphot=7.69−0.23+0.33 LBG with magF125W = 26.41 at 223 projected proper kpc (pkpc) from the quasar. We find no HST counterpart to one [C ii] emitter previously found with the Atacama Large millimeter/submillimeter Array (ALMA) at 27 projected pkpc and z [C II]=7.5341 ± 0.0009 (Venemans et al. ). We estimate the completeness of our LBG candidates using results from Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey/GOODS deep blank field searches sharing a similar filter setup. We find that >50% of the z ∼ 7.5 Lyman-break galaxies (LBGs) with magF125W > 25.5 are missed due to the absence of a filter redward of the Lyman break in F105W, hindering the UV color accuracy of the candidates. We conduct a QSO-LBG clustering analysis revealing a low LBG excess of 0.46−0.08+1.52 in this quasar field, consistent with an average or low-density field. Consequently, this result does not present strong evidence of an LBG overdensity around ULAS J1342+0928. Furthermore, we identify two LBG candidates with a z phot matching a confirmed z = 6.84 absorber along the line of sight to the quasar. All these galaxy candidates are excellent targets for follow-up observations with JWST and/or ALMA to confirm their redshift and physical properties.

Investigations of Massive Filaments and Star Formation (INFANT). I. Core Identification and Core Mass Function

Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (Investigations of Massive Filaments and Sar Formation) survey, a multiscale, multiwavelength survey of massive filamentary clouds with Atacama Large Millimeter/submillimeter Array (ALMA) band 3/band 6 and Very Large Array K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of eight high-mass star-forming filaments. We covered each target with an approximately rectangular mosaic field of view with two 12 m array configurations, achieving an angular resolution of ∼0.″6 (2700 au at 4.5 kpc) and a continuum rms of ∼0.1 mJy beam−1 (∼0.06 M ⊙ in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15–36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after the formation of protostars. For the high-mass end ( Mcore > 1.5 M ⊙) of the core mass function (CMF) we derive a power-law index of −1.15 ± 0.12 for the whole sample and −1.70 ± 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution (−0.89 ± 0.15 for the young group versus. −1.44 ± 0.25 for the evolved group) and discuss its implications for cluster formation.

Gas-dynamical Mass Measurements of the Supermassive Black Holes in the Early-type Galaxies NGC 4786 and NGC 5193 from ALMA and HST Observations* *Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. These observations are associated

We present molecular gas-dynamical mass measurements of the central black holes in the giant elliptical galaxies NGC 4786 and NGC 5193, based on CO (2−1) observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope near-infrared imaging. The central region in each galaxy contains a circumnuclear disk that exhibits orderly rotation with projected line-of-sight velocities of ∼270 km s−1. We build gas-dynamical models for the rotating disk in each galaxy and fit them directly to the ALMA data cubes. At 0.″31 resolution, the ALMA observations do not fully resolve the black hole sphere of influence (SOI), and neither galaxy exhibits a central rise in rotation speed, indicating that emission from deep within the SOI is not detected. As a result, our models do not tightly constrain the central black hole mass in either galaxy, but they prefer the presence of a central massive object in both galaxies. We measure the black hole mass to be (MBH/108M⊙)=5.0±0.2[1σstatistical]−1.3+1.4[systematic] in NGC 4786 and (MBH/108M⊙)=1.4±0.03[1σstatistical]−0.1+1.5[systematic] in NGC 5193. The largest component of each measurement’s error budget is from the systematic uncertainty associated with the extinction correction in the host galaxy models. This underscores the importance of assessing the impact of dust attenuation on the inferred M BH.

Disk Wind Feedback from High-mass Protostars. III. Synthetic CO Line Emission

Abstract To test theoretical models of massive star formation it is important to compare their predictions with observed systems. To this end, we conduct CO molecular line radiative transfer post-processing of 3D magnetohydrodynamic simulations of various stages in the evolutionary sequence of a massive protostellar core, including its infall envelope and disk wind outflow. Synthetic position–position–velocity cubes of various transitions of 12CO, 13CO, and C18O emission are generated. We also carry out simulated Atacama Large Millimeter/submillimeter Array (ALMA) observations of this emission. We compare the mass, momentum, and kinetic energy estimates obtained from molecular lines to the true values, finding that the mass and momentum estimates can have uncertainties of up to a factor of 4. However, the kinetic energy estimated from molecular lines is more significantly underestimated. Additionally, we compare the mass outflow rate and momentum outflow rate obtained from the synthetic spectra with the true values. Finally, we compare the synthetic spectra with real examples of ALMA-observed protostars and determine the best-fitting protostellar masses and outflow inclination angles. We then calculate the mass outflow rate and momentum outflow rate for these sources, finding that both rates agree with theoretical protostellar evolutionary tracks.

The youngest of hot jupiters in action: episodic accretion outbursts in Gaia20eae

Abstract Recent imaging observations with ALMA and other telescopes found widespread signatures of planet presence in protoplanetary discs at tens of au separations from their host stars. Here we point out that the presence of very massive planets at 0.1 au sized orbits can be deduced for protostars accreting gas at very high rates, when their discs display powerful Thermal Instability bursts. Earlier work showed that a massive planet modifies the nature of this instability, with outbursts triggered at the outer edge of the deep gap opened by the planet. We present simulations of this effect, finding two types of TI outbursts: downstream and upstream of the planet, which may or may not be causally connected. We apply our model to the outburst in Gaia20eae. We find that the agreement between the data and our disc thermal instability model is improved if there is a planet of 6 Jupiter masses orbiting the star at 0.062 au separation. Gaia20eae thus becomes the second episodically erupting star, after FU Ori, where the presence of a massive planet is strongly suspected. Future observations of similar systems will constrain the mode and the frequency of planet formation in such an early epoch.

Small and Large Dust Cavities in Disks around Mid-M Stars in Taurus

High angular resolution imaging by Atacama Large Millimeter/submillimeter Array (ALMA) has revealed the near universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M dwarfs in the Galaxy. Here we present high-resolution (∼50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities, and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ∼10 au for 2M0434, 2M0436, and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared spectral energy distributions to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet–disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ∼10 au cavities under current observations and theories. Currently, disks around mid- to late M stars in Taurus show a higher detection frequency of cavities than earlier-type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass.

Into the thick of it: ALMA 0.45 mm observations of HL Tau at a resolution of 2 au

To comprehend the efficiency of dust evolution within protoplanetary disks, it is crucial to conduct studies of these disks using high-resolution observations at multiple wavelengths with the Atacama Large Millimeter/submillimeter Array (ALMA). In this work, we present high-frequency ALMA observations of the HL Tau disk using its Band 9 centered at a wavelength of 0.45 mm. These observations achieve the highest angular resolution in a protoplanetary disk to date, 12 milliarcseconds (mas), allowing the study of the dust emission at scales of 2 au. We used these data to extend the previously published multiwavelength analysis of the HL Tau disk, constraining the dust temperature, dust surface density, and maximum grain size throughout the disk. We performed this modeling for compact solid dust particles as well as for porous particles. Our new 0.45 mm data mainly trace optically thick emission, providing a tight constraint to the dust temperature profile. We derive maximum particle sizes of $ cm from the inner disk to $ au. Beyond this radius, we find particles between 300 mu m and 1 mm. The total dust mass of the disk is 2.1 J $ with compact grains, and it increases to 6.3 M$_ J $ assuming porous particles. Moreover, an intriguing asymmetry is observed at 32 au in the northeast inner part of the HL Tau disk at 0.45 mm. We propose that this asymmetry is the outcome of a combination of factors, including the optically thick nature of the emission, the orientation of the disk, and a relatively large dust scale height of the grains that is preferentially traced at 0.45 mm. To validate this, we conducted a series of radiative transfer models using the software RADMC-3D. Our models varying dust masses and scale heights successfully replicate the observed asymmetry in the HL Tau disk. If this scenario is correct, our measured dust mass within 32 au would suggest a dust scale height H/R$&gt;0.08$ for the inner disk. Finally, the unprecedented resolution allowed us to probe the dust emission down to scales of a few au for the first time. We observed an increase in brightness temperature inside the estimated water snowline, and we speculate whether this might indicate a traffic-jam effect in the inner disk. Our results show that 0.45 mm observations of protoplanetary disks can be used to robustly constrain the radial profile of their dust temperature. Additionally, the higher optical depths at this wavelength can be used to constrain the vertical scale height of the dust. Finally, these higher frequencies allow us to reach higher spatial resolutions, which have the potential to resolve the region within the water snowline in disks.

Molecular line polarisation from the circumstellar envelopes of asymptotic giant branch stars

Polarisation observations of masers in the circumstellar envelopes (CSEs) around asymptotic giant branch (AGB) stars have revealed strong magnetic fields. However, masers probe only specific lines of sight through the CSE. Non-masing molecular line polarisation observations can more directly reveal the large-scale magnetic field morphology and hence probe the effect of the magnetic field on AGB mass loss and the shaping of the AGB wind. Observations and models of CSE molecular line polarisation can now be used to describe the magnetic field morphology and estimate its strength throughout the entire CSE. We used observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) of molecular line polarisation in the envelope of two AGB stars: CW Leo and R Leo. We modelled the observations using the multi-dimensional polarised radiative transfer tool PORTAL. We found linearly polarised emission, with maximum fractional polarisation on the order of a few percent, in several molecular lines towards both stars. Towards R Leo, we also found a high level of linear polarisation (up to $ for one of the SiO $v=1$ maser transitions. We can explain the observed differences in polarisation structure between the different molecular lines by alignment of the molecules through a combination of the Goldreich-Kylafis effect and radiative alignment effects. We specifically show that the polarisation of CO traces the morphology of the magnetic field. Competition between the alignment mechanisms allowed us to describe the behaviour of the magnetic field strength as a function of the radius throughout the circumstellar envelope of CW Leo. The magnetic field strength derived using this method is inconsistent with the magnetic field strength derived using a structure-function analysis of the CO polarisation and the strength previously derived using CN Zeeman observations. In contrast with CW Leo, the magnetic field in the outer envelope of R Leo appears to be advected outwards by the stellar wind. The ALMA observations and our polarised radiative transfer models show the power of using multiple molecular species to trace the magnetic field behaviour throughout the circumstellar envelope. While the observations appear to confirm the existence of a large-scale magnetic field, further observations and modelling are needed to understand the apparent inconsistency of the magnetic field strength derived with different methods in the envelope of CW Leo.

JWST Observations of Young protoStars (JOYS): Linked accretion and ejection in a Class I protobinary system

Accretion and ejection dictate the outcomes of star and planet formation processes. The mid-infrared (MIR) wavelength range offers key tracers of processes that have been difficult to detect and spatially resolve in protostars until now. We aim to characterize the interplay between accretion and ejection in the low-mass Class I protobinary system TMC1, comprising two young stellar objects: TMC1-W and TMC1-E at a 85 au separation. Using the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) observations in 5 -- 28 mu m range, we measured the intensities of emission lines of H$_2$, atoms, and ions, for instance, the\ Fe II and Ne II and HI recombination lines. We analyzed the spatial distribution of the different species using the MIRI Medium Resolution Spectrometer (MRS) capabilities to spatially resolve emission at 0 scales. We compared these results are with the corresponding Atacama Large Millimeter/submillimeter Array (ALMA) maps tracing cold gas and dust. We detected H$_2$ outflow coming from TMC1-E, with no significant H$_2$ emission from TMC1-W. The H$_2$ emission from TMC1-E outflow appears narrow and extends to wider opening angles with decreasing from S(8) to S(1) rotational transitions, indicating the disk wind as its origin. The outflow from TMC1-E protostar shows spatially extended emission lines of Ne II Ne III Ar II and Ar III with their line ratios consistent with UV radiation as a source of ionization. With ALMA, we detected an accretion streamer infalling from $&gt;$ 1000 au scaling onto the TMC1-E component. The TMC1-W protostar powers a collimated jet, detected with Fe II and Ni II making it consistent with energetic flow. A much weaker ionized jet is observed from TMC1-E, and both jets appear strikingly parallel to each other, indicating that the disks are co-planar. TMC1-W is associated with strong emission from hydrogen recombination lines, tracing the accretion onto the young star. MIRI-MRS observations provide an unprecedented view of protostellar accretion and ejection processes on 20 au scales. Observations of a binary Class I protostellar system show that the two processes are clearly intertwined, with accretion from the envelope onto the disk influencing a wide-angle wind ejected on disk scales. Finally, the accretion from the protostellar disk onto the protostar is associated with the source launching a collimated high-velocity jet within the innermost regions of the disk.