Continuum Data Research Articles

Deeply Buried Nuclei in the Infrared-luminous Galaxies NGC 4418 and Arp 220. I. ALMA Observations at λ = 1.4–0.4 mm and Continuum Analysis

Abstract We observed with Atacama Large Millimeter/submillimeter Array three deeply buried nuclei in two galaxies, NGC 4418 and Arp 220, at ∼0.″2 resolution over a total bandwidth of 67 GHz in f rest = 215–697 GHz. Here we (1) introduce our program, (2) describe our data reduction method for wide-band, high-resolution imaging spectroscopy, (3) analyze in visibilities the compact nuclei with line forests, (4) develop a continuum-based estimation method of dust opacity and gas column density in heavily obscured nuclei, which uses the buried galactic nuclei (BGN) model and is sensitive to log ( N H 2 / cm − 2 ) ∼ 25–26 at λ ∼ 1 mm, and (5) present the continuum data and diagnosis of our targets. The three continuum nuclei have major-axis FWHMs of ∼0.″1–0.″3 (20–140 pc) aligned to their rotating nuclear disks of molecular gas. However, each nucleus is described better with two or three concentric components than with a single Gaussian. The innermost cores have sizes of 0.″05–0.″10 (8–40 pc), peak brightness temperatures of ∼100–500 K at 350 GHz, and more fractional flux at lower frequencies. The intermediate components correspond to the nuclear disks. They have axial ratios of ≈0.5 and hence inclinations ≳60°. The outermost elements include the bipolar outflow from Arp 220W. We estimate 1 mm dust opacity of τ d,1 mm ≈ 2.2, 1.2, and ≲0.4, respectively, for NGC 4418, Arp 220W, and Arp 220E. The first two correspond to log ( N H / cm − 2 ) ∼ 26 for conventional dust-opacity laws, and hence the nuclei are highly Compton thick.

Detection of Substructures in Young Transition Disk WL 17

Abstract WL 17 is a young transition disk in the Ophiuchus L1688 molecular cloud complex. Even though WL 17 is among the brightest disks in L1688 and massive enough to expect dust self-scattering, it was undetected in polarization down to Atacama Large Millimeter/submillimeter Array’s (ALMA’s) instrument sensitivity limit. Such low polarization fractions could indicate unresolved polarization within the beam or optically thin dust emission. We test the latter case by combining the high-sensitivity 233 GHz Stokes I data from the polarization observations with previous ALMA data at 345 and 100 GHz. We use simple geometric shapes to fit the observed disk visibilities in each band. Using our simple models and assumed dust temperature profiles, we estimate the optical depth in all three bands. The optical depth at 233 GHz peaks at τ 233 ∼ 0.3, which suggests the dust emission may not be optically thick enough for dust self-scattering to be efficient. We also find the higher-sensitivity 233 GHz data show substructure in the disk for the first time. The substructure appears as brighter lobes along the major axis, on either side of the star. We attempt to fit the lobes with a simple geometric model, but they are unresolved in the 233 GHz data. We propose that the disk may be flared at 1 mm such that there is a higher column of dust along the major axis than the minor axis when viewed at an inclination. These observations highlight the strength of high-sensitivity continuum data from dust polarization observations to study disk structures.

The SPECFIND V3.0 catalog of radio continuum cross-identifications and spectra: Reaching lower frequencies

Context. Many radio continuum catalogs with different sensitivity limits and spatial resolutions are published via the VizieR database at the Centre de Données astronomiques de Strasbourg. The diversity of spatial resolutions of different catalogs makes the cross-identification of different flux density measurements of individual sources complex. The SPECFIND tool is able to handle radio surveys at different frequencies from different instruments with different resolutions. Aims. Since the former version of the SPECFIND catalog was released ten years ago, hundreds of new radio continuum catalogs have been published. We upgraded the SPECFIND tool to reach a wider frequency range, especially the lower-frequency radio regime, as well as to have better spatial sky coverage. Methods. We adapted selection criteria and applied them to all of the radio catalogs listed in the VizieR database to define a final sample of new catalogs. We unified the new catalogs and implemented them in the SPECFIND tool. The new SPECFIND V3.0 radio cross-identification catalog was constructed using 204 input tables from 160 VizieR radio continuum catalogs to cross-identify flux density measurements of individual sources and fit their spectral slopes. We discuss the frequency and sky coverage of all processed catalogs and compare the results to the previous version. Furthermore, we present and investigate peaked spectrum (PS) sources with spectral breaks around 1.4 GHz and 325 MHz. Results. By increasing the number of input catalog tables that were implemented in SPECFIND from 115 to 204 (89 new catalog tables and two updates), we improved the number of resulting spectra from ∼107 500 to ∼340 000 and increased the number of cross-identified sources from ∼600 000 to ∼1.6 million. The final SPECFIND V3.0 catalog is publicly available via VizieR. By applying SPECFIND to two subsamples of the catalogs with frequency cuts at 325 MHz and 1.4 GHz, spectral break and PS source candidates could be identified. We encourage follow-up observations of these candidates to confirm their nature because the population we identify has a relatively low reliability. Conclusions. The SPECFIND V3.0 catalog is a very useful resource and a powerful open access tool, reachable via VizieR. By tripling the resulting spectra and including many radio continuum surveys from the last 50 years, we provide a significantly extended catalog of cross-identified radio continuum sources. Furthermore, the SIMBAD database will be updated using the SPECFIND V3.0 catalog and will contain more radio continuum data, serving the needs of future projects.

MIGHTEE: total intensity radio continuum imaging and the COSMOS/XMM-LSS Early Science fields

ABSTRACT MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectropolarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image ∼20 deg2 over the COSMOS, E-CDFS, ELAIS-S1, and XMM-Newton Large Scale Structure field (XMM-LSS) extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength data sets from numerous existing and forthcoming observational campaigns. Here, we describe and validate the data processing strategy for the total intensity continuum aspect of MIGHTEE, using a single deep pointing in COSMOS (1.6 deg2) and a three-pointing mosaic in XMM-LSS (3.5 deg2). The processing includes the correction of direction-dependent effects, and results in thermal noise levels below 2 $\mathrm{\mu }$Jy beam−1 in both fields, limited in the central regions by classical confusion at ∼8 arcsec angular resolution, and meeting the survey specifications. We also produce images at ∼5 arcsec resolution that are ∼3 times shallower. The resulting image products form the basis of the Early Science continuum data release for MIGHTEE. From these images we extract catalogues containing 9896 and 20 274 radio components in COSMOS and XMM-LSS, respectively. We also process a close-packed mosaic of 14 additional pointings in COSMOS and use these in conjunction with the Early Science pointing to investigate methods for primary beam correction of broad-band radio images, an analysis that is of relevance to all full-band MeerKAT continuum observations, and wide-field interferometric imaging in general. A public release of the MIGHTEE Early Science continuum data products accompanies this article.

A Decade of SCUBA-2: A Comprehensive Guide to Calibrating 450 μm and 850 μm Continuum Data at the JCMT

Abstract The Submillimetre Common User Bolometer Array 2 (SCUBA-2) is the James Clerk Maxwell Telescope’s continuum imager, operating simultaneously at 450 and 850 μm. SCUBA-2 was commissioned in 2009–2011, and since that time, regular observations of point-like standard sources have been performed whenever the instrument is in use. Expanding the calibrator observation sample by an order of magnitude compared to previous work, in this paper we derive updated opacity relations at each wavelength for a new atmospheric extinction correction, analyze the Flux Conversion Factors used to convert instrumental units to physical flux units as a function of date and observation time, present information on the beam profiles for each wavelength, and update secondary calibrator source fluxes. Between 07:00 and 17:00 UTC, the portion of the night that is most stable to temperature gradients that cause dish deformation, the total flux uncertainty and the peak flux uncertainty measured at 450 μm are found to be 14% and 17%, respectively. Measured at 850 μm, the total flux and peak flux uncertainties are 6% and 7%, respectively. The analysis presented in this work is applicable to all SCUBA-2 projects observed since 2011.

Accurate dust temperature determination in a z = 7.13 galaxy

ABSTRACT We report ALMA Band 9 continuum observations of the normal, dusty star-forming galaxy A1689-zD1 at z = 7.13, resulting in a ∼4.6 σ detection at 702 GHz. For the first time, these observations probe the far-infrared spectrum shortward of the emission peak of a galaxy in the Epoch of Reionization (EoR). Together with ancillary data from earlier works, we derive the dust temperature, Td, and mass, Md, of A1689-zD1 using both traditional modified blackbody spectral energy density fitting, and a new method that relies only on the [C ii] 158 μm line and underlying continuum data. The two methods give $T_{\rm d} = (42^{+13}_{-7}, 40^{+13}_{-7}$) K, and $M_{\rm d} = (1.7^{+1.3}_{-0.7}, 2.0^{+1.8}_{-1.0})\, \times {}\, 10^{7} \, \mathrm{ M}_{\odot }$. Band 9 observations improve the accuracy of the dust temperature (mass) estimate by ∼50 per cent (6 times). The derived temperatures confirm the reported increasing Td-redshift trend between z = 0 and 8; the dust mass is consistent with a supernova origin. Although A1689-zD1 is a normal UV-selected galaxy, our results, implying that ∼85 per cent of its star-formation rate is obscured, underline the non-negligible effects of dust in EoR galaxies.

Probing protoplanetary disk evolution in the Chamaeleon II region

Context. Characterizing the evolution of protoplanetary disks is necessary to improve our understanding of planet formation. Constraints on both dust and gas are needed to determine the dominant disk dissipation mechanisms. Aims. We aim to compare the disk dust masses in the Chamaeleon II (Cha II) star-forming region with other regions with ages between 1 and 10 Myr. Methods. We use ALMA band 6 observations (1.3 mm) to survey 29 protoplanetary disks in Cha II. Dust mass estimates are derived from the continuum data. Results. Out of our initial sample of 29 disks, we detect 22 sources in the continuum, 10 in 12CO, 3 in 13CO, and none in C18O (J = 2−1). Additionally, we detect two companion candidates in the continuum and 12CO emission. Most disk dust masses are lower than 10 M⊕, assuming thermal emission from optically thin dust. Including non-detections, we derive a median dust mass of 4.5 ± 1.5 M⊕ from survival analysis. We compare consistent estimations of the distributions of the disk dust mass and the disk-to-stellar mass ratios in Cha II with six other low mass and isolated star-forming regions in the age range of 1–10 Myr: Upper Sco, CrA, IC 348, Cha I, Lupus, and Taurus. When comparing the dust-to-stellar mass ratio, we find that the masses of disks in Cha II are statistically different from those in Upper Sco and Taurus, and we confirm that disks in Upper Sco, the oldest region of the sample, are statistically less massive than in all other regions. Performing a second statistical test of the dust mass distributions from similar mass bins, we find no statistical differences between these regions and Cha II. Conclusions. We interpret these trends, most simply, as a sign of decline in the disk dust masses with time or dust evolution. Different global initial conditions in star-forming regions may also play a role, but their impact on the properties of a disk population is difficult to isolate in star-forming regions lacking nearby massive stars.

Observation of the Apoapsis of S62 in 2019 with NIRC2 and SINFONI

Abstract Given the increased attention toward the detection of faint stars in the Galactic center, we would like to address the detectability of S62 in its apoapsis with the Spectrograph for INtegral Field Observations in the Near Infrared at the Very Large Telescope (SINFONI; VLT) and the Near-Infrared Camera 2 (NIRC2; KECK) in this work. Because of the nearby stars and the chance of confusion, we are using Lucy–Richardson deconvolved images to detect S62 on its Keplerian orbit around Sgr A* with a period of less than 10 yr. We use the same data set as for S62 to trace additionally the S-cluster star S29 at the expected position based on the orbital elements presented in this work. To verify the results of the filtering technique, we are analyzing K-band continuum data of the same epoch independently observed with NIRC2/KECK. Based on the well-derived orbit of S62, we find the star in projection at the expected position in 2019.4 and 2019.5. By analyzing the SINFONI data of 2019.5, we confirm the 16.1 ± 0.2 mag for S62 that was formerly derived with NAOS-CONICA (NACO; VLT). We base our NACO imaging analysis on the robust data set that was previously used to investigate the Schwarzschild precision of S2. We also present a critical discussion of the elsewhere proposed linear trajectory of S62 and its disputed identification with a 19 mag star found with GRAVITY mounted at the VLT Interferometer.

Census of High- and Medium-mass Protostars. V. CO Abundance and the Galactic X CO Factor

Abstract We present the second dust continuum data release in the Census of High- and Medium-mass Protostars (CHaMP), expanding the methodology trialed in Pitts et al. to the entire CHaMP survey area (280° < ℓ < 300°, − 4° < b < + 2°). This release includes maps of dust temperature (T d), H2 column density ( N H 2 ), gas-phase CO abundance, and temperature–density plots for every prestellar clump with Herschel coverage, showing no evidence of internal heating for most clumps in our sample. We show that CO abundance is a strong function of T d and can be fit with a second-order polynomial in log-space, with a typical dispersion of a factor of 2–3. The CO abundance peaks at 20.0 − 1.0 + 0.4 K with a value of 7.4 − 0.3 + 0.2 × 10 − 5 per H2; the low T d at which this maximal abundance occurs relative to laboratory results is likely due to interstellar UV bombardment in the largest survey fields. Finally, we show that, as predicted by theoretical literature and hinted at in previous studies of individual clouds, the conversion factor from integrated 12CO line intensity ( I 12 CO ) to N H 2 , the X CO factor, varies as a broken power law in I 12 CO with a transition zone between 70 and 90 K km s−1. The X CO function we propose has N H 2 ∝ I 12 CO 0.51 for I 12 CO ≲ 70 K km s−1 and N H 2 ∝ I 12 CO 2.3 for I 12 CO ≳ 90 K km s−1. The high- I 12 CO side should be generalizable with known adjustments for metallicity, but the influence of interstellar UV fields on the low- I 12 CO side may be sample specific. We discuss how these results expand on previous works in the CHaMP series and help tie together observational, theoretical, and laboratory studies on CO over the past decade.

No Impact of Core-scale Magnetic Field, Turbulence, or Velocity Gradient on Sizes of Protostellar Disks in Orion A

Abstract We compared the sizes and fluxes of a sample of protostellar disks in Orion A measured with the Atacama Large Millimeter/submillimeter Array 0.87 mm continuum data from the VLA/ALMA Nascent Disk and Multiplicity survey with the physical properties of their ambient environments on the core scale of 0.6 pc estimated with the Green Bank Ammonia Survey data and the legacy catalog datasets of the polarimeter for the Submillimeter Common-User Bolometer Array. We did not find any significant dependence of the disk radii and continuum fluxes on a single parameter on the core scale, such as nonthermal line width, magnetic field orientation and strength, or magnitude and orientation of the velocity gradient. Among these parameters, we only found a positive correlation between the magnitude of the velocity gradient and the nonthermal line width. Thus, the observed velocity gradients are more likely related to turbulent motion but not large-scale rotation. Our results of no clear dependence of the disk radii on these parameters are more consistent with the expectation from nonideal MHD simulations of disk formation in collapsing cores, where the disk size is self-regulated by magnetic braking and diffusion, compared to other simulations that only include turbulence and/or a magnetic field misaligned with the rotational axis. Therefore, our results could hint that the nonideal MHD effects play a more important role in the disk formation. Nevertheless, we cannot exclude the influences on the observed disk size distribution by dynamical interaction in a stellar cluster or amounts of angular momentum on the core scale, which cannot be probed with the current data.

A global view on star formation: The GLOSTAR Galactic plane survey

Context. The properties of the population of Galactic supernova remnants (SNRs) are essential to our understanding of the dynamics of the interstellar medium (ISM) in the Milky Way. However, the completeness of the catalog of Galactic SNRs is expected to be only ~30%, with on order 700 SNRs yet to be detected. Deep interferometric radio continuum surveys of the Galactic plane help in rectifying this apparent deficiency by identifying low surface brightness SNRs and compact SNRs that have not been detected in previous surveys. However, SNRs are routinely confused with H II regions, which can have similar radio morphologies. Radio spectral index, polarization, and emission at mid-infrared (MIR) wavelengths can help distinguish between SNRs and H II regions. Aims. We aim to identify SNR candidates using continuum images from the Karl G. Jansky Very Large Array GLObal view of the STAR formation in the Milky Way (GLOSTAR) survey. Methods. GLOSTAR is a C-band (4–8 GHz) radio wavelength survey of the Galactic plane covering 358° ≤ l ≤ 60°, |b|≤ 1°. The continuum images from this survey, which resulted from observations with the most compact configuration of the array, have an angular resolution of 18″. We searched for SNRs in these images to identify known SNRs, previously identified SNR candidates, and new SNR candidates. We study these objects in MIR surveys and the GLOSTAR polarization data to classify their emission as thermal or nonthermal. Results. We identify 157 SNR candidates, of which 80 are new. Polarization measurements provide evidence of nonthermal emission from nine of these candidates. We find that two previously identified candidates are filaments. We also detect emission from 91 of the 94 known SNRs in the survey region. Four of these are reclassified as H II regions following detection in MIR surveys. Conclusions. The better sensitivity and resolution of the GLOSTAR data have led to the identification of 157 SNR candidates, along with the reclassification of several misidentified objects. We show that the polarization measurements can identify nonthermal emission, despite the diffuse Galactic synchrotron emission. These results underscore the importance of higher resolution and higher sensitivity radio continuum data in identifying and confirming SNRs.

Identification of pre-stellar cores in high-mass star forming clumps via H2D+ observations with ALMA

Context. The different theoretical models concerning the formation of high-mass stars make distinct predictions regarding their progenitors, which are the high-mass pre-stellar cores. However, no conclusive observation of such objects has been made to date. Aims. We aim to study the very early stages of high-mass star formation in two infrared-dark massive clumps. Our goal is to identify the core population that they harbour and to investigate their physical and chemical properties at high spatial resolution. Methods. We obtained Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 6 observations of continuum emission at 0.8 mm and of the ortho-H2D+ transition at 372 GHz towards the two clumps. We used the SCIMES algorithm to identify substructures (i.e. cores) in the position-position-velocity space, finding 16 cores. We modelled their observed spectra using a Bayesian fitting approach in the approximation of local thermodynamic equilibrium. We derived the centroid velocity, the line width, and the molecular column density maps. We also studied the correlation between the continuum and molecular data, which in general do not present the same structure. Results. We report, for the first time, the detection of ortho-H2D+ in high-mass star-forming regions performed with an interferometer. The molecular emission shows narrow and subsonic lines, suggesting that locally, the temperature of the gas is below 10 K. From the continuum emission, we estimated the cores’ total masses and compare them with the respective virial masses. We also computed the volume density values, which are found to be higher than 106 cm−3. Conclusions. Our data confirm that ortho-H2D+ is an ideal tracer of cold and dense gas. Interestingly, almost all the H2D+-identified cores are less massive than ≈13 M⊙, with the exception of one core in AG354, which could be as massive as 39 M⊙ under the assumption of low dust temperature (5 K). Furthermore, most of them are sub-virial and larger than their Jeans masses. These results are hence difficult to explain in the context of the turbulent accretion models, which predict massive and virialised pre-stellar cores. However, we cannot exclude that the cores are still in the process of accreting mass and that magnetic fields are providing enough support for the virialisation. ALMA could also be seeing only the innermost parts of the cores, and hence the cores’ total masses could be higher than inferred in this work. Furthermore, we note that the total masses of the investigated clumps are below the average for typical high-mass clumps, and thus studies of more massive sources are needed.

A Semi-Automatic Method to Measure the Rotation of Sunspots

Sunspots have been observed to undergo rotation about their umbral centre. This is typically a slow rotation, with even the fastest sunspot rotations only reaching angular velocities of a few degrees per hour. This rotation may inject magnetic energy into the Sun’s atmosphere, which can be stored in the coronal magnetic field and later released in eruptive events such as solar flares and coronal mass ejections. To usefully investigate rotating sunspots long periods of data need to be analysed, often of the order of several days, to build up a bulk rotation profile for the sunspot over time. This article outlines a semi-automated approach for analysing series of solar continuum data to extract the rotation profile of a sunspot as it transits across the solar disc. Moving towards an automated approach is vital for generating large, unbiased statistical samples of rotating sunspots in order to understand their contribution to solar activity. Existing methods typically focus on sunspots near disc centre for short time periods, neglecting much of the rotation history of the sunspot. The method is tested on six sunspots observed in continuum data from the Helioseismic and Magnetic Imager (HMI) instrument on board the Solar Dynamics Observatory (SDO). These have been chosen to test the method for a range of different types of sunspots, including well-behaved sunspots, shape-changing sunspots, fast rotators, non-rotators, and interacting sunspots. The rotation profiles are compared by eye to animations of the sunspot from the data and are in acceptable visual agreement with the observed bulk rotation of the sunspot for all of the cases, except for the one which contains two sunspots in a shared penumbra. The method is also tested against sunspot rotations in active region (AR) 11158 that have been reported in the literature. While the results compare to some degree, the method outlined in this article reports lower rotations than those reported in the literature. Some of this discrepancy can be attributed to selection bias by the approaches in the literature, where only features that undergo larger rotation are tracked in sunspots that exhibit non-uniform rotation. The method also provides uncertainties on the calculated rotation profile which can be broken down to allow the principal sources of error to be identified. For the test sunspots in this article, the dominant source of uncertainty is the resolution of the SDO/HMI instrument.

Racial/Ethnic Heterogeneity and Rural-Urban Disparity of COVID-19 Case Fatality Ratio in the USA: a Negative Binomial and GIS-Based Analysis.

The 2019 coronavirus disease (COVID-19) has exacerbated inequality in the United States of America (USA). Black, indigenous, and people of color (BIPOC) are disproportionately affected by the pandemic. This study examines determinants of COVID-19 case fatality ratio (CFR) based on publicly sourced data from January 1 to December 18, 2020, and sociodemographic and rural-urban continuum data from the US Census Bureau. Nonspatial negative binomial Poisson regression and geographically weighted Poisson regression were applied to estimate the global and local relationships between the CFR and predictors—rural-urban continuum, political inclination, and race/ethnicity in 2407 rural counties. The mean COVID-19 CFR among rural counties was 1.79 (standard deviation (SD) = 1.07; 95% CI 1.73-1.84) higher than the total US counties (M = 1.69, SD = 1.18; 95% CI: 1.65-1.73). Based on the global NB model, CFR was positively associated with counties classified as “completely rural” (incidence rate ratio (IRR) = 1.24; 95% CI: 1.12-1.39) and “mostly rural” (IRR = 1.26; 95% CI: 1.15-1.38) relative to “mostly urban” counties. Nonspatial regression indicates that COVID-19 CFR increases by a factor of 8.62, 5.87, 2.61, and 1.36 for one unit increase in county-level percent Blacks, Hispanics, American Indians, and Asian/Pacific Islanders, respectively. Local spatial regression shows CFR was significantly higher in rural counties with a higher share of BIPOC in the Northeast and Midwest regions, and political inclination predicted COVID-19 CFR in rural counties in the Midwest region. In conclusion, spatial and racial/ethnic disparities exist for COVID-19 CFR across the US rural counties, and findings from this study have implications for public health.Supplementary InformationThe online version contains supplementary material available at 10.1007/s40615-021-01006-7.

A proto-pseudobulge in ESO 320-G030 fed by a massive molecular inflow driven by a nuclear bar

Galaxies with nuclear bars are believed to efficiently drive gas inward, generating a nuclear starburst and possibly an active galactic nucleus. We confirm this scenario for the isolated, double-barred, luminous infrared galaxy ESO 320-G030 based on an analysis of Herschel and ALMA spectroscopic observations. Herschel/PACS and SPIRE observations of ESO 320-G030 show absorption or emission in 18 lines of H2O, which we combine with the ALMA H2O 423 − 330 448 GHz line (Eupper ∼ 400 K) and continuum images to study the physical properties of the nuclear region. Radiative transfer models indicate that three nuclear components are required to account for the multi-transition H2O and continuum data. An envelope, with radius R ∼ 130 − 150 pc, dust temperature Tdust ≈ 50 K, and NH2 ∼ 2 × 1023 cm−2, surrounds a nuclear disk with R ∼ 40 pc that is optically thick in the far-infrared (τ100 μm ∼ 1.5 − 3, NH2 ∼ 2 × 1024 cm−2). In addition, an extremely compact (R ∼ 12 pc), warm (≈100 K), and buried (τ100 μm > 5, NH2 ≳ 5 × 1024 cm−2) core component is required to account for the very high-lying H2O absorption lines. The three nuclear components account for 70% of the galaxy luminosity (SFR ∼ 16 − 18 M⊙ yr−1). The nucleus is fed by a molecular inflow observed in CO 2-1 with ALMA, which is associated with the nuclear bar. With decreasing radius (r = 450 − 225 pc), the mass inflow rate increases up to Ṁinf ∼ 20 Ṁ yr−1, which is similar to the nuclear star formation rate (SFR), indicating that the starburst is sustained by the inflow. At lower r, ∼100 − 150 pc, the inflow is best probed by the far-infrared OH ground-state doublets, with an estimated Ṁinf ∼ 30 Ṁ yr−1. The inferred short timescale of ∼20 Myr for nuclear gas replenishment indicates quick secular evolution, and indicates that we are witnessing an intermediate stage (< 100 Myr) proto-pseudobulge fed by a massive inflow that is driven by a strong nuclear bar. We also apply the H2O model to the Herschel far-infrared spectroscopic observations of H218O, OH, 18OH, OH+, H2O+, H3O+, NH, NH2, NH3, CH, CH+, 13CH+, HF, SH, and C3, and we estimate their abundances.

Where’s the Dust?: The Deepening Anomaly of Microwave Emission in NGC 4725 B

Abstract We present new Atacama Large Millimeter Array (ALMA) observations toward NGC 4725 B, a discrete, compact, optically faint region within the star-forming disk of the nearby galaxy NGC 4725 that exhibits strong anomalous microwave emission (AME). These new ALMA data include continuum observations centered at 92, 133, 203, and 221 GHz accompanied by spectral observations of the 12CO ( ) line. NGC 4725 B is detected in the continuum at all frequencies, although the detection at 203 GHz is marginal. While molecular gas is not detected at the exact location of NGC 4725 B, there is molecular gas in the immediate vicinity (i.e., ≲100 pc) along with associated diffuse 8 μm emission. When combined with existing Very Large Array continuum data at 1.5, 3, 5.5, 9, 14, 22, 33, and 44 GHz, the spectrum is best fit by a combination of AME, synchrotron, and free–free emission that is free–free absorbed below ∼6 GHz. Given the strength of the AME, there is surprisingly no indication of millimeter dust emission associated with NGC 4725 B on ≲6″ spatial scales at the sensitivity of the ALMA interferometric data. Based on the properties of the nearest molecular gas complex and the inferred star formation rate, NGC 4725 B is consistent with being an extremely young (∼3–5 Myr) massive (≲105 M ⊙) cluster that is undergoing active cluster feedback. However, the lack of millimeter thermal dust emission is difficult to reconcile with a spinning dust origin of the 30 GHz emission. On the other hand, modeling NGC 4725 B as a new class of background radio galaxy is also unsatisfactory.

Dust dynamics in planet-driven spirals

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

Ammonia observations towards the Aquila Rift cloud complex

We surveyed the Aquila Rift complex including the Serpens South and W 40 regions in the NH3 (1,1) and (2,2) transitions making use of the Nanshan 26-m telescope. Our observations cover an area of ~ 1.5° × 2.2° (11.4 pc × 16.7 pc). The kinetic temperatures of the dense gas in the Aquila Rift complex obtained from NH3 (2,2)/(1,1) ratios range from 8.9 to 35.0 K with an average of 15.3 ± 6.1 K (errors are standard deviations of the mean). Low gas temperatures are associated with Serpens South ranging from 8.9 to 16.8 K with an average of 12.3 ± 1.7 K, while dense gas in the W 40 region shows higher temperatures ranging from 17.7 to 35.0 K with an average of 25.1 ± 4.9 K. A comparison of kinetic temperatures derived from para-NH3 (2,2)/(1,1) against HiGal dust temperatures indicates that the gas and dust temperatures are in agreement in the low-mass-star formation region of Serpens South. In the high-mass-star formation region W 40, the measured gas kinetic temperatures are higher than those of the dust. The turbulent component of the velocity dispersion of NH3 (1,1) is found to be positively correlated with the gas kinetic temperature, which indicates that the dense gas may be heated by dissipation of turbulent energy. For the fractional total-NH3 (para+ortho) abundance obtained by a comparison with Herschel infrared continuum data representing dust emission, we find values from 0.1 ×10−8 to 2.1 ×10−7 with an average of 6.9 (±4.5) × 10−8. Serpens South also shows a fractional total-NH3 (para+ortho) abundance ranging from 0.2 ×10−8 to 2.1 ×10−7 with an average of 8.6 (±3.8) × 10−8. In W 40, values are lower, between 0.1 and 4.3 ×10−8 with an average of 1.6 (±1.4) × 10−8. Weak velocity gradients demonstrate that the rotational energy is a negligible fraction of the gravitational energy. In W 40, gas and dust temperatures are not strongly dependent on the projected distance to the recently formed massive stars. Overall, the morphology of the mapped region is ring-like, with strong emission at lower and weak emission at higher Galactic longitudes. However, the presence of a physical connection between the two parts remains questionable.

The ALPINE-ALMA [CII] survey

We present ALMA observations of a merging system atz ∼ 4.57, observed as a part of the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) survey. Combining ALMA [CII]158 μm and far-infrared continuum data with multi-wavelength ancillary data, we find that the system is composed of two massive (M⋆ ≳ 1010 M⊙) star-forming galaxies experiencing a major merger (stellar mass ratiormass ≳ 0.9) at close spatial (∼13 kpc; projected) and velocity (Δv < 300 km s−1) separations, and two additional faint narrow [CII]-emitting satellites. The overall system belongs to a larger scale protocluster environment and is coincident to one of its overdensity peaks. Additionally, ALMA reveals the presence of [CII] emission arising from a circumgalactic gas structure, extending up to a diameter-scale of ∼30 kpc. Our morpho-spectral decomposition analysis shows that about 50% of the total flux resides between the individual galaxy components, in a metal-enriched gaseous envelope characterised by a disturbed morphology and complex kinematics. Similarly to observations of shock-excited [CII] emitted from tidal tails in local groups, our results can be interpreted as a possible signature of interstellar gas stripped by strong gravitational interactions, with a possible contribution from material ejected by galactic outflows and emission triggered by star formation in small faint satellites. Our findings suggest that mergers could be an efficient mechanism of gas mixing in the circumgalactic medium around high-zgalaxies, and thus play a key role in the galaxy baryon cycle at early epochs.

The ALPINE-ALMA [CII] survey: Data processing, catalogs, and statistical source properties

The Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) targets the [CII] 158 μm line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z = 4.4 and z = 5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details regarding the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates (SFRs), we measured the conversion factor from the ALMA 158 μm rest-frame dust continuum luminosity to the total infrared luminosity (LIR) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median LIR of 4.4 × 1011 L⊙. We also detected 57 additional continuum sources in our ALMA pointings. They are at a lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5 ± 0.2. We measured the 850 μm number counts between 0.35 and 3.5 mJy, thus improving the current interferometric constraints in this flux density range. We found a slope break in the number counts around 3 mJy with a shallower slope below this value. More than 40% of the cosmic infrared background is emitted by sources brighter than 0.35 mJy. Finally, we detected the [CII] line in 75 of our targets. Their median [CII] luminosity is 4.8 × 108 L⊙ and their median full width at half maximum is 252 km s−1. After measuring the mean obscured SFR in various [CII] luminosity bins by stacking ALPINE continuum data, we find a good agreement between our data and the local and predicted SFR–L[CII] relations.