How would observers differentiate Beacons from pulsars or other exotic sources, in light of likely Beacon observables? Bandwidth, pulse width and frequency may be distinguishing features. Such transients could be evidence of civilizations slightly higher than ourselves on the Kardashev scale.

Context. Global MHD simulations show Kelvin-Helmholtz (KH) instabilities at the contact surface of two merging neutron stars. That region has been identified as the site of efficient amplification of magnetic fields. However, these global simulations, due to numerical limitations, were unable to determine the saturation level of the field strength, and thus the possible back-reaction of the magnetic field onto the flow. Aims. We investigate the amplification of initially weak fields in KH unstable shear flows, and the back-reaction of the field onto the flow. Methods. We use a high-resolution ideal MHD code to perform 2D and 3D local simulations of shear flows. Results. In 2D, the magnetic field is amplified in less than 0.01ms until it reaches locally equipartition with the kinetic energy. Subsequently, it saturates due to resistive instabilities that disrupt the KH vortex and decelerate the shear flow on a secular time scale. We determine scaling laws of the field amplification with the initial field strength and the grid resolution. In 3D, this hydromagnetic mechanism may be dominated by purely hydrodynamic instabilities limiting the amplification. We find maximum magnetic fields of 10^16 G locally, and r.m.s. maxima within the box of 10^15 G. However, such strong fields exist only for a short period. In the saturated state, the magnetic field is mainly oriented parallel to the shear flow for strong initial fields, while weaker initial fields tend to lead to a more balanced distribution of the field energy. In all models the flow shows small-scale features. The magnetic field is at most in equipartition with the decaying shear flow. (abridged)

Recent observations of the Crab pulsar show no evidence for a spectral break in the infrared regime. It is argued that the observations are consistent with a power-law spectrum in the whole observable infrared - optical range. This is taken as the starting point for an evaluation of how self-consistent incoherent synchrotron models fare in a comparison with observations. Inclusion of synchrotron self-absorption proves important as does the restriction on the observed size of the emission region imposed by the relativistic beaming thought to define the pulse profile. It is shown that the observations can be used to derive two independent constraints on the distance from the neutron star to the emission region; in addition to a direct lower limit, an indirect measure is obtained from an upper limit to the magnetic field strength. Both of these limits indicate that the emission region is located at a distance considerably larger than the light cylinder radius. The implications of this result are discussed and it is emphasized that, in order for standard incoherent synchrotron models to fit inside the light cylinder, rather special physical conditions need to be invoked.

We use the growing data sets of very-metal-poor stars to study the impact of stellar winds of fast rotating massive stars on the chemical enrichment of the early Galaxy. We use an inhomogeneous chemical evolution model for the Galactic halo to predict both the mean trend and scatter of C/O and N/O. In one set of models, we assume that massive stars enrich the interstellar medium during both the stellar wind and supernovae phases. In the second set, we consider that in the earliest phases (Z <10^-8), stars with masses above 40 Msun only enrich the interstellar medium via stellar winds, collapsing directly into black holes. We predict a larger scatter in the C/O and N/O ratios at low metallicities when allowing the more massive fast-rotating stars to contribute to the chemical enrichment only via stellar winds. The latter assumption, combined with the stochasticity in the star formation process in the primordial Galactic halo can explain the wide spread observed in the N/O and C/O ratios in normal very-metal-poor stars. For chemical elements with stellar yields that depend strongly on initial mass (and rotation) such as C, N, and neutron capture elements, within the range of massive stars, a large scatter is expected once the stochastic enrichment of the early interstellar medium is taken into account. We also find that stellar winds of fast rotators mixed with interstellar medium gas are not enough to explain the large CNO enhancements found in most of the carbon-enhanced very-metal-poor stars. In particular, this is the case of the most metal-poor star known to date, HE 1327-2326, for which our models predict lower N enhancements than observed when assuming a mixture of stellar winds and interstellar medium. We suggest that these carbon-enhanced very metal-poor stars were formed from almost pure stellar wind material, without dilution with the pristine interstellar medium.

We have re-analyzed the stability of pulse arrival times from pulsars and white dwarfs using several analysis tools for measuring the noise characteristics of sampled time and frequency data. We show that the best terrestrial artificial clocks substantially exceed the performance of astronomical sources as time-keepers in terms of accuracy (as defined by cesium primary frequency standards) and stability. The superiority in stability is demonstrated over timescales up to 2 years. Beyond 2 years there is a deficiency of data for clock/clock comparisons and both terrestrial and astronomical clocks are equally limited by the quality of the time dissemination systems used to make the comparisons. Nonetheless, we show that detailed accuracy evaluations of modern terrestrial clocks imply that these clocks are likely to have a stability better than any astronomical source up to timescales of at least hundreds of years. This is in conflict with many claims in the literature and so we believe it is crucial to rectify this misunderstanding so that there may be a correct appreciation of the relative merits of natural and artificial clocks. The use of natural clocks as tests of physics under the most extreme conditions is entirely appropriate; however, the contention that these natural clocks, particularly white dwarfs, can compete as timekeepers against devices constructed by mankind is shown to be doubtful.

The work presented here examines populations of double compact binary systems and tidally enhanced collapsars. We make use of BINPOP and BINKIN, two components of a recently developed population synthesis package. Results focus on correlations of both binary and spatial evolutionary population characteristics. Pulsar and long duration gamma-ray burst observations are used in concert with our models to draw the conclusions that: double neutron star binaries can merge rapidly on timescales of a few million years (much less than that found for the observed double neutron star population), common envelope evolution within these models is a very important phase in double neutron star formation, and observations of long gamma-ray burst projected distances are more centrally concentrated than our simulated coalescing double neutron star and collapsar Galactic populations. Better agreement is found with dwarf galaxy models although the outcome is strongly linked to the assumed birth radial distribution. The birth rate of the double neutron star population in our models range from 4-160 Myr^-1 and the merger rate ranges from 3-150 Myr^-1. The upper and lower limits of the rates results from including electron capture supernova kicks to neutron stars and decreasing the common envelope efficiency respectively. Our double black hole merger rates suggest that black holes should receive an asymmetric kick at birth.

We review X-ray properties of the Supergiant Fast X-ray Transients following from their observations with INTEGRAL and show that a compact object in these systems is a neutron star with strong magnetic field accreting from the stellar wind of a donor star. We show that presence of a centrifugal barrier at the magnetospheric boundary of the neutron star may be a key to understanding of abrupt short X-ray outbursts of these transients and long intervals of their quiescence.

We discuss new limits on masses and radii of compact stars and we conclude that they can be interpreted as an indication of the existence of two classes of stars: "normal" compact stars and "ultra-compact" stars. We estimate the critical mass at which the first configuration collapses into the second.

(abbrev.) Linear acceleration emission occurs when a charged particle is accelerated parallel to its velocity. We evaluate the spectral and angular distribution of this radiation for several special cases, including constant acceleration (hyperbolic motion) of finite duration. Based on these results, we find the following general properties of the emission from an electron in a linear accelerator that can be characterized by an electric field E acting over a distance L: (i) the spectrum extends to a cut-off photon energy ~ LE^2 MeV, where E is in units of the Schwinger critical field and L in units of the Compton wavelength of the electron. (ii) the total energy emitted by a particle traversing the accelerator is in agreement with the standard Larmor formula (iii) the low frequency spectrum is flat for hyperbolic trajectories, but in general depends on the details of the accelerator. We also show that linear acceleration emission complements curvature radiation in the strongly magnetized pair formation regions in pulsar magnetospheres. It dominates when the length L of the accelerator is less than the formation length of curvature photons, which is given by the ratio of the radius of curvature of the magnetic field lines to the Lorentz factor of the particle. In standard static models of pair creating regions linear acceleration emission is negligible, but it is important in more realistic dynamical models in which the accelerating field fluctuates on a short length-scale.

We investigate the evolution of Type Ib/c supernova (SN Ib/c) progenitors in close binary systems, using new evolutionary models that include the effects of rotation, with initial masses of 12 - 25 Msun for the primary components, and of single helium stars with initial masses of 2.8 - 20 Msun. We find that, despite the impact of tidal interaction on the rotation of primary stars, the amount of angular momentum retained in the core at the presupernova stage in different binary model sequences converge to a value similar to those found in previous single star models. This amount is large enough to produce millisecond pulsars, but too small to produce magnetars or long gamma-ray bursts. We employ the most up-to-date estimate for the Wolf-Rayet mass loss rate, and its implications for SN Ib/c progenitors are discussed in detail. In terms of stellar structure, SN Ib/c progenitors in binary systems are predicted to have a wide range of final masses even up to 10 Msun, with helium envelopes of 0.16 - 1.7 Msun. Our results indicate that, if the lack of helium lines in the spectra of SNe Ic were due to small amounts of helium, the distribution of both initial and final masses of SN Ic progenitors should be bimodal. Furthermore, we find that a thin hydrogen layer (0.001 - 0.01 Msun) is expected to be present in many SN Ib progenitors at the presupernova stage. We show that the presence of hydrogen, together with a rather thick helium envelope, can lead to a significant expansion of some SN Ib/c progenitors by the time of supernova explosion. This may have important consequences for the shock break-out and supernova light curve. We also argue that some SN progenitors with thin hydrogen layers produced via Case AB/B transfer might be related to Type IIb supernova progenitors with relatively small radii of about 10 Rsun.

We compute the contribution of kinks on cosmic string loops to stochastic background of gravitational waves (SBGW).We find that kinks contribute at the same order as cusps to the SBGW.We discuss the accessibility of the total background due to kinks as well as cusps to current and planned gravitational wave detectors, as well as to the big bang nucleosynthesis (BBN), the cosmic microwave background (CMB), and pulsar timing constraints. As in the case of cusps, we find that current data from interferometric gravitational wave detectors, such as LIGO, are sensitive to areas of parameter space of cosmic string models complementary to those accessible to pulsar, BBN, and CMB bounds.

We calculate the screening of the ion-ion potential due to electrons in the presence of a large background magnetic field, at densities of relevance to neutron star crusts. Using the standard approach to incorporate electron screening through the one-loop polarization function, we show that the magnetic field produces important corrections both at short and long distances. In extreme fields, realized in highly magnetized neutron stars called magnetars, electrons occupy only the lowest Landau levels in the relatively low density region of the crust. Here our results show that the screening length for Coulomb interactions between ions can be smaller than the inter-ion spacing. More interestingly, we find that the screening is anisotropic and the screened potential between two static charges exhibits long range Friedel oscillations parallel to the magnetic field. This long-range oscillatory behavior is likely to affect the lattice structure of ions, and can possibly create rod-like structures in the magnetar crusts. We also calculate the imaginary part of the electron polarization function which determines the spectrum of electron-hole excitations and plays a role in damping lattice phonon excitations. We demonstrate that even for modest magnetic fields this damping is highly anisotropic and will likely lead to anisotropic phonon heat transport in the outer neutron star crust.

Recent cosmic ray, gamma ray, and microwave signals observed by Fermi, PAMELA, and WMAP indicate an unexpected primary source of e+e- at 10-1000 GeV. We fit these data to "standard backgrounds" plus a new source, assumed to be a separable function of position and energy. For the spatial part, we consider three cases: annihilating dark matter, decaying dark matter, and pulsars. In each case, we use GALPROP to inject energy in log-spaced energy bins and compute the expected cosmic-ray and photon signals for each bin. We then fit a linear combination of energy bins, plus backgrounds, to the data. We use a non-parametric fit, with no prior constraints on the spectrum except smoothness and non-negativity. In addition, we consider arbitrary modifications to the energy spectrum of the "ordinary" primary source function, fixing its spatial part, finding this alone to be inadequate to explain the PAMELA or WMAP signals. We explore variations in the fits due to choice of magnetic field, primary electron injection index, spatial profiles, propagation parameters, and fit regularization method. Dark matter annihilation fits well, where our fit finds a mass of ~1 TeV and a boost factor times energy fraction of ~70. While it is possible for dark matter decay and pulsars to fit the data, unconventionally high magnetic fields and radiation densities are required near the Galactic Center to counter the relative shallowness of the assumed spatial profiles. We also fit to linear combinations of these three scenarios, though the fit is much less constrained.

We discuss outer gap closure mechanism in the trans-field direction with the magnetic pair-creation process near the stellar surface. The gap closure by the magnetic pair-creation is possible if some fraction of the pairs are produced with an outgoing momentum. By assuming that multiple magnetic field will affect the local field near the stellar surface, we show a specific magnetic field geometry near the stellar surface resulting in the outflow of the pairs. Together with the fact that the electric field is weak below null charge surface, the characteristic curvature photon energy emitted by incoming particles, which were accelerated in the outer gap, decreases drastically to $\sim 100$MeV near the stellar surface. We estimate the height measured from the last-open field line, above which 100MeV photons is converted into pairs by the magnetic pair-creation. We also show the resultant multiplicity due to the magnetic pair-creation process could acquire $M_{e^{\pm}}\sim 10^4-10^5$. In this model the fractional outer gap size is proportional to $P^{-1/2}$. The predicted gamma-ray luminosity ($L_{\gamma}$) and the characteristic curvature photon energy ($E_c$) emitted from the outer gap are proportional to $B^2P^{-5/2}$ and $B^{3/4}P^{-1}$ respectively. This model also predicts that $L_{\gamma}$ and $E_c$ are related to the spin down power ($L_{sd}$) or the spin down age of pulsars ($\tau$) as $L_{\gamma} \propto L_{sd}^{5/8}$ or $L_{\gamma} \propto \tau^{-5/4}$, and $E_c \propto L_{sd}^{1/4}$ or $E_c \propto \tau^{-1/2}$ respectively.

We report results from a deep high-frequency search for pulsars within the central parsec of Sgr A* using the Green Bank Telescope. The observing frequency of 15 GHz was chosen to maximize the likelihood of detecting normal pulsars (i.e. with periods of $\sim 500$\,ms and spectral indices of $\sim -1.7$) close to Sgr A*, that might be used as probes of gravity in the strong-field regime; this is the highest frequency used for such pulsar searches of the Galactic Center to date. No convincing candidate was detected in the survey, with a $10\sigma$ detection threshold of $\sim 10 \mu$Jy achieved in two separate observing sessions. This survey represents a significant improvement over previous searches for pulsars at the Galactic Center and would have detected a significant fraction ($\gtrsim 5%) of the pulsars around Sgr A*, if they had properties similar to those of the known population. Using our best current knowledge of the properties of the Galactic pulsar population and the scattering material toward Sgr A*, we estimate an upper limit of 90 normal pulsars in orbit within the central parsec of Sgr A*.

Swift observations suggest that the central compact objects of some gamma-ray bursts (GRBs) could be newly born millisecond magnetars. Therefore, by considering the spin evolution of the magnetars against r-mode instability, we investigate the role of the magnetars in GRB X-ray afterglow emission. Besides modifying the conventional energy injection model, we pay particular attention to the internal X-ray afterglow emission, whose luminosity is assumed to track the magnetic dipole luminosity of the magentars with a certain fraction. Following a comparison between the model and some selected observational samples, we suggest that some so-called "canonical" X-ray afterglows including the shallow decay, normal decay, and steeper-than-normal decay phases could be internally produced by the magnetars (possibly through some internal dissipations of the magnetar winds), while the (energized) external shocks are associated with another type of X-ray afterglows. If this is true, from those internal X-ray afterglows, we can further determine the magnetic field strengths and the initial spin periods of the corresponding magnetars.

We study the magnetosphere of a slowly rotating magnetized neutron star subject to toroidal oscillations in the relativistic regime. Under the assumption of a zero inclination angle between the magnetic moment and the angular momentum of the star, we analyze the Goldreich-Julian charge density and derive a second-order differential equation for the electrostatic potential. The analytical solution of this equation in the polar cap region of the magnetosphere shows the modification induced by stellar toroidal oscillations on the accelerating electric field and on the charge density. We also find that, after decomposing the oscillation velocity in terms of spherical harmonics, the first few modes with $m=0,1$ are responsible for energy losses that are almost linearly dependent on the amplitude of the oscillation and that, for the mode $(l,m)=(2,1)$, can be a factor $\sim8$ larger than the rotational energy losses, even for a velocity oscillation amplitude at the star surface as small as $\eta=0.05 \ \Omega \ R$. The results obtained in this paper clarify the extent to which stellar oscillations are reflected in the time variation of the physical properties at the surface of the rotating neutron star, mainly by showing the existence of a relation between $P\dot{P}$ and the oscillation amplitude. Finally, we propose a qualitative model for the explanation of the phenomenology of intermittent pulsars in terms of stellar oscillations that are periodically excited by star glitches.

Supergiant Fast X-ray Transients (SFXTs) are a new class of High Mass X-ray Binaries, discovered by the INTEGRAL satellite, which display flares lasting from minutes to hours, with peak luminosity of 1E36-1E37 erg/s. Outside the bright outbursts, they show a frequent long-term flaring activity reaching an X-ray luminosity level of 1E33-1E34 erg/s, as recently observed with the Swift satellite. Since a few persistent High Mass X-ray Binaries (HMXBs) with supergiant donors show flares with properties similar to those observed in SFXTs, it has been suggested that the flaring activity in both classes could be produced by the same mechanism, probably the accretion of clumps composing the supergiant wind. We have developed a new clumpy wind model for OB supergiants with both a spherical and a non spherical symmetry for the outflow. We have investigated the effects of the accretion of a clumpy wind onto a neutron star in both classes of persistent and transient HMXBs.

We analyzed two XMM-Newton observations of the bright atoll source 4U 1705-44, which can be considered a prototype of the class of the persistent NS LMXBs showing both hard and soft states. The first observation was performed when the source was in a hard low flux state, the second during a soft, high-flux state. Both the spectra show broad iron emission lines. We fit the spectra using a two-component model, together with a reflection model specifically suited to the case of a neutron star, where the incident spectrum has a blackbody shape. In the soft state, the reflection model, convolved with a relativistic smearing component, consistently describes the broad features present in the spectrum, and we find a clear relation between the temperature of the incident flux and the temperature of the harder X-ray component that we interpret as the boundary layer emission. In this state we find converging evidence that the boundary layer outer radius is ~ 2 times the neutron star radius. In the low flux state, we observe a change in the continuum shape of the spectrum with respect to the soft state. Still, the broad local emission features can be associated with a disk reflecting matter, but in a lower ionization state, and possibly produced in an accretion disk truncated at greater distance. Our analysis provides strong evidence that the reflection component in soft states of LMXBs comes from to hard X-ray thermal irradiation, which we identify with the boundary layer emission, also present in the continuum model. In the hard state, the broad iron line if also produced by reflection, and the continuum disk emission can be self-consistently accounted if the disk is truncated at a greater distance than the soft state.

Context. Supergiant fast X-ray transients are a subclass of high mass X-ray binaries which host a neutron star accreting mass from the wind of its OB supergiant companion. They are characterized by an extremely pronounced and rapid variability in X-rays, which still lacks an unambiguous interpretation. Aims. Recently, a number of deep pointed observations with XMM-Newton have been carried out in order to study the quiescent emission of these sources and gain further insight on the mechanism that causes their X-ray variability. Methods. In this paper, we analyzed the first three deep-pointed XMM-Newton observations of the two supergiant fast X-ray transient XTE J1739-302 and IGR J08408-4503 in quiescence. Results. We found that the quiescent emission of these sources is characterized by a complex timing and spectral variability, with multiple small flares occurring sporadically after periods of lower X-ray emission. Some evidence is found in the XMM-Newton spectra for the presence of a soft component below ~2 keV, similar to that observed in the two supergiant fast X-ray transients AX J1845.0-0433 and IGR J16207-5129 and in many other high mass X-ray binaries. Conclusions. We suggest some possible interpretations of the timing and spectral properties of the quiescent emission of XTE J1739- 302 and IGR J08408-4503 in the context of the different theoretical models proposed to interpret the behavior of the supergiant fast X-ray transients.

Deep observations of the Small Magellanic Cloud (SMC) and region were carried out in the hard X-ray band by the INTEGRAL observatory in 2008-2009. The field of view of the instrument permitted simultaneous coverage of the entire SMC and the eastern end of the Magellanic Bridge. In total, INTEGRAL detected seven sources in the SMC and five in the Magellanic Bridge; the majority of the sources were previously unknown systems. Several of the new sources were detected undergoing bright X- ray outbursts and all the sources exhibited transient behaviour except the supergiant system SMC X-1. They are all thought to be High Mass X-ray Binary (HMXB) systems in which the compact object is a neutron star.

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 investigate a model for the shallow decay phases of Gamma-ray Burst (GRB) afterglows discovered by Swift/XRT in the first hours following a GRB event. In the context of the fireball scenario, we consider the possibility that long-lived energy injection from a millisecond spinning, ultramagnetic neutron star (magnetar) powers afterglow emission during this phase. We consider the energy evolution in a relativistic shock subject to both radiative losses and energy injection from a spinning down magnetar in spherical symmetry. We model the energy injection term through magnetic dipole losses and discuss an approximate treatment for the dynamical evolution of the blastwave. We obtain an analytic solution for the energy evolution in the shock and associated lightcurves. To fully illustrate the potential of our solution we calculate lightcurves for a few selected X-ray afterglows observed by Swift and fit them using our theoretical lightcurves. Our solution naturally describes in a single picture the properties of the shallow decay phase and the transition to the so-called normal decay phase. In particular, we obtain remarkably good fits to X-ray afterglows for plausible parameters of the magnetar. Even though approximate, our treatment provides a step forward with respect to previously adopted approximations and provides additional support to the idea that a millisecond spinning (1-3 ms), ultramagnetic (B$\sim 10^{14}-10^{15}$ G) neutron star loosing spin energy through magnetic dipole radiation can explain the luminosity, durations and shapes of X-ray GRB afterglows.

The Crab pulsar is well-known for its anomalous giant radio pulse emission. Past studies have concentrated only on the very bright pulses or were insensitive to the faint end of the giant pulse luminosity distribution. With our new instrumentation offering a large bandwidth and high time resolution combined with the narrow radio beam of the Westerbork Synthesis Radio Telescope (WSRT), we seek to probe the weak giant pulse emission regime. The WSRT was used in a phased array mode, resolving a large fraction of the Crab nebula. The resulting pulsar signal was recorded using the PuMa II pulsar backend and then coherently dedispersed and searched for giant pulse emission. After careful flux calibration, the data were analysed to study the giant pulse properties. The analysis includes the distributions of the measured pulse widths, intensities, energies, and scattering times. The weak giant pulses are shown to form a separate part of the intensity distribution. The large number of giant pulses detected were used to analyse scattering and scintillation in giant pulses. We report for the first time the detection of giant pulse emission at both the main- and interpulse phases within a single rotation period. The rate of detection is consistent with the appearance of pulses at either pulse phase as being independent. These pulse pairs were used to examine the scintillation timescales within a single pulse period.

We present broadband (radio, optical, and X-ray) light curves and spectra of the afterglows of four long-duration gamma-ray bursts (GRBs 090323, 090328, 090902B, and 090926A) detected by the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) instruments on the Fermi satellite. With its wide spectral bandpass, extending to GeV energies, Fermi is sensitive to GRBs with very large isotropic energy releases (10e54 erg). Although rare, these events are particularly important for testing GRB central-engine models. When combined with spectroscopic redshifts, our afterglow data for these four events are able to constrain jet collimation angles, the density structure of the circumburst medium, and both the true radiated energy release and the kinetic energy of the outflows. In agreement with our earlier work, we find that the relativistic energy budget of at least one of these events (GRB 090926A) exceeds the canonical value of 10e51 erg by an order of magnitude. Such energies pose a severe challenge for models in which the GRB is powered by a magnetar or neutrino-driven collapsar, but remain compatible with theoretical expectations for magneto-hydrodynamical (MHD) collapsar models. Our jet opening angles (theta) are similar to those found for pre-Fermi GRBs, but the large initial Lorentz factors (Gamma_0) inferred from the detection of GeV photons imply theta Gamma_0 ~ 70-90, values which are above those predicted in MHD models of jet acceleration. Finally, we find that these Fermi-LAT events preferentially occur in a low-density circumburst environment, and we speculate that this might result from the lower mass-loss rates of their lower-metallicity progenitor stars. Future studies of Fermi-LAT afterglows in the radio with the order-of-magnitude improvement in sensitivity offered by the EVLA should definitively establish the relativistic energy budgets of these events.

We present observations of HESS J1640-465 with the Fermi-LAT. The source is detected with high confidence as an emitter of high-energy gamma-rays. The spectrum lacks any evidence for the characteristic cutoff associated with emission from pulsars, indicating that the emission arises primarily from the pulsar wind nebula. Broadband modeling implies an evolved nebula with a low magnetic field resulting in a high gamma-ray to X-ray flux ratio. The Fermi emission exceeds predictions of the broadband model, and has a steeper spectrum, possibly resulting from a distinct excess of low energy electrons similar to what is inferred for both the Vela X and Crab pulsar wind nebulae.

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are $\sim 10^{8}$G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The $\sim 10^{8}$ G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have $< 10^{11}$ G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit $10^{11}$ G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict $\sim 10^{8}$ G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted $\sim 10^{8}$ G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields $\gg 10^{8}$ G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields $> 10^{8}$ G with their magnetic dipole emission being balanced by low level accretion from their companions.

We study the spectral evolution of PWNe taking into account the energy injected when they are young. We model the evolution of the magnetic field inside a uniformly expanding PWN. Considering time dependent injection from the pulsar and coolings by radiative and adiabatic losses, we solve the evolution of the particle distribution function. The model is calibrated by fitting the calculated spectrum to the observations of the Crab Nebula at an age of a thousand years. The spectral evolution of the Crab Nebula in our model shows that the flux ratio of TeV gamma-rays to X-rays increases with time, which implies that old PWNe are faint in X-rays, but not in TeV gamma-rays. The increase of this ratio is because the magnetic field decreases with time and is not because the X-ray emitting particles are cooled more rapidly than the TeV gamma-ray emitting particles. Our spectral evolution model matches the observed rate of the radio flux decrease of the Crab Nebula. This result implies that our magnetic field evolution model is close to the reality. Finally, from the viewpoint of the spectral evolution, only a small fraction of the injected energy from the Crab Pulsar needs to go to the magnetic field, which is consistent with previous studies.

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, we use the 3D self-consistent annular gap model to study light curves for both young and millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap can generate high energy emission for short-period pulsars. The annular gap regions are so large that they have enough electric potential drop to accelerate charged particles to produce gamma-ray photons. For young pulsars, the emission region is from the neutron star surface to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null charge surface. The emission region for the millisecond pulsars is located much lower than that of the young pulsars. The higher energy gamma-ray emission comes from higher altitudes in the magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J0218+4232, PSR J0437-3715) using the annular gap model. Our simulations can reproduce the main properties of observed light curves.

We show that several features reminiscent of short-hard GRBs arise naturally when Quark-Novae occur in post-accretion low-mass X-ray binaries with a circumbinary disk. Post-accretion conditions in a neutron star-white dwarf binary are just right for the conversion of the neutron star to a quark star (Quark-Nova). In our model, the subsequent interaction of material from the neutron star's ejected crust with the circumbinary disk explains the duration, variability and near-universal nature of the prompt emission in short-hard GRBs. We also describe a statistical approach to ejecta break-up and collision to obtain the photon spectrum in our model, which turns out remarkably similar to the empirical Band function (Band et al. 1993). We apply the model to the fluence and spectrum of GRB 000727, GRB 000218, and GRB980706A obtaining excellent fits. Extended emission (spectrum and duration) is explained by shock-heating and ablation of the white dwarf by the highly energetic ejecta. Depending on the orbital separation when the Quark-Nova occurs, we isolate interesting regimes within our model when both prompt and extended emission can occur. We find that the spectrum can carry signatures typical of Type Ib/c SNe, thus providing an alternative to the collapsar scenario. Late X-ray activity is due to accretion onto the quark star as well as its spin-down luminosity, while afterglows arise from the expanding shell of material from the shock-heated expanding circumbinary disk. We find a correlation between the duration and spectrum of short-hard GRBs as well as modest hard-to-soft time evolution of the peak energy.

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) exhibit characteristic X-ray luminosities (both soft and hard) of the order of $10^{35}$ erg s$^{-1}$ and characteristic power-law hard X-ray spectra extending to about 200 keV. Two AXPs exhibit also pulsed radio emission. Assuming that AXPs and SGRs accrete matter from a fall-back disk, we attempt to explain both the soft and the hard X-ray emission as a result of the accretion process. We also attempt to explain their radio emission or the lack of it. We propose that the power-law, hard X-ray spectra are produced in the accretion flow mainly by bulk motion Comptonization of soft photons emitted at the neutron star surface. Unlike normal X-ray pulsars, for which the accretion rate is highly super-Eddington, in AXPs and SGRs the accretion rate is approximately Eddington and thus the bulk motion Comptonization operates efficiently. Fall-back disk models invoke surface dipole magnetic fields of $10^{12} - 10^{13}$ G and this is what we assume here. As an illustrative example we reproduce extremely well both the hard and the soft X-ray spectrum of AXP 4U 0142+61 using the XSPEC package compTB. Our model seems to explain in a natural way both the hard and the soft X-ray spectra of AXPs and SGRs as well as their radio emission or the lack of it. It can also possibly explain the short bursts observed in these sources. On the other hand, it cannot explain the giant X-ray outbursts observed in SGRs, which may be due to the conversion of magnetic energy in local multipole fields.

Efforts to detect gravitational waves by timing an array of pulsars have focused traditionally on stationary gravitational waves: e.g., stochastic or periodic signals. Gravitational wave bursts --- signals whose duration is much shorter than the observation period --- will also arise in the pulsar timing array waveband. Sources that give rise to detectable bursts include the formation or coalescence of supermassive black holes (SMBHs), the periapsis passage of compact objects in highly elliptic or unbound orbits about a SMBH, or cusps on cosmic strings. Here we describe how pulsar timing array data may be analyzed to detect and characterize these bursts. Our analysis addresses, in a mutually consistent manner, a hierarchy of three questions: \emph{i}) What are the odds that a dataset includes the signal from a gravitational wave burst? \emph{ii}) Assuming the presence of a burst, what is the direction to its source? and \emph{iii}) Assuming the burst propagation direction, what is the burst waveform's time dependence in each of its polarization states? Applying our analysis to synthetic data sets we find that we can \emph{detect} gravitational waves even when the radiation is too weak to either localize the source of infer the waveform, and \emph{detect} and \emph{localize} sources even when the radiation amplitude is too weak to permit the waveform to be determined. While the context of our discussion is gravitational wave detection via pulsar timing arrays, the analysis itself is directly applicable to gravitational wave detection using either ground or space-based detector data.

Recent X-ray observations have proved to be very effective in detecting previously unknown supernova remnant shells around pulsar wind nebulae (PWNe), and in these cases the characteristics of the shell provide further clues on the evolutionary stage of the embedded PWN. However, it is not clear why some PWNe are still "naked". We carried out an X-ray observational campaign targeted at the PWN G54.1+0.3, the "close cousin" of the Crab, with the aim to detect the associated SNR shell. We analyzed an XMM-Newton and Suzaku observations of G54.1+0.3 and we model out the contribution of dust scattering halo. We detected an intrinsic faint diffuse X-ray emission surrounding a hard spectrum, which can be modeled either with a power-law (gamma= 2.9) or with a thermal plasma model (kT=2.0 keV.). If the shell is thermal, we derive an explosion energy E=0.5-1.6x10^51 erg, a pre-shock ISM density of 0.2 cm^-3 and an age of about 2000 yr. Using these results in the MHD model of PWN-SNR evolution, we obtain an excellent agreement between the predicted and observed location of the shell and PWN shock.

Pulsars are very stable clocks in space which have many applications to problems in physics and astrophysics. Observations of double-neutron-star binary systems have given the first observational evidence for the existence of gravitational waves (GWs) and shown that Einstein's general theory of relativity is an accurate description of gravitational interactions in the regime of strong gravity. Observations of a large sample of pulsars spread across the celestial sphere forming a "Pulsar Timing Array" (PTA), can in principle enable a positive detection of the GW background in the Galaxy. The Parkes Pulsar Timing Array (PPTA) is making precise timing measurements of 20 millisecond pulsars at three radio frequencies and is approaching the level of timing precision and data spans which are needed for GW detection. These observations will also allow us to establish a "Pulsar Timescale" and to detect or limit errors in the Solar System ephemerides used in pulsar timing analyses. Combination of PPTA data with that of other groups to form an International Pulsar Timing Array (IPTA) will enhance the sensitivity to GWs and facilitate reaching other PTA goals. The principal source of GWs at the nanoHertz frequencies to which PTAs are sensitive is believed to be super-massive binary black holes in the cores of distant galaxies. Current results do not signficantly limit models for formation of such black-hole binary systems, but in a few years we expect that PTAs will either detect GWs or seriously constrain current ideas about black-hole formation and galaxy mergers. Future instruments such as the Square Kilometre Array (SKA) should not only detect GWs from astrophysical sources but also enable detailed studies of the sources and the gravitational theories used to account for the GW emission.

We present results on the Soft Gamma Repeater (SGR) 0501+4516, discovered by the SWIFT Burst Alert Telescope (BAT) on 2008 August 22. More than 50 bursts were identified from this source, out of which 18 bursts had enough counts to carry out spectral analysis. We performed time-averaged spectral analysis on these 18 bursts using 8 models, among which the cut-off powerlaw and the two-blackbody models provided the best fit in the 15-150 keV energy range. The cut-off powerlaw model fit yields a mean photon index Gamma_{CPL} = 0.54+/-0.11 and a cut-off energy E_C = 19.1+/-1.8 keV for the bursts. The mean hard and soft blackbody temperatures are found to be kT_{BB_h} = 12.8+/-0.7 keV and kT_{BB_s} = 4.6+/-0.5 keV, respectively, and are anti-correlated with the square of the radii of the hard and soft emitting regions (R_{BB_h} and R_{BB_s}) as R_{BB_h}^2 $\propto$ kT^{-5.8} and R_{BB_s}^2 $\propto$ kT^{-2.7}, respectively. The soft and hard component temperatures with different indices support the idea of two distinct emitting regions with the hard component corresponding to a smaller radius and the soft component corresponding to a larger radius, which further corroborate the idea of the propagation of extraordinary (E) and ordinary (O) mode photons across the photosphere, as predicted in the magnetar model. We notice strong burst fluence-duration correlation as well as hardness ratio-duration and hardness ratio-fluence anti-correlations for the SGR 0501+4516 bursts. The burst fluences range from ~ 4.4\times10^{-9} ergs/cm^-2 to ~ 2.7\times10^{-6} ergs/cm^{-2}, consistent with those observed for typical short SGR bursts.

We study time evolutions of superfluid neutron stars, focussing on the nature of the oscillation spectrum, the effect of mutual friction force on the oscillations and the hydrodynamical spin-up phase of pulsar glitches. We linearise the dynamical equations of a Newtonian two-fluid model for rapidly rotating backgrounds. In the axisymmetric equilibrium configurations, the two fluid components corotate and are in beta-equilibrium. We use analytical equations of state that generate stratified and non-stratified stellar models, which enable us to study the coupling between the dynamical degrees of freedom of the system. By means of time evolutions of the linearised dynamical equations, we determine the spectrum of axisymmetric and non-axisymmetric oscillation modes, accounting for the contribution of the gravitational potential perturbations, i.e. without adopting the Cowling approximation. We study the mutual friction damping of the superfluid oscillations and consider the effects of the non-dissipative part of the mutual friction force on the mode frequencies. We also provide technical details and relevant tests for the hydrodynamical model of pulsar glitches discussed by Sidery, Passamonti and Andersson (2010). In particular, we describe the method used to generate the initial data that mimic the pre-glitch state, and derive the equations that are used to extract the gravitational-wave signal.

Context: Chemical element abundances for distant Galactic globular clusters (GCs) hold important clues to the origin of the Milky Way halo and its substructures. Aims: We study the chemical composition of red giant stars in Pal 4 - one of the most remote GCs in the Milky Way - and compare our abundance measurements to those for both low surface brightness dwarf galaxies, and GCs in the inner and the outer halo. Methods: By co-adding high-resolution, low-S/N Keck/HIRES spectra of 19 stars along the red giant branch, we estimate chemical abundance ratios of 20 alpha-, iron peak-, and neutron-capture elements. Our method gives total uncertainties on most element-to-iron ratios of typically 0.2 dex. Results: We measure [Fe/H] = -1.41 +- 0.04 (statistical) +- 0.17 (systematic) and an alpha-enhancement of [alpha/Fe] = +0.38 +- 0.11 dex, which is consistent with the canonical value of +0.4 dex found for Galactic halo field stars and most halo GCs at this metallicity. Although Pal 4 has higher enhancements in the heavier elements with respect to the halo, the majority of the element ratios are, within the measurement errors, consistent with those for local halo field stars. We find, however, evidence for a lower [Mg/Ca] ratio than in other halo clusters. Conclusions: Based on the available evidence, we conclude that the material from which Pal 4 and the Galactic halo formed experienced similar enrichment processes, despite the apparently younger age of this cluster. Within the limitations of our methodology we find no significant indication of an iron spread, as is typical of genuine GCs of the Milky Way. However, abundance ratios for individual stars in Pal 4 and other distant satellites are urgently needed to understand the relationship, if any, between remote GCs and other halo substructures (i.e., luminous and ultra-faint dwarf spheroidal galaxies).

The pulse profiles of the transient Be/X-ray binary EXO 2030+375 show strong dependence on energy, as well as on its luminosity state, and are asymmetric in shape. We want to identify the emission components of the two magnetic poles in the pulsed emission to understand the geometry of the neutron star and its beam pattern. We utilize a pulse-profile decomposition method that enables us to find two symmetric pulse profiles from the magnetic poles of the neutron star. The symmetry characteristics of these single-pole pulse profiles give information about the position of the magnetic poles of the neutron star relative to its rotation axis. We find a possible geometry for the neutron star in EXO 2030+375 through the decomposition of the pulse profiles, which suggests that one pole gets closer to the line of sight than the other and that, during the revolution of the neutron star, both poles disappear behind the horizon for a short period of time. A considerable fraction of the emission arises from a halo while the pole is facing the observer and from the accretion stream of the other pole while it is behind the neutron star, but the gravitational line bending makes the emission visible to us.

We study evolution of isolated neutron stars on long time scale and calculate distribution of these sources in the main evolutionary stages: Ejector, Propeller, Accretor, and Georotator. We compare different initial magnetic field distributions taking into account a possibility of magnetic field decay, and include in our calculations the stage of subsonic Propeller. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0\ga 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.

Coalescing binary systems, consisting of two collapsed objects, are among the most promising sources of high frequency gravitational waves signals detectable, in principle, by ground-based interferometers. Binary systems of Neutron Star or Black Hole/Neutron Star mergers should also give rise to short Gamma Ray Bursts, a subclass of Gamma Ray Bursts. Short-hard-Gamma Ray Bursts might thus provide a powerful way to infer the merger rate of two-collapsed object binaries. Under the hypothesis that most short Gamma Ray Bursts originate from binaries of Neutron Star or Black Hole/Neutron Star mergers, we outline here the possibility to associate short Gamma Ray Bursts as electromagnetic counterpart of coalescing binary systems.

A precise electromagnetic measurement of the sky coordinates and redshift of a coalescing black hole binary holds the key for using its gravitational wave (GW) signal to constrain cosmological parameters and to test general relativity. Here we show that the merger of ~10^{6-8}M_sun black holes is generically followed over a period of years by multiple electromagnetic flares from tidally disrupted stars. The sudden recoil imparted to the merged black hole by GW emission results promptly in a tidal disruption rate of stars as high as ~0.1-1 per year. The sequential disruption of stars within a single galaxy over a short period provides a unique electromagnetic flag of a recent black hole coalescence event, and can be used on its own to calibrate the expected rate of GW sources for pulsar timing arrays or the proposed Laser Interferometer Space Antenna (LISA).

Thermal emission during X-ray bursts is a powerful tool to determine neutron star masses and radii, if the Eddington flux and the apparent radius in the cooling tail can be measured accurately, and distances to the sources are known. We propose here an improved method of determining the basic stellar parameters using the data from the cooling phase of long, photospheric radius expansion bursts covering a large range of luminosities. For this purpose, we computed a large set of atmosphere models for burst luminosities varying by two orders of magnitude and for various chemical compositions and surface gravities. We show that the variation of the inverse square root of the apparent blackbody radius with the flux, observed during the photospheric radius expansion burst from 4U 1724-307 located in globular cluster Terzan 2, is entirely consistent with the theoretical expectations of the color-correction factor evolution. Our method allows us to determine both the Eddington flux and the ratio of the stellar apparent radius to the distance much more reliably. We then find a lower limit on the neutron star radius of 13 km, independently of the chemical composition. These results suggest that the matter inside neutron stars is characterized by a stiff equation of state. We also find evidences in favor of hydrogen rich accreting matter and obtain an upper limit to the distance of 7 kpc. Our approach improves the old way of distances determination to X-ray bursters using Eddington fluxes.

Although they are but a small fraction of the mass ejected in core-collapse supernovae, neutrino-driven winds (NDWs) from nascent proto-neutron stars (PNSs) have the potential to contribute significantly to supernova nucleosynthesis. In previous works, the NDW has been implicated as a possible source of r-process and light p-process isotopes. In this paper we present time-dependent hydrodynamic calculations of nucleosynthesis in the NDW which include accurate weak interaction physics coupled to a full nuclear reaction network. Using two published models of PNS neutrino luminosities, we predict the contribution of the NDW to the integrated nucleosynthetic yield of the entire supernova. For the neutrino luminosity histories considered, no true r-process occurs in the most basic scenario. The wind driven from an older $1.4 M_\odot$ model for a PNS is moderately neutron-rich at late times however, and produces $^{87}$Rb, $^{88}$Sr, $^{89}$Y, and $^{90}$Zr in near solar proportions relative to oxygen. The wind from a more recently studied $1.27 M_\odot$ PNS is proton-rich throughout its entire evolution and does not contribute significantly to the abundance of any element. It thus seems very unlikely that the simplest model of the NDW can produce the r-process. At most, it contributes to the production of the N = 50 closed shell elements and some light p-nuclei. In doing so, it may have left a distinctive signature on the abundances in metal poor stars, but the results are sensitive to both uncertain models for the explosion and the masses of the neutron stars involved.

Gravitational waves from the final stages of inspiralling binary neutron stars are expected to be one of the most important sources for ground-based gravitational wave detectors. The masses of the components are determinable from the orbital and chirp frequencies during the early part of the evolution during which tidal effects provide small correction; however, during this phase the signal is relatively clean. The accumulated phase shift due to tidal corrections is characterized by a single quantity, the Love number, which is sensitive to the compactness parameter M/R and the star's internal structure, and its determination could constrain the star's radius. We show that the Love number of normal neutron stars are much different from those of self-bound strange quark matter stars and could therefore provide an important way to distinguish between these two classes of stars.

We discuss a possible generation of radio bursts preceding final stages of binary neutron star mergings which can be accompanied by short gamma-ray bursts. Detection of such bursts appear to be advantageous in the low-frequency radio band due to a time delay of ten to several hundred seconds required for radio signal to propagate in the ionized intergalactic medium. This delay makes it possible to use short gamma-ray burst alerts to promptly monitor specific regions on the sky by low-frequency radio facilities, especially by LOFAR. To estimate the strength of the radio signal, we assume a power-law dependence of the radio luminosity on the total energy release in a magnetically dominated outflow, as found in millisecond pulsars. Based on the planned LOFAR sensitivity at 120 MHz, we estimate that the LOFAR detection rate of such radio transients could be about several events per month from redshifts up to $z\sim1.3$ in the most optimistic scenario. The LOFAR ability to detect such events would crucially depend on exact efficiency of low-frequency radio emission mechanism.

X-ray Flash (XRF) 100316D, a nearby super-long under-luminous burst with a peak energy E_p \sim 20 keV, was detected by Swift and was found to be associated with an energetic supernova SN 2010bh. Both the spectral and the temporal behavior of this burst are rather similar to that of XRF 060218, except that the latter was associated with a "less energetic" SN 2006aj and had a prominent soft thermal emission component in the spectrum. We analyze the spectral and temporal properties of this burst, and interpret the prompt gamma-ray emission and the early X-ray plateau emission as synchrotron emission from a dissipating Poynting-flux-dominated outflow, probably powered by a magnetar with a spin period of $P \sim 10$ ms and the polar cap magnetic field $B_{\rm p} \sim 3\times 10^{15}$ G. The energetic supernova SN 2010bh associated with this burst is, however, difficult to interpret within the slow magnetar model, which implies that the nascent magnetar may spin much faster with an initial rotation period $\sim 1$ ms, and thus suggests a delay between the core collapse and the emergence of the relativistic magnetar wind from the star. The diverse behaviors of low-luminosity GRBs and their associated SNe may be understood within a unified picture that invokes different initial powers of the central engine and different delay times between the core collapse and the emergence of the relativistic jet from the star.

We present near-infrared (NIR) broadband (0.80--2.42 $\mu$m) spectroscopy of two low mass X-ray binaries: V404 Cyg and Cen X-4. One important parameter required in the determination of the mass of the compact objects in these systems is the binary inclination. We can determine the inclination by modeling the ellipsoidal modulations of the Roche-lobe filling donor star, but the contamination of the donor star light from other components of the binary, particularly the accretion disk, must be taken into account. To this end, we determined the donor star contribution to the infrared flux by comparing the spectra of V404 Cyg and Cen X-4 to those of various field K-stars of known spectral type. For V404 Cyg, we determined that the donor star has a spectral type of K3 III. We determined the fractional donor contribution to the NIR flux in the H- and K-bands as $0.98 \pm .05$ and $0.97 \pm .09$, respectively. We remodeled the H-band light curve from \citet{sanwal1996} after correcting for the donor star contribution to obtain a new value for the binary inclination. From this, we determined the mass of the black hole in V404 Cyg to be $M_{BH}= 9.0^{+.2}_{-.6}M_{\odot}$. We performed the same spectral analysis for Cen X-4 and found the spectral type of the donor star to be in the range K5 -- M1V. The donor star contribution in Cen X-4 is $0.94\pm.14$ in the H-band while in the K-band, the accretion disk can contribute up to 10% of the infrared flux. We remodeled the H-band light curve from \citet{shahbaz1993}, again correcting for the fractional contribution of the donor star to obtain the inclination. From this, we determined the mass of the neutron star as $M_{NS}= 1.5^{+.1}_{-.4}M_{\odot}$. However, the masses obtained for both systems should be viewed with some caution since contemporaneous light curve and spectral data are required to obtain definitive masses.