We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae.

The light curves of Type Ia supernovae are routinely used to constrain cosmology models. Driven by radioactive decay of 56Ni, the light curves steadily decline over time, but >150 days past explosion, the near-infrared portion is poorly characterized. We report a year-long plateau in the near-infrared light curve at 150-500 days, followed by a second decline phase accompanied by a possible appearance of [Fe I] emission lines. This near-infrared plateau contrasts sharply with Type IIP plateaus and requires a new physical mechanism. We suggest a such as masking of the "near-infrared catastrophe," a predicted yet unobserved sharp light-curve decline, by scattering of ultraviolet photons to longer wavelengths. The transition off the plateau could be due to a change in the dominant ionization state of the supernova ejecta. Our results shed new light on the complex radiative transfer processes that take place in Type Ia supernovae and enhance their use as "standard candles."

The late-time spectra of Type Ia supernovae (SNe Ia) are powerful probes of the underlying physics of their explosions. We investigate the late-time optical and near-infrared spectra of seven SNe Ia obtained at the VLT with XShooter at $>$200 d after explosion. At these epochs, the inner Fe-rich ejecta can be studied. We use a line-fitting analysis to determine the relative line fluxes, velocity shifts, and line widths of prominent features contributing to the spectra ([Fe II], [Ni II], and [Co III]). By focussing on [Fe II] and [Ni II] emission lines in the ~7000-7500 \AA\ region of the spectrum, we find that the ratio of stable [Ni II] to mainly radioactively-produced [Fe II] for most SNe Ia in the sample is consistent with Chandrasekhar-mass delayed-detonation explosion models, as well as sub-Chandrasekhar mass explosions that have metallicity values above solar. The mean measured Ni/Fe abundance of our sample is consistent with the solar value. The more highly ionised [Co III] emission lines are found to be more centrally located in the ejecta and have broader lines than the [Fe II] and [Ni II] features. Our analysis also strengthens previous results that SNe Ia with higher Si II velocities at maximum light preferentially display blueshifted [Fe II] 7155 \AA\ lines at late times. Our combined results lead us to speculate that the majority of normal SN Ia explosions produce ejecta distributions that deviate significantly from spherical symmetry.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present a chemical abundance analysis of a metal-poor star, ROA 276, in the stellar system omega Centauri. We confirm that this star has an unusually high [Sr/Ba] abundance ratio. Additionally, ROA 276 exhibits remarkably high abundance ratios, [X/Fe], for all elements from Cu to Mo along with normal abundance ratios for the elements from Ba to Pb. The chemical abundance pattern of ROA 276, relative to a primordial omega Cen star ROA 46, is best fit by a fast-rotating low-metallicity massive stellar model of 20 Msun, [Fe/H] = -1.8, and an initial rotation 0.4 times the critical value; no other nucleosynthetic source can match the neutron-capture element distribution. ROA 276 arguably offers the most definitive proof to date that fast-rotating massive stars contributed to the production of heavy elements in the early Universe.

There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $\omega$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $\omega$ Centauri.

We present a high resolution spectroscopic analysis of 62 red giants in the Milky Way globular cluster NGC5286. We have determined abundances of representative light proton-capture, alpha, Fe-peak and neutron-capture element groups, and combined them with photometry of multiple sequences observed along the colour-magnitude diagram. Our principal results are: (i) a broad, bimodal distribution in s-process element abundance ratios, with two main groups, the s-poor and s-rich groups; (ii) substantial star-to-star Fe variations, with the s-rich stars having higher Fe, e.g. <[Fe/H]>_s-rich - <[Fe/H]>_s-poor ~ 0.2~dex; and (iii) the presence of O-Na-Al (anti-)correlations in both stellar groups. We have defined a new photometric index, c_{BVI}=(B-V)-(V-I), to maximise the separation in the colour-magnitude diagram between the two stellar groups with different Fe and s-element content, and this index is not significantly affected by variations in light elements (such as the O-Na anticorrelation). The variations in the overall metallicity present in NGC5286 add this object to the class of "anomalous" GCs. Furthermore, the chemical abundance pattern of NGC5286 resembles that observed in some of the anomalous GCs, e.g. M22, NGC1851, M2, and the more extreme Omega Centauri, that also show internal variations in s-elements, and in light elements within stars with different Fe and s-elements content. In view of the common variations in s-elements, we propose the term s-Fe-anomalous GCs to describe this sub-class of objects. The similarities in chemical abundance ratios between these objects strongly suggest similar formation and evolution histories, possibly associated with an origin in tidally disrupted dwarf satellites.

We investigate the enrichment in elements produced by the slow neutron-capture process ($s$-process) in the globular clusters M4 (NGC 6121) and M22 (NGC 6656). Stars in M4 have homogeneous abundances of Fe and neutron-capture elements, but the entire cluster is enhanced in $s$-process elements (Sr, Y, Ba, Pb) relative to other clusters with a similar metallicity. In M22, two stellar groups exhibit different abundances of Fe and $s$-process elements. By subtracting the mean abundances of $s$-poor from $s$-rich stars, we derive $s$-process residuals or empirical $s$-process distributions for M4 and M22. We find that the $s$-process distribution in M22 is more weighted toward the heavy $s$-peak (Ba, La, Ce) and Pb than M4, which has been enriched mostly with light $s$-peak elements (Sr, Y, Zr). We construct simple chemical evolution models using yields from massive star models that include rotation, which dramatically increases $s$-process production at low metallicity. We show that our massive star models with rotation rates of up to 50\% of the critical (break-up) velocity and changes to the preferred $^{17}$O($\alpha$,$\gamma$)$^{21}$Ne rate produce insufficient heavy $s$-elements and Pb to match the empirical distributions. For models that incorporate AGB yields, we find that intermediate-mass yields (with a $^{22}$Ne neutron source) alone do not reproduce the light-to-heavy $s$-element ratios for M4 and M22, and that a small contribution from models with a $^{13}$C pocket is required. With our assumption that $^{13}$C pockets form for initial masses below a transition range between 3.0 and 3.5 M$_\odot$, we match the light-to-heavy s-element ratio in the s-process residual of M22 and predict a minimum enrichment timescale of between 240 and 360 Myr. Our predicted value is consistent with the 300 Myr upper limit age difference between the two groups derived from isochrone fitting.

We present CN and CH indices and CaII triplet metallicities for 34 giant stars and chemical abundances for 33 elements in 14 giants in the globular cluster M2. Assuming the program stars are cluster members, our analysis reveals (i) an extreme variation in CN and CH line strengths, (ii) a metallicity dispersion with a dominant peak at [Fe/H] = -1.7 and smaller peaks at -1.5 and -1.0, (iii) star-to-star abundance variations and correlations for the light elements O, Na, Al and Si and (iv) a large (and possibly bimodal) distribution in the abundances of all elements produced mainly via the s-process in solar system material. Following Roederer et al. (2011), we define two groups of stars, "r+s" and "r-only", and subtract the average abundances of the latter from the former group to obtain a "s-process residual". This s-process residual is remarkably similar to that found in M22 and in M4 despite the range in metallicity covered by these three systems. With recent studies identifying a double subgiant branch in M2 and a dispersion in Sr and Ba abundances, our spectroscopic analysis confirms that this globular cluster has experienced a complex formation history with similarities to M22, NGC 1851 and omega Centauri.

Augmented Reality consists of merging live images with virtual layers of information. The rapid growth in the popularity of smartphones and tablets over recent years has provided a large base of potential users of Augmented Reality technology, and virtual layers of information can now be attached to a wide variety of physical objects. In this article, we explore the potential of Augmented Reality for astrophysical research with two distinct experiments: (1) Augmented Posters and (2) Augmented Articles. We demonstrate that the emerging technology of Augmented Reality can already be used and implemented without expert knowledge using currently available apps. Our experiments highlight the potential of Augmented Reality to improve the communication of scientific results in the field of astrophysics. We also present feedback gathered from the Australian astrophysics community that reveals evidence of some interest in this technology by astronomers who experimented with Augmented Posters. In addition, we discuss possible future trends for Augmented Reality applications in astrophysics, and explore the current limitations associated with the technology. This Augmented Article, the first of its kind, is designed to allow the reader to directly experiment with this technology.

Motivated by unexplained observations of low sulphur abundances in planetary nebulae (PNe) and the PG1159 class of post asymptotic giant branch (AGB) stars, we investigate the possibility that sulphur may be destroyed by nucleosynthetic processes in low-to-intermediate mass stars during stellar evolution. We use a 3 Msun, Z=0.01 evolutionary sequence to examine the consequences of high and low reaction rate estimates of neutron captures onto sulphur and neighbouring elements. In addition, we have tested high and low rates for the neutron producing reactions C13(alpha,n)O16 and Ne22(alpha,n)Mg25. We vary the mass width of a partially mixed zone (PMZ), which is responsible for the formation of a C13 pocket and is the site of the C13(alpha,n)O16 neutron source. We test PMZ masses from zero up to an extreme upper limit of the entire He-intershell mass at 10^-2 Msun. We find that the alternative reaction rates and variations to the partially mixed zone have almost no effect on surface sulphur abundances and do not reproduce the anomaly. To understand the effect of initial mass on our conclusions, 1.8 Msun and 6 Msun evolutionary sequences are also tested with similar results for sulphur abundances. We are able to set a constraint on the size of the PMZ, as PMZ sizes that are greater than half of the He-intershell mass (in the 3 Msun model) are excluded by comparison with neon abundances in planetary nebulae. We compare the 1.8 Msun model's intershell abundances with observations of PG1159-035, whose surface abundances are thought to reflect the intershell composition of a progenitor AGB star. We find general agreement between the patterns of F, Ne, Si, P, and Fe abundances and a very large discrepancy for sulphur where our model predicts abundances that are 30-40 times higher than is observed in the star.