An objective classification of 147 globular clusters (GCs) in the inner region of the giant elliptical galaxy M87 is carried out with the help of two methods of multivariate analysis. First, independent component analysis (ICA) is used to determine a set of independent variables that are linear combinations of various observed parameters (mostly Lick indices) of the GCs. Next, K-means cluster analysis (CA) is applied on the independent components (ICs), to find the optimum number of homogeneous groups having an underlying structure. The properties of the four groups of GCs thus uncovered are used to explain the formation mechanism of the host galaxy. It is suggested that M87 formed in two successive phases. First a monolithic collapse, which gave rise to an inner group of metal-rich clusters with little systematic rotation and an outer group of metal-poor clusters in eccentric orbits. In a second phase, the galaxy accreted low-mass satellites in a dissipationless fashion, from the gas of which the two other groups of GCs formed. Evidence is given for a blue stellar population in the more metal rich clusters, which we interpret by Helium enrichment. Finally, it is found that the clusters of M87 differ in some of their chemical properties (NaD, TiO1, light-element abundances) from GCs in our Galaxy and M31.

The regular or chaotic dynamics of an analytical realistic three dimensional model composed of a spherically symmetric central nucleus, a bar and a flat disk is investigated. For describing the properties of the bar, we introduce a new simple dynamical model and we explore the influence on the character of orbits of all the involved parameters of it, such as the mass and the scale length of the bar, the major semi-axis and the angular velocity of the bar, as well as the energy. Regions of phase space with ordered and chaotic motion are identified in dependence on these parameters and for breaking the rotational symmetry. First, we study in detail the dynamics in the invariant plane z = pz = 0 using the Poincaré map as a basic tool and then study the full three-dimensional case using the Smaller Alignment index method as principal tool for distinguishing between order and chaos. We also present strong evidence obtained through the numerical simulations that our new bar model can realistically describe the formation and the evolution of the observed twin spiral structure in barred galaxies.

We estimated iron and metallicity gradients in the radial and vertical directions with the F and G type dwarfs taken from the Radial Velocity Experiment Data Release 4 database. The sample defined by the constraints Zmax ⩽ 825 pc and ep ⩽ 0.10 consists of stars with metal abundances and space velocity components agreeable with the thin-disc stars. The radial iron and metallicity gradients estimated for the vertical distance intervals 0 < Zmax ⩽ 500 and 500 < Zmax ⩽ 800 pc are $\text{d}[\text{Fe}/\text{H}]/\text{d}R_{\rm m}=-0.083\pm 0.030$ and $\text{d}[\text{Fe}/\text{H}]/\text{d}R_{\rm m}=-0.048\pm 0.037$ dex kpc−1, and $\text{d}[\text{M}/\text{H}]/\text{d}R_{\rm m}=-0.063\pm 0.011$ and $\text{d}[\text{M}/\text{H}]/\text{d}R_{\rm m}=-0.028\pm 0.057$ dex kpc−1, respectively, where Rm is the mean Galactocentric distance. The iron and metallicity gradients for less number of stars at further vertical distances, 800 < Zmax ⩽ 1500 pc, are mostly positive. Compatible iron and metallicity gradients could be estimated with guiding radius (Rg) for the same vertical distance intervals 0 < Zmax ⩽ 500 and 500 < Zmax ⩽ 800 pc, i.e. $\text{d}[\text{Fe}/\text{H}]/\text{d}R_{\rm g}=-0.083\pm 0.030$ and $\text{d}[\text{Fe}/\text{H}]/\text{d}R_{\rm g}=-0.065\pm 0.039$ dex kpc−1; $\text{d}[\text{M}/\text{H}]/\text{d}R_{\rm g}=-0.062\pm 0.018$ and $\text{d}[\text{M}/\text{H}]/\text{d}R_{\rm g}=-0.055\pm 0.045$ dex kpc−1. F and G type dwarfs on elongated orbits show a complicated radial iron and metallicity gradient distribution in different vertical distance intervals. Significant radial iron and metallicity gradients could be derived neither for the sub-sample stars with Rm ⩽ 8 kpc, nor for the ones at larger distances, Rm > 8 kpc. The range of the iron and metallicity abundance for the F and G type dwarfs on elongated orbits, [−0.13, −0.01), is similar to the thin-disc stars, while at least half of their space velocity components agree better with those of the thick-disc stars. The vertical iron gradients estimated for the F and G type dwarfs on circular orbits are $\text{d}[\text{Fe}/\text{H}]/\text{d}Z_{{\rm max}}=-0.176\pm 0.039$ dex kpc−1 and $\text{d}[\text{Fe}/\text{H}]/\text{d}Z_{{\rm max}}=-0.119\pm 0.036$ dex kpc−1 for the intervals Zmax ⩽ 825 and Zmax ⩽ 1500 pc, respectively.

The total magnification due to a point lens has been of particular interest as the theorem that gravitational lensing results in light amplification for all observers appears to contradict the conservation of photon number. This has been discussed several times, and various resolutions have been offered. In this note, we use a kinematic approach to provide a formula for the magnification factor for the primary image accurate to first order and valid for rays leaving the source at any trajectory. We thus determine the magnification over a sphere surrounding the system. A new result found is that while the magnification dips below unity far from the optical axis as noted by others, it returns to unity directly behind the source.

Determining when and how the first galaxies reionised the intergalactic medium promises to shed light on both the nature of the first objects and the cosmic history of baryons. Towards this goal, quasar absorption lines play a unique role by probing the properties of diffuse gas on galactic and intergalactic scales. In this review, we examine the multiple ways in which absorption lines trace the connection between galaxies and the intergalactic medium near the reionisation epoch. We first describe how the Ly α forest is used to determine the intensity of the ionising ultraviolet background and the global ionising emissivity budget. Critically, these measurements reflect the escaping ionising radiation from all galaxies, including those too faint to detect directly. We then discuss insights from metal absorption lines into reionisation-era galaxies and their surroundings. Current observations suggest a buildup of metals in the circumgalactic environments of galaxies over z ~ 6 to 5, although changes in ionisation will also affect the evolution of metal line properties. A substantial fraction of metal absorbers at these redshifts may trace relatively low-mass galaxies. Finally, we review constraints from the Ly α forest and quasar near zones on the timing of reionisation. Along with other probes of the high-redshift Universe, absorption line data are consistent with a relatively late end to reionisation (5.5 ≲ z ≲ 7); however, the constraints are still fairly week. Significant progress is expected to come through improved analysis techniques, increases in the number of known high-redshift quasars from optical and infrared sky surveys, large gains in sensitivity from next-generation observing facilities, and synergies with other probes of the reionisation era.

The first observations by a worldwide network of advanced interferometric gravitational wave detectors offer a unique opportunity for the astronomical community. At design sensitivity, these facilities will be able to detect coalescing binary neutron stars to distances approaching 400 Mpc, and neutron star–black hole systems to 1 Gpc. Both of these sources are associated with gamma-ray bursts which are known to emit across the entire electromagnetic spectrum. Gravitational wave detections provide the opportunity for ‘multi-messenger’ observations, combining gravitational wave with electromagnetic, cosmic ray, or neutrino observations. This review provides an overview of how Australian astronomical facilities and collaborations with the gravitational wave community can contribute to this new era of discovery, via contemporaneous follow-up observations from the radio to the optical and high energy. We discuss some of the frontier discoveries that will be made possible when this new window to the Universe is opened.

Theoretical models of high-mass star formation lie between two extreme scenarios. At one extreme, all the mass comes from an initially gravitationally bound core. At the other extreme, the majority of the mass comes from cluster scale gas, which lies far outside the initial core boundary. One way to unambiguously show high-mass stars can assemble their gas through the former route would be to find a high-mass star forming in isolation. Making use of recently available CORNISH and ATLASGAL Galactic plane survey data, we develop sample selection criteria to try and find such an object. From an initial list of approximately 200 sources, we identify the high-mass star-forming region G13.384 + 0.064 as the most promising candidate. The region contains a strong radio continuum source, that is powered by an early B-type star. The bolometric luminosity, derived from infrared measurements, is consistent with this. However, sub-millimetre continuum emission, measured in ATLASGAL, as well as dense gas tracers, such as HCO+(3–2) and N2H+(3–2) indicate that there is less than ~ 100 M⊙ of material surrounding this star. We conclude that this region is indeed a promising candidate for a high-mass star forming in isolation.

We review existing smoothed particle hydrodynamics setup methods and outline their advantages, limitations, and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tessellations. Particles are placed and iteratively moved based on their proximity to neighbouring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements.