Weibel Turbulence in Laboratory Experiments and GRB/SN Shocks

It has recently been realized that the Weibel instability plays a major role in the formation and dynamics of astrophysical shocks of gamma-ray bursts and supernovae. Thanks to technological advances in the recent years, experimental studies of the Weibel instability are now possible in laser-plasma interaction devices. We, thus, have a unique opportunity to model and study astrophysical conditions in laboratory experiments – a key goal of the Laboratory Astrophysics program. Here we briefly review the theory of strong non-magnetized collisionless GRB and SN shocks, emphasizing the crucial role of the Weibel instability and discuss the properties of radiation emitted by (isotropic) electrons moving through the Weibel-generated magnetic fields, which is referred to as the jitter radiation. We demonstrate that the jitter radiation field is anisotropic with respect to the direction of the Weibel current filaments and that its spectral and polarization characteristics are determined by microphysical plasma parameters. We stress that the spectral analysis can be used for accurate diagnostics of the plasma conditions in laboratory experiments and in astrophysical GRB and SN shocks.     

Irregular Satellites of Jupiter: a study of the capture direction

The irregular satellites of Jupiter are believed to be captured asteroids or planetesimals. In the present work is studied the direction of capture of these objects as a function of their orbital inclination. We performed numerical simulations of the restricted three-body problem, Sun-Jupiter-particle, taking into account the growth of Jupiter. The integration was made backward in time. Initially, the particles have orbits as satellites of Jupiter, which has its present mass. Then, the system evolved with Jupiter losing mass and the satellites escaping from the planet. The reverse of the escape direction corresponds to the capture direction. The results show that the Lagrangian points L1 and L2 mainly guide the direction of capture. Prograde satellites are captured through these two gates with very narrow amplitude angles. In the case of retrograde satellites, these two gates are wider. The capture region increases as the orbital inclination increases. In the case of planar retrograde satellites the directions of capture cover the whole 360° around Jupiter. We also verified that prograde satellites are captured earlier in actual time than retrograde ones. This paper was presented at the Asteriods, Comets and Meteors meeting held at Búzios, Rio de Janeiro, Brazil in August 2005 and could not be included in the special issue related to that conference.     

Ion acceleration and Alfvén wave generation at the Earth’s bow shock

The processes of ion acceleration and Alfvén wave generation by accelerated particles at the Earth’s bow shock are studied within a quasi-linear approach. Steady-state ion and wave spectra are shown to be established in a time of 0.3–4 h, depending on the background level of Alfvénic turbulence in the solar wind. The Alfvén waves produced by accelerated ions are confined within the frequency range 10^?2–1 Hz and their spectral peak with a wave amplitude βB ～ B comparable to the interplanetary magnetic field strength B corresponds to the frequency v = (2–3) × 10^?2 Hz. The high-frequency part of the wave spectrum (v > 0.2 Hz) undergoes damping by thermal ions. The calculated spectra of the accelerated ions and the Alfvén waves generated by them reproduce the main features observed in experiments.     

A Two—Temperature Supernova Fallback Disk Model for Anomalous X—ray Pulsars

We present a case study of the relevance of the radially pulsational instability of a two-temperature accretion disk around a neutron star to anomalous X-ray pulsars (AXPs). Our estimates are based on the approximation that such a neutron star disk with mass in the range of 10^-6-10^-5M⊙ is formed by supernova fallback. We derive several peculiar properties of the accretion disk instability: a narrow interval of X-ray pulse periods; lower X-ray luminosities; a period derivative and an evolution time scale. All these results are in good agreement with the observations of the AXPs.     

Validating the Flash Code: Vortex-Dominated Flows

As a component of the Flash Center’s validation program, we compare FLASH simulation results with experimental results from Los Alamos National Laboratory. The flow of interest involves the lateral interaction between a planar M _a = 1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SF₆) in air, and in particular the development of primary and secondary instabilities after the passage of the shock. While the overall evolution of the flow is comparable in the simulations and experiments, small-scale features are difficult to match. We focus on the sensitivity of numerical results to simulation parameters.     

Studying Hydrodynamic Instability Using Shock-Tube Experiments

The hydrodynamic instability, which develops on the contact surface between two fluids, has great importance in astrophysical phenomena such as the inhomogeneous density distribution following a supernova event. In this event acceleration waves pass across a material interface and initiate and enhance unstable conditions in which small perturbations grow dramatically. In the present study, an experimental technique aimed at investigating the above-mentioned hydrodynamic instability is presented. The experimental investigation is based on a shock-tube apparatus by which a shock wave is generated and initiates the instability that develops on the contact surface between two gases. The flexibility of the system enables one to vary the initial shape of the contact surface, the shock-wave Mach number, and the density ratio across the contact surface. Three selected sets of shock-tube experiments are presented in order to demonstrate the system capabilities: (1) large-initial amplitudes with low-Mach-number incident shock waves; (2) small-initial amplitudes with moderate-Mach-number incident shock waves; and (3) shock bubble interaction. In the large-amplitude experiments a reduction of the initial velocity with respect to the linear growth prediction was measured. The results were compared to those predicted by a vorticity-deposition model and to previous experiments with moderate- and high-Mach number incident shock waves that were conducted by others. In this case, a reduction of the initial velocity was noted. However, at late times the growth rate had a 1/t behavior as in the small-amplitude low-Mach number case. In the small-amplitude moderate-Mach number shock experiments a reduction from the impulsive theory was noted at the late stages. The passage of a shock wave through a spherical bubble results in the formation of a vortex ring. Simple dimensional analysis shows that the circulation depends linearly on the speed of sound of the surrounding material and on the initial bubble radius.     

Homoclinic and heteroclinic orbits in the photogravitational restricted three-body problem

We study numerically the asymptotic homoclinic and heteroclinic orbits around the hyperbolic Lyapunov periodic orbits which emanate from Euler's critical points L ₁ and L ₂, in the photogravitational restricted plane circular three-body problem. The invariant stable-unstable manifolds associated to these Lyapunov orbits, are also presented. Poincaré surface of sections of these manifolds on appropriate planes and several homoclinic and heteroclinic orbits for the gravitational case as well as for varying radiation factor q ₁, are displayed. Homoclinic-homoclinic and homoclinic-heteroclinic-homoclinic chains which link the interior with the exterior Hill's regions, are illustrated. We adopt the Sun-Jupiter system and assume that only the larger primary radiates. It is found that for small deviations of its value from the gravitational case (q ₁ = 1), the radiation pressure exerts a significant impact on the Hill's regions and on these asymptotic orbits.     

NIGHTGLOW: an instrument to measure the Earth’s nighttime ultraviolet glow—results from the first engineering flight

We have designed and built an instrument to measure and monitor the “nightglow” of the Earth’s atmosphere in the near ultraviolet (NUV). In this paper we describe the design of this instrument, called NIGHTGLOW. NIGHTGLOW is designed to be flown from a high altitude research balloon, and circumnavigate the globe. NIGHTGLOW is a NASA, University of Utah, and New Mexico State University project. A test flight took place from Palestine, Texas on July 5, 2000, lasting about 8 h. The instrument performed well and landed safely in Stiles, Texas with little damage. The resulting measurements of the NUV nightglow are compared with previous measurements from sounding rockets and balloons.     

Numerical Stabilization of Orbital Motion

Mainly, the author focuses on Baumgarte's method and its applications in satellite, asteroid, stellar and planetary problems. In the paper arguments are given for the use of energy relations for stabilization in the elliptical two-body problem. Stabilizing properties of Baumgarte's equations and others are discussed. A simple approach is proposed for stabilizing the equations of almost circular motion. By using Baumgarte's technique, the author derives stabilized equations of perturbed restricted three-body problem. It is shown experimentally that stabilization in the problems mentioned above can raise the accuracy of numerical integration by several orders.     

A method for averaging the perturbing function in Hill’s problem

A modified method for averaging the perturbing function in Hill’s problem is suggested. The averaging is performed in the revolution period of the satellite over the mean anomaly of its motion with a full allowance for a variation in the position of the perturbing body. At its fixed position, the semimajor axis of the satellite orbit during the revolution of the satellite is constant in view of the evolution equations, while the remaining orbital elements undergo secular and long-period perturbations. Therefore, when the motion of the perturbing body is taken into account, the semimajor axis of the satellite orbit undergoes the strongest perturbations. The suggested approach generalizes the averaging method in which only the linear (in time) term is included in the perturbing function. This method requires no expansion in powers of time. The described method is illustrated by calculating the perturbations of the semimajor axes for two distant satellites of Saturn, S/2000 S 1 and S/2000 S5. An approximate analytic solution is compared with the results of numerical integration of the averaged system of equations of motion for these satellites.