Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

In a hot bubble: why does superbubble feedback work, but isolated supernovae do not?

Using idealized one-dimensional Eulerian hydrodynamic simulations, we contrast the behaviour of isolated supernovae with the superbubbles driven by multiple, collocated supernovae. Continuous energy injection via successive supernovae exploding within the hot/dilute bubble maintains a strong termination shock. This strong shock keeps the superbubble over-pressured and drives the outer shock well after it becomes radiative. Isolated supernovae, in contrast, with no further energy injection, become radiative quite early (less than or similar to 0.1Myr, tens of pc), and stall at scales less than or similar to 100 pc. We show that isolated supernovae lose almost all of their mechanical energy by 1 Myr, but superbubbles can retain up to similar to 40 per cent of the input energy in the form of mechanical energy over the lifetime of the star cluster (a few tens of Myr). These conclusions hold even in the presence of realistic magnetic fields and thermal conduction. We also compare various methods for implementing supernova feedback in numerical simulations. For various feedback prescriptions, we derive the spatial scale below which the energy needs to be deposited in order for it to couple to the interstellar medium. We show that a steady thermal wind within the superbubble appears only for a large number (greater than or similar to 10(4)) of supernovae. For smaller clusters, we expect multiple internal shocks instead of a smooth, dense thermalized wind.

Angular momentum transport and evolution of lopsided galaxies

The surface brightness distribution in the majority of stellar galactic discs falls off exponentially. Often what lies beyond such a stellar disc is the neutral hydrogen gas whose distribution also follows a nearly exponential profile at least for a number of nearby disc galaxies. Both the stars and gas are commonly known to host lopsided asymmetry especially in the outer parts of a galaxy. The role of such asymmetry in the dynamical evolution of a galaxy has not been explored so far. Following Lindblad's original idea of kinematic density waves, we show that the outer part of an exponential disc is ideally suitable for hosting lopsided asymmetry. Further, we compute the transport of angular momentum in the combined stars and gas disc embedded in a dark matter halo. We show that in a pure star and gas disc, there is a transition point where the free precession frequency of a lopsided mode, Omega - kappa, changes from retrograde to prograde and this in turn reverses the direction of angular momentum flow in the disc leading to an unphysical behaviour. We show that this problem is overcome in the presence of a dark matter halo, which sets the angular momentum flow outwards as required for disc evolution, provided the lopsidedness is leading in nature. This, plus the well-known angular momentum transport in the inner parts due to spiral arms, can facilitate an inflow of gas from outside perhaps through the cosmic filaments.

Deep searches for decametre-wavelength pulsed emission from radio-quiet gamma-ray pulsars

We report the results of extensive follow-up observations of the gamma-ray pulsar J1732-3131, which has recently been detected at decametre wavelengths, and the results of deep searches for the counterparts of nine other radio-quiet gamma-ray pulsars at 34 MHz, using the Gauribidanur radio telescope. No periodic signal from J1732-3131 could be detected above a detection threshold of 8 sigma, even with an effective integration time of more than 40 h. However, the average profile obtained by combining data from several epochs, at a dispersion measure of 15.44 pc cm(-3), is found to be consistent with that from the earlier detection of this pulsar at a confidence level of 99.2 per cent. We present this consistency between the two profiles as evidence that J1732-3131 is a faint radio pulsar with an average flux density of 200-400 mJy at 34 MHz. Despite the extremely bright sky background at such low frequencies, the detection sensitivity of our deep searches is generally comparable to that of higher frequency searches for these pulsars, when scaled using reasonable assumptions about the underlying pulsar spectrum. We provide details of our deep searches, and put stringent upper limits on the decametre-wavelength flux densities of several radio-quiet gamma-ray pulsars.

Long way to go: how outflows from large galaxies propagate through the hot halo gas

Using hydrodynamic simulations, we study the mass-loss due to supernova-driven outflows from Milky Way type disc galaxies, paying particular attention to the effect of the extended hot halo gas. We find that the total mass-loss at inner radii scales roughly linearly with total mass of stars formed, and that the mass loading factor at the virial radius can be several times its value at inner radii because of the swept up hot halo gas. The temperature distribution of the outflowing material in the inner region (similar to 10 kpc) is bimodal in nature, peaking at 10(5) K and 10(6.5) K, responsible for optical and X-ray emission, respectively. The contribution of cold/warm gas with temperature <= 10(5.5) K to the outflow rate within 10 kpc is approximate to 0.3-0.5. The warm mass loading factor, eta(3e5) (T <= 3 x 10(5) K) is related to the mass loading factor at the virial radius (eta(v)) as eta(v) approximate to 25 eta(3e5) (SFR/M-circle dot yr(-1))(-0.15) for a baryon fraction of 0.1 and a starburst period of 50 Myr. We also discuss the effect of multiple bursts that are separated by both short and long periods. The outflow speed at the virial radius is close to the sound speed in the hot halo, less than or similar to 200 km s(-1). We identify two `sequences' of outflowing cold gas at small scales: a fast (approximate to 500 km s(-1)) sequence, driven by the unshocked free-wind; and a slow sequence (approximate to +/- 100 km s(-1)) at the conical interface of the superwind and the hot halo.

Formation of Prestellar Cores via Non-Isothermal Gas Fragmentation

Sheet-like clouds are common in turbulent gas and perhaps form via collisions between turbulent gas flows. Having examined the evolution of an isothermal shocked slab in an earlier contribution, in this work we follow the evolution of a sheet-like cloud confined by (thermal) pressure and gas in it is allowed to cool. The extant purpose of this endeavour is to study the early phases of core-formation. The observed evolution of this cloud supports the conjecture that molecular clouds themselves are three-phase media (comprising viz. a stable cold and warm medium, and a third thermally unstable medium), though it appears, clouds may evolve in this manner irrespective of whether they are gravitationally bound. We report, this sheet fragments initially due to the growth of the thermal instability (TI) and some fragments are elongated, filament-like. Subsequently, relatively large fragments become gravitationally unstable and sub-fragment into smaller cores. The formation of cores appears to be a three stage process: first, growth of the TI leads to rapid fragmentation of the slab; second, relatively small fragments acquire mass via gas-accretion and/or merger and third, sufficiently massive fragments become susceptible to the gravitational instability and sub-fragment to form smaller cores. We investigate typical properties of clumps (and smaller cores) resulting from this fragmentation process. Findings of this work support the suggestion that the weak velocity field usually observed in dense clumps and smaller cores is likely seeded by the growth of dynamic instabilities. Simulations were performed using the smooth particle hydrodynamics algorithm.

A Suzaku view of IGR J16393-4643

The pulsar IGR J16393-4643 belongs to a class of highly absorbed supergiant high-mass X-ray binaries (HMXBs), characterized by a very high column density of absorbing matter. We present the results of simultaneous broad-band pulsation and spectrum analysis from a 44-ks Suzaku observation of the source. The orbital intensity profile created with the Swift Burst Alert Telescope (Swift-BAT) light curve shows an indication of IGR J16393-4643 being an eclipsing system with a short eclipse semi-angle theta(E) similar to 17 degrees. For a supergiant companion star with a 20-R-circle dot radius, this implies an inclination of the orbital plane in the range 39 degrees-57 degrees, whereas for a main-sequence B star as the companion with a 10-R-circle dot radius, the inclination of the orbital plane is in the range 60 degrees-77 degrees. Pulse profiles created for different energy bands have complex morphology, which shows some energy dependence and increases in pulse fraction with energy. We have also investigated broad-band spectral characteristics, phase-averaged spectra and resolving the pulse phase into peak and trough phases. The phase-averaged spectrum has a very high N-H(similar to 3 x 10(23) cm(-2)) and is described by a power law (Gamma similar to 0.9) with a high-energy cut-off above 20 keV. We find a change in the spectral index in the peak and trough phases, implying an underlying change in the source spectrum.

The cold mode: a phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a non-linear phenomenological model for the evolution of density perturbations in the intracluster medium (ICM) due to local thermal instability and gravity. We have analysed and extended a model for the evolution of an overdense blob in the ICM. We find two regimes in which the overdense blobs can cool to thermally stable low temperatures. One for large t(cool)/t(ff) (t(cool) is the cooling time and t(ff) is the free-fall time), where a large initial overdensity is required for thermal runaway to occur; this is the regime which was previously analysed in detail. We discover a second regime for t(cool)/t(ff) less than or similar to 1 (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if t(cool)/t(ff) less than or similar to 10.

Stability of Filaments in Star-Forming Clouds and the Formation of Prestellar Cores in Them

The exact process(es) that generate(s) dense filaments which then form prestellar cores within them is unclear. Here we study the formation of a dense filament using a relatively simple set-up of a pressure-confined, uniform-density cylinder. We examine if its propensity to form a dense filament and further, to the formation of prestellar cores along this filament, bears on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised by the approximate equality between the mass line-density for this volume and its maximum value). Though self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then fragments along its length due to the growth of a Jeans-like instability to form prestellar cores. We suggest dense filaments in typical star-forming clouds classified as gravitationally super-critical under the assumption of: (i) isothermality when in fact, they are not, and (ii) extended radial profiles as against pressure-truncated, that significantly over-estimates their mass line-density, are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent surveys of nearby star-forming regions.

Role of gas in supporting grand spiral structure

The density wave theory for the grand-design two-armed spiral pattern in galaxies is successful in explaining several observed features. However, the long-term persistence of this spiral structure is a serious problem since the group transport would destroy it within about a billion years as shown in a classic paper by Toomre. In this paper, we include the low-velocity dispersion component, namely gas, on an equal footing with stars in the formulation of the density wave theory, and obtain the dispersion relation for this coupled system. We show that the inclusion of gas makes the group transport slower by a factor of few, thus allowing the pattern to persist longer - for several billion years. Though still less than the Hubble time, this helps in making the spiral structure more long-lived. Further we show that addition of gas is essential to get a stable wave for the observed pattern speed for the Galaxy, which otherwise is not possible for a one-component stellar disc.

Narrow escape: how ionizing photons escape from disc galaxies

In this paper we calculate the escape fraction (f(esc)) of ionizing photons from starburst galaxies. Using 2D axisymmetric hydrodynamic simulations, we study superbubbles created by overlapping supernovae in OB associations. We calculate the escape fraction of ionizing photons from the centre of the disc along different angles through the superbubble and the gas disc. After convolving with the luminosity function of OB associations, we show that the ionizing photons escape within a cone of similar to 40 degrees, consistent with observations of nearby galaxies. The evolution of the escape fraction with time shows that it falls initially as cold gas is accumulated in a dense shell. After the shell crosses a few scaleheights and fragments, the escape fraction through the polar regions rises again. The angle-averaged escape fraction cannot exceed similar to1 - cos (1 rad)] = 0.5 from geometrical considerations (using the emission cone opening angle). We calculate the dependence of the time-and angle-averaged escape fraction on the mid-plane disc gas density (in the range n(0) = 0.15-50 cm(-3)) and the disc scaleheight (between z(0) = 10 and 600 pc). We find that the escape fraction is related to the disc parameters (the mid-plane disc density and scaleheight) roughly so that f(esc)(alpha)n(0)(2)z(0)(3) (with alpha approximate to 2.2) is a constant. For discs with a given warm neutral medium temperature, massive discs have lower escape fraction than low-mass galaxies. For Milky Way ISM parameters, we find f(esc) similar to 5 per cent, and it increases to approximate to 10 per cent for a galaxy 10 times less massive. We discuss the possible effects of clumpiness of the ISM on the estimate of the escape fraction and the implications of our results for the reionization of the Universe.

Variations in the cyclotron resonant scattering features during 2011 outburst of 4U 0115+63

We study the variations in the Cyclotron Resonant Scattering Feature (CRSF) during 2011 outburst of the high mass X-ray binary 4U 0115+63 using observations performed with Suzaku, RXTE, Swift and INTEGRAL satellites. The wide-band spectral data with low-energy coverage allowed us to characterize the broad-band continuum and detect the CRSFs. We find that the broad-band continuum is adequately described by a combination of a low temperature (kT similar to 0.8 keV) blackbody and a power law with high energy cutoff (E-cut similar to 5.4 keV) without the need for a broad Gaussian at similar to 10 keV as used in some earlier studies. Though winds from the companion can affect the emission from the neutron star at low energies (<3 keV), the blackbody component shows a significant presence in our continuum model. We report evidence for the possible presence of two independent sets of CRSFs with fundamentals at similar to 11 and similar to 15 keV. These two sets of CRSFs could arise from spatially distinct emitting regions. We also find evidence for variations in the line equivalent widths, with the 11 keV CRSF weakening and the 15 keV line strengthening with decreasing luminosity. Finally, we propose that the reason for the earlier observed anticorrelation of line energy with luminosity could be due to modelling of these two independent line sets (similar to 11 and similar to 15 keV) as a single CRSF.

GRMHD formulation of highly super-Chandrasekhar rotating magnetized white dwarfs: stable configurations of non-spherical white dwarfs

Here we extend the exploration of significantly super-Chandrasekhar magnetized white dwarfs by numerically computing axisymmetric stationary equilibria of differentially rotating magnetized polytropic compact stars in general relativity (GR), within the ideal magnetohydrodynamic regime. We use a general relativistic magnetohydrodynamic (GRMHD) framework that describes rotating and magnetized axisymmetric white dwarfs, choosing appropriate rotation laws and magnetic field profiles (toroidal and poloidal). The numerical procedure for finding solutions in this framework uses the 3 + 1 formalism of numerical relativity, implemented in the open source XNS code. We construct equilibrium sequences by varying different physical quantities in turn, and highlight the plausible existence of super-Chandrasekhar white dwarfs, with masses in the range of 2-3 solar mass, with central (deep interior) magnetic fields of the order of 10(14) G and differential rotation with surface time periods of about 1-10 s. We note that such white dwarfs are candidates for the progenitors of peculiar, overluminous Type Ia supernovae, to which observational evidence ascribes mass in the range 2.1-2.8 solar mass. We also present some interesting results related to the structure of such white dwarfs, especially the existence of polar hollows in special cases.

Multiwavelength features of Fermi bubbles as signatures of a Galactic wind

Using hydrodynamical simulations, we show for the first time that an episode of star formation in the centre of the Milky Way, with a star formation rate (SFR) similar to 0.5 M-circle dot yr(-1) for similar to 30 Myr, can produce bubbles that resemble the Fermi bubbles (FBs), when viewed from the solar position. The morphology, extent and multiwavelength observations of FBs, especially X-rays, constrain various physical parameters such as SFR, age, and the circumgalactic medium (CGM) density. We show that the interaction of the CGM with the Galactic wind driven by star formation in the central region can explain the observed surface brightness and morphological features of X-rays associated with the FBs. Furthermore, assuming that cosmic ray electrons are accelerated in situ by shocks and/or turbulence, the brightness and morphology of gamma-ray emission and the microwave haze can be explained. The kinematics of the cold and warm clumps in our model also matches with recent observations of absorption lines through the bubbles.

Effect of dark matter halo on global spiral modes in galaxies

Low surface brightness (LSB) galaxies form a major class of galaxies, and are characterized by low disc surface density and low star formation rate. These are known to be dominated by dark matter halo from the innermost regions. Here, we study the role of the dark matter halo on the grand-design, m = 2, spiral modes in a galactic disc by carrying out a global mode analysis in the WKB approximation. The Bohr-Sommerfeld quantization rule is used to determine how many discrete global spiral modes are permitted. First, a typical superthin, LSB galaxy UGC 7321 is studied by taking only the galactic disc, modelled as a fluid; and then the disc embedded in a dark matter halo. We find that both cases permit the existence of global spiral modes. This is in contrast to earlier results where the inclusion of dark matter halo was shown to nearly fully suppress local, swing-amplified spiral features. Although technically global modes are permitted in the fluid model as shown here, we argue that due to lack of tidal interactions, these are not triggered in LSB galaxies. For comparison, we carried out a similar analysis for the Galaxy, for which the dark matter halo does not dominate in the inner regions. We show that here too the dark matter halo has little effect, hence the disc embedded in a halo is also able to support global modes. The derived pattern speed of the global mode agrees fairly well with the observed value for the Galaxy.