ALMA Super-resolution Imaging of T Tau: r = 12 au Gap in the Compact Dust Disk around T Tau N

Based on Atacama Large Millimeter/submillimeter Array (ALMA) observations, compact protoplanetary disks with dust radii of r ≲ 20–40 au were found to be dominant in nearby low-mass star formation regions. However, their substructures have not been investigated because of the limited spatial resolution achieved so far. We apply a newly developed super-resolution imaging technique utilizing sparse modeling (SpM) to explore several au-scale structures in such compact disks. SpM imaging can directly solve for the incomplete sampling of visibilities in the spatial frequency and potentially improve the fidelity and effective spatial resolution of ALMA images. Here we present the results of the application to the T Tau system. We use the ALMA 1.3 mm continuum data and achieve an effective spatial resolution of ∼30% (5 au) compared with the conventional CLEAN beam size at a resolution of 17 au. The reconstructed image reveals a new annular gap structure at r = 12 au in the T Tau N compact disk, with a dust radius of 24 au, and resolves the T Tau Sa and Sb binary into two sources. If the observed gap structure in the T Tau N disk is caused by an embedded planet, we estimate a Saturn-mass planet when the viscous parameter of the disk is 10−3. Ultimately, ALMA observations with enough angular resolution and sensitivity should be able to verify the consistency of the super-resolution imaging and definitely confirm the existence of this disk substructure.

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

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

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.

Accurate dust temperature determination in a z = 7.13 galaxy

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 (FIR) 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} \, 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.

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.

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.

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.

The ALPINE-ALMA [CII] survey

We present ALMA observations of a merging system at z ∼ 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 ratio rmass ≳ 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-z galaxies, 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.