The extreme colliding-wind system Apep : resolved imagery of the central binary and dust plume in the infrared

The recent discovery of a spectacular dust plume in the system 2XMM J160050.7–514245 (referred to as ‘Apep’) suggested a physical origin in a colliding-wind binary by way of the ‘Pinwheel’ mechanism. Observational data pointed to a hierarchical triple-star system, however, several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf–Rayet binary: a system that has so far remained spatially unresolved. Here, we present an updated proper motion study deriving the expansion speed of Apep’s dust plume over a 2-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high angular resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf–Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf–Rayet binary and lending further support to the anisotropic wind model

Variability of M giant stars based on Kepler photometry: general characteristics

M giants are among the longest-period pulsating stars which is why their studies were traditionally restricted to analyses of low-precision visual observations, and more recently, accurate ground-based data. Here we present an overview of M giant variability on a wide range of time-scales (hours to years), based on analysis of thirteen quarters of Kepler long-cadence observations (one point per every 29.4 minutes), with a total time-span of over 1000 days. About two-thirds of the sample stars have been selected from the ASAS-North survey of the Kepler field, with the rest supplemented from a randomly chosen M giant control sample. We first describe the correction of the light curves from different quarters, which was found to be essential. We use Fourier analysis to calculate multiple frequencies for all stars in the sample. Over 50 stars show a relatively strong signal with a period equal to the Kepler-year and a characteristic phase dependence across the whole field-of-view. We interpret this as a so far unidentified systematic effect in the Kepler data. We discuss the presence of regular patterns in the distribution of multiple periodicities and amplitudes. In the period-amplitude plane we find that it is possible to distinguish between solar-like oscillations and larger amplitude pulsations which are characteristic for Mira/SR stars. This may indicate the region of the transition between two types of oscillations as we move upward along the giant branch.Comment: 12 pages, 13 figures, accepted for publication in MNRAS. The normalized light curves are available upon reques

Calibration and Stokes Imaging with Full Embedded Element Primary Beam Model for the Murchison Widefield Array

15 pages, 11 figures. Accepted for publication in PASA. © Astronomical Society of Australia 2017The Murchison Widefield Array (MWA), located in Western Australia, is one of the low-frequency precursors of the international Square Kilometre Array (SKA) project. In addition to pursuing its own ambitious science program, it is also a testbed for wide range of future SKA activities ranging from hardware, software to data analysis. The key science programs for the MWA and SKA require very high dynamic ranges, which challenges calibration and imaging systems. Correct calibration of the instrument and accurate measurements of source flux densities and polarisations require precise characterisation of the telescope's primary beam. Recent results from the MWA GaLactic Extragalactic All-sky MWA (GLEAM) survey show that the previously implemented Average Embedded Element (AEE) model still leaves residual polarisations errors of up to 10-20 % in Stokes Q. We present a new simulation-based Full Embedded Element (FEE) model which is the most rigorous realisation yet of the MWA's primary beam model. It enables efficient calculation of the MWA beam response in arbitrary directions without necessity of spatial interpolation. In the new model, every dipole in the MWA tile (4 x 4 bow-tie dipoles) is simulated separately, taking into account all mutual coupling, ground screen and soil effects, and therefore accounts for the different properties of the individual dipoles within a tile. We have applied the FEE beam model to GLEAM observations at 200 - 231 MHz and used false Stokes parameter leakage as a metric to compare the models. We have determined that the FEE model reduced the magnitude and declination-dependent behaviour of false polarisation in Stokes Q and V while retaining low levels of false polarisation in Stokes U.Peer reviewedFinal Accepted Versio

The Spectral Energy Distribution of Powerful Starburst Galaxies I: Modelling the Radio Continuum

We have acquired radio continuum data between 70\,MHz and 48\,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500\,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500\,MHz the radio-continuum at low frequency ($\nu < 200$\,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $\alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution

The Spectral Energy Distribution of Powerful Starburst Galaxies I : Modelling the Radio Continuum

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. © 2018 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.We have acquired radio-continuum data between 70MHz and 48 GHz for a sample of 19 southern starburst galaxies at moderate redshifts (0.067 < z < 0.227) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework, we find the radio continuum is rarely characterized well by a single power law, instead often exhibiting lowfrequency turnovers below 500 MHz, steepening at mid to high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 and 500MHz the radio continuum at low frequency (v < 200 MHz) could be overestimated by upwards of a factor of 12 if a simple power-law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be α = -1.06, which is steeper than the canonical value of -0.8 for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.Peer reviewe

Murchison Widefield Array and XMM-Newton observations of the Galactic supernova remnant G5.9+3.1

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$.Comment: Accepted for publication in A&

The SAMI Galaxy Survey: The Low-Redshift Stellar Mass Tully-Fisher Relation

We investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use 'simulated slits' to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis

Signatures of star-planet interactions across the electromagnetic spectrum

Stars and planetary system

The Murchison Widefield Array Transients Survey (MWATS). A search for low frequency variability in a bright Southern hemisphere sample

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year.Comment: 19 pages, 11 figure

Low Radio Frequency Observations and Spectral Modelling of the Remnant of Supernova 1987A

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. ©: 2016 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.We present Murchison Widefield Array observations of the supernova remnant (SNR) 1987A between 72 and 230 MHz, representing the lowest frequency observations of the source to date. This large lever arm in frequency space constrains the properties of the circumstellar medium created by the progenitor of SNR 1987A when it was in its red supergiant phase. As of late-2013, the radio spectrum of SNR 1987A between 72 MHz and 8.64 GHz does not show any deviation from a non-thermal power-law with a spectral index of $-0.74 \pm 0.02$. This spectral index is consistent with that derived at higher frequencies, beneath 100 GHz, and with a shock in its adiabatic phase. A spectral turnover due to free-free absorption by the circumstellar medium has to occur below 72 MHz, which places upper limits on the optical depth of $\leq$ 0.1 at a reference frequency of 72 MHz, emission measure of $\lesssim$ 13,000 cm$^{-6}$ pc, and an electron density of $\lesssim$ 110 cm$^{-3}$. This upper limit on the electron density is consistent with the detection of prompt radio emission and models of the X-ray emission from the supernova. The electron density upper limit implies that some hydrodynamic simulations derived a red supergiant mass loss rate that is too high, or a wind velocity that is too low. The mass loss rate of $\sim 5 \times 10^{-6}$ $M_{\odot}$ yr$^{-1}$ and wind velocity of 10 km s$^{-1}$ obtained from optical observations are consistent with our upper limits, predicting a current turnover frequency due to free-free absorption between 5 and 60 MHz.Peer reviewedFinal Published versio