We present the Pilot Survey Phase 2 data release for the Wide-field ASKAP L-band Legacy All-sky Blind surveY (WALLABY), carried-out using the Australian SKA Pathfinder (ASKAP). We present 1760 H i detections (with a default spatial resolution of 30′′) from three pilot fields including the NGC 5044 and NGC 4808 groups as well as the Vela field, covering a total of $\sim 180$ deg$^2$ of the sky and spanning a redshift up to $z \simeq 0.09$. This release also includes kinematic models for over 126 spatially resolved galaxies. The observed median rms noise in the image cubes is 1.7 mJy per 30′′ beam and 18.5 kHz channel. This corresponds to a 5$\sigma$ H i column density sensitivity of $\sim 9.1\times10^{19}(1 + z)^4$ cm$^{-2}$ per 30′′ beam and $\sim 20$ km s$^{-1}$ channel and a 5$\sigma$ H i mass sensitivity of $\sim 5.5\times10^8 (D/100$ Mpc)$^{2}$ M$_{\odot}$ for point sources. Furthermore, we also present for the first time 12′′ high-resolution images (“cut-outs”) and catalogues for a sub-sample of 80 sources from the Pilot Survey Phase 2 fields. While we are able to recover sources with lower signal-to-noise ratio compared to sources in the Public Data Release 1, we do note that some data quality issues still persist, notably, flux discrepancies that are linked to the impact of side lobes associated with the dirty beams due to inadequate deconvolution. However, in spite of these limitations, the WALLABY Pilot Survey Phase 2 has already produced roughly a third of the number of HIPASS sources, making this the largest spatially resolved H i sample from a single survey to date.

We investigated the kinematics and dynamics of gas structures on galaxy-cloud scales in two spiral galaxies NGC5236 (M83) and NGC4321 (M100) using CO (2$-$1) line. We utilised the FILFINDER algorithm on integrated intensity maps for the identification of filaments in two galaxies. Clear fluctuations in velocity and density were observed along these filaments, enabling the fitting of velocity gradients around intensity peaks. The variations in velocity gradient across different scales suggest a gradual and consistent increase in velocity gradient from large to small scales, indicative of gravitational collapse, something also revealed by the correlation between velocity dispersion and column density of gas structures. Gas structures at different scales in the galaxy may be organised into hierarchical systems through gravitational coupling. All the features of gas kinematics on galaxy-cloud scale are very similar to that on cloud-clump and clump-core scales studied in previous works. Thus, the interstellar medium from galaxy to dense core scales presents multi-scale/hierarchical hub-filament structures. Like dense core as the hub in clump, clump as the hub in molecular cloud, now we verify that cloud or cloud complex can be the hub in spiral galaxies. Although the scaling relations and the measured velocity gradients support the gravitational collapse of gas structures on galaxy-cloud scales, the collapse is much slower than a pure free-fall gravitational collapse.

After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass loading factor for every galaxy. Building on the model’s resolved angular momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus–age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z = 0$ stellar and H i mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and python, is publicly available at https://github.com/arhstevens/DarkSage.

Cold, neutral interstellar gas, the reservoir for star formation, is traced through the absorption of the 21-cm continuum radiation by neutral hydrogen (H i). Although detected in one hundred cases in the host galaxies of distant radio sources, only recently have column densities approaching the maximum value observed in Lyman-$\alpha$ absorption systems ($N_{{\textrm{H}\,\scriptsize{\textrm{I}}}}\sim 10^{22}$ $\textrm{cm}^{-2}$) been found. Here, we explore the implications these have for the hypothesis that the detection rate of H i absorption is dominated by photo-ionisation from the active galactic nucleus (AGN). We find, with the addition all of the current searches for H i absorption at $z\geq0.1$, a strong correlation between the H i absorption strength and the ionising photon rate, with the maximum value at which H i is detected remaining close to the theoretical value in which all of the neutral gas would be ionised in a large spiral galaxy ($Q_{{\textrm{H}\,\scriptsize{\textrm{I}}}} = 2.9\times10^{56}$ ionising photons s$^{-1}$). We also rule out other effects (excitation by the radio continuum and changing gas properties) as the dominant cause for the decrease in the detection rate with redshift. Furthermore, from the maximum theoretical column density, we find that the five high column density systems have spin temperatures close to those of the Milky Way ($T_{\textrm{spin}} \lesssim 300$ K), whereas, from our model of a gaseous galactic disc, the H i detection at $Q_{{\textrm{H}\,\scriptsize{\textrm{I}}}} =2.9\times10^{56}$ s$^{-1}$ yields $T_{\textrm{spin}}\sim10\,000$ K, consistent with the gas being highly ionised.

We here present 0.02–0.04″ resolution ALMA observation of the compact obscured nucleus (CON) of IRAS 17578-0400. A dusty torus within the nucleus, approximately 4 pc in radius, has been uncovered, exhibiting a usually flat spectral index at ALMA band 3, likely due to the millimeter corona emission from the central supermassive black hole (SMBH). The dense gas disk, traced by 13CO(1-0), spans 7 pc in radius and suggests an outflow driven by a disk wind due to its asymmetrical structure along the minor axis. Collimated molecular outflows (CMO), traced by the low-velocity components of the HCN(3-2) and HCO+(3-2) lines, align with the minor axis gas disk. Examination of position-velocity plots of HCN(3-2) and HCO+(3-2) reveals a flared dense gas disk extended a radius of ∼60 pc, infalling and rotating at speeds of about 200 km s−1 and 300 km s−1 respectively. A centrifugal barrier, located around 4 pc from the dynamical center, implies an SMBH mass of approximately 108 Mȯ, consistent with millimeter corona emission estimates. The CMO maintains a steady rotation speed of 200 km s−1 over the 100 pc scale along the minor axis. The projected speed of the CMO is about 80 km s−1, corresponding to around ∼500 km s−1, assuming an inclination angle of 80°. Such a kinematics structure of disk-driven collimated rotating molecular outflow with gas supplies from a falling rotating disk indicates that the feedback of the compact obscured nucleus is likely regulated by the momentum transfer of the molecular gas that connects to both the feeding of the nuclear starburst and supermassive black hole.

Galaxy gas kinematics are sensitive to the physical processes that contribute to a galaxy’s evolution. It is expected that external processes will cause more significant kinematic disturbances in the outer regions, while internal processes will cause more disturbances for the inner regions. Using a subsample of 47 galaxies ($0.27<z<0.36$) from the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, we conduct a study into the source of kinematic disturbances by measuring the asymmetry present in the ionised gas line-of-sight velocity maps at the $0.5R_e$ (inner regions) and $1.5R_e$ (outer regions) elliptical annuli. By comparing the inner and outer kinematic asymmetries, we aim to better understand what physical processes are driving the asymmetries in galaxies. We find the local environment plays a role in kinematic disturbance, in agreement with other integral field spectroscopy studies of the local universe, with most asymmetric systems being in close proximity to a more massive neighbour. We do not find evidence suggesting that hosting an Active Galactic Nucleus contributes to asymmetry within the inner regions, with some caveats due to emission line modelling. In contrast to previous studies, we do not find evidence that processes leading to asymmetry also enhance star formation in MAGPI galaxies. Finally, we find a weak anti-correlation between stellar mass and asymmetry (i.e., high stellar mass galaxies are less asymmetric). We conclude by discussing possible sources driving the asymmetry in the ionised gas, such as disturbances being present in the colder gas phase (either molecular or atomic) prior to the gas being ionised, and non-axisymmetric features (e.g., a bar) being present in the galactic disk. Our results highlight the complex interplay between ionised gas kinematic disturbances and physical processes involved in galaxy evolution.

We investigate the diversity in the sizes and average surface densities of the neutral atomic hydrogen (H i) gas discs in $\sim$280 nearby galaxies detected by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). We combine the uniformly observed, interferometric H i data from pilot observations of the Hydra cluster and NGC 4636 group fields with photometry measured from ultraviolet, optical, and near-infrared imaging surveys to investigate the interplay between stellar structure, star formation, and H i structural parameters. We quantify the H i structure by the size of the H i relative to the optical disc and the average H i surface density measured using effective and isodensity radii. For galaxies resolved by $>$$1.3$ beams, we find that galaxies with higher stellar masses and stellar surface densities tend to have less extended H i discs and lower H i surface densities: the isodensity H i structural parameters show a weak negative dependence on stellar mass and stellar mass surface density. These trends strengthen when we limit our sample to galaxies resolved by $>$2 beams. We find that galaxies with higher H i surface densities and more extended H i discs tend to be more star forming: the isodensity H i structural parameters have stronger correlations with star formation. Normalising the H i disc size by the optical effective radius (instead of the isophotal radius) produces positive correlations with stellar masses and stellar surface densities and removes the correlations with star formation. This is due to the effective and isodensity H i radii increasing with mass at similar rates while, in the optical, the effective radius increases slower than the isophotal radius. Our results are in qualitative agreement with previous studies and demonstrate that with WALLABY we can begin to bridge the gap between small galaxy samples with high spatial resolution H i data and large, statistical studies using spatially unresolved, single-dish data.

The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (H i) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival H i data, we study the H i and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($\Sigma_{\star}$), H i ($\Sigma_{\text{H}\,{\small\text{I}}}$), and molecular gas ($\Sigma_{\text{mol}}$) surface densities. Adopting H i deficiency as a measure of environmental impact, we find evidence that, in addition to removing the H i at large radii, the cluster processes also lower the average $\Sigma_{\text{H}\,{\small\text{I}}}$ of the remaining gas even in the central $1.2\,$kpc. The impact on molecular gas is comparatively weaker than on the H i, and we show that the lower $\Sigma_{\text{mol}}$ gas is removed first. In the most H i-deficient galaxies, however, we find evidence that environmental processes reduce the typical $\Sigma_{\text{mol}}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $\Sigma_{\text{H}\,{\small\text{I}}}$ or $\Sigma_{\text{mol}}$ in H i-deficient galaxies. Using the ratio of $\Sigma_{\text{mol}}$-to-$\Sigma_{\text{H}\,{\small\text{I}}}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results.

We present a new high-resolution neutral atomic hydrogen (Hi) survey of ring galaxies using the Australia Telescope Compact Array (ATCA). We target a sample of 24 ring galaxies from the Buta (1995) Southern Ring Galaxy Survey Catalogue in order to study the origin of resonance-, collisional- and interaction-driven ring galaxies. In this work, we present an overview of the sample and study their global and resolved Hi properties. In addition, we also probe their star formation properties by measuring their star formation rates (SFR) and their resolved SFR surface density profiles. We find that a majority of the barred galaxies in our sample are Hi-deficient, alluding to the effects of the bar in driving their Hi deficiency. Furthermore, for the secularly evolving barred ring galaxies in our sample, we apply Lindblad’s resonance theory to predict the location of the resonance rings and find very good agreement between predictions and observations. We identify rings of Hi gas and/or star formation co-located at one or the other major resonances. Lastly, we measure the bar pattern speed ($\Omega_{\textrm{bar}}$) for a sub-sample of our galaxies and find that the values range from 10–90 $\textrm{km s}^{-1}$kpc$^{-1}$, in good agreement with previous studies.

We present observations of the four $^2 \Pi _{3/2}\,J=3/2$ ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitation temperature ($T_{\rm ex}$) and optical depth, while our ATCA observations give optical depth only. We perform Gaussian decomposition using the Automated Molecular Excitation Bayesian line-fitting Algorithm ‘Amoeba’ (Petzler, Dawson, & Wardle 2021, ApJ, 923, 261) fitting all four transitions simultaneously with shared centroid velocity and width. We identify 109 features across 38 sightlines (including 58 detections along 27 sightlines with excitation temperature measurements). While the main lines at 1665 and 1667 MHz tend to have similar excitation temperatures (median $|\Delta T_{\rm ex}({\rm main})|=0.6\,$K, 84% show $|\Delta T_{\rm ex}({\rm main})|<2\,$K), large differences in the 1612 and 1720 MHz satellite line excitation temperatures show that the gas is generally not in LTE. For a selection of sightlines, we compare our OH features to associated (on-sky and in velocity) Hi cold gas components (CNM) identified by Nguyen et al. (2019, ApJ, 880, 141) and find no strong correlations. We speculate that this may indicate an effective decoupling of the molecular gas from the CNM once it accumulates.

Understanding properties of galaxies in the epoch of reionization (EoR) is a frontier in the modern astronomy. ALMA observations have demonstrated that i) some [O iii] 88 μm emitters have matured stellar populations at z>6, implying early star formation activity at z>10, and that ii) high-z star-forming galaxies typically have very high [O iii] 88 μm-to-[C ii] 158 μm luminosity ratios ranging from 3 to 12 or higher, indicating interstellar media of high-z galaxies could be highly ionized. We discuss initial results of a medium-sized JWST GO1 program that targets a sample of 12 z=6–8 ALMA [O iii] 88 μm emitters with NIRCam and NIRSPec IFU modes (GO-1840). Our JWST GO1 program, in conjunction with ALMA data, will characterize the stellar, nebular, and dust properties of these [O iii] 88 μm emitters and place this galaxies in the context of reionization.

Investigating the thermal and non-thermal processes in galaxies is vital to understand their evolution over cosmic time. This can best be studied by combining the radio and optical/near-infrared observations of galaxies. The JWST can resolve the evolution of the thermal processes by mapping ionized gas and dust in distant galaxies. This information combined with the upcoming surveys with the Square Kilometer Array (SKA) will make a major breakthrough in mapping the non-thermal processes and understanding their role in the evolution of galaxies. Our simulations show that SKA surveys will be able to trace the evolution history of spiral galaxies such as M 51 and NGC 6946 back to a redshift of 3 already in its first phase of construction. This study indicates the important role of the non-thermal pressure inserted by cosmic rays and magnetic fields in deriving winds and outflows at cosmic noon as deduced by a flat synchrotron spectrum in star forming galaxies.

The shape of emission lines in the optical spectra of star-forming galaxies reveals the kinematics of the diffuse gaseous component. We analyse the shape of prominent emission lines in a sample of $\sim$ 53000 star-forming galaxies from the Sloan Digital Sky Survey, focusing on departures from gaussianity. Departures from a single gaussian profile allow us to probe the motion of gas and to assess the role of outflows. The sample is divided into groups according to their stellar velocity dispersion and star formation rate (SFR). The spectra within each group are stacked to improve the signal-to-noise ratio of the emission lines, to remove individual signatures, and to enhance the effect of SFR on the shapes of the emission lines. The moments of the emission lines, including kurtosis and skewness, are determined. We find that most of the emission lines in strong star-forming systems unequivocally feature negative kurtosis. This signature is present in $\mathrm{H}\unicode{x03B2}$ , $\mathrm{H}\unicode{x03B1}$ , [N ii], and [S ii] in massive galaxies with high SFRs. We attribute it as evidence of radial outflows of ionised gas driven by the star formation of the galaxies. Also, most of the emission lines in low-mass systems with high SFRs feature negative skewness, and we interpret it as evidence of dust obscuration in the galactic disk. These signatures are however absent in the [O iii] line, which is believed to trace a different gas component. The observed trend is significantly stronger in face-on galaxies, indicating that star formation drives the outflows along the galactic rotation axis, presumably the path of least resistance. The data suggest that outflows driven by star formation exert accumulated impacts on the interstellar medium, and the outflow signature is more evident in older galaxies as they have experienced a longer total duration of star formation.

Nebular Heii emission implies the presence of energetic photons (E≽54 eV). Despite the great deal of effort dedicated to understanding Heii ionization, its origin has remained mysterious, particularly in metal-deficient star-forming galaxies. Unfolding Heii-emitting, metal-poor starbursts at z∼0 can yield insight into the powerful ionization processes occurring in the primordial universe. Here we present a study on the origin of the extended nebular Heii emission in SBS 0335-052E, one of the most metal-poor (Z ∼ 3% Z⊙ Heii-emitter starbursts known locally. Based on optical VLT/MUSE spectroscopic and Chandra X-ray observations, and current stellar models we found that the Heii-ionization budget of SBS 0335-052E can only be produced by peculiar, nearly metal-free ionizing stars (called here “PopIII-like” stars) with a top-heavy initial mass function. This result is in line with recent simulations for PopIII star formation down to z=0.

We describe the scientific goals and survey design of the First Large Absorption Survey in H i (FLASH), a wide field survey for 21-cm line absorption in neutral atomic hydrogen (H i) at intermediate cosmological redshifts. FLASH will be carried out with the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope and is planned to cover the sky south of $\delta \approx +40\,\deg$ at frequencies between 711.5 and 999.5 MHz. At redshifts between $z = 0.4$ and $1.0$ (look-back times of 4 – 8 Gyr), the H i content of the Universe has been poorly explored due to the difficulty of carrying out radio surveys for faint 21-cm line emission and, at ultra-violet wavelengths, space-borne searches for Damped Lyman- $\alpha$ absorption in quasar spectra. The ASKAP wide field of view and large spectral bandwidth, in combination with a radio-quiet site, will enable a search for absorption lines in the radio spectra of bright continuum sources over 80% of the sky. This survey is expected to detect at least several hundred intervening 21-cm absorbers and will produce an H i-absorption-selected catalogue of galaxies rich in cool, star-forming gas, some of which may be concealed from optical surveys. Likewise, at least several hundred associated 21-cm absorbers are expected to be detected within the host galaxies of radio sources at $0.4 < z < 1.0$ , providing valuable kinematical information for models of gas accretion and jet-driven feedback in radio-loud active galactic nuclei. FLASH will also detect OH 18-cm absorbers in diffuse molecular gas, megamaser OH emission, radio recombination lines, and stacked H i emission.

Even though galactic winds are common in galaxies with starbursts or active galactic nuclei (AGN), the role of such gas flows in galaxy evolution remains uncertain. Here we examine how winds vary along a likely evolutionary sequence connecting starburst to post-starburst to quiescent galaxies. To detect the interstellar medium and measure its bulk flows, we examine the residual Na D absorption line doublet after the stellar contribution has been removed from each galaxy’s spectrum. We discover that outflows diminish along this sequence, i.e., as star formation ends. We then focus on the wind behavior within the post-starburst sample, for which we have measured the time elapsed since the starburst ended (post-burst age) via detailed modeling of their star formation histories (French et al.2018). Even within our post-starburst sample, the fraction of galaxies with significant winds and the average wind velocities decrease with post-burst age after controlling for stellar mass.

The Atacama Large Millimetre/Sub-millimetre Array (ALMA) is obtaining the deepest observations of early galaxies ever achieved at (sub-)millimetre wavelengths, and detecting the dust emission of young galaxies in the first billion years of cosmic history, well in the epoch of reionization. Here I review some of the latest results from these observations, with special focus on the REBELS large programme, which targets a sample of 40 star-forming galaxies at z ⋍ 7. ALMA detects significant amounts of dust in very young galaxies, and this dust might have different properties to dust in lower-redshift galaxies. I describe the evidence for this, and discuss theoretical/modelling efforts to explain the dust properties of these young galaxies. Finally, I describe two additional surprising results to come out of the REBELS survey: (i) a new population of completely dust-obscured galaxies at z ⋍ 7, and (ii) the prevalence of spatial offsets between the ultraviolet and infrared emission of UV-bright, high-redshift star-forming galaxies.

Star-formation is one of the main processes that shape galaxies, defining its stellar population and metallicity production and enrichment. It is nowadays known that this process is ruled by a set of relations that connect three parameters: the molecular gas mass, the stellar mass and the star-formation rate itself. These relations are fulfilled at a wide range of scales in galaxies, from galaxy wide to kpc-scales. At which scales they are broken, and how universal they are (i.e., if they change at different scales or for different galaxy types) it is still an open question. We explore here how those relations compare at different scales using as proxy the new analysis done using Integral Field Spectroscopy data and CO observations data from the EDGE-CALIFA survey and the AMUSSING++ compilation.

Remarkable progress has been made in the last few years in understanding the global properties of galaxies and how they evolve through cosmic time. Major focus has been given to studies of how the availability of molecular gas regulates star-forming activity and galaxy growth, the eventual quenching of star formation, and how these mechanisms evolve through cosmic time. Most of these advances have been made thanks to ALMA and the upgraded capabilities of NOEMA. In this contribution, I briey review the latest constraints on the molecular gas content based on dierent tracers of the interstellar medium (ISM; dust continuum and CO, [CI] and [CII] line emission), including recent determinations of the molecular gas fraction, gas depletion timescales, and molecular gas cosmic density provided by the recent ALMA programs out to z ∼ 7. Finally, I concentrate on recent and ongoing studies aiming to spatially and kinematically resolve the cold ISM and star formation activity down to kpc scales in galaxies out to z ∼ 6 – 7, which represent an unprecedented view of the galaxy assembly and feedback processes in the early universe.

Whether the star formation efficiency (SFE) in the bar region is lower than those in the other regions in a barred galaxy has recently been debated. We statistically investigate the SFEs along the bars in nearby gas-rich massive star-forming barred galaxies by distinguishing the center, bar-end, and bar regions for the first time. The molecular gas surface density is derived from archival CO(1–0) and/or CO(2–1) data and the star formation rate surface density is derived from a linear combination of far-ultraviolet and mid-infrared intensities. To distinguish the three regions, we targeted 18 galaxies with a large apparent bar length (≥ 75"). The resulting SFE in the bars is about 0.6 – 0.8 times lower than that in the disks, which suggests the star formation in the bars tends to be systematically suppressed.