White dwarfs are routinely observed to have polluted atmospheres, and sometimes significant infrared excesses, that indicate ongoing accretion of circumstellar dust and rocky debris. Typically this debris is assumed to be in the form of a (circular) disc, and to originate from asteroids that passed close enough to the white dwarf to be pulled apart by tides. However, theoretical considerations suggest that the circularisation of the debris, which initially occupies highly eccentric orbits, is very slow. We therefore hypothesise that the observations may be readily explained by the debris remaining on highly eccentric orbits, and we explore the properties of such debris. For the generic case of an asteroid originating at several au from the white dwarf, we find that all of the tidal debris is always bound to the white dwarf and that the orbital energy distribution of the debris is narrow enough that it executes similar elliptical orbits with only a narrow spread. Assuming that the tidal field of the white dwarf is sufficient to minimise the effects of self-gravity and collisions within the debris, we estimate the time over which the debris spreads into a single elliptical ring, and we generate toy spectra and lightcurves from the initial disruption to late times when the debris distribution is essentially time steady. Finally we speculate on the connection between these simple considerations and the observed properties of these systems, and on additional physical processes that may change this simple picture.

We have made high precision polarimetric observations of the polluted white dwarf G29-38 with the HIgh Precision Polarimetric Instrument 2. The observations were made at two different observatories -- using the 8.1-m Gemini North Telescope and the 3.9-m Anglo AustralianTelescope -- and are consistent with each other. After allowing for a small amount of interstellar polarization, the intrinsic linear polarization of the system is found to be 275.3 +/- 31.9 parts-per-million at a position angle of 90.8 +/- 3.8 degrees in the SDSS g' band. We compare the observed polarization with the predictions of circumstellar disc models. The measured polarization is small in the context of the models we develop which only allows us to place limits on disc inclination and Bond albedo for optically thin disc geometries. In this case either the inclination is near face-on or the albedo is small -- likely in the range 0.05 to 0.15 -- which is in line with other debris disc measurements. A preliminary search for the effects of G29-38's pulsations in the polarization signal produced inconsistent results. This may be caused by beating effects, indicate a clumpy dust distribution, or be a consequence of measurement systematics.

In order to provide an explanation for the unexpected radial brightness distribution of the steady accretion discs seen in nova-like variables, Nixon & Pringle (2019) proposed that the accretion energy is redistributed outwards by means of strong, magnetically driven, surface flows. In this paper we note that the "powerful, rotating disc winds" observed in the soft states of black hole X-ray binaries, and also in the disc around a magnetized neutron star in Her X-1, have the properties of the outflows postulated by Nixon & Pringle to exist in the nova-like variable accretion discs around white dwarfs. The relevant properties are that the flows are not winds, but are, instead, bound flows (traveling at less than the escape velocity) and that the mass fluxes in the flows are a substantial fraction of the accretion rate in the disc.

We present the results of a multi-wavelength follow up campaign for the luminous nuclear transient Gaia16aax, which was first identified in January 2016. The transient is spatially consistent with the nucleus of an active galaxy at z=0.25, hosting a black hole of mass $\rm \sim6\times10^8M_\odot$. The nucleus brightened by more than 1 magnitude in the Gaia G-band over a timescale of less than one year, before fading back to its pre-outburst state over the following three years. The optical spectra of the source show broad Balmer lines similar to the ones present in a pre-outburst spectrum. During the outburst, the $\rm H\alpha$ and $\rm H\beta$ emission lines develop a secondary peak. We also report on the discovery of two transients with similar light curve evolution and spectra: Gaia16aka and Gaia16ajq. We consider possible scenarios to explain the observed outbursts. We exclude that the transient event could be caused by a microlensing event, variable dust absorption or a tidal encounter between a neutron star and a stellar mass black hole in the accretion disk. We consider variability in the accretion flow in the inner part of the disk, or a tidal disruption event of a star $\geq 1 M_{\odot}$ by a rapidly spinning supermassive black hole as the most plausible scenarios. We note that the similarity between the light curves of the three Gaia transients may be a function of the Gaia alerts selection criteria.

In spiral galaxies, the pitch angle, $\alpha$, of the spiral arms is often proposed as a discriminator between theories for the formation of the spiral structure. In Lin-Shu density wave theory, $\alpha$ stays constant in time, being simply a property of the underlying galaxy. In other theories (e.g tidal interaction, self-gravity) it is expected that the arms wind up in time, so that to a first approximation $\cot \alpha \propto t$. For these theories, it would be expected that a sample of galaxies observed at random times should show a uniform distribution of $\cot \alpha$. We show that a recent set of measurements of spiral pitch angles (Yu & Ho 2018) is broadly consistent with this expectation.

In a standard, steady, thin accretion disc, the radial distribution of the dissipation of the accretion energy is determined simply by energy considerations. Here we draw attention to the fact that while the (quasi-)steady discs in dwarf novae in outburst are in agreement with the expected emission distribution, the steady discs in the nova-like variables are not. We note that essentially the only difference between these two sets of discs is the time for which they have been in the high viscosity, high accretion rate state. In such discs, the major process by which angular momentum is transported outwards is MHD turbulence. We speculate that such turbulence gives rise to corona-like structures (here called magnetically controlled zones, or MCZs) which are also able to provide non-negligible angular momentum transport, the magnitude of which depends on the spatial scale $L$ of the magnetic field structures in such zones. For short-lived, high accretion rate discs (such as those in dwarf novae) we expect $L \sim H$ and the MCZ to have little effect. But, with time (such as in the nova-like variables) an inverse cascade in the MHD turbulence enables $L$, and the net effect of the MCZ, to grow. We present a simple toy model which demonstrates that such ideas can provide an explanation for the difference between the dwarf novae and the nova-like variable discs.

The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(\alpha < H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R<\alpha \ll 1$ it evolves diffusively. Here $\alpha$ is the viscosity parameter and $H/R$ the disk aspect ratio. Currently there is no simple set of equations which describes the evolution in both regimes. In this paper we describe a possible solution to this problem and introduce a set of one-dimensional equations that describe the evolution of a warped disk that are applicable in both high and low viscosity regimes for arbitrary tilts, but small warps.

We use well-established observational evidence to draw conclusions about the fundamental nature of the viscosity in accretion discs. To do this, we first summarise the observational evidence for the value of the dimensionless accretion disc viscosity parameter $\alpha$, defined by Shakura & Sunyaev (1973, 1976). We find that, for fully ionized discs, the value of $\alpha$ is readily amenable to reliable estimation and that the observations are consistent with the hypothesis that $\alpha \sim 0.2 - 0.3$. In contrast in discs that are not fully ionized, estimates of the value of $\alpha$ are generally less direct and the values obtained are generally $ < 0.01$ and often $ \ll 0.01$. We conclude that this gives us crucial information about the nature of viscosity in accretion discs. First, in fully ionized discs the strength of the turbulence is always limited by being at most trans-sonic. This implies that it is necessary that credible models of the turbulence reflect this fact. Second, the smaller values of $\alpha$ found for less ionized, and therefore less strongly conducting, discs imply that magnetism plays a dominant role. This provides important observational support for the concept of magneto-rotational instability (MRI) driven hydromagnetic turbulence.

Magnetic fields are observed in star forming regions. However simulations of the late stages of star formation that do not include magnetic fields provide a good fit to the properties of young stars including the initial mass function (IMF) and the multiplicity. We argue here that the simulations that do include magnetic fields are unable to capture the correct physics, in particular the high value of the magnetic Prandtl number, and the low value of the magnetic diffusivity. The artificially high (numerical and uncontrolled) magnetic diffusivity leads to a large magnetic flux pervading the star forming region. We argue further that in reality the dynamics of high magnetic Prandtl number turbulence may lead to local regions of magnetic energy dissipation through reconnection, meaning that the regions of molecular clouds which are forming stars might be essentially free of magnetic fields. Thus the simulations that ignore magnetic fields on the scales on which the properties of stellar masses, stellar multiplicities and planet-forming discs are determined, may be closer to reality than those which include magnetic fields, but can only do so in an unrealistic parameter regime.

We perform simulations of isolated galaxies in order to investigate the likely origin of the spiral structure in M33. In our models, we find that gravitational instabilities in the stars and gas are able to reproduce the observed spiral pattern and velocity field of M33, as seen in HI, and no interaction is required. We also find that the optimum models have high levels of stellar feedback which create large holes similar to those observed in M33, whilst lower levels of feedback tend to produce a large amount of small scale structure, and undisturbed long filaments of high surface density gas, hardly detected in the M33 disc. The gas component appears to have a significant role in producing the structure, so if there is little feedback, both the gas and stars organise into clear spiral arms, likely due to a lower combined $Q$ (using gas and stars), and the ready ability of cold gas to undergo spiral shocks. By contrast models with higher feedback have weaker spiral structure, especially in the stellar component, compared to grand design galaxies. We did not see a large difference in the behaviour of $Q_{stars}$ with most of these models, however, because $Q_{stars}$ stayed relatively constant unless the disc was more strongly unstable. Our models suggest that although the stars produce some underlying spiral structure, this is relatively weak, and the gas physics has a considerable role in producing the large scale structure of the ISM in flocculent spirals.

Supermassive black holes at the centre of galactic nuclei mostly grow in mass through gas accretion over cosmic time. This process also modifies the angular momentum (or spin) of black holes, both in magnitude and in orientation. Despite being often neglected in galaxy formation simulations, spin plays a crucial role in modulating accretion power, driving jet feedback, and determining recoil velocity of coalescing black hole binaries. We present a new accretion model for the moving-mesh code {\sc arepo} that incorporates (i) mass accretion through a thin $\alpha$-disc, and (ii) spin evolution through the Bardeen-Petterson effect. We use a diverse suite of idealised simulations to explore the physical connection between spin evolution and larger scale environment. We find that black holes with mass $\lesssim 10^{7}$ M$_{\odot}$ experience quick alignment with the accretion disc. This favours prolonged phases of spin-up, and the spin direction evolves according to the gas inflow on timescales as short as $\lesssim 100$ Myr, which might explain the observed jet direction distribution in Seyfert galaxies. Heavier black holes ($\gtrsim 10^{8}$ M$_{\odot}$) are instead more sensitive to the local gas kinematic. Here we find a wider distribution in spin magnitudes: spin-ups are favoured if gas inflow maintains a preferential direction, and spin-downs occur for nearly isotropic infall, while the spin direction does not change much over short timescales $\sim 100$ Myr. We therefore conclude that supermassive black holes with masses $\gtrsim 5 \times 10^{8}$ M$_{\odot}$ may be the ideal testbed to determine the main mode of black hole fuelling over cosmic time.

Current planet formation theories provide successful frameworks with which to interpret the array of new observational data in this field. However, each of the two main theories (core accretion, gravitational instability) is unable to explain some key aspects. In many planet formation calculations, it is usual to treat the initial properties of the planet forming disc (mass, radius, etc.) as free parameters. In this paper, we stress the importance of setting the formation of planet forming discs within the context of the formation of the central stars. By exploring the early stages of disc formation, we introduce the concept of the Maximum Mass Solar Nebula (MMSN), as opposed to the oft-used Minimum Mass Solar Nebula (here mmsn). It is evident that almost all protoplanetary discs start their evolution in a strongly self-gravitating state. In agreement with almost all previous work in this area, we conclude that on the scales relevant to planet formation these discs are not gravitationally unstable to gas fragmentation, but instead form strong, transient spiral arms. These spiral arms can act as efficient dust traps allowing the accumulation and subsequent fragmentation of the dust (but not the gas). This phase is likely to populate the disc with relatively large planetesimals on short timescales while the disc is still veiled by a dusty-gas envelope. Crucially, the early formation of large planetesimals overcomes the main barriers remaining within the core accretion model. A prediction of this picture is that essentially all observable protoplanetary discs are already planet hosting.

The large-scale magnetic fields observed in spiral disc galaxies are often thought to result from dynamo action in the disc plane. However, the increasing importance of Faraday depolarization along any line of sight towards the galactic plane suggests that the strongest polarization signal may come from well above (~0.3-1 kpc) this plane, from the vicinity of the warm interstellar medium (WIM)/halo interface. We propose (see also Henriksen & Irwin 2016) that the observed spiral fields (polarization patterns) result from the action of vertical shear on an initially poloidal field. We show that this simple model accounts for the main observed properties of large-scale fields. We speculate as to how current models of optical spiral structure may generate the observed arm/interarm spiral polarization patterns.

Accretion discs are generally warped. If a warp in a disc is too large, the disc can `break' apart into two or more distinct planes, with only tenuous connections between them. Further if an initially planar disc is subject to a strong differential precession, then it can be torn apart into discrete annuli that precess effectively independently. In previous investigations, torque-balance formulae have been used to predict where and when the disc breaks into distinct parts. In this work, focusing on discs with Keplerian rotation and where the shearing motions driving the radial communication of the warp are damped locally by turbulence (the `diffusive' regime), we investigate the stability of warped discs to determine the precise criterion for an isolated warped disc to break. We find and solve the dispersion relation, which in general yields three roots. We provide a comprehensive analysis of this viscous-warp instability and the emergent growth rates and their dependence on disc parameters. The physics of the instability can be understood as a combination of (1) a term which would generally encapsulate the classical Lightman-Eardley instability in planar discs (given by $\partial(\nu\Sigma)/\partial\Sigma < 0$) but is here modified by the warp to include $\partial(\nu_1|\psi|)/\partial|\psi| < 0$ and (2) a similar condition acting on the diffusion of the warp amplitude given in simplified form by $\partial(\nu_2|\psi|)/\partial|\psi| < 0$. We discuss our findings in the context of discs with an imposed precession, and comment on the implications for different astrophysical systems.

The rapidly evolving dust and gas extinction observed towards WD 1145+017 has opened a real-time window onto the mechanisms for destruction-accretion of planetary bodies onto white dwarf stars, and has served to underline the importance of considering the dynamics of dust particles around such objects. Here it is argued that the interaction between (charged) dust grains and the stellar magnetic field is an important ingredient in understanding the physical distribution of infrared emitting particles in the vicinity of such white dwarfs. These ideas are used to suggest a possible model for WD 1145+017 in which the unusual transit shapes are caused by opaque clouds of dust trapped in the stellar magnetosphere. The model can account for the observed transit periodicities if the stellar rotation is near 4.5 h, as the clouds of trapped dust are then located near or within the co-rotation radius. The model requires the surface magnetic field to be at least around some tens of kG. In contrast to the eccentric orbits expected for large planetesimals undergoing tidal disintegration, the orbits of magnetospherically-trapped dust clouds are essentially circular, consistent with the observations.

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of $\beta$ Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain-grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionisation fraction and population levels of carbon and oxygen are followed with the photodissociation region model Cloudy, which is coupled to a dynamical viscous $\alpha$ model. We present new gas observations of $\beta$ Pic, of C I observed with APEX and O I observed with Herschel, and show that these along with published C II and CO observations can all be explained with this new model. Our model requires a viscosity $\alpha$ > 0.1, similar to that found in sufficiently ionised discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionised and dilute medium. This new model can be tested from its predictions for high resolution ALMA observations of C I. We also constrain the water content of the planetesimals in $\beta$ Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.

Laming (1990) predicted that the narrow Balmer line core of the ~3000 km/s shock in the SN 1006 remnant would be significantly polarized due to electron and proton impact polarization. Here, based on deep spectrally resolved polarimetry obtained with the European Southern Observatory (ESO)'s Very Large Telescope (VLT), we report the discovery of polarized line emission of polarization degree approx 1.3 percent with position angle orthogonal to the SNR filament. Correcting for an unpolarized broad line component, the implied narrow line polarization is approx 2.0 percent, close to the predictions of Laming (1990). The predicted polarization is primarily sensitive to shock velocity and post-shock temperature equilibration. By measuring polarization for the SN1006 remnant, we validate and enable a new diagnostic that has important applications in a wide variety of astrophysical situations, such as shocks, intense radiation fields, high energy particle streams and conductive interfaces.

We investigate cloud-cloud collisions, and GMC evolution, in hydrodynamic simulations of isolated galaxies. The simulations include heating and cooling of the ISM, self--gravity and stellar feedback. Over timescales $<5$ Myr most clouds undergo no change, and mergers and splits are found to be typically two body processes, but evolution over longer timescales is more complex and involves a greater fraction of intercloud material. We find that mergers, or collisions, occur every 8-10 Myr (1/15th of an orbit) in a simulation with spiral arms, and once every 28 Myr (1/5th of an orbit) with no imposed spiral arms. Both figures are higher than expected from analytic estimates, as clouds are not uniformly distributed in the galaxy. Thus clouds can be expected to undergo between zero and a few collisions over their lifetime. We present specific examples of cloud--cloud interactions in our results, including synthetic CO maps. We would expect cloud--cloud interactions to be observable, but find they appear to have little or no impact on the ISM. Due to a combination of the clouds' typical geometries, and moderate velocity dispersions, cloud--cloud interactions often better resemble a smaller cloud nudging a larger cloud. Our findings are consistent with the view that spiral arms make little difference to overall star formation rates in galaxies, and we see no evidence that collisions likely produce massive clusters. However, to confirm the outcome of such massive cloud collisions we ideally need higher resolution simulations.

This paper explores how the stochastic accretion of planetesimals onto white dwarfs would be manifested in observations of their atmospheric pollution. Archival observations of pollution levels for unbiased samples of DA and non-DA white dwarfs are used to derive the distribution of accretion rates, confirming that rates become systematically lower as sinking time is decreased, with no discernable dependence on cooling age. The accretion rates expected from planetesimals that are all the same mass (ie, a mono-mass distribution) are explored both analytically and using a Monte Carlo model, quantifying how measured accretion rates inevitably depend on sinking time, since different sinking times probe different times since the last accretion event. However, that dependence is so dramatic that a mono-mass distribution can be excluded. Consideration of accretion from a broad distribution of planetesimal masses uncovers an important conceptual difference: accretion is continuous (rather than stochastic) for planetesimals below a certain mass, and the accretion of such planetesimals determines the rate typically inferred from observations; smaller planetesimals dominate the rates for shorter sinking times. A reasonable fit to the observationally inferred accretion rate distributions is found with model parameters consistent with a collisionally evolved mass distribution up to Pluto-mass, and an underlying accretion rate distribution consistent with that expected from descendants of debris discs of main sequence A stars. With these parameters, while both DA and non-DA white dwarfs accrete from the same broad planetesimal distribution, this model predicts that the pollution seen in DAs is dominated by the continuous accretion of <35km objects, and that in non-DAs by >35km objects (though the dominant size varies between stars by around an order of magnitude from this value)...

The nature of the interaction between low-excitation gas filaments at ~10^4 K, seen in optical line emission, and diffuse X-ray emitting coronal gas at ~10^7 K in the centers of galaxy clusters remains a puzzle. The presence of a strong, empirical correlation between the two gas phases is indicative of a fundamental relationship between them, though as yet of undetermined cause. The cooler filaments, originally thought to have condensed from the hot gas, could also arise from a merger or the disturbance of cool circumnuclear gas by nuclear activity. Here, we have searched for intrinsic line emission polarization in cool core galaxy clusters as a diagnostic of fundamental transport processes. Drawing on developments in solar astrophysics, direct energetic particle impact induced polarization holds the promise to definitively determine the role of collisional processes such as thermal conduction in the ISM physics of galaxy clusters, while providing insight into other highly anisotropic excitation mechanisms such as shocks, intense radiation fields and suprathermal particles. Under certain physical conditions, theoretical calculations predict of order ten percent polarization. Our observations of the filaments in four nearby cool core clusters place stringent upper limits (<0.1%) on the presence of emission line polarization, requiring that if thermal conduction is operative, the thermal gradients are not in the saturated regime. This limit is consistent with theoretical models of the thermal structure of filament interfaces.

In this Letter, we analyse the distributions of stellar ages in Giant Molecular Clouds (GMCs) in spiral arms, inter-arm spurs, and at large galactic radii, where the spiral arms are relatively weak. We use the results of numerical simulations of galaxies, which follow the evolution of GMCs and include star particles where star formation events occur. We find that GMCs in spiral arms tend to have predominantly young (< 10 Myr) stars. By contrast, clouds which are the remainders of spiral arm GMAs that have been sheared into inter-arm GMCs, contain fewer young (< 10 Myr) stars, and more ~20 Myr stars. We also show that clouds which form in the absence of spiral arms, due to local gravitational and thermal instabilities, contain preferentially young stars. We propose the age distributions of stars in GMCs will be a useful diagnostic to test different cloud evolution scenarios, the origin of spiral arms, and the success of numerical models of galactic star formation. We discuss the implications of our results in the context of Galactic and extragalactic molecular clouds.

We present a detailed study of the evolution of GMCs in a galactic disc simulation. We follow individual GMCs (defined in our simulations by a total column density criterion), including their level of star formation, from their formation to dispersal. We find the evolution of GMCs is highly complex. GMCs often form from a combination of smaller clouds and ambient ISM, and similarly disperse by splitting into a number of smaller clouds and ambient ISM. However some clouds emerge as the result of the disruption of a more massive GMC, rather than from the assembly of smaller clouds. Likewise in some cases, clouds accrete onto more massive clouds rather than disperse. Because of the difficulty of determining a precursor or successor of a given GMC, determining GMC histories and lifetimes is highly non-trivial. Using a definition relating to the continuous evolution of a cloud, we obtain lifetimes typically of 4-25 Myr for >10^5 M$_{\odot}$ GMCs, over which time the star formation efficiency is about 1 %. We also relate the lifetime of GMCs to their crossing time. We find that the crossing time is a reasonable measure of the actual lifetime of the cloud, although there is considerable scatter. The scatter is found to be unavoidable because of the complex and varied shapes and dynamics of the clouds. We study cloud dispersal in detail and find both stellar feedback and shear contribute to cloud disruption. We also demonstrate that GMCs do not behave as ridge clouds, rather massive spiral arm GMCs evolve into smaller clouds in inter-arm spurs.

Standard, planar accretion discs operate through a dissipative mechanism, usually thought to be turbulent, and often modelled as a viscosity. This acts to take energy from the radial shear, enabling the flow of mass and angular momentum in the radial direction. In a previous paper we discussed observational evidence for the magnitude of this viscosity, and pointed out discrepancies between these values and those obtained in numerical simulations. In this paper we discuss the observational evidence for the magnitude of the dissipative effects which act in non--planar discs, both to transfer and to eliminate the non--planarity. Estimates based on the model by Ogilvie (1999), which assumes a small--scale, isotropic viscosity, give alignment timescales for fully ionized discs which are apparently too short by a factor of a few compared with observations, although we emphasize that more detailed computations as well as tighter observational constraints are required to verify this conclusion. For discs with low temperature and conductivity, we find that the timescales for disc alignment based on isotropic viscosity are too short by around two orders of magnitude. This large discrepancy suggests that our understanding of viscosity in quiescent discs is currently inadequate.

Motivated by recent observations of short-timescale variations in the infrared emission of circumstellar disks, we propose that coronal mass ejections can remove dust grains on timescales as short as a few days. Continuous monitoring of stellar activity, coupled to infrared observations, can place meaningful constraints on the proposed mechanism.

We analyse the results of four simulations of isolated galaxies: two with a rigid spiral potential of fixed pattern speed, but with different degrees of star-formation induced feedback, one with an axisymmetric galactic potential and one with a `live' self-gravitating stellar component. Since we use a Lagrangian method we are able to select gas that lies within giant molecular clouds (GMCs) at a particular timeframe, and to then study the properties of this gas at earlier and later times. We find that gas which forms GMCs is not typical of the interstellar medium at least 50 Myr before the clouds form and reaches mean densities within an order of magnitude of mean cloud densities by around 10 Myr before. The gas in GMCs takes at least 50 Myr to return to typical ISM gas after dispersal by stellar feedback, and in some cases the gas is never fully recycled. We also present a study of the two-dimensional, vertically-averaged velocity fields within the ISM. We show that the velocity fields corresponding to the shortest timescales (that is, those timescales closest to the immediate formation and dissipation of the clouds) can be readily understood in terms of the various cloud formation and dissipation mechanisms. Properties of the flow patterns can be used to distinguish the processes which drive converging flows (e.g.\ spiral shocks, supernovae) and thus molecular cloud formation, and we note that such properties may be detectable with future observations of nearby galaxies.

A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points to historical, time-averaged rates significantly higher than the currently observed episodes at their DAZ counterparts. The difference between the ongoing, instantaneous accretion rates witnessed at DAZ white dwarfs, which often exceed 1e8 g/s, and those inferred over the past 1e5-1e6 yr for the DBZ stars can be a few orders of magnitude, and therefore must result from high-rate episodes of tens to hundreds of years so they remain undetected to date. This paper explores the likelihood that such brief, intense accretion episodes of gas-phase material can account for existing data. For reasonable assumptions about the circumstellar gas, accretion rates approaching or exceeding 1e15 g/s are possible, similar to rates observed in quiescent cataclysmic variables, and potentially detectable with future x-ray missions or wide-field monitoring facilities. Gaseous debris that is prone to such rapid accretion may be abundant immediately following a tidal disruption event via collisions and sublimation, or if additional bodies impinge upon an extant disk. Particulate disk matter accretes at or near the Poynting-Robertson drag rate for long periods between gas-producing events, consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ stars without infrared excesses have rates consistent with accretion from particulate disks that remain undetected. This overall picture has implications for quasi-steady state models of accretion and the derived chemical composition of asteroidal debris in DBZ white dwarfs.

We perform global time-dependent simulations of an accretion disc around a young stellar object with a dead zone (a region where the magneto-rotational instability cannot drive turbulence because the material is not sufficiently ionised). For infall accretion rates on to the disc of around 10^-7 Msun/yr, dead zones occur if the critical magnetic Reynolds number is larger than about 10^4. We model the collapse of a molecular gas cloud. At early times when the infall accretion rate is high, the disc is thermally ionised and fully turbulent. However, as the infall accretion rate drops, a dead zone may form if the critical magnetic Reynolds number is sufficiently large, otherwise the disc remains fully turbulent. With a dead zone the disc can become unstable to the gravo-magneto instability. The mass of the star grows in large accretion outbursts that may explain FU Orionis events. At late times there is not sufficient mass in the disc for outbursts to occur but the dead zone becomes even more prominent as the disc cools. Large inner dead zones in the later stages of disc evolution may help to explain observations of transition discs with an inner hole.

The physical relationship between low-excitation gas filaments at ~10^4 K, seen in optical line emission, and diffuse X-ray emitting coronal gas at ~10^7 K in the centers of many galaxy clusters is not understood. It is unclear whether the ~10^4 K filaments have cooled and condensed from the ambient hot (~10^7 K) medium or have some other origin such as the infall of cold gas in a merger, or the disturbance of an internal cool reservoir of gas by nuclear activity. Observations of gas at intermediate temperatures (~10^5-10^6 K) can potentially reveal whether the central massive galaxies are gaining cool gas through condensation or losing it through conductive evaporation and hence identify plausible scenarios for transport processes in galaxy cluster gas. Here we present spectroscopic detection of ~10^5 K gas spatially associated with the H-alpha filaments in a central cluster galaxy, M87 in the Virgo Cluster. The measured emission-line fluxes from triply ionized carbon (CIV 1549 A) and singly ionized helium (HeII 1640 A) are consistent with a model in which thermal conduction determines the interaction between hot and cold phases.

Angular momentum is transported outwards through an accretion disc by magnetohydrodynamical (MHD) turbulence thus allowing material to accrete on to the central object. The magneto-rotational instability (MRI) requires a minimum ionisation fraction to drive turbulence in a disc. The inner parts of the disc around a young stellar object are sufficiently hot to be thermally ionised. Further out, cosmic rays ionise the surface layers and a dead zone forms at the mid-plane where the disc is too cool for the MRI to operate. The surface density in the turbulent active layer is often assumed to be constant with radius because the cosmic rays penetrate a constant layer. However, if a critical magnetic Reynolds number, Re_{M,crit}, is used to determine the extent of the dead zone, the surface density in the layer generally increases with radius. For small critical magnetic Reynolds number of order 1, the constant layer approximation may be a reasonable fit. However, MHD simulations suggest the critical magnetic Reynolds number may be much larger, of order 10^4. Analytical fits for the surface density in the magnetic active layer show that \Sigma_m \propto Re_{M,crit}^{-2} R^{9/2} T^{\,2}$, at temperature T and radius R, are a good fit for higher critical magnetic Reynolds number. For the metallicity variation between our galaxy, the LMC and the SMC, there should be no significant difference in the extent of the dead zone. Observations suggest an increase in the lifetime of the disc with decreasing metallicity that cannot be explained by the dead zone structure (ignoring possible differences in dust abundances).

We propose two new means of identifying the main class of progenitors of Type Ia supernovae--single or double degenerate: (i) If the range of supernova properties is significantly determined by the range of viewing angles of non-spherically symmetric explosions, then the nature of the correlation between polarization and another property (for example, the velocity gradient) can be used to determine the geometry of the asymmetry and hence the nature of the progenitor, and (ii) in the double- but not in the single-degenerate case, the range in the observed properties (e.g., velocity gradients) is likely to increase with the amount of carbon seen in the ejecta.

We show that for young stars which are still accreting and for which measurements of stellar age, disc mass and accretion rate are available, nominal disc age (Disc Age = Disc Mass / Accretion Rate) is approximately equal to the stellar age, at least within the considerable observational scatter. We then consider theoretical models of proto-stellar discs through analytic and numerical models. A variety of viscosity prescriptions including empirical power laws, magnetohydrodynamic turbulence and gravitational instability were considered within models describing the disc phenomena of dead zones, photoevaporation and planet formation. These models are generally poor fits to the observational data, showing values of 'Disc Age' which are too high by factors of 3 - 10. We then ask whether a systematic error in the measurement of one of the observational quantities might provide a reasonable explanation for this discrepancy. We show that for the observed systems only disc mass shows a systematic dependence on the value of 'Disc Age / Stellar Age' and we note that a systematic underestimate of the value of disc mass by a factor of around 3 - 5, would account for the discrepancy between theory and observations.

We study the stability of disc accretion in the recurrent nova RS Ophiuchi. We construct a one-dimensional time-dependent model of the binary-disc system, which includes viscous heating and radiative cooling and a self-consistent treatment of the binary potential. We find that the extended accretion disc in this system is always unstable to the thermal-viscous instability, and undergoes repeated disc outbursts on ~10-20yr time-scales. This is similar to the recurrence time-scale of observed outbursts in the RS Oph system, but we show that the disc's accretion luminosity during outburst is insufficient to explain the observed outbursts. We explore a range of models, and find that in most cases the accretion rate during outbursts reaches or exceeds the critical accretion rate for stable nuclear burning on the white dwarf surface. Consequently we suggest that a surface nuclear burning triggered by disc instability may be responsible for the observed outbursts. This allows the white dwarf mass to grow over time, and we suggest that disc instability in RS Oph and similar systems may represent a path to Type Ia supernovae.

We consider the interaction between a binary system (e.g. two supermassive black holes or two stars) and an external accretion disc with misaligned angular momentum. This situation occurs in galaxy merger events involving supermassive black holes, and in the formation of stellar--mass binaries in star clusters. We work out the gravitational torque between the binary and disc, and show that their angular momenta J_b, J_d stably counteralign if their initial orientation is sufficiently retrograde, specifically if the angle theta between them obeys cos(theta) < -J_d/2J_b, on a time short compared with the mass gain time of the central accretor(s). The magnitude J_b remains unchanged in this process. Counteralignment can promote the rapid merger of supermassive black hole binaries, and possibly the formation of coplanar but retrograde planets around stars in binary systems.

We perform calculations of isolated disc galaxies to investigate how the properties of the ISM, the nature of molecular clouds, and the global star formation rate depend on the level of stellar feedback. We adopt a simple physical model, which includes a galactic potential, a standard cooling and heating prescription of the ISM, and self gravity of the gas. Stellar feedback is implemented by injecting energy into dense, gravitationally collapsing gas, but is independent of the Schmidt-Kennicutt relation. We obtain fractions of gas, and filling factors for different phases of the ISM in reasonable ageement with observations. Supernovae are found to be vital to reproduce the scale heights of the different components of the ISM, and velocity dispersions. The GMCs formed in the simulations display mass spectra similar to the observations, their normalisation dependent on the level of feedback. We find ~40 per cent of the clouds exhibit retrograde rotation, induced by cloud-cloud collisions. The star formation rates we obtain are in good agreement with the observed Schmidt-Kennicutt relation, and are not strongly dependent on the star formation efficiency we assume, being largely self regulated by the feedback. We also investigate the effect of spiral structure by comparing calculations with and without the spiral component of the potential. The main difference with a spiral potential is that more massive GMCs are able to accumulate in the spiral arms. Thus we are able to reproduce massive GMCs, and the spurs seen in many grand design galaxies, even with stellar feedback. The presence of the spiral potential does not have an explicit effect on the star formation rate, but can increase the star formation rate indirectly by enabling the formation of long-lived, strongly bound clouds.

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed edge on. Transitions between the two implies that the discs may be warped and precessing. Type II X-ray outbursts are thought to occur when the warped disc interacts with the fast stellar wind. We suggest that tides from a misaligned companion neutron star can cause the observed effects. We make numerical models of a Be star decretion disc in which the spin of the Be star is misaligned with the orbital axis of a neutron star companion. Tidal torques from the neutron star truncate the disc at a radius small enough that the neutron star orbit does not intersect the disc unless the eccentricity or misalignment is very large. A magnetic torque from the Be star that is largest at the equator, where the rotation is fastest, is approximated by an inner boundary condition. There are large oscillations in the mass and inclination of the disc as it moves towards a steady state. These large variations may explain the observed changes from Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find the tidal timescale on which the disc warps, precesses and reaches a steady state to be around a year up to a few hundred years. If present, the oscillations in mass and disc inclination occur on a fraction of this timescale depending on the orbital parameters of the binary. The timescales associated with the tidal torque for observed Be star binaries suggest that these effects are important in all but the longest period binaries.

The most recent observational evidence seems to indicate that giant molecular clouds are predominantly gravitationally unbound objects. In this paper we show that this is a natural consequence of a scenario in which cloud-cloud collisions and stellar feedback regulate the internal velocity dispersion of the gas, and so prevent global gravitational forces from becoming dominant. Thus, while the molecular gas is for the most part gravitationally unbound, local regions within the denser parts of the gas (within the clouds) do become bound and are able to form stars. We find that the observations, in terms of distributions of virial parameters and cloud structures, can be well modelled provided that the star formation efficiency in these bound regions is of order 5 - 10 percent. We also find that in this picture the constituent gas of individual molecular clouds changes over relatively short time scales, typically a few Myr.

We investigate whether a circumbinary gas disc can coalesce a supermassive black hole binary system in the centre of a galaxy. This is known to be problematic for a prograde disc. We show that in contrast, interaction with a retrograde circumbinary disc is considerably more effective in shrinking the binary because there are no orbital resonances. The binary directly absorbs negative angular momentum from the circumbinary disc by capturing gas into a disc around the secondary black hole, or discs around both holes if the binary mass ratio is close to unity. In many cases the binary orbit becomes eccentric, shortening the pericentre distance as the eccentricity grows. In all cases the binary coalesces once it has absorbed the angular momentum of a gas mass comparable to that of the secondary black hole. Importantly, this conclusion is unaffected even if the gas inflow rate through the disc is formally super--Eddington for either hole. The coalescence timescale is therefore always $\sim M_2/\dot M$, where $M_2$ is the secondary black hole mass and $\dot M$ the inflow rate through the circumbinary disc.

The stellar origin of gamma-ray bursts can be explained by the rapid release of energy in a highly collimated, extremely relativistic jet. This in turn appears to require a rapidly spinning highly magnetised stellar core that collapses into a magnetic neutron star or a black hole within a relatively massive envelope. They appear to be associated with type Ib/c supernovae but, with a birthrate of around 10^{-6}-10^{-5} per year per galaxy, they are considerably rarer than such supernovae in general. To satisfy all these requirements we hypothesize a binary star model that ends with the merging of an oxygen neon white dwarf with the carbon-oxygen core of a naked helium star during a common envelope phase of evolution. The rapid spin and high magnetic field are natural consequences of such a merging. The evolution that leads to these progenitors is convoluted and so naturally occurs only very rarely. To test the hypothesis we evolve a population of progenitors and find that the rate is as required. At low metallicity we calculate that a similar fraction of stars evolve to this point and so would expect the gamma-ray burst rate to correlate with the star formation rate in any galaxy. This too is consistent with observations. These progenitors, being of intermediate mass, differ radically from the usually postulated high-mass stars. Thus we can reconcile observations that the bursts occur close to but not within massive star associations.

We consider models of gas flow in spiral galaxies in which the spiral structure has been excited by various possible mechanisms: a global steady density wave, self-gravity of the stellar disc and an external tidal interaction, as well as the case of a galaxy with a central rotating bar. In each model we estimate in a simple manner the likely current positions of star clusters of a variety of ages, ranging from ~ 2 Myr to around 130 Myr, depending on the model. We find that the spatial distribution of cluster of different ages varies markedly depending on the model, and propose that observations of the locations of age-dated stellar clusters is a possible discriminant between excitation mechanisms for spiral structure in an individual galaxy.

The nearest accretion disc to us in space if not time was the protosolar nebula. Remnants of this nebula thus potentially offer unique insight into how discs work. In particular the existence of chondrules, which must have formed in the disc as small molten droplets, requires strong and intermittent heating of disc material. We argue that this places important constraints on the way gravitational energy is released in accretion discs, which are not met by current shearing--box simulations of MRI--driven dynamos. A deeper understanding of accretion energy release in discs may require a better model for these dynamos.

We consider the response of a protostellar disc to a tidally induced warp and the resultant changes in the spectral energy distribution (SED). We argue that for typical protostellar disc parameters the warp is communicated through the disc in a wave-like fashion. We find that the main effects of the warp tend to be at large radii (greater than 30 AU) and, for sufficiently small viscosity, can be quite long-lived. This can result in non-uniform illumination of the disc at these radii and can induce significant changes to the SED at wavelengths greater than 100 microns.

We present hydrodynamical models of the grand design spiral M51 (NGC 5194), and its interaction with its companion NGC 5195. Despite the simplicity of our models, our simulations capture the present day spiral structure of M51 remarkably well, and even reproduce details such as a kink along one spiral arm, and spiral arm bifurcations. We investigate the offset between the stellar and gaseous spiral arms, and find at most times (including the present day) there is no offset between the stars and gas to within our error bars. We also compare our simulations with recent observational analysis of M51. We compute the pattern speed versus radius, and like the observations, find no single global pattern speed. We also show that the spiral arms cannot be fitted well by logarithmic spirals. We interpret these findings as evidence that M51 does not exhibit a quasi-steady density wave, as would be predicted by density wave theory. The internal structure of M51 derives from the complicated and dynamical interaction with its companion, resulting in spiral arms showing considerable structure in the form of short-lived kinks and bifurcations. Rather than trying to model such galaxies in terms of global spiral modes with fixed pattern speeds, it is more realistic to start from a picture in which the spiral arms, while not being simple material arms, are the result of tidally induced kinematic density `waves' or density patterns, which wind up slowly over time.

The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.

The low-mass X-ray binary microquasar GRO J1655-40 is observed to have a misalignment between the jets and the binary orbital plane. Since the current black hole spin axis is likely to be parallel to the jets, this implies a misalignment between the spin axis of the black hole and the binary orbital plane. It is likely the black holes formed with an asymmetric supernova which caused the orbital axis to misalign with the spin of the stars. We ask whether the null hypothesis that the supernova explosion did not affect the spin axis of the black hole can be ruled out by what can be deduced about the properties of the explosion from the known system parameters. We find that this null hypothesis cannot be disproved but we find that the most likely requirements to form the system include a small natal black hole kick (of a few tens of km/s) and a relatively wide pre-supernova binary. In such cases the observed close binary system could have formed by tidal circularisation without a common envelope phase.

We investigate the evolution of the relative angle between the stellar rotation axis and the circumstellar disc axis of a star that forms in a stellar cluster from the collapse of a turbulent molecular cloud. This is an inherently chaotic environment with variable accretion, both in terms of rate and the angular momentum of the material, and dynamical interactions between stars. We find that the final stellar rotation axis and disc spin axis can be strongly misaligned, but this occurs primarily when the disc is truncated by a dynamical encounter so that the final disc rotation axis depends simply on what fell in last. This may lead to planetary systems with orbits that are misaligned with the stellar rotation axis, but only if the final disc contains enough mass to form planets. We also investigate the time variability of the inner disc spin axis, which is likely to determine the direction of a protostellar jet. We find that the jet direction varies more strongly for lighter discs, such as those that have been truncated by dynamical interactions or have suffered a period of rapid accretion. Finally, we note that variability of the angular momentum of the material accreting by a star implies that the internal velocity field of such stars may be more complicated than that of aligned differential rotation.

Gas at intermediate temperature between the hot X-ray emitting coronal gas in galaxies at the centers of galaxy clusters, and the much cooler optical line emitting filaments, yields information on transport processes and plausible scenarios for the relationship between X-ray cool cores and other galactic phenomena such as mergers or the onset of an active galactic nucleus. Hitherto, detection of intermediate temperature gas has proven elusive. Here, we present FUV imaging of the "low excitation" emission filaments of M87 and show strong evidence for the presence of CIV 1549 A emission which arises in gas at temperature ~10^5K co-located with Halpha+[NII] emission from cooler ~10^4K gas. We infer that the hot and cool phases are in thermal communication, and show that quantitatively the emission strength is consistent with thermal conduction, which in turn may account for many of the observed characteristics of cool core galaxy clusters.

The M51 galaxy is one of best examples of grand design structure in galaxies. The most likely origin of the spiral arms in M51 is the ongoing interaction with NGC 5195. Here we report recent calculations which model the orbit of M51 and NGC 5195. These calculations, which for the first time focus on the gas dynamics, are able to reproduce even detailed features in M51, as well as demonstrate the nature of the spiral structure of M51.

We model the overall shape of an accretion disc in a semi-detached binary system in which mass is transfered on to a spinning black hole the spin axis of which is misaligned with the orbital rotation axis. We assume the disc is in a steady state. Its outer regions are subject to differential precession caused by tidal torques of the companion star. These tend to align the outer parts of the disc with the orbital plane. Its inner regions are subject to differential precession caused by the Lense-Thirring effect. These tend to align the inner parts of the disc with the spin of the black hole. We give full numerical solutions for the shape of the disc for some particular disc parameters. We then show how an analytic approximation to these solutions can be obtained for the case when the disc surface density varies as a power law with radius. These analytic solutions for the shape of the disc are reasonably accurate even for large misalignments and can be simply applied for general disc parameters. They are particularly useful when the numerical solutions would be slow.

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre--existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework. The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high--photon--energy variability on the relatively long orbital timescale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary. We also find that the disc can drive the binary to merger on a reasonable timescale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e. $a_0 \lesssim 0.05$ pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is $\la 0.001$. Gas discs thus appear not to be an effective solution to the `last parsec' problem for a significant class of mergers.

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are thought to be caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion, $\sigma_k \approx 15\rm km s^{-1}$, is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically $i \ge 25^\circ$. Alternatively a higher velocity kick distribution, $\sigma_k = 265 \rm km s^{-1}$, is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularisation of the binary orbits has taken place. We discuss briefly how this might be achieved.

Sokoloski et al (2008) have recently reported evidence that the recurrent nova RS Ophiuchi produced a pair of highly collimated radio jets within days of its 2006 outburst. This suggests that an accretion disc must be present during the outburst. However in the standard picture of recurrent novae as thermonuclear events, any such disc must be expelled from the white dwarf vicinity, as the nuclear energy yield greatly exceeds its binding energy. We suggest instead that the outbursts of RS Oph are thermal--viscous instabilities in a disc irradiated by the central accreting white dwarf. The distinctive feature of RS Oph is the very large size of its accretion disc. Given this, it fits naturally into a consistent picture of systems with unstable accretion discs. This picture explains the presence and speed of the jets, the brightness and duration of the outburst, and its rise time and linear decay, as well as the faintness of the quiescence. By contrast, the hitherto standard picture of recurrent thermonuclear explosions has a number of severe difficulties. These include the presence of jets, the faintness of quiescence, and the fact the the accretion disc must be unstable unless it is far smaller than any reasonable estimate.

We investigate the relationship between the star formation rate per unit area and the surface density of the ISM (the local Kennicutt-Schmitt law) using a simplified model of the ISM and a simple estimate of the star formation rate based on the mass of gas in bound clumps, the local dynamical timescales of the clumps, and an efficiency parameter of around 5 per cent. Despite the simplicity of the approach, we are able to reproduce the observed linear relation between star formation rate and surface density of dense (molecular) gas. We use a simple model for the dependence of H_2 fraction on total surface density to argue why neither total surface density nor the HI surface density are good local indicators of star formation rate. We also investigate the dependence of the star formation rate on the depth of the spiral potential. Our model indicates that the mean star formation rate does not depend significantly on the strength of the spiral potential, but that a stronger spiral potential, for a given mean surface density, does result in more of the star formation occurring close to the spiral arms. This agrees with the observation that grand design galaxies do not appear to show a larger degree of star formation compared to their flocculent counterparts.

In this paper we revisit the arguments for the basis of the time evolution of the flares expected to arise when a star is disrupted by a supermassive black hole. We present a simple analytic model relating the lightcurve to the internal density structure of the star. We thus show that the standard lightcurve proportional to $t^{-5/3}$ only holds at late times. Close to the peak luminosity the lightcurve is shallower, deviating more strongly from $t^{-5/3}$ for more centrally concentrated (e.g. solar--type) stars. We test our model numerically by simulating the tidal disruption of several stellar models, described by simple polytropic spheres with index $\gamma$. The simulations agree with the analytical model given two considerations. First, the stars are somewhat inflated on reaching pericentre because of the effective reduction of gravity in the tidal field of the black hole. This is well described by a homologous expansion by a factor which becomes smaller as the polytropic index becomes larger. Second, for large polytropic indices wings appear in the tails of the energy distribution, indicating that some material is pushed further away from parabolic orbits by shocks in the tidal tails. In all our simulations, the $t^{-5/3}$ lightcurve is achieved only at late stages. In particular we predict that for solar type stars, this happens only after the luminosity has dropped by at least two magnitudes from the peak. We discuss our results in the light of recent observations of flares in otherwise quiescent galaxies and note the dependence of these results on further parameters, such as the star/hole mass ratio and the stellar orbit.

In the microquasar V4641 Sgr the spin of the black hole is thought to be misaligned with the binary orbital axis. The accretion disc aligns with the black hole spin by the Lense-Thirring effect near to the black hole and further out becomes aligned with the binary orbital axis. The inclination of the radio jets and the Fe$K\alpha$ line profile have both been used to determine the inclination of the inner accretion disc but the measurements are inconsistent. Using a steady state analytical warped disc model for V4641 Sgr we find that the inner disc region is flat and aligned with the black hole up to about $900 R_{\rm g}$. Thus if both the radio jet and fluorescent emission originates in the same inner region then the measurements of the inner disc inclination should be the same.

White dwarfs with surface magnetic fields in excess of $1 $MG are found as isolated single stars and relatively more often in magnetic cataclysmic variables. Some 1,253 white dwarfs with a detached low-mass main-sequence companion are identified in the Sloan Digital Sky Survey but none of these is observed to show evidence for Zeeman splitting of hydrogen lines associated with a magnetic field in excess of 1MG. If such high magnetic fields on white dwarfs result from the isolated evolution of a single star then there should be the same fraction of high field white dwarfs among this SDSS binary sample as among single stars. Thus we deduce that the origin of such high magnetic fields must be intimately tied to the formation of cataclysmic variables. CVs emerge from common envelope evolution as very close but detached binary stars that are then brought together by magnetic braking or gravitational radiation. We propose that the smaller the orbital separation at the end of the common envelope phase, the stronger the magnetic field. The magnetic cataclysmic variables originate from those common envelope systems that almost merge. We propose further that those common envelope systems that do merge are the progenitors of the single high field white dwarfs. Thus all highly magnetic white dwarfs, be they single stars or the components of MCVs, have a binary origin. This hypothesis also accounts for the relative dearth of single white dwarfs with fields of 10,000 - 1,000,000G. Such intermediate-field white dwarfs are found preferentially in cataclysmic variables. In addition the bias towards higher masses for highly magnetic white dwarfs is expected if a fraction of these form when two degenerate cores merge in a common envelope. Similar scenarios may account for very high field neutron stars.

The microquasar GRO J1655-40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady state disc warped by Lense-Thirring precession and find the timescale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically fifty times higher in the latter case but we find that, in both cases, the lifetime of the mass transfer state is at most a few times the alignment timescale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter-aligned after its supernova kick but that it is most likely to be close to alignment rather than counteralignment now.

We investigate the evolution of accretion luminosity $L_{\rm acc}$ and stellar luminosity ${L_\ast}$ in pre-mainsequence stars. We make the assumption that when the star appears as a Class II object, the major phase of accretion is long past, and the accretion disc has entered its asymptotic phase. We use an approximate stellar evolution scheme for accreting pre-mainsequence stars based on Hartmann, Cassen & Kenyon, 1997. We show that the observed range of values $k = L_{\rm acc}/L_\ast$ between 0.01 and 1 can be reproduced if the values of the disc mass fraction $M_{\rm disc}/M_*$ at the start of the T Tauri phase lie in the range 0.01 -- 0.2, independent of stellar mass. We also show that the observed upper bound of $L_{\rm acc} \sim L_\ast$ is a generic feature of such disc accretion. We conclude that as long as the data uniformly fills the region between this upper bound and observational detection thresholds, then the degeneracies between age, mass and accretion history severely limit the use of this data for constraining possible scalings between disc properties and stellar mass.

We argue that supermassive black hole growth in AGN occurs via sequences of randomly--oriented accretion discs with angular momentum limited by self--gravity. These stably co-- or counter--align with the black hole spin with almost equal frequency. Accretion from these discs very rapidly adjusts the hole's spin parameter to average values $\bar a \sim 0.1-0.3$ (the precise range depending slightly on the disc vertical viscosity coefficient $\alpha_2$) from any initial conditions, but with significant fluctuations ($\Delta a\sim \pm 0.2$) about these. We conclude (a) AGN black holes should on average spin moderately, with the mean value $\bar a$ decreasing slowly as the mass increases; (b) SMBH coalescences leave little long--term effect on $\bar a$; (c) SMBH coalescence products in general have modest recoil velocities, so that there is little likelihood of their being ejected from the host galaxy; (d) black holes can grow even from stellar masses to $\sim 5\times 10^9 \msun$ at high redshift $z\sim 6$; (e) jets produced in successive accretion episodes can have similar directions, but after several episodes the jet direction deviates significantly. Rare examples of massive holes with larger spin parameters could result from prograde coalescences with SMBH of similar mass, and are most likely to be found in giant ellipticals. We compare these results with observation. (abridged)

We consider the shape of an accretion disc whose outer regions are misaligned with the spin axis of a central black hole and calculate the steady state form of the warped disc in the case where the viscosity and surface densities are power laws in the distance from the central black hole. We discuss the shape of the resulting disc in both the frame of the black hole and that of the outer disc. We note that some parts of the disc and also any companion star maybe shadowed from the central regions by the warp. We compute the torque on the black hole caused by the Lense-Thirring precession and hence compute the alignment and precession timescales. We generalise the case with viscosity and hence surface density independent of radius to more realistic density distributions for which the surface density is a decreasing function of radius. We find that the alignment timescale does not change greatly but the precession timescale is more sensitive. We also determine the effect on this timescale if we truncate the disc. For a given truncation radius, the the timescales are less affected for more sharply falling density distributions.

We observed the brightest central galaxy (BCG) in the nearby (z=0.0821) cool core galaxy cluster Abell 2597 with the IRAC and MIPS instruments on board the Spitzer Space Telescope. The BCG was clearly detected in all Spitzer bandpasses, including the 70 and 160 micron wavebands. We report aperture photometry of the BCG. The spectral energy distribution exhibits a clear excess in the FIR over a Rayleigh-Jeans stellar tail, indicating a star formation rate of ~4-5 solar masses per year, consistent with the estimates from the UV and its H-alpha luminosity. This large FIR luminosity is consistent with that of a starburst or a Luminous Infrared Galaxy (LIRG), but together with a very massive and old population of stars that dominate the energy output of the galaxy. If the dust is at one temperature, the ratio of 70 to 160 micron fluxes indicate that the dust emitting mid-IR in this source is somewhat hotter than the dust emitting mid-IR in two BCGs at higher-redshift (z~0.2-0.3) and higher FIR luminosities observed earlier by Spitzer, in clusters Abell 1835 and Zwicky 3146.

In this paper we explore numerically the evolution of a warped accretion disc. Here, we focus here on the regime where the warp evolves diffusively. By comparing the numerical results to a simple diffusion model, we are able to determine the diffusion coefficient of the warp, $\alpha_2$, as a function of the relevant disc parameters. We find that while in general the disc behaviour is well reproduced by the diffusion model and for relatively large viscosities the warp diffusion is well described by the linear theory (in particular confirming that the warp diffusion coefficient is inversely proportional to viscosity), significant non-linear effects are present as the viscosity becomes smaller, but still dominates over wave-propagation effects. In particular, we find that the inverse dependence of the diffusion coefficient on viscosity breaks down at low viscosities, so that $\alpha_2$ never becomes larger than a saturation value $\alpha_{\rm max}$ of order unity. This can have major consequences in the evolution of systems where a warped disc is present. In particular, it affects the location of the warp radius in the Bardeen-Petterson effect and therefore the spin up (or spin down) of supermassive black holes in the nuclei of galaxies. Additionally, we also find that while the rate of warp diffusion does not depend significantly on the detailed viscosity formulation, the rate of internal precession generated by the warp is strongly affected by it. Such effects should be considered with care when modeling the evolution of warped discs. This emphasises the need to test the above results using different numerical schemes, and with higher resolution, in order to investigate the degree to which numerical simulations are able to provide accurate modeling of the complex fluid dynamics of warped discs. (Abridged)

There is evidence for the existence of massive planets at orbital radii of several hundred AU from their parent stars where the timescale for planet formation by core accretion is longer than the disc lifetime. These planets could have formed close to their star and then migrated outwards. We consider how the transfer of angular momentum by viscous disc interactions from a massive inner planet could cause significant outward migration of a smaller outer planet. We find that it is in principle possible for planets to migrate to large radii. We note, however, a number of effects which may render the process somewhat problematic.

We consider observational and theoretical estimates of the accretion disc viscosity parameter $\alpha$. We find that in thin, fully-ionized discs, the best observational evidence suggests a typical range $\alpha \sim 0.1 - 0.4$, whereas the relevant numerical simulations tend to derive estimates for $\alpha$ which are an order of magnitude smaller. We discuss possible reasons for this apparent discrepancy.

We suggest that most nearby active galactic nuclei are fed by a series of small--scale, randomly--oriented accretion events. Outside a certain radius these events promote rapid star formation, while within it they fuel the supermassive black hole. We show that the events have a characteristic time evolution. This picture agrees with several observational facts. The expected luminosity function is broadly in agreement with that observed for moderate--mass black holes. The spin of the black hole is low, and aligns with the inner disc in each individual feeding event. This implies radio jets aligned with the axis of the obscuring torus, and uncorrelated with the large--scale structure of the host galaxy. The ring of young stars observed about the Galactic Centre are close to where our picture predicts that star formation should occur.

We present an extension of the King et al. (2004) model for the flickering of black hole accretion discs by taking proper account for the thermal properties of the disc. First we develop a one-dimensional, vertically averaged, one-zone model for an optically thick accretion disc and study the temporal evolution. This limits the current model to the so-called high/soft state, where the X-Ray spectrum is dominated by a thermal black-body component. Then we couple this disc model to the flickering process as described in King et al. (2004). Thus we consider the evolution of a poloidal magnetic field subject to a magnetic dynamo. By comparing to observations of X-Ray binaries in the high-soft state, we can constrain the strength of the energy density of the poloidal magnetic field to a few percent of the energy density of the intrinsic disc magnetic field.

We present time-dependent simulations of a two-phase accretion flow around a black hole. The accretion flow initially is composed of an optically thick and cool disc close to the midplane, while on top and below the disc there is a hot and optically thin corona. We consider several interaction mechanisms as heating of the disc by the corona and Compton cooling of the corona by the soft photons of the disc. Mass and energy can be exchanged between the disc and the corona due to thermal conduction. For the course of this more exploratory work, we limit ourselves to one particular model for a stellar mass black hole accreting at a low accretion rate. We confirm earlier both theoretical and observational results which show that at low accretion rates the disc close to the black hole cannot survive and is evaporated. Given the framework of this model, we now can follow through this phase of disc evaporation time dependently.

Simulations and analytic arguments suggest that the turbulence driven by magnetorotational instability (MRI) in accretion discs can amplify the toroidal (azimuthal) component of the magnetic field to a point at which magnetic pressure exceeds the combined gas + radiation pressure in the disc. Arguing from the recent analysis by Pessah and Psaltis, and other MRI results in the literature, we conjecture that the limiting field strength for a thin disc is such that the Alfven speed roughly equals the geometric mean of the Keplerian speed and the gas sound speed. We examine the properties of such magnetically-dominated discs, and show that they resolve a number of outstanding problems in accretion disc theory. The discs would be thicker than standard (Shakura-Sunyaev) discs at the same radius and accretion rate, and would tend to have higher colour temperatures. If they transport angular momentum according to an alpha-prescription, they would be stable against the thermal and viscous instabilities that are found in standard disc models. In discs fuelling active galactic nuclei, magnetic pressure support could also alleviate the restriction on accretion rate imposed by disc self-gravity.

We consider the problem of growing the largest supermassive black holes from stellar--mass seeds at high redshift. Rapid growth without violating the Eddington limit requires that most mass is gained while the hole has a low spin and thus a low radiative accretion efficiency. If, as was formerly thought, the black--hole spin aligns very rapidly with the accretion flow, even a randomly--oriented sequence of accretion events would all spin up the hole and prevent rapid mass growth. However, using a recent result that the Bardeen--Petterson effect causes {\it counter}alignment of hole and disc spins under certain conditions, we show that holes can grow rapidly in mass if they acquire most of it in a sequence of randomly oriented accretion episodes whose angular momenta $J_d$ are no larger than the hole's angular momentum $J_h$. Ultimately the hole has total angular momentum comparable with the last accretion episode. This points to a picture in which the accretion is chaotic on a lengthscale of order the disc size, that is $\la 0.1$ pc.

We present Smoothed Particle Hydrodynamics (SPH) simulations of molecular cloud formation in spiral galaxies. These simulations model the response of a non-self-gravitating gaseous disk to a galactic potential. The spiral shock induces high densities in the gas, and considerable structure in the spiral arms, which we identify as molecular clouds. We regard the formation of these structures as due to the dynamics of clumpy shocks, which perturb the flow of gas through the spiral arms. In addition, the spiral shocks induce a large velocity dispersion in the spiral arms, comparable with the magnitude of the velocity dispersion observed in molecular clouds. We estimate the formation of molecular hydrogen, by post-processing our results and assuming the gas is isothermal. Provided the gas is cold ($T\le100$ K), the gas is compressed sufficiently in the spiral shock for molecular hydrogen formation to occur in the dense spiral arm clumps. These molecular clouds are largely confined to the spiral arms, since most molecular gas is photodissociated to atomic hydrogen upon leaving the arms.

An unsolved issue in the standard core accretion model for gaseous planet formation is how kilometre-sized planetesimals form from, initially, micron-sized dust grains. Solid growth beyond metre sizes can be difficult both because the sticking efficiency becomes very small, and because these particles should rapidly migrate into the central star. We consider here how metre-sized particles evolve in self-gravitating accretion discs using simulations in which the gravitational influence of the solid particles is also included. Metre-sized particles become strongly concentrated in the spiral structures present in the disc and, if the solid to gas density ratio is sufficiently high, can fragment due to their own self-gravity to form planetesimals directly. This result suggests that planetesimal formation may occur very early in the star formation process while discs are still massive enough to be self-gravitating. The dependence of this process on the surface density of the solids is also consistent with the observation that extrasolar planets are preferentially found around high metallicity stars.

The periodic precession (162--day) and nodding (6.3--day) motions of the jets in SS433 are driven in the outer regions of the disc, whereas the jets themselves, being relativistic, are launched near the black hole at the disc centre. Given that the nutation period is comparable to the dynamical timescale in the outer regions of the disc, it seems unlikely that these periods can be communicated efficiently to the disc centre. Here we propose that the massive outflow observed in SS433 is launched at large radii in the disc, about 1/10 of the outer disc edge, and that it is this outflow which responds to the oscillations of the outer disc and determines the direction of the jets. The massive outflow is launched at large radius because the mass transfer rate is hyper-Eddington. This implies not only that the total luminosity of SS433 exceeds $\le$ by a considerable factor, but also that the radiative output is collimated along the outflow. We thus suggest that SS433 is an ultraluminous X-ray source (ULX) viewed `from the side'. We also suggest that the obscured {\it INTEGRAL} sources may be SS433--like objects, but with slightly lower mass transfer rates.

We examine the recent data and analysis of Natta et al. concerning the accretion rate on to young stars as a function of stellar mass, and conclude that the apparently steep dependence of accretion rate on mass is strongly driven by selection/detection thresholds. We argue that a convincing demonstration of a physical relationship between accretion and stellar mass requires further studies which, as is the case for Natta et al., include information on upper limits, and which quantify the possible incompleteness of the sample, at both low and high accretion rates. We point out that the distribution of detections in the ($M, \dot{M}$)-plane can in principle be used to test conventional accretion disc evolutionary models, and that higher sensitivity observations might be able to test the hypothesis of accelerated disc clearing at late times.

We present a new model for protoplanetary disc evolution. This model combines viscous evolution with photoevaporation of the disc, in a manner similar to Clarke, Gendrin & Sotomayor (2001). However in a companion paper (Alexander, Clarke & Pringle 2006a) we have shown that at late times such models must consider the effect of stellar radiation directly incident on the inner disc edge, and here we model the observational implications of this process. We find that the entire disc is dispersed on a time-scale of order $10^5$yr after a disc lifetime of a few Myr, consistent with observations of T Tauri (TT) stars. We use a simple prescription to model the spectral energy distribution of the evolving disc, and demonstrate that the model is consistent with observational data across a wide range of wavelengths. We note also that the model predicts a short ``inner hole'' phase in the evolution of all TT discs, and make predictions for future observations at mid-infrared and millimetre wavelengths.

In this paper we consider the effect of the direct ionizing stellar radiation field on the evolution of protoplanetary discs subject to photoevaporative winds. We suggest that models which combine viscous evolution with photoevaporation of the disc (e.g. Clarke, Gendrin & Sotomayor 2001) incorrectly neglect the direct field after the inner disc has drained, at late times in the evolution. We construct models of the photoevaporative wind produced by the direct field, first using simple analytic arguments and later using detailed numerical hydrodynamics. We find that the wind produced by the direct field at late times is much larger than has previously been assumed, and we show that the mass-loss rate scales as $R_{in}^{1/2}$ (where $R_{in}$ is the radius of the instantaneous inner disc edge). We suggest that this result has important consequences for theories of disc evolution, and go on to consider the effects of this result on disc evolution in detail in a companion paper (Alexander, Clarke & Pringle 2006b).

In this paper we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in Active Galactic Nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole-disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius $R_{\rm w}$ and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines, and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.

We extend the model of King et al. (2004) for variability in black hole accretion discs, by taking proper account of the thermal properties of the disc. Because the degree of variability in the King et al. (2004) model depends sensitively on the ratio of disc thickness to radius, H/R, it is important to follow the time-dependence of the local disc structure as the variability proceeds. In common with previous authors, we develop a one-zone model for the local disc structure. We agree that radial heat advection plays an important role in determining the inner disc structure, and also find limit-cycle behaviour. When the stochastic magnetic dynamo model of King et al. (2004) is added to these models, we find similar variability behaviour to before. We are now better placed to put physical constraints on model parameters. In particular, we find that in order to be consistent with the low degree of variability seen in the thermal disc component of black hole binaries, we need to limit the energy density of the poloidal field that can be produced by local dynamo cells in the disc to less than a few percent of the energy density of the dynamo field within the disc itself.

We present numerical simulations of the passage of clumpy gas through a galactic spiral shock and the subsequent formation of giant molecular clouds (GMCs). The spiral shock forms dense clouds while dissipating kinetic energy, producing regions that are locally gravitationally bound and collapse to form stars. The effect of the clumpiness of gas as it passes through the shock is to generate chaotic internal motions in the gas. The kinematics of these motions are found to agree with the observed velocity-dispersion/size relation found in star-forming regions. In contrast to the standard picture where continuously driven turbulence generates the density inhomogeneities in star-forming clouds, we find here that it is the clumpiness of the interstellar gas that produces the chaotic motions as it passes through the spiral shock and initiates the star formation process. The velocity dispersion can be understood as being due to the random mass loading of clumps as they converge in the spiral shock. In this model there is no need for any internal or external continuous driving mechanism for the 'turbulence'. The coupling of the clouds' internal kinematics to their externally triggered formation removes the need for the clouds to be self-gravitating. Indeed, while clearly some parts of the clouds are self-gravitating and able to form stars, most of the molecular material remains gravitationally unbound. This can provide a simple explanation for the low efficiency of star formation.

We predict the existence of dusty disks and possibly CO planets around massive white dwarfs. We show that the thermal emission from these disks should be detectable in the infrared. The planets may also be detectable either by direct IR imaging, spectroscopy, or using the pulsations of the white dwarfs.

We consider the alignment torque between a spinning black hole and an accretion disc whose angular momenta are misaligned. This situation must hold initially in almost all gas accretion events on to supermassive black holes, and may occur in binaries where the black hole receives a natal supernova kick. We show that the torque always acts to align the hole's spin with the total angular momentum without changing its magnitude. The torque acts dissipatively on the disc, reducing its angular momentum, and aligning it with the hole if and only if the angle theta between the angular momenta J_d of the disc and J_h of the hole satisfies the inequality cos theta > -J_d / 2 J_h. If this condition fails, which requires both theta > pi/2 and J_d < 2 J_h, the disc counteraligns.

We calculate the radial diffusion coefficient for a passive contaminant in an accretion disc which is turbulent due to the action of the magnetorotational instability. Numerical MHD simulations are used to follow the evolution of a local patch of the disc using the shearing box formalism. A separate continuity equation for the mass fraction of contaminant is integrated along with the MHD system, and radial profiles of this fraction are obtained as a function of time. Solutions of a linear diffusion equation are fitted to the numerical measured profiles of the contaminant, treating the diffusion coefficient D as the fitting parameter. At early times, the value of D is found to vary, however once the contaminant is spread over scales comparable to the box size, it saturates at a steady value. The ratio of D to the transport coefficient of angular momentum due to shear stress is small. If D can be used as a proxy for the turbulent magnetic diffusivity, the effective magnetic Prandtl number P_eff=\nu/D (where \nu is the coefficient of "efective viscosity" associated with shear stress) would be large.

We present the results of an analysis of ultraviolet observations of T Tauri Stars (TTS). By analysing emission measures taken from the literature we derive rates of ionizing photons from the chromospheres of 5 classical TTS in the range ~10^41-10^44 photons/s, although these values are subject to large uncertainties. We propose that the HeII/CIV line ratio can be used as a reddening-independent indicator of the hardness of the ultraviolet spectrum emitted by TTS. By studying this line ratio in a much larger sample of objects we find evidence for an ionizing flux which does not decrease, and may even increase, as TTS evolve. This implies that a significant fraction of the ionizing flux from TTS is not powered by the accretion of disc material onto the central object, and we discuss the significance of this result and its implications for models of disc evolution. The presence of a significant ionizing flux in the later stages of circumstellar disc evolution provides an important new constraint on disc photoevaporation models.

We consider the formation and evolution of the planetary system PSR B1620-26 in the globular cluster M4. We propose that as M4 is a very-low metallicity environment the standard model of planet formation around main-sequence stars through the accretion of gas onto metallic rocky cores should not be applied. Consequently the previously suggested methods for formation are unlikely. We propose that the planet formed through the interaction of a passing star with a circumbinary disc during the common-envelope phase of the inner binary's evolution. This formation route is favoured by dense stellar systems such as globular clusters.

In this paper we consider the evolution of small planetesimals in marginally stable, self-gravitating protoplanetary discs. The drag force between the disc gas and the embedded planetesimals generally causes the planetesimals to drift inwards through the disc at a rate that depends on the particle size. In a marginally stable, self-gravitating disc, however, the planetesimals are significantly influenced by the non-axisymmetric spiral structures resulting from the growth of the gravitational instability. For small particles, that are strongly coupled to the disc gas, and for large particles, that have essentially decoupled from the disc gas, the effect is not particularly significant. Intermediate sized particles, which would generally have the largest radial drift rates, do, however, become significantly concentrated at the peaks of the spiral arms. These high density regions may persist for, of order, an orbital period and may attain densities comparable to that of the disc gas. Although at the end of the simulation only ~ 25 % of the planetesimal particles lie in regions of enhanced density, during the course of the simulation at least 75 % of the planetesimal particles have at some stage been in a such a region. We find that the concentration of particles in the spiral arms results in an increased collision rate, an effect that could significantly accelerate planetesimal growth. The density enhancements may also be sufficient for the growth of planetesimals through direct gravitational collapse. The interaction between small planetesimals and self-gravitating spiral structures may therefore play an important role in the formation of large planetesimals that will ultimately coagulate to form terrestrial planets or the cores of gas/ice giant planets.

Most mechanisms proposed for the formation of planets are modified versions of the mechanism proposed for the solar system. Here we argue that, in terms of those planetary systems which have been observed, the case for the solar system being a typical planetary system has yet to be established. We consider the possibility that most observed planetary systems have been formed in some quite different way. If so, it may be that none of the observed planetary systems is likely to harbour an earth-like planet.

We present the results of a theoretical study investigating the effects of photoionization and heating by X-rays on discs around low-mass stars. In particular we address the question of whether or not X-rays can drive a disc wind. First, we construct a 1-dimensional ``quasi-hydrostatic'' model, which solves for the vertical structure introduced by X-ray heating. We consider uniform X-ray illumination of the disc, but the X-ray fluxes required to heat the disc significantly are much greater than those seen by recent observations. When the model is extended to consider heating from a central X-ray source we find that the 1-dimensional model is only valid very close to the star. We extend our analysis to consider a simple 2-dimensional model, treating the disc as a two-layered structure and solving for its density profile self-consistently. For T Tauri stars we are able to set a crude upper limit on the mass-loss rate that can be driven by X-ray photoevaporation, with a value of ~10^-13 g/cm^2/s. Our model is designed to maximise this value, and most likely over-estimates it significantly. However we still find a mass-loss rate which is less than that found in studies of ultraviolet photoevaporation. We conclude that in the presence of a significant UV field, X-ray driven disc winds are unlikely to play a significant role in the evolution of discs around low-mass stars.

We show how the viscous evolution of Keplerian accretion discs can be understood in terms of simple kinetic theory. Although standard physics texts give a simple derivation of momentum transfer in a linear shear flow using kinetic theory, many authors, as detailed by Hayashi & Matsuda 2001, have had difficulties applying the same considerations to a circular shear flow. We show here how this may be done, and note that the essential ingredients are to take proper account of, first, isotropy locally in the frame of the fluid and, second, the geometry of the mean flow.

We address the issue of the production of Lyman continuum photons by T Tauri stars, in an attempt to provide constraints on theoretical models of disc photoionisation. By treating the accretion shock as a hotspot on the stellar surface we show that Lyman continuum photons are produced at a rate approximately three orders of magnitude lower than that produced by a corresponding black body, and that a strong Lyman continuum is only emitted for high mass accretion rates. When our models are extended to include a column of material accreting on to the hotspot we find that the accretion column is extremely optically thick to Lyman continuum photons. Further, we find that radiative recombination of hydrogen atoms within the column is not an efficient means of producing photons with energies greater than 13.6eV, and find that an accretion column of any conceivable height suppresses the emission of Lyman continuum photons to a level below or comparable to that expected from the stellar photosphere. The photospheric Lyman continuum is itself much too weak to affect disc evolution significantly, and we find that the Lyman continuum emitted by an accretion shock is similarly unable to influence disc evolution significantly. This result has important consequences for models which use photoionisation as a mechanism to drive the dispersal of circumstellar discs, essentially proving that an additional source of Lyman continuum photons must exist if disc photoionisation is to be significant.

Observations of accreting systems often show significant variability (10-20 percent of accretion luminosity) on timescales much longer than expected for the disc regions releasing most of the luminosity. We propose an explicit physical model for disc variability, consistent with Lyubarskii's (1997) general scheme for solving this problem. We suggest that local dynamo processes can affect the evolution of an accretion disc by driving angular momentum loss in the form of an outflow (a wind or jet). We model the dynamo as a small-scale stochastic phenomenon, operating on roughly the local dynamical timescale. We argue that large-scale outflow can only occur when the small-scale random processes in neighbouring disc annuli give rise by chance to a coherent large-scale magnetic field. This occurs on much longer timescales, and causes a bright large-amplitude flare as a wide range of disc radii evolve in a coherent fashion. Most of the time, dynamo action instead produces small-amplitude flickering. We reproduce power spectra similar to those observed, including a 1/f power spectrum below a break frequency given by the magnetic alignment timescale at the inner disc edge. However the relation between the black hole mass and the value of the break frequency is less straightforward than often assumed in the literature. The effect of an outer disc edge is to flatten the spectrum below the magnetic alignment frequency there. We also find a correlation between the variability amplitude and luminosity, similar to that found in some AGN.

Recent observations show a clear correlation between the probability of hosting a planet and the metallicity of the parent star. Since radial velocity surveys are biased, however, towards detecting planets with short orbital periods, the probability-metallicity correlation could merely reflect a dependence of migration rates on metallicity. We investigated the possibility, but find no basis to suggest that the migration process is sensitive to the metallicity. The indication is, therefore, that a higher metallicity results in a higher probability for planet formation.

We present the results of a Hubble Space Telescope (HST) survey of extended [OIII] emission in a sample of 60 nearby Seyfert galaxies (22 Seyfert 1's and 38 Seyfert 2's), selected by mostly isotropic properties. The comparison between the semi major axis size of their [OIII] emitting regions (R_Maj) shows that Seyfert 1's and Seyfert 2's have similar distributions, which seems to contradict Unified Model predictions. We discuss possible ways to explain this result, which could be due either to observational limitations or the models used for the comparison with our data. We show that Seyfert 1 Narrow Line Regions (NLR's) are more circular and concentrated than Seyfert 2's, which can be attributed to foreshortening in the former. We find a good correlation between the NLR size and luminosity, following the relation R_Maj propto L([OIII])^0.33, which is flatter than a previous one found for QSO's and Seyfert 2's. We discuss possible reasons for the different results, and their implications to photoionization models. We confirm previous results which show that the [OIII] and radio emission are well aligned, and also find no correlation between the orientation of the extended [OIII] emission and the host galaxy major axis. This agrees with results showing that the torus axis and radio jet are not aligned with the host galaxy rotation axis, indicating that the orientation of the gas in the torus, and not the spin of the black hole, determine the orientation of the accretion disk, and consequently the orientation of the radio jet.

We present an analytical method for constraining the substellar initial mass function in globular clusters, based on the observed frequency of transit events. Globular clusters typically have very high stellar densities where close encounters are relatively common, and thus tidal capture can occur to form close binary systems. Encounters between main sequence stars and lower-mass objects can result in tidal capture if the mass ratio is > 0.01. If brown dwarfs exist in significant numbers, they too will be found in close binaries, and some fraction of their number should be revealed as they transit their stellar companions. We calculate the rate of tidal capture of brown dwarfs in both segregated and unsegregated clusters, and find that the tidal capture is more likely to occur over an initial relaxation time before equipartition occurs. The lack of any such transits in recent HST monitoring of 47 Tuc implies an upper limit on the frequency of brown dwarfs (< 15 % relative to stars) which is significantly below that measured in the galactic field and young clusters.

We have investigated the dynamics of the merging process in the minor merger hypothesis for active galactic nuclei. We find that for a satellite galaxy to be able to merge directly with the nucleus of the host galaxy (for example, to give rise to the compact dust discs which are seen in early type active galaxies) requires the initial orbit of the satellite to be well aimed. For the case of the host galaxy being a disc galaxy, if the initial orbits of the satellites are randomly oriented with respect to the host galaxy, then the orbits of those which reach the host nuclear regions in a reasonable time, are also fairly randomly oriented once they reach the nucleus. We note that this result might be able to provide an explanation of why the jet directions in the nuclei of Seyfert galaxies are apparently unrelated to the plane of the galaxy discs.

We propose a new model for the disk-jet connection in black hole x-ray transients and active galactic nuclei. In our model, the inner part of the accretion disk around the central black hole switches between two states. In one state, the accretion energy is dissipated locally (some within the disk, some within a corona) to produce the observed disk luminosity. In the second state, the accretion energy is deposited into the bulk flow of a relativistic jet. We associate the transition between the two states with the generation of a global, poloidal magnetic field. We show that this model can explain the observed behavior of black hole accretors.

We examine the hypothesis that there exists a simple scaling between the observed velocities of jets found in Young Stellar Objects (YSOs) and jets found in Active Galactic Nuclei (AGN). We employ a very simplified physical model of the jet acceleration process. We use time-dependent, spherically symmetric wind models in Newtonian and relativistic gravitational fields to ask whether the energy input rates required to produce the jet velocities observed in YSOs (of about 2 $\times$ the escape velocity from the central object) can also produce AGN jet velocities (Lorentz factors of $\gamma \sim$ 10). Such a scaling would be expected if there is a common production mechanism for such jets. We demonstrate that such a scaling does exist, provided that the energy input process takes place sufficiently deep in the gravitational potential well, enabling physical use to be made of the speed of light as a limiting velocity, and provided that the energy released in the accretion process is imparted to a small fraction of the available accreting material.

We study the viscous effects on the interaction between the coplanar Be-star disc and the neutron star in Be/X-ray binaries, using a three-dimensional, smoothed particle hydrodynamics code. For simplicity, we assume the Be disc to be isothermal at the temperature of half the stellar effective temperature. In order to mimic the gas ejection process from the Be star, we inject particles with the Keplerian rotation velocity at a radius just outside the star. Both Be star and neutron star are treated as point masses. We find that the Be-star disc is effectively truncated if the Shakura-Sunyaev viscosity parameter alpha_SS >> 1, which confirms the previous semi-analytical result. In the truncated disc, the material decreted from the Be star accumulates, so that the disc becomes denser more rapidly than if around an isolated Be star. The resonant truncation of the Be disc results in a significant reduction of the amount of gas captured by the neutron star and a strong dependence of the mass capture rate on the orbital phase. We also find that an eccentric mode is excited in the Be disc through direct driving due to a one-armed bar potential of the binary. The strength of the mode becomes greater in the case of a smaller viscosity. In a high-resolution simulation with alpha_SS=0.1, the eccentric mode is found to precess in a prograde sense. The mass capture rate by the neutron star modulates as the mode precesses.

We consider the evolution of a warped disc around a Kerr black hole, under conditions such that the warp propagates in a wavelike manner. This occurs when the dimensionless effective viscosity, alpha, that damps the warp is less than the characteristic angular semi-thickness, H/R, of the disc. We adopt linearized equations that are valid for warps of sufficiently small amplitude in a Newtonian disc, but also account for the apsidal and nodal precession that occur in the Kerr metric. Through analytical and time-dependent studies, we confirm the results of Demianski & Ivanov, and of Ivanov & Illarionov, that such a disc takes on a characteristic warped shape. The inner part of the disc is not necessarily aligned with the equator of the hole, even in the presence of dissipation. We draw attention to the fact that this might have important implications for the directionality of jets emanating from discs around rotating black holes.

In a recent paper (astro-ph/0207561, posted on July 25, 2002), Fromang & Balbus state that ``a uniformly rotating barotropic fluid in an external potential attains a true energy minimum if and only if the rotation profile is everywhere subsonic. If regions of supersonic rotation are present, fluid variations exist that could take the system to states of lower energy.'' We present a simple counter-example which demonstrates that the stablity criterion presented by Fromang & Balbus is incorrect.

We use HST broad band images, VLA and VLBI continuum data to study the three dimensional orientation of jets relative to nuclear dust disks in 20 radio galaxies. The comparison between the position angles of the jets with those of the dust disks major axes shows a wide distribution, suggesting that they are not aligned preferentially perpendicular to each other. We use a statistical technique to determine the 3 dimensional distribution of angles between jets and dust disks rotation axes. This analysis shows that the observations are consistent with jets homogeneously distributed over a large region, extending over polar caps of 55 degrees to 77 degrees, but seem to avoid lying close to the plane of the dust disks. We argue that the lack of close alignment between jets and dust disks axes is not likely to be caused by feeding the nucleus with gas from mergers originated from random directions. We suggest that the misalignment can be due by a warping mechanism in the accretion disk, like self-irradiation instability or the Bardeen-Petterson effect, or that the gravitational potential in the inner regions of the galaxy is misaligned with respect to that of the dust disk.

We investigate the migration of massive extrasolar planets due to gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 AU at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log (a), where a is the orbital radius. Between 1 AU and 3 AU, the expected number of planets per logarithmic interval in orbital radius roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, about 15 percent of stars targeted by existing radial velocity searches are predicted to possess planets with masses 0.3 M_Jupiter < M_p sin (i) < 10 M_Jupiter, and radii 0.1 AU < a < 5 AU. A third of these planets (around 5 percent of the target stars) lie at the radii most amenable to detection via microlensing. A further 5-10 percent of stars could have planets at radii of 5 AU < a < 8 AU that have migrated outwards. We discuss the probability of forming a system (akin to the Solar System) in which significant radial migration of the most massive planet does not occur. About 10-15 percent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low mass planets than high mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.

Using current observational data and simple dynamical modeling, we demonstrate that Omega Cen is not special among the Galactic globular clusters in its ability to produce and retain the heavy elements dispersed in the AGB phase of stellar evolution. The multiple stellar populations observed in Omega Cen cannot be explained if it had formed as an isolated star cluster. The formation within a progenitor galaxy of the Milky Way is more likely, although the unique properties of Omega Cen still remain a mystery.

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic ($\gamma = 5/3$). Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance ${\cal M}_{\rm r} \ga 0.01$ undergo shocks within a distance of order the resonance radius.

We present a series of statistical tests using homogeneous data and measurements for a sample of Seyfert galaxies. These galaxies were selected from mostly isotropic properties, their far infrared fluxes and warm infrared colors, which provide a considerable advantage over the criteria used by most investigators in the past, like ultraviolet excess. Our results provide strong support for a Unified Model in which Seyferts 2's contain a torus seen more edge-on than in Seyferts 1's and show that previous results showing the opposite were most likely due to selection effects.

Withdrawn in current form.

In a recent paper, Elmegreen (2000) has made a cogent case, from an observational point of view, that the lifetimes of molecular clouds are comparable to their dynamical timescales. If so, this has important implications for the mechanisms by which molecular clouds form. In particular we consider the hypothesis that molecular clouds may form not by {\it in situ} cooling of atomic gas, but rather by the agglomeration of the dense phase of the interstellar medium (ISM), much, if not most, of which is already in molecular form.

We present a series of statistical tests done to a sample of 29 Seyfert 1 and 59 Seyfert 2 galaxies selected from mostly isotropic properties, their far infrared fluxes and warm infrared colors. Such selection criteria provide a profound advantage over the criteria used by most investigators in the past, such as ultraviolet excess. These tests were done using ground based high resolution VLA A-configuration 3.6 cm radio and optical B and I imaging data. From the relative number of Seyfert 1's and Seyfert 2's we calculate that the torus half opening angle is 48deg. We show that, as seen in previous papers, there is a lack of edge-on Seyfert 1 galaxies, suggesting dust and gas along the host galaxy disk probably play an important role in hiding some nuclei from direct view. We find that there is no statistically significant difference in the distribution of host galaxy morphological types and radio luminosities of Seyfert 1's and Seyfert 2's, suggesting that previous results showing the opposite may have been due to selection effects. The average extension of the radio emission of Seyfert 1's is smaller than that of Seyfert 2's by a factor of ~2-3, as predicted by the Unified Model. A search for galaxies around our Seyferts allows us to put a lower and an upper limit on the possible number of companions around these galaxies of 19% and 28%, respectively, with no significant difference in the number of companion galaxies between Seyfert 1's and Seyfert 2's. We also show that there is no preference for the radio jets to be aligned closer to the host galaxy disk axis in late type Seyferts, unlike results claimed by previous papers. These results, taken together, provide strong support for a Unified Model in which type 2 Seyferts contain a torus seen more edge-on than the torus in type 1 Seyferts.

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 present a simple physical mechanism that can account for the observed stellar mass spectrum for masses $\ms \simgreat 0.5 \solm$. The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a $\rho\propto R^{-2}$ distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of $\gamma=-3/2$ (where Salpeter is $\gamma=-2.35$). Once the stars dominate the potential and are virialised, which occurs first in the cluster core, the accretion radius is the Bondi-Hoyle radius. The resultant mass spectrum has an asymptotic limit of $\gamma=-2$ with slightly steeper slopes ($\gamma\approx-2.5$) if the stars are already mass-segr egated. Simulations of accretion onto clusters containing 1000 stars show that as expected, the low-mass stars accumulate the majority of their masses during the gas dominated phase whereas the high-mass stars accumulate the majority of their massed during the stellar dominated phase. This results in a mass spectrum with a relatively shallow $\gamma\approx 3/2$ power-law for low-mass stars and a steeper, power-law for high-mass stars $ -2.5\simless\gamma\le -2$. This competitive accretion model also results in a mass segregated cluster.

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity that results means that the tidal radius is smaller than the Bondi-Hoyle radius in these systems. In contrast, when the stars dominate the potential and are virialised, the Bondi-Hoyle radius is smaller than the tidal radius and thus Bondi-Hoyle accretion is a better fit to the accretion rates.

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 AU such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion onto the protostar occurs in bursts with an accretion rate of 10^{-5} solar masses / yr, separated by quiescent intervals where the accretion rate is 10^{-8} solar masses / yr. Such bursts could lead to repeated episodes of strong mass outflow in Young Stellar Objects. The transition to this episodic mode of accretion occurs at an early epoch (t < 1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at around 1 AU. The predicted rate of low mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1-3 AU.

On the basis of the current observational evidence, we put forward the case that the merger of two CO white dwarfs produces both a Type Ia supernova explosion and a stellar remnant, the latter in the form of a magnetar. The estimated occurrence rates raise the possibility that many, if not most, SNe Ia might result from white dwarf mergers.

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow to previous hydrodynamical simulations. The net mass accretion rate is small compared to the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared to the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

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.

We present the study of the relative angle between the accretion disk (or radio jet) and the galaxy disk for a sample of Seyfert galaxies selected from a mostly isotropic property, the 60$\mu$m flux. For each galaxy we have a pair of points ($i$,$\delta$), which are the inclination of the galaxy relative to the line of sight and the angle between the jet projected into the plane of the sky and the host galaxy major axis, respectively. This data is combined with a statistical technique to determine the distribution of $\beta$ angles {\it in 3 dimensions}, the angle between the jet and the host galaxy plane axis. The analysis of the data, not differentiating between Seyfert 1's and 2's, showed that the observed distribution of $i$ and $\delta$ values can be well represented by a homogeneous $\beta-$distribution in the range $0^{\circ}\leq\beta\leq90^{\circ}$. However, when we distinguished between Seyfert 1's and 2's, the models could not represent Seyfert 1's adequately. It was necessary to introduce viewing angle restrictions, that a galaxy can only be recognized as a Seyfert 1 if the angle between the jet and the line of sight ($|\phi|$) is smaller than a given angle $\phi_c$ and that the galaxy inclination $i$ is smaller than an angle $i_c$, in order to have statistically acceptable models. This is a direct and independent confirmation of the underlying concepts of the Unified Model. We discuss several ways to explain the misalignment between the accretion disk axis and the host galaxy disk axis.

We demonstrate that, by combining optical, radio and X-ray observations of a Seyfert, it is possible to provide a direct measurement of the angle $\beta$ between the direction of the radio jet and the normal to the plane of the spiral host galaxy. To do so, we make the assumptions that the inner radio jet is perpendicular to the X-ray observed inner accretion disk, and that the observed jet (or the stronger component, if the jet is two-sided) is physically closer to Earth than the plane of the galaxy. We draw attention to the possibility of measurement producing a result which is not self-consistent, in which case for that galaxy, one of the assumptions must fail.

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows which cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function \nu. We study the flows generated as the amplitude and form of \nu are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by, and a large fraction of the mass participating in, the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. Near the equatorial plane, the radial profiles of the time-averaged variables are power-laws with an index that depends on the radial scaling of the shear stress. We find that regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.

Precessing accretion discs have long been suggested as explanations for the long periods observed in a variety of X-ray binaries, most notably Her X-1/HZ Her. We show that an instability of the disc's response to the radiation reaction force from the illumination by the central source can cause the disc to tilt out of the orbital plane and precess in something like the required manner. The rate of precession and disc tilt obtained for realistic values of system parameters compare favourably with the known body of data on X-ray binaries with long periods. We explore other possible types of behaviour than steadily precessing discs that might be observable in systems with somewhat different parameters. At high luminosities, the inner disc tilts through more than 90 degrees, i.e. it rotates counter to the usual direction, which may explain the torque reversals in systems such as 4U 1626-67.

We critically assess the role of the Blandford-Znajek mechanism in the powering of outflows from accretion disk-fed black holes. We argue that there is no reason to suppose that the magnetic field threading the central spinning black hole differs significantly in strength from that threading the central regions of the disk. In this case, we show that the electromagnetic output from the inner disk regions is expected in general to dominate over that from the hole. Thus the spin (or not) of the hole is probably irrelevant to the expected electromagnetic power output from the system. We also point out that the strength of the poloidal field in the center of a standard accretion disk has been generally overestimated, and discuss scenarios which might give rise to more significant central poloidal fields.

The inner parts of an accretion disk around a spinning black hole are forced to align with the spin of the hole by the Bardeen-Petterson effect. Assuming that any jet produced by such a system is aligned with the angular momentum of either the hole or the inner disk, this can, in principle provide a mechanism for producing steady jets in AGN whose direction is independent of the angular momentum of the accreted material. However, the torque which aligns the inner disk with the hole, also, by Newton's third law, tends to align the spin of the hole with the outer accretion disk. In this letter, we calculate this alignment timescale for a black hole powering an AGN, and show that it is relatively short. This timescale is typically much less than the derived ages for jets in radio loud AGN, and implies that the jet directions are not in general controlled by the spin of the black hole. We speculate that the jet directions are most likely controlled either by the angular momentum of the accreted material or by the gravitational potential of the host galaxy.

Enigmatic transitions between spin-up and spin-down have been observed in several X-ray pulsars accreting matter via an accretion disk. In these transitions, the torque changes sign but remains at nearly the same magnitude. It has been noted previously that alternating prograde and retrograde disk flows would explain many features of the torque reversals, although it has been unclear how a stable retrograde disk could be formed. We suggest that the reversals may be related to the disk at times being warped to such an extent that the inner region becomes tilted by more than 90 degrees. This region would thus become retrograde, leading to a negative torque. Accretion disk models can show such behavior, if account is taken of a warping instability due to irradiation. The resulting `flip-overs' of the inner parts of the disk can reproduce most characteristics of the observations, although it remains unclear what sets the timescale on which the phenomenon occurs. If this model were correct, it would have a number of ramifications, for instance that in the spin-down state the X-ray source would mostly be observed through the accretion disk.

We present a new method for the statistical investigation into the distributions of the angle beta between the radio axis and the normal to the galactic disk for a sample of Seyfert galaxies. We discuss how further observations of the sample galaxies can strengthen the conclusions. Our data are consistent with the hypothesis that AGN jets are oriented randomly in space, independent of the position of the plane of the galaxy. By making the simple assumption that the Standard Model of AGN holds, with a universal opening angle of the thick torus of phi_c, we demonstrate a statistical method to obtain an estimate of phi_c. Our data are not consistent with the simple-minded idea that Seyfert 1s and Seyfert 2s are differentiated solely by whether or not our line of sight lies within some fixed angle of the jet axis. Our result is significant on the 2 sigma level and can thus be considered only suggestive, not conclusive. A complete sample of Seyfert galaxies selected on an isotropic property is required to obtain a conclusive result.

We examine the consequences of radiatively driven warping of accretion disks surrounding pre-main-sequence stars. These disks are stable against warping if the luminosity arises from a steady accretion flow, but are unstable at late times when the intrinsic luminosity of the star overwhelms that provided by the disk. Warps can be excited for stars with luminosities of around 10 solar luminosities or greater, with larger and more severe warps in the more luminous systems. A twisted inner disk may lead to high extinction towards stars often viewed through their disks. After the disk at all radii becomes optically thin, the warp decays gradually on the local viscous timescale, which is likely to be long. We suggest that radiation induced warping may account for the origin of the warped dust disk seen in Beta Pictoris, if the star is only around 10-20 Myr old, and could lead to non-coplanar planetary systems around higher mass stars.

A geometrically thin, optically thick, warped accretion disk with a central source of luminosity is subject to non-axisymmetric forces due to radiation pressure; the resulting torque acts to modify the warp. In a recent paper, \cite{pri96} used a local analysis to show that initially planar accretion disks are unstable to warping driven by radiation torque. Here we extend this work with a global analysis of the stable and unstable modes. We confirm Pringle's conclusion that thin centrally-illuminated accretion disks are generically unstable to warping via this mechanism; we discuss the time-evolution and likely steady-state of such systems and show specifically that this mechanism can explain the warping of the disk of water masers in NGC 4258 and the 164-day precession period of the accretion disk in SS 433. Radiation-driven warping and precession provides a robust mechanism for producing warped, precessing accretion disks in active galactic nuclei and X-ray binary systems.

We explore the consequences of a magnetic dynamo origin for the viscosity in accretion discs, for the structure and evolution of discs in dwarf nova systems. We propose that the rapid cooling that sets in at the end of a dwarf nova eruption acts to inhibit the Balbus-Hawley instability, and thereby to quench dynamo action and so reduce disc viscosity. We demonstrate that a modified disc instability model can reproduce the basic properties of dwarf nova eruptions, as well as some properties of quiescent discs. We also discuss some observational consequences of our model.

We investigate the response of a circumstellar accretion disc to the fly-by of a perturbing mass on a parabolic orbit. The energy and angular momentum transferred during the encounter are calculated using a reduced three-body method. In almost all close encounters the energy and angular momentum transfer is dominated by disc material becoming unbound from the system, with the contributions from close disc particle -- star encounters being significant. For more distant encounters with some prograde element to the motion the disc material loses energy and angular momentum to the perturber's orbit through a resonance feature. The magnitude of the energy transfer calculated in our simulations is greater than that of the binding energy of material exterior to periastron by a factor of two in the prograde case, and up to a factor of five in the case of the retrograde encounter. The destructive nature of the encounters indicates that a non-linear treatment is essential in all but the most distant encounters.

We examine how accretion on to T Tauri stars may be modulated by a time-dependent `magnetic gate' where the inner edge of the accretion disc is disrupted by a varying stellar field. We show that magnetic field variations on time-scales shorter than 10^5 yr can modulate the accretion flow, thus providing a possible mechanism both for the marked photometric variability of T Tauri stars and for the possible conversion of T Tauri stars between classical and weak line status. We thus suggest that archival data relating to the spectrophotometric variability of T Tauri stars may provide an indirect record of magnetic activity cycles in low-mass pre-main-sequence stars.