Inverse compton emission of gamma rays near the pulsar surface

The physical conditions near pulsar surface that might give rise to gamma ray emission from Crab and Vela pulsars are not yet well understood. Here I suggest that, in the context of the vacuum discharge mechanism proposed by Ruderman and Sutherland (1975), gamma rays are produced by inverse Compton scattering of secondary electrons with the thermal radiation of the star surface as well as for curvature and synchrotron radiation. It is found that inverse Compton scattering is relevant if the neutron star surface temperature is greater than 10⁶K or if the polar cap temperature is of the order of 5×10⁶K. Inverse Compton scattering in anisotropic photon fields and Klein-Nishina regime is here carefully considered.     

Hydrogen molecular ion in the magnetic field of a neutron star

A two-dimensional potential energy surface of an H ₂ ⁺ molecular ion is calculated for the case of the strong magnetic field of the neutron starB=10¹¹–10¹³ G. It is shown that the dependence of the potential energy from the angle between the magnetic field direction and the internuclear axis becomes very sharp as the magnetic field increases. The obtained potential energy surfaces can be used for studying the vibrational-rotational structure of the H ₂ ⁺ spectrum in a strong magnetic field and the development of the observational methods for the determination of the magnetic field of a neutron star.     

Influence of MHD waves on thermonuclear burning in surface layers of neutron stars: The nature of gamma-bursts

A possibility for gamma-ray bursts to arise due to thermonuclear flashes in the surface layers of accreting neutron stars is discussed. The principal difference of the sources of gamma-ray bursts from bursters is supposed to result from the existence of strong magnetic fields (10¹²–10¹³G) on the neutron star surface. It is shown that the thermonuclear energy released may be rapidly and effectively transported to the outer layers by MHD waves (in particular, by Alfvén waves). A very short growth time and rapid variations of some gamma-ray bursts may be easily explained in this case.     

Tunnel effect in molecules in strong magnetic fields of neutron stars

It is shown that, in the strong magnetic field of the neutron starB=10¹²–10¹³ G, the probability of the tunnel effect in the molecules increases significantly. It is quite probable that this effect can catalyze nuclear reactions at the neutron star surface.     

The drift model of “magnetars”

It is shown that the drift waves near the light cylinder can cause the modulation of emission with periods of order several seconds. These periods explain the intervals between successive pulses observed in AXPs, SGRs and radio pulsars with long periods. The model under consideration gives the possibility to calculate real rotation periods P of host neutron stars. It is shown that P≤1 s for the investigated objects. The magnetic fields at the surface of the neutron star are of order 10¹¹–10¹³ G and equal to the fields usual for the known radio pulsars.      

Hydrogen molecular ion in the magnetic field of a neutron star

The paper discusses the structure and energies of rotational-vibrational levels of a molecular ion H ₂ ⁺ in a strong magnetic field typical of neutron stars,B=10¹²–10¹³ G. The study is based on the calculations of the potential energy surface of this molecular ion presented in the paper by Khersonskii (1984a). The energies of the rovibrational levels are calculated with the aid of the perturbation theory. The number of levels in a H ₂ ⁺ potential hole is determined at different magnetic intensities. In particular, it is ascertained that the number of levels decreases as the strength of the magnetic field increases. The effect of nuclear spins on the structure of the rotational levels is considered.     

Relaxation time of the superfluid core of a neutron star

The superconducting proton condensate in the “npe” phase of a neutron star is considered. It is shown to be a type II superconductor in the outer layer of the “npe” phase and a type I superconductor in the inner layer. Relaxation times are found for elastic scattering of normal relativistic electrons from the magnetic fields of proton vortex clusters in the case of a type II superconductor and elastic scattering from the magnetic field at the center of a neutron vortex in the case of a superconductor of the first kind. The dynamical relaxation times obtained for the angular velocity of the pulsar PSR 0833—45 vary, as a function of the density of the layers taking part in the relaxation process, within a fairly wide range: from several hours to l0⁹ years. This means that the characteristic times of variation of pulsar angular velocity may be observed to lie in the indicated time range. Translated from Astrofizika, Vol. 40, No. 4, op. 497–506, October–December, 1997.     

One-photon and two-photon annihilation of positronium in strong magnetic fields

In the presence of a strong magnetic field (such as those believed to be characteristic of neutron stars:B-10¹² Gauss) positronium may annihilate through the emission of a single photon, the magnetic field providing the photon momentum. We report on calculations of the one-photon and two-photon annihilation rates for the ground state of positronium, for magnetic fields in the range (1–44)×10¹² Gauss, and give, in the two-photon case, the minimum energy half-width of the emission line due to the momentum contributions from the magnetic field. We find that unless neutron stars have magnetic fields in excess of 10¹³ Gauss, it is unlikely that the one-photon process will be observable.Research supported in part by the National Research Council of Canada.     

Evolution of superhigh magnetic fields of magnetars

In this paper, we consider the effect of Landau levels on the decay of superhigh magnetic fields of magnetars. Applying ³ P ₂ anisotropic neutron superfluid theory yield a second-order differential equation for a superhigh magnetic field B and its evolutionary timescale t. The superhigh magnetic fields may evolve on timescales ∼(10⁶–10⁷) yrs for common magnetars. According to our model, the activity of a magnetar may originate from instability caused by the high electron Fermi energy.     

Spectral Constraints for Millisecond Pulsars Due to General Relativistic Frame Dragging

We develop a numerical code for simulating the magnetospheres of millisecond pulsars, which are expected to have unscreened electric potentials due to the lack of magnetic pair production. We incorporate General Relativistic (GR) expressions for the electric field and charge density and include curvature radiation (CR) due to primary electrons accelerated above the stellar surface, whereas inverse Compton scattering (ICS) of thermal X-ray photons by these electrons are neglected as a second-order effect. We apply the model to PSR J0437-4715, a prime candidate for testing the GR-Electrodynamic theory, and find that the curvature radiation spectrum cuts off at energies below 15 GeV, which are well below the threshold of the H.E.S.S. telescope, whereas Classical Electrodynamics predict a much higher cutoff near 100 GeV, which should be visible for H.E.S.S., if standard assumed Classical Electrodynamics apply. GR theory also predicts a relatively narrow pulse (2φ _L ∼ 0.2 phase width) centered on the magnetic axis, which sets the beaming solid angle to ∼0.5 sr per polar cap (PC) for a magnetic inclination angle of 35^∘ relative to the spin axis, given an observer which sweeps close to the magnetic axis. We also find that EGRET observations above 100 MeV of this pulsar constrain the polar magnetic field strength to B _pc < 4× 10⁸ G for a pulsar radius of 10 km and moment of inertia of 10⁴⁵ g cm². The field strength constraint becomes even tighter for a larger radius and moment of inertia. Furthermore, a reanalysis of the full EGRET data set of this pulsar, assuming the predicted pulse shape and position, should lead to even tighter constraints on neutron star and GR parameters, up to the point where the GR-derived potential and polar cap current may be questioned.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.     

XMM-Newton observations of the isolated neutron star 1RXS J214303.7+065419/RBS1774

We report XMM-Newton observations of the isolated neutron star RBS1774 and confirm its membership as an XDINS. The X-ray spectrum is best fit with an absorbed blackbody with temperature kT=101 eV and absorption edge at 0.7 keV. No power law component is required. An absorption feature in the RGS data at 0.4 keV is not evident in the EPIC data, but it is not possible to resolve this inconsistency. The star is not seen in the UV OM data to m _AB ∼21. There is a sinusoidal variation in the X-ray flux at a period of 9.437 s with an amplitude of 4%. The age as determined from cooling and magnetic field decay arguments is 10⁵–10⁶ yr for a neutron star mass of 1.35–1.5 M_⊙.      

Vertical structure of the outer accretion disk in persistent low-mass X-ray binaries

We have investigated the influence of X-ray irradiation on the vertical structure of the outer accretion disk in low-mass X-ray binaries by performing a self-consistent calculation of the vertical structure and X-ray radiation transfer in the disk. Penetrating deep into the disk, the field of scattered X-ray photons with energy E ≳ 10 keV exerts a significant influence on the vertical structure of the accretion disk at a distance R ≳ 10¹⁰ cm from the neutron star. At a distance R ∼ 10¹¹ cm, where the total surface density in the disk reaches Σ₀ ∼ 20 g cm⁻², X-ray heating affects all layers of an optically thick disk. The X-ray heating effect is enhanced significantly in the presence of an extended atmospheric layer with a temperature T _atm ≈ (2–3) × 10⁶ K above the accretion disk. We have derived simple analytic formulas for the disk heating by scattered X-ray photons using an approximate solution of the transfer equation by the Sobolev method. This approximation has a ≲10% accuracy in the range of X-ray photon energies E < 20 keV.     

Ohm’s law for plasma in general relativity and Cowling’s theorem

The general-relativistic Ohm’s law for a two-component plasma which includes the gravitomagnetic force terms even in the case of quasi-neutrality has been derived. The equations that describe the electromagnetic processes in a plasma surrounding a neutron star are obtained by using the general relativistic form of Maxwell equations in a geometry of slow rotating gravitational object. In addition to the general-relativistic effect first discussed by Khanna and Camenzind (Astron. Astrophys. 307:665, 1996) we predict a mechanism of the generation of azimuthal current under the general relativistic effect of dragging of inertial frames on radial current in a plasma around neutron star. The azimuthal current being proportional to the angular velocity ω of the dragging of inertial frames can give valuable contribution on the evolution of the stellar magnetic field if ω exceeds 2.7×10¹⁷(n/σ) s⁻¹ (n is the number density of the charged particles, σ is the conductivity of plasma). Thus in general relativity a rotating neutron star, embedded in plasma, can in principle generate axial-symmetric magnetic fields even in axisymmetry. However, classical Cowling’s antidynamo theorem, according to which a stationary axial-symmetric magnetic field can not be sustained against ohmic diffusion, has to be hold in the general-relativistic case for the typical plasma being responsible for the rotating neutron star.     

On the differential rotation of the polar caps of neutron stars

The model of a magnetized rotating neutron star with an electric current in the region of its fluid polar magnetic caps is considered. The presence of an electric current leads to differential rotation of the magnetic caps. The rotation structure is determined by the electric current density distribution over the surface. In the simplest axisymmetric configuration, the current flows in one direction near the polar cap center and in the opposite direction in the outer ring (the total current is zero for the neutron star charge conservation). In this case, two rings with opposite directions of rotation appear on the neutron star surface, with the inner ring always lagging behind the star’s main rotation. The differential rotation velocity is directly proportional to the electric current density gradient along the polar cap radius. At a width of the region of change in the electric current from 1 to 10² cm and a period ～1 s and a magnetic field B ～ 10¹² G typical of radio pulsars, the linear differential rotation velocity is ～10^?2–10^?4 cm s^?1 (corresponding to a revolution time of ～0.1–10 yr).     

Nuclear reactions on accreting neutron star surface

The production of X-rays and gamma-rays in bursts is believed to be due to the rapid burning of matter accreted onto a neutron star surface from its companion, most probably a giant star. The accreted matter consists mainly of hydrogen and helium and a very small amount of heavy elements. Due to the infall of matter, the temperature at the bottom layers is raised to a value of the order of 10⁸ K. The neutron star surface density is>10⁷ g cm^–3. As hydrogen burning is a slow process under any temperature and density conditions, we consider the helium-burning reactions as the source of gamma-rays in the neutron star surface. Under high-density conditions the ordinary laboratory reaction rates should become modified. At high-density conditions, the strong screening effect due to the polarising cloud of electrons around the ions become important and enhances the reaction rates considerably. The helium-burning reactions are calculated under such conditions. The abundances of helium-burning products such as¹²C, 1¹⁶O, and²⁰Ne, etc., are computed. Under high-density and temperature conditions carbon is found to be more abundant than oxygen. Neon is completely absent in almost all the relevant physical conditions in which a strong screening effect is operative. It is suggested that explosive burning of accreted helium of 10^–13 M will account for the observed energy of gamm-ray burst.     

MHD Simulations of the ISM: The Importance of the Galactic Magnetic Field on the ISM “Phases”

We have carried out 1.25 pc resolution MHD simulations of the ISM, on a Cartesian grid of 0 ≤ (x, y) ≤ 1 kpc size in the galactic plane and ?10 ≤ z ≤ 10 kpc into the halo, thus being able to fully trace the time-dependent evolution of the galactic fountain. The simulations show that large scale gas streams emerge, driven by SN explosions, which are responsible for the formation and destruction of shocked compressed layers. The shocked gas can have densities as high as 800 cm^?3 and lifetimes up to 15 Myr. The cold gas is distributed into filaments which tend to show a preferred orientation due to the anisotropy of the flow induced by the galactic magnetic field. Ram pressure dominates the flow in the unstable branch 10² < T ≤ 10^3.9 K, whereas for T ≤ 100 K (stable branch) magnetic pressure takes over. Near supernovae thermal and ram pressures determine the dynamics of the flow. Up to 80% of the mass in the disk is concentrated in the thermally unstable regime 10² < T ≤ 10^3.9 K with ～30% of the disk mass enclosed in the T ≤ 10³ K gas. The hot gas in contrast is controlled by the thermal pressure, since magnetic field lines are swept towards the dense compressed walls.     

On the slow-down of the millisecond pulsar PSR 1937+214

There are indications that less than 10^–3 of the spin-down energy of the millisecond pulsar PSR 1937+214 emerges as electromagnetic radiation. The implications of this result are discussed. The surface magnetic field would then be 10⁷ G, making the pulsar optically undetectable, and casting aspersions on the accretion disc spin-up neutron star models for the pulsar. The pulsar should have an equatorial ellipticity 10^–9, which can be accounted for if the equatorial magnetic field departs from axisymmetry by one part in 10³.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

Tenma observation of the X-ray pulsar 4U1626-67

Energy-dependent pulse profiles from the X-ray pulsar 4U1626-67, observed with Tenma in May 1983, were studied using an anisotropic radiation transfer model in an accretion column on a strongly magnetized neutron star. The result indicates that the strength of the magnetic field is about 8×10¹² G at the surface of the neutron star.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.