Neutron star magnetic fields may have polar caps (PC) that are offset from the dipole axis, through field-line sweepback near the light cylinder or non-symmetric currents within the star. The effects of such offsets on electron-positron pair cascades are investigated, using simple models of dipole magnetic fields with small distortions that shift the PCs by different amounts or directions. Using a Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, multiplicity and energy flux across the PC, as well as the trends in pair flux and pair energy flux with spin-down luminosity, L_{sd}. We also give an estimate of the distribution of heating flux from returning positrons on the PC for different offsets. We find that even modest offsets can produce significant increases in pair multiplicity, especially for pulsars that are near or beyond the pair death lines for centered PCs, primarily because of higher accelerating fields. Pair spectra cover several decades in energy, with the spectral range of millisecond pulsars (MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets allow significant extension of all spectra to lower pair energies. We find that the total PC pair luminosity L_{pair} is proportional to L_{sd}, with L_{pair} ~ 10^{-3} L_{sd} for normal pulsars and L_{pair} ~ 10^{-2} L_{sd} for MSPs. Remarkably, the total PC heating luminosity for even large offsets increases by less than a factor of two, even though the PC area increases by much larger factors, because most of the heating occurs near the magnetic axis.

We investigate the effect of a distorted neutron star dipole magnetic field on pulsar pair cascade multiplicity and pair death lines. Using a simple model for a distorted dipole field that produces an offset polar cap, we derive the accelerating electric field above the polar cap in space charge limited flow. We find that even a modest azimuthally asymmetric distortion can significantly increase the accelerating electric field on one side of the polar cap and, combined with a smaller field line radius of curvature, leads to larger pair multiplicity. The death line for producing pairs by curvature radiation moves downward in the P-Pdot diagram, allowing high pair multiplicities in a larger percentage of the radio pulsar population. These results could have important implications for the radio pulsar population, high energy pulsed emission and the pulsar contribution to cosmic ray positrons.

We propose a simple analytic model for the innermost (within the light cylinder of canonical radius) structure of open-magnetic-field lines of a rotating neutron star (NS) with relativistic outflow of charged particles (electrons/positrons) and arbitrary angle between the NS spin and magnetic axes. We present the self-consistent solution of Maxwell's equations for the magnetic field and electric current in the pair-starved regime where the density of electron-positron plasma generated above the pulsar polar cap is not sufficient to completely screen the accelerating electric field and thus establish the E \cdot B = 0 condition above the pair-formation front up to very high altitudes within the light cylinder. The proposed model may provide a theoretical framework for developing a refined model of the global pair-starved pulsar magnetosphere.

A number of rotation-powered millisecond pulsars are powerful sources of X-ray emission. We present predictions for the spectral characteristics of these sources at gamma-ray energies, using a model for acceleration and pair cascades on open magnetic field lines above the polar caps. Since these pulsars have low surface magnetic fields, the majority do not produce sufficient pairs to completely screen the accelerating electric field allowing particle acceleration to high altitude. The resulting emission above 1 GeV comes from curvature radiation by primary electrons with radiation-reaction-limited Lorentz factors. The spectra are very hard power-laws with exponential cutoffs between 1 and 50 GeV, and the spectral power peaks near the cutoff energy. Millisecond pulsars are thus ideal targets for air-Cherenkov detectors that are able to reach energy thresholds below 50 GeV.

We propose an analytical 3-D model of the open field-line region of a neutron star (NS) magnetosphere. We construct an explicit analytic solution for arbitrary obliquity (angle between the rotation and magnetic axes) incorporating the effects of magnetospheric rotation, relativistic flow of charges (e.g. primary electron beam) along the open field lines, and E X B drift of these charges. Our solution employs the space-charge-limited longitudinal current calculated in the electrodynamic model of Muslimov & Tsygan (1992) and is valid up to very high altitudes nearly approaching the light cylinder. We assume that in the innermost magnetosphere, the NS magnetic field can be well represented by a static magnetic dipole configuration. At high altitudes the open magnetic field lines significantly deviate from those of a static dipole and tend to focus into a cylindrical bundle, swept back in the direction opposite to the rotation, and also bent towards the rotational equator. We briefly discuss some implications of our study to spin-powered pulsars.

The X-ray and gamma-ray spectrum of rotation-powered millisecond pulsars is investigated in a model for acceleration and pair cascades on open field lines above the polar caps. Although these pulsars have low surface magnetic fields, their short periods allow them to have large magnetospheric potential drops, but the majority do not produce sufficient pairs to completely screen the accelerating electric field. The accelerating particles maintain high Lorentz factors and undergo cyclotron resonant absorption of radio emission, that produces and maintains a large pitch angle, resulting in a strong synchrotron component. The resulting spectra consist of several distinct components: curvature radiation from primary electrons dominating from 1 - 100 GeV, synchrotron radiation from primary and secondary electrons dominating up to about 100 MeV, and much weaker inverse-Compton radiation from primary electrons at 0.1 - 1 TeV. We find that the relative size of these components depends on pulsar period, period derivative, and neutron star mass and radius with the level of the synchrotron component also depending sensitively on the radio emission properties. This model is successful in describing the observed X-ray and gamma-ray spectrum of PSR J0218+4232 as synchrotron radiation, peaking around 100 MeV and extending up to a turnover around several GeV. The predicted curvature radiation components from a number of millisecond pulsars, as well as the collective emission from the millisecond pulsars in globular clusters, should be detectable with AGILE and GLAST. We also discuss a hidden population of X-ray-quiet and radio-quiet millisecond pulsars which have evolved below the pair death line, some of which may be detectable by telescopes sensitive above 1 GeV.

A new picture of pulsar high-energy emission is proposed that is different from both the traditional polar cap and outer gap models, but combines elements of each. The slot gap model is based on electron acceleration along the edge of the open field region from the neutron star surface to near the light cylinder. Along the last open field line, the pair formation front rises to very high altitude forming a slot gap, where the accelerating electric field is unscreened by pairs. The resulting radiation features both hollow cones from the lower-altitude pair cascades, seen at small viewing angles, as well as caustic emission on the trailing-edge field lines at high altitude, seen from both poles at large viewing angle. The combination of the small solid angle of slot gap emission (<< 1 sr) with a high probability of viewing the emission predicts that more gamma-ray pulsars could be detected at larger distances. In this picture, many of the positional coincidences of radio pulsars with unidentified EGRET sources become plausible as real associations, as the flux predicted by the slot gap model for many of the pulsars would provide the observed EGRET source flux. The expected probability of seeing radio-quiet gamma-ray pulsars in this model will also be discussed.

We investigate the physical situation above the pulsar polar cap (PC) where the accelerating primaries (electrons) are not capable of producing sufficient numbers of electron-positron pairs at low altitudes (within 1-2 stellar radii above the PC surface) to screen the accelerating electric field, and continue accelerating up to, at least, very high altitudes nearly approaching the light cylinder. We derive an analytic solution for the parallel electric field valid at high altitudes. The solution is based on the physical condition of asymptotic vanishing of the rotationally induced transverse electric field within the magnetic flux tube. This condition constrains the asymptotic value of the effective space charge that determines the distribution of the parallel electric field within the magnetic tube. Our estimates of low- to high-altitude values of the parallel electric field imply the occurrence of a regime of primary acceleration (with the characteristic Lorentz factor up to 1-2 X 10^7) all the way from the PC to the light cylinder limited by curvature-radiation reaction. In this model the primary outflow becomes asymptotically force-free, and may turn into a relativistic wind beyond the light cylinder. Such a solution will apply to both older pulsars producing only inverse Compton scattering pairs and younger very high B pulsars (magnetars). We suggest that pulsars, which are lying below the pair death line, may be radio-quiet gamma-ray sources.

We explore the pulsar slot gap electrodynamics up to very high altitudes, where for most relatively rapidly rotating pulsars both the standard small-angle approximation and the assumption that the magnetic field lines are ideal stream lines break down. We address the importance of the electrodynamic conditions at the slot gap boundaries and the occurrence of a steady-state drift of charged particles across the slot gap field lines at very high altitudes. These boundary conditions and the deviation of particle trajectories from stream lines determine the asymptotic behavior of the scalar potential at all radii from the polar cap to near the light cylinder. As a result, we demonstrate that the steady-state accelerating electric field must approach a small and constant value at high altitude above the polar cap. This parallel electric field is capable of maintaining electrons moving with high Lorentz factors (a few times 10^7) and emitting curvature gamma-ray photons up to nearly the light cylinder. By numerical simulations, we show that primary electrons accelerating from the polar cap surface to high altitude in the slot gap along the outer edge of the open field region will form caustic emission patterns on the trailing dipole field lines. Acceleration and emission in such an extended slot gap may form the physical basis of a model that can successfully reproduce some pulsar high-energy light curves.

The characteristics of the high-energy emission from polar cap accelerators will be discussed. Particles accelerated in the "slot gap" near the polar cap rim will reach altitudes of several stellar radii before initiating pair cascades, producing a wide hollow cone of emission in young pulsars and some millisecond pulsars. Model X-ray and gamma-ray spectra and pulse profiles, based on Monte-Carlo simulations of polar cap pair cascades, will be presented.

We revise the physics of primary electron acceleration in the "slot gap" (SG) above the pulsar polar caps (PCs), a regime originally proposed by Arons and Scharlemann (1979) in their electrodynamic model of pulsar PCs. We employ the standard definition of the SG as a pair-free space between the last open field lines and the boundary of the pair plasma column which is expected to develop above the bulk of the PC. The rationale for our revision is that the proper treatment of primary acceleration within the pulsar SGs should take into account the effect of the narrow geometry of the gap on the electrodynamics within the gap and also to include the effect of inertial frame dragging on the particle acceleration. The combination of the effects of frame dragging and geometrical screening in the gap region naturally gives rise to a regime of extended acceleration, that is not limited to "favorably curved" field lines as in earlier models, and the possibility of multiple-pair production by curvature photons at very high altitudes, up to several stellar radii. The estimated theoretical high-energy luminosities of the SG cascade radiation are in good agreement with the corresponding empirical relationships for gamma-ray pulsars. We illustrate the results of our modeling of the pair cascades and gamma-ray emission from the high altitudes in the SG for the Crab pulsar. The combination of the frame-dragging field and high-altitude SG emission enables both acceleration at the smaller inclination angles and a larger emission beam, both necessary to produce widely-spaced double-peaked profiles.

We investigate the conditions required for the production of electron-positron pairs above a pulsar polar cap (PC) and the influence of pair production on the energetics of the primary particle acceleration. Assuming space-charge limited flow acceleration including the inertial frame-dragging effect, we allow both one-photon and two-photon pair production by either curvature radiation (CR) photons or photons resulting from inverse-Compton scattering of thermal photons from the PC by primary electrons. We find that, while only the younger pulsars can produce pairs through CR, nearly all known radio pulsars are capable of producing pairs through non-resonant inverse-Compton scatterings. The effect of the neutron star equations of state on the pair death lines is explored. We show that pair production is facilitated in more compact stars and more massive stars. Therefore accretion of mass by pulsars in binary systems may allow pair production in most of the millisecond pulsar population. We also find that two-photon pair production may be important in millisecond pulsars if their surface temperatures are above about three million degrees K. Pulsars that produce pairs through CR will have their primary acceleration limited by the effect of screening of the electric field. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}^{1/2}$ dependence. The acceleration voltage drop in pulsars that produce pairs only through inverse-Compton emission will not be limited by electric field screening. In this regime, the high-energy luminosity should follow a $L_{HE} \propto \dot E_{rot}$ dependence. Thus, older pulsars will have significantly lower gamma-ray luminosity.

We investigate the production of electron-positron pairs by inverse Compton scattered (ICS) photons above a pulsar polar cap (PC) and surface heating by returning positrons. This paper is a continuation of our self-consistent treatment of acceleration, pair dynamics and electric field screening above pulsar PCs. We calculate the altitude of the inverse Compton pair formation fronts, the flux of returning positrons and present the heating efficiencies and X-ray luminosities. We revise pulsar death lines implying cessation of pair formation, and present them in surface magnetic field-period space. We find that virtually all known radio pulsars are capable of producing pairs by resonant and non-resonant ICS photons radiated by particles accelerated above the PC in a pure star-centered dipole field, so that our ICS pair death line coincides with empirical radio pulsar death. Our calculations show that ICS pairs are able to screen the accelerating electric field only for high neutron star surface temperatures and magnetic fields. We argue that such screening at ICS pair fronts occurs locally, slowing but not turning off acceleration of particles until screening can occur at a curvature radiation (CR) pair front at higher altitude. In the case where no screening occurs above the PC surface, we anticipate that the pulsar gamma-ray luminosity will be a substantial fraction of its spin-down luminosity. The X-ray luminosity resulting from PC heating by ICS pair fronts is significantly lower than the PC heating luminosity from CR pair fronts, which dominates for most pulsars. PC heating from ICS pair fronts is highest in millisecond pulsars, which cannot produce CR pairs, and may account for observed thermal X-ray components in the spectra of these old pulsars.

We investigate the effect of pulsar polar cap (PC) heating produced by positrons returning from the upper pair formation front. Our calculations are based on a self-consistent treatment of the pair dynamics and the effect of electric field screening by the returning positrons. We calculate the resultant X-ray luminosities, and discuss the dependence of the PC heating efficiencies on pulsar parameters, such as characteristic spin-down age, spin period, and surface magnetic field strength. In this study we concentrate on the regime where the pairs are produced in a magnetic field by curvature photons emitted by accelerating electrons. Our theoretical results are not in conflict with the available observational X-ray data and suggest that the effect of PC heating should significantly contribute to the thermal X-ray fluxes from middle-aged and old pulsars. The implications for current and future X-ray observations of pulsars are briefly outlined.

We reinvestigate the radio pulsar ``death lines'' within the framework of two different types of polar cap acceleration models, i.e., the vacuum gap model and the space-charge-limited flow model, with either curvature radiation or inverse Compton scattering photons as the source of pairs. General relativistic frame-dragging is taken into account in both models. We find that the inverse Compton scattering induced space-charge-limited flow model can sustain strong pair production in some long-period pulsars, which allows the newly detected 8.5s pulsar PSR J2144-3933 to be radio loud, without assuming a special neutron star equation-of-state or ad hoc magnetic field configurations.

We investigate self-consistent particle acceleration near a pulsar polar cap (PC) by the electrostatic field due to the effect of inertial frame dragging. Test particles gain energy from the electric field parallel to the open magnetic field lines and lose energy by both curvature radiation (CR) and resonant and non-resonant inverse Compton scattering (ICS) with soft thermal X-rays from the neutron star (NS) surface. Gamma-rays radiated by electrons accelerated from the stellar surface produce pairs in the strong magnetic field, which screen the electric field beyond a pair formation front (PFF). Some of the created positrons can be accelerated back toward the surface and produce gamma-rays and pairs that create another PFF above the surface. We find that ICS photons control PFF formation near the surface, but due to the different angles at which the electron and positron scatter the soft photons, positron initiated cascades develop above the surface and screen the accelerating electric field. Stable acceleration from the NS surface is therefore not possible in the presence of dominant ICS energy losses. However, we find that stable acceleration zones may occur at some distance above the surface, where CR dominates the electron and positron energy losses, and there is up-down symmetry between the electron and positron PFFs. We examine the dependence of CR-controlled acceleration zone voltage, width and height above the surface on parameters of the pulsar and its soft X-ray emission. For most pulsars, we find that acceleration will start at a height of 0.5 - 1 stellar radii above the NS surface.

The study of physical processes associated with particle acceleration in the open field line region above the polar cap (PC) of an isolated neutron star (NS) plays a fundamental role in our understanding and interpretation of high-energy emission from pulsars. The systematic study of particle acceleration and formation of electron-positron pair fronts above the PCs of NSs was initiated two decades ago. The detailed analysis of these processes is now possible with the development of pair cascade codes that enables us to calculate the spectra and pulse profiles of high-energy emission from pulsars. The calculation of pair formation and gamma-ray production is being improved to include new results on the PC physics. We briefly outline the current status of the PC model for pulsar high-energy emission, focusing on some of our most recent results on the theoretical modeling of the PC acceleration and gamma-ray emission.

We model the thermal X-ray profiles of Geminga, Vela and PSR 0656+14, which have also been detected as gamma-ray pulsars, to constrain the phase space of obliquity and observer angles required to reproduce the observed X-ray pulsed fractions and pulse widths. These geometrical constraints derived from the X-ray light curves are explored for various assumptions about surface temperature distribution and flux anisotropy caused by the magnetized atmosphere. We include curved spacetime effects on photon trajectories and magnetic field. The observed gamma-ray pulse profiles are double peaked with phase separations of 0.4 - 0.5 between the peaks. Assuming that the gamma-ray profiles are due to emission in a hollow cone centered on the magnetic pole, we derive the constraints on the phase space of obliquity and observer angles, for different gamma-ray beam sizes, required to produce the observed gamma-ray peak phase separations. We compare the constraints from the X-ray emission to those derived from the observed gamma-ray pulse profiles, and find that the overlapping phase space requires both obliquity and observer angles to be smaller than 20-30 degrees, implying gamma-ray beam opening angles of at most 30-35 degrees.

We explain the millisecond variability detected by Rossi X-ray Timing Explorer (RXTE) in the X-ray emission from a number of low mass X-ray binary systems (Sco X-1, 4U1728-34, 4U1608-522, 4U1636-536, 4U0614+091, 4U1735-44, 4U1820-30, GX5-1 and etc) in terms of dynamics of the centrifugal barrier, a hot boundary region surrounding a neutron star. We demonstrate that this region may experience the relaxation oscillations, and that the displacements of a gas element both in radial and vertical directions occur at the same main frequency, of order of the local Keplerian frequency. We show the importance of the effect of a splitting of the main frequency produced by the Coriolis force in a rotating disk for the interpretation of a spacing between the QPO peaks. We estimate a magnitude of the splitting effect and present a simple formula for the whole spectrum of the split frequencies. It is interesting that the first three lowest-order overtones fall in the range of 200-1200 Hz and match the kHz-QPO frequencies observed by RXTE. Similar phenomena should also occur in Black Hole (BH) systems, but, since the QPO frequency is inversely proportional to the mass of a compact object, the frequency of the centrifugal-barrier oscillations in the BH systems should be a factor of 5-10 lower than that for the NS systems. The X-ray spectrum formed in this region is a result of upscattering of a soft radiation (from a disk and a NS surface) off relatively hot electrons in the boundary layer. We also briefly discuss some alternative QPO models, including a possibility of acoustic oscillations in the boundary layer, the proper stellar rotation, and g-mode disk oscillations.

We explain the peak spacing in the power spectra of millisecond variability detected by the Rossi X-ray Timing Explorer (RXTE) in the X-ray emission from the LMXBs (Sco X-1, 4U1728-34, 4U1608-522, 4U1636-536, etc) in terms of the rotational splitting of an intrinsic frequency (which is of order of the local Keplerian frequency) caused by an accretion disk. We calculate this effect and demonstrate that there is a striking agreement with the observational data. We show that the observed discrete frequencies ranging from 200 to 1200 Hz can be described by a whole set of overtones. For higher overtones (lower frequencies, $\lsim $ 100 Hz) the ratio between frequencies is determined by the quantum numbers alone. We suggest that a similar phenomena should be observed in Black Hole (BH) systems for which the QPO (quasi-periodic oscillation) frequency should be inversely proportional to the mass of the compact object. For BH systems the characteristic frequency of oscillations should therefore be a factor of 5-10 less than for a neutron star system. keywords

After Michel (1994) introduced a phenomenological picture of `rapid magnetization' of newly born neutron stars (NSs), Muslimov & Page (1995) suggested that the physical conditions accompanying the formation of a NS may result in the surface magnetic field of the NS having a low value (~ 10^8 - 10^9 G) while the bulk of the flux is submerged under the stellar surface. The subsequent growth of the surface magnetic field occurs owing to ohmic diffusion of a strong internal field (10^12-10^13 G) and proceeds over a time scale of 10^2-10^3 yrs, depending on the early thermal history of the NS, the initial distribution of magnetic flux and the electrical resistivity in the stellar crust, as well as on the details of the stellar structure. Referring to this suggestion, we perform numerical computations to demonstrate its possible relevance to young radio pulsars. In our calculations, we use different NS structures based on a model of dense matter presenting a phase transition to kaon condensation which implies both slow and fast neutrino cooling, depending on the stellar mass. We present sequences of secular evolution of the surface magnetic field, spin-down luminosity, spin period, `braking index', and spin-down age best matching the corresponding measured quantities for PSR 0531+21 (Crab), PSR 0540-69, and PSR 1509-58. We conclude that the effect under consideration reproduces remarkably well both the observed rotational characteristics and derived values of the surface magnetic field strength for these radio pulsars, the only ones having reliable measurements of the braking indices. In addition, our analysis suggests that the NSs in the Crab and PSR 0540-69 have experienced fast neutrino cooling while the NS in PSR 1509-58 has undergone slow cooling and hence has a mass below the critical mass.