A numerical simulation of the evolution and fate of a FRI jet. The case of 3C 31

The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 % of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, in agreement with recent X-ray observations that show the presence of bow shocks in FRI sources. The jet expands until it becomes underpressured with respect to the ambient medium, and then recollimates. Subsequent oscillations around pressure equilibrium and generation of standing shocks lead to the mass loading and disruption of the jet flow. We derive an estimate for the minimum age of the source of $t>1. 10^8 \rm{yrs}$, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak CSS sources may be the young counterparts of FRIs. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario.Comment: Accepted for publication in MNRAS. For better quality figures, please check http://www.mpifr-bonn.mpg.de/staff/mperucho/Research.htm

Limit to the radio emission from a putative central compact source in SN1993J

SN1993J in M81 is the best studied young radio-luminous supernova in the Northern Hemisphere. We recently reported results from the analysis of a complete set of VLBI observations of this supernova at 1.7, 2.3, 5.0, and 8.4 GHz, covering a time baseline of more than one decade. Those reported results were focused on the kinematics of the expanding shock, the particulars of its evolving non-thermal emission, the density profile of the circumstellar medium, and the evolving free-free opacity by the supernova ejecta. In the present paper, we complete our analysis by performing a search for any possible signal from a compact source (i.e., a stellar-mass black hole or a young pulsar nebula) at the center of the expanding shell. We have performed a stacking of all our VLBI images at each frequency, after subtraction of our best-fit shell model at each epoch, and measured the peak intensity in the stacked residual image. Given the large amount of available global VLBI observations, the stacking of all the residual images allows us to put upper limits to the eventual emission of a putative compact central source at the level of $\sim102$ $\mu$Jy at 5 GHz (or, more conservatively, $\sim192$ $\mu$Jy, if we make a further correction for the ejecta opacity) and somewhat larger at other wavelengths.Comment: 4 pages, 3 figures. Accepted for publication in A&

Simulations of the relativistic parsec-scale jet in 3C273

We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets.Comment: 4 pages, 3 figures; to be published in Proceedings of the workshop "Multiband Approach to AGN", held on Sep.30-Oct.2 in Bonn. Publication: Memorie della Societa Astronomica Italiana, v. 26, No.1 (2005). Reduced figure resolution! Version with original figures is availavble at http://www.mpifr-bonn.mpg.de/bonn04/proceedings/perucho.pd

RAM: A Relativistic Adaptive Mesh Refinement Hydrodynamics Code

We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite difference, weighted essentially non-oscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth and fifth order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods and physics modules. In addition to WENO we have implemented a finite volume module with the piecewise parabolic method (PPM) for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that under-resolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two and three dimensions and in Cartesian, cylindrical and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.Comment: ApJS in press, 21 pages including 18 figures (6 color figures

Front dynamics in turbulent media

A study of a stable front propagating in a turbulent medium is presented. The front is generated through a reaction-diffusion equation, and the turbulent medium is statistically modeled using a Langevin equation. Numerical simulations indicate the presence of two different dynamical regimes. These regimes appear when the turbulent flow either wrinkles a still rather sharp propagating interfase or broadens it. Specific dependences of the propagating velocities on stirring intensities appropriate to each case are found and fitted when possible according to theoretically predicted laws. Different turbulent spectra are considered.Comment: 8 pages, REVTEX, 6 postscript figures included. To appear in Phys. Fluids (1997

The gamma-ray emitting microquasar LSI+61303

LS I +61 303 is one of the most studied X-ray binary systems because of its two peculiarities: On the one hand being the probable counterpart of the variable gamma ray source 2CG 135+01 (Gregory and Taylor 1978; Tavani et al. 1998) and on the other hand being a periodic radio source (Taylor and Gregory 1982). The recent discovery of a radio emitting jet extending ca. 200 AU at both sides of a central core (Massi et al. 2004) in all evidence has shown the occurrence of accretion/ejection processes in this system. However, the radio outbursts do not occur at periastron passage, where the accretion is at its maximum, but several days later. In addition, when the gamma-ray emission of 2CG 135+01 is examined along the orbital phase of LS I +61 303 one sees that this emission seems to peak at periastron passage (Massi 2004). Here in detail we analyse the trend of the gamma-ray data versus orbital phase and discuss the delay between peaks at gamma-rays and in the radio band within the framework of a two-peak accretion/ejection model proposed by Taylor et al. (1992) and further developed by Marti' and Paredes (1995).Comment: To be published in the proceedings of the Symposium on High-Energy Gamma-Ray Astronomy, Heidelberg, July 26-30, 2004 (AIP Proceedings Series