The MeerKAT telescope represents an outstanding opportunity for radio pulsar timing science with its unique combination of a large collecting area and aperture efficiency (effective area $\sim$7500 m$^2$), system temperature ($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and ability to form up to four sub-arrays. The MeerTime project is a five-year program on the MeerKAT array by an international consortium that will regularly time over 1000 radio pulsars to perform tests of relativistic gravity, search for the gravitational wave signature induced by supermassive black hole binaries in the timing residuals of millisecond pulsars, explore the interiors of neutron stars through a pulsar glitch monitoring programme, explore the origin and evolution of binary pulsars, monitor the swarms of pulsars that inhabit globular clusters and monitor radio magnetars. In addition to these primary programmes, over 1000 pulsars will have their arrival times monitored and the data made immediately public. The MeerTime pulsar backend comprises two server-class machines each of which possess four Graphics Processing Units. Up to four pulsars can be coherently dedispersed simultaneously up to dispersion measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits format. The MeerTime backend will be capable of producing coherently dedispersed filterbank data for timing multiple pulsars in the cores of globular clusters that is useful for pulsar searches of tied array beams. All MeerTime data will ultimately be made available for public use, and any published results will include the arrival times and profiles used in the results.

This paper gives a complete characterization of resonant orbits in a Kerr spacetime. A resonant orbit is defined as a geodesic for which the longitudinal and radial orbital frequencies are commensurate. Our analysis is based on expressing the resonance condition in its most symmetric form using Carlson's integrals. We provide a number of concise formulae for the dependence of resonances on the system parameters. Resonant effects may be observable during the in-spiral of a compact object into a super-massive black hole. When the slowly evolving orbital frequencies pass through a series of low-order resonances, rapid changes in the orbital parameters could produce measurable phase shifts in the emitted gravitational radiation (GW). Resonant orbits may also capture dust leading to electromagnetic emission. The KAM theorem indicates that, low order resonant orbits demarcate the regions where the onset of chaos could occur around a perturbed black-hole. We find that the 1/2 and 2/3 resonances occur at ~4 and 5.4 Schwarzschild radii (Rs) from the event horizon. For compact object in-spirals around super-massive black holes, this region lies within the sensitivity band of space-based GW detectors. For Sgr A*, length scales of ~41 and 55 microarcseconds and timescales of 50 and 79 min respectively should be associated with resonant effects, if Sgr A* is non-spinning. Spin decreases these values by up to ~32% and ~28%. These length-scales are potentially resolvable with VLBI measurements. We find that all low-order resonances are localized to the strong field region r < 50 Rs. This fact guarantees the validity of using approximations based on averaging to model the frequency evolution of a test object in region 50 Rs <r <1000 Rs. The systematic determination of the multipole moments of the central object by observing the orbit of a pulsar, free of chaotic effects, is thus possible.

We compute the length and timescales associated with resonant orbits in the Kerr Metric for all orbital and spin parameters. Resonance induced effects are potentially observable when the Event Horizon telescope resolves the inner structure of Sgr A*, space-based gravitational wave detectors record phase-shifts in the waveform during the resonant passage of a neutron star as it spirals in to the black hole and in the frequencies of quasi periodic oscillations for accreting black holes. The onset of geodesic chaos for non-Kerr spacetimes should occur at the resonance locations quantified here.

When one splits spacetime into space plus time, the spacetime curvature (Weyl tensor) gets split into an "electric" part E_{jk} that describes tidal gravity and a "magnetic" part B_{jk} that describes differential dragging of inertial frames. We introduce tools for visualizing B_{jk} (frame-drag vortex lines, their vorticity, and vortexes) and E_{jk} (tidal tendex lines, their tendicity, and tendexes), and also visualizations of a black-hole horizon's (scalar) vorticity and tendicity. We use these tools to elucidate the nonlinear dynamics of curved spacetime in merging black-hole binaries.

We explore the prospects for Advanced LIGO to detect gravitational waves from neutron stars and stellar mass black holes spiraling into intermediate-mass ($M\sim 50 M_\odot$ to $350 M_\odot$) black holes. We estimate an event rate for such \emph{intermediate-mass-ratio inspirals} (IMRIs) of up to $\sim 10$--$30 \mathrm{yr}^{-1}$. Our numerical simulations show that if the central body is not a black hole but its metric is stationary, axisymmetric, reflection symmetric and asymptotically flat then the waves will likely be tri-periodic, as for a black hole. We report generalizations of a theorem due to Ryan (1995) which suggest that the evolutions of the waves' three fundamental frequencies and of the complex amplitudes of their spectral components encode (in principle) a full map of the central body's metric, full details of the energy and angular momentum exchange between the central body and the orbit, and the time-evolving orbital elements. We estimate that Advanced LIGO can measure or constrain deviations of the central body from a Kerr black hole with modest but interesting accuracy.

We describe briefly the properties of the recently completed Southern African Large Telescope (SALT), along with its first light imager SALTICAM. Using this instrument, we present 4.3 hr of high speed unfiltered photometric observations of the eclipsing polar SDSSJ015543.40+002807.2 with time resolution as short as 112 ms, the highest quality observations of this kind of any polar to date. The system was observed during its high luminosity state. Two accreting poles are clearly seen in the eclipse light curve. The binary system parameters have been constrained: the white dwarf mass is at the low end of the range expected for cataclysmic variables. Correlations between the positions of the accretion regions on or near the surface of the white dwarf and the binary system parameters were established. The sizes of the accretion regions and their relative movement from eclipse to eclipse were estimated: they are typically 4-7 deg depending on the mass of the white dwarf. The potential of these observations will only fully be realised when low state data of the same kind are obtained and the contact phases of the eclipse of the white dwarf are measured.

Two mechanisms for nonlinear mode saturation of the r-mode in neutron stars have been suggested: the parametric instability mechanism involving a small number of modes and the formation of a nearly continuous Kolmogorov-type cascade. Using a network of oscillators constructed from the eigenmodes of a perfect fluid incompressible star, we investigate the transition between the two regimes numerically. Our network includes the 4995 inertial modes up to n<= 30 with 146,998 direct couplings to the r-mode and 1,306,999 couplings with detuning< 0.002 (out of a total of approximately 10^9 possible couplings). The lowest parametric instability thresholds for a range of temperatures are calculated and it is found that the r-mode becomes unstable to modes with 13<n<15. In the undriven, undamped, Hamiltonian version of the network the rate to achieve equipartition is found to be amplitude dependent, reminiscent of the Fermi-Pasta-Ulam problem. More realistic models driven unstable by gravitational radiation and damped by shear viscosity are explored next. A range of damping rates, corresponding to temperatures 10^6K to 10^9K, is considered. Exponential growth of the r-mode is found to cease at small amplitudes, approximately 10^-4. For strongly damped, low temperature models, a few modes dominate the dynamics. The behavior of the r-mode is complicated, but its amplitude is still no larger than about 10^-4 on average. For high temperature, weakly damped models the r-mode feeds energy into a sea of oscillators that achieve approximate equipartition. In this case the r-mode amplitude settles to a value for which the rate to achieve equipartition is approximately the linear instability growth rate.

R-modes of a rotating neutron star are unstable because of the emission of gravitational radiation. We explore the saturation amplitudes of these modes determined by nonlinear mode-mode coupling. Modelling the star as incompressible allows the analytic computation of the coupling coefficients. All couplings up to n=30 are obtained, and analytic values for the shear damping and mode normalization are presented. In a subsequent paper we perform numerical simulations of a large set of coupled modes.

Non-linear interactions among the inertial modes of a rotating fluid can be described by a network of coupled oscillators. We use such a description for an incompressible fluid to study the development of the r-mode instability of rotating neutron stars. A previous hydrodynamical simulation of the r-mode reported the catastrophic decay of large amplitude r-modes. We explain the dynamics and timescale of this decay analytically by means of a single three mode coupling. We argue that at realistic driving and damping rates such large amplitudes will never actually be reached. By numerically integrating a network of nearly 5000 coupled modes, we find that the linear growth of the r-mode ceases before it reaches an amplitude of around 10^(-4). The lowest parametric instability thresholds for the r-mode are calculated and it is found that the r-mode becomes unstable to modes with 13<n<15 if modes up to n=30 are included. Using the network of coupled oscillators, integration times of 10^6 rotational periods are attainable for realistic values of driving and damping rates. Complicated dynamics of the modal amplitudes are observed. The initial development is governed by the three mode coupling with the lowest parametric instability. Subsequently a large number of modes are excited, which greatly decreases the linear growth rate of the r-mode.