Search

The nature of the gas in CO-rich debris discs remains poorly understood, as it could either be a remnant from the earlier Class II phase or of secondary origin, driven by the destruction of icy planetesimals. The aim of this paper was to elucidate the origin of the gas content in the debris discs via various simple molecules that are often detected in the less-evolved Class II discs. We present millimetre molecular line observations of nine circumstellar discs around A-type stars: four CO-rich debris discs (HD 21997, HD 121617, HD 131488, HD 131835) and five old Herbig Ae protoplanetary discs (HD 139614, HD 141569, HD 142666, HD 145718, HD 100453). The sources were observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in Bands 5 and 6 with 1-2'' resolution. The Herbig Ae discs are detected in the CO isotopologues, CN, HCN, HCO$^+$, C2H, and CS lines. In contrast, only CO isotopologues are detected in the debris discs, showing a similar amount of CO to that found in the Herbig Ae protoplanetary discs. Using chemical and radiative transfer modelling, we show that the abundances of molecules other than CO in debris discs are expected to be very low. We consider multiple sets of initial elemental abundances with various degrees of H2 depletion. We find that the HCO+ lines should be the second brightest after the CO lines and that their intensities strongly depend on the overall CO/H$_2$ ratio of the gas. However, even in the ISM-like scenario, the simulated HCO+ emission remains weak as required by our non-detections., 15 pages, 9 figures, 7 tables Accepted for publication in MNRAS

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of $\mathrm{^{13}CO(J=1-0)}$ line and 104 GHz continuum emission from NGC 604, a giant HII region (GHR) in the nearby spiral galaxy M33. Our high spatial resolution images ( 3.2"$\times$ 2.4", corresponding to $13 \times 10$ pc physical scale) allow us to detect fifteen molecular clouds. We find spatial offsets between the $^{13}CO$ and 104 GHz continuum emission and also detect continuum emission near the centre of the GHR. The identified molecular clouds have sizes ranging from 5-21 pc, linewidths of 0.3-3.0 $\mathrm{kms^{-1}}$ and luminosity-derived masses of (0.4-80.5) $\times 10^3$ M$_{\bigodot}$. These molecular clouds are in near virial equilibrium, with a spearman correlation coefficient of 0.98. The linewidth-size relationship for these clouds is offset from the corresponding relations for the Milky Way and for NGC 300, although this may be an artefact of the dendrogram process., The article has 9 pages, 7 figures and 1 table

The interstellar traveler, 2I/Borisov, is the first clearly active extrasolar comet, ever detected in our Solar system. We obtained high-resolution interferometric observations of 2I/Borisov with the Atacama Large Millimeter/submillimeter Array (ALMA), and multi-color optical observations with the Very Large Telescope (VLT) to gain a comprehensive understanding of the dust properties of this comet. We found that the dust coma of 2I/Borisov consists of compact "pebbles" of radii exceeding ~1 mm, suggesting that the dust particles have experienced compaction through mutual impacts during the bouncing collision phase in the protoplanetary disk. We derived a dust mass loss rate of >= 200 kg/s and a dust-to-gas ratio >=3. Our long term monitoring of 2I/Borisov with VLT indicates a steady dust mass loss with no significant dust fragmentation and/or sublimation occurring in the coma. We also detected emissions from carbon monoxide gas (CO) with ALMA and derived the gas production rate of Q(CO) (3.3+/-0.8)x10^{26} mole/s. We found that the CO/H$_2$O mixing ratio of 2I/Borisov changed drastically before and after perihelion, indicating the heterogeneity of the cometary nucleus, with components formed at different locations beyond the volatile snow-line with different chemical abundances. Our observations suggest that 2I/Borisov's home system, much like our own system, experienced efficient radial mixing from the innermost parts of its protoplanetary disk to beyond the frost line of CO., 12 pages, 4 figures

We investigated the radio spectra of two magnetars, PSR J1622$-$4950 and 1E 1547.0$-$5408, using observations from the Australia Telescope Compact Array and the Atacama Large Millimeter/submillimeter Array taken in 2017. Our observations of PSR J1622$-$4950 show a steep spectrum with a spectral index of $-$1.3 $\pm$ 0.2 in the range of 5.5-45 GHz during its re-activating X-ray outburst in 2017. By comparing the data taken at different epochs, we found significant enhancement in the radio flux density. The spectrum of 1E 1547.0$-$5408 was inverted in the range of 43-95 GHz, suggesting a spectral peak at a few hundred gigahertz. Moreover, we obtained the X-ray and radio data of radio magnetars, PSR J1622$-$4950 and SGR J1745$-$2900, from literature and found two interesting properties. First, radio emission is known to be associated with X-ray outburst but has different evolution. We further found that the rising time of the radio emission is much longer than that of the X-ray during the outburst. Second, the radio magnetars may have double peak spectra at a few GHz and a few hundred GHz. This could indicate that the emission mechanism is different in the cm and the sub-mm bands. These two phenomenons could provide a hint to understand the origin of radio emission and its connection with the X-ray properties., 7 pages, 4 figures, to be published in MNRAS

Supermassive black hole (SMBH) masses can be measured by resolving the dynamical influences of the SMBHs on tracers of the central potentials. Modern long-baseline interferometers have enabled the use of molecular gas as such a tracer. We present here Atacama Large Millimeter/submillimeter Array observations of the elliptical galaxy NGC 7052 at 0${^{\prime\prime}_{.}}$11 ($37\,$pc) resolution in the 12CO(2-1) line and $1.3\,$ mm continuum emission. This resolution is sufficient to resolve the region in which the potential is dominated by the SMBH. We forward model these observations, using a multi-Gaussian expansion of a Hubble Space Telescope F814W image and a spatially constant mass-to-light ratio to model the stellar mass distribution. We infer an SMBH mass of $2.5\pm 0.3\times 10^{9}\, \mathrm{M_\odot }$ and a stellar I-band mass-to-light ratio of $4.6\pm 0.2\, \mathrm{M_\odot /L_{\odot ,I}}$ (3σ confidence intervals). This SMBH mass is significantly larger than that derived using ionized gas kinematics, which however appears significantly more kinematically disturbed than the molecular gas. We also show that a central molecular gas deficit is likely to be the result of tidal disruption of molecular gas clouds due to the strong gradient in the central gravitational potential.

We present the highest resolution CO (2-1) observations obtained to date (0.25") of NGC 3256 and use them to determine the detailed properties of the molecular interstellar medium in the central 6 kpc of this merger. Distributions of physical quantities are reported from pixel-by-pixel measurements at 55 and 120 pc scales and compared to disc galaxies observed by PHANGS-ALMA. Mass surface densities range from 8 to 5500 M$_{\odot}$ pc$^{-2}$ and velocity dispersions from 10 to 200 km s$^{-1}$. Peak brightness temperatures as large as 37 K are measured, indicating the gas in NGC 3256 may be hotter than all regions in nearby disc galaxies measured by PHANGS-ALMA. Brightness temperatures even surpass those in the overlap region of NGC 4038/9 at the same scales. The majority of the gas appears unbound with median virial parameters of 7 to 19, although external pressure may bind some of the gas. High internal turbulent pressures of 10$^{5}$ to 10$^{10}$ K cm$^{-3}$ are found. Given the lack of significant trends in surface density, brightness temperature, and velocity dispersion with physical scale we argue the molecular gas is made up of a smooth medium down to 55 pc scales, unlike the more structured medium found in the PHANGS-ALMA disc galaxies., accepted to MNRAS, 21 pages, 8 figures

Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe., 19 pages main journal, 97 pages Supplementary Information

Numerical simulations have predicted that substructures such as spiral arms can be produced through a gravitationally unstable disk around high-mass young stellar objects (HMYSOs). Recent high-resolution observations from the Atacama Large Millimeter/submillimeter Array have investigated these substructures at a spatial resolution of ∼100 au. An accretion burst, which is a manifestation of an increase in the accretion rate caused by a gravitational instability in the disk, can result in luminosity outbursting phenomena - as has been seen in several HMYSOs. However, no clear relationship between the accretion bursts and disk substructures has been established. Here we report the detections of three new molecular maser species, HDO, HNCO and 13CH3OH, from the direction of the HMYSO G358.93-0.03 during a 6.7 GHz methanol maser flaring event. High-quality imaging of the three new maser species exhibits consistent observational evidence that these masers closely trace the spiral-arm substructures around this HMYSO. The rapid decay of the spectral lines emitted from these molecules suggests that these are transient phenomena (for only ∼1 month), probably associated with rapid changes in radiation field due to an accretion burst. Therefore, these new maser species provide evidence linking the spiral-arm substructure with an accretion burst, both expected from massive disk instabilities.

We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring., 32 pages, 19 figures, 6 tables, accepted for publication in MNRAS

We report on the detection of the [CII] 157.7 $\mu$m emission from the Lyman break galaxy (LBG) MACS0416_Y1 at z = 8.3113, by using the Atacama Large Millimeter/submillimeter Array (ALMA). The luminosity ratio of [OIII] 88 $\mu$m (from previous campaigns) to [CII] is 9.31 $\pm$ 2.6, indicative of hard interstellar radiation fields and/or a low covering fraction of photo-dissociation regions. The emission of [CII] is cospatial to the 850 $\mu$m dust emission (90 $\mu$m rest-frame, from previous campaigns), however the peak [CII] emission does not agree with the peak [OIII] emission, suggesting that the lines originate from different conditions in the interstellar medium. We fail to detect continuum emission at 1.5 mm (160 $\mu$m rest-frame) down to 18 $\mu$Jy (3$\sigma$). This nondetection places a strong limit on the dust spectrum, considering the 137 $\pm$ 26 $\mu$Jy continuum emission at 850 $\mu$m. This suggests an unusually warm dust component (T $>$ 80 K, 90% confidence limit), and/or a steep dust-emissivity index ($\beta_{\rm dust}$ $>$ 2), compared to galaxy-wide dust emission found at lower redshifts (typically T $\sim$ 30 - 50 K, $\beta_{\rm dust}$ $\sim$ 1 - 2). If such temperatures are common, this would reduce the required dust mass and relax the dust production problem at the highest redshifts. We therefore warn against the use of only single-wavelength information to derive physical properties, recommend a more thorough examination of dust temperatures in the early Universe, and stress the need for instrumentation that probes the peak of warm dust in the Epoch of Reionization., Comment: 15 pages, 8 figures, Accepted by MNRAS

The gas content of galaxies is a key factor for their growth, starting from star formation and black hole accretion to galaxy mergers. Thus, characterising its properties via observations of tracers like the CO emission line is of big importance in order to understand the bigger picture of galaxy evolution. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum, CO(5-4) and CO(8-7) line emission in the quasar--star-forming companion system SMM J04135+10277 (z=2.84). Earlier low-$J$ CO studies of this system found a huge molecular gas reservoir associated to the companion galaxy, while the quasar appeared gas-poor. Our CO observations revealed that the host galaxy of the quasar is also gas-rich, with an estimated molecular gas mass of $\sim(0.7-2.3)\times10^{10}$ M$_{\odot}$. The CO line profiles of the companion galaxy are broad ($\sim1000$ km s$^{-1}$), and show signs of rotation of a compact, massive system. In contrast to previous far-infrared observations, we resolve the continuum emission and detect both sources, with the companion galaxy dominating the dust continuum and the quasar having a $\sim25\%$ contribution to the total dust emission. By fitting the infrared spectral energy distribution of the sources with \textsc{MR-MOOSE} and empirical templates, the infrared luminosities of the quasar and the companion are in the range of $L_{\rm IR, QSO}\sim(2.1-9.6)\times10^{12}$ L$_{\odot}$ and $L_{\rm IR, Comp.}\sim(2.4-24)\times10^{12}$ L$_{\odot}$, while the estimated star formation rates are $\sim210-960$ M$_{\odot}$ yr$^{-1}$ and $\sim240-2400$ M$_{\odot}$ yr$^{-1}$, respectively. Our results demonstrate that non-detection of low-$J$ CO transition lines in similar sources does not necessarily imply the absence of massive molecular gas reservoir but that the excitation conditions favour the excitation of high-$J$ transitions., Comment: Accepted by MNRAS; 14 pages, 9 figures

We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope. We established a responsivity model that converts frequency responses of the MKIDs to line-of-sight brightness temperature. We estimated two parameters of the responsivity model using a set of skydip data taken under various precipitable water vapor (PWV, 0.4-3.0 mm) conditions for each MKID. The line-of-sight brightness temperature of sky is estimated using an atmospheric transmission model and the PWVs. As a result, we obtain an average temperature calibration uncertainty of $1\sigma=4$%, which is smaller than other photometric biases. In addition, the average forward efficiency of 0.88 in our responsivity model is consistent with the value expected from the geometrical support structure of the telescope. We also estimate line-of-sight PWVs of each skydip observation using the frequency response of MKIDs, and confirm the consistency with PWVs reported by the Atacama Large Millimeter/submillimeter Array., Comment: Accepted for publication in the Journal of Low Temperature Physics for LTD-18 special issue

The ever increasing challenges posed by the science projects in astronomy have skyrocketed the complexity of the new generation telescopes. Due to the climate and sky requirements, these high-precision instruments are generally located in remote areas, suffering from the harsh environments around it. These modern telescopes not only produce massive amounts of scientific data, but they also generate an enormous amount of operational information. The Atacama Large Millimeter/submillimeter Array (ALMA) is one of these unique instruments, generating more than 50 Gb of operational data every day while functioning in conditions of extreme dryness and altitude. To maintain the array working under extreme conditions, the engineering teams must check over 130,000 monitoring points, combing through the massive datasets produced every day. To make this possible, predictive tools are needed to identify, hopefully beforehand, the occurrence of failures in all the different subsystems. This work presents a novel fault detection scheme for one of these subsystems, the Intermediate Frequency Processors (IFP). This subsystem is critical to process the information gathered by each antenna and communicate it, reliably, to the correlator for processing. Our approach is based on echo state networks, a configuration of artificial neural networks, used to learn and predict the signal patterns. These patterns are later compared to the actual signal, to identify failure modes. Additional preprocessing techniques were also added since the signal-to-noise ratio of the data used was very low. The proposed scheme was tested in over seven years of data from 132 IFPs at ALMA, showing an accuracy of over 70%. Furthermore, the detection was done several months earlier, on average, when compared to what human operators did. These results help the maintenance procedures, increasing reliability while reducing humans’ exposure to the harsh environment where the antennas are. Although applied to a specific fault, this technique is broad enough to be applied to other types of faults and settings.

The distributions of temperature and density in the sunspot atmosphere are crucial for an understanding of the process of energy transfer from the photosphere upward to the corona in strong magnetic fields. The joint analysis of millimeter observations at ALMA (Atacama Large Millimeter/Submillimeter Array) and observations performed in the microwave range (centimeter wavelengths) with the RATAN-600 radio telescope of the Russian Academy of Sciences provides new data on the temperature distribution and physical processes at different altitudes above a spot. The NOAA 12470 active region was observed and mapped at 1.3 and 3 mm (ALMA) and 2–10 cm (RATAN-600) in December 2015. These observations are analyzed, and the results are compared to models of sunspot atmospheres. The fundamental problems arising in the study of atmospheres of active regions in the millimeter and centimeter ranges are determined, and the importance of the 3–18 mm range in the physics of the generation and transfer of energy for corona heating is demonstrated.

Traditionally, corrugated conical horn antennas have been the main choice for use in astronomical receivers in the range of millimeter and submillimeter waves. They present low cross-polar level and high coupling efficiency into the fundamental Gaussian mode. However, this type of antenna is difficult to manufacture, inevitably increasing its price and extending the production process. In this article, we present two kinds of feed horn antennas, aimed for use in a frequency range equivalent to atacama large millimeter/submillimeter array (ALMA) Band 6 (211–275 GHz), which can be fabricated in a much simpler way with the conventional machining tools. Specifically, we present the design and performance comparison of smooth-walled spline-profile horns in two geometries, diagonal, and conical. Optimization of the designs has been made by means of an algorithm that allowed us to obtain models whose electrical and mechanical characteristics make them competitive when compared with corrugated horns. In particular, they are $\text{40}\%$ shorter than the conventional corrugated horns suited for this band, representing an advantage given the stringent space constraints of most astronomical receivers. We also demonstrate that they can be coupled efficiently to an astronomical-grade optical system, using ALMA Band-6 receiver as an example. Furthermore, we have constructed the diagonal horn and characterized it thoroughly. Experimental results of the radiation pattern at room temperature show a good cross-polar performance with levels below −20 dB and Gaussicity above $\text{96}\%$ . Our calculations show a good antenna-efficiency level with losses less than $\text{1}\%$ . All these properties demonstrate the feasibility of this type of horns to become the main option at the time of choosing a feed system for cutting-edge astronomical applications.

The physical processes behind the transfer of mass from parsec-scale clumps to massive-star-forming cores remain elusive. We investigate the relation between the clump morphology and the mass fraction that ends up in its most massive core (MMC) as a function of infrared brightness, i.e. a clump evolutionary tracer. Using ALMA 12 m and ACA we surveyed 6 infrared-dark hubs in 2.9mm continuum at $\sim$3" resolution. To put our sample into context, we also re-analysed published ALMA data from a sample of 29 high mass-surface density ATLASGAL sources. We characterise the size, mass, morphology, and infrared brightness of the clumps using Herschel and Spitzer data. Within the 6 newly observed hubs, we identify 67 cores, and find that the MMCs have masses between 15-911 $\mathrm{M}_{\odot}$ within a radius of 0.018-0.156 pc. The MMC of each hub contains 3-24% of the clump mass ($f_\mathrm{MMC}$), becoming 5-36% once core masses are normalised to the median core radius. Across the 35 clumps, we find no significant difference in the median $f_\mathrm{MMC}$ values of hub and non-hub systems, likely the consequence of a sample bias. However, we find that $f_\mathrm{MMC}$ is $\sim$7.9 times larger for infrared-dark clumps compared to infrared-bright ones. This factor increases up to $\sim$14.5 when comparing our sample of 6 infrared-dark hubs to infrared-bright clumps. We speculate that hub-filament systems efficiently concentrate mass within their MMC early on during its evolution. As clumps evolve, they grow in mass, but such growth does not lead to the formation of more massive MMCs., 15 pages, 9 figures

We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster., 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) and Band 3 (3 mm) at spatial resolutions 0.5 - 3" for 16 protostellar sources. The protostellar envelope is well traced by C$^{18}$O, DCO$^+$ and N$_2$D$^+$, with the freeze-out of CO governing the chemistry at envelope scales. Molecular outflows are seen in classical shock tracers like SiO and SO, but ice-mantle products such as CH$_3$OH and HNCO released with the shock are also observed. The molecular jet is prominent not only in SiO and SO but also occasionally in H$_2$CO. The cavity walls show tracers of UV-irradiation such as C$_2$H c-C$_3$H$_2$ and CN. The hot inner envelope, apart from showing emission from complex organic molecules (COMs), also presents compact emission from small molecules like H$_2$S, SO, OCS and H$^{13}$CN, most likely related to ice sublimation and high-temperature chemistry. Sub-arcsecond millimeter-wave observations allow to identify those (simple) molecules that best trace each of the physical components of a protostellar system. COMs are found both in the hot inner envelope (high excitation lines) and in the outflows (lower-excitation lines) with comparable abundances. COMs can coexist with hydrocarbons in the same protostellar sources, but they trace different components. In the near future, mid-IR observations with JWST-MIRI will provide complementary information about the hottest gas and the ice mantle content, at unprecedented sensitivity and at resolutions comparable to ALMA for the same sources., Comment: 38 pages, 33 figures, accepted for publication to A\&A

The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location of the "dust lines". This is a new concept in the field, linking the disk size at different continuum wavelengths with the radial distribution of grain populations of different sizes. We aim to use a dust evolution model to test the dependence of the dust line location on disk gas mass. In particular, we are interested in the reliability of the method for disks showing radial substructures, as recent high-resolution observations revealed. Our models show that the determination of the dust line location is a promising approach to the mass estimate of protoplanetay disks, but the exact relation between the dust line location and disk mass depends on the structure of the particular disk. We calibrated the relation for disks without evidence of radial structures, while for more complex structures we ran a simple dust evolution model. However, this method fails when there is evidence of strong dust traps. It is possible to reveal when dust evolution is dominated by traps, providing the necessary information for when the method should be applied with caution., 14 pages, 18 figures, accepted by Astronomy & Astrophysics

We report observations of polarized line and continuum emission from the disk of TW~Hya using the Atacama Large Millimeter/submillimeter Array. We target three emission lines, $^{12}$CO (3-2), $^{13}$CO (3-2) and CS (7-6), to search for linear polarization due to the Goldreich-Kylafis effect, while simultaneously tracing the continuum polarization morphology at 332\,GHz (900\,\micron{}), achieving a spatial resolution of 0.5\arcsec{} (30~au). We detect linear polarization in the dust continuum emission; the polarization position angles show an azimuthal morphology, and the median polarization fraction is $\sim$\,0.2\%, comparable to previous, lower frequency observations. Adopting a `shift-and-stack' technique to boost the sensitivity of the data, combined with a linear combination of the $Q$ and $U$ components to account for their azimuthal dependence, we detect weak linear polarization of $^{12}$CO and $^{13}$CO line emission at a $\sim 10\sigma$ and $\sim 5\sigma$ significance, respectively. The polarization was detected in the line wings, reaching a peak polarization fraction of $\sim 5\%$ and $\sim 3\%$ for the two molecules between disk radii of 0.5" and 1". The sign of the polarization was found to flip from the blue-shifted side of the emission to the red-shifted side, suggesting a complex, asymmetric polarization morphology. Polarization is not robustly detected for the CS emission; however, a tentative signal, comparable in morphology to that found for the $^{12}$CO and $^{13}$CO emission, is found at a $\lesssim 3\sigma$ significance. We are able to reconstruct a polarization morphology, consistent with the azimuthally averaged profiles, under the assumption that this is also azimuthally symmetric, which can be compared with future higher-sensitivity observations., Comment: Accepted by ApJ

Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. High-resolution observations of the CO emission around AGB stars now have the brightness temperature sensitivity to detect possible weak CO maser emission. We used high angular resolution observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the small-scale structure of CO $J=3-2$ emission around the oxygen-rich AGB star W Hya. We find CO maser emission amplifying the stellar continuum with an optical depth $\tau\approx-0.55$. The maser predominantly amplifies the limb of the star because CO $J=3-2$ absorption from the extended stellar atmosphere is strongest towards the centre of the star. The CO maser velocity corresponds to a previously observed variable component of high-frequency H$_2$O masers and with the OH maser that was identified as the amplified stellar image. This implies that the maser originates beyond the acceleration region and constrains the velocity profile since we find the population inversion primarily in the inner circumstellar envelope. We find that inversion can be explained by the radiation field at 4.6 $\mu$m and that the existence of CO maser emission is consistent with the estimated mass-loss rates for W Hya. However, the pumping mechanism requires a complex interplay between absorption and emission lines in the extended atmosphere. Excess from dust in the circumstellar envelope of W Hya is not sufficient to contribute significantly to the required radiation field at 4.6 $\mu$m. The interplay between molecular lines that cause the pumping can be constrained by future multi-level CO observations., Comment: Accepted for publication in A\&A, 14 pages, 14 figures

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

We present observations of a giant Lyman-alpha blob in the SSA22 proto-cluster at z=3.1, SSA22-LAB1, taken with the Atacama Large Millimeter/submillimeter Array (ALMA). Dust continuum, along with [C II]158um, and CO(4-3) line emission have been detected in LAB1, showing complex morphology and kinematics across a ~100 kpc central region. Seven galaxies at z=3.0987-3.1016 in the surroundings are identified in [C II] and dust continuum emission, with two of them potential companions or tidal structures associated with the most massive galaxies. Spatially resolved [C II] and infrared luminosity ratios for the widely distributed media (L[C II]/LIR~0.01-0.001) suggest that the observed extended interstellar media are likely to have originated from star-formation activity and the contribution from shocked gas is probably not dominant. LAB1 is found to harbour a total molecular gas mass Mmol=(8.7+/-2.0)e+10 Msun, concentrated in the core region of the Ly-alpha-emitting area. While (primarily obscured) star-formation activity in the LAB1 core is one of the most plausible power sources for the Ly-alpha emission, multiple major-mergers found in the core may also play a role in making LAB1 exceptionally bright and extended in Ly-alpha as a result of cooling radiation induced by gravitational interactions., 23 pages, 11 figures, accepted for publication in ApJ

The detection of phosphine (PH3) in the atmosphere of Venus has been recently reported based on millimeter-wave radio observations (Greaves et al. 2020), and its re-analyses (Greaves et al. 2021a/b). In this Matters Arising we perform an independent reanalysis, identifying several issues in the interpretation of the spectroscopic data. As a result, we determine sensitive upper-limits for PH3 in Venus' atmosphere (>75 km, above the cloud decks) that are discrepant with the findings in G2020 and G2021a/b. The measurements target the fundamental first rotational transition of PH3 (J=1-0) at 266.944513 GHz, which was observed with the James Clerk Maxwell Telescope (JCMT) in June 2017 and with the Atacama Large Millimeter/submillimeter Array (ALMA) in March 2019. This line's center is near the SO2 (J=309,21-318,24) transition at 266.943329 GHz (only 1.3 km/s away from the PH3 line) which represents a potential source of contamination. The JCMT and ALMA data, as presented in G2020, are at spectral resolutions comparable to the frequency separation of the two lines. Moreover, the spectral features identified are several km/s in width, and therefore do not permit distinct spectroscopic separation of the candidate spectral lines of PH3 and SO2. We present the radiative transfer modelling we have performed and then discuss the ALMA and JCMT analyses in turn.

Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star., 21 pages, 23 figures

We carried out new CO($J$ = 2-1) observations toward the mixed-morphology supernova remnant (SNR) W49B with the Atacama Large Millimeter/submillimeter Array (ALMA). We found that CO clouds at $\sim$10 km s$^{-1}$ show a good spatial correspondence with synchrotron radio continuum as well as an X-ray deformed shell. The bulk mass of molecular clouds accounts for the western part of the shell, not for the eastern shell where near-infrared H$_2$ emission is detected. The molecular clouds at $\sim$10 km s$^{-1}$ show higher kinetic temperature of $\sim$20-60 K, suggesting that modest shock-heating occurred. The expanding motion of the clouds with $\Delta V \sim$6 km s$^{-1}$ was formed by strong winds from the progenitor system. We argue that the barrel-like structure of Fe rich ejecta was possibly formed not only by an asymmetric explosion, but also by interactions with dense molecular clouds. We also found a negative correlation between the CO intensity and the electron temperature of recombining plasma, implying that the origin of the high-temperature recombining plasma in W49B can be understood as the thermal conduction model. The total energy of accelerated cosmic-ray protons $W_\mathrm{p}$ is estimated to be $\sim$$2\times 10^{49}$ erg by adopting an averaged gas density of $\sim$$650\pm200$ cm$^{-3}$. The SNR age-$W_\mathrm{p}$ diagram indicates that W49B shows one of the highest in-situ values of $W_\mathrm{p}$ in the gamma-ray bright SNRs., Comment: 19 pages, 10 figures, 2 tables, published in The Astrophysical Journal (ApJ)

We report the discovery of a very dense jet-like fast molecular outflow surrounded by a wide-angle wind in a massive young stellar object (MYSO) G18.88MME (stellar mass $\sim$8 M$_{\odot}$) powering an Extended Green Object G18.89$-$0.47. Four cores MM1-4 are identified in the Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum map (resolution $\sim$0.$"$8) toward G18.88MME, and are seen at the center of the emission structure (extent $\sim$0.3 pc $\times$ 0.2 pc) detected in the ALMA map. G18.88MME is embedded in the core MM1 (mass $\sim$13-18 M$_{\odot}$), where no radio continuum emission is detected. The molecular outflow centered at MM1 is investigated in the SiO(5-4), HC$_{3}$N(24-23) and $^{13}$CO(2-1) lines. The detection of HC$_{3}$N in the outflow is rare in MYSOs and indicates its very high density. The position-velocity diagrams display a fast narrow outflow (extent $\sim$28000 AU) and a slower wide-angle more extended outflow toward MM1, and both of these components show a transverse velocity gradient indicative of a possible rotation. All these observed features together make G18.88MME as a unique object for studying the unification of the jet-driven and wind-driven scenarios of molecular outflows in MYSOs., 5 pages, 6 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

The origin of the mysterious multiphase filamentary structures surrounding Brightest Cluster Galaxies (BCGs) remains unknown. We present Atacama Large Millimeter/submillimeter Array (ALMA) and Multi Unit Spectroscopic Explorer (MUSE) observations for a sample of 15 BCGs to investigate the origin and life-cycle of the gas. Those observations show clumpy and massive molecular filaments, preferentially located around the radio bubbles inflated by the active galactic nuclei (AGN). We investigate where the cold gas condenses from the intra-cluster medium, by comparing the radial extent of the filaments with predictions from numerical simulations.

We present a 4.7σ detection of the [OIII] 88 μm line in a gravitationally-lensed Lyman break galaxy, RXC J2248-ID3, using the Atacama Large Millimeter/submillimeter Array (ALMA). We did not detect [CII] 158 μm and rest-frame 90 μm dust continuum emission, suggesting that the bulk of the interstellar medium (ISM) is ionized. Our two-component SED model combining the previous Hubble Space Telescope (HST) data and new photometry obtained from Very Large Telescope (VLT), Spitzer and ALMA suggests the presence of young (∼2 Myr) and mature (∼600 Myr) stellar components with the metallicity of Z = 0.2Z⊙. Our findings are in contrast with previous results claiming a very young, metal-poor stellar component.

Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability., Published at MNRAS, link: http://adsabs.harvard.edu/abs/2019MNRAS.487.4473M Previous version (2019arXiv190602015M) was identical but had an incomplete title

The advent of the Atacama Large Millimeter/submillimeter Array (ALMA) has enabled a new era for studies of the formation and assembly of distant galaxies. Cosmological deep field surveys with ALMA and other interferometers have flourished in the last few years covering wide ranges of galaxy properties and redshift, and allowing us to gain critical insights into the physical mechanisms behind the galaxy growth. Here, we present a brief review of recent studies that aim to characterize the interstellar medium properties of galaxies at high redshift ($z>1$), focusing on blank-field ALMA surveys of dust continuum and molecular line emission. In particular, we show recent results from the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) large program., 10 pages, 7 figures. Proceeding of IAU Symposium No. 352, 2019 "Uncovering early galaxy evolution in the ALMA and JWST era"

We present our ongoing work of using two independent tracers to estimate the supermassive black hole mass in the nearby early-type galaxy NGC 6958; namely integrated stellar and molecular gas kinematics. We used data from the Atacama Large Millimeter/submillimeter Array (ALMA), and the adaptive-optics assisted Multi-Unit Spectroscopic Explorer (MUSE) and constructed state-of-the-art dynamical models. The different methods provide black hole masses of $(2.89\pm 2.05) \times 10^8M_{\odot}$ from stellar kinematics and $(1.35\pm 0.09) \times 10^8M_{\odot}$ from molecular gas kinematics which are consistent within their $3\sigma$ uncertainties. Compared to recent M$_{\rm BH}$ - $\sigma_{\rm e}$ scaling relations, we derive a slightly over-massive black hole. Our results also confirm previous findings that gas-based methods tend to provide lower black hole masses than stellar-based methods. More black hole mass measurements and an extensive analysis of the method-dependent systematics are needed in the future to understand this noticeable discrepancy., Comment: 5 pages, 3 figures, Proceeding of IAU Symposium 353, Galactic Dynamics in the Era of Large Surveys, ed. M. Valluri & J. A. Sellwood, Cambridge Univ. Press, in press; submitted September, 2019

Two laboratory emission spectrometers have been designed and described previously. Here, we present a follow-up study with special focus on absolute intensity calibration of the new SURFER-spectrometer (SUbmillimeter Receiver For Emission spectroscopy of Rotational transitions), operational between 300 and 400~GHz and mostly coincident with ALMA (Atacama Large Millimeter/submillimeter Array) Band 7. Furthermore, we present a feasibility study to extend the detection frequencies up to 2~THz. First results have been obtained using the SOFIA (Stratospheric Observatory for IR Astronomy) upGREAT laboratory setup at the University of Cologne. Pure rotational spectra of the complex molecule vinyl cyanide have been obtained and are used to give an estimate on the sensitivity to record ro-vibrational transitions of molecules with astrophysical importance at 2~THz., 4 pages, 3 Figures, conference proceedings IAUS 350, accepted amuscript

We present a polarimetric investigation of the protoplanetary nebula Frosty Leo performed with the Submillimeter Array. We were able to detect, in the low continuum level (peak at 14.4 mJy beam$^{-1}$), a marginal polarization at $\sim2.6\sigma$. The molecular line investigation based on the CO $J=3\rightarrow2$ emission shows a peak emission of 68.1 Jy beam$^{-1}$ km s$^{-1}$ and the polarization detection in this CO line is also marginal, with a peak at $\sim3.8\sigma$. In both cases, it was therefore not possible to use the electric vector maps (B-field) to accurately trace the magnetic field (B-field) within the PPN. The spatio-kinematic modelling realised with the different velocity channel maps indicates three main structures: a distorted torus accompanied by a bipolar outflow or jet aligned with its axis and a flattened spherical "cap". The comparison of the CO polarization segments with our model suggests that the polarized emission probably arises in the first two components., Comment: 8 pages, Accepted for publication in Monthly Notices of the Royal Astronomical Society

We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.

Validation of global climate models (GCMs) for planets in our solar system requires observational data, but observations from the orbit of Mars and its surface are limited in number and are constrained by their orbit or landing site. Ground-based observations of Mars can help by providing data across the entire Martian hemisphere, yet historically, ground-based observations at submillimeter wavelengths have been limited to disk-average, or at best, a few resolution elements across Mars. We used Atacama Large Millimeter/submillimeter Array (ALMA) observations of Mars to determine the spatial distribution of carbon monoxide in the Martian atmosphere, which can be related to the atmospheric temperature. ALMA’s comparably high spatial and spectral resolutions in the submillimeter wavelengths could allow the mapping of abundances and temperature profiles, and the comparison of these data to simulations generated by the Laboratoire de Météorologie Dynamique (LMD) Mars GCM. However, the long baselines associated with the high spatial resolution of ALMA introduced systematic errors that resulted in radiative transfer modeling degeneracies. We serve to provide insight to facilitate proposed ALMA observations of Mars in the future so that the systematic errors encountered within these observations might be avoided.

Ambipolar diffusion likely plays a pivotal role in the formation and evolution of dense cores in weakly-ionized molecular clouds. Linear analyses show that the evolutionary times and fragmentation scales are significantly greater than the hydrodynamic (Jeans) values even for clouds with mildly supercritical mass-to-flux ratio. We utilize values of fragmentation scales and growth times that correspond to typical ionization fractions within a molecular cloud, and apply to the context of the observed estimated lifetime of prestellar cores as well as the observed number of such embedded cores forming in a parent clump. By varying a single parameter, the mass-to-flux ratio, over the range of observationally measured densities, we fit the range of estimated prestellar core lifetimes ($\sim 0.1$ to a few Myr) identified with Herschel as well as the number of embedded cores formed in a parent clump measured in Perseus with the Submillimeter Array (SMA). Our model suggests that the prestellar cores are formed with a transcritical mass-to-flux ratio and higher densities correspond to somewhat higher mass-to-flux ratio but the normalized mass-to-flux ratio $\mu$ remains in the range $1 \lesssim \mu \lesssim 2$. Our best-fit model exhibits $B \propto n^{0.43}$ for prestellar cores, due to partial flux-freezing as a consequence of ambipolar diffusion., Comment: Accepted For Publication in Astronomy & Astrophysics Letters, 9 pages, 5 figures

We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration and its derivative in the pulsar search, and also search in segments of the entire observation to compensate for potential orbital motion of the pulsar. While no new pulsar is found, our observations yield the polarization profile of the GC magnetar PSR J1745-2900 at mm-wavelength for the first time, which turns out to be nearly 100 % linearly polarized. Additionally, we estimate the survey sensitivity placed by both system and red noise, and evaluate its capability of finding pulsars in orbital motion with either Sgr A* or a binary companion. We show that the survey is sensitive to only the most luminous pulsars in the known population, and future observations with ALMA in Band-1 will deliver significantly deeper survey sensitivity on the GC pulsar population., 16 pages, 11 figures, accepted for publication in ApJ

Observations of the Galactic Center supermassive black hole Sagittarius A∗ (Sgr A∗) with very long-baseline interferometry (VLBI) are affected by interstellar scattering along our line of sight. At long radio observing wavelengths (≲1 cm), the scattering heavily dominates image morphology. At 3.5 mm (86 GHz), the intrinsic source structure is no longer sub-dominant to scattering, and thus the intrinsic emission from Sgr A∗ is resolvable with the Global Millimeter VLBI Array (GMVA). Long-baseline detections to the phased Atacama Large Millimeter/submillimeter Array (ALMA) in 2017 provided new constraints on the intrinsic and scattering properties of Sgr A∗, but the stochastic nature of the scattering Requires multiple observing epochs to reliably estimate its statistical properties. We present new observations with the GMVA+ALMA, taken in 2018, which confirm non-Gaussian structure in the scattered image seen in 2017. In particular, the ALMA-GBT baseline shows more flux density than expected for an anistropic Gaussian model, providing a tight constraint on the source size and an upper limit on the dissipation scale of interstellar turbulence. We find an intrinsic source extent along the minor axis of ∼100 μas both via extrapolation of longer wavelength scattering constraints and direct modeling of the 3.5 mm observations. Simultaneously fitting for the scattering parameters, we find an at-most modestly asymmetrical (major-to-minor axis ratio of 1.5 0.2) intrinsic source morphology for Sgr A∗. © 2021. The American Astronomical Society. All rights reserved., This work is supported by the ERC Synergy grant "BlackHoleCam: Imaging the Event Horizon of Black Holes", grant 610058. We thank the National Science Foundation (AST-1716536, AST-1440254, AST-1935980, AST-2034306) and the Gordon and Betty Moore Foundation (GBMF-5278) for financial support of this work. This work was supported in part by the black hole Initiative, which is funded by grants from the John Templeton Foundation and the Gordon and Betty Moore Foundation to Harvard University. A.C. is supported by Hubble Fellowship grant HST-HF2-51431.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. I.C. is supported by the National Research Foundation of Korea (NRF) via a Global PhD Fellowship grant (NRF-2015H1A2A1033752). L.L. acknowledges the financial support of DGAPA, UNAM (projects IN112417 and IN112820), and CONACyT, Mexico (projects 275201-Agencia Espacial Mexicana; and 263356-Ciencia de Frontera). M.K. was supported by JSPS KAKENHI grant Nos. JP18K03656 and JP18H03721. This paper makes use of the following ALMA data: ADS/JAO.ALMA2017.1.00795.V. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The ALMA data required non-standard processing by the VLBI QA2 team (C. Goddi, I. Marti-Vidal, G. B. Crew, H. Rottmann, and H. Messias). This research has made use of data obtained with the Global Millimeter VLBI Array (GMVA), coordinated by the VLBI group at the Max-Planck-Institut fur Radioastronomie (MPIfR). The GMVA consists of telescopes operated by MPIfR, IRAM, Onsala, Metsahovi, Yebes, the Korean VLBI Network, the Green Bank Observatory and the Very Long Baseline Array (VLBA). The VLBA is a facility of the National Science Foundation under cooperative agreement by Associated Universities, Inc. The data were correlated at the DiFX correlator of the MPIfR in Bonn, Germany. This work made use of the Swinburne University of Technology software correlator (Deller et al. 2011), developed as part of the Australian Major National Research Facilities Programme and operated under licence., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

Remote and in-situ observations of cometary gases have revealed the presence of a wealth of complex organic molecules, including carbon chains, alcohols, imines and the amino acid glycine. Such chemical complexity in cometary material implies that impacts by comets could have supplied reagents for prebiotic chemistry to young planetary surfaces. However, the assumption that some of the molecules observed in cometary comae at millimetre wavelengths originate from ices stored inside the nucleus has not yet been proven. In fact, the comae of moderately-active comets reach sufficient densities within a few thousand kilometers of the nucleus for an active (solar radiation-driven) photochemistry to ensue. Here we present results from our latest chemical-hydrodynamic models incorporating an updated reaction network, and show that the commonly-observed HC3N (cyanoacetylene) and NH2CHO (formamide) molecules can be efficiently produced in cometary comae as a result of two-body, neutral-neutral, gas-phase reactions involving well-known coma gases. In the presence of a near-nucleus distributed source of CN (similar to that observed by the Rosetta spacecraft at comet 67P), we find that sufficient HC$_3$N and NH2CHO can be synthesized to match the abundances of these molecules in previous observations of Oort Cloud comets. The precise origin of these (and other) complex organic molecules in cometary comae can be verified through interferometric mapping observations, for example, using the Atacama Large Millimeter/submillimeter Array (ALMA)., Accepted for publication in MNRAS, April 15th 2021

We present multi-scale and multi-wavelength data of the Galactic HII region G25.4-0.14 (hereafter G25.4NW, distance ~5.7 kpc). The SHARC-II 350 micron continuum map displays a hub-filament configuration containing five parsec scale filaments and a central compact hub. Through the 5 GHz radio continuum map, four ionized clumps (i.e., Ia-Id) are identified toward the central hub, and are powered by massive OB-stars. The Herschel temperature map depicts the warm dust emission (i.e., Td ~23-39 K) toward the hub. High resolution Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum map (resolution ~0".82 X 0".58) reveals three cores (c1-c3; mass ~80-130 Msun) toward the ionized clumps Ia, and another one (c4; mass ~70 Msun) toward the ionized clump Ib. A compact near-infrared (NIR) emission feature (extent ~0.2 pc) is investigated toward the ionized clump Ia excited by an O8V-type star, and contains at least three embedded K-band stars. In the direction of the ionized clump Ia, the ALMA map also shows an elongated feature (extent ~0.2 pc) hosting the cores c1-c3. All these findings together illustrate the existence of a small cluster of massive stars in the central hub. Considering the detection of the hub-filament morphology and the spatial locations of the mm cores, a global non-isotropic collapse (GNIC) scenario appears to be applicable in G25.4NW, which includes the basic ingredients of the global hierarchical collapse and clump-fed accretion models. Overall, the GNIC scenario explains the birth of massive stars in G25.4NW., 12 pages, 6 figures, 1 Table, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS) Journal

We present detailed radio observations of the tidal disruption event (TDE) AT2019dsg, obtained with the Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), and spanning $55-560$ days post-disruption. We find that the peak brightness of the radio emission increases until ~200 days and subsequently begins to decrease steadily. Using the standard equipartition analysis, including the effects of synchrotron cooling as determined by the joint VLA-ALMA spectral energy distributions, we find that the outflow powering the radio emission is in roughly free expansion with a velocity of $\approx 0.07c$, while its kinetic energy increases by a factor of about 5 from 55 to 200 days and plateaus at $\approx 5\times 10^{48}$ erg thereafter. The ambient density traced by the outflow declines as $\approx R^{-1.6}$ on a scale of $\approx (1-4)\times 10^{16}$ cm ($\approx 6300-25000$ $R_s$), followed by a steeper decline to $\approx 6\times 10^{16}$ cm ($\approx 37500$ $R_s$). Allowing for a collimated geometry, we find that to reach even mildly relativistic velocities ($\Gamma=2$) the outflow requires an opening angle of $\theta_j\approx 2^\circ$, which is narrow even by the standards of GRB jets; a truly relativistic outflow requires an unphysically narrow jet. The outflow velocity and kinetic energy in AT2019dsg are typical of previous non-relativistic TDEs, and comparable to those from Type Ib/c supernovae, raising doubts about the claimed association with a high-energy neutrino event., Comment: Accepted to the Astrophysical Journal

We report an analysis of the dust disk around DM~Tau, newly observed with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm. The ALMA observations with high sensitivity (8.4~$\mu$Jy/beam) and high angular resolution (35~mas, 5.1~au) detect two asymmetries on the ring at $r\sim$20~au. They could be two vortices in early evolution, the destruction of a large scale vortex, or double continuum emission peaks with different dust sizes. We also found millimeter emissions with $\sim$50~$\mu$Jy (a lower limit dust mass of 0.3~$M_{\rm Moon}$) inside the 3-au ring. To characterize these emissions, we modeled the spectral energy distribution (SED) of DM~Tau using a Monte Carlo radiative transfer code. We found that an additional ring at $r=$ 1~au could explain both the DM~Tau SED and the central point source. The disk midplane temperature at the 1-au ring calculated in our modeling is less than the typical water sublimation temperature of 150~K, prompting the possibility of forming small icy planets there., Comment: 26 pages, 17 figures

We present the first linear-polarization mosaicked observations performed by the Atacama Large Millimeter/submillimeter Array (ALMA). We mapped the Orion-KLeinmann-Low (Orion-KL) nebula using super-sampled mosaics at 3.1 and 1.3 mm as part of the ALMA Extension and Optimization of Capabilities (EOC) program. We derive the magnetic field morphology in the plane of the sky by assuming that dust grains are aligned with respect to the ambient magnetic field. At the center of the nebula, we find a quasi-radial magnetic field pattern that is aligned with the explosive CO outflow up to a radius of approximately 12 arc-seconds (~ 5000 au), beyond which the pattern smoothly transitions into a quasi-hourglass shape resembling the morphology seen in larger-scale observations by the James-Clerk-Maxwell Telescope (JCMT). We estimate an average magnetic field strength $\langle B\rangle = 9.4$ mG and a total magnetic energy of 2 x 10^45 ergs, which is three orders of magnitude less than the energy in the explosive CO outflow. We conclude that the field has been overwhelmed by the outflow and that a shock is propagating from the center of the nebula, where the shock front is seen in the magnetic field lines at a distance of ~ 5000 au from the explosion center., Comment: Accepted for publication in ApJ

We present a method to mitigate the atmospheric effects (residual atmospheric lines) in single-dish radio spectroscopy caused by the elevation difference between the target and reference positions. The method is developed as a script using the Atmospheric Transmission at Microwaves (ATM) library built into the Common Astronomy Software Applications (CASA) package. We apply the method to the data taken with the Total Power Array of the Atacama Large Millimeter/submillimeter Array. The intensities of the residual atmospheric (mostly O3) lines are suppressed by, typically, an order of magnitude for the tested cases. The parameters for the ATM model can be optimized to minimize the residual line and, for a specific O3 line at 231.28 GHz, a seasonal dependence of a best-fitting model parameter is demonstrated. The method will be provided as a task within the CASA package in the near future. The atmospheric removal method we developed can be used by any radio/millimeter/submillimeter observatory to improve the quality of its spectroscopic measurements., 21 pages, 12 figures, accepted for publication in PASP

Debris disks are tenuous, dusty belts surrounding main sequence stars generated by collisions between planetesimals. HD 206893 is one of only two stars known to host a directly imaged brown dwarf orbiting interior to its debris ring, in this case at a projected separation of 10.4 au. Here we resolve structure in the debris disk around HD 206893 at an angular resolution of 0.6" (24 au) and wavelength of 1.3 mm with the Atacama Large Millimeter/submillimeter Array (ALMA). We observe a broad disk extending from a radius of, 24 pages, 12 figures, Accepted for publication in ApJ

We observed the high-mass star-forming region G335.579-0.292 with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz with an angular resolution of 0.3'' ($\sim 1000$ au resolution at the source distance). G335.579-0.292 hosts one of the most massive cores in the Galaxy (G335-MM1). The continuum emission shows that G335-MM1 fragments into at least five sources, while molecular line emission is detected in two of the continuum sources (ALMA1 and ALMA3). We found evidence of large and small scale infall in ALMA1 revealed by an inverse P-Cygni profile and the presence of a blue-shifted spot at the center of the first moment map of the CH$_3$CN emission. In addition, hot gas expansion in the innermost region is unveiled by a red-shifted spot in the first moment map of HDCO and (CH$_3$)$_2$CO (both with $E_u > 1100$ K). Our modeling reveals that this expansion motion originates close to the central source, likely due to reversal of the accretion flow induced by the expansion of the HII region, while infall and rotation motions originate in the outer regions. ALMA3 shows clear signs of rotation, with a rotation axis inclination with respect to the line of sight close to $90^\circ$, and a system mass (disk + star) in the range of 10-30 M$_\odot$., 15 pages, 12 figures, 3 tables. Accepted for publication in ApJ

We have observed the very low-mass Class 0 protostar IRAS 15398−3359 at scales ranging from 50 to 1800 au, as part of the Atacama Large Millimeter/Submillimeter Array Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source in a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398−3359, by 1200 au. The arc-like structure is blueshifted with respect to the systemic velocity. A velocity gradient of 1.2 km s−1 over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398−3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution. The authors thank the anonymous reviewer for invaluable comments. They are also grateful to Hauyu Baobab Liu for useful discussions. This paper makes use of the following ALMA data set: ADS/JAO.ALMA#2018.1.01205.L (PI: Satoshi Yamamoto). ALMA is a partnership of the ESO (representing its member states), the NSF (USA) and NINS (Japan), together with the NRC (Canada) and the NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by the ESO, the AUI/NRAO, and the NAOJ. The authors thank to the ALMA staff for their excellent support. This work is supported by the projects PRIN-INAF-MAIN-STREAM 2017 Protoplanetary disks seen through the eyes of new-generation instruments, ARIEL and the astrochemical link between circumstellar disks and planets, and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programmes: The Dawn of Organic Chemistry (DOC), grant agreement No. 741002, and Astro-Chemistry Origins (ACO), grant No. 811312. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. I.J.-S. has received partial support from the Spanish FEDER (project number ESP2017-86582-C4-1-R) and the State Research Agency (AEI; project number PID2019-105552RB-C41). D.J. is supported by NRC Canada and by an NSERC Discovery Grant. This project is also supported by a Grant-in-Aid from Japan Society for the Promotion of Science (KAKENHI: Nos. 18H05222, 19H05069, 19K14753). Yuki Okoda thanks the Advanced Leading Graduate Course for Photon Science (ALPS) and Japan Society for the Promotion of Science (JSPS) for financial support. Peerreview