Source author record

J. E. Pringle

J. E. Pringle appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2024arXiv

The Morphology of Asteroidal Dust Around White Dwarf Stars: Optical and Near-infrared Pulsations in G29-38

More than 36 years have passed since the discovery of the infrared excess from circumstellar dust orbiting the white dwarf G29-38, which at 17.5 pc it is the nearest and brightest of its class. The precise morphology of the orbiting dust remains only marginally constrained by existing data, subject to model-dependent inferences, and thus fundamental questions of its dynamical origin and evolution persist. This study presents a means to constrain the geometric distribution of the emitting dust using stellar pulsations measured at optical wavelengths as a variable illumination source of the dust, which re-radiates primarily in the infrared. By combining optical photometry from the Whole Earth Telescope with 0.7-2.5 micron spectroscopy obtained with SpeX at NASA's Infrared Telescope Facility, we detect luminosity variations at all observed wavelengths, with variations at most wavelengths corresponding to the behavior of the pulsating stellar photosphere, but towards the longest wavelengths the light curves probe the corresponding time-variability of the circumstellar dust. In addition to developing methodology, we find pulsation amplitudes decrease with increasing wavelength for principal pulsation modes, yet increase beyond approximately 2 microns for nonlinear combination frequencies. We interpret these results as combination modes deriving from principal modes of identical l values and discuss the implications for the morphology of the warm dust. We also draw attention to some discrepancies between our findings and theoretical expectations for the results of the non-linearity imposed by the surface convection zone on mode--mode interactions and on the behavior of the first harmonic of the highest-amplitude pulsation mode.

preprint2020arXiv

Emission from elliptical streams of dusty debris around white dwarfs

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.

preprint2020arXiv

Extreme variability in an active galactic nucleus: Gaia16aax

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α$ and $\rm Hβ$ 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.

preprint2020arXiv

Polarization measurements of the polluted white dwarf G29-38

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.

preprint2020arXiv

Short gamma-ray bursts and the decompression of neutron star matter in tidal streams

Short gamma-ray bursts (sGRBs) are generally thought to result from the merger of two neutron stars or the merger of a neutron star with a black hole. It is becoming standard practise to model these mergers with hydrodynamical simulations that employ equations of state that are derived, for example, for determining the behaviour of matter in core-collapse supernovae (CCSNe), and which therefore make use of the assumption that the matter is hot and in nuclear statistical equilibrium (NSE). In this Letter we draw attention to the fact that the hydrodynamical timescale (roughly the gravitational timescale of the neutron star) may be several orders of magnitude shorter than the timescale on which such equilibrium can be re-established in the tidal debris ejected during a sGRB, and that on the initial decompression timescales the unshocked tidal ejecta may remain sufficiently cool that the employed equations of state are not appropriate for modelling the dynamics of this part of the flow. On timescales short compared with the timescale on which NSE can be (re)established, the equation of state can remain relatively stiff and thus the stream of tidal debris can remain narrow and vulnerable to gravitational instability, as has recently been suggested. These findings suggest that estimates of the type and abundances of heavy elements formed in short gamma-ray bursts need to be revisited. We suggest that the most direct method of testing the physical and dynamical properties of tidal ejecta in sGRBs will come from modelling of their light curves, which provides the cleanest source of information on the system dynamics.

preprint2020arXiv

Strong surface outflows on accretion discs

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.

preprint2014arXiv

Stochastic accretion of planetesimals onto white dwarfs: constraints on the mass distribution of accreted material from atmospheric pollution

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

preprint2014arXiv

The frequency and nature of `cloud-cloud collisions' in galaxies

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.

preprint2013arXiv

Accretion disc viscosity: what do warped discs tell us?

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.

preprint2013arXiv

The dependence of stellar age distributions on GMC environment

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.

preprint2013arXiv

The Exciting Lives of Giant 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.

preprint2012arXiv

Dead Zones around Young Stellar Objects: FU Orionis Outbursts and Transition Discs

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.

preprint2012arXiv

Giant Molecular clouds: what are they made from, and how do they get there?

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.

preprint2012arXiv

Hundred Thousand Degree Gas in the Virgo Cluster of Galaxies

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.

preprint2012arXiv

Scars of Intense Accretion Episodes at Metal-Rich White Dwarfs

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.

preprint2011arXiv

Dead Zones around Young Stellar Objects: Dependence on Physical Parameters

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 Σ_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).

preprint2011arXiv

Disc instability in RS Ophiuchi: a path to Type Ia supernovae?

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.

preprint2011arXiv

On Identifying the Progenitors of Type Ia Supernovae

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.

preprint2011arXiv

The properties of the ISM in disc galaxies with stellar feedback

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.

preprint2011arXiv

Tidal Warping of Be Star Decretion Discs

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.

preprint2011arXiv

Why are most molecular clouds not gravitationally bound?

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.

preprint2010arXiv

A Common Envelope Binary Star Origin of Long Gamma-ray Bursts

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.

preprint2010arXiv

Retrograde Accretion and Merging Supermassive Black Holes

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.

preprint2010arXiv

The Accretion Disc Dynamo in the Solar Nebula

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.

preprint2010arXiv

The observable effects of tidally induced warps in protostellar discs

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.

preprint2009arXiv

A Simple Model for the Relationship Between Star Formation and Surface Density

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.

preprint2009arXiv

Black hole mergers: can gas discs solve the `final parsec' problem?

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.

preprint2009arXiv

Black hole mergers: the first light

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.

preprint2009arXiv

Chaotic star formation and the alignment of stellar rotation with disc and planetary orbital axes

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.

preprint2009arXiv

Discovery of CIV Emission Filaments in M87

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.

preprint2009arXiv

RS Ophiuchi: Thermonuclear Explosion or Disc Instability?

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.

preprint2009arXiv

Supernova Kicks and Misaligned Be Star Binaries

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, $σ_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, $σ_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.

preprint2009arXiv

Supernova Kicks and Misaligned Microquasars

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.

preprint2009arXiv

The Shape of an Accretion Disc in a Misaligned Black Hole Binary

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.

preprint2005arXiv

Spiral shocks, triggering of star formation and the velocity dispersion in Giant Molecular Clouds

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.

preprint2001arXiv

Competitive accretion in embedded stellar 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.

preprint2001arXiv

The Formation of Molecular Clouds

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.

preprint2000arXiv

Jet directions in Seyfert galaxies

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$μ$m flux. For each galaxy we have a pair of points ($i$,$δ$), 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 $β$ 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 $δ$ values can be well represented by a homogeneous $β-$distribution in the range $0^{\circ}\leqβ\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 ($|ϕ|$) is smaller than a given angle $ϕ_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.

preprint1995arXiv

Dynamo driven accretion discs and dwarf nova eruptions

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.

preprint1995arXiv

Energetics of star-disc encounters in the non-linear regime

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.

preprint1995arXiv

Magnetically modulated accretion in T Tauri stars

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.