Cosmic rays from collapsing magnetized stars

We examine the acceleration of cosmic rays in the magnetospheres of collapsing stars with initial dipole magnetic fields and various initial energy distributions of charged particles in their magnetospheres (the exponential, relativistic Maxwellian, and Boltzmann distributions were considered). When a magnetized star contracts at the gravitational collapse stage, its magnetic field grows considerably. Such a variable magnetic field generates an eddy electric field. Our calculations suggest that this electric field can accelerate charged particles to relativistic energies. In this way collapsing stars can be sources of high-energy cosmic rays in our Galaxy as well as in other galaxies.     

Active galactic nuclei and pulsars as cosmic ray sources

Relativistic e^± particles and cosmic rays are accelerated in the magnetospheres of supermassive black holes and neutron stars. The possibility of synchrotron radiation with extremely high intensity inside the deepest regions of magnetospheres is investigated. Very high brightness temperatures are expected for such radiation by relativistic protons, which can be made even higher in the presence of non-stationary conditions, Doppler boosting and coherent processes. The main parameters for models of such high-brightness-temperature radiation are determined. Two types of active galactic nuclei (AGNs) are expected. One type is associated with the acceleration and ejection of relativistic e^± particles only (probably non-IDV sources and FR-I radio galaxies). The second type of AGN is also associated with e^± acceleration, but is dominated by the contribution of relativistic protons (probably IDV sources and FR-II radio galaxies). Analogous objects for pulsars are plerion and shell supernova remnants with neutron stars or pulsars without synchrotron nebulae, respectively.     

Energization of charged particles in planetary magnetospheres

A model is presented to describe the energization of charged particles in planetary magnetospheres. The model is based on the stochastic acceleration produced by a random electric field that is induced by the magnetic field fluctuations measured within the magnetospheres. The stochastic behavior of the electric field is simulated through a Monte Carlo method. We solve the equation of motion for a single charged particle—which comprises the stochastic acceleration due to the stochastic electric field, the Lorentz acceleration (containing the local magnetic field and the corotational electric field) and the gravitational planetary acceleration of the particle—under several initial conditions. The initial conditions include the ion species and the velocity distribution of the particles which depends on the sources they come from (solar wind, ionospheres, rings and satellites). We applied this model to Saturn’s inner magnetosphere using a sample of particles (H⁺, H₂O⁺, N⁺, O⁺ and OH⁺) initially located on Saturn’s north pole, above the C-Ring, on the south pole of Enceladus, in the north pole of Dione and above the E-Ring. The results show that the particles tend to increase the value of their energy with time reaching several eV in a few seconds and the large energization is observed far from the planet. We can distinguish three main energization regions within Saturn’s inner magnetosphere: minimum (Saturn’s ionosphere), intermediate (Dione) and high-energy (Enceladus and the E-ring). The resulting energy spectrum follows a power-law distribution (>1 keV), a logistic, an exponential decay or an asymmetric sigmoidal (<1 keV).     

Numerical construction of magnetosphere with relativistic two-fluid plasma flows

We present a numerical model in which a cold pair plasma is ejected with relativistic speed through a polar cap region and flows almost radially outside the light cylinder. Stationary axisymmetric structures of electromagnetic fields and plasma flows are self-consistently calculated. In our model, motions of positively and negatively charged particles are assumed to be determined by electromagnetic forces and inertial terms, without pair creation and annihilation or radiation loss. The global electromagnetic fields are calculated by the Maxwell's equations for the plasma density and velocity, without using ideal magnetohydrodynamic condition. Numerical result demonstrates the acceleration and deceleration of plasma due to parallel component of the electric fields. Numerical model is successfully constructed for weak magnetic fields or highly relativistic fluid velocity, i.e. kinetic energy dominated outflow. It is found that appropriate choices of boundary conditions and plasma injection model at the polar cap should be explored in order to extend present method to more realistic pulsar magnetosphere, in which the Poynting flux is dominated.     

The Converter Mechanism of Particle Acceleration and Its Applications to the Unidentified Egret Sources

We discuss the properties of gamma-ray radiation accompanying the acceleration of cosmic rays via the converter mechanism. The mechanism exploits multiple photon-induced conversions of high-energy particles from charged into neutral state (namely, protons to neutrons and electrons to photons) and back. Because a particle in the neutral state can freely cross the magnetic field lines, this allows to avoid both particle losses downstream and reduction in the energy gain factor, which normally takes place due to highly collimated distribution of accelerated particles. The converter mechanism efficiently operates in relativistic outflows under the conditions typical for Active Galactic Nuclei, Gamma-Ray Bursts, and microquasars, where it outperforms the standard diffusive shock acceleration. The accompanying radiation has a number of distinctive features, such as an increase of the maximum energy of synchrotron photons and peculiar radiation beam-pattern, whose opening angle is much wider at larger photon energies. This provides an opportunity to observe off-axis relativistic jets in GeV–TeV energy range. One of the implications is the possibility to explain high-latitude unidentified EGRET sources as off-axis but otherwise typical relativistic-jet sources, such as blazars.     

The X-Ray Emission of the Centaurus A Jet

The extended nonthermal X-ray emission of extragalactic jets like Centaurus A can only be explained by in situ particle acceleration. The only energy source in the entire jet region is the magnetic field. Magnetic reconnection can convert the free energy stored in the helical configuration to particle kinetic energy. In the collisionless magnetized jet plasma, the inertia-driven reconnection is operating in a highly filamentary magnetic flux rope, and this results in a continuously charged particle acceleration. The synchrotron radiation of these particles can cause the observed X-ray emission in Centaurus A.     

Stochastic Acceleration of Energetic Particles in the Magnetosphere of Saturn

Voyager's plasma probe observations suggest that there are at least three fundamentally different plasma regimes in Saturn: the hot outer magnetosphere, the extended plasma sheet, and the inner plasma torus. At the outer regions of the inner torus some ions have been accelerated to reach energies of the order of 43 keV. We develop a model that calculates the acceleration of charged particles in the Saturn's magnetosphere. We propose that the stochastic electric field associated to the observed magnetic field fluctuations is responsible of such acceleration. A random electric field is derived from the fluctuating magnetic field – via a Monte Carlo simulation – which then is applied to the momentum equation of charged particles seeded in the magnetosphere. Taking different initial conditions, like the source of charged particles and the distribution function of their velocities, we find that particles injected with very low energies ranging from 0.129 eV to 5.659 keV can be strongly accelerated to reach much higher energies ranging from 22.220 eV to 9.711 keV as a result of 125,000 hitting events (the latter are used in the numerical code to produce the particle acceleration over a predetermined distance).     

Stability of Accretion Models

A critical analysis of standard accretion models is presented. We consider the stability of models in the theories of disc accretion onto black holes and spherical/disc accretion onto a magnetosphere. We take into account realistic physics processes and geometry (inner magnetic field in the accreted plasma, finite conductivity, finite length of the field lines, finite rotation of the accreted object, and magnetic shear on the boundary between the magnetosphere and accreted plasma). The influence of these factors leads to radical changes of both the accretion as whole and the energy release in the accreting system. Strong current-sheet and Z-pinch-like structures should arise over the polar region of the accreting object. Particle acceleration in the electric fields of current discharges in these regions may be a source of efficient conversion of energy into nonthermal particles and of the emission observed from many accreting objects.     

Scale transformations, tree-level perturbation theory and the cosmological matter bispectrum

Soon after the discovery of radio pulsars in 1967, the pulsars are identified as strongly magnetic (typically 10¹² G) rapidly rotating (∼10²− 0.1 Hz) neutron stars. However, the mechanism of particle acceleration in the pulsar magnetosphere has been a longstanding problem. The central problem is why the rotation power manifests itself in both gamma-ray beams and a highly relativistic wind of electron–positron plasmas, which excites surrounding nebulae observed in X-ray. Here we show with a three-dimensional particle simulation for the global axisymmetric magnetosphere that a steady outflow of electron–positron pairs is formed with associated pair sources, which are the gamma-ray emitting regions within the light cylinder. The magnetic field is assumed to be a dipole, and to be consistent, the pair creation rate is taken to be small, so that the model might be applicable to old pulsars such as Geminga. The pair sources are charge-deficient regions around the null surface, and we identify them as the outer gap. The wind mechanism is the electromagnetic induction which brings about fast azimuthal motion and eventually trans-field drift by radiation drag in the close vicinity of the light cylinder and beyond. The wind causes loss of particles from the system. This maintains charge deficiency in the outer gap and pair creation. The model is thus in a steady state, balancing loss and supply of particles. Our simulation implies how the wind coexists with the gamma-ray emitting regions in the pulsar magnetosphere.     

The example of effective plasma acceleration in a magnetosphere

The problem of effective transform of Poynting flux energy into the kinetic energy of relativistic plasma outflow in a magnetosphere is considered. In this article we present an example of such acceleration. In order to perform it, we use the approach of ideal axisymmetric magnetohydrodynamics (MHD). For highly magnetized plasma outflow we show that a linear growth of Lorentz factor with a cylindrical distance from the rotational axis is a general result for any field configuration in the sub-magnetosonic flow. In the far region the full magnetohydrodynamics problem for one-dimensional flow is considered. It turns out that the effective plasma outflow acceleration is possible in the paraboloidal magnetic field. It is shown that such an acceleration is due to the drift of charged particles in the crossed electric and magnetic field. The clear explanation of the absence of acceleration in the monopole magnetic field if given.      

Relativistic anisotropic pair plasmas

The properties of waves able to propagate in a relativistic pair plasma are at the basis of the interpretation of several astrophysical observations. For instance, they are invoked in relation to radio emission processes in pulsar magnetospheres and to radiation mechanisms for relativistic radio jets. In such physical environments, pair plasma particles probably have relativistic, or even ultrarelativistic, temperatures. Besides, the presence of an extremely strong magnetic field in the emission region constrains the particles to one-dimensional motion: all the charged particles strictly move along magnetic field lines.
We take anisotropic effects and relativistic effects into account by choosing one-dimensional relativistic Jűttner–Synge distribution functions to characterize the distribution of electrons and/or positrons in a relativistic, anisotropic pair plasma. The dielectric tensor, from which the dispersion relation associated with plane wave perturbations of such a pair plasma is derived, involves specific coefficients that depend on the distribution function of particles. A precise determination of these coefficients, using the relativistic one-dimensional Jűttner–Synge distribution function, allows us to obtain the appropriate dispersion relation. The properties of waves able to propagate in anisotropic relativistic pair plasmas are deduced from this dispersion relation. The conditions in which a beam and a plasma, both ultrarelativistic, may interact and trigger off a two-stream instability are obtained from this same dispersion relation. Two astrophysical applications are discussed.     

The Environment of YSO Jets

It is commonly accepted that stars form in molecular clouds by the gravitational collapse of dense gas. However, it is precisely not the infalling but the outflowing material that is primarily observed. Outflow motions prevail around both low and high mass young stellar objects. We present here results from a family of self-similar models that could possibly help to understand this paradox. The models take into account the heating of the central protostar for the deflection and acceleration of the gas. The models make room for all the ingredients observed around the central objects, i.e. molecular outflows, fast jets, accretion disks and infalling envelopes. We suggest that radiative heating and magnetic field may ultimately be the main energy sources driving outflows for both low and high mass stars. The models show that the ambient medium surrounding the jet is unhomogeneous in density, velocity, magnetic field. Consequently, we suggest that jets and outflows have a prehistory that is inprinted in their environment, and that this should have direct consequences on the setting of jet numerical simulations.     

Influence of relativistic effects and vacuum polarization on the transfer of gyroresonance radiation and the stability of the atmospheres of compact stars

We consider the transfer of radiation and calculate the force of its pressure in the electron gyroresonance line in the atmospheres of magnetic degenerate stars. We specify the atmospheric parameters for which an outflow of plasma is possible under radiation pressure in the cyclotron line. We show that the permittivity tensor of a mildly relativistic plasma in a strong magnetic field found by applying relativistic corrections to the cyclotron resonance condition and by taking into account the vacuum polarization and recoil effects during photon scattering should be used to obtain proper results. We have determined the real and imaginary parts of the refractive indices and the polarization coefficients for normal electromagnetic waves when scattering dominates over absorption. Relativistic effects, which change greatly the dispersion and resonant absorption of waves propagating almost perpendicular to the magnetic field, and vacuum polarization have been found to change qualitatively the gyroresonance radiation spectrum and pressure for a wide range of parameters of stellar magnetospheres.     

Studying millisecond pulsars in X-rays

Millisecond pulsars represent an evolutionarily distinct group among rotation-powered pulsars. Outside the radio band, the soft X-ray range (～0.1–10 keV) is most suitable for studying radiative mechanisms operating in these fascinating objects. X-ray observations revealed diverse properties of emission from millisecond pulsars. For the most of them, the bulk of radiation is of a thermal origin, emitted from small spots (polar caps) on the neutron star surface heated by relativistic particles produced in pulsar acceleration zones. On the other hand, a few other very fast rotating pulsars exhibit almost pure nonthermal emission generated, most probably, in pulsar magnetospheres. There are also examples of nonthermal emission detected from X-ray nebulae powered by millisecond pulsars, as well as from pulsar winds shocked in binary systems with millisecond pulsars as companions. These and other most important results obtained from X-ray observations of millisecond pulsars are reviewed in this paper, as well as results from the search for millisecond pulsations in X-ray flux of the radio-quite neutron star RX J1856.5-3754.     

Formation of Hydrogen Emission Lines in the Magnetospheres of Young Stars

Astronomy Letters - The formation of hydrogen emission lines in the magnetospheres of young stars is considered. The magnetosphere is assumed to be formed by a dipolar magnetic field whose axis is...     

Plasma magnetosphere of rotating magnetized neutron star in the braneworld

Plasma magnetosphere surrounding rotating magnetized neutron star in the braneworld has been studied. For the simplicity of calculations Goldreich-Julian charge density is analyzed for the aligned neutron star with zero inclination between magnetic field and rotation axis. From the system of Maxwell equations in spacetime of slowly rotating star in braneworld, second-order differential equation for electrostatic potential is derived. Analytical solution of this equation indicates the general relativistic modification of an accelerating electric field and charge density along the open field lines by brane tension. The implication of this effect to the magnetospheric energy loss problem is underlined. It was found that for initially zero potential and field on the surface of a neutron star, the amplitude of the plasma mode created by Goldreich-Julian charge density will increase in the presence of the negative brane charge. Finally we derive the equations of motion of test particles in magnetosphere of slowly rotating star in the braneworld. Then we analyze particle motion in the polar cap and show that brane tension can significantly change conditions for particle acceleration in the polar cap region of the neutron star.     

Outline of a theory of solar wind interaction with the magnetosphere

Some new ideas on the interaction of the solar wind with the magnetosphere are brought forward. The mechanism of reflection of charged particles at the magnetopause is examined. It is shown that in general the reflection is not specular but that a component of momentum of the particle parallel to the magnetopause changes. A critical angle is derived such that particles whose trajectories make an angle less than it with the magnetopause enter the magnetosphere freely, so transferring their forward momentum to it. Spatially or temporally non-uniform entry of charged particles into the magnetosphere causes electric fields parallel to the magnetopause which either allow the free passage of solar wind across it or vacuum reconnection to the interplanetary magnetic field depending on the direction of the latter. These electric fields can be discharged in the ionosphere and so account qualitatively for the dayside agitation of the geomagnetic field observed on the polar caps. The solar wind wind plasma which enters the magnetosphere creates (1) a dawn-dusk electric field across the tail (2) enough force to account for the geomagnetic tail and (3) enough current during disturbed times to account for the auroral electrojets. The entry of solar wind plasma across the magnetosphere and connection of the geomagnetic to interplanetary field can be assisted by wind generated electric field in the ionosphere transferred by the good conductivity along the geomagnetic field to the magnetopause. This may account for some of the observed correlations between phenomena in the lower atmosphere and a component of magnetic disturbance.     

3-D Rpic Simulations of Relativistic Jets: Particle Acceleration, Magnetic Field Generation, and Emission

We have applied numerical simulations and modeling to the particle acceleration, magnetic field generation, and emission from relativistic shocks. We investigate the nonlinear stage of theWeibel instability and compare our simulations with the observed gamma-ray burst emission. In collisionless shocks, plasma waves and their associated instabilities (e.g., the Weibel, Buneman and other two-stream instabilities) are responsible for particle (electron, positron, and ion) acceleration and magnetic field generation. 3-D relativistic electromagnetic particle (REMP) simulations with three different electron-positron jet velocity distributions and also with an electron-ion plasma have been performed and show shock processes including spatial and temporal evolution of shocks in unmagnetized ambient plasmas. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-positron shocks are larger for smaller jet Lorentz factor. This comes from the fact that the growth time of the Weibel instability is proportional to the square of the jet Lorentz factor. We have performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which is assumed to be created by photon annihilation. Simulation results with this broad distribution show that the Weibel instability is excited continuously by the wide-range of jet Lorentz factor from lower to higher values. In all simulations the Weibel instability is responsible for generating and amplifying magnetic fields perpendicular to the jet propagation direction, and contributes to the electron’s (positron’s) transverse deflection behind the jet head. This small scale magnetic field structure contributes to the generation of “jitter” radiation from deflected electrons (positrons), which is different from synchrotron radiation in uniform magnetic fields. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks. The detailed studies of shock microscopic process evolution may provide some insights into early and later GRB afterglows.     

Magnetorotational and Tayler Instabilities in the Pulsar Magnetosphere

The magnetospheres around neutron stars should be very particular because of their strong magnetic field and rapid rotation. A study of the pulsar magnetospheres is of crucial importance since it is the key issue to understand how energy outflow to the exterior is produced. In this paper, we discuss magnetohydrodynamic processes in the pulsar magnetosphere. We consider in detail the properties of magnetohydrodynamic waves that can exist in the magnetosphere and their instabilities. These instabilities lead to formation of magnetic structures and can be responsible for short-term variability of the pulsar emission.     

Radiation pressure supported stars in Einstein gravity: eternally collapsing objects

Even when we consider Newtonian stars, that is, stars with surface gravitational redshift   z ≪ 1  , it is well known that, theoretically, it is possible to have stars supported against self-gravity almost entirely by radiation pressure. However, such Newtonian stars must necessarily be supermassive. We point out that this requirement for excessively large M in the Newtonian case is a consequence of the occurrence of low   z ≪ 1  . However, if we remove such restrictions, and allow for the possible occurrence of a highly general relativistic regime,   z ≫ 1  , we show that it is possible to have radiation pressure supported stars (RPSSs) at an arbitrary value of M . Since RPSSs necessarily radiate at the Eddington limit, in Einstein gravity, they are never in strict hydrodynamical equilibrium. Further, it is believed that sufficiently massive or dense objects undergo continued gravitational collapse to the black hole (BH) stage characterized by   z =∞  . Thus, late stages of BH formation, by definition, would have   z ≫ 1  , and hence would be examples of quasi-stable general relativistic RPSSs. It is shown that the observed duration of such Eddington limited radiation pressure dominated states is   t ≈ 5 × 10⁸ (1 + z ) yr  . Thus,   t →∞  as BH formation  ( z →∞)  takes place. Consequently, such radiation pressure dominated extreme general relativistic stars become eternally collapsing objects (ECOs) and the BH state is preceded by such an ECO phase. This result is also supported by our previous finding that trapped surfaces are not formed in gravitational collapse and the value of the integration constant in the vacuum Schwarzschild solution is zero. Hence the supposed observed BHs are actually ECOs.