We have embarked on a survey for pulsars and fast transients using the 13-beam Multibeam receiver on the Parkes radio telescope. Installation of a digital backend allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels and sampled every 64 us. Limits of the receiver package restrict us to a 340 MHz observing band centred at 1352 MHz. The factor of eight improvement in frequency resolution over previous multibeam surveys allows us to probe deeper into the Galactic plane for short duration signals such as the pulses from millisecond pulsars. We plan to survey the entire southern sky in 42641 pointings, split into low, mid and high Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond pulsars. With ~30% of the mid-latitude survey complete, we have re-detected 223 previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. The newly discovered millisecond pulsars tend to have larger dispersion measures than those discovered in previous surveys, as expected from the improved time and frequency resolution of our instrument.

The S stars orbiting the Galactic center black hole reach speeds of up to a few percent the speed of light during pericenter passage. This makes, for example, S2 at pericenter much more relativistic than known binary pulsars, and opens up new possibilities for testing general relativity. This paper develops a technique for fitting nearly-Keplerian orbits with perturbations from Schwarzschild curvature, frame dragging, and spin-induced torque, to redshift measurements distributed along the orbit but concentrated around pericenter. Both orbital and light-path effects are taken into account. It turns out that absolute calibration of rest-frame frequency is not required. Hence, if pulsars on orbits similar to the S stars are discovered, the technique described here can be applied without change, allowing the much greater accuracies of pulsar timing to be taken advantage of. For example, pulse timing of 3 microsec over one hour amounts to an effective redshift precision of 30 cm/s, enough to measure frame dragging and the quadrupole moment from an S2-like orbit, provided problems like the Newtonian "foreground" due to other masses can be overcome. On the other hand, if stars with orbital periods of order a month are discovered, the same could be accomplished with stellar spectroscopy from the E-ELT at the level of 1 km/s.

We report the detection of a large glitch in the pulsar B0919+06 (J0922+0638). The glitch occurred in 2009 November 5 (MJD 55140) and was characterized by a fractional increase in the rotation frequency of Deltanu/nu=1.3x10^{-6}. A large glitch happens in the pulsar whose rotation has unstable character. We present the results of the analysis of the rotation behavior of this pulsar over the 30-year time span from 1979 to 2009. These results show that the pulsar's rotation frequency underwent continuous, slow oscillations which look like glitch-like events. During the 1991-2009 interval, the pulsar experienced a continuous sequence of 12 slow glitches with a fractional increase in the rotation frequency Deltanu/nu=1.5x10^{-9}. All the slow glitches observed have a similar signature related to a slow increase in the rotation frequency during 200 days and the subsequent relaxation back to the pre-glitch value during 400 days. We show that a continuous sequence of such slow glitches is characterized by practically identical amplitudes equal to Deltanu=3.5x10^{-9} Hz and identical time intervals between glitches of about 600 days and is well described by a periodic sawtooth-like function. The detection of two different phenomena, such as a large glitch and a sequence of slow glitches, indicates the presence of two types of discontinuities in the rotation frequency of the pulsar B0919+06. These discontinuities can be classified as normal and slow glitches.

In late 2008, the quasi-persistent neutron star X-ray transient and eclipsing binary EXO 0748-676 started a transition from outburst to quiescence, after it had been actively accreting for more than 24 years. In a previous work, we discussed Chandra and Swift observations obtained during the first five months after this transition. Here, we report on further X-ray observations of EXO 0748-676, extending the quiescent monitoring to 1.6 years. Chandra and XMM-Newton data reveal quiescent X-ray spectra composed of a soft, thermal component that is well-fitted by a neutron star atmosphere model. An additional hard powerlaw tail is detected that changes non-monotonically over time, contributing between 4 and 20 percent to the total unabsorbed 0.5-10 keV flux. The combined set of Chandra, XMM-Newton and Swift data reveals that the thermal bolometric luminosity fades from ~1E34 to 6E33 (d/7.4 kpc)^2 erg/s, whereas the inferred neutron star effective temperature decreases from ~124 to 109 eV. We interpret the observed decay as cooling of the neutron star crust and show that that the quiescent lightcurve of EXO 0748-676 is markedly shallower than that observed for three other neutron star X-ray binaries that underwent prolonged accretion outbursts.

We report on continued monitoring observations of the Galactic center carried out by the X-ray telescope aboard the Swift satellite in 2008 and 2009. This campaign revealed activity of the five known X-ray transients AX J1745.6-2901, CXOGC J174535.5-290124, GRS 1741-2853, XMM J174457-2850.3 and CXOGC J174538.0-290022. All these sources are known to undergo very faint X-ray outbursts with 2-10 keV peak luminosities of Lx~1E34-1E36 erg/s, although the two confirmed neutron star low-mass X-ray binaries AX J1745.6-2901 and GRS 1741-2853 can also become brighter (Lx~1E36-1E37 erg/s). We discuss the observed long-term lightcurves and X-ray spectra of these five enigmatic transients. In 2008, AX J1745.6-2901 returned to quiescence following an unusually long accretion outburst of ~1.5 years. GRS 1741-2853 was active in 2009 and displayed the brightest outburst ever recorded for this source, reaching up to a 2-10 keV luminosity of Lx~1E37 (D/7.2 kpc)^2 erg/s. This system appears to undergo recurrent accretion outbursts every ~2 years. Furthermore, we find that the unclassified transient XMM J174457-2850.3 becomes bright only during short episodes (days) and is often found active in between quiescence (Lx~1E32 erg/s) and its maximum outburst luminosity of Lx~1E36 erg/s. CXOGC J174535.5-290124 and CXOGC J174538.0-290022, as well as three other very-faint X-ray transients that were detected by Swift monitoring observations in 2006, have very low time-averaged mass-accretion rates of ~< 2E-12 Msun/yr. Despite having obtained two years of new data in 2008 and 2009, no new X-ray transients were detected.

Recent studies of black hole and neutron star low mass X-ray binaries (LMXBs) show a positive correlation between the X-ray flux at which the low/hard(LH)-to-high/soft(HS) state transition occurs and the peak flux of the following HS state. By analyzing the data from the All Sky Monitor (ASM) onboard the Rossi X-ray Timing Explorer (RXTE), we show that the HS state flux after the source reaches its HS flux peak still correlates with the transition flux during soft X-ray transient (SXT) outbursts. By studying large outbursts or flares of GX 339-4, Aql X-1 and 4U 1705-44, we have found that the correlation holds up to 250, 40, and 50 d after the LH-to-HS state transition, respectively. These time scales correspond to the viscous time scale in a standard accretion disk around a stellar mass black hole or a neutron star at a radius of ~104-5 Rg, indicating that the mass accretion rates in the accretion flow either correlate over a large range of radii at a given time or correlate over a long period of time at a given radius. If the accretion geometry is a two-flow geometry composed of a sub-Keplerian inflow or outflow and a disk flow in the LH state, the disk flow with a radius up to ~105 Rg would have contributed to the nearly instantaneous non-thermal radiation directly or indirectly, and therefore affects the time when the state transition occurs.

We use a simple two-layer outer gap model, whose accelerator consists of a primary region and a screening region, to discuss $\gamma$-ray spectrum of mature pulsars detected by $Fermi$. By solving the Poisson equation with an assumed simple step-function distributions of the charge density in these two regions, the distribution of the electric field and the curvature radiation process of the accelerated particles can be calculated. In the our model, the properties of the phase-averaged spectrum can be completely specified by three gap parameters, i.e. the fractional gap size in the outer magnetosphere, the gap current in the primary region and the gap size ratio between the primary region and the total gap size. We discuss how these parameters affect the spectral properties. We argue that although the radiation mechanism in the outer gap is curvature radiation process, the observed gamma-ray spectrum can substantially deviate from the simple curvature spectrum because the overall spectrum consists of two components, i.e. the primary region and screening region. In some pulsars the radiation from the screening region is so strong that the photon index from 100MeV to several GeV can be as flat as $\sim 2$. We show the fitting fractional gap thickness of the canonical pulsars increases with the spin down age. We find that the total gap current is about 50~\% of the Goldreich-Julian value and the thickness of the screening region is a few percent of the total gap thickness. We also find that the predicted \gamma$-ray luminosity is less dependent on the spin down power ($L_{sd}$) for the pulsars with $L_{sd}\ga 10^{36}$~erg/s, while the $\gamma$-ray luminosity decreases with the spin down power for the pulsars with $L_{sd}\la 10^{36}$~erg/s.

Recently, the neutron star X-ray binary EXO 0748-676 underwent a transition to quiescence. We analyzed an XMM-Newton observation of this source in quiescence, where we fitted the spectrum with two different neutron-star atmosphere models. From the fits we constrained the allowed parameter space in the mass-radius diagram for this source for an assumed range of distances to the system. Comparing the results with different neutron-star equations of state, we constrained the distance to EXO 0748-676. We found that the EOS model 'SQM1' is rejected by the atmosphere model fits for the known distance, and the 'AP3' and 'MS1' is fully consistent with the known distance.

Magnetic reconnection, a fundamental plasma process associated with a rapid dissipation of magnetic energy, is believed to power many disruptive phenomena in laboratory plasma devices, the Earth magnetosphere, and the solar corona. Traditional reconnection research, geared towards these rather tenuous environments, has justifiably ignored the effects of radiation on the reconnection process. However, in many reconnecting systems in high-energy astrophysics (e.g., accretion-disk coronae, relativistic jets, magnetar flares) and, potentially, in powerful laser plasma and z-pinch experiments, the energy density is so high that radiation, in particular radiative cooling, may start to play an important role. This observation motivates the development of a theory of high-energy-density radiative magnetic reconnection. As a first step towards this goal, we present in this paper a simple Sweet--Parker-like theory of non-relativistic resistive-MHD reconnection with strong radiative cooling. First, we show how, in the absence of a guide magnetic field, intense cooling leads to a strong compression of the plasma in the reconnection layer, resulting in a higher reconnection rate. The compression ratio and the layer temperature are determined by the balance between ohmic heating and radiative cooling. The lower temperature in the radiatively-cooled layer leads to a higher Spitzer resistivity and hence to an extra enhancement of the reconnection rate. We then apply our general theory to several specific astrophysically important radiative processes (bremsstrahlung, cyclotron, and inverse-Compton) in the optically thin regime, for both the zero- and strong-guide-field cases. We derive specific expressions for key reconnection parameters, including the reconnection rate. We also discuss the limitations and conditions for applicability of our theory.

We discuss results from general-relativistic magnetohydrodynamical simulations of magnetised neutron star models (magnetars) including the effects of an elastic crust. The simulations reveal three distinct regimes: (a) a weak-field limit for magnetic field strengths $B<5\times10^{13}\,$G where purely crustal shear oscillations are recovered, (b) a strong-field limit $ B > 10^{15}$G where the magnetic field dominates the dynamics and the resulting quasi-periodic oscillations (QPOs) agree qualitatively with previous work, and (c) an intermediate regime, where purely crustal modes are damped rapidly with increasing magnetic field strength. Due to the presence of a solid crust a polar region exists where the standing-wave condition is significantly modified. As a result, strong QPOs are localised at a substantial angular distance from the pole. The boundary conditions at the base of the crust lead to a reversal in the order of the various families of QPOs. Pure crustal oscillations are strongly absorbed by the Alfv\'en continuum even for relatively low values of the poloidal magnetic field strength. This excludes torsional, axisymmetric shear modes of the crust as a viable interpretation of observed long-lived QPOs in giant flares of soft-gamma repeaters, if magnetic fields in magnetars are dominated by an axisymmetric dipolar component.

(abridged) In this paper, we express the relativistic propagational delay of light in the space-time of a binary system (commonly known as the "Shapiro delay") as a sum of harmonics of the orbital period of the system. We do this first for near-circular orbits as a natural expansion of an existing orbital model for low-eccentricity binary systems. The amplitudes of the 3rd and higher harmonics can be described by two new post-Keplerian (PK) parameters proportional to the amplitudes of the third and fourth harmonics (h_3, h_4). For high orbital inclinations we use a PK parameter proportional to the ratio of amplitudes of successive harmonics (sigma) instead of h_4. The new PK parameters are much less correlated with each other than r and s and provide a superior description of the constraints introduced by the Shapiro delay on the orbital inclination and the masses of the components of the binary (...). We extend the h_3,sigma parameterisation to eccentric binaries with high orbital inclinations. For some such binaries we can measure extra PK parameters and test general relativity using the Shapiro delay parameters. In this case we can use the measurement of h_3 as a test of general relativity. We show that this new test is not only more stringent than the r test, but it is even more stringent than the previous s test. Until now this new parametric test could only be derived statistically from an analysis of a probabilistic chi2 map.

The Advanced X-ray Timing Array (AXTAR) is a mission concept for X-ray timing of compact objects that combines very large collecting area, broadband spectral coverage, high time resolution, highly flexible scheduling, and an ability to respond promptly to time-critical targets of opportunity. It is optimized for submillisecond timing of bright Galactic X-ray sources in order to study phenomena at the natural time scales of neutron star surfaces and black hole event horizons, thus probing the physics of ultradense matter, strongly curved spacetimes, and intense magnetic fields. AXTAR's main instrument, the Large Area Timing Array (LATA) is a collimated instrument with 2-50 keV coverage and over 3 square meters effective area. The LATA is made up of an array of supermodules that house 2-mm thick silicon pixel detectors. AXTAR will provide a significant improvement in effective area (a factor of 7 at 4 keV and a factor of 36 at 30 keV) over the RXTE PCA. AXTAR will also carry a sensitive Sky Monitor (SM) that acts as a trigger for pointed observations of X-ray transients in addition to providing high duty cycle monitoring of the X-ray sky. We review the science goals and technical concept for AXTAR and present results from a preliminary mission design study.

As part of a survey for radio pulsars with the Parkes 64-m telescope we have discovered PSR J1622-4950, a pulsar with a 4.3-s rotation period. Follow-up observations show that the pulsar has the highest inferred surface magnetic field of the known radio pulsars (B ~ 3e14 G), exhibits significant timing noise and appears to have an inverted spectrum. Unlike the vast majority of the known pulsar population, PSR J1622-4950 appears to switch off for many hundreds of days and even in its on-state exhibits extreme variability in its flux density. Furthermore, the integrated pulse profile changes shape with epoch. All of these properties are remarkably similar to the only two magnetars previously known to emit radio pulsations. The position of PSR J1622-4950 is coincident with an X-ray source that, unlike the other radio pulsating magnetars, was found to be in quiescence. We conclude that our newly discovered pulsar is a magnetar - the first to be discovered via its radio emission.

Measuring the spin of Accreting Neutron Stars is important because it can provide constraints on the Equation of State of ultra-dense matter. Particularly crucial to our physical understanding is the discovery of sub-millisecond pulsars, because this will immediately rule out many proposed models for the ground state of dense matter. So far, it has been impossible to accomplish this because, for still unknown reasons, only a small amount of Accreting Neutron Stars exhibit coherent pulsations. An intriguing explanation for the lack of pulsations is that they form only on neutron stars accreting with a very low average mass accretion rate. I have searched pulsations in the faintest persistent X-ray source known to date and I found no evidence for pulsations. The implications for accretion theory are very stringent, clearly showing that our understanding of the pulse formation process is not complete. I discuss which sources are optimal to continue the search of sub-ms pulsars and which are the new constraints that theoretical models need to explain to provide a complete description of these systems

We present results from multi-epoch neutral hydrogen (HI) absorption observations of six bright pulsars with the Arecibo telescope. Moving through the interstellar medium (ISM) with transverse velocities of 10--150 AU/yr, these pulsars have swept across 1--200 AU over the course of our experiment, allowing us to probe the existence and properties of the tiny scale atomic structure (TSAS) in the cold neutral medium (CNM). While most of the observed pulsars show no significant change in their HI absorption spectra, we have identified at least two clear TSAS-induced opacity variations in the direction of B1929+10. These observations require strong spatial inhomogeneities in either the TSAS clouds' physical properties themselves or else in the clouds' galactic distribution. While TSAS is occasionally detected on spatial scales down to 10 AU, it is too rare to be characterized by a spectrum of turbulent CNM fluctuations on scales of 10-1000 AU, as previously suggested by some work. In the direction of B1929+10, an apparent correlation between TSAS and interstellar clouds inside the warm Local Bubble (LB) indicates that TSAS may be tracing the fragmentation of the LB wall via hydrodynamic instabilities. While similar fragmentation events occur frequently throughout the ISM, the warm medium surrounding these cold cloudlets induces a natural selection effect wherein small TSAS clouds evaporate quickly and are rare, while large clouds survive longer and become a general property of the ISM.

We report on the \textit{Fermi}-LAT observations of the Geminga pulsar, the second brightest non-variable GeV source in the $\gamma$-ray sky and the first example of a radio-quiet $\gamma$-ray pulsar. The observations cover one year, from the launch of the $Fermi$ satellite through 2009 June 15. A data sample of over 60,000 photons enabled us to build a timing solution based solely on $\gamma$ rays. Timing analysis shows two prominent peaks, separated by $\Delta \phi$ = 0.497 $\pm$ 0.004 in phase, which narrow with increasing energy. Pulsed $\gamma$ rays are observed beyond 18 GeV, precluding emission below 2.7 stellar radii because of magnetic absorption. The phase-averaged spectrum was fitted with a power law with exponential cut-off of spectral index $\Gamma$ = (1.30 $\pm$ 0.01 $\pm$ 0.04), cut-off energy $E_{0}$ = (2.46 $\pm$ 0.04 $\pm$ 0.17) GeV and an integral photon flux above 0.1 GeV of (4.14 $\pm$ 0.02 $\pm$ 0.32) $\times$ 10$^{-6}$ cm$^{-2}$ s$^{-1}$. The first uncertainties are statistical and the second are systematic. The phase-resolved spectroscopy shows a clear evolution of the spectral parameters, with the spectral index reaching a minimum value just before the leading peak and the cut-off energy having maxima around the peaks. Phase-resolved spectroscopy reveals that pulsar emission is present at all rotational phases. The spectral shape, broad pulse profile, and maximum photon energy favor the outer magnetospheric emission scenarios.

Seven years of pulse time-of-arrival measurements have been collected from observations of the young pulsar PSR B2334+61 using the Nanshan radio telescope of Urumqi Observatory. A phase-connected solution has been obtained over the whole data span, 2002 August to 2009 August. This includes a very large glitch that occurred between 2005 August 26 and September 8 (MJDs 53608 and 53621). The relative increase in rotational frequency for this glitch, $\Delta\nu_{g}/\nu~\sim~20.5\times10^{-6}$, is the largest ever seen. Although accounting for less than 1\% of the glitch, there were two well-defined exponential decay terms with time constants of 21 and 147 days respectively. There was also a large long-term increase in the spindown rate with $\Delta\dot\nu_p/\dot\nu \sim 0.011$ at the time of the glitch. A highly significant oscillation with a period of close to one year is seen in the post-glitch residuals. It is very unlikely that this can be accounted for by a pulsar position error or proper motion -- it appears to result from effects interior to the neutron star. Implications of these results for pulsar glitch models are discussed.

Using observations made with the XMM-Newton Observatory, we report the probable X-ray detection of the high-magnetic-field radio pulsar PSR J1734-3333. This pulsar has an inferred surface dipole magnetic field of B = 5.2e13 G, just below that of one anomalous X-ray pulsar (AXP). We find that the pulsar has an absorbed 0.5-2.0 keV flux of (5-15)e-15 erg/s/cm^2 and that its X-ray luminosity L_X is well below its spin down luminosity E_dot, with L_X < 0.1E_dot. No pulsations were detected in these data although our derived upper limit is unconstraining. Like most of the other high-B pulsars, PSR J1734-3333 is X-ray faint with no sign of magnetar activity. We collect and tabulate the properties of this and all other known high-B radio pulsars with measured X-ray luminosities or luminosity upper limits and plot L_X versus B for them all.

If the cosmological QCD phase transition is strongly first order and lasts sufficiently long, it generates a background of gravitational waves which may be detected via pulsar timing experiments. We estimate the amplitude and the spectral shape of such a background and we discuss its detectability prospects.

We present first comparisons between Light Element Primary Process (LEPP) abundances observed in ultra metal poor (UMP) stars and nucleosynthesis calculations based on long-time hydrodynamical simulations of core-collapse supernovae and their neutrino-driven wind. Observations indicate that r-process elements have at least two components: heavy r-process nuclei (A > 130) that are synthesized by rapid neutron capture in a yet unknown site and LEPP elements (mainly Sr, Y, Zr). We show that our neutrino-driven wind simulations can explain the observed LEPP pattern. We explore in detail the sensitivity of the calculated abundances to the electron fraction, which is a key nucleosynthesis parameter but poorly known due to uncertainties in neutrino interactions and transport. Our results show that the observed LEPP pattern can also be realized in proton-rich winds, which are obtained in the most recent supernova simulations. However, a small amount of neutron-rich matter from supernovae is necessary to account for the expected LEPP contribution to the solar system.

We calculate light curves produced by a hot spot of a rapidly rotating neutron star, assuming that the spot is perturbed by a core $r$-mode, which is destabilized by emitting gravitational waves. To calculate light curves, we take account of relativistic effects such as the Doppler boost due to the rapid rotation and light bending assuming the Schwarzschild metric around the neutron star. We assume that the core $r$-modes penetrate to the surface fluid ocean to have sufficiently large amplitudes to disturb the spot. For a $l'=m$ core $r$-mode, the oscillation frequency $\omega\approx2m\Omega/[l'(l'+1)]$ defined in the co-rotating frame of the star will be detected by a distant observer, where $l'$ and $m$ are respectively the spherical harmonic degree and the azimuthal wave number of the mode, and $\Omega$ is the spin frequency of the star. In a linear theory of oscillation, using a parameter $A$ we parametrize the mode amplitudes such that ${\rm max}\left(|\xi_\theta|,|\xi_\phi|\right)/R=A$ at the surface, where $\xi_\theta$ and $\xi_\phi$ are the $\theta$ and $\phi$ components of the displacement vector of the mode and $R$ is the radius of the star. For the $l'=m=2$ $r$-mode with $\omega=2\Omega/3$, we find that the fractional Fourier amplitudes at $\omega=2\Omega/3$ in light curves depend on the angular distance $\theta_s$ of the spot centre measured from the rotation axis and become comparable to or even larger than $A\sim0.001$ for small values of $\theta_s$.

The Fermi Large Area Telescope, in collaboration with several groups from the radio community, have had marvellous success at uncovering new gamma-ray millisecond pulsars (MSPs). In fact, MSPs now make up a sizable fraction of the total number of known gamma-ray pulsars. The MSP population is characterized by a variety of pulse profile shapes, peak separations, and radio-to-gamma phase lags, with some members exhibiting nearly phase-aligned radio and gamma-ray light curves (LCs). The MSPs' short spin periods underline the importance of including special relativistic effects in LC calculations, even for emission originating from near the stellar surface. We present results on modelling and classification of MSP LCs using standard pulsar model geometries.

The energy release due to neutralino WIMP self-annihilation in the thermalization volume inside a compact object is shown to be comparable to the energy needed to create a long-lived lump of strange quark matter, or strangelet, for WIMP masses above a few GeV. Since strange matter is the most stable state of matter, accretion of self-annihilating dark matter onto neutron stars provides a mechanism to seed compact objects with lumps of strange quark matter and this effect may trigger a conversion of most of the star into a strange star. Using an energy estimate based on the Fermi gas model combined with the MIT bag model for the long-lived strangelet, a new limit on the possible values of the WIMP mass can be set that is competitive with those from direct searches. Our limit is especially important for subdominant species of massive neutralinos.

X-ray observations of EXO~0748-676 during thermonuclear bursts revealed a set of narrow (\Delta \lambda /\lambda = 0.018) absorption lines that potentially originate from the stellar photosphere. The identification of these lines with particular atomic transitions led to the measurement of the surface gravitational redshift of the neutron star and to constraints on its mass and radius. However, the recent detection of 552 Hz oscillations at 15% rms amplitude revealed the spin frequency of the neutron star and brought into question the consistency of such a rapid spin with the narrow width of the absorption lines. Here, we calculate the amplitudes of burst oscillations and the width of absorption lines emerging from the surface of a rapidly rotating neutron star for a wide range of model parameters. We show that no combination of neutron-star and geometric parameters can simultaneously reproduce the narrowness of the absorption lines, the high amplitude of the oscillations, and the observed flux at the time the oscillations were detected. We, therefore, conclude that the observed absorption lines are unlikely to originate from the surface of this neutron star.

Many of the characteristic properties of the millisecond pulsars found in globular clusters are markedly different from those in the Galactic disc. We find that one such physical parameter is the surface magnetic field strength. Even though the average spin-periods do not differ much the average surface magnetic field is 2-5 times larger in the globular cluster pulsars. This effect could be apparent, arising due to one or more of several biases. Alternatively, if future observations confirm this effect to be real, then this could be interpreted as a preferential recycling of pulsars in tight binaries where the mass transfer takes place at high accretion rates.

Astronomy depends on ever increasing computing power. Processor clock-rates have plateaued, and increased performance is now appearing in the form of additional processor cores on a single chip. This poses significant challenges to the astronomy software community. Graphics Processing Units (GPUs), now capable of general-purpose computation, exemplify both the difficult learning-curve and the significant speedups exhibited by massively-parallel hardware architectures. We present a generalised approach to tackling this paradigm shift, based on the analysis of algorithms. We describe a small collection of foundation algorithms relevant to astronomy and explain how they may be used to ease the transition to massively-parallel computing architectures. We demonstrate the effectiveness of our approach by applying it to four well-known astronomy problems: Hogbom CLEAN, inverse ray-shooting for gravitational lensing, pulsar dedispersion and volume rendering. Algorithms with well-defined memory access patterns and high arithmetic intensity stand to receive the greatest performance boost from massively-parallel architectures, while those that involve a significant amount of decision-making may struggle to take advantage of the available processing power.

Almost all sources of high energy particles and photons are associated with jet phenomena. Prominent sources of such highly relativistic outflows are pulsar winds and Active Galactic Nuclei. The current understanding of these jets assumes diluted plasmas which are best described as kinetic phenomena. In this kinetic description particle acceleration to ultra-relativistic speeds can occur in completely unmagnetized and neutral plasmas through insetting effects of instabilities. Even though the morphology and nature of particle spectra are understood to a certain extent, the composition of the jets is not known yet. While Poynting-flux dominated jets are certainly composed of electron-positron plasmas, the understanding of the governing physics in AGN jets is mostly unclear. In this article we investigate how the constituting elements of an electron-positron-proton plasma behave differently under the variation of the fundamental mass-ratio m_p/m_e. We studied initially unmagnetized counterstreaming plasmas using fully relativistic three-dimensional particle-in-cell simulations to investigate the influence of the mass-ratio on particle acceleration and magnetic field generation in electron-positron-proton plasmas. We covered a range of mass-ratios m_p/m_e between 1 and 100 with a particle number composition of n_{p^+}/n_{e^+} of 1 in one stream, only protons are injected in the other, whereas electrons are present in both to guarantee charge neutrality in the simulation box. We find that with increasing proton mass the instability takes longer to develop and for mass-ratios > 20 the particles seem to be accelerated in two phases which can be accounted to the individual instabilities of the different species. This means that for high mass ratios the coupling between electrons/positrons and the heavier protons, which occurs in low mass-ratios, disappears.

We present a single pulse study of pulsar B1944+17, whose non-random nulls dominate nearly 70% of its pulses and usually occur at mode boundaries. When not in the null state, this pulsar displays four bright modes of emission, three of which exhibit drifting subpulses. B1944+17 displays a weak interpulse whose position relative to the main pulse we find to be frequency independent. Its emission is nearly 100% polarized, its polarization-angle traverse is very shallow and opposite in direction to that of the main pulse, and it nulls approximately two-thirds of the time. Geometric modeling indicates that this pulsar is a nearly aligned rotator whose alpha value is hardly 2 degrees--i.e., its magnetic axis is so closely aligned with its rotation axis that its sightline orbit remains within its conal beam. The star's nulls appear to be of two distinct types: those with lengths less than about 8 rotation periods appear to be pseudonulls--that is, produced by "empty" sightline traverses through the conal beam system; whereas the longer nulls appear to represent actual cessations of the pulsar's emission engine.

HESS J1632-478 is an extended and still unidentified TeV source in the galactic plane. In order to identify the source of the very high energy emission and to constrain its spectral energy distribution, we used a deep observation of the field obtained with XMM-Newton together with data from Molonglo, Spitzer and Fermi to detect counterparts at other wavelengths. The flux density emitted by HESS J1632-478 peaks at very high energies and is more than 20 times weaker at all other wavelengths probed. The source spectrum features two large prominent bumps with the synchrotron emission peaking in the ultraviolet and the external inverse Compton emission peaking in the TeV. HESS J1632-478 is an energetic pulsar wind nebula with an age of the order of 10^4 years. Its bolometric (mostly GeV-TeV) luminosity reaches 10% of the current pulsar spin down power. The synchrotron nebula has a size of 1 pc and contains an unresolved point-like X-ray source, probably the pulsar with its wind termination shock.

The launch of the Fermi Gamma-ray Space Telescope has heralded a new era in the study of gamma-ray pulsars. The population of confirmed gamma-ray pulsars has gone from 6-7 to more than 60, and the superb sensitivity of the Large Area Telescope (LAT) on Fermi has allowed the detailed study of their spectra and light curves. Twenty-four of these pulsars were discovered in blind searches of the gamma-ray data, and twenty-one of these are, at present, radio quiet, despite deep radio follow-up observations. In addition, millisecond pulsars have been confirmed as a class of gamma-ray emitters, both individually and collectively in globular clusters. Recently, radio searches in the direction of LAT sources with no likely counterparts have been highly productive, leading to the discovery of a large number of new millisecond pulsars. Taken together, these discoveries promise a great improvement in the understanding of the gamma-ray emission properties and Galactic population of pulsars. We summarize some of the results stemming from these newly-detected pulsars and their timing and multi-wavelength follow-up observations.

The majority of short gamma-ray bursts (SGRBs) are thought to originate from the merger of compact binary systems collapsing directly to form a black hole. However, it has been proposed that both SGRBs and long gamma-ray bursts (LGRBs) may, on rare occasions, form an unstable millisecond pulsar (magnetar) prior to final collapse. GRB 090515, detected by the Swift satellite was extremely short, with a T_90 of 0.036 +/- 0.016 s, and had a very low fluence of 2 x 10^-8 erg cm^-2 and faint optical afterglow. Despite this, the 0.3 - 10 keV flux in the first 200 s was the highest observed for a SGRB by the Swift X-ray Telescope (XRT). The X-ray light curve showed an unusual plateau and steep decay, becoming undetectable after ~500 s. This behaviour is similar to that observed in some long bursts proposed to have magnetars contributing to their emission. In this paper, we present the Swift observations of GRB 090515 and compare it to other gamma-ray bursts (GRBs) in the Swift sample. Additionally, we present optical observations from Gemini, which detected an afterglow of magnitude 26.4 +/- 0.1 at T+ 1.7 hours after the burst. We discuss potential causes of the unusual 0.3 - 10 keV emission and suggest it might be energy injection from an unstable millisecond pulsar. Using the duration and flux of the plateau of GRB 090515, we place constraints on the millisecond pulsar spin period and magnetic field.

Using the Chandra X-ray Observatory, we have pinpointed the location of a faint X-ray point source (CXOUJ182913.1-125113) and an associated diffuse nebula in the composite supernova remnant G18.95-1.1. These objects appear to be the long-sought pulsar and its wind nebula. The X-ray spectrum of the point source is best described by an absorbed powerlaw model with Gamma=1.6 and an N_H of ~1x10^(22) cm^(-2). This model predicts a relatively low unabsorbed X-ray luminosity of about L_X (0.5-8.0keV) = 4.1x10^(31)D_2^2 erg s^(-1), where D_2 is the distance in units of 2kpc. The best-fitted model of the diffuse nebula is a combination of thermal (kT = 0.48keV) and non-thermal (1.4 < Gamma < 1.9) emission. The unabsorbed X-ray luminosity of L_X = 5.4x10^(33)D_2^2 erg s^(-1) in the 0.5-8keV energy band seems to be largely dominated by the thermal component from the SNR, providing 87% of L_X in this band. No radio or X-ray pulsations have been reported for CXOUJ182913.1-125113. If we assume an age of ~5300yr for G18.95-1.1 and use the X-ray luminosity for the pulsar and the wind nebula together with the relationship between spin-down luminosity (via magnetic dipole radiation) and period, we estimate the pulsar's period to be P = 0.4s. Compared to other rotation-powered pulsars, a magnetic field of 2.2x10^(13)G is implied by its location in the P-Pdot diagram, a value which is close to that of the quantum critical field.

This paper reports the spectral and timing analyses of the quiescent low-mass X-ray binary U24 observed during five archived Chandra-ACIS exposures of the nearby globular cluster NGC 6397, for a total of 350 ksec. We find that the X-ray flux and the parameters of the hydrogen atmosphere spectral model are consistent with those previously published. Following the timing analysis, we find no evidence of short or long-term intensity variability. We also report the improved neutron star physical radius measurements, with statistical accuracy of the order of ~10%: R_ns = 8.9(+0.9)(-0.6) km for M_ns = 1.4 Msun. Alternatively, we provide the best-fit projected radius R_infinity= 11.9(+2.2)(-2.5)km, as seen by an observer at infinity. The best-fit effective temperature, kTeff = 80(+4)(-5) eV, is used to estimate the neutron star core temperature which falls in the range T_core = (3.0 - 9.8) x10 7 K, depending on the atmosphere model considered. This makes U24 the fourth most precisely measured neutron star radius among qLMXBs, after those in OmCen, in M13 and the qLMXB 47Tuc X7.

Since the last phase coherent timing solution of the nearby radio-quiet isolated neutron star RX J0720.4-3125 six new XMM-Newton and three Chandra observations were carried out. The phase coherent timing solutions from previous authors were performed without restricting to a fixed energy band. However, we recently showed that the phase residuals are energy dependent, and thus phase coherent solutions must be computed referring always to the same energy band. We updated the phase coherent timing solution for RX J0720.4-3125 by including the recent XMM-Newton EPIC-pn, MOS1, MOS2 and Chandra ACIS data in the energy range 400-1000~eV. Altogether these observations cover a time span of almost 10~yrs. A further timing solution was obtained including the ROSAT pointed data. In this case, observations cover a time span of $\approx$16~yrs. To illustrate the timing differences between the soft band (120-400~eV) and the hard band (400-1000~eV) a timing solution for the soft band is also presented and the results are verified using a $\mathrm{Z_{n}^{2}}$ test. In contrast to previous work, we obtain almost identical solutions whether or not we include the ROSAT or Chandra data. Thanks to the restriction to the hard band, the data points from EPIC-pn are in better agreement with those from MOS1, MOS2 and Chandra than in previous works. In general the phase residuals are still large and vary with time. In particular, the latest XMM-Newton and Chandra data show that the phase residuals have attained relatively large and negative values.

In the past ten years of regular operations, a new generation of Cherenkov telescopes have established binary systems as a new class of Very High Energy gamma-ray (VHE) emitters. Particle acceleration in these systems may occur either in an accretion-powered jet (microquasar) or in the shock between a pulsar wind and a stellar wind (wind-wind). This paper describes the phenomenology of the three VHE binaries PSR~B1259-63, LS 5039 and LS I +61 303. Two other objects may belong to this new class: HESS J0632+057 is a point-like variable VHE source whose multiwavelength behaviour resembles that of the other binaries, whereas Cyg X-1 is a well-known accreting system which may have been detected in VHE during a flaring episode. The paper concludes with a review of the latest searches for other binaries with Cherenkov telescopes, with special emphasis on Cyg X-3.

We show that equations governing pulsations of superfluid neutron stars can be splitted into two sets of weakly coupled equations, one describing the superfluid modes and another one -- the normal modes. The coupling parameter s is small, |s| ~ 0.01-0.05, for realistic equations of state. Already an approximation s=0 is sufficient to calculate the pulsation spectrum within the accuracy of a few percents. Our results indicate, in particular, that emission of gravitational waves from superfluid pulsation modes is suppressed in comparison to that from normal modes. The proposed approach allows to drastically simplify modeling of pulsations of superfluid neutron stars.

The nearby neutron star low-mass X-ray binary, Cen X-4, has been in a quiescent state since its last outburst in 1979. Typically, quiescent emission from these objects consists of thermal emission (presumably from the neutron star surface) with an additional hard power-law tail of unknown nature. Variability has been observed during quiescence in Cen X-4 on both timescales as short as hundreds of seconds and as long as years. However, the nature of this variability is still unknown. Early observations seemed to show it was all due to a variable hard X-ray tail. Here, we present new and archival observations that contradict this. The most recent Suzaku observation of Cen X-4 finds it in a historically low state, a factor of 4.4 fainter than the brightest quiescent observation. As the spectrum during the brightest observation was comprised of approximately 60% from the thermal component and 40% from the power-law component, such a large change cannot be explained by just power-law variability. Spectral fits with a variable thermal component fit the data well, while spectral fits allowing both the column density and the power-law to vary do not, leading to the conclusion that the thermal component must be variable. Interestingly, we also find that the thermal fraction remains consistent between all epochs, implying that the thermal and power-law fluxes vary by approximately the same amount. If the emitting area remains unchanged between observations, then the effective surface temperature must change. Alternatively, if the temperature remains constant, then the emitting area must change. The nature of this thermal variability is unclear, but may be explained by variable low-level accretion.

We report results of Swift observations for the high mass Be/X-ray binary system 1A 1118-615, during an outburst stage in January, 2009 and at a flaring stage in March, 2009. Using the epoch-folding method, we successfully detected a pulsed period of 407.69(2) sec in the outburst of January and of 407.26(1) sec after the flare detection in March. We find that the spectral detection for the source during outburst can be described by a blackbody model with a high temperature (kT ~ 1-3 keV) and a small radius (R ~ 1 km), indicating that the emission results from the polar cap of the neutron star. On the other hand, the spectra obtained after the outburst can further be described by adding an additional component with a lower temperature (kT ~ 0.1-0.2 keV) and a larger emission radius (R ~ 10-500 km), which indicates the emission from around the inner region of an accretion disk. We find that the thermal emission from the hot spot of the accreting neutron star dominates the radiation in outburst; the existence of both this X-ray contribution and the additional soft component suggest that the polar cap and the accretion disk emission might co-exist after the outburst. Because the two-blackbody signature at the flaring stage is a unique feature of 1A 1118-615, our spectral results may provide a new insight to interpret the X-ray emission for the accreting neutron star. The time separation between the three main outbursts of this system is ~17 years and it might be related to the orbital period. We derive and discuss the associated physical properties by assuming the elongated orbit for this specific Be/X-ray transient.

The ~800 yr-old pulsar J1846-0258 is a unique transition object between rotation-powered pulsars and magnetars: though behaving like a rotation-powered pulsar most of the time, in 2006 it exhibited a distinctly magnetar-like outburst accompanied by a large glitch. Here we present X-ray timing observations taken with the Rossi X-ray Timing Explorer over a 2.2-yr period after the X-ray outburst and glitch had recovered. We observe that the braking index of the pulsar, previously measured to be n=2.65+/-0.01, is now n=2.16+/-0.13, a decrease of 20+/-5%. We also note a persistent increase in the timing noise relative to the pre-outburst level. Despite the timing changes, a 2009 Chandra X-ray Observatory observation shows that the X-ray flux and spectrum of the pulsar and its wind nebula are consistent with the quiescent levels observed in 2000.

PSRB1055-52 is a middle-aged (~535 kyr) radio, X-ray, and gamma-ray pulsar showing X-ray thermal emission from the neutron star (NS) surface. A candidate optical counterpart to PSRB1055-52 was proposed by Mignani and coworkers based on Hubble Space Telescope (HST) observations performed in 1996, in one spectral band only. We report on HST observations of this field carried out in 2008, in four spectral bands. The astrometric and photometric analyses of these data confirm the identification of the proposed candidate as the pulsar's optical counterpart. Similarly to other middle-aged pulsars, its optical-UV spectrum can be described by the sum of a power-law (PLO) component, presumably emitted from the pulsar magnetosphere, and a Rayleigh-Jeans (RJ) component emitted from the NS surface. The spectral index of the PLO component, alpha_O=1.05+/-0.34, is larger than for other pulsars with optical counterparts. The RJ component, with the brightness temperature TO=(0.66+/-0.10) d_350**2 R_O,13**-2 MK (where d_350 and R_O,13 are the distance to the pulsar in units of 350 pc and the radius of the emitting area in units of 13 km), shows a factor of 4 excess with respect to the extrapolation of the X-ray thermal component into the UV-optical. This hints that the RJ component is emitted from a larger, colder area, and suggests that the distance to the pulsar is smaller than previously thought. From the absolute astrometry of the HST images we measured the pulsar coordinates with a position accuracy of 0.15". From the comparison with previous observations we measured the pulsar proper motion, mu = 42+/-5 mas/yr, which corresponds to a transverse velocity V_t = (70+/-8) d_350 km/s.

During the search for counterparts of very-high-energy gamma-ray sources, we serendipitously discovered large, extended, low surface brightness emission from PWNe around pulsars with the ages up to ~100 kyrs, a discovery made possible by the low and stable background of the Suzaku X-ray satellite. A systematic study of a sample of 8 of these PWNe, together with Chandra datasets, has revealed us that the nebulae keep expanding up to for ~100 kyrs, although time scale of the synchrotron X-ray emission is only ~60 yr for typical magnetic fields of 100 microG. Our result suggests that the accelerated electrons up to ~80 TeV can escape from the PWNe without losing most energies. Moreover, in order to explain the observed correlation between the X-ray size and the pulsar spindwon age, the magnetic field strength in the PWNe must decrease with time.

We argue that small pitch angle synchrotron emission provides an important dissipation mechanism which has to be taken into account in the models of formation of relativistic magnetized gamma-ray burst (GRB) outflows from the newborn black holes and/or magnetars. We show that if the GRB outflow is proton loaded, spectral energy distribution of this emission is expected to sharply peak in 0.1-1~MeV energy band. If the small pitch angle synchrotron emission efficiently cools relativistic particles of the outflow, its spectrum below the peak energy is a powerlaw with spectral index alpha ~ -1, close to the typical spectral index of the time-resolved GRB spectra. Otherwise, the low energy spectral index can be as hard as alpha ~ 0, as observed at the beginning of the GRB pulses. We make a conjecture that small pitch angle synchrotron emission from proton-loaded magnetized GRB outflow could significantly contribute to the Band component of the prompt emission of GRBs while electromagnetic cascade initiated by the protons could be responsible for the GeV component.

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

Set of neutron star observational results is used to test some selected equations of state of dense nuclear matter. The first observational result comes from the mass--baryon number relation for pulsar B of the double pulsar system J 0737--3039. The second one is based on the mass--radius relation coming from observation of the thermal radiation of the neutron star RX J 1856.35--3754. The third one follows the population analysis of isolated neutron star thermal radiation sources. The last one is the test of maximum mass. The equation of state of asymmetric nuclear matter is given by the parameterized form of the relativistic Brueckner-Hartree-Fock mean field, and we test selected parameterizations that represent fits of full relativistic mean field calculation. We show that only one of them is capable to pass the observational tests. This equation of state represents the first equation of state that is able to explain all the mentioned observational tests, especially the very accurate test given by the double pulsar even if no mass loss is assumed.

Pulsars provide a wealth of information about General Relativity, the equation of state of superdense matter, relativistic particle acceleration in high magnetic fields, the Galaxy's interstellar medium and magnetic field, stellar and binary evolution, celestial mechanics, planetary physics and even cosmology. The wide variety of physical applications currently being investigated through studies of radio pulsars rely on: (i) finding interesting objects to study via large-scale and targeted surveys; (ii) high-precision timing measurements which exploit their remarkable clock-like stability. We review current surveys and the principles of pulsar timing and highlight progress made in the rotating radio transients, intermittent pulsars, tests of relativity, understanding pulsar evolution, measuring neutron star masses and the pulsar timing array.

The arrival directions of ultrahigh energy extensive air showers (EAS) by Yakutsk, AGASA and P. Auger data are considered. It is found that the arrival directions of EAS with a deficit muons by Yakutsk data are not isotropy. Majority of these EAS form doublets which have maximum at side anticenter of Galaxy. It is shown that some EAS by data of Yakutsk, AGASA, P.Auger correlate with pulsars. These pulsars are situated near Input and Output Local arm Galaxy Orion. The majority of these pulsars have a short period rotate around of their axes. The problem of cosmic ray origin is discussed.

1E 1207.4-5209, the peculiar Central Compact object in the G296.5+10.0 supernova remnant, has been proposed to be an "anti-magnetar" - a young neutron star born with a weak dipole field. Accretion, possibly of supernova fallback material, has also been invoked to explain a large surface temperature anisotropy as well as the generation of peculiar cyclotron absorption features superimposed to its thermal spectrum. Interestingly enough, a faint optical/infrared source was proposed as a possible counterpart to 1E 1207.4-5209, but later questioned, based on coarse positional coincidence. Considering the large offset of 1E 1207.4-5209 with respect to the center of its host supernova remnant, the source should move at ~70 mas/yr. Thus, we tested the association by measuring the proper motion of the proposed optical counterpart. Using HST observations spanning 3.75 years, we computed a 3 sigma upper limit of 7 mas/yr. Absolute astrometry on the same HST data set also places the optical source significantly off the 99% confidence Chandra position. This allows us to safely rule out the association. Using the HST data set, coupled to ground-based observations collected at the ESO/VLT, we set the deepest limits ever obtained to the optical/infrared emission from 1E 1207.4-5209. By combining such limits to the constraints derived from X-ray timing, we rule out accretion as the source of the thermal anisotropy of the neutron star.

Pulsar Wind Nebulae (PWNe) act as calorimeters for the relativistic pair winds emanating from within the pulsar light cylinder. Their radiative dissipation in various wavebands is significantly different from that of their pulsar central engines: the broadband spectra of PWNe possess characteristics distinct from those of pulsars, thereby demanding a site of lepton acceleration remote from the pulsar magnetosphere. A principal candidate for this locale is the pulsar wind termination shock, a putatively highly-oblique, ultra-relativistic MHD discontinuity. This paper summarizes key characteristics of relativistic shock acceleration germane to PWNe, using predominantly Monte Carlo simulation techniques that compare well with semi-analytic solutions of the diffusion-convection equation. The array of potential spectral indices for the pair distribution function is explored, defining how these depend critically on the parameters of the turbulent plasma in the shock environs. Injection efficiencies into the acceleration process are also addressed. Informative constraints on the frequency of particle scattering and the level of field turbulence are identified using the multiwavelength observations of selected PWNe. These suggest that the termination shock can be comfortably invoked as a principal injector of energetic leptons into PWNe without resorting to unrealistic properties for the shock layer turbulence or MHD structure.

There have been suggestions that the abundance of Extremely Metal-Poor (EMP) stars can be reproduced by Hypernovae (HNe), not by normal supernovae (SNe). However, recently it was also suggested that if the innermost neutron-rich or proton-rich matter is ejected, the abundance patterns of ejected matter are changed, and normal SNe may also reproduce the observations of EMP stars. In this letter, we calculate explosive nucleosynthesis with various Ye and entropy, and investigate whether normal SNe with this innermost matter, which we call "hot-bubble" component, can reproduce the abundance of EMP stars. We find that neutron-rich (Ye = 0.45-0.49) and proton-rich (Ye = 0.51-0.55) matter can increase Zn/Fe and Co/Fe ratios as observed, but tend to overproduce other Fe-peak elements. In addition to it, we find that if slightly proton-rich matter with 0.50 <= Ye < 0.501 with s/kb ~ 15-40 is ejected as much as ~ 0.06 Msolar, even normal SNe can reproduce the abundance of EMP stars, though it requires fine-tuning of Ye. On the other hand, HNe can more easily reproduce the observations of EMP stars without fine-tuning. Our results imply that HNe are the most possible origin of the abundance pattern of EMP stars.

The merger of a binary system composed of a black hole and a neutron star may leave behind a torus of hot, dense matter orbiting around the black hole. While numerical-relativity simulations are necessary to simulate this process accurately, they are also computationally expensive and unable at present to cover the large space of possible parameters, which include the relative mass ratio, the stellar compactness and the black hole spin. To mitigate this and provide a first reasonable coverage of the space of parameters, we have developed a method for estimating the mass of the remnant torus from black hole-neutron star mergers. The toy model makes use of the relativistic affine model to describe the tidal deformations of an extended tri-axial ellipsoid orbiting around a Kerr black hole and measures the mass of the remnant torus by considering which of the fluid particles composing the star are on bound orbits at the time of the tidal disruption. We tune the toy model by using the results of fully general-relativistic simulations obtaining relative precisions of few a percent and use it to investigate the space of parameters extensively. In this way we find that the torus mass is largest for systems with highly spinning black holes, small stellar compactnesses and large mass ratios. As an example, tori as massive as $M_{b,\text{tor}} \simeq 1.33\,M_{\odot}$ can be produced for a very extended star with compactness $C\simeq 0.1$ inspiralling around a black hole with dimensionless spin $a=0.85$ and mass ratio $q\simeq 0.3$. However, for a more astrophysically reasonable mass ratio $q \simeq 0.14$ and a canonical value of the stellar compactness $C\simeq 0.145$, the toy model sets a considerably smaller upper limit of $M_{b,\text{tor}} \lesssim 0.34\,M_{\odot}$.

The spin of black holes in black hole-neutron star (BHNS) binaries can have a strong influence on the merger dynamics and the postmerger state; a wide variety of spin magnitudes and orientations are expected to occur in nature. In this paper, we report the first simulations in full general relativity of BHNS mergers with misaligned black hole spin. We vary the spin magnitude from a/m=0 to a/m=0.9 for aligned cases, and we vary the misalignment angle from 0 to 80 degrees for a/m=0.5. We restrict our study to 3:1 mass ratio systems and use a simple Gamma-law equation of state. We find that the misalignment angle has a strong effect on the mass of the postmerger accretion disk, but only for angles greater than ~ 40 degrees. Although the disk mass varies significantly with spin magnitude and misalignment angle, we find that all disks have very similar lifetimes ~ 100ms. Their thermal and rotational profiles are also very similar. For a misaligned merger, the disk is tilted with respect to the final black hole's spin axis. This will cause the disk to precess, but on a timescale longer than the accretion time. In all cases, we find promising setups for gamma-ray burst production: the disks are hot, thick, and hyperaccreting, and a baryon-clear region exists above the black hole.

A two-component hydrodynamic model is constructed of the global superfluid flow induced by two-component Ekman pumping during the recovery stage of a glitch. The model successfully accounts for the quasi-exponential recovery observed in pulsars like Vela and the "overshoot" observed in pulsars like the Crab. By fitting the model to high-resolution timing data, three important constitutive coefficients in bulk nuclear matter can be extracted: the shear viscosity, the mutual friction parameter, and the charged fluid fraction. The fitted coefficients for the Crab and Vela are compared with theoretical predictions for several equations of state, including the color-flavor locked and two-flavor color superconductor phases of quark matter.

U 2206+54 is a wind-fed high mass X-ray binary with a main-sequence donor star. The nature of its compact object was recently identified as a slow-pulsation magnetized neutron star. INTEGRAL/IBIS observations have a long-term hard X-ray monitoring of 4U 2206+54 and detected a hard X-ray outburst around 15 December 2005 combined with the RXTE/ASM data.The hard X-ray outburst had a double-flare feature with a duration of $\sim$ 2 days. The first flare showed a fast rise and long time decaying light curve about 15 hours with a peak luminosity of $\sim 4\times 10^{36}$ erg s$^{-1}$ from 1.5 -- 12 keV and a hard spectrum (only significantly seen above 5 keV). The second one had the mean hard X-ray luminosity of $1.3\times 10^{36}$ erg s$^{-1}$ from 20 -- 150 keV with a modulation period at $\sim 5550$ s which is the pulse period of the neutron star in 4U 2206+54; its hard X-ray spectrum from 20 -- 300 keV can be fitted with a broken power-law model with the photon indexes $\Gamma_1 \sim 2.3,\ \Gamma_2 \sim 3.3$, and the break energy is $E_b \sim 31$ keV or a bremsstrahlung model of $kT\sim 23$ keV. We suggest that the hard X-ray flare could be induced by suddenly enhanced accretion dense materials from stellar winds hitting the polar cap region of the neutron star. This hard X-ray outburst may be a link to supergiant fast X-ray transients though 4U 2206+54 has a different type of companion.

The deconfinement phase transition from hadronic matter to quark matter in the interior of compact stars is investigated. The hadronic phase is described in the framework of relativistic mean-field (RMF) theory, when also the scalar-isovector delta-meson effective field is taken into account. The MIT bag model for describing a quark phase is used. The changes of the parameters of phase transition caused by the presence of delta-meson field are explored. Finally, alterations in the integral and structural parameters of hybrid stars due to both deconfinement phase transition and inclusion of delta-meson field are discussed.

The cooling rate of young neutron stars gives direct insight into their internal makeup. Although the temperatures of several young neutron stars have been measured, until now a young neutron star has never been observed to decrease in temperature over time. We fit 9 years of archival Chandra ACIS spectra of the likely neutron star in the ~330 years old Cassiopeia A supernova remnant with our non-magnetic carbon atmosphere model. Our fits show a relative decline in the surface temperature by 4% (5.4 sigma, from 2.12+-0.01*10^6 K in 2000 to 2.04+-0.01*10^6 K in 2009) and observed flux (by 21%). Using a simple model for neutron star cooling, we show that this temperature decline could indicate that the neutron star became isothermal sometime between 1965 and 1980, and constrains some combinations of neutrino emission mechanisms and envelope compositions. However, the neutron star is likely to have become isothermal soon after formation, in which case the temperature history suggests episodes of additional heating or more rapid cooling. Observations over the next few years will allow us to test possible explanations for the temperature evolution.

Pulsars, or more generally rotation powered neutron stars, are excellent factories of antimatter in the Galaxy, in the form of pairs of electrons and positrons. Electrons are initially extracted from the surface of the star by the intense rotation induced electric fields and later transformed into electron-positron pairs through electromagnetic cascading. Observations of Pulsar Wind Nebulae (PWNe) show that cascades in the pulsar magnetosphere must ensure pair multiplicities of order $10^{4}-10^{5}$. These pairs finally end up as part of the relativistic magnetized wind emanating from the pulsar. The wind is slowed down, from its highly relativistic bulk motion, at a termination shock, which represents the reverse shock due to its interaction with the surrounding ejecta of the progenitor supernova. At the (relativistic) termination shock, acceleration of the pairs occurs, as part of the dissipation process, so that the cold wind is transformed into a plasma of relativistic non-thermal particles, plus a potential thermal component, which however has never been observed. As long as the pulsar wind is embedded in the supernova remnant these pairs are forced to escavate a bubble and lose energy adiabatically (because of the expansion) and radiatively (because of magnetic and radiation fields). We discuss here the observational constraints on the energy and number content of such pairs and discuss the scenarios that may allow for the pairs to escape in the interstellar medium and possibly contribute to the positron excess that has recently been detected by the PAMELA satellite. Special attention is dedicated to the case of Pulsar Bow Shock Nebulae. The pairs produced in these objects may be effectively carried out of the Supernova Remnant and released in the Interstellar Medium. As a result, Bow Shock Pulsar Wind Nebulae might be the main contributors to the positron excess in the Galaxy.

We have discovered an accreting black hole (BH) in a spectroscopically confirmed globular cluster (GC) in NGC 1399 through monitoring of its X-ray activity. The source, with a peak luminosity of L_x=2x10^39 ergs/s, reveals an order of magnitude change in the count rate within ~10 ks in a Chandra observation. The BH resides in a metal-rich [Fe/H]~0.2 globular cluster. After RZ2109 in NGC 4472 this is only the second black-hole X-ray source in a GC confirmed via rapid X-ray variability. Unlike RZ2109, the X-ray spectrum of this BH source did not change during the period of rapid variability. In addition to the short-term variability the source also exhibits long-term variability. After being bright for at least a decade since 1993 within a span of 2 years it became progressively fainter, and eventually undetectable, or marginally detectable, in deep Chandra and XMM observations. The source also became harder as it faded. The characteristics of the long term variability in itself provide sufficient evidence to identify the source as a BH. The long term decline in the luminosity of this object was likely not recognized in previous studies because the rapid variability within the bright epoch suppressed the average luminosity in that integration. The hardening of the spectrum accompanying the fading would also make this black hole source indistinguishable from an accreting neutron star in some epochs. Therefore some low mass X-ray binaries identified as NS accretors in snapshot studies of nearby galaxies may also be BHs. Thus the discovery of the second confirmed BH in an extragalactic GC through rapid variability at the very least suggests that accreting BHs in GCs are not exceedingly rare occurences.

Compact stars made of quark matter rather than confined hadronic matter, are expected to form a color superconductor. This superconductor ought to be threaded with rotational vortex lines within which the star's interior magnetic field is confined. The vortices (and thus magnetic flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at the surface of the star and the prolific production of thermal energy. In this Letter, we show that this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star enters the color superconducting phase, which can be used to estimate the density at which quarks deconfine. We find this density to be five times that of nuclear saturation.

Traditionally runaway stars are O and B type stars with large peculiar velocities.We want to extend this definition to young stars (up to ~50 Myr) of any spectral type and identify those present in the Hipparcos catalogue applying different selection criteria such as peculiar space velocities or peculiar one-dimensional velocities. Runaway stars are important to study the evolution of multiple star systems or star clusters as well as to identify origins of neutron stars. We compile distances, proper motions, spectral types, luminosity classes, V magnitudes and B-V colours and utilise evolutionary models from different authors to obtain star ages and study a sample of 7663 young Hipparcos stars within 3 kpc from the Sun. Radial velocities are obtained from the literature. We investigate the distributions of the peculiar spatial velocity, the peculiar radial velocity as well as the peculiar tangential velocity and its one-dimensional components and obtain runaway star probabilities for each star in the sample. In addition, we look for stars that are situated outside any OB association or OB cluster and the Galactic plane as well as stars of which the velocity vector points away from the median velocity vector of neighbouring stars or the surrounding local OB association/ cluster although the absolute velocity might be small. We find a total of 2547 runaway star candidates (with a contamination of normal Population I stars of 20 per cent at most). Thus, after subtraction of those 20 per cent, the runaway frequency among young stars is about 27 per cent. We compile a catalogue of runaway stars which will be available via VizieR.

The high-frequency quasi-periodic oscillations (HF QPOs) that appear in the X-ray fluxes of low-mass X-ray binaries remain an unexplained phenomenon. Among other ideas, it has been suggested that a non-linear resonance between two oscillation modes in an accretion disc orbiting either a black hole or a neutron star plays a role in exciting the observed modulation. Several possible resonances have been discussed. A particular model assumes resonances in which the disc-oscillation modes have the eigenfrequencies equal to the radial and vertical epicyclic frequencies of geodesic orbital motion. This model has been discussed for black hole microquasar sources as well as for a group of neutron star sources. Assuming several neutron (strange) star equations of state and Hartle-Thorne geometry of rotating stars, we briefly compare the frequencies expected from the model to those observed. Our comparison implies that the inferred neutron star radius "RNS" is larger than the related radius of the marginally stable circular orbit "rms" for nuclear matter equations of state and spin frequencies up to 800Hz. For the same range of spin and a strange star (MIT) equation of state, the inferrred radius RNS is roughly equal to rms. The Paczynski modulation mechanism considered within the model requires that RNS < rms. However, we find this condition to be fulfilled only for the strange matter equation of state, masses below one solar mass, and spin frequencies above 800Hz. This result most likely falsifies the postulation of the neutron star 3:2 resonant eigenfrequencies being equal to the frequencies of geodesic radial and vertical epicyclic modes. We suggest that the 3:2 epicyclic modes could stay among the possible choices only if a fairly non-geodesic accretion flow is assumed, or if a different modulation mechanism operates.

Being fast rotating objects, Isolated Neutron Stars (INSs) are natural targets for high-time resolution observations across the whole electromagnetic spectrum. With the number of objects detected at optical (plus ultraviolet and infrared) wavelengths now increased to 24, high-time resolution observations of INSs at these wavelengths are becoming more and more important. While classical rotation-powered radio pulsars, like the Crab and Vela pulsars, have been the first INSs studied at high-time resolution in the optical domain, observations performed in the last two decades have unveiled potential targets in other types of INSs which are not rotation powered, although their periodic variability is still related to the neutron star rotation. In this paper I review the current status of high-time resolution observations of INSs in the optical domain for different classes of objects: rotation-powered pulsars, magnetars, thermally emitting neutron stars, and rapid radio transients, I describe their timing properties, and I outline the scientific potentials of their optical timing studies.

We describe a Phased Array Feed (PAF) system, called Apertif, which will be installed in the Westerbork Synthesis Radio Telescope (WSRT). The aim of Apertif is, at frequencies from 1.0 to 1.7 GHz, to increase the instantaneous field of view of the WSRT 8 deg^2 and its observing bandwidth to 300 MHz with high spectral resolution. This system will turn the WSRT into an effective survey telescope with scientific applications ranging from deep surveys of the northern sky of HI and OH emission and polarised continuum to efficient searches for pulsars and transients. We present results obtained with a prototype PAF installed in one of the WSRT dishes. These results demonstrate that at decimetre wavelengths PAFs have excellent performance and that even for a single beam on the sky they outperform single feed radio dishes. PAFs turn radio telescopes into very effective survey instruments. Apertif is now fully funded and the community is invited to express their interest in using Apertif (this http URL )

Aims. We report results of an X-ray study of the supernova remnant (SNR) G344.7-0.1 and the point-like X-ray source located at the geometrical center of the SNR radio structure. Methods. The morphology and spectral properties of the remnant and the central X-ray point-like source were studied using data from the XMM-Newton and Chandra satellites. Archival radio data and infrared Spitzer observations at 8 and 24 $\mu$m were used to compare and study its multi-band properties at different wavelengths. Results. The XMM-Newton and Chandra observations reveal that the overall X-ray emission of G344.7-0.1 is extended and correlates very well with regions of bright radio and infrared emission. The X-ray spectrum is dominated by prominent atomic emission lines. These characteristics suggest that the X-ray emission originated in a thin thermal plasma, whose radiation is represented well by a plane-parallel shock plasma model (PSHOCK). Our study favors the scenario in which G344.7-0.1 is a 6 x 10^3 year old SNR expanding in a medium with a high density gradient and is most likely encountering a molecular cloud on the western side. In addition, we report the discovery of a soft point-like X-ray source located at the geometrical center of the radio SNR structure. The object presents some characteristics of the so-called compact central objects (CCO). However, its neutral hydrogen absorption column (N_{H}) is inconsistent with that of the SNR. Coincident with the position of the source, we found infrared and optical objects with typical early-K star characteristics. The X-ray source may be a foreground star or the CCO associated with the SNR. If this latter possibility were confirmed, the point-like source would be the farthest CCO detected so far and the eighth member of the new population of isolated and weakly magnetized neutron stars.

We investigate the thermonuclear bursting behaviour of IGR J17473-2721, an X-ray transient that in 2008 underwent a six month long outburst, starting (unusually) with an X-ray burst. We detected a total of 57 thermonuclear bursts throughout the outburst with AGILE, Swift, RXTE, and INTEGRAL. The wide range of inferred accretion rates (between <1% and about 20% of the Eddington accretion rate m-dot_Edd) spanned during the outburst allows us to study changes in the nuclear burning processes and to identify up to seven different phases. The burst rate increased gradually with the accretion rate until it dropped (at a persistent flux corresponding to about 15% of m-dot_Edd) a few days before the outburst peak, after which bursts were not detected for a month. As the persistent emission subsequently decreased, the bursting activity resumed with a much lower rate than during the outburst rise. This hysteresis may arise from the thermal effect of the accretion on the surface nuclear burning processes, and the timescale is roughly consistent with that expected for the neutron star crust thermal response. On the other hand, an undetected superburst, occurring within a data gap near the outburst peak, could have produced a similar quenching of burst activity.