A novel look at the pulsar force-free magnetosphere

The stationary axisymmetric force-free magnetosphere of a pulsar is considered. We present an exact dipolar solution of the pulsar equation, construct the magnetospheric model on its basis and examine its observational support. The new model has toroidal rather than common cylindrical geometry, in line with that of the plasma outflow observed directly as the pulsar wind nebula at much larger spatial scale. In its new configuration, the axisymmetric magnetosphere consumes the neutron star rotational energy much more efficiently, implying re-estimation of the stellar magnetic field, \(B_{\mathrm{new}}^{0}=3.3\times10^{-4}B/P\), where \(P\) is the pulsar period. Then the 7-order scatter of the magnetic field derived from the rotational characteristics of the pulsars observed appears consistent with the \(\cot\chi\)-law, where \(\chi\) is a random quantity uniformly distributed in the interval \([0,\pi/2]\). Our result is suggestive of a unique actual magnetic field strength of the neutron stars along with a random angle between the magnetic and rotational axes and gives insight into the neutron star unification on the geometrical basis.     

The magnetohydrodynamic rotational model of supernova explosion

The calculations of supernova explosion are made, using the one-dimensional nonstationary equations of magnetic hydrodynamics for the case of cylindrical symmetry. The energy source is supposed to be the rotational energy of the system (the neutron star in the centre and the surrounding envelope). The magnetic field plays the role of a mechanism of the transfer of rotational momentum. The calculations show that the envelope split up during the dynamical evolution of the system, the main part of the envelope joins the neutron star and becomes uniformly rotating with it, the outer part of the envelope (10% mass) expands with large velocity, carrying out a considerable part of rotational energy and rotational momentum.These results correspond qualitatively with the observational picture of supernovae explosions.     

Stellar pulsation: (I) analysis of stability of envelope

The life-time of the star on AGB is approximately 6 × 10⁴ yr. We divide it into front half and back half of AGB (including to optical Mira variable and OH/IR star) according to their evolution character. The observations show that the star has non-pulsation, but constant mass loss rate ( 5 × 10^–7 M yr^–1) on front half of AGB. Its circumstellar envelope is formed. When the star has pulsation on back half of AGB, its mass loss rate is relative with time, and increases gradually. In this time the star is on the stage of optical Mira variable. When the mass loss rate reaches the value of 3 × 10^–6 M yr^–1, the star evoluted from the stage of optical variable into the stage of radio bright OH/IR star. On the end of AGB the mass loss rate reaches 10^–4 M yr^–1. (Band and Habing 1983, Hermen and Habing 1985).     

Rapidly rotating supermassive stars in the first post-Newtonian approximation to general relativity

The structure and stability of rapidly uniformly rotating supermassive stars is investigated using the full post-Newtonian equations of hydrodynamics. The standard model of a supermassive star, a polytrope of index three, is adopted. All rotation terms up to and including those of order ⁴, where is the angular velocity, are retained. The effects of rotation and post-Newtonian gravitation on the classical configuration are explicitly evaluated and shown to be very small. The dynamical stability of the model is treated by using the binding energy approach. The most massive objects are found to be dynamically unstable when =1/c ².p _c / _c 2.2 × 10^–3, wherep _c and _c are the central pressure and density, respectively. Hence, the higher-order terms considered in this analysis do not appreciably alter the previously known stability limits.The maximum mass that can be stabilized by uniform rotation in the hydrogen-burning phase is found to be 2.9×10⁶ M , whereM is the solar mass. The corresponding nuclear-generated luminosity of 6×10⁴⁴ erg/sec^–1 is too small for the model to be applicable to the quasi-stellar objects. The maximum kinetic energy of a uniformly rotating supermassive star is found to be 3×10^–5 Mc ², whereM is the mass of the star. Masses in excess of 10¹⁰ M are required if an adequate store of kinetic energy is to be made available to a pulsar like QSO. However such large masses have rotation periods in excess of 100 yr and thus could not account for any short term periodic variability. It is concluded then that the uniformly rotating supermassive star does not provide a suitable base for a model of a QSO.     

Pulsar radio emission

A new version of the theory of pulsar radio emission is developed for the case of a coaxial rotator. It is based on the electric field that we established [G. S. Sahakian, Astrofizika, 37, 97 (1994)] for the radiation channel (the channel of open magnetic field lines) and on convenient approximations for the electron energy obtained in [G. S. Sahakian and É. S. Chubarian, Astrofizika, 37, 255 (1994)]. It is shown that, owing to the emission of photons of curvature radiation by particles, e e+_c', and photon annihilation, _c e⁺e^– in the lower part of the radiation channel, a special region (the magnetic funnel) is formed in which vigorous cascade multiplication of particles occurs. The height of the magnetic funnel is h 6R^0.2, where R is the radius of the neutron star and is its angular rotation rate. As a result of supersaturation of the plasma density in the magnetic funnel, a discharge occurs after each time intervalt5·10^–7–0.8B ₁₂ ^–1.4 R ₆ ^–0.2 , i.e., the longitudinal electric field disappears (B is the magnetic induction in the star). During the active radiative processes in the magnetic funnel, two main fluxes of particles with high ultrarelativistic energies are formed: an upward flux of electrons and a positron flux falling onto the star's magnetic cap. These fluxes are accompanied by narrow strips of positron and electron fluxes, respectively, of considerably lower energy, which are fairly powerful, coherent radio sources. The pulsar's radio luminosity is calculated to be L7.4·10^223.8 ₃₀ ³ R ₆ ^–2 erg/sec, where =BR 3/2 is the star's magnetic moment. Comparing this result with observations, we conclude that the magnetic moment and hence the mass of the neutron star evidently must be considerably smaller, on the average, for fast pulsars than for slow ones. It is shown that the magnetic moment of the neutron star can be determined from the intervals between micropulses in the pulse profiles. The problem of the origin of the macrostructure of the radio pulse is discussed.Translated from Astrofizika, Vol. 38, No. 1, pp. 141–185, January – March, 1995.     

Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

Search for the sources of cosmic rays with energies above 10<Superscript>20</Superscript> eV

The sources of ultrahigh energy cosmic rays (UHECRs, E >10¹⁸ eV) are still unknown, mainly due to the loss of the direction to the source after the deflection of cosmic rays’ (CRs) trajectories in the galactic and extragalactic magnetic fields. With the increase in CR energy (rigidity), the influence of the magnetic field weakens; therefore, the most promising approach is to search for the sources of events with the highest energy. In our work, we expand the existing UHECR (E > 10²⁰ eV) sample from 33 to 42 events by calibrating the AUGER events. The sample is characterized by the presence of an event triplet in a circle of radius 3°. The highest-energy event is still the shower (E = 3.2 × 10²⁰ eV) detected with the Fly’s Eye fluorescent detector (FE-event) in 1993. The possible sources of the triplet and the FE-event are analyzed. Taking into account the deflection of CR trajectories in the extragalactic and galactic magnetic fields, it is shown that transient sources of the FE-event and the triplet may be galaxies with active star formation, where CRs are accelerated by newborn millisecond pulsars. Among the galactic sources, the potential candidates are young pulsars that might have had millisecond periods at birth and giant magnetar flares.     

Two X-ray sources model of Cyg X-3

A model of Cyg X-3, as a binary cocooned star system with two sources of X-rays, one above the polar caps of the neutron star — the usual pulsar radiation — and the other around the equatorial plane of the magneto-bounding surface formed due to the interaction of the infalling plasma and the magnetic field of the neutron star, is made. The X-ray, -ray, and IR radiation light curves are considered from the shadow effect. An upper limit on the mass of the neutron star is estimated from the consideration of periodic derivative purely due to mass loss. A comparison is made with the results of Elsneret al. (1980) and Ghoshet al. (1981), which they derived from the consideration of period derivative purely from apsidal motion.     

Relativistic supermassive stars

A plausible scenario for the formation of a stable supermassive star in the relativistic regime has been discussed. The onset of the negativity of the `distribution function' in the stable sequences of the star clusters [the stability of star clusters is assured by using the variational method (Chandrasekhar, 1964a,b) for equivalent gas spheres] described by Tolman's type VII solution with vanishing surface density has been regarded as an indication of the conversion of the cluster structure into a supermassive star. For the critical values of the `adiabatic index', (4/3) < _crit (5/3) (forwhich a supermassive star represent neutrally stable system), the mass, and the size of this object comes out to be 6.87 × 10⁷ M M 1.7 × 10⁹M, and 2.74× 10¹⁴ cm a 1.43× 10¹⁵cm respectively, for the central temperature,T₀ = 6× 10⁷°K, which is sufficient for the release ofnuclear energy. The total energy released during their evolution rangesfrom 2.46× 10⁶⁰ - 3.18× 10⁶² erg, which is sufficient to power these objects at least for a period of 10⁶ - 10⁷years. These figures agrees quite well with those cited for Quasi Stellar Objects (QSOs) in the literature.     

Generation of γ-bursts by old magnetic neutron stars

A model of -bursts is considered that treats the flares of neutron stars as a result of convectiveoscillation instability associated with the stars having strong internal magnetic fields ( 10¹³ to 10¹⁴ G). In the context of this model only sufficiently old (10⁴ to 10⁷ yr), drastically cooled-down neutron stars may be sources of -bursts. The paper shows that major characteristics of a -burster in the Supernova N 49 remnant (energy release during burst up to 10⁴⁴ erg, age 10⁴ yr, burst-to-burst interval (I to 3)×10⁶s; rotation period P=8 s) may be explained under the assumption that the mass of the neutron star is about 0.14M _· while its mean magnetic field strength is 1.5×10¹⁴ G abd 10¹³ G within the star and on its surface, respectively. The observational tests of the model discussed conclude the paper.     

Mechanisms of optical,X-ray andγ-radiation from Crab pulsar

The synchrotron mechanism of radiation from the Crab pulsar has been investigated on the assumption that the mechanism acts in a source moving with relativistic velocity round a neutron star. A detailed matching has been made of the theoretical spectra of synchrotron radiation from relativistic electrons with the results of measurements of the radiation flux from the Crab pulsar in the infrared, optical and X-ray ranges. The parameters of the radiating region (intensity of the magnetic field, source dimensions, density and lifetime of radiating electrons) have been found. They are expressed through the ratio of the energy density of the magnetic field in the source to that of radiating electrons. The level of Compton-radiation in this region is estimated. Possible values of at which the level will correspond to the available results of measurements of the-radiation flux from the Crab pulsar are given. An estimate is presented for the surface magnetic field of the neutron star which does not contradict those obtained from considerations of the magnetic flux conservation when compressing the object up to the neutron star dimensions.     

Flux-Vortex Pinning and Neutron Star Evolution

G. Srinivasan et al. (1990) proposed a simple and elegant explanation for the reduction of the neutron star magnetic dipole moment during binary evolution leading to low mass X-ray binaries and eventually to millisecond pulsars: Quantized vortex lines in the neutron star core superfluid will pin against the quantized flux lines of the proton superconductor. As the neutron star spins down in the wind accretion phase of binary evolution, outward motion of vortex lines will reduce the dipole magnetic moment in proportion to the rotation rate. The presence of a toroidal array of flux lines makes this mechanism inevitable and independent of the angle between the rotation and magnetic axes. The incompressibility of the flux-line array (Abrikosov lattice) determines the epoch when the mechanism will be effective throughout the neutron star. Flux vortex pinning will not be effective during the initial young radio pulsar phase. It will, however, be effective and reduce the dipole moment in proportion with the rotation rate during the epoch of spindown by wind accretion as proposed by Srinivasan et al. The mechanism operates also in the presence of vortex creep.     

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³.     

Microwave radiation from magnetic stars

Cyclotron microwave emission from magnetic stars is considered, assuming that they have coronae with the temperatureT10⁷ K and the emission measureEM10⁵⁴ cm^–3. It has been shown that the cyclotron radiation from a star with a dipole magnetic field has a specific spectrum with a maximum in the frequency rangesv _o/2 >v >sv _o/2 (s being the number of cyclotron harmonic, andv _o the gyrofrequency corresponding to the polar magnetic field) and radiation flux decreasing towards lower frequencies asv _4/3. The frequency of the spectrum maximum depends on the angle between the line-of-sight and the magnetic axis of the star. The observed radiation from a rotating magnetic star can be modulated with a modulation depth of about 0.2 at frequencies near maximum. The radiation is partially circularly-polarized in the sense of an extraordinary mode. The degree of polarization is almost constant at frequenciesv >sv _o/2 and increases with frequency atv >sv _o/2. The estimation of cyclotron radio fluxes of the nearest magnetic stars shows that they are observable in microwaves by means of modern radio astronomy.     

NuSTAR observations of the X-ray pulsar LMC X-4: A constraint on the magnetic field and tomography of the system in the fluorescent iron line

We present the results of the spectral and timing analysis of the X-ray pulsar LMC X-4 based on data from the NuSTAR observatory in the broad X-ray energy range 3–79 keV. Along with a detailed analysis of the source’s averaged spectrum, high-precision spectra corresponding to different phases of the neutron star spin cycle have been obtained for the first time. The Comptonization model is shown to describe best the source’s spectrum, and the evolution of its parameters as a function of the pulse phase has been traced. For all spectra (the averaged and phase-resolved ones) in the energy range 5–55 keV we have searched for the cyclotron absorption line. The derived upper limit on the optical depth of the cyclotron line τ ~ 0.15 (3σ) points to the absence of this feature in the given energy range, which provides a constraint on the magnetic field of the neutron star: B <3 × 10¹¹ or >6.5 × 10¹² G. The latter constraint is consistent with the magnetic field estimate obtained by analyzing the pulsar’s power spectrum, B ? 3 × 10¹³ G. Based on our analysis of the phase-resolved spectra, we have determined the delay between the emission peaks and the equivalent width of the fluorescent iron line. This delay depends on the orbital phase and is apparently associated with the travel time of photons between the emitting regions in the vicinity of the neutron star and the region where the flux is reflected (presumably in the inflowing stream or at the place of interaction between the stream and the outer edge of the accretion disk).     

Vector magnetic field evolution,energy storage,and associated photospheric velocity shear within a flare-productive active region

The evolution of vector photospheric magnetic fields has been studied in concert with photospheric spot motions for a flare-productive active region. Over a three-day period (5–7 April, 1980), sheared photospheric velocity fields inferred from spot motions are compared both with changes in the orientation of transverse magnetic fields and with the flare history of the region. Rapid spot motions and high inferred velocity shear coincide with increased field alignment along the B _L= 0 line and with increased flare activity; a later decrease in velocity shear precedes a more relaxed magnetic configuration and decrease in flare activity. Crude energy estimates show that magnetic reconfiguration produced by the relative velocities of the spots could cause storage of 10³² erg day^–1, while the flares occurring during this time expended 10³¹ erg day^–1.Maps of vertical current density suggest that parallel (as contrasted with antiparallel) currents flow along the stressed magnetic loops. For the active region, a constant-, force-free magnetic field (J = B) at the photosphere is ruled out by the observations.Presently located at NASA/MSFC, Huntsville, Ala. 35812, U.S.A.     

Magnetic Model of HD 2453

A model is constructed for the magnetic field of the star HD 2453, which has a very long rotation period (P=521^d). It is found that the structure of the field corresponds to the model of a dipole shifted by r=0.09R from the center. The angle of inclination of the axis of the dipole to the axis of rotation, =5°; that is, the star is viewed almost from its equator of rotation and magnetic equator. This explains the low amplitude of the phase dependence of the magnetic field, Be(P), and the low amplitude of the photometric variability. The field at the magnetic poles is equal to Bp=+4400 and -7660 G. The magnetic field parameters turn out to be close to those obtained by Landstreet and Mathys assuming a dipole-quadrupole-octupole model. A Mercator map of the magnetic field distribution of HD 2453 is produced.     

Structure of the Magnetic Field of a Neutron Star

The magnetic field distribution in the superfluid, spherical, hadronic core of a rotating neutron star, which consists of vortex and vortex-free zones, is investigated. Due to the effect of entrainment of superconducting protons by rotating superfluid neutrons, a nonuniform magnetic field, the average value of which is constant, is formed in the vortex zone of the neutron star, directed parallel to the star's axis of rotation. It is shown that at the stellar surface, near the equatorial plane, there is a vortex-free zone of macroscopic size in which there is no magnetic field. The magnetic field near the boundaries of the vortex-free zone falls off exponentially with depth into the interior of this zone. This result essentially alters earlier concepts about the magnetic field distribution in the superfluid hadronic core of a neutron star. Outside the hadronic core the magnetic field has a dipole character with a magnetic moment on the order of 10³⁰ g×cm³.     

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.     

Plasmon neutrinos emission in a strong magnetic field

The decay of a longitudinal plasmon into two neutrinos is studied in the presence of a strong magnetic field. Contrary to the transverse case, for longitudinal plasmons the existence of a new mode, entirely dependent on the magnetic field, greatly enhances the energy loss at high densities. Denoting byQ _HandQ ₀the neutrino energy losses with and without magnetic field respectively, the situation is as follows: atH10¹¹ G andT10⁹K,Q ₀10⁵ Q _Hfor <10¹¹g cm^–3, WhileQ _H10¹⁰ Q ₀for >10¹¹g cm^–3. A second physically interesting feature is the anisotropic character of the energy loss which is highly peaked along the field lines, giving rise to a shorter cooling time in that direction than in any other.