The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

We present a Wide-field ASKAP L-Band Legacy All-sky Blind surveY (WALLABY) study of the nearby ($v_{\rm sys}$ = 915 km s$^{-1}$) spiral galaxy IC 5201 using the Australian Square Kilometre Array Pathfinder (ASKAP). IC 5201 is a blue, barred spiral galaxy that follows the known scaling relations between stellar mass, SFR, HI mass and diameter. We create a four-beam mosaicked HI image cube, from 175 hours of observations made with a 12-antenna sub-array. The RMS noise level of the cube is 1.7 mJy beam$^{-1}$ per channel, equivalent to a column density of $N_{\rm HI}$ = 1.4 $\times$ 10$^{20}$ cm$^{-2}$ over 25 km s$^{-1}$. We report 9 extragalactic HI detections $-$ 5 new HI detections including the first velocity measurements for 2 galaxies. These sources are IC 5201, 3 dwarf satellite galaxies, 2 galaxies and a tidal feature belonging to the NGC 7232/3 triplet and 2 potential infalling galaxies to the triplet. There is evidence of a previous tidal interaction between IC 5201 and the irregular satellite AM 2220$-$460. A close fly-by is likely responsible for the asymmetric optical morphology of IC 5201 and warping its disc, resulting in the irregular morphology of AM 2220$-$460. We quantify the HI kinematics of IC 5201, presenting its rotation curve as well as showing that the warp starts at 14 kpc along the major axis, increasing as a function of radius with a maximum difference in position angle of 20$^\circ$. There is no evidence of stripped HI, triggered or quenched star formation in the system as measured using DECam optical and $GALEX$ UV photometry.

This paper reports on the atomic hydrogen gas (HI) observations of the spiral galaxy NGC 1566 using the newly commissioned Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope. We measure an integrated HI flux density of $180.2$ Jy km s$^{-1}$ emanating from this galaxy, which translates to an HI mass of $1.94\times10^{10}$M$_\circ$ at an assumed distance of $21.3$ Mpc. Our observations show that NGC 1566 has an asymmetric and mildly warped HI disc. The HI-to-stellar mass fraction of NGC 1566 is $0.29$, which is high in comparison with galaxies that have the same stellar mass ($10^{10.8}$M$_\circ$). We also derive the rotation curve of this galaxy to a radius of $50$ kpc and fit different mass models to it. The NFW, Burkert and pseudo-isothermal dark matter halo profiles fit the observed rotation curve reasonably well and recover dark matter fractions of $0.62$, $0.58$ and $0.66$, respectively. Down to the column density sensitivity of our observations ($N_{HI} = 3.7\times10^{19}$ cm$^{-2}$), we detect no HI clouds connected to, or in the nearby vicinity of, the HI disc of NGC 1566 nor nearby interacting systems. We conclude that, based on a simple analytic model, ram pressure interactions with the IGM can affect the HI disc of NGC 1566 and is possibly the reason for the asymmetries seen in the HI morphology of NGC 1566.

We present Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY) early science results from the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the NGC 7162 galaxy group. We use archival HIPASS and Australia Telescope Compact Array (ATCA) observations of this group to validate the new ASKAP data and the data reduction pipeline ASKAPsoft. We detect six galaxies in the neutral hydrogen (HI) 21-cm line, expanding the NGC 7162 group membership from four to seven galaxies. Two of the new detections are also the first HI detections of the dwarf galaxies, AM 2159-434 and GALEXASC J220338.65-431128.7, for which we have measured velocities of $cz=2558$ and $cz=2727$ km s$^{-1}$, respectively. We confirm that there is extended HI emission around NGC 7162 possibly due to past interactions in the group as indicated by the $40^{\circ}$ offset between the kinematic and morphological major axes for NGC 7162A, and its HI richness. Taking advantage of the increased resolution (factor of $\sim1.5$) of the ASKAP data over archival ATCA observations, we fit a tilted ring model and use envelope tracing to determine the galaxies' rotation curves. Using these we estimate the dynamical masses and find, as expected, high dark matter fractions of $f_{\mathrm{DM}}\sim0.81-0.95$ for all group members. The ASKAP data are publicly available.

We obtained ALMA spectroscopy and imaging to investigate the origin of the unexpected sub-mm emission toward the most distant quiescent galaxy known to date, ZF-COSMOS-20115 at z=3.717. We show here that this sub-mm emission is produced by another massive, compact and extremely obscured galaxy, located only 3.1 kpc away from the quiescent galaxy. We dub the quiescent and dusty galaxies Jekyll and Hyde, respectively. No dust emission is detected at the location of Jekyll, implying SFR < 13 Msun/yr, which is the most stringent upper limit ever obtained for a quiescent galaxy at these redshifts. The two sources are confirmed to lie at the same redshift thanks to the detection of [CII]158 in Hyde, which provides one the few robust redshifts for an "H-dropout" galaxy. The line has a rotating-disk velocity profile blueshifted from Jekyll by 549+/-60 km/s, demonstrating that it is produced by another galaxy. Careful de-blending of the Spitzer imaging confirms the existence of Hyde, and its non-detection with Hubble requires extreme attenuation by dust. Modeling the photometry of both galaxies shows that Jekyll has fully quenched >200 Myr prior to observation and still presents a challenge for models, while Hyde only harbors moderate star-formation (SFR<120 Msun/yr) and is located at least a factor 1.4 below the z~4 main sequence. Hyde could also have stopped forming stars <200 Myr before being observed, which would be consistent with its hight compactness similar to z~4 quiescent galaxies and its low [CII]/FIR ratio, but significant SF cannot be ruled out. Finally, we show that Hyde hosts a dense reservoir of gas comparable to that of extreme starbursts, suggesting that its SFR was reduced without expelling the gas outside of the galaxy. We argue that Jekyll and Hyde can be seen as two stages of quenching, and provide a unique laboratory to study this poorly understood phenomenon. [abridged]

We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($\sigma=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$\alpha$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result.

In this letter we study the mean sizes of Halpha clumps in turbulent disk galaxies relative to kinematics, gas fractions, and Toomre Q. We use 100~pc resolution HST images, IFU kinematics, and gas fractions of a sample of rare, nearby turbulent disks with properties closely matched to z~1.5-2 main-sequence galaxies (the DYNAMO sample). We find linear correlations of normalized mean clump sizes with both the gas fraction and the velocity dispersion-to-rotation velocity ratio of the host galaxy. We show that these correlations are consistent with predictions derived from a model of instabilities in a self-gravitating disk (the so-called "violent disk instability model"). We also observe, using a two-fluid model for Q, a correlation between the size of clumps and self-gravity driven unstable regions. These results are most consistent with the hypothesis that massive star forming clumps in turbulent disks are the result of instabilities in self-gravitating gas-rich disks, and therefore provide a direct connection between resolved clump sizes and this in situ mechanism.

The gas accretion and star-formation histories of galaxies like the Milky Way remain an outstanding problem in astrophysics. Observations show that 8 billion years ago, the progenitors to Milky Way-mass galaxies were forming stars 30 times faster than today and predicted to be rich in molecular gas, in contrast with low present-day gas fractions ($<$10%). Here we show detections of molecular gas from the CO(J=3-2) emission (rest-frame 345.8 GHz) in galaxies at redshifts z=1.2-1.3, selected to have the stellar mass and star-formation rate of the progenitors of today's Milky Way-mass galaxies. The CO emission reveals large molecular gas masses, comparable to or exceeding the galaxy stellar masses, and implying most of the baryons are in cold gas, not stars. The galaxies' total luminosities from star formation and CO luminosities yield long gas-consumption timescales. Compared to local spiral galaxies, the star-formation efficiency, estimated from the ratio of total IR luminosity to CO emission,} has remained nearly constant since redshift z=1.2, despite the order of magnitude decrease in gas fraction, consistent with results for other galaxies at this epoch. Therefore the physical processes that determine the rate at which gas cools to form stars in distant galaxies appear to be similar to that in local galaxies.

With large-scale Integral Field Spectroscopy (IFS) surveys of thousands of galaxies currently under-way or planned, the astronomical community is in need of methods, techniques and tools that will allow the analysis of huge amounts of data. We focus on the kinematic modelling of disc galaxies and investigate the potential use of massively parallel architectures, such as the Graphics Processing Unit (GPU), as an accelerator for the computationally expensive model-fitting procedure. We review the algorithms involved in model-fitting and evaluate their suitability for GPU implementation. We employ different optimization techniques, including the Levenberg-Marquardt and Nested Sampling algorithms, but also a naive brute-force approach based on Nested Grids. We find that the GPU can accelerate the model-fitting procedure up to a factor of ~100 when compared to a single-threaded CPU, and up to a factor of ~10 when compared to a multi-threaded dual CPU configuration. Our method's accuracy, precision and robustness are assessed by successfully recovering the kinematic properties of simulated data, and also by verifying the kinematic modelling results of galaxies from the GHASP and DYNAMO surveys as found in the literature. The resulting GBKFIT code is available for download from: http://supercomputing.swin.edu.au/gbkfit.