Search

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., Comment: 21 pages, 17 figure - accepted for publication in MNRAS main journal

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., 17 pages, 17 figures, accepted for publication in MNRAS

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 $, 13 pages, 15 figures, accepted for publication in MNRAS

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 M⊙ yr −1 , 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.

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 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.

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., Comment: 7 pages, 9 figures. Accepted by MNRAS

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.

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., 11 pages, 10 figures, MNRAS accepted

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), Comment: MNRAS, in press

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., Accepted for publication in MNRAS

We study the contamination of the B-mode of the cosmic microwave background polarization (CMBP) by Galactic synchrotron in the lowest emission regions of the sky. The 22.8-GHz polarization map of the 3-yr Wilkinson Microwave Anisotropy Probe (WMAP) data release is used to identify and analyse such regions. Two areas are selected with signal-to-noise ratio S/N 0.3). Such regions allow T/S similar to 10(-3) to be measured directly which approximately corresponds to the limit imposed by using a sky coverage of 15 per cent. This opens interesting perspectives on the investigation of the inflationary model space in lowest emission regions.

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., 13 pages, 16 figures, accepted for publication in MNRAS

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., Comment: 8 pages, 2 figures, MNRAS accepted

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., 6 pages, 3 figures, MNRAS accepted

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., 15 pages, including 20 figures and 1 table, accepted for publication in MNRAS

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., 6 pages, no figures, accepted on 20 June 2001 for publication in MNRAS

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 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., MNRAS, in press

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., Comment: submitted to MNRAS, 17-Dec-97, revised submit 2-Nov-98. 15 pages LaTeX, 11 postscript figures in-text

We consider the thermal and dynamical balance in dense molecular cloud cores. We point out that because such cores are heated throughout their volumes by cosmic rays, but cool radiatively by means of optically thick lines, they are likely to be dynamically unstable and at best in only a form of quasi-static equilibrium. We also draw attention to the possibility that if such cores undergo an external impulse that is large enough to produce a velocity difference ~v;::S c., where Cs is the thermal sound speed, then the cooling rate rises suddenly and pressure support drops, which leads to gravitational collapse and fragmentation.

We consider a thin accretion disc of half-thickness H, vertically threaded by a magnetic field. The field is due to contributions from both the disc current and an external current (giving rise to a uniform external field). We derive an integro-differential equation for the evolution of the magnetic field, subject to magnetic diffusivity η and disc accretion with radial velocity v r . The evolution equation is solved numerically, and a steady state is reached. The evolution equation depends upon a single, dimensionless parameter D=2η/(3H|v r |)=(R/H)(η/ν), where the latter equality holds for a viscous disc having viscosity ν. At the disc surface, field lines are bent by angle i from the vertical, such that tan i=1.52D −1

Ultraviolet spectra of nova-like variables appear to be cooler than those of dwarf novae in outburst. This suggests that, once a steady-state accretion disc becomes established, it releases less of its energy thermally and more in other forms such as mass loss or magnetic flux. Here we use archival IUE data for the anti-dwarf nova TT Ari to show that its spectrum also is cooler at the end of a period in a high state than it is just after a minimum. Models such as that proposed by Tout & Pringle for magnetic viscosity in the disc could account for this kind of behaviour

We consider the linear stability of an isothermal, pressure-bounded, self-gravitating gas slab. Such a configuration is unstable at sufficiently long wavelengths and grows at rate ∼ √Gρ, where ρ is the density of the gas slab, but at high external pressure the nature of this instability is quite different from that of the usual Jeans instability. We develop an analytical model that reproduces the behaviour of the instability found numerically by Elmegreen & Elmegreen. They found, surprisingly, that the critical wavenumber for the onset of the instability is always of the order of the layer thickness, independently of the level of self- gravity and the external pressure

We investigate the effect on an accretion disc around one star of a parabolic fly-by of another. In the prograde, coplanar encounter, the disc is tidally stripped down to about one half of the periastron radius, and a substantial fraction of the stripped material is captured by the intruder. In the retrograde, coplanar encounter, the disc is essentially unaffected within periastron. We also consider an encounter in which the disc and orbit are initially orthogonal. Here the disc outside periastron is only partially tidally stripped, essentially no material is captured by the intruder, and the disc which remains is twisted. We comment briefly on the significance of the above results to models of binary star formation

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., Comment: MNRAS, accepted

A simple set of equations is introduced which governs the time evolution of a twisted accretion disc. The time evolution is governed by two «viscosities», one governing shear within the plane of the disc and the other governing shear perpendicular to the disc (brought about by non-planarity of the disc). It is shown that these equations can be put in numerical difference form so that angular momentum is locally conserved. Numerical simulations are given for two simple cases, the evolution of a simple twist in a steady Keplerian accretion disc, and the effect of precession about an axis which is not aligned with the disc on the inner regions of the disc (the Bardeen-Petterson effect)

A model is constructed for the spin-down of rapidly rotating, convective stars. It is assumed that rotation and convection cause a star to rotate differentially, that differential rotation and convection generate a magnetic dynamo, that a magnetic dynamo gives rise to mass loss and a magnetically controlled stellar wind, and that such a wind results in angular momentum loss and hence stellar spin-down. Using the model, it is shown that a protostar accreting from a circumstellar disc reaches an equilibrium rotation period significantly below break-up

We use the properties of OB runaway stars and, assuming that they are ejected from their parent clusters by dynamical binary-binary interactions, deduce certain necessary characteristics of high-mass star formation. We deduce that massive stars form in small groups of binaries that are biased towards unit mass ratio, and that within these groups the IMF is severely under-represented in low-mass stars compared with the field star IMF

The possibility of a large, massive protostellar accretion disc playing a role in the formation of binary stars by enabling the capture of a passing star within a dense star-forming region is analyzed. It is found that capture rates are too low to play a major role in all known star-forming environments, particularly when the probability of prior disc dispersal by the more frequent high-velocity interactions is taken into account

The dynamical evolution of a small number of protostars formed by the fragmentation of a dense molecular cloud core, taking into account in a simplified way the dissipative encounters between the protostars caused by the presence of their concomitant circumstellar discs, is analysed. For an initial system of two stars, the resulting binary parameters are determined essentially by the initial angular momentum and to a lesser extent by the rate at which the accretion discs evolve. For an initial system of three stars, however, a narrow range of initial conditions leads to a broad spread of final configurations

The properties of external accretion discs (discs with a central source of angular momentum) are explored both analytically and numerically. An illustrative example of the effect of a disc of material around a binary star on the stellar separation is considered. A Greens-function-type solution is considered, in which an initial ring of matter is put in orbit around the central binary.