Considering the upcoming OWFA, we use simulations of the foregrounds and the $z = 3.35$ HI 21-cm intensity mapping signal to identify the $(k_{\perp},k_{\parallel})$ modes where the expected 21-cm power spectrum $P(k_{\perp},k_{\parallel})$ is substantially larger than the predicted foreground contribution. Only these uncontaminated $k$-modes are used for measuring $P(k_{\perp},k_{\parallel})$ in the "Foreground Avoidance" technique. Though the foregrounds are largely localised within a wedge. we find that the small leakage beyond the wedge surpasses the 21-cm signal across a significant part of the $(k_{\perp},k_{\parallel})$ plane. The extent of foreground leakage is extremely sensitive to the frequency window function used to estimate $P(k_{\perp},k_{\parallel})$. It is possible to reduce the leakage by making the window function narrower, however this comes at the expense of losing a larger fraction of the 21-cm signal. It is necessary to balance these competing effects to identify an optimal window function. Considering a broad class of cosine window functions, we identify a six term window function as optimal for 21-cm power spectrum estimation with OWFA. Considering only the $k$-modes where the expected 21-cm power spectrum exceeds the predicted foregrounds by a factor of $100$ or larger, a $5\,\sigma$ detection of the binned power spectrum is possible in the $k$ ranges $0.18 \leq k \leq 0.3 \, {\rm Mpc}^{-1}$ and $0.18 \le k \le 0.8 \, {\rm Mpc}^{-1}$ with $1,000-2,000$ hours and $10^4$ hours of observation respectively.

Pulsar scintillation allows a glimpse into small-scale plasma structures in the interstellar medium, if we can infer their properties from the scintillation pattern. With Very Long Baseline Interferometry and working in delay-delay rate space, where the contributions of pairs of images to the interference pattern become localized, the scattering geometry and distribution of scattered images on the sky can be determined if a single, highly-anisotropic scattering screen is responsible for the scintillation. However, many pulsars are subject to much more complex scattering environments where this method cannot be used. We present a novel technique to reconstruct the scattered flux of the pulsar and solve for the scattering geometry in these cases by combining interferometric visibilities with cross-correlations of single-station intensities. This takes advantage of the fact that, considering a single image pair in delay-delay rate space, the visibilities are sensitive to the sum of the image angular displacements, while the cross-correlated intensities are sensitive to the difference, so that their combination can be used to localize both images of the pair. We show that this technique is able to reconstruct the published scattering geometry of PSR B0834+06, then apply it to simulations of more complicated scattering systems, where we find that it can distinguish features from different scattering screens even when the presence of multiple screens is not obvious in the Fourier transform of the dynamic spectrum. This technique will allow us to both better understand the distribution of scattering within the interstellar medium and to apply current scintillometry techniques, such as modelling scintillation and constraining the location of pulsar emission, to sources for which a current lack of understanding of the scattering environment precludes the use of these techniques. (abridged)

Simulations play a vital role in testing and validating HI 21-cm power spectrum estimation techniques. Conventional methods use techniques like N-body simulations to simulate the sky signal which is then passed through a model of the instrument. This makes it necessary to simulate the HI distribution in a large cosmological volume, and incorporate both the light-cone effect and the telescope's chromatic response. The computational requirements may be particularly large if one wishes to simulate many realizations of the signal. In this paper we present an analytical method to simulate the HI visibility signal. This is particularly efficient if one wishes to simulate a large number of realizations of the signal. Our method is based on theoretical predictions of the visibility correlation which incorporate both the light-cone effect and the telescope's chromatic response. We have demonstrated this method by applying it to simulate the HI visibility signal for the upcoming Ooty Wide Field Array Phase I.

Foreground removal is the most important step in detecting the large-scale redshifted HI 21-cm signal. Modelling foreground spectra is challenging and is further complicated by the chromatic response of the telescope. We present a multi-frequency angular power spectrum (MAPS) estimator for use in a survey for redshifted HI 21-cm emission from z~3.35, and demonstrate its ability to accurately characterize the foregrounds. This survey will be carried out with the two wide-field interferometer modes of the upgraded Ooty Radio Telescope, called the Ooty Wide Field Array (OWFA), at 326.5 MHz. We have tailored the two-visibility correlation for OWFA to estimate the MAPS and test it with simulated foregrounds. In the process, we describe a software model that encodes the geometry and the details of the telescope, and simulates a realistic model for the bright radio sky. This article presents simulations which include the full chromatic response of the telescope, in addition to the frequency dependence intrinsic to the foregrounds. We find that the visibility correlation MAPS estimator recovers the input angular power spectrum accurately, and that the instrument response to the foregrounds dominates the systematic errors in the recovered foreground power spectra.

The upcoming Ooty Wide Field Array (OWFA) will operate at $326.5 \, {\rm MHz}$ which corresponds to the redshifted 21-cm signal from neutral hydrogen (HI) at z = 3.35. We present two different prescriptions to simulate this signal and calculate the visibilities expected in radio-interferometric observations with OWFA. In the first method we use an input model for the expected 21-cm power spectrum to directly simulate different random realizations of the brightness temperature fluctuations and calculate the visibilities. This method, which models the HI signal entirely as a diffuse radiation, is completely oblivious to the discrete nature of the astrophysical sources which host the HI. While each discrete source subtends an angle that is much smaller than the angular resolution of OWFA, the velocity structure of the HI inside the individual sources is well within reach of OWFA's frequency resolution and this is expected to have an impact on the observed HI signal. The second prescription is based on cosmological N-body simulations. Here we identify each simulation particle with a source that hosts the HI, and we have the freedom to implement any desired line profile for the HI emission from the individual sources. Implementing a simple model for the line profile, we have generated several random realizations of the complex visibilities. Correlations between the visibilities measured at different baselines and channels provides an unique method to quantify the statistical properties of the \HI signal. We have used this to quantify the results of our simulations, and explore the relation between the expected visibility correlations and the underlying HI power spectrum.

One of the scientific objectives of the Ooty Wide Field Array(OWFA) is to observe the redshifted HI emission from z ~ 3.35. Although predictions spell out optimistic outcomes in reasonable integration times, these studies were based purely on analytical assumptions, without accounting for limiting systematics. A software model for OWFA has been developed with a view to understanding the instrument-induced systematics, by describing a complete software model for the instrument. This model has been implemented through a suite of programs, together called Prowess, which has been conceived with a dual role of an emulator as well as observatory data analysis software. The programming philosophy followed in building Prowess enables any user to define an own set of functions and add new functionality. This paper describes a co-ordinate system suitable for OWFA in which the baselines are defined. The foregrounds are simulated from their angular power spectra. The visibilities are then computed from the foregrounds. These visibilities are then used for further processing, such as calibration and power spectrum estimation. The package allows for rich visualisation features in multiple output formats in an interactive fashion, giving the user an intuitive feel for the data. Prowess has been extensively used for numerical predictions of the foregrounds for the OWFA HI experiment.

For radio interferometric arrays with a sufficient number of redundant spacings the multiplicity of measurements of the same sky visibility can be used to determine both the antenna gains as well as the true visibilities. Many of the earlier approaches to this problem focused on linearized versions of the relation between the measured and the true visibilities. Here we propose to use a standard non-linear minimization algorithm to solve for both the antenna gains as well as the true visibilities. We show through simulations done in the context of the ongoing upgrade to the Ooty Radio Telescope that the non-linear minimization algorithm is fast compared to the earlier approaches. Further, unlike the most straightforward linearized approach, which works with the logarithms of the visibilities and the gains, the non-linear minimization algorithm leads to unbiased solutions. Finally we present error estimates for the estimated gains and visibilities. Monte-Carlo simulations establish that the estimator is indeed statistically efficient, achieving the Cramer-Rao bound.

We describe the serendipitous discovery of a radio point source in a 618 MHz image of the supernova remnant(SNR) G76.9+1.0. The SNR has a bipolar structure and the point source is located near a faint bridge of emission joining the two lobes of emission. The point source was also detected in follow-up higher frequency(1170 MHz) observations. The spectral index for the point source obtained from the GMRT observations is alpha = -2.1. The steep spectrum, as well as the location of the point source near the centre of the SNR establish the fact that it is indeed the pulsar J2022+3842 associated with this SNR. Consistent with this, subsequent analysis of archival Chandra X-ray data shows a point source coincident with the radio point source, as well as diffuse extended X-ray emission surrounding the unresolved source. However, no pulsed emission was detected despite deep searches at both 610 MHz and 1160 MHz although pulsed emission has been seen at 2 GHz with the GBT. It appears that the most likely reason for not detecting the pulsed signal at the GMRT is temporal broadening: for the estimated DM towards this SNR, the pulse broadening time could be as large as tens of milliseconds. The diffuse X-ray emission is elongated along the same direction as the bipolar structure seen in the radio. We interpret the radio lobes as having been formed from an equatorial wind. Although direct detection of pulsed signal has not been possible, we show convincingly that sensitive, high-resolution, radio imaging at multiple frequencies is a useful method to search for pulsar candidates.

The OH molecule, found abundantly in the Milky Way, has four transitions at the ground state rotational level(J = 3/2) at cm wavelengths. These are E1 transitions between the F+ and F- hyperfine levels of the Lambda doublet of the J=3/2 state. There are also forbidden M1 transitions between the hyperfine levels within each of the doublet states occuring at frequencies 53.171 MHz and 55.128 MHz. These are extremely weak and hence difficult to detect. However there is a possibility that the level populations giving rise to these lines are inverted under special conditions, in which case it may be possible to detect them through their maser emission. We describe the observational diagnostics for determining when the hyperfine levels are inverted, and identify a region near W44 where these conditions are satisfied. A high-velocity-resolution search for these hyperfine OH lines using the low frequency feeds on four antennas of the GMRT and the new GMRT Software Backend(GSB) was performed on this target near W44. We place a 3-sigma upper limit of ~17.3 Jy (at 1 km/s velocity resolution) for the 55 MHz line from this region. This corresponds to an upper limit of 3 X 10^8 for the amplification of the Galactic synchrotron emission providing the background.