Precision pulsar timing with the ORT and the GMRT and its applications in pulsar astrophysics

Radio pulsars show remarkable clock-like stability, which make them useful astronomy tools in experiments to test equation of state of neutron stars and detecting gravitational waves using pulsar timing techniques. A brief review of relevant astrophysical experiments is provided in this paper highlighting the current state-of-the-art of these experiments. A program to monitor frequently glitching pulsars with Indian radio telescopes using high cadence observations is presented, with illustrations of glitches detected in this program, including the largest ever glitch in PSR B0531+21. An Indian initiative to discover sub-\(\mu \)Hz gravitational waves, called Indian Pulsar Timing Array (InPTA), is also described briefly, where time-of-arrival uncertainties and post-fit residuals of the order of \(\mu \)s are already achievable, comparable to other international pulsar timing array experiments. While timing the glitches and their recoveries are likely to provide constraints on the structure of neutron stars, InPTA will provide upper limits on sub-\(\mu \)Hz gravitational waves apart from auxiliary pulsar science. Future directions for these experiments are outlined.     

Spectral Properties of Black Holes in Gamma Rays

Black holes are the most compact objects in the universe. Therefore, matter accreting onto them is likely to radiate photons of energy comparable to very high gravitational potential energy. We discuss the nature of the emitted radiation in X-rays and gamma-rays from black hole candidates. We present theoretical solutions, which comprise both Keplerian and sub-Keplerian components and suggest that shocks in accretion and outflows may play a major role in producing these spectra.     

Estimation of Sea Ice Freeboard from SARAL/AltiKa Data

An attempt has been made to derive sea ice freeboard from Ka-band Altimeter (SARAL/AltiKa) over Arctic region for 15 March–15 April 2013 (spring) and 15 September–15 October 2013 (autumn). A waveform template matching technique is employed for classification of leads and floe pixels. The estimated sea ice freeboards were found in close agreement with “Operation IceBridge quick look” freeboards (RMSD = 0.30 m). The differences between the two freeboards were largely due to snow layer over sea ice (R = 0.8). The estimated freeboards were of the order of 0.08–0.15 m during the two seasons.     

Turbulent unsteady flow profiles over an adverse slope

When an unsteady free surface flow encounters an adverse slope, it results in a decelerating flow up the adverse slope. The time dependent turbulent flow is treated here by appropriately reducing the two-dimensional Reynolds averaged Navier-Stokes equation along with the equation of continuity considering turbulence closure. With suitable choice of parameters, the resulting differential equations are numerically solved to compute free surface and streamwise velocity profiles with time. It is found that initially the advancing free surface is convex upwards for a short time, followed by a jump of the free surface with a negative streamwise velocity that is a backwater rolling breaker due to deceleration of flow. At later time, however, the velocity becomes positive, that is, the breakers roll forward. This dual feature of motion, that is a surge followed by rolling breakers, is repeated for sometime before the jumps stop. The theoretical analysis presented here is motivated by tidal bores propagating upstream in an estuarine river.     

Role of interannual equatorial forcing on the subsurface temperature dipole in the Bay of Bengal during IOD and ENSO events

Role of equatorial forcing on the thermocline variability in the Bay of Bengal (BoB) during positive and negative phases of the Indian Ocean Dipole (IOD) and El Niño Southern Oscillation (ENSO) was investigated using the Regional Ocean Modeling System (ROMS) simulations during 1988 to 2015. Two numerical experiments were carried out for (i) the Indian Ocean Model (IOM) with interannual open boundary conditions and (ii) the BoB Model (BoBM) with climatological boundary conditions. The first mode of Sea Surface Height Anomalies (SSHA) variability showed a west-east dipole nature in both IOM and altimetry observations around 11°N, which was absent in the BoBM. The vertical section of temperature along the same latitude showed a sharp subsurface temperature dipole with a core at ~ 100 m depth. The positive (negative) subsurface temperature anomalies were observed over the whole northeastern BoB during NIOD (PIOD) and LN (EN) composites due to stronger (weaker) second downwelling Kelvin Waves. During the negative phases of IOD and ENSO, the cyclonic eddy on the southwestern BoB strengthened due to intensified southward coastal current along the western BoB and local wind stress. The subsurface temperature dipole was at its peak during October–December (OND) with 1-month lag from IOD and was evident from the Argo observations and other reanalysis datasets as well. A new BoB dipole index (BDI) was defined as the normalized difference of 100-m temperature anomaly and found to be closely related to the frequency of cyclones and the surface chlorophyll-a concentration in the BoB.

     

Smoothed particle hydrodynamic simulations of viscous accretion discs around black holes

Viscous Keplerian discs become sub-Keplerian close to a black hole since they pass through sonic points before entering into it. We study the time evolution of polytropic viscous accretion discs (both in one- and two-dimensional flows) using smoothed particle hydrodynamics. We discover that for a large region of the parameter space spanned by energy, angular momentum and polytropic index, when the flow viscosity parameter is less than a critical value, standing shock waves are formed. If the viscosity is very high then the shock wave disappears. In the intermediate viscosity, the disc oscillates very significantly in the viscous time-scale. Our simulations indicate that these centrifugally supported high density regions close to a black hole play an active role in the flow dynamics, and consequently, the radiation dynamics.     

Seismic fragility of base-isolated water storage tanks under non-stationary earthquakes

Seismic fragility curves for fixed-base and base-isolated liquid storage tanks are developed under non-stationary earthquakes, and their seismic performance are compared. The correlation between different earthquake intensity measure (IM) parameters and peak response quantities of the base-isolated liquid storage tanks are investigated. The failure criteria are chosen based on (1) the elastic buckling strength of the tank wall, which is defined in terms of critical base shear and critical overturning moment, and (2) in terms of the critical isolation displacement. The uncertainty involved is considered in the earthquake characteristics. Non-stationary earthquake ground motions are generated using Monte Carlo (MC) simulation. Influence of the isolator characteristic parameters and modeling approaches on the seismic fragility of the base-isolated liquid storage tanks is also investigated. Peak ground acceleration is found to be the well correlated IM parameter with the peak response quantities of the base-isolated liquid storage tanks. Substantial decrease in the seismic fragility of the base-isolated liquid storage tanks is observed as compared to the fixed-base tanks. Significant influence of the isolator characteristic parameters on the seismic fragility of the base-isolated liquid storage tanks are reported in the present study.