Numerical studies of the transport of solar protons in interplanetary space

Numerical solutions of the Fokker-Planck equation governing the transport of solar protons are obtained using the Crank-Nicholson technique with the diffusion coefficient represented by K_r=K₀r^b where r is radial distance from the Sun and b can take on positive or negative values. As b ranges from +1 to ?3, the time to the observation of peak flux decreases by a factor of 5 for 1 MeV protons when $\text{V}\text{K}{}_0\text{=}\text{3 AU}{}^{\mathrm{b?1}}$ where V is the solar wind speed. The time to peak flux is found to be very insensitive to assumptions concerning the solar and outer scattering boundary conditions and the presence of exponential time decay in the flux does not depend on the existence of an outer boundary. At $\text{V}\text{K}{}_0\text{? 15 AU}{}^{\mathrm{b?1}}$, 1 MeV particles come from the Sun by an almost entirely convective process and suffer large adiabatic deceleration at b?0 but for b=+1, large Fermi acceleration is possible at all reasonable $\text{V}\text{K}{}_0$ values. Implications of this result for the calculation and measurement of particle diffusion coefficients is discussed. At $\text{b?0}$, the pure diffusion approximation to transport overestimates by a factor 2 or more the time to peak flux but as b becomes more negative, the additional effects of convection and energy loss become less important.     

Electron runaway in turbulent astrophysical plasmas

An astrophysical electron acceleration process is described which involves turbulent plasma effects: the acceleration mechanism will operate in ‘collision free’ magnetoactive astrophysical plasmas when ion-acoustic turbulence is generated by an electric field which acts parallel to the ambient magnetic lines of force. The role of ‘anomalous’ (ion-sound) resistivity is crucial in maintaining the parallel electric field. It is shown that, in spite of the turbulence, a small fraction of the electron population can accelerate freely, i.e. runaway, in the high parallel electric potential. The number density n^(B) of the runaway electron component is of order $\text{n}{}^{\mathrm{(B)}}\text{?}\text{n}\text{2}\text{(}\text{c}{}_{\mathrm{s}}\text{U}{}_{\mathrm{\parallel }}{}^{\mathrm{?}}\text{)}{}^2$, where n = background electron number density, c_s = ion-sound speed and U_∥^? = relative drift velocity between the electron and ion populations. The runaway mechanism and the number density n^(B) do not depend critically on the details of the non-linear saturation of the ion-sound instability.     

Total current to cylindrical collectors in collisionless plasma flow

A theory is presented for charged-particle collection by a cylindrical conducting object, such as a spacecraft or an electrostatic probe, which is moving transversely through a collisionless plasma, such as those in the upper atmosphere and space. The calculation is approximate, using symmetric potential profiles which are exact for the infinite-cylinder stationary case. Theoretical current predictions are presented for ratios of collector potential to electron thermal energy eφ_c/kT_e from 0 to ?25, for ion-to-electron temperature ratios T_i/T_c = 1 and 0.5, ratio of collector radius to electron Debye length r_c/λ_D from 0 to 100, and ratio of flow speed to ion thermal speed $\text{S}{}_{\mathrm{i}}\text{= U/(2kT}{}_{\mathrm{i}}\text{/m}{}_{\mathrm{i}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}$ from 0 to 10. Comparisons with existing exact calculations by other authors show that none of these fulfil all of the requirements for nontrivial comparison. Appropriate parameter ranges for future exact calculations are thereby suggested. These are as follows: (a) r_c/λ_D should be large enough that the collector not be in or near orbit-limited conditions; (b) the ratio $\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,\; i}}\text{¦}$ of ion directed energy to potential energy change in the sheath, should be close to unity or if

$\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,i}}\text{¦? 1,}\text{then}\text{S}{}_{\mathrm{i}}\text{? 1}$

.     

On the stability of the inner planets,satellites and close binaries

I consider the range of Hill stability in the restricted circular problem of three bodies when the larger one of the two principal bodies has a finite oblateness. I show that the range r satisfies the equation

$\text{r =}\text{1? μ}\text{C}{}_{\mathrm{cr}}\text{? 3μ + μr}{}^2\text{? v (1? μ)r}{}^{\mathrm{?3}}\text{±(2 + 3v)}\text{(1 ? μ}\text{1 +}\text{3v}\text{r}{}^2\text{r}\text{,}$

where μ is the mass parameter and v is an oblateness parameter. This result is applied to the solar system, the Earth-Moon system and binary star systems. It is then shown that, all the inner planets of the solar system, the great majority of asteroids and some short-period comets are Hill stable, that direct artificial satellites of the Earth are more stable than retrograde ones, and that contact binaries possess cores between which no mass exchange takes place.     

The locating of hydromagnetic whistler sources and determination of their generating proton spectra

Results are given of the calculations of the group delay time propagating τ(ω, φ₀) of hydromagnetic whistlers, using outer ionospheric models closely resembling actual conditions. The τ(ω, φ₀) dependencies were compared with the experimental data of τ_exp(ω, φ₀) obtained from sonagrams. The sonagrams were recorded in the frequency range $\text{? ? (0.5?2.5) Hz}$ at observation points located at geomagnetic latitudes φ₀ = (53?66)° and in the vicinity of the geomagnetic poles. This investigation has led us to new and important conclusions.The wave packets (W.P.) forming hydromagnetic whistlers (H.W.) are mainly generated in the plasma regions at L = 3.5?4.0. This is not consistent with ideas already expressed in the literature that their generation region is $\text{L ? 3?10}$. The overwhelming majority of the τ_exp values differ considerably from the times at which wave packets would, in theory, propagate along the magnetic field lines corresponding to those of the geomagnetic latitudes φ₀ of the observation points. The second important fact is that the W.P. frequency ω is less than $\text{Ω}{}_{\mathrm{H}}$ everywhere along its propagation trajectory, including the apogee of the magnetic force line ($\text{Ω}{}_{\mathrm{H}}$ is the proton gyrofrequency). Proton flux spectra $\text{E ? (30?120)}$ keV, responsible for H.W. generation, were determined. Comparison of the Explorer-45 and OGO-3 measurements published in the literature, with our data, showed that the proton flux density energy responsible for the H.W. excitation $\text{N}{}_{\mathrm{p}}\text{(}\text{MV}{}_6{}^2\text{2}\text{) ? (5 × 10}{}^{\mathrm{?3}}\text{?10}{}^{\mathrm{?1}}\text{)}\text{H}{}_{\mathrm{a}}{}^2\text{8π}$ where H_a is the magnetic field force in the generation region of these W.P. The electron concentration is $\text{N}{}_{\mathrm{a}}\text{? (10}{}^2\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The values given in the literature are $\text{N}{}_{\mathrm{a}}\text{? (10}{}^{\mathrm{?10}}\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The e data considered also leads to the conclusion that the generating mechanism of the W.P. studied probably always co-exists with the mechanism of their amplification.     

A note on the asymptotic density-potential relation of a spiral galaxy with a finite thickness

It is shown that the asymptotic σ₁(r) and ψ₁(r) relations can be derived very simply by using the method of double series expansion, where σ₁, ψ₁(r,0) and ψ₁ are the surface density perturbation, the gravitational potential perturbation at the symmetric plane Z=0 and the average potential perturbation respectively. The results are accurate to the order of both m²(kr)^?2 and k〈∣z∣〉, where m is the number of spiral arms, k is the radial wave number, r is the distance from the centre of the galaxy, and 〈∣z∣〉 is the average vertical distance of a star from the Symmetrie plane Z=0. Such an accuracy is usually sufficient for the discussion of spiral modes in a spiral galaxy of small but finite disk thickness. It is pointed out that ψ₁(r,0)～(σ₁(r) relation can be expressed in a unified form for different vertical density profiles if 〈∣z∣〉 is adopted as the thickness scale, and that ψ₁(r,0)～(σ₁(r) can be expressed in a unified form for different vertical density profiles if 〈∣z?z∣〉 the average vertical separation between two stars, is adopted as the thickness scale. Only the value of the ratio $\text{〈|z?z′|〉z}\text{〈|z|〉}$ is a functional of the vertical density profile. However, for the usual physically meaningful profiles, these values are very close to each other: It is $\text{2}$ for the Gaussian profile, $\text{1}\text{Ln}\text{2}\text{= 1.443}$ for the $\text{rmsech}{}^2\text{(}\text{z}\text{z}{}_1\text{(r)}\text{)}$ profile, and 1.5 for the $\text{exp}\text{[}\text{?|z|}\text{z}{}_1\text{(r)}\text{]}$ profile.     

Love numbers of the giant planets

We calculate the Love numbers k_n for n = 2 to 10, and determine the “gravitational noise” from tides. The new values k₂ for Jupiter, Saturn, and Uranus yield new estimates for the planetary dissipation functions: $\text{Q}{}_{\mathrm{<mtext>J</mtext>}}\text{? 2.5 × 10}{}^4$, $\text{Q}{}_{\mathrm{<mtext>S</mtext>}}\text{? 1.4 × 10}{}^4\text{, Q}{}_{\mathrm{<mtext>U</mtext>}}\text{? 5 × 10}{}^3$.     

New high-speed photometry of EH Lib and a binary interpretation of its period change

Six times of maxima of the ultrashort-period cepheid variable EH Librae were measured in 1980 May to June and in 1981 January, with a three-channel photocounting high-speed photoelectric photometer. These, together with all the photoelectric times of maxima over the past 30 years, are used to re-examine the nature of the change of the period. We found that we can fix the times of maxima by the following formula

$\text{T}{}_{\max }\text{= T}{}_0\text{+P}{}_0\text{E+}\text{1}\text{2}\text{βE}{}^2\text{+A}\text{sin}\text{2π}\text{EP}{}_0\text{E}{}_0$

where T₀ = HJD 2433438.6088 and P₀ = 0.0884132445 d represent the initial maximum epoch and the pulsation period, β = ?2.8 × 10^?8/yr; A = 0.0015 d, P₀ = 6251 d = 17.1 yr are the semi-amplitude and the period of the sine curve, and E is the number of periods elapsed since T₀, and (E₀ = 70700).If we interpret this 17.1 year periodicity as a modulation of the phase of maximum by binary motion, then the semi-amplitude of the orbital radial velocity variation is K = 2πasini/E₀ = 0.45 km/s and the mass function is

$\text{f(m)=}\text{m}{}^3{}_2\text{sin}{}^3\text{i}\text{(m}{}_1\text{m}{}_2\text{)}{}^2\text{=}\text{(a}\text{sin}\text{i)}{}^3\text{E}{}^2{}_0\text{=6 x 10}{}^{\mathrm{?5}}\text{M}{}_{\mathrm{\odot }}$

     

Electric fields and potential patterns in the high-latitude ionosphere for different situations in interplanetary space

The potential ? of the electric field at high latitudes has been obtained by solving numerically the second order differential equation in spherical coordinates:

$\text{?}\text{1}\text{2}\text{(rσ}{}^{\mathrm{H}}\text{?}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\theta }}\text{+}\text{1}\text{r}\text{(σ}{}^{\mathrm{H}}\text{?λ)}{}_{\mathrm{\lambda }}\text{+}\text{1}\text{r}\text{(σ}{}^{\mathrm{P}}\text{?}{}_{\mathrm{\lambda }}\text{)}{}_{\mathrm{\theta }}\text{?(σ}{}^{\mathrm{P}}\text{?}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\lambda }}\text{=}\text{1}\text{r}\text{(rψ}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\theta }}\text{+}\text{1}\text{r}{}^2\text{ψ}{}_{\mathrm{\lambda \lambda }}$

, where θ is colatitude, λ is longitude, σ^H and σ^P are the height-integrated Hall and Perdersen ionospheric conductivities, r = sinθ, and ψ is the current function. The boundary condition is ? = 0 on the geomagnetic parallel θ = 34°. Values of ψ are determined from geomagnetic field variations at the Earth's surface from geomagnetic field variations at the Earth's surface for various conditions in interplanetary space. σ^P and σ^H are taken to vary with season, local time, tilt of the geomagnetic dipole axis (UT), and intensity of corpuscular precipitation (the model proposed by Wallis and Budzinski, 1981). The model distributions of ?^M and E^M = -▽?^m so obtained are compared with observational results. The feasibility has been demonstrated of interpreting the statistical results and individual measurement data in terms of a unified dynamic model of ionospheric electric fields. The model makes allowance for the changes of electromagnetic “weather” in interplanetary space.     

On fourier-analysis of geomagnetic pulsations pc-1, “two-fold” and “three-fold” pulsations pc-1 and fundamental frequencies of the magnetospheric cavity—I

The paper gives the results of detailed studies of the frequency spectra $\text{S}{}_{\mathrm{s}}\text{(?)}$ of the chain of the wave packets F_s(t) of geomagnetic pulsations PC-1 recorded at the Novolazarevskaya station. The bulk of the energy of F_s(t) is concentrated in the vicinity of the central frequencies $\text{?}{}_{\mathrm{s0}}$ of spectra—the carrier frequencies of the signals. The velocity $\text{V}{}_0\text{≌ 6.10}{}^3\text{km s}{}^{\mathrm{?1}}$ of the flux of protons generating these signals correspond to them. The spectra of the signals have oscillations—“satellites” irregularly distributed in frequency. These satellites, as the authors believe, testify to the presence of the individual groups of protons of low concentration whose velocities vary within 10³–10⁴ km s^?1.Their energy is only of the order of 10^?2–10^?3 of the energy of the main proton flux. Clearly pronounced maxima on double and triple frequencies $\text{? = 2?}{}_{\mathrm{s0}}\text{and}\text{3?}{}_{\mathrm{s0}}$ are detected. They show that the generation of pulsations PC-1 is accompanied by the generation on the overtones of wave packets called in this paper “two-fold” and “three-fold” pulsations PC-1. Intensive symmetrical satellites of a modulation character have been discovered on frequencies $\text{?}{}^{\mathrm{\pm }}{}_{\mathrm{sK}}$. Frequency differences $\text{Δ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{=}\text{¦?}{}_{\mathrm{s0}}\text{? ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{¦}\text{= (0.011,0.022}\text{and}\text{0.035)}\text{Hz}$ correspond to them. The authors believe that the values of Δ?^±_sK are resonance frequencies of the magnetospheric cavity in which geomagnetic pulsations PC-1 are generated. It is established that the values of Δ?^±_sK coincide closely with the carrier frequencies of geomagnetic pulsations PC-3 and PC-4 generated in the magnetosphere. This leads to the conclusion that the resonance oscillations of the magnetospheric cavity are their source. Thus, the generation of geomagnetic pulsations of different types and resonance oscillations in the magnetosphere are integrated into a unified process. The importance of the results obtained and the necessity to check further their trustworthiness and universality, using experimental data gathered in different conditions, is stressed.     

On the photodissociation of water vapour in the mesosphere

The photodissociation of water vapour in the mesosphere depends on the absorption of solar radiation in the region (175–200 nm) of the O₂ Schumann-Runge band system and also at H-Lyman alpha. The photodissociation products are OH + H, OH + H, O + 2H and H₂ + O at Lyman alpha; the percentages for these four channels are 70, 8, 12 and 10%, respectively, but OH + H is the only channel between 175 and 200 nm. Such proportions lead to a production of H atoms corresponding to practically the total photodissociation of H₂O, while the production of H₂ molecules is only 10% of the H₂O photodissociation by Lyman alpha.The photodissociation frequency (s^?1) at Lyman alpha can be expressed by a simple formula

$\text{J}{}_{\mathrm{<mtext>Ly</mtext><mtext>\alpha </mtext>}}\text{H}{}_2\text{O=4.5 ×10}{}^{\mathrm{?6}}\text{1+0.2}\text{F}{}_{10.7}\text{?65}\text{100}\text{exp}\text{[?4.4 ×10}{}^{\mathrm{?19}}\text{N}{}^{0.917}\text{]}$

where F_{10.7 cm} is the solar radioflux at 10.7 cm and N the total number of O₂ molecules (cm^?2), and when the following conventional value is accepted for the Lyman alpha solar irradiance at the top of the Earth's atmosphere ($\text{Δλ = 3.5}\text{A}\text{?}$) q_Lyα,∞ = 3 × 10¹¹ photons cm^?2 s^1?.The photodissociation frequency for the Schumann-Runge band region is also given for mesospheric conditions by a simple formula

$\text{J}{}_{\mathrm{<mtext>SRB</mtext>}}\text{(}\text{H}{}_2\text{O}\text{) = J}{}_{\mathrm{<mtext>SRB</mtext><mtext>,\infty </mtext>}}\text{(}\text{H}{}_2\text{O}\text{) exp [?10}{}^{\mathrm{?7}}\text{N}{}^{0.35}\text{]}$

where J_SRB,∞(H₂O) = 1.2 × 10^?6 and 1.4 × 10^?6 s^?1 for quiet and active sun conditions, respectively.The precision of both formulae is good, with an uncertainty less than 10%, but their accuracy depends on the accuracy of observational and experimental parameters such as the absolute solar irradiances, the variable transmittance of O₂ and the H₂O effective absorption cross sections. The various uncertainties are discussed. As an example, the absolute values deduced from the above formulae could be decreased by about 25-20% if the possible minimum values of the solar irradiances were used.     

The excitation of atomic oxygen to the O(1S) level by energy transfer from N2(A3∑u+) molecules in aurora

The part that the energy transfer reaction $\text{N}{}_2\text{(A}{}^3\text{∑}{}_{\mathrm{u}}{}^{\mathrm{+}}\text{) + O(}{}^3\text{P) → N}{}_2\text{(X}{}^1\text{∑}{}_{\mathrm{g}}{}^{\mathrm{+}}\text{) + O(}{}^1\text{S)}$ plays in the excitation of the auroral green line has been investigated. The contribution is estimated to be 40 per cent in this case, containing pulsating aurora in class IBC 1. Due to greater quenching of the A³∑_u⁺ state, the centroid of the VK emission is displaced 10 km upwards of the green line height, which is centred at 110 km.     

Characteristics of localized resonance coupling oscillations of the slow magnetosonic wave in a non-uniform plasma

We analyze linear resonance oscillations in a non-uniform one-fluid finite-β plasma, which is oversimplified to understand easily fundamental characteristics of the resonance oscillations. A linear resonance oscillation of localized slow magnetosonic mode $\text{[ω}{}^2{}_{\mathrm{s}}\text{=}\text{ω}{}^2{}_{\mathrm{A}}\text{(1 +}\text{V}{}^2{}_{\mathrm{A}}\text{V}{}^2{}_{\mathrm{s}}\text{)}\text{]}$, which has the diamagnetic property in a uniform plasma, is newly found to be excited in the radially non-uniform plasma. The localized slow resonance indicates a radially polarized compressional oscillation (δB_∥ ? δB_H ? δB_D). The sense of the Alfvénic polarizations in the H-D plane near the resonant point is a function of both the propagation in the azimuthal direction and the slope of wave amplitude in the radial direction, whereas the sense of the resonant slow magnetosonic polarizations changes in accordance only with the switch in the azimuthal propagation direction. Further multi-satellite studies are necessary to establish the resonant structures of the slow magnetosonic waves in the magnetosphere.     

Analysis of HI content in galaxies

We have collected data on 241 galaxies from 13 sources and made a statistical analysis after reduction to a uniform system. We found that the Hubble sequence is one of increasing $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{〉}$, these mean values increasing monotonically from .0016 and .024 at E to .084 and .83 at Im, but the dispersion is large.The HI content in barred spiral is greater than that in ordinary spirals, and this is consistent with their statistics of angular momentum and colour.The HI content is related to colour; it is greater in bluer systems. The large dispersion suggests that it also depends on some other factors, but these are smoothed out when averaged over each type, resulting in a linear relation between $\text{〈log(}\text{M}{}_{\mathrm{<mtext>M</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and 〈(B ? V^O_T)〉. Unlike the colour-colour diagram, the large dispersion on the $\text{log}\text{(}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{) ? (B ? V}{}^0{}_{\mathrm{<mtext>T</mtext>}}\text{)}$ is not related to peculiar galaxies.     

Comment on the group velocity of surface waves on the plasma sheet boundary

In the recent estimation by Maltsev and Lyatsky (1984) of the group velocity of surface waves on the inner boundary of the plasma sheet, the effect of the curvature of the field lines of the ambient magnetic field of the Earth on the spectrum has been assessed. The authors have not accounted for the fact, however, that the group velocity of the compressional surface magnetohydrodynamic waves itself is nonzero transverse to the magnetic field—a characteristic which has been omitted in the spectrum of Chen and Hasegawa (1974), being used by Maltsev and Lyatsky.This characteristic of compressional surface MHD waves is inherent for the spectrum $\text{ω = (}\text{k}{}_6\text{k}\text{)V}{}_{\mathrm{A}}\text{(k}{}^2{}_6\text{+ 2k}{}^2{}_{\mathrm{\perp }}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, obtained by Nenovski (1978) in the cold plasma limit V_A ? V_S(V_A is Alfvén velocity, and V_S, sound velocity). A comment has been made on the restrictions, proceeding from the approximation, used by Maltsev and Lyatsky. The estimation of the velocities for movements of auroral riometer absorption bays have been reviewed.     

Interplanetary energy flux associated with magnetospheric substorms

It is shown that the interplanetary quantity ε(t), obtained by Perreault and Akasofu (1978), for intense geomagnetic storms, also correlates well with individual magnetospheric substonns. This quantity is given by $\text{ε(t) = VB}{}^2\text{sin}{}^4\text{(}\text{θ}\text{2}\text{)l}{}_{\mathrm{o}}{}^2$, where V and B denote the solar wind speed and the magnitude of the interplanetary magnetic field (IMF), respectively, and θ denotes the polar angle of the IMF; l_o is a constant ? 7 Earth radii. The AE index is used in this correlation study. The correlation is good enough to predict both the occurrence and intensity of magnetospheric substonns observed in the auroral zone, by monitoring the quantity ε(t) upstream of the solar wind.     

Thinning of the near-earth (10∼15 RE) plasma sheet preceding the substorm expansion phase

The timing of the plasma-sheet thinning relative to the onset of the expansion phase of substorms is examined by the analysis of the OGO 5 electron (79 ± 23 keV) and proton (100～150 keV) data with the aid of simultaneous magnetic field observations. It is found that the timing of the thinning is significantly dependent on the distance. At $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ the thinning often starts before the onset, while at $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ it tends to occur after the onset, where x and y refer to solar magnetospheric coordinates. The thinning that precedes the expansion-phase onset has been found to reduce the thickness to ～1 R_E, and further thinning may occur in a spatially limited region. Hence it is conceivable that the formation of the neutral line characterizing the substorm expansion phase is the consequence of the thinning of the plasma sheet in the near-Earth region.     

An indirect mechanism for the production of O(1S) in the aurora

Recent laboratory measurements of the deactivation rate constants for O(¹S) have suggested that the dominant production mechanism for the green line in the nightglow is a two-step process. A similar mechanism involving energy transfer from an excited state of molecular oxygen is considered as a potential source of the OI (5577 Å) emission in the aurora. It is shown that the mechanism, $\text{O}{}_2\text{+ e → O}{}_2{}^1\text{+ e O}{}_2{}^1\text{+ O → O}{}_2\text{+ O(}{}^1\text{S}\text{)}$, is consistent with auroral observations; the intermediate excited state has been tentatively identified as the O₂(c¹∑^?_u) state. For the proposed energy transfer mechanism to be the primary source of the auroral green line, the peak electron impact cross-section for O₂¹ production must be approximately 2 × 10^?17 cm².     

The influence of a flow field on the stability of the force-free magnetic field

In order to analyse the convective instability of the force-free magnetic field, an exact solution of the MHD equation for the magnetic field (1) together with the flow field (2) of constant speed V₀ making an angle θ with the magnetic field, was chosen as the unperturbed state. The stability of the fields between two parallel conducting walls of seperation d was studied by a linear perturbation method, which led to the eigenvalue problem (12), X being given by (13). It was shown by an approximate variational method that instability will set in by the flow field if V₀ is greater than 1/ $\text{3}$ times Alfven velocity V_A. For , the stability of the force-free field (1) is not influenced by the flow field, which may still be significant in other respects. Perturbations transverse to the magnetic field were found to be the most unstable modes.