In this paper, we introduce Nicil: Non-Ideal magnetohydrodynamics Coefficients and Ionisation Library. Nicil is a stand-alone Fortran90 module that calculates the ionisation values and the coefficients of the non-ideal magnetohydrodynamics terms of Ohmic resistivity, the Hall effect, and ambipolar diffusion. The module is fully parameterised such that the user can decide which processes to include and decide upon the values of the free parameters, making this a versatile and customisable code. The module includes both cosmic ray and thermal ionisation; the former includes two ion species and three species of dust grains (positively charged, negatively charged, and neutral), and the latter includes five elements which can be doubly ionised. We demonstrate tests of the module, and then describe how to implement it into an existing numerical code.

We describe the performance of the Boolardy Engineering Test Array, the prototype for the Australian Square Kilometre Array Pathfinder telescope. Boolardy Engineering Test Array is the first aperture synthesis radio telescope to use phased array feed technology, giving it the ability to electronically form up to nine dual-polarisation beams. We report the methods developed for forming and measuring the beams, and the adaptations that have been made to the traditional calibration and imaging procedures in order to allow BETA to function as a multi-beam aperture synthesis telescope. We describe the commissioning of the instrument and present details of Boolardy Engineering Test Array’s performance: sensitivity, beam characteristics, polarimetric properties, and image quality. We summarise the astronomical science that it has produced and draw lessons from operating Boolardy Engineering Test Array that will be relevant to the commissioning and operation of the final Australian Square Kilometre Array Path telescope.

We aim to discover the accuracy of photometric mass ratios (qph) determined for eclipsing binary stars, in the case of the system having at least one ‘flat bottom’ as a minimum profile, as well as the accuracy of data used in that sense. Within this context, we present the results of two-dimensional grid search (q – i) for some W UMa-type eclipsing binaries showing total eclipses, based on the high precision photometric data provided by the KEPLER Mission. The radial velocity data obtained for KIC10618253 in this study, enables us to compare both qph and the corresponding spectroscopic mass ratio (qsp) values. The results indicate that the high precision photometric data for overcontact eclipsing binaries showing total eclipses allow us to obtain the photometric mass ratios as accurate as the spectroscopic values.

We investigate the Milky Way Galaxy’s radial and vertical metallicity gradients using a sample of 47 406 red clump stars from the RAdial Velocity Experiment Data Release 4. Distances are calculated by adopting Ks-band absolute magnitude as −1.54±0.04 mag for the sample. The metallicity gradients are calculated with their current orbital positions (Rgc and Z) and with their orbital properties (Rm and zmax): d[Fe/H]/dRgc = −0.047±0.003 dex kpc−1 for |Z| ≤ 0.5 kpc and d[Fe/H]/dRm = −0.025±0.002 dex kpc−1 for zmax ≤ 0.5 kpc. This reaffirms the radial metallicity gradient in the thin disc but highlights that gradients are sensitive to the selection effects caused by the difference between Rgc and Rm. The radial gradient is flat in the distance interval 0.5-1 kpc from the plane and then becomes positive greater than 1 kpc from the plane. The radial metallicity gradients are also eccentricity dependent. We showed that d[Fe/H]/dRm = −0.089±0.010, −0.073±0.007, −0.053±0.004 and −0.044±0.002 dex kpc−1 for ep ≤ 0.05, ep ≤ 0.07, ep ≤ 0.10 and ep ≤ 0.20 sub-samples, respectively, in the distance interval zmax ≤ 0.5 kpc. Similar trend is found for vertical metallicity gradients. Both the radial and vertical metallicity gradients are found to become shallower as the eccentricity of the sample increases. These findings can be used to constrain different formation scenarios of the thick and thin discs.

Here, we present a catalogue of known Fast Radio Burst sources in the form of an online catalogue, FRBCAT. The catalogue includes information about the instrumentation used for the observations for each detected burst, the measured quantities from each observation, and model-dependent quantities derived from observed quantities. To aid in consistent comparisons of burst properties such as width and signal-to-noise ratios, we have re-processed all the bursts for which we have access to the raw data, with software which we make available. The originally derived properties are also listed for comparison. The catalogue is hosted online as a Mysql database which can also be downloaded in tabular or plain text format for off-line use. This database will be maintained for use by the community for studies of the Fast Radio Burst population as it grows.

An early-type, massive, short-period ($P_{\text{orb}}=2^d.310951$) eclipsing spectroscopic binary DN Cas has been re-visited with new spectral and photometric data. The masses and radii of the components have been obtained as $M_1=19.04\pm 0.07\,\text{M}_\odot$, $M_2=13.73\pm 0.05\,\text{M}_\odot$ and $R_1=7.22\pm 0.06\,\text{R}_\odot$, $R_2=5.79\pm 0.06\,\text{R}_\odot$, respectively. Both components present synchronous rotation ($V_{\text{rot}1}=160\,\text{km } \text{s}^{-1}$, $V_{\text{rot}2}=130\ \text{km} \,\text{s}^{-1}$) with their orbit. Orbital period analysis yielded a physically bound additional component in the system with a minimum mass of $M_3=0.88\,\text{M}_\odot$ orbiting in an eccentric orbit (e = 0.37 ± 0.2) with an orbital period of P12 = 42 ± 9 yr. High precision absolute parameters of the system allowed us to derive a distance to DN Cas as 1.7 ± 0.2 kpc which locates the system within the borders of the Cas OB6 association (d = 1.8 kpc). The space velocities and the age of DN Cas are in agreement with those of Cas OB6. The age of DN Cas (τ = 3–5 Myr) is found to be 1–2 Myr older than the embedded clusters (IC 1795, IC 1805, and IC 1848) in the Cas OB6 association, which implies a sequential star formation in the association.

The high Antarctic plateau provides exceptional conditions for infrared observations on account of the cold, dry and stable atmosphere above the ice surface. This paper describes the scientific goals behind the first program to examine the time-varying universe in the infrared from Antarctica — the Kunlun Infrared Sky Survey (KISS). This will employ a 50cm telescope to monitor the southern skies in the 2.4μmKdark window from China's Kunlun station at Dome A, on the summit of the Antarctic plateau, through the uninterrupted 4-month period of winter darkness. An earlier paper discussed optimisation of the Kdark filter for sensitivity (Li et al. 2016). This paper examines the scientific program for KISS. We calculate the sensitivity of the camera for the extrema of observing conditions that will be encountered. We present the parameters for sample surveys that could then be carried out for a range of cadences and sensitivities. We then discuss several science programs that could be conducted with these capabilities, involving star formation, brown dwarfs and hot Jupiters, exoplanets around M dwarfs, the terminal phases of stellar evolution, fast transients, embedded supernova searches, reverberation mapping of AGN, gamma ray bursts and the detection of the cosmic infrared background.

Models of neutrino-driven core-collapse supernova explosions have matured considerably in recent years. Explosions of low-mass progenitors can routinely be simulated in 1D, 2D, and 3D. Nucleosynthesis calculations indicate that these supernovae could be contributors of some lighter neutron-rich elements beyond iron. The explosion mechanism of more massive stars remains under investigation, although first 3D models of neutrino-driven explosions employing multi-group neutrino transport have become available. Together with earlier 2D models and more simplified 3D simulations, these have elucidated the interplay between neutrino heating and hydrodynamic instabilities in the post-shock region that is essential for shock revival. However, some physical ingredients may still need to be added/improved before simulations can robustly explain supernova explosions over a wide range of progenitors. Solutions recently suggested in the literature include uncertainties in the neutrino rates, rotation, and seed perturbations from convective shell burning. We review the implications of 3D simulations of shell burning in supernova progenitors for the ‘perturbations-aided neutrino-driven mechanism,’ whose efficacy is illustrated by the first successful multi-group neutrino hydrodynamics simulation of an 18 solar mass progenitor with 3D initial conditions. We conclude with speculations about the impact of 3D effects on the structure of massive stars through convective boundary mixing.

The Dark Energy Survey is undertaking an observational programme imaging 1/4 of the southern hemisphere sky with unprecedented photometric accuracy. In the process of observing millions of faint stars and galaxies to constrain the parameters of the dark energy equation of state, the Dark Energy Survey will obtain pre-discovery images of the regions surrounding an estimated 100 gamma-ray bursts over 5 yr. Once gamma-ray bursts are detected by, e.g., the Swift satellite, the DES data will be extremely useful for follow-up observations by the transient astronomy community. We describe a recently-commissioned suite of software that listens continuously for automated notices of gamma-ray burst activity, collates information from archival DES data, and disseminates relevant data products back to the community in near-real-time. Of particular importance are the opportunities that non-public DES data provide for relative photometry of the optical counterparts of gamma-ray bursts, as well as for identifying key characteristics (e.g., photometric redshifts) of potential gamma-ray burst host galaxies. We provide the functional details of the DESAlert software, and its data products, and we show sample results from the application of DESAlert to numerous previously detected gamma-ray bursts, including the possible identification of several heretofore unknown gamma-ray burst hosts.

Wepresent and evaluate several strategies to search for prompt, low-frequency radio emission associated with gravitational wave transients using the Murchison Widefield Array. As we are able to repoint the Murchison Widefield Array on timescales of tens of seconds, we can search for the dispersed radio signal that has been predicted to originate along with or shortly after a neutron star-neutron star merger. We find that given the large, 600 deg2 instantaneous field of view of the Murchison Widefield Array, we can cover a significant fraction of the predicted gravitational wave error region, although due to the complicated geometry of the latter, we only cover > 50% of the error region for approximately 5% of events, and roughly 15% of events will be located < 10° from the Murchison Widefield Array pointing centre such that they will be covered in the radio images. For optimal conditions, our limiting flux density for a 10-s long transient would be 0.1 Jy, increasing to about 1 Jy for a wider range of events. This corresponds to luminosity limits of 1038−39 erg s−1 based on expectations for the distances of the gravitational wave transients, which should be sufficient to detect or significantly constrain a range of models for prompt emission.

The detection of quasars at z > 6 unveils the presence of supermassive black holes of a few billion solar masses. The rapid formation process of these extreme objects remains a fascinating and open issue. Such discovery implies that seed black holes must have formed early on, and grown via either rapid accretion or BH/galaxy mergers. In this theoretical review, we discuss in detail various BH seed formation mechanisms and the physical processes at play during their assembly. We discuss the three most popular BH formation scenarios, involving the (i) core-collapse of massive stars, (ii) dynamical evolution of dense nuclear star clusters, (iii) collapse of a protogalactic metal free gas cloud. This article aims at giving a broad introduction and an overview of the most advanced research in the field.

This automated catalogue combines all the largest published optical, radio, and X-ray sky catalogues to find probable radio/X-ray associations to optical objects, plus double radio lobes, using uniform processing against all input data. The total count is 1 002 855 optical objects so presented. Each object is displayed with J2000 astrometry, optical and radio/X-ray identifiers, red and blue photometry, and calculated probabilities and optical field solutions of the associations. This is the third and final edition of this method.

The Protoplanetary Discussions conference—held in Edinburgh, UK, from 2016 March 7th–11th—included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multi-physics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer, and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical, and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on community-driven input, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

Observational constraints on the birth and early evolution of massive black holes come from two extreme regimes. At high redshift, quasars signal the rapid growth of billion-solar-mass black holes and indicate that these objects began remarkably heavy and/or accreted mass at rates above the Eddington limit. At low redshift, the smallest nuclear black holes known are found in dwarf galaxies and provide the most concrete limits on the mass of black hole seeds. Here, we review current observational work in these fields that together are critical for our understanding of the origin of massive black holes in the Universe.

This paper presents the first major data release and survey description for the ANU WiFeS SuperNovA Programme. ANU WiFeS SuperNovA Programme is an ongoing supernova spectroscopy campaign utilising the Wide Field Spectrograph on the Australian National University 2.3-m telescope. The first and primary data release of this programme (AWSNAP-DR1) releases 357 spectra of 175 unique objects collected over 82 equivalent full nights of observing from 2012 July to 2015 August. These spectra have been made publicly available via the WISEREP supernova spectroscopy repository.

We analyse the ANU WiFeS SuperNovA Programme sample of Type Ia supernova spectra, including measurements of narrow sodium absorption features afforded by the high spectral resolution of the Wide Field Spectrograph instrument. In some cases, we were able to use the integral-field nature of the Wide Field Spectrograph instrument to measure the rotation velocity of the SN host galaxy near the SN location in order to obtain precision sodium absorption velocities. We also present an extensive time series of SN 2012dn, including a near-nebular spectrum which both confirms its ‘super-Chandrasekhar’ status and enables measurement of the sub-solar host metallicity at the SN site.

We have used the Australia Telescope Compact Array to search for a number of centimetre wavelength methanol transitions which are predicted to show weak maser emission towards star formation regions. Sensitive, high spatial, and spectral resolution observations towards four high-mass star formation regions which show emission in a large number of class II methanol maser transitions did not result in any detections. From these observations, we are able to place an upper limit of ≲ 1300 K on the brightness temperature of any emission from the 31A+–31A−, 17−2–18−3 E (vt = 1), 124–133 A−, 124–133 A+, and 41A+–41A− transitions of methanol in these sources on angular scales of 2 arcsec. This upper limit is consistent with current models for class II methanol masers in high-mass star formation regions and better constraints than those provided here will likely require observations with next-generation radio telescopes.

Correlations are identified between the strength of the λ10830 He I triplet line and the following tracers of stellar activity amongst FGK dwarfs with colours of (B − V) > 0.47: coronal soft X-ray emission, emission in the λ1549 C IV and λ1335 C II lines originating from the transition region, and Ca II H and K emission from the chromosphere. No such correlations are present amongst dwarfs with spectral type earlier than F6. In addition, G and K dwarfs with strong triplet lines show evidence of excess flux in the GALEX FUV band compared to weak-triplet-line dwarfs. The X-ray spectra of late-F, G, and K dwarfs with He I triplets stronger than 160 mÅ have greater values of the ROSAT hardness ratio HR1 than are typical of weak-triplet dwarfs in the same range of spectral type. In other words, dwarfs later than F7V with strong He I triplet lines tend towards harder 0.1–2.0 keV X-ray spectra than weak-triplet dwarfs, although values of HR1 ~ −0.2 to +0.1 can still be encountered amongst a minority of weak-He-triplet stars. As regards, FGK main sequence stars the observational data on the λ10830 triplet line remains sparse. Progress could be made through spectroscopy of high resolution for samples of hundreds of stars, selected on the basis of having other measures of chromospheric and coronal activity available.

A novel sky-subtraction method based on non-negative matrix factorisation with sparsity is proposed in this paper. The proposed non-negative matrix factorisation with sparsity method is redesigned for sky-subtraction considering the characteristics of the skylights. It has two constraint terms, one for sparsity and the other for homogeneity. Different from the standard sky-subtraction techniques, such as the B-spline curve fitting methods and the Principal Components Analysis approaches, sky-subtraction based on non-negative matrix factorisation with sparsity method has higher accuracy and flexibility. The non-negative matrix factorisation with sparsity method has research value for the sky-subtraction on multi-object fibre spectroscopic telescope surveys. To demonstrate the effectiveness and superiority of the proposed algorithm, experiments are performed on Large Sky Area Multi-Object Fiber Spectroscopic Telescope data, as the mechanisms of the multi-object fibre spectroscopic telescopes are similar.

Recent progress in telescope development has brought us different ways to observe protoplanetary disks: interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together to understand disk dispersal and planet formation.

This paper was born at the ‘Protoplanetary Discussions’ meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations converge in unveiling disk structure and evolution, and where they challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions from limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized, specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how their connections and apparent contradictions can help us to disentangle the disk physics and to learn about disk evolution.

Multi-frequency, multi-epoch ATCA observations of a sample of AGN resulted in the identification of nine new candidate Giga-hertz Peaked Spectrum sources. Here, we present Long Baseline Array observations at 4.8 GHz of the four candidates with no previously published VLBI image, and consider these together with previously published VLBI images of the other five sources. We find core-jet or compact double morphologies dominate, with further observations required to distinguish between these two possibilities for some sources. One of the nine candidates, PKS 1831–711, displays appreciable variability, suggesting its GPS spectrum is more ephemeral in nature. We focus in particular on the apparent relationship between a narrow spectral width and ‘compact double’ parsec-scale morphology, finding further examples, but also exceptions to this trend. An examination of the VLBI morphologies high-redshift (z > 3) sub-class of GPS sources suggests that core-jet morphologies predominate in this class.