Iron line profiles including emission from within the innermost stable orbit of a black hole accretion disc

We propose a model for the source of the X-ray background (XRB) in which low-luminosity active nuclei ( L  ∼ 10⁴³ erg s⁻¹) are obscured ( N  ∼ 10²³ cm⁻²) by nuclear starbursts within the inner ∼ 100 pc. The obscuring material covers most of the sky as seen from the central source, rather than being distributed in a toroidal structure, and hardens the averaged X-ray spectrum by photoelectric absorption. The gas is turbulent with velocity dispersion ∼ few × 100 km s⁻¹ and cloud–cloud collisions lead to copious star formation. Although supernovae tend to produce outflows, most of the gas is trapped in the gravity field of the star-forming cluster itself and the central black hole. A hot ( T  ∼ 10⁶ − 10⁷ K) virialized phase of this gas, comprising a few per cent of the total obscuring material, feeds the central engine of ∼ 10⁷ M⊙ through Bondi accretion, at a sub-Eddington rate appropriate for the luminosity of these objects. If starburst-obscured objects give rise to the residual XRB, then only 10 per cent of the accretion in active galaxies occurs close to the Eddington limit in unabsorbed objects.     

On the MHD effects on the force-free monopole outflow

The stationary axisymmetric outflow from a rotating sphere with a (split) monopole magnetic field is considered. The stream equation describing the outflow is linearized in terms of the Michel magnetization parameter σ⁻¹ ≪ 1, which allows a self-consistent analysis of the direct problem. It is shown that for a finite σ the fast magnetosonic surface is located at a finite distance ∼ σ^1/3 R _L ( R _L =  c /Ω_F is the light cylinder). We have also found that the particle energy at the fast surface is just equal to the Michel value γ ∼ 1/3σ. The particle acceleration and magnetic field collimation are shown to become ineffective outside the fast magnetosonic surface.     

ROSAT ASCA observations of X-ray luminous starburst galaxies: NGC 3310 and 3690

We present ROSAT [High Resolution Imager (HRI) and Position Sensitive Proportional Counter (PSPC)] and ASCA observations of the two luminous ( L _x ∼ 10⁴¹⁻⁴² erg s⁻¹) star-forming galaxies NGC 3310 and 3690. The HRI shows clearly that the sources are extended with the X-ray emission in NGC 3690 coming from at least three regions. The combined 0.1–10 keV spectrum of NGC 3310 can be described by two components, a Raymond–Smith plasma with temperature kT  = 0.81^+0.09_−0.12 keV and a hard power law, Γ = 1.44^−0.20_−0.11 (or alternatively a harder Raymond–Smith plasma with kT  ∼ 15 keV), while there is no substantial excess absorption above the Galactic column value. The soft component emission is probably a super wind while the nature of the hard emission is more uncertain with the likely origins being X-ray binaries, inverse Compton scattering of infrared photons, an active galactic nucleus or a very hot gas component (∼10⁸ K). The spectrum of NGC 3690 is similar, with kT  = 0.83^+0.02_−0.04 keV and Γ = 1.56^+0.11_−0.11. We also employ more complicated models such as a multi-temperature thermal plasma, a non-equilibrium ionization code or the addition of a third softer component, which improve the fit but not at a statistically significant level (2σ). These results are similar to recent results on the archetypal star-forming galaxies M82 and NGC 253.     

A tenuous X-ray corona enveloping AE Aquarii

In this paper we propose that the observed unpulsed X-ray emission in AE Aquarii is the result of a very tenuous hot corona associated with the secondary star, which is pumped magnetohydrodynamically by the propeller action of the fast rotating white dwarf. It is shown that the closed coronal field of the secondary star envelops a substantial portion of the binary system, including the fast rotating magnetized white dwarf. This implies that the propeller outflow of material in AE Aquarii is initiated inside an enveloping magnetic cavity. The outflow crossing the secondary dead-zone field constitutes a  β_gen= (8πρ v ²_esc/ B ²) ≫ 1  plasma, acting as a magnetohydrodynamic generator resulting in the induction of field-aligned currents in these closed magnetospheric circuits where  β_cir= (8π nkT / B ²) ≪ 1  . The Ohmic heating of the coronal circuit can readily account for a   T _x≥ 10⁷  K  plasma in the coronal flux tubes connecting the generator and the stellar surface. Further, the bremsstrahlung losses of the thermal electrons in the coronal circuit can readily drive the observed unpulsed X-ray luminosity of   L _x∼ 10³¹  erg s ⁻¹  , which correlates with the luminosity and relatively large source implied by recent XMM–Newton observations.     

Self-similar spherical collapse with non-radial motions

We derive the asymptotic mass profile near the collapse centre of an initial spherical density perturbation, δ ∝ M ^{− ε} , of collisionless particles with non-radial motions. We show that angular momenta introduced at the initial time do not affect the mass profile. Alternatively, we consider a scheme in which a particle moves on a radial orbit until it reaches its turnaround radius, r ∗. At turnaround the particle acquires an angular momentum L =ℒ√ GM _* r _* per unit mass, where M ∗ is the mass interior to r ∗. In this scheme, the mass profile is M ∝ r ^{3/(1+3 ε )} for all ε >0 , in the region r / r _t ≪ℒ , where r _t is the current turnaround radius. If ℒ≪1 then the profile in the region ℒ≪ r / r _t ≪1 is M ∝ r for ε <2/3 , and remains M ∝ r ^{3/(1+3 ε )} for ε ≥2/3 . The derivation relies on a general property of non-radial orbits which is that the ratio of the pericentre to apocentre is constant in a force field k ( t ) r ⁿ with k ( t ) varying adiabatically.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Testing fundamental physics with distant star clusters: theoretical models for pressure-supported stellar systems

We investigate the mean velocity dispersion and the velocity dispersion profile of stellar systems in modified Newtonian dynamics (MOND), using the N -body code n-mody , which is a particle-mesh-based code with a numerical MOND potential solver developed by Ciotti, Londrillo & Nipoti. We have calculated mean velocity dispersions for stellar systems following Plummer density distributions with masses in the range of 10⁴ to  10⁹ M_⊙  and which are either isolated or immersed in an external field. Our integrations reproduce previous analytic estimates for stellar velocities in systems in the deep MOND regime  ( a _i, a _e≪ a ₀)  , where the motion of stars is either dominated by internal accelerations  ( a _i≫ a _e)  or constant external accelerations  ( a _e≫ a _i)  . In addition, we derive for the first time analytic formulae for the line-of-sight velocity dispersion in the intermediate regime  ( a _i∼ a _e∼ a ₀)  . This allows for a much-improved comparison of MOND with observed velocity dispersions of stellar systems. We finally derive the velocity dispersion of the globular cluster Pal 14 as one of the outer Milky Way halo globular clusters that have recently been proposed as a differentiator between Newtonian and MONDian dynamics.     

A Chandra observation of the distant radio galaxy B2 0902+343: a powerful obscured active galaxy

We report a Chandra observation of the   z =3.395  radio galaxy B2 0902+343. The unresolved X-ray source is centred on the active nucleus. The spectrum is well fitted by a flat power law of photon index of  Γ∼1.1  with intrinsic absorption of  8×10²² cm^-2  , and an intrinsic  2–10 keV  luminosity of  3.3×10⁴⁵ erg s^-1  . More complex models that allow for a steeper spectral index cause the column density and intrinsic luminosity to increase. The data limit any thermal luminosity of the hot magnetized medium, assumed responsible for high Faraday rotation measures seen in the radio source, to less than ∼10⁴⁵ erg s⁻¹.     

The 'Carina Flare' supershell: probing the atomic and molecular ISM in a Galactic chimney

The 'Carina Flare' supershell, GSH 287+04−17, is a molecular supershell originally discovered in  ¹²CO( J = 1–0)  with the NANTEN 4 m telescope. We present the first study of the shell's atomic ISM, using H  i 21-cm line data from the Parkes 64-m telescope Southern Galactic Plane Survey. The data reveal a gently expanding,  ∼230 × 360  pc H  i supershell that shows strong evidence of Galactic Plane blowout, with a break in its main body at   z ∼ 280  pc and a capped high-latitude extension reaching   z ∼ 450  pc. The molecular clouds form comoving parts of the atomic shell, and the morphology of the two phases reflects the supershell's influence on the structure of the ISM. We also report the first discovery of an ionized component of the supershell, in the form of delicate, streamer-like filaments aligned with the proposed direction of blowout. The distance estimate to the shell is re-examined, and we find strong evidence to support the original suggestion that it is located in the Carina Arm at a distance of  2.6 ± 0.4 kpc  . Associated H  i and H₂ masses are estimated as   M _H I≈ 7 ± 3 × 10⁵ M_⊙  and     , and the kinetic energy of the expanding shell as   E _K ∼ 1 × 10⁵¹  erg. We examine the results of analytical and numerical models to estimate a required formation energy of several 10⁵¹ to  ∼10⁵²  erg, and an age of  ∼10⁷ yr  . This age is compatible with molecular cloud formation time-scales, and we briefly consider the viability of a supershell-triggered origin for the molecular component.     

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

Radiation pressure acceleration by X-rays in active galactic nuclei

We present calculations of the dynamics of highly ionized gas (HIG) clouds that are confined by external pressure, and are photoionized by AGN continuum. We focus on the gas that is seen, in absorption, in the X-ray spectrum of many AGN and show that such gas can reach hydrostatic equilibrium under various conditions. The principal conclusion is that the clouds can be accelerated to high velocities by the central X-ray source. The dynamical problem can be reduced to the calculation of a single parameter, the average force multiplier, 〈 M 〉. The typical value of 〈 M 〉 is ∼10 suggesting that radiation pressure acceleration by X-rays is efficient for L / L _Edd≳0.1 . The terminal velocity scales with the escape velocity at the base of the flow and can exceed it by a large factor. The typical velocity for a HIG flow that originates at R =10¹⁷ cm in a source with L _x =10⁴⁴ erg s⁻¹ is ∼1000 km s⁻¹, i.e. similar to the velocities observed in several X-ray and UV absorption systems.
Highly ionized AGN clouds are driven mainly by bound–free absorption, and bound–bound processes are less important unless the lines are significantly broadened or the column density is very small. Pressure laws that result in constant or outward decreasing ionization parameters are most effective in accelerating the flow.     

Episodic accretion in magnetically layered protoplanetary discs

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 on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ 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 r ∼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.     

The census of nuclear activity of late-type galaxies in the Virgo cluster

The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M _dyn > 10^8.5 M_⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M _dyn≳ 10¹⁰ M_⊙  . Their frequency increases steeply with the dynamical mass from zero at   M _dyn≈ 10^9.5 M_⊙  to virtually 1 at   M _dyn > 10^11.5 M_⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼10⁶ M_⊙  are found active in our sample.     

Explosive reconnection in magnetars

The X-ray activity of anomalous X-ray pulsars and soft γ-ray repeaters may result from the heating of their magnetic corona by direct currents dissipated by magnetic reconnection. We investigate the possibility that X-ray flares and bursts observed from anomalous X-ray pulsars and soft γ-ray repeaters result from magnetospheric reconnection events initiated by development of the tearing mode in magnetically dominated relativistic plasma. We formulate equations of resistive force-free electrodynamics, discuss the relation of the latter to ideal electrodynamics, and give examples of both ideal and resistive equilibria. Resistive force-free current layers are unstable towards the development of small-scale current sheets where resistive effects become important. Thin current sheets are found to be unstable due to the development of the resistive force-free tearing mode. The growth rate of the tearing mode is intermediate between the short Alfvén time-scale  τ_A  and a long resistive time-scale  τ_R: Γ∼ 1/(τ_Rτ_A)^1/2  , similar to the case of non-relativistic non-force-free plasma. We propose that growth of the tearing mode is related to the typical rise time of flares, ∼10 ms. Finally, we discuss how reconnection may explain other magnetar phenomena and ways to test the model.     

The high-velocity molecular stream in W51B: a ring structure with compact H ii regions

We present ¹³ CO J  = 1 − 0 line observations of the H  ii region complex W51B located in the high-velocity (HV) stream. These observations reveal a filamentary and clumpy structure in the molecular gas. The mean local standard of rest (LSR) velocity ∼ + 65 km s⁻¹ of the molecular gas in this region is greater than the maximum velocities allowed by kinematic Galactic rotation curves. The size and mass of the molecular clouds are ∼ 48 × 17 pc² and ∼ 2.4 × 10⁵ M⊙ respectively. In a position–velocity diagram, molecular gas in the southern part comprises a redshifted ring structure with v _LSR=+ 60 to +73 km s⁻¹. The velocity gradient of this ring is ∼ 0.5 km s⁻¹ pc⁻¹, and the mass is ∼ 6.2 × 10⁴ M⊙. If we assume that the ring is expanding with a uniform velocity, the expansion velocity, radius and kinetic energy are ∼ 7 km s⁻¹, ∼ 13 pc and ∼ 3.0 × 10 ⁴⁹ erg respectively. The kinetic energy and mass spectrum of the ring could be explained by an expanding cylindrical cloud with a centrally condensed mass distribution. The locations of two compact H  ii regions, G49.0−0.3 and G48.9−0.3, coincide with the two molecular clumps in this ring. We discuss star formation, and the mechanism that produced the ring structure.     

The X-ray luminosity function of AGN at z∼ 3

We combine Lyman-break colour selection with ultradeep (≳200 ks) Chandra X-ray imaging over a survey area of ∼0.35 deg² to select high-redshift active galactic nuclei (AGN). Applying careful corrections for both the optical and X-ray selection functions, the data allow us to make the most accurate determination to date of the faint end of the X-ray luminosity function (XLF) at   z ∼ 3  . Our methodology recovers a number density of X-ray sources at this redshift which is at least as high as previous surveys, demonstrating that it is an effective way of selecting high z AGN. Comparing to results at   z = 1  , we find no evidence that the faint slope of the XLF flattens at high z , but we do find significant (factor ∼3.6) negative evolution of the space density of low luminosity AGN. Combining with bright end data from very wide surveys we also see marginal evidence for continued positive evolution of the characteristic break luminosity   L _*  . Our data therefore support models of luminosity-dependent density evolution between   z = 1  and   z = 3  . A sharp upturn in the the XLF is seen at the very lowest luminosities  ( L _X≲ 10^42.5 erg s⁻¹)  , most likely due to the contribution of pure X-ray starburst galaxies at very faint fluxes.     

A unified picture for the γ-ray and prompt optical emissions of GRB 990123

The prompt optical emission of GRB 990123 was uncorrelated to the γ-ray light curve and exhibited temporal properties similar to those of the steeply decaying, early X-ray emission observed by Swift at the end of many bursts. These facts suggest that the optical counterpart of GRB 990123 was the large-angle emission released during (the second pulse of) the burst. If the optical and γ-ray emissions of GRB 990123 have, indeed, the same origin then their properties require that (i) the optical counterpart was synchrotron emission and γ-rays arose from inverse-Compton scatterings (the 'synchrotron self-Compton model'), (ii) the peak energy of the optical-synchrotron component was at ∼20 eV and (iii) the burst emission was produced by a relativistic outflow moving at Lorentz factor  ≳450  and at a radius  ≳10¹⁵  cm, which is comparable to the outflow deceleration radius. Because the spectrum of GRB 990123 was optically thin above 2 keV, the magnetic field behind the shock must have decayed on a length-scale of  ≲1  per cent  of the thickness of the shocked gas, which corresponds to  10⁶–10⁷  plasma skin depths. Consistency of the optical counterpart decay rate and its spectral slope (or that of the burst, if they represent different spectral components) with the expectations for the large-angle burst emission represents the most direct test of the unifying picture proposed here for GRB 990123.     

On the evolution of the black hole: spheroid mass ratio

We present the results of a study which uses the 3C RR sample of radio-loud active galactic nuclei to investigate the evolution of the black hole:spheroid mass ratio in the most massive early-type galaxies from  0 < z < 2  . Radio-loud unification is exploited to obtain virial (linewidth) black hole mass estimates from the 3C RR quasars, and stellar mass estimates from the 3C RR radio galaxies, thereby providing black hole and stellar mass estimates for a single population of early-type galaxies. At low redshift  ( z ≲ 1)  , the 3C RR sample is consistent with a black hole:spheroid mass ratio of   M _bh/ M _sph≃ 0.002  , in good agreement with that observed locally for quiescent galaxies of similar stellar mass  ( M _sph≃ 5 × 10¹¹ M_⊙)  . However, over the redshift interval  0 < z < 2  the 3C RR black hole:spheroid mass ratio is found to evolve as   M _bh/ M _sph∝ (1 + z )^2.07±0.76  , reaching   M _bh/ M _sph≃ 0.008  by redshift   z ≃ 2  . This evolution is found to be inconsistent with the local black hole:spheroid mass ratio remaining constant at a moderately significant level (98 per cent). If confirmed, the detection of evolution in the 3C RR black hole:spheroid mass ratio further strengthens the evidence that, at least for massive early-type galaxies, the growth of the central supermassive black hole may be completed before that of the host spheroid.