For studies of galaxy formation and evolution, one of the major benefits of the James Webb Space Telescope is that space-based IFUs like those on its NIRSpec and MIRI instruments will enable spatially resolved spectroscopy of distant galaxies, including spectroscopy at the scale of individual star-forming regions in galaxies that have been gravitationally lensed. In the meantime, there is only a very small subset of lensed sources where work like this is possible even with the Hubble Space Telescope's Wide Field Camera 3 infrared channel grisms. We examine two of these sources, SDSS J1723+3411 and SDSS J2340+2947, using HST WFC3/IR grism data and supporting spatially-unresolved spectroscopy from several ground-based instruments to explore the size of spatial variations in observed strong emission line ratios like O32, R23, which are sensitive to ionization parameter and metallicity, and the Balmer decrement as an indicator of reddening. We find significant spatial variation in the reddening and the reddening-corrected O32 and R23 values which correspond to spreads of a few tenths of a dex in ionization parameter and metallicity. We also find clear evidence of a negative radial gradient in star formation in SDSS J2340+2947 and tentative evidence of one in SDSS J1723+3411, though its star formation is quite asymmetric. Finally, we find that reddening can vary enough spatially to make spatially-resolved reddening corrections necessary in order to characterize gradients in line ratios and the physical conditions inferred from them, necessitating the use of space-based IFUs for future work on larger, more statistically robust samples.

Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements.

Our ability to study the properties of the interstellar medium (ISM) in the earliest galaxies will rely on emission line diagnostics at rest-frame ultraviolet (UV) wavelengths. In this work, we identify metallicity-sensitive diagnostics using UV emission lines. We compare UV-derived metallicities with standard, well-established optical metallicities using a sample of galaxies with rest-frame UV and optical spectroscopy. We find that the He2-O3C3 diagnostic (He II 1640 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, particularly at low metallicity (12+log(O/H) < 8), where stellar contributions are minimal. We find that the Si3-O3C3 diagnostic (Si III 1883 / C III 1906,1909 vs. O III 1666 / C III 1906,1909) is a reliable metallicity tracer, though with large scatter (0.2-0.3 dex), which we suggest is driven by variations in gas-phase abundances. We find that the C4-O3C3 diagnostic (C IV 1548,1550 / O III 1666 vs. O III 1666 / C III 1906,1909) correlates poorly with optically-derived metallicities. We discuss possible explanations for these discrepant metallicity determinations, including the hardness of the ionizing spectrum, contribution from stellar wind emission, and non-solar-scaled gas-phase abundances. Finally, we provide two new UV oxygen abundance diagnostics, calculated from polynomial fits to the model grid surface in the He2-O3C3 and Si3-O3C3 diagrams.

We derive new self-consistent theoretical UV, optical, and IR diagnostics for the ISM pressure and electron density in the ionized nebulae of star-forming galaxies. Our UV diagnostics utilize the inter-combination, forbidden and resonance lines of silicon, carbon, aluminum, neon, and nitrogen. We also calibrate the optical and IR forbidden lines of oxygen, argon, nitrogen and sulfur. We show that line ratios used as ISM pressure diagnostics depend on the gas-phase metallicity with a residual dependence on the ionization parameter of the gas. In addition, the traditional electron density diagnostic [S II] {\lambda}6731/[S II] {\lambda}6717 is strongly dependent on the gas-phase metallicity. We show how different emission-line ratios are produced in different ionization zones in our theoretical nebulae. The [S II] and [O II] ratios are produced in different zones, and should not be used interchangeably to measure the electron density of the gas unless the electron temperature is known to be constant. We review the temperature and density distributions observed within H II regions and discuss the implications of these distributions on measuring the electron density of the gas. Many H II regions contain radial variations in density. We suggest that the ISM pressure is a more meaningful quantity to measure in H II regions or galaxies. Specific combinations of line ratios can cover the full range of ISM pressures (4 < log(P/k) < 9). As H II regions become resolved at increasingly high redshift through the next generation telescopes, we anticipate that these diagnostics will be important for understanding the conditions around the young, hot stars from the early universe to the present day.

We examine the diagnostic power of rest-frame ultraviolet (UV) nebular emission lines, and compare them to more commonly used rest-frame optical emission lines, using the test case of a single star-forming knot of the bright lensed galaxy RCSGA 032727-132609 at redshift z~1.7. This galaxy has complete coverage of all the major rest-frame UV and optical emission lines from Magellan/MagE and Keck/NIRSPEC. Using the full suite of diagnostic lines, we infer the physical properties: nebular electron temperature (T_e), electron density (n_e), oxygen abundance (log(O/H)), ionisation parameter (log(q)) and interstellar medium (ISM) pressure (log(P/k)). We examine the effectiveness of the different UV, optical and joint UV-optical spectra in constraining the physical conditions. Using UV lines alone we can reliably estimate log(q), but the same is difficult for log(O/H). UV lines yield a higher (~1.5 dex) log(P/k) than the optical lines, as the former probes a further inner nebular region than the latter. For this comparison, we extend the existing Bayesian inference code IZI, adding to it the capability to infer ISM pressure simultaneously with metallicity and ionisation parameter. This work anticipates future rest-frame UV spectral datasets from the James Webb Space Telescope (JWST) at high redshift and from the Extremely Large Telescope (ELT) at moderate redshift.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

We stack the rest-frame ultraviolet spectra of N=14 highly magnified gravitationally lensed galaxies at redshifts 1.6<z<3.6. The resulting new composite spans $900< \lambda_{rest} < 3000$ \AA, with a peak signal-to-noise ratio of 103 per spectral resolution element ($\sim$100 km/s). It is the highest signal-to-noise ratio, highest spectral resolution composite spectrum of $z\sim2$--3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of $z\sim0$ COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at $z\sim 0$ and $z\sim2$, suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorption-weighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least $-2000$ km/s . We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at $z \sim 2$. We publicly release our composite spectra.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We report the detection of extended Lyman-$\alpha$ emission from the host galaxy of SDSS~J2222+2745, a strongly lensed quasar at $z = 2.8$. Spectroscopic follow-up clearly reveals extended Lyman-$\alpha$ in emission between two images of the central active galactic nucleus (AGN). We reconstruct the lensed quasar host galaxy in the source plane by applying a strong lens model to HST imaging, and resolve spatial scales as small as $\sim$200 parsecs. In the source plane we recover the host galaxy morphology to within a few hundred parsecs of the central AGN, and map the extended Lyman-$\alpha$ emission to its physical origin on one side of the host galaxy at radii $\sim$0.5-2 kpc from the central AGN. There are clear morphological differences between the Lyman-$\alpha$ and rest-frame ultraviolet stellar continuum emission from the quasar host galaxy. Furthermore, the relative velocity profiles of quasar Lyman-$\alpha$, host galaxy Lyman-$\alpha$, and metal lines in outflowing gas reveal differences in the absorbing material affecting the AGN and host galaxy. These data indicate the presence of patchy local intervening gas in front of the central quasar and its host galaxy. This interpretation is consistent with the central luminous quasar being obscured across a substantial fraction of its surrounding solid angle, resulting in strong anisotropy in the exposure of the host galaxy to ionizing radiation from the AGN. This work demonstrates the power of strong lensing-assisted studies to probe spatial scales that are currently inaccessible by other means.

The Epoch of Reionization (EoR) 21-cm signal is expected to become increasingly non-Gaussian as reionization proceeds. We have used semi-numerical simulations to study how this affects the error predictions for the EoR 21-cm power spectrum. We expect $SNR=\sqrt{N_k}$ for a Gaussian random field where $N_k$ is the number of Fourier modes in each $k$ bin. We find that non-Gaussianity is important at high $SNR$ where it imposes an upper limit $[SNR]_l$. For a fixed volume $V$, it is not possible to achieve $SNR > [SNR]_l$ even if $N_k$ is increased. The value of $[SNR]_l$ falls as reionization proceeds, dropping from $\sim 500$ at $\bar{x}_{HI} = 0.8-0.9$ to $\sim 10$ at $\bar{x}_{HI} = 0.15 $ for a $[150.08\, {\rm Mpc}]^3$ simulation. We show that it is possible to interpret $[SNR]_l$ in terms of the trispectrum, and we expect $[SNR]_l \propto \sqrt{V}$ if the volume is increased. For $SNR \ll [SNR]_l$ we find $SNR = \sqrt{N_k}/A $ with $A \sim 0.95 - 1.75$, roughly consistent with the Gaussian prediction. We present a fitting formula for the $SNR$ as a function of $N_k$, with two parameters $A$ and $[SNR]_l$ that have to be determined using simulations. Our results are relevant for predicting the sensitivity of different instruments to measure the EoR 21-cm power spectrum, which till date have been largely based on the Gaussian assumption.

The magnetospheric locations of pulsar radio emission region are not well known. The actual form of the so--called radius--to--frequency mapping should be reflected in the aberration--retardation (A/R) effects that shift and/or delay the photons depending on the emission height in the magnetosphere. Recent studies suggest that in a handful of pulsars the A/R effect can be discerned w.r.t the peak of the central core emission region. To verify these effects in an ensemble of pulsars we launched a project analysing multi--frequency total intensity pulsar profiles obtained from the new observations from the Giant Meterwave Radio Telescope (GMRT), Arecibo Observatory (AO) and archival European Pulsar Network (EPN) data. For all these profiles we measure the shift of the outer cone components with respect to the core component which is necessary for establishing the A/R effect. Within our sample of 23 pulsars 7 show the A/R effects, 12 of them (doubtful cases) show a tendency towards this effect, while the remaining 4 are obvious counter examples. The counter--examples and doubtful cases may arise from uncertainties in determination of the location of the meridional plane and/or the core emission component. It hence appears that the A/R effects are likely to operate in most pulsars from our sample. We conclude that in cases where those effects are present the core emission has to originate below the conal emission region.

We report the discovery of a young pulsar associated with the supernova remnant G21.5-0.9, using the Giant Metrewave Radio Telescope (GMRT) located near Pune, India. Discovered at a frequency of 610 MHz, J1833-1034 has a period of 61.86 ms and a period derivative of $2.0 \times 10^{-13}$, making it similar to other known young pulsars. The characteristic age of the pulsar is $\approx 4900$ yr, somewhat higher than estimates for the age of the remnant, but not incompatible with it. The pulsar has a spin-down luminosity of $3.3 \times 10^{37}$ erg s$^{-1}$, which is the second highest amongst all the known Galactic pulsars.