We introduce pinta, a pipeline for reducing the uGMRT raw pulsar timing data developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. Furthermore, the results of a calibration experiment carried out to determine the relative time offsets between different back-end modes and the correct interpretation of the observatory frequency settings at the uGMRT, which are crucial for performing precision pulsar timing, are also discussed.

The observed values of the time-derivatives of the spin or orbital frequency of pulsars are affected by their dynamical properties. We derive thorough analytical expressions for such dynamical contributions in terms of the Galactic coordinates, the proper motion, the pulsar distance and the radial velocity. We see that although for most of the cases, the effects of the dynamical terms in the second-derivative of frequencies or parameters that depend on such second derivatives, e.g., braking index, are negligible, there might be unique pulsars for which the effects of the dynamical terms would be large enough. We emphasize the fact that our expressions provide more accurate results than pre-existing approximate ones that exclude some of the terms.

XTE J1810-197 (PSR J1809-1943) was the first ever magnetar which was found to emit transient radio emission. It has recently undergone another radio and high-energy outburst. This is only the second radio outburst that has been observed from this source. We observed J1810-197 soon after its recent radio outburst at low radio frequencies using the Giant Metrewave Radio Telescope. We present the 650 MHz flux density evolution of the source in the early phases of the outburst, and its radio spectrum down to frequencies as low as 300 MHz. The magnetar also exhibits radio emission in the form of strong, narrow bursts. We show that the bursts have a characteristic intrinsic width of the order of 0.5-0.7 ms, and discuss their properties in the context of giant pulses and giant micropulses from other pulsars. We also show that the bursts exhibit spectral structures which cannot be explained by interstellar propagation effects. These structures might indicate a phenomenological link with the repeating fast radio bursts which also show interesting, more detailed frequency structures. While the spectral structures are particularly noticeable in the early phases of the outburst, these seem to be less prominent as well as less frequent in the later phases, suggesting an evolution of the underlying cause of these spectral structures.

Recycled pulsars are old ($\gtrsim10^{8}$ yr) neutron stars that are descendants from close, interacting stellar systems. In order to understand their evolution and population, we must find and study the largest number possible of recycled pulsars in a way that is as unbiased as possible. In this work, we present the discovery and timing solutions of five recycled pulsars in binary systems (PSRs J0509$+$0856, J0709$+$0458, J0732$+$2314, J0824$+$0028, J2204$+$2700) and one isolated millisecond pulsar (PSR J0154$+$1833). These were found in data from the Arecibo 327-MHz Drift-Scan Pulsar Survey (AO327). All these pulsars have a low dispersion measure (DM) ($\lesssim 45 \, \rm{pc}\, cm^{-3}$), and have a DM-determined distance of $\lesssim$ 3 kpc. Their timing solutions, have data spans ranging from 1 to $\sim$ 7 years, include precise estimates of their spin and astrometric parameters, and for the binaries, precise estimates of their Keplerian binary parameters. Their orbital periods range from about 4 to 815 days and the minimum companion masses (assuming a pulsar mass of 1.4 $\rm{M_{\odot}}$) range from $\sim$ 0.06--1.11 $\rm{M_{\odot}}$. For two of the binaries we detect post-Keplerian parameters; in the case of PSR~J0709$+$0458 we measure the component masses but with a low precision, in the not too distant future the measurement of the rate of advance of periastron and the Shapiro delay will allow very precise mass measurements for this system. Like several other systems found in the AO327 data, PSRs J0509$+$0854, J0709$+$0458 and J0732$+$2314 are now part of the NANOGrav timing array for gravitational wave detection.

For rotation powered pulsars, the rate of change of the spin period is expected to be positive. On the other hand, for a clean binary system (where gravity is the only force acting on the pulsar) the rate of change of the orbital period is expected to be negative. There are, however, some pulsars with the measurements of the rate of change of the spin period as negative quantities. Likewise, there are some pulsars with positive measurements of the rate of change of the orbital period. We investigate these cases by eliminating the external dynamical effects. We also look at those cases where at first, the measured values possess the correct sign, but on subtracting the dynamical effects, the sign changes. We use `GalDynPsr' package to perform these tasks. Moreover, we investigate possible reasons for such anomalies.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

The observed values of the rate of change of the orbital and the spin periods of pulsars are affected by different dynamical effects, for example, the line-of-sight acceleration and the proper motion of the pulsar relative to the sun. We explore these dynamical effects thoroughly and point out the drawbacks of popular methods. We introduce a package, `GalDynPsr', that evaluates different dynamical effects following traditional as well as improved methods based on the model of the Galactic potential provided in a publicly available package called `galpy'. We argue that the improved methods introduced in this paper should be used for pulsars located 1 kpc or farther away from the solar system, especially when precise values of the rate of change of the periods are required, e.g., while placing limits on alternative theories of gravity, calculating the spin-down limit of the continuous gravitational waves emitted from a rotationally deformed neutron star, understanding pulsar `death-line', etc. GalDynPsr is available online and open for contributions.

The Lense-Thirring effect from spinning neutron stars in double neutron star binaries contribute to the periastron advance of the orbit. This extra term involves the moment of inertia of the neutron stars. Moment of inertia, on the other hand, depends on the mass and spin of the neutron star as well as the equation of state of the matter. If at least one member of the double neutron star binary (better the faster one) is a radio pulsar, then accurate timing analysis might lead to the estimation of the contribution of the Lense-Thirring effect to the periastron advance, which will lead to the measurement of the moment of inertia of the pulsar. Combination of the knowledge on the values of the moment of inertia, the mass, and the spin of the pulsar, will give a new constraint on the equation of state. Pulsars in double neutron star binaries are the best for this purpose as short orbits and moderately high eccentricities make the Lense-Thirring effect substantial, whereas tidal effects are negligible (unlike pulsars with main sequence or white-dwarf binaries). The most promising pulsars are PSR J0737-3039A and PSR J1757-1854. The spin-precession of pulsars due to the misalignment between the spin and the orbital angular momentum vectors affect the contribution of the Lense-Thirring effect to the periastron advance. This effect has been explored for both PSR J0737-3039A and PSR J1757-1854, and as the misalignment angles for both of these pulsars are small, the variation in the Lense-Thirring term is not much. However, to extract the Lense-Thirring effect from the observed rate of the periastron advance, more accurate timing solutions including precise proper motion and distance measurements are essential.

In this work, we report the discovery and characterization of PSR J1411+2551, a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our timing observations of the radio pulsar in the system span a period of about 2.5 years. This timing campaign allowed a precise measurement of its spin period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\, s^{-1}$; from these, we derive a characteristic age of $\sim 10\,$Gyr and a surface magnetic field strength of 2.5 $\times 10^{9}$ G. These numbers indicate that this pulsar was mildly recycled by accretion of matter from the progenitor of the companion star. The system has an eccentric ($e\, = \, 0.17$) 2.61 day orbit. This eccentricity allows a highly significant measurement of the rate of advance of periastron, $\dot{\omega} = 0.07686 \pm 0.00046 ^{\circ}~{\rm yr}^{-1}$. Assuming general relativity accurately models the orbital motion, this implies a total system mass M = $2.538 \pm 0.022 M_{\odot}$. The minimum companion mass is $0.92\, M_{\odot}$ and the maximum pulsar mass is $1.62\, M_{\odot}$. The large companion mass and the orbital eccentricity suggest that PSR J1411+2551 is a double neutron star system; the lightest known to date including the DNS merger GW 170817. Furthermore, the relatively low orbital eccentricity and small proper motion limits suggest that the second supernova had a relatively small associated kick; this and the low system mass suggest that it was an ultra-stripped supernova.

The model that fast radio bursts are caused by plunges of asteroids onto neutron stars can explain both repetitive and non-repetitive bursts. If a neutron star passes through an asteroid belt around another star, there would be a series of bursts caused by series of asteroid impacts. Moreover, the neutron star would cross the same belt repetitively if it is in a binary with the star hosting the asteroid belt, leading to a repeated series of bursts. I explore the properties of neutron star binaries which could lead to the only known repetitive fast radio burst so far (FRB121102). In this model, the next two epochs of bursts are expected around 27-February-2017 and 18-December-2017. On the other hand, if the asteroid belt is located around the neutron star itself, then a chance fall of an asteroid from that belt onto the neutron star would lead to a non-repetitive burst. Even a neutron star grazing an asteroid belt can lead to a non-repetitive burst caused by just one asteroid plunge during the grazing. This is possible even when the neutron star is in a binary with the asteroid hosting star, if the belt and the neutron star orbit are non-coplanar.

It is an exceptionally opportune time for Astrophysics when a number of next-generation mega-instruments are poised to observe the universe across the entire electromagnetic spectrum with unprecedented data quality. The Square Kilometre Array (SKA) is undoubtedly one of the major components of this scenario. In particular, the SKA is expected to discover tens of thousands of new neutron stars giving a major fillip to a wide range of scientific investigations. India has a sizeable community of scientists working on different aspects of neutron star physics with immediate access to both the uGMRT (an SKA pathfinder) and the recently launched X-ray observatory Astrosat. The current interests of the community largely centre around studies of - a) the generation of neutron stars and the SNe connection}, b) the neutron star population and evolutionary pathways}, c) the evolution of neutron stars in binaries and the magnetic fields}, d) the neutron star equation of state}, e) the radio pulsar emission mechanism}, and, f) the radio pulsars as probes of gravitational physics}. Most of these studies are the main goals of the SKA first phase, which is likely to be operational in the next four years. This article summarises the science goals of the Indian neutron star community in the SKA era, with significant focus on coordinated efforts among the SKA and other existing/upcoming instruments.

The Square Kilometre Array (SKA), when it becomes functional, is expected to enrich neutron star (NS) catalogues by at least an order of magnitude over their current state. This includes the discovery of new NS objects leading to better sampling of under-represented NS categories, precision measurements of intrinsic properties such as spin period and magnetic field, as also data on related phenomena such as microstructure, nulling, glitching, etc. This will present a unique opportunity to seek answers to interesting and fundamental questions about the extreme physics underlying these exotic objects in the universe. In this paper, we first present a meta-analysis (from a methodological viewpoint) of statistical analyses performed using existing NS data, with a two-fold goal: First, this should bring out how statistical models and methods are shaped and dictated by the science problem being addressed. Second, it is hoped that these analyses will provide useful starting points for deeper analyses involving richer data from SKA whenever it becomes available. We also describe a few other areas of NS science which we believe will benefit from SKA which are of interest to the Indian NS community.

We present Clusterrank, a new algorithm for identifying dispersed astrophysical pulses. Such pulses are commonly detected from Galactic pulsars and rotating radio transients (RRATs), which are neutron stars with sporadic radio emission. More recently, isolated, highly dispersed pulses dubbed fast radio bursts (FRBs) have been identified as the potential signature of an extragalactic cataclysmic radio source distinct from pulsars and RRATs. Clusterrank helped us discover 14 pulsars and 8 RRATs in data from the Arecibo 327 MHz Drift Pulsar Survey (AO327). The new RRATs have DMs in the range $23.5 - 86.6$ pc cm$^{-3}$ and periods in the range $0.172 - 3.901$ s. The new pulsars have DMs in the range $23.6 - 133.3$ pc cm$^{-3}$ and periods in the range $1.249 - 5.012$ s, and include two nullers and a mode-switching object. We estimate an upper limit on the all-sky FRB rate of $10^5$ day$^{-1}$ for bursts with a width of 10 ms and flux density $\gtrsim 83$ mJy. The DMs of all new discoveries are consistent with a Galactic origin. In comparing statistics of the new RRATs with sources from the RRATalog, we find that both sets are drawn from the same period distribution. In contrast, we find that the period distribution of the new pulsars is different from the period distributions of canonical pulsars in the ATNF catalog or pulsars found in AO327 data by a periodicity search. This indicates that Clusterrank is a powerful complement to periodicity searches and uncovers a subset of the pulsar population that has so far been underrepresented in survey results and therefore in Galactic pulsar population models.

To understand the nature of supernovae and neutron star (NS) formation, as well as binary stellar evolution and their interactions, it is important to probe the distribution of NS masses. Until now, all double NS (DNS) systems have been measured to have a mass ratio close to unity (q $\geq$ 0.91). Here we report the measurement of the individual masses of the 4.07-day binary pulsar J0453+1559 from measurements of the rate of advance of periastron and Shapiro delay: The mass of the pulsar is 1.559(5) $M_{\odot}$ and that of its companion is 1.174(4) $M_{\odot}$; q = 0.75. If this companion is also a neutron star (NS), as indicated by the orbital eccentricity of the system (e=0.11), then its mass is the smallest precisely measured for any such object. The pulsar has a spin period of 45.7 ms and a spin derivative of 1.8616(7) x$10^-19$; from these we derive a characteristic age of ~ 4.1 x $10^9$ years and a magnetic field of ~ 2.9 x $10^9$ G,i.e, this pulsar was mildly recycled by accretion of matter from the progenitor of the companion star. This suggests that it was formed with (very approximately) its current mass. Thus NSs form with a wide range of masses, which is important for understanding their formation in supernovae. It is also important for the search for gravitational waves released during a NS-NS merger: it is now evident that we should not assume all DNS systems are symmetric.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

Pulsars in binary systems have been very successful to test the validity of general relativity in the strong field regime. So far, such binaries include neutron star-white dwarf (NS-WD) and neutron star-neutron star (NS-NS) systems. It is commonly believed that a neutron star-black hole (NS-BH) binary will be much superior for this purpose. But in what sense is this true? Does it apply to all possible deviations?

Studies of binary pulsars provide insight into various theories of physics. Detection of such systems is challenging due to the Doppler modulation of the pulsed signal caused by the orbital motion of the pulsar. We investigated the loss of sensitivity in eccentric binary systems for different types of companions. This reduction of sensitivity should be considered in future population synthesis models for binary pulsars. This loss can be recovered partially by employing the `acceleration search' algorithm and even better by using the `acceleration-jerk search' algorithm.

We report initial results from AO327, a drift survey for pulsars with the Arecibo telescope at 327 MHz. The first phase of AO327 will cover the sky at declinations of -1 to 28 degrees, excluding the region within 5 degrees of the Galactic plane, where high scattering and dispersion make low-frequency surveys sub-optimal. We record data from a 57 MHz bandwidth with 1024 channels and 125 us sampling time. The 60 s transit time through the AO327 beam means that the survey is sensitive to very tight relativistic binaries even with no acceleration searches. To date we have detected 44 known pulsars with periods ranging from 3 ms to 2.21 s and discovered 24 new pulsars. The new discoveries include three millisecond pulsars, three objects with periods of a few tens of milliseconds typical of young as well as mildly recycled pulsars, a nuller, and a rotating radio transient. Five of the new discoveries are in binary systems. The second phase of AO327 will cover the sky at declinations of 28 to 38 degrees. We compare the sensitivity and search volume of AO327 to the Green Bank North Celestial Cap survey and the GBT350 drift survey, both of which operate at 350 MHz.

Luminosity is an intrinsic property of radio pulsars related to the properties of the magnetospheric plasma and the beam geometry, and inversely proportional to the observing frequency. In traditional models, luminosity has been considered as a function of the spin parameters of pulsars. On the other hand, parameter independent models like power law and lognormal have been also used to fit the observed luminosities. Some of the older studies on pulsar luminosities neglected observational biases, but all of the recent studies tried to model observational effects as accurately as possible. Luminosities of pulsars in globular clusters and in the Galactic disk have been studied separately. Older studies concluded that these two categories of pulsars have different luminosity distributions, but the most recent study concluded that those are the same. This article reviews all significant works on pulsar luminosities and discusses open questions.

By generalizing earlier work of Johnston & Kulkarni, we present a detailed description of the reduction in the signal-to-noise ratio for observations of binary pulsars. We present analytical expressions, and provide software, to calculate the sensitivity reduction for orbits of arbitrary eccentricity. We find that this reduction can be quite significant, especially in the case of a massive companion like another neutron star or a black hole. On the other hand, the reduction is less for highly eccentric orbits. We also demonstrate that this loss of sensitivity can be recovered by employing "acceleration search" or "acceleration-jerk search" algorithms.

Studies of the Galactic population of radio pulsars have shown that their luminosity distribution appears to be log-normal in form. We investigate some of the consequences that occur when one applies this functional form to populations of pulsars in globular clusters. We use Bayesian methods to explore constraints on the mean and standard deviation of the luminosity function, as well as the total number of pulsars, given an observed sample of pulsars down to some limiting flux density, accounting for measurements of flux densities of individual pulsars as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 5, 47 Tucanae and M 28, and demonstrate, under reasonable assumptions, that the number of potentially observable pulsars should be within 95% credible intervals of $147^{+112}_{-65}$, $83^{+54}_{-35}$ and $100^{+91}_{-52}$, respectively. Beaming considerations would increase the true population size by approximately a factor of two. Using non-informative priors, however, the constraints are not tight due to the paucity and quality of flux density measurements. Future cluster pulsar discoveries and improved flux density measurements would allow this method to be used to more accurately constrain the luminosity function, and to compare the luminosity function between different clusters.

The luminosity distribution of Galactic radio pulsars is believed to be log-normal in form. Applying this functional form to populations of pulsars in globular clusters, we employ Bayesian methods to explore constraints on the mean and standard deviation of the function, as well as the total number of pulsars in the cluster. Our analysis is based on an observed number of pulsars down to some limiting flux density, measurements of flux densities of individual pulsars, as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 5 and demonstrate, under reasonable assumptions, that the number of potentially observable pulsars is in a 95.45% credible interval of 133$^{+101}_{-58}$. Beaming considerations would increase the true population size by approximately a factor of two.

As neutron star - black hole binaries are expected to be discovered through future pulsar surveys using upcoming facilities, it is necessary to understand various observable properties of such systems. In the present work, we study the advance of the periastron of such binaries under the post-Newtonian formalism over a wide range of parameters. We find that the first and second order post-Newtonian effects and the leading order spin-orbit coupling effects are significant for such binaries but higher order effects can be neglected.

A number of different classes of potentially extra-terrestrial bursts of radio emission have been observed in surveys with the Parkes 64m radio telescope, including "Rotating Radio Transients", the "Lorimer burst" and "perytons". Rotating Radio Transients are radio pulsars which are best detectable in single-pulse searches. The Lorimer burst is a highly dispersed isolated radio burst with properties suggestive of extragalactic origin. Perytons share the frequency-swept nature of the Rotating Radio Transients and Lorimer burst, but unlike these events appear in all thirteen beams of the Parkes Multibeam receiver and are probably a form of peculiar radio frequency interference. In order to constrain these and other radio source populations further, we searched the archival Parkes Multibeam Pulsar Survey data for events similar to any of these. We did not find any new Rotating Radio Transients or bursts like the Lorimer burst. We did, however, discover four peryton-like events. Similar to the perytons, these four bursts are highly dispersed, detected in all thirteen beams of the Parkes multibeam receiver, and have pulse widths between 20--30 ms. Unlike perytons, these bursts are not associated with atmospheric events like rain or lightning. These facts may indicate that lightning was not responsible for the peryton phenomenon. Moreover, the lack of highly dispersed celestial signals is the evidence that the Lorimer burst is unlikely to belong to a cosmological source population.

Using Monte Carlo simulations, we model the luminosity distribution of recycled pulsars in globular clusters as the brighter, observable part of an intrinsic distribution and find that the observed luminosities can be reproduced using either log-normal or power-law distributions as the underlying luminosity function. For both distributions, a wide range of model parameters provide an acceptable match to the observed sample, with the log-normal function providing statistically better agreement in general than the power-law models. Moreover, the power-law models predict a parent population size that is a factor of between two and ten times higher than for the log-normal models. We note that the log-normal luminosity distribution found for the normal pulsar population by Faucher-Gigu\`ere and Kaspi is consistent with the observed luminosities of globular cluster pulsars. For Terzan~5, our simulations show that the sample of detectable radio pulsars, and the diffuse radio flux measurement, can be explained using the log-normal luminosity law with a parent population of $\sim 150$ pulsars. Measurements of diffuse gamma-ray fluxes for several clusters can be explained by both power-law and log-normal models, with the log-normal distributions again providing a better match in general. In contrast to previous studies, we do not find any strong evidence for a correlation between the number of pulsars inferred in globular clusters and globular cluster parameters including metallicity and stellar encounter rate.

Millisecond pulsars are believed to descend from low mass x-ray binaries. Observable parameters of binary millisecond pulsars e.g. mass of the pulsar, mass of the companion, spin period of the pulsar, orbital period, orbital eccentricity etc are used to probe the past accretion history of the millisecond pulsars. But unfortunately in these studies, binding energy of the neutron star is not considered commonly. We show that the effect of the binding energy is significant in the estimation of the amount of mass accretion and thus should be incorporated in models for binary evolutions. Moreover, different Equations of State for dense matter give different values for the accreted mass for the same amount of increase in the gravitational mass of the neutron star implying the need of constraining dense matter Equations of State even to understand the spin-up procedure properly.

Orbital parameters of binary radio pulsars reveal the history of the pulsars' formation and evolution including dynamic interactions with other objects. Advanced technology has enabled us to determine these orbital parameters accurately in most of the cases. Determination of post-Keplerian parameters of double neutron star binaries (especially of the double pulsar) provide clean tests of GTR and in the future may lead us to constrain the dense matter EoS. For binary pulsars with MS or WD companions, knowledge about the values of the orbital parameters as well as of the spin periods and the masses of the pulsars and the companions might be useful to understand the evolutionary history of the systems. As accreting neutron star binaries lead to orbit circularization due to the tidal coupling during accretion, their descendants i.e. binary MSPs are expected to be in circular orbits. On the other hand, dense stellar environments inside globular clusters (GCs) cause different types of interactions of single stars with pulsar binaries. These interactions can impart high eccentricities to the pulsar binaries. So it is quite common to get eccentric millisecond pulsar binaries in GCs and we find that "fly-by" causes intermediate values of eccentricities while "exchange" or "merger" causes high values of eccentricities. We also show that "ionization" is not much effective in the present stage of GCs. Even in the absence of such kinds of stellar interactions, a millisecond pulsar can have an eccentric orbit as a result of Kozai resonance if the pulsar binary is a member of a hierarchical triple system. PSR J1903+0327 is the only one eccentric millisecond pulsar binary in the galactic disk where stellar interactions are negligible. The possibility of this system to be a member of a hierarchical triple system or past association of a GC have been studied and found to be less likely.

We study luminosities of millisecond pulsars in globular clusters by fitting the observed luminosity distribution with single and double power laws. We use simulations to model the observed distribution as the brighter part of some parent distribution for Terzan 5 and try to find a model which simultaneously agrees with the observed diffuse radio flux, total predicted number of pulsars and observed luminosity distribution. We find that wide ranges of parameters for log-normal and power-law distributions give such good models. No clear difference between the luminosity distributions of millisecond pulsars in globular clusters and normal disk pulsars was seen.

The observed distribution of globular cluster binary radio pulsars in the eccentricity versus orbital period plane can be explained as a result of binary-single star interactions. Our numerical and analytical study hints that the highest eccentricity binaries in clusters are likely to be from exchange and/or merger of a single star with a binary component, while the intermediate eccentricity systems are probably results of fly-by interactions.

We investigate the observational signatures associated with one of the proposed formation scenario for the recently discovered highly eccentric binary millisecond pulsar (MSP) PSR J1903+0327 in the galactic plane. The scenario requires that the MSP to be part of a hierarchical triple (HT), consisting of inner and outer binaries, experiencing the Kozai resonance. Numerical modeling of a bound point mass HT, while incorporating the effects due to the quadrupolar interactions between the binary orbits and dominant contributions to the general relativistic periastron precession in the inner binary, reveals that, at the present epoch, the orbital eccentricity of the binary MSP should decrease for reasonable ranges in the HT parameters. The estimated decrements in the orbital eccentricity of the inner binary are few parts in $10^{5}$, substantially higher than the reported accuracies in the estimation of the orbital eccentricity of the binary MSP, while employing various general relativistic timing models for isolated binary pulsars. For wide ranges in the allowed orbital parameters, the estimated rate of change in the eccentricity of the inner binary is orders of magnitude higher than the value recently measured by the pulsar timing analysis. Therefore, we rule out the scenario that the MSP is part of a HT undergoing the Kozai oscillations. The origin of this system in a typical globular cluster is also shown to be less likely than inferred in the discovery paper.

I study stellar structures $i.e.$ the mass, the radius, the moment of inertia and the oblateness parameter at different spin frequencies for strange stars and neutron stars in a comparative manner. I also calculate the values of the radii of the marginally stable orbits and Keplerian orbital frequencies. By equating kHz QPO frequencies to Keplerian orbital frequencies, I find corresponding orbital radii. Knowledge about these parameters might be useful in further modeling of the observed features from LMXBs with advanced and improved future techniques for observations and data analysis.

High sensitivity searches of globular clusters (GC) for radio pulsars by improved pulsar search algorithms and sustained pulsar timing observations have so far yielded some 140 pulsars in more than two dozen GCs. The observed distribution of orbital eccentricity and period of binary radio pulsars in GCs have imprints of the past interaction between single pulsars and binary systems or of binary pulsars and single passing non-compact stars. It is seen that GCs have different groups of pulsars. These may have arisen out of exchange or merger of a component of the binary with the incoming star or a "fly-by" in which the original binary remains intact but undergoes a change of eccentricity and orbital period. We consider the genesis of the distribution of pulsars using analytical and computational tools such as STARLAB, which performs numerical scattering experiments with direct N-body integration. Cluster pulsars with intermediate eccentricities can mostly be accounted for by fly-bys whereas those with high eccentricities are likely to be the result of exchanges and/or mergers of single stars with the binary companion of the pulsar, although there are a few objects which do not easily fit into this description. The corresponding distribution for galactic field pulsars shows notable differences from the GC pulsar orbital period and eccentricity distribution. The long orbital period pulsars in the galactic field with frozen out low eccentricities are largely missing from the globular clusters, and we show that ionization of these systems in GCs cannot alone account for the peculiarities.

We analyze the observed distribution of the orbital eccentricity and period of binary radio pulsars in globular clusters using computational tools to simulate binary-single star interactions. Globular clusters have different groups of pulsars arising from separate interaction scenarios. Intermediate eccentricities of cluster pulsars can be mostly accounted by fly-bys although locally lower stellar densities at pulsar positions may alter the situation. Very high eccentricities are likely to be results of exchanges and/or mergers of single stars with the binary companion of the pulsar.

We study the dependence of the surface magnetic fields of radio pulsars on the choice of Equations of State, pulsar masses and the values of the angle between the magnetic axis and the spin axis of the pulsars within simple dipole model. We show that the values of the surface magnetic field can be even order of magnitude different from its canonical values. This difference will effect any magnetosphere related model to explain observational features of radio pulsars and magnetars. We find a significant difference of the value of the surface magnetic field from the commonly quoted value for the faster member of the double pulsar system, i.e. PSR J0737-3039A as here both the mass of the pulsar and the angle between the magnetic axis and the spin axis are known. Our study reveals the importance of constraining the dense matter Equations of State in pulsar astrophysics as well as hints an alternative way to constrain these by independent determination of the pulsar magnetic field.

High energy density ($\eps$) and temperature (T) links general relativity and hydrodynamics leading to a lower bound for the ratio of shear viscosity ($\eta$) and entropy density ($s$). We get the interesting result that the bound is saturated in the simple model for quark matter that we use for strange stars at the surface for $T \sim 80 MeV$. At this $T$ we have the possibility of cosmic separation of phases. At the surface of the star where the pressure is zero - the density $\eps$ has a fixed value for all stars of various masses with correspondingly varying central energy density $\eps_c$. Inside the star where this density is higher, the ratio of $\eta/s$ is larger and are like the known results found for perturbative QCD. This serves as a check of our calculation. The deconfined quarks at the surface of the strange star at $T = 80 MeV$ seem to constitute the most perfect interacting fluid permitted by nature.

One interesting method of constraining the dense matter Equations of State is to measure the advancement of the periastron of the orbit of a binary radio pulsar (when it belongs to a double neutron star system). There is a great deal of interest on applicability of this procedure to the double pulsar system PSR J0737-3039 (A/B). Although the above method can be applied to PSR A in future within some limitations, for PSR B this method can not be applied. On the other hand, the study of genesis of PSR B might be useful in this connection and its low mass might be an indication that it could be a strange star.

We have studied strange star properties both at zero temperature and at finite temperatures and searched signatures of strange stars in gamma-ray, x-ray and radio astronomy. We have a set of Equations of State (EoS) for strange quark matter (SQM) and solving the TOV equations, we get the structure of strange stars. The maximum mass for a strange star decreases with the increase of temperature, because at high temperatures, the EoS become softer. One important aspect of strange star is that, surface tension depends on the size and structure of the star and is significantly larger than the conventional values. Moment of inertia is another important parameter for compact stars as by comparing theoretical values with observed estimate, it is possible to constrain the dense matter Equation of State. We hope that this approach will help us to decide whether the members of the double pulsar system PSR J0737-3039 are neutron stars or strange stars.

We examine neutron star properties based on a model of dense matter composed of B=1 skyrmions immersed in a mesonic mean field background. The model realizes spontaneous chiral symmetry breaking non-linearly and incorporates scale-breaking of QCD through a dilaton VEV that also affects the mean fields. Quartic self-interactions among the vector mesons are introduced on grounds of naturalness in the corresponding effective field theory. Within a plausible range of the quartic couplings, the model generates neutron star masses and radii that are consistent with a preponderance of observational constraints, including recent ones that point to the existence of relatively massive neutron stars with mass M 1.7 Msun and radius R (12-14) km. If the existence of neutron stars with such dimensions is confirmed, matter at supra-nuclear density is stiffer than extrapolations of most microscopic models suggest.

The effects of temperature on strange stars are studied and it is found that the maximum mass of the star decreases with the increase of temperature since at high temperatures the equations of state become softer. Moreover, if the temperature of a strange star increases, keeping its baryon number fixed, its gravitational mass increases and radius decreases. This leads to a limiting temperature where it turns into a black hole. These features are due to a combined effect of the change of gluon mass and the quark distribution with temperature. We report a new kind of radial oscillations of strange stars driven by what we call chromo-thermal instability. We also discuss the relevance of our findings in the astrophysics of core collapse supernovae and gamma ray bursts.

We propose a new model within the ``Quark-nova'' scenario to interpret the recent observations of early afterglows of long Gamma-Ray Bursts (GRB) with the Swift satellite. This is a three-stage model within the context of a core-collapse supernova. Stage 1 is an accreting (proto-) neutron star leading to a possible delay between the core collapse and the GRB. Stage 2 is an accreting quark-star, generating the prompt GRB. Stage 3, which occurs only if the quark-star collapses to form a black-hole, consists of an accreting black-hole. The jet launched in this accretion process interacts with the ejecta from stage 2, and could generate the flaring activity frequently seen in X-ray afterglows. This model may be able to account for both the energies and the timescales of GRBs, in addition to the newly discovered early X-ray afterglow features.

At finite temperature (T) there is a link with general relativity and hydrodynamics that leads to a lower bound for the ratio of shear viscosity and entropy density (\eta/s). We find that the bound is saturated in the simple model for quark matter that we use for strange stars at T = 80 MeV, at the surface of a strange star. At this T we have the possibility of cosmic separation of phases. We find that, although strongly correlated, the quark matter at the surface of strange stars constitute the most perfect interacting fluid permitted by nature. At the centre of the star, however, the density is higher and conditions are more like the results found for perturbative QCD.

The pions decouple from quarks while their coupling to nucleons $g_{\pi NN}$ remains unchanged at high density. This will have consequences for some recent papers.

Richardson potential is an phenomenological interquark interaction taking care of two aspects of QCD, namely the asymptotic freedom and the confinement. The original potential has a scale parameter having value around 400 MeV and is well tested in hadronic property calculations. This potential was then used in strange star calculation. Strange stars are very compact stars composed of strange quark matter i.e. a very high density strange quark phase consisting of deconfined u, d and s quarks. Here the value of the scale parameter was taken as 100 MeV. The argument was that for a deconfined quark system like a strange star, the scale parameter may have a value quite different from that used in hadronic sector. To remove this discrepancy we introduced two scale parameters in the potential, one for the asymptotic freedom part and the other for the confining part. With suitable values of the parameters, this modified potential has been successfully used in both baryonic property and strange star calculations. The Equation of States obtained with the modified potential are also used to obtain mass-radius relations for the strange stars.

From recent reports on terrestrial heavy ion collision experiments it appears that one may not obtain information about the existence of asymptotic freedom (AF) and chiral symmetry restoration (CSR) for quarks of QCD at high density. This information may still be obtained from compact stars - if they are made up of strange quark matter. Very high gravitational redshift lines (GRL), seen from some compact stars, seem to suggest high ratios of mass and radius (M/R) for them. This is suggestive of strange stars (SS) and can in fact be fitted very well with SQM equation of state deduced with built in AF and CSR. In some other stars broad absorption bands appear at about ~ 0.3 keV and multiples thereof, that may fit in very well with resonance with harmonic compressional breathing mode frequencies of these SS. Emission at these frequencies are also observed in six stars. If these two features of large GRL and BAB were observed together in a single star, it would strengthen the possibility for the existence of SS in nature and would vindicate the current dogma of AF and CSR that we believe in QCD. Recently, in 4U 1700-24, both features appear to be detected, which may well be interpreted as observation of SS - although the group that analyzed the data did not observe this possibility. We predict that if the shifted lines, that has been observed, are from neon with GRL shift z = 0.4 - then the compact object emitting it is a SS of mass 1.2 M_sun and radius 7 km. In addition the fit to the spectrum leaves a residual with broad dips at 0.35 keV and multiples thereof, as in 1E1207-5209 which is again suggestive of SS.

We calculate strange star properties, using large N_c approximation with built-in chiral symmetry restoration (CSM). We used a relativistic Hartree Fock mean field approximation method, using a modified Richardson potential with two scale parameters \Lambda and \Lambda^\prime, to find a new set of equation of states for strange quark matter. We take the effect of temperature (T) on gluon mass, in addition to the usual density dependence, and find that the transition T from hadronic matter to strange matter is 80 MeV. Therefore formation of strange stars may be the only signal for formation of QGP with asymptotic freedom and CSM.

The possible existence of strange stars in the universe will help in the understanding of various properties of quantum chromodynamics, like asymptotic freedom and chiral symmetry restoration, which is otherwise very difficult to prove in laboratory experiments. Strange star properties were calculated using large $N_c$ approximation with built-in chiral symmetry restoration. A relativistic Hartree Fock calculation was performed using the Richardson potential as an interquark interaction. This potential has the asymptotic freedom and a confinement-deconfinement mechanism built into it and the present calculation employs an application of this potential with modified two scale parameters $\Lambda$ and $\Lambda^{\prime}$, to find a new set of equations of state for strange quark matter. The linear confinement string tension from lattice calculations is 350 $MeV$ and the Coulomb -like part has the parameter 100 $MeV$ from deep inelastic scattering experiments. We also consider the effect of temperature, $T$, on gluon mass in a simple way, in addition to the usual density dependence, and find that the transition $T$ from hadronic matter to strange matter is at 80 MeV, close to the 100 MeV estimated in litarature. Therefore formation of strange stars may be the only signal for formation of quark-gluon plasma with asymptotic freedom and chiral symmetry restoration and this may be observable through many processes -such as for example through delayed $\gamma$ ray afterglow.

It appears that there is a genuine shortage of radio pulsars with surface magnetic fields significantly smaller than $\sim 10^8$ Gauss. We propose that the pulsars with very low magnetic fields are actually strange stars locked in a state of minimum free energy and therefore at a limiting value of the magnetic field which can not be lowered by the system spontaneously.

Surface tension (S) is due to the inward force experienced by particles at the surface and usually gravitation does not play an important role in this force. But in compact stars the gravitational force on the particles is very large and S is found to depend not only on the interactions in the strange quark matter, but also on the structure of the star, i.e. on its mass and radius. Indeed, it has been claimed recently that 511 keV photons observed by the space probe INTEGRAL from the galactic bulge may be due to electron-positron annihilation, and their source may be the positron cloud outside of an antiquark star. Such stars, if they exist, may also go a long way towards explaining away the antibaryon deficit of the universe. For that to happen S must be high enough to allow for survival of quark/antiquark stars born in early stages of the formation of the universe. High value of S may also assist explanation of delayed gamma-ray burst after a supernova explosion, as conversion from normal matter to strange matter takes place. The possibility of some implications from formation of surface waves are also discussed.

Chiral symmetry is restored at high density, quarks become nearly massless and pion, the Goldstone of the symmetry breaking decouples from the quarks. What happens at high density is important for finding the density dependence of Strange Quark Matter (SQM), - which in turn is relevant for understanding the structure of compact stars.

Strange stars calculated from a realistic equation of state (EOS), that incorporate chiral symmetry restoration as well as deconfinement at high density show compact objects in the mass radius curve. We compare our calculations of incompressibility for this EOS with that of nuclear matter. One of the nuclear matter EOS has a continuous transition to ud-matter at about five times normal density. Another nuclear matter EOS incorporates density dependent coupling constants. From a look at the consequent velocity of sound, it is found that the transition to ud-matter seems necessary.