Warped discs and the directional stability of jets in Active Galactic Nuclei

Warped accretion discs in Active Galactic Nuclei (AGN) exert a torque on the black hole that tends to align the rotation axis with the angular momentum of the outer disc. We compute the magnitude of this torque by solving numerically for the steady state shape of the warped disc, and verify that the analytic solution of Scheuer and Feiler (1996) provides an excellent approximation. We generalise these results for discs with strong warps and arbitrary surface density profiles, and calculate the timescale on which the black hole becomes aligned with the angular momentum in the outer disc. For massive holes and accretion rates of the order of the Eddington limit the alignment timescale is always short (less than a Myr), so that jets accelerated from the inner disc region provide a prompt tracer of the angular momentum of gas at large radii in the disc. Longer timescales are predicted for low luminosity systems, depending on the degree of anisotropy in the disc's hydrodynamic response to shear and warp, and for the final decay of modest warps at large radii in the disc that are potentially observable via VLBI. We discuss the implications of this for the inferred accretion history of those Active Galactic Nuclei whose jet directions appear to be stable over long timescales. The large energy deposition rate at modest disc radii during rapid realignment episodes should make such objects transiently bright at optical and infrared wavelengths.Comment: MNRAS, in press. Revised to match accepted version, with one new figure showing alignment timescale as a function of black hole mas

The alignment of disk and black hole spins in active galactic nuclei

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

Alignment Timescale of the Microquasar GRO J1655-40

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

Lense-Thirring precession of accretion disks around compact objects

Misaligned accretion disks surrounding rotating compact objects experience a torque due to the Lense-Thirring effect, which leads to precession of the inner disk. It has been suggested that this effect could be responsible for some low frequency Quasi-Periodic Oscillations observed in the X-ray lightcurves of neutron star and galactic black hole systems. We investigate this possibility via time-dependent calculations of the response of the inner disk to impulsive perturbations for both Newtonian point mass and Paczynski-Wiita potentials, and compare the results to the predictions of the linearized twisted accretion disk equations. For most of a wide range of disk models that we have considered, the combination of differential precession and viscosity causes the warps to decay extremely rapidly. Moreover, at least for relatively slowly rotating objects, linear calculations in a Newtonian point mass potential provide a good measure of the damping rate, provided only that the timescale for precession is much shorter than the viscous time in the inner disk. The typically rapid decay rates suggest that coherent precession of a fluid disk would not be observable, though it remains possible that the damping rate of warp in the disk could be low enough to permit weakly coherent signals from Lense-Thirring precession.Comment: ApJ, in press. Minor revisions to match accepted version. Animations showing warp evolution are available at http://www.cita.utoronto.ca/~armitage/lense_thirring.htm

Evaluation of a global aerosol microphysics model against size-resolved particle statistics in the marine atmosphere

A statistical synthesis of marine aerosol measurements from experiments in four different oceans is used to evaluate a global aerosol microphysics model (GLOMAP). We compare the model against observed size resolved particle concentrations, probability distributions, and the temporal persistence of different size particles. We attempt to explain the observed sub-micrometre size distributions in terms of sulfate and sea spray and quantify the possible contributions of anthropogenic sulfate and carbonaceous material to the number and mass distribution. The model predicts a bimodal size distribution that agrees well with observations as a grand average over all regions, but there are large regional differences. Notably, observed Aitken mode number concentrations are more than a factor 10 higher than in the model for the N Atlantic but a factor 7 lower than the model in the NW Pacific. We also find that modelled Aitken mode and accumulation mode geometric mean diameters are generally smaller in the model by 10–30%. Comparison with observed free tropospheric Aitken mode distributions suggests that the model underpredicts growth of these particles during descent to the marine boundary layer (MBL). Recent observations of a substantial organic component of free tropospheric aerosol could explain this discrepancy. We find that anthropogenic continental material makes a substantial contribution to N Atlantic MBL aerosol, with typically 60–90% of sulfate across the particle size range coming from anthropogenic sources, even if we analyse air that has spent an average of >120 h away from land. However, anthropogenic primary black carbon and organic carbon particles (at the emission size and quantity assumed here) do not explain the large discrepancies in Aitken mode number. Several explanations for the discrepancy are suggested. The lack of lower atmospheric particle formation in the model may explain low N Atlantic particle concentrations. However, the observed and modelled particle persistence at Cape Grim in the Southern Ocean, does not reveal a diurnal cycle consistent with a photochemically driven local particle source. We also show that a physically based cloud drop activation scheme better explains the observed change in accumulation mode geometric mean diameter with particle number

Theoretical Modeling of the Thermal State of Accreting White Dwarfs Undergoing Classical Novae

White dwarfs experience a thermal renaissance when they receive mass from a stellar companion in a binary. For accretion rates < 10^-8 Msun/yr, the freshly accumulated hydrogen/helium envelope ignites in a thermally unstable manner that results in a classical novae (CN) outburst and ejection of material. We have undertaken a theoretical study of the impact of the accumulating envelope on the thermal state of the underlying white dwarf (WD). This has allowed us to find the equilibrium WD core temperatures (T_c), the classical nova ignition masses (M_ign) and the thermal luminosities for WDs accreting at rates of 10^-11 - 10^-8 Msun/yr. These accretion rates are most appropriate to WDs in cataclysmic variables (CVs) of P_orb <~ 7 hr, many of which accrete sporadically as dwarf novae. We have included ^3He in the accreted material at levels appropriate for CVs and find that it significantly modifies the CN ignition mass. We compare our results with several others from the CN literature and find that the inclusion of ^3He leads to lower M_ign for >~ 10^-10 Msun/yr, and that for below this the particular author's assumption concerning T_c, which we calculate consistently, is a determining factor. Initial comparisons of our CN ignition masses with measured ejected masses find reasonable agreement and point to ejection of material comparable to that accreted.Comment: 14 pages, 11 figures; uses emulateapj; accepted by the Astrophysical Journal; revised for clarity, added short discussion of diffusio

Radio Continuum Evidence for Outflow and Absorption in the Seyfert 1 Galaxy Markarian 231

The VLBA and the VLA have been used to image the continuum radio emission from Mrk 231, a Seyfert 1 galaxy and the brightest infrared galaxy in the local universe. The smallest scales reveal a double source less than 2 pc in extent. The components of this central source have minimum brightness temperatures of 10^9 to 10^{10} K, spectral turnovers between 2 and 10 GHz, and appear to define the galaxy nucleus plus the inner regions of a jet. The components may be free-free absorbed or synchtrotron self-absorbed. On larger scales, the images confirm a previously known north-south triple source extending 40 pc and elongated perpendicular to a 350-pc starburst disk. Both lobes show evidence for free-free absorption near 2 GHz, probably due to ionized gas with a density of 1-2 X 10^3 cm^{-3} in the innermost parts of the starburst disk. The absorbing gas may be ionized by the active nucleus or by local regions of enhanced star formation. The elongation of the 40-pc triple differs by 65 deg from that of the 2-pc source. The different symmetry axes on different scales imply strong curvature in the inner part of the radio jet. The radio continuum from the 350-pc disk has a spectral index near -0.4 above 1.4 GHz and is plausibly energized by a massive burst of star formation. On VLA scales, asymmetric and diffuse emission extends for more than 25 kpc. This emission has a steep spectrum, linear polarization exceeding 50% at some locations, and shares the symmetry axis of the 40-pc triple. The diffuse radio source is probably generated by energy deposition from a slow-moving nuclear jet, which conceivably could help energize the off-nuclear starburst as well.Comment: 34 pages, 7 Postscript figures, LaTeX file in AASTeX format, accepted in ApJ, Vol. 516, May 1, 199

Megamaser Disks in Active Galactic Nuclei

Recent spectroscopic and VLBI-imaging observations of bright extragalactic water maser sources have revealed that the megamaser emission often originates in thin circumnuclear disks near the centers of active galactic nuclei (AGNs). Using general radiative and kinematic considerations and taking account of the observed flux variability, we argue that the maser emission regions are clumpy, a conclusion that is independent of the detailed mechanism (X-ray heating, shocks, etc.) driving the collisionally pumped masers. We examine scenarios in which the clumps represent discrete gas condensations (i.e., clouds) and do not merely correspond to velocity irregularities in the disk. We show that even two clouds that overlap within the velocity coherence length along the line of sight could account (through self-amplification) for the entire maser flux of a high-velocity ``satellite'' feature in sources like NGC 4258 and NGC 1068, and we suggest that cloud self-amplification likely contributes also to the flux of the background-amplifying ``systemic'' features in these objects. Analogous interpretations have previously been proposed for water maser sources in Galactic star-forming regions. We argue that this picture provides a natural explanation of the time-variability characteristics of extragalactic megamaser sources and of their apparent association with Seyfert 2-like galaxies. We also show that the requisite cloud space densities and internal densities are consistent with the typical values of nuclear (broad emission-line region-type) clouds.Comment: 55 pages, 7 figures, AASTeX4.0, to appear in The Astrophysical Journal (1999 March 1 issue

A Radio Study of the Seyfert galaxy Markarian 6: Implications for Seyfert life-cycles

We have carried out an extensive radio study with the Very Large Array on the Seyfert 1.5 galaxy Mrk 6 and imaged a spectacular radio structure in the source. The radio emission occurs on three different spatial scales, from ~7.5 kpc bubbles to ~1.5 kpc bubbles lying nearly orthogonal to them and a ~1 kpc radio jet lying orthogonal to the kpc-scale bubble. To explain the complex morphology, we first consider a scenario in which the radio structures are the result of superwinds ejected by a nuclear starburst. However, recent Spitzer observations of Mrk 6 provide an upper limit to the star formation rate (SFR) of ~5.5 M_sun/yr, an estimate much lower than the SFR of ~33 M_sun/yr derived assuming that the bubbles are a result of starburst winds energized by supernovae explosions. Thus, a starburst alone cannot meet the energy requirements for the creation of the bubbles in Mrk 6. We show that a single plasmon model is energetically infeasible, and we argue that a jet-driven bubble model while energetically feasible does not produce the complex radio morphologies. Finally, we consider a model in which the complex radio structure is a result of an episodically-powered precessing jet that changes its orientation. This model is the most attractive as it can naturally explain the complex radio morphology, and is consistent with the energetics, the spectral index and the polarization structure. Radio emission in this scenario is a short-lived phenomenon in the lifetime of a Seyfert galaxy which results due to an accretion event.Comment: Accepted for publication in Ap

Relative periodic orbits in transitional pipe flow

A dynamical system description of the transition process in shear flows with no linear instability starts with a knowledge of exact coherent solutions, among them travelling waves (TWs) and relative periodic orbits (RPOs). We describe a numerical method to find such solutions in pipe flow and apply it in the vicinity of a Hopf bifurcation from a TW which looks to be especially relevant for transition. The dominant structural feature of the RPO solution is the presence of weakly modulated streaks. This RPO, like the TW from which it bifurcates, sits on the laminar-turbulent boundary separating initial conditions which lead to turbulence from those which immediately relaminarise