The role of charge-exchange cross-section for pickup protons and neutrals in the inner heliosheath

The process of deceleration of the solar wind downstream of the termination shock is studied on the basis of a one-dimensional multi-component model. It is assumed that the solar wind consists of thermal protons, electrons and interstellar pickup protons. The protons interact with interstellar hydrogen atoms by charge-exchange. Two cases are considered. In the first one, the charge-exchange cross-section for thermal protons and hydrogen atoms is the same as for pickup protons and atoms. Under this condition, there is a strong dependence of the solar wind velocity on the downstream temperature of pickup protons. When the proton temperature is close to 10 keV, the change in the velocity with the distance from the termination shock is similar to that measured on the Voyager 1 spacecraft: linear velocity decrease is accompanied by an extended transition region with near-zero velocity. However, with a more careful approach to the choice of the charge-exchange cross-section, the situation changes dramatically. The strong dependence of the solar wind speed on the pickup proton temperature disappears and the transition region in the heliosheath disappears as well, at least at reasonable distances from the TS.     

Motion of solar cosmic rays in the coronal magnetic field

Trajectories of solar cosmic rays have been calculated in a static ninth-order coronal magnetic field. It is found that as a result of field curvature and gradients, protons drift across the field lines at a rate of up to 200 ² deg hr^–1. These drift rates are of the same order as, but somewhat smaller than, empirically derived rates. Localized enhancements of magnetic field have been inserted into the ninth-order field in order to model (in a highly idealized manner) the effects of the small-scale magnetic features which give rise to X-ray bright points. The motions of the particles in the presence of these scattering centers can be parameterized approximately by a cross-field diffusion coefficient. Our estimates of this coefficient, although crude, overlap with empirical values which have been deduced over a wide range of energies.We propose that coronal propagation of solar cosmic rays has two components. One is independent of particle velocity, and is associated with dynamic field phenomena (such as an expanding magnetic bottle): this is the only component which is important in flares which occur close to the foot-point of the Sun-Earth field line. The second component is velocity dependent, but is independent of mass, and is associated with scattering off (relatively static) magnetic inhomogeneities with scale sizes of at least 500 km: the second component contributes to coronal propagation if the flare occurs more than about 50–60 deg away from the Sun-Earth field line.     

Energetic protons from the solar flare of March 24, 1966

Ten to 100 meV protons from the solar flare of March 24, 1966 were observed on the University of California scintillation counter on OGO-I. The short rise and decay times observed in the count rates of the 32 channels of pulse-height analysis show that scattering of the protons by the interplanetary field was much less important in this event than in previously observed proton flares. A diffusion theory in which D = M _r is found to be inadequate to account for the time behavior of the count rates of this event. Small fluctuations of the otherwise smooth decay phase may be due to flare protons reflected from the back of a shock front, which passed the earth on March 23.     

Photoelectrons and solar wind/lunar limb interaction

It is suggested that boundary conditions for solar wind/lunar limb interactions are active. The whole-Moon limb does not evoke a shock cone because warm (13 eV/electron) solar wind electrons are replaced by cool (2 eV/electron) photoelectrons that are ejected from the generally smooth areas of the lunar terminator illuminated at glazing angles by the Sun. A localized volume of low thermal pressure is created in the solar wind by these cool photoelectrons. The solar wind expands into this turbulence-suppressive volume without shock production. Conversely, directly illuminated highland areas exchange hot photoelectrons (> 20 eV/electron) for warm solar wind electrons. The hot electrons generate a localized pressure increase (p) in the adjacent solar wind flow which evokes a shock streamer in the solar wind. Shock streamers are identifiable by a coincident increase in the magnitude (B p) of the solar wind magnetic field immediately external to the lunar wake. Shock occurrence is controlled by lunar topography, solar activity in the hard ultraviolet (> 20 eV), solar wind electron density and thermal velocity, and the intensity of the solar wind magnetic field.Paper dedicated to Professor Harold C. Urey on the occasion of his 80th birthday on 29 April 1973.The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR 09-051-001 with the National Aeronautics and Space Administration.     

A model of the magnetized solar wind

A steady state, inviscid, single fluid model of the solar win d in the equatorial plane is developed using magneto-hydrodynamics and including the heat equation wit h thermal conduction but no non-thermal heating (i.e. a conduction model). The effects of solar rotation and magnetic field are included enabling both radial and azimuthal components of the velocity and magnetic fields to be found in a conduction model for the first time.The magnetic field cuts off the thermal conduction far from the sun and leads to an increased temperature at 1 AU and relatively small changes to the radial velocity and density. Models have been found which fit the experimental electron densities in 2 R < r < 16 R . These models predict at 1 AU a radial velocity of 300–380 km·sec^-1 and a density of 8 protons·cm^-3. The latter velocity corresponds to a density profile obtained by Blackwell and Petford (1966) during the last sunspot minimum, and is about 100 km·sec^-1 above that found in previous conduction models which fit the coronal electron densities. The radial velocities are now consistent with the mean quiet solar wind, as are the densities when the experimental values are averaged over a magnetic sector. However, the azimuthal velocity at 1 AU is only 1–2 km·sec^-1 which is low compared to the experimental values, as found by previous authors.     

Two-dimensional guiding-center model of the solar wind-moon interaction

This paper presents a continued study of the two-dimensional guiding-center model of the solar wind interaction with the Moon. The characteristics theory and the computational method are discussed. The magnetic permeability of plasma is (1 + /2)^–1 in the solar wind flow upstream of the Moon, and it changes to 1 in the void region of the lunar wake. The gradual change of the magnetic permeability in the penumbral region from the interplanetary condition to the void condition is explained as the source of field perturbations in the lunar wake. Perturbations of the magnetic field propagate as magnetoacoustic waves in a frame of reference moving with the plasma flow. Computer solutions were obtained to show that (i) the two principal perturbations of the magnetic field in the lunar wake (the umbral increase and the penumbral decrease) are confined to a region bounded by a Mach cone tangent to the lunar body, and (ii) the penumbral increases occur outside the lunar Mach cone. Computer solutions are also used to identify the source of field perturbations and to simulate the solar wind-moon interaction under varying interplanetary conditions.     

Collisional heating of interplanetary gas: Fokker-Planck treatment

Interstellar gas streaming through the solar system undergoes both elastic collisions with solar wind ions and destructive, ionizing processes. The Boltzmann equation is set up, with linear Fokker-Planck terms describing the glancing elastic collisions. Solutions combining the dynamical effects of the central force field and the diffusion in velocity space are derived, appropriate to cool gas.For the He component of the streaming gas, if initially at 100 K, the collisional heating dominates inside 2 a.u. upstream and 5 a.u. downstream. A modified formula is given for the density in the downstream wake, as enhanced by gravitational focussing. Calculations of the helium resonant radiation backscatter require substantial modification.     

Anisotropy characteristics of low energy cosmic ray population of solar origin

Comparison of Explorer 34 observations on solar protons in the energy range 0.7–55 MeV with similar observations from other spacecrafts show that the large field aligned anisotropies which are observed during the rise time of a flare event change to an equilibrium anisotropy coming radially from the sunward direction due to the convective removal of the solar particles. At very late times during the decay (T 4 days) the anisotropy is observed to be from a direction 45° E of the satellite-Sun line which is interpreted as indicative of positive density gradient of solar cosmic ray population. The dependence of both types of equilibrium anisotropies on the energy and the velocity of the particles and on plasma velocity are shown to be in agreement with the theoretical predictions. The amplitude of the large field aligned anisotropies observed earlier in the event is found to be independent of the rise time of the event and to vary as (Vt)^–1.Interplanetary magnetic sector crossings during a flare event, cause abrupt changes in both the amplitude and phase of the non-equilibrium anisotropy whereas they have no significant effect on equilibrium anisotropy. The effect of azimuthal density gradients on the decay time constants of flare enhancements are also examined in an attempt to understand the complicated structures, often observed, in the time intensity profiles at low energies.Part of this work was done while the author was at the University of Texas as Dallas, U.S.A.Now at the National Academy of Sciences, Washington, D.C.     

A note on the instantaneous rotational velocity of Mercury

The fluctuation in the angular velocity of the present rotation of Mercury is investigated. The instantaneous rotational rate in terms of orbital mean motion at different positions along Mercury's orbit is given. At aphelion the rotational velocity decreases substantially because the solar gravitational torque on the two thermal bulges on Mercury's surface tends to retard the rotation of Mercury before aphelion passage. It is found that the difference between the rotational periods derived from the motions at perihelion and aphelion is 4.68 min and that the maximum rate of rotation occurs atf=/4 andf=7/4, wheref is the true anomaly.     

Long-Term Radio Scintillation Variations in the Circumsolar Plasma

Long-term scintillation measurements of the solar wind formation zone at solar elongations ranging from 1°–8° (Sun impact parameters: 4–30 R ) were recorded using the water maser source IRC-20431 at the wavelength =1.35 cm during its annual solar occultations in December 1981–1998. Dramatic changes in the spatial dependence of the scintillation index were recorded over the course of the 11-year solar cycle. Markedly diminished scattering, attributed to a pronounced heliolatitude effect, was observed at the closest solar approach distances in the years around solar activity minimum. From parallel investigations of the solar magnetic field structure it was determined that the field strength at the source of the solar wind streamlines is the governing factor for the solar wind acceleration process. Particularly apparent in the scintillation data during solar activity minimum is the increasing role of the polar coronal holes with their associated open magnetic field structure. The dependence of the solar scattering intensity on heliolatitude fades in the years of high solar activity as the level of scintillations increases at polar latitudes.     

A dominant role for protons at the onset of solar flares

The energetics of the onset of the impulsive phase of solar flares are examined on the premise that a single acceleration mechanism is operating in the corona. From considerations of recent observations of plasma turbulence and upflows, and nuclear gamma-rays it is concluded that a model where the bulk of the energy resides in a non-thermal electron beam with a low energy cut-off at 20–25 keV is incompatible with many of the observations. Conversely, a model where the bulk of the energy resides in non-thermal protons is consistent with the majority, if not all, of the observations. It is suggested that the bulk of the energy in the impulsive phase is initially transferred to 10²–10³ keV protons. Acceleration by a series of small shocks is an energy transfer mechanism which gives particles increments in velocity rather than energy and would naturally favour protons over electrons. An important consequence of this result is that the hard X-ray burst must be thermal. At this time the precise mechanism for thermal X-ray production is unclear; however recent theoretical plasma physics results have indicated promising avenues of research in this context.     

Recurrency and the origin of the vertical component of the interplanetary magnetic field

The autocorrelation functions of the solar wind velocity and of the IMF components as well as of the geomagnetic activity indices are studied for two periods: August–December, 1965 and January–May, 1974. The vertical component of the IMF is shown to exhibit a rather definite recurrency relatively independent of the recurrency of the solar wind velocity.The daily mean values of the Z-component of the IMF are shown to correlate ( = -0.5) with the intensity of the meridional component of the large scale solar magnetic field with time delay of about 5 days with respect to the latter. This result is interpreted as an evidence for the Z- component of the IMF to be carried away by the solar wind from the Sun.     

A new outlook on the `Differential Theory' of the solar Quadrupole Moment and Oblateness

The modeling of the quadrupole moment J ₂ and of the oblateness , two key solar parameters, derives from the development in successive spherical harmonics of the gravitational potential. These harmonics are representative of the shape of the Sun, by studying the local distortion of the internal layers, under their distribution of mass and velocity. The first aim of this paper is to study, over the radius r and the colatitude , the structure of the internal layers of the Sun through a geometrical approach, considering J ₂ and under a differential form. The second aim is to determine their theoretical values, after integration over r and , taking the best available models of density and rotation into account constrained by helioseismic data. The novelty of our approach lies in the analysis of the profiles of the two above-mentioned solar parameters, under differential form, from the core to the surface. This analysis allows us to comply with the physical processes located in the transition regions, namely the tachocline and maybe a new subsurface layer which could be called the leptocline. The profiles of tildeJ ₂ show two sharp decreases, which are directly connected to the shear layers located at 0.7 R and beneath the surface. The profiles of tilde show five changes of curvature, which seem to be connected to solar processes, such as the matter circulation flows, seismic events or the storage of the magnetic field, phenomena taking place in the transition regions. These sets of profiles allow us to propose a configuration scenario composed of a double layer. Moreover, as a result of the theoretical determination of tildeJ ₂ and tilde, the values at the surface of the quadrupole moment and of the oblateness can be deduced, which are 1.60×10^–7 and 8.77×10^–6, respectively. As a result of an analysis of available data, we may admit J ₂=(2.0±0.4)×10^–7. The theoretical computations of J ₂ and at the surface will be confronted in the near future with the values measured in space by means of the Picard microsatellite. This mission should permit one to measure at the same time both the total solar irradiance and the latitudinal diameters in any position angle (after removing the passing spots or faculae at the limb).     

Dynamics of a cloud of fast electrons travelling through the plasma

A simulation of normal type III radio bursts has been made in a whole frequency range of about 200 MHz to 30 kHz by the usage of the semi-analytical method as developed in previous papers for the plasma waves excited by a cloud of fast electrons. Three-dimensional plasma waves are computed, though the velocities of fast electrons are assumed to be one-dimensional. Many basic problems about type III radio bursts and associated solar electrons have been solved showing the following striking or unexpected results.Induced scattering of plasma waves, by thermal ions, into the plasma waves with opposite wave vectors is efficient even for a solar electron cloud of rather low number density. Therefore, the second harmonic radio emission as attributed to the coalescence of two plasma waves predominates in a whole range from meter waves to km waves. Fundamental radio emission as ascribed to the scattering of plasma waves by thermal ions is negligibly small almost in the whole range. On the other hand, third harmonic radio emission can be strong enough to be observed in a limited frequency range.If, however, the time integral of electron flux is, for example, 2 × 10¹³ cm^–2 (>5 keV) or more at the height of 4.3 × 10¹⁰ cm ( _p = 40 MHz) above the photosphere, the fundamental may be comparable with or greater than the second harmonic, but an effective area of cross-section of the electron beam is required to be very small, 10¹⁷ cm² or less, and hence much larger sizes of the observed radio sources must be attributed to the scattering alone of radio waves.The radio flux density expected at the Earth for the second harmonic can increase with decreasing frequencies giving high flux densities at low frequencies as observed, if x-dependence of the cross-sectional area of the electron beam is x ^1.5 or less instead of x ², at least at x 2 × 10¹² cm.The second harmonic radio waves are emitted predominantly into forward direction at first, but the direction of emission may reverse a few times in a course of a single burst showing a greater backward emission at the low frequencies.In a standard low frequency model, a total number of solar electrons above 18 keV arriving at the Earth orbit reduces to 12% of the initial value due mainly to the collisional decay of plasma waves before the waves are reabsorbed by the beam electrons arriving later. However, no deceleration of the apparent velocity of exciter appears. A change in the apparent velocity, if any, results from a change in growth rate of the plasma waves instead of the deceleration of individual electrons.Near the Earth, the peak of second harmonic radio flux as emitted from the local plasma appears well after the passage of a whole solar electron cloud through this layer. This is ascribed to the secondary and the third plasma waves as caused in non-resonant regions by the induced scattering of primary plasma waves in a resonant region.     

Energy Spectra of low-Flux Protons in the Inner Heliosphere under Quiet Solar Conditions

The energy spectra of protons at energies in the range of about 1–100 MeV are investigated during time periods of low solar activity using data sets from near Earth spacecraft. These populations pose a tough experimental and theoretical problem that remains unsolved up to now. We attempt to provide a consistent definition of low-flux quiet-time periods relevant to low solar activity as well as quasi-stationary periods useful at higher levels of solar activity. Using statistical methods, the possible instrumental contribution to the lowest observed proton fluxes for various detectors is estimated. We suggest and prove that there exists a low-flux population of charged particles in the energy range of about 1–10 MeV, which is present in the inner heliosphere even during the quietest conditions at lowest solar activity. The dynamics of the variations of proton spectra over the solar cycle is investigated. A series of low-flux periods is examined in detail and energy spectra of protons are approximated in the form of J(E)=AE ^–+CE. By determining the best fitting parameters to the energy spectra correlations are made among them as well as with monthly sunspot numbers characterizing solar activity. It has been demonstrated that the value of the energy minimum of proton spectrum E _min that `divides' the two populations – `solar/heliospheric' and `galactic' – is shifted towards higher values with increasing solar activity. Protons have been argued to be predominantly of solar origin up to several MeV near the solar cycle minimum and up to 20–30 MeV at maximum. The slope of the lower spectrum branch (parameter ) slightly decreases with increasing solar activity. The minimum fluxes observed during the last 3 minima of solar activity are compared; the lowest fluxes were those during the 1985–1987 period.     

The acceleration of minor ion species in the solar wind

This paper provides a comprehensive analysis of the dynamics of the flow of minor ion species in the solar wind under the combined influences of gravity, Coulomb friction (with protons), rotational forces (arising from the Sun's rotation and the interplanetary spiral magnetic field) and wave forces (induced in the minor ion flow by Alfvén waves propagating in the solar wind). It is assumed that the solar wind can be considered as a proton-electron plasma which is, to a first approximation, unaffected by the presence of minor ions. In the dense hot region near the Sun Coulomb friction accelerates minor ions outwards against the gravitational force, part of which is cancelled by the charge-separation electric field. Once the initial acceleration has been achieved, wave and rotational forces assist Coulomb friction in further increasing the minor ion speed so that it becomes comparable with, or perhaps even exceeds, the solar wind speed. A characteristic feature of the non-resonant wave force is that it tends to bring the minor ion flow into an equilibrium where the radial speed matches the Alfvén speed relative to the solar wind speed, whereas Coulomb friction and rotational forces tend to bring the flow into an equilibrium where the radial speed of the minor ions equals the solar wind speed. Therefore, provided that there is sufficient wave energy and Coulomb friction is weak, the minor ion speed can be trapped between these two speeds. This inteststing result is in qualitative agreement with observational findings to the effect that the differential flow speed between helium ions and protons is controlled by the ratio of the solar wind expansion time to the ion-proton collision time. If the thermal speeds of the protons and minor ions are small compared to the Alfvén speed, two stable equilibrium speeds can exist because the rapid decrease in the Coulomb cross-section with increasing differential flow speed allows the non-resonant wave force to balance Coulomb friction at more than one ion speed. However, it must be emphasized that resonant wave acceleration and/or strong ion partial pressure gradients are required to achieve radial speeds of minor ions in excess of the proton speed, since, as is shown in Section 4, the non-resonant wave acceleration on protons and minor ions are identical when their radial speeds are the same, with the result that, in the solar wind, non-resonant wave acceleration tends (asymptotically) to equalize minor ion and proton speeds.     

Thermal and tidal effect on the rotation of Mercury

It is shown that the influences of the thermal and tidal effects on Mercury's libration are in equilibrium with the periods of rotation and revolution of Mercury locked in the 32 resonant state. The suggestion by Liu that the solar gravitational couple on the thermal bulges accelerates Mercury's rotation is investigated and the production of mechanical energy to balance the dissipation of the bodily tides is discussed. It is possible for Mercury to rotate with two bulges as a solar thermal engine; the tidal effect causes this engine to function and its maximum power is close to 10¹⁶ ergs per sec.     

Mass motions in a flare spray

High velocity H ejections were observed in association with an important solar flare on March 12, 1969, and simultaneously with Type II followed by Type IV radio emission detectable to 3–4 solar radii (R ) from the center of the Sun. From a sequence of H coronagraph photographs, trajectories and velocity determinations were made for fragments of the flare spray which was visible to a distance of 2 R . The temporal and spatial relationship between the optical and radio events is discussed. The mass motions appear to be controlled by the gravitational field while the fragments move in the direction of the open magnetic field lines.     

Analysis of the H ring velocity distribution function near the Martian and Venusian bow shocks by an analytical model

We have studied the H⁺ velocity distribution function at Mars and Venus near the bow shock both in the solar wind and in the magnetosheath by a simple analytical one-dimensional model. We found that over half of the ions in the ring velocity distribution which moved towards the magnetosheath were scattered back into the bow shock. The original ring distribution is destroyed in less than an ion gyro period. Ions contained in the magnetosphere which hit the bow shock bounce back into the solar wind with a maximum energy exceeding twice the energy of solar wind protons. The ions finite gyroradius causes an asymmetric flow in the magnetosheath with respect to the direction of the convective electric field, which can be observed already a few ion gyroradius downstream of the bow shock.     

Distribution of Planetesimals around a Protoplanet in the Nebula Gas

We have developed a semi-analytic method of calculating the changes in heliocentric Keplerian orbital elements due to gravitational scattering by a protoplanet as a three-body problem. In encounters with high incident velocities, either the gravity of the protoplanet or the solar gravity can be regarded as perturbation force. In close encounters, by taking into account the solar gravity as a perturbation, we modified the two-body gravitational scattering. On the other hand, in slightly distant encounters, we apply the perturbing force of the protoplanet to the heliocentric Keplerian orbit of planetesimals. As a result, as for high-velocity encounters, the three-body problem is semi-analytically solvable. Our semi-analytic method can reproduce the numerical result of the orbital changes of individual planetesimals for the broad range of high-energy encounters with surprising high accuracy. We found that our method is valid under the conditions (i)b₀? 2 and (ii) (e²₀+i²₀− b²₀)^1/2? 4, wheree₀andi₀are eccentricity and inclination of relative motion normalized by the reduced Hill radius andb₀is the difference between semimajor axes normalized by the Hill radius. Though our method needs some numerical procedure, its cpu time is negligibly short compared with that of the direct orbital integration. In simulation of orbital evolution of planetesimals around a protoplanet in the gas, which we will perform in the subsequent paper, most encounters can be calculated by the semi-analytic method. This makes it possible to perform the long term (∼10⁵years) orbital calculation of ∼10^3–4planetesimals.