Studies of accretion flows around rotating black holes – II. Standing shocks in the pseudo-Kerr geometry

Standing, propagating or oscillating shock waves are common in accretion and winds around compact objects. We study the topology of all possible solutions using the pseudo-Kerr geometry. We present the parameter space spanned by the specific energy and angular momentum and compare it with that obtained from the full general relativity to show that the potential can work satisfactorily in fluid dynamics also, provided the polytropic index is suitably modified. We then divide the parameter space depending on the nature of the solution topology. We specifically study the nature of the standing Rankine–Hugoniot shocks. We also show that as the Kerr parameter is increased, the shock location generally moves closer to the black hole. In future, these solutions can be used as guidelines to test numerical simulations around compact objects.     

General relativistic electromagnetic fields of a slowly rotating magnetized neutron star – I. Formulation of the equations

We present analytic solutions of Maxwell equations in the internal and external background space–time of a slowly rotating magnetized neutron star. The star is considered isolated and in vacuum, with a dipolar magnetic field not aligned with the axis of rotation. With respect to a flat space–time solution, general relativity introduces corrections related both to the monopolar and the dipolar parts of the gravitational field. In particular, we show that in the case of infinite electrical conductivity general relativistic corrections resulting from the dragging of reference frames are present, but only in the expression for the electric field. In the case of finite electrical conductivity, however, corrections resulting from both the space–time curvature and the dragging of reference frames are shown to be present in the induction equation. These corrections could be relevant for the evolution of the magnetic fields of pulsars and magnetars. The solutions found, while obtained through some simplifying assumption, reflect a rather general physical configuration and could therefore be used in a variety of astrophysical situations.     

Topology of magnetars external field – I. Axially symmetric fields

There is an increasing theoretical and observational evidence that the external magnetic field of magnetars may contain a toroidal component, likely of the same order of the poloidal one. Such 'twisted magnetospheres' are threaded by currents flowing along the closed field lines which can efficiently interact with soft thermal photons via resonant cyclotron scatterings (RCS). Actually, RCS spectral models proved quite successful in explaining the persistent ∼1–10 keV emission from the magnetar candidates, the soft γ-ray repeaters (SGRs) and the anomalous X-ray pulsars (AXPs). Moreover, it has been proposed that, in the presence of highly relativistic electrons, the same process can give rise to the observed hard X-ray spectral tails extending up to  ∼200 keV  . Spectral calculations have been restricted up to now to the case of a globally twisted dipolar magnetosphere, although there are indications that the twist may be confined only to a portion of the magnetosphere, and/or that the large-scale field is more complex than a simple dipole. In this paper, we investigate multipolar, force–free magnetospheres of ultramagnetized neutron stars. We first discuss a general method to generate multipolar solutions of the Grad-Schlüter-Shafranov (GSS) equation, and analyse in detail dipolar, quadrupolar and octupolar fields. The spectra and lightcurves for these multipolar, globally twisted fields are then computed using a Monte Carlo code and compared with those of a purely dipolar configuration. Finally, the phase-resolved spectra and energy-dependent lightcurves obtained with a simple model of a locally sheared field are confronted with the International Gamma-Ray Astrophysics Laboratory ( INTEGRAL ) observations of the AXPs 1RXS J1708−4009 and 4U 0142+61. Results support a picture in which the field in these two sources is not globally twisted.     

General relativistic electromagnetic fields of a slowly rotating magnetized neutron star – II. Solution of the induction equations

We have solved numerically the general relativistic induction equations in the interior background space–time of a slowly rotating magnetized neutron star. The analytic form of these equations was discussed recently (Paper I), where corrections due to both the space–time curvature and the dragging of reference frames were shown to be present. Through a number of calculations we have investigated the evolution of the magnetic field with different rates of stellar rotation, different inclination angles between the magnetic moment and the rotation axis, as well as different values of the electrical conductivity. All of these calculations have been performed for a constant-temperature relativistic polytropic star and make use of a consistent solution of the initial-value problem which avoids the use of artificial analytic functions. Our results show that there exist general relativistic effects introduced by the rotation of the space–time which tend to decrease the decay rate of the magnetic field. The rotation-induced corrections are however generally hidden by the high electrical conductivity of the neutron star matter, and when realistic values for the electrical conductivity are considered, these corrections become negligible even for the fastest known pulsar.     

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.     

Newtonian hydrodynamics of the coalescence of black holes with neutron stars – II. Tidally locked binaries with a soft equation of state

We present a numerical study of the hydrodynamics in the final stages of inspiral of a black hole–neutron star binary, when the binary separation becomes comparable to the stellar radius. We use a Newtonian three-dimensional smooth particle hydrodynamics (SPH) code, and model the neutron star with a soft (adiabatic index Γ=5/3) polytropic equation of state, and the black hole as a Newtonian point mass that accretes matter via an absorbing boundary at the Schwarzschild radius. Our initial conditions correspond to tidally locked binaries in equilibrium, and we have explored configurations with different values of the mass ratio q M _NS M _BH, ranging from q =1 to 0.1. The dynamical evolution is followed for approximately 23 ms, and in every case studied here we find that the neutron star is tidally disrupted on a dynamical time-scale, forming a dense torus around the black hole that contains a few tenths of a solar mass. A nearly baryon-free axis is present in the system throughout the coalescence, and only modest beaming of a fireball that could give rise to a gamma-ray burst would be sufficient to avoid excessive baryon contamination. We find that some mass (of the order of 10⁻³–10⁻² M_⊙) may be dynamically ejected from the system, and could thus contribute substantially to the amount of observed r-process material in the galaxy. We calculate the gravitational radiation waveforms and luminosity emitted during the coalescence in the quadrupole approximation.     

Dual black holes in merger remnants – I. Linking accretion to dynamics

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or corotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an 'orbital angular momentum flip' due to the gas-dynamical friction, i.e. it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi–Hoyle–Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly 'lights up' at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ≲10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is corotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double active galactic nuclei at separations of ≲100 pc.     

Newtonian hydrodynamics of the coalescence of black holes with neutron stars – IV. Irrotational binaries with a soft equation of state

We present the results of three-dimensional hydrodynamical simulations of the final stages of in-spiral in a black hole–neutron star binary, when the separation is comparable to the stellar radius. We use a Newtonian smooth particle hydrodynamics (SPH) code to model the evolution of the system, and take the neutron star to be a polytrope with a soft (adiabatic indices     and     equation of state and the black hole to be a Newtonian point mass. The only non-Newtonian effect we include is a gravitational radiation back reaction force, computed in the quadrupole approximation for point masses. We use irrotational binaries as initial conditions for our dynamical simulations, which are begun when the system is on the verge of initiating mass transfer and followed for approximately 23 ms. For all the cases studied we find that the star is disrupted on a dynamical time-scale, and forms a massive     accretion torus around the spinning (Kerr) black hole. The rotation axis is clear of baryons (less than 10⁻⁵ M_⊙ within 10°) to an extent that would not preclude the formation of a relativistic fireball capable of powering a cosmological gamma-ray burst. Some mass (the specific amount is sensitive to the stiffness of the equation of state) may be dynamically ejected from the system during the coalescence and could undergo r-process nucleosynthesis. We calculate the waveforms, luminosities and energy spectra of the gravitational radiation signal, and show how they reflect the global outcome of the coalescence process.     

One-dimensional electric field structure of an outer gap accelerator – I. γ-ray production resulting from curvature radiation

We study the structure of a stationary and axisymmetric charge-deficient region (or potential gap) in the outer magnetosphere of a spinning neutron star. Assuming the existence of global current flow patterns in the magnetosphere, the charge depletion causes a large electric field along the magnetic field lines. This longitudinal electric field accelerates migratory electrons and/or positrons to ultrarelativistic energies. These relativistic electrons/positrons radiate γ -ray photons by curvature radiation. These γ -rays, in turn, produce yet more radiating particles by colliding with ambient X-ray photons, leading to a pair production cascade in the gap. The replenished charges partially screen the longitudinal electric field, which is self-consistently solved together with the distribution of e^± and γ -ray photons. We find the voltage drop in the gap as a function of the soft photon luminosity. It is demonstrated that the voltage drop is less than 3×10¹³ V when the background X-ray radiation is as luminous as Vela . However, this value increases with decreasing X-ray luminosity and attains 3×10¹⁵ V when the X-ray luminosity is as low as L _X=10³¹ erg s⁻¹.     

Mean-field dynamo in partially ionized plasmas – I

There are several astrophysical situations where one needs to study the dynamics of magnetic flux in partially ionized turbulent plasmas. In a partially ionized plasma, the magnetic induction is subjected to the ambipolar diffusion and the Hall effect in addition to the usual resistive dissipation. In this paper, we initiate the study of the kinematic dynamo in a partially ionized turbulent plasma. The Hall effect arises from the treatment of the electrons and the ions as two separate fluids and the ambipolar diffusion due to the inclusion of neutrals as the third fluid. It is shown that these non-ideal effects modify the so-called α effect and the turbulent diffusion coefficient β in a rather substantial way. The Hall effect may enhance or quench the dynamo action altogether. The ambipolar diffusion brings in an α which depends on the mean magnetic field. The new correlations embodying the coupling of the charged fluids and the neutral fluid appear in a decisive manner. The turbulence is necessarily magnetohydrodynamic with new spatial and time-scales. The nature of the new correlations is demonstrated by taking the Alfvénic turbulence as an example.     

Alfvén polar oscillations of relativistic stars

We study polar Alfvén oscillations of relativistic stars endowed with a strong global poloidal dipole magnetic field. Here, we focus only on the axisymmetric oscillations which are studied by numerically evolving the two-dimensional perturbation equations. Our study shows that the spectrum of the polar Alfvén oscillations is discrete in contrast to the spectrum of axial Alfvén oscillations which is continuous. We also show that the typical fluid modes, such as the f and p modes, are not significantly affected by the presence of the strong magnetic field.     

Testing the locality of transport in self-gravitating accretion discs – II. The massive disc case

We present hydrogen emission-line profile models of magnetospheric accretion on to classical T Tauri stars. The models are computed under the Sobolev approximation using the three-dimensional Monte Carlo radiative-transfer code torus . We have calculated four illustrative models in which the accretion flows are confined to azimuthal curtains – a geometry predicted by magnetohydrodynamic simulations. Properties of the line profile variability of our models are discussed, with reference to dynamic spectra and cross-correlation images. We find that some gross characteristics of the observed line profile variability are reproduced by our models, although in general the level of variability predicted is larger than that observed. We conclude that this excessive variability probably excludes dynamical simulations that predict accretion flows with low degrees of axisymmetry.     

Studies of accretion flows around rotating black holes – III. Shock oscillations and an estimation of the spin parameter from QPO frequencies

In the present communication of our series of papers dealing with the accretion flows in the pseudo-Kerr geometry, we discuss the effects of viscosity on the accretion flow around a rotating black hole. We find the solution topologies and give special attention to the solutions containing shocks. We draw the parameter space where standing shocks are possible and where the shocks could be oscillating and could produce quasi-periodic oscillations (QPOs) of X-rays observed from black hole candidates. In this model, the extreme locations of the shocks give the upper limits of the QPO frequencies  (ν_QPO)  which could be observed. We show that both the viscosity of the flow and the spin of the black hole a increase the QPO frequency while, as expected, the black hole mass reduces the QPO frequencies. Our major conclusion is that the highest observed frequency gives a strict lower limit of the spin. For instance, a black hole exhibiting  ν_QPO∼ 400  and  700 Hz  must have the spin parameters of   a > 0.25  and  >0.75  , respectively, provided viscosity of the flow is small. We discuss the implications of our results in the light of observations of QPOs from black hole candidates.