J. M. Ibáñez

Hydrodynamical Models of Superluminal Sources

We present numerical simulations of the generation, evolution, and radio emission of superluminal components in relativistic jets. We perform the fluid dynamical calculations using a relativistic time-dependent code based on a high-resolution shock-capturing scheme, and then we calculate the radio emission by integrating the transfer equations for synchrotron radiation. These simulations show that a temporary increase in the flow velocity at the base of the jet produces a moving perturbation that contains both a forward and a reverse shock and is trailed by a rarefaction. The perturbation appears in the simulated maps as a region of enhanced emission moving downstream at a superluminal apparent velocity. Interactions of the perturbation with the underlying steady jet result in changes in the internal brightness distribution of the superluminal component, which are manifested as low-level fluctuations about the long-term evolution of both the apparent velocity and the exponential decay of the light curves.

GENESIS: A High-Resolution Code for Three-dimensional Relativistic Hydrodynamics

The main features of a three-dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations that demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of the performance of the code on single and multiprocessor machines are given. Simulations (in double precision) with ≤7×106 computational cells require less than 1 Gbyte of RAM memory and ≈ 7×10-5 CPU s per zone and time step (on a SCI Cray-Origin 2000 with a R10000 processor). Currently, a version of the numerical code is under development, which is suited for massively parallel computers with distributed memory architecture (such as, e.g., Cray T3E).

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Radio Emission from Three-dimensional Relativistic Hydrodynamic Jets: Observational Evidence of Jet Stratification

We present the first radio emission simulations from high-resolution three-dimensional relativistic hydrodynamic jets; these simulations allow us to study the observational implications of the interaction between the jet and the external medium. This interaction gives rise to a stratification of the jet in which a fast spine is surrounded by a slow high-energy shear layer. The stratification (in particular, the large specific internal energy and slow flow in the shear layer) largely determines the emission from the jet. If the magnetic field in the shear layer becomes helical (e.g., resulting from an initial toroidal field and an aligned field component generated by shear), the emission shows a cross section asymmetry, in which either the top or the bottom of the jet dominates the emission. This, as well as limb or spine brightening, is a function of the viewing angle and flow velocity, and the top/bottom jet emission predominance can be reversed if the jet changes direction with respect to the observer or if it presents a change in velocity. The asymmetry is more prominent in the polarized flux because of field cancellation (or amplification) along the line of sight. Recent observations of jet cross section emission asymmetries in the blazar 1055+018 can be explained by assuming the existence of a shear layer with a helical magnetic field.

Stability analysis of relativistic jets from collapsars and its implications on the short-term variability of gamma-ray bursts

We consider the transverse structure and stability properties of relativistic jets formed in the course of the collapse of a massive progenitor. Our numerical simulations show the presence of a strong shear in the bulk velocity of such jets. This shear can be responsible for a very rapid shear-driven instability that arises for any velocity profile. This conclusion has been confirmed both by numerical simulations and theoretical analysis. The instability leads to rapid fluctuations of the main hydrodynamical parameters (density, pressure, Lorentz factor, etc.). However, the perturbations of the density are eectively decoupled from those of the pressure because the beam of the jet is radiation-dominated. The characteristic growth time of instability is much shorter than the life time of the jet and, therefore, may lead to a complete turbulent beam. In the course of the non-linear evolution, these fluctuations may yield to internal shocks which can be randomly distributed in the jet. In the case that internal shocks in a ultrarelativistic outflow are responsible for the observed phenomenology of gamma-ray bursts, the proposed instability can well account for the short-term variability of gamma-ray light curves down to milliseconds.We consider the transverse structure and stability properties of relativistic jets formed in the course of the collapse of a massive progenitor. Our numerical simulations show the presence of a strong shear in the bulk velocity of such jets. This shear can be responsible for a very rapid shear--driven instability that arises for any velocity profile. This conclusion has been confirmed both by numerical simulations and theoretical analysis. The instability leads to rapid fluctuations of the main hydrodynamical parameters (density, pressure, Lorentz factor, etc.). However, the perturbations of the density are effectively decoupled from those of the pressure because the beam of the jet is radiation--dominated. The characteristic growth time of instability is much shorter than the life time of the jet and, therefore, may lead to a complete turbulent beam. In the course of the non-linear evolution, these fluctuations may yield to internal shocks which can be randomly distributed in the jet. In the case that internal shocks in a ultrarelativistic outflow are responsible for the observed phenomenology of gamma-ray bursts, the proposed instability can well account for the short-term variability of gamma-ray light curves down to milliseconds.

CFC+: Improved dynamics and gravitational waveforms from relativistic core collapse simulations

Received date / Accepted date Abstract. Core collapse supernovae are a promising source of detectable gravitational waves. Most of the existing (multidimensional) numerical simulations of core collapse in general relativity have been done using approxima- tions of the Einstein field equations. As recently shown by Dimmelmeier et al. (2002a,b), one of the most inter- esting such approximation is the so-called conformal flatness condition (CFC) of Isenberg, Wilson and Mathews. Building on this previous work we present here new results from numerical simulations of relativistic rotational core collapse in axisymmetry, aiming at improving the dynamics and the gravitational waveforms. The computer code used for these simulations evolves the coupled system of metric and fluid equations using the 3+1 formalism, specialized to a new framework for the gravitational field equations which we call CFC+. In this approach we add new degrees of freedom to the original CFC equations, which extend them by terms of second post-Newtonian or- der. The resulting metric equations are still of elliptic type but the number of equations is significantly augmented in comparison to the original CFC approach. The hydrodynamics evolution and the CFC spacetime metric are calculated using the code developed by Dimmelmeier et al. (2002a), which has been conveniently extended to account for the additional CFC+ equations. The corrections for CFC+ are computed solving a system of elliptic linear equations. The new formalism is assessed with time evolutions of both rotating neutron stars in equilibrium and gravitational core collapse of rotating polytropes. Gravitational wave signals for a comprehensive sample of collapse models are extracted using either the quadrupole formula or directly from the metric. We discuss our results on the dynamics and the gravitational wave emission through a detailed comparison between CFC and CFC+ simulations. The main conclusion is that, for the neutron star spacetimes analyzed in the present work, no significant differences are found among CFC, CFC+, and full general relativity, which highlights the suitability of the former.

Relativistic simulations of superluminal sources

Abstract We present numerical simulations of the radio emission from hydrodynamical relativistic jets. The quiescent-state jet emission consists of quasi-periodic knots of high emission, associated with internal recollimation shocks. Superluminal components can be reproduced by introducing a square-wave perturbation in the injection velocity of the jet. Strong interactions of the resulting moving shock and the standing recollimations result in a “drag” and increase in emission of the latter.

A Numerical Study of Relativistic Jets

Numerical simulations of supersonic jets are able to explain the structures observed in many VLA images of radio sources. The improvements achieved in classical simulations (see Hardee, these proceedings) are in contrast with the almost complete lack of relativistic simulations the reason being that numerical difficulties arise from the highly relativistic flows typical of extragalactic jets. For our study, we have developed a two-dimensional code which is based on (i) an explicit conservative differencing of the special relativistic hydrodynamics (SRH) equations and (ii) the use of an approximate Riemann solver (see Marti et al. 1995a,b and references therein).

Anomalous dynamics triggered by a non-convex equation of state in relativistic flows

The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density (

Upwind Relativistic MHD Code for Astrophysical Applications

n_0 \approx 0.16\,