We present three dimensional smoothed particle hydrodynamics (SPH) calculations of irradiation-driven warping of accretion discs. Initially unwarped planar discs are unstable to the radiation reaction when the disc is illuminated by a central radiation source. The disc warps and tilts and precesses slowly in a retrograde direction; its shape continuously flexes in response to the changing orientation of the Roche potential. We simulate ten systems: eight X-ray binaries, one cataclysmic variable (CV), and a `generic' low mass X-ray binary (LMXB). We adopt system parameters from observations and tune a single parameter: our model X-ray luminosity ($L_{*}$) to reproduce the observed or inferred super-orbital periods. Without exception, across a wide range of parameter space, we find an astonishingly good match between the observed $L_{X}$ and the model $L_{*}$. We conclude irradiation-driven warping is the mechanism underlying the long periods in X-ray binaries. Our Her X-1 simulation simultaneously reproduces the observed $L_{X}$, the "main-" and "short-high" X-ray states and the orbital inclination. Our simulations of SS 433 give a maximum warp angle of $18.6^{\circ}$, a good match to the cone traced by the jets, but this angle is reached only in the outer disc. In all cases, the overall disc tilt is less than $\degrees{13}$ and the maximum disc warp is less than and or equal to $\degrees{21}$.

(Abridged) We use 3D SPH calculations with higher resolution, as well as with more realistic viscosity and sound-speed prescriptions than previous work to examine the eccentric instability which underlies the superhump phenomenon in semi-detached binaries. We illustrate the importance of the two-armed spiral mode in the generation of superhumps. Differential motions in the fluid disc cause converging flows which lead to strong spiral shocks once each superhump cycle. The dissipation associated with these shocks powers the superhump. We compare 2D and 3D results, and conclude that 3D simulations are necessary to faithfully simulate the disc dynamics. We ran our simulations for unprecedented durations, so that an eccentric equilibrium is established except at high mass ratios where the growth rate of the instability is very low. Our improved simulations give a closer match to the observed relationship between superhump period excess and binary mass ratio than previous numerical work. The observed black hole X-ray transient superhumpers appear to have systematically lower disc precession rates than the cataclysmic variables. This could be due to higher disc temperatures and thicknesses. The modulation in total viscous dissipation on the superhump period is overwhelmingly from the region of the disc within the 3:1 resonance radius. As the eccentric instability develops, the viscous torques are enhanced, and the disc consequently adjusts to a new equilibrium state, as suggested in the thermal-tidal instability model. We quantify this enhancement in the viscosity, which is ~10 per cent for q=0.08. We characterise the eccentricity distributions in our accretion discs, and show that the entire body of the disc partakes in the eccentricity.

We present three dimensional smoothed particle hydrodynamics (SPH) calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and in general the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios q ~ 0.1 a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (c.f. X 1916-053). For less extreme mass ratios the entire disc tilts out of the orbital plane (c.f. Her X-1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction we obtain both positive and negative superhumps simultaneously in the dissipation light curve (c.f. V603 Aql).

We present results from three XMM-Newton observations of the M31 low mass X-ray binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-hr period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anti-correlated with source intensity. However, our new observations do not favour a strict intensity dependance, but rather suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio <~ 0.3 (period <~ 4 hr), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of ~11 P_orb, and prograde disc precession on a period of ~29 P_orb. This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on a 81+/-3 hr cycle.

The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that the dip depth was anticorrelated with the X-ray luminosity. However, we present results from three new XMM-Newton observations that suggest that the evolution of dipping is instead due to precession of the accretion disc. Such precession is expected in neutron star LMXBs with mass ratios <0.3 (i.e. with orbital periods <4 hr), such as the Galactic dipping LMXB 4U 1916-053. We simulated the accretion disc of Bo 158 using cutting-edge 3D smoothed particle hydrodynamics (SPH), and using the observed parameters. Our results show disc variability on two time-scales. The disc precesses in a prograde direction on a period of 81+/-3 hr. Also, a radiatively-driven disc warp is present in the inner disc, which undergoes retrograde precesson on a \~31 hr period. From the system geometry, we conclude that the dipping evolution is driven by the disc precession. Hence we predict that the dipping behaviour repeats on a ~81 hr cycle.

In the proceedings of this, and of several recent close binary conferences, there have been several contributions describing smoothed particle hydrodynamics simulations of accretion disks. It is apposite therefore to review the numerical scheme itself with emphasis on its advantages for disk modelling, and the methods used for modelling viscous processes.

Two dimensional SPH simulations of a precessing accretion disc in a q=0.1 binary system (such as XTE J1118+480) reveal complex and continuously varying shape, kinematics, and dissipation. The stream-disc impact region and disc spiral density waves are prominent sources of energy dissipation.The dissipated energy is modulated on the period P_{sh} = ({P_{orb}}^{-1}-{P_{prec}}^{-1}^{-1} with which the orientation of the disc relative to the mass donor repeats. This superhump modulation in dissipation energy has a variation in amplitude of ~10% relative to the total dissipation energy and evolves, repeating exactly only after a full disc precession cycle. A sharp component in the light curve is associated with centrifugally expelled material falling back and impacting the disc. Synthetic trailed spectrograms reveal two distinct "S-wave" features, produced respectively by the stream gas and the disc gas at the stream-disc impact shock. These S-waves are non-sinusoidal, and evolve with disc precession phase. We identify the spiral density wave emission in the trailed spectrogram. Instantaneous Doppler maps show how the stream impact moves in velocity space during an orbit. In our maximum entropy Doppler tomogram the stream impact region emission is distorted, and the spiral density wave emission is uppressed. A significant radial velocity modulation of the whole line profile occurs on the disc precession period. We compare our SPH simulation with a simple 3D model: the former is appropriate for comparison with emission lines while the latter is preferable for skewed absorption lines from precessing discs.