Optical Flash of GRB 990123: Constraints on the Physical Parameters of the Reverse Shock

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best-fit physical parameters for GRB 990123 and the assumption that the parameters in the optical flash are the same as in the afterglow, we show that: 1) the shell is thick rather than thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3×10~-4 of the observed. In order to remove this discrepancy, the physical parameters, the electron energy and magnetic ratios, εe and εB, should be 0.61 and 0.39, which are very different from their values for the late afterglow.     

General relativistic self-similar solutions in cosmology

We present general relativistic solutions for self-similar spherical perturbations in an expanding cosmological background of cold pressure-less gas. We focus on solutions having shock discontinuities propagating in the surrounding cold gas. The pressure, p , and energy density, μ, in the shock-heated matter are assumed to obey   p = w μ  , where w is a positive constant. Consistent solutions are found for shocks propagating from the symmetry centre of a region of a positive density excess over the background. In these solutions, shocks exist outside the radius marking the event horizon of the black hole which would be present in a shock-less collapse. For large jumps in the energy density at the shock, a black hole is avoided altogether and the solutions are regular at the centre. The shock-heated gas does not contain any sonic points, provided the motion of the cold gas ahead of the shock deviates significantly from the Hubble flow. For shocks propagating in the uniform background, sonic points always appear for small jumps in the energy density. We also discuss self-similar solutions without shocks in fluids with   w < −1/3  .     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

Partial Reduction in the N-Body Planetary Problem using the Angular Momentum Integral

We present a new set of variables for the reduction of the planetary n-body problem, associated to the angular momentum integral, which can be of any use for perturbation theory. The construction of these variables is performed in two steps. A first reduction, called partial is based only on the fixed direction of the angular momentum. The reduction can then be completed using the norm of the angular momentum. In fact, the partial reduction presents many advantages. In particular, we keep some symmetries in the equations of motion (d'Alembert relations). Moreover, in the reduced secular system, we can construct a Birkhoff normal form at any order. Finally, the topology of this problem remains the same as for the non-reduced system, contrarily to Jacobi's reduction where a singularity is present for zero inclinations. For three bodies, these reductions can be done in a very simple way in Poincaré's rectangular variables. In the general n-body case, the reduction can be performed up to a fixed degree in eccentricities and inclinations, using computer algebra expansions. As an example, we provide the truncated expressions for the change of variable in the 4-body case, obtained using the computer algebra system TRIP.     

Analysis of the Space Debris Impacts Risk on the International Space Station

Öpik's analytical expressions relate in a simple way the semimajor axis, eccentricity and inclination of the projectile orbit to the magnitude and direction of the relative velocity vector at impact on a given target on circular orbit. These interesting quantities, along with the impact probability of any given projectile, can be all represented on a suitable projection giving a comprehensive picture of the impact risk on the selected target. By means of this theory a complete analysis of the impact risk on the International Space Station (ISS) is performed. It is found that the large majority of the debris population is on orbits such that a correlation exists between their impact velocity on the ISS and the angle between the velocity vector of the impactor and that of the ISS. The impactor population also is separated in terms of nature of the projectiles, with most of the low-medium velocity ones being particles related to solid rocket motor slag condensates. On the other hand, the highest velocity projectiles are composed mainly by fragments of past in-orbit explosions.The flux of projectiles on the ISS has been calculated for the planned operative lifetime of the Station, by assuming a realistic scenario of the future debris environment evolution and the actual planned altitude profile for the ISS mission. There is a factor 2-3 variation of the flux due to the changing ISS altitude. The most dangerous part of the mission appears to be the central one, when the ISS will orbit at about 450 km above the Earth.     

Theoretical estimates of intrinsic galaxy alignment

It has recently been argued that the observed ellipticities of galaxies may be determined at least in part by the primordial tidal gravitational field in which the galaxy formed. Long-range correlations in the tidal field could thus lead to an ellipticity–ellipticity correlation for widely separated galaxies. We present a new model relating ellipticity to angular momentum, which can be calculated in linear theory. We use this model to calculate the angular power spectrum of intrinsic galaxy shape correlations. We show that, for low-redshift galaxy surveys, our model predicts that intrinsic correlations will dominate correlations induced by weak lensing, in good agreement with previous theoretical work and observations. We find that our model produces ' E -mode' correlations enhanced by a factor of 3.5 over B -modes on small scales, making it harder to disentangle intrinsic correlations from those induced by weak gravitational lensing.     

Virialization of cosmological structures in models with time-varying equation of state

We study the virialization of the cosmic structures in the framework of flat cosmological models where the dark energy component plays an important role in the global dynamics of the Universe. In particular, our analysis focuses on the study of the spherical matter perturbations, as the latter decouple from the background expansion, start to 'turn around' and finally collapse. We generalize this procedure, taking into account models with an equation of state which vary with time, and provide a complete formulation of the cluster virialization attempting to address the non-linear regime of structure formation. In particular, assuming that clusters have collapsed prior to the epoch of z _f≃ 1.4, in which the most distant cluster has been found, we show that the behaviour of the spherical collapse model depends on the functional form of the equation of state.