Investigation of shock waves and tangential discontinuities of the cosmic plasma by the radio interference technique

The possibility of investigation of the cosmic plasma dynamics by the radio interference technique based on a finite time of radio wave propagation between the sounding and responding stations is shown. By locating the sounding station on a spacecraft the greatest contribution to the phase difference ΔΦ(t)or the phase difference growth rate $\text{Δ?}$ between the sounding and response signals are caused by disturbances in close proximity to the spacecraft. This method permits interplanetary shock waves and tangential discontinuities to be registered and the velocities and plasma density changes on their fronts to be determined. By using experimental data of ΔΦ(t) or $\text{Δ?(t)}$ one can also obtain information about plasma concentration jump, location and motion of bow shock wave and magnetopause and plasmapause. Available experimental data about different disturbances of cosmic plasma were analysed and the requirements on frequency stability of spacecraft-borne and groundbased radio equipment were estimated to register those disturbances. In most cases relative stability 10^?11–10^?13 provided by present atomic frequency standards is sufficient.     

Ring current simulation in connection with interplanetary space conditions

A ring current model has been obtained which permits calculations ofD_st variations on the Earth's surface during magnetic storms. The changes in D_st are described by the equation

$\text{d}\text{d}\text{t}\text{D}{}_{\mathrm{sto}}\text{= F(EM)?}\text{D}{}_{\mathrm{sto}}\text{tau;}$

where $\text{D}{}_{\mathrm{sto}}\text{= D}{}_{\mathrm{st}}\text{-bp}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{+~tc}$; p = mnv² is solar wind pressure; F(EM) is the function, controlled by the electromagnetic parameters of interplanetary medium, of injection into ring current ; τ is the constant of ring current decay. $\text{C}\text{= C}\text{u}\text{τ?=18}\text{nT}$, where C is the level of the D_st-variation field measurements; ? is the injection function characterizing the quasisteady-state injection of energy into the ring-current region. The constant Ç is determined from the condition that the change of the ring current energy from magnetic storm commencement to end should equal the difference between the injected and dissipated energy throughout the storm. The values of the factors b and τ were found by the method of minimizing the sum of the quadratic deviations of the calculated D_st from the values observed throughout the storm : b = 0.23 nT/(eV cm^?3)^{$\text{1}\text{2}$}, τ = 8.2 h at D_st? ? 55 nT and τ = 5.8 h at -120 ? D_st ? — 55 nT. The injection function F(EM) is of the form F(EM) = d(E_y? A) at the values of the azimuthal component of the solar wind electric field E_y ? A, and F(EM) =0 at A?E_y.d = ? 1.2 × 10^?3 Ts^?1 (mV/m)^?1 and A = ? 0.9 mV m^?1 . Thus, the injection to ring current is possible at the northward B_z component of the IMF.     

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.     

Determination of magnetic field direction in a comet's head

The Stokes parameters of resonance radiation scattered by a Na atom with the angular momentum F aligned by directed unpolarized radiation in a magnetic field H ～ 10^?5?10^?1 Oe are presented. An influence of the orientation of the magnetic field on these parameters are studied; the intensity ratio $\text{I(D}{}_2\text{)}\text{I(D}{}_1\text{)}$ changes within ±5%, and the polarization degree P(D₂) within ±25%. Measurements of $\text{I(D}{}_2\text{)}\text{I(D}{}_1\text{)}$ and P(D₂), if the geometry of scattering is known, may give information on the direction of the magnetic field in the sodium atmospheres of comets, as well as Io's sodium cloud or man-made cosmic clouds.     

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.     

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.     

Venus: Cloud optical depth and surface albedo from Venera 8

We have used the measurements of the solar flux obtained by the Venera 8 spacecraft inside the atmosphere of Venus and the values of the Venus spherical albedo to deduce the characteristics of the clouds and of the ground. The method used is the exponential kernel approximation and the results have been tested by exact computations with the spherical harmonics method.A cloud layer with an optical thickness $\text{τ}{}_1\text{? 144}$, an albedo for single scattering $\text{ω}{}_0\text{= 0.9998}$ in the rear infrared, above a Rayleigh layer between 0 and 32 km and a ground of reflectivity ? = 0.4, gives a good agreement with the experimental results. A model with two cloud layers is also discussed.     

Application of methane band-model parameters to the visible and near-infrared spectrum of Uranus

Laboratory band-model absorption coefficients of CH₄ are used to calculate the Uranus spectrum from 5400 to 10,400 Å. A good fit of both strong and weak bands for the Uranus spectrum over the entire wavelength interval is achieved for the first time. Three different atmospheric models are employed: a reflecting layer model, a homogeneous scattering layer model, and a clear atmosphere sandwiched between two scattering layers. The spectrum for the reflecting layer model exhibits serious discrepancies but shows that large amounts of CH₄ (5–10 km-am) are necessary to reproduce the Uranus spectrum. Both scattering models give reasonably good fits. The homogeneous model requires a particle scattering albedo $\text{(}\text{g}\text{?}\text{w}{}_{\mathrm{<mtext>p</mtext>}}\text{) ? 0.998}$ and an abundance per scattering mean free path $\text{(}\text{a}\text{?}\text{)}\text{of}\text{a}\text{?}\text{1}\text{km-am}$. The parameters derived from the sandwich layer model are: forsb the upper scattering layer a continuum single scattering albedo ($\text{g}\text{?}\text{w}{}_0$) of 0.995 and a scattering optical depth variable with wavelength consistent with Rayleigh scattering; for the clear layer they are a CH₄ abundance (a) of 2.2 km-am and an effective pressure (p) ? 0.1 atm; for the lower cloud deck a Lambert reflectivity (L) of 0.9 resulted. A severe depletion of CH₄ in the upper scattering layer is required. An enrichment of CH₄/H₂ over the solar ratio by a factor of 4–14 in the lower atmosphere is, however, indicated.     

Residual mass from atmospheric ablation of small meteoroids

Numerical solutions of the equations of meteor ablation in the Earth's atmosphere have been obtained using a variable step size Runge-Kutta technique in order to determine the size of the residual mass resulting from atmospheric flight. The equations used include effects of meteoroid heat capacity and thermal radiation, and a realistic atmospheric density profile. Results were obtained for initial masses in the range 10^?7–10^?2 g, and for initial velocities less than 24 km s^?1 (results indicated no appreciable residual mass for meteors with velocities above 24 km s^?1 in this mass range). The following function has been obtained to provide the logarithm of the ratio of the residual mass following atmospheric ablation to the original preatmospheric mass

$\text{log r = 4.7 ?0.33v}{}_{\mathrm{\infty }}\text{?0.013v}{}_{\mathrm{\infty }}{}^2\text{+ 1.2 log m}{}_{\mathrm{\infty }}\text{+ 0.08 log}{}^2\text{m}{}_{\mathrm{\infty }}\text{?0.083v}{}_{\mathrm{\infty }}\text{log m}{}_{\mathrm{\infty }}\text{M}$

The pre-atmospheric mass and velocity are represented by m_∞ and v_∞.When the results are expressed in terms of the size of the residual mass following atmospheric ablation as a function of the initial mass and velocity, it is found that the final residual mass is almost independent of the original mass of the meteoroid, but very strongly dependent on the original velocity. For example, the residual mass is very nearly 10^?7 g for a meteoroid with velocity 18 kms^?1 for initial masses from 10^?7 to 10^?3 g. On the other hand, a slight change in the initial velocity to 20 km s^?1 will shift the residual mass to approx. 10^?8 g. This strong velocity dependence coupled with the weak dependence on the original mass has important consequences for the sampling of ablation product micrometeorites.     

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.     

Substorm energy

It is shown that the area A_k(× 10⁶km²) covered by brightest auroras and the area A_q bounded by the auroral oval have a simple relation given by

$\text{A}{}_{\mathrm{k}}\text{= 0.05(A}{}_{\mathrm{q}}\text{? A}{}_0\text{)}{}^2$

, where A₀ denotes the area of the minimum size oval and the quantity (A_q ? A₀)² is proportional to the energy ε_q which is stored in the magnetotail and is available for substorms. Following the definition of the intensity of solar flares, A_k may be chosen as a measure of the intensity of substorms. It is also found that the joule heat energy produced by the auroral electrojet is also proportional to (A_q ? A₀)². Thus, it may be concluded that the intensity of substorms is proportional to the energy ε_q stored in the magnetotail.     

Radiation forces on small particles in the solar system

We present a new and more accurate expression for the radiation pressure and Poynting-Robertson drag forces; it is more complete than previous ones, which considered only perfectly absorbing particles or artificial scattering laws. Using a simple heuristic derivation, the equation of motion for a particle of mass m and geometrical cross section A, moving with velocity v through a radiation field of energy flux density S, is found to be (to terms of order $\text{v}\text{c}$)

$\text{m}\text{v}\text{?}\text{= (}\text{SA}\text{c}\text{)Q}{}_{\mathrm{<mtext>pr</mtext>}}\text{[(1 ?}\text{r}\text{?}\text{c}\text{)}\text{S}\text{?}\text{?}\text{v}\text{c}\text{]}$

, where ? is a unit vector in the direction of the incident radiation, $\text{r}\text{?}$ is the particle's radial velocity, and c is the speed of light; the radiation pressure efficiency factor Q_pr ≡ Q_abs + Q_sca(1 ? 〈cos α〉), where Q_abs and Q_sca are the efficiency factors for absorption and scattering, and 〈cos α〉 accounts for the asymmetry of the scattered radiation. This result is confirmed by a new formal derivation applying special relativistic transformations for the incoming and outgoing energy and momentum as seen in the particle and solar frames of reference. Q_pr is evaluated from Mie theory for small spherical particles with measured optical properties, irradiated by the actual solar spectrum. Of the eight materials studied, only for iron, magnetite , and graphite grains does the radiation pressure force exceed gravity and then just for sizes around 0.1 μm; very small particles are not easily blown out of the solar system nor are they rapidly dragged into the Sun by the Poynting-Robertson effect. The solar wind counterpart of the Poynting-Robertson drag may be effective, however, for these particles. The orbital consequences of these radiation forces-including ejection from the solar system by relatively small radiation pressures-and of the Poynting-Robertson drag are considered both for heliocentric and planetocentric orbiting particles. We discuss the coupling between the dynamics of particles and their sizes (which diminish due to sputtering and sublimation). A qualitative derivation is given for the differential Doppler effect, which occurs because the light received by an orbiting particle is slightly red-shifted by the solar rotation velocity when coming from the eastern hemisphere of the Sun but blue-shifted when from the western hemisphere; the ratio of this force to the Poynting-Robertson force is $\text{(}\text{R}{}_{\mathrm{\odot }}\text{r}\text{)}{}^2\text{[(}\text{w}{}_{\mathrm{\odot }}\text{n}\text{) ? 1]}$, where R_⊙ and w_⊙ are the solar radius and spin rate, and n is the particle's mean motion. The Yarkovsky effect, caused by the asymmetry in the reradiated thermal emission of a rotating body, is also developed relying on new physical arguments. Throughout the paper, representative calculations use the physical and orbital properties of interplanetary dust, as known from various recent measurements.     

Pitch angle dependence of the charge-exchange lifetime of ring current ions

Previous work has parameterized the pitch angle dependence of the charge-exchange lifetime τ of ring current ions in terms of γ, the power of the cosine of the mirror latitude λ_m of the particles, such that $\text{τ(λ}{}_{\mathrm{m}}\text{)}\text{τ(0) ≌}\text{cos}{}^{\mathrm{\gamma }}\text{λ}{}_{\mathrm{m}}$ at given L. Using the atomic hydrogen density model of Johnson and Fish, previous authors have suggested values of γ = 5 or 6. We here evaluate γ as a function of λ_m and L using the more recent Chamberlain density models, and show that γ = 3?4 is more appropriate over most of the pitch angle and L range. Consequently, ion distributions in the ring current decay phase are expected to become rather less anisotropic in pitch angle due to chargeexchange than previously believed. We have also investigated the use of several other simple approximate analytic forms for $\text{τ(λ}{}_{\mathrm{m}}\text{)}\text{τ(0)}$, one of which gives far better agreement with the numerical results than the cos^γ λ_m, variation, and should hence be used in future studies.     

The roles of the north-south component of the interplanetary magnetic field on large-scale auroral dynamics observed by the DMSP satellite

An extensive study of DMSP photographs and the simultaneous interplanetary magnetic field data suggests that the quantity defined by

$\text{S=∫}{}^{\mathrm{\tau }}{}_0\text{(Ф}{}_{\mathrm{D}}\text{? Ф}{}_{\mathrm{N}}\text{)dt}$

has a fundamental importance in substorm processes, where Φ_D and Φ_N denote the production rate of merged (or open) field lines along the dayside X-line and of reconnected (or closed) field lines along the nightside X-line, respectively; t = 0 is measured from the time when the B_z component begins to decrease after a prolonged period of a large positive B_z value. It is shown, first of all, that substorms occur so long as S > 0, regardless of the sign of the B_z component and its changes (namely, the southward and northward turnings) and of its time derivative as well. Secondly, the intensity of substorms is proportional to S². By introducing the quantity S, the recent confusion of the problem of the roles of the north-south component of the interplanetary magnetic field on substorm processes can be removed.Since S is equal to the amount of the open magnetic fluxes at a time reckoned from t = 0, it is proportional to (A₁ ? A₀), where A₀ denotes the minimum polar cap area (namely, the area bounded by the minimum auroral oval) and A₁ the polar cap area at an arbitrary time t. Therefore, substorms can occur whenever the auroral oval is larger than its minimum size. Further, an intense substorm tends to occur along a large oval.The quantity S can also be considered as an excess flux, and thus the substorm can be considered as a process by which the magnetosphere tends to remove sporadically the excess energy associated with S.     

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.     

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.     

Cosmic ray ionization of the Jovian atmosphere

An approximate form of the Boltzmann equation has been used to obtain local ionization rates due to the absorption of galactic cosmic rays in the Jovian atmosphere. It is shown that the muon flux component of the cosmic ray-induced cascade may be especially importannt in ionizing the atmosphere at levels where the total number density exceeds 10¹⁹ cm^?3 (well below the ionospheric layers produced by solar euv). A model containing both positive and negative ion reactions has been employed to compute equilibrium electron and ion number densities. Peak electron number densities on the order of 10³ cm^?3 may be expected even at relatively low magnetic latitudes. The dominant positive ions are NH₄⁺ and C_nH_m⁺ cluster ions, with n ? 2; it is suggested that the absorption of galactic cosmic ray energy at such relatively high pressures in the Jovian atmosphere $\text{(M ? 10}{}^{18}\text{to}\text{10}{}^{20}\text{cm}{}^{\mathrm{?3}}\text{)}$ and the subsequent chemical reactions may be instrumental in the local formation of complex hydrocarbons.     

A possible mechanism of formation of faculae

We propose a new heating mechanism of faculae. We think that the formation of faculae is a result of the Joule dissipation of the Hall current generated by the interaction of the convection field of granules in an active region and the inter-granular magnetic field. For a region to generate effectively Hall current, its characteristic length must be such that the magnetic Reynolds number is less than 1. The equation of energy balance in the facula region is

$\text{16σT}{}^3{}_{\mathrm{p}}\text{(T}{}_{\mathrm{l}}\text{? T}{}_{\mathrm{p}}\text{)}{}^{\mathrm{n}}\text{HP}{}_{\mathrm{s}}{}^{\mathrm{a}}\text{H}{}^{\mathrm{?}}\text{=}\text{Qn}{}_{\mathrm{s}}\text{m}{}_{\mathrm{i}}\text{u}{}_{\mathrm{x}}{}^2\text{(τ}{}^2{}_{\mathrm{in}}\text{ω}{}_{\mathrm{i}}\text{)}$

.For five observational models of faculae, we calculated the corresponding velocity fields, and the results are in basic agreement with the observed fields. The present mechanism explains the dependence of the facula brightness on the magnetic and velocity fields, the apparent distribution of the faculae on the solar disk and suggest a possible interpretation of the five structures of faculae.