Model dependence of transonic properties of accretion flows around black holes

We analytically study how the behaviour of accretion flows changes when the flow model is varied. We study the transonic properties of the conical flow, a flow of constant height and a flow in vertical equilibrium, and show that all these models are basically identical, provided that the polytropic constant is suitably changed from one model to another. We show that this behaviour is extendible even when standing shocks are produced in the flow. The parameter space where shocks are produced remains roughly identical in all these models when the same transformation among the polytropic indices is used. We present applications of these findings.     

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.     

Chaotic zone boundary for low free eccentricity particles near an eccentric planet

We consider particles with low free or proper eccentricity that are orbiting near planets on eccentric orbits. Through collisionless particle integration, we numerically find the location of the boundary of the chaotic zone in the planet's corotation region. We find that the distance in semimajor axis between the planet and boundary depends on the planet mass to the 2/7 power and is independent of the planet eccentricity, at least for planet eccentricities below 0.3. Our integrations reveal a similarity between the dynamics of particles at zero eccentricity near a planet in a circular orbit and with zero free eccentricity particles near an eccentric planet. The 2/7th law has been previously explained by estimating the semimajor at which the first-order mean motion resonances are large enough to overlap. Orbital dynamics near an eccentric planet could differ due to first-order corotation resonances that have strength proportional to the planet's eccentricity. However, we find that the corotation resonance width at low free eccentricity is small; also the first-order resonance width at zero free eccentricity is the same as that for a zero-eccentricity particle near a planet in a circular orbit. This accounts for insensitivity of the chaotic zone width to planet eccentricity. Particles at zero free eccentricity near an eccentric planet have similar dynamics to those at zero eccentricity near a planet in a circular orbit.     

Superhumps: confronting theory with observation

We review the theory and observations related to the 'superhump' precession of eccentric accretion discs in close binary systems. We agree with earlier work, although for different reasons, that the discrepancy between observation and dynamical theory implies that the effect of pressure in the disc cannot be neglected. We extend earlier work that investigates this effect to include the correct expression for the radius at which resonant orbits occur. Using analytic expressions for the accretion disc structure, we derive a relationship between the period excess and mass ratio with the pressure effects included. This is compared to the observed data, recently derived results for detailed integration of the disc equations and the equivalent empirically derived relations and used to predict values for the mass ratio based on measured values of the period excess for 88 systems.     

Radio and X-ray properties of relativistic beaming models for ultraluminous X-ray sources

We calculate the broad-band radio–X-ray spectra predicted by microblazar and microquasar models for ultraluminous X-ray sources (ULXs), exploring the possibility that their dominant power-law component is produced by a relativistic jet, even at near-Eddington mass accretion rates. We do this by first constructing a generalized disc–jet theoretical framework in which some fraction of the total accretion power, P _a, is efficiently removed from the accretion disc by a magnetic torque responsible for jet formation. Thus, for different black hole masses, mass accretion rates and magnetic coupling strength, we self-consistently calculate the relative importance of the modified disc spectrum, as well as the overall jet emission due to synchrotron and Compton processes. In general, transferring accretion power to a jet makes the disc fainter and cooler than a standard disc at the same mass accretion rate; this may explain why the soft spectral component appears less prominent than the dominant power-law component in most bright ULXs. We show that the apparent X-ray luminosity and spectrum predicted by the microquasar model are consistent with the observed properties of most ULXs. We predict that the radio synchrotron jet emission is too faint to be detected at the typical threshold of radio surveys to date. This is consistent with the high rate of non-detections over detections in radio counterpart searches. Conversely, we conclude that the observed radio emission found associated with a few ULXs cannot be due to beamed synchrotron emission from a relativistic jet.     

The origin of the initial mass function and its dependence on the mean Jeans mass in molecular clouds

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of 9 and 3, respectively.
We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 au), one of which is a wide binary brown dwarf system.
Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high masses may depend on the dispersion in the accretion rates of protostars.     

Physics of accretion in the millisecond pulsar XTE J1751 − 305

We have undertaken an extensive study of X-ray data from the accreting millisecond pulsar XTE J1751 − 305 observed by RXTE and XMM–Newton during its 2002 outburst. In all aspects this source is similar to the prototypical millisecond pulsar SAX J1808.4 − 3658, except for the higher peak luminosity of 13 per cent of Eddington, and the optical depth of the hard X-ray source, which is larger by a factor ∼2. Its broad-band X-ray spectrum can be modelled by three components. We interpret the two soft components as thermal emission from a colder  ( kT ∼ 0.6 keV)  accretion disc and a hotter (∼1 keV) spot on the neutron star surface. We interpret the hard component as thermal Comptonization in plasma of temperature ∼40 keV and optical depth ∼1.5 in a slab geometry. The plasma is heated by the accretion shock as the material collimated by the magnetic field impacts on to the surface. The seed photons for Comptonization are provided by the hotspot, not by the disc. The Compton reflection is weak and the disc is probably truncated into an optically thin flow above the magnetospheric radius. Rotation of the emission region with the star creates an almost sinusoidal pulse profile with an rms amplitude of 3.3 per cent. The energy-dependent soft phase lags can be modelled by two pulsating components shifted in phase, which is naturally explained by a different character of emission of the optically thick spot and optically thin shock combined with the action of the Doppler boosting. The observed variability amplitude constrains the hotspot to lie within 3°–4° of the rotational pole. We estimate the inner radius of the optically thick accreting disc to be about 40 km. In that case, the absence of emission from the antipodal spot, which can be blocked by the accretion disc, gives the inclination of the system as ≳70°.