Exact relativistic force free fields with cylindrical symmetry are explored. Such fields are generated in the interstellar gas via their connection to pulsar magnetospheres both inside and outside their light cylinders. The possibility of much enhanced interstellar fields wound on cylinders of Solar system dimensions is discussed but these are most likely unstable.

Our goal is to understand the specificities of highly absorbed sgHMXB and in particular of the companion stellar wind, thought to be responsible for the strong absorption. We have monitored IGR J17252-3616, a highly absorbed system featuring eclipses, with XMM-Newton to study the vari- ability of the column density and of the Fe K{\alpha} emission line along the orbit and during the eclipses. We also built a 3D model of the structure of the stellar wind to reproduce the observed variability. We first derived a refined orbital solution built from INTEGRAL, RXTE and XMM data. The XMM monitoring campaign revealed significant variation of intrinsic absorbing column density along the orbit and of the Fe K{\alpha} line equivalent width around the eclipses. The origin of the soft X-ray absorption is modeled with an dense and extended hydrodynamical tail, trailing the neutron star. This structure extends along most of the orbit, indicating that the stellar wind is strongly disrupted by the neutron star. The variability of the absorbing column density suggests that the terminal velocity of the wind is smaller (~400 km/s) than observed in classical systems. This can also explain the much stronger density perturbation inferred from the observations. Most of the Fe K{\alpha} emission is generated in the most inner region of the hydrodynamical tail. This region, that extends over a few accretion radii, is ionized and does not contribute to the soft X-ray absorption. We have built a qualitative model of the stellar wind of IGR J17252-3616 that can represent the observations and suggest that highly absorbed systems have a lower wind velocity than classical sgHMXB. This proposal could be tested with de- tailed numerical simulations and high-resolution infrared/optical observations. If confirmed, it may turn out that half of the persistent sgHMXB have low stellar wind speeds.

We consider isolated compact remnants (ICoRs), i.e. neutrons stars and black holes that do not reside in binary systems and therefore cannot be detected as X-ray binaries. ICoRs may represent $\sim\,5$ percent of the stellar mass budget of the Galaxy, but they are very hard to detect. Here we explore the possibility of using microlensing to identify ICoRs. In a previous paper we described a simulation of neutron star evolution in phase space in the Galaxy, taking into account the distribution of the progenitors and the kick at formation. Here we first reconsider the evolution and distribution of neutron stars and black holes adding a bulge component. From the new distributions we calculate the microlensing optical depth, event rate and distribution of event time scales, comparing and contrasting the case of ICoRs and "normal stars". We find that the contribution of remnants to optical depth is slightly lower than without kinematics, owing to the evaporation from the Galaxy. On the other hand, the relative contribution to the rate of events is a factor $\sim\,5$ higher. In all, $\sim\,6-7$ percent of the events are likely related to ICoRs. In particular, $\sim\,30-40$ percent of the events with duration $>\,100$ days are possibly related to black holes. It seems therefore that microlensing observations are a suitable tool to probe the population of Galactic ICoRs.

We present the results of new Agile observations of PSR B1509-58 performed over a period of 2.5 years following the detection obtained with a subset of the present data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier by Comptel up to 30 MeV: a broad asymmetric first peak reaching its maximum 0.39 +/- 0.02 cycles after the radio peak plus a second peak at 0.94 +/- 0.03. The gamma-ray spectral energy distribution of the pulsed flux detected by Comptel and Agile is well described by a power-law (photon index alpha=1.87+/-0.09) with a remarkable cutoff at E_c=81 +/- 20 MeV, representing the softest spectrum observed among gamma-ray pulsars so far. The pulsar luminosity at E > 1 MeV is $L_{\gamma}=4.2^{+0.5}_{-0.2} \times10^{35}$ erg/s, assuming a distance of 5.2 kpc, which implies a spin-down conversion efficiency to gamma-rays of $\sim 0.03$. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In this interpretation our spectrum constrains the magnetic altitude of the emission point(s) at 3 km above the neutron star surface, implying that the attenuation may not be as strong as formerly suggested because pair production can substitute photon splitting in regions of the magnetosphere where the magnetic field becomes too low to sustain photon splitting. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce Agile lightcurves and spectra match the polarization measurements.

Gamma-ray burst X-ray flares are believed to mark the late time activity of the central engine. We compute the temporal evolution of the average flare luminosity $< L >$ in the common rest frame energy band of 44 GRBs taken from the large \emph{Swift} 5-years data base. Our work highlights the importance of a proper consideration of the threshold of detection of flares against the contemporaneous continuous X-ray emission. In the time interval $30 \rm{s}<t<1000\,\rm{s}$ we find $< L >\propto t^{-2.7\pm 0.1}$; this implies that the flare isotropic energy scaling is $E_{\rm{iso,flare}}\propto t^{-1.7}$. The decay of the continuum underlying the flare emission closely tracks the average flare luminosity evolution, with a typical flare to steep-decay luminosity ratio which is $L_{\rm{flare}}/L_{\rm{steep}}=4.7$: this suggests that flares and continuum emission are deeply related to one another. We infer on the progenitor properties considering different models. According to the hyper-accreting black hole scenario, the average flare luminosity scaling can be obtained in the case of rapid accretion ($t_{\rm{acc}}\ll t$) or when the last $\sim 0.5 M_{\sun}$ of the original $14 M_{\sun}$ progenitor star are accreted. Alternatively, the steep $\propto t^{-2.7}$ behaviour could be triggered by a rapid outward expansion of an accretion shock in the material feeding a convective disk. If instead we assume the engine to be a rapidly spinning magnetar, then its rotational energy can be extracted to power a jet whose luminosity is likely to be between the monopole ($L\propto e^{-2t}$) and dipole ($L\propto t^{-2}$) cases. In both scenarios we suggest the variability, which is the main signature of the flaring activity, to be established as a consequence of different kinds of instabilities.

A design study is currently in progress for a third generation gravitational-wave (GW) detector called Einstein Telescope (ET). An important kind of source for ET will be the inspiral and merger of binary neutron stars (BNS) up to $z \sim 2$. If BNS mergers are the progenitors of short-hard $\gamma$-ray bursts, then some fraction of them will be seen both electromagnetically and through GW, so that the luminosity distance and the redshift of the source can be determined separately. An important property of these `standard sirens' is that they are \emph{self-calibrating}: the luminosity distance can be inferred directly from the GW signal, with no need for a cosmic distance ladder. Thus, standard sirens will provide a powerful independent check of the $\Lambda$CDM model. In previous work, estimates were made of how well ET would be able to measure a subset of the cosmological parameters (such as the dark energy parameter $w_0$) it will have access to, assuming that the others had been determined to great accuracy by alternative means. Here we perform a more careful analysis by explicitly using the potential Planck CMB data as prior information for these other parameters. We find that ET will be able to constrain $w_0$ and $w_a$ with accuracies $\Delta w_0 = 0.096$ and $\Delta w_a = 0.296$, respectively. These results are compared with projected accuracies for the JDEM Baryon Acoustic Oscillations (BAO) project and the SNAP Type Ia supernovae (SNIa) observations. Comparing with the combination of the future CMB(Planck)+BAO(JDEM)+SNIa(SNAP) projects, the contribution of GW standard sirens can decrease the uncertainties on $w_0$ and $w_a$ by $\sim 6%$.

We report on a multi-wavelength study of the compact object candidate 1RXS J141256.0+792204 (Calvera). Calvera was observed in the X-rays with XMM/EPIC twice for a total exposure time of ~50 ks. The source spectrum is thermal and well reproduced by a two component model composed of either two hydrogen atmosphere models, or two blackbodies (kT_1~ 55/150 eV, kT_2~ 80/250 eV, respectively, as measured at infinity). Evidence was found for an absorption feature at ~0.65 keV; no power-law high-energy tail is statistically required. Using pn and MOS data we discovered pulsations in the X-ray emission at a period P=59.2 ms. The detection is highly significant (> 11 sigma), and unambiguously confirms the neutron star nature of Calvera. The pulse profile is nearly sinusoidal, with a pulsed fraction of ~18%. We looked for the timing signature of Calvera in the Fermi Large Area Telescope (LAT) database and found a significant (~5 sigma) pulsed signal at a period coincident with the X-ray value. The gamma-ray timing analysis yielded a tight upper limit on the period derivative, dP/dt < 5E-18 s/s (dE_rot/dt <1E33 erg/s, B<5E10 G for magneto- dipolar spin-down). Radio searches at 1.36 GHz with the 100-m Effelsberg radio telescope yielded negative results, with a deep upper limit on the pulsed flux of 0.05 mJy. Diffuse, soft (< 1 keV) X-ray emission about 13' west of the Calvera position is present both in our pointed observations and in archive ROSAT all-sky survey images, but is unlikely associated with the X-ray pulsar. Its spectrum is compatible with an old supernova remnant (SNR); no evidence for diffuse emission in the radio and optical bands was found. The most likely interpretations are that Calvera is either a central compact object escaped from a SNR or a mildly recycled pulsar; in both cases the source would be the first ever member of the class detected at gamma-ray energies.

The stellar origin of gamma-ray bursts can be explained by the rapid release of energy in a highly collimated, extremely relativistic jet. This in turn appears to require a rapidly spinning highly magnetised stellar core that collapses into a magnetic neutron star or a black hole within a relatively massive envelope. They appear to be associated with type Ib/c supernovae but, with a birthrate of around 10^{-6}-10^{-5} per year per galaxy, they are considerably rarer than such supernovae in general. To satisfy all these requirements we hypothesize a binary star model that ends with the merging of an oxygen neon white dwarf with the carbon-oxygen core of a naked helium star during a common envelope phase of evolution. The rapid spin and high magnetic field are natural consequences of such a merging. The evolution that leads to these progenitors is convoluted and so naturally occurs only very rarely. To test the hypothesis we evolve a population of progenitors and find that the rate is as required. At low metallicity we calculate that a similar fraction of stars evolve to this point and so would expect the gamma-ray burst rate to correlate with the star formation rate in any galaxy. This too is consistent with observations. These progenitors, being of intermediate mass, differ radically from the usually postulated high-mass stars. Thus we can reconcile observations that the bursts occur close to but not within massive star associations.

The pulsar current, in the $P$--$ {\dot P}$ plane where $P$ is the pulsar period and ${\dot P}$ is the period derivative, is used as a supposedly ``model free'' way to estimate the pulsar birthrate from statistical data on pulsars. We reconsider the derivation of the kinetic equation on which this is based, and argue that the interpretation of the pulsar current is strongly model dependent, being sensitive to the form of the assumed evolution law for pulsars. We discuss the case where the trajectory of a pulsar is assumed to be of the form ${\dot P}=KP^{2-n}$ with $K$ and $n$ constant, and show that (except for $n=2$) one needs to introduce a pseudo source term in order to infer the birthrate from the pulsar current. We illustrate the effect of this pseudo source term using pulsar data to estimate the birthrate for different choices of $n$. We define and discuss an alternative ``potential'' class of evolution laws for which this complication is avoided due to the pseudo source term being identically zero.

We report the results of a Suzaku observation of the anomalous X-ray pulsar (AXP) 1E 1841-045 at a center of the supernova remnant Kes 73. We confirmed that the energy-dependent spectral models obtained by the previous separate observations were also satisfied over a wide energy range from 0.4 to ~70 keV, simultaneously. Here, the models below ~10 keV were a combination of blackbody (BB) and power-law (PL) functions or of two BBs wit h different temperatures at 0.6 - 7.0 keV (Morii et al. 2003), and that above ~20 keV was a PL function (Kuiper Hermsen Mendez 2004). The combination BB + PL + PL was found to best represent the phase-averaged spectrum. Phase-resolved spectroscopy indicated the existence of two emission regions, one with a thermal and the other with a non-thermal nature. The combination BB + BB + PL was also found to represent the phase-averaged spectrum well. However, we found that this model is physically unacceptable due to an excessively large area of the emission region of the blackbody. Nonetheless, we found that the temperatures and radii of the two blackbody components showed moderate correlations in the phase-resolved spectra. The fact that the same correlations have been observed between the phase-averaged spectra of various magnetars (Nakagawa et al. 2009) suggests that a self-similar function can approximate the intrinsic energy spectra of magnetars below ~10 keV.

We present the Second Palermo Swift-BAT hard X-ray catalogue obtained by analysing data acquired in the first 54 months of the Swift mission. Using our software dedicated to the analysis of data from coded mask telescopes, we analysed the BAT survey data in three energy bands (15-30 keV, 15-70 keV, 15-150 keV), obtaining a list of 1256 detections above a significance threshold of 4.8 standard deviations. The identification of the source counterparts is pursued using two strategies: the analysis of field observations of soft X-ray instruments and cross-correlation of our catalogue with source databases.The survey covers 50% of the sky to a 15--150 keV flux limit of 1.0 x 10^-11 erg s^-1 cm^-2 and 9.2 x 10^-12 erg s^-1 cm^-2 for |b|< 10 degrees and |b|> 10 degrees, respectively. The Second Palermo Swift-BAT hard X-ray catalogue includes 1079 (86%) hard X-ray sources with an associated counterpart (26 with a double association and 2 with a triple association) and 177 BAT excesses (14%) that still lack a counterpart. The distribution of the BAT sources among the different object classes consists of 19% Galactic sources, 57% extragalactic sources, and 10% sources with a counterpart at softer energies whose nature has not yet been determined. About half of the BAT associated sources lack a counterpart in the ROSAT catalogues. This suggests that either moderate or strong absorption may be preventing their detection in the ROSAT energy band. The comparison of our BAT catalogue with the Fermi Large Area Telescope First Source Catalogue identifies 59 BAT/Fermi correspondences: 48 blazars, 3 Seyfert galaxies, 1 interacting galaxy, 3 high mass X-ray binaries, and 4 pulsars/supernova remnants. This small number of correspondences indicates that different populations make the sky shine in these two different energy bands.

It has recently been suggested that neutron stars inside the shells of young supernova remnants (SNR) are the sources of PeV cosmic rays and that the interaction of the particles with the radiation field in the SNR causes electron pair production, which has relevance to recent observations of 'high' positron fluxes. Furthermore, the character of the interaction is such that the well-known knee in the cosmic ray energy spectrum can be explained. Our examination of the mechanism leads us to believe that the required parameters of SN and pulses are so uncommon that the knee and positron fraction can only be explained if a single, local and recent SN - and associated pulsar - are concerned.

Progenitors of long GRBs, and core-collapse supernovae in general, may have two separate mechanisms driving the outflows: quasi-isotropic neutrino-driven supernova explosions followed by a highly collimated relativistic outflow driven by the GRB central engine, a black hole or a magnetar. We consider the dynamics of the second GRB-driven explosion propagating through expanding envelope generated by the passage of the primary supernova shock. Beyond the central core, in the region of steep density gradient created by the SN shock breakout, the accelerating secondary quasi-spherical GRB shock become unstable to corrugation and under certain conditions may form a highly collimated jet, a "chimney", when a flow expands almost exclusively along a nearly cylindrically collimated channel. Thus, weakly non-spherical driving and/or non-spherical initial conditions of the wind cavity may produce highly non-spherical, jetted outflows. For a constant luminosity GRB central engine, this occurs for density gradient in the envelope \rho ~ r^{-\omega} steeper than \omega >4.

Our high-time-resolution observations reveal that individual main pulses from the Crab pulsar contain one or more short-lived microbursts. Both the energy and duration of bursts measured above 1 GHz can vary dramatically in less than a millisecond. These fluctuations are too rapid to be caused by propagation through turbulence in the Crab Nebula or the interstellar medium; they must be intrinsic to the radio emission process in the pulsar. The mean duration of a burst varies with frequency as $\nu^{-2}$, significantly different from the broadening caused by interstellar scattering. We compare the properties of the bursts to some simple models of microstructure in the radio emission region.

Methods are developed for constructing spectral representations of cold (barotropic) neutron-star equations of state. These representations are faithful in the sense that every physical equation of state has a representation of this type, and conversely every such representation satisfies the minimal thermodynamic stability criteria required of any physical equation of state. These spectral representations are also efficient, in the sense that only a few spectral coefficients are generally required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency is illustrated by constructing spectral fits to a large collection of "realistic" neutron-star equations of state.

Pulsars are rapidly-rotating, highly-magnetized neutron stars emitting radiation across the electromagnetic spectrum. Although there are more than 1800 known radio pulsars, until recently, only seven were observed to pulse in gamma rays and these were all discovered at other wavelengths. The Fermi Large Area Telescope makes it possible to pinpoint neutron stars through their gamma-ray pulsations. We report the detection of 16 gamma-ray pulsars in blind frequency searches using the LAT. Most of these pulsars are coincident with previously unidentified gamma-ray sources, and many are associated with supernova remnants. Direct detection of gamma-ray pulsars enables studies of emission mechanisms, population statistics and the energetics of pulsar wind nebulae and supernova remnants.

We examine electron-capture supernovae (ECSNe) as sources of elements heavier than iron in the solar system and in Galactic halo stars. Nucleosynthesis calculations are performed on the basis of thermodynamic histories of mass elements from a fully self-consistent, two-dimensional (2D) hydrodynamic explosion model of an ECSN. We find that neutron-rich convective lumps with an electron fraction down to Ye,min=0.40, which are absent in the one-dimensional (1D) counterpart, allow for interesting production of elements between the iron group and N=50 nuclei (from Zn to Zr, with little Ga) in nuclear statistical equilibrium and by the alpha-process. Our models yield very good agreement with the Ge, Sr, Y, and Zr abundances of r-process deficient Galactic halo stars and constrain the occurrence of ECSNe to ~4% of all stellar core-collapse events. If tiny amounts of additional material with slightly lower Ye,min down to ~0.30-0.35 were also ejected - which presently cannot be excluded because of the limitations of resolution and two-dimensionality of the model -, a weak r-process can yield elements beyond N=50 up to Pd, Ag, and Cd as observed in the r-process deficient stars.

Neutron stars are the densest objects known in the Universe. Being the final product of stellar evolution, their internal composition and structure is rather poorly constrained by measurements.
It is the purpose of this paper to put some constrains on the mass and moment of inertia of neutron stars based on the interpretation of kHz quasi-periodic oscillations observed in low mass X-ray binaries.
We use observations of high-frequency quasi-periodic observations (HF-QPOs) in low mass X-ray binaries (LMXBs) to look for the average mass and moment of inertia of neutron stars. This is done by applying our parametric resonance model to discriminate between slow and fast rotators.
We fit our model to data from ten LMXBs for which HF-QPOs have been seen and the spin of the enclosed accreting neutron star is known. For a simplified analysis we assume that all neutron stars possess the same properties (same mass $M_*$ and same moment of inertia $I_*$). We find an average mass $M_* \approx 2.0-2.2\, M_{\odot}$. The corresponding average moment of inertia is then $I_* \approx 1-3 \times 10^{38}\;{\rm kg\,m^2} \approx 0.5-1.5 \, (10\;\textrm{ km})^2 \, M_\odot$ which equals to dimensionless spin parameter $\tilde{a} \approx 0.05-0.15$ for slow rotators (neutron stars with a spin frequency roughly about 300~Hz) respectively $\tilde{a} \approx 0.1-0.3$ for fast rotators (neutron stars with the spin frequency roughly about 600~Hz).

Most of the spectra of neutron star low mass X-ray binaries (NS LMXBs), being them persistent or transient, are characterized by the presence of a strong thermal Comptonization bump, thought to originate in the transition layer (TL) between the accretion disk and the NS surface. The observable quantities which characterize this component dominating the emission below 30 keV, are the spectral index alpha and the rollover energy, both related to the electron temperature and optical depth of the plasma. Starting from observational results on a sample of NS LMXBs in different spectral states, we formulate the problem of X-ray spectral formation in the TL of these sources. We predict a stability of the thermal Comptonization spectral index in different spectral states if the energy release in the TL is much higher than the intercepted flux coming from the accretion disk. We use an equation for the energy balance and the radiative transfer diffusion equation for a slab geometry in the TL, to derive a formula for the thermal Comptonization index alpha. We show that in this approximation the TL electron temperature kTe and optical depth tau_0 can be written as a function of the energy flux from the disk intercepted by the corona (TL) and that in the corona itself Qdisk/Qcor, in turn leading to a relation alpha=f(Qdisk/Qcor), with alpha ~ 1 when Qdisk/Qcor <<1. We show that the observed spectral index alpha for the sample of sources here considered lies in a belt around 1 +/- 0.2 a part for the case of GX 354--0. Comparing our theoretical predictions with observations, we claim that this result, which is consistent with the condition Qdisk/Qcor <<1, can give us constraints on the accretion geometry of these systems, an issue that seems difficult to be solved using only the spectral analysis method.

We suggest that white dwarf (WD) pulsars can compete with neutron star (NS) pulsars for producing the excesses of cosmic ray electrons/positrons observed by the PAMELA, ATIC/PPB-BETS, Fermi and HESS experiments. A merger of two WDs leads to a rapidly spinning WD with a rotational energy comparable to the NS case. The birth rate is also similar, providing the right energy budget for the cosmic ray electrons/positrons. Applying the NS theory, we suggest that the WD pulsars can in principle produce electrons/positrons up to 10 TeV. In contrast to the NS model, the adiabatic and radiative energy losses of electrons/positrons are negligible since their injection continues after the expansion of the pulsar wind nebula, and hence it is enough that a fraction 1% of WDs are magnetized as observed. The long activity also increases the number of nearby sources, which reduces the Poisson fluctuation in the flux. The WD pulsars could dominate the quickly cooling electrons/positrons above TeV energy as a second spectral bump or even surpass the NS pulsars in the observing energy range 10 GeV - 1 TeV, providing a background for the dark matter signals and a nice target for the future AMS-02, CALET and CTA experiment.

One of the most essential but uncertain processes for producing cosmic-rays (CRs) and their spectra is how accelerated particles escape into the interstellar space. We propose that the CR electron spectra at >~TeV energy can provide more direct and powerful probe of the CR escape than the CR nuclei and gamma-rays. We calculate the electron spectra from a young pulsar embedded in the supernova remnant (SNR), like Vela, taking into account the energy-dependent CR escape. Since SNRs would accelerate and hence confine particles with energy up to ~10^{15.5}eV, only energetic particles can escape first, while the lower energy particles are confined and released later. Then the observed electron spectrum should have a low energy cutoff whose position marks the age of the pulsar/SNR. The low energy cutoff is observable in the >~ TeV energy window, where other contaminating sources are expected to be small due to the fast cooling of electrons. The spectrum looks similar to a dark matter annihilation line if the low energy cutoff is close to the high energy intrinsic or cooling break. The future experiments such as CALET and CTA are capable of directly detecting the CR escape features toward revealing the origin of CRs.

The Galactic Black hole candidate XTE J1752-223 was observed during the decay of its 2009 outburst with the Suzaku and XMM-Newton observatories. The observed spectra are consistent with the source being in the ''intermediate`` and ''low-hard state`` respectively. The presence of a strong, relativistic iron emission line is clearly detected in both observations and the line profiles are found to be remarkably consistent and robust to a variety of continuum models. This strongly points to the compact object in \j\ being a stellar-mass black hole accretor and not a neutron star. Physically-motivated and self-consistent reflection models for the Fe-\ka\ emission-line profile and disk reflection spectrum rule out either a non-rotating, Schwarzchild black hole or a maximally rotating, Kerr black hole at greater than 3sigma level of confidence. Using a fully relativistic line function in which the black hole spin parameter is a variable, we have formally constrained the spin parameter to be $0.52\pm0.11 (1\sigma)$. Furthermore, we show that the source in the low--hard state still requires an optically--thick disk component having a luminosity which is consistent with the $L\propto T^4$ relation expected for a thin disk extending down to the inner--most stable circular orbit. Our result is in contrast to the prevailing paradigm that the disk is truncated in the low-hard state.

Two faint X-ray pulsars, AX J1749.2-2725 and AX J1749.1-2733, located in the direction to the Galactic Center, were studied in detail using data of INTEGRAL, XMM-Newton and Chandra observatories in X-rays, the SOFI/NTT instrument in infrared and the RTT150 telescope in optics. X-ray positions of both sources were determined with the uncertainty better than ~1 arcsec, that allowed us to identify their infrared counterparts. From the subsequent analysis of infrared and optical data we conclude that counterparts of both pulsars are likely massive stars of B0-B3 classes located behind the Galactic Center at distances of 12-20 kpc, depending on the type, probably in further parts of galactic spiral arms. In addition, we investigated the extinction law towards the galactic bulge and found that it is significantly different from standard one.

Two classes of high energy sources in our galaxy are believed to host magnetars, neutron stars whose emission results from the dissipation of their magnetic field. The extremely high magnetic field of magnetars distorts their shape, and causes the emission of a conspicuous gravitational waves signal if rotation is fast and takes place around a different axis than the symmetry axis of the magnetic distortion. Based on a numerical model of the cosmic star formation history, we derive the cosmological background of gravitational waves produced by magnetars, when they are very young and fast spinning. We adopt different models for the configuration and strength of the internal magnetic field (which determines the distortion) as well as different values of the external dipole field strength (which governs the spin evolution of magnetars over a wide range of parameters). We find that the expected gravitational wave background differs considerably from one model to another. The strongest signals are generated for magnetars with very intense toroidal internal fields ($\sim 10^{16}$ G range) and external dipole fields of $\sim 10^{14}$, as envisaged in models aimed at explaining the properties of the Dec 2004 giant flare from SGR 1806-20. Such signals should be easily detectable with third generation ground based interferometers such as the Einstein Telescope.

Until 2004, Anomalous X-ray Pulsars (AXPs) were known as strong emitters of soft X-rays only (< 10 keV). The discovery of hard X-ray component from AXPs provided important insight about their emission properties while it posed a serious challenge to explain its origin. The physical mechanism of the hard emission component has still not been fully resolved. We investigate the high-energy gamma-ray properties of the brightest AXP, 4U 0142+61 using data collected with the Large Area Telescope (LAT) on board Fermi Gamma-ray Space Telescope to establish the spectral behavior of the source on a very broad energy span and search for pulsed emission. Here, we present our results of detailed search for the persistent and pulsed high-energy gamma-ray emission from 4U 0142+61 which result in no significant detection. However, we obtain upper limits to the persistent high-energy gamma-ray emission flux which helps us to constrain existing physical models.

Motivated by suggestions that binaries with almost equal-mass components ("twins") play an important role in the formation of double neutron stars and may be rather abundant among binaries, we study the stability of synchronized close and contact binaries with identical components in circular orbits. In particular, we investigate the dependency of the innermost stable circular orbit on the core mass, and we study the coalescence of the binary that occurs at smaller separations. For twin binaries composed of convective main-sequence stars, subgiants, or giants with low mass cores (M_c <~0.15M, where M is the mass of a component), a secular instability is reached during the contact phase, accompanied by a dynamical mass transfer instability at the same or at a slightly smaller orbital separation. Binaries that come inside this instability limit transfer mass gradually from one component to the other and then coalesce quickly as mass is lost through the outer Lagrangian points. For twin giant binaries with moderate to massive cores (M_c >~0.15M), we find that stable contact configurations exist at all separations down to the Roche limit, when mass shedding through the outer Lagrangian points triggers a coalescence of the envelopes and leaves the cores orbiting in a central tight binary. We discuss the implications of our results to the formation of binary neutron stars.

Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (this http URL).

Hyperaccreting disks around neutron stars or magnetars cooled via neutrino emission can be the potential central engine of GRBs. The neutron-star disk can cool more efficiently, produce much higher neutrino luminosity and neutrino annihilation luminosity than its black hole counterpart with the same accretion rate. The neutron star surface boundary layer could increase the annihilation luminosity as well. An ultra relativistic jet via neutrino annihilation can be produced along the stellar poles. Moreover, we investigate the effects of strong fields on the disks around magnetars. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity; the neutrino annihilation mechanism and the magnetically-driven pulsar wind which extracts the stellar rotational energy can work together to generate and feed an even stronger ultra-relativistic jet along the stellar magnetic poles.

We present the results of new AGILE observations of PSR B1509-58 performed over a period of $\sim$2.5 years following the detection obtained with preliminary data. The modulation significance of the lightcurve above 30 MeV is at a 5$\sigma$ confidence level and the lightcurve is similar to those found earlier up to 30 MeV by COMPTEL: a broad asymmetric first peak reaching its maximum $0.39 \pm 0.02$ cycles after the radio peak plus a second peak at $0.94 \pm 0.03$. The $\gamma$-ray spectral energy distribution of pulsed flux is well described by a power-law (photon index $\alpha=1.87\pm0.09$) with a remarkable cutoff at $E_c=81\pm 20$ MeV, representing the softest spectrum observed among $\gamma$-ray pulsars so far. The unusual soft break in the spectrum of PSR B1509-58 has been interpreted in the framework of polar cap models as a signature of the exotic photon splitting process in the strong magnetic field of this pulsar. In the case of an outer-gap scenario, or the two pole caustic model, better constraints on the geometry of the emission would be needed from the radio band in order to establish whether the conditions required by the models to reproduce AGILE lightcurves and spectra match the polarization measurements.

Globular clusters (GCs) with their ages of the order of several billion years contain many final products of evolution of stars such as: neutron stars, white dwarfs and probably also black holes. These compact objects can be at present responsible for the acceleration of particles to relativistic energies. Therefore, gamma-ray emission is expected from GCs as a result of radiation processes occurring either in the inner magnetosperes of millisecond pulsars or in the vicinity of accreting neutron stars and white dwarfs or as a result of interaction of particles leaving the compact objects with the strong radiation field within the GC. Recently, GeV gamma-ray emission has been detected from several GCs by the new satellite observatory Fermi. Also Cherenkov telescopes reported interesting upper limits at the TeV energies which start to constrain the content of GCs. We review the results of these gamma-ray observations in the context of recent scenarios for their origin.

This contribution to the proceedings of "A New Golden Age for Radio Astronomy" is simply intended to give some of the highlights from pulsar observations with LOFAR at the time of its official opening: June 12th, 2010. These observations illustrate that, though LOFAR is still under construction and astronomical commissioning, it is already starting to deliver on its promise to revolutionize radio astronomy in the low-frequency regime. These observations also demonstrate how LOFAR has many "next-generation" capabilities, such as wide-field multi-beaming, that will be vital to open a new Golden Age in radio astronomy through the Square Kilometer Array and its precursors.

Radio astronomy has benefited greatly from advances in technology and will continue to do so in the future. In fact, we are experiencing a revolution in the way radio astronomy is conducted as our instruments allow us now to directly "digitize" our photons. This has enormous consequences, since we can greatly benefit from the continuing advances in digital electronics, telecommunication and computing. The results are dramatic increase in observable bandwidths, FoVs, frequency coverage and collecting area. The global efforts will culminate in the construction of the SKA as the world's largest and most powerful telescope. On the way projects like LOFAR, LEAP and others will revolutionize many areas of astrophysics and fundamental physics. Observations of pulsars will play a central role in these scientific endeavours. We briefly summarize here some recent scientific developments that help us in defining our expectations for the the new generation of radio telescopes for pulsar astrophysics.

Gamma-ray binaries are systems containing a massive star and a compact object that have been detected up to TeV energies. The high energy emission could result from particle acceleration in the region where the stellar wind from the massive star interacts with the relativistic wind from a young pulsar. LS 5039 has the most compact orbit amongst gamma-ray binaries and its X-ray lightcurve shows a stable modulation synchronized with the orbital period. Photoelectric absorption of X-rays in the O star wind and occultation of the X-ray emitting region by the massive star can alter the X-ray lightcurve and spectrum along the orbit. Yet, the X-ray spectrum and lightcurve of LS 5039 do not show intrinsic absorption or X-ray eclipses. We study these effects in the framework of the pulsar wind scenario as a function of the binary inclination angle, the stellar wind mass-loss rate and the size of the X-ray emitter. An extended X-ray emission region >~ 3R_star appears necessary to reconcile the pulsar wind scenario with observations.

In this work, we perform the detailed analysis of absorption features in spectra of magnetar candidates observed by XMM-Newton satellite. No significant line-like feature has been found. This negative result may indicate the possible presence of smoothing out the absorption features mechanisms.

LS I +61 303 is an exceptionally rare example of a high mass X-ray binary (HMXB) that also exhibits MeV-TeV emission, making it one of only a handful of "gamma-ray binaries". Here we present H-alpha spectra that show strong variability during the 26.5 day orbital period and over decadal time scales. We detect evidence of a spiral density wave in the Be circumstellar disk over part of the orbit. The H-alpha line profile also exhibits a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the line profile, as the compact source moves almost directly away from the observer. We investigate several possible origins for this red shoulder, including an accretion disk, mass transfer stream, and a compact pulsar wind nebula that forms via a shock between the Be star's wind and the relativistic pulsar wind.

We present Chandra HETGS spectroscopy of the Be/X-ray binary 1A 0535+262 obtained during the 2009/2010 giant outburst. These are the first CCD grating spectra of this type of system during a giant outburst. Our spectra reveal a number of lines including a narrow Fe K_alpha emission line with a FWHM of ~ 5000 km s^-1. For the first time, we detect the presence of a highly ionized outflow in a Be/X-ray binary. Assuming that the line is He-like Fe XXV, fits with a simple Gaussian imply an outflow velocity of ~ 1500 km s^-1. However, self-consistent photoionization modeling with XSTAR suggests that Fe XXIII-XXIV must also contribute. In this case, an outflow velocity of ~ 3000 km s^-1 is implied. These results are discussed in the context of the accretion flow in Be-star, neutron star, and black hole X-ray binaries.

Fermi has detected gamma-ray emission from eight globular clusters. We suggest that the gamma-ray emission from globular clusters may result from the inverse Compton scattering between relativistic electrons/positrons in the pulsar wind of MSPs in the globular clusters and background soft photons including cosmic microwave/relic photons, background star lights in the clusters, the galactic infrared photons and the galactic star lights. We show that the gamma-ray spectrum from 47 Tuc can be explained equally well by upward scattering of either the relic photons, the galactic infrared photons or the galactic star lights whereas the gamma-ray spectra from other seven globular clusters are best fitted by the upward scattering of either the galactic infrared photons or the galactic star lights. We also find that the observed gamma-ray luminosity is correlated better with the combined factor of the encounter rate and the background soft photon energy density. Therefore the inverse Compton scattering may also contribute to the observed gamma-ray emission from globular clusters detected by Fermi in addition to the standard curvature radiation process. Furthermore, we find that the emission region of high energy photons from globular cluster produced by inverse Compton scattering is substantially larger than the core of globular cluster with a radius >10pc. The diffuse radio and X-rays emitted from globular clusters can also be produced by synchrotron radiation and inverse Compton scattering respectively. We suggest that future observations including radio, X-rays, and gamma-rays with energy higher than 10 GeV and better angular resolution can provide better constraints for the models.

A simple description of superfluid hydrodynamics in the inner crust of a neutron star is given. Particular attention is paid to the effect of the lattice of nuclei on the properties of the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei.

2009 has been an extraordinary year for gamma-ray pulsar astronomy and 2010 promises to be equally good. Not only have we registered an extraordinary increase in the number of pulsars detected in gamma rays, but we have also witnessed the birth of new sub-families: first of all, the radio-quiet gamma pulsars and later an ever growing number of millisecond pulsars, a real surprise. We started with a sample of 7 gamma-ray emitting neutron stars (6 radio pulsars and Geminga) and now the Fermi-LAT harvest encompasses 24 "Geminga-like" new gamma-ray pulsars, a dozen millisecond pulsars and about thirty radio pulsars. Moreover, radio searches targeted to LAT unidentified sources yielded 18 new radio millisecond pulsars, several of which have been already detected also in gamma rays. Thus, currently the family of gamma-ray emitting neutron stars seems to be evenly divided between classical radio pulsars, millisecond pulsars and radio quiet neutron stars.

We study $m=1$ oscillations and instabilities of magnetised neutron stars, by numerical time-evolution of linear perturbations of the system. The background stars are stationary equilibrium configurations with purely toroidal magnetic fields. We find that an $m=1$ instability of toroidal magnetic fields, already known from local analyses, may also be found in our relatively low-resolution global study. We present quantitative results for the instability growth rate and its suppression by rotation. The instability is discussed as a possible trigger mechanism for Soft Gamma Repeater (SGR) flares. Although our primary focus is evolutions of magnetised stars, we also consider perturbations about unmagnetised background stars in order to study $m=1$ inertial modes. We track these modes up to break-up frequency, extending known slow-rotation results.

Motivated by the recent proposal that one can obtain quasi-periodic oscillations (QPOs) by photon echoes manifesting as non-trivial features in the autocorrelation function (ACF), we study the ACFs of the light curves of three accreting black hole candidates and a neutron star already known to exhibit QPOs namely, GRS 1915+105, XTE J1550-564, XTE J1859+226 and Cygnus X-2. We compute and focus on the form of the ACFs in search of systematics or specific temporal properties at the time scales associated with the known QPO frequencies in comparison with the corresponding PDS. Even within our small object sample we find both similarities as well as significant and subtle differences in the form of the ACFs both amongst black holes and between black holes and neutron stars to warrant a closer look at the QPO phenomenon in the time domain: The QPO features manifest as an oscillatory behavior of the ACF at lags near zero; the oscillation damps exponentially on time scales equal to the inverse QPO width to a level of a percent or so. In black holes this oscillatory behavior is preserved and easily discerned at much longer lags while this is not the case for the neutron star system Cyg X-2. The ACF of GRS 1915+105 provides an exception to this general behavior in that its decay is linear in time indicating an undamped oscillation of coherent phase. We present simple ad hoc models that reproduce these diverse time domain behaviors and we speculate that their origin is the phase coherence of the underlying oscillation. It appears plausible that time domain analyses, complementary to the more common frequency domain ones, could impose tighter constraints and provide clues for the driving mechanisms behind the QPO phenomenon.

I outline a new model of particle acceleration in the current sheet separating the closed from the open field lines in the force-free model of pulsar magnetospheres, based on reconnection at the light cylinder and "auroral" acceleration occurring in the return current channel that connects the light cylinder to the neutron star surface. I discuss recent studies of Pulsar Wind Nebulae, which find that pair outflow rates in excess of those predicted by existing theories of pair creation occur, and use those results to point out that dissipation of the magnetic field in a pulsar's wind upstream of the termination shock is restored to life as a viable model for the solution of the "$\sigma$" problem as a consequence of the lower wind 4-velocity implied by the larger mass loading.

Broad-band (0.8-70 keV) spectra of the persistent X-ray emission from 9 magnetars were obtained with Suzaku, including 3 objects in apparent outburst. The soft X-ray component was detected from all of them, with a typical blackbody temperature of kT ~ 0.5 keV, while the hard-tail component, dominating above ~10 keV, was detected at ~1 mCrab intensity from 7 of them. Therefore, the spectrum composed of a soft emission and a hard-tail component may be considered to be a common property of magnetars, both in their active and quiescent states. Wide-band spectral analyses revealed that the hard-tail component has a 1-60 keV flux, Fh, comparable to or even higher than that carried by the 1-60 keV soft component, Fs. The hardness ratio of these objects, defined as xi=Fh/Fs, was found to be tightly anti-correlated with their characteristic age tau as xi=(3.3+/-0.3)x(tau/1 kyr)^(-0.67+/-0.04) with a correlation coefficient of -0.989, over the range from xi~10 to xi~0.1. Magnetars in outburst states were found to lie on the same correlation as relatively quiescent ones. This hardness ratio is also positively correlated with their surface magnetic fields with a correlation coefficient of 0.873. In addition, the hard-tail component becomes harder towards sources with older characteristic ages, with the photon index changing from ~1.7 to ~0.4.

Context. The Crab nebula has been used as a celestial calibration source of the X-ray flux and spectral shape for many years by X-ray astronomy missions. However, the object is often too bright for current and future missions equipped with instruments with improved sensitivity. Aims. We use G21.5-0.9 as a viable, fainter substitute to the Crab, which is another pulsar-wind nebula with a time-constant powerlaw spectrum with a flux of a few milli Crab in the X-ray band. Using this source, we conduct a cross-calibration study of the instruments onboard currently active observatories: Chandra ACIS, Suzaku XIS, Swift XRT, XMM-Newton EPIC (MOS and pn) for the soft-band, and INTEGRAL IBIS-ISGRI, RXTE PCA, and Suzaku HXD-PIN for the hard band. Methods. We extract spectra from all the instruments and fit them under the same astrophysical assumptions. We compare the spectral parameters of the G21.5-0.9 model: power-law photon index, H-equivalent column density of the interstellar photoelectric absorption, flux in the soft (2-8 keV) or hard (15-50 keV) energy band. Results. We identify the systematic differences in the best-fit parameter values unattributable to the statistical scatter of the data alone. We interpret these differences as due to residual cross-calibration problems. The differences can be as large as 20% and 9% for the soft-band flux and power-law index, respectively, and 46% for the hard-band flux. The results are plotted and tabulated as a useful reference for future calibration and scientific studies using multiple missions.

Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; E>100 GeV) {\gamma}-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discovery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to further investigate the nature of this {\gamma}-ray source. Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spectroscopy. Results. A new, hard spectrum VHE {\gamma}-ray source, HESSJ1026-582, was discovered with a statistical significance of 7{\sigma}. It is positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/{\gamma}-ray characteristics of the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar PSRJ1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar and/or the open cluster remains a possibility.

Though pulsars spin regularly, the differences between the observed and predicted ToA (time of arrival), known as "timing noise", can still reach a few milliseconds or more. We try to understand the noise in this paper. As proposed by Xu & Qiao in 2001, both dipole radiation and particle emission would result in pulsar braking. Accordingly, possible fluctuation of particle current flow is suggested here to contribute significant ToA variation of pulsars. We find that the particle emission fluctuation could lead to timing noise which can't be eliminated in timing process, and that a longer period fluctuation would arouse a stronger noise. The simulated timing noise profile and amplitude are in accord with the observed timing behaviors on the timescale of years.

About 25 isolated neutron stars (INSs) are now detected in the optical domain, mainly thanks to the HST and to VLT-class telescopes. The European Extremely Large Telescope (E-ELT) will yield ~100 new identifications, many of which from the follow-up of SKA, IXO, and Fermi observations. Moreover, the E-ELT will allow to carry out, on a much larger sample, INS observations which still challenge VLT-class telescopes, enabling studies on the structure and composition of the NS interior, of its atmosphere and magnetosphere, as well as to search for debris discs. In this contribution, I outline future perspectives for NS optical astronomy with the E-ELT.

Significant researches of compact stars are currently focused on two kinds of enigma sources: anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs). Although AXPs/SGRs are popularly thought to be magnetars, other models (e.g. accretion model) to understand the observations can still not be ruled out. It is worth noting that a null result of Fermi/LAT observation of AXP 4U 0142+61 has been reported recently by Mus and Gogus. We propose here that Fermi/LAT observation may result in distinguishing magnetar model and accretion model for AXPs and SGRs. We address that this null result of AXP 4U 0142+61 would prefer the accretion model.

Using the population synthesis binary evolution code StarTrack, we present theoretical rates and delay times of Type Ia supernovae arising from various formation channels. These channels include binaries in which the exploding white dwarf reaches the Chandrasekhar mass limit (DDS, SDS, and helium-rich donor scenario) as well as the sub-Chandrasekhar mass scenario, in which a white dwarf accretes from a helium-rich companion and explodes as a SN Ia before reaching the Chandrasekhar mass limit. We find that using a common envelope parameterization employing energy balance with alpha=1 and lambda=1, the supernova rates per unit mass (born in stars) of sub-Chandrasekhar mass SNe Ia exceed those of all other progenitor channels at epochs t=0.7 - 4 Gyr for a burst of star formation at t=0. Additionally, the delay time distribution of the sub-Chandrasekhar model can be divided in to two distinct evolutionary channels: the `prompt' helium-star channel with delay times < 500 Myr, and the `delayed' double white dwarf channel with delay times > 800 Myr spanning up to a Hubble time. These findings are in agreement with recent observationally-derived delay time distributions which predict that a large number of SNe Ia have delay times < 1 Gyr, with a significant fraction having delay times < 500 Myr. We find that the DDS channel is also able to account for the observed rates of SNe Ia. However, detailed simulations of white dwarf mergers have shown that most of these mergers will not lead to SNe Ia but rather to the formation of a neutron star via accretion-induced collapse. If this is true, our standard population synthesis model predicts that the only progenitor channel which can account for the rates of SNe Ia is the sub-Chandrasekhar mass scenario, and none of the other progenitors considered can fully account for the observed rates.

(Abridged) We present an analysis of X-ray observations made of the Galactic supernova remnants (SNRs) HB21 (G89.0+4.7) and CTB 1 (G116.9+0.2), two well-known mixed-morphology (MM) SNRs. We find a marked contrast between the X-ray properties of these SNRs: for HB21, the extracted ASCA spectra of the northwest and southeast regions of the X-ray emitting plasma can be fit with a single thermal model with marginally enhanced silicon and sulfur abundances. For both of these regions, the derived column density and temperature are N_H~0.3x10^22 cm^-2 and kT~0.7 keV, respectively. No significant spatial differences in temperature or elemental abundances between the two regions are detected and the X-ray-emitting plasma in both regions is close to ionization equilibrium. Our Chandra spectral analysis of CTB 1 reveals that this source is likely an oxygen-rich SNR with enhanced abundances of oxygen and neon. The extracted ASCA spectra for the southwestern and northeastern regions of CTB 1 cannot be fit with a single thermal component. Based on our fits to these spectra, we derive a column density N_H~0.6x10^22 cm^-2 and a temperature for the soft thermal component of kT_soft~0.28 keV. The hard emission from the southwest may be modeled with either a thermal component (kT_hard~3 keV) or by a power law component (Gamma~2-3) while the hard emission from the northeast may be modeled with a power law component (Gamma~1.4). We have also extracted ASCA GIS spectra of the discrete X-ray source 1WGA J0001.4+6229 which is seen in projection toward CTB 1. These spectra are best fit using a power-law model with a photon index Gamma=2.2^{+0.5}_{-1.2} which is typical for featureless power-law continua produced by rotation-powered pulsars. This source may be a neutron star associated with CTB 1.

For Crab-like pulsars we consider the synchrotron mechanism influenced by relativistic effects of rotation to study the production of the very high energy (VHE) pulsed radiation. The process of quasi-linear diffusion (QLD) is applied to prevent the damping of the synchrotron emission due to extremely strong magnetic field. By examining the kinetic equation governing the QLD, apart from the synchrotron radiative force, we taken into account the the so-called reaction force, that is responsible for corotation and influences plasma processes in the nearby zone of the light cylinder (LC) surface. We have found that the relativistic effects of rotation significantly change efficiency of the quasi-linear diffusion. In particular, examining magnetospheric parameters typical for Crab-like pulsars, it has been shown that unlike the situation, where relativistic effects of rotation are not important, on the LC surface, the relativistic electrons via the synchrotron mechanism may produce photons even in the TeV domain. It is shown that the VHE radiation is strongly correlated with the relatively low frequency emission.

A stochastic gravitational wave background causes the apparent positions of distant sources to fluctuate, with angular deflections of order the characteristic strain amplitude of the gravitational waves. These fluctuations may be detectable with high precision astrometry, as first suggested by Braginsky et al. in 1990. Several researchers have made order of magnitude estimates of the upper limits obtainable on the gravitational wave spectrum \Omega_gw(f), at frequencies of order f ~ 1 yr^-1, both for the future space-based optical interferometry missions GAIA and SIM, and for VLBI interferometry in radio wavelengths with the SKA. For GAIA, tracking N ~ 10^6 sources over a time of T ~ 1 yr with an angular accuracy of \Delta \theta ~ 10 \mu as would yield a sensitivity level of \Omega_gw ~ (\Delta \theta)^2/(N T^2 H_0^2) ~ 10^-6, which would be comparable with pulsar timing. In this paper we take a first step toward firming up these estimates by computing in detail the statistical properties of the angular deflections caused by a stochastic background. We compute analytically the two point correlation function of the deflections on the sphere, and the spectrum as a function of frequency and angular scale. The fluctuations are concentrated at low frequencies (for a scale invariant stochastic background), and at large angular scales, starting with the quadrupole. The magnetic-type and electric-type pieces of the fluctuations have equal amounts of power.

It is still a matter of debate to understand the equation of state of cold supra-nuclear matter in compact stars because of unknown on-perturbative strong interaction between quarks. Nevertheless, it is speculated from an astrophysical view point that quark clusters could form in cold quark matter due to strong coupling at realistic baryon densities. Although it is hard to calculate this conjectured matter from first principles, one can expect the inter-cluster interaction to share some general features to nucleon-nucleon interaction. We adopt a two-Gaussian component soft-core potential with these general features and show that quark clusters can form stable simple cubic crystal structure if we assume Gaussian form wave function. With this parameterizing, Tolman-Oppenheimer-Volkoff equation is solved with reasonable constrained parameter space to give mass-radius relation of crystalline solid quark star. With baryon densities truncated at 2 times nuclear density at surface and range of interaction fixed at 2fm we can reproduce similar mass-radius relation to that obtained with bag model equations of state. The maximum mass ranges from about 0.5 to 3 solar mass. Observed maximum pulsar mass (about 2 solar mass) is then used to constrain parameters of this simple interaction potential.

Sedimentation of the neutron rich isotope $^{22}$Ne may be an important source of gravitational energy during the cooling of white dwarf stars. This depends on the diffusion constant for $^{22}$Ne in strongly coupled plasma mixtures. We calculate self-diffusion constants $D_i$ from molecular dynamics simulations of carbon, oxygen, and neon mixtures. We find that $D_i$ in a mixture does not differ greatly from earlier one component plasma results. For strong coupling (coulomb parameter $\Gamma>$ few), $D_i$ has a modest dependence on the charge $Z_i$ of the ion species, $D_i \propto Z_i^{-2/3}$. However $D_i$ depends more strongly on $Z_i$ for weak coupling (smaller $\Gamma$). We conclude that the self-diffusion constant $D_{\rm Ne}$ for $^{22}$Ne in carbon, oxygen, and neon plasma mixtures is accurately known so that uncertainties in $D_{\rm Ne}$ should be unimportant for simulations of white dwarf cooling.

The Collapsar model, in which a massive star (greater than 20 solar masses) fails to produce a SN and forms a BH, provides the main framework for understanding long Gamma-Ray Bursts (GRB) and the accompanying hypernovae (HN). However, single massive-star models that explain the population of pulsars, predict cores that rotate too slowly to produce GRBs/HNe. We present a model of binary evolution that allows the formation of Kerr black holes (BH) where the spin of the BH can be estimated from the pre-collapse orbit, and use the Blandford-Znajek (BZ) mechanism to estimate the available energy for a GRB/HN. A population synthesis study shows that this model can account for both, the long GRB and the subluminous GRB populations.

X-ray charge-coupled devices (CCDs) are the workhorse detectors of modern X-ray astronomy. Typically covering the 0.3-10.0 keV energy range, CCDs are able to detect photoelectric absorption edges and K shell lines from most abundant metals. New CCDs also offer resolutions of 30-50 (E/dE), which is sufficient to detect lines in hot plasmas and to resolve many lines shaped by dynamical processes in accretion flows. The spectral capabilities of X-ray CCDs have been particularly important in detecting relativistic emission lines from the inner disks around accreting neutron stars and black holes. One drawback of X-ray CCDs is that spectra can be distorted by photon "pile-up", wherein two or more photons may be registered as a single event during one frame time. We have conducted a large number of simulations using a statistical model of photon pile-up to assess its impacts on relativistic disk line and continuum spectra from stellar-mass black holes and neutron stars. The simulations cover the range of current X-ray CCD spectrometers and operational modes typically used to observe neutron stars and black holes in X-ray binaries. Our results suggest that severe photon pile-up acts to falsely narrow emission lines, leading to falsely large disk radii and falsely low spin values. In contrast, our simulations suggest that disk continua affected by severe pile-up are measured to have falsely low flux values, leading to falsely small radii and falsely high spin values. The results of these simulations and existing data appear to suggest that relativistic disk spectroscopy is generally robust against pile-up when this effect is modest.

We show that a pair of thermal, antipodal hot-spots on the neutron star surface is able to fully account for the pulsar's double blackbody spectrum and energy-dependent pulse profile, including the observed 180 degree phase reversal at approximately 1.2 keV. By comparing the observed pulse modulation and phase to the model predictions, we strongly constrain the hot-spot pole (xi) and the line-of-sight (psi) angles with respect to the spin axis. For a nominal radius of R = 12 km and distance D = 2.2 kpc, we find (xi,psi) = (86d,6d), with 1-sigma error ellipse of (2d,1d); this solution is degenerate in the two angles. The best-fit spectral model for this geometry requires that the temperatures of the two emission spots differ by a factor of 2 and their areas by a factor of ~ 20. Including a cosine-beamed pattern for the emitted intensity modifies the result, decreasing the angles to (84d,3d); however this model is not statistically distinguishable from the isotropic emission case. We also present a new upper limit on the period derivative of Pdot < 3.5E-16 (2-sigma), which limits the global dipole magnetic field to B_s < 2.0E11 G, confirming PSR J0821-4300 as an "anti-magnetar." We discuss the results in the context of observations and theories of nonuniform surface temperature on isolated NSs of both weak and strong magnetic field. To explain the nonuniform temperature of PSR J0821-4300 may require a crustal field that is much stronger than the external, global dipole field.

To better characterize the abundance patterns produced by the r-process, we have derived new abundances or upper limits for the heavy elements zinc (Zn), yttrium (Y), lanthanum (La), europium (Eu), and lead (Pb). Our sample of 161 metal-poor stars includes new measurements from 88 high resolution and high signal-to-noise spectra obtained with the Tull Spectrograph on the 2.7m Smith Telescope at McDonald Observatory, and other abundances are adopted from the literature. We use models of the s-process in AGB stars to characterize the high Pb/Eu ratios produced in the s-process at low metallicity, and our new observations then allow us to identify a sample of stars with no detectable s-process material. In these stars, we find no significant increase in the Pb/Eu ratios with increasing metallicity. This suggests that s-process material was not widely dispersed until the overall Galactic metallicity grew considerably, perhaps even as high as [Fe/H]=-1.4. We identify a dispersion of at least 0.5 dex in [La/Eu] in metal-poor stars with [Eu/Fe]<+0.6 attributable to the r-process, suggesting that there is no unique "pure" r-process elemental ratio among pairs of rare earth elements. We confirm earlier detections of an anti-correlation between Y/Eu and Eu/Fe bookended by stars strongly enriched in the r-process (e.g., CS 22892-052) and those with deficiencies of the heavy elements (e.g., HD 122563). We can reproduce the range of Y/Eu ratios using simulations of high-entropy neutrino winds of core-collapse supernovae that include charged-particle and neutron-capture components of r-process nucleosynthesis. The heavy element abundance patterns in most metal-poor stars do not resemble that of CS 22892-052, but the presence of heavy elements such as Ba in nearly all metal-poor stars without s-process enrichment suggests that the r-process is a common phenomenon.

Anomalous X-ray Pulsars and Soft-Gamma Repeaters groups are magnetar candidates featuring low characteristic ages ($\tau = {P\over{2 {\dot P}}}$). At least some of them they should still be associated with the remnants of the explosive events in which they were born, giving clues to the type of events leading to their birth and the physics behind the apparent high value of the magnetar magnetic fields. To explain the high values of $B$, a self-consistent picture of field growth also suggests that energy injection into the SNR is large and unavoidable, in contrast with the evolution of {\it conventional} SNR. This modified dynamics, in turn, has important implications for the proposed associations. We show that this scenario yields low ages for the new candidates CXOU J171405.7-381031/CTB 37B and XMMU J173203.3-344518/G353.6-0.7, and predicted values agree with recently found ${\dot P}$, giving support to the overall picture.

We investigate the physical and chemical evolution of population II stars with initial masses in the range 6.5-8 Msun, which undergo an off centre carbon ignition under partially degenerate conditions, followed by a series of thermal pulses, and supported energetically by a CNO burning shell, above a O-Ne degenerate core. In agreement with the results by other research groups, we find that the O-Ne core is formed via the formation of a convective flame that proceeds to the centre of the star. The evolution which follows is strongly determined by the description of the mass loss mechanism. Use of the traditional formalism with the super-wind phase favours a long evolution with many thermal pulses, and the achievement of an advanced nucleosynthesis, due the large temperatures reached by the bottom of the external mantle. Use of a mass loss recipe with a strong dependence on the luminosity favours an early consumption of the stellar envelope, so that the extent of the nucleosynthesis, and thus the chemical composition of the ejecta, is less extreme. The implications for the multiple populations in globular clusters are discussed. If the "extreme" populations present in the most massive clusters are a result of direct formation from the super-AGB ejecta, their abundances may constitute a powerful way of calibrating the mass loss rate of this phase. This calibration will also provide informations on the fraction of super-AGBs exploding as single e-capture supernova, leaving a neutron star remnant in the cluster.

We present abundances of Mn, Cu, Zn and various light and heavy elements for a sample of barium and normal giant stars, and present correlations between abundances contributed to different degrees by the weak-s, main-s, and r-processes of neutron capture, between Fe-peak elements and heavy elements. Data from the literature are also considered in order to better study the abundance pattern of peculiar stars. The stellar spectra were observed with FEROS/ESO, and the stellar atmospheric parameters of the 8 barium giant stars and 6 normal giants analyzed by us lie in the range 4300 < T_eff/K < 5300, -0.7 < [Fe/H] <= 0.12 and 1.5 <= log g < 2.9. Carbon and nitrogen abundances were derived by the spectral synthesis of molecular bands of C_2, CH and CN. For all the other elements, atomic lines were used to perform the spectral synthesis. A very large scatter was found when data from the literature were considered, mainly for the Mn abundances. We found that [Zn/Fe] correlates well with the heavy element excesses, its abundance clearly increasing as the heavy element excesses increase, a trend not shown by the [Cu/Fe] and [Mn/Fe] ratios. Also, the excess of Mn, Cu and Zn relative to heavy elements usually show an increasing trend toward higher metallicities. Our results suggest that a larger fraction of the synthesis of Zn is owed to massive stars than is the case of Cu, and that the contribution of the main-s process to the synthesis of both elements is small. We also conclude that Mn is mostly synthesized by SN Ia, and that a non-negligible fraction of the synthesis of Mn, Cu and Zn is owed to the weak s-process.

Westerlund 1 is one of the most massive young clusters known in the Local Group, with an age of 3-5 Myr. It contains an assortment of rare evolved massive stars, such as blue, yellow and red supergiants, Wolf-Rayet stars, a luminous blue variable, and a magnetar, as well as 4 massive eclipsing binary systems (Wddeb, Wd13, Wd36, WR77o, see Bonanos 2007). The eclipsing binaries present a rare opportunity to constrain evolutionary models of massive stars, the distance to the cluster and furthermore, to determine a dynamical lower limit for the mass of a magnetar progenitor. Wddeb, being a detached system, is of great interest as it allows determination of the masses of 2 of the most massive unevolved stars in the cluster. We have analyzed spectra of all 4 eclipsing binaries, taken in 2007-2008 with the 6.5 meter Magellan telescope at Las Campanas Observatory, Chile, and present fundamental parameters (masses, radii) for their component stars.

The heating associated with the deposition of $\gamma$-rays in an accretion disk is proposed as a mechanism to facilitate the transformation of a low mass X-ray binary to the radio millisecond pulsar phase. The $\gamma$-ray emission produced in the outer gap accelerator in the pulsar magnetosphere likely irradiates the surrounding disk, resulting in its heating and to the possible escape of matter from the system. We apply the model to PSR J1023+0038, which has recently been discovered as a newly born rotation powered millisecond pulsar. The predicted $\gamma$-ray luminosity $\sim 6 \times 10^{34}~\mathrm{erg~s^{-1}}$ can be sufficient to explain the disappearance of the truncated disk existing during the 8~month$\sim 2$~yr period prior to the 2002 observations of J1023+0038 and the energy input required for the anomalously bright optical emission of its companion star.

We simulate head-on collisions from rest at large separation of binary white dwarf -- neutron stars (WDNSs) in full general relativity. Our study serves as a prelude to our analysis of the circular binary WDNS problem. We focus on compact binaries whose total mass exceeds the maximum mass that a cold degenerate star can support, and our goal is to determine the fate of such systems. A fully general relativistic hydrodynamic computation of a realistic WDNS head-on collision is prohibitive due to the large range of dynamical time scales and length scales involved. For this reason, we construct an equation of state (EOS) which captures the main physical features of NSs while, at the same time, scales down the size of WDs. We call these scaled-down WD models "pseudo-WDs (pWDs)". Using pWDs, we can study these systems via a sequence of simulations where the size of the pWD gradually increases toward the realistic case. We perform two sets of simulations; One set studies the effects of the NS mass on the final outcome, when the pWD is kept fixed. The other set studies the effect of the pWD compaction on the final outcome, when the pWD mass and the NS are kept fixed. All simulations show that 14%-18% of the initial total rest mass escapes to infinity. All remnant masses still exceed the maximum rest mass that our cold EOS can support (1.92 solar masses), but no case leads to prompt collapse to a black hole. This outcome arises because the final configurations are hot. All cases settle into spherical, quasiequilibrium configurations consisting of a cold NS core surrounded by a hot mantle, resembling Thorne-Zytkow objects. Extrapolating our results to realistic WD compactions, we predict that the likely outcome of a head-on collision of a realistic, massive WDNS system will be the formation of a quasiequilibrium Thorne-Zytkow-like object.

The emission-line Be star HD 215227 lies within the positional error circle of the newly identified gamma-ray source AGL J2241+4454. We present new blue spectra of the star, and we point out the morphological and variability similarities to other Be binaries. An analysis of the available optical photometry indicates a variation with a period of 60.37 +/- 0.04 d, which may correspond to an orbital modulation of the flux from the disk surrounding the Be star. The distance to the star of 2.6 kpc and its relatively large Galactic latitude suggest that the binary was ejected from the plane by a supernova explosion that created the neutron star or black hole companion. The binary and runaway properties of HD 215227 make it an attractive candidate as the optical counterpart of AGL J2241+4454 and as a new member of the small class of gamma-ray emitting binaries.

The effects of pinning between fluxoids and vortices in the core of a neutron star, on the dynamics of the core neutron superfluid are considered. The pinning impedes, but does not absolutely block, any radial as well as {\em azimuthal} motion of the neutron vortices with respect to the lattice of fluxoids. The time scale for the coupling of rotation of the core superfluid to the rest of the star is calculated, allowing for the effect of the finite frictional force on the neutron vortices due to their pinning with the fluxoids. This turns out to be the dominant mechanism for the coupling of the core of a neutron star to its crust, as compared to the role of electron scattering, for most cases of interest. Furthermore, different behaviors for the post-glitch response of the core superfluid are distinguished that might be tested against the relevant observational data. Also, a conceptually important case (and controversial too, in the earlier studies on the role of the crustal superfluid) is realized where a superfluid may remain decoupled in spite of an spinning up of its vortices.

We analyzed XMM-Newton observations of two center-filled supernova remnants, G27.8+0.6 and G28.8+1.5, to search for pulsars/pulsar wind nebulae (PWNe) and other types of neutron stars associated with these remnants. We discovered a PWN candidate within the extent of the centrally-peaked radio emission of G27.8+0.6. The X-ray morphology of the PWN candidate and its position with respect to the host supernova remnant suggest that the alleged pulsar might be moving away from the supernova remnant's center with a transverse velocity of around 100-200 km/s. The majority of the detected X-ray point sources in both fields are classified as main-sequence stars, based on their bright optical counterparts and relatively soft X-ray spectra. The remaining medium-hard and hard sources are most probably either AGNs or cataclysmic variables, although we cannot completely rule out the possibility that some of these sources are neutron stars.

The impact of magnetic field decay on radiative activity of quaking neutron star undergoing Lorentz-force-driven torsional seismic vibrations about axis of its dipole magnetic moment is studied. We found that monotonic depletion of internal magnetic field pressure is accompanied by the loss of vibration energy of the star that causes its vibration period to lengthen at a rate proportional to the rate of magnetic field decay. Particular attention is given to the magnetic-field-decay induced conversion of the energy of differentially rotational Alfven vibrations into the energy of oscillating magneto-dipole radiation. A set of representative examples of magnetic field decay illustrating the vibration energy powered emission with elongating periods produced by quaking neutron star are considered and discussed in the context of theory of magnetars.

We have simulated a Galactic population of young pulsars and compared with the Fermi LAT sample, constraining the birth properties, beaming and evolution of these spin-powered objects. Using quantitative tests of agreement with the distributions of observed spin and pulse properties, we find that short birth periods P_0 ~ 50ms and gamma-ray beams arising in the outer magnetosphere, dominated by a single pole, are strongly preferred. The modeled relative numbers of radio-detected and radio-quiet objects agree well with the data. Although the sample is local, extrapolation to the full Galaxy implies a gamma-ray pulsar birthrate 1/(59 yr). This is shown to be in good agreement with the estimated Galactic core collapse rate and with the local density of OB star progenitors. We give predictions for the numbers of expected young pulsar detections if Fermi LAT observations continue 10 years and estimate the (small) contributions of unresolved young pulsars to the Galactic gamma-ray background.

PSR B1259-63/LS 2883 is a binary system where a 48 ms pulsar orbits a massive Be star with a highly eccentric orbit (e=0.87) with a period of 3.4 years. The system exhibits variable, non-thermal radiation visible from radio to very high energies (VHE) around periastron passage. This radiation is thought to originate from particles accelerated in the shock region between the pulsar wind (PW) and stellar outflows. The consistency of the H.E.S.S. data with the inverse Compton (IC) scenario is studied in the context of dominant orbital phase dependent adiabatic losses. The dependence of the observed TeV flux with the separation distance is analyzed. Model calculations based on IC scattering of shock accelerated PW electrons and UV photons are performed. Different non-radiative cooling profiles are suggested for the primary particle population to account for the variable TeV flux. The TeV fluxes obtained with H.E.S.S. in the years 2004 and 2007 seem to be only dependent on the binary separation. The presented results hint at a peculiar non-radiative cooling profile around periastron dominating the VHE emission in PSR B1259-63. The location of the stellar disc derived from this non-radiative cooling profile is in good agreement with that inferred from radio observations.

We have detected an Halpha bow shock nebula around PSR J1741-2054, a pulsar discovered through its GeV gamma-ray pulsations. The pulsar is only ~1.5" behind the leading edge of the shock. Optical spectroscopy shows that the nebula is non-radiative, dominated by Balmer emission. The Halpha images and spectra suggest that the pulsar wind momentum is equatorially concentrated and implies a pulsar space velocity ~150km/s, directed 15+/-10deg out of the plane of the sky. The complex Halpha profile indicates that different portions of the post-shock flow dominate line emission as gas moves along the nebula and provide an opportunity to study the structure of this unusual slow non-radiative shock under a variety of conditions. CXO ACIS observations reveal an X-ray PWN within this nebula, with a compact ~2.5" equatorial structure and a trail extending several arcmin behind. Together these data support a close (<0.5kpc) distance, a spin geometry viewed edge-on and highly efficient gamma-ray production for this unusual, energetic pulsar.

An unusual Eddington-limited thermonuclear X-ray burst was detected from the accreting neutron star in 2S 0918-549 with the Rossi X-ray Timing Explorer. The burst commenced with a brief (40 ms) precursor and maintained near-Eddington fluxes during the initial 77 s. These characteristics are indicative of a nova-like expulsion of a shell from the neutron star surface. Starting 122 s into the burst, the burst shows strong (87 +/- 1% peak-to-peak amplitude) achromatic fluctuations for 60 s. We speculate that the fluctuations are due to Thompson scattering by fully-ionized inhomogeneities in a resettling accretion disk that was disrupted by the effects of super-Eddington fluxes. An expanding shell may be the necessary prerequisite for the fluctuations.

The existence of stable magnetic configurations in white dwarfs, neutron stars and various non-convective stellar {regions} is now well recognized. It has recently been shown numerically that various families of equilibria, including axisymmetric mixed poloidal-toroidal configurations, are stable. Here we test the stability of an analytically-derived non force-free magnetic equilibrium, using three-dimensional magnetohydrodynamic simulations: the mixed configuration is compared with the dynamical evolution of its purely poloidal and purely toroidal components, both known to be unstable. The mixed equilibrium shows no sign of instability under white noise perturbations. {This configuration therefore provides a good description of magnetic equilibrium topology inside non-convective stellar objects and will be useful to initialize magneto-rotational transport in stellar evolution codes.

The most massive neutron stars constrain the behavior of ultra-dense matter, with larger masses possible only for increasingly stiff equations of state. Here, we present evidence that the black widow pulsar, PSR B1957+20, has a high mass. We took spectra of its strongly irradiated companion and found an observed radial-velocity amplitude of K_obs=324+/-3 km/s. Correcting this for the fact that, due to the irradiation, the center of light lies inward relative to the center of mass, we infer a true radial-velocity amplitude of K_2=353+/-4 km/s and a mass ratio q=M_PSR/M_2=69.2+/-0.8. Combined with the inclination i=65+/-2 deg inferred from models of the lightcurve, our best-fit pulsar mass is M_PSR=2.40+/-0.12 M_sun. We discuss possible systematic uncertainties, in particular in the lightcurve modeling. Taking an upper limit of i<85 deg based on the absence of radio eclipses at high frequency, combined with a conservative lower-limit to the motion of the center of mass, K_2>343 km/s (q>67.3), we infer a lower limit to the pulsar mass of M_PSR>1.66 M_sun.

Fundamental constants are a cornerstone of our physical laws. Any constant varying in space and/or time would reflect the existence of an almost massless field that couples to matter. This will induce a violation of the universality of free fall. It is thus of utmost importance for our understanding of gravity and of the domain of validity of general relativity to test for their constancy. We thus detail the relations between the constants, the tests of the local position invariance and of the universality of free fall. We then review the main experimental and observational constraints that have been obtained from atomic clocks, the Oklo phenomenon, Solar system observations, meteorites dating, quasar absorption spectra, stellar physics, pulsar timing, the cosmic microwave background and big bang nucleosynthesis. At each step we describe the basics of each system, its dependence with respect to the constants, the known systematic effects and the most recent constraints that have been obtained. We then describe the main theoretical frameworks in which the low-energy constants may actually be varying and we focus on the unification mechanisms and the relations between the variation of different constants. To finish, we discuss the more speculative possibility of understanding their numerical values and the apparent fine-tuning that they confront us with.

Identification of gamma-ray-emitting Galactic sources is a long-standing problem in astrophysics. One such source, 1AGL J2022+4032, coincident with the interior of the radio shell of the supernova remnant Gamma Cygni (SNR G78.2+2.1) in the Cygnus Region, has recently been identified by Fermi as a gamma-ray pulsar, LAT PSR J2021+4026. We present long-term observations of 1AGL J2022+4032 with the AGILE gamma-ray telescope, measuring its flux and light curve. We compare the light curve of 1AGL J2022+4032 with that of 1AGL J2021+3652 (PSR J2021+3651), showing that the flux variability of 1AGL J2022+4032 appears to be greater than the level predicted from statistical and systematic effects and producing detailed simulations to estimate the probability of the apparent observed variability. We evaluate the possibility that the gamma-ray emission may be due to the superposition of two or more point sources, some of which may be variable, considering a number of possible counterparts. We consider the possibility of a nearby X-ray quiet microquasar contributing to the flux of 1AGL J2022+4032 to be more likely than the hypotheses of a background blazar or intrinsic gamma-ray variabilty of LAT PSR J2021+4026.

A survey of the Galactic plane in the region $-60\degree \leq l \leq 30\degree$, $|b| \leq 0.25\degree$ was carried out using the seven-beam Parkes Methanol Multibeam (MMB) receiver, which operates at a frequency of 6.5 GHz. Three pulsars were discovered, and 16 previously known pulsars detected. In this paper we present two previously-unpublished discoveries, both with extremely high dispersion measures, one of which is very close, in angular distance, to the Galactic centre. The survey data also contain the first known detection, at radio frequencies, of the radio magnetar PSR J1550-5418. Our survey observation was made 46 days prior to that previously published and places constraints on the beginning of pulsed radio emission from the source.
The detection of only three previously undiscovered pulsars argues that there are few pulsars in the direction of the inner Galaxy whose flux density spectrum is governed by a flat power law. However, these pulsars would be likely to remain undetected at lower frequencies due to the large amount of scatter broadening which affects pulsars with high values of dispersion measure. Surveys with future telescopes at high observing frequencies will, therefore, play an important role in the discovery of pulsars at the Galactic centre. By simulating pulsar surveys of the Galaxy with Phase 1 SKA at frequencies of 1.4 GHz and 10 GHz, we find that high-frequency observations are the only way to discover and observe the Galactic-centre pulsar population.

We perform population synthesis studies of different types of neutron stars taking into account the magnetic field decay. For the first time, we confront our results with observations using {\it simultaneously} the Log N -- Log S distribution for nearby isolated neutron stars, the Log N -- Log L distribution for magnetars, and the distribution of radio pulsars in the $P$ -- $\dot P$ diagram. We find that our theoretical model is consistent with all sets of data if the initial magnetic field distribution function follows a log-normal law with $<\log (B_0/[G]) > \sim 13.25$ and $\sigma_{\log B_0}\sim 0.6$. The typical scenario includes about 10% of neutron stars born as magnetars, significant magnetic field decay during the first million years of a NS life. Evolutionary links between different subclasses may exist, although robust conclusions are not yet possible.
We apply the obtained field distribution and the model of decay to study long-term evolution of neuton stars till the stage of accretion from the interstellar medium. It is shown that though the subsonic propeller stage can be relatively long, initially highly magnetized neutron stars ($B_0 > \sim 10^{13}$ G) reach the accretion regime within the Galactic lifetime if their kick velocities are not too large. The fact that in previous studies made $>$10 years ago, such objects were not considered results in a slight increase of the Accretor fraction in comparison with earlier conclusions. Most of the neutron stars similar to the Magnificent seven are expected to become accreting from the interstellar medium after few billion years of their evolution. They are the main predecestors of accreting isolated neutron stars.