I review the evolutionary connection between low-mass X-ray binaries (LMXBs) and pulsars with binary companions (bPSRs) from a stellar binary evolution perspective. I focus on the evolution of stellar binaries with end-states consisting of a pulsar with a low-mass (<1.0 solar mass) companion, starting at the point the companion's progenitor first initiates mass transfer onto the neutron star. Whether this mass transfer is stable and the physics driving ongoing mass transfer partitions the phase space of the companions's initial mass and initial orbital period into five regions. The qualitative nature of the mass-transfer process and the binary's final end-state differ between systems in each region; four of these regions each produce a particular class of LMXBs. I compare the theoretical expectations to the populations of galactic field LMXBs with companion-mass constraints and field bPSRs. I show that the population of accreting millisecond pulsars are all identified with only two of the four LMXB classes and that these systems do not have readily identifiable progeny in the bPSR population. I discuss which sub-populations of bPSRs can be explained by binary evolution theory and those that currently are not. Finally I discuss some outstanding questions in this field.

We investigate the effects of neutrino-nucleus interactions (the nu-process) on the production of iron-peak elements in Population III core-collapse supernovae. The nu-process and the following proton and neutron capture reactions produce odd-Z iron-peak elements in complete and incomplete Si burning region. This reaction sequence enhances the abundances of Sc, Mn, and Co in the supernova ejecta. The supernova explosion models of 15 M_sol and 25 M_sol stars with the nu-process well reproduce the averaged Mn/Fe ratio observed in extremely metal-poor halo stars. In order to reproduce the observed Mn/Fe ratio, the total neutrino energy in the supernovae should be 3 - 9 x 10^{53} ergs. Stronger neutrino irradiation and other production sites are necessary to reproduce the observed Sc/Fe and Co/Fe ratios, although these ratios increase by the nu-process.

The INTEGRAL/IBIS survey was performed collecting all the GPS and GCDE data together with all the available public data . The second catalogue, published in 2006 by Bird et al., is dominated by detection of 113 X-ray binaries, with 38 being high-mass and 67 low-mass. In most systems the compact object is a neutron star, but the sample also contains 4 confirmed Black Holes and 6 LMXB black hole candidates (BHC). There are also, in additional, 6 tentative associations as BHCs based simply on spectral and timing properties. In the sample of 12 sources (BHC and tentatively associated BHC) there are 7 transient sources that went into outbursts during the INTEGRAL survey observations. We present here the monitoring of the time and spectral evolution of these 7 outbursts.

We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter when moving through the envelope of a giant companion.

Combined X-ray synchrotron and inverse-Compton gamma-ray observations of pulsar wind nebulae (PWN) may help to elucidate the processes of acceleration and energy loss in these systems. In particular, such observations provide constraints on the particle injection history and the magnetic field strength in these objects. The newly discovered TeV gamma-ray source HESS J1718-385 has been proposed as the likely PWN of the high spin-down luminosity pulsar PSR J1718-3825. The absence of previous sensitive X-ray measurements of this pulsar, and the unusual energy spectrum of the TeV source, motivated observations of this region with XMM-Newton. The data obtained reveal a hard spectrum X-ray source at the position of PSR 1718-3825 and evidence for diffuse emission in the vicinity of the pulsar. We derive limits on the keV emission from the centroid of HESS J1718-385 and discuss the implications of these findings for the PWN nature of this object.

We discuss three specific modes of accretion disks around rotating magnetized neutron stars which may explain the separations of the kilo Hertz quasi periodic oscillations (QPO) seen in low mass X-ray binaries. The existence of these modes requires that there be a maximum in the angular velocity of the accreting material, and that the fluid is in stable, nearly circular motion near this maximum rather than moving rapidly towards the star or out of the disk plane into funnel flows. It is presently not known if these conditions occur, but we are exploring this with 3D magnetohydrodynamic simulations and will report the results elsewhere. The first mode is a corotation mode which is radially trapped in the vicinity of the maximum of the disk rotation rate and is unstable. The second mode, relevant to relatively slowly rotating stars, is a magnetically driven eccentric ($m=1$) oscillation of the disk excited at a Lindblad radius in the vicinity of the maximum of the disk rotation. The third mode, relevant to rapidly rotating stars, is a magnetically coupled eccentric ($m=1$) and an axisymmetric ($m=0$) radial disk perturbation which has an inner Lindblad radius also in the vicinity of the maximum of the disk rotation. We suggest that the first mode is associated with the upper QPO frequency, $\nu_u$, the second with the lower QPO frequency, $\nu_\ell =\nu_u-\nu_*$, and the third with the lower QPO frequency, $\nu_\ell=\nu_u-\nu_*/2$, where $\nu_*$ is the star's rotation rate.

We investigate the number and type of pulsars that will be discovered with the low-frequency radio telescope LOFAR. We consider different search strategies for the Galaxy, for globular clusters and for galaxies other than our own. We show an all-sky Galactic survey can be optimally carried out by incoherently combining the LOFAR stations. In a 60-day all-sky Galactic survey LOFAR can find over a thousand pulsars, probing the local pulsar population to a very deep luminosity limit. For targets of smaller angular size, globular clusters and galaxies, the LOFAR stations can be combined coherently, making use of the full sensitivity. Searches of nearby northern-sky globular clusters can find large numbers of low luminosity millisecond pulsars (eg. over 10 new millisecond pulsars in a 10-hour observation of M15). If the pulsar population in nearby galaxies is similar to that of the Milky Way, a 10-hour observation could find the 10 brightest pulsars in M33, or pulsars in other galaxies out to a distance of 1.2Mpc.

The star HDE 226868 known as an optical counterpart of the black hole candidate Cyg X-1 has been observed in H_alpha region using spectrograph at Ondrejov 2-m telescope. The orbital parameters are determined from HeI-line by means of the author's method of Fourier disentangling. Preliminary results are also presented of disentangling the H_alpha-line into a P-Cyg profile of the (optical) primary and an emission profile of the circumstellar matter (and a telluric component).

We present a complete set of diagnostic tools aimed at reproducing synthetic non-thermal (synchrotron and/or Inverse Compton, IC) emissivity, integrated flux energy, polarization and spectral index simulated maps in comparison to observations. The time dependent relativistic magnetohydrodynamic (RMHD) equations are solved with a shock capturing code together with the evolution of the maximum particles energy. Applications to Pulsar Wind Nebulae (PWNe) are shown.

Context: Many thermally emitting isolated neutron stars have magnetic fields larger than 10^{13}G. A realistic cooling model should be reconsidered including the presence of high magnetic fields. Aims: We investigate the effects of anisotropic temperature distribution and Joule heating on the cooling of magnetized neutron stars. Methods: The 2D heat transfer equation with anisotropic thermal conductivity tensor and including all relevant neutrino emission processes is solved for realistic models of the neutron star interior and crust. Results: The presence of the magnetic field affects significantly the thermal surface distribution and the cooling history during both, the early neutrino cooling era and the late photon cooling era. Conclusions: There is a huge effect of the Joule heating on the thermal evolution of strongly magnetized neutron stars. Magnetic fields and Joule heating play a key role in maintaining magnetars warm for a long time. Moreover, this effect is also important for intermediate field neutron stars and should be considered in radio-quiet isolated neutron stars or high magnetic field radio-pulsars.

At the time the third EGRET catalog was published, unidentified sources accounted for a substantial fraction of the detections. To this day, the vast majority of these sources have not yet been associated with low-energy counterparts. In addition to known classes of gamma-ray emitters such as pulsars, supernova remnants, and blazars, a number of theoretically motivated candidate emitters have been suggested as the origin of these detections. We take a new approach to evaluate the plausibility of a Galactic population accounting for some or all of the unidentified EGRET sources. Rather than focusing on the properties of a specific candidate emitter, we constrain the abundance and spatial distribution of any objects of Galactic origin which may be among the EGRET unidentified sources by making the simple assumption that galaxies similar to the Milky Way host comparable populations of gamma-ray emitters. We find that it is highly improbable that the unidentified EGRET sources contain more than ~20 members of a Galactic halo population (e.g., annihilating dark matter clumps or intermediate mass black holes), but that current observations are consistent with all of these sources being Galactic objects if they reside entirely in the disk and bulge. However, upcoming observations by GLAST have the potential to exclude association of a large number of the unidentified sources with any Galactic source class. We discuss the additional constraints and new insights into the nature of Galactic gamma-ray emitting populations that GLAST is expected to provide.

Using a population synthesis approach, we compute the total merger rate in the local Universe for double neutron stars, double black holes, and black hole -- neutron star binaries. These compact binaries are the prime source candidates for gravitational-wave detection by LIGO and VIRGO. We account for mergers originating from field populations and from dense stellar clusters. For both populations we use the same treatment of stellar evolution. Our results indicate that the merger rates of double neutron stars and black hole -- neutron star binaries are strongly dominated by field populations, while merging black hole binaries are formed much more effectively in dense stellar clusters. The overall merger rate of double compact objects depends sensitively on the (largely unknown) initial mass fraction contained in dense clusters f_cl. For f_cl < 0.0001, the Advanced LIGO detection rate will be dominated by field populations of double neutron star mergers, with a small but significant number of detections ~20 yr^-1. However for a higher mass fraction in clusters, f_cl > 0.001, the detection rate will be dominated by numerous mergers of double black holes originating from dense clusters, and it will be considerably higher, ~25 - 300 yr^-1. In addition, we show that, once mergers of double black holes are detected, it is easy to differentiate between systems formed in the field and in dense clusters, since the chirp mass distributions are strikingly different. Finally, we point out that there may exist a population of merging black hole binaries in intergalactic space.

A key parameter of the polar-cap accelerator is a=j/(c rho_GJ), where j is the electric current and rho_GJ is the corotation charge density. The customary model assumed a=1+-epsilon and found epsilon << 1, steady acceleration of particles, and e+- creation. The fine-tuning of a to 1 is not, however, favored by the magnetosphere, and the accelerator with |a-1| > epsilon ~ 0.001 is qualitatively different. For 0<a<1-epsilon, a charge-separated outflow forms at a low voltage, and no e+- pairs are created. For a>1+epsilon, the accelerator is unstable and a high voltage is generated by electrostatic and induction effects. We argue that the pulsar activity is generated by an unsteady gap near the surface that forms if a>1 or a<0. The unsteady voltage implies unsteady rotation of the open magnetic tube and generation of Alfven waves, which are ducted along the tube. Their energy flux can directly feed radio emission as Alfven waves convert to electromagnetic waves.

The transient X-ray binary pulsar A0535+262 was observed with Suzaku on 2005 September 14 when the source was in the declining phase of the August-September minor outburst. The ~103 s X-ray pulse profile was strongly energy dependent, a double peaked profile at soft X-ray energy band (<3 keV) and a single peaked smooth profile at hard X-rays. The width of the primary dip is found to be increasing with energy. The broad-band energy spectrum of the pulsar is well described with a Negative and Positive power-law with EXponential (NPEX) continuum model along with a blackbody component for soft excess. A weak iron K_alpha emission line with an equivalent width ~25 eV was detected in the source spectrum. The blackbody component is found to be pulsating over the pulse phase implying the accretion column and/or the inner edge of the accretion disk may be the possible emission site of the soft excess in A0535+262. The higher value of the column density is believed to be the cause of the secondary dip at the soft X-ray energy band. The iron line equivalent width is found to be constant (within errors) over the pulse phase. However, a sinusoidal type of flux variation of iron emission line, in phase with the hard X-ray flux suggests that the inner accretion disk is the possible emission region of the iron fluorescence line.

Gamma-ray binaries, composed of a massive star and compact object, have been established as a new class of sources of very high energy (VHE) photons. The gamma-rays are produced by inverse Compton scattering of the stellar light by VHE electrons accelerated in the vicinity of the compact object. The VHE emission from LS 5039 displays an orbital modulation.
The inverse Compton spectrum depends on the angle between the incoming and outgoing photon in the electron rest frame. Since the angle at which an observer sees the star and electrons changes with the orbit, a phase dependence of the spectrum is expected. The phase-dependent spectrum of LS 5039 is calculated, assuming a continuous injection of electrons. The shape of the electron distribution depends on the injected power-law and on the magnetic field intensity.
Anisotropic scattering produces hard emission at inferior conjunction, when attenuation due to pair production of the VHE gamma-rays on star light is minimum. The computed lightcurve and spectra provide good fits to the HESS and EGRET observations, except at phases of maximum attenuation where pair cascade emission may be significant for HESS. Detailed predictions are made for a modulation in the GLAST energy range. The magnetic field intensity at periastron is 0.8+-0.2 G.
Anisotropic inverse Compton scattering plays a major role in LS 5039. The derived magnetic field intensity, injection energy and slope suggest a rotation-powered pulsar wind nebula. Gamma-ray binaries are promising sources to study the environment of pulsars on small scales.

We consider the electron-positron plasma generation processes in the magnetospheres of magnetars - neutron stars with strong surface magnetic fields, B = 10^(14) - 10^(15) G. We show that the photon splitting in a magnetic field, which is effective at large field strengths, does not lead to the suppression of plasma multiplication, but manifests itself in a high polarization of gamma-ray photons. A high magnetic field strength does not give rise to the second generation of particles produced by synchrotron photons. However, the density of the first-generation particles produced by curvature photons in the magnetospheres of magnetars can exceed the density of the same particles in the magnetospheres of ordinary radio pulsars. The plasma generation inefficiency can be attributed only to slow magnetar rotation, which causes the energy range of the produced particles to narrow. We have found a boundary in the P - Pdot diagram that defines the plasma generation threshold in a magnetar magnetosphere.

The evolution of single stars at low metallicity has attracted a large interest, while the effect of metallicity on binary evolution remains still relatively unexplored. We study the effect of metallicity on the number of binary systems that undergo different cases of mass transfer. We find that binaries at low metallicity are more likely to start transferring mass after the onset of central helium burning, often referred to as case C mass transfer. In other words, the donor star in a metal poor binary is more likely to have formed a massive CO core before the onset of mass transfer.
At solar metallicity the range of initial binary separations that result in case C evolution is very small for massive stars, because they do not expand much after the ignition of helium and because mass loss from the system by stellar winds causes the orbit to widen, preventing the primary star to fill its Roche lobe. This effect is likely to have important consequences for the metallicity dependence of the formation rate of various objects through binary evolution channels, such as long GRBs, double neutron stars and double white dwarfs.

We test statistically the hypothesis that radio pulsar glitches result from an avalanche process, in which angular momentum is transferred erratically from the flywheel-like superfluid in the star to the slowly decelerating, solid crust via spatially connected chains of local, impulsive, threshold-activated events, so that the system fluctuates around a self-organised critical state. Analysis of the glitch population (currently 285 events from 101 pulsars) demonstrates that the size distribution in individual pulsars is consistent with being scale invariant, as expected for an avalanche process. The waiting-time distribution is consistent with being exponential in seven out of nine pulsars where it can be measured reliably, after adjusting for observational limits on the minimum waiting time, as for a constant-rate Poisson process. PSR J0537$-$6910 and PSR J0835$-$4510 are the exceptions; their waiting-time distributions show evidence of quasiperiodicity. In each object, stationarity requires that the rate $\lambda$ equals $- \epsilon \dot{\nu} / <\Delta\nu>$, where $\dot{\nu}$ is the angular acceleration of the crust, $<\Delta\nu>$ is the mean glitch size, and $\epsilon\dot{\nu}$ is the relative angular acceleration of the crust and superfluid. There is no evidence that $\lambda$ changes monotonically with spin-down age. The rate distribution itself is fitted reasonably well by an exponential for $\lambda \geq 0.25 {\rm yr^{-1}}$. For $\lambda < 0.25 {\rm yr^{-1}}$, its exact form is unknown; the exponential overestimates the number of glitching pulsars observed at low $\lambda$, where the limited total observation time exercises a selection bias.

We investigate torsional Alfv\'en oscillations of relativistic stars with a global dipole magnetic field, via two-dimensional numerical simulations. We find that a) there exist two families of quasi-periodic oscillations (QPOs) with harmonics at integer multiples of the fundamental frequency, b) the lower-frequency QPO is related to the region of closed field lines, near the equator, while the higher-frequency QPO is generated near the magnetic axis, c) the QPOs are long-lived, d) for the chosen form of dipolar magnetic field, the frequency ratio of the lower to upper fundamental QPOs is ~0.6, independent of the equilibrium model or of the strength of the magnetic field, and e) within a representative sample of equations of state and of various magnetar masses, the Alfv\'en QPO frequencies are given by accurate empirical relations that depend only on the compactness of the star and on the magnetic field strength. The lower and upper QPOs can be interpreted as corresponding to the edges or turning points of an Alfv\'en continuum, according to the model proposed by Levin (2007). Several of the low-frequency QPOs observed in the X-ray tail of SGR 1806-20 can readily be identified with the Alfv\'en QPOs we compute. In particular, one could identify the 18Hz and 30Hz observed frequencies with the fundamental lower and upper QPOs, correspondingly, while the observed frequencies of 92Hz and 150Hz are then integer multiples of the fundamental upper QPO frequency (three times and five times, correspondingly). With this identification, we obtain an upper limit on the strength of magnetic field of SGR 1806-20 (if is dominated by a dipolar component) between ~3 and $7\times 10^{15}$G. Furthermore, we discuss the implications for the high-density EOS of compact stars. (Abridged)

We present an XMM-Newton detection of two low radio surface brightness SNRs, G85.4+0.7 and G85.9-0.6, discovered with the Canadian Galactic Plane Survey (CGPS). High-resolution XMM-Newton images revealing the morphology of the diffuse emission, as well as discrete point sources, are presented and correlated with radio and Chandra images. The new data also permit a spectroscopic analysis of the diffuse emission regions, and a spectroscopic and timing analysis of the point sources. Distances have been determined from HI and CO data to be 3.5 +/- 1.0 kpc for SNR G85.4+0.7 and 4.8 +/- 1.6 kpc for SNR G85.9-0.6. The SNR G85.4+0.7 is found to have a temperature of ~12-13 MK and a 0.5-2.5 keV luminosity of ~1-4 x 10^33 D(3.5)^2 erg/s (where D(3.5) is the distance in units of 3.5 kpc), with an electron density n_e of ~0.07-0.16(fD(3.5))^-1/2 cm^-3 (where f is the volume filling factor), and a shock age of ~9-49(fD(3.5))^1/2 kyr. The SNR G85.9-0.6 is found to have a temperature of ~15-19 MK and a 0.5-2.5 keV luminosity of ~1-4 x 10^34 D(4.8)^2 erg/s (where D(4.8) is the distance in units of 4.8 kpc), with an electron density n_e of ~0.04-0.10(fD(4.8))^-1/2 cm^-3 and a shock age of ~12-42(fD(4.8))^1/2 kyr. Based on the data presented here, none of the point sources appears to be the neutron star associated with either SNR.

We show that: (i) the long-term X-ray outburst light curve of the transient AXP XTE J1810-197 can be accounted for by a fallback disk that is evolving towards quiescence through a disk instability after having been heated by a soft gamma-ray burst, (ii) the spin-frequency evolution of this source in the same period can also be explained by the disk torque acting on the magnetosphere of the neutron star, (iii) most significantly, recently observed pulsed-radio emission from this source coincides with the epoch of minimum X-ray luminosity. This is natural in terms of a fallback disk model, as the accretion power becomes so low that it is not sufficient to suppress the beamed radio emission from XTE J1810-197.

We perform a timing analysis on RXTE data of the accreting millisecond pulsar XTE J1751-305 observed during the April 2002 outburst. After having corrected for Doppler effects on the pulse phases due to the orbital motion of the source, we performed a timing analysis on the phase delays, which gives, for the first time for this source, an estimate of the average spin frequency derivative <nu_dot> = (3.7 +/- 1.0)E-13 Hz/s. We discuss the torque resulting from the spin-up of the neutron star deriving a dynamical estimate of the mass accretion rate and comparing it with the one obtained from X-ray flux. Constraints on the distance to the source are discussed, leading to a lower limit of \sim 6.7 kpc.

We present the results of a detailed study of the high mass X-ray binary 2S 0114+650 made with the pointed instruments onboard the Rossi X-ray Timing Explorer. The spectral and temporal behaviour of this source was examined over the pulse, orbital, and super-orbital timescales, covering 2 cycles of the 30.7 d super-orbital modulation. Marginal evidence for variability of the power law photon index over the pulse period was identified, similar to that observed from other X-ray pulsars. If this variability is real it could be attributed to a varying viewing geometry of the accretion region with the spin of the neutron star. Variability of the neutral hydrogen column density over the orbital period was observed, which we attribute to the line of sight motion of the neutron star through the dense circumstellar environment. A hardening of the spectrum was observed during the orbital maximum, which we attribute to absorption effects as the neutron star undergoes partial eclipse. No significant variability of the column density was observed over the super-orbital period, indicating that variable obscuration by a precessing warp in an accretion disc is not the mechanism behind the super-orbital modulation. In contrast, a significant softening of the spectrum was observed during the super-orbital minimum. We conclude that the observed super-orbital modulation is tied to variability in the mass accretion rate due to some as yet unidentified mechanism.

We model the interaction between the wind from a newly formed rapidly rotating magnetar and the surrounding progenitor. In the first few seconds after core collapse the magnetar inflates a bubble of plasma and magnetic fields behind the supernova shock, which expands asymmetrically because of the pinching effect of the toroidal magnetic field, as in PWNe, even if the host star is spherically symmetric. The degree of asymmetry depends on the ratio of the magnetic energy to the total energy in the bubble. We assume that the wind by newly formed magnetars inflating these bubbles is more magnetized than for PWNe. We show that for a magnetic to total power supplied by the central magnetar $\sim 0.1$ the bubble expands relatively spherically while for values greater than 0.3, most of the pressure in the bubble is exerted close to the rotation axis, driving a collimated outflow out through the host star. This can account for the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). Given that the wind magnetization increases in time, we thus suggest that the magnetar-driven bubble initially expands relatively spherically (enhancing the energy of the associated supernova) while at late times it becomes progressivelymore collimated (producing the GRB). Similar processes may operate in more modestly rotating neutron stars to produce asymmetric supernovae and lower energy transients such as X-ray flashes.

Gamma ray bursts have been divided into two classes, long-soft gamma ray burst and short-hard gamma ray burst according to the bimodal distribution in duration time. Due to the harder spectrum and the lack of afterglows of short-hard bursts in optical and radio observations, different progenitors for short-hard bursts and long-soft bursts have been suggested. Based on the X-ray afterglow observation and the cumulative red-shift distribution of short-hard bursts, Nakar et al. (2006) found that the progenitors of short-hard bursts are consistent with old populations, such as mergers of binary neutron stars. Recently, the existence of two subclasses in long-soft bursts has been suggested after considering multiple characteristics of gamma-ray bursts, including fluences and the duration time. In this work, we extended the analysis of cumulative red-shift distribution to two possible subclasses in L-GRBs. We found that two possible subclass GRBs show different red-shift distributions, especially for red-shifts z > 1. Our results indicate that the accumulative red-shift distribution can be used as a tool to constrain the progenitor characteristics of possible subclasses in L-GRBs.

We present the general relativistic calculation of the energy release associated with a first order phase transition (PT) at the center of a rotating neutron star (NS). The energy release, E_rel, is equal to the difference in mass-energies between the initial (normal) phase configuration and the final configuration containing a superdense matter core, assuming constant total baryon number and the angular momentum. The calculations are performed with the use of precise pseudo-spectral 2-D numerical code; the polytropic equations of state (EOS) as well as realistic EOSs (Skyrme interactions, Mean Field Theory kaon condensate) are used. The results are obtained for a broad range of metastability of initial configuration and size of the new superdense phase core in the final configuration. For a fixed ``overpressure'', dP, defined as the relative excess of central pressure of a collapsing metastable star over the pressure of the equilibrium first-order PT, the energy release up to numerical accuracy does not depend on the stellar angular momentum and coincides with that for nonrotating stars with the same dP. When the equatorial radius of the superdense phase core is much smaller than the equatorial radius of the star, analytical expressions for the E_rel can be obtained: E_rel is proportional to dP^2.5 for small dP. At higher dP, the results of 1-D calculations of E_rel(dP) for non-rotating stars reproduce with very high precision exact 2-D results for fast-rotating stars. The energy release-angular momentum independence for a given overpressure holds also for the so-called ``strong'' PTs (that destabilise the star against the axi-symmetric perturbations), as well as for PTs with ``jumping'' over the energy barrier.

We present our new advanced model for population synthesis of close-by cooling NSs. Detailed treatment of the initial spatial distribution of NS progenitors and a detailed ISM structure up to 3 kpc give us an opportunity to discuss the strategy to look for new isolated cooling NSs. Our main results in this respect are the following: new candidates are expected to be identified behind the Gould Belt, in directions to rich OB associations, in particular in the Cygnus-Cepheus region; new candidates, on average, are expected to be hotter than the known population of cooling NS. Besides the usual approach (looking for soft X-ray sources), the search in 'empty' $\gamma$-ray error boxes or among run-away OB stars may yield new X-ray thermally emitting NS candidates.

The accreting millisecond pulsar SAX J1808.4-3658 may be a transition object between accreting X-ray binaries and millisecond radio pulsars. We have constrained the thermal radiation from its surface through XMM-Newton X-ray observations, providing strong evidence for neutrino cooling processes from the neutron star core. We have also undertaken simultaneous X-ray and optical (Gemini) observations, shedding light on whether the strong heating of the companion star in quiescence may be due to X-ray irradiation, or to a radio pulsar turning on when accretion stops.

Prior to the incineration of a white dwarf (WD) that makes a Type Ia supernova (SN Ia), the star "simmers" for ~1000 years in a convecting, carbon burning region. We have found that weak interactions during this time increase the neutron excess by an amount that depends on the total quantity of carbon burned prior to the explosion. This contribution is in addition to the metallicity (Z) dependent neutronization through the 22Ne abundance (as studied by Timmes, Brown, & Truran). The main consequence is that we expect a floor to the level of neutronization that dominates over the metallicity contribution when Z/Z_\odot<2/3, and it can be important for even larger metallicities if substantial energy is lost to neutrinos via the convective Urca process. This would mask any correlations between SN Ia properties and galactic environments at low metallicities. In addition, we show that recent observations of the dependences of SNe Ia on galactic environments make it clear that metallicity alone cannot provide for the full observed diversity of events.

A dissipative Lorentz-covariant Ohm's law which uses only the electromagnetic degrees of freedom is proposed. For large conductivity, Maxwell equations equipped with this Ohm's law reduce to the equations of Force-Free Electrodynamics (FFE) with small dissipative corrections, but only in the regions where the ideal FFE 4-current is space-like. This might indicate that the pulsar emission comes primarily from the magnetic separartrix.

Our high time resolution observations of individual pulses from the Crab pulsar show that the main pulse and interpulse differ in temporal behavior, spectral behavior, polarization and dispersion. The main pulse properties are consistent with one current model of pulsar radio emission, namely, soliton collapse in strong plasma turbulence. The high-frequency interpulse is quite another story. Its dynamic spectrum cannot easily be explained by any current emission model; its excess dispersion must come from propagation through the star's magnetosphere. We suspect the high-frequency interpulse does not follow the ``standard model'', but rather comes from some unexpected region within the star's magnetosphere. Similar observations of other pulsars will reveal whether the radio emission mechanisms operating in the Crab pulsar are unique to that star, or can be identified in the general population.

We present an analysis of five X-ray Multi-Mirror Mission (XMM) observations of the Anomalous X-ray Pulsar (AXP) 1E 2259+586 taken in 2004 and 2005 during its relaxation following its 2002 outburst. We compare these data with those of five previous XMM observations taken in 2002 and 2003, and find the observed flux decay is well described by a power-law of index -0.69+/-0.03. As of mid-2005, the source may still have been brighter than pre-outburst, and was certainly hotter. We find a strong correlation between hardness and flux, as seen in other AXP outbursts. We discuss the implications of these results for the magnetar model.

Inspiral signals from binary compact objects (black holes and neutron stars) are primary targets of the ongoing searches by a number of ground-based gravitational-wave interferometers (LIGO, Virgo, GEO-600 and TAMA-300). Detection of such inspirals and ensuing mergers is expected to provide us with important physical information about the properties of the sources, bearing on outstanding issues in compact-object astrophysics, including the progenitors of short gamma-ray bursts. Compact-object spin effects add to the challenges associated with searches and anticipated detections, but on the other hand they provide some interesting possibilities for extracting astrophysical information. We present parameter-estimation simulations for inspirals of black-hole binaries with neutron-star companions using Markov-Chain Monte-Carlo methods. We specifically highlight the potential for measurements of masses, spins, source sky location and distance of such objects with just one or two gravitational-wave detectors.

In this paper, we present a detailed hydrodynamical study of the properties of the flow produced by the collision of a pulsar wind with the surrounding in a binary system. This work is the first attempt to simulate interaction of the ultrarelativistic flow (pulsar wind) with the nonrelativistic stellar wind. Obtained results show that the wind collision could result in the formation of an "unclosed" (at spatial scales comparable to the binary system size) pulsar wind termination shock even when the stellar wind ram pressure exceeds significantly the pulsar wind kinetical pressure. Moreover, the post-shock flow propagates in a rather narrow region, with very high bulk Lorentz factor ($\gamma\sim100$). This flow acceleration is related to adiabatical losses, which are purely hydrodynamical effects. Interestingly, in this particular case, no magnetic field is required for formation of the ultrarelativistic bulk outflow. The obtained results provide a new interpretation for the orbital variability of radio, X-ray and gamma-ray signals detected from binary pulsar system PSR 1259-63/SS2883.

We propose that the strong millisecond extragalactic radio burst (mERB) discovered by Lorimer et al. (2007) may be related to a hyperflare from an extragalactic soft gamma-ray repeater. The expected rate of such hyperflares, $\sim$ 20 - 50 d$^{-1}$ Gpc$^{-3}$, is in good correspondence with the value estimated by Lorimer et al. The possible mechanism of radio emission can be related to the tearing mode instability in the magnetar magnetosphere as discussed by Lyutikov (2002), and can produce the radio flux corresponding to the observed $\sim$ 30 Jy from the mERB using a simple scaling of the burst energy.

We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and superfluid properties of neutron star cores by comparing the cooling theory with observations of thermal radiation from isolated neutron stars.

Rotating Radio Transients (RRATs) are a new class of neutron stars discovered through the emission of radio bursts. Eleven sources are known up to now, but population studies predict these objects to be more numerous than the normal radio pulsar population. Multiwavelength observations of these peculiar objects are in progress to disentangle their spectral energy distribution, and then study in detail their nature. In this review I report on the current state of the art on these objects, and in particular on the results of new X-ray observations.

We present new determination of the birth rate of AXPs and SGRS and their associated SNRs. We find a high birth rate of 1/(500 yr) for AXPs/SGRs and 1/(1700 yr) for associated SNRs. These high rates suggest that all massive stars (greater than ~ 25 M_sun) give rise to remnants with magnetar-like fields. Recent observations indicate that fossil fields cannot explain such high fields in the progenitor stars. Dynamo mechanisms during the birth of the neutron stars require spin rates much faster than either observations or theory indicate. Here, we propose the neutron stars form with normal (~ 10^{12} G) magnetic fields, which are then amplified to 10^{14}-10^{15} G after a delay of a few hundred years. The amplification is a consequence of color ferromagnetism and occurs after the neutron star core reaches quark-deconfinement density. This delayed amplification alleviates many difficulties in interpreting simultaneously the high birth rate and high magnetic fields of AXPs/SGRs and their link to massive stars.

We have searched for pulsed radio emission from magnetar 4U 0142+61 at the frequency of 111 MHz. No pulsed signal was detected from this source. Upper limits for mean flux density are 0.9 - 9 mJy depending on assumed duty cycle (.05 - .5) of the pulsar.

We have studied non-axisymmetric standing accretion shock instability, or SASI, by 3D hydrodynamical simulations. This is an extention of our previous study on axisymmetric SASI. We have prepared a spherically symmetric and steady accretion flow through a standing shock wave onto a proto-neutron star, taking into account a realistic equation of state and neutrino heating and cooling. This unperturbed model is supposed to represent approximately the typical post-bounce phase of core-collapse supernovae. We then have added a small perturbation (~1%) to the radial velocity and computed the ensuing evolutions. Not only axisymmetric but non-axisymmetric perturbations have been also imposed. We have applied mode analysis to the non-spherical deformation of the shock surface, using the spherical harmonics. We have found that (1) the growth rates of SASI are degenerate with respect to the azimuthal index m of the spherical harmonics Y_l^m, just as expected for a spherically symmetric background, (2) nonlinear mode couplings produce only m=0 modes for the axisymmetric perturbations, whereas m=!0 modes are also generated in the non-axisymmetric cases according to the selection rule for the quadratic couplings, (3) the nonlinear saturation level of each mode is lower in general for 3D than for 2D because a larger number of modes are contributing to turbulence in 3D, (4) low l modes are dominant in the nonlinear phase, (5) the equi-partition is nearly established among different m modes in the nonlinear phase, (6) the spectra with respect to l obey power laws with a slope slightly steeper for 3D, and (7) although these features are common to the models with and without a shock revival at the end of simulation, the dominance of low l modes is more remarkable in the models with a shock revival.

We have detected 523 sources in a survey of the Small Magellanic Cloud (SMC) Wing with Chandra. By cross-correlating the X-ray data with optical and near-infrared catalogues we have found 300 matches. Using a technique that combines X-ray colours and X-ray to optical flux ratios we have been able to assign preliminary classifications to 265 of the objects. Our identifications include four pulsars, one high-mass X-ray binary (HMXB) candidate, 34 stars and 185 active galactic nuclei (AGNs). In addition, we have classified 32 sources as 'hard' AGNs which are likely absorbed by local gas and dust, and nine 'soft' AGNs whose nature is still unclear. Considering the abundance of HMXBs discovered so far in the Bar of the SMC the number that we have detected in the Wing is low.

The deeper and more extended survey of the central parts of the Galactic Plane by H.E.S.S. during 2005-2007 has revealed a number of new point-like, as well as, extended sources. Two point-like sources can be associated to two remarkable objects around "Crab-like" young and energetic pulsars in our Galaxy: G21.5-0.9 and Kes 75. The characteristics of each of the sources are presented and possible interpretations are briefly discussed.

We present measurements of the neutron-capture elements Rb and Pb for bright giants in the globular clusters M4 and M5. The clusters are of similar metallicity ([Fe/H] = -1.2) but M4 is decidedly s-process enriched relative to M5: [Ba/Fe] = +0.6 for M4 but 0.0 for M5. The Rb and Pb abundances were derived by comparing synthetic spectra with high-resolution, high signal-to-noise ratio spectra obtained with MIKE on the Magellan telescope. Abundances of Y, Zr, La, and Eu were also obtained. In M4, the mean abundances from 12 giants are [Rb/Fe] = 0.39 +/- 0.02 (sigma = 0.07), [Rb/Zr] = 0.17 +/- 0.03 (sigma = 0.08), and [Pb/Fe] = 0.30 +/- 0.02 (sigma = 0.07). In M5, the mean abundances from two giants are [Rb/Fe] = 0.00 +/- 0.05 (sigma = 0.06), [Rb/Zr] = 0.08 +/- 0.08 (sigma = 0.11), and [Pb/Fe] = -0.35 +/- 0.02 (sigma = 0.04). Within the measurement uncertainties, the abundance ratios [Rb/Fe], [Pb/Fe] and [Rb/X] for X = Y, Zr, La are constant from star-to-star in each cluster and none of these ratios are correlated with O or Na abundances. While M4 has a higher Rb abundance than M5, the ratios [Rb/X] are similar in both clusters indicating that the nature of the s-products are very similar for each cluster but the gas from which M4's stars formed had a higher concentration of these products.

Experimental evidence is reviewed for the existence of superfluid turbulence in a differentially rotating, spherical shell at high Reynolds numbers ($\Rey\gsim 10^3$), such as the outer core of a neutron star. It is shown that torque variability increases with $\Rey$, suggesting that glitch activity in radio pulsars may be a function of $\Rey$ as well. The $\Rey$ distribution of the 67 glitching radio pulsars with characteristic ages $\tau_c \leq 10^6$ {\rm yr} is constructed from radio timing data and cooling curves and compared with the $\Rey$ distribution of all 348 known pulsars with $\tau_c \leq 10^6$ {\rm yr}. The two distributions are different, with a Kolmogorov-Smirnov probability $\geq 1 - 3.9 \times 10^{-3}$. The conclusion holds for (modified) Urca and nonstandard cooling, and for Newtonian and superfluid viscosities.

The giant flares produced by highly magnetized neutron stars, "magnetars," are the brightest sources of high energy radiation outside our solar system. Serendipitous observations with NASA's Rossi X-ray Timing Explorer (RXTE) of the two most recent flares resulted in the discovery of high frequency oscillations in their X-ray fluxes. The frequencies of these oscillations range from about 20 Hz to as high as 1800 Hz, and may represent the first detection of global oscillation modes of neutron stars. Here I will present an observational and theoretical overview of these oscillations and discuss how they might allow us to probe neutron star interiors and dense matter physics.

INTEGRAL, the European Space Agency's gamma-ray observatory, tripled the number of super-giant high-mass X-ray binaries (sgHMXB) known in the Galaxy by revealing absorbed and fast transient (SFXT) systems. In these sources, quantitative constraints on the wind clumping of the massive stars could be obtained from the study of the hard X-ray variability of the compact accreting object. Hard X-ray flares and quiescent emission of SFXT systems have been characterized and used to derive wind clump parameters. A large fraction of the hard X-ray emission is emitted in the form of flares with a typical duration of 3 ks, frequency of 7 days and luminosity of 1E36 erg/s. Such flares are most probably emitted by the interaction of a compact object orbiting at ~10 R* with wind clumps (1E(22-23) g) representing a large fraction of the stellar mass-loss rate. The density ratio between the clumps and the inter-clump medium is 1E(2-4) in SFXT systems. The parameters of the clumps and of the inter-clump medium, derived from the SFXT flaring behavior, are in good agreement with macro-clumping scenario and line driven instability simulations. SFXT have probably a larger orbital radius than classical sgHMXB.

We observed four southern AXPs in 1999 near 1400 MHz with the Parkes 64-m radio telescope to search for periodic radio emission. No Fourier candidates were discovered in the initial analysis, but the recent radio activity observed for the AXP XTE J1810-197 has prompted us to revisit these data to search for single radio pulses and bursts. The data were searched for both persistent and bursting radio emission at a wide range of dispersion measures, but no detections of either kind were made. These results further weaken the proposed link between rotating radio transient sources and magnetars. However, continued radio searches of these and other AXPs at different epochs are warranted given the transient nature of the radio emission seen from XTE J1810-197, which until very recently was the only known radio-emitting AXP.

We report the detection of 10 new X-ray filaments using the data from the {\sl Chandra} X-ray satellite for the inner $6^{\prime}$ ($\sim 15$ parsec) around the Galactic center (GC). All these X-ray filaments are characterized by non-thermal energy spectra, and most of them have point-like features at their heads that point inward. Fitted with the simple absorbed power-law model, the measured X-ray flux from an individual filament in the 2-10 keV band is $\sim 2.8\times10^{-14}$ to $10^{-13}$ ergs cm$^{-2}$ s$^{-1}$ and the absorption-corrected X-ray luminosity is $\sim 10^{32}-10^{33}$ ergs s$^{-1}$ at a presumed distance of 8 kpc to the GC. We speculate the origin(s) of these filaments by morphologies and by comparing their X-ray images with the corresponding radio and infrared images. On the basis of combined information available, we suspect that these X-ray filaments might be pulsar wind nebulae (PWNe) associated with pulsars of age $10^3 \sim 3\times 10^5$ yr. The fact that most of the filament tails point outward may further suggest a high velocity wind blowing away form the GC.

In December 2004, the soft gamma-ray repeater SGR 1806-20 emitted the most powerful giant flare ever observed. This probably involved a large-scale rearrangement of the magnetosphere leading to observable variations in the properties of its X-ray emission. Here we present the results of the first Suzaku observation of SGR 1806-20, together with almost simultaneous observations with XMM-Newton and INTEGRAL. The source seems to have reached a state characterized by a flux close to the pre-flare level and by a relatively soft spectrum. Despite this, SGR 1806-20 remained quite active also after the giant flare, allowing us to study several short bursts observed by Suzaku in the 1-100 keV range. We discuss the broad-band spectral properties of SGR 1806-20, covering both persistent and bursting emission, in the context of the magnetar model, and consider its recent theoretical developments.

SGRs denote ``soft $\gamma$-ray repeaters'', a small class of slowly spinning neutron stars with strong magnetic fields. On 27 December 2004, a giant flare was detected from magnetar SGR 1806-20. The initial spike was followed by a hard-X-ray tail persisting for 380 s with a modulation period of 7.56 s. This event has received considerable attention, particularly in the astrophysics area. Its relevance to the geophysics community lies in the importance of investigating the effects of such an event on the near-earth electromagnetic environment. However, the signature of a magnetar flare on the geomagnetic field has not previously been investigated. Here, by applying wavelet analysis to the high-resolution magnetic data provided by the CHAMP satellite, a modulated signal with a period of 7.5 s over the duration of the giant flare appears in the observed data. Moreover, this event was detected by the energetic ion counters onboard the DEMETER satellite.

Neutron stars contain persistent, ordered magnetic fields that are the strongest known in the Universe. However, their magnetic fluxes are similar to those in magnetic A and B stars and white dwarfs, suggesting that flux conservation during gravitational collapse may play an important role in establishing the field, although it might also be modified substantially by early convection, differential rotation, and magnetic instabilities. The equilibrium field configuration, established within hours (at most) of the formation of the star, is likely to be roughly axisymmetric, involving both poloidal and toroidal components. The stable stratification of the neutron star matter (due to its radial composition gradient) probably plays a crucial role in holding this magnetic structure inside the star. The field can evolve on long time scales by processes that overcome the stable stratification, such as weak interactions changing the relative abundances and ambipolar diffusion of charged particles with respect to neutrons. These processes become more effective for stronger magnetic fields, thus naturally explaining the magnetic energy dissipation expected in magnetars, at the same time as the longer-lived, weaker fields in classical and millisecond pulsars.

We report on our initial analysis of a deep 510 ks observation of the Galactic oxygen-rich supernova remnant (SNR) G292.0+1.8 with the {\it Chandra X-ray Observatory}. Our new {\it Chandra} ACIS-I observation has a larger field of view and an order of magnitude deeper exposure than the previous {\it Chandra} observation, which allows us to cover the entire SNR and to detect new metal-rich ejecta features. We find a highly non-uniform distribution of thermodynamic conditions of the X-ray emitting hot gas that correlates well with the optical [O {\small III}] emission, suggesting the possibility that the originating supernova explosion of G292.0+1.8 was itself asymmetric. We also reveal spectacular substructures of a torus, a jet, and an extended central compact nebula all associated with the embedded pulsar J1124$-$5916.

We report on the discovery of three new pulsars in the first blind survey of the north Galactic plane (45 < l < 135 ; |b| < 1) with the Giant Meterwave Radio telescope (GMRT) at an intermediate frequency of 610 MHz. The timing parameters, obtained in follow up observations with the Lovell Telescope at Jodrell Bank Observatory and the GMRT, are presented.

A perturbative analysis is used to investigate the effect of rotation on the instability of a steady accretion shock (SASI) in a simple toy-model, in view of better understanding supernova explosions in which the collapsing core contains angular momentum. A cylindrical geometry is chosen for the sake of simplicity. Even when the centrifugal force is very small, rotation can have a strong effect on the non-axisymmetric modes of SASI by increasing the growth rate of the spiral modes rotating in the same direction as the steady flow. Counter-rotating spiral modes are significantly damped, while axisymmetric modes are hardly affected by rotation. The growth rates of spiral modes have a nearly linear dependence on the specific angular momentum of the flow. The fundamental one-armed spiral mode (m=1) is favoured for small rotation rates, whereas stronger rotation rates favour the mode m=2. A WKB analysis of higher harmonics indicates that the efficiency of the advective-acoustic cycles associated to spiral modes is strongly affected by rotation in the same manner as low frequency modes, whereas the purely acoustic cycles are stable. These results suggest that the linear phase of SASI in rotating core-collapse supernovae naturally selects a spiral mode rotating in the same direction of the flow, as observed in the 3D numerical simulations of Blondin & Mezzacappa (2007). This emphasizes the need for a 3D approach of core-collapse, before conclusions on the explosion mechanisms and pulsar kicks can be drawn.

The Parity doublet model containing the SU(2) multiplets including the baryons identified as the chiral partners of the nucleons is applied for neutron star matter. The chiral restoration is analyzed and the maximum mass of the star is calculated.

Short-hard Gamma Ray Bursts (SGRBs) are currently thought to arise from gravitational wave driven coalescences of double neutron star systems forming either in the field or dynamically in globular clusters. For both channels we fit the peak flux distribution of BATSE SGRBs to derive the local burst formation rate and luminosity function. We then compare the resulting redshift distribution with the Swift data, showing that both formation channels are needed in order to reproduce the observations. Double neutron stars forming in globular clusters are found to dominate the distribution at z<0.3, whereas the field population from primordial binaries can account for the high-z SGRBs. This result is not in contradiction with the observed host galaxy type of SGRBs.

We study CP-violation effects when neutrinos are present in dense matter, such as outside the proto-neutron star formed in a core-collapse supernova. General arguments based on the Standard Model predict that there are no CP-violating effects at the tree level in core-collapse supernovae, justifying the assumption of ignoring the Dirac delta phase made in most of the literature. Physics beyond the Standard Model, especially flavor-changing interactions might induce CP-violating effects. We analyze this possibility.

We present a study of the globular cluster (GC) systems of nearby elliptical and S0 galaxies at a variety of wavelengths from the X-ray to the infrared. Our analysis shows that roughly half of the low mass X-ray binaries (LMXBs), that are the luminous tracers of accreting neutron star or black hole systems, are in clusters. There is a surprisingly strong correlation between the LMXB frequency and the metallicity of the GCs, with metal-rich GCs hosting three times as many LMXBs as metal-poor ones, and no convincing evidence of a correlation with GC age so far. In some of the galaxies the LMXB formation rate varies with GC color even within the red peak of the typical bimodal cluster color distribution, providing some of the strongest evidence to date that there are metallicity variations within the metal-rich GC peak as is expected in a hierarchical galaxy formation scenario. We also note that any analysis of subtler variations in GC color distributions must carefully account for both statistical and systematic errors. We caution that some published GC correlations, such as the apparent 'blue-tilt' or mass-metallicity effect might not have a physical origin and may be caused by systematic observational biases.

During accretion, a neutron star (NS) is spun up as angular momentum is transported through its liquid surface layers. We study the resulting differentially rotating profile, focusing on the impact this has for type I X-ray bursts. The viscous heating is found to be negligible, but turbulent mixing can be activated. Mixing has the greatest impact when the buoyancy at the compositional discontinuity between accreted matter and ashes is overcome. This occurs preferentially at high accretion rates or low spin frequencies and may depend on the ash composition from the previous burst. We then find two new regimes of burning. The first is ignition in a layer containing a mixture of heavier elements with recurrence times as short as ~5-30 minutes, similar to short recurrence time bursts. When mixing is sufficiently strong, a second regime is found where accreted helium mixes deep enough to burn stably, quenching X-ray bursts altogether. The carbon-rich material produced by stable helium burning would be important for triggering and fueling superbursts.

Very high energy (VHE) gamma-ray observations have proven to be very successful in localizing Galactic acceleration sites of VHE particles. Observations of shell-type supernova remnants have confirmed that particles are accelerated to VHE energies in supernova blast waves; the interpretation of the gamma-ray data in terms of hadronic or leptonic particle components in these objects relies nevertheless strongly on input from X-ray observations. The largest identified Galactic VHE source class consists of pulsar wind nebulae, as detected in X-rays. Many of the remaining VHE sources remain however unidentified until now. With X-ray observations of these enigmatic "dark" objects one hopes to solve the following questions: What is the astrophysical nature of these sources? Are they predominantly electron or hadron accelerators? And what is their contribution to the overall cosmic ray energy budget? The paper aims to provide an overview over the identification status of the Galactic VHE source population.

The focusing of the radiation generated by a polarization current with a superluminally rotating distribution pattern is of a higher order in the plane of rotation than in other directions. Consequently, our previously published asymptotic approximation to the value of this field outside the equatorial plane breaks down as the line of sight approaches a direction normal to the rotation axis, i.e., is nonuniform with respect to the polar angle. Here we employ an alternative asymptotic expansion to show that, though having a rate of decay with frequency (mu) that is by a factor of order mu^(2/3) slower, the equatorial radiation field has the same dependence on distance as the nonspherically decaying component of the generated field in other directions: it, too, diminishes as the inverse square root of the distance from its source. We also briefly discuss the relevance of these results to the giant pulses received from pulsars: the focused, nonspherically decaying pulses that arise from a superluminal polarization current in a highly magnetized plasma have a power-law spectrum (i.e., a flux density proportional to mu^alpha) whose index (alpha) is given by one of the values -2/3, -2, -8/3, or -4.

Using a numerical simulation, we study the effects of ambipolar diffusion and ohmic diffusion on the magnetic field evolution in the interior of an isolated neutron star. We are interested in the behavior of the magnetic field on a long time scale, over which all Alfven and sound waves have been damped. We model the stellar interior as an electrically neutral plasma composed of neutrons, protons and electrons, which can interact with each other through collisions and electromagnetic forces. Weak interactions convert neutrons and charged particles into each other, erasing chemical imbalances. As a first step, we assume that the magnetic field points in one fixed Cartesian direction but can vary along an orthogonal direction. We start with a uniform-density background threaded by a homogeneous magnetic field and study the evolution of a magnetic perturbation as well as the density fluctuations it induces in the particles. We show that the system evolves through different quasi-equilibrium states and estimate the characteristic time scales on which these quasi-equilibria occur.

We performed a timing analysis of the 2003 outburst of the accreting X-ray millisecond pulsar XTE J1807-294 observed by RXTE. Using recently refined orbital parameters we report for the first time a precise estimate of the spin frequency and of the spin frequency derivative. The phase delays of the pulse profile show a strong erratic behavior superposed to what appears as a global spin-up trend. The erratic behavior of the pulse phases is strongly related to rapid variations of the light curve, making it very difficult to fit these phase delays with a simple law. As in previous cases, we have therefore analyzed separately the phase delays of the first harmonic and of the second harmonic of the spin frequency, finding that the phases of the second harmonic are far less affected by the erratic behavior. In the hypothesis that the second harmonic pulse phase delays are a good tracer of the spin frequency evolution we give for the first time a estimation of the spin frequency derivative in this source. The source shows a clear spin-up of $\dot \nu = 2.5^{+1.0}_{-1.1} \times 10^{-14}$ Hz sec$^{-1}$ (90% confidence level). The largest source of uncertainty in the value of the spin-up rate is given by the uncertainties on the source position in the sky. We discuss this systematics on the spin frequency and its derivative.

Forty years after the discovery of rotation-powered pulsars, we still do not understand many aspects of their pulsed emission. In the last few years there have been some fundamental developments in acceleration and emission models. I will review both the basic physics of the models as well as the latest developments in understanding the high-energy emission of rotation-powered pulsars, with particular emphasis on the polar-cap and slot-gap models. Special and general relativistic effects play important roles in pulsar emission, from inertial frame-dragging near the stellar surface to aberration, time-of-flight and retardation of the magnetic field near the light cylinder. Understanding how these effects determine what we observe at different wavelengths is critical to unraveling the emission physics. I will discuss how current and future X-ray and gamma-ray detectors can test the predictions of these models.

Globular clusters produce orders of magnitude more millisecond pulsars per unit mass than the Galactic disk. Since the first cluster pulsar was uncovered twenty years ago, at least 138 have been identified - most of which are binary millisecond pulsars. Because of their origins involving stellar encounters, many of these systems are exotic objects that would never be observed in the Galactic disk. Examples include pulsar-main sequence binaries, extremely rapid rotators (including the current record holder), and millisecond pulsars in highly eccentric orbits. These systems are allowing new probes of the interstellar medium, the equation of state of material at supra-nuclear density, the mass distribution of neutron stars, and the dynamics of globular clusters.

The Newtonian solid-mechanical theory of nodeless spheroidal and torsional seismic elastic vibrations trapped in the crust of quaking neutron star is outlined. The spectral equations for the frequency of these modes are obtained and applied to the modal classification of the quasi-periodic oscillations of X-ray luminosity in the aftermath of giant flares in SGR 1806-20 and SGR 1900+14. The presented analysis is heavily relied on the currently accepted identification of the QPOs frequency from the range [30-200] Hz with those for torsional nodeless vibrations. Based on this identification, which is used to fix the input parameters entering the obtained spectral formulae, we compute frequency spectrum of nodeless spheroidal elastic vibrations. Focus is placed on the low-frequency QPOs in the data for SGR 1806-20 whose physical origin has been called into question. Our calculations suggest that QPOs with frequencies 18 Hz and 26 Hz are due to dipole spheroidal and dipole torsional shear vibrations of the crust against unperturbed by starquake core, respectively. The uncertainties in association of the high-frequency QPOs with nodeless elastic shear vibrations are briefly discussed.

We report results of spectral studies of steady emission and short bursts for magnetar candidates using data obtained by the XMM-Newton, Chandra and Swift spacecrafts. The spectra of the steady emission and the short bursts for most magnetar candidates are represented by a two blackbody function with an absorption (2BB). Three AXPs (4U 0142+614, 1RXS J170849.0-400910 and 1E 2259+586) seem to have an excess above ~7 keV in their soft X-ray spectra. The excess could be related to a non-thermal hard component discovered by the International Gamma-Ray Astrophysics Laboratory. For the magnetar candidates of which the spectra are well reproduced by 2BB, there is a strong linear correlation between lower blackbody temperatures (kT_LT) and higher blackbody temperatures (kT_HT), and between squares of lower blackbody radii (R_LT^2) and squares of higher blackbody radii (R_HT^2). The ratio of kT_LT/kT_HT ~ 0.4 is almost constant independently with magnetar candidates and/or emission types (the steady emission and the short bursts). This would imply a common emission mechanism among the magnetar candidates. The relationship between the steady emission and the short bursts might be analogous to a relationship between microflares and solar flares of the Sun.

We present here X-ray spectra of the HMXB SMC X-1 obtained in an observation with the XMM observatory beginning before eclipse and ending near the end of eclipse. With the Reflection Grating Spectrometers (RGS) on board XMM, we observe emission lines from hydrogen-like and helium-like ions of nitrogen, oxygen, neon, magnesium, and silicon. Though the resolution of the RGS is sufficient to resolve the helium-like n=2->1 emission into three line components, only one of these components, the intercombination line, is detected in our data. The lack of flux in the forbidden lines of the helium-like triplets is explained by pumping by ultraviolet photons from the B0 star and, from this, we set an upper limit on the distance of the emitting ions from the star. The lack of observable flux in the resonance lines of the helium-like triplets indicate a lack of enhancement due to resonance line scattering and, from this, we derive a new observational constraint on the distribution of the wind in SMC X-1 in velocity and coordinate space. We find that the solid angle subtended by the volume containing the helium-like ions at the neutron star multiplied by the velocity dispersion of the helium-like ions must be less than 4pi steradians km/s. This constraint will be satisfied if the helium-like ions are located primarily in clumps distributed throughout the wind or in a thin layer along the surface of the B0 star.

We argue that giant flares in SGRs can be associated to the core conversion of an isolated neutron star having a subcritical magnetic field $\sim 10^{12}$ G and a fallback disk around it. We show that, in a timescale of $\lesssim 10^5$ yrs, accretion from the fallback disk can increase the mass of the central object up to the critical mass for the conversion of the core of the star into quark matter. A small fraction of the neutrino-antineutrino emission from the just-converted quark-matter hot core annihilates into $e^+e^-$ pairs above the neutron star surface originating the gamma emission of the spike while the further cooling of the heated neutron star envelope originates the tail of the burst. We show that several characteristics of the giant flare of the SGR 1806-20 of 27 December 2004 (spike and tail energies, timescales, and spectra) can be explained by this mechanism.

The scattering of electromagnetic radiation by the particle gyrating in an external magnetic field is considered. Particular attention is paid to the low-frequency case, when the frequencies of incident radiation are much less than the electron gyrofrequency. The spectral and polarization features of the scattering cross-section are analyzed in detail. It is found that the scattering transfers the low-frequency photons to high harmonics of the gyrofrequency, into the range of the synchrotron emission of the electron. The total scattering cross-section appears much larger than that for the particle at rest. The problem studied is directly applicable to the radio wave scattering in the magnetosphere of a pulsar. The particles acquire relativistic rotational energies as a result of resonant absorption of the high-frequency radio waves and concurrently scatter the low-frequency radio waves, which are still below the resonance. It is shown that the scattering can affect the radio intensity and polarization at the lowest frequencies and can compete with the resonant absorption in contributing to the low-frequency turnover in the pulsar spectrum. Moreover, the scattering can be an efficient mechanism of the pulsar high-energy emission, in addition to the synchrotron re-emission of the particles. Other astrophysical applications of the scattering by gyrating particles are pointed out as well.

The pulse sequence is interpreted as a realization of a random electron discharge process in a vacuum gap over the polar cap (PC) - the open magnetic force line region on the neutron star surface. This point of view is illustrated by an example based on the dice. The generators of the random numbers are a cube and a coin. Throwing of dice and coin tossing determine the discharge places on the "light" and "dark" sides of PC, and correspondingly - the "shape" of the individual pulses and their statistical properties The physical mechanism giving such discharge scheme is shortly discussed. It may be a charge drained down from a sharp top of surface waves in a parallel electric field on the liquid PC surface of a neutron star.

The majority of X-ray burst sources do not display a burst rate that increases with luminosity as expected, but this is seen in the two clocked bursters XB1323-619 and GS1826-24. We present a detailed investigation of these two sources which in the case of the first source, spans 18 years. Based on measurements of the burst rate, X-ray luminosity, the alpha-parameter and the two time constants generally present in the burst decays, we demonstrate the importance of the rp nuclear burning process. A detailed comparison with theory shows that although the burst rate in each source agrees well with the theoretical value, there is a difference of more than a factor of 5 in the burst rate at a given luminosity between the sources. We show that the main reason for this is that the two sources have substantially different emitting areas on the neutron star in non-burst emission, a factor often neglected. Variation of this area may explain the inverse relation of burst rate with luminosity in the majority of burst sources.

Based on the results of applying the extended ADC emission model to three Z-track sources: GX340+0, GX5-1 and CygX-2, we propose an explanation of the Z-track sources in which the Normal and Horizontal Branches are dominated by the increasing radiation pressure of the neutron star. The emitted flux becomes several times super-Eddington at the Hard Apex and Horizontal Branch and we suggest that the inner accretion disk is disrupted by this and that part of the accretion flow is diverted vertically. This position on the Z-track is exactly the position where radio emission is detected showing the presence of jets. We thus propose that high radiation pressure is a necessary condition for the launching of jets. We also show that flaring must consist of unstable nuclear burning and that the mass accretion rate per unit emitting area of the neutron star mdot at the onset of flaring agrees well with the critical theoretical value at which burning becomes unstable.

We have applied the torus fitting procedure described in Ng & Romani (2004) to PWNe observations in the Chandra data archive. This study provides quantitative measurement of the PWN geometry and we characterize the uncertainties in the fits, with statistical errors coming from the fit uncertainties and systematic errors estimated by varying the assumed fitting model. The symmetry axis $\Psi$ of the PWN are generally well determined, and highly model-independent. We often derive a robust value for the spin inclination $\zeta$. We briefly discuss the utility of these results in comparison with new radio and high energy pulse measurements

We construct for the first time, the sequences of stable neutron star (NS) models capable of explaining simultaneously, the glitch healing parameters, $Q$, of both the pulsars, the Crab ($Q \geq 0.7$) and the Vela ($Q \leq 0.2$), on the basis of starquake mechanism of glitch generation. These models yield an upper bound on surface redshift of NSs, $z_R \simeq 0.77$. If the lower limit of the observational constraint of (i) $Q \geq 0.7$ for the Crab pulsar and (ii) the recently evaluated value of the moment of inertia for the Crab pulsar, $I_{\rm Crab,45} \geq 3.04$ (where $I_{45}=I/10^{45} {\rm g.cm}^2$), both are imposed together on these models, the models yield the value of matching density, $E_b = 7.0794 \times 10^{14}{\rm g cm}^{-3}$ at the core-envelope boundary. This value of matching density yields a model-independent upper bound on neutron star masses, $M_{\rm max} \leq 2.59 M_\odot$, and the strong lower bounds on surface redshift $z_R \simeq 0.6232$ and mass $M \simeq 2.455 M_\odot$ for the Crab ($Q \simeq 0.7$) and the strong upper bound on surface redshift $z_R \simeq 0.2016 $ and mass $M \simeq 1.142 M_\odot$ for the Vela ($Q \simeq 0.2$) pulsar. However, for the observational constraint of the `central' weighted mean value $Q \approx 0.72$, and $I_{\rm Crab,45} > 3.04$, for the Crab pulsar, the minimum surface redshift and mass of the Crab pulsar are slightly increased to the values $z_R \simeq 0.655$ and $M \simeq 2.5 M_\odot$ respectively, whereas for the central weighted mean value $Q \approx 0.12$ for the Vela pulsar, the maximum surface redshift and mass of the Vela pulsar are slightly decreased to the values $z_R \simeq 0.1645$ and $M \simeq 0.9635 M_\odot$ respectively.

We study the effect of the neutron star spin -- kick velocity alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios of the neutron star formation are considered: when the kick is always present and when it is small or absent if a neutron star is formed in a binary system due to electron-capture degenerate core collapse. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the expected galactic rate of binary neutron star coalescences compared to calculations with randomly directed kicks. The spin-kick correlation also leads to a much narrower NS spin-orbit misalignment.

We investigate the birth and evolution of isolated radio pulsars using a population synthesis method, modeling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam (MB) surveys. Together, the Parkes and Swinburne MB surveys have detected nearly 2/3 of the known pulsars and provide a remarkably homogeneous sample to compare with simulations. New proper motion measurements and an improved model of the distribution of free electrons in the interstellar medium, NE2001, also make revisiting these issues particularly worthwhile. We present a simple population model that reproduces the actual observations well, and consider others that fail. We conclude that: pulsars are born in the spiral arms, with the birthrate of 2.8+/-0.5 pulsars/century peaking at a distance ~3 kpc from the Galactic centre, and with mean initial speed of 380^{+40}_{-60} km/s; the birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality, implying that characteristic ages overestimate the true ages of the pulsars by a median factor >2 for true ages <30,000 yr; models in which the radio luminosities of the pulsars are random generically fail to reproduce the observed P-Pdot diagram, suggesting a relation between intrinsic radio luminosity and (P, Pdot); radio luminosities L Edot^0.5 provide a good match to the observed P-Pdot diagram; for this favored radio luminosity model, we find no evidence for significant magnetic field decay over the lifetime of the pulsars as radio sources ~100 Myr.

We present the analysis and results of 12.5 hours of high-energy gamma-ray observations of the EGRET-detected pulsar PSR B1951+32 using the Solar Tower Atmospheric Cherenkov Effect Experiment (STACEE). STACEE is an atmospheric Cherenkov detector, in Albuquerque, New Mexico, that detects cosmic gamma rays using the shower-front-sampling technique. STACEE's sensitivity to astrophysical sources at energies around 100 GeV allows it to investigate emission from gamma-ray pulsars with expected pulsed emission cutoffs below 100 GeV. We discuss the observations and analysis of STACEE's PSR 1951+32 data, accumulated during the 2005 and 2006 observing seasons.

We present 2D simulations of the cooling of neutron stars with strong magnetic fields (B \geq 10^{13} G). We solve the diffusion equation in axial symmetry including the state of the art microphysics that controls the cooling such as slow/fast neutrino processes, superfluidity, as well as possible heating mechanisms. We study how the cooling curves depend on the the magnetic field strength and geometry. Special attention is given to discuss the influence of magnetic field decay. We show that Joule heating effects are very large and in some cases control the thermal evolution. We characterize the temperature anisotropy induced by the magnetic field for the early and late stages of the evolution of isolated neutron stars.

We calculate axisymmetric toroidal modes of magnetized neutron stars with a solid crust. We assume the interior of the star is threaded by a poloidal magnetic field that is continuous at the surface with the outside dipole field whose strength $B_p$ at the magnetic pole is $B_p\sim 10^{16}$G. Since separation of variables is not possible for oscillations of magnetized stars, we employ finite series expansions of the perturbations using spherical harmonic functions to represent the angular dependence of the oscillation modes. For $B_p\sim 10^{16}$G, we find distinct mode sequences, in each of which the oscillation frequency of the toroidal mode slowly increases as the number of radial nodes of the eigenfunction increases. The frequency spectrum of the toroidal modes for $B_p\sim 10^{16}$G is largely different from that of the crustal toroidal modes of the non-magnetized model, although the frequency ranges are overlapped each other. This suggests that an interpretation of the observed QPOs based on the magnetic toroidal modes may be possible if the field strength of the star is as strong as $B_p\sim 10^{16}$G.

Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz. Such rotational frequency would be the first for which the effects of rotation are significant. We study the consequences of very fast rotating neutron stars for the potentially observable quantities as stellar mass and pulsar period.

Detection and study of gravitational waves from astrophysical sources is a major goal of current astrophysics. Ground-based laser-interferometer systems such as LIGO and VIRGO are sensitive to gravitational waves with frequencies of order 100 Hz, whereas space-based systems such as LISA are sensitive in the millihertz regime. Precise timing observations of a sample of millisecond pulsars widely distributed on the sky have the potential to detect gravitational waves at nanohertz frequencies. Potential sources of such waves include binary super-massive black holes in the cores of galaxies, relic radiation from the inflationary era and oscillations of cosmic strings. The Parkes Pulsar Timing Array (PPTA) is an implementation of such a system in which 20 millisecond pulsars have been observed using the Parkes radio telescope at three frequencies at intervals of two -- three weeks for more than two years. Analysis of these data has been used to limit the gravitational wave background in our Galaxy and to constrain some models for its generation. The data have also been used to investigate fluctuations in the interstellar and Solar-wind electron density and have the potential to investigate the stability of terrestrial time standards and the accuracy of solar-system ephemerides.

The group of 7 thermally emitting and radio-quiet isolated neutron stars (INSs) discovered by ROSAT constitutes a nearby population which locally appears to be as numerous as that of the classical radio pulsars. So far, attempts to enlarge this particular group of INSs finding more remote objects failed to confirm any candidate. We found in the 2XMMp catalogue a handful of sources with no catalogued counterparts and with X-ray spectra similar to those of the ROSAT discovered INSs, but seen at larger distances and thus undergoing higher interstellar absorptions. In order to rule out alternative identifications such as an AGN or a CV, we obtained deep ESO-VLT and SOAR optical imaging for the X-ray brightest candidates. We report here on the current status of our search and discuss the possible nature of our candidates. We focus particularly on the X-ray brightest source of our sample, 2XMM J104608.7-594306, observed serendipitously over more than four years by the XMM-Newton Observatory. A lower limit on the X-ray to optical flux ratio of ~ 300 together with a stable flux and soft X-ray spectrum make it the most promising thermally emitting INS candidate. Beyond the finding of new members, our study aims at constraining the space density of this population at large distances and at determining whether their apparently high local density is an anomaly or not.

I outline a mechanism, akin to Weibel instabilities of interpenetrating beams, in which the neighboring current sheets in a striped wind from an oblique rotator interact through a two stream-like mechanism (a Weibel instability in flatland), to create an anomalous resistivity that heats the sheets and causes the magnetic field to diffusively annihilate in the wind upstream of the termination shock. The heating has consequences for observable unpulsed emission from pulsars.

I share a few reminiscences and observations of 40 years of Pulsars.

While isolated neutron stars (INSs) are among the brightest gamma-ray sources, they are among the faintest ones in the optical, and their study is a challenging task which require the most powerful telescopes. HST has lead neutron star optical astronomy yielding nearly all the identifications achieved since the early 1990s. Here, the major HST contributions in the optical studies of INSs and their relevance for neutron stars' astronomy are reviewed.

We investigate the structure of protoneutron stars (PNS) formed by hadronic and quark matter in $\beta$-equilibrium described by appropriate equations of state (EOS). For the hadronic matter, we use a finite temperature EOS based on the Brueckner-Bethe-Goldstone many-body theory, with realistic two- and three-body forces. For the quark sector, we employ the Nambu--Jona-Lasinio model. We find that the maximum allowed masses are comprised in a narrow range around 1.8 solar masses, with a slight dependence on the temperature. Metastable hybrid protoneutron stars are not found.

We study X-ray spectra of Cyg X-3 from BeppoSAX, taking into account absorption and emission in the strong stellar wind of its companion. We find the intrinsic X-ray spectra are well modelled by disc blackbody emission, its upscattering by hot electrons with a hybrid distribution, and by Compton reflection. These spectra are strongly modified by absorption and reprocessing in the stellar wind, which we model using the photoionization code cloudy. The form of the observed spectra implies the wind is composed of two phases. A hot tenuous plasma containing most of the wind mass is required to account for the observed features of very strongly ionized Fe. Small dense cool clumps filling <0.01 of the volume are required to absorb the soft X-ray excess, which is emitted by the hot phase but not present in the data. The total mass-loss rate is found to be (0.6--1.6) 10^-5 solar masses per year. We also discuss the feasibility of the continuum model dominated by Compton reflection, which we find to best describe our data. The intrinsic luminosities of our models are compatible with the compact object being either a black hole or a neutron star.

During the analysis of the INTEGRAL observatory archival data we found a powerful X-ray burst, registered by JEM-X and IBIS/ISGRI telescopes on April 16, 2005 from a weak and poorly known source AX J1754.2-2754. Analysis of the burst profiles and spectrum shows, that it was a type I burst, which result from thermonuclear explosion on the surface of nutron star. It means that we can consider AX J1754.2-2754 as an X-ray burster. Certain features of burst profile at its initial stage witness of a radiation presure driven strong expansion and a corresponding cooling of the nutron star photosphere. Assuming, that the luminosity of the source at this phase was close to the Eddington limit, we estimated the distance to the burst source d=6.6+/-0.3 kpc (for hidrogen atmosphere of the neutron star) and d=9.2+/-0.4 kpc (for helium atmosphere).

Fusion reactions in the crust of an accreting neutron star are an important source of heat, and the depth at which these reactions occur is important for determining the temperature profile of the star. Fusion reactions depend strongly on the nuclear charge $Z$. Nuclei with $Z\le 6$ can fuse at low densities in a liquid ocean. However, nuclei with $Z=8$ or 10 may not burn until higher densities where the crust is solid and electron capture has made the nuclei neutron rich. We calculate the $S$ factor for fusion reactions of neutron rich nuclei including $^{24}$O + $^{24}$O and $^{28}$Ne + $^{28}$Ne. We use a simple barrier penetration model. The $S$ factor could be further enhanced by dynamical effects involving the neutron rich skin. This possible enhancement in $S$ should be studied in the laboratory with neutron rich radioactive beams. We model the structure of the crust with molecular dynamics simulations. We find that the crust of accreting neutron stars may contain micro-crystals or regions of phase separation. Nevertheless, the screening factors that we determine for the enhancement of the rate of thermonuclear reactions are insensitive to these features. Finally, we calculate the rate of thermonuclear $^{24}$O + $^{24}$O fusion and find that $^{24}$O should burn at densities near $10^{11}$ g/cm$^3$. The energy released from this and similar reactions may be important for the temperature profile of the star.

We communicate the detection of soft (20--200 keV) gamma-rays from the pulsar and pulsar wind nebula of PSR J1846-0258 and aim to identify the component of the system which is responsible for the gamma-ray emission. To pinpoint the source of gamma-ray emission we combine spectral information from the INTEGRAL gamma-ray mission with archival data from the Chandra X-ray Observatory. Our analysis shows that the soft gamma-rays detected from PSR J1846-0258 include emission from both the pulsar and the pulsar wind nebula, but the measured spectral shape is dominated by the pulsar wind nebula. We further discuss PSR J1846-0258 in the context of rotation-powered pulsars with high magnetic field strengths and review the anomalously high fraction of spin-down luminosity converted into X- and gamma-ray emission in light of a possible overestimate of the distance to this pulsar.

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars yield iron cores at collapse endowed with high angular momentum. It is thought that high angular momentum and black hole formation are critical ingredients of the collapsar model of long-soft GRBs. Here, we present 2D multi-group, flux-limited-diffusion MHD simulations of the collapse, bounce, and immediate post-bounce phases of a 35Msun collapsar-candidate model of Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the differentially-rotating surface layers of the millisecond-period proto-neutron star (PNS), a magnetically-driven explosion ensues during the PNS phase, in the form of a baryon-loaded non-relativistic jet, and that a black hole, central to the collapsar model, does not form. Paradoxically, and although much uncertainty surrounds mass loss, angular momentum transport, magnetic fields, and the MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, explosion in the immediate post-bounce phase that does not lead to black hole formation. We surmise that fast rotation in the iron core may inhibit, rather than enable, collapsar formation, which requires a large angular momentum not in the core but above it. Variations in the angular momentum distribution of massive stars at core collapse might provide an explanation both for the diversity of Type Ic supernovae/hypernovae and for when they are associated with a GRB. A corollary is that, rather than the progenitor mass, the angular momentum distribution, through its effect on magnetic field amplification, distinguishes these outcomes.

Magnetar's persistent emission above 10 keV was recently discovered thanks to the imaging capabilities of the IBIS coded mask telescope on board the INTEGRAL satellite. The only two sources that show some degree of long term variability are SGR 1806-20 and 1RXS J170849.0-400910. We find some indications that variability of these hard tails could be the driver of the spectral variability measured in these sources below 10 keV. In addition we report for the first time the detection at 2.8 sigma level of pulsations in the hard X-ray tail of SGR 1806-20.

We present the first results of our VLT observation campaign of the Central Compact Objects (CCOs) in SNRs RX J085201.4-461753 (Vela Jr), 1E 1648-5051 (RCW 103) and RX J171328.4-394955 (G347.3-0.5). For Vela Jr., we found that the source is embedded in a compact optical nebulosity, possibly a bow-shock or a photo-ionization nebula, and we identified a candidate IR counterpart to the CCO. For RCW 103, we found no convincing evidence neither for 6 hrs IR modulation nor for variability on any time scale from the proposed counterpart, as well as for the other candidates close to the revised Chandra position. For G347.3-0.5, we identified few possible IR counterparts but none of them is apparently associated with the CCO.

We report on our Spitzer observations of the anomalous X-ray pulsar 4U 0142+61, made following a large X-ray burst that occurred on 2007 February 7. To search for mid-infrared flux variations, four imaging observations were carried out at 4.5 and 8.0 $\mu$m with the Infrared Array Camera from February 14 to 21. No significant flux variations were detected, and the average fluxes were 32.1$\pm$2.0 $\mu$Jy at 4.5 $\mu$m and 59.8$\pm8.5$ $\mu$Jy at 8.0 $\mu$m, consistent with those obtained in 2005. The non-detection of variability is interesting in light of reported rapid variability from this source in the near-infrared, but consistent with the fact that the source already went back to its quiescent state before our observations began, as indicated by contemporaneous X-ray observations. In order to understand the origin of the near-infrared variability, frequent, simultaneous multi-wavelength observations are needed.

The problem of reconstruction of the pulsar magnetospheres nearby the light cylinder surface is studied. It is shownthat on the basis of of the Euler, continuity and induction equations, there is the possibility of parametrically excited rotational energy pumping process into drift modes is shown. As a result a toroidal component of magnetic field increases very rapidly, the increment of which is analyzed for typical magnetospheric plasma parameters. The presented hydrodynamic instability of the generated mode is specific in a sense that the feedback of the instability on particles is insignificant. We analytically investigate dynamics of the reconstruction of a magnetosphere, from generation of the toroidal component of magnetic field up to transformation of field lines into a shape of the Archimedes spiral, when plasma particles do not experience any forces caused by pulsar magnetic field and motion of particles is characterized by the so called the force-free regime. As a result generation of the toroidal component comes to the end and the quasi stable state of the pulsar wind is formed.

We present phase resolved optical photometry and spectroscopy of the accreting millisecond pulsar HETE J1900.1-2455. Our R-band light curves exhibit a sinusoidal modulation, at close to the orbital period, which we initially attributed to X-ray heating of the irradiated face of the secondary star. However, further analysis reveals that the source of the modulation is more likely due to superhumps caused by a precessing accretion disc. Doppler tomography of a broad Halpha emission line reveals an emission ring, consistent with that expected from an accretion disc. Using the velocity of the emission ring as an estimate for the projected outer disc velocity, we constrain the maximum projected velocity of the secondary to be 200 km/s, placing a lower limit of 0.05 Msun on the secondary mass. For a 1.4 Msun primary, this implies that the orbital inclination is low, < 20 degrees. Utilising the observed relationship between the secondary mass and orbital period in short period cataclysmic variables, we estimate the secondary mass to be ~0.085 Msun, which implies an upper limit of ~2.4 Msun for the primary mass.

The soft gamma-ray repeater SGR 1806-20 has been attracting a lot of attention owing to the fact that in December 2004 it emitted the most powerful giant flare ever observed. Here we present the results of the first Suzaku observation of SGR 1806-20, that seems to have reached a state characterized by a flux close to the pre-flare level and by a relatively soft spectrum. Despite this, the source remained quite active, as testified by several short bursts observed by Suzaku. We discuss the broadband spectral properties of SGR 1806-20 in the context of the magnetar model, considering its recent theoretical developments.

The deconfinement phase transition will lead to the release of latent heat during spins down of neutron stars if the transition is the first-order one.We have investigated the thermal evolution of neutron stars undergoing such deconfinement phase transition. The results show that neutron stars may be heated to higher temperature.This feature could be particularly interesting for high temperature of low-magnetic field millisecond pulsar at late stage.

In 2003 a previously unpulsed Einstein and ROSAT source cataloged as soft and dim (Lx of few 10^33 ergs) thermal emitting object, namely XTE J1810-197, was identified as the first unambiguous transient Anomalous X-ray Pulsar. Two years later this source was also found to be a bright highly polarized transient radio pulsar, a unique property among both AXPs and radio pulsars. In September 2006 the Swift Burst Alert Telescope (BAT) detected an intense burst from the candidate AXP CXOU J164710.2-455216, which entered in an outburst state reaching a peak emission of at least a factor of 300 higher than quiescence. Here, we briefly outline the recent results concerning the outburst phenomena observed in these two AXPs. In particular, XTE J1810-197 has probed to be a unique laboratory to monitor the timing and spectral properties of a cooling/fading AXP, while new important information have been inferred from X-ray and radio band simultaneous observations. CXOU J164710.2-455216 has been monitored in X-rays and radio bands since the very beginning of its outbursting state allowing us to cover the first phases of the outburst and to study the timing and spectral behavior during the first months.

Many astrophysicists believe that Anomalous X-Ray Pulsars (AXP), Soft Gamma-Ray Repeaters (SGR), Rotational Radio Transients (RRAT), Compact Central Objects (CCO), and X-Ray Dim Isolated Neutron Stars (XDINS) belong to different classes of anomalous objects with neutron stars as the central bodies inducing all their observable peculiarities. We have shown earlier (I.F.Malov and G.Z.Machabeli, Astron. Astrophys. Trans. 25, 7, 2006) that AXPs and SGRs could be described by the drift model in the framework of preposition on usual properties of the central neutron star (rotation periods P ~ 0.1 - 1 sec, surface magnetic fields B ~ 10^11 - 10^13 G). Here we shall try to show that some differences of considered sources will be explained by their geometry (particularly, by the angle BETA between their rotation and magnetic axes). If BETA < 10 deg. (the aligned rotator) the drift waves at the outer layers of the neutron star magnetosphere should play a key role in the observable periodicity. For large values of BETA (the case of the nearly orthogonal rotator) an accretion from the surrounding medium (for example, from the relic disk) can cause some modulation and transient events in received radiation.

We present a spatially resolved spectroscopic study of the thermal composite supernova remnant Kes 27 with Chandra. The X-ray spectrum of Kes 27 is characterized by K lines from Mg, Si, S, Ar, and Ca. The X-ray emitting gas is found to be enriched in sulphur and calcium. The broadband and tri-color images show two incomplete shell-like features in the northeastern half and brightness fading with increasing radius in the southwest. There are over 30 unresolved sources within the remnant. None show characteristics typical of young neutron stars. The maximum diffuse X-ray intensity coincides with a radio bright region along the eastern border. In general, gas in the inner region is at higher temperature and emission is brighter than from the outer region. The gas in the remnant appears to approach ionization equilibrium. The overall morphology can be explained by the evolution of the remnant in an ambient medium with a density enhancement from west to east. We suggest that the remnant was born in a pre-existing cavity and that the inner bright emission is due to the reflection of the initial shock from the dense cavity wall. This scenario may provide a new candidate mechanism for the X-ray morphology of other thermal composite supernova remnants.

We report new results obtained from multi-frequency observations of PSR B0826-34 with the Giant Metrewave Radio Telescope (GMRT). (1) We find no evidence of weak emission during the typical long null state of this pulsar, simultaneously at 303 and 610 MHz, as well as individually at 157, 325, 610 and 1060 MHz at separate epochs. Our limit of non-detection is at ~ 1% or better of the peak of the active state profile, and corresponds to ~ 2 mJy at 610 MHz. (2) Significant correlation in the total intensity of the individual pulses between 303 and 610 MHz is reported from the simultaneous dual frequency observations, which is indicative of the broadband nature of the emission. We also report correlation between total energy in the main pulse and inter-pulse region from the high sensitivity single frequency observations at 610 and 1060 MHz. (3) Though we find the drift pattern to be very similar in the simultaneous 303 and 610 MHz data, we observe that the drift band separation (P2) evolves significantly between these two frequencies, and in a manner opposite to the average profile evolution. In addition, we confirm the dependence of P2 on pulse longitude at 303 MHz and find indications for the same at 610 MHz. We also present results for subpulse width at different frequencies, and as well as a function of pulse longitude. (4) As a natural out-come of the simultaneous dual frequency observations, we obtain an accurate DM value, equal to 52.2(6) pc/cc, for this pulsar.

In the standard supernova picture, type Ib/c and type II supernovae are powered by the potential energy released in the collapse of the core of a massive star. In studying supernovae, we primarily focus on the ejecta that makes it beyond the potential well of the collapsed core. But, as we shall show in this paper, in most supernova explosions, a tenth of a solar mass or more of the ejecta is decelerated enough that it does not escape the potential well of that compact object. This material falls back onto the proto-neutron star within the first 10-15 seconds after the launch of the explosion, releasing more than 1e52erg of additional potential energy. Most of this energy is emitted in the form of neutrinos and we must understand this fallback neutrino emission if we are to use neutrino observations to study the behavior of matter at high densities. Here we present both a 1-dimensional study of fallback using energy-injected, supernova explosions and a first study of neutrino emission from fallback using a suite of 2-dimensional simulations.

A number of X-ray binaries (neutron stars or black holes accreting from a companion star) have such short orbital periods that ordinary, hydrogen rich, stars do not fit in. Instead the mass-losing star must be a compact, evolved star, leading to the transfer of hydrogen deficient material to the neutron star. I discuss the current knowledge of these objects, with focus on optical spectroscopy.

We present spectral measurements made in the soft (20-100 keV) gamma-ray band of the region containing the composite supernova remnant G11.2-0.3 and its associated pulsar PSR J1811-1925. Analysis of INTEGRAL/IBIS data allows characterisation of the system above 10 keV. The IBIS spectrum is best fitted by a power law having photon index of 1.8^{+0.4}_{-0.3} and a 20-100 keV flux of 1.5E{-11} erg/cm^2/s. Analysis of archival Chandra data over different energy bands rules out the supernova shell as the site of the soft gamma-ray emission while broad band (1-200 keV) spectral analysis strongly indicates that the INTEGRAL/IBIS photons originate in the central zone of the system which contains both the pulsar and its nebula. The composite X-ray and soft gamma-ray spectrum indicates that the pulsar provides around half of the emission seen in the soft gamma-ray domain; its spectrum is hard with no sign of a cut off up to at least 80 keV. The other half of the emission above 10 keV comes from the PWN; with a power law slope of 1.7 its spectrum is softer than that of the pulsar. From the IBIS/ISGRI mosaics we are able to derive 2 sigma upper limits for the 20-100 keV flux from the location of the nearby TeV source HESS J1809-193 to be 4.8E{-12} erg/cm^2/s. We have also examined the likelihood of an association between PSR J1811-1925 and HESS J1809-193. Although PSR J1811-1925 is the most energetic pulsar in the region, the only one detected above 10 keV and thus a possible source of energy to fuel the TeV fluxes, there is no morphological evidence to support this pairing, making it an unlikely counterpart.

We analyzed the available LIGO data coincident with GRB 070201, a short duration hard spectrum gamma-ray burst whose electromagnetically determined sky position is coincident with the spiral arms of the Andromeda galaxy (M31). Possible progenitors of such short hard GRBs include mergers of neutron stars or a neutron star and black hole, or soft gamma-ray repeater (SGR) flares. These events can be accompanied by gravitational-wave emission. No plausible gravitational wave candidates were found within a 180 s long window around the time of GRB 070201. This result implies that a compact binary progenitor of GRB 070201, with masses in the range 1 M_sun < m_1 < 3 M_sun and 1 M_sun < m_2 < 40 M_sun, located in M31 is excluded at >99% confidence. Indeed, if GRB 070201 were caused by a binary neutron star merger, we find that D < 3.5 Mpc is excluded, assuming random inclination, at 90% confidence. The result also implies that an unmodeled gravitational wave burst from GRB 070201 most probably emitted less than 4.4 x 10^(-4) M_sun c^2 (7.9 x 10^(50) ergs) in any 100 ms long period within the signal region if the source was in M31 and radiated isotropically at the same frequency as LIGO's peak sensitivity (f ~ 150 Hz). This upper limit does not exclude current models of SGRs at the M31 distance.

Six glitches have been recently observed in the rotational frequency of the young pulsar PSR B1737$-$30 (J1740$-$3015) using the 25-m Nanshan telescope of Urumqi Observatory. With a total of 20 glitches in 20 years, it is one of the most frequently glitching pulsars of the about 1750 known pulsars. Glitch amplitudes are very variable with fractional increases in rotation rate ranging from 10^{-9} to 10^{-6}. Inter-glitch intervals are also very variable, but no relationship is observed between interval and the size of the preceding glitch. There is a persistent increase in |\dot\nu|, opposite in sign to that expected from slowdown with a positive braking index, which may result from changes in the effective magnetic dipole moment of the star during the glitch.

A new method is suggested for search of the axions constituting Dark Matter that utilizes observations of neutron stars (NS) in radio-frequency region. It uses the conversion of axions to photons in strong magnetic fields of NS. The observations of Magnificent Seven objects are proposed. Whether the conversion takes place, the radio spectrum of the object would have a very distinctive feature -- a narrow spike at a frequency corresponding to the rest mass of the axions. For example, if the coupling constant of the photon-axion interaction is M=10^10 GeV, the density of DM axions is \rho=10^(-24)*g*cm^(-3) and the NS is located at a distance of 300 pc from the Solar system, then the flux density of excess signal for axions with the rest mass of 1 \mu eV will be as large as several mJy at the frequency 240 MHz in the bandwidth 0.5 MHz.

From a uniform analysis of a large (8.5 Ms) Rossi X-ray Timing Explorer data set of Low Mass X-ray Binaries, we present a complete identification of all the variability components in the power spectra of black holes in their canonical states. It is based on gradual frequency shifts of the components observed between states, and uses a previous identification in the black hole low hard state as a starting point. It is supported by correlations between the frequencies in agreement with those previously found to hold for black hole and neutron stars. Similar variability components are observed in neutron stars and black holes (only the component observed at the highest frequencies is different) which therefore cannot depend on source-specific characteristics such as the magnetic field or surface of the neutron star or spin of the black hole. As the same variability components are also observed across the jet-line the X-ray variability cannot originate from the outer-jet but is most likely produced in either the disk or the corona. We use the identification to directly compare the difference in strength of the black hole and neutron star variability and find these can be attributed to differences in frequency and strength of high frequency features, and do not require the absence of any components. Black holes attain their highest frequencies (in the hard-intermediate and very-high states) at a level a factor ~6 below the highest frequencies attained by the corresponding neutron star components, which can be related to the mass difference between the compact objects in these systems.

Detection of paired kilohertz quasi-periodic oscillations (kHz QPOs) in the X-ray emission of a compact object is compelling evidence that the object is an accreting neutron star. In many neutron stars, the stellar spin rate is equal or roughly equal to Delta-nu, the frequency separation of the QPO pair, or to 2Delta-nu. Hence, if the mechanism that produces the kilohertz QPOs is similar in all stars, measurement of Delta-nu can provide an estimate of the star's spin rate. The involvement of the stellar spin in producing Delta-nu indicates that the magnetic fields of these stars are dynamically important.
We focus here on the implications of the paired kHz QPOs recently discovered in the low-mass X-ray binary (LMXB) system Cir X-1 (Boutloukos et al. 2006). The kHz QPOs discovered in Cir X-1 are generally similar to those seen in other stars, establishing that the compact object in the Cir X-1 system is a neutron star. However, the frequency nu-u of its upper kHz QPO is up to a factor of three smaller than is typical, and Delta-nu varies by about a factor 2 (167 Hz, the largest variation so far observed). Periodic oscillations have not yet been detected from Cir X-1, so its spin rate has not yet been measured directly. The low values of nu-u and the large variation of Delta-nu challenge current models of the generation of kHz QPOs. Improving our understanding of Cir X-1 will improve our knowledge of the spin rates and magnetic fields of all neutron stars.

During the first 40 s after their birth, proto-neutron stars are expected to be subject to at least two types of instability: the convective instability and the neutron-finger one. Both instabilities involve convective motions and hence can trigger dynamo actions which may be responsible for the large magnetic fields in neutron stars and magnetars. We have solved the mean-field induction equation in a simplified one-dimensional model of both the convective and the neutron-finger instability zones. Although very idealized, the model includes the nonlinearities introduced by the feedback processes which tend to saturate the growth of the magnetic field (alpha-quenching) and suppress its turbulent diffusion (eta-quenching). The possibility of a dynamo action is studied within a dynamical model of turbulent diffusivity where the boundary of the unstable zone is allowed to move. We show that the dynamo action can be operative and that the amplification of the magnetic field can still be very effective. Furthermore, we confirm the existence of a critical spin-period, below which the dynamo is always excited independently of the degree of differential rotation, and whose value is related to the size of the neutron-finger instability zone. Finally we provide a relation for the intensity of the final field as a function of the spin of the star and of its differential rotation. Although they were obtained by using a toy model, we expect that our results are able to capture the qualitative and asymptotic behaviour of a mean-field dynamo action developing in the neutron-finger instability zone. Overall, we find that such a dynamo is very efficient in producing magnetic fields well above equipartition and thus that it could represent a possible explanation for the large surface magnetic fields observed in neutron stars.

We show that if early proton-rich supernova wind ejecta are responsible for the solar production of Mo92 and Mo94 then the proton separation energy for Rh93 is S_p = 1.64+/-0.1 MeV. The current experimentally-based best estimate for this proton separation energy is S_p = 2.05 MeV, with an uncertainty of 0.5 MeV (Audi, Wapstra & Thibault 2003), a factor of five larger than the uncertainty inferred here. A more accurate experimental determination of this separation energy could provide strong circumstantial evidence that neutrino-irradiated winds from a nascent neutron star are responsible for producing the "light p" nuclei, which include Mo92 and Mo94.

We present an analysis of the X-ray variability of three symbiotic X-ray binaries, GX 1+4, 4U 1700+24, and 4U 1954+31, using observations made with the Swift Burst Alert Telescope (BAT) and the Rossi X-ray Timing Explorer (RXTE) All-Sky Monitor (ASM). Observations of 4U 1954+31 with the Swift BAT show modulation at a period near 5 hours. Models to explain this modulation are discussed including the presence of an exceptionally slow X-ray pulsar in the system and accretion instabilities. We conclude that the most likely interpretation is that 4U 1954+31 contains one of the slowest known X-ray pulsars. Unlike 4U 1954+31, neither GX 1+4 nor 4U 1700+24 show any evidence for modulation on a timescale of hours. An analysis of the RXTE ASM light curves of GX 1+4, 4U 1700+24, and 4U 1954+31 does not show the presence of periodic modulation in any source, although there is considerable variability on long timescales for all three sources. There is no modulation in GX 1+4 on either the optical 1161 day orbital period or a previously reported 304 day X-ray period. For 4U 1700+24 we do not confirm the 404 day X-ray period previously proposed for this source from a shorter duration ASM light curve. We conclude that all three sources have substantial low-frequency noise in their power spectra that may give the appearance of periodic modulation if this noise is not properly accounted for, particularly if short duration light curves are examined.

Black holes are popping up all over the place: in compact binary X-ray sources and GRBs, in quasars, AGNs and the cores of all bulge galaxies, in binary black holes and binary black hole-neutron stars, and maybe even in the LHC! Black holes are strong-field objects governed by Einstein's equations of general relativity. Hence general relativistic, numerical simulations of dynamical phenomena involving black holes may help reveal ways in which black holes can form, grow and be detected in the universe. To convey the state-of-the art, we summarize several representative simulations here, including the collapse of a hypermassive neutron star to a black hole following the merger of a binary neutron star, the magnetorotational collapse of a massive star to a black hole, and the formation and growth of supermassive black hole seeds by relativistic MHD accretion in the early universe.

Our new study of the two central compact object pulsars, PSR J1210-5226 (P = 424 ms) and PSR J1852+0040 (P = 105 ms), leads us to conclude that a weak natal magnetic field shaped their unique observational properties. In the dipole spin-down formalism, the 2-sigma upper limits on their period derivatives, < 2E-16 for both pulsars, implies surface magnetic field strengths of B_s < 3E11 G and spin periods at birth equal to their present periods to three significant digits. Their X-ray luminosities exceed their respective spin-down luminosities, implying that their thermal spectra are derived from residual cooling and perhaps partly from accretion of supernova debris. For sufficiently weak magnetic fields an accretion disk can penetrate the light cylinder and interact with the magnetosphere while resulting torques on the neutron star remain within the observed limits. We propose the following as the origin of radio-quiet CCOs: the magnetic field, derived from a turbulent dynamo, is weaker if the NS is formed spinning slowly, which enables it to accrete SN debris. Accretion excludes neutron stars born with both B_s < 1E11 G and P > 0.1 s from radio pulsar surveys, where such weak fields are not encountered except among very old (> 40 Myr) or recycled pulsars. We predict that these birth properties are common, and may be attributes of the youngest detected neutron star, the CCO in Cassiopeia A, as well as an undetected infant neutron star in the SN 1987A remnant. In view of the far-infrared light echo discovered around Cas A and attributed to an SGR-like outburst, it is especially important to determine via timing whether Cas A hosts a magnetar or not. If not a magnetar, the Cas A NS may instead have undergone a one-time phase transition (corequake) that powered the light echo.

We describe in this paper a new binary millisecond pulsar, PSR J1738+0333. Using Arecibo, we have achieved good timing accuracy for this object, about 220 ns for 1-hour integrations over 100 MHz. This allowed us to measure a precise proper motion, parallax and orbital parameters for this system. We highlight the system's potential for constraining alternative theories of gravitation.

PSR J0514-4002A is a 5-ms pulsar is located in the globular cluster NGC 1851; it belongs to a highly eccentric (e = 0.888) binary system. It is one of the earliest known examples of a numerous and fast-growing class of eccentric binary MSPs recently discovered in globular clusters. Using the GBT, we have obtained a phase-coherent timing solution for the pulsar, which includes a measurement of the rate of advance of periastron: 0.01289(4) degrees per year, which if due completely to general relativity, implies a total system mass of 2.453(14) solar masses. We also derive m_p < 1.5 solar masses and m_c > 0.96 solar masses. The companion is likely to be a massive white dwarf star.

We report the discovery of five new millisecond pulsars in the globular cluster NGC 6440 and three new ones in NGC 6441; each cluster has one previously known pulsar. Four of the new pulsars are found in binary systems. One of the new pulsars, PSR J1748-2021B in NGC 6440, is notable for its eccentric (e = 0.57) and wide (P_b = 20.5 days) orbit. If the rate of advance of periastron is due solely to general relativity, we can estimate of the total mass of this binary system: 2.92 +/- 0.20 solar masses. This would imply an anomalously large mass for this pulsar, which could introduce important constraints in the study of the equation of state for cold neutron matter.

We report on optical and infrared observations of the anomalous X-ray pulsar (AXP) 1E 1048.1-5937, made during its ongoing X-ray flare which started in 2007 March. We detected the source in the optical I and near-infrared K_s bands in two ground-based observations and obtained deep flux upper limits from four observations, including one with the Spitzer Space Telescope at 4.5 and 8.0 microns. The detections indicate that the source was approximately 1.3--1.6 magnitudes brighter than in 2003 April/June, when it was at the tail of a previous similar X-ray flare. Comparing the optical/near-infrared fluxes to the X-ray flux at the two epochs, we find that they are correlated. Similar related flux variations have been seen in two other AXPs during their X-ray outbursts, suggesting common behavior for large X-ray flux variation events in AXPs. The Spitzer flux limits are sufficiently deep that we can exclude mid-infrared emission similar to that from the AXP 4U 0142+61, which has been interpreted as arising from a dust disk around the AXP. The optical/near-infrared emission probably arises from the magnetosphere of 1E 1048.1-5937, given the similar flux spectra from the source and 4U 0142+61 and the magnetospheric origin suggested for the latter. The flux spectrum similarity challenges the dust disk model proposed for 4U 0142+61, since in the model, $K_s$ band emission primarily arises from the disk, not from the magnetosphere.

In recent years, X-ray observations with Chandra and XMM-Newton have significantly increased our understanding of rotation-powered (radio) millisecond pulsars (MSPs). Deep Chandra studies of several globular clusters have detected X-ray counterparts to a host of MSPs, including 19 in 47 Tuc alone. These surveys have revealed that most MSPs exhibit thermal emission from their heated magnetic polar caps. Realistic models of this thermal X-ray emission have provided important insight into the basic physics of pulsars and neutron stars. In addition, intrabinary shock X-ray radiation observed in ``black-widow'' and peculiar globular cluster ``exchanged'' binary MSPs give interesting insight into MSP winds and relativistic shock. Thus, the X-ray band contains valuable information regarding the basic properties of MSPs that are not accesible by radio timing observations.

Recent astrophysical observations of neutron stars and heavy-ion data are confronted with our present understanding of the equation of state of dense hadronic matter. Emphasis is put on the possible role of the presence of hyperons in the interior of compact stars. We argue that data from low-mass pulsars provide an important cross-check between high-density astrophysics and heavy-ion physics.

(abridged) The physics of the pulsar magnetosphere remains poorly constrained by observations. Little is known about their emission mechanism. Large vacuum gaps probably exist, and a non-neutral plasma partially fills the neutron star surroundings to form an electrosphere. We showed that the differentially rotating equatorial disk in the pulsar's electrosphere is diocotron unstable and that it tends to stabilise when relativistic effects are included. However, when approaching the light cylinder, particle inertia becomes significant and the electric drift approximation is violated. In this paper, we study the most general instability, i.e. by including particle inertia effects, as well as relativistic motions. This general non-neutral plasma instability is called the magnetron instability. We linearise the coupled relativistic cold-fluid and Maxwell equations. The non-linear eigenvalue problem for the perturbed azimuthal electric field component is solved numerically. The spectrum of the magnetron instability in a non-neutral plasma column confined between two cylindrically conducting walls is computed for several cylindrical configurations. For a pulsar electrosphere, no outer wall exists. In this case, we allow for electromagnetic wave emission propagating to infinity. When the self-field induced by the plasma becomes significant, it can first increase the growth rate of the magnetron instability. However, equilibrium solutions are only possible when the self-electric field, measured by the parameter $s_{\rm e}$ and tending to disrupt the plasma configuration, is bounded to an upper limit, $s_{\rm e,max}$. For $s_{\rm e}$ close to but smaller than this value $s_{\rm e,max}$, the instability becomes weaker or can be suppressed as was the case in the diocotron regime.

For a variety of fully relativistic polytropic neutron star models we calculate the star's tidal Love number k2. Most realistic equations of state for neutron stars can be approximated as a polytrope with an effective index n~0.5-1.0. The equilibrium stellar model is obtained by numerical integration of the Tolman-Oppenheimer-Volkhov equations. We calculate the linear l=2 static perturbations to the Schwarzschild spacetime following the method of Thorne and Campolattaro. Combining the perturbed Einstein equations into a single second order differential equation for the perturbation to the metric coefficient g_tt, and matching the exterior solution to the asymptotic expansion of the metric in the star's local asymptotic rest frame gives the Love number. Our results agree well with the Newtonian results in the weak field limit. The fully relativistic values differ from the Newtonian values by up to ~24%. The Love number is potentially measurable in gravitational wave signals from inspiralling binary neutron stars.

We revise the distance estimate to the supernova remnant (SNR) G109.1-1.0 (CTB 109) and its associated anomalous X-ray pulsar (AXP) 1E 2259+586 by re-analyzing both 408 and 1420 MHz continuum images, 21 cm HI absorption/emission spectra and $^{12}$CO emission spectra of CTB 109 and an adjacent compact HII region Sh 152, using the data from the Canadian Galactic Plane Survey. The HI images show that CTB 109 is morphologically located inside a shell in HI emission in the velocity range -61 to -70 km/s. Such a shell is likely blown by a stellar wind. The relation is supported by newly detected HI absorption features between $\sim$ -61 km/s and -66 km/s, which are in the direction of the bright radio continuum emission within the north shell of CTB 109. We therefore give a distance of $\sim$ 6 kpc to CTB 109, which is consistent with the distance of AXP 1E 2259+586 derived by the red clump stars method. The distance of $\sim$ 6 kpc to the SNR/AXP system argues against a SNR/CO cloud interaction and instead favors a SNR/HI cloud interaction. The large distance argues for an anomalously large explosion energy which might be expected by a magnetar (AXP 1E 2259+586) harboring at the center of CTB 109.

We have identified three possible ways in which future XMM-Newton observations can provide significant constraints on the equation of state of neutron stars. First, using a long observation of the neutron star X-ray transient CenX-4 in quiescence one can use the RGS spectrum to constrain the interstellar extinction to the source. This removes this parameter from the X-ray spectral fitting of the pn and MOS spectra and allows us to investigate whether the variability observed in the quiescent X-ray spectrum of this source is due to variations in the soft thermal spectral component or variations in the power law spectral component coupled with variations in N_H. This will test whether the soft thermal spectral component can indeed be due to the hot thermal glow of the neutron star. Potentially such an observation could also reveal redshifted spectral lines from the neutron star surface. Second, XMM-Newton observations of radius expansion type I X-ray bursts might reveal redshifted absorption lines from the surface of the neutron star. Third, XMM-Newton observations of eclipsing quiescent low-mass X-ray binaries provide the eclipse duration. With this the system inclination can be determined accurately. The inclination determined from the X-ray eclipse duration in quiescence, the rotational velocity of the companion star and the semi-amplitude of the radial velocity curve determined through optical spectroscopy, yield the neutron star mass.

Chandra or XMM-Newton observations of quiescent low-mass X-ray binaries can provide important constraints on the equation of state of neutron stars. The mass and radius of the neutron star can potentially be determined from fitting a neutron star atmosphere model to the observed X-ray spectrum. For a radius measurement it is of critical importance that the distance to the source is well constrained since the fractional uncertainty in the radius is at least as large as the fractional uncertainty in the distance. Uncertainties in modelling the neutron star atmosphere remain. At this stage it is not yet clear if the soft thermal component in the spectra of many quiescent X-ray binaries is variable on timescales too short to be accommodated by the cooling neutron star scenario. This can be tested with a long XMM-Newton observation of the neutron star X-ray transient CenX-4 in quiescence. With such an observation one can use the Reflection Grating Spectrometer spectrum to constrain the interstellar extinction to the source. This removes this parameter from the X-ray spectral fitting of the EPIC pn and MOS spectra and allows one to investigate whether the variability observed in the quiescent X-ray spectrum of this source is due to variations in the soft thermal spectral component or variations in the power law spectral component coupled with variations in N_H. This will test whether the soft thermal component can indeed be due to the hot thermal glow of the neutron star. Irrespective of the outcome of such a study, the observed cooling in quiescence in sources for which the crust is significantly out of thermal equilibrium with the core due to a prolonged outburst, such as KS 1731-260, seem excellent candidates for mass and radius determinations through modelling the observed X-rays with a neutron star atmosphere model.

We present an updated timing solution and an analysis of the profile evolution - including precession and beam shape - of the young, relativistic binary pulsar J1906+0746. The 144-ms pulsar, in a 3.98-hour orbit with eccentricity 0.085 (Lorimer et al. 2006), was initially discovered during the early stages of the ALFA (Arecibo L-band Feed Array) pulsar survey (Cordes et al. 2006) using the 305-metre Arecibo telescope and was subsequently found in archival Parkes Multibeam Survey data. We have since been regularly monitoring the system using the Arecibo and Green Bank telescopes, and include data from the Jodrell Bank, Parkes, Nancay and Westerbork telescopes. The nature of the binary companion will also be discussed based on improved estimates of the total and companion masses obtained from the updated timing solution.

The hard X-ray detector (HXD) on board the X-ray satellite Suzaku is designed to have a good timing capability with a 61 $\mu$s time resolution. In addition to detailed descriptions of the HXD timing system, results of in-orbit timing calibration and performance of the HXD are summarized. The relative accuracy of time measurements of the HXD event was confirmed to have an accuracy of $1.9\times 10^{-9}$ s s$^{-1}$ per day, and the absolute timing was confirmed to be accurate to 360 $\mu$s or better. The results were achieved mainly through observations of the Crab pulsar, including simultaneous ones with RXTE, INTEGRAL, and Swift.

We use a modified outer gap model to study the multi-frequency phase-resolved spectra of the Crab pulsar. The emissions from both poles contribute to the light curve and the phase-resolved spectra. Using the synchrotron self-Compton mechanism and by considering the incomplete conversion of curvature photons into secondary pairs, the observed phase-averaged spectrum from 100 eV - 10 GeV can be explained very well. The predicted phase-resolved spectra can match the observed data reasonably well, too. We find that the emission from the north pole mainly contributes to Leading Wing 1. The emissions in the remaining phases are mainly dominated by the south pole. The widening of the azimuthal extension of the outer gap explains Trailing Wing 2. The complicated phase-resolved spectra for the phases between the two peaks, namely Trailing Wing 1, Bridge and Leading Wing 2, strongly suggest that there are at least two well-separated emission regions with multiple emission mechanisms, i.e. synchrotron radiation, inverse Compton scattering and curvature radiation. Our best fit results indicate that there may exist some asymmetry between the south and the north poles. Our model predictions can be examined by GLAST.

We present a detailed analysis of observations of the high mass X-ray binary Cen X-3 spanning two consecutive binary orbits performed with the RXTE satellite in early March 1997. The PCA and HEXTE light curves both show a clear reduction in count rate after mid-orbit for both binary revolutions. We therefore analyze two broad band spectra for each orbit, before and after mid-orbit. Consistent with earlier observations these four joint PCA and HEXTE spectra can be well described using a phenomenological pulsar continuum model, including an iron emission line and a cyclotron resonance scattering feature. While no strong spectral variations were detected, the second half of orbit 2 shows a tendency toward being softer and more strongly absorbed. In order to follow the orbital phase-dependent evolution of the spectrum in greater detail, we model spectra for shorter exposures, confirming that most spectral parameters show either a gradual or sudden change for the second half of the second orbit. A comparison with a simple wind model indicates the existence of an accretion wake in this system. We also present and discuss high resolution pulse profiles for several different energy bands, as well as their hardness ratios. PCA and HEXTE spectra were created for 24 phase bins and fitted using the same model as in the phase averaged case. Systematic pulse phase-dependent variations of several continuum and cyclotron line parameters were detected, most notably a significant increase of the cyclotron line energy during the early rise of the main peak, followed by a gradual decrease. We show that applying a simple dipole model for the magnetic field is not sufficient to describe our data.

Sterile neutrinos with keV masses can constitute all or part of the cosmological dark matter. The electroweak-singlet fermions, which are usually introduced to explain the masses of active neutrinos, need not be heavier than the electroweak scale; if one of them has a keV-scale mass, it can be the dark-matter particle, and it can also explain the observed pulsar kicks. The relic sterile neutrinos could be produced by several different mechanisms. If they originate primarily from the Higgs decays at temperatures of the order of 100 GeV, the resulting dark matter is much ``colder'' than the warm dark matter produced in neutrino oscillations. The signature of this form of dark matter is the spectral line from the two-body decay, which can be detected by the X-ray telescopes. The same X-rays can have other observable manifestations, in particular, though their effects on the formation of the first stars.

In this contribution we discuss how neutron stars are produced and retained in globular clusters, outlining the most important dynamical channels and evolutionary events that affect thepopulation of mass-transferring binaries with neutron stars and result in the formation of recycled pulsars. We confirm the importance of electron-capture supernovae in globular clusters as the major supplier of retained neutron stars.By comparing the observed millisecond pulsar population and the results obtained from simulations, we discuss several constraints on the evolution of mass-transferring systems.In particular, we find that in our cluster model the following mass-gaining events create populations of MSPs that do not match the observations (with respect to binary periods and companion masses or the number of produced systems) and therefore likely do not lead to NSs spun up to millisecond periods: (i) accretion during a common envelope event with a NS formed through accretion-induced collapse, and (ii) mass transfer from a WD donor. By restricting ourselves to the evolutionary and dynamical paths that most likely lead to neutron star recycling, we obtain good agreement between our models and the numbers and characteristics of observed millisecond pulsars in the clusters Terzan 5 and 47 Tuc.

We present a three-stage model for a long GRB inner engine to explain the prompt gamma ray emission, and interpret recent Swift satellite observations of early X-ray afterglow plateaus followed by a sharp drop off or a shallow power law decay. The three stages involves a neutron star phase, a quark star (QS) and a black hole phase as described in Staff et al. (2007). We find that the QS stage allows for more energy to be extracted from neutron star to QS conversion as well as from ensuing accretion onto the QS. The QS accretion phase naturally extends the engine activity and can account for both the prompt emission and irregular early X-ray afterglow activity. Following the accretion phase, the QS can spin-down by emission of a baryon-free outflow. The magnetar-like magnetic field strengths resulting from the NS to QS transition provide enough spin-down energy, for the correct amount of time, to account for the plateau in the X-ray afterglow. In our model, a sharp drop-off following the plateau occurs when the QS collapses to a BH during the spin-down, thus shutting-off the secondary outflow. We applied our model to GRB 070110 and GRB 060607A and found that we can consistently account for the energetics and duration during the prompt and plateau phases.

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical open issues. Here, after scrutinizing the ATNF Catalog searching for pulsar distances and proper motions, we verify that for an ATNF sample of 212 galactic run away pulsars (RAPs) that were flung at very high speed and undergo large displacements, some gravitational-wave (GW) signals produced by such large accelerations appear to be detectable, after calibration against the Advanced LIGO (LII). Motivated by this insight, we address the pulsar kick at birth or {\sl short rise fling} with the theory for phenomena of emission of GW with memory. We show that during the short rise fling each run away pulsar (RAP) should have generated a GW signal with characteristic amplitude and frequency potentially detectable by current GW interferometers. An efficiency parameter quantifies the use of the rise time kinetic energy, which is estimated from the linear momentum conservation law applied to the supernova explosion that flings the pulsar. The remaining energy is supposed to be used to make the star to spin. The spin of ATNF pulsars with velocity in the interval 400-500 km s$^{-1}$ is compared. The resulting difference suggests that other mechanisms (like differential rotation, magnetic breaking or magneto-rotational instability) should dissipate part of that energy to produce the observed pulsar spin periods. Meanwhile, the kick phenomenon may also occur in globular and open star clusters at the formation or disruption of very short period compact binary systems where abrupt velocity and acceleration similar to RAPs can be imparted. In this case, pulsar astrometry from micro- to nano-arsec scales may be of much help. In case of a supernova, the RAP GW signal could be a benchmark for the GW signal from the core collapse.

Dynamical interactions that occur between objects in dense stellar systems are particularly important for the question of formation of X-ray binaries. We present results of numerical simulations of 70 globular clusters with different dynamical properties and a total stellar mass of 2*10^7 Msun. We find that in order to retain enough neutron stars to match observations we must assume that NSs can be formed via electron-capture supernovae. Our simulations explain the observed dependence of the number of LMXBs on ``collision number'' as well as the large scatter observed between different globular clusters. For millisecond pulsars, we obtain good agreement between our models and the numbers and characteristics of observed pulsars in the clusters Terzan 5 and 47 Tuc

White dwarfs represent the endpoint of stellar evolution for stars with initial masses between approximately 0.07 msun and 8-10 msun, where msun is the mass of the Sun (more massive stars end their life as either black holes or neutron stars). The theory of stellar evolution predicts that the majority of white dwarfs have a core made of carbon and oxygen, which itself is surrounded by a helium layer and, for ~80 per cent of known white dwarfs, by an additional hydrogen layer. All white dwarfs therefore have been traditionally found to belong to one of two categories: those with a hydrogen-rich atmosphere (the DA spectral type) and those with a helium-rich atmosphere (the non-DAs). Here we report the discovery of several white dwarfs with atmospheres primarily composed of carbon, with little or no trace of hydrogen or helium. Our analysis shows that the atmospheric parameters found for these stars do not fit satisfactorily in any of the currently known theories of post-asymptotic giant branch evolution, although these objects might be the cooler counterpart of the unique and extensively studied PG1159 star H1504+65. These stars, together with H1504+65, might accordingly form a new evolutionary sequence that follow the asymptotic giant branch.

We investigate when the energy that pins a superfluid vortex to the lattice of nuclei in the inner crust of neutron stars can be approximated by the energy that binds the vortex to a single nucleus. Indeed, although the pinning energy is the quantity relevant to the theory of pulsar glitches, so far full quantum calculations have been possible only for the binding energy. Physically, the presence of nearby nuclei can be neglected if the lattice is dilute, namely with nuclei sufficiently distant from each other. We find that the dilute limit is reached only for quite large Wigner-Seitz cells, with radii > 55 fm; these are found only in the outermost low-density regions of the inner crust. We conclude that present quantum calculations do not correspond to the pinning energies in almost the entire inner crust and thus their results are not predictive for the theory of glitches.

White dwarfs represent the endpoint of stellar evolution for stars with initial masses between approximately 0.07 msun and 8-10 msun, where msun is the mass of the Sun (more massive stars end their life as either black holes or neutron stars). The theory of stellar evolution predicts that the majority of white dwarfs have a core made of carbon and oxygen, which itself is surrounded by a helium layer and, for ~80 per cent of known white dwarfs, by an additional hydrogen layer. All white dwarfs therefore have been traditionally found to belong to one of two categories: those with a hydrogen-rich atmosphere (the DA spectral type) and those with a helium-rich atmosphere (the non-DAs). Here we report the discovery of several white dwarfs with atmospheres primarily composed of carbon, with little or no trace of hydrogen or helium. Our analysis shows that the atmospheric parameters found for these stars do not fit satisfactorily in any of the currently known theories of post-asymptotic giant branch evolution, although these objects might be the cooler counterpart of the unique and extensively studied PG1159 star H1504+65. These stars, together with H1504+65, might accordingly form a new evolutionary sequence that follow the asymptotic giant branch.

We investigate when the energy that pins a superfluid vortex to the lattice of nuclei in the inner crust of neutron stars can be approximated by the energy that binds the vortex to a single nucleus. Indeed, although the pinning energy is the quantity relevant to the theory of pulsar glitches, so far full quantum calculations have been possible only for the binding energy. Physically, the presence of nearby nuclei can be neglected if the lattice is dilute, namely with nuclei sufficiently distant from each other. We find that the dilute limit is reached only for quite large Wigner-Seitz cells, with radii > 55 fm; these are found only in the outermost low-density regions of the inner crust. We conclude that present quantum calculations do not correspond to the pinning energies in almost the entire inner crust and thus their results are not predictive for the theory of glitches.

Possible origins of the magnetic fields of neutron stars include inheritance from the main sequence progenitor and dynamo action at some stage of evolution of progenitor. Inheritance is not sufficient to explain the fields of magnetars. Energetic considerations point to differential rotation in the final stages of core collapse process as the most likely source of field generation, at least for magnetars. A runaway phase of exponential growth is needed to achieve sufficient field amplification during relevant phase of core collapse; it can probably be provided by a some form of magnetorotational instability. Once formed in core collapse, the field is in danger of decaying again by magnetic instabilities. The evolution of a magnetic field in a newly formed neutron star is discussed, with emphasis on the existence of stable equilibrium configurations as end products of this evolution, and the role of magnetic helicity in their existence.

We discuss the influence of the cosmological constant $\Lambda$ on the gravitational equations of motion of bodies with arbitrary masses and eventually solve the two-body problem. Observational constraints are derived from measurements of the periastron advance in stellar systems, in particular binary pulsars and the solar system. For the latter we consider also the change in the mean motion due to $\Lambda$. Up to now, Earth and Mars data give the best constraint, $\Lambda \sim 10^{-36} \mathrm{km}^{-2}$. If properly accounting for the gravito-magnetic effect, this upper limit on $\Lambda$ could greatly improve in the near future thanks to new data from planned or already operating space-missions. Dark matter or modifications of the Newtonian inverse-square law in the solar system are discussed as well. Variations in the $1/r^2$ behavior are considered in the form of either a possible Yukawa-like interaction or a modification of gravity of MOND type.

Following our discovery of radio pulsations from the newly recognized Anomalous X-ray Pulsar (AXP) 1E 1547.0-5408, we initiated X-ray monitoring with the Swift X-ray Telescope, and obtained a single target-of-opportunity observation with the Newton X-ray Multi-Mirror Mission (XMM-Newton). In comparison with its historic minimum flux of 3e-13 ergs cm^-2 s^-1, the source was found to be in a record high state, f_X(1-8 keV) = 5e-12 ergs cm^-2 s^-1, or L_X = 1.7e35(d/9 kpc)^2 ergs s^-1, and declining by 25% in 1 month. Extrapolating the decay, we bound the total energy in this outburst to 1e42 < E < 1e43 ergs. The spectra (fitted with a Comptonized blackbody) show that an increase in the temperature and area of a hot region, to 0.5 keV and ~16% of the surface area of the neutron star, respectively, are primarily responsible for its increase in luminosity. The energy, spectrum, and timescale of decay are consistent with a deep crustal heating event, similar to an interpretation of the X-ray turn-on of the transient AXP XTE J1810-197. Simultaneous with the 4.6 hour XMM-Newton observation, we observed at 6.4 GHz with the Parkes telescope, measuring the phase relationship of the radio and X-ray pulse. The X-ray pulsed fraction of 1E 1547.0-5408 is only ~7%, while its radio pulse is relatively broad for such a slow pulsar, which may indicate a nearly aligned rotator. As also inferred from the transient behavior of XTE J1810-197, the only other AXP known to emit in the radio, the magnetic field rearrangement responsible for this X-ray outburst of 1E 1547.0-5408 is probably the cause of its radio turn-on.

After our tentative detection of an optical counterpart to CXOU J010043.1-721134 from archival Hubble Space Telescope (HST) imaging, we have followed up with further images in four bands. Unfortunately, the source originally identified is not confirmed. We provide deep photometric limits in four bands and accurate photometry of field stars around the location of the magnetar.

We present long-term optical and RXTE data of two X-ray binary pulsars in the Small Magellanic Cloud, SXP46.6 and SXP6.85. The optical light curves of both sources show substantial (~0.5-0.8 mag) changes over the time span of the observations. While the optical data for SXP6.85 do not reveal any periodic behaviour, by detrending the optical measurements for SXP46.6 we find an orbital period of ~137 days, consistent with results from the X-ray data. The detection of Type I X-ray outbursts from SXP46.6, combined with the fact that we also see optical outbursts at these times, implies that SXP46.6 is a high orbital eccentricity system. Using contemporaneous optical spectra of SXP46.6 we find that the equivalent width of the H_alpha emission line changes over time indicating that the size of the circumstellar disc varies. By studying the history of the colour variations for SXP6.85 we find that the source gets redder as it brightens which can also be attributed to changes in the circumstellar disc. We do not find any correlation between the X-ray and optical data for SXP6.85. The results for SXP6.85 suggest that it is a low eccentricity binary and that the optical modulations are due to the Be phenomenon.

We report the discovery of burst oscillations at 414.7 Hz during a thermonuclear X-ray burst from the low mass X-ray binary (LMXB) 4U 0614+091 with the Burst Alert Telescope (BAT) onboard SWIFT. In a search of the BAT archive, we found two burst triggers consistent with the position of 4U 0614+091. We searched both bursts for high frequency timing signatures, and found a significant detection at 414.7 Hz during a 5 s interval in the cooling tail of the brighter burst. This result establishes the spin frequency of the neutron star in 4U 0614+091 as 415 Hz. The oscillation had an average amplitude (rms) of 14%, These results are consistent with those known for burst oscillations seen in other LMXBs. The inferred ratio of the frequency difference between the twin kHz QPOs, and the spin frequency in this source is strongly inconsistent with either 0.5 or 1, and tends to support the recent suggestions by Yin et al., and Mendez & Belloni, that the kHz QPO frequency difference may not have a strong connection to the neutron star spin frequency.

Here we consider the nonlinear evolution of Alfven waves that have been excited by gravitational waves from merging binary pulsars. We derive a wave equation for strongly nonlinear and dispersive Alfven waves. Due to the weak dispersion of the Alfven waves, significant wave steepening can occur, which in turn implies strong harmonic generation. We find that the harmonic generation is saturated due to dispersive effects, and use this to estimate the resulting spectrum. Finally we discuss the possibility of observing the above process.

The supernova remnant (SNR) Kes 75/PSR J1846-0258 association can be regarded as certain due to the accurate location of young PSR J1846-0258 at the center of Kes 75 and the detected bright radio/X-ray synchrotron nebula surrounding the pulsar. We provide a new distance estimate to the SNR/pulsar system by analyzing the HI and 13CO maps, the HI emission and absorption spectra, and the 13CO emission spectrum of Kes 75. No absorption features at negative velocities strongly argue against the widely-used large distance of 19 to 21 kpc for Kes 75, and show that Kes 75 is within the Solar circle, i.e. a distance d< 13.2 kpc. Kes 75 is likely at distance of 5.1 to 7.5 kpc because the highest HI absorption velocity is at 95 km/s and no absorption is associated with a nearby HI emission peak at 102 km/s in the direction of Kes 75. This distance to Kes 75 is consistent with the DM distance of pulsar PSR J1846-0258, gives a reasonable luminosity of PSR J1846-0258 and its PWN, and also leads to a much smaller radius for Kes 75. So the age of the SNR is consistent with the spin-down age of PSR J1846-0258, confirming this pulsar as the second-youngest in the Galaxy.

The majority of pulsar population synthesis studies performed to date have focused on isolated pulsar evolution. Those that have incorporated pulsar evolution within binary systems have tended to either treat binary evolution poorly of evolve the pulsar population in an ad-hoc manner. Here we present the first model of the Galactic field pulsar population that includes a comprehensive treatment of both binary and pulsar evolution. Synthetic observational surveys mimicking a variety of radio telescopes are then performed on this population. As such, a complete and direct comparison of model data with observations of the pulsar population within the Galactic disk is now possible. The tool used for completing this work is a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Here we give a brief overview of the method and assumptions involved with each component. Some preliminary results are also presented as well as plans for future applications of the code.

We present a review of observed parameters of giant radio pulses, based on the observations conducted by our group during recent years. The observations cover a broad frequency range of about 3 octaves, concentrating between 600 and 4850 MHz. Giant pulses of both the Crab pulsar and the millisecond pulsar B1937+21 were studied with the 70-m Tidbinbilla, the 100-m GBT, 64-m Kalyazin and Westerbork radio telescopes. We discuss pulse energy distribution, dependence of peak flux density from the pulse width, peculiarities of radio spectra, and polarization properties of giant radio pulses.

We report on two Chandra observations of the 3-Myr pulsar B1929+10, which reveal a faint compact (~9"x4") nebula elongated in the direction perpendicular to the pulsar's proper motion, two patchy wings, and a possible short (~3") jet emerging from the pulsar. In addition, we detect a tail extending up to at least 4' in the direction opposite to the pulsar's proper motion, aligned with the 15'-long tail detected in ROSAT and XMM-Newton observations. The overall morphology of the nebula suggests that the shocked pulsar wind is confined by the ram pressure due to the pulsar's supersonic speed. The shape of the compact nebula in the immediate vicinity of the pulsar seems to be consistent with the current MHD models. However, since these models do not account yet for the change of the flow velocity at larger distances from the pulsar, they are not able to constrain the extent of the long pulsar tail. The luminosity of the whole nebula as seen by Chandra is ~10^30 ergs/s in the 0.3-8 keV band, for the distance of 361 pc. Using the Chandra and XMM-Newton data, we found that the pulsar spectrum is comprised of non-thermal (magnetospheric) and thermal components. The non-thermal component can be described by a power-law model with photon index ~1.7 and luminosity 1.7x10^30 ergs/s in the 0.3-10 keV band. The blackbody fit for the thermal component, which presumably emerges from hot polar caps, gives the temperature kT~0.3 keV and projected emitting area 3x10^3 m^2, corresponding to the bolometric luminosity ~(1-2)x10^30 ergs/s.

Suzaku TOO observation of the anomalous X-ray pulsar CXOU J164710.2-455216 was performed on 2006 September 23--24 for a net exposure of 38.8 ks. During the observation, the XIS was operated in 1/8 window option to achieve a time resolution of 1 s. Pulsations are clearly detected in the XIS light curves with a barycenter corrected pulse period of 10.61063(2) s. The XIS pulse profile shows 3 peaks of different amplitudes with RMS fractional amplitude of ~11% in 0.2--6.0 keV energy band. Though the source was observed with the HXD of Suzaku, the data is highly contaminated by the nearby bright X-ray source GX 340+0 which was in the HXD field of view. The 1-10 keV XIS spectra are well fitted by two blackbody components. The temperatures of two blackbody components are found to be 0.61+/-0.01 keV and 1.22+/-0.06 keV and the value of the absorption column density is 1.73+/-0.03 x 10^{22} atoms cm^{-2}. The observed source flux in 1-10 keV energy range is calculated to be 2.6 x 10^{-11} ergs cm^{-2} s^{-1} with significant contribution from the soft blackbody component (kT = 0.61 keV). Pulse phase resolved spectroscopy of XIS data shows that the flux of the soft blackbody component consists of three narrow peaks, whereas the flux of the other component shows a single peak over the pulse period of the AXP. The blackbody radii changes between 2.2-2.7 km and 0.28-0.38 km (assuming the source distance to be 5 kpc) over pulse phases for the soft and hard components, respectively. The details of the results obtained from the timing and spectral analysis is presented.

We report the analysis of the first superburst from a transiently accreting neutron star system with the All-Sky Monitor (ASM) on the Rossi X-ray Timing Explorer. The superburst occurred 55 days after the onset of an accretion outburst in 4U 1608-522. During that time interval, the accretion rate was at least 7% of the Eddington limit. The peak flux of the superburst is 22 to 45% of the Eddington limit, and its radiation energy output is between 4e41 and 9e41 erg for a distance of 3.2 kpc. Fits of cooling models to the superburst light curve indicate an ignition column depth between 1.5e12 and 4.1e12 g/cm2. Extrapolating the accretion history observed by the ASM, we derive that this column was accreted over a period of 26 to 72 years. The superburst characteristics are consistent with those seen in other superbursting low-mass X-ray binaries. However, the transient nature of the hosting binary presents significant challenges for superburst theory, requiring additional ingredients for the models. The carbon that fuels the superburst is thought to be produced mostly during the accretion outbursts and destroyed in the frequent type-I X-ray bursts. Mixing and sedimentation of the elements in the neutron star envelope may significantly influence the balance between the creation and destruction of carbon. Furthermore, predictions for the temperature of the neutron star crust fail to reach the values required for the ignition of carbon at the inferred column depth.

Some pulsars have their maximum observable energy output in the gamma-ray band, offering the possibility of using these high-energy photons as probes of the particle acceleration and interaction processes in pulsar magnetospheres. After an extended hiatus between satellite missions, the recently-launched AGILE mission and the upcoming Gamma-ray Large Area Space Telescope (GLAST) Large Area Telescope (LAT) will allow gamma-ray tests of the theoretical models developed based on past discoveries. With its greatly improved sensitivity, better angular resolution, and larger energy reach than older instruments, GLAST LAT should detect dozens to hundreds of new gamma-ray pulsars and measure luminosities, light curves, and phase-resolved spectra with unprecedented resolution. It will also have the potential to find radio-quiet pulsars like Geminga, using blind search techniques. Cooperation with radio and X-ray pulsar astronomers is an important aspect of the LAT team's planning for pulsar studies.

X-ray timing of neutron stars in low-mass X-ray binaries (LMXBs) with RXTE has since 1996 revealed several distinct high-frequency phenomena. Among these are oscillations during thermonuclear (type-I) bursts, which (in addition to persistent X-ray pulsations) are thought to trace the neutron star spin. Recent discoveries bring the total number of measured LMXB spin rates to 22. An open question is why the majority of the ~100 known neutron stars in LMXBs show neither pulsations nor burst oscillations.
Recent observations suggest that persistent pulsations may be more common than previously thought, although detectable intermittently, and in some cases at very low duty cycles. For example, the 377.3 Hz pulsations in HETE J1900.1-2455 were only present in the first few months of it's outburst, and have been absent since (although X-ray activity continues). Intermittent (persistent) pulsations have since been detected in a further two sources. In two of these three systems the pulsations appear to be related to the thermonuclear burst activity, but in the third (Aql X-1) they are not. This phenomenon offers new opportunities for spin measurements in known systems.
Such measurements can constrain the poorly-known neutron star equation of state, and neutron stars in LMXBs offer observational advantages over rotation-powered pulsars which make the detection of more rapidly-spinning examples more likely. Even so, spin rates of at least 50% faster than the present maximum appear necessary to give constraints stringent enough to discriminate between the various models. Although the future prospects for such rapidly-spinning objects do not appear optimistic, several additional observational approaches are possible for LMXBs.

Motivated by the string gas cosmological model, which predicts a blue tilt of the primordial gravitational wave spectrum, we examine the constraints imposed by current and planned observations on a blue tilted tensor spectrum. Starting from an expression for the primordial gravitational wave spectrum normalized using cosmic microwave background observations, pulsar timing, direct detection and nucleosynthesis bounds are examined. If we assume a tensor to scalar ratio on scales of the CMB which equals the current observational upper bound, we obtain from these current observations constraints on the tensor spectral index of $n_{T} \lesssim 0.79$, $n_{T} \lesssim 0.53$, and $n_{T} \lesssim 0.15$ respectively.

We report on a comprehensive analysis of the kilo-Hz (>~600 Hz) quasi-periodic oscillations (kHz QPOs) detected from the neutron star X-ray transient Aquila X-1 (Aql X-1) with the Rossi X-ray Timing Explorer, between 1997 and 2007. Among kHz QPO sources, Aql X-1 is peculiar because so far only one kHz QPO has been reported, whereas in most sources, two kHz QPOs are usually detected (a lower and an upper kHz QPO). The identification of the QPOs reported so far has therefore been ambiguous, although it has been proposed that they were likely to be the lower QPO. Following up on previous work, we confirm the identification of the QPOs previously reported as lower QPOs, because of their high quality factors and the quality factor versus frequency dependency, which are similar to those observed in other sources. Combining all segments of data containing a lower QPO, we detect for the first time an upper kHz QPO. As in other sources for which the neutron star spin frequency is larger than 400 Hz (550.25 Hz in Aql X-1), the frequency difference between the two kHz QPOs is close to half the spin frequency. Based on this result, we re-examine the link between the neutron star spin and the frequency of the kHz QPOs, to show that a model in which the separation of the lower and upper QPOs relates to the neutron star spin frequency is still as good as any comparably simple model.

Hot spots residing on the surface of an accretion disc have been considered as a model of short-term variability of active galactic nuclei. In this paper we apply the theory of random point processes to model the observed signal from an ensemble of randomly generated spots. The influence of general relativistic effects near a black hole is taken into account and it is shown that typical features of power spectral density can be reproduced. Connection among spots is also discussed in terms of Hawkes' process, which produces more power at low frequencies. We derive a semi-analytical way to approximate the resulting power-spectral density.

Millihertz quasi-periodic oscillations reported in three neutron-star low mass X-ray binaries have been suggested to be a mode of marginally stable nuclear burning on the neutron star surface. In this Letter, we show that close to the transition between the island and the banana state, 4U~1636--53 shows mHz QPOs whose frequency systematically decreases with time until the oscillations disappear and a Type I X-ray burst occurs. There is a strong correlation between the QPO frequency $\nu$ and the occurrence of X-ray bursts: when $\nu\gtrsim9$ mHz no bursts occur, while $\nu\lesssim9$ mHz does allow the occurrence of bursts. The mHz QPO frequency constitutes the first identified observable that can be used to predict the occurrence of X-ray bursts. If a systematic frequency drift occurs, then a burst happens within a few kilo-seconds after $\nu$ drops below 9 mHz. This observational result confirms that the mHz QPO phenomenon is intimately related with the processes that lead to a thermonuclear burst.

The double pulsar system PSR J0737-3039A/B is a highly relativistic double neutron star (DNS) binary, with a 2.4-hour orbital period. The low mass of the second-formed NS, as well the low system eccentricity and proper motion, point to a different evolutionary scenario compared to other known DNS systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to have a low misalignment between its spin and orbital angular momentum axes, with a 68.3% upper limit of 6.1 degrees, assuming emission from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

Anomalous X-ray Pulsars (AXPs) belong to a class of neutron stars believed to harbor the strongest magnetic fields in the universe, as indicated by their energetic bursts and their rapid spindowns. We have developed a theoretical model that takes into account processes in the atmospheres and magnetospheres of ultramagnetic neutron stars, as well as the effects of their strong gravitational fields on the observable properties. Using this model, we have analyzed the X-ray spectra of a number of AXPs. We find that in all cases, the X-ray spectra are described very well with this emission model. The spectroscopically measured magnetic field strengths of these sources are in close agreement with the values inferred from their spindown properties and provide independent evidence for their magnetar nature. The analysis of spectral data using this physical model also sheds light on the long-term evolution of AXPs.

We present here the results of Arecibo timing of PSR B1516+02B, a 7.95-ms pulsar in a binary system with a ~0.17 solar mass companion and an orbital period of 6.85 days located in the globular cluster M5. The eccentricity of the orbit (e = 0.14) has allowed a measurement of the rate of advance of periastron: (0.0136 +/- 0.0007) degrees per year. It is very likely that the periastron advance is due to the effects of general relativity; the total mass of the binary system is (2.14 +/-0.16) solar masses. The small measured mass function implies, in a statistical sense, that a very large fraction of this total mass is contained in the pulsar: (1.94+0.17 -0.19) solar masses (1-sigma); there is a 5% probability that the mass of this object is below 1.59 solar masses. With the possible exception of PSR J1748-2021B, this is the largest neutron star mass measured to date. When combined with similar measurements made previously for Terzan 5 I and J, we can exclude, in a statistical sense, the ``soft'' equations of state for dense neutron matter, implying that matter at the center of a neutron star is highly incompressible. There is also some evidence for a bimodal distribution of MSP masses, the reasons for that are not clear.

We use a model with two sterile neutrinos obtained by fits to the MiniBoone and LSND experiments. Using formulations with neutrinos created by URCA Processes in a strong magnetic field, so the lowest Landau level has a sizable probability, we find that with known paramenters the assymetric sterile neutrino emissivity might account for large pulsar kicks.

The unprecedented spatial resolution of the Chandra observatory opens the possibility to detect with relatively high accuracy proper motions at X-ray wavelengths. We have conducted an astrometric study of three of the "Magnificent Seven", the thermally emitting radio quiet isolated neutron stars (INSs) discovered by ROSAT. These three INSs (RX J0420.0-5022, RX J0806.4-4123 and RX J1308.6+2127/RBS 1223) either lack an optical counterpart or have one too faint to be used for astrometric purposes. We obtained ACIS observations 3 to 5 years apart to constrain or measure the displacement of the sources on the X-ray sky using as reference the background of extragalactic or remote galactic X-ray sources. Upper limits of 138 mas/yr and 76 mas/yr on the proper motion of RX J0420.0-5022 and RX J0806.4-4123, respectively, have already been presented in Motch et al. (2007). Here we report the very significant measurement (~ 10 sigma) of the proper motion of the third INS of our program, RX J1308.6+2127/RBS1223. Comparing observations obtained in 2002 and 2007 reveals a displacement of 1.1 arcsec implying a yearly proper motion of 223 mas, the second fastest measured for the ROSAT discovered INSs. The source is rapidly moving away from the galactic plane at a speed which precludes any significant accretion of matter from the interstellar medium. Its transverse velocity of ~ 740 (d/700pc) km/s might be the largest of the "Magnificent Seven" and among the fastest recorded for neutron stars. RX J1308.6+2127/RBS1223 is thus a young high velocity cooling neutron star. The source may have its origin in the closest part of the Scutum OB2 association about 0.8 Myr ago, an age consistent with that expected from cooling curves, but significantly younger than inferred from pulse timing measurements (1.5 Myr).

Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically symmetric or axisymmetric. In Newtonian or fixed-background relativistic approximations it is common practice to use spherical polar coordinates for computational grids; however, these coordinates have singularities and are difficult to use in fully relativistic models. We present, in this series of papers, a numerical technique which is able to use effectively spherical grids by employing multiple patches. We provide detailed instructions on how to implement such a scheme, and present a number of code tests for the fixed background case, including an accretion torus around a black hole.

Eddington-limited X-ray bursts from neutron stars can be used in conjunction with other spectroscopic observations to measure neutron star masses, radii, and distances. In order to quantify some of the uncertainties in the determination of the Eddington limit, we analysed a large sample of photospheric radius-expansion thermonuclear bursts observed with the Rossi X-ray Timing Explorer. We identified the instant at which the expanded photosphere "touches down" back onto the surface of the neutron star and compared the corresponding touchdown flux to the peak flux of each burst. We found that for the majority of sources, the ratio of these fluxes is smaller than 1.6, which is the maximum value expected from the changing gravitational redshift during the radius expansion episodes (for a 2M_sun neutron star). The only sources for which this ratio is larger than 1.6 are high inclination sources that include dippers and Cyg X-2. We discuss two possible geometric interpretations of this effect and show that the inferred masses and radii of neutron stars are not affected by this bias. On the other hand, systematic uncertainties as large as ~50% may be introduced to the distance determination.

It is usually proposed that hyperaccretion disks surrounding stellar-mass black holes at an accretion rate of a fraction of one solar mass per second, which are produced during the mergers of double compact stars or the collapses of massive stars, are central engines of gamma-ray bursts (GRBs). In some origin/afterglow models, however, newborn compact objects are invoked to be neutron stars rather than black holes. Thus, hyperaccretion disks around neutron stars seem to exist in some GRBs. Such disks may also occur in type-II supernovae. In this paper we study the structure of a hyperaccretion disk around a neutron star. We consider a steady-state disk model and divide the disk into two regions, called inner and outer disks. The inner disk satisfies an adiabatic self-similar structure and the outer disk is similar to the outer region of a hyperaccretion disk around a black hole. By using analytical and numerical methods, we explore the size of the inner disk, the radial distributions of the density, temperature and pressure of the whole disk, the mechanisms of energy heating and cooling, and the efficiency of neutrino cooling. We find that, compared with a black-hole disk, the hyperaccretion disk around a neutron star can cool more efficiently and produce a much higher neutrino luminosity.

(shortened) We study the Gamma-Ray Burst (GRB) 060218: a particularly close source at z=0.033 with an extremely long duration, namely T_{90} ~ 2000 s, related to SN 2006aj. [...] I present the fitting time consuming procedure. In order to show its sensitivity I also present two examples of fits with the same value of B and different value of E_{e^\pm}^{tot}. We fit the X- and \gamma-ray observations by Swift of GRB 060218 in the 0.1-150 keV energy band during the entire time of observations from 0 all the way to 10^6 s within a unified theoretical model. The free parameters of our theory are only three, namely the total energy E_{e\pm}^{tot} of the e^\pm plasma, its baryon loading B \equiv M_Bc^2/E_{e\pm}^{tot}, as well as the CircumBurst Medium (CBM) distribution. We justify the extremely long duration of this GRB by a total energy E_{e\pm}^{tot} = 2.32\times 10^{50} erg, a very high value of the baryon loading B=1.0\times 10^{-2} and the effective CircumBurst Medium (CBM) density which shows a radial dependence n_{cbm} \propto r^{-\alpha} with 1.0 \leq \alpha \leq 1.7 and monotonically decreases from 1 to $10^{-6}$ particles/cm$^3$. We recall that this value of the $B$ parameter is the highest among the sources we have analyzed and it is very close to its absolute upper limit expected. [...] We also think that the smallest possible black hole, formed by the gravitational collapse of a neutron star in a binary system, is consistent with the especially low energetics of the class of GRBs associated with SNe Ib/c.

Galactic globular clusters harbour binary systems that are detected as faint X-ray sources. These close binaries are thought to play an important role in the stability of the clusters by liberating energy and delaying the inevitable core collapse of globular clusters. The inventory of close binaries and their identification is therefore essential. We present XMM-Newton observations of two Galactic globular clusters: NGC 2808 and NGC 4372. We use X-ray spectral and variability analysis combined with ultra-violet observations made with the XMM-Newton optical monitor and published data from the Hubble Space Telescope to identify sources associated with the clusters. We compare the results of our observations with estimates from population synthesis models. Five sources out of 96 are likely to be related to NGC 2808. Nine sources are found in the field of view of NGC 4372, none being located inside its half-mass radius. We find one quiescent neutron star low mass X-ray binary candidate in the core of NGC 2808, and propose that the majority of the central sources in NGC 2808 are cataclysmic variables. An estimation leads to ~20+/-10 cataclysmic variables with luminosity above 4.25 x 10^31 erg s^-1. Millisecond pulsars could also be present in the core of NGC 2808, and some sources outside of the half-mass radius could possibly be linked to the cluster.

Aims: We performed deep optical observations of the field of an old, fast-moving radio pulsar PSR B1133+16 in an attempt to detect its optical counterpart and a bow shock nebula.
Methods: The observations were carried out using the direct imaging mode of
FORS1 at the ESO VLT/UT1 telescope in the B, R, and H_alpha bands. We also used archival images of the same field obtained with the VLT in the B band and with the Chandra/ACIS in X-rays.
Results: In the B band we detected a faint (B=28.1+/-0.3) source that may be the optical counterpart of PSR B1133+16, as it is positionally consistent with the radio pulsar and with the X-ray counterpart candidate published earlier. Its upper limit in the R band implies a color index B-R <0.5, which is compatible with the index values for most pulsars identified in the optical range. The derived optical luminosity and its ratio to the X-ray luminosity of the candidate are consistent with expected values derived from a sample of pulsars detected in both spectral domains. No Balmer bow shock was detected, implying a low density of ambient matter around the pulsar. However, in the X-ray and H_alpha images we found the signature of a trail extending ~4"-5" behind the pulsar and coinciding with the direction of its proper motion. If confirmed by deeper studies, this is the first time such a trail has been seen in the optical and X-ray wavelengths.
Conclusions: Further observations at later epochs are necessary to confirm the identification of the pulsar by the candidate's proper motion measurements.

Context: Since its launch in 2002, \integral has discovered many new hard X-ray sources. A lot of them still lack sufficient positional accuracy, preventing counterparts at other wavelengths from being found. Their true nature is, therefore, still unknown. Aims: The goal of this study is to give an accurate X-ray position for 12 of these sources with the view to further identify their counterpart at optical, infrared and radio wavelength, and unveil their true nature. We also make use of the X-ray spectral parameters to tentatively discriminate between the various possible types. Methods: We made use of public X-ray observations with the X-ray telescope on-board the Swift observatory to refine the X-ray position to 3-5\arcsec accuracy, and performed 0.1--10 keV spectral analysis. We then searched the online catalogues (e.g. NED, SIMBAD, 2MASS, 2MASX and NVSS) to search for counterparts at other wavelengths. Results:
For all sources, we give a refined X-ray position, provide X-ray spectral parameters and identify infrared counterparts. We confirm the nature of five sources formerly suspected to be AGN (IGR J02343+3229, J13149+4422, J14579$-$4308, J16385$-$2057, J19378$-$0617) while we possibly reject the AGN nature for one (IGR J18559+1535), and suggest that it is a Galactic source instead. All other sources may be Galactic sources, and in that case their spectral shape may suggest that they are X-ray Binaries. In one case (IGR J19308+0530), the Galactic nature is confirmed through the identification of an F8 star as the counterpart. We favour a distance to the source not higher than 1 kpc. The source is likely to be a neutron star XRB or a CV. We also report the discovery of six serendipitous sources of unknown nature.

Mergers of double neutron stars are considered the most likely progenitors for short gamma-ray bursts. Indeed such a merger can produce a black hole with a transient accreting torus of nuclear matter and the conversion of the torus mass-energy to radiation can power a gamma-ray burst. Using available binary pulsar observations supported by our extensive evolutionary calculations of double neutron star formation, we demonstrate that the fraction of mergers that can form a black hole -- torus system depends very sensitively on the (largely unknown) maximum neutron star mass. We show that the available observations and models put a very stringent constraint on this maximum mass under the assumption that a majority of short gamma-ray bursts originate in double neutron star mergers. Specifically, we find that the maximum neutron star mass must be within 2--2.5 Msun. Moreover, a single unambiguous measurement of a neutron star mass above 2.5 Msun would exclude double neutron star mergers as short gamma-ray burst progenitors.

Traditionally, studies aimed at inferring the distribution of birth periods of neutron stars are based on radio surveys. Here we propose an independent method to constrain the pulsar spin periods at birth based on their X-ray luminosities. In particular, the observed luminosity distribution of supernovae poses a constraint on the initial rotational energy of the embedded pulsars, via the L_X-dot{E}_{rot} correlation found for radio pulsars, and under the assumption that this relation continues to hold beyond the observed range. We have extracted X-ray luminosities (or limits) for a large sample of historical SNe observed with Chandra, XMM and Swift, that have been firmly classified as core-collapse supernovae. We have then compared these observational limits with the results of Monte Carlo simulations of the pulsar X-ray luminosity distribution, for a range of values of the birth parameters. We find that a pulsar population dominated by millisecond periods at birth is ruled out by the data.

We study light curves and spectra (equivalent widths of the iron line and some other spectral characteristics) which arise by reprocessing on the surface of an accretion disc, following its illumination by a primary off-axis source - an X-ray 'flare', assumed to be a point-like source just above the accretion disc. We consider all general relativity effects (energy shifts, light bending, time delays, delay amplification due to the spot motion) near a rotating black hole. For some sets of parameters the reflected flux exceeds the flux from the primary component. We show that the orbit-induced variations of the equivalent width with respect to its mean value can be as high as 30% for an observer's inclination of 30 degrees, and much more at higher inclinations. We calculate the ratio of the reflected flux to the primary flux and the hardness ratio which we find to vary significantly with the spot phase mainly for small orbital radii. This offers the chance to estimate the lower limit of the black hole spin if the flare arises close to the black hole. We show the results for different values of the flare orbital radius.

Evidence is growing for a class of gamma-ray bursts (GRBs) characterized by an initial ~0.1-1 s spike of hard radiation followed, after a ~10 s lull in emission, by a softer period of extended emission lasting ~10-100 s. In a few well-studied cases, these ``short GRBs with extended emission'' show no evidence for a bright associated supernova (SN). We propose that these events are produced by the formation and early evolution of a highly magnetized, rapidly rotating neutron star (a ``proto-magnetar'') which is formed from the accretion-induced collapse (AIC) of a white dwarf (WD), the merger and collapse of a WD-WD binary, or, perhaps, the merger of a double neutron star binary. The initial emission spike is powered by accretion onto the proto-magnetar from a small disk that is formed during the AIC or merger event. The extended emission is produced by a relativistic wind that extracts the rotational energy of the proto-magnetar on a timescale ~10-100 s. The ~10 s delay between the prompt and extended emission is the time required for the newly-formed proto-magnetar to cool sufficiently that the neutrino-heated wind from its surface becomes ultra-relativistic. Because a proto-magnetar ejects little or no Ni56 (< 1e-3 M_sun), these events should not produce a bright SN-like transient. We model the extended emission from GRB060614 using spin-down calculations of a cooling proto-magnetar, finding reasonable agreement with observations for a magnetar with an initial rotation period of ~1 ms and a surface dipole field of ~3e15 G. If GRBs are indeed produced by AIC or WD-WD mergers, they should occur within a mixture of both early and late-type galaxies and should not produce strong gravitational wave emission. An additional consequence of our model is the existence of X-ray flashes unaccompanied by a bright SN.

We analyze 1050 ks of INTEGRAL data of the high mass X-ray binary pulsar 1E 1145.1-6141 to study its properties over a long time baseline, from June 2003 to June 2004, with wide spectral coverage.
We study three high luminosity episodes, two of them at the system apoastron, three brightening with lower intensity, two at the periastron, and one extended period of intermediate luminosity spanning one orbital cycle. We perform timing analysis to determine the pulse period and pulse profiles at different energy ranges. We also analyze the broad band phase average spectrum of different luminosity states and perform phase resolved spectroscopy for the first flare.
From the timing analysis, we find a pulse period of ~297 s around MJD 53000 with a significant scatter around the mean value. From the spectral analysis we find that the source emission can be described by an absorbed bremsstrahlung model in which the electron temperature varies between ~25 and ~37 keV, without any correlation to luminosity, and the intrinsic absorbing column is constantly of the order of 10^23 cm^-2. Phase resolved spectral analysis evidences a different temperature of the plasma in the ascending and descending edges of the pulse during the first flare. This justifies the pulse maximum shift by ~0.4 phase units between 20 and 100 keV observed in the pulse profiles.
The comparison with the previous period measurements reveals that the source is currently spinning-down, in contrast to the long term secular trend observed so far indicating that at least a temporary accretion disk is formed. The study of the spectral property variations with respect to time and spin phase suggests the presence of two emitting components at different temperatures whose relative intensity varies with time.

New GMRT observations of the five-component pulsar B1857--26 provide detailed insight into its pulse-sequence modulation phenomena for the first time. The outer conal components exhibit a 7.4-rotation-period, longitude-stationary modulation. Several lines of evidence indicate a carousel circulation time $\P3hat$ of about 147 stellar rotations, characteristic of a pattern with 20 beamlets. The pulsar nulls some 20% of the time, usually for only a single pulse, and these nulls show no discernible order or periodicity. Finally, the pulsar's polarization-angle traverse raises interesting issues: if most of its emission is comprised of a single polarization mode, the full traverse exceeds 180\degr; or if both polarization modes are present, then the leading and the trailing halves of the profiles exhibit two different modes. In either case the rotating vector model fails to fit the polarization-angle traverse of the core component.

The impact of strong magnetic fields B>10e13 G on the thermal evolution of neutron stars is investigated, including crustal heating by magnetic field decay. For this purpose, we perform 2D cooling simulations with anisotropic thermal conductivity considering all relevant neutrino emission processes for realistic neutron stars. The standard cooling models of neutron stars are called into question by showing that the magnetic field has relevant (and in many cases dominant) effects on the thermal evolution. The presence of the magnetic field significantly affects the thermal surface distribution and the cooling history of these objects during both, the early neutrino cooling era and the late photon cooling era. The minimal cooling scenario is thus more complex than generally assumed. A consistent magneto-thermal evolution of magnetized neutron stars is needed to explain the observations.

We report about very high energy (VHE) gamma-ray observations of the Crab Nebula with the MAGIC telescope. The gamma-ray flux from the nebula was measured between 60 GeV and 9 TeV. The energy spectrum can be described with a curved power law dF/dE=f_0 (E/300GeV)^(a+b log10(E/300GeV)) with a flux normalization f_0 of (6.0+-0.2stat)*10^-10 cm^-2 s^-1 TeV^-1, a=-2.31+-0.06stat and b=-0.26+-0.07stat. The position of the IC-peak is determined at 77+-47 GeV. Within the observation time and the experimental resolution of the telescope, the gamma-ray emission is steady and pointlike. The emission's center of gravity coincides with the position of the pulsar. Pulsed gamma-ray emission from the pulsar could not be detected. We constrain the cutoff energy of the spectrum to be less than ~30 GeV, assuming that the differential energy spectrum has an exponential cutoff. For a super-exponential shape, the cutoff energy can be as high as ~60GeV.

We present results of a timing analysis of various isolated pulsars using ESA's \emph{XMM-Newton} observatory. Isolated pulsars are useful for calibration purposes because of their stable emission. We have analyzed six pulsars with different pulse profiles in a range of periods between 15 and 200 ms. All observations were made using the \emph{EPIC-pn camera} in its faster modes (Small window, Timing and Burst modes). We investigate the relative timing accuracy of the camera by comparing the pulse periods determined from the \emph{EPIC-pn camera} observations with those from radio observations. As a result of our analysis we conclude that the relative timing accuracy of the \emph{EPIC-pn camera} is of the order of $1\times 10^{-8}$.

LS I +61 303 is a puzzling object detected from radio up to high-energy gamma-rays. Variability has recently been observed in its high-energy emission. The object is a binary system, with a compact object and a Be star as primary. The nature of the secondary and the origin of the gamma-ray emission are not clearly established at present. Recent VLBA radio data have been used to claim that the system is a Be/neutron star colliding wind binary, instead of a microquasar. We review the main views on the nature of LS I +61 303 and present results of 3D SPH simulations that can shed some light on the nature of the system. Our results support an accretion powered source, compatible with a microquasar interpretation.

We describe a model within the ``Quark-nova'' scenario to interpret the recent observations of early X-ray afterglows of long Gamma-Ray Bursts (GRB) with the Swift satellite. This is a three-stage model within the context of a core-collapse supernova. STAGE 1 is an accreting (proto-) neutron star leading to a possible delay between the core collapse and the GRB. STAGE 2 is accretion onto a quark-star, launching an ultrarelativistic jet generating the prompt GRB. This jet also creates the afterglow as the jet interacts with the surrounding medium creating an external shock. Slower shells ejected from the quark star (during accretion), can re-energize the external shock leading to a flatter segment in the X-ray afterglow. STAGE 3, which occurs only if the quark-star collapses to form a black-hole, consists of an accreting black-hole. The jet launched in this accretion process interacts with the preceding quark star jet, and could generate the flaring activity frequently seen in early X-ray afterglows. Alternatively, a STAGE 2b can occur in our model if the quark star does not collapse to a black hole. The quark star in this case can then spin down due to magnetic braking, and the spin down energy may lead to flattening in the X-ray afterglow as well. This model seems to account for both the energies and the timescales of GRBs, in addition to the newly discovered early X-ray afterglow features.

The ultra-compact Low Mass X-ray Binary (LMXB) X1916-053, composed of a neutron star and a semi-degenerated white dwarf, exhibits periodic X-ray dips with variable width and depth. We have developed new methods to parameterize the dip to systematically study its variations. This helps to further understand binary and accretion disk behaviors. The RXTE 1998 observations clearly show a 4.87d periodic variation of the dip width. This is probably due to the nodal precession of the accretion disk, although there are no significant sidebands in the spectrum from the epoch folding search. From the negative superhump model (Larwood et. al. 1996), the mass ratio can be estimated as q = 0.045. Combined with more than 24 years of historical data, we found an orbital period derivative of $\dot{P}_{orb}/P_{orb}=(1.62 \pm 0.48)\times 10^{-7} yr^{-1}$ and established a quadratic ephemeris for the X-ray dips. The period derivative seems inconsistent with the prediction of the standard model of binary orbital evolution proposed by Rappaport et. al. (1987). On the other hand, the radiation-driven model (Tavani et. al. 1991) may properly interpret the period derivative even though the large mass outflow predicted by this model has never been observed in this system. With the best ephemeris, we obtained that the standard deviation of primary dips are smaller than that of secondary dips. This means that the primary dips are more stable than the secondary dips. Thus, we conclude that the primary dips of X1916-053 occur from the bulge at the rim instead of the ring of the disk proposed by Frank et. al. (1987).

I present polarization modelling of Active Galactic Nuclei in the optical/UV range. The modelling is conducted using the Monte-Carlo radiative transfer code Stokes, which self-consistently models the polarization signature of a complex model arrangement for an active nucleus. In this work I include three different scattering regions around the central source: an equatorial electron scattering disk, an equatorial obscuring dusty torus, and polar electron scattering cones. I investigate the resulting dependencies of the V-band polarization for different optical depths of the scattering cones, different dust compositions inside the torus, and various half-opening angles of the torus/polar cones. The observed polarization dichotomy can be successfully reproduced by the model.

INTEGRAL tripled the number of super-giant high-mass X-ray binaries (sgHMXB) known in the Galaxy by revealing absorbed and fast transient (SFXT) systems. Quantitative constraints on the wind clumping of massive stars can be obtained from the study of the hard X-ray variability of SFXT. A large fraction of the hard X-ray emission is emitted in the form of flares with a typical duration of 3 ksec, frequency of 7 days and luminosity of 1E36 ergs/s. Such flares are most probably emitted by the interaction of a compact object orbiting at ~ 10 R* with wind clumps (1E(22-23) g) representing a large fraction of the stellar mass-loss rate. The density ratio between the clumps and the inter-clump medium is 1E(2-4) . The parameters of the clumps and of the inter-clump medium, derived from the SFXT flaring behavior, are in good agreement with macro-clumping scenario and line-driven instability simulations. SFXT are likely to have larger orbital radius than classical sgHMXB.

We perform a series of two-dimensional magnetohydrodynamic core-collapse simulations of rapidly rotating and strongly magnetized massive stars. To study the properties of magnetic explosions for a longer time stretch of postbounce evolution, we develop a new code under the framework of special relativity including a realistic equation of state with a multiflavor neutrino leakage scheme. Our results show the generation of the magnetically-dominated jets in the two ways. One is launched just after the core-bounce in a prompt way and another is launched at $ \sim 100 $ ms after the stall of the prompt shock. We find that the shock-revival occurs when the magnetic pressure becomes strong, due to the field wrapping, enough to overwhelm the ram pressure of the accreting matter. The critical toroidal magnetic fields for the magnetic shock-revival are found to be universal of $\sim 10^{15}\mathrm{G}$ behind the jets. We point out that the time difference before the shock-revival has a strong correlation with the explosions energies. Our results suggest that the magnetically dominated jets are accompanied by the formation of the magnetars. Since the jets are mildly relativistic, we speculate that they might be the origin of some observed X-ray flashes.

IGR J16493-4348 was one of the first new sources to be detected by the INTEGRAL gamma-ray telescope in the 18-100 keV energy band. Based upon spatial coincidence the source was originally associated with the free radio pulsar PSR J1649-4349. Presented here are the results of 2.8 Ms of observations made by the INTEGRAL mission and a 5.6 ks observation with the Swift X-ray Telescope. Spectral analysis indicates that the source is best modeled by an absorbed power law with a high energy cut-off at E$_{cut}$~15 keV and a hydrogen absorbing column of NH=5.4$^{+1.3}_{-1}$ x 10$^{22}$ cm$^{-2}$. Analysis of the light curves indicates that the source is a weak, persistent gamma-ray emitter showing indications of variability in the 2-9 and 22-100 keV bands. The average source flux is ~1.1 x 10^{-10} erg cm$^{-2}$ s$^{-1}$ in the 1-100 keV energy band. No coherent timing signal is identified at any timescale in the INTEGRAL or Swift data.
The refined source location and positional uncertainty of IGR J16493-4348 places PSR J1649-4349 outside of the 90% error circle. We conclude that IGR J16493-4348 is not associated with PSR J1649-4349. Combining the INTEGRAL observations with Swift/XRT data and information gathered by RXTE and Chandra we suggest that IGR J16493-4348 is an X-ray binary; and that the source characteristics favour a high mass X-ray binary although an LMXB nature cannot be ruled out.

We find numerical solutions of the coupled system of Einstein-Maxwell's equations with a linear approach, in which the magnetic field acts as a perturbation of a spherical neutron star. In our study, magnetic fields having both poloidal and toroidal components are considered, and higher order multipoles are also included. We evaluate the deformations induced by different field configurations, paying special attention to those for which the star has a prolate shape. We also explore the dependence of the stellar deformation on the particular choice of the equation of state and on the mass of the star. Our results show that, for neutron stars with mass M = 1.4 Msun and surface magnetic fields of the order of 10^15 G, a quadrupole ellipticity of the order of 10^(-6) - 10^(-5) should be expected. Low mass neutron stars are in principle subject to larger deformations (quadrupole ellipticities up to 10^(-3) in the most extreme case). The effect of quadrupolar magnetic fields is comparable to that of dipolar components. A magnetic field permeating the whole star is normally needed to obtain negative quadrupole ellipticities, while fields confined to the crust typically produce positive quadrupole ellipticities.

Central Compact Objects (CCOs) are a handful of soft X-ray sources located close to the centers of Supernova Remnants and supposed to be young, radio-quiet Isolated Neutron Stars (INSs). A clear understanding of their physics would be crucial in order to complete our view of the birth properties of INSs. We will review the phenomenologies of CCOs, underlining the most important, recent results, and we will discuss the possible relationships of such sources with other classes of INSs.

We have proposed the new explanation of some magnetic chemically peculiar (MCP) stars anomalies, which is based on assumption that such stars can be the close binary systems with a secondary component being neutron star. Within this hypothesis one can naturally explain the main anomalous features of MCP stars: first of all, an existence of the short-lived radioactive isotopes detected in some stars (like Przybylski's star and HR465), and some others peculiarities (e.g. the behavior of CU Vir in radio range, the phenomenon of the roAp stars).

Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss our new fully general relativistic calculations of merging BHNS binaries, which use high-accuracy, low-eccentricity, conformal thin-sandwich configurations as initial data. Our evolutions are performed using the moving puncture method and include a fully relativistic, high-resolution shock-capturing hydrodynamics treatment. Focusing on systems in which the neutron star is irrotational and the black hole is nonspinning with a 3:1 mass ratio, we investigate the inspiral, merger, and disk formation in the system. We find that the vast majority of material is promptly accreted and no more than 3% of the neutron star's rest mass is ejected into a tenuous, gravitationally bound disk. We find similar results for mass ratios of 2:1 and 1:1, even when we reduce the NS compaction in the 2:1 mass ratio case. These ambient disks reach temperatures suitable for triggering SGRBs, but their masses may be too small to produce the required total energy output. We measure gravitational waveforms and compute the effective strain in frequency space, finding measurable differences between our waveforms and those produced by binary black hole mergers within the advanced LIGO band. These differences appear at frequencies corresponding to the emission that occurs when the NS is tidally disrupted and accreted by the black hole. The resulting information about the radius of the neutron star may be used to constrain the neutron star equation of state.

We have derived isotopic fractions of europium, samarium, and neodymium in two metal-poor giants with differing neutron-capture nucleosynthetic histories. These isotopic fractions were measured from new high resolution (R ~ 120,000), high signal-to-noise (S/N ~ 160-1000) spectra obtained with the 2dCoude spectrograph of McDonald Observatory's 2.7m Smith telescope. Synthetic spectra were generated using recent high-precision laboratory measurements of hyperfine and isotopic subcomponents of several transitions of these elements and matched quantitatively to the observed spectra. We interpret our isotopic fractions by the nucleosynthesis predictions of the stellar model, which reproduces s-process nucleosynthesis from the physical conditions expected in low-mass, thermally-pulsing stars on the AGB, and the classical method, which approximates s-process nucleosynthesis by a steady neutron flux impinging upon Fe-peak seed nuclei. Our Eu isotopic fraction in HD 175305 is consistent with an r-process origin by the classical method and is consistent with either an r- or an s-process origin by the stellar model. Our Sm isotopic fraction in HD 175305 suggests a predominantly r-process origin, and our Sm isotopic fraction in HD 196944 is consistent with an s-process origin. The Nd isotopic fractions, while consistent with either r-process or s-process origins, have very little ability to distinguish between any physical values for the isotopic fraction in either star. This study for the first time extends the n-capture origin of multiple rare earths in metal-poor stars from elemental abundances to the isotopic level, strengthening the r-process interpretation for HD 175305 and the s-process interpretation for HD196944.

One of the primary goals when studying stellar systems with neutron stars has been to reveal the physical properties of progenitors and understand how neutron star spins and birth kicks are determined. Over the years a consensus understanding had been developed, but recently some of the basic elements of this understanding are being challenged by current observations of some binary systems and their theoretical interpretation. In what follows we review such recent developments and highlight how they are interconnected; we particularly emphasize some of the assumptions and caveats of theoretical interpretations and examine their validity (e.g., in connection to the unknown radial velocities of pulsars or the nuances of multi-dimensional statistical analysis). The emerging picture does not erase our earlier understanding; instead it broadens it as it reveals additional pathways for neutron star formation and evolution: neutron stars probably form at the end of both core collapse of Fe cores of massive stars and electron-capture supernovae of ONeMg cores of lower-mass stars; birth kicks are required to be high (well in excess of 100 km/s) for some neutron stars and low (< 100 km/s) for others depending on the formation process; and spin up may occur not just through Roche-lobe overflow but also through wind accretion or phases of hypercritical accretion during common envelope evolution.

We present the results of X-ray and near-IR observations of the anomalous X-ray pulsar 1E 1048.1-5937, believed to be a magnetar. This AXP underwent a period of extreme variability during 2001-2004, but subsequently entered an extended and unexpected quiescence in 2004-2006, during which we monitored it with RXTE, CXO, and HST. Its timing properties were stable for >3 years throughout the quiescent period. 1E 1048.1-5937 again went into outburst in March 2007, which saw a factor of >7 total X-ray flux increase which was anti-correlated with a pulsed fraction decrease, and correlated with spectral hardening, among other effects. The near-IR counterpart also brightened following the 2007 event. We discuss our findings in the context of the magnetar and other models.

We consider a situation in which a pulsar is formed inside or close to a high density region of a molecular cloud. Right after birth, the pulsar was very active and accelerated hadrons and leptons to very high energies. Hadrons diffuse through the supernova remnant (SNR) and some of them are trapped in the nearby cloud interacting with the matter. We extend a recent time-dependent model for the gamma-radiation of pulsar wind nebulae (PWNe) to describe this more complicated astrophysical scenario. The example calculations have been performed for two objects, IC443 and W41, which have recently been discovered as sources of TeV gamma-rays. In this model the low energy TeV emission should be correlated with the birth place of the pulsar and the region of dense soft radiation rather than with its present position, provided that the injection rate of relativistic particles into the PWNa has been much more efficient at early times. The high energy TeV emission should be correlated with the location of dense clouds which were able to capture high energy hadrons due to their strong magnetic fields.

We discuss the flavor conversion of neutrinos from core collapse supernovae that have oxygen-neon-magnesium (ONeMg) cores. Using the numerically calculated evolution of the star up to 650 ms post bounce, we find that, for the normal mass hierarchy, the electron neutrino flux in a detector shows signatures of two typical features of an ONeMg-core supernova: a sharp step in the density profile at the base of the He shell and a faster shock wave propagation compared to iron core supernovae. Before the shock hits the density step (t ~ 150 ms), the survival probability of electron neutrinos is about 0.68, in contrast to values of 0.32 or less for an iron core supernova. The passage of the shock through the step and its subsequent propagation cause a decrease of the survival probability and a decrease of the amplitude of oscillations in the Earth, reflecting the transition to a more adiabatic propagation inside the star. These changes affect the lower energy neutrinos first; they are faster and more sizable for larger theta_13. They are unique of ONeMg-core supernovae, and give the possibility to test the speed of the shock wave. The time modulation of the Earth effect and its negative sign at the neutronization peak are the most robust signatures in a detector.

We report here on recent progress in understanding the birth conditions of neutron stars and the way how supernovae explode. More sophisticated numerical models have led to the discovery of new phenomena in the supernova core, for example a generic hydrodynamic instability of the stagnant supernova shock against low-mode nonradial deformation and the excitation of gravity-wave activity in the surface and core of the nascent neutron star. Both can have supportive or decisive influence on the inauguration of the explosion, the former by improving the conditions for energy deposition by neutrino heating in the postshock gas, the latter by supplying the developing blast with a flux of acoustic power that adds to the energy transfer by neutrinos. While recent two-dimensional models suggest that the neutrino-driven mechanism may be viable for stars from about 8 solar masses to at least 15 solar masses, acoustic energy input has been advocated as an alternative if neutrino heating fails. Magnetohydrodynamic effects constitute another way to trigger explosions in connection with the collapse of sufficiently rapidly rotating stellar cores, perhaps linked to the birth of magnetars. The global explosion asymmetries seen in the recent simulations offer an explanation of even the highest measured kick velocities of young neutron stars.

With LOFAR beginning operation in 2008 there is huge potential for studying pulsars with high signal to noise at low frequencies. We present results of observations made with the Westerbork Synthesis Radio Telescope to revisit, with modern technology, this frequency range. Coherently dedispersed profiles of millisecond pulsars obtained simultaneously between 115-175 MHz are presented. We consider the detections and non-detections of 14 MSPs in light of previous observations and the fluxes, dispersion measures and spectral indices of these pulsars. The excellent prospects for LOFAR finding new MSPs and studying the existing systems are then discussed in light of these results.

We investigate the hyperon bulk viscosity due to the non-leptonic process $n + p \rightleftharpoons p + \Lambda $ in $K^-$ condensed matter and its effect on the r-mode instability in neutron stars. We find that the hyperon bulk viscosity coefficient in the condensate phase is significantly suppressed than that in the hadronic phase. The suppressed hyperon bulk viscosity in the superconducting phase is still an efficient mechanism to damp the r-mode instability in neutron stars.

RX J1856.5-3754 is the X-ray brightest among the nearby isolated neutron stars. Its X-ray spectrum is thermal, and is reproduced remarkably well by a black-body, but its interpretation has remained puzzling. One reason is that the source did not exhibit pulsations, and hence a magnetic field strength--vital input to atmosphere models--could not be estimated. Recently, however, very weak pulsations were discovered. Here, we analyze these in detail, using all available data from the XMM-Newton and Chandra X-ray observatories. From frequency measurements, we set a 2-sigma upper limit to the frequency derivative of \dot\nu<1.3e-14 Hz/s. Trying possible phase-connected timing solutions, we find that one solution is far more likely than the others, and we infer a most probable value of \dot\nu=(-5.98+/-0.14)e-16 Hz/s. The inferred magnetic field strength is 1.5e13 G, comparable to what was found for similar neutron stars. From models, the field seems too strong to be consistent with the absence of spectral features for non-condensed atmospheres. It is sufficiently strong, however, that the surface could be condensed, but only if it is consists of heavy elements like iron. Our measurements imply a characteristic age of about 4 Myr. This is longer than the cooling and kinematic ages, as was found for similar objects, but at almost a factor ten, the discrepancy is more extreme. A puzzle raised by our measurement is that the implied rotational energy loss rate of about 3e30 erg/s is orders of magnitude smaller than what was inferred from the H-alpha nebula surrounding the source.

A detailed high-resolution spectroscopic analysis is presented for the carbon-rich low metallicity Galactic halo object CS 22964-161. We have discovered that CS 22964-161 is a double-lined spectroscopic binary, and have derived accurate orbital components for the system. From a model atmosphere analysis we show that both components are near the metal-poor main-sequence turnoff. Both stars are very enriched in carbon and in neutron-capture elements that can be created in the s-process, including lead. The primary star also possesses an abundance of lithium close to the value of the ``Spite-Plateau''. The simplest interpretation is that the binary members seen today were the recipients of these anomalous abundances from a third star that was losing mass as part of its AGB evolution. We compare the observed CS 22964-161 abundance set with nucleosynthesis predictions of AGB stars, and discuss issues of envelope stability in the observed stars under mass transfer conditions, and consider the dynamical stability of the alleged original triple star. Finally, we consider the circumstances that permit survival of lithium, whatever its origin, in the spectrum of this extraordinary system.

We investigate the effect of including a significant ``binary twin'' population (binaries with almost equal mass stars, q = M2/M1 > 0.95) for the production of double compact objects and some resulting consequences, including LIGO inspiral rate and some properties of short-hard gamma-ray bursts. We employ very optimistic assumptions on the twin fraction (50%) among all binaries, and therefore our calculations place an upper limits on the influence of twins on double compact object populations. We show that for LIGO the effect of including twins is relatively minor: although the merger rates does indeed increase when twins are considered, the rate increase is fairly small (1.5). Also, chirp mass distribution for double compact objects formed with or without twins are almost indistinguishable. If double compact object are short-hard GRB progenitors, including twins in population synthesis calculations does not alter significantly the earlier rate predictions for the event rate. However, for one channel of binary evolution, introducing twins more than doubles the rate of ``very prompt'' NS-NS mergers (time to merger less than 1 Myr) compared to models with the ``flat'' q distribution. In that case, 70% of all NS-NS binaries merge within 100 Myr after their formation, indicating a possibility of a very significant population of ``prompt'' short-hard gamma-ray bursts, associated with star forming galaxies. We also point out that, independent of assumptions, fraction of such prompt neutron star mergers is always high, 35--70%. We note that recent observations (e.g., Berger et al.) indicate that fraction of short-hard GRBs found in young hosts is at least 40% and possibly even 80%.

Anomalous X-ray Pulsars (AXPs), thought to be magnetars, exhibit poorly understood deviations from a simple spin-down called "timing noise". AXP timing noise has strong low-frequency components which pose significant challenges for quantification. We describe a procedure for extracting two quantities of interest, the intensity and power spectral index of timing noise. We apply this procedure to timing data from three sources: a monitoring campaign of five AXPs, observations of five young pulsars, and the stable rotator PSR B1937+21.

Nearby masses can have a high probability of lensing stars in a distant background field. High-probability lensing, or mesolensing, can therefore be used to dramatically increase our knowledge of dark and dim objects in the solar neighborhood, where it can discover and study members of the local dark population (free-floating planets, low-mass dwarfs, white dwarfs, neutron stars, and stellar mass black holes). We can measure the mass and transverse velocity of those objects discovered (or already known), and determine whether or not they are in binaries with dim companions. We explore these and other applications of mesolensing, including the study of forms of matter that have been hypothesized but not discovered, such as intermediate-mass black holes, dark matter objects free-streaming through the Galactic disk, and planets in the outermost regions of the solar system. In each case we discuss the feasibility of deriving results based on present-day monitoring systems, and also consider the vistas that will open with the advent of all-sky monitoring in the era of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), and the Large Synoptic Survey Telescope (LSST).

The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures.

Internal flows inside gravitationally stable astrophysical objects, such as the Sun, stars and compact stars are compressed and extremely subsonic. Such low Mach number flows are usually encountered when studying for example dynamo action in stars, planets, the hydro-thermodynamics of X-ray bursts on neutron stars and dwarf novae. Treating such flows is numerically complicated and challenging task. We aim to present a robust numerical tool that enables modeling the time-evolution or quasi-stationary of stratified low Mach number flows under astrophysical conditions. It is argued that astrophysical low Mach number flows cannot be considered as an asymptotic limit of incompressible flows, but rather as highly compressed flows with extremely stiff pressure terms. Unlike the pseudo-pressure in incompressible fluids, a Possion-like treatment for the pressure would smooth unnecessarily the physically induced acoustic perturbations, thereby violating the conservation character of the compressible equations. Moreover, classical dimensional splitting techniques, such as ADI or Line-Gauss-Seidel methods are found to be unsuited for modeling compressible flows with low Mach numbers. In this paper we present a nonlinear Newton-type solver that is based on the defect-correction iteration procedure and in which the Approximate Factorization Method (AFM) is used as a preconditioner. This solver is found to be sufficiently robust and is capable of capturing stationary solutions for viscous rotating flows with Mach number as small as $\mcal{M} \approx 10^{-3},$ i.e., near the incompressibility limit.

We report the results of 18 years of Arecibo timing of two pulsars in the globular cluster NGC 5904 (M5), PSR B1516+02A and PSR B1516+02B. This has allowed the measurement of the proper motions of these pulsars and of the cluster. PSR B1516+02B is a 7.95-ms pulsar in a binary system with a ~0.2 solar-mass companion and an orbital period of 6.86 days. In deep HST images, no optical counterpart is detected at the position of the pulsar, implying the companion is either a white dwarf or a low-mass MS star. The eccentricity of the orbit (e = 0.14) has allowed a measurement of the rate of advance of periastron: 0.0136 +/ 0.0007 degrees per year. It is very likely that the periastron advance is due to the effects of general relativity; the total mass of the binary system is then 2.14 +/- 0.16 solar masses. The small measured mass function implies, in a statistical sense, that a very large fraction of this total mass is contained in the pulsar: 1.94 +0.17/-0.19 solar masses (1 sigma$); there is a 5 % probability that the mass of this object is smaller than 1.59 solar masses and a 1.3% probability that it is between 1.2 and 1.44 solar masses. With the possible exception of PSR J1748-2021B, this is the largest neutron star mass measured to date. When combined with similar measurements made previously for Terzan 5 I and J, we conclude that there is a 99 % probability that at least one of these MSPs is more massive than 1.72 solar masses. Confirmation of these mass measurements would exclude most of the ``soft'' equations of state for dense neutron matter, implying that matter at the center of a neutron star is highly incompressible. Furthermore, we see evidence for a bi-modal MSP mass distribution, but the reasons for this are not clear.

<Context>. We report on near-infrared (IR) observations of the three anomalous X-ray pulsars XTE J1810-197, 1RXS J1708-4009, 1E 1841-045 and the soft gamma-ray repeater SGR 1900+14, taken with the ESO-VLT, the Gemini, and the CFHT telescopes. <Aims>. This work is aimed at identifying and/or confirming the IR counterparts of these magnetars, as well as at measuring their possible IR variability. <Methods>. In order to perform photometry of objects as faint as Ks~20, we have used data taken with the largest telescopes, equipped with the most advanced IR detectors and in most of the cases with Adaptive Optics devices. The latter are critical to achieve the sharp spatial accuracy required to pinpoint faint objects in crowded fields. <Results>. We confirm with high confidence the identification of the IR counterpart to XTE J1810-197, and its IR variability. For 1E 1841-045 and SGR 1900+14 we propose two candidate IR counterparts based on the detection of IR variability. For 1RXS J1708-4009 we show that none of the potential counterparts within the source X-ray error circle can be yet convincingly associated with this AXP. <Conclusions>. The IR variability of the AXP XTE J1810-197 does not follow the same monotonic decrease of its post-outburst X-ray emission. Instead, the IR variability appears more similar to the one observed in radio band, although simultaneous IR and radio observations are crucial to draw any conclusion in this respect. For 1E 1841-045 and SGR 1900+14, follow-up observations are needed to confirm our proposed candidates with higher confidence.

After 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March that lasted several months. During the active phase, several bursts were detected, and many aspects of the X-ray emission changed. We report on the discovery of six X-ray bursts, the first ever seen from this AXP in ~10 years of Rossi X-ray Timing Explorer (RXTE) monitoring. All the bursts occurred in the interval between 2006 April 6 and 2007 February 7. The bursts had the canonical fast rise slow decay profiles characteristic of SGR/AXP bursts. The burst durations ranged from 8-3x10^3 s as characterized by T90,these are very long durations even when compared to the broad T90 distributions of other bursts from SGRs and AXPs. The first five burst spectra are well modeled by simple blackbodies, with temperature kT ~2-6 keV. However, the sixth burst had a complicated spectrum consisting of at least three emission lines with possible additional emission and absorption lines. The most significant feature was at ~14 keV. Similar 14-keV spectral features were seen in bursts from AXPs 1E 1048.1-5937 and XTE J1810-197. If this feature is interpreted as a proton cyclotron line, then it supports the existence of a magnetar-strength field for these AXPs. Several of the bursts were accompanied by a short-term pulsed flux enhancement. We discuss these events in the context of the magnetar model.

We propose an explanation for the origin of hyperfast neutron stars (e.g. PSR B1508+55, PSR B2224+65, RX J0822-4300) based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star (or its helium core) which attained its peculiar velocity (similar to that of the neutron star) in the course of a strong three- or four-body dynamical encounter in the core of a young massive star cluster. This hypothesis implies that the dense cores of star clusters (located either in the Galactic disk or near the Galactic centre) could also produce the so-called hypervelocity stars -- the ordinary stars moving with a speed of ~1000 km/s.

The recent discovery of burst oscillation at 1122Hz in the X-ray transient XTE J1739-285 supports the suggestion that it contains a submillisecond pulsar\cite{1}. We here find for the first time the enormous dissipation effect in the transition boundary layer between quark matter core and hadron matter envelope. Just combining the estimation with previous dissipation mechanism together, we show that XTE J1739-285 can be uniquely restricted to a quark star masquerading as a neutron star (hybrid star) that contains a pure quark matter or mixed quark-hadron matter core from synthesizing both gravitational wave radiation (r-mode) instability and Keplerian motion constraints at 1122Hz lever. Such constraints allow the radii in the range $9{\rm km}\leq R\leq 12{\rm km}$ and the masses in the range $1.2M_\odot\leq M\leq 2.0M_\odot$. The normal neutron stars, hyperon stars and strange stars within the mass-radius limits are excluded.

Using the archival RXTE/ASM and SWIFT/BAT observations, the new orbital phases of Type I outbursts of EXO 2030+375 are estimated. A possible correlation between the Type II outburst and optical brightness variations is investigated. In order to estimate the phases of Type I outbursts, we fitted Gaussian profiles to the RXTE/ASM and SWIFT/BAT light curves. The time corresponding to the maximum value of the profiles is treated as the arrival time of Type I outburst. We used differential magnitudes in the time-series analysis of the optical light curve. MIDAS and its suitable packages were used to reduce and analyze the spectra. Prior to the Type II outburst, orbital phases of Type I outbursts were delayed for 6 days after the periastron passage, which is consistent with findings of Wilson et al., (2002, 2005). After the giant Type II outburst, the phase of Type I outbursts underwent a sudden shift of 13 days after the periastron passage. The amplitudes of Type I outbursts were increased between MJD 52500 and 53500. These amplitudes then decreased for 10 orbital cycles until the Type II outburst was triggered. If the change of outburst amplitudes correlated with the mass accretion, then during the decrease of these amplitudes mass should be deposited in a disk around neutron star temporarily. The release of this stored mass may ignite the Type II outburst. We report that the optical light curve became fainter by 0.4 mag during the decrease of amplitude of the Type I outbursts. The observed H$\alpha$ profiles and their equivalent widths during the decay and after the giant outburst are consistent with previous observations of the system.

We investigate the effect of antikaon condensed matter on bulk viscosity in rotating neutron stars. We use relativistic field theoretical models to construct the equation of state of neutron stars with the condensate, where the phase transition from nucleonic to $K^-$ condensed phase is assumed to be of first order. We calculate the coefficient of bulk viscosity due to the non-leptonic weak interaction n --> p + K^-. The influence of antikaon bulk viscosity on the gravitational radiation reaction driven instability in the r-modes is investigated. We compare our results with the previously studied non-leptonic weak interaction $n + p --> p + \Lambda$ involving hyperons on the damping of the r-mode oscillations.
We find that the bulk viscosity coefficient due to the non-leptonic weak process involving the condensate is suppressed by several orders of magnitude in comparison with the non-superfluid hyperon bulk viscosity coefficient. Consequently, the antikaon bulk viscosity may not be able to damp the r-mode instability, while hyperon bulk viscosity can effectively suppress r-mode oscillations at low temperatures. Hence neutron stars containing $K^-$ condensate in their core could be possible sources of gravitational waves.

The recent discovery of burst oscillation at 1122Hz in the X-ray transient XTE J1739-285 supports the suggestion that it contains a submillisecond pulsar\cite{1}. We here find for the first time the enormous dissipation effect in the transition boundary layer between quark matter core and hadron matter envelope. Just combining the estimation with previous dissipation mechanism together, we show that XTE J1739-285 can be uniquely restricted to a quark star masquerading as a neutron star (hybrid star) that contains a pure quark matter or mixed quark-hadron matter core from synthesizing both gravitational wave radiation (r-mode) instability and Keplerian motion constraints at 1122Hz lever. Such constraints allow the radii in the range $9{\rm km}\leq R\leq 12{\rm km}$ and the masses in the range $1.2M_\odot\leq M\leq 2.0M_\odot$. The normal neutron stars, hyperon stars and strange stars within the mass-radius limits are excluded.

Using the archival RXTE/ASM and SWIFT/BAT observations, the new orbital phases of Type I outbursts of EXO 2030+375 are estimated. A possible correlation between the Type II outburst and optical brightness variations is investigated. In order to estimate the phases of Type I outbursts, we fitted Gaussian profiles to the RXTE/ASM and SWIFT/BAT light curves. The time corresponding to the maximum value of the profiles is treated as the arrival time of Type I outburst. We used differential magnitudes in the time-series analysis of the optical light curve. MIDAS and its suitable packages were used to reduce and analyze the spectra. Prior to the Type II outburst, orbital phases of Type I outbursts were delayed for 6 days after the periastron passage, which is consistent with findings of Wilson et al., (2002, 2005). After the giant Type II outburst, the phase of Type I outbursts underwent a sudden shift of 13 days after the periastron passage. The amplitudes of Type I outbursts were increased between MJD 52500 and 53500. These amplitudes then decreased for 10 orbital cycles until the Type II outburst was triggered. If the change of outburst amplitudes correlated with the mass accretion, then during the decrease of these amplitudes mass should be deposited in a disk around neutron star temporarily. The release of this stored mass may ignite the Type II outburst. We report that the optical light curve became fainter by 0.4 mag during the decrease of amplitude of the Type I outbursts. The observed H$\alpha$ profiles and their equivalent widths during the decay and after the giant outburst are consistent with previous observations of the system.

We investigate the effect of antikaon condensed matter on bulk viscosity in rotating neutron stars. We use relativistic field theoretical models to construct the equation of state of neutron stars with the condensate, where the phase transition from nucleonic to $K^-$ condensed phase is assumed to be of first order. We calculate the coefficient of bulk viscosity due to the non-leptonic weak interaction n --> p + K^-. The influence of antikaon bulk viscosity on the gravitational radiation reaction driven instability in the r-modes is investigated. We compare our results with the previously studied non-leptonic weak interaction $n + p --> p + \Lambda$ involving hyperons on the damping of the r-mode oscillations.
We find that the bulk viscosity coefficient due to the non-leptonic weak process involving the condensate is suppressed by several orders of magnitude in comparison with the non-superfluid hyperon bulk viscosity coefficient. Consequently, the antikaon bulk viscosity may not be able to damp the r-mode instability, while hyperon bulk viscosity can effectively suppress r-mode oscillations at low temperatures. Hence neutron stars containing $K^-$ condensate in their core could be possible sources of gravitational waves.

We have observed the supernova remnant MSH 15-52 (G320.4-1.2), which contains the gamma-ray pulsar PSR B1509-58, using the CANGAROO-III imaging atmospheric Cherenkov telescope array from April to June in 2006. We detected gamma rays above 810 GeV at the 7 sigma level during a total effective exposure of 48.4 hours. We obtained a differential gamma-ray flux at 2.35 TeV of (7.9+/-1.5_{stat}+/-1.7_{sys}) \times 10^{-13} cm^{-2}s^{-1}TeV^{-1} with a photon index of 2.21+/-0.39_{stat}+/-0.40_{sys}, which is compatible with that of the H.E.S.S. observation in 2004. The morphology shows extended emission compared to our Point Spread Function. We consider the plausible origin of the high energy emission based on a multi-wavelength spectral analysis and energetics arguments.

Anomalous X-ray pulsars (AXPs) are thought to be magnetars, which are neutron stars with ultra strong magnetic field of $10^{14}$-- $10^{15}$ G. Their energy spectra below $\sim$10 keV are modeled well by two components consisting of a blackbody (BB) ($\sim$0.4 keV) and rather steep power-law (POW) function (photon index $\sim$2-4). Kuiper et al.(2004) discovered hard X-ray component above $\sim$10 keV from some AXPs. Here, we present the Suzaku observation of the AXP 1E 1841-045 at the center of supernova remnant Kes 73. By this observation, we could analyze the spectrum from 0.4 to 50 keV at the same time. Then, we could test whether the spectral model above was valid or not in this wide energy range. We found that there were residual in the spectral fits when fit by the model of BB + POW. Fits were improved by adding another BB or POW component. But the meaning of each component became ambiguous in the phase-resolved spectroscopy. Alternatively we found that NPEX model fit well for both phase-averaged spectrum and phase-resolved spectra. In this case, the photon indices were constant during all phase, and spectral variation seemed to be very clear. This fact suggests somewhat fundamental meaning for the emission from magnetars.

The model of pulsar radio emission is discussed in which a coherent radio emis-sion is excited in a vacuum gap above polar cap of neutron star. Pulsar X and gamma radiation are considered as the result of low-frequency radio emission inverse Comp-ton scattering on ultra relativistic electrons accelerated in the gap. The influence of the pulsar magnetic field on Compton scattering is taken into account. The relation of radio and gamma radiation spectra has been found in the framework of the model.

We explore the abundance of light clusters in core-collapse supernovae at post-bounce stage in a quantum statistical approach. Adopting the profile of a supernova core from detailed numerical simulations, we study the distribution of light bound clusters up to alpha particles (A=2-4) as well as heavy nuclei (A > 4) in dense matter at finite temperature. Within the frame of a cluster-mean field approach, the abundances of light clusters are evaluated accounting for self-energy, Pauli blocking and effects of continuum correlations. We find that deuterons and tritons, in addition to 3He and 4He, appear abundantly in a wide region from the surface of the proto-neutron star to the position of the shock wave. The appearance of light clusters may modify the neutrino emission in the cooling region and the neutrino absorption in the heating region, and, thereby, influence the supernova mechanism.

Over the last decade, considerable effort has been made to measure the proper motions of the pulsars B1757-24 and B1951+32 in order to establish or refute associations with nearby supernova remnants and to understand better the complicated geometries of their surrounding nebulae. We present proper motion measurements of both pulsars with the Very Large Array, increasing the time baselines of the measurements from 3.9 yr to 6.5 yr and from 12.0 yr to 14.5 yr, respectively, compared to previous observations. We confirm the non-detection of proper motion of PSR B1757-24, and our measurement of (mu_a, mu_d) = (-11 +/- 9, -1 +/- 15) mas yr^{-1} confirms that the association of PSR B1757-24 with SNR G5.4-1.2 is unlikely for the pulsar characteristic age of 15.5 kyr, although an association can not be excluded for a significantly larger age. For PSR B1951+32, we measure a proper motion of (mu_a, mu_d) = (-28.8 +/- 0.9, -14.7 +/- 0.9) mas yr^{-1}, reducing the uncertainty in the proper motion by a factor of two compared to previous results. After correcting to the local standard of rest, the proper motion indicates a kinetic age of ~51 kyr for the pulsar, assuming it was born near the geometric center of the supernova remnant. The radio-bright arc of emission along the pulsar proper motion vector shows time-variable structure, but moves with the pulsar at an approximately constant separation ~2.5", lending weight to its interpretation as a shock structure driven by the pulsar.

There is accumulating evidence that (fundamental) scalar fields may exist in Nature. The gravitational collapse of such a boson cloud would lead to a boson star (BS) as a new type of a compact object. Similarly as for white dwarfs and neutron stars, there exists a limiting mass, below which a BS is stable against complete gravitational collapse to a black hole. According to the form of the self-interaction of the basic constituents and the spacetime symmetry, we can distinguish mini-, axidilaton, soliton, charged, oscillating and rotating BSs. Their compactness prevents a Newtonian approximation, however, modifications of general relativity, as in the case of Jordan-Brans-Dicke theory as a low energy limit of strings, would provide them with gravitational memory. In general, a BS is a compact, completely regular configuration with structured layers due to the anisotropy of scalar matter, an exponentially decreasing 'halo', a critical mass inversely proportional to constituent mass, an effective radius, and a large particle number. Due to the Heisenberg principle, there exists a completely stable branch, and as a coherent state, it allows for rotating solutions with quantised angular momentum. In this review, we concentrate on the fascinating possibilities of detecting the various subtypes of (excited) BSs: Possible signals include gravitational redshift and (micro-)lensing, emission of gravitational waves, or, in the case of a giant BS, its dark matter contribution to the rotation curves of galactic halos.

The discovery of pulsars in 1968 heralded an era where the temporal characteristics of detectors had to be reassessed. Up to this point detector integration times would normally be measured in minutes rather seconds and definitely not on sub-second time scales. At the start of the 21st century pulsar observations are still pushing the limits of detector telescope capabilities. Flux variations on times scales less than 1 nsec have been observed during giant radio pulses. Pulsar studies over the next 10 to 20 years will require instruments with time resolutions down to microseconds and below, high-quantum quantum efficiency, reasonable energy resolution and sensitive to circular and linear polarisation of stochastic signals. This chapter is review of temporally resolved optical observations of pulsars. It concludes with estimates of the observability of pulsars with both existing telescopes and into the ELT era.

We report observations of Crab giant pulses made with the Australia Telescope Compact Array and a baseband recorder system, made simultaneously at two frequencies, 1300 and 1470 MHz. These observations were sensitive to pulses with amplitudes \ga 3 kJy and widths \ga 0.5 $\mu$s. Our analysis led to the detection of more than 700 such bright giant pulses over 3 hours, and using this large sample we investigate their amplitude, width, arrival time and energy distributions. The brightest pulse detected in our data has a peak amplitude of $\sim$ 45 kJy and a width of $\sim$ 0.5 $\mu$s, and therefore an inferred brightness temperature of $\sim 10^{35}$ K. The duration of giant-pulse emission is typically $\sim$1 $\mu$s, however it can also be as long as 10 $\mu$s. The pulse shape at a high time resolution (128 ns) shows rich diversity and complexity in structure and is marked by an unusually low degree of scattering. We discuss possible implications for scattering due to the nebula, and for underlying structures and electron densities.

The thermal evolution of neutron stars is coupled to their spin down and the resulting changes in structure and chemical composition. This coupling correlates stellar surface temperatures with rotational state as well as time. We report an extensive investigation of the coupling between spin down and cooling for hybrid stars which undergo a phase transition to deconfined quark matter at the high densities present in stars at low rotation frequencies. The thermal balance of neutron stars is reanalyzed to incorporate phase transitions and the related latent heat self-consistently, and numerical calculations are undertaken to simultaneously evolve the stellar structure and temperature distribution. We find that the changes in stellar structure and chemical composition with the introduction of a pure quark matter phase in the core delay the cooling and produce a period of increasing surface temperature for strongly superfluid stars of strong and intermediate magnetic field strength. The latent heat of deconfinement is found to reinforce this signature but is itself relatively less significant.

Forty years after their discovery, and in spite of a very large body of observations, the operation of the 'engine' responsible for long-duration Gamma-Ray Bursts (GRBs) and X-ray flashes --as well as the mechanisms generating their radiation-- are still the subject of debate and study. In this respect a recent event, XRF 060218, associated with SN 2006aj, is particularly significant. It has been argued that, for the first time, the break-out of the shock involved in the supernova explosion has been observed, thanks to the detection of a thermal component in the event's radiation; that this XRF was not a GRB seen 'off-axis', but a member of a new class of energetically feeble GRBs; and that its 'continued engine activity' may have been driven by a remnant highly-magnetized neutron star, a magnetar. I argue, on grounds based on observations and on limpid verified hypothesis, that there is a common, simpler alternative to these views, with no thermal component, no new feeble GRBs, and no steady engine activity.

There is ${}^3P_2$ neutron superfluid region in NS (neutron star) interior. For a rotating NS, the ${}^3P_2$ superfluid region is like a system of rotating magnetic dipoles. It will give out electromagnetic radiation, which may provides a new heating mechanism of NSs. This heating mechanism plus some cooling agent may give sound explanation to NS glitches.