We present a detailed reconstruction of the star-formation history of the Constellation III region in the Large Magellanic Cloud, to constrain the formation mechanism of this enigmatic feature. Star formation in Constellation III seems to have taken place during two distinct epochs: there is the 8–15 Myr epoch that had previously been recognized, but we also see strong evidence for a separate ‘burst’ of star formation 25–30 Myr ago. The ‘super-supernova' or GRB blast wave model for the formation of Constellation III is difficult to reconcile with such an extended, two-epoch star formation history, because the shock wave should have induced star formation throughout the structure simultaneously, and any unconsumed gas would quickly be dissipated, leaving nothing from which to form a subsequent burst of activity. We propose a ‘truly stochastic’ self-propagating star formation model, distinct from the canonical model in which star formation proceeds in a radially directed wave from the center of Constellation III to its perimeter. As others have noted, and we now confirm, the bulk age gradients demanded by such a model are simply not present in Constellation III. In our scenario, the prestellar gas is somehow pushed into these large-scale arc structures, without simultaneously triggering immediate and violent star formation throughout the structure. Rather, star formation proceeds in the arc according to the local physical conditions of the gas. Self-propagating star formation is certainly possible, but in a truly stochastic manner, without a directed, large scale pattern.

The spin-down power of a pulsar is thought to be carried away in an MHD wind in which, at least close to the star, the energy transport is dominated by Poynting flux. The pulsar drives a low frequency wave in this wind, consisting of stripes of toroidal magnetic field of alternating polarity, propagating in a region around the equatorial plane. The current implied by this configuration falls off more slowly with radius than the number of charged particles available to carry it, so that the MHD picture must, at some point, fail. Recently, magnetic reconnection in such a structure has been shown to accelerate the wind significantly. This reduces the magnetic field in the comoving frame and, consequently, the required current, enabling the solution to extend to much larger radius. This scenario is discussed and, for the Crab Nebula, the range of validity of the MHD solution is compared with the radius at which the flow appears to terminate. For sufficiently high particle densities, it is shown that a low frequency entropy wave can propagate out to the termination point. In this case, the ‘termination shock’ itself must be responsible for dissipating the wave.

This paper is dedicated to Don Melrose on his 60th birthday.

We present an overview of the SkyMapper optical follow-up programme for gravitational-wave event triggers from the LIGO/Virgo observatories, which aims at identifying early GW170817-like kilonovae out to $\sim200\,\mathrm{Mpc}$ distance. We describe our robotic facility for rapid transient follow-up, which can target most of the sky at $\delta<+10\deg $ to a depth of $i_\mathrm{AB}\approx 20\,\mathrm{mag}$ . We have implemented a new software pipeline to receive LIGO/Virgo alerts, schedule observations and examine the incoming real-time data stream for transient candidates. We adopt a real-bogus classifier using ensemble-based machine learning techniques, attaining high completeness ( $\sim98\%$ ) and purity ( $\sim91\%$ ) over our whole magnitude range. Applying further filtering to remove common image artefacts and known sources of transients, such as asteroids and variable stars, reduces the number of candidates by a factor of more than 10. We demonstrate the system performance with data obtained for GW190425, a binary neutron star merger detected during the LIGO/Virgo O3 observing campaign. In time for the LIGO/Virgo O4 run, we will have deeper reference images allowing transient detection to $i_\mathrm{AB}\approx 21\,\mathrm{mag}$ .

Four massive, early-type stars, three of which are confirmed binaries, have been observed with the Australia Telescope Compact Array at 1.4, 2.4, 4.8, and 8.6 GHz. The earliest star cataloged so far, HD 93129A, was also observed at 17.8 and 24.5 GHz. Here we present an analysis of the spectra as well as the structure of the stellar systems. All four spectra show clear evidence of non-thermal emission, indicative of a binary system with a colliding wind region. We discuss the magnetic field of the emitting region of HD 93129A and make predictions on the radiation at high energies. Archive X-ray observations towards the target sources are also investigated and interpreted in the light of the non-thermal emission detected.

The low-redshift Lyα forest of absorption lines provides a probe of large-scale baryonic structures in the intergalactic medium, some of which may be remnants of physical conditions set up during the epoch of galaxy formation. We discuss our recent Hubble Space Telescope (HST) observations and interpretation of low-z Lyα clouds toward nearby Seyferts and QSOs, including their frequency, space density, estimated mass, association with galaxies, and contribution to Ωb. Our HST/GHRS detections of ∼ 70 Lyα absorbers with Nhi ≥ 1012·6 cm−2 along 11 sightlines covering pathlength Δ(cz) = 114,000 km s−1 show f (>Nhi) α Nhi−0·63±0·04 and a line frequency dN/dz = 200 ± 40 for Nhi > 1012·6 cm−2 (one every 1500 km s−1 of redshift). A group of strong absorbers toward PKS 2155–304 may be associated with gas (400–800) kpc from four large galaxies, with low metallicity (≤0·003 solar) and D/H ≤ 2 × 10−4. At low-z, we derive a metagalactic ionising radiation field from AGN of J0 = × 10−23 erg cm−2 s−1 Hz−1 sr−1 and a Lyα-forest baryon density Ωb =(0·008 ± 0·004)[J−23N14b100]½ for clouds of characteristic size b = (100 kpc)b100.

We present a new procedure for photometric parallax estimation. The data for 1236 stars provide calibrations between the absolute magnitude offset from the Hyades main-sequence and the ultraviolet-excess for eight different (B–V)0 colour-index intervals, (0.3 0.4), (0.4 0.5), (0.5 0.6), (0.6 0.7), (0.7 0.8), (0.8 0.9), (0.9 1.0) and (1.0 1.1). The mean difference between the original and estimated absolute magnitudes and the corresponding standard deviation are rather small, +0.0002 and ±0.0613 mag. The procedure has been adapted to the Sloan photometry by means of colour equations and applied to a set of artificial stars with different metallicities. The comparison of the absolute magnitudes estimated by the new procedure and the canonical one indicates that a single colour–magnitude diagram does not supply reliable absolute magnitudes for stars with large range of metallicity.

The Australia Telescope was used in March–April 2005 to observe the 1.384 and 2.368-GHz emissions from the RS CVn binary HR 1099 in two sessions, each of 9-h duration and 11 days apart. Two intervals of highly polarised emission, each lasting 2–3 h, were recorded. During this coherent emission we employed a recently installed facility to sample the data at 78-ms intervals to measure the fine temporal structure and, in addition, all the data were used to search for fine spectral structure. We present the following observational results: (1) ∼100% left-hand circularly polarised emission was seen at both 1.384 and 2.368 GHz during separate epochs; (2) the intervals of highly polarised emission lasted for 2–3 h on each occasion; (3) three 22-min integrations made at 78-ms time resolution showed that the modulation index of the Stokes V parameter increased monotonically as the integration time was decreased and was still increasing at our resolution limit; (4) the extremely fine temporal structure strongly indicates that the highly polarised emission is due to an electron-cyclotron maser operating in the corona of one of the binary components; (5) the first episode of what we believe is ECME (electron-cyclotron maser emission) at 1.384 GHz contained a regular frequency structure of bursts with FWHM ∼48 MHz, which drifted across the spectrum at ∼0.7 MHz min−1. Our second episode of ECME at 2.368 GHz contained wider-band frequency structure, which did not permit us to estimate an accurate bandwidth or direction of drift; (6) the two ECME events reported in this paper agree with six others reported in the literature in occurring in the binary orbital phase range 0.5–0.7; (7) in one event of 8-h duration, two independent maser sources were operating simultaneously at 1.384 and 2.368 GHz.

We discuss two kinds of maser sources that may be responsible for driving the observed events that we believe are powered by ECME. One is based on the widely reported ‘loss-cone anisotropy', the second on an auroral analogue, which is driven by an unstable ‘horseshoe' distribution of fast-electron velocities with respect to the magnetic field direction. Generally, we favour the latter, because of its higher growth rate and the possibility of the escape of radiation which has been emitted at the fundamental electron cyclotron frequency. If the auroral analogue is operating, the magnetic field in the source cavity is ∼500 G at 1.384 GHz and ∼850 G at 2.368 GHz; the source brightness temperatures are of the order TB ∼ 1015 K.

We suggest that the ECME source may be an aurora-like phenomenon due to the transfer of plasma from the K2 subgiant to the G5 dwarf in a strong stellar wind, an idea that is based on VLBA maps showing the establishment of an 8.4 GHz source near the G5 dwarf at times of enhanced radio activity in HR 1099.

It is well known that thermally pulsing Asymptotic Giant Branch stars with low mass play a relevant role in the chemical evolution. They have synthesized about 30% of the galactic carbon and provide an important contribution to the nucleosynthesis of heavy elements (A > 80). The relevant nucleosynthesis site is the He-rich intermediate zone (less than 10−2 M⊙), where α(2α,γ)12C reactions and slow neutron captures on seed nuclei (essentially iron) take place. A key ingredient is the interplay between nuclear processes and convective mixing. It is the partial overlap of internal and external convective zones that allows the dredge-up of the material enriched in C and heavy elements. We review the progresses made in the last 50 years in the comprehension of the s process in AGB stars, with special attention to the identification of the main neutron sources and to the particular physical conditions allowing this important nucleosynthesis.

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5-m optical/infrared telescope to be located at Dome C on the Antarctic plateau. Conditions at Dome C are known to be exceptional for astronomy. The seeing (above ∼30 m height), coherence time, and isoplanatic angle are all twice as good as at typical mid-latitude sites, while the water-vapour column, and the atmosphere and telescope thermal emission are all an order of magnitude better. These conditions enable a unique scientific capability for PILOT, which is addressed in this series of papers. The current paper presents an overview of the optical and instrumentation suite for PILOT and its expected performance, a summary of the key science goals and observational approach for the facility, a discussion of the synergies between the science goals for PILOT and other telescopes, and a discussion of the future of Antarctic astronomy. Paper II and Paper III present details of the science projects divided, respectively, between the distant Universe (i.e. studies of first light, and the assembly and evolution of structure) and the nearby Universe (i.e. studies of Local Group galaxies, the Milky Way, and the Solar System).

The spectral line datacubes obtained from the Square Kilometre Array (SKA) and its precursors, such as the Australian SKA Pathfinder (ASKAP), will be sufficiently large to necessitate automated detection and parametrisation of sources. Matched filtering is widely acknowledged as the best possible method for the automated detection of sources. This paper presents the Characterised Noise Hi (cnhi) source finder, which employs a novel implementation of matched filtering. This implementation is optimised for the 3-D nature of the Hi spectral line observations of the planned Wide-field ASKAP Legacy L-band All-sky Blind surveY (WALLABY). The CNHI source finder also employs a novel sparse representation of 3-D objects, with a high compression rate, to implement the Lutz one-pass algorithm on datacubes that are too large to process in a single pass. WALLABY will use ASKAP's phenomenal 30 square degree field of view to image ∼70% of the sky. It is expected that WALLABY will find 500 000 Hi galaxies out to z ∼0.2.

A radio continuum survey of X-ray-identified weak-lined T Tauri stars (WTTs) in the newly-discovered η Chamaeleontis cluster has been completed using the Australia Telescope Compact Array (ATCA). The 10 known WTTs in the cluster form a unique sample of codistant late-type pre-main-sequence stars with ages of ~8 Myr and masses ranging from 0·1–1·0 M⊙. Our survey detected none of the 10 X-ray-emitting WTTs with 3σ sensitivity limits at 4·8 and 8·6 GHz (6·2 and 3·5 cm) of typically 0·4 mJy, corresponding to a radio luminosity of 4·5 ×1015 erg Hz−1 s−1. Rotation periods for these stars indicate that they are not, as a group, fast-rotating stars. The non-detection in the radio bands supports the findings of other radio surveys of inhomogeneous samples of young stars, where radio emission is fairly common (10–30%) among very young T Tauri stars across all late spectral types, but confined to rapidly-rotating F-G-K stars amongst older zero-age main sequence stars. Rotation, more than youth, appears to be the key to radio emission in young stars.

Barred spirals are well-known to include several ‘distinct components’ or morphological structures in the light distribution (Sandage 1961; de Vaucouleurs 1959; Kormendy 1979). In addition to the bar, spheroid, and disk components, one often observes ring-like enhancements at discrete relative positions within the disk. These inner rings (r), outer rings (R), and nuclear rings (nr) have been the subject of recent, fairly intensive observational (Buta 1984a; Kormendy J979; de Vaucouleurs and Buta 1980a,b; Athanassoula et al. 1982) and theoretical (Schwarz 1979,1981,1984a,b,c; Simkin, Su, and Schwarz 1980) research, and it now seems well-established that most rings are linked to dynamical resonances with the bar.

Extreme scattering events (ESEs) are dramatic variations of the flux density at gigahertz frequencies caused by ray path distortions within an isolated inhomogeneity (‘plasma lens’) in the interstellar medium. These events are characterised by a deep flux density minimum in the light curve with, in some cases, surrounding maxima. The variability time scales range from weeks to months. These phenomena show a strong frequency dependence, in which the variability amplitudes increase with wavelength. During an intraday variability (IDV) monitoring project (March 2000), a feature resembling an ESE-like event appeared in the variable light curve of 0954+658, however with a time scale of less than two days. We will discuss this effect and its implications for a better description of the interstellar medium.

In 1975, Hoessel, Elias, Wade and Huchra commenced a near infrared survey of 80 fields in the northern Milky Way with the Palomar 1.2 m Schmidt telescope, (Hoessel et al. 1979). This has now been issued as an atlas reproduced in the form of photographic paper prints. In 1977, the SRC 1.2 m Schmidt telescope at Siding Spring was authorized to commence a complementary survey of the southern Milky Way, consisting of the 151 ESO/SRC survey fields which have centres within 10° of the galactic plane and negative declinations (see Fig. 1). A further 12 fields have subsequently been added to the survey to permit coverage of the Large and Small Magellanic Clouds.

During the solar maximum of 1989–91 an unprecedented sequence of 13 cosmic ray ground-level enhancements (GLEs) was observed by the world-wide neutron monitor network. Of particular interest were two GLEs observed by the Australian network. The 1989 September 29 event was the largest GLE in the space era while the October 22 GLE included an highly anisotropic precursor peak.

Analysis of both these GLEs, taking into account disturbed geomagnetic conditions, shows that the particle arrivals at the earth were unusual. The September 29 GLE had significant particle propagation in the reverse direction and as the particle flux decreased following the peak the spectrum also softened. In contrast, the 1989 October 22 precursor exhibited extreme anisotropy while the particles involved in the main GLE showed a complex temporal structure possibly indicating multiple particle injection at the solar acceleration region.

A low-cost differential image motion monitor (DIMM), consisting of an 11-inch Celestron, an SBIG ST-4 autoguiding CCD camera and a PC, is described. Two such systems were used during June–July and November–December 1993 to make near-simultaneous seeing measurements at Freeling Heights in the Northern Flinders Ranges and at Siding Spring Observatory. The results of these campaigns show that the seeing-distribution is generally similar at both sites, with the most common seeing value being ~l·2″. Siding Spring does, however, have slightly more bad seeing (>2″) than Freeling Heights. Weather records from Arkaroola Resort (15 km south of Freeling Heights) indicate that there is ~15% less cloud cover at Freeling Heights than at Siding Spring. Episodes of rapid seeing deterioration at Siding Spring in winter are shown to coincide with warm air masses crossing the mountain.

Pulse arrival time measurements allow the determination of accurate pulsar periods, period derivatives and, provided the data span is at least one year, precise pulsar positions. If observations are frequent and reasonably regular, irregularities in the period can also be investigated. To minimize the effect of possible variations in dispersion measure, it is important that these observations be made at a relatively high frequency, preferably above 1 GHz. To eliminate pulse shape variations due to variable ionospheric Faradayrotation, the pulse total intensity or one of the circular polarizations must be recorded.

Most if not all metre-wave Type III bursts appear to be associated with solar active regions (e.g. Duncan 1977). There is a well-known association of Type III bursts and flares within a few minutes of the onset of the brightenings seen in Hα (Wild et al. 1954). Early descriptions of these observations were given by Malville (1961) (or see Kundu (1965) for review). More recently Jackson et al. (1978) have reported the occurrence of a broad maximum of Type III bursts 10 to 5 h prior to mass ejection transients seen by the white light coronagraph on Skylab.

Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Monte Carlo techniques such as Markov Chain Monte Carlo (MCMC) and nested sampling. However, these techniques are time-consuming and their convergence to the posterior could be difficult to determine. In this study, we advocate variational inference as an alternative to solve the above problems, and demonstrate its usefulness for parameter estimation and model selection in astrophysics. Variational inference converts the inference problem into an optimisation problem by approximating the posterior from a known family of distributions and using Kullback–Leibler divergence to characterise the difference. It takes advantage of fast optimisation techniques, which make it ideal to deal with large datasets and makes it trivial to parallelise on a multicore platform. We also derive a new approximate evidence estimation based on variational posterior, and importance sampling technique called posterior-weighted importance sampling for the calculation of evidence, which is useful to perform Bayesian model selection. As a proof of principle, we apply variational inference to five different problems in astrophysics, where Monte Carlo techniques were previously used. These include assessment of significance of annual modulation in the COSINE-100 dark matter experiment, measuring exoplanet orbital parameters from radial velocity data, tests of periodicities in measurements of Newton’s constant G, assessing the significance of a turnover in the spectral lag data of GRB 160625B, and estimating the mass of a galaxy cluster using weak gravitational lensing. We find that variational inference is much faster than MCMC and nested sampling techniques for most of these problems while providing competitive results. All our analysis codes have been made publicly available.

Although the Sloan Digital Sky Survey (SDSS) was primarily envisioned as a tool for understanding the nature of the ‘high redshift’ universe, significant discoveries have already been made at lower redshift, z ∼ 0, through studies of stars in the Milky Way galaxy. We have begun to explore the nature of the Milky Way by detailed investigation of the publicly accessible SDSS archive, using spectroscopically targeted stars of special interest (e.g. field horizontal-branch stars, carbon-enhanced stars, and F- and G-type turnoff stars), as well as the stars originally selected as photometric and reddening standards. The first step is to use the SDSS data (which includes independently calibrated five-band photometry and spectrophotometry of individual stars) to derive reliable estimates of the stellar physical parameters, such as Teff, log g, and [Fe/H], for stars that have been observed to date. Of particular interest, at present, are the stars that are apparently associated with the Monoceros Stream (also known as the SDSS ‘Ring around the Galaxy’), for which we report derived metallicities. The techniques we have developed for derivation of the physical parameters for these stars are presently being applied to other stars in the SDSS database, including the Early Data Release (EDR), as well as the first official public database, DR-1. Here we report on the progress made to date, and comment on what might be explored in the near future from a dedicated extension of the SDSS survey (SEGUE) that specifically targets stars in the Milky Way.