The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

Millisecond x-ray pulsars have weak magnetic dipole moments of $\sim 10^{16}$\,T\,m$^3$ compared to ordinary X-ray pulsars with dipole moments of $10^{20}$\,T\,m$^3$. For this reason a surrounding accretion disc can extend closer to the neutron star, and thus reach a higher temperature, at which the opacity is dominated by electron scattering and radiation pressure is strong. We compute the self-similar structure of such a geometrically thin axisymmetric accretion disc with an internal dynamo. Such models produce significantly stronger torques on the neutron star than models without dynamos, and can explain the strong spin variations in some millisecond X-ray pulsars.

Aims: We determine the steady-state of an axisymmetric thin accretion disc with an internal dynamo around a magnetised star. Methods: Starting from the vertically integrated equations of magnetohydrodynamics we derive a single ordinary differential equation for a thin accretion disc around a massive magnetic dipole and integrate this equation numerically from the outside inwards. Results: Our numerical solution shows that the torque between the star and the accretion disc is dominated by the contribution from the dynamo in the disc. The location of the inner edge of the accretion disc varies between $R_{\rm A}$ and $10R_{\rm A}$ depending mainly on the strength and direction of the magnetic field generated by the dynamo in the disc

BP Cru is a well known high-mass X-ray binary composed of a late B hypergiant (Wray 977) and a neutron star, also observed as the X-ray pulsar GX 301-2. No information about emission from BP Cru in other bands than X-rays and optical has been reported to date in the literature, though massive X-ray binaries containing black holes can have radio emission from a jet. In order to assess the presence of a radio jet, we searched for radio emission towards BP Cru using the Australia Compact Array Telescope during a survey for radio emission from Be/X-ray transients. We probed the 41.5d orbit of BP Cru with the Australia Telescope Compact Array not only close to periastron but also close to apastron. BP Cru was clearly detected in our data on 4, possibly 6, of 12 occasions at 4.8 and 8.6 GHz. Our data suggest that the spectral index of the radio emission is modulated either by the X-ray flux or the orbital phase of the system. We propose that the radio emission of BP Cru probably arises from two components: a persistent component, coming from the mass donor Wray 977, and a periodic component connected to the accretion onto the neutron star, possibly coming from a (weak and short lived) jet.

Context. The high-frequency quasi-periodic oscillations (HF QPOs) in neutron star and stellar-mass black hole X-ray binaries may be the result of a resonance between the radial and vertical epicyclic oscillations in strong gravity. Aims. In this paper we investigate the resonant coupling between the epicyclic modes in a torus in a strong gravitational field. Methods. We perform numerical simulations of axisymmetric constant angular momentum tori in the pseudo-Newtonian potential. The epicyclic motion is excited by adding a constant radial velocity to the torus. Results. We verify that slender tori perform epicyclic motions at the frequencies of free particles, but the epicyclic frequencies decrease as the tori grow thicker. More importantly, and in contrast to previous numerical studies, we do not find a coupling between the radial and vertical epicyclic motions. The appearance of other modes than the radial epicyclic motion in our simulations is rather due to small numerical deviations from exact equilibrium in the initial state of our torus. Conclusions. We find that there is no pressure coupling between the two axisymmetric epicyclic modes as long as the torus is symmetric with respect to the equatorial plane. However we also find that there are other modes in the disc that may be more attractive for explaining the HF QPOs.

We study epicyclic oscillations of fluid tori around black holes (in the Paczynski-Wiita potential) and derive exact analytic expressions for their radial and vertical eigenfrequencies nu_r and nu_z to second-order accuracy in the width of the torus. We prove that pressure effects make the eigenfrequencies smaller than those for free particles. However, the particular ratio nu_z/nu_r=3/2, which is important for the theory of high-frequency quasi-periodic oscillations (QPOs), occurs when the fluid tori epicyclic frequencies nu_r and nu_z are about 15% higher than the ones corresponding to free particles. Our results therefore suggest that previous estimates of black hole spins from QPOs have produced values that are too high.

We study the evolution of non-linear circularly polarised Alfven waves by solving numerically the time-dependent equations of magnetohydrodynamics (MHD) in one dimension. We examine the behaviour of the waves and find that different physical mechanisms are relevant in different ranges of beta. In a low beta plasma the wave may undergo a parametric decay. This is because the wave excites a density enhancement that travels slower than the wave itself and thus interacts with the wave. When beta is greater or equal to one the density enhancement does not interact with the wave and no decay takes place, instead the Alfven wave is reflected against the density enhancement. The reflection zone propagates with the speed 1/n of the Alfven speed. Because of that the magnetic flux is conserved which results in an amplification of the oscillating magnetic field by a factor 1/n. We find that n depends on beta, and that in particular it is less or equal to one for values of beta ~ 1 and is greater or equal to one for beta >> 1. We discuss the relevance of these mechanisms to the acceleration of the solar wind, and the triggering of MHD turbulence in the polar wind region. In particular these simulations can explain the presence of inward propagating Alfven waves in the solar corona.

We study propagating Alfven waves by solving the time-dependent equations of magnetohydrodynamics (MHD) in one dimension numerically. In a homogeneous medium the circularly polarized Alfven wave is an exact solution of the ideal MHD equations, and therefore it does not suffer from any dissipation. A high-amplitude linearly polarized Alfven wave, on the other hand, steepens and form current sheets, in which the Poynting flux is lost. In a stratified medium, however, a high-amplitude circularly polarized Alfven wave can also lose a significant fraction of its Poynting flux.

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $\alpha_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.

I show in this paper that two types of magnetic torques can appear in the interaction between an accretion disc and a magnetic accretor. There is the well-known torque resulting from the difference in angular velocity between the accretion disc and the star, but in addition there is a torque coming from the interaction between the stellar magnetic field and the disc's own magnetic field. The latter form of magnetic torque decreases in strength more slowly with increasing radius, and will therefore dominate at large radii. The direction of the disc field is not determined by the difference in angular velocity between the star and the disc as in the Ghosh & Lamb model, but is rather a free parameter. The magnetic torque may therefore either spin up or spin down the star, and the torque changes sign if the magnetic field in the disc reverses. I suggest that this mechanism can explain the torque reversals that have been observed in some disc-fed X-ray pulsars.

We present an one-dimensional numerical study of Alfven waves propagating along a radial magnetic field. Neglecting losses, any spherical Alfven wave, no matter how small its initial amplitude is, becomes nonlinear at sufficiently large radii. From previous simulations of Alfven waves in plane parallel atmospheres we did expect the waves to steepen and produce current sheets in the nonlinear regime, which was confirmed by our new calculations. On the other hand we did find that even the least nonlinear waves were damped out almost completely before 10 solar radii. A damping of that kind is required by models of Alfven wave-driven winds from old low-mass stars as these winds are mainly accelerated within a few stellar radii.