[Abridged] Superluminous Supernovae (SN2006gy, SN2005gj, SN2005ap, SN2008fz, SN2003ma) have been a challenge to explain by standard models. We present an alternative scenario involving a quark-nova (QN), an explosive transition of the newly born neutron star to a quark star in which a second explosion (delayed) occurs inside the already expanding ejecta of a normal SN. The reheated SN ejecta can radiate at higher levels for longer periods of time primarily due to reduced adiabatic expansion losses, unlike the standard SN case. Our model is successfully applied to SN2006gy, SN2005gj, SN2005ap, SN2008fz, SN2003ma with encouraging fits to the lightcurves. There are four predictions in our model: (i) superluminous SNe optical lightcurves should show a double-hump with the SN hump at weaker magnitudes occurring days to weeks before the QN; (ii) Two shock breakouts should be observed vis-a-vis one for a normal SN. Depending on the time delay, this would manifest as two distinct spikes in the X-ray region or a broadening of the first spike for extremely short delays; (iii) The QN deposits heavy elements of mass number A> 130 at the base of the preceeding SN ejecta. These QN r-processed elements should be visible in the late spectrum (few days-weeks in case of strong ejecta mixing) of the superluminous SN; (iv) The QN yield will also contain lighter elements (Hydrogen and Helium). We expect the late spectra to include H_alpha emission lines that should be distinct in their velocity signature from standard H_alpha emission.

(Abridged) The hyperaccreting neutron star or magnetar disks cooled via neutrino emission can be a candidate of gamma-ray burst (GRB) central engines. The strong field $\geq10^{15}-10^{16}$ G of the magnetar can play a significant role in affecting the disk properties and even lead to the funnel accretion process. We investigate the effects of strong fields on the disks around magnetars, and discuss implications of such accreting magnetar systems for GRB and GRB-like events. We discuss quantum effects of the strong fields on the disk, and use the MHD conservation equations to describe the behavior of the disk flow coupled with a large scale field, which is generated by the star-disk interaction. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity, and change the electron fraction and degeneracy state significantly. A magnetized disk is always viscously stable outside the Alfv\'{e}n radius, but will be thermally unstable near the Alfv\'{e}n radius where the magnetic field plays a more important role in transferring the angular momentum and heating the disk than the viscous stress. The funnel accretion process will be only important for an extremely strong field, which creates a magnetosphere inside the Alfv\'{e}n radius and truncates the plane disk. Because of higher temperature and more concentrated neutrino emission of the magnetar surface ring-like belt region covered by funnel accretion, the neutrino annihilation rate from the accreting magnetars can be much higher than that from accreting neutron stars without fields. Furthermore, the neutrino annihilation mechanism and the magnetically-driven pulsar wind from the magnetar surface can work together to generate and feed an ultra-relativistic jet along the stellar magnetic poles.

The purpose of this review paper is to summarise the pulsar timing method, to provide an overview of recent research into the spin-down of pulsars over decadal timescales and to highlight the science that can be achieved using high-precision timing of millisecond pulsars.

On 2009 June 5, the Gamma-ray Burst Monitor (GBM) onboard the Fermi Gamma-ray Space Telescope triggered on two short, and relatively dim bursts with spectral properties similar to Soft Gamma Repeater (SGR) bursts. Independent localizations of the bursts by triangulation with the Konus-RF and with the Swift satellite, confirmed their origin from the same, previously unknown, source. The subsequent discovery of X-ray pulsations with the Rossi X-ray Timing Explorer (RXTE), confirmed the magnetar nature of the new source, SGR J0418+5729. We describe here the Fermi/GBM observations, the discovery and the localization of this new SGR, and our infrared and Chandra X-ray observations. We also present a detailed temporal and spectral study of the two GBM bursts. SGR J0418+5729 is the second source discovered in the same region of the sky in the last year, the other one being SGR J0501+4516. Both sources lie in the direction of the galactic anti-center and presumably at the nearby distance of ~2 kpc (assuming they reside in the Perseus arm of our galaxy). The near-threshold GBM detection of bursts from SGR J0418+5729 suggests that there may be more such dim SGRs throughout our galaxy, possibly exceeding the population of bright SGRs. Finally, using sample statistics, we conclude that the implications of the new SGR discovery on the number of observable active magnetars in our galaxy at any given time is <10, in agreement with our earlier estimates.

We investigate strongly interacting dense matter and neutron stars using a flavor-SU(3) approach based on a non-linear realization of chiral symmetry as well as a hadronic flavor-SU(2) parity-doublet model. We study chiral symmetry restoration and the equation of state of stellar matter and determine neutron star properties using different sets of degrees of freedom. Finally, we include quarks in the model approach. We show the resulting phase diagram as well as hybrid star solutions for this model.

We report an indication (3.22 sigma) of ~ 1860 Hz quasi-periodic oscillations from a neutron star low-mass X-ray binary 4U 1636-536. If confirmed, this will be by far the highest frequency feature observed from an accreting neutron star system, and hence could be very useful to understand such systems. This plausible timing feature was observed simultaneously with lower (~ 585 Hz) and upper (~ 904 Hz) kilohertz quasi-periodic oscillations. The two kilohertz quasi-periodic oscillation frequencies had the ratio of ~ 1.5, and the frequency of the alleged ~ 1860 Hz feature was close to the triple and the double of these frequencies. This can be useful to constrain the models of all the three features. In particular, the ~ 1860 Hz feature could be (1) from a new and heretofore unknown class of quasi-periodic oscillations, or (2) the first observed overtone of lower or upper kilohertz quasi-periodic oscillations. Finally we note that, although the relatively low significance of the ~ 1860 Hz feature argues for caution, even a 3.22 sigma feature at such a uniquely high frequency should be interesting enough to spur a systematic search in the archival data, as well as to scientifically motivate sufficiently large timing instruments for the next generation X-ray missions.

We argue that the recent groundbreaking discovery by Badenes et al. (2009) of a nearby (~50 pc) white dwarf-neutron star (or black hole) binary (SDSS 1257+5428) with a merger timescale ~500 Myr implies that such systems are common; we estimate that there are of order 10^6 in the Galaxy. Although subject to large uncertainties, the nominal derived merger rate is ~5 x 10^-4 per yr in the Milky Way, just ~3-6 and ~20-40 times less than the Type Ia and core-collapse supernova (SN) rates, respectively. This implies that the merger rate is ~0.5-1 x 10^4 per Gpc^3 per yr in the local universe, ~5000-10000 times more than the observed (beaming-uncorrected) long-duration gamma-ray burst (GRB) rate. We estimate the lower limit on the rate in the Galaxy to be >2.5 x 10^-5 per yr at 95% confidence. We briefly discuss the implications of this finding for the census of long- and short-duration GRBs and their progenitors, the frequency of tight binary companions to Type Ib/c SN progenitors, the origin of ultra-high energy cosmic rays (UHECRs), the formation of rapidly rotating neutron stars and ~2-3 M_sun black holes, the census of faint Ia-like SNe, as well as for upcoming and current transient surveys (e.g., LOSS, PTF, LSST), and for high- (LIGO) and low-frequency (LISA) gravitational wave searches.

The hadron-quark phase transition in the interior of compact stars is investigated, when the transition proceeds through a mixed phase. The hadronic phase is described in the framework of relativistic mean-field theory, when also the scalar-isovector delta-meson mean-field is taken into account. The changes of the parameters of phase transition caused by the presence of delta-meson field are explored. The results of calculation of structure of the mixed phase(Glendenning construction) are compared with the results of usual first-order phase transition (Maxwell construction).

It depends on the state of matter at supra-nuclear density to model pulsar's structure, which is unfortunately not certain due to the difficulties in physics. In cold quark matter at realistic baryon densities of compact stars (with an average value of $\sim 2-3\rho_0$), the interaction between quarks is so strong that they would condensate in position space to form quark-clusters. We argue that quarks in quark stars are grouped in clusters, then we apply two phenomenological models for quark stars, the polytropic model and Lennard-Jones model. Both of the two models have stiffer EoS, and larger maximum mass for quark stars (larger than 2 $M_\odot$). The gravitational energy releases during the AIQ process could explain the observed energy of three supergiant flares from soft gamma-ray repeaters ($\sim 10^{47}$ ergs).

We describe a focal plane array (FPA) system, called Apertif, that is being developed for the Westerbork Synthesis Radio Telescope (WSRT). The aim of Apertif is to increase the instantaneous field of view of the WSRT by a factor of 37 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 such as deep imaging surveys of the northern sky of HI and OH emission, of the polarised continuum and efficient searches for pulsars and transients. Such surveys will detect the HI in more than 100,000 galaxies out to z = 0.4, will allow to determine the detailed structure of the magnetic field of the Galaxy, and will discover more than 1,000 pulsars. We present experimental results obtained with a prototype FPA installed in one of the WSRT dishes. These results demonstrate that FPAs do have the performance that is required to make all these surveys possible.

A new accreting millisecond X-ray pulsar, IGR J17511-3057, was discovered in outburst on 2009 September 12 during the INTEGRAL Galactic bulge monitoring programme. To study the evolution of the source X-ray flux and spectral properties during the outburst, we requested a Swift monitoring of IGRJ17511-3057. In this paper we report on the results of the first two weeks of monitoring the source. The persistent emission of IGR J17511-3057 during the outburst is modeled well with an absorbed blackbody (kT~0.9 keV) and a power-law component (photon index~1-2), similar to what has been observed from other previously known millisecond pulsars. Swift also detected three type-I Xray bursts from this source. By assuming that the peak luminosity of these bursts is equal to the Eddington value for a pure helium type-I X-ray burst, we derived an upper limit to the source distance of ~10 kpc. The theoretical, expected recurrence time of the bursts according to the helium burst hypothesis is 0.2-0.9 days, in agreement with the observations.

In past decades a lot of progress has been made towards understanding the main s-process component that takes place in thermally pulsing Asymptotic Giant Branch (AGB) stars. During this process about half of the heavy elements, mainly between 90<=A<=209 are synthesized. Improvements were made in stellar modeling as well as in measuring relevant nuclear data for a better description of the main s process. The weak s process, which contributes to the production of lighter nuclei in the mass range 56<=A<=90 operates in massive stars (M>=8Msolar) and is much less understood. A better characterization of the weak s component would help disentangle the various contributions to element production in this region. For this purpose, a series of measurements of neutron-capture cross sections have been performed on medium-mass nuclei at the 3.7-MV Van de Graaff accelerator at FZK using the activation method. Also, neutron captures on abundant light elements with A<56 play an important role for s-process nucleosynthesis, since they act as neutron poisons and affect the stellar neutron balance. New results are presented for the (n,g) cross sections of 41K and 45Sc, and revisions are reported for a number of cross sections based on improved spectroscopic information.

3D Magnetohydrodynamic simulations show that when matter accretes onto neutron stars, in particular if the misalignment angle is small, it does not constantly fall at a fixed spot. Instead, the location at which matter reaches the star moves. These moving hot spots can be produced both during stable accretion, where matter falls near the magnetic poles of the star, and unstable accretion, characterized by the presence of several tongues of matter which fall on the star near the equator, due to Rayleigh-Taylor instabilities. Precise modeling with Monte Carlo simulations shows that those movements could be observed as high frequency Quasi Periodic Oscillations. We performed a number of new simulation runs with a much wider set of parameters, focusing on neutron stars with a small misalignment angle. In most cases we observe oscillations whose frequency is correlated with the mass accretion rate $\dot{M}$. Moreover, in some cases double QPOs appear, each of them showing the same correlation with $\dot{M}$.

The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August/September 2005, which reached a maximum X-ray flux of 400mCrab in the 5-100keV range. The outburst was observed by INTEGRAL and RXTE. The energy of the fundamental cyclotron line has been measured with INTEGRAL and RXTE at ~45keV. Flaring activity was observed during the rise to the peak of the outburst. RXTE observations during one of these flares found the energy of the fundamental cyclotron line shifted to a significantly higher position than during the rest of the outburst, where it remains constant. Also, the energy-dependent pulse profiles during the flare differ significantly from the rest of the outburst. These differences have been interpreted with the presence of magnetospheric instabilities at the onset of the accretion.
A decomposition method is applied to A 0535+26 pulse profiles. Basic assumptions of the method are that the asymmetry observed in the pulse profiles is caused by a distorted magnetic dipole field, and that the emission regions have axisymmetric beam patterns. Using pulse profiles obtained from RXTE observations, the contribution of the two emission regions has been disentangled. Constraints on geometry of the pulsar and a possible solution of the beam pattern are given. First results of the comparison of the reconstructed beam pattern with a geometrical model that includes relativistic light deflection are presented.

We report on the first Suzaku observation of IGR J16318-4848, the most extreme example of a new group of highly absorbed X-ray binaries that have recently been discovered by the International Gamma-Ray Astrophysics Laboratory INTEGRAL. The Suzaku observation was carried out between 2006 August 14 and 17, with a net exposure time of 97 ks.
The average X-ray spectrum of the source can be well described with a continuum model typical for neutron stars i.e., a strongly absorbed power law continuum with a photon index of 0.676(42) and an exponential cutoff at 20.5(6) keV. The absorbing column is 1.95(3)x10e24 cm-2. Consistent with earlier work, strong fluorescent emission lines of Fe Kalpha, Fe Kbeta, and Ni Kalpha are observed. Despite the large absorbing column, no Compton shoulder is seen in the lines, arguing for a non-spherical and inhomogeneous absorber.
Seen at an average 5-60 keV absorbed flux of 3.4x10e-10 erg cm-2 s-1, the source exhibits significant variability on timescales of hours.

An 8.8 solar mass electron-capture supernova (SN) was simulated in spherical symmetry consistently from collapse through explosion to nearly complete deleptonization of the forming neutron star. The evolution time of about 9 s is short because of nucleon-nucleon correlations in the neutrino opacities. After a brief phase of accretion-enhanced luminosities (~200 ms), luminosity equipartition among all species becomes almost perfect and the spectra of electron antineutrinos and muon/tau antineutrinos very similar. We discuss consequences for the neutrino-driven wind as a nucleosynthesis site and for flavor oscillations of SN neutrinos.

(Abridged) We studied an exceptional period of activity of the anomalous X-ray pulsar 1E 1547.0-5408 in January 2009, during which about 200 bursts were detected by INTEGRAL. The major activity episode happened when the source was outside the field of view of all the INTEGRAL instruments. But we were still able to study the properties of 84 bursts detected simultaneously by the anti-coincidence shield of the spectrometer SPI and by the detector of the imager ISGRI. We find that the luminosity of the 22 January 2009 bursts of 1E 1547.0-5408 was > 1e42 erg/s. This luminosity is comparable to that of the bursts of soft gamma repeaters (SGR) and is at least two orders of magnitude larger than the luminosity of the previously reported bursts from AXPs. Similarly to the SGR bursts, the brightest bursts of 1E 1547.0-5408 consist of a short spike of ~100 ms duration with a hard spectrum, followed by a softer extended tail of 1-10 s duration, which occasionally exhibits pulsations with the source spin period of ~2 s. The observation of AXP bursts with luminosities comparable to the one of SGR bursts strengthens the conjecture that AXPs and SGRs are different representatives of one and the same source type.

The GeV light curve of a pulsar is an important probe to detect acceleration regions in its magnetosphere. Motivated by the recent reports on the observations of pulsars by {\it Fermi} Large Area Telescope (LAT), we restudy the two-pole caustic model and revise it to investigate the properties of the light curves in the GeV band. In the revised model, although acceleration gaps can extend from the star surface to the light cylinder along near the last open field lines, the extension of the gaps along the azimuthal direction is limited because of photon-photon pair production process. In such gaps, high-energy photons are emitted uniformly and tangentially to the field lines but cannot be efficiently produced along these field lines where the distances to the null charge surface are larger than $\sim0.9$ times of the distance of the light cylinder, and the effective azimuth extension of the gaps is about $230^\circ$. The model is applied to the four pulsars Vela, PSR J1028-5819, PSR J0205+6449, and PSR J2021+3651 whose light curves obtained with {\it Fermi} have been recently released. The model is successful in reproducing the general feature of the light curves for the four pulsars, and the radial distances of the radio pulse for the four pulsars are estimated.

IGR J01583+6713 is a new transient source discovered by the INTEGRAL/IBIS hard X-ray surveys. Optical observations suggested that it should be a high mass X-ray binary, but its nature is still unclear. We used the INTEGRAL/IBIS data to study the nature of the transient hard X-ray source IGR J01583+6713 during its outburst. Temporal profiles and spectral properties of IGR J01583+6713 around its outburst on 2005 Dec 6 were obtained. During the outburst, the mean X-ray luminosity reached around 4\times 10^{35} erg s^{-1} in the energy range of 20 -- 100 keV. The continuum spectrum can be fitted by a bremsstrahlung model of kT\sim 45 keV or a power-law model of \Gamma\sim 2.1, and the electron resonant cyclotron absorption lines at \sim 35 keV and possible at \sim 63 keV were detected, suggesting that a magnetic neutron star of B\sim 4\times 10^{12} G is located in IGR J01583+6713. Timing analysis of the outburst light curve found a pulsation period of \sim 5.47 hr. Then we identify IGR J01583+6713 as a transient slow-pulsation X-ray pulsar with a magnetic neutron star. After the outburst, the flux of IGR J01583+6713 was decreasing and it could not be detected by IBIS after 2005 Dec 10. This property of a highly magnetized neutron star with a long pulsation period suggested that the transient X-ray pulsar IGR J01583+6713 could be a wind-fed accretion system.

We analyze pulsar fluxes at 1400 MHz ($S_{1400}$) and distances ($d$) extracted from the Parkes Multibeam Survey. Under the assumption that distribution of pulsar luminosities is distance-independent, we find that either (a) pulsar fluxes diminish with distance according to a non-standard power law, due, we suggest, to the presence of a component with $S_{1400} \propto 1/d$, or (b) that there are very significant (i.e. order of magnitude) errors in the dispersion-measure method for estimating pulsar distances. The former conclusion (a) supports a model for pulsar emission that has also successfully explained the frequency spectrum of the Crab and 8 other pulsars over 16 orders of magnitude of frequency, whilst alternative (b) would necessitate a radical re-evaluation of both the dispersion-measure method and current ideas about the distribution of free electrons within our Galaxy.

Determining radiation location observationally plays a very important role in testing the pulsar radiation models. One-photon pair production in the strong magnetic field, $\gamma-e^{+}e^{1}$, is one of the important physical processes in pulsar radiation mechanisms. Photons near pulsar surface with sufficient energy will be absorbed in the magnetosphere and the absorption optical depth for these GeV $\gamma$-ray photons is usually large. In this paper, we include the aberrational, rotational and general relativistic effects and calculate the $\gamma$-B optical depth for $\gamma$-ray photons. Then we use the derived optical depth to determine the radiation altitude lower bounds for photons with given energies. As a case study, we calculate the lower bounds of radiation altitudes of Crab pulsar for photons with energy from 5 GeV to 1 TeV.

We discuss different ways that neutron stars can generate gravitational waves, describe recent improvements in modelling the relevant scenarios in the context of improving detector sensitivity, and show how observations are beginning to test our understanding of fundamental physics. The main purpose of the discussion is to establish promising science goals for third-generation ground-based detectors, like the Einstein Telescope, and identify the various challenges that need to be met if we want to use gravitational-wave data to probe neutron star physics.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

This note gives a correction to the standard analysis of the delay pattern in the radio signals from a pulsar in a binary system; the same coordinate frame should be used for the transmission of the signal as for the motion of the pulsar in the field of its companion.

We report the results of a search for an emission line from radiatively decaying dark matter in the Chandra X-ray Observatory spectrum of the ultra-faint dwarf spheroidal galaxy Willman 1. 99% confidence line flux upper limits over the 0.4-7 keV Chandra bandpass are derived and mapped to an allowed region in the sterile neutrino mass-mixing angle plane that is consistent with recent constraints from Suzaku X-ray Observatory and Chandra observations of the Ursa Minor and Draco dwarf spheroidals. A significant excess to the continuum, detected by fitting the particle-background-subtracted source spectrum, indicates the presence of a narrow emission feature with energy 2.51 +/- 0.07 (0.11) keV and flux [3.53 +/- 1.95 (2.77)] X 10^(-6) photons/cm^2/s at 68% (90%) confidence. Interpreting this as an emission line from sterile neutrino radiative decay, we derive the corresponding allowed range of sterile neutrino mass and mixing angle using two approaches. The first assumes that dark matter is solely composed of sterile neutrinos, and the second relaxes that requirement. The detection is consistent with the sterile neutrino mass of 5.0 +/- 0.2 keV and a mixing angle in a narrow range for which neutrino oscillations can produce all of the dark matter and for which sterile neutrino emission from the cooling neutron stars can explain pulsar kicks, thus bolstering both the statistical and physical significance of our measurement.

The high accuracy and long time span of photoelectric observations allow them to be used for multifactor analyses and refining some of the "fine" photometric effects in the light curves of close binary systems. The results obtained can be subsequently interpreted in terms of the model of mass flow from the optical component of the close binary systems onto the accretion disk of the neutron star, which can explain satisfactorily the irregularities of the gaseous flow, the "hot spot", and the presence of splashes moving in individual Keplerian trajectories about the outer parts of the accretion disk of the neutron star Her X-1.

This paper reports new observations of pulsars B0943+10 and B1822--09 carried out with the Arecibo Observatory (AO) and the Giant Metrewave Radio Telescope (GMRT), respectively. Both stars exhibit two stable emission modes. We report the discovery in B0943+10 of a highly linearly polarized precursor component that occurs primarily in only one mode. This emission feature closely resembles B1822-09's precursor which also occurs brightly in only one mode. B0943+10's other mode is well known for its highly regular drifting subpulses that are apparently produced by a rotating carousel system of 20 beamlets. Similary, B1822-09 exhibits subpulse-modulation behavior only in the mode where its precursor is absent. We survey our 18 hours of B0943+10 observations and find that the sideband-modulation features, from which the carousel-rotation time can be directly determined, occur rarely--less than 5% of the time--but always indicating 20 beamlets. We present an analysis of B1822-09's modal modulation characteristics at 325-MHz and compare them in detail with B0943+10. The pulsar never seems to null, and we find a 43-rotation-period feature in the star's Q mode that modulates the interpulse as well as the conal features in the main pulse. We conclude that B1822-09 must have a nearly orthogonal geometry and that its carousel circulation time is long compared to the modal sub-sequences available in our observations, and the mainpulse/interpulse separation is almost exactly 180 degrees. We conclude the precursors for both stars are incompatible with core-cone emission. We assess the interesting suggestion by Dyks et al. that downward-going radiation produces B1822-09's precursor emission.

The one-electron Li-containing Coulomb systems of atomic type $(li, e)$ and molecular type $(li, li, e)$, $(li, \alpha, e)$ and $(li, p, e)$ are studied in the presence of a strong magnetic field $B \leq 10^{7}$ a.u. in the non-relativistic framework. They are considered at the Born-Oppenheimer approximation of zero order (infinitely massive centers) within the parallel configuration (molecular axis parallel to the magnetic field). The variational and Lagrange-mesh methods are employed in complement to each other. It is demonstrated that the molecular systems ${\rm LiH}^{3+}$, ${\rm LiHe}^{4+}$ and ${\rm Li}_{2}^{5+}$ can exist for sufficiently strong magnetic fields $B \gtrsim 10^{4}$ a.u. and that ${\rm Li}_{2}^{5+}$ can even be stable at magnetic fields typical of magnetars.

As the ground-based gravitational-wave telescopes LIGO, Virgo, and GEO 600 approach the era of first detections, we review the current knowledge of the coalescence rates and the mass and spin distributions of merging neutron-star and black-hole binaries. We emphasize the bi-directional connection between gravitational-wave astronomy and conventional astrophysics. Astrophysical input will make possible informed decisions about optimal detector configurations and search techniques. Meanwhile, rate upper limits, detected merger rates, and the distribution of masses and spins measured by gravitational-wave searches will constrain astrophysical parameters through comparisons with astrophysical models. Future developments necessary to the success of gravitational-wave astronomy are discussed.

The results from a systematic study of eleven pulsar wind nebulae with a torus structure observed with the Chandra X-ray observatory are presented. A significant observational correlation is found between the radius of the tori, r, and the spin-down luminosity of the pulsars, Edot. A logarithmic linear fit between the two parameters yields log r = (0.57 +- 0.22) log Edot -22.3 +- 8.0 with a correlation coefficient of 0.82, where the units of r and Edot are pc and ergs s^-1, respectively. The value obtained for the Edot dependency of r is consistent with a square root law, which is theoretically expected. This is the first observational evidence of this dependency, and provides a useful tool to estimate the spin-down energies of pulsars without direct detections of pulsation. Applications of this dependency to some other samples are also shown.

We observed the neutron star (NS) ultra-compact X-ray binary 4U0614+091 quasi-simultaneously in the radio band (VLA), mid-IR/IR (Spitzer/MIPS and IRAC), near-IR/optical (SMARTS), optical-UV (Swift/UVOT), soft and hard X-rays (Swift/XRT and RXTE). The source was steadily in its `hard state'. We detected the source in the whole range, for the first time in the radio band at 4.86 and 8.46 GHz and in the mid-IR at 24 um, up to 100 keV. The optically thick synchrotron spectrum of the jet is consistent with being flat from the radio to the mid-IR band. The flat jet spectrum breaks in the range (1-4)x10^(13) Hz to an optically-thin power-law synchrotron spectrum with spectral index ~-0.5. These observations allow us to estimate a lower limit on the jet radiative power of ~3x10^(32) erg/s and a total jet power Lj~10^(34) u_(0.05)^(-1) Ec^(0.53) erg/s (where Ec is the high-energy cutoff of the synchrotron spectrum in eV and u_(0.05) is the radiative efficiency in units of 0.05). The contemporaneous detection of the optically thin part of the compact jet and the X-ray tail above 30 keV allows us to assess the contribution of the jet to the hard X-ray tail by synchrotron self-Compton (SSC) processes. We conclude that, for realistic jet size, boosting, viewing angle and energy partition, the SSC emission alone, from the post-shock, accelerated, non-thermal population in the jet, is not a viable mechanism to explain the observed hard X-ray tail of the neutron star 4U0614+091.

We report on archival Hubble Space Telescope (HST) observations of the globular cluster M71 (NGC 6838). These observations, covering the core of the globular cluster, were performed by the Advanced Camera for Surveys (ACS) and the Wide Field Planetary Camera 2 (WFPC2). Inside the half-mass radius (r_h = 1.65') of M71, we find 33 candidate optical counterparts to 25 out of 29 Chandra X-ray sources while outside the half-mass radius, 6 possible optical counterparts to 4 X-ray sources are found. Based on the X-ray and optical properties of the identifications, we find 1 certain and 7 candidate cataclysmic variables (CVs). We also classify 2 and 12 X-ray sources as certain and potential chromospherically active binaries (ABs), respectively. The only star in the error circle of the known millisecond pulsar (MSP) is inconsistent with being the optical counterpart. The number of X-ray faint sources with L_x>4x10^{30} ergs/s (0.5-6.0 keV) found in M71 is higher than extrapolations from other clusters on the basis of either collision frequency or mass. Since the core density of M71 is relatively low, we suggest that those CVs and ABs are primordial in origin.

We investigate the possibilities that pulsars act as the lens in a gravitational microlensing event towards the galactic bulge or a spiral arm, especially the proper motion of pulsars is taken into consideration. Our estimation is based on expectant observations of FAST (Five hundred meter Aperture Spherical Telescope) and SKA (Square Kilometer Array), and two different models of pulsar distribution are used. We find that the lensing rate is > 1 event/decade, being high enough to search the real events, and then suggest that microlensing observations focusing on those pulsars identified by FAST or SKA in the near future are meaningful. Certainly, as an independent determination of pulsar mass, a future detection event of microlensing pulsars should be significant in the history of studying pulsar-like stars, and constrain the state of matter (either hadronic or quark matter) at supra-nuclear densities.

With the development of high contrast imaging techniques and infrared detectors, vast efforts have been devoted during the past decade to detect and characterize lighter, cooler and closer companions to nearby stars, and ultimately image new planetary systems. Complementary to other observing techniques (radial velocity, transit, micro-lensing, pulsar-timing), this approach has opened a new astrophysical window to study the physical properties and the formation mechanisms of brown dwarfs and planets. I here will briefly present the observing challenge, the different observing techniques, strategies and samples of current exoplanet imaging searches that have been selected in the context of the LAOG-Planet Imaging Surveys. I will finally describe the most recent results that led to the discovery of giant planets probably formed like the ones of our solar system, offering exciting and attractive perspectives for the future generation of deep imaging instruments.

We present here a phenomenological solid quark star pulsar model to interpret the observed thermal X-ray emission of isolated pulsars. The heat capacity for solid quark stars was found to be quite small, so that the residual internal stellar heat gained at the birth of the star could be dissipated in an extremely short timescale. However, the bombardment induced by backflowing plasma at the poles of solid quark stars would get the stars be reheated, so that long term soft X-ray emission can be sustained. Such a scenario could be used for those X-ray pulsars with significant magnetospheric activities, and their cooling processes would thus be established. Dim X-ray isolated neutron stars (XDINs) as well as compact central objects (CCOs) have been observed with dominant soft X-ray radiation combined with little magnetospheric manifestations. Such sources could be solid quark stars accreting in the propeller regime.

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.

Low mass neutron stars may be uniquely strong sources of gravitational waves (GW). The neutron star crust can support large deformations for low mass stars. This is because of the star's weaker gravity. We find maximum ellipticities $\epsilon$ (fractional difference in moments of inertia) that are 1000 times larger, and maximum quadrupole moments $Q_{22}$ over 100 times larger, for low mass stars than for 1.4 $M_\odot$ neutron stars. Indeed, we calculate that the crust can support an $\epsilon$ as large as 0.01 for a minimum mass neutron star. A 0.12 $M_\odot$ star, that is maximally strained and rotating at 100 Hz, will produce a characteristic gravitational wave strain of $h_0=4.2\times 10^{-24}$ at a distance of 1 kpc. The GW detector Advanced LIGO should be sensitive to such objects through out the Milky Way Galaxy.

We present an analysis of archival X-ray observations of the Type IIL supernova SN 1979C. We find that its X-ray luminosity is remarkably constant at (6.5 +/- 0.1) x 10^{38} erg/s. The high and steady luminosity is evidence for a stellar-mass (~5-10M_{sun}) black hole accreting material from either a supernova fallback disk or possibly from a binary companion. We find that the bright and steady X-ray light curve is not consistent with either a model for a supernova powered by magnetic braking of a rapidly rotating magnetar, or a model where the blast wave is expanding into a dense circumstellar wind.

Short duration lensing events tend to be generated by low-mass lenses or by lenses with high transverse velocities. Furthermore, for any given lens mass and speed, events of short duration are preferentially caused by nearby lenses (mesolenses) that can be studied in detail, or else by lenses so close to the source star that finite-source-size effects may be detected, yielding information about both the Einstein ring radius and the surface of the lensed star. Planets causing short-duration events may be in orbits with any orientation, and may have semimajor axes smaller than an AU, or they may reach the outer limits of their planetary systems, in the region corresponding to the Solar System's Oort Cloud. They can have masses larger than Jupiter's or smaller than Pluto's. Lensing therefore has a unique potential to expand our understanding of planetary systems. A particular advantage of lensing is that it can provide precision measurements of system parameters, including the masses of and projected separation between star and planet. We demonstrate how the parameters can be extracted and show that a great deal can be learned. For example, it is remarkable that the gravitational mass of nearby free-floating planet-mass lenses can be measured by complementing observations of a photometric event with deep images that detect the planet itself. A fraction of short events may be caused by high-velocity stars located within a kpc. Many high-velocity lenses are likely to be neutron stars that received large natal kicks. Other high-speed stars may be members of the halo population. Still others may be hypervelocity stars that have been ejected from the Galactic Center, or runaway stars escaped from close binaries, possibly including the progenitor binaries of Type Ia supernovae.

This paper defines the mathematical convention adopted to describe an electromagnetic wave and its polarisation state, as implemented in the PSRCHIVE software and represented in the PSRFITS definition. Contrast is made between the convention that has been widely accepted by pulsar astronomers and the IAU/IEEE definitions of the Stokes parameters. The former is adopted as the PSR/IEEE convention, and a set of useful parameters are presented for describing the differences between the PSR/IEEE standard and the conventions (either implicit or explicit) that form part of the design of observatory instrumentation. To aid in the empirical determination of instrumental convention parameters, well-calibrated average polarisation profiles of PSR J0304+1932 and PSR J0742-2822 are presented at radio wavelengths of approximately 10, 20, and 40 cm.

We report the results from archival XMM-Newton and INTEGRAL observations of the Supergiant Fast X-ray Transient (SFXT) IGR J18483-0311 in quiescence. The 18-60 keV hard X-ray behaviour of the source is presented here for the first time, it is characterized by a spectral shape ($\Gamma$ about 2.5) similar to that during outburst activity and the lowest measured luminosity level is about 10^34 erg s^-1. The 0.5-10 keV luminosity state, measured by XMM-Newton during the apastron passage, is about one order of magnitude lower and it is reasonably fitted by an absorbed black body model yielding parameters consistent with previous measurements. In addition, we find evidence (about 3.5 sigma significance) of an emission-like feature at about 3.3 keV in the quiescent 0.5-10 keV source spectrum. The absence of any known or found systematic effects, which could artificially introduce the observed feature, give us confidence about its non-instrumental nature. We show that its physical explanation in terms of atomic emission line appears unlikely and conversely we attempt to ascribe it to an electron cyclotron emission line which would imply a neutron star magnetic field of the order of about 3x10^11 G. Importantly, such direct estimation is in very good agreement with that independently inferred by us in the framework of accretion from a spherically symmetric stellar wind. If firmly confirmed by future longer X-ray observations, this would be the first detection ever of a cyclotron feature in the X-ray spectrum of a SFXT, with important implications on theoretical models.

We present a new method to estimate the numerical viscosity in simulations of astrophysical ob jects, which is based in the damping of fluid oscillations. We apply the method to general relativistic hydrodynamic simulations using a spherical coordinates. We perform 1D-spherical and 2D- axisymmetric simulations of radial oscillations in spherical systems. We calibrate first the method with simulations with added bulk viscosity and study the differences between different numerical schemes. We apply the method to radial oscillations of neutron stars and we conclude that the main source of numerical viscosity in this case is the surface of the star. We expect that this method could be useful to compute the resolution requirements and limitations of the numerical simulations in different astrophysical scenarios in the future.

The nine millisecond pulsars (MSPs) that have now been detected by Fermi-LAT are providing an excellent opportunity to probe the emission geometry of these ancient compact objects. As they are radio-loud, one may use the relative phase lags across wavebands to obtain constraints on the orientation, size, and location of their radio and gamma-ray beams. We model the gamma-ray light curves using geometric outer gap (OG) and two-pole caustic (TPC) models, in addition to a pair-starved polar cap (PSPC) model which incorporates the full General Relativistic E-field. We find that most MSP light curves are fit by OG and TPC models, while PSPC is more appropriate for two others. The light curves of the newest discovery, PSR J0034-0534, are best modeled using outer magnetosphere OG / TPC models of limited extension for both radio and gamma-ray beams. We model the radio emission of the other eight MSPs using a fixed-altitude conal model at lower altitude. We lastly deduce values for inclination and observer angles (alpha and zeta), as well as the flux correction factor, in each case.

Suzaku observed the region including HESS J1809-193, one of the TeV unidentified (unID) sources, and confirmed existence of the extended hard X-ray emission previously reported by ASCA, as well as hard X-ray emission from the pulsar PSR J1809-1917 in the region. One-dimensional profile of the diffuse emission is represented with a Gaussian model with the best-fit sigma of 7+-1 arcmin. The diffuse emission extends for at least 21 pc (at the 3sigma level, assuming the distance of 3.5 kpc), and has a hard spectrum with the photon index of Gamma ~1.7. The hard spectrum suggests the pulsar wind nebula (PWN) origin, which is also strengthened by the hard X-ray emission from PSR J1809-1917 itself. Thanks to the low background of Suzaku XIS, we were able to investigate spatial variation of the energy spectrum, but no systematic spectral change in the extended emission is found. These results imply that the X-ray emitting pulsar wind electrons can travel up to 21 pc from the pulsar without noticeable energy loss via synchrotron emission.

We present the results from a study of wide profile pulsars using high sensitivity multifrequency observations with the GMRT. Since the line of sight samples a large region of the polar cap in case of the wide profile pulsars, presence of simultaneous multiple drift regions is quite probable (as seen in PSR B0826-34 and PSR B0818-41). We solve the aliasing problem of PSR B0818-41 using the observed phase relationship of the drift regions, and determine its pattern rotation period P4 to be ~ 10s, which makes it the fastest known carousel. We find that, for all the pulsars showing drifting in multiple rings of emission, the drift pattern from the rings are phase locked. This can constraint the theoretical models of pulsar emission as it favors a pan magnetospeheric radiation mechanism.

The origin and composition of Ultra-High Energy Cosmic Ray Events (UHECRs) are under debate. Here we improve constraints on the source population(s) and compositions of UHECRs by accounting for UHECR deflections within existing Galactic magnetic field models (GMFs). We used Monte Carlo simulations for UHECRs detected by the Pierre Auger Observatory and AGASA in order to determine their outside-the-Galaxy arrival directions, and compared these with Galactic and extragalactic sources. The simulations, which used UHECR compositions from protons to Iron and seven models of the ordered GMF, include uncertainties in the GMF and a turbulent magnetic field. The correlation between UHECRs and nearby extended radiogalaxies (Nagar & Matulich 2008) remains valid, even strengthened, within several GMF models. Both the nearest radiogalaxy CenA, and the nearest radio-extended BL Lac, CGCG 413-019, are likely sources of multiple UHECRs. The correlation appears to be linked to the presence of the extended radio source rather than a tracer of an underlying population. It is possible, but unlikely, that all UHECRs originate in the nearby radiogalaxy CenA. For light UHECRs about a third of UHECRs can be "matched" to nearby galaxies with extended radio jets. The remaining UHECRs could also be explained as originating in extended radiogalaxies if one has at least one of: a large UHECR mean free path, a high cluster and/or intergalactic magnetic field, a heavy composition for two-thirds of the detected UHECRs. Several UHECRs have trajectories which pass close to Galactic magnetars and/or microquasars.
If extended radiogalaxies are, or trace, UHECR sources, the most consistent models for the ordered GMF are the BS-S and BS-A models; the GMF models of Sun et al. 2008 are acceptable if a dipole component is added.

An updated grid of stellar yields for low to intermediate-mass thermally-pulsing Asymptotic Giant Branch (AGB) stars are presented. The models cover a range in metallicity Z = 0.02, 0.008, 0.004, and 0.0001, and masses between 1Msun to 6Msun. New intermediate-mass Z = 0.0001 AGB models are also presented, along with a finer mass grid than used in previous studies. The yields are computed using an updated reaction rate network that includes the latest NeNa and MgAl proton capture rates, with the main result that between ~6 to 30 times less Na is produced by intermediate-mass models with hot bottom burning. In low-mass AGB models we investigate the effect on the production of light elements of including some partial mixing of protons into the intershell region during the deepest extent of each third dredge-up episode. The protons are captured by the abundant 12C to form a 13C pocket. The 13C pocket increases the yields of 19F, 23Na, the neutron-rich Mg and Si isotopes, 60Fe, and 31P. The increase in 31P is by factors of ~4 to 20, depending on the metallicity. Any structural changes caused by the addition of the 13C pocket into the He-intershell are ignored. However, the models considered are of low mass and any such feedback is likely to be small. Further study is required to test the accuracy of the yields from the partial-mixing models. For each mass and metallicity, the yields are presented in a tabular form suitable for use in galactic chemical evolution studies or for comparison to the composition of planetary nebulae.

Neutrino emission due to the pair breaking and formation processes in the bulk triplet superfluid in neutron stars is investigated with taking into account of anomalous weak interactions. We consider the problem in the BCS approximation discarding Fermi-liquid effects. In this approach we derive self-consistent equations for anomalous vector and axial-vector vertices of weak interactions taking into account the $^{3}P_{2}- ^{3}F_{2}$ mixing. Further we simplify the problem and consider the pure $^{3}P_{2}$ pairing with $m_{j}=0$, as is adopted in the minimal cooling paradigm. As was expected because of current conservation we have obtained a large suppression of the neutrino emissivity in the vector channel. More exactly, the neutrino emission through the vector channel vanishes in the nonrelativistic limit $V_F=0$. The axial channel is also found to be moderately suppressed. The total neutrino emissivity is suppressed by a factor of $1.9\times10^{-1}$ relative to original estimates using bare weak vertices.

Forty years have gone by since the optical identification of the Crab pulsar (PSR B0531+21), and 25 isolated neutron stars of different types have been now identified. Observations with ESO telescopes historically played a pivotal role in the optical studies of INSs, first with the 3.6m telescope and the New Technology Telescope (NTT), at the La Silla Observatory, and after 1998 with the Very Large Telescope (VLT), at the Paranal Observatory. In this review I summarise some of the most important results obtained in ten years of VLT observations of isolated neutron stars.

We investigate the strong electric current sheet that develops at the tip of the pulsar closed line region through time dependent three-dimensional numerical simulations of a rotating magnetic dipole. We show that curvature radiation from relativistic electrons and positrons in the current sheet may naturally account for several features of the high-energy pulsar emission. We obtain light curves and polarization profiles for the complete range of magnetic field inclination angles and observer orientations, and compare our results to recent observations from the Fermi gamma-ray telescope.

We present the new open-source spherically-symmetric general-relativistic (GR) hydrodynamics code GR1D. It is based on the Eulerian formulation of GR hydrodynamics (GRHD) put forth by Romero-Ibanez-Gourgoulhon and employs radial-gauge, polar-slicing coordinates in which the 3+1 equations simplify substantially. We discretize the GRHD equations with a finite-volume scheme, employing piecewise-parabolic reconstruction and an approximate Riemann solver. GR1D is intended for the simulation of stellar collapse to neutron stars and black holes and will also serve as a testbed for modeling technology to be incorporated in multi-D GR codes. Its GRHD part is coupled to various finite-temperature microphysical equations of state in tabulated form that we make available with GR1D. An approximate deleptonization scheme for the collapse phase and a neutrino-leakage/heating scheme for the postbounce epoch are included and described. We also derive the equations for effective rotation in 1D and implement them in GR1D. We present an array of standard test calculations and also show for the first time how simple analytic equations of state in combination with presupernova models from stellar evolutionary calculations can be used to study the qualitative aspects of black hole formation in failing rotating core-collapse supernovae. Going beyond this, we present a simulation with microphysical EOS and neutrino leakage/heating of a failing core-collapse supernova and black hole formation in a 40 solar mass star. We find favorably close agreement on the time of black hole formation and last stable protoneutron star mass with predictions from simulations with full Boltzmann neutrino radiation hydrodynamics.

Context: When a rotating neutron star loses angular momentum, the reduction in the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. This mechanism is named rotochemical heating and has previously been studied for neutron stars of nonsuperfluid matter, finding that they reach a quasi-steady configuration in which the rate at which the spin-down modifies the equilibrium concentrations is the same at which neutrino reactions restore the equilibrium. Aims: We describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons \Delta_n and protons \Delta_p on the rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. By this, we may determine the amplitude of the superfluid energy gaps for the neutron and protons needed to produce different thermal evolution of MSPs. Results: We find that the chemical imbalances in the star grow up to the threshold value \Delta_{thr}= min(\Delta_n+ 3\Delta_p, 3\Delta_n+\Delta_p), which is higher than the quasi-steady state achieved in absence of superfluidity. Therefore, the superfluid MSPs will take longer to reach the quasi-steady state than their nonsuperfluid counterparts, and they will have a higher a luminosity in this state, given by L_\gamma ~ (1-4) 10^{32}\Delta_{thr}/MeV \dot{P}_{-20}/P_{ms}^3 erg s^-1. We can explain the UV emission of the PSR J0437-4715 for 0.05 MeV<\Delta_{thr}<0.45 MeV.

Observations of core-collapse supernovae can be used to probe the quark-hadron phase transition at high baryon densities. Recent simulations of stellar core-collapse that include descriptions of quark matter predict a sharp burst of anti-nu_e several hundred milliseconds after the prompt nu_e neutronization burst. We study the observational signatures of the anti-nu_e burst at the current neutrino detectors IceCube and Super-Kamiokande, and find that a QCD phase transition in a Galactic core-collapse can clearly be detected, regardless of the neutrino oscillation scenario. The detection would constitute direct evidence of quark matter in the neutron star and would have significant implications for our understanding of supernovae and the state of matter at extreme conditions.

The Hall effect is an important nonlinear mechanism affecting the evolution of magnetic fields in neutron stars. Studies of the governing equation, both theoretical and numerical, have shown that the Hall effect proceeds in a turbulent cascade of energy from large to small scales. We investigate the small-scale Hall instability conjectured to exist from the linear stability analysis of Rheinhardt and Geppert. Identical linear stability analyses are performed to find a suitable background field to model Rheinhardt and Geppert's ideas. The nonlinear evolution of this field is then modelled using a three-dimensional pseudospectral numerical MHD code. Combined with the background field, energy was injected at the ten specific eigenmodes with the greatest positive eigenvalues as inferred by the linear stability analysis. Energy is transferred to different scales in the system, but not into small scales to any extent that could be interpreted as a Hall instability. Any instabilities are overwhelmed by a late-onset turbulent Hall cascade, initially avoided by the choice of background field, but soon generated by nonlinear interactions between the growing eigenmodes. The Hall cascade is shown here, and by several authors elsewhere, to be the dominant mechanism in this system.

The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers.

The aim of this work is to investigate in a physical and quantitative way the spectral evolution of bright Neutron Star Low-Mass X-ray Binaries (NS LMXBs), with special regard to the transient hard X-ray tails. We analyzed INTEGRAL data for five sources (GX 5-1, GX 349+2, GX 13+1, GX 3+1, GX 9+1) and built broad-band X-ray spectra from JEM-X1 and IBIS/ISGRI data. For each source, X-ray spectra from different states were fitted with the recently proposed model compTB. The spectra have been fit with a two-compTB model. In all cases the first compTB describes the dominant part of the spectrum that we interpret as thermal Comptonization of soft seed photons (< 1 keV), likely from the accretion disk, by a 3-5 keV corona. In all cases, this component does not evolve much in terms of Comptonization efficiency, with the system converging to thermal equilibrium for increasing accretion rate. The second compTB varies more dramatically spanning from bulk plus thermal Comptonization of blackbody seed photons to the blackbody emission alone. These seed photons (R < 12 km, kT_s > 1 keV), likely from the neutron star and the innermost part of the system, the Transition Layer, are Comptonized by matter in a converging flow. The presence and nature of this second compTB component (be it a pure blackbody or Comptonized) are related to the inner local accretion rate which can influence the transient behaviour of the hard tail: high values of accretion rates correspond to an efficient Bulk Comptonization process (bulk parameter delta > 0) while even higher values of accretion rates suppress the Comptonization, resulting in simple blackbody emission (delta=0).

Neutron stars are sensitive laboratories for testing general relativity, especially when considering deviations where velocities are relativistic and gravitational fields are strong. One such deviation is described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modified through the addition of the gravitational parity-violating Pontryagin density coupled to a field. This four-dimensional effective theory arises naturally both in perturbative and non-perturbative string theory, loop quantum gravity, and generic effective field theory expansions. We calculate here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neutron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of the metric to leading order, thus introducing modifications to the moment of inertia but not to the mass-radius relation. We show that an observational determination of the moment of inertia to an accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a constraint on the Chern-Simons coupling constant of \xi^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

We present results of optical spectroscopic observations of the mass donor star in SS 433 with Subaru and Gemini, with an aim to best constrain the mass of the compact object. Subaru/FOCAS observations were performed on October 6-8 and 10, 2007, covering the orbital phase of phi=0.96-0.26. We first calculate cross correlation function of these spectra with that of the reference star HD 9233 in the wavelength range of 4740-4840 Angstrom. This region is selected to avoid 'strong' absorption lines accompanied with contaminating emission components. The same analysis is applied to archive data of Gemini/GMOS taken at phi=0.84-0.30 by Hillwig & Gies (2008). From the Subaru and Gemini CCF results, the amplitude of radial velocity curve of the donor star is determined to be 58.3+/-3.8 km s-1 with a systemic velocity of 59.2+/-2.5 km s-1. Together with the radial velocity curve of the compact object, we derive the mass of the donor star and compact object to be M_O=12.4+/-1.9 M_sun and M_X=4.3+/-0.6 M_sun, respectively. We conclude, however, that these values should be taken as upper limits. From the analysis of the averaged absorption line profiles of strong lines and weak lines observed with Subaru, we find evidence for heating effects from the compact object. Using a simple model, we find that the true radial velocity amplitude of the donor star could be as low as 40+/-5 km s-1 in order to produce the observed absorption-line profiles. Taking into account the heating of the donor star may lower the derived masses to M_O=10.4 +2.3/-1.9 M_sun and M_X=2.5 +0.7/-0.6 M_sun. Our final constraint, 1.9 M_sun< M_X <4.9 M_sun, indicates that the compact object in SS 433 is most likely a low mass black hole, although the possibility of a massive neutron star cannot be firmly excluded.

Pulsars are known to power winds of relativistic particles that can produce bright nebulae by interacting with the surrounding medium. These pulsar wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some cases, also at TeV energies, but the lack of information in the gamma-ray band prevents from drawing a comprehensive multiwavelength picture of their phenomenology and emission mechanisms. Using data from the AGILE satellite, we detected the Vela pulsar wind nebula in the energy range from 100 MeV to 3 GeV. This result constrains the particle population responsible for the GeV emission, probing multivavelength PWN models, and establishes a class of gamma-ray emitters that could account for a fraction of the unidentified Galactic gamma-ray sources.

Forty years passed since the optical identification of the first isolated neutron star (INS), the Crab pulsar. 25 INSs have been now identified in the optical (O), near-ultraviolet (nUV), or near-infrared (nIR), hereafter UVOIR, including rotation-powered pulsars (RPPs), magnetars, and X-ray-dim INSs (XDINSs), while deep investigations have been carried out for compact central objects (CCOs), Rotating RAdio transients (RRATs), and high-magnetic field radio pulsars (HBRPs). In this review I describe the status of UVOIR observations of INSs, their emission properties, and I present the results from recent observations.

The Small Magellanic Cloud (SMC) Be/X-ray binary pulsar SXP6.85 = XTE J0103-728 underwent a large Type II outburst beginning on 2008 August 10. The source was consistently seen for the following 20 weeks (MJD = 54688 - 54830). We present X-ray timing and spectroscopic analysis of the source as part of our ongoing Rossi X-ray Timing Explorer (RXTE) monitoring campaign and INTEGRAL key programme monitoring the SMC and 47 Tuc. A comparison with the Optical Gravitational Lensing Experiment (OGLE) III light curve of the Be counterpart shows the X-ray outbursts from this source coincide with times of optical maximum. We attribute this to the circumstellar disk increasing in size, causing mass accretion onto the neutron star. Ground based IR photometry and H-alpha spectroscopy obtained during the outburst are used as a measure of the size of the circumstellar disk and lend support to this picture. In addition, folded RXTE light curves seem to indicate complex changes in the geometry of the accretion regions on the surface of the neutron star, which may be indicative of an inhomogeneous density distribution in the circumstellar material causing a variable accretion rate onto the neutron star. Finally, the assumed inclination of the system and H-alpha equivalent width measurements are used to make a simplistic estimate of the size of the circumstellar disk.

We suggest that nulling, mode changing and may be RRAT phenomena could be due to changes in the configuration of pulsar magnetosphere. It could be changes in the magnetosphere geometry or/and redistribution of currents flowing in the magnetosphere. This alters pulsar emission beam and depending on the orientation of the line of sight relative to the emission beam radiopulsar either changes mode of gets invisible causing 'nulls', what in an extreme case makes pulsar appearing as a RRAT. Changes in the magnetosphere configuration result in non-negligible changes of the pulsar spindown rate. We estimate that for different magnetosphere configurations of the same pulsar the dependence between the pulsar spindown rate W and the opening angle of the emission beam \alpha can be as strong as W\propto\alpha^4. We speculate about physical mechanisms which may cause reconfiguration of the magnetosphere.

We present results from two observations of the wind-accreting X-ray pulsar 4U 1907+09 using the Suzaku observatory. The broadband time-averaged spectrum allows us to examine the continuum emission of the source and the cyclotron resonance scattering feature at ~19 keV. Additionally, using the narrow CCD response of Suzaku near 6 keV allows us to study in detail the Fe K bandpass and to quantify the Fe K beta line for this source for the first time. The source is absorbed by fully-covering material along the line of sight with a column density of NH ~2e22 /cm^2, consistent with a wind accreting geometry, and a high Fe abundance (~3-4 x solar). Time and phase-resolved analyses allow us to study variations in the source spectrum. In particular, dips found in the 2006 observation which are consistent with earlier observations occur in the hard X-ray bandpass, implying a variation of the whole continuum rather than occultation by intervening material, while a dip near the end of the 2007 observation occurs mainly in the lower energies implying an increase in NH along the line of sight, perhaps indicating clumpiness in the stellar wind.

The thermal evolution of neutron stars (NSs) is investigated by coupling with the evolution of $\textit{r}$-mode instability that is described by a second order model.The heating effect due to shear viscous damping of the $\textit{r}$-modes enables us to understand the high temperature of two young pulsars (i.e., PSR B0531+21 and RX J0822-4300) in the framework of the simple $npe$ NS model, without superfluidity or exotic particles.Moreover, the light curves predicted by the model within an acceptable parameter regime may probably cover all of the young and middle-aged pulsars in the $\lg T_s^{\infty}-\lg t$ panel, and an artificially strong $p$ superfluidity invoked in some early works is not needed here. Additionally, by considering the radiative viscous damping of the $\textit{r}$-modes, a surprising extra cooling effect is found, which can even exceed the heating effect sometimes although plays an ignorable role in the thermal history.

We study non-linear effects of radiative viscosity of $npe$ matter in neutron stars for both direct Urca process and modified Urca process, and find that non-linear effects will decrease the ratio of radiative viscosity to bulk viscosity from 1.5 to 0.5 (for direct Urca process) and 0.375 (for modified Urca process). Which means that for small oscillations of neutron star, the large fraction of oscillation energy is emitted as neutrinos; but for large enough ones, bulk viscous dissipation dominates.

We present the results of our study of X-ray pulsar Cen X-3 using XMM-Newton observations. The light curve and the spectrum for this observations were built and Fe triplet within 6.5-7 keV region was detected. The geometric model of relativistic accretion disk for iron emission lines Fe I K alpha, Fe XXV and Fe XXVI in 6.4-7.0 keV region was applied. The values of disc inclination, inner and outer radii of the disc and mass of the central compact object (neutron star) were obtained. Intensity variations of these lines during orbital motion were also detected. The largest variation was detected for Fe I K alpha line, that agrees with the results of other authors. These results conform the model in which Fe I K alpha line forms in hot plasma of accretion disc and highly ionized iron lines form in outer regions of binary system. Probably the most interesting feature of Cen X-3 spectrum is Fe XXV triplet which was found by Iaria et al. (2005) from Chandra data analysis. We did not find this triplet in our analysis of XMM-Newton data and explain its absence by the insufficient energy resolution of XMM-Newton instruments.

We investigate the oscillation spectrum of rotating Newtonian neutron stars endowed with purely toroidal magnetic fields, using a time evolution code to evolve linear perturbations in the Cowling approximation. The background star is generated by numerically solving the MHD equilibrium equations and may be nonspherical by virtue of both rotation and magnetic effects; hence our perturbations and background are fully consistent. Whilst the background field is purely toroidal, the perturbed field is mixed poloidal-toroidal. From Fourier analysis of the perturbations we are able to identify a number of magnetically-restored Alfv\'en (or $a$-) modes. We show that in a rotating star pure inertial and $a$-modes are replaced by hybrid magneto-inertial modes, which reduce to $a$-modes in the nonrotating limit and inertial modes in the nonmagnetic limit. We show that the $r$-mode instability is suppressed by magnetic fields in sufficiently slowly rotating stars. In addition, we determine magnetic frequency shifts in the $f$-mode. We discuss the astrophysical relevance of our results, in particular for magnetar oscillations.

Recent X-ray observations by Fermi/GBM discovered a new torque reversal of 4U 1626-67 after 18 years of steady spinning down. Using Swift/BAT observations we were able to center this new torque reversal on Feb 4 2008, lasting approximately 150 days. From 2004 up to the end of 2007, the spin-down rate averaged at a mean rate of ~dnu/dt=-4.8e-13 Hz s-1 until the torque reversal reported here. Since then it has been following a steady spin-up at a mean rate of ~dnu/dt= 4e-13 Hz s-1. The properties of this torque reversal, as well as the lack of correlation between the X-ray flux and the torque applied to the neutron star before this transition, challenges our understanding of the physical mechanisms operating in this system.

The Vela pulsar is one of the most exciting gamma-ray sources and has been at the forefront of high-energy pulsar science since the detection of gamma-ray pulsations at the radio period by SAS-2 in 1975. With the unprecedented angular resolution, effective area, field of view, and timing resolution, in the GeV band, of the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope, the light curve of the Vela pulsar can be studied in greater detail than ever before. Using a timing solution derived solely from the LAT data, phase aligned with the radio emission, the spectrum of the Vela pulsar has been fit in intervals as small as 0.0016 in phase. Significant variation is seen in the cutoff energy and photon index across the light curve, strongly supporting curvature radiation as the source of the high-energy gamma-rays from the Vela pulsar.

We show that kicks generated by topological currents may be responsible for the large velocities seen in a number of pulsars. The majority of the kick builds up within the first second of the star's birth and generates a force about two orders of magnitude larger than a neutrino kick in the same temperature and magnetic field regime. Because of the nature of the topological currents the star's cooling is not affected until it reaches 10^9 K; thereafter the current replaces neutrino emission as the dominant cooling process. A requirement for the kick to occur is a suitably thin crust on the star; this leads us to speculate that pulsars with large kicks are quark stars and those with small kicks are neutron stars. If true this would be an elegant way to distinguish quark stars from neutron stars.

In the last couple of years the Magic air Cherenkov telescope has made significant contributions to very high energy $\gamma$-ray astronomy. These include the detection of the galactic binary system LSI +61 303 and the observation of pulsed emission from the Crab pulsar. Extragalactic objects like the famous FSRQ 3C 279 and the LBL S5 0716+714 have both been detected during optical high states, and the radio galaxy M87 could be observed during an unexpected strong $\gamma$-ray outburst. Given its low energy trigger threshold (~50 GeV) and fast repositioning time of less than 30s the Magic air Cherenkov telescope is particularly well suited for the observation of fast transient objects like AGN or GRBs. So far no GRB could be detected with Magic, however. In this paper we present selected highlights from recent MAGIC observations of extragalactic objects.

We demonstrate a method for quantitatively comparing gamma-ray pulsar light curves with magnetosphere beaming models. With the Fermi LAT providing many pulsar discoveries and high quality pulsar light curves for the brighter objects, such comparison allows greatly improved constraints on the emission zone geometry and the magnetospheric physics. Here we apply the method to Fermi LAT light curves of a set of bright pulsars known since EGRET or before. We test three approximate models for the magnetosphere structure and two popular schemes for the location of the emission zone, the Two Pole Caustic (TPC) model and the Outer Gap (OG) model. We find that OG models and relatively physical B fields approximating force-free dipole magnetospheres are preferred at high statistical significance. An application to the full LAT pulsar sample will allow us to follow the emission zone's evolution with pulsar spindown.

Using the Gamma ray Burst Monitor (GBM) on Fermi we are monitoring accreting pulsar systems. We use the rates from GBM's 12 NaI detectors in the 8-50 keV range to detect and monitor pulsations with periods between 0.5 and 1000 seconds. After discussing our analysis approach we present results for individual sources from the first year of monitoring. Updated figures for these and other sources are available at this http URL .

The Large Area Telescope (LAT), Fermi's main instrument, is providing a new view of the local energetic pulsar population. In addition to identifying a pulsar origin of a large fraction of the bright unidentified Galactic EGRET sources, the LAT results provide a great opportunity to study a sizable population of high-energy pulsars. Correlations of their physical properties, such as the trend of the luminosity versus the rotational energy loss rate, help identify global features of the gamma-ray pulsar population. Several lines of evidence, including the light curve and spectral features, suggest that gamma-ray emission from the brightest pulsars arises largely in the outer magnetosphere.

The Fermi Large Area Telescope (LAT) has provided the measurement of the high energy cosmic ray electrons plus positrons (CRE) spectrum with unprecedented accuracy form 20 GeV to 1 TeV. Recently this range has been extended down to ~ 7 GeV. The spectrum shows no prominent features and it is significantly harder than that inferred from several previous experiments. While the reported Fermi-LAT data alone may be interpreted in terms of a single (electron dominated) Galactic component, when combined with other complementary experimental results, specifically the CRE spectrum measured by H.E.S.S., and especially the positron fraction measured by PAMELA, an additional electron and positron component seems to be required. We show that the acceleration of electron-positron pairs in Galactic pulsars may offer a natural (though not unique) interpretation of all those results.

To constrain the giant pulse (GP) emission mechanism and test the model of Lyutikov (2007) of GP emission, we are carrying out a campaign of simultaneous observations of the Crab pulsar between gamma-rays (Fermi) and radio wavelengths. The correlation between times of arrival of radio GPs and high-energy photons, whether it exists or not, will allow us to choose between different origins of GP emission and further constrain the emission physics. Our foremost goal was testing whether radio GPs are due to changes in the coherence of the radio emission mechanism, variations in the pair creation rate in the pulsar magnetosphere, or changes in the beaming direction. Accomplishing this goal requires an enormous number of simultaneous radio GPs and gamma-photons. Thus, we organized a radio observations campaign using the 42-ft telescope at the Jodrell Bank Observatory (UK), the 140-ft telescope, and the 100-m Richard C. Byrd Green Bank Telescope (GBT) at the Green Bank Observatory (WV). While the observations with the two first ones are ongoing, here we present the preliminary results of 20 hrs of observations with the GBT at the high frequency of 8.9 GHz. These particular observations were aimed to probe the model of GP emission by Lyutikov (2007) which predicts that GPs at frequencies > 4 GHz should be accompanied by gamma -ray photons of energies of 1-100 GeV.

Since its launch in June 2008, the Large Area Telescope (LAT), onboard the \emph{Fermi} Gamma-ray Space Telescope, has greatly added to our understanding of gamma-ray pulsars. Its fine point spread function and large effective area, combined with the time-differencing method, make it the first gamma-ray instrument capable of discovering a new population of gamma-ray pulsars. We will present the recent discovery of the youngest (~4600 yr) radio-quiet gamma-ray pulsar discovered in a blind frequency search so far: PSR J1022-5746, a pulsar associated with an extended TeV source. We also present multiwavelength observations of the source, including X-ray observations.

The study of short-duration gamma-ray bursts (GRBs) has undergone a revolution in recent years thanks to the discovery of the first afterglows and host galaxies in May 2005. In this review we summarize our current knowledge of the galactic and sub-galactic environments of short GRBs, and the implications for the progenitor population. The most crucial results are: (i) some short GRBs occur in elliptical galaxies; (ii) the majority of short GRBs occur in star forming galaxies; (iii) the star forming hosts of short GRBs are distinct from the host galaxies of long GRBs in terms of star formation rates, luminosities, and metallicities, and instead appear to be drawn from the general field galaxy population; (iv) the physical offsets of short GRBs relative to their host galaxy centers are significantly larger than for long GRBs; (v) the observed offset distribution agrees well with predictions for the locations of NS-NS binary mergers; and (vi) unlike long GRBs, which tend to occur in the brightest regions of their hosts, the environments of short GRBs generally under-represent the light distribution of their host galaxies. Taken together, these observations suggest that short GRB progenitors have a wide age distribution and generally track stellar mass rather than star formation activity. These results are fully consistent with NS-NS binary mergers, but partial contribution from prompt or delayed magnetar formation is also consistent with the data.

We investigate the shear viscosity of neutron star matter in the presence of an antikaon condensate. The electron and muon number densities are reduced due to the appearance of a $K^-$ condensate in neutron star matter, whereas the proton number density increases. Consequently the shear viscosity due to scatterings of electrons and muons with themselves and protons is lowered compared to the case without the condensate. On the other hand, the contribution of proton-proton collisions to the proton shear viscosity through electromagnetic and strong interactions, becomes important and comparable to the neutron shear viscosity.

Supernova remnants (SNRs) are among the strongest candidates to explain the flux of cosmic rays below the knee around 10^15 eV. Pulsar wind nebulae (PWNe), synchrotron nebulae powered by the spin-down of energetic young pulsars, comprise one of the most populous VHE gamma-ray source classes. Gamma-ray studies in the GeV and TeV bands probe the nature (ions vs. electrons), production, and diffusion of high-energy particles in SNRs and PWNe. For sources that are visible across both the GeV and TeV bands, such as IC 443, the spatial and spectral distribution of gamma rays can be studied over an unprecedented energy range. This presentation will review recent VERITAS results, including studies of Cassiopeia A, IC 443, PSR J1930+1852, and the SNR G106.3+2.7/Boomerang region, and discuss prospects for complementary studies of SNRs and PWNe in the Fermi and VHE gamma-ray bands.

The VERITAS IACT observatory has carried out an extensive survey of the Cygnus region between 67 and 82 degrees in galactic longitude and between -1 and 4 degrees in galactic latitude. This region is a natural choice for a Very High Energy (VHE) gamma-ray survey in the Northern Hemisphere, as it contains a substantial number of potential VHE gamma-ray emitters such as supernova remnants, pulsar wind nebulae, high-mass X-ray binaries, and massive star clusters, in addition to a few previously detected VHE gamma-ray sources. It is also home to a number of GeV gamma-ray sources, including no less than four new high-significance sources detected in the first six months of Fermi data. The VERITAS survey, comprising more than 140 hours of observations, reaches an average VHE point-source flux sensitivity of better than 4% of the Crab Nebula flux at energies above 200 GeV. Here we report on preliminary results from this survey, including two source detections, and discuss the prospects for further studies that would exploit the joint coverage provided by VERITAS and Fermi data in this region.

We provide an analysis of timing irregularities observed for 366 pulsars. Observations were obtained using the 76-m Lovell radio telescope at the Jodrell Bank Observatory over the past 36 years. These data sets have allowed us to carry out the first large-scale analysis of pulsar timing noise over time scales of > 10yr, with multiple observing frequencies and for a large sample of pulsars. Our sample includes both normal and recycled pulsars. The timing residuals for the pulsars with the smallest characteristic ages are shown to be dominated by the recovery from glitch events, whereas the timing irregularities seen for older pulsars are quasi-periodic. We emphasise that previous models that explained timing residuals as a low-frequency noise process are not consistent with observation.

G54.1+0.3 is a young pulsar wind nebula (PWN), closely resembling the Crab, for which no thermal shell emission has been detected in X-rays. Recent Spitzer observations revealed an infrared (IR) shell containing a dozen point sources arranged in a ring-like structure, previously proposed to be young stellar objects. An extended knot of emission located in the NW part of the shell appears to be aligned with the pulsar's X-ray jet, suggesting a possible interaction with the shell material. Surprisingly, the IRS spectrum of the knot resembles the spectrum of freshly formed dust in Cas A, and is dominated by an unidentified dust emission feature at 21 microns. The spectra of the shell also contain various emission lines and show that some are significantly broadened, suggesting that they originate in rapidly expanding supernova (SN) ejecta. We present the first evidence that the PWN is driving shocks into expanding SN ejecta and we propose an alternative explanation for the origin of the IR emission in which the shell is composed entirely of SN ejecta. In this scenario, the freshly formed SN dust is being heated by early-type stars belonging to a cluster in which the SN exploded. Simple dust models show that this interpretation can give rise to the observed shell emission and the IR point sources.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

When it comes to identifying or to characterizing gamma-ray sources, X-ray observations are of paramount importance. Correlated X-and-gamma-ray flux variations are a powerful identification tool, if the gamma-ray source is unidentified, or an important diagnostic tool to understand the behaviour of an identified source. Moreover, X-ray observations of non-variable unidentified gamma-ray sources, both galactic and extragalactic, can unveil interesting candidate counterparts, narrowing down the search space and improving significantly the chances for a successful identification. Swift observations of Fermi gamma-ray-selected pulsar error boxes provide accurate positions of likely counterparts. This makes it possible both to confirm pulsations and to improve timing solution, opening up a new synergy between X and gamma ray observations.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

When it comes to identifying or to characterizing gamma-ray sources, X-ray observations are of paramount importance. Correlated X-and-gamma-ray flux variations are a powerful identification tool, if the gamma-ray source is unidentified, or an important diagnostic tool to understand the behaviour of an identified source. Moreover, X-ray observations of non-variable unidentified gamma-ray sources, both galactic and extragalactic, can unveil interesting candidate counterparts, narrowing down the search space and improving significantly the chances for a successful identification. Swift observations of Fermi gamma-ray-selected pulsar error boxes provide accurate positions of likely counterparts. This makes it possible both to confirm pulsations and to improve timing solution, opening up a new synergy between X and gamma ray observations.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

When it comes to identifying or to characterizing gamma-ray sources, X-ray observations are of paramount importance. Correlated X-and-gamma-ray flux variations are a powerful identification tool, if the gamma-ray source is unidentified, or an important diagnostic tool to understand the behaviour of an identified source. Moreover, X-ray observations of non-variable unidentified gamma-ray sources, both galactic and extragalactic, can unveil interesting candidate counterparts, narrowing down the search space and improving significantly the chances for a successful identification. Swift observations of Fermi gamma-ray-selected pulsar error boxes provide accurate positions of likely counterparts. This makes it possible both to confirm pulsations and to improve timing solution, opening up a new synergy between X and gamma ray observations.

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

When it comes to identifying or to characterizing gamma-ray sources, X-ray observations are of paramount importance. Correlated X-and-gamma-ray flux variations are a powerful identification tool, if the gamma-ray source is unidentified, or an important diagnostic tool to understand the behaviour of an identified source. Moreover, X-ray observations of non-variable unidentified gamma-ray sources, both galactic and extragalactic, can unveil interesting candidate counterparts, narrowing down the search space and improving significantly the chances for a successful identification. Swift observations of Fermi gamma-ray-selected pulsar error boxes provide accurate positions of likely counterparts. This makes it possible both to confirm pulsations and to improve timing solution, opening up a new synergy between X and gamma ray observations.

We identify two candidate magnetars in archival X-ray observations of HESS detected shell-type SNRs. X-ray point sources in CTB 37B coincident with HESS J1713-381 and in G353.6-0.7 coincident with HESS J1731-347 both have AXP-like spectra, much softer than those of ordinary, rotation powered pulsars, and no optical/IR counterparts. The spectrum of CXOU J171405.7-381031 in CTB 37B has a hard excess above 6 keV, which may be similar to such components seen in some AXPs. A new Chandra observation of this object reveals a highly significant pulsed signal at P = 3.82 s with pulsed fraction f_p = 0.31. Analysis of an XMM-Newton observation of the second candidate, XMMU J173203.3-344518 in G353.6-0.7, yields only marginal evidence for a 1 s period. If it is not a magnetar, then it could be a weakly magnetized central compact object (CCO). Considering that these HESS sources previously attributed to the SNR shells are possibly centrally peaked, we hypothesize that their pulsars may contribute to diffuse TeV emission. These identifications potentially double the number of magnetar/SNR associations in the Galaxy, and can be used to investigate the energetics and asymmetries of the supernovae that give rise to magnetars.

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to explore the mechanism by which it was ejected from the open cluster of its birth. Here we present recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic binary and physical parameters of the O star primary. We also present the first XMM-Newton observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place limits on the emission from a proposed compact companion in the system, and we rule out a quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved secondary companion is also not detected in our optical spectra of the binary, and it is difficult to conclude that a gravitational interaction could have ejected this runaway binary with a low mass optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent black hole with conditions unfavorable for accretion at the time of our observations.

The population of stellar black holes (SBHs) in the Galaxy and galaxies generally is poorly known in both number and distribution. SBHs are the fossil record of the massive stars in galaxy evolution and may have produced some (if not all) of the intermediate mass (\gsim100\Msun) black holes (IMBHs) and, in turn, the central supermassive black holes (SMBHs) in galactic nuclei. For the first time, a Galaxy-wide census of accreting black holes, and their more readily recognizable tracer population, accreting neutron stars (NSs), could be measured with a wide-field hard X-ray imaging survey and soft X-ray and optical/IR prompt followup -- as proposed for the EXIST mission.

We discuss possible radio bursts which can be generated during binary neutron stars mergers associated with short gamma-ray bursts. Low-frequency radio band appear to be advantageous due to the time delay of a radio signal propagating in the intergalactic medium, which makes it possible to use short gamma-ray burst alerts to search for specific regions on the sky by low-frequency radio instruments, especially LOFAR. The LOFAR sensitivity will allow significant detection of such bursts up to redshifts z~0.3 with a rate of ~25 events per year.

We consider the motion of charged particles in the vacuum magnetospheres of rotating neutron stars with a strong surface magnetic field, B>10^(12) G. The electrons and positrons falling into the magnetosphere or produced in it are shown to be captured by the force-free surface EB=0. Using the Dirac-Lorentz equation, we investigate the dynamics of particle capture and subsequent motion near the force-free surface. The particle energy far from the force-free surface has been found to be determined by the balance between the power of the forces of an accelerating electric field and the intensity of curvature radiation. When captured, the particles perform adiabatic oscillations along the magnetic field lines and simultaneously move along the force-free surface. We have found the oscillation parameters and trajectories of the captured particles. We have calculated the characteristic capture times and energy losses of the particles through the emission of both bremsstrahlung and curvature photons by them. The capture of particles is shown to lead to a monotonic increase in the thickness of the layer of charged plasma accumulating near the force-free surface. The time it takes for a vacuum magnetosphere to be filled with plasma has been estimated.

The gamma-ray survey of the sky by the Fermi Gamma-ray Space Telescope offers both opportunities and challenges for multiwavelength and multi-messenger studies. Gamma-ray bursts, pulsars, binary sources, flaring Active Galactic Nuclei, and Galactic transient sources are all phenomena that can best be studied with a wide variety of instruments simultaneously or contemporaneously. Identification of newly-discovered gamma-ray sources is largely a multiwavelength effort. From the gamma-ray side, a principal challenge is the latency from the time of an astrophysical event to the recognition of this event in the data. Obtaining quick and complete multiwavelength coverage of gamma-ray sources can be difficult both in terms of logistics and in terms of generating scientific interest. The Fermi LAT team continues to welcome cooperative efforts aimed at maximizing the scientific return from the mission through multiwavelength studies.

Geminga is the second brightest persistent source in the GeV gamma-ray sky. Discovered in 1975 by SAS-2 mission, it was identified as a pulsar only in the 90s, when ROSAT detected the 237 ms X-ray periodicity, that was later also found by EGRET in gamma rays. Even though Geminga has been one of the most intensively studied isolated neutron star during the last 30 years, its interest remains intact especially at gamma-ray energies, where instruments like the Large Area Telescope (LAT) aboard the Fermi mission will provide an unprecedented view of this pulsars. We will report on the preliminary results obtained on the analysis of the first year of observations. We have been able to do precise timing of Geminga using solely gamma rays, producing a timing solution and allowing a deep study of the evolution of the light curve with energy. We have also measured and studied the high-energy cutoff in the phase-averaged spectrum and produced a detailed study of the spectral evolution with phase.

Following an extremely interesting idea \cite{R1}, published long ago, the work function at the outer crust region of a strongly magnetized neutron star is obtained using relativistic version of Thomas-Fermi type model. In the present scenario, the work function becomes anisotropic; the longitudinal part is an increasing function of magnetic field strength, whereas the transverse part diverges. An approximate estimate of the electron density in the magnetosphere due to field emission and photo emission current, from the polar cap region are obtained.

Here we briefly report on first results of self-consistent simulation of non-stationary electron-positron cascades in the polar cap of pulsar using specially developed hybrid PIC/Monte-Carlo numerical code. We consider the case of Ruderman-Sutherland cascade -- when particles cannot be extracted from the surface of the neutrons star.