A three-ring circuit model of the magnetosphere

We have modeled the magnetosphere by superimposing a dipole field, a uniform field and a perturbation field due to a simple current system. This current system consists of a ring current in the neutral line of the dipole plus uniform fields, together with vertical currents representing field-aligned currents to the neutral line. The current circuit is closed by two additional ring currents above and below the equatorial plane representing distributed adiabatic perpendicular currents. This system produces many magnetospheric features including a magnetopause, bending of magnetic field lines in the anti-solar direction, a magnetotail, and cusps on the day-side of the Earth. Our aim is to demonstrate that it is not necessary to think of the magnetic field topology as being caused by the flowing plasma carrying field lines. The fundamental physical problem is to derive the current system from the self-consistent interaction of the solar-wind and magnetospheric plasmas and fields.     

Inductive electric field of a time-dependent ring current

The inductive electric field generated by a time-dependent symmetric ring current has been investigated. The symmetric ring current was modelled by a population of protons drifting in a magnetic dipole field. The interaction of these protons with ion-cyclotron waves was assumed to be the dominant energy loss process for the ring current protons, at least under certain conditions. The calculation, with spectral densities for the ion-cyclotron waves that are based on experimental data, showed that an azimuthal inductive electric field of as much as 0.25 mV/m can be produced by this mechanism. Furthermore there is evidence that if the spectral density of the waves is substantially larger than the one adopted here, the electric field might increase to the order of 1.0 mV/m or more.     

The effect of an electric field induced by a time-dependent ring current on the particle drift motion

The effect of an electric field induced by a rapidly decaying ring current on the motion of charged particles in the magnetosphere has been investigated using Euler potentials. For a model consisting of the earth dipole and the symmetric ring current, the electric field satisfies the condition E . B = 0.

Under this circumstance, the E × B drift of the particle can be identified as the motion of the magnetic field lines and vice versa. The time dependent electric field induced can be evaluated in a Spherical polar coordinate system by the formula

where and β are Euler potentials.

A model calculation on the particle drift velocity v_D = E × B/B² shows that the radial component of the drift velocity is in good agreement with those deduced from whistler duct studies.     

A mathematical magnetospheric field model with independent physical parameters

A quantitative magnetospheric magnetic field model has been calculated in three dimensions. The model is based on an analytical solution of the Chapman-Ferraro problem. For this solution, the magnetopause was assumed to be an infinitesimally thin discontinuity with given geometry. The shape of the dayside magnetopause is in agreement with measurements derived from spacecraft boundary crossings.The magnetic field of the magnetopause currents can be derived from scalar potentials. The scalar potentials result from solutions of Laplace's equation with Neumann's boundary conditions. The boundary values and the magnetic flux through the magnetopause are determined by all magnetic sources which are located inside and outside the magnetospheric cavity. They include the Earth's dipole field, the fields of the equatorial ring current and tail current systems, and the homogeneous interplanetary magnetic field. In addition, the flux through the magnetopause depends on two constants of interconnection which provide the possibility of calculating static interconnection between magnetospheric and interplanetary field lines. Realistic numerical values for both constants have been derived empirically from observed displacements of the polar cusps which are due to changes in the orientation of the interplanetary field. The transition from a closed to an open magnetosphere and vice versa can be computed in terms of a change of the magnetic boundary conditions on the magnetopause. The magnetic field configuration of the closed magnetosphere is independent of the amount and orientation of the interplanetary field. In contrast, the configuration of the open magnetosphere confirms the observational finding that field line interconnection occurs primarily in the polar cusp and high latitude tail regions.The tail current system reflects explicitly the effect of dayside magnetospheric compression which is caused by the solar wind. In addition, the position of the plasma sheet relative to the ecliptic plane depends explicitly on the tilt angle of the Earth's dipole. Near the tail axis, the tail field is approximately in a self-consistent equilibrium with the tail currents and the isotropic thermal plasma.The models for the equatorial ring current depend on the D_st-parameter. They are self-consistent with respect to measured energy distributions of ring current protons and the axially symmetric part of the magnetospheric field.     

Formation of whistler ducts

A theory of whistler duct formation is presented. By means of order of magnitude calculations it is shown that, when the ring current overlaps the outer plasmasphere, irregularities will cause field-aligned currents to flow, which are below the threshold sensitivity of satellite-borne magnetometers. These currents must be continuous with horizontal ionospheric currents, which produce horizontal electric fields. These fields map up to the equatorial plane and are large enough to produce flux tube interchange and hence the formation of whistler ducts in the outer plasmasphere.     

Nonthermal electrons in the thick-target reverse-current model for hard X-ray bremsstrahlung

The behaviour of the accelerated electrons escaping from a high-temperature source of primary energy in a solar flare is investigated. The direct current of fast electrons is supposed to be balanced by the reverse current of thermal electrons in the ambient colder plasma inside flare loops. The self-consistent kinetic problem is formulated; and the reverse-current electric field and the fast electron distribution function are found from its solution. The X-ray bremsstrahlung polarization is then calculated from the distribution function. The difference of results from those in the case of thermal runaway electrons (Diakonov and Somov, 1988) is discussed. The solutions with and without account of the affect of a reverse-current electric field are also compared.     

On current continuity at the Harang discontinuity

Although the Harang discontinuity has so far been identified in terms of various phenomena (such as ground magnetic fields, ionospheric currents, auroral features, and electric fields), the loci defined by those different phenomena do not always coincide. It is suggested that the Harang discontinuity may not be a line boundary across which the electric field changes its direction simply from poleward to equatorward, but that the field gradually rotates counterclockwise in a narrow region; thus the westward electric field dominates there. In such a case, no field-aligned current is necessarily required to flow from or into the discontinuity region. This view may be contrasted with the conventional view that an intense upward field-aligned current should flow from the Harang discontinuity. A model is presented in which the poleward ionospheric current (the Hall current resulting from the westward electric field) in the Harang discontinuity region connects the eastward electrojet and the westward electrojet.     

Thermal electrons runaway from a hot plasma during a flare in the reverse-current model and their X-ray bremsstrahlung

The behaviour of the thermal electrons escaping from a hot plasma to a cold one during a solar flare is investigated. We suppose that the direct current of fast electrons is compensated by the reverse current of the thermal electrons in ambient plasma. It is shown that the direct current strength is determined only by the regular energy losses due to Coulomb collisions. The reverse-current electric field and the distribution function of fast electrons are found in the form of an approximate analytical solution to the self-consistent kinetic problem of the dynamics of a beam of escaping thermal electrons and its associated reverse current.The reverse-current electric field in solar flares leads to a significant reduction of the convective heat flux carried by fast electrons escaping from the high-temperature plasma to the cold one. The spectrum and polarization of hard X-ray bremsstrahlung, and its spatial distribution along flare loops are calculated and can be used for diagnostics of flare plasmas and escaping electrons.Send offprint requests to B. V. Somov.     

Model for the onset of a magnetospheric substorm

In view of observations which show that a substorm often begins in a small local time sector, a model is assumed in which the neutral sheet current is diverted around a small region we call a bubble. The simplest assumption is that of a linear variation of current with distance from the centre of the bubble in the x-direction in a SM coordinate system, with the diverted current being channelled within narrow paths of width δ_y on the dawn and dusk sides of the bubble. This assumption leads to vector potential integrals that can be evaluated analytically. The addition of this current loop into the magnetotail results in a magnetic field structure where new neutral lines of X- and 0-type can be observed; these are connected to each other as a continuous neutral ring in the xy equatorial plane. The magnetic and electric field components around the neutral regions are calculated, and the time dependent evolution of the neutral ring is studied. Comparison with some published satellite observations shows good agreement. Taking typical values for the various quantities on the basis of actual observations within the magnetotail, we show that the induced electric field is at least comparable to the average cross-tail electrostatic field, and it may well be one or two orders of magnitude greater. The response of the plasma to the induction field is discussed qualitatively. It is concluded that field aligned currents may be produced due to inertial forces of the expanding disturbance. Interpretation of the ground based precipitation patterns of energized particles during auroral breakup is given.     

Electrostatic Potentials in Supernova Remnant Shocks

Recent advances in the understanding of the properties of supernova remnant shocks have been precipitated by theChandra and XMM X-ray Observatories, and the HESS Atmospheric Čerenkov Telescope in the TeV band. A critical problem for this field is the understanding of the relative degree of dissipative heating/energization of electrons and ions in the shock layer. This impacts the interpretation of X-ray observations, and moreover influences the efficiency of injection into the acceleration process, which in turn feeds back into the thermal shock layer energetics and dynamics. This paper outlines the first stages of our exploration of the role of charge separation potentials in non-relativistic electron-ion shocks where the inertial gyro-scales are widely disparate, using results from a Monte Carlo simulation. Charge density spatial profiles were obtained in the linear regime, sampling the inertial scales for both ions and electrons, for different magnetic field obliquities. These were readily integrated to acquire electric field profiles in the absence of self-consistent, spatial readjustments between the electrons and the ions. It was found that while diffusion plays little role in modulating the linear field structure in highly oblique and perpendicular shocks, in quasi-parallel shocks, where charge separations induced by gyrations are small, and shock-layer electric fields are predominantly generated on diffusive scales.     

A model for the formation of spokes in Saturn's ring

We suggest that spokes consist of charged micron-sized dust particles elevated from the rings by radially moving dense plasma columns created by meteor impacts on the ring. Dense plasma causes electrostatic wall-sheaths at the ring and charging of the ring with electric fields strong enough to overcome the gravitational force on small dust particles. Under “ordinary” conditions only very few dust particles will be elevated as the probability of a dust particle having at least one excess electronic charge is very low. Dense plasma raises this probability significantly. The radial motion of the plasma column is due to an azimuthal polarization electric field built up by the relative motion between the corotating plasma and the negatively charged dust particles which move with a Keplerian speed.     

Self-consistent electromagnetic field in a nearly co-rotating magnetosphere

An example of the self-consistent solution which belongs to the non-trivial solution, obtained in a previous paper (Kaburaki, 1985), is found in a nearly co-rotating inner magnetosphere. Though the stellar wind is neglected there compared with the co-rotatinal velocity, drift motion around the magnetic axis, which is a manifestation of inertial effects, is determined self-consistently with the electromagnetic field. In this process, explicit expressions for the energy integral in the rotating frame and for the density distribution are also obtained. These expressions contain a fundamental length, which is to be evaluated according to physical conditions of a magnetosphere and determines the asymptotic-kinetic energy of a plasma particle at infinity. The electric current associated with the drift motion is too small to alter the original magnetic field, but the electric field is modified by the inertial effects even in the inner magnetosphere. The integrated Ohm's law is used to describe a force balance in the rotating frame, in the limits of weak and strong magnetic field.     

A numerical study of polar ionospheric currents

Numerical calculations for the electric current in the polar ionosphere have been made by assuming some realistic distributions of the electric field and conductivity. Two dynamo actions are taken into account; one of which is induced by ionospheric winds and the other by the solar wind. For the solar wind dynamo action, it is found that the secondary polarization field caused by non-uniform distribution of ionospheric conductivity is much larger than the primary field induced by the solar wind, suggesting its important effect on charged particles in the magnetosphere, and that the irrotational current having a source and sink is of the same order of magnitude as the solenoidal current closing its circuit in the ionosphere. It is also found that the solar wind is, in general, more effective than the ionospheric winds in producing polar current systems such as DP 1 and 2, but in some cases the ionospheric winds have a significant effect on the current distribution.     

Collisionless Magnetic Reconnection in the Presence of Initial Guide Field

By using the method of 2-dimensional, 3-component full particle simulation, collisionless magnetic reconnection in the presence of various initial guide fields and the Harris current sheet with 1-dimensional initial state are studied. The results show that strong guide fields with B_z0 > 0.5B₀ can evidently alter not only the trajectory of the particles, but also the structure of the electric and velocity fields in the vicinity of the reconnection region, thereby affecting the rate of reconnection and the acceleration of electrons. The generalized Ohm's law is employed to interpret the structural characteristics of the electric fields with various guide fields. Also, via the tracing of the electron beam near he diffusion region, it is revealed that in the 2-D model, for both strong and weak guide fields, the induced electric field perpendicular to the simulation plane at the center of the diffusion region plays the major role in the acceleration of electrons. The contribution of the planar electric field outside the diffusion region is very small.     

Potential and inductive electric fields in the magnetosphere during auroras

During quiescent auroras the large-scale electric field is essentially irrotational. The volume formed by the plasma sheet and its extension into the auroral oval is connected to an external source by electric currents, which enter and leave the volume at different electric potentials and which supply sufficient energy to support the auroral activity. The location of the actual acceleration of particles depends on the internal distribution of electric fields and currents. One important feature is the energization of the carriers of the cross-tail current and another is the acceleration of electrons precipitated through relatively low-altitude magnetic-field-aligned potential drops.Substorm auroras depend on rapid and (especially initially) localized release of energy that can only be supplied by tapping stored magnetic energy. The energy is transmitted to the charged particle via electric inductive fields.The primary electric field due to changing electric currents is redistributed in a complicated way—but never extinguished—by polarization of charges. As a consequence, any tendency of the plasma to suppress magnetic-field-aligned components of the electric fields leads to a corresponding enhancement of the transverse component.     

Pulsar magnetosphere

Pulsars are presently believed to be rotating neutron stars with large frozen-in magnetic fields normally assumed to be dipole fields. It has been shown that such a star must possess a magnetosphere if it rotates sufficiently rapidly. By assuming that the magnetic field is dipolar, and unaffected by the trapped particles in the magnetosphere, and that the field dipole axis is parallel to the rotation axis, Goldreich and Julian determined many of the properties of the magnetosphere. In this paper is given a self-consistent model of the closed field lines of a pulsar magnetosphere. Using this model, it is shown that, close to the star, the above assumptions of Goldreich and Julian are justified. Their results are extended to the oblique rotator as well as to stars with magnetic multipoles of arbitrary order and arbitrary orientation.Supported in part by the U.S. Atomic Energy Commission under Grant 2171T.     

Saturn's sweeper moons predicted

The recent observation of the absorption of radiation belts in the vicinity of Saturn's bright rings and historical observations of the ring system make the following related results apparent:

- The gaps in the rings are caused by the presence of at least 6 small, extremely dense and probably electrically charged ‘sweeper’ moons which effectively sweep the ring matter clean from the gaps. This is known due to the fading of the inner ring edges whereas the outer edges are well defined. Their orbital periods will differ from the expected Keplerian periods if the moons and Saturn do possess electric fields.

- Absorption of radiation belts near the rings (of Jupiter also) implies that the ring particles themselves are not absorbing the radiation but the small moons are. This is consistent with the observed radiation belt absorption near the outer Saturnian moons.

- If electric fields of the sweeper moons cause the ring edge fading as observed (and not simply gravitational), then Saturn itself must maintain an electric field in its vicinity by way of a sizeable proton wind to affect the uneven ring edge fading and will be surrounded by an H⁺ cloud at least to approximately the A-ring. this is consistent with the detection of an H⁺ cloud surrounding Saturn (Weiseret al., 1977, p. 755). The other possibility is that these moons are extremely dense and have very large internal magnetic fields.

- Because of their location, these moons must be captured and if very dense as believed, may be core remnants of a nova.

     

激光驱动亥姆霍兹电容线圈靶磁重联实验的数值模拟

激光驱动亥姆霍兹电容线圈靶的磁重联实验已经提出并进行了多年.当实验中的金属板被强激光照射时产生自由电子,这些自由电子的运动在连接两金属板的两个平行线圈中产生电流,由两个平行线圈内部电流产生的磁场之间随即发生重联.该实验不同于其他直接由Biermann电池效应所产生高β(等离子体热压与磁压的比值)环境下的磁重联实验.对该类实验进行了3维磁流体动力学数值模拟,首次展示了亥姆霍兹电容器线圈靶如何驱动磁重联的过程.数值模拟结果清楚地表明,磁重联的出流等离子体在线圈周围发生与实验结果相一致的堆积现象.线圈电流产生的磁场可高达100 T,使得磁重联区域周围的等离子体β值达到10^-2.与实验室结果进行比较,数值模拟重复了实验展示的大多数特征,可有助于深入认识和理解实验结果背后的物理学原理.