Atacama Large Millimeter/submillimeter Array Research Articles

New constraints on the presence of debris disks around G 196-3 B and VHS J125601.92–125723.9 b

Context. The existence of warm (protoplanetary) disks around very young isolated planetary and brown dwarf mass objects is known based on near- and mid-infrared flux excesses and millimeter observations. These disks may later evolve into debris disks or rings, although none have been observed or confirmed so far. Little is known about circum(sub)stellar and debris disks around substellar objects. Aims. We aim to investigate the presence of debris disks around two of the closest (~20 pc), young substellar companions, namely G196-3 B and VHS J125601.92–125723.9 b (VHS J1256–1257 b), whose masses straddle the borderline between planets and brown dwarfs. Both are companions at wide orbits (≥100 au) of M-type dwarfs and their ages (50–100 Myr and 150–300 Myr, respectively) are thought to be adequate for the detection of second-generation disks. Methods. We obtained deep images of G196-3 B and VHS J1256–1257 b with the NOrthern Extended Millimeter Array (NOEMA) at 1.3 mm. These data were combined with recently published Atacama Large Millimeter Array (ALMA) and Very Large Array (VLA) data of VHS J1256–1257 b at 0.87 mm and 0.9 cm, respectively. Results. Neither G196-3 B nor VHS J1256–1257 b were detected in the NOEMA, ALMA, and VLA data. At 1.3 mm, we imposed flux upper limits of 0.108 mJy (G196-3 B) and 0.153 mJy (VHS J1256–1257 b) with a 3-σ confidence. Using the flux upper limits at the millimeter and radio wavelength regimes, we derived maximum values of 1.38×10−2 MEarth and 5.46 × 10−3 MEarth for the mass of any cold dust that might be surrounding G196-3 B and VHS J1256–1257 b, respectively. Conclusions. We put our results in the context of other deep millimeter observations of free-floating and companion objects with substellar masses smaller than 20 MJup and ages between approximately one and a few hundred million years. Only two very young (2–5.4 Myr) objects are detected out of a few tens of them. This implies that the disks around these very low-mass objects must have small masses, and possibly reduced sizes, in agreement with findings by other groups. If debris disks around substellar objects scale down (in mass and size) in a similar manner as protoplanetary disks do, millimeter observations of moderately young brown dwarfs and planets must be at least two orders of magnitude deeper to be able to detect and characterize their surrounding debris disks.

ATOMIUM: Probing the inner wind of evolved O-rich stars with new, highly excited H2O and OH lines

Context. Water (H2O) and the hydroxyl radical (OH) are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are rotationally or vibrationally highly excited (energies ≳4000 K) have been observed in both H2O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution in the inner wind close to the stellar photosphere with the Atacama Large Millimeter/submillimeter Array (ALMA). Aims. Our goals are: (1) to identify and map the emission and absorption of H2O in several vibrational states, and of OH in Λ-doubling transitions with similar excitation energies; and (2) to determine the physical conditions and kinematics in gas layers close to the extended atmosphere in a sample of asymptotic giant branch stars (AGBs) and red supergiants (RSGs). Methods. Spectra and maps of H2O and OH lines observed in a 27 GHz aggregated bandwidth and with an angular resolution of ~0."02−1."0 were obtained at two epochs with the main ALMA array. Additional observations with the Atacama Compact Array (ACA) were used to check for time variability of water transitions. Radiative transfer models of H2O were revisited to characterize masing conditions. Up-to-date chemical models were used for comparison with the observed OH/H2O abundance ratio. Results. Ten rotational transitions of H2O with excitation energies ~4000–9000 K were observed in vibrational states up to (υ1,υ2,υ3) = (0,1,1). All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split Λ-doubling transitions in υ = 0, J = 27/2 and 29/2 levels of the 2Π3/2 state, as well as J = 33/2 and 35/2 of the 2Π1/2 state of OH with excitation energies of ~4780–8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the υ = 0 and υ = 1 Λ-doubling frequencies have been improved. Our H2O maps show compact emission toward the central star and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the υ2 = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited in AGBs and RSGs. The widespread but weaker 262.898 GHz water line in the υ2 = 1 state also shows signs of maser emission. The OH emission is weak and quasithermally excited. Emission and absorption features of H2O and OH reveal an infall of matter and complex kinematics influenced by binarity. From the OH and H2O column densities derived with nonmasing transitions in a few sources, we obtain OH/H2O abundance ratios of ~(0.7–2.8) × 10−2.

Neutral outflows in high-z QSOs

The OH+(11 − 10) absorption line is a powerful tracer of inflowing and outflowing gas in the predominantly atomic diffuse and turbulent halo surrounding galaxies. In this Letter, we present observations of OH+(11 − 10), CO(9-8) and the underlying dust continuum in five strongly lensed z ∼ 2 − 4 quasi-stellar objects (QSOs), using the Atacama Large Millimeter/submillimeter Array (ALMA) to detect outflowing neutral gas. Blue-shifted OH+(11 − 10) absorption is detected in three out of five QSOs and tentatively detected in a fourth. Absorption at systemic velocities is also detected in one source also displaying blue-shifted absorption. OH+(11 − 10) emission is observed in three out of five QSOs at systemic velocities and the OH+(21 − 10) transition is also detected in one source. CO(9-8) is detected in all five QSOs at high S/N, providing information on the dense molecular gas within the host galaxy. We compare our sample to high-z far-infrared (FIR) luminous star-forming and active galaxies from the literature. We find no difference in OH+ absorption line properties between active and star-forming galaxies with both samples roughly following the same optical depth-dust temperature relation. This suggests that these observables are driven by the same mechanism in both samples. Similarly, star-forming and active galaxies both follow the same OH+ emission–FIR relation. Obscured QSOs display broader (> 800 km s−1) emission than the unobscured QSOs and all but one of the high-z star-forming galaxies (likely caused by the warm molecular gas reservoir obscuring the accreting nucleus). Broader CO(9-8) emission (> 500 km s−1) is found in obscured versus unobscured QSOs, but overall they cover a similar range in line widths as the star-forming galaxies and follow the CO(9-8)–FIR luminosity relation found in low-z galaxies. We find that outflows traced by OH+ are only detected in extreme star-forming galaxies (indicated by broad CO(9-8) emission) and in both types of QSOs, which, in turn, display no red-shifted absorption. This suggests that diffuse neutral outflows in galaxy halos may be associated with the most energetic evolutionary phases leading up to and following the obscured QSO phase.

Investigating the Time Evolution of the Thermal Emission from the Putative Neutron Star in SN 1987A for 50+ Years

Observations collected with the Atacama Large Millimeter/submillimeter Array (ALMA) and analysis of broadband X-ray spectra have recently suggested the presence of a central compact object (CCO) in SN 1987A. However, no direct evidence of the CCO has been found yet. Here we analyze Chandra X-ray observations of SN 1987A collected in 2007 and 2018, and synthesize 2027 Chandra and 2037 Lynx spectra of the faint inner region of SN 1987A. We estimate the temporal evolution of the upper limits of the intrinsic luminosity of the putative CCO in three epochs (2018, 2027, and 2037). We find that these upper limits are higher for higher neutron star (NS) kick velocities due to increased absorption from the surrounding cold ejecta. We compare NS cooling models with both the intrinsic luminosity limits obtained from the X-ray spectra and the ALMA constraints with the assumption that the observed blob of SN 1987A is primarily heated by thermal emission. We find that the synthetic Lynx spectra are crucial to constrain the physical properties of the CCO, which will be confirmed by future observations in the 2040s. We draw our conclusions based on two scenarios, namely the nondetection and detection of the NS by Lynx. If the NS is not detected, its kick velocity should be ≃700 km s−1. Furthermore, nondetection of the NS would suggest rapid cooling processes at the age of 40 yr, implying strong crust superfluidity. Conversely, in the case of NS detection, the mass of the NS envelope must be high.

Challenging interferometric imaging: Machine learning-based source localization from uv-plane observations

Context. Rising interest in radio astronomy and upcoming projects in the field is expected to produce petabytes of data per day, questioning the applicability of traditional radio astronomy data analysis approaches under the new large-scale conditions. This requires new, intelligent, fast, and efficient methods that potentially involve less input from the domain expert. Aims. In our work, we examine, for the first time, the possibility of fast and efficient source localization directly from the uv-observations, omitting the recovering of the dirty or clean images. Methods. We propose a deep neural network-based framework that takes as its input a low-dimensional vector of sampled uv-data and outputs source positions on the sky. We investigated a representation of the complex-valued input uv-data via the real and imaginary and the magnitude and phase components. We provided a comparison of the efficiency of the proposed framework with the traditional source localization pipeline based on the state-of-the-art Python Blob Detection and Source Finder (PyBDSF) method. The investigation was performed on a data set of 9164 sky models simulated using the Common Astronomy Software Applications (CASA) tool for the Atacama Large Millimeter Array (ALMA) Cycle 5.3 antenna configuration. Results. We investigated two scenarios: (i) noise-free as an ideal case and (ii) sky simulations including noise representative of typical extra-galactic millimeter observations. In the noise-free case, the proposed localization framework demonstrates the same high performance as the state-of-the-art PyBDSF method. For noisy data, however, our new method demonstrates significantly better performance, achieving a completeness level that is three times higher for sources with uniform signal-to-noise ratios (S/N) between 1 and 10, and a high increase in completeness in the low S/N regime. Furthermore, the execution time of the proposed framework is significantly reduced (by factors ~30) as compared to traditional methods that include image reconstructions from the uv-plane and subsequent source detections. Conclusions. The proposed framework for obtaining fast and efficient source localization directly from uv-plane observations shows very encouraging results, which could open new horizons for interferometric imaging with existing and future facilities.

Highly structured turbulence in high-mass star formation: An evolved infrared-dark cloud G35.20-0.74 N

Context. Massive stars are generally believed to form in environments characterized by supersonic turbulence. However, recent observations challenge this traditional view. High-spatial- and spectral-resolution observations of the Orion Molecular Cloud (OMC, the closest massive star formation region) and an infrared-dark cloud (IRDC) G35.39 (a typical distant massive star formation region) show a resolution-dependent turbulence, and that high-mass stars are forming exclusively in subsonic to transonic cores in those clouds. These studies demand a re-evaluation of the role of turbulence in massive star formation. Aims. We aim to study the turbulence in a typical massive-star-forming region G35.20-0.74 N (G35.20 in short) with sufficient spatial resolution to resolve the thermal Jeans length, and sufficient spectral resolution to resolve the thermal line width. Methods. We use the Atacama Large Millimeter/submillimeter Array (ALMA) dust continuum emission to resolve fragmentation, the Karl G. Jansky Very Large Array (JVLA) 1.2 cm continuum to trace ionized gas, and JVLA NH3 (1,1) to (7,7) inversion transition lines to trace line width, temperature, and dynamics. We fit those lines and remove line broadening due to channel width, thermal pressure, and velocity gradient to obtain a clean map of intrinsic turbulence. Results. We find that (1) the turbulence in G35.20 is overall supersonic, with mean and median Mach numbers 3.7 and 2.8, respectively. (2) Mach number decreases from 6–7 at a 0.1 pc scale to less than 3 toward the central cores at a 0.01 pc scale. (3) The central ALMA cores appear to be decoupled from the host filament, which is made evident by an opposite velocity gradient and significantly reduced turbulence. Because of intense star-formation activity in G35.20 (as compared to the relatively young and quiescent IRDC G35.39), the supersonic turbulence is likely replenished by protostellar outflows. G35.20 is therefore representative of an evolved form of IRDC G35.39. More observations of a sample of IRDCs are highly demanded in order to further investigate the role of turbulence in the initial conditions required for massive-star formation.

Synergizing 21 cm and submillimetre surveys during reionization: new empirical insights

ABSTRACT We use the latest results from Atacama Large Millimetre/submillimetre Array (ALMA) surveys targeting the ionized carbon [C ii] 158 μm and oxygen [O iii] 88 μm lines, in combination with data-driven predictions for the evolution of neutral hydrogen (H i), to illustrate the prospects for intensity mapping cross-correlations between 21 cm and submillimetre surveys over z ∼ 5 – 7. We work with a data set including the ALPINE and REBELS surveys for [C ii] over z ∼ 4.5 – 7, and ALMA [O iii] detections over z ∼ 6 – 9. The resultant evolution of the [C ii] luminosity–halo mass relation is well described by a double power law at high redshifts, with the best-fitting parameters in good agreement with the results of simulations. The data favour secure detections of the autopower spectrum of [C ii] at all redshifts with an enhanced Fred Young Submillimetre Telescope (FYST)-like configuration. Such an experiment, along with the Murchison Widefield Array (MWA) will be able to measure the 21 cm - [C ii] cross-correlation power with a signal-to-noise ratio of a few tens to a few hundreds. We find that a balloon-borne experiment improving upon the Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) should be able to detect the 21 cm - [O iii] cross-correlation with the MWA and the Square Kilometre Array (SKA)-LOW out to z ∼ 7. Our results have implications for constraining the evolution of luminous sources during the mid-to-end stages of reionization.

Clusters, clouds, and correlations: relating young clusters to giant molecular clouds in M33 and M31

ABSTRACT We use young clusters and giant molecular clouds (GMCs) in the galaxies M33 and M31 to constrain temporal and spatial scales in the star formation process. In M33, we compare the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) catalogue of 1214 clusters with ages measured via colour–magnitude diagram (CMD) fitting to 444 GMCs identified from a new 35 pc resolution Atacama Large Millimeter/submillimeter Array (ALMA) 12CO(2–1) survey. In M31, we compare the Panchromatic Hubble Andromeda Treasury (PHAT) catalogue of 1249 clusters to 251 GMCs measured from a Combined Array for Research in Millimeter-wave Astronomy (CARMA) 12CO(1–0) survey with 20 pc resolution. Through two-point correlation analysis, we find that young clusters have a high probability of being near other young clusters, but correlation between GMCs is suppressed by the cloud identification algorithm. By comparing the positions, we find that younger clusters are closer to GMCs than older clusters. Through cross-correlation analysis of the M33 cluster data, we find that clusters are statistically associated when they are ≤10 Myr old. Utilizing the high precision ages of the clusters, we find that clusters older than ≈18 Myr are uncorrelated with the molecular interstellar medium (ISM). Using the spatial coincidence of the youngest clusters and GMCs in M33, we estimate that clusters spend ≈4–6 Myr inside their parent GMC. Through similar analysis, we find that the GMCs in M33 have a total lifetime of ≈11–15 Myr. We also develop a drift model and show that the above correlations can be explained if the clusters in M33 have a 5–10 km s−1 velocity dispersion relative to the molecular ISM.

The ALMA REBELS Survey: discovery of a massive, highly star-forming, and morphologically complex ULIRG at z = 7.31

ABSTRACT We present Atacama Large Millimeter/Submillimeter Array (ALMA) [C ii] and ∼158 $\rm \mu m$ continuum observations of REBELS-25, a massive, morphologically complex ultra-luminous infrared galaxy (ULIRG; LIR = $1.5^{+0.8}_{-0.5} \times 10^{12}$ L⊙) at z = 7.31, spectroscopically confirmed by the Reionization Era Bright Emission Line Survey (REBELS) ALMA Large Programme. REBELS-25 has a significant stellar mass of $M_{*} = \mbox{$ 8^{+4}_{-2} \times 10^{9} $}{} ~\mbox{M$_\odot $}{}$. From dust-continuum and ultraviolet observations, we determine a total obscured + unobscured star formation rate of SFR $= \mbox{$199^{+101}_{-63}$}{} ~ \mbox{M$_\odot $}~ \mbox{${\rm yr}$}^{-1}$. This is about four times the SFR estimated from an extrapolated main sequence. We also infer a [C ii]-based molecular gas mass of $M_{{\rm H}_{2}} = \mbox{$5.1^{+5.1}_{-2.6} \times 10^{10}$}{} ~\mbox{M$_\odot $}{}$, implying a molecular gas depletion time of $t_{\rm depl, {\rm H}_{2}} = \mbox{$0.3^{+0.3}_{-0.2} $}{}$ Gyr. We observe a [C ii] velocity gradient consistent with disc rotation, but given the current resolution we cannot rule out a more complex velocity structure such as a merger. The spectrum exhibits excess [C ii] emission at large positive velocities (∼500 km s−1), which we interpret as either a merging companion or an outflow. In the outflow scenario, we derive a lower limit of the mass outflow rate of 200 $\mbox{M$_\odot $}~ \mbox{${\rm yr}$}^{-1}$, which is consistent with expectations for a star-formation-driven outflow. Given its large stellar mass, SFR, and molecular gas reservoir ∼700 Myr after the big bang, we explore the future evolution of REBELS-25. Considering a simple, conservative model assuming an exponentially declining star formation history, constant star formation efficiency, and no additional gas inflow, we find that REBELS-25 has the potential to evolve into a galaxy consistent with the properties of high-mass quiescent galaxies recently observed at z ∼ 4.

WISDOM project – XIV. SMBH mass in the early-type galaxies NGC 0612, NGC 1574, and NGC 4261 from CO dynamical modelling

ABSTRACT We present a CO dynamical estimate of the mass of the supermassive black hole (SMBH) in three nearby early-type galaxies: NGC 0612, NGC 1574, and NGC 4261. Our analysis is based on Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 3–6 observations of the 12CO(2–1) emission line with spatial resolutions of 14–58 pc (0.01″–0.26″). We detect disc-like CO distributions on scales from ≲ 200 pc (NGC 1574 and NGC 4261) to ≈10 kpc (NGC 0612). In NGC 0612 and NGC 1574 the bulk of the gas is regularly rotating. The data also provide evidence for the presence of a massive dark object at the centre of NGC 1574, allowing us to obtain the first measure of its mass, MBH = (1.0 ± 0.2) × 108 M⊙ (1σ uncertainty). In NGC 4261, the CO kinematics is clearly dominated by the SMBH gravitational influence, allowing us to determine an accurate black hole mass of (1.62 ± 0.04) × 109 M⊙ (1σ uncertainty). This is fully consistent with a previous CO dynamical estimate obtained using a different modelling technique. Signs of non-circular gas motions (likely outflow) are also identified in the inner regions of NGC 4261. In NGC 0612, we are only able to obtain a (conservative) upper limit of MBH ≲ 3.2 × 109 M⊙. This has likely to be ascribed to the presence of a central CO hole (with a radius much larger than that of the SMBH sphere of influence), combined with the inability of obtaining a robust prediction for the CO velocity curve. The three SMBH mass estimates are overall in agreement with predictions from the MBH − σ* relation.

Comparison of chromospheric diagnostics in a 3D model atmosphere

Context. The Hα line, one of the most studied chromospheric diagnostics, is a tracer of magnetic field structures, while the intensity of its line core provides an estimate of the mass density. The interpretation of Hα observations is complicated by deviations from local thermodynamic equilibrium (LTE) or instantaneous statistical equilibrium conditions. Meanwhile, millimetre (mm) continuum radiation is formed in LTE, and therefore the brightness temperatures from Atacama Large Millimetre-submillimetre Array (ALMA) observations provide a complementary view of the activity and the thermal structure of stellar atmospheres. These two diagnostics together can provide insights into the physical properties of stellar atmospheres, such as their temperature stratification, magnetic structures, and mass density distribution. Aims. In this paper, we present a comparative study between synthetic continuum brightness temperature maps at mm wavelengths (0.3 mm to 8.5 mm) and the width of the Hα 6565 Å line. Methods. We used the 3D radiative-transfer codes Multi3D and Advanced Radiative Transfer (ART) to calculate synthetic spectra for the Hα line and the mm continua, respectively, from an enhanced network atmosphere model with non-equilibrium hydrogen ionisation generated with the state-of-the-art 3D radiation magnetohydrodynamics (rMHD) code Bifrost. We use a Gaussian point spread function (PSF) to simulate the effect of ALMA’s limited spatial resolution and calculate the Hα versus mm continuum correlations and slopes of scatter plots for the original and degraded resolution of the whole box, quiet sun, and enhanced network patches separately. Results. The Hα linewidth and mm brightness temperatures are highly correlated and the correlation is highest at a wavelength of 0.8 mm, that is, in ALMA Band 7. The correlation systematically increases with decreasing resolution. On the other hand, the slopes decrease with increasing wavelength. The degradation of resolution does not have a significant impact on the calculated slopes. Conclusions. With decreasing spatial resolution, the standard deviations of the observables, Hα linewidth, and brightness temperatures decrease and the correlations between them increase, but the slopes do not change significantly. These relations may therefore prove useful in calibrating the mm continuum maps observed with ALMA.

High-angular-resolution NIR view of the Orion Bar revealed by Keck/NIRC2

Context. Nearby photo-dissociation regions (PDRs), where the gas and dust are heated by the far-ultraviolet (FUV) irradiation emitted from stars, are ideal templates with which to study the main stellar feedback processes. Aims. With this study, we aim to probe the detailed structures at the interfaces between ionized, atomic, and molecular gas in the Orion Bar. This nearby prototypical strongly irradiated PDR are among the first targets of the James Webb Space Telescope (JWST) within the framework of the PDRs4All Early Release Science program. Methods. We employed the subarcsecond resolution accessible with Keck-II NIRC2 and its adaptive optics system to obtain images of the vibrationally excited line H2 1−0 S(1) at 2.12 µm that are more detailed and complete than ever before. H2 1−0 S(1) traces the dissociation front (DF), and the [FeII] and Brγ lines, at 1.64 and 2.16 µm, respectively, trace the ionization front (IF). The former is a powerful tracer of the FUV radiation field strength and gas density distribution at the PDR edge, while the last two trace the temperature and density distribution from the ionized gas to the PDR. We obtained narrow-band filter images in these key gas line diagnostics over ~40″ at spatial scales of ~0.1″ (~0.0002 pc or ~40 AU at 414 pc). Results. The Keck/Near Infrared Camera 2 (NIRC2) observations spatially resolve a plethora of irradiated substructures such as ridges, filaments, globules, and proplyds. This portends what JWST should accomplish and how it will complement the highest resolution Atacama Large Millimeter/submillimeter Array (ALMA) maps of the molecular cloud. We observe a remarkable spatial coincidence between the H2 1−0 S(1) vibrational and HCO+ J = 4−3 rotational emission previously obtained with ALMA. This likely indicates the intimate link between these two molecular species and highlights that in high-pressure PDRs, the H/H2 and C+/C/CO transitions zones come closer than in a typical layered structure of a constant density PDR. The H/H2 dissociation front appears as a highly structured region containing substructures with a typical thickness of a few ~10−3 pc.

Another X-ray UFO without a momentum-boosted molecular outflow

We present Atacama Large Millimetre/submillimetre Array (ALMA) CO(1–0) observations of the nearby infrared luminous (LIRG) galaxy pair IRAS 05054+1718 (also known as CGCG 468-002), as well as a new analysis of X-ray data of this source collected between 2012 and 2021 using the Nuclear Spectroscopic Telescope Array (NuSTAR), Swift, and the XMM-Newton satellites. The western component of the pair, NED01, hosts a Seyfert 1.9 nucleus that is responsible for launching a powerful X-ray ultra-fast outflow (UFO). Our X-ray spectral analysis suggests that the UFO could be variable or multi-component in velocity, ranging from v/c ∼ −0.12 (as seen in Swift) to v/c ∼ −0.23 (as seen in NuSTAR), and constrains its momentum flux to be ṗoutX−ray ∼ (4 ± 2) × 1034 g cm s−2. The ALMA CO(1–0) observations, obtained with an angular resolution of 2.2″, although targeting mainly NED01, also include the eastern component of the pair, NED02, a less-studied LIRG with no clear evidence of an active galactic nucleus (AGN). We study the CO(1–0) kinematics in the two galaxies using the 3D-BAROLO code. In both sources we can model the bulk of the CO(1–0) emission with rotating disks and, after subtracting the best-fit models, we detect compact residual emission at S/N = 15 within ∼3 kpc of the centre. A molecular outflow in NED01, if present, cannot be brighter than such residuals, implying an upper limit on its outflow rate of Ṁoutmol ≲ 19 ± 14 M⊙ yr−1 and on its momentum rate of ṗoutmol ≲ (2.7 ± 2.4) × 1034 g cm s−1. Combined with the revised energetics of the X-ray wind, we derive an upper limit on the momentum rate ratio of ṗoutmol/ṗoutX−ray < 0.67. We discuss these results in the context of the expectations of AGN feedback models, and we propose that the X-ray disk wind in NED01 has not significantly impacted the molecular gas reservoir (yet), and we can constrain its effect to be much smaller than expectations of AGN ‘energy-driven’ feedback models. We also consider and discuss the hypothesis of asymmetries of the molecular disk not properly captured by the 3D-BAROLO code. Our results highlight the challenges in testing the predictions of popular AGN disk-wind feedback theories, even in the presence of good-quality multi-wavelength observations.

A major asymmetric ice trap in a planet-forming disk

Context. Most well-resolved disks observed with the Atacama Large Millimeter/submillimeter Array (ALMA) show signs of dust traps. These dust traps set the chemical composition of the planet-forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust at smaller radii of volatiles. Aims. In this work, we analyse the first detection of nitric oxide (NO) in a protoplanetary disk. We aim to constrain the nitrogen chemistry and the gas-phase C/O ratio in the highly asymmetric dust trap in the Oph-IRS 48 disk. Methods. We used ALMA observations of NO, CN, C2H, and related molecules in the Oph-IRS 48 disk. We modeled the effect of the increased dust-to-gas ratio in the dust trap on the physical and chemical structure using a dedicated nitrogen chemistry network in the thermochemical code DALI. Furthermore, we explored how ice sublimation contributes to the observed emission lines. Finally, we used the model to put constraints on the nitrogen-bearing ices. Results. Nitric oxide (NO) is only observed at the location of the dust trap, but CN and C2H are not detected in the Oph-IRS 48 disk. This results in an CN/NO column density ratio of <0.05 and thus a low C/O ratio at the location of the dust trap. Models show that the dust trap cools the disk midplane down to ~30 K, just above the NO sublimation temperature of ~25 K. The main gas-phase formation pathways to NO though OH and NH in the fiducial model predict NO emission that is an order of magnitude lower than what has been observed. The gaseous NO column density can be increased by factors ranging from 2.8 to 10 when the H2O and NH3 gas abundances are significantly boosted by ice sublimation. However, these models are inconsistent with the upper limits on the H2O and OH column densities derived from Herschel PACS observations and the upper limit on CN derived from ALMA observations. As the models require an additional source of NO to explain its detection, the NO seen in the observations is likely the photodissociation product of a larger molecule sublimating from the ices. The non-detection of CN provides a tighter constraint on the disk C/O ratio than the C2H upper limit. Conclusions. We propose that the NO emission in the Oph-IRS 48 disk is closely related to the nitrogen-bearing ices sublimating in the dust trap. The non-detection of CN constrains the C/O ratio both inside and outside the dust trap to be <1 if all nitrogen initially starts as N2 and ≤ 0.6, consistent with the Solar value, if (at least part of) the nitrogen initially starts as N or NH3.

Detection of complex nitrogen-bearing molecule ethyl cyanide towards the hot molecular core G10.47 + 0.03

The studies of the complex organic molecular lines towards the hot molecular cores at millimeter and submillimeter wavelengths provide instructive knowledge about the chemical complexity in the interstellar medium (ISM). We present the detection of the rotational emission lines of the complex nitrogen-bearing molecule ethyl cyanide (C2H5CN) towards the chemically rich hot molecular core G10.47 + 0.03 using the Atacama Large Millimeter/Submillimeter Array (ALMA) band 4 observations. The estimated column density of C2H5CN towards the G10.47 + 0.03 is (7.7 ± 0.5) × $10^{16}\text{ cm}^{-2}$ with the high rotational temperature of 352.9 ± 66.8 K. The estimated fractional abundance of C2H5CN with respect to H2 towards the G10.47 + 0.03 is 5.70 × $10^{-9}$ . We observe that the estimated fractional abundance of C2H5CN is similar to the existing three-phase warm-up chemical modelling abundance of C2H5CN. We also discuss the possible formation mechanism of C2H5CN towards the hot molecular cores, and we claim the barrierless and exothermic radical-radical reaction between CH2 and CH2CN is responsible for the production of low abundance of C2H5CN ( ${\sim} 10^{-9}$ ) in the grain surface of G10.47 + 0.03.

JWST and ALMA Imaging of Dust-obscured, Massive Substructures in a Typical z ∼ 3 Star-forming Disk Galaxy

We present an identification of dust-attenuated star-forming galactic-disk substructures in a typical star-forming galaxy (SFG), UDF2, at z = 2.696. To date, substructures containing significant buildup of stellar mass and actively forming stars have yet to be found in typical (i.e., main-sequence) SFGs at z > 2. This is due to the strong dust attenuation common in massive galaxies at the epoch and the scarcity of high-resolution, high-sensitivity extinction-independent imaging. To search for disk substructures, we subtracted the central stellar-mass disk from the JWST/NIRCam rest-frame 1.2 μm image (0.″13 resolution) and subtracted, in the visibility plane, the central starburst disk from Atacama Large Millimeter/submillimeter Array (ALMA) rest-frame 240 μm observations (0.″03 resolution). The residual images revealed substructures at rest-frame 1.2 μm colocated with those found at rest-frame 240 μm, ≃2 kpc away from the galactic center. The largest substructure contains ≃20% of the total stellar mass and ≃5% of the total star formation rate of the galaxy. While UDF2 exhibits a kinematically ordered velocity field of molecular gas consistent with a secularly evolving disk, more sensitive observations are required to characterize the nature and the origin of this substructure (spiral arms, minor merger, or other types of disk instabilities). UDF2 resides in an overdense region (N ≥ 4 massive galaxies within 70 kpc projected distance at z = 2.690–2.697) and the substructures may be associated with interaction-induced instabilities. Importantly, a statistical sample of such substructures identified with JWST and ALMA could play a key role in bridging the gap between the bulge-forming starburst and the rest of the galaxy.

Evolution of the Outflow in the Water Fountain Source IRAS 18043–2116* *Released on 2021 March 1.

We present the spectral and spatial evolution of H2O masers associated with the water fountain source IRAS 18043−2116, found in observations with the Nobeyama 45 m Telescope and the Australia Telescope Compact Array. We have found new highest-velocity components of the H2O masers (at the redshifted side V LSR ≃ 376 km s−1 and at the blueshifted side V LSR ≃ −165 km s−1), and the resulting velocity spread of ≃540 km s−1 breaks the speed record of fast jets/outflows in this type of sources. The locations of those components have offsets from the axis joining the two major maser clusters, indicating a large opening angle of the outflow (∼60°). The evolution of the maser cluster separation of ∼2.9 mas yr−1 and the compact (∼0.″2) CO emission source mapped with the Atacama Large Millimeter-submillimeter Array suggest a very short (∼30 yr) timescale of the outflow. We also confirmed an increase in the flux density of the 22 GHz continuum source. The properties of the jet and the continuum sources and their possible evolution in the transition to the planetary nebula phase are further discussed.

The ALMA REBELS Survey: the first infrared luminosity function measurement at z ∼ 7

ABSTRACT We present the first observational infrared luminosity function (IRLF) measurement in the Epoch of Reionization (EoR) based on a ultraviolet (UV)-selected galaxy sample with the Atacama Large Millimeter Array (ALMA) spectroscopic observations. Our analysis is based on the ALMA large program Reionization Era Bright Emission Line Survey (REBELS), which targets 42 galaxies at z = 6.4–7.7 with [C ii] 158 $\rm{\mu m}$ line scans. 16 sources exhibit dust detection, 15 of which are also spectroscopically confirmed through the [C ii] line. The infrared (IR) luminosities of the sample range from log LIR/L⊙ = 11.4 to 12.2. Using the UV luminosity function as a proxy to derive the effective volume for each of our target sources, we derive IRLF estimates, both for detections and for the full sample including IR luminosity upper limits. The resulting IRLFs are well reproduced by a Schechter function with the characteristic luminosity of $\log L_{*}/\mathrm{ L}_\odot =11.6^{+0.2}_{-0.1}$ . Our observational results are in broad agreement with the average of predicted IRLFs from simulations at z ∼ 7. Conversely, our IRLFs lie significantly below lower redshift estimates, suggesting a rapid evolution from z ∼ 4 to z ∼ 7, into the reionization epoch. The IR obscured contribution to the cosmic star formation rate density at z ∼ 7 amounts to $\mathrm{log(SFRD/{\rm M}_{\odot }\,yr^{-1}\,Mpc^{-3}) = -2.66^{+0.17}_{-0.14} }$ that is at least ∼10 per cent of UV-based estimates. We conclude that the presence of dust is already abundant in the EoR and discuss the possibility of unveiling larger samples of dusty galaxies with future ALMA and JWST observations.

The HH 24 Complex: Jets, Multiple Star Formation, and Orphaned Protostars

The HH 24 complex harbors five collimated jets emanating from a small protostellar multiple system. We have carried out a multiwavelength study of the jets, their driving sources, and the cloud core hosting the embedded stellar system, based on data from the Hubble Space Telescope, Gemini, Subaru, Apache Point Observatory 3.5 m, Karl G. Jansky Very Large Array, and Atacama Large Millimeter/submillimeter Array (ALMA) telescopes. The data show that the multiple system, SSV 63, contains at least 7 sources, ranging in mass from the hydrogen-burning limit to proto-Herbig Ae stars. The stars are in an unstable nonhierarchical configuration, and one member, a borderline brown dwarf, is moving away from the protostellar system with 25 km s−1, after being ejected ∼5800 yr ago as an orphaned protostar. Five of the embedded sources are surrounded by small, possibly truncated, disks resolved at 1.3 mm with ALMA. Proper motions and radial velocities imply jet speeds of 200–300 km s−1. The two main HH 24 jets, E and C, form a bipolar jet system that traces the innermost portions of parsec-scale chains of Herbig–Haro and H2 shocks with a total extent of at least 3 pc. H2CO and C18O observations show that the core has been churned and continuously fed by an infalling streamer. 13CO and 12CO trace compact, low-velocity, cavity walls carved by the jets and an ultracompact molecular outflow from the most embedded object. Chaotic N-body dynamics likely will eject several more of these objects. The ejection of stars from their feeding zones sets their masses. Dynamical decay of nonhierarchical systems can thus be a major contributor to establishing the initial mass function.

FAUST – VIII. The protostellar disc of VLA 1623–2417W and its streamers imaged by ALMA

ABSTRACT More than 50 per cent of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C18O (2–1) observations of the VLA 1623–2417 protostellar system at 50 au angular resolution as part of the ALMA (Atacama Large Millimeter/submillimeter Array) Large Program FAUST. The 1.3 mm continuum probes the discs of VLA 1623A, B, and W, and the circumbinary disc of the A1 + A2 binary. The C18O emission reveals, for the first time, the gas in the disc envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, Mdyn = 0.45 ± 0.08 M⊙, and the mass of its disc, Mdisc ∼ 6 × 10−3 M⊙. C18O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1 + A2 + B system with the disc of VLA 1623W. The presence of the streamers, as well as the spatial (∼1300 au) and velocity (∼2.2 km s−1) offsets of VLA 1623W, suggests that either sources W and A + B formed in different cores, interacting between them, or source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB on to VLA 1623W, thus concurring to set its final mass and chemical content.