Theoretical determination of HI vertical scale heights in the dwarf galaxies DDO 154, Ho II, IC 2574 and NGC 2366

In this paper, we model dwarf galaxies as a two-component system of gravitationally coupled stars and atomic hydrogen gas in the external force field of a pseudo-isothermal dark matter halo, and numerically obtain the radial distribution of HI vertical scale heights. This is done for a group of four dwarf galaxies (DDO 154, Ho II, IC 2574 and NGC 2366) for which most necessary input parameters are available from observations. The formulation of the equations takes into account the rising rotation curves generally observed in dwarf galaxies. The inclusion of self-gravity of the gas into the model at par with that of the stars results in scale heights that are smaller than what was obtained by previous authors. This is important as the gas scale height is often used for deriving other physical quantities. The inclusion of gas self-gravity is particularly relevant in the case of dwarf galaxies where the gas cannot be considered a minor perturbation to the mass distribution of the stars. We find that three out of four galaxies studied show a flaring of their HI discs with increasing radius, by a factor of a few within several disc scale lengths. The fourth galaxy has a thick HI disc throughout. This flaring arises as a result of the gas velocity dispersion remaining constant or decreasing only slightly while the disc mass distribution declines exponentially as a function of radius.

Non-Gaussianity and direction-dependent systematics in HST Key Project data

Two new statistics, namely Delta(chi 2) and Delta(chi), based on the extreme value theory, were derived by Gupta et al. We use these statistics to study the direction dependence in the HST Key Project data, which provides one of the most precise measurements of the Hubble constant. We also study the non-Gaussianity in this data set using these statistics. Our results for Delta(chi 2) show that the significance of direction-dependent systematics is restricted to well below the 1 sigma confidence limit; however, the presence of non-Gaussian features is subtle. On the other hand, the Delta(chi). statistic, which is more sensitive to direction dependence, shows direction dependence systematics to be at a slightly higher confidence level, and the presence of non-Gaussian features at a level similar to the Delta(chi 2) statistic.

Automatic Compilation of MATLAB Programs for Synergistic Execution on Heterogeneous Processors

MATLAB is an array language, initially popular for rapid prototyping, but is now being increasingly used to develop production code for numerical and scientific applications. Typical MATLAB programs have abundant data parallelism. These programs also have control flow dominated scalar regions that have an impact on the program's execution time. Today's computer systems have tremendous computing power in the form of traditional CPU cores and throughput oriented accelerators such as graphics processing units(GPUs). Thus, an approach that maps the control flow dominated regions to the CPU and the data parallel regions to the GPU can significantly improve program performance. In this paper, we present the design and implementation of MEGHA, a compiler that automatically compiles MATLAB programs to enable synergistic execution on heterogeneous processors. Our solution is fully automated and does not require programmer input for identifying data parallel regions. We propose a set of compiler optimizations tailored for MATLAB. Our compiler identifies data parallel regions of the program and composes them into kernels. The problem of combining statements into kernels is formulated as a constrained graph clustering problem. Heuristics are presented to map identified kernels to either the CPU or GPU so that kernel execution on the CPU and the GPU happens synergistically and the amount of data transfer needed is minimized. In order to ensure required data movement for dependencies across basic blocks, we propose a data flow analysis and edge splitting strategy. Thus our compiler automatically handles composition of kernels, mapping of kernels to CPU and GPU, scheduling and insertion of required data transfer. The proposed compiler was implemented and experimental evaluation using a set of MATLAB benchmarks shows that our approach achieves a geometric mean speedup of 19.8X for data parallel benchmarks over native execution of MATLAB.

The spectral index-luminosity relationship for steep-spectrum cores in extra-galactic radio-sources

The spectral index-luminosity relationship for steep-spectrum cores in galaxies and quasars has been investigated, and it is found that the sample of galaxies supports earlier suggestions of a strong correlation, while there is weak evidence for a similar relationship for the quasars. It is shown that a strong spectral index-luminosity correlation can be used to set an upper limit to the velocities of the radio-emitting material which is expelled from the nucleus in the form of collimated beams or jets having relativistic bulk velocities. The data on cores in galaxies indicate that the Lorentz factors of the radiating material are less than about 2.

Null Dereference Verification via Over-approximated Weakest Pre-conditions Analysis

Null dereferences are a bane of programming in languages such as Java. In this paper we propose a sound, demand-driven, inter-procedurally context-sensitive dataflow analysis technique to verify a given dereference as safe or potentially unsafe. Our analysis uses an abstract lattice of formulas to find a pre-condition at the entry of the program such that a null-dereference can occur only if the initial state of the program satisfies this pre-condition. We use a simplified domain of formulas, abstracting out integer arithmetic, as well as unbounded access paths due to recursive data structures. For the sake of precision we model aliasing relationships explicitly in our abstract lattice, enable strong updates, and use a limited notion of path sensitivity. For the sake of scalability we prune formulas continually as they get propagated, reducing to true conjuncts that are less likely to be useful in validating or invalidating the formula. We have implemented our approach, and present an evaluation of it on a set of ten real Java programs. Our results show that the set of design features we have incorporated enable the analysis to (a) explore long, inter-procedural paths to verify each dereference, with (b) reasonable accuracy, and (c) very quick response time per dereference, making it suitable for use in desktop development environments.

Timing and broad-band spectroscopy of 1A 1118-61 with Suzaku

We present a timing and broad-band pulse-phase-resolved spectral analysis of the transient Be X-ray binary pulsar 1A 1118-61 observed during its outburst in 2009 January using Suzaku observations. The Suzaku observations were made twice, once at the peak of the outburst, and the other 13 d later at its declining phase. Pulse profiles from both observations exhibit strong energy dependence with several peaks at low energies and a single peak above similar to 10 keV. A weak, narrow peak is detected at the main dip of the pulse profiles from both observations in the energy bands below 3 keV, indicating the presence of a phase-dependent soft excess in the source continuum. The broad-band energy spectrum of the pulsar could be fitted well with a partial covering cut-off power-law model and a narrow iron fluorescence line. We also detect a broad cyclotron feature at similar to 50 keV from both observations which is a feature common for accretion-powered pulsars with high magnetic field strength. The pulse-phase-resolved spectral analysis shows an increase in the absorption column density of the partial covering component, as well as variation in the covering fraction at the dips of the pulse profiles, which naturally explains energy dependence of the same. The cyclotron line parameters also show significant variation with pulse phase with an similar to 10 keV variation in the cyclotron line energy and a variation in depth by a factor of 3. This can be explained either as the effect of different viewing angles of the dipole field at different pulse phases, or due to a more complex underlying magnetic field geometry.

Thermal instability and the feedback regulation of hot haloes in clusters, groups and galaxies

We present global multidimensional numerical simulations of the plasma that pervades the dark matter haloes of clusters, groups and massive galaxies (the intracluster medium; ICM). Observations of clusters and groups imply that such haloes are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally motivated heating prescriptions, we show that local thermal instability (TI) can produce a multiphase medium with similar to 104 K cold filaments condensing out of the hot ICM only when the ratio of the TI time-scale in the hot plasma (tTI) to the free-fall time-scale (tff) satisfies tTI/tff? 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of similar to 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling flow is similar to 10-3 for clusters and decreases for lower mass haloes; supernova heating may be energetically sufficient to balance cooling in galactic haloes. We further argue that the ICM self-adjusts so that tTI/tff? 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives tTI/tff? 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminositytemperature () relation. The high-velocity clouds observed in the Galactic halo can be due to local TI producing multiphase gas close to the virial radius if the density of the hot plasma in the Galactic halo is >rsim 10-5 cm-3 at large radii.

Unstable m=1 modes of counter-rotating Keplerian discs

We study the linear m= 1 counter-rotating instability in a two-component, nearly Keplerian disc. Our goal is to understand these slow modes in discs orbiting massive black holes in galactic nuclei. They are of interest not only because they are of large spatial scale and can hence dominate observations but also because they can be growing modes that are readily excited by accretion events. Self-gravity being non-local, the eigenvalue problem results in a pair of coupled integral equations, which we derive for a two-component softened gravity disc. We solve this integral eigenvalue problem numerically for various values of mass fraction in the counter-rotating component. The eigenvalues are in general complex, being real only in the absence of the counter-rotating component, or imaginary when both components have identical surface density profiles. Our main results are as follows: (i) the pattern speed appears to be non-negative, with the growth (or damping) rate being larger for larger values of the pattern speed; (ii) for a given value of the pattern speed, the growth (or damping) rate increases as the mass in the counter-rotating component increases; (iii) the number of nodes of the eigenfunctions decreases with increasing pattern speed and growth rate. Observations of lopsided brightness distributions would then be dominated by modes with the least number of nodes, which also possess the largest pattern speeds and growth rates.

Pulsed radio emission from the Fermi-LAT pulsar J1732-3131: Search and a possible detection at 34.5 MHz

We report our search for and a possible detection of periodic radio pulses at 34.5 MHz from the Fermi Large Area Telescope pulsar J1732-3131. The candidate detection has been possible in only one of the many sessions of observations made with the low-frequency array at Gauribidanur, India, when the otherwise radio weak pulsar may have apparently brightened many folds. The candidate dispersion measure along the sight line, based on the broad periodic profiles from �20min of data, is estimated to be 15.44 ± 0.32 pccc -1. We present the details of our periodic and single-pulse search, and discuss the results and their implications relevant to both, the pulsar and the intervening medium. © 2012 RAS.

Energy-dependent orbital modulation of X-rays and constraints on emission of the jet in Cyg X-3

We study the orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and RXTE. Using the wealth of data presently available and an improved averaging method, we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. We find that above similar to 5?keV the modulation depth decreases with increasing energy, which is consistent with the modulation being caused by both boundfree absorption and Compton scattering in the stellar wind of the donor, with minima corresponding to the highest optical depth, which occurs around the superior conjunction. We find a decrease of the depth below similar to 3?keV, which appears to be due to re-emission of the absorbed continuum by the wind in soft X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet in Cyg X-3, which emits ?-rays detected at energies >0.1?GeV in the soft spectral states, is found to be minor up to similar to 100?keV. This implies the presence of a rather sharp low-energy break in the jet MeV-range spectrum. We also calculate phase-resolved RXTE X-ray spectra and show that the difference between the spectra corresponding to phases around superior and inferior conjunctions can indeed be accounted for by the combined effect of boundfree absorption in an ionized medium and Compton scattering.

Model Pseudoconvex Domains and Bumping

The Levi geometry at weakly pseudoconvex boundary points of domains in C-n, n >= 3, is sufficiently complicated that there are no universal model domains with which to compare a general domain. Good models may be constructed by bumping outward a pseudoconvex, finite- type Omega subset of C-3 in such a way that: (i) pseudoconvexity is preserved, (ii) the (locally) larger domain has a simpler defining function, and (iii) the lowest possible orders of contact of the bumped domain with partial derivative Omega, at the site of the bumping, are realized. When Omega subset of C-n, n >= 3, it is, in general, hard to meet the last two requirements. Such well-controlled bumping is possible when Omega is h-extendible/semiregular. We examine a family of domains in C-3 that is strictly larger than the family of h-extendible/semiregular domains and construct explicit models for these domains by bumping.

Exploiting the Structure of the Constraint Graph for Efficient Points-to Analysis

Points-to analysis is a key compiler analysis. Several memory related optimizations use points-to information to improve their effectiveness. Points-to analysis is performed by building a constraint graph of pointer variables and dynamically updating it to propagate more and more points-to information across its subset edges. So far, the structure of the constraint graph has been only trivially exploited for efficient propagation of information, e.g., in identifying cyclic components or to propagate information in topological order. We perform a careful study of its structure and propose a new inclusion-based flow-insensitive context-sensitive points-to analysis algorithm based on the notion of dominant pointers. We also propose a new kind of pointer-equivalence based on dominant pointers which provides significantly more opportunities for reducing the number of pointers tracked during the analysis. Based on this hitherto unexplored form of pointer-equivalence, we develop a new context-sensitive flow-insensitive points-to analysis algorithm which uses incremental dominator update to efficiently compute points-to information. Using a large suite of programs consisting of SPEC 2000 benchmarks and five large open source programs we show that our points-to analysis is 88% faster than BDD-based Lazy Cycle Detection and 2x faster than Deep Propagation. We argue that our approach of detecting dominator-based pointer-equivalence is a key to improve points-to analysis efficiency.

Why are some galaxy discs extremely thin?

Some low-surface-brightness galaxies are known to have extremely thin stellar discs with the vertical-to-planar axes ratio 0.1 or less, often referred to as superthin galaxies. Although their existence is now known for over three decades, the physical origin of the superthin discs is still not understood. We model the vertical thickness of the stellar disc using our model of a two-component (gravitationally coupled stars and gas) disc embedded in a dark matter halo, for a bulgeless, superthin galaxy UGC 7321 which has a dense, compact halo, and is compare with a typical dwarf irregular galaxy Holmberg II which has a low-density, non-compact halo. We show that while the presence of gas does constrain the stellar disc thickness and hence its axial ratio, it is the compact dark matter halo which plays the decisive role in determining the mean distribution of stars in the vertical direction in low-luminosity bulgeless galaxies like UGC 7321, and causes the stellar disc to be superthin. Thus, the compactness of the dark matter halo significantly affects the disc structure and this could be important for the early evolution of galaxies.

Dynamics of prolate spheroidal mass distributions with varying eccentricity

In this paper we calculate the potential for a prolate spheroidal distribution as in a dark matter halo with a radially varying eccentricity. This is obtained by summing up the shell-by-shell contributions of isodensity surfaces, which are taken to be concentric and with a common polar axis and with an axis ratio that varies with radius. Interestingly, the constancy of potential inside a shell is shown to be a good approximation even when the isodensity contours are dissimilar spheroids, as long as the radial variation in eccentricity is small as seen in realistic systems. We consider three cases where the isodensity contours are more prolate at large radii, or are less prolate or have a constant eccentricity. Other relevant physical quantities like the rotation velocity, the net orbital and vertical frequency due to the halo and an exponential disc of finite thickness embedded in it are obtained. We apply this to the kinematical origin of Galactic warp, and show that a prolate-shaped halo is not conducive to making long-lived warps - contrary to what has been proposed in the literature. The results for a prolate mass distribution with a variable axis ratio obtained are general, and can be applied to other astrophysical systems, such as prolate bars, for a more realistic treatment.

Radiatively inefficient accretion flow simulations with cooling: implications for black hole transients

We study the effects of optically thin radiative cooling on the structure of radiatively inefficient accretion flows (RIAFs). The flow structure is geometrically thick, and independent of the gas density and cooling, if the cooling time is longer than the viscous time-scale (i.e. t(cool) greater than or similar to t(visc)). For higher densities, the gas can cool before it can accrete and forms the standard geometrically thin, optically thick Shakura-Sunyaev disc. For usual cooling processes (such as bremsstrahlung), we expect an inner hot flow and an outer thin disc. For a short cooling time the accretion flow separates into two phases: a radiatively inefficient hot coronal phase and a cold thin disc. We argue that there is an upper limit on the density of the hot corona corresponding to a critical value of t(cool)/t(ff)( similar to 10-100), the ratio of the cooling time and the free-fall time. Based on our simulations, we have developed a model for transients observed in black hole X-ray binaries (XRBs). An XRB in a quiescent hot RIAF state can transition to a cold blackbody-dominated state because of an increase in the mass accretion rate. The transition from a thin disc to a RIAF happens because of mass exhaustion due to accretion; the transition happens when the cooling time becomes longer than the viscous time at inner radii. Since the viscous time-scale for a geometrically thin disc is quite long, the high-soft state is expected to be long-lived. The different time-scales in black hole transients correspond to different physical processes such as viscous evolution, cooling and free fall. Our model captures the overall features of observed state transitions in XRBs.