When a rotating neutron star loses angular momentum, the reduction in the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium that enhance the neutrino emissivity and produce thermal energy. This mechanism is named rotochemical heating and has previously been studied for neutron stars of non-superfluid matter, finding that they reach a quasi-steady state in which the rate that the spin-down modifies the equilibrium concentrations is the same to that of the neutrino reactions restoring the equilibrium. On the other hand, the neutron star interior is believed to contain superfluid nucleons, which affect the thermal evolution of the star by suppressing the neutrino reactions and the specific heat, and opening new Cooper pairing reactions.
In this work we describe the thermal effects of Cooper pairing with spatially uniform energy gaps of neutrons and protons on rotochemical heating in millisecond pulsars (MSPs) when only modified Urca reactions are allowed. We find that the chemical imbalances grow up to a value close to the energy gaps, which is higher than the one of the nonsuperfluid case. Therefore, the surface temperatures predicted with Cooper pairing are higher and explain the recent measurement of MSP J0437-4715.

The thermal emission detected from the millisecond pulsar J0437-4715 is not explained by standard cooling models of neutron stars without a heating mechanism. We investigated three heating mechanisms controlled by the rotational braking of the pulsar: breaking of the solid crust, superfluid vortex creep, and non-equilibrium reactions ('rotochemical heating'). We find that the crust cracking mechanism does not produce detectable heating. Given the dependence of the heating mechanisms on spin-down parameters, which leads to different temperatures for different pulsars, we study the thermal evolution for two types of pulsars: young, slowly rotating 'classical' pulsars and old, fast rotating millisecond pulsars (MSPs). We find that the rotochemical heating and vortex creep mechanism can be important both for classical pulsars and MSPs.

The concept of "self-organized criticality" (SOC) has been introduced by Bak, Tang, and Wiesenfeld (1987) to describe the statistics of avalanches on the surface of a sandpile with a critical slope, which produces a scale-free powerlaw size distribution of avalanches. In the meantime, SOC behavior has been identified in many nonlinear dissipative systems that are driven to a critical state. On a most general level, SOC is the statistics of coherent nonlinear processes, in contrast to the Poisson statistics of incoherent random processes. The SOC concept has been applied to laboratory experiments (of rice or sand piles), to human activities (population growth, language, economy, traffic jams, wars), to biophysics, geophysics (earthquakes, landslides, forest fires), magnetospheric physics, solar physics (flares), stellar physics (flares, cataclysmic variables, accretion disks, black holes, pulsar glitches, gamma ray bursts), and to galactic physics and cosmology.

As neutron stars spin-down and contract, the deconfinement phase transition can continue to occur, resulting in energy release(so-called deconfinement heating) in case of the first-order phase transition. The thermal evolution of neutron stars is investigated to combine phase transition and the related energy release self-consistently. We find that the appearance of deconfinement heating during spin-down result in not only the cooling delay but also the increase of surface temperature of stars. For stars characterized by intermediate and weak magnetic field strength, a period of increasing surface temperature could exist. Especially, a sharp jump in surface temperature can be produced as soon as quark matter appears in the core of stars with a weak magnetic field. We think that this may serve as evidence for the existence of deconfinement quark matter. The results show that deconfinement heating facilitates the emergence of such characteristic signature during the thermal evolution process of neutron stars.

We have timed four millisecond pulses, PSRs J1721-2457, J1745-0952, J1810-2005, and J1918-0642, for up to a total of 10.5 years each using multiple telescopes in the European Pulsar Timing Array network: the Westerbork Synthesis Radio Telescope in The Netherlands, the Nancay Radio Telescope in France and the Lovell telescope at Jodrell Bank in the UK. The long time span has enabled us to measure the proper motions of J1745-0952 and J1918-0642, indicating that they have transverse velocities of 200(50) and 54(7) km/s respectively. We have obtained upper limits on the proper motion of J1721-2457 and J1810-2005, which imply that they have transverse velocities less than 140 and 400 km/s respectively. In all cases, the velocities lie in the range typical of millisecond pulsars. We present pulse profiles for each pulsar taken from observations at multiple frequencies in the range of 350 to 2600 MHz, and show that J1810-2005 shows significant profile evolution in this range. Using our multi-frequency observations, we measured the spectral indices for all four pulsars, and for J1810-2005 it appears to be very flat. The flux density of J1918-0642 shows extensive modulation which we attribute to the combined effects of refractive and diffractive scintillation. We discuss the possible use of including J1721-2457 or J1918-0642 in a pulsar timing array, and find that J1918-0642 will be useful to include when the timing precision of this pulsar is improved over the next few years. We have searched archival optical observations to detect companions of the binary pulsars, but none were detected. However, we provide lower limits on the masses of the white dwarf companions of PSRs J1745-0952 and J1918-0642.

The observation of massive black hole binaries (MBHBs) with Pulsar Timing Arrays (PTAs) is one of the goals of gravitational wave astronomy in the coming years. Massive (>10^8 solar masses) and low-redshift (< 1.5) sources are expected to be individually resolved by up-coming PTAs, and our ability to use them as astrophysical probes will depend on the accuracy with which their parameters can be measured. In this paper we estimate the precision of such measurements using the Fisher-information-matrix formalism. We restrict to "monochromatic" sources. In this approximation, the system is described by seven parameters and we determine their expected statistical errors as a function of the number of pulsars in the array, the array sky coverage, and the signal-to-noise ratio (SNR) of the signal. At fixed SNR, the gravitational wave astronomy capability of a PTA is achieved with ~20 pulsars; adding more pulsars (up to 1000) to the array reduces the source error-box in the sky \Delta\Omega by a factor ~5 and has negligible consequences on the statistical errors on the other parameters. \Delta\Omega improves as 1/SNR^2 and the other parameters as 1/SNR. For a fiducial PTA of 100 pulsars uniformly distributed in the sky and a coherent SNR = 10, we find \Delta\Omega~40 deg^2, a fractional error on the signal amplitude of ~30% (which constraints only very poorly the chirp mass - luminosity distance combination M_c^{5/3}/D_L), and the source inclination and polarization angles are recovered at the ~0.3 rad level. The ongoing Parkes PTA is particularly sensitive to systems located in the southern hemisphere, where at SNR = 10 the source position can be determined with \Delta\Omega ~10 deg^2, but has poorer performance for sources in the northern hemisphere. (Abridged)

Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if they have a neutral point in the spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations in full general relativity, determine the approximate locations of the neutral points, and derive limits on the observable frequency band available to the instability in these sequences. We find that general relativity enhances the detectability of a CFS-unstable neutron star substantially, both by widening the instability window and enlarging the band into the optimal range for interferometric detectors like LIGO, VIRGO, and GEO-600.

The study of the younger, and brighter, pulsars is important to understand the optical emission properties of isolated neutron stars. PSRB0540-69, the second brightest (V~22) optical pulsar, is obviously a very interesting target for these investigations. The aim of this work is threefold: constraining the pulsar proper motion and its velocity on the plane of the sky through optical astrometry, obtaining a more precise characterisation of the pulsar optical spectral energy distribution (SED) through a consistent set of multi-band, high-resolution, imaging photometry observations, measuring the pulsar optical phase-averaged linear polarisation, for which only a preliminary and uncertain measurement was obtained so far from ground-based observations. We performed high-resolution observations of PSRB0540-69 with the WFPC2 aboard the HST, in both direct imaging and polarimetry modes. From multi-epoch astrometry we set a 3sigma upper limit of 1 mas/yr on the pulsar proper motion, implying a transverse velocity <250 km/s at the 50 kpc LMC distance. Moreover, we determined the pulsar absolute position with an unprecedented accuracy of 70 mas. From multi-band photometry we characterised the pulsar power-law spectrum and we derived the most accurate measurement of the spectral index (0.70+/-0.07) which indicates a spectral turnover between the optical and X-ray bands. Finally, from polarimetry we obtained a new measurement of the pulsar phase-averaged polarisation degree (16+/-4%),consistent with magnetosphere models depending on the actual intrinsic polarisation degree and depolarisation factor, and we found that the polarisation vector (22+/-12deg position angle) is possibly aligned with the semi-major axis of the pulsar-wind nebula and with the apparent proper motion direction of its bright emission knot.

Some of radiopulsars have anomalous braking index values $n = \Omega \ddot{\Omega} / \dot{\Omega}^2 \sim \pm (10^3 \div 10^4) $. It is shown that such values may be related with nondipolar magnetic field. The precession of neutron star lead to rotation (in reference frame related with neutron star) of vector of angular velocity $\vec{\Omega}$ along direction of neutron star magnetic dipole moment $\vec{m}$ with angular velocity $\vec{\Omega}_{p}$. This process may cause the altering of electric current flow through inner gap and consequently the current losses with the same time scale as precession period $T_{p} = 2\pi / \Omega_{p}$. It occurs because of electric current in inner gaps is determined by Goldreich-Julian charge density $\rho_{GJ} = -\frac{\vec{\Omega} \cdot \vec{B}}{2\pi c}$, which are depend on angle between direction of small scale magnetic field and angular velocity $\vec{\Omega}$. It is essential that pulsar tubes nearby neutron star surface are curved. In current paper it is considered the only inner gaps with steady, electron charge limited flow regime.

Observations using X-ray telescopes can help understand the origin of the electron and positron signals reported by ATIC, PAMELA, PPB-BETS, and Fermi. It remains unclear whether the observed high-energy electrons and positrons are produced by the relic particles, or by some astrophysical sources. To distinguish between the two possibilities, one can compare the electron population in the local neighborhood with that in the dwarf spheroidal galaxies, which are not expected to host as many pulsars and other astrophysical sources. This can be accomplished using the X-ray observations of the dwarf spheroidal galaxies. Assuming the Fermi signal comes from dark matter and using the inferred dark matter profile of the Draco dwarf spheroidal galaxy, we calculate the spectrum of X-rays produced by electrons via inverse Compton scattering. The next generation of X-ray telescopes may be able to detect such a signal.

The objects known as anomalous X-ray pulsars and soft gamma repeaters are commonly identified with magnetars, neutron stars with ultrastrong magnetic fields. The rotational history of these objects has, so far, revealed no evidence of free precession. At the same time these objects do not generally appear to have magnetic axes nearly parallel or orthogonal to their spin axes. In this paper we show that the combination of these two observations, together with simple rigid-body dynamics, leads to non-trivial predictions about the interior properties of magnetars: either (i) elastic stresses in magnetar crusts are close to the theoretical upper limit above which the crustal matter yields or (ii) there is a "pinned" superfluid component in the magnetar interior. As a potentially observable consequence of these ideas we point out that, in the case of no pinned superfluidity, magnetars of stronger magnetic field strength than those currently observed would have to be nearly aligned/orthogonal rotators.

A recent research shows that particles with a spectrum of a relativistic Maxwellian plus a high-energy tail can be accelerated by relativistic collisionless shocks. We investigate the possibility of the high-energy particles with this new spectrum injected in pulsar wind nebulae (PWNe) from the terminate shock based on the study of multiwavelength emission from PWNe.} {The dynamics of a supernova remnant (SNR) and multiband nonthermal emission from the PWN inside the remnant are investigated using a dynamical model with electrons/positrons injected with the new spectrum. In this model, the dynamical and radiative evolution of a pulsar wind nebula in a non-radiative supernova remnant can be self-consistently described.} {This model is applied to the three composite SNRs, G0.9+0.1, MSH 15-52, G338.3-0.0, and the multiband observed emission from the three PWNe can be well reproduced.} {Our studies on the three remnant provide evidence for the new spectrum of the particles, which are accelerated by the terminate shock, injected into a PWN.

We present an X-ray morphological and spectroscopic study of the pulsar B2224+65 and its apparent jet-like X-ray features based on two epoch Chandra observations. The main X-ray feature, which shows a large directional offset from the ram-pressure confined pulsar wind nebula (Guitar Nebula), is broader in apparent width and more luminous in the second epoch than the first. Furthermore, the sharp leading edge is found to have a proper motion consistent with that of the pulsar (~180 mas/yr). The combined data set also provides evidence for the presence of a counter feature, albeit substantially fainter and shorter than the main one. These results are consistent with a simple model of relativistic jet outflow originating from the pulsar and ram-pressure confined by the unusually rapid motion of the pulsar.

A stalled spherical accretion shock, such as that arising in core-collapse supernovae, is unstable to non-spherical perturbations. In three dimensions, this Standing Accretion Shock Instability (SASI) has been observed to develop spiral modes that can impart a sizable torque to the protoneutron star. Previous studies have arrived at diverging conclusions regarding the need for rotation in the collapsing progenitor to trigger these modes. Here we combine linear stability analysis and three-dimensional, time-dependent hydrodynamic simulations with Zeus-MP to study these non-axisymmetric modes in the linear regime. We do not impose any rotation on the accretion flow, and use simplified microphysics with no neutrino heating or nuclear dissociation. The three-dimensional structure of the linear eigenmodes is constructed, and their evolution is compared with time-dependent simulations. We show that spiral modes are most easily understood as a superposition of sloshing modes out of phase. In the absence of rotation, there is no preferred direction in the system, hence multiple sloshing modes with different relative phases can be combined, resulting in a larger class of spiral-like modes that require less specific conditions for their excitation than the standard cases.

We report on the discovery of a molecular cavity in the Norma near arm in the general direction of Westerlund 1 (Wd1), but not associated with it. The cavity has a mean radial velocity of -91.5 kms^{-1}, which differs by as much as ~40 kms^{-1} from the mean radial velocity of the Wd1 stars. The cavity is surrounded by a fragmented molecular shell of an outer diameter of about 100 pc and 10^{6} M_odot, which is expanding at velocities of 6 to 8 kms^{-1}. The amount of kinetic energy involved in the expanding shell is ~10^{51} erg. Inside this cavity the atomic HI gas surface density is also the lowest. Structure of the extended Very High Energetic (VHE) gamma-ray emission, recently reported by the H.E.S.S. collaboration Ohm et al. 2009, coincides with the cavity. The observed morphology suggests that the inner wall of the molecular shell is the zone of the gamma-ray emission, and not the dense gas surrounding massive stars of Wd1 as had been speculated by the H.E.S.S. collaboration. A likely candidate responsible for creating the observed cavity and the gamma-ray emission is the pulsar PSR J1648-4611.

When a rotating neutron star loses angular momentum, the reduction of the centrifugal force makes it contract. This perturbs each fluid element, raising the local pressure and originating deviations from beta equilibrium, inducing reactions that release heat (rotochemical heating). This effect has previously been studied by Fern\'andez and Reisenegger for neutron stars of non-superfluid matter and by Petrovich and Reisenegger for superfluid matter, finding that the system in both cases reaches a quasi-steady state, corresponding to a partial equilibration between compression, due to the loss of angular momentum, and reactions that try to restore the equilibrium. However, Petrovich and Reisenegger assumes a constant value of the superfluid energy gap, whereas theoretical models predict density-dependent gap amplitudes, and therefore gaps that depend on the location in the star. In this work, we try to discriminate between several proposed gap models, comparing predicted surface temperatures to the value measured for the nearest millisecond pulsar, J0437-4715.

Anomalous X-ray Pulsars (AXPs) are now established to exhibit significant X-ray variability and be prolific glitchers, with some glitches being accompanied by large radiative changes. An open issue is whether AXP glitches are generically accompanied by radiative changes, relevant for understanding magnetar physical properties. Here we report on an analysis of archival X-ray data from the AXP 1E~1841$-$045, obtained between 1993 and 2007. This AXP, located in the center of SNR Kes~73, has exhibited three glitches between 2002 and 2007, as determined by {\it RXTE} monitoring since 1999. We have searched for evidence of phase-averaged flux variability that could be present if glitches in AXPs are usually accompanied by radiative changes. We find no evidence for glitch-correlated flux changes from this source, arguing that such behavior is not generic to AXPs.

We report on the first near-infrared observations obtained to date for Rotating RAdio Transients (RRATs). Using adaptive optics devices mounted on the ESO Very Large Telescope (VLT), we observed two objects of this class: RRAT J1819-1458, and RRAT J1317-5759. These observations have been performed in 2006 and 2008, in the J, H and Ks bands. We found no candidate infrared counterpart to RRAT J1317-5759, down to a limiting magnitude of Ks ~ 21. On the other hand, we found a possible candidate counterpart for RRAT J1819-1458, having a magnitude of Ks=20.96+/-0.10 . In particular, this is the only source within a 1 sigma error circle around the source's accurate X-ray position, although given the crowded field we cannot exclude that this is due to a chance coincidence. The infrared flux of the putative counterpart to the highly magnetic RRAT J1819-1458, is higher than expected from a normal radio pulsar, but consistent with that seen from magnetars. We also searched for the near-infrared counterpart to the X-ray diffuse emission recently discovered around RRAT J1819-1458, but we did not detect this component in the near-infrared band. We discuss the luminosity of the putative counterpart to RRAT J1819-1458, in comparison with the near-infrared emission of all isolated neutron stars detected to date in this band (5 pulsars and 7 magnetars).

We construct general relativistic models of stationary, strongly magnetized neutron stars. The magnetic field configuration, obtained by solving the relativistic Grad-Shafranov equation, is a generalization of the twisted torus model recently proposed in the literature; the stellar deformations induced by the magnetic field are computed by solving the perturbed Einstein's equations; stellar matter is modeled using realistic equations of state. We find that in these configurations the poloidal field dominates over the toroidal field and that, if the magnetic field is sufficiently strong during the first phases of the stellar life, it can produce large deformations.

SDSS 1257+5428 is a white dwarf in a close orbit with a companion that has been suggested to be a neutron star. If so, it hosts the closest known neutron star, and its existence implies a great abundance of similar systems and a rate of white-dwarf neutron-star mergers similar to that of the type Ia supernova rate. Here, we present high signal-to-noise spectra of SDSS 1257+5428, which confirm an independent finding that the system is in fact composed of two white dwarfs, one relatively cool and with low mass, and the other hotter and more massive. With this, the demographics and merger rate are no longer puzzling (various factors combine to lower the latter by more than two orders of magnitude). We show that the spectra are fit well with a combination of two hydrogen model atmospheres, as long as the lines of the higher-gravity component are broadened significantly relative to what is expected from just pressure broadening. Interpreting this additional broadening as due to rotation, the inferred spin period is short, about 1 minute. Similarly rapid rotation is only seen in accreting white dwarfs that are magnetic; empirically, it appears that in non-magnetized white dwarfs, accreted angular momentum is lost by nova explosions before it can be transferred to the white dwarf. This suggests that the massive white dwarf in SDSS 1257+5428 is magnetic as well, with B~10^5 G. Alternatively, the broadening seen in the spectral lines could be due to a stronger magnetic field, of ~10^6 G. The two models could be distinguished by further observations.

Explosive astrophysical systems, such as supernovae or compact star binary mergers, provide conditions where strange quark matter can appear. The high degree of isospin asymmetry and temperatures of several MeV in such systems may cause a transition to the quark phase already around saturation density. Observable signals from the appearance of quark matter can be predicted and studied in astrophysical simulations. As input in such simulations, an equation of state with an integrated quark matter phase transition for a large temperature, density and proton fraction range is required. Additionally, restrictions from heavy ion data and pulsar observation must be considered. In this work we present such an approach. We implement a quark matter phase transition in a hadronic equation of state widely used for astrophysical simulations and discuss its compatibility with heavy ion collisions and pulsar data. Furthermore, we review the recently studied implications of the QCD phase transition during the early post-bounce evolution of core-collapse supernovae and introduce the effects from strong interactions to increase the maximum mass of hybrid stars. In the MIT bag model, together with the strange quark mass and the bag constant, the strong coupling constant $\alpha_s$ provides a parameter to set the beginning and extension of the quark phase and with this the mass and radius of hybrid stars.

Massive stars are of interest as progenitors of super novae, i.e. neutron stars and black holes, which can be sources of gravitational waves. Recent population synthesis models can predict neutron star and gravitational wave observations but deal with a fixed super nova rate or an assumed initial mass function for the population of massive stars. Here we investigate those massive stars, which are supernova progenitors, i.e. with O and early B type stars, and also all super giants within 3kpc. We restrict our sample to those massive stars detected both in 2MASS and observed by Hipparcos, i.e. only those stars with parallax and precise photometry. To determine the luminosities we calculated the extinctions from published multi-colour photometry, spectral types, luminosity class, all corrected for multiplicity and recently revised Hipparcos distances. We use luminosities and temperatures to estimate the masses and ages of these stars using different models from different authors. Having estimated the luminosities of all our stars within 3kpc, in particular for all O- and early B-type stars, we have determined the median and mean luminosities for all spectral types for luminosity classes I, III, and V. Our luminosity values for super giants deviate from earlier results: Previous work generally overestimates distances and luminosities compared to our data, this is likely due to Hipparcos parallaxes (generally more accurate and larger than previous ground-based data) and the fact that many massive stars have recently been resolved into multiples of lower masses and luminosities. From luminosities and effective temperatures we derived masses and ages using mass tracks and isochrones from different authors. From masses and ages we estimated lifetimes and derived a lower limit for the supernova rate of ~20 events/Myr averaged over the next 10 Myrs within 600 pc from the sun. These data are then used to search for areas in the sky with higher likelihood for a supernova or gravitational wave event (like OB associations).

We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.

In order to examine the rotational effect around neutron star in tensor-vector-scalar (TeVeS) theory, we consider the slowly rotating relativistic stars with a uniform angular velocity. As a result, we find that similar to the case in general relativity (GR), the angular momentum is proportional to the angular velocity. Additionally, as the value of coupling constant $K$ becomes higher, the frame dragging in TeVeS becomes quite different distribution from that in GR, where we can also see the deviation even in the interior of star. While with smaller value of $K$, although the frame dragging approaches to that expected in GR, the induced vector field due to the rotation does not vanish and still exists. Thus, through the observations associated with relativistic object, one could be possible to distinguish the gravitational theory in strong field regime even in the case that the value of coupling constant $K$ is quite small.

We recently presented a series of dark energy theorems that place constraints on the equation of state of dark energy ($\wdark$), the ime-variation of Newton's constant ($\dot G$), and the violation of energy conditions in theories with extra dimensions. In this paper, we explore how current and future measurements of $\wdark$ and $\dot G$ can be used to place tight limits on large classes of these theories (including some of the most well-motivated examples) independent of the size of the extra dimensions. As an example, we show that models with conformally Ricci-flat metrics obeying the null energy condition (a common ansatz for Kaluza-Klein and string constructions) are highly constrained by current ata and may be ruled out entirely by future dark energy and pulsar observations.

From a deep multi-epoch Chandra observation of the elliptical galaxy NGC 3379 we report the spectral properties of nine luminous LMXBs (LX>1E38 erg/s). We also present a set of spectral simulations, produced to aid the interpretation of low-count single-component spectral modeling. These simulations demonstrate that it is possible to infer the spectral states of X-ray binaries from these simple models and thereby constrain the properties of the source. Of the nine LMXBs studied, three reside within globular clusters, and one is a confirmed field source. Due to the nature of the luminosity cut all sources are either neutron star binaries emitting at or above the Eddington luminosity or black hole binaries. The spectra from these sources are well described by single-component models, with parameters consistent with Galactic LMXB observations, where hard-state sources have a range in photon index of 1.4-1.9 and thermally dominated sources have inner disc temperatures between ~0.7-1.55 keV.
The large variability observed in the brightest globular cluster source (LX>4E38 erg/s) suggests the presence of a black hole binary. At its most luminous this source is observed in a thermally dominated state with kT=1.5 keV, consistent with a black hole mass of 10 Msol. This observation provides further evidence that globular clusters are able to retain such massive binaries. We also observed a source transitioning from a bright state (LX~1E39 erg/s), with prominent thermal and non-thermal components, to a less luminous hard state (LX=4.5E38 erg/s, Gamma=1.85). In its high flux emission this source exhibits a cool-disc component of ~0.14 keV, similar to spectra observed in some ultraluminous X-ray sources. Such a similarity indicates a possible link between `normal' stellar mass black holes in a high accretion state and ULXs.

The stability of the color flavor locked phase in the presence of a strong magnetic field is investigated within the phenomenological MIT bag model, taking into account the variation of the strange quark mass, the baryon density, the magnetic field, as well as the bag and gap parameters. It is found that the minimum value of the energy per baryon in a color flavor locked state at vanishing pressure is lower than the corresponding one for unpaired magnetized strange quark matter and, as the magnetic field increases, the energy per baryon decreases. This implies that magnetized color flavor locked matter is more stable and could become the ground state inside neutron stars. The mass-radius relation for such stars is also studied.

Galaxies are missing most of their baryons, and many models predict these baryons lie in a hot halo around galaxies. We establish observationally motivated constraints on the mass and radii of these haloes using a variety of independent arguments. First, the observed dispersion measure of pulsars in the Large Magellanic Cloud allows us to constrain the hot halo around the Milky Way: if it obeys the standard NFW profile, it must contain less than 4-5% of the missing baryons from the Galaxy. This is similar to other upper limits on the Galactic hot halo, such as the soft X-ray background and the pressure around high velocity clouds. Second, we note that the X-ray surface brightness of hot haloes with NFW profiles around large isolated galaxies is high enough that such emission should be observed, unless their haloes contain less than 10-25% of their missing baryons. Third, we place constraints on the column density of hot haloes using nondetections of OVII absorption along AGN sightlines: in general they must contain less than 70% of the missing baryons or extend to no more than 40 kpc. Flattening the density profile of galactic hot haloes weakens the surface brightness constraint so that a typical L$_*$ galaxy may hold half its missing baryons in its halo, but the OVII constraint remains unchanged, and around the Milky Way a flattened profile may only hold $6-13%$ of the missing baryons from the Galaxy ($2-4 \times 10^{10} M_{\odot}$). We also show that AGN and supernovae at low to moderate redshift - the theoretical sources of winds responsible for driving out the missing baryons - do not produce the expected correlations with the baryonic Tully-Fisher relationship and so are insufficient to explain the missing baryons from galaxies. We conclude that most of missing baryons from galaxies do not lie in hot haloes around the galaxies, and that the missing baryons never fell into the potential wells of protogalaxies in the first place. They may have been expelled from the galaxies as part of the process of galaxy formation.

The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.

We present Rossi X-Ray Timing Explorer and Swift observations made during the final three weeks of the 2006-2007 outburst of the super-Eddington neutron star transient XTE J1701-462, as well as Chandra and XMM-Newton observations covering the first ~800 days of the subsequent quiescent phase. The source transitioned quickly from active accretion to quiescence, with the luminosity dropping by over three orders of magnitude in ~13 days. The spectra obtained during quiescence exhibit both a thermal component, presumed to originate in emission from the neutron star surface, and a non-thermal component of uncertain origin, which has shown large and irregular variability. We interpret the observed decay of the inferred effective surface temperature of the neutron star in quiescence as the cooling of the neutron star crust after having been heated and brought out of thermal equilibrium with the core during the outburst. The interpretation of the data is complicated by an apparent temporary increase in temperature ~220 days into quiescence, possibly due to an additional spurt of accretion. We derive an exponential decay timescale of ~120 (+30/-20) days for the inferred temperature (excluding observations affected by the temporary increase). This short timescale indicates a highly conductive neutron star crust. Further observations are needed to confirm whether the crust is still slowly cooling or has already reached thermal equilibrium with the core at a surface temperature of ~125 eV. The latter would imply a high equilibrium bolometric thermal luminosity of ~5x10^{33} erg/s for an assumed distance of 8.8 kpc.

We report on a detailed abundance analysis of the strongly r-process enhanced giant star, HE 2327-5642 ([Fe/H] = -2.78, [r/Fe] = +0.99). Determination of stellar parameters and element abundances was based on analysis of high-quality VLT/UVES spectra. The surface gravity was calculated from the NLTE ionization balance between Fe I and Fe II, and Ca I and Ca II. Accurate abundances for a total of 40 elements and for 23 neutron-capture elements beyond Sr and up to Th were determined. The heavy element abundance pattern of HE 2327-5642 is in excellent agreement with those previously derived for other strongly r-process enhanced stars. Elements in the range from Ba to Hf match the scaled Solar r-process pattern very well. No firm conclusion can be drawn with respect to a relationship between the fisrt neutron-capture peak elements, Sr to Pd, in HE 2327-5642 and the Solar r-process, due to the uncertainty of the latter. A clear distinction in Sr/Eu abundance ratios was found between the halo stars with different europium enhancement. The strongly r-process enhanced stars reveal a low Sr/Eu abundance ratio at [Sr/Eu] = -0.92+-0.13, while the stars with 0 < [Eu/Fe] < 1 and [Eu/Fe] < 0 have 0.36 dex and 0.93 dex larger Sr/Eu values, respectively. Radioactive dating for HE 2327-5642 with the observed thorium and rare-earth element abundance pairs results in an average age of 13.3 Gyr, when based on the high-entropy wind calculations, and 5.9 Gyr, when using the Solar r-residuals. HE 2327-5642 is suspected to be radial-velocity variable based on our high-resolution spectra, covering ~4.3 years.

Globular clusters with their large populations of millisecond pulsars (MSPs) are believed to be potential emitters of high-energy gamma-ray emission. The observation of this emission provides a powerful tool to assess the millisecond pulsar population of a cluster, is essential for understanding the importance of binary systems for the evolution of globular clusters, and provides complementary insights into magnetospheric emission processes. Our goal is to constrain the millisecond pulsar populations in globular clusters from analysis of gamma-ray observations. We use 546 days of continuous sky-survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope to study the gamma-ray emission towards 13 globular clusters. Steady point-like high-energy gamma-ray emission has been significantly detected towards 8 globular clusters. Five of them (47 Tucanae, Omega Cen, NGC 6388, Terzan 5, and M 28) show hard spectral power indices $(0.7 < \Gamma <1.4)$ and clear evidence for an exponential cut-off in the range 1.0-2.6 GeV, which is the characteristic signature of magnetospheric emission from MSPs. Three of them (M 62, NGC 6440 and NGC 6652) also show hard spectral indices $(1.0 < \Gamma < 1.7)$, however the presence of an exponential cut-off can not be unambiguously established. Three of them (Omega Cen, NGC 6388, NGC 6652) have no known radio or X-ray MSPs yet still exhibit MSP spectral properties. From the observed gamma-ray luminosities, we estimate the total number of MSPs that is expected to be present in these globular clusters. We show that our estimates of the MSP population correlate with the stellar encounter rate and we estimate 2600-4700 MSPs in Galactic globular clusters, commensurate with previous estimates. The observation of high-energy gamma-ray emission from a globular cluster thus provides a reliable independent method to assess their millisecond pulsar populations that can be used to make constraints on the original neutron star X-ray binary population, essential for understanding the importance of binary systems in slowing the inevitable core collapse of globular clusters.

We report the detection of high energy gamma-ray emission from the young and energetic pulsar PSR B1509$-$58 and its pulsar wind nebula (PWN) in the composite supernova remnant SNR G320.4-1.2 (aka MSH 15-52). Using 1 year of survey data with the Fermi-Large Area Telescope (LAT), we detected pulsations from PSR B1509-58 up to 1 GeV and extended gamma-ray emission above 1 GeV spatially coincident with the PWN. The pulsar light curve presents two peaks offset from the radio peak by phases 0.96 $\pm$ 0.01 and 0.33 $\pm$ 0.02. New constraining upper limits on the pulsar emission are derived below 1 GeV and confirm a severe spectral break at a few tens of MeV. The nebular spectrum in the 1 - 100 GeV energy range is well described by a power-law with a spectral index of (1.57 $\pm$ 0.17 $\pm$ 0.13) and a flux above 1 GeV of (2.91 $\pm$ 0.79 $\pm$ 1.35) 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. The LAT spectrum of the nebula connects nicely with Cherenkov observations, and indicates a spectral break between GeV and TeV energies.

Using new and archival observations made with the Swift satellite and other facilities, we examine 147 X-ray sources selected from the ROSAT All-Sky-Survey Bright Source Catalog (RASS/BSC) to produce a new limit on the number of isolated neutron stars (INSs) in the RASS/BSC, the most constraining such limit to-date. Independent of X-ray spectrum and variability, the number of INSs is <=48 (90% confidence). Restricting attention to soft (having an effective temperature of < 200 eV), non-variable X-ray sources -- as in a previous study -- yields an all-sky limit of <=31 INSs. In the course of our analysis, we identify five new high-quality INS candidates for targeted follow-up observations. A future all-sky X-ray survey with eROSITA, or another mission with similar capabilities, can be expected to increase the detected population of X-ray-discovered INSs from the 8 to 50 in the BSC, to (for a disk population) 240 to 1500, which will enable a more detailed study of neutron star population models.

The cumulative luminosity distribution functions (CLFs) of radio millisecond pulsars (MSPs) in globular clusters (GCs) and in the Galactic field at a frequency of 1.4 GHz have been examined. Assuming a functional form, $N \propto L^q$ where $N$ is the number of MSPs and $L$ is the luminosity at 1.4 GHz, it is found that the CLFs significantly differ with a steeper slope, $q=-0.83 \pm 0.05$, in GCs than in the Galactic field ($q=-0.48 \pm 0.04$), suggesting a different formation or evolutionary history of MSPs in these two regions of the Galaxy. To probe the production mechanism of MSPs in clusters, a search of the possible relationships between the MSP population and cluster properties was carried out. The results of an investigation of 9 GCs indicate positive correlations between the MSP population and the stellar encounter rate and metallicity. This provides additional evidence suggesting that stellar dynamical interactions are important in the formation of the MSP population in GCs.

We present for the first time a coherent model of the polarized Galactic synchrotron and thermal dust emissions which are the main diffuse foreground for the measurement of the polarized power spectra of the CMB fluctuations with the Planck satellite mission. We produce 3D models of the Galactic magnetic field including regular and turbulent components, and of the distribution of matter in the Galaxy, relativistic electrons and dust grains. By integrating along the line of sight we construct maps of the polarized Galactic synchrotron and thermal dust emission for each of these models and compare them to currently available data. We consider the 408 MHz all-sky continuum survey, the 23 GHz band of the Wilkinson Microwave Anisotropy Probe and the 353 GHz Archeops data.}{The best-fit parameters obtained are consistent with previous estimates in the literature based only on synchrotron emission and pulsar rotation measurements. They allows us to reproduce the large scale structures observed on the data. Poorly understood local Galactic structures and turbulence make difficult an accurate reconstruction of the observations in the Galactic plane. Finally, using the best-fit model we are able to estimate the expected polarized foreground contamination at the Planck frequency bands. For the CMB bands, 70, 100, 143 and 217 GHz, at high Galactic latitudes although the CMB signal dominates in general, a significant foreground contribution is expected at large angular scales. In particular, this contribution will dominate the CMB signal for the B modes expected from realistic models of a background of primordial gravitational waves.

We report on submillimetre bolometer observations of the isolated neutron star RX J1856.5--3754 using the LABOCA bolometer array on the Atacama Pathfinder Experiment (APEX) Telescope. No cold dust continuum emission peak at the position of RX J1856.5--3754 was detected. The 3 sigma flux density upper limit of 5 mJy translates into a cold dust mass limit of a few earth masses. We use the new submillimetre limit, together with a previously obtained H-band limit, to constrain the presence of a gaseous, circumpulsar disc. Adopting a simple irradiated-disc model, we obtain a mass accretion limit of dM/dt less than 10^{14} g/s, and a maximum outer disc radius of around 10^{14} cm. By examining the projected proper motion of RX J1856.5--3754, we speculate about a possible encounter of the neutron star with a dense fragment of the CrA molecular cloud a few thousand years ago.

We report the detection of Cd I (Z = 48), Lu II (Z = 71), and Os II (Z = 76) in the metal-poor star BD+173248. These abundances are derived from an ultraviolet spectrum obtained with the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. This is the first detection of these neutron-capture species in a metal-poor star enriched by the r-process. We supplement these measurements with new abundances of Mo I, Ru I, and Rh I derived from an optical spectrum obtained with the High Resolution Echelle Spectrograph on Keck. Combined with previous abundance derivations, 32 neutron-capture elements have been detected in BD+173248, the most complete neutron-capture abundance pattern in any metal-poor star to date. The light neutron-capture elements (38 <= Z <= 48) show a more pronounced even-odd effect than expected from current Solar system r-process abundance predictions. The age for BD+173248 derived from the Th II/Os II chronometer is in better agreement with the age derived from other chronometers than the age derived from Th II/Os I. New Hf II abundance derivations from transitions in the ultraviolet are lower than those derived from transitions in the optical, and the lower Hf abundance is in better agreement with the scaled Solar system r-process distribution.

We have developed a method to compute the possible distribution of radio emission regions in a typical pulsar magnetosphere, taking into account the viewing geometry and rotational effects of the neutron star. Our method can estimate the emission altitude and the radius of curvature of particle trajectory as a function of rotation phase for a given inclination angle, impact angle, spin-period, Lorentz factor, field line constant and the observation frequency. Further, using curvature radiation as the basic emission mechanism, we simulate the radio intensity profiles that would be observed from a given distribution of emission regions, for different values of radio frequency and Lorentz factor. We show clearly that rotation effects can introduce significant asymmetries into the observed radio profiles. We investigate the dependency of profile features on various pulsar parameters. We find that the radiation from a given ring of field lines can be seen over a large range of pulse longitudes, originating at different altitudes, with varying spectral intensity. Preferred heights of emission along discrete sets of field lines are required to reproduce realistic pulsar profiles, and we illustrate this for a known pulsar. Finally, we show how our model provides feasible explanations for the origin of core emission, and also for one-sided cones which have been observed in some pulsars.

We study the mean profiles of the multi--component pulsars PSRs B1839+09, B1916+14 and B2111+46. We estimate the emission height of the core components, and hence find the absolute emission altitudes corresponding to the conal components. By fitting Gaussians to the emission components, we determine the phase location of the component peaks. Our findings indicate that the emission beams of these pulsars have the nested core--cone structures. Based on the phase location of the component peaks, we estimate the aberration--retardation (A/R) phase shifts in the profiles. Due to the A/R phase shift, the peak of the core component in the intensity profile and the inflection point of the polarization angle swing are found to be symmetrically shifted in the opposite directions with respect to the meridional plane in such a way that the core shifts towards the leading side and the polarization angle inflection point towards the trailing side. We have been able to locate the phase location of the meridional plane and to estimate the absolute emission altitude of both the core and the conal components relative to the neutron star centre, using the exact expression for the A/R phase shift given by Gangadhara (2005).

Some of X-ray dim isolated neutron stars (XDINS) and central compact objects in supernova remnants (CCO) show absorption features in their thermal soft X-ray spectra. There have been suggestions in the literature that these features could be due to the periodic peaks in free-free absorption opacities, caused either by Landau quantization of electron motion in magnetic fields B<10^{11} G or analogous quantization of ion motion in magnetic fields B>10^{13} G. I review the physics behind cyclotron quantum harmonics in free-free photoabsorption, discuss different approximations for their calculation, and explain why the ion cyclotron harmonics (beyond the fundamental) cannot be observed.

Millisecond pulsars are rapidly rotating neutron stars where general relativity plays a strong role in the propagation of light from the neutron star to observer. The observed X-ray pulse shapes carry information on the mass, radius and surface shape of the neutron star. Comparison of theoretical calculations of pulse shapes with observed pulse shapes can give useful constraints on neutron star properties. Then comparison with calculated properties giving an assumed equation of state (EOS) can confirm or rule out the assumed EOS.

We study the effect of the non-linear process of ambipolar diffusion (joint transport of magnetic flux and charged particles relative to neutral particles) on the long-term behavior of a non-uniform magnetic field in a one-dimensional geometry. Our main focus is the dissipation of magnetic energy inside neutron stars(particularly magnetars), but our results have a wider application, particularly to the interstellar medium and the loss of magnetic flux from collapsing molecular cloud cores. Our system is a weakly ionized plasma in which neutral and charged particles can be converted into each other through nuclear beta decays (or ionization-recombination processes). In the "weak-coupling" limit of infrequent inter-particle interactions, the evolution of the magnetic field is controlled by the beta decay rate and can be described by a non-linear partial integro-differential equation. In the opposite, "strong-coupling" regime, the evolution is controlled by the inter-particle collisions and can be modelled through a non-linear diffusion equation. We show numerically that, in both regimes, ambipolar diffusion tends to spread out the magnetic flux, but, contrary to the normal Ohmic diffusion, it produces sharp magnetic field gradients with associated current sheets around those regions where the magnetic field is weak. In the case of the non-linear diffusion equation, this leads to magnetic reconnection even in the absence of resistivity.

We report the first detection of a linear correlation between rms variability amplitude and flux in the Ultraluminous X-ray source NGC 5408 X-1. The rms-flux relation has previously been observed in several Galactic black hole X-ray binaries (BHBs), several Active Galactic Nuclei (AGN) and at least one neutron star X-ray binary. This result supports the hypothesis that a linear rms-flux relation is common to all luminous black hole accretion and perhaps even a fundamental property of accretion flows about compact objects. We also show for the first time the cross-spectral properties of the variability of this ULX, comparing variations below and above 1 keV. The coherence and time delays are poorly constrained but consistent with high coherence between the two bands, over most of the observable frequency range, and a significant time delay (with hard leading soft variations). The magnitude and frequency dependence of the lags are broadly consistent with those commonly observed in BHBs, but the direction of the lag is reversed. These results indicate that ULX variability studies, using long X-ray observations, hold great promise for constraining the processes driving ULXs behaviour, and the position of ULXs in the scheme of black hole accretion from BHBs to AGN.

In order to establish whether the unstable r-modes in a rotating neutron star provide a detectable source of gravitational waves, we need to understand the details of the many dissipative processes that tend to counteract the instability. It has been established that the bulk viscosity due to exotic particles, like hyperons, may be particularly important in this respect. However, the effects of hyperon superfluidity have so far not been fully accounted for. While the associated suppression of the reaction rates that give rise to the bulk viscosity has been estimated, superfluid aspects of the fluid dynamics have not been considered. In this paper we determine the r-mode instability window for a neutron star with a $\Sigma^{-}$ hyperon core, using the appropriate multifluid formalism including, for the first time, the effect of the "superfluid" bulk viscosity coefficients. We demonstrate that, even though the extra terms may increase the bulk viscosity damping somewhat, their presence does not affect the qualitative features of the r-mode instability window.