Computer simulation of the particle acceleration in pulsar magnetospheres

In the spherically-symmetric case, a computer simulation of the electron acceleration inside the outflow channel of the pulsar magnetosphere is produced. The stationary motion of electrons is shown to be unstable in the case of > _c, where is a parameter describing inhomogeneity of the background charge, and _c is its critical value. The arising non-stationary motion of electrons leads to a formation of electron bunches, which move chaotically. The mean electron energy appears to be much greater at the non-stationary motion, than at the stationary one. The time-averaged parameters of the non-stationary electron flow and their dependence upon have been investigated. Distributions of the mean values of parameters (charge density, electron velocity, electric field energy density, pressure, and internal energy of the gas composed of the electron bunches) over the magnetosphere altitude have been investigated. The mean spectra of the charge density have been obtained. The results of numerical investigation of the spherically-symmetric model are used for estimation of the electron energy and of the electron flux in the case of the more realistic model. The radioemission loss is estimated, and is shown to be large enough for explaining the radiopulsar phenomenon as a thermal radioemission of the pulsar magnetosphere. In particular, such common properties of the pulsar radioemission as the high bright temperature, the sharp radioemission directivity, and the characteristic turn-over of the radioemission spectrum at the frequency of the order 10⁸ Hz are found a natural explanation in frames of this model.     

Hemispherical asymmetries in sunspot areas and auroral frequencies

One thousand and fifty-two aurorae boreales and 554 aurorae australes recorded during the nineteenth century at medium latitudes 55° N or 55° S are compared statistically with the known hemispherical asymmetry of the sums of the areas of sunspots. According to the present study, the solar hemispherical asymmetry may be accompanied by an analogous pattern of the hemispherical frequency of auroral days. For the number of auroral days in each hemisphere beyond the two auroral ovals, a remarkable degree of phase equality with the sunspot areas during the second Gleissberg cycle can be demonstrated.     

Diffuse extragalactic gamma-radiation above ‘black-body cutoff’

The possible contribution of gamma-rays predicted within the universal cosmic-ray (CR) hypothesis to the energy range of CR spectrum above black-body cutoff is calculated. These gamma-rays arise from the relativistic electromagnetic cascade generated in the field of microwave background radiation (MBR). The ultra-high energy photons and electrons that initiate the cascade are produced at the decay of -mesons created in interactions of photons with the MBR. Simple analytic expression for cascade gamma-ray spectrum is obtained from the solution of kinetic equations for electrons and photons as well as for protons propagating in the MBR field. It is shown that at certain values of magnetic field and radio-wave density in the intergalactic space te flux of cascade gamma-rays may at least partly mask the black-body cutoff in the CR spectrum.Deceased, August 13, 1989.     

Containment of an adiabatic plasma on magnetic lines of force by a self-generated electrostatic field

One component of a three-fluid adiabatic plasma is under certain conditions contained in a restricted region of space by a large-scale electrostatic field generated within the plasma. The containment is discussed here for plasma consisting of ions and two populations of electrons characterized by different pitch angle distribution functions.The bouncing motion of electrons along open fieldlines between a magnetic mirror and an electrostatic mirror produces a velocity distribution function similar to that generated by bouncing particles on closed fieldlines.     

The global structure of the pulsar magnetospheres

The model of the pulsar magnetosphere filled with massless charged particles (rest massm=0) is considered. The gas of charged massless particles can be found in two different phases: (1) dynamical phase (DP), when the particles move with nonvanishing energy along some base line, determined by the electromagnetic field only, (2) statical phase (SP), when the particles have vanishing energy =0. The pulsar magnetosphere occurs to be divided into regions of different types: (a) the accelerating regions (DP-regions) containing only DP, (b) the capture regions, containing only SP, (3) leaky capture regions, where DP moves through SP. The leaky capture regions are the active regions, which are responsible for the pulsar radioemission. In the case of oblique magnetic moment the equation for the capture surface has been obtained. The capture region is formed around this surface. Expressions for jumps of the electromagnetic field, the current density and the charge density on the capture region boundary have been obtained. The problem of the pulsar magnetosphere is stated mathematically in the case of oblique magnetic moment and ejection of only electrons.     

An examination of directional discontinuities and magnetic polarity changes around interplanetary sector boundaries usingE > 2 keV electrons

Past studies of interplanetary magnetic sector boundaries have been based on the assumption that one can determine the field polarities by comparing the field directions with those of the nominal Parker spiral angles. Previous investigators have found evidence for decreases of B, the magnitude of the magnetic fieldB, and increases of , the angle betweenB and the ecliptic plane, at sector boundaries. Others have argued that the characteristic thickness of sector boundaries exceeds that of tangential discontinuities, making sector boundaries a separate class of structures.We use a simple technique for inferring the polarities of interplanetary magnetic fields based on the assumption thatE > 2 keV electrons are always flowing along the magnetic field away from the Sun. Electron data from the UC Berkeley experiment on the ISEE-3 spacecraft are used to examine periods around several apparent sector boundaries in 1978 and 1979. We compare properties of (a) boundaries with field polarity changes and (b) large-angle ( > 60°) directional discontinuities with no field polarity changes. We find no significant differences between the sector boundaries and the directional discontinuities in terms of associated decreases in B or of values of . These results suggest no significant difference between sector boundaries and directional discontinuities other than the change in field polarities. Within limited statistics we find that about half the polarity changes would not have been identified using a requirement that > 90° and that half of the > 120° discontinuities would have been misidentified as polarity changes.     

The spectra of extended radio sources with a nonuniform magnetic field and hard injection of particles

On the basis of an analytical solution of the diffusion-type kinetic equation for electrons, electron distributions and radiation spectra have been found which result from a hard injection of particles in sources of the core halo type, characterized by spatially nonuniform magnetic fields and diffusion parameters. Such radio sources are shown to possess nonlinear radiation spectra containing universal (=0.5) and diffusion-controlled power-law sections shaped by synchrotron losses, spatial diffusion and radiation conditions of the electrons. The diffusion-controlled sections can be described by spectral indices 0.5<1, if the magnetic field decreases towards the source edge, and by <0.5 where the magnetic field increases.     

Complete determination of the magnetic field vector and of the electron density in 14 prominences from linear polarizaton measurements in the HeI D3 and Hα lines

The present paper is devoted to the interpretation of linear polarization data obtained in 14 quiescent prominences with the Pic-du-Midi coronagraph-polarimeter by J. L. Leroy, in the two lines Hei D₃ andH quasi-simultaneously. The linear polarization of the lines is due to scattering of the anisotropic photospheric radiation, modified by the Hanle effect due to the local magnetic field. The interpretation of the polarization data in the two lines is able to provide the 3 components of the magnetic field vector, and one extra parameter, namely the electron density, because the linear polarization of H is also sensitive to the depolarizing effect of collisions with the electrons and protons of the medium. Moreover, by using two lines with different optical thicknesses, namely Hei D₃, which is optically thin, and H, which is optically thick ( = 1), it is possible to solve the fundamental ambiguity, each line providing two field vector solutions that are symmetrical in direction with respect to the line of sight in the case of the optically thin line, and which have a different symmetry in the case of the optically thick line.It is then possible to determine without ambiguity the polarity of the prominence magnetic field with respect to that of the photospheric field: 12 prominences are found to be Inverse polarity prominences, whereas 2 prominences are found to be Normal polarity prominences. It must be noticed that in 12 of the 14 cases, the line-of-sight component of the magnetic field vector has a Normal polarity (to the extent that the notion of polarity of a vector component is meaningful; no polarity can be derived in the 2 remaining cases); this may explain the controversy between the results obtained with methods based on the Hanle effect with results obtained through the Zeeman effect. A dip of the magnetic field lines across the prominence has been assumed, to which the optically thick H line is sensitive, and the optically thin Hei D₃ line is insensitive.For the Inverse prominences, the average field strength is 7.5±1.2 G, the average angle,, between the field vector and the prominence long axis is 36° ± 15°, the average angle, , between the outgoing field lines and the solar surface at the prominence boundary is 29° ± 20°, and the average electron density is 2.1 × 10¹⁰ ± 0.7 × 10¹⁰ cm^–3. For the Normal prominences, the average field strength is 13.2±2.0 G, the average angle,, between the field vector and the prominence long axis is 53° ± 15°, the average angle, , between the outgoing field lines and the solar surface at the prominence boundary is 0° ± 20° (horizontal field), and the average electron density is 8.7 × 10⁹ ± 3.0 × 10⁹ cm^–3.     

PROBLEMS AND PROGRESS IN COMPUTING THREE-DIMENSIONAL CORONAL ACTIVE REGION MAGNETIC FIELDS FROM BOUNDARY DATA

An overview of the whole process of reconstructing the coronal magnetic field from boundary data measured at the photosphere is presented. We discuss the errors and uncertainties in the data and in the data reduction process. The problems include noise in the magnetograph measurements, uncertainties in the interpretation of polarization signals, the 180° ambiguity in the transverse field, and the fact that the photosphere is not force-free. Methods for computing the three-dimensional structure of coronal active region magnetic fields, under the force-free assumption, from these boundary data, are then discussed. The methods fall into three classes: the extrapolation technique, which seeks to integrate upwards from the photosphere using only local values at the boundary; the current-field iteration technique, which propagates currents measured at the boundary along field lines, then iteratively recomputes the magnetic field due to this current distribution; and the evolutionary technique, which simulates the evolution of the coronal field, under quasi-physical resistive magnetohydrodynamic equations, as currents injected at the boundary are driven towards the observed values. The extrapolation method is mathematically ill-posed, and must be heavily smoothed to avoid exponential divergence. It may be useful for tracing low-lying field lines, but appears incapable of reconstructing the magnetic field higher in the corona. The original formulation of the current-field iteration method had problems achieving convergence, but a recent reformulation appears promising. Evolutionary methods have been applied to several real datasets, with apparent success.     

High-resolution photography of the solar chromosphere

A method is described for reconstructing the true geometry of a solar loop observed on the disk which takes account of tilt in its own plane. Reconstructions of three H loops yield small tilt angles (14°) and provide further evidence that H loops show a close correspondence to the field lines of a magnetic dipole. The method offers new opportunities for exploring the physics of individual solar loops.     

Basic concepts in the electrodynamics of the auroral ionosphere

Equations governing the large-scale electrodynamic processes in the auroral ionosphere are systematically discussed and the limits to drawing conclusions from incomplete sets of equations are evaluated. The vectors of electic current density,j, and electric field,E, are expressed as explicit functions of the densities, pressures and velocities of the constituents of the ionosphere.The equation div (·E)=0 is an identity satisfied by any solution of the full set of equations governing the problem and cannot be treated as a differential equation forE in which the components of the conductivity tensor are given parameters. The concept of the height-integrated conductivities and the conclusions based on it are inconsistent with the equations of momentum balance for the ionospheric constituents.The global structure of the auroral ionosphere is determined by the state of equilibrium between the pressure gradients, the inertial forces and thej×B-force associated with the auronal electrojets flowing along the auroral oval. The time-averaged, global, electric field is directed across the auroral oval. Its value is substantially affected by the motions of neutral particles. The velocity vector of the neutrals has a substantial component directed across the oval.     

Nonpotential waves in the auroral electrojet

An ionospheric plasma instability in the auroral electrojet region believed to be responsible for electromagnetic waves with angular frequency smaller thanv _i is discussed. The threshold current velocityU _t and oscillation growth rate are found for the realistic physical situation when a nonuniform plasma density and an ionization function are considered. It is found that the gradient of plasma density n may be neglected in the expressions forU _t and and that the spectral density of fluctuations in the density of the plasma agrees well with the experimental data.     

Pulsar radio emission

Radio emission by pulsars is calculated from first principles. In an almost current-free magnetosphere, the two charged components (of the unsteadily escaping pair plasma) have different (and varying) bulk Lorentz factors. Curvature radiation emitted by the more energetic component is thus locally coherent, (so-called antenna mechanism). Strong enough seed signals cause the relativistically streaming charges to enhance their radiation, via an induced drift that can largely exceed the curvature drift. This amplification mechanism is similar to - but different from - that of a maser; we call it a MAIDER. Maximal amplification occurs at an (emission) altitude where the two components have sufficiently separated in energy though not yet separated too strongly in space.     

Simulation for electron acceleration by DC electric field in the presence of ion sound waves and associated hard X-ray emission

Time-dependent Fokker-Planck equation was numerically solved to demonstrate the dynamics of electrons in a uniform coronal loop with an applied axial DC electric field in the presence of ion-sound waves. This electric field is attributed to an anomalous resistivity due to the ion-sound turbulence caused by an initially given critical current density.The electron momentum distribution becomes a steady state in the whole turbulent region in a short time for which some electrons can be accelerated to the maximum electric potential K _c. The steady energy distribution of electrons flowing out the end of the turbulent region has a very hard power-law-like spectrum with an index of about 0.75. The associated hard X-rays from a thick target also show a hard spectrum with a photon spectral index of 1.3. In order for to be much greater as observed in impulsive X-ray bursts, it is required that the source is a sum of many elementary loops with a power-law-like distribution in K _c with an index = – + 2.5.     

Magnetic field and electric current structure in the chromosphere

Three dimensional vector magnetic field structure throughout the chromosphere above an active region is deduced by combining high resolution H filtergrams with a simultaneous digital magnetogram. An analog model of the field is made with 400 metal wires representing fieldlines which are assumed to outline the H structure. The height extent of the field is determined from vertical field gradient observations around sunspots, from observed fibril heights and from an assumption that the sources of the field should be largely local. After digitization the magnetic field H matrix is retrieved. Electric current densities j are computed from j=curl H. The currents (typically 10 mA m^–2) flow in patterns not similar to observed features and not parallel to magnetic fields. Lorentz forces are computed from {ie0323-01}. The force structures correspond to observed solar features and a series of observed dynamics may be expected: downward motion in bipolar areas in lower chromosphere, an outflow of the outer chromosphere into the corona with radially outward flow above bipolar plage regions (where coronal streamers are observed) and motions of arch filament systems. Observed current structure and magnitude agree well with previous vector magnetograph observations but disagree with theoretical current-free or force-free concepts. A dynamic chromosphere with electromagnetic forces in action is thus inferred from observations.     

Magnetic reconnection in high-temperature plasma of solar flares

Quasi-steady high-temperature current sheets are an energy source during the main or hot phase of solar flares. Such sheets are shown to be stabilized with respect to the tearing instability by a small transverse component of magnetic field existing in the sheets.     

Signatures of solar wind interaction with the nightside ionosphere of Venus

Plasma and field relationships observed across the nightside of Venus evidence a chaotic variety of interactions between the ionosphere and the combined effect of the solar wind and interplanetary magnetic field draped about the planet. Close examination of these data reveal within the chaos a number of repeatable signatures key to understanding fundamental field-plasma interactions. Observed from the Pioneer Venus Orbiter, (PVO), nightside conditions range from extensive, full-up ionospheres with little evidence of dynamic or energetic perturbations, to an almost full depletion, sometimes described as disappearing ionospheres. Between these extremes, the ionospheric structure is often irregular, sometimes exhibiting well-defined density troughs, at other times complex intervals of either abundant or minimal plasma concentration. Consistently, large B-fields (typically exceeding 5–10 nanoteslas) coincide with plasma decreases, whereas stable, abundant plasma distributions are associated with very low-level field. We examine hundreds of nightside orbits, identifying close correlations between regions of elevated magnetic fields featuring polarity reversals, and (a) exclusive low-frequency or distinctive broadband noise, or both, in the electric field data, (b) turbulent, superthermal behavior of the the ions and electrons. We review extensive studies of nightside fields to show that the correlations observed are consistent with theoretical arguments that the presence of strong magnetic fields within normal ionospheric heights indicates the intrusion of magnetosheath fields and plasma within such regions. We find abundant evidence that the ionosphere is frequently disrupted by such events, exhibiting a chaotic, auroral-like complexity appearing over a wide range of altitude and local time. We show that field-plasma disturbances, widely suggested to be similar to conditions in the Earth's auroral regions, are tightly linked to the electric field noise otherwise attributed to lightning. Owing to the coincidence inherent in this relationship, we suggest that natural, predictable plasma instabilities associated with the plasma gradients and current sheets evident within these events produce the E-field noise. The data relationships argue for a more detailed investigation of solar wind induced E-field noise mechanisms as the appropriate scientific procedure for invoking sources for the noise previously attributed to lightning. Consistent with these views, we note that independent analyses have offered alternative explanations of the noise as arising from ionospheric disturbances, that repeated searches for optical evidence of lightning have found no such evidence, and that no accepted theoretical work has yet surfaced to support the inference of lightning at Venus.     

A cometary ionosphere model for Io

The ionosphere of Jupiter's satellite Io, discovered by the Pioneer 10 radio-occultation experiment, cannot easily be understood in terms of a model of a gravitationally bound, Earth-like ionosphere. Io's gravitational field is so weak that a gravitationally bound ionosphere would probably be blown away by the ram force of the Jovian magnetospheric wind — i.e., the plasma corotating in the Jovian magnetosphere. We propose here a model in which the material for Io's atmosphere and ionosphere is drawn from the ionosphere of Jupiter through a Birkeland current system that is driven by the potential induced across Io by the Jovian corotation electric field. We argue that the ionization near Io is caused by a comet-like interaction between the corotating plasma and Io's atmosphere. The initial interaction employs the critical velocity phenomenon proposed many years ago by Alfvén. Further ionization is produced by the impact of Jovian trapped energetic electrons, and the ionization thus created is swept out ahead of Io in its orbit. Thus, we suggest that what has been reported as a day-night ionospheric asymmetry is in fact an upstream-downstream asymmetry caused by the Jovian magnetospheric wind.Paper dedicated to Professor Hannes Alfvén on the occasion of his 70th birthday, 30th May, 1978.     

Circulating subbeam systemsand the physics of pulsar emission

The purpose of this paper is to suggest how detailed single-pulse observations of slow radio pulsars may be utilized to construct an empirical model for their emission. It links the observational synthesis developed in a series of papers by Rankin in the 1980s and 90s to the more recent empirical feedback model of Wright (2003a) by regarding the entire pulsar magnetosphere as a non-steady, non-linear interactive system with a natural built-in delay. It is argued that the enhanced role of the outer gap in such a system indicates an evolutionary link to younger pulsars, in which this region is thought to be highly active, and that pulsar magnetospheres should no longer be seen as being driven by events on the neutron stars polar cap, but as having more in common with planetary magnetospheres and auroral phenomena.Received: 8 May 2003, Published online: 14 November 2003 Correspondence to: Joanna M. Rankin. On leave from: Physics Department, University of Vermont, Burlington, VT 05405, USA     

Chaos and secondary resonances in the mimas–tethys system

We have investigated the role of the 200yr period discovered by Vienne and Duriez (1992) on the tidal evolution of the Mimas–Tethys system through the 2:4 ii present resonance. Three terms are found to generate this period. We present a perturbedpendulum model in which these terms bring about a perturbation to the ideal ii resonance pendulum, which is in a direct ratio to the eccentricity e of Tethys. Although e is now very small, it is shown that this quantity could have been much greater in the past. We also show, thanks to this model, that these terms may have brought about a stochastic layer of noticeable width at the time of capture in the ii resonance, with the consequence that the possible values of the inclination i of Mimas before capture range from 0.4° to 0.6° (these uncertainties arise from the present uncertainties on e). The role of each one of the three terms is examined in the appearance of chaos. A capture into the 1/1 secondary resonance (between the libration period of the primary ii resonance and the period of about 200yr) is found possible. It means that the system could have experienced several captures in the primary resonance, instead of a single one, and that i could have been, with this assumption, much lower than 0.4°. A probability of capture into this secondary resonance as a function of the eccentricity of Tethys on encounter is derived, using Malhotra's method (Malhotra, 1990). Allan's values of i = 0.42° and e 0 (Allan, 1969) are therefore called into question, and taking e 0 is shown to be absolutely necessary if we want to understand the phenomena at work in the Mimas–Tethys system.