Interfacial instabilities driven by self-gravity: first numerical results

The evolution of the interstellar medium (ISM) is driven by a variety of phenomena, including turbulence, shearing flows, magnetic fields and the thermal properties of the gas. Among the most important forces at work is self-gravity, which ultimately drives protostellar collapse. As part of an ongoing study of instabilities in the ISM, Hunter, Whitaker & Lovelace have discovered another process driven by self-gravity: the instability of an interface of discontinuous density. Theory predicts that this self-gravity driven interfacial instability persists in the static limit and in the absence of a constant background acceleration. Disturbances to a density interface are found to grow on a time-scale of the order of the free-fall time, even when the perturbation wavelength is much less than the Jeans length. Here we present the first numerical simulations of this instability. The theoretical growth rate is confirmed and the non-linear morphology displayed. The self-gravity interfacial instability is shown to be fundamentally different from the Rayleigh–Taylor instability, although both exhibit similar morphologies under the condition of a high density contrast, such as is commonly found in the ISM. Such instabilities are a possible mechanism by which observed features, such as the pillars of gas seen near the boundaries of interstellar clouds, are formed.     

The role of multipolar magnetic fields in pulsar magnetospheres

We explore the role of complex multipolar magnetic fields in determining physical processes near the surface of rotation powered pulsars. We model the actual magnetic field as the sum of global dipolar and star-centred multipolar fields. In configurations involving axisymmetric and uniform multipolar fields, 'neutral points' and 'neutral lines' exist close to the stellar surface. Also, the curvature radii of magnetic field lines near the stellar surface can never be smaller than the stellar radius, even for very high-order multipoles. Consequently, such configurations are unable to provide an efficient pair-creation process above pulsar polar caps, necessary for plasma mechanisms of generation of pulsar radiation. In configurations involving axisymmetric and non-uniform multipoles, the periphery of the pulsar polar cap becomes fragmented into symmetrically distributed narrow subregions where curvature radii of complex magnetic field lines are less than the radius of the star. The pair-production process is only possible just above these 'favourable' subregions. As a result, the pair plasma flow is confined within narrow filaments regularly distributed around the margin of the open magnetic flux tube. Such a magnetic topology allows us to model the system of 20 isolated subbeams observed in PSR B0943+10 by Deshpande & Rankin. We suggest a physical mechanism for the generation of pulsar radio emission in the ensemble of finite subbeams, based on specific instabilities. We propose an explanation for the subpulse drift phenomenon observed in some long-period pulsars.     

Location and parameters of a microwave millisecond spike event

A typical microwave millisecond spike event on November 2, 1997 was observed by the radio spectrograph of National Astronomical Observatories (NAOs) at 2.6–3.8 GHz with high time and frequency resolution. This event was also recorded by Nobeyama Radio Polarimeters (NoRP) at 1–35 GHz and Radio Heliograph (NoRH) at 17 GHz. The source at 17 GHz is located in one foot-point of a small bright coronal loop of YOHKOH SXT and SOHO EIT images with strong photospheric magnetic field in SOHO MDI magnetograph. It is assumed that the electron cyclotron maser instability and gyro-resonance absorption dominate, respectively, the rising and decay phase of the spike event. For different harmonic number of gyro-frequency or magnetic field strength, a fitting program with free plasma parameters is used to minimize the difference between the observational and theoretical values of the exponential growth and decay rates for a given spike. The plasma parameters at third harmonic number are more comparable to their typical values in solar corona. Hence, it is able to provide a diagnosis for the source parameters (magnetic field, density, and temperature), the properties of radiations (wave vector and propagation angle), and the properties of non-thermal electrons (density, pitch angle, and energy). The results are also comparable with the diagnosis of the gyro-synchrotron radiation model, the frequency drift rates and a dipole magnetic field model, as well as the YOHKOH SXT and SOHO MDI data. This study is supported by the NFSC project nos. 10333030 and 10273025, and “973” program with no. G2000078403.     

Tayler instability with Hall effect in young neutron stars

Collapse calculations indicate that the hot young neutron stars rotate differentially so that strong toroidal magnetic field components should exist in the outer shell where also the Hall effect appears to be important when the Hall parameter = ω_Bτ exceeds unity. The amplitudes of the induced toroidal magnetic fields are limited by the current‐induced Tayler instability. An important characteristics of the Hall effect is its distinct dependence on the sign of the magnetic field. We find for fast rotation that positive (negative) Hall parameters essentially reduce (increase) the stability domain. It is thus concluded that the toroidal field belts in young neutron stars induced by their differential rotation should have different amplitudes in both hemispheres which later are frozen in. Due to the effect of magnetic suppression of the heat conductivity also the brightness of the two hemispheres should be different. As a possible example for our scenario the isolated neutron star RBS 1223 is considered which has been found to exhibit different X‐ray brightness at both hemispheres (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Do superfluid instabilities prevent neutron star precession?

We discuss short wavelength (inertial wave) instabilities present in the standard two-fluid neutron star model when there is sufficient relative flow along the superfluid neutron vortex array. We demonstrate that these instabilities may be triggered in precessing neutron stars, since the angular velocity vectors of the neutron and proton fluids are misaligned during precession. Our results suggest that the standard (Eulerian) slow precession that results for weak drag between the vortices and the charged fluid (protons and electrons) is not seriously affected by the instability. In contrast, the fast precession, which results when vortices are strongly coupled to the charged component, is generally unstable. The presence of this instability renders the standard (solid body) rotation model for free precession inconsistent and makes unsafe conclusions that have recently been drawn regarding neutron star interiors based on observations of precession in radio pulsars.     

Quasi-viscous accretion flow – I. Equilibrium conditions and asymptotic behaviour

In a novel approach to studying viscous accretion flows, viscosity has been introduced as a perturbative effect, involving a first-order correction in the α-viscosity parameter. This method reduces the problem of solving a second-order non-linear differential equation (Navier–Stokes equation) to that of an effective first-order equation. Viscosity breaks down the invariance of the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes accretion from wind. Under a dynamical systems classification, the only feasible critical points of this 'quasi-viscous' flow are saddle points and spirals. On large spatial scales of the disc, where a linearized and radially propagating time-dependent perturbation is known to cause a secular instability, the velocity evolution equation of the quasi-viscous flow has been transformed to bear a formal closeness with Schrödinger's equation with a repulsive potential. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length-scale has been defined for the full spatial extent over which the accretion process would be viable.     

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 time-scale 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 time-scale 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.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 ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

Warping modes in discs around accreting neutron stars

The origin and stability of a thin sheet of plasma in the magnetosphere of an accreting neutron star are investigated. First, the radial extension of such a magnetospheric disc is explored. Then a mechanism for magnetospheric accretion is proposed, reconsidering the bending wave explored by Agapitou, Papaloizou & Terquem, that was found to be stable in ideal magnetohydrodynamics. We show that this warping becomes unstable and can reach high amplitudes, in a variant of Pringle's radiation-driven model for the warping of active galactic nucleus accretion discs. Finally, we discuss how this mechanism might give a clue to explain the observed X-ray kilohertz quasi-periodic oscillation of neutron star binaries.     

Acceleration of a spherical neutral shell produced by an ionization-shock front in an inhomogeneous interstellar medium

We numerically model the formation and acceleration of a neutral gas shell as an ionization-shock front propagates in a spherical cloud by taking into account the photoionization and radiative heating of the gas, the spectral radiative transfer. We suggest and implement an approximation of the cooling function that allows calculations to be performed in a wide range of gas ionization fractions and temperatures. The total mass, average velocity, and thickness of the shell have been determined. The results are compared with approximate formulas known in the literature. Based on the parameters of the shell found, we estimate its acceleration, characteristic scales, and the growth times of unstable perturbations. We analyze the influence of the cloud particle density, cloud radius, stellar temperature, and radiation spectrum on the integrated characteristics of the neutral gas in the layer between the ionization and shock fronts. The distribution of matter in the shell and its thickness are shown to differ significantly from those used in approximate models.