Line shape of the non-thermal 6300 Å O(1D) emission

The line shape of the non-thermal O(¹D) 6300 Å emission is calculated using the two population model of Schmitt, Abreu and Hays (Planet. Space Sci.29, 1095, 1981). The calculated line shapes simulate observations made from a space platform at different zenith angles and altitudes. The non-thermal line shapes observed at zenith angles other than the local vertical have been obtained by using the Addition theorem for spherical harmonics of a Legendre polynomial expansion of the non-thermal population distribution function.     

On the Origin of Pulsations of Sub-THz Emission from Solar Flares

We propose a model to explain fast pulsations in sub-THz emission from solar flares. The model is based on the approach of a flaring loop as an equivalent electric circuit and explains the pulse-repetition rate, the high-quality factor, Q≥10³, low modulation depth, pulse synchronism at different frequencies, and the dependence of the pulse-repetition rate on the emission flux, observed by Kaufmann et al. (Astrophys. J. 697, 420, 2009). We solved the nonlinear equation for electric current oscillations using a Van der Pol method and found the steady-state value for the amplitude of the current oscillations. Using the pulse rate variation during the flare on 4 November 2003, we found a decrease of the electric current from 1.7×10¹² A in the flare maximum to 4×10¹⁰ A just after the burst. Our model is consistent with the plasma mechanism of sub-THz emission suggested recently by Zaitsev, Stepanov, and Melnikov (Astron. Lett. 39, 650, 2013).     

Anisotropic universe and observational constraint on the dark energy models with the cosmological redshift drift

This study set out to examine the effect of anisotropy on the various dark energy models by using the observational data, including the Sandage-Loeb test, Strongly gravitationally lensing, observational Hubble data, and Baryon Acoustic Oscillations data. In particular, we consider three cases of dark energy models: the cosmological constant model, which is most favored by current observations, the wCDM model where dark energy is introduced with constant w equation of state parameter and in Chevalier-Polarski-Linder parametrization where ω is allowed to evolve with redshift. With an anisotropy framework, a maximum likelihood method to constrain the cosmological parameters was implemented. With an anisotropic universe, we also study the behavior of different cosmological parameters such as Hubble parameter, EoS parameter, and deceleration parameter of dark energy models mentioned. The results indicate that the Bianchi type I model for the dark energy models are consistent with the combined observational data.     

Ionosphere-plasmasheet field-aligned currents and parallel electric fields

Two kinetic models for the auroral topside ionosphere are compared. The collisionless plasma distributed along an auroral magnetic field line behaves like a non-Ohmic conducting medium with highly non-linear characteristic curves relating the parallel current density to the potential difference between the cold ionosphere and the hot plasmasheet region. The (zero-electric current) potential difference, required to balance the current carried by the precipitating plasmasheet particles and the current transported by the outflowing ionospheric particles, depends on the ratio n_ps.e/n_th.e and T_ps.e/T_th.e of the plasmasheet and ionospheric electron densities and temperatures. When in the E-region the magnetic field lines are interconnected by a high conductivity plasma the resulting field-aligned currents driven by the magnetospheric potential distribution are limited by the integrated Pedersen conductivity of the ionospheric layers. These currents are not related to the parallel electric field intensity as they would be in Ohmic materials. The parallel electric field intensity is necessarily determined by the local quasi-neutrality of the plasma.     

A model of cometary gas and dust production and nongravitational forces with application to P/Halley

A program that computes gas and dust production rates and idealized nongravitational force components has been developed and applied to the case of Comet Halley. We use a modified form of our earlier comet model (F.P. Fanale and J.R. Salvali[(1984) Icarus 60, 476–511] to which coma effects and a section on nongravitational forces have been added. The possibility of grain cohesion is also included. These models are used together with observations from 1910 and semiempirically derived data to investigate the effects of obliquity and thermal conductivity of the near thermal conductivity of the nucleus on gas and dust production. The results indicate that the thermal conductivity of the nucleus is of the order of 10⁵ ergs/cm-s-°K, which implies that the ice near the surface is in the crystalline form. A general method is presented for calculating the radii of cometary nuclei using theoretically derived and semiempirically derived nongravitational force components. This method is used to calculate possible radii for Comet Halley that depend on the model variation chosen. The method used and the results presented herein should have greater significance and value when the observational data from Halley's current perihelion passage become available.     

Three-dimensional current system in different phases of a substorm

It is assumed that the three-dimensional current system of a substorm passes three successive stages. (1) When a dawn-to-dusk magnetospheric electric field appears, a current system with field-aligned currents at the poleward boundary of the auroral zone arises. An equivalent ionospheric current system calculated, taking into account a day-night asymmetry of ionospheric conductivity, looks like the well-known DP-2 system including an eastward low-latitude current and a greater magnitude of the dusk vortex in comparison with the dawn one. (2) An electric drift of plasma towards the Earth leads to the appearance of a westward partial ring current increasing in time. This current is closed by field-aligned currents at the equatorward boundary of the auroral zone. The calculated equivalent current system is similar to the well-known one of the precursory phase. (3) An increase of the auroral ionospheric conductivity during the expansive phase produces an increase of all currents and a turning of field-aligned currents at the equatorward boundary of the auroral zone relative to those at the poleward one. The calculated equivalent current system is similar to the DP-1 system.     

A study of the dynamics of a discrete auroral arc

High resolution electric field and particle data, obtained by the S23L1 rocket crossing over a discrete prebreakup arc in January 1979, are studied in coordination with ground observations (Scandinavian Magnetometer Array—SMA, TV and all-sky cameras) in order to clarify the electrodynamics of the arc and its surroundings. Height-integrated conductivities have been calculated from the particle data, including the ionization effects of precipitating protons and assuming a steady state balance between ion production and recombination losses. High resolution optical information of arc location relative to the rocket permitted a check of the validity of this assumption for each flux tube passed by the rocket. Another check was provided by a comparison between calculated (equilibrium values) and observed electron densities along the rocket trajectory. A way to compensate for the finite precipitation time when calculating the electron densities is outlined. The height-integrated HalI-Pedersen conductivity ratio is typically 1.4 within the arc and about 1 at the arc edges, indicative of a relatively softer energy spectrum there. The height-integrated conductivities combined with the DC electric field measurements permitted calculation of the horizontal ionospheric current vectors (J_⊥), Birkeland currents (from div J_⊥) and energy dissipation through Joule heating (Σ_pE²). An eastward current of typically 1 A m^?1 was found to be concentrated mainly to the arc region and equatorward of it. A comparison has been made with the equivalent current system deduced from ground based magnetometer data (SMA) showing a generally good agreement with the rocket results. An intense Pedersen current peak (1.2 A m^?1) was found at the southern arc edge. This edge constituted a division line between a very intense (> 10 μA m^?1) and localized (~ 6 km) downward current sheet to the south, probably carried by upward flowing cold ionospheric electrons and a more extended upward current sheet (> 10 μA m^?2) over the arc carried by measured precipitating electrons. Joule and particle heating across the arc were anticorrelated, consistent with the findings of Evans et al. (1977) with a total value of about 100mW m^?2.     

Dissipation of magnetic flux and magnetic energy in partially ionized gases

Macroscopic equations of motion are used to derive several forms of the generalized Ohm's law for partially ionized ternary gases in magnetic fields, and a conductivity σ is defined that is independent of the magnetic field. A flux theorem is derived using a velocityu _H that can be defined to be the velocity of magnetic field lines;u _H is only slightly different from the velocity of the electron component of the gas. It is shown that σ is the conductivity relevant to the decay of magnetic flux through any surface moving everywhere with velocityu _H . The rate of increase of the thermal energy density of the gas arising through collisions between particles of different species can be resolved into Joule heating at the ratej ²/σ, wherej is the current density, and heating associated with ambipolar drift. The latter, contrary to what has been claimed by some authors, is not necessarily fully compensated by a decrease in the energy of the electromagnetic field. In many applications such compensation does occur, but it may not in interstellar clouds where large amounts of gravitational energy can be made available by collapse, and then both heating and an increase in electromagnetic field energy may occur.     

An electrodynamic model of electric currents and magnetic fields in the dayside ionosphere of Venus

The electric current configuration induced in the ionosphere of Venus by the interaction of the solar wind has been calculated in previous papers (Cloutier and Daniell, Planet. Space Sci. 21, 463, 1973; Daniell and Cloutier. Planet. Space Sci.25, 621, 1977; Cloutier and Daniell, Planet. Space Sci.27, 1111, 1979) for average steady-state solar wind conditions and interplanetary magnetic field. This model is generalized to include the effects of (a) plasma depletion and magnetic field enhancement near the ionopause, (b) velocity-shear-induced MHD instabilities of the Kelvin-Helmholtz type within the ionosphere, and (c) variations in solar wind parameters and interplanetary magnetic field. It is shown that the magnetic field configuration resulting from the model varies in response to changes in solar wind and interplanetary field conditions, and that these variations produce magnetic field profiles in excellent agreement with those seen by the PIONEER-VENUS Orbiter. The formation of “flux-ropes” by the Kelvin-Helmholtz instability is shown to be a natural consequence of the model, with the spatial distribution and size of the flux-ropes determined by the magnetic Reynolds number.     

On high-latitude convection field inhomogeneities,parallel electric fields and inverted-V precipitation events

Assuming a certain horizontal distribution of the convection field at a certain altitude above the ionosphere, the associated electric field and current distributions in a vertical plane are calculated using a model with finite current-dependent conductivity along the magnetic field lines. It is seen that given the kind of horizontal distribution of E₆ commonly observed by polar-orbiting satellites at inverted-V electron precipitation events, the calculated distribution of E_∥ is able to reproduce the basic spatial structure of the precipitation. It is also seen that the combined effect of a locally increased ionization within auroral forms and a large potential difference (ΔV)_∥ along the magnetic field lines at higher altitudes is a strong reduction of E₆ within the auroral forms. From the basic features of the electric field, it is concluded that an interpretation of auroral precipitation in terms of a static E_∥ may require a mechanism that can support a large (ΔV)_∥ even at relatively weak current densities and at the same time allow local enhancements of the parallel conductivity within the region of non-zero E_∥. It is suggested that the magnetic mirroring combined with gyro-resonant wave-particle interactions may be a suitable mechanism.     

Observations of the 22 April 1982 stellar occultation by Uranus and the rings

The 22 April 1982 stellar occultation of KME 14 by Uranus was observed from Tenerife, Canary Islands, using the Teide Observatory 1.5-m telescope. From model fits to the immersion and emersion ring profiles, we obtained (1) midtimes of the ring events with a typical uncertainty of 0.01 sec; (2) ring widths for rings 4, α, β, γ, δ, and ϵ with a typical uncertainty of a few tenths of a kilometer; and (3) equivalent widths and normal optical depths of all nine rings. The immersion planetary occultation was clouded out, but emersion was successfully observed, and the stratospheric temperature profile was obtained by numerical inversion. The profile shows a temperature maximum near the 8-μbar pressure level, characterized by T(8 μbar) = 141°K and T(8 μbar)–T(13 μbar) = 5°K for the sampled suboccultation latitude of −11°.9. Both the mean temperature and the temperature variations are consistent with the latitude-dependent atmospheric structure found by B. Sicardy et al. (1985, Icarus 64, 88–106) from widely separated observations of the same event.     

Magnetic field and field-aligned current variations in the polar cusp controlled by interplanetary medium parameters

The distribution of large-scale field-aligned currents in the dayside sector of the auroral oval has been presented for different situations in the interplanetary space. The j_∥ distribution has been calculated on the basis of a model, each part of which is controlled by a corresponding parameter of the interplanetary space. It has been shown that the field-aligned current models, proposed by Iijima and Potemra and by McDiarmid et al. describe the planetary j_∥ distribution for only particular situations in the interplanetary space and represent some particular cases of a more general model.     

Numerical simulation of cometary nuclei: III. Internal temperatures of cometary nuclei

A one-dimensional heat-diffusion model was used to calculate internal temperatures in cometary nuclei composed of either crystalline or amorphous ice, and for a range of orbits. It was found that the final central temperature, T_c, was a complex function of the comet's orbital semimajor axis, a, and eccentricity, e, as well as the functional form of the thermal conductivity. For cometary nuclei with identical thermal properties, T_c was found to decrease with eccentricity for a short-period orbit with a = 3 AU. For an intermediate-period orbit with a = 20 AU, T_c initially increased with eccentricity but then declined at large values of e for a crystalline ice nucleus, while for amorphous ice T_c increased monotonically. In addition, it was found that for conductivities of similar magnitude, crystalline ice (for which the conductivity varies inversely proportional to temperature) reached the final central temperature twice as fast as amorphouslike ice (for which the conductivity is proportional to temperature). T_c also depended on the magnitude of the conductivity. A four- to fivefold decrease in the conductivity resulted in a 3–4°K decrease in T_c at large eccentricities, while at small eccentricities T_c was only weakly dependent on the conductivity. Finally, the numerical results are compared to the analytical solutions of J. Klinger (1981, Icarus 47, 320–324) and C. P. McKay, S. W. Squyres, and R. T. Reynolds (1986, Icarus, 66, 625–629), and a numerical correction factor is derived for the McKay et al. expression for the central temperature.     

Field-aligned and ionospheric currents

Zmuda and Armstrong (1974) showed that the field-aligned currents consist of two pairs; one is located in the morning sector and the other in the evening sector. Our analysis of magnetic records from the TRIAD satellite suggests that in each pair the poleward field-aligned current is more intense than the equatorward current, a typical ratio being 2:1. This difference has a fundamental importance in understanding the coupling between the magnetosphere and the ionosphere. We demonstrate this importance by computing the ionospheric current distribution by solving the continuity equation $\text{▽ .}\text{I}\text{= j}{}_{\mathrm{\parallel }}$ using the “observed” distribution of j_∥ for several models of the ionosphere with a high conductive annular ring (simulating the auroral oval).It is shown that the actual field-aligned and ionospheric current system is neither a simple Birkeland type, Boström type nor Zmuda-Armstrong type, but is a complicated combination of them. The relative importance among them varies considerably, depending on the conductivity distribution, the location of the peak of the field-aligned currents, etc. Further, it is found that the north-south segment of ionospheric current which connects the pair of the field-aligned currents in the morning sector does not close in the same meridian and has a large westward deflection. Thus, it has an appreciable contribution to the westward electrojet. One of the model calculations shows that the entire north-south closure current contributes to the westward electrojet.     

Evolution of Planetesimal Velocities Based on Three-Body Orbital Integrations and Growth of Protoplanets

We obtain the viscous stirring and dynamical friction rates of planetesimals with a Rayleigh distribution of eccentricities and inclinations, using three-body orbital integration and the procedure described by Ohtsuki (1999, Icarus137, 152), who evaluated these rates for ring particles. We find that these rates based on orbital integrations agree quite well with the analytic results of Stewart and Ida (2000, Icarus 143, 28) in high-velocity cases. In low-velocity cases where Kepler shear dominates the relative velocity, however, the three-body calculations show significant deviation from the formulas of Stewart and Ida, who did not investigate the rates for low velocities in detail but just presented a simple interpolation formula between their high-velocity formula and the numerical results for circular orbits. We calculate evolution of root mean square eccentricities and inclinations using the above stirring rates based on orbital integrations, and find excellent agreement with N-body simulations for both one- and two-component systems, even in the low-velocity cases. We derive semi-analytic formulas for the stirring and dynamical friction rates based on our numerical results, and confirm that they reproduce the results of N-body simulations with sufficient accuracy. Using these formulas, we calculate equilibrium velocities of planetesimals with given size distributions. At a stage before the onset of runaway growth of large bodies, the velocity distribution calculated by our new formulas are found to agree quite well with those obtained by using the formulas of Stewart and Ida or Wetherill and Stewart (1993, Icarus106, 190). However, at later stages, we find that the inclinations of small collisional fragments calculated by our new formulas can be much smaller than those calculated by the previously obtained formulas, so that they are more easily accreted by larger bodies in our case. The results essentially support the previous results such as runaway growth of protoplanets, but they could enhance their growth rate by 10-30% after early runaway growth, where those fragments with low random velocities can significantly contribute to rapid growth of runaway bodies.     

A quasi-self-consistent axially symmetric model for the growth of a ring current through earthward motion from a pre-storm configuration

A quasi-self-consistent axially symmetric model of a storm-time ring current that has evolved quasi-adiabatically through inward motion from a prestorm state is presented in which the disturbance field as a function of geocentric distance, r, in the equatorial plane (including the value, H, at r = 0 from which Dst can be found), the beta of the plasma as a function of r, and the ring current magnetic moment are all given in terms of analytic expressions, having only two independent ring current parameters: the geocentric distance to the inner edge of the ring current, R, and the distance at which the beta of the plasma is unity, k—a constant of the model. The model is used to find H as a function of R at constant k, which corresponds to the growth of Dst as the ring current moves earthward, and to illustrate how the radial profile of the disturbance field varies also as the ring current moves earthward. Comparison with published S³-A storm-time, particle and field data shows reasonable agreement. The model is also able to fit the empirically determined relationship between Dst and the equatorward extension of quiet auroral arcs. This result applied to the great magnetic storm of September 1859 predicts a Dst of the order of 2000 nT, a value for which the sparse magnetometer data for that storm offer some support. Because the model (quasi) self-consistently determines the particle distribution function by allowing it to (quasi) adiabatically evolve in a field of its own making, the model avoids difficulties encountered by previous ring current models in which the particle distribution function had to be assumed. For example, the ring current magnetic moment is not limited to values less than the geomagnetic dipole moment, and null points in the field do not develop even in very high beta situations.     

Atmospheric odd oxygen production due to the photodissociation of ordinary and isotopic molecular oxygen

Through a line by line calculation, the contributions of the Schumann-Runge bands of the ordinary and isotopic oxygen to the photodissociation of these molecules at different altitudes have been calculated. The photodissociation rates are expressed analytically. Contribution of the satellite lines has been taken into account. Due to the broadening of the SR lines, this contribution is insignificant. Similarly, it is shown that the first and higher vibrational states of the initial molecular states contribute insignificantly to the dissociation rates. It is also shown that the main contribution to the odd oxygen production in the important ozone producing altitudes is from the low vibrational and high rotational quantum numbers. The effect of the temperature on dissociation rates has similarly been studied.Due to its selective absorption, the isotopic oxygen ¹⁶O¹⁸O produces at 70 km 10 times as much odd oxygen as would be produced if the isotope did not have selective absorption. At this altitude 6% of the odd oxygen produced is due to this isotope. Also, 1.45% of the odd oxygen produced per second in an atmospheric column is due to ¹⁶O¹⁸O. However, the excess odd oxygen produced is not enough to explain the excess amount of ozone observed in the atmosphere which cannot be accounted for in the photochemical models.The calculated dissociation rates for the isotope are in moderate agreement with similar rates obtained by Blake et al. (1984, J. geophys. Res.89, 7277), but are by an order of magnitude smaller than similar rates given by Cicerone and McCrumb (1980, Geophys. Res. Lett.7, 251).     

Unipolar induction in the moon and a lunar limb shock mechanism

The unipolar induction mechanism is employed to calculate electric field profiles in the interior of a chemically homogeneous Moon possessing a steep radial thermal gradient characteristic of long-term radioactive heating. The thermal models used are those of Fricker, Reynolds, and Summers. From the magnetic field, the magnetic back pressure upon the solar wind is found. The electric field profile is shown to depend only upon the activation energy,E _o, of the geological material and the radial gradient of the reciprocal temperature. The current is additionally dependent upon the coefficient of the electrical conductivity function but only by a scale factor. Since the Moon is experimentally known to correspond to the case of weak interaction with the solar wind, the magnetic back pressure is calculated without the need for an iterative procedure. The results indicate that a hot Moon can yield sufficient current flow so that the magnetic back pressure is observable as a vestigial limb shock wave using an activation energy of about 2/3 eV together with a conductivity coefficient of about 10³ mhos/m. Such matter is approximated by diabase-like composition, although the result that both the activation energy and coefficient enter into the current determination does not rule out the possibility of a match with other similar substances. The calculations are entirely consistent with earlier results which indicated a model where the unipolar current density is dominated by a high impedance surface layer and a strong shock wave is inhibited. In addition to the magnetic back pressure, the integration of the current continuity equation permits current densities and joule heating rates to be calculated, though the magnitude of the latter for present solar wind conditions is not thermally important.On leave from NASA Ames Research Center     

Measurement of Impact Ejecta from Regolith Targets in Oblique Impacts

This paper reports the results of experiments on projectile impact into regolith targets at various impact angles. Copper projectiles of 240 mg are accelerated to 197 to 272 m s⁻¹ using an electromagnetic gun. The ejecta are detected by thin Al foil targets as secondary targets, and the resulting holes on the foil are measured to derive the spatial distribution of the ejecta. The ejecta that penetrated the foil are concentrated toward the downrange azimuths of impacting projectiles in oblique impacts. In order to investigate the ejecta velocity distribution, the nondimensional volume of ejecta with velocities higher than a given value is calculated from the spatial distribution. In the case of the vertical impact of the projectile, most ejecta have velocities lower than 24% of the projectile speed (∼50 m s⁻¹), and there are only several ejecta with velocities higher than 72 m s⁻¹. This result confirms the existence of an upper limit to the ejection velocity in the ejecta velocity distribution (Hartmann cutoff velocity) (W. K. Hartmann, 1985, Icarus63, 69-98). On the other hand, it is found that, in the oblique impacts, there are a large number of ejecta with velocities higher than the Hartmann cutoff velocity. The relative quantity of ejecta above the Hartmann cutoff velocity increases as the projectile impact angle decreases. Taking these results with the results of S. Yamamoto and A. M. Nakamura (1997, Icarus128, 160-170) from impact experiments using an impact angle of 30°, it can be concluded that the ejecta from these regolith targets exhibit a bimodal velocity distribution. Below a few tens of m s⁻¹, we see the expected velocity distribution of ejecta, but above this velocity we see a separate group of high-velocity ejecta.