Showing 1–50 of 52 results for author: Meliani, Z

Hydrodynamical simulations of wind interaction in spider systems : A step toward understanding transitional millisecond pulsars

Authors: C. Guerra, Z. Meliani, G. Voisin

Abstract: The detected population of "spiders" has significantly grown in the past decade thanks to multiwavelength follow-up investigations of unidentified Fermi sources. These systems consist of low-mass stellar companions orbiting rotation-powered millisecond pulsars in short periods of a few hours up to day. Among them, a subset of intriguing objects called transitional millisecond pulsars (tMSPs) has b… ▽ More The detected population of "spiders" has significantly grown in the past decade thanks to multiwavelength follow-up investigations of unidentified Fermi sources. These systems consist of low-mass stellar companions orbiting rotation-powered millisecond pulsars in short periods of a few hours up to day. Among them, a subset of intriguing objects called transitional millisecond pulsars (tMSPs) has been shown to exhibit a remarkable behavior, transitioning between pulsar-binary and faint low-mass X-ray binary states over a span of a few years. Our objective is to study the interaction of stellar winds in tMSPs in order to understand their observational properties. To this end we focus on the parameter range that places the system near Roche-lobe overflow. Employing the adaptative mesh refinement (AMR) AMRVAC 2.0 code, we performed 2D hydrodynamical (HD) simulations of the interaction between the flows from both stars, accounting for the effects of gravity and orbital motion. By studying the mass loss and launch speed of the winds, we successfully recreated two phenomenologically distinct regimes: the accretion stream and the radio pulsar state. We also identified the tipping point that marks the sharp transition between these two states. In the pulsar state, we reconstructed the corresponding X-ray light curves of the system that produces the characteristic double-peak pattern of these systems. The position of the peaks is shifted due to orbital motion and the leading peak is weaker due to eclipsing by the companion. We suggest that a smaller leading peak in X-rays is indicative of a nearly edge-on system. This study highlights the importance of gravity and orbital motion in the interaction between the companion and pulsar winds. Our setup allows the study of the complex interaction between the pulsar wind and an accretion stream during mass transfer. △ Less

Submitted 23 July, 2024; v1 submitted 20 July, 2024; originally announced July 2024.

Comments: A&A accepted

The galactic bubbles of starburst galaxies The influence of galactic large-scale magnetic fields

Authors: Z. Meliani, P. Cristofari, A. Rodriguez-Gonzalez, G. Fichet de Clairfontaine, E. Proust, E. Peretti

Abstract: Context. The galactic winds of starburst galaxies (SBGs) give rise to remarkable structures on kiloparsec scales. However, the evolution and shape of these giant wind bubbles, as well as the properties of the shocks they develop, are not yet fully understood. Aims. We aim to understand what shapes the galactic winds of SBGs, with a particular focus on the role of large-scale magnetic fields in the… ▽ More Context. The galactic winds of starburst galaxies (SBGs) give rise to remarkable structures on kiloparsec scales. However, the evolution and shape of these giant wind bubbles, as well as the properties of the shocks they develop, are not yet fully understood. Aims. We aim to understand what shapes the galactic winds of SBGs, with a particular focus on the role of large-scale magnetic fields in the dynamical evolution of galactic wind-inflated bubbles. In addition, we aim to explore where the conditions for efficient particle acceleration are met in these systems. Methods. We performed magnetohydrodynamic simulations with the AMRVAC code (Adaptive Mesh Refinement Versatile Advection Code) with various configurations of the galactic medium density profile and magnetization. Results. We observe that the large-scale magnetic field, in which galactic winds expand, can impact the structure and evolution of inflated bubbles. However, the typical structures observed in starburst galaxies, such as M82, cannot be solely explained by the magnetic field structures that have been considered. This highlights the importance of other factors, such as the galactic disk, in shaping the galactic bubble. Furthermore, in all the magnetized cases we investigated, the forward wave resulting from the expanding bubbles only results in compression waves, whereas the wind termination shock features high Mach numbers, making it a promising site for diffusive shock acceleration up to $\sim 10^{2}$ PeV. The synthetic X-ray images generated from our models reveal an envelope surrounding the bubbles that extends up to 2 kpc, which could correspond to the polarized emission observed from planar geometry in M82, as well as a large structure inside the bubble corresponding to the shocked galactic wind. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 17 pages

Pulsar wind nebulae of runaway massive stars

Authors: D. M. -A. Meyer, Z. Meliani

Abstract: A significant fraction of massive stars move at speed through the interstellar medium of galaxies. After their death as core collapse supernovae, a possible final evolutionary state is that of a fast rotating magnetised neutron star, shaping its circumstellar medium into a pulsar wind nebula. Understanding the properties of pulsar wind nebulae requires knowledge of the evolutionary history of thei… ▽ More A significant fraction of massive stars move at speed through the interstellar medium of galaxies. After their death as core collapse supernovae, a possible final evolutionary state is that of a fast rotating magnetised neutron star, shaping its circumstellar medium into a pulsar wind nebula. Understanding the properties of pulsar wind nebulae requires knowledge of the evolutionary history of their massive progenitors. Using 2.5D magnetohydrodynamical simulations, we demonstrate that, in the context of a runaway high mass red supergiant supernova progenitor, the morphology of its subsequent pulsar wind nebula is strongly affected by the wind of the defunct progenitor star preshaping the stellar surroundings throughout its entire past life. In particular, pulsar wind nebulae of obscured runaway massive stars harbour asymmetries function of the morphology of the progenitors wind blown cavity, inducing projected asymmetric up down synchrotron emission. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: Accepted at MNRAS Letters (5 pages, 3 figures)

Wide-angle protostellar outflows driven by narrow jets in stratified cores

Authors: M. Rabenanahary, S. Cabrit, Z. Meliani, G. Pineau des Forêts

Abstract: Most simulations of outflow feedback on star formation are based on the assumption that outflows are driven by a wide angle "X-wind," rather than a narrow jet. However, the arguments initially raised against pure jet-driven flows were based on steady ejection in a uniform medium, a notion that is no longer supported based on recent observations. We aim to determine whether a pulsed narrow jet laun… ▽ More Most simulations of outflow feedback on star formation are based on the assumption that outflows are driven by a wide angle "X-wind," rather than a narrow jet. However, the arguments initially raised against pure jet-driven flows were based on steady ejection in a uniform medium, a notion that is no longer supported based on recent observations. We aim to determine whether a pulsed narrow jet launched in a density-stratified, self-gravitating core could reproduce typical molecular outflow properties, without the help of a wide-angle wind component. We performed axisymmetric hydrodynamic simulations using the MPI-AMRVAC code with optically thin radiative cooling on timescales up to 10000 yrs. Then we computed and compared the predicted properties with observational data. First, the jet-driven shell expands faster and wider through a core with steeply decreasing density than through an uniform core. Second, when blown into the same singular flattened core, a jet-driven shell has a similar width as a wide-angle wind-driven shell in the first few hundred years, but a decelerating expansion on long timescales. The flow adopts a conical shape and a base opening angle reaching up to $90\unicode{xb0}$. Third, after $\sim$ 10000 yrs, a pulsed jet-driven shell shows fitting features and a qualitative resemblance with recent observations of protostellar outflows with the Atacama Large Millimeter Array (ALMA), such as HH46-47 and CARMA-7. In particular, similarities are seen in the shell widths, opening angles, position-velocity diagrams, and mass-velocity distribution, with some showing a closer resemblance than in simulations based on a wide-angle "X-wind" model. Therefore, a realistic ambient density stratification in addition to millenia-long integration times are equally essential to reliably predict the properties of outflows driven by a pulsed jet and to confront them with the observations. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 20 pages, 13 figures, accepted for publication in A&A, Reproduced with permission from Astronomy & Astrophysics, ©ESO

Journal ref: A&A 664, A118 (2022)

Flare echos from relaxation shocks in perturbed relativistic jets

Authors: Gaëtan Fichet de Clairfontaine, Zakaria Meliani, Andreas Zech

Abstract: One of the main scenarios to account for the multi-wavelength flux variability observed in relativistic jets of AGNs is based on diffusive shock acceleration of a population of relativistic electrons on internal shocks of various origins. Any complete AGN emission scenario has to be able to explain the wide range of observed variability time scales between the radio and gamma-ray band. In addition… ▽ More One of the main scenarios to account for the multi-wavelength flux variability observed in relativistic jets of AGNs is based on diffusive shock acceleration of a population of relativistic electrons on internal shocks of various origins. Any complete AGN emission scenario has to be able to explain the wide range of observed variability time scales between the radio and gamma-ray band. In addition constraints are also provided by very-long-baseline interferometry (VLBI), which shows a large variety of moving and standing emission zones with distinct behaviors. We aim to characterize the evolution of stationary and moving emission zones in the jet and to study their multi-wavelength signatures through emission maps and light curves. We focus our study on flare events that occur during strong interactions between moving ejecta and stationary recollimation shocks. We simulate relativistic jets with the magneto-hydrodynamic code MPI-AMRVAC and inject non-thermal particle distributions of electrons in shock regions. We follow the propagation of a moving shock and its interactions with a structure of standing re-collimation shocks in the jet. Synchrotron emission and radiative transfer are calculated in the post-processing code RIPTIDE for given observation angles and frequencies. We demonstrate the appearance of trailing components behind the leading moving shock. The latter destabilizes the jet, causing the emergence of oscillating standing shocks and relaxation shocks. Emissions from these regions can dominate the overall flux or lead to flare echos in the light curve. Another observational marker of relaxation shocks appears in time-distance plots of bright VLBI components of the jet. Our scenario provides a plausible explanation for radio VLBI observations of the radio-galaxy 3C 111 where trailing components have been observed during a radio outburst event in 1997. △ Less

Submitted 5 March, 2022; originally announced March 2022.

Comments: 13 pages, 12 figures

Journal ref: A&A 661, A54 (2022)

Two-wind interactions in binaries with two orbiting giant stars

Authors: A. Castellanos-Ramírez, A. Rodríguez-González, Z. Meliani, P. R. Rivera-Ortiz, A. C. Raga, J. Cantó

Abstract: Some red giant envelopes present spiral structures (seen either in dust scattered stellar continuum or in molecular line emission), the most striking example probably being AFGL 3068. This object has been modeled (both analytically and numerically) in terms of a wind ejected from a star in orbit around a binary companion. We revisit both analytical models and 3D simulations of a wind from an orbit… ▽ More Some red giant envelopes present spiral structures (seen either in dust scattered stellar continuum or in molecular line emission), the most striking example probably being AFGL 3068. This object has been modeled (both analytically and numerically) in terms of a wind ejected from a star in orbit around a binary companion. We revisit both analytical models and 3D simulations of a wind from an orbiting red giant star, and extend the numerical simulations to the case of a binary with two red giants with strong winds. We find that most two-wind models on the orbital plane show a "double spiral" structure close to the binary source, and that these two arms merge into a single spiral structure at larger distances. However, for the case of a binary with two identical winds the two spiral arms are still present at large distances from the binary source. We also find that for models of two (not identical) dynamically important winds, a region close to the orbital plane has material from both winds. Also, an approximately conical region centered on the orbital axis is filled exclusively by the wind with larger momentum rate. These two structures lead to morphologies reminiscent of the so-called "hour glass" planetary nebulae. Finally, we find that increasing wind velocity disparities lead to the formation of clumpy structures along the spiral amrs. Observations of "clumpy spirals" are therefore likely to indicate the presence of two strong winds from the stars in the central binary system. △ Less

Submitted 14 August, 2021; originally announced August 2021.

Comments: 9 pages, 10 figures

Flux variability from ejectas in structured relativistic jets with large-scale magnetic fields

Authors: Gaëtan Fichet de Clairfontaine, Zakaria Meliani, Andreas Zech, Olivier Hervet

Abstract: Using the relativistic MHD code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks an… ▽ More Using the relativistic MHD code MPI-AMRVAC and a radiative transfer code in post-processing, we explore the influence of the magnetic-field configuration and transverse stratification of an over-pressured jet on its morphology, on the moving shock dynamics, and on the emitted radio light curve. First, we investigate different large-scale magnetic fields with their effects on the standing shocks and on the stratified jet morphology. Secondly, we study the interaction of a moving shock wave with the standing shocks. We calculate the synthetic synchrotron maps and radio light curves and analyse the variability at two frequencies 1 and 15.3 GHz and for several observation angles. Finally, we compare the characteristics of our simulated light curves with radio flares observed from the blazar 3C 273 with OVRO and VLBA in the MOJAVE survey between 2008 and 2019. We find that, in a structured, over-pressured relativistic jet, the presence of the large-scale magnetic field structure changes the properties of the standing shock waves and leads to an opening of the jet. When crossing such standing shocks, moving shock waves accompanying overdensities injected in the base of the jet are causing very luminous radio flares. The observation of the temporal structure of these flares under different viewing angles probes the jet at different optical depths. At 1 GHz and for small angles, the self-absorption caused by the moving shock wave becomes more important and leads to a drop in the observed flux after it interacts with the brightest standing knot. A weak asymmetry is seen in the shape of the simulated flares, resulting from the remnant emission of the shocked standing shocks. The characteristics of the simulated flares and the correlation of peaks in the light curve with the crossing of moving and standing shocks favor this scenario as an explanation of the observed radio flares of 3C 273. △ Less

Submitted 25 January, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: 19 pages, 14 figures

Journal ref: A&A 647, A77 (2021)

Multi-scale simulations of particle acceleration in astrophysical systems

Authors: A. Marcowith, G. Ferrand, M. Grech, Z. Meliani, I. Plotnikov, R. Walder

Abstract: This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a n… ▽ More This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a non-thermal distribution in both Newtonian and relativistic astrophysical flows, namely the first and second order Fermi acceleration processes. We also discuss shock drift and surfing acceleration, two processes important in the context of particle injection in shock acceleration. We analyze with some details the particle-in-cell (PIC) approach used to describe particle kinetics. We review the main results obtained with PIC simulations in the recent years concerning particle acceleration at shocks and in reconnection events. The review discusses the solution of Fokker-Planck problems with application to the study of particle acceleration at shocks but also in hot coronal plasmas surrounding compact objects. We continue by considering large scale physics. We describe recent developments in magnetohydrodynamic (MHD) simulations. We give a special emphasize on the way energetic particle dynamics can be coupled to MHD solutions either using a multi-fluid calculation or directly coupling kinetic and fluid calculations. This aspect is mandatory to investigate the acceleration of particles in the deep relativistic regimes to explain the highest Cosmic Ray energies. △ Less

Submitted 25 February, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: 202 pages, 43 figures, to be published in the Living Reviews of Computational Astrophysics

Report number: RIKEN-iTHEMS-Report-20

Numerical models for the dust in RCW 120

Authors: A. Rodríguez-González, Z. Meliani, M. Sánchez-Cruces, P. R. Rivera-Ortiz, A. Castellanos-Ramírez

Abstract: The interstellar bubble RCW 120 seen around a type O runaway star is driven by the stellar wind and the ionising radiation emitted by the star. The boundary between the stellar wind and interstellar medium (ISM) is associated with the arc-shaped mid-infrared dust emission around the star within the HII region. The interstellar bubble RCW 120 seen around a type O runaway star is driven by the stellar wind and the ionising radiation emitted by the star. The boundary between the stellar wind and interstellar medium (ISM) is associated with the arc-shaped mid-infrared dust emission around the star within the HII region. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 6 pages, 5 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 631, A170 (2019)

Comparing timelike geodesics around a Kerr black hole and a boson star

Authors: M. Grould, Z. Meliani, F. H. Vincent, P. Grandclément, E. Gourgoulhon

Abstract: The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $μ$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are su… ▽ More The second-generation beam combiner at the Very Large Telescope (VLT), GRAVITY, observes the stars orbiting the compact object located at the center of our galaxy, with an unprecedented astrometric accuracy of 10 $μ$as. The nature of this compact source is still unknown since black holes are not the only candidates explaining the four million solar masses at the Galactic center. Boson stars are such an alternative model to black holes. This paper focuses on the study of trajectories of stars orbiting a boson star and a Kerr black hole. We put in light strong differences between orbits obtained in both metrics when considering stars with sufficiently close pericenters to the compact object, typically $\lesssim 30~M$. Discovery of closer stars to the Galactic center than the S2 star by the GRAVITY instrument would thus be a powerful tool to possibly constrain the nature of the central source. △ Less

Submitted 18 September, 2017; originally announced September 2017.

Comments: 30 pages, 20 figures

Rotation and toroidal magnetic field effects on the stability of two-component jets

Authors: Dimitrios Millas, Rony Keppens, Zakaria Meliani

Abstract: Several observations of astrophysical jets show evidence of a structure in the direction perpendicular to the jet axis, leading to the development of "spine & sheath" models of jets. Most studies focus on a two-component jet consisting of a highly relativistic inner jet and a slower - but still relativistic - outer jet surrounded by an unmagnetized environment. These jets are believed to be suscep… ▽ More Several observations of astrophysical jets show evidence of a structure in the direction perpendicular to the jet axis, leading to the development of "spine & sheath" models of jets. Most studies focus on a two-component jet consisting of a highly relativistic inner jet and a slower - but still relativistic - outer jet surrounded by an unmagnetized environment. These jets are believed to be susceptible to a relativistic Rayleigh-Taylor-type instability, depending on the effective inertia ratio of the two components. We extend previous studies by taking into account the presence of a non-zero toroidal magnetic field. Different values of magnetization are examined, to detect possible differences in the evolution and stability of the jet. We find that the toroidal field, above a certain level of magnetization $σ$, roughly equal to 0.01, can stabilize the jet against the previously mentioned instabilities and that there is a clear trend in the behaviour of the average Lorentz factor and the effective radius of the jet when we continuously increase the magnetization. The simulations are performed using the relativistic MHD module from the open source, parallel, grid adaptive, MPI-AMRVAC code. △ Less

Submitted 21 June, 2017; originally announced June 2017.

Comments: 15 pages, 17 figures

Shocks in relativistic transverse stratified jets, a new paradigm for radio-loud AGN

Authors: O. Hervet, Z. Meliani, A. Zech, C. Boisson, V. Cayatte, C. Sauty, H. Sol

Abstract: The transverse stratification of active galactic nuclei (AGN) jets is suggested by observations and theoretical arguments, as a consequence of intrinsic properties of the central engine (accretion disc + black hole) and external medium. On the other hand, the one-component jet approaches are heavily challenged by the various observed properties of plasmoids in radio jets (knots), often associated… ▽ More The transverse stratification of active galactic nuclei (AGN) jets is suggested by observations and theoretical arguments, as a consequence of intrinsic properties of the central engine (accretion disc + black hole) and external medium. On the other hand, the one-component jet approaches are heavily challenged by the various observed properties of plasmoids in radio jets (knots), often associated with internal shocks. Given that such a transverse stratification plays an important role on the jets acceleration, stability, and interaction with the external medium, it should also induce internal shocks with various strengths and configurations, able to describe the observed knots behaviours. By establishing a relation between the transverse stratification of the jets, the internal shock properties, and the multiple observed AGN jet morphologies and behaviours, our aim is to provide a consistent global scheme of the various AGN jet structures. Working on a large sample of AGN radio jets monitored in very long baseline interferometry (VLBI) by the MOJAVE collaboration, we determined the consistency of a systematic association of the multiple knots with successive re-collimation shocks. We then investigated the re-collimation shock formation and the influence of different transverse stratified structures by parametrically exploring the two relativistic outflow components with the specific relativistic hydrodynamic (SRHD) code AMRVAC. We were able to link the different spectral classes of AGN with specific stratified jet characteristics, in good accordance with their VLBI radio properties and their accretion regimes. △ Less

Submitted 30 May, 2017; originally announced May 2017.

Comments: 16 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 606, A103 (2017)

Impact of Schwarzschild black hole's gravity upon the Rossby wave instability

Authors: F. Casse, P. Varniere, Z. Meliani

Abstract: In an early work the Rossby Wave Instability (RWI) has been proposed to explain variability thought to originate in the close vicinity of black-holes but this was done in the pseudo-Newtonian approach. Here we present the first general relativistic hydrodynamics simulations of this instability not only proving its theorized existence in a full general relativistic (GR) environment but also studyin… ▽ More In an early work the Rossby Wave Instability (RWI) has been proposed to explain variability thought to originate in the close vicinity of black-holes but this was done in the pseudo-Newtonian approach. Here we present the first general relativistic hydrodynamics simulations of this instability not only proving its theorized existence in a full general relativistic (GR) environment but also studying the effect of the strong gravity on the instability. To that end we performed a set of simulations increasingly closer to the black hole with our new GR version of the MPI-AMRVAC code.This allows us to study the minute changes in the behaviour of the instability.We found that the pseudo-Newtonian approach gives adequate results provided that time shifting induced by the black hole gravity is taken into account. Hence, to view the disc as a distant observer would a full GR ray-tracing post treatment of the simulations is a must. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 10 pages, 11 figures, MNRAS accepted

Simulations of recoiling black holes: adaptive mesh refinement and radiative transfer

Authors: Zakaria Meliani, Yosuke Mizuno, Hector Olivares, Oliver Porth, Luciano Rezzolla, Ziri Younsi

Abstract: (Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accur… ▽ More (Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of AMR techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to compute accurately the electromagnetic emissions from such accretion flows. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry and performed either in 2D or 3D. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black hole binary interacting with the surrounding circumbinary disc. In this way we can present, for the first time, ray-traced images of the shocked fluid and the light-curve resulting from consistent general-relativistic radiation-transport calculations from this process. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to deal accurately and self-consistently with accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are presently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars. △ Less

Submitted 10 November, 2016; v1 submitted 27 June, 2016; originally announced June 2016.

Comments: 20 pages, 20 figures, accepted for publication in A&A

Journal ref: A&A 598, A38 (2017)

Pinwheels in the sky, with dust: 3D modeling of the Wolf-Rayet 98a environment

Authors: Tom Hendrix, Rony Keppens, Allard Jan van Marle, Peter Camps, Maarten Baes, Zakaria Meliani

Abstract: The Wolf-Rayet 98a (WR 98a) system is a prime target for interferometric surveys, since its identification as a "rotating pinwheel nebulae", where infrared images display a spiral dust lane revolving with a 1.4 year periodicity. WR 98a hosts a WC9+OB star, and the presence of dust is puzzling given the extreme luminosities of Wolf-Rayet stars. We present 3D hydrodynamic models for WR 98a, where du… ▽ More The Wolf-Rayet 98a (WR 98a) system is a prime target for interferometric surveys, since its identification as a "rotating pinwheel nebulae", where infrared images display a spiral dust lane revolving with a 1.4 year periodicity. WR 98a hosts a WC9+OB star, and the presence of dust is puzzling given the extreme luminosities of Wolf-Rayet stars. We present 3D hydrodynamic models for WR 98a, where dust creation and redistribution are self-consistently incorporated. Our grid-adaptive simulations resolve details in the wind collision region at scales below one percent of the orbital separation (~4 AU), while simulating up to 1300 AU. We cover several orbital periods under conditions where the gas component alone behaves adiabatic, or is subject to effective radiative cooling. In the adiabatic case, mixing between stellar winds is effective in a well-defined spiral pattern, where optimal conditions for dust creation are met. When radiative cooling is incorporated, the interaction gets dominated by thermal instabilities along the wind collision region, and dust concentrates in clumps and filaments in a volume-filling fashion, so WR 98a must obey close to adiabatic evolutions to demonstrate the rotating pinwheel structure. We mimic Keck, ALMA or future E-ELT observations and confront photometric long-term monitoring. We predict an asymmetry in the dust distribution between leading and trailing edge of the spiral, show that ALMA and E-ELT would be able to detect fine-structure in the spiral indicative of Kelvin-Helmholtz development, and confirm the variation in photometry due to the orientation. Historic Keck images are reproduced, but their resolution is insufficient to detect the details we predict. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: Accepted for publication in mnras

Circular geodesics and thick tori around rotating boson stars

Authors: Z. Meliani, F. H. Vincent, P. Grandclément, E. Gourgoulhon, R. Monceau-Baroux, O. Straub

Abstract: Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of i… ▽ More Accretion disks play an important role in the evolution of their relativistic inner compact objects. The emergence of a new generation of interferometers will allow to resolve these accretion disks and provide more information about the properties of the central gravitating object. Due to this instrumental leap forward it is crucial to investigate the accretion disk physics near various types of inner compact objects now to deduce later constraints on the central objects from observations. A possible candidate for the inner object is the boson star. Here, we will try to analyze the differences between accretion structures surrounding boson stars and black holes. We aim at analysing the physics of circular geodesics around boson stars and study simple thick accretion tori (so-called Polish doughnuts) in the vicinity of these stars. We realize a detailed study of the properties of circular geodesics around boson stars. We then perform a parameter study of thick tori with constant angular momentum surrounding boson stars. This is done using the boson star models computed by a code constructed with the spectral solver library KADATH. We demonstrate that all the circular stable orbits are bound. In the case of a constant angular momentum torus, a cusp in the torus surface exists only for boson stars with a strong gravitational scalar field. Moreover, for each inner radius of the disk, the allowed specific angular momentum values lie within a constrained range which depends on the boson star considered. We show that the accretion tori around boson stars have different characteristics than in the vicinity of a black hole. With future instruments it could be possible to use these differences to constrain the nature of compact objects. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Comments: Accepted for publication in CQG

Imaging a boson star at the Galactic center

Authors: F. H. Vincent, Z. Meliani, P. Grandclement, E. Gourgoulhon, O. Straub

Abstract: Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an… ▽ More Millimeter very long baseline interferometry will soon produce accurate images of the closest surroundings of the supermassive compact object at the center of the Galaxy, Sgr A*. These images may reveal the existence of a central faint region, the so-called shadow, which is often interpreted as the observable consequence of the event horizon of a black hole. In this paper, we compute images of an accretion torus around Sgr A* assuming this compact object is a boson star, i.e. an alternative to black holes within general relativity, with no event horizon and no hard surface. We show that very relativistic rotating boson stars produce images extremely similar to Kerr black holes, showing in particular shadow-like and photon-ring-like structures. This result highlights the extreme difficulty of unambiguously telling the existence of an event horizon from strong-field images. △ Less

Submitted 16 March, 2016; v1 submitted 14 October, 2015; originally announced October 2015.

Comments: 21 pages, 9 figures, accepted in CQG; main difference wrt previous version is the last paragraph of the conclusion

Shape and evolution of wind-blown bubbles of massive stars: on the effect of the interstellar magnetic field

Authors: Allard Jan van Marle, Zakaria Meliani, Alexandre Marcowith

Abstract: The winds of massive stars create large (>10 pc) bubbles around their progenitors. As these bubbles expand they encounter the interstellar coherent magnetic field which, depending on its strength, can influence the shape of the bubble. We wish to investigate if, and how much, the interstellar magnetic field can contribute to the shape of an expanding circumstellar bubble around a massive star. We… ▽ More The winds of massive stars create large (>10 pc) bubbles around their progenitors. As these bubbles expand they encounter the interstellar coherent magnetic field which, depending on its strength, can influence the shape of the bubble. We wish to investigate if, and how much, the interstellar magnetic field can contribute to the shape of an expanding circumstellar bubble around a massive star. We use the MPI-AMRVAC code to make magneto-hydrodynamical simulations of bubbles, using a single star model, combined with several different field strengths: B = 5, 10, and 20 muG for the interstellar magnetic field. This covers the typical field strengths of the interstellar magnetic fields found in the galactic disk and bulge. Furthermore, we present two simulations that include both a 5 muG interstellar magnetic field and a 10,000 K interstellar medium and two different ISM densities to demonstrate how the magnetic field can combine with other external factors to influence the morphology of the circumstellar bubbles. Our results show that low magnetic fields, as found in the galactic disk, inhibit the growth of the circumstellar bubbles in the direction perpendicular to the field. As a result, the bubbles become ovoid, rather than spherical. Strong interstellar fields, such as observed for the galactic bulge, can completely stop the expansion of the bubble in the direction perpendicular to the field, leading to the formation of a tube-like bubble. When combined with a warm, high-density ISM the bubble is greatly reduced in size, causing a dramatic change in the evolution of temporary features inside the bubble. The magnetic field of the interstellar medium can affect the shape of circumstellar bubbles. This effect may have consequences for the shape and evolution of circumstellar nebulae and supernova remnants, which are formed within the main wind-blown bubble. △ Less

Submitted 17 September, 2015; v1 submitted 1 September, 2015; originally announced September 2015.

Comments: Proposed for acceptance for publication in Astronomy & Astrophysics. The published version will contain animations of each simulation

Journal ref: A&A 584, A49 (2015)

Using numerical models of bow shocks to investigate the circumstellar medium of massive stars

Authors: Allard Jan van Marle, Leen Decin, Nick Cox, Zakaria Meliani

Abstract: Many massive stars travel through the interstellar medium at supersonic speeds. As a result they form bow shocks at the interface between the stellar wind. We use numerical hydrodynamics to reproduce such bow shocks numerically, creating models that can be compared to observations. In this paper we discuss the influence of two physical phenomena, interstellar magnetic fields and the presence of in… ▽ More Many massive stars travel through the interstellar medium at supersonic speeds. As a result they form bow shocks at the interface between the stellar wind. We use numerical hydrodynamics to reproduce such bow shocks numerically, creating models that can be compared to observations. In this paper we discuss the influence of two physical phenomena, interstellar magnetic fields and the presence of interstellar dust grains on the observable shape of the bow shocks of massive stars. We find that the interstellar magnetic field, though too weak to restrict the general shape of the bow shock, reduces the size of the instabilities that would otherwise be observed in the bow shock of a red supergiant. The interstellar dust grains, due to their inertia can penetrate deep into the bow shock structure of a main sequence O-supergiant, crossing over from the ISM into the stellar wind. Therefore, the dust distribution may not always reflect the morphology of the gas. This is an important consideration for infrared observations, which are dominated by dust emission. Our models clearly show, that the bow shocks of massive stars are useful diagnostic tools that can used to investigate the properties of both the stellar wind as well as the interstellar medium. △ Less

Submitted 7 July, 2014; originally announced July 2014.

Comments: 7 pages, 4 figures, to be published in the Journal of Physics: Conference Series (JPCS) as part of the proceedings of the 13th Annual International Astrophysics Conference

Can the magnetic field in the Orion arm inhibit the growth of instabilities in the bow shock of Betelgeuse?

Authors: Allard Jan van Marle, Leen Decin, Zakaria Meliani

Abstract: Many evolved stars travel through space at supersonic velocities, which leads to the formation of bow shocks ahead of the star where the stellar wind collides with the interstellar medium (ISM). Herschel observations of the bow shock of $α$-Orionis show that the shock is almost free of instabilities, despite being, at least in theory, subject to both Kelvin-Helmholtz and Rayleigh-Taylor instabilit… ▽ More Many evolved stars travel through space at supersonic velocities, which leads to the formation of bow shocks ahead of the star where the stellar wind collides with the interstellar medium (ISM). Herschel observations of the bow shock of $α$-Orionis show that the shock is almost free of instabilities, despite being, at least in theory, subject to both Kelvin-Helmholtz and Rayleigh-Taylor instabilities. A possible explanation for the lack of instabilities lies in the presence of an interstellar magnetic field. We wish to investigate whether the magnetic field of the interstellar medium (ISM) in the Orion arm can inhibit the growth of instabilities in the bow shock of $α$-Orionis. We used the code MPI-AMRVAC to make magneto-hydrodynamic simulations of a circumstellar bow shock, using the wind parameters derived for $α$-Orionis and interstellar magnetic field strengths of $B\,=\,1.4,\, 3.0$, and $5.0\, μ$G, which fall within the boundaries of the observed magnetic field strength in the Orion arm of the Milky Way. Our results show that even a relatively weak magnetic field in the interstellar medium can suppress the growth of Rayleigh-Taylor and Kelvin-Helmholtz instabilities, which occur along the contact discontinuity between the shocked wind and the shocked ISM. The presence of even a weak magnetic field in the ISM effectively inhibits the growth of instabilities in the bow shock. This may explain the absence of such instabilities in the Herschel observations of $α$-Orionis. △ Less

Submitted 20 December, 2013; originally announced December 2013.

Comments: 5 pages, including 7 figures. The published version will include 4 animations. Accepted for publication in A&A

Relativistic 3D jet simulations for the X-ray binary SS433

Authors: Remi Monceau-Baroux, Rony Keppens, Zakaria Meliani, Oliver Porth

Abstract: Context. Modern high resolution observations allow to view closer into the objects powering relativistic jets. This is especially the case for SS433, an X-ray binary from which a precessing jet is observed down to the sub-parsec scale. Aims. We want to study full 3D dynamics of relativistic jets associated with AGN or XRB. We study the precessing motion of a jet as a model for the jet associated… ▽ More Context. Modern high resolution observations allow to view closer into the objects powering relativistic jets. This is especially the case for SS433, an X-ray binary from which a precessing jet is observed down to the sub-parsec scale. Aims. We want to study full 3D dynamics of relativistic jets associated with AGN or XRB. We study the precessing motion of a jet as a model for the jet associated with the XRB SS433. Our study of the jet dynamics in this system focuses on the sub-parsec scales. We investigate the impact of jet precession and the variation of the Lorentz factor of the injected matter on the general 3D jet dynamics and its energy transfer to the surrounding medium. We realize synthetic radio mapping of the data, to compare our results with observations. Methods. For our study we use the code MPI-AMRVAC with SRHD model of a baryonic jet. We use a AMR scheme and an inner time-dependent boundary prescription to inject the jets. Parameters extracted from observations were used. 3D jet realizations that match the SS433 jet are intercompared. We track the energy content, as deposited in different regions of the domain affected by the jet. Our code also follows a population of particles injected with the jet. This evolving energy spectrum of accelerated electrons, allows to obtain the radio emission from our simulation. Results. we obtain meaningful observations. We find increased energy transfer for a precessing jet compared to a standing jets. We obtain synthetic radio maps for all jets. Conclusions. The synthetic radio map matches best for a model using the canonical kinematic model. Overdense precessing jets experience significant deceleration in their propagation through the ISM, and while the overall jet is of helical shape. This argument show that the kinematic model for SS433 assuming ballistic propagation has to be corrected for this decelretaion. △ Less

Submitted 2 December, 2013; v1 submitted 29 November, 2013; originally announced November 2013.

Comments: Accepted for publication in Astrophysic and Astronomy

Relativistic AGN jets I. The delicate interplay between jet structure, cocoon morphology and jet-head propagation

Authors: S. Walg, A. Achterberg, S. Markoff, R. Keppens, Z. Meliani

Abstract: Current observations have shown that astrophysical jets reveal strong signs of radial structure. They suggest that the inner region of the jet, the jet spine, consists of a low-density, fast-moving gas, while the outer region of the jet consists of a more dense and slower moving gas, called the jet sheath. Moreover, if jets carry angular momentum, the resultant centrifugal forces lead to a radial… ▽ More Current observations have shown that astrophysical jets reveal strong signs of radial structure. They suggest that the inner region of the jet, the jet spine, consists of a low-density, fast-moving gas, while the outer region of the jet consists of a more dense and slower moving gas, called the jet sheath. Moreover, if jets carry angular momentum, the resultant centrifugal forces lead to a radial stratification. Current observations are not able to fully resolve the radial structure, so little is known about its actual profile. We present three AGN jet models in $2.5D$ of which two have been given a radial structure. The first model is a homogeneous jet, the only model that doesn't carry angular momentum; the second model is a spine-sheath jet with an isothermal equation of state; and the third jet model is a (piecewise) isochoric spine-sheath jet, with constant but different densities for jet spine and jet sheath. In this paper, we look at the effects of radial stratification on jet integrity, mixing between the different jet components and global morphology of the jet-head and surrounding cocoon. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: 27 pages, 9 figures

Effect of angular opening on the dynamics of relativistic hydro jets

Authors: Remi Monceau-Baroux, Rony Keppens, Zakaria Meliani

Abstract: Context. Relativistic jets emerging from AGN cores transfer energy from the core to their surrounding ISM/IGM. Because jets are observed to have finite opening angles, one needs to quantify the role of conical versus cylindrical jet propagation in this energy transfer. Aims. We use FR-II AGN jets parameter with finite opening angles. We study the effect of the variation of the opening angle on the… ▽ More Context. Relativistic jets emerging from AGN cores transfer energy from the core to their surrounding ISM/IGM. Because jets are observed to have finite opening angles, one needs to quantify the role of conical versus cylindrical jet propagation in this energy transfer. Aims. We use FR-II AGN jets parameter with finite opening angles. We study the effect of the variation of the opening angle on the dynamics and energy transfer of the jet. We also point out how the characteristics of this external medium, such as its density profile, play a role in the dynamics. Methods. This study exploits our parallel AMR code MPI-AMRVAC with its special relativistic hydrodynamic model, incorporating an equation of state with varying effective polytropic index. We studied mildly under-dense jets up to opening angles of 10 degrees, at Lorentz factors of about 10, inspired by observations. Instantaneous quantification of the various ISM volumes and their energy content allows one to quantify the role of mixing versus shock-heated cocoon regions over the time intervals. Results. We show that a wider opening angle jet results in a faster deceleration of the jet and leads to a wider cocoon dominated by Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The energy transfer mainly occurs in the shocked ISM region by both the frontal bow shock and cocoon-traversing shock waves, in a roughly 3 to 1 ratio to the energy transfer of the mixing zone, for a 5 degree opening angle jet. A rarefaction wave induces a dynamically formed layered structure of the jet beam. Conclusions. Finite opening angle jets can efficiently transfer significant fractions (25 % up to 70 %) of their injected energy over a growing region of shocked ISM matter. The role of the ISM stratification is prominent for determining the overall volume that is affected by relativistic jet injection. △ Less

Submitted 7 November, 2012; originally announced November 2012.

Journal ref: A&A 545 A62 (2012)

Dust-trapping Rossby vortices in protoplanetary disks

Authors: H. Meheut, Z. Meliani, P. Varniere, W. Benz

Abstract: One of the most challenging steps in planet formation theory is the one leading to the formation of planetesimals of kilometre size. A promising scenario involves the existence of vortices able to concentrate a large amount of dust and grains in their centres. Up to now this scenario has been studied mostly in 2D razor thin disks. A 3D study including, simultaneously, the formation and resulting d… ▽ More One of the most challenging steps in planet formation theory is the one leading to the formation of planetesimals of kilometre size. A promising scenario involves the existence of vortices able to concentrate a large amount of dust and grains in their centres. Up to now this scenario has been studied mostly in 2D razor thin disks. A 3D study including, simultaneously, the formation and resulting dust concentration of the vortices with vertical settling, was still missing. The Rossby wave instability self-consistently forms 3D vortices, which have the unique quality of presenting a large scale vertical velocity in their centre. Here we aim to study how this newly discovered effect can alter the dynamic evolution of the dust. We perform global 3D simulations of the RWI in a radially and vertically stratified disk using the code MPI-AMRVAC. After the growth phase of the instability, the gas and solid phases are modelled by a bi-fluid approach, where the dust is considered as a fluid without pressure. Both the drag force of the gas on the dust and the back-reaction of the dust on the gas are included. Multiple grain sizes from 1mm to 5cm are used with a constant density distribution. We obtain in a short timescale a high concentration of the largest grains in the vortices. Indeed, in 3 rotations the dust-to-gas density ratio grows from 10^-2 to unity leading to a concentration of mass up to that of Mars in one vortex. The presence of the radial drift is also at the origin of a dust pile-up at the radius of the vortices. Lastly, the vertical velocity of the gas in the vortex causes the sedimentation process to be reversed, the mm size dust is lifted and higher concentrations are obtained in the upper layer than in the mid-plane. △ Less

Submitted 24 August, 2012; originally announced August 2012.

Comments: Accepted for publication in Astronomy and Astrophysics

Practical flux-prescriptions for gamma-ray burst afterglows, from early to late times

Authors: Konstantinos Leventis, Hendrik J. van Eerten, Zakaria Meliani, Ralph A. M. J. Wijers

Abstract: We present analytic flux prescriptions for broadband spectra of self-absorbed and optically thin synchrotron radiation from gamma-ray burst afterglows, based on one-dimensional relativistic hydrodynamic simulations. By treating the evolution of critical spectrum parameters as a power-law break between the ultrarelativistic and non-relativistic asymptotic solutions, we generalize the prescriptions… ▽ More We present analytic flux prescriptions for broadband spectra of self-absorbed and optically thin synchrotron radiation from gamma-ray burst afterglows, based on one-dimensional relativistic hydrodynamic simulations. By treating the evolution of critical spectrum parameters as a power-law break between the ultrarelativistic and non-relativistic asymptotic solutions, we generalize the prescriptions to any observer time. Our aim is to provide a set of formulas that constitutes a useful tool for accurate fitting of model-parameters to observational data, regardless of the dynamical phase of the outflow. The applicability range is not confined to gamma-ray burst afterglows, but includes all spherical outflows (also jets before the jet-break) that produce synchrotron radiation as they adiabatically decelerate in a cold, power-law medium. We test the accuracy of the prescriptions and show that numerical evidence suggests that typical relative errors in the derivation of physical quantities are about 10 per cent. A software implementation of the presented flux prescriptions combined with a fitting code is freely available on request and on-line. Together they can be used in order to directly fit model parameters to data. △ Less

Submitted 12 November, 2012; v1 submitted 13 June, 2012; originally announced June 2012.

Comments: 17 pages, 10 figures, submitted to MNRAS

A hydrodynamical model of the circumstellar bubble created by two massive stars

Authors: Allard Jan van Marle, Zakaria Meliani, Alexandre Marcowith

Abstract: Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by several stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass… ▽ More Numerical models of the wind-blown bubble of massive stars usually only account for the wind of a single star. However, since massive stars are usually formed in clusters, it would be more realistic to follow the evolution of a bubble created by several stars. We develope a two-dimensional (2D) model of the circumstellar bubble created by two massive stars, a 40 solar mass star and a 25 solar mass star, and follow its evolution. The stars are separated by approximately 16 pc and surrounded by a cold medium with a density of 20 particles per cubic cm. We use the MPI-AMRVAC hydrodynamics code to solve the conservation equations of hydrodynamics on a 2D cylindrical grid using time-dependent models for the wind parameters of the two stars. At the end of the stellar evolution (4.5 and 7.0 million years for the 40 and 25 solar mass stars, respectively), we simulate the supernova explosion of each star. Each star initially creates its own bubble. However, as the bubbles expand they merge, creating a combined, aspherical bubble. The combined bubble evolves over time, influenced by the stellar winds and supernova explosions. The evolution of a wind-blown bubble created by two stars deviates from that of the bubbles around single stars. In particular, once one of the stars has exploded, the bubble is too large for the wind of the remaining star to maintain and the outer shell starts to disintegrate. The lack of thermal pressure inside the bubble also changes the behavior of circumstellar features close to the remaining star. The supernovae are contained inside the bubble, which reflects part of the energy back into the circumstellar medium. △ Less

Submitted 18 April, 2012; v1 submitted 10 April, 2012; originally announced April 2012.

Comments: Accepted for publication in A&A. Six .avi files to be published online (uploaded to ArXiv DC and available as ancillary files) (updated after language corrections)

A far-infrared survey of bow shocks and detached shells around AGB stars and red supergiants

Authors: N. L. J. Cox, F. Kerschbaum, A. -J. van Marle, L. Decin, D. Ladjal, A. Mayer, M. A. T. Groenewegen, S. van Eck, P. Royer, R. Ottensamer, T. Ueta, A. Jorissen, M. Mecina, Z. Meliani, A. Luntzer, J. A. D. L. Blommaert, Th. Posch, B. Vandenbussche, C. Waelkens

Abstract: Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. "fermat… ▽ More Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. "fermata", "eyes", "irregular", and "rings") are identified and basic parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. De-projected stand-off distances (R_0) -- defined as the distance between the central star and the nearest point of the interaction region -- of the detected bow shocks ("fermata" and "eyes") are derived from the PACS images and compared to previous results, model predictions and the simulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used together to independently derive the local ISM density. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar properties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection of a bow shock -- for the distance limited sample -- appears to be governed by its physical size as set roughly by the stand-off distance. In particular the star's peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or detached rings. Tentatively, the "eyes" class objects are associated to (visual) binaries, while the "rings" generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse. △ Less

Submitted 25 October, 2011; originally announced October 2011.

Comments: 29 pages, 15 figures. Accepted for publication in A&A

Dust distribution in circumstellar shells

Authors: Allard Jan van Marle, Zakaria Meliani, Rony Keppens, Leen Decin

Abstract: We present numerical simulations of the hydrodynamical interactions that produce circumstellar shells. These simulations include several scenarios, such as wind-wind interaction and wind-ISM collisions. In our calculations we have taken into account the presence of dust in the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to the gas, large dust grains are… ▽ More We present numerical simulations of the hydrodynamical interactions that produce circumstellar shells. These simulations include several scenarios, such as wind-wind interaction and wind-ISM collisions. In our calculations we have taken into account the presence of dust in the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to the gas, large dust grains are only weakly coupled. As a result, the distribution of the large dust grains is not representative of the gas distribution. Combining these results with observations may give us a new way of validating hydrodynamical models of the circumstellar medium. △ Less

Submitted 14 October, 2011; originally announced October 2011.

Comments: 2 pagers, 2 figures, to appear in the proceedings of: "Planetary nebulae, an Eye to the Future", Tenerife, Spain, July 2011

Nonradial and nonpolytropic astrophysical outflows IX. Modeling T Tauri jets with a low mass-accretion rate

Authors: C. Sauty, Z. Meliani, J. J. G. Lima, K. Tsinganos, V. Cayatte, N. Globus

Abstract: Context: A large sample of T Tauri stars exhibits optical jets, approximately half of which rotate slowly, only at ten per cent of their breakup velocity. The disk-locking mechanism has been shown to be inefficient to explain this observational fact. Aims: We show that low mass accreting T Tauri stars may have a strong stellar jet component that can effectively brake the star to the observed rot… ▽ More Context: A large sample of T Tauri stars exhibits optical jets, approximately half of which rotate slowly, only at ten per cent of their breakup velocity. The disk-locking mechanism has been shown to be inefficient to explain this observational fact. Aims: We show that low mass accreting T Tauri stars may have a strong stellar jet component that can effectively brake the star to the observed rotation speed. Methods: By means of a nonlinear separation of the variables in the full set of the MHD equations we construct semi- analytical solutions describing the dynamics and topology of the stellar component of the jet that emerges from the corona of the star. Results: We analyze two typical solutions with the same mass loss rate but different magnetic lever arms and jet radii. The first solution with a long lever arm and a wide jet radius effectively brakes the star and can be applied to the visible jets of T Tauri stars, such as RY Tau. The second solution with a shorter lever arm and a very narrow jet radius may explain why similar stars, either Weak line T Tauri Stars (WTTS) or Classical T Tauri Stars (CTTS) do not all have visible jets. For instance, RY Tau itself seems to have different phases that probably depend on the activity of the star. Conclusions: First, stellar jets seem to be able to brake pre-main sequence stars with a low mass accreting rate. Second, jets may be visible only part time owing to changes in their boundary conditions. We also suggest a possible scenario for explaining the dichotomy between CTTS and WTTS, which rotate faster and do not have visible jets. △ Less

Submitted 13 July, 2011; originally announced July 2011.

Synchrotron radiation of self-collimating relativistic MHD jets

Authors: Oliver Porth, Christian Fendt, Zakaria Meliani, Bhargav Vaidya

Abstract: The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launchi… ▽ More The goal of this paper is to derive signatures of synchrotron radiation from state-of-the-art simulation models of collimating relativistic magnetohydrodynamic (MHD) jets featuring a large-scale helical magnetic field. We perform axisymmetric special relativistic MHD simulations of the jet acceleration region using the PLUTO code. The computational domain extends from the slow magnetosonic launching surface of the disk up to 6000^2 Schwarzschild radii allowing to reach highly relativistic Lorentz factors. The Poynting dominated disk wind develops into a jet with Lorentz factors of 8 and is collimated to 1 degree. In addition to the disk jet, we evolve a thermally driven spine jet, emanating from a hypothetical black hole corona. Solving the linearly polarized synchrotron radiation transport within the jet, we derive VLBI radio and (sub-) mm diagnostics such as core shift, polarization structure, intensity maps, spectra and Faraday rotation measure (RM), directly from the Stokes parameters. We also investigate depolarization and the detectability of a lambda^2-law RM depending on beam resolution and observing frequency. We find non-monotonic intrinsic RM profiles which could be detected at a resolution of 100 Schwarzschild radii. In our collimating jet geometry, the strict bi-modality in polarization direction (as predicted by Pariev et al.) can be circumvented. Due to relativistic aberration, asymmetries in the polarization vectors across the jet can hint to the spin direction of the central engine. △ Less

Submitted 21 May, 2011; originally announced May 2011.

Comments: Submitted to ApJ

Computing the dust distribution in the bowshock of a fast moving, evolved star

Authors: Allard Jan van Marle, Zakaria Meliani, Rony Keppens, Leen Decin

Abstract: We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dus… ▽ More We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind, and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grainsize, for which we use ten representative grainsize bins. Our simulations allow to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 micro-meters) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate finestructure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer, and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grainsize, this should leave a clear imprint in infrared observations of bowshocks of red supergiants and other evolved stars. △ Less

Submitted 12 May, 2011; originally announced May 2011.

Comments: Accepted for publication in ApJL, 4 figures

Two-shell collisions in the GRB afterglow phase

Authors: A. Vlasis, H. van Eerten, Z. Meliani, R. Keppens

Abstract: Strong optical and radio flares often appear in the afterglow phase of Gamma-Ray Bursts (GRBs). It has been proposed that colliding ultra-relativistic shells can produce these flares. Such consecutive shells can be formed due to the variability in the central source of a GRB. We perform high resolution 1D numerical simulations of late collisions between two ultra-relativistic shells in order to ex… ▽ More Strong optical and radio flares often appear in the afterglow phase of Gamma-Ray Bursts (GRBs). It has been proposed that colliding ultra-relativistic shells can produce these flares. Such consecutive shells can be formed due to the variability in the central source of a GRB. We perform high resolution 1D numerical simulations of late collisions between two ultra-relativistic shells in order to explore these events. We examine the case where a cold uniform shell collides with a self-similar Blandford and McKee shell in a constant density environment and consider cases with different Lorentz factor and energy for the uniform shell. We produce the corresponding on-axis light curves and emission images for the afterglow phase and examine the occurrence of optical and radio flares assuming a spherical explosion and a hard-edged jet scenario. For our simulations we use the Adaptive Mesh Refinement version of the Versatile Advection Code (AMRVAC) coupled to a linear radiative transfer code to calculate synchrotron emission. We find steeply rising flare like behavior for small jet opening angles and more gradual rebrightenings for large opening angles. Synchrotron self-absorption is found to strongly influence the onset and shape of the radio flare. △ Less

Submitted 15 March, 2011; originally announced March 2011.

Comments: Accepted for publication in MNRAS

3-D simulations of shells around massive stars

Authors: Allard Jan van Marle, Rony Keppens, Zakaria Meliani

Abstract: As massive stars evolve, their winds change. This causes a series of hydrodynamical interactions in the surrounding medium. Whenever a fast wind follows a slow wind phase, the fast wind sweeps up the slow wind in a shell, which can be observed as a circumstellar nebula. One of the most striking examples of such an interaction is when a massive star changes from a red supergiant into a Wolf-Rayet… ▽ More As massive stars evolve, their winds change. This causes a series of hydrodynamical interactions in the surrounding medium. Whenever a fast wind follows a slow wind phase, the fast wind sweeps up the slow wind in a shell, which can be observed as a circumstellar nebula. One of the most striking examples of such an interaction is when a massive star changes from a red supergiant into a Wolf-Rayet star. Nebulae resulting from such a transition have been observed around many Wolf-Rayet stars and show detailed, complicated structures owing to local instabilities in the swept-up shells. Shells also form in the case of massive binary stars, where the winds of two stars collide with one another. Along the collision front gas piles up, forming a shell that rotates along with the orbital motion of the binary stars. In this case the shell follows the surface along which the ram pressure of the two colliding winds is in balance. Using the MPI-AMRVAC hydrodynamics code we have made multi-dimensional simulations of these interactions in order to model the formation and evolution of these circumstellar nebulae and explore whether full 3D simulations are necessary to obtain accurate models of such nebulae. △ Less

Submitted 1 February, 2011; originally announced February 2011.

Comments: 5 Pages, 4 figures, Proceedings of the 39th Liege Astrophysical Colloquium, held in Liege 12-16 July 2010

Journal ref: Societe Royale des Sciences de Liege, Bulletin, July 2010, published

Thin shell morphology in the circumstellar medium of massive binaries

Authors: Allard Jan van Marle, Rony Keppens, Zakaria Meliani

Abstract: We investigate the morphology of the collision front between the stellar winds of binary components in two long-period binary systems, one consisting of a hydrogen rich Wolf-Rayet star (WNL) and an O-star and the other of a Luminous Blue Variable (LBV) and an O-star. Specifically, we follow the development and evolution of instabilities that form in such a shell, if it is sufficiently compressed,… ▽ More We investigate the morphology of the collision front between the stellar winds of binary components in two long-period binary systems, one consisting of a hydrogen rich Wolf-Rayet star (WNL) and an O-star and the other of a Luminous Blue Variable (LBV) and an O-star. Specifically, we follow the development and evolution of instabilities that form in such a shell, if it is sufficiently compressed, due to both the wind interaction and the orbital motion. We use MPI-AMRVAC to time-integrate the equations of hydrodynamics, combined with optically thin radiative cooling, on an adaptive mesh 3D grid. Using parameters for generic binary systems, we simulate the interaction between the winds of the two stars. The WNL+O star binary shows a typical example of an adiabatic wind collision. The resulting shell is thick and smooth, showing no instabilities. On the other hand, the shell created by the collision of the O star wind with the LBV wind, combined with the orbital motion of the binary components, is susceptible to thin shell instabilities, which create a highly structured morphology. We identify the nature of the instabilities as both linear and non-linear thin-shell instabilities, with distinct differences between the leading and the trailing parts of the collision front. We also find that for binaries containing a star with a (relatively) slow wind, the global shape of the shell is determined more by the slow wind velocity and the orbital motion of the binary, than the ram pressure balance between the two winds. The interaction between massive binary winds needs further parametric exploration, to identify the role and dynamical importance of multiple instabilities at the collision front, as shown here for an LBV+O star system. △ Less

Submitted 8 November, 2010; originally announced November 2010.

Comments: 10 pages, 13 figures. Accepted for publication in A&A

Dynamics and stability of relativistic GRB blast waves

Authors: Z. Meliani, R. Keppens

Abstract: In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprecedented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov phase. Our goal is to determine the stability properties of its frontal shock. We carried out a grid-adaptive relativistic 2D hydro-simulation at extr… ▽ More In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprecedented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov phase. Our goal is to determine the stability properties of its frontal shock. We carried out a grid-adaptive relativistic 2D hydro-simulation at extreme resolving power, following the GRB jet during the entire afterglow phase. We investigate the effect of the finite initial jet opening angle on the deceleration of the blast wave, and identify the growth of various instabilities throughout the coasting shock front. We find that during the relativistic phase, the blast wave is subject to pressure-ram pressure instabilities that ripple and fragment the frontal shock. These instabilities manifest themselves in the ultra-relativistic phase alone, remain in full agreement with causality arguments, and decay slowly to finally disappear in the near-Newtonian phase as the shell Lorentz factor drops below 3. From then on, the compression rate decreases to levels predicted to be stable by a linear analysis of the Sedov phase. Our simulations confirm previous findings that the shell also spreads laterally because a rarefaction wave slowly propagates to the jet axis, inducing a clear shell deformation from its initial spherical shape. The blast front becomes meridionally stratified, with decreasing speed from axis to jet edge. △ Less

Submitted 7 September, 2010; originally announced September 2010.

Comments: accepted for publication in A&A

Relativistic spine jets from Schwarzschild black holes: "Application to AGN radioloud sources"

Authors: Z. Meliani, C. Sauty, K. Tsinganos, E. Trussoni, V. Cayatte

Abstract: The two types of Fanaroff-Riley radio loud galaxies, FRI and FRII, exhibit strong jets but with different properties. These differences may be associated to the central engine and/or the external medium. Aims: The AGN classification FRI and FRII can be linked to the rate of electromagnetic Poynting flux extraction from the inner corona of the central engine by the jet. The collimation results from… ▽ More The two types of Fanaroff-Riley radio loud galaxies, FRI and FRII, exhibit strong jets but with different properties. These differences may be associated to the central engine and/or the external medium. Aims: The AGN classification FRI and FRII can be linked to the rate of electromagnetic Poynting flux extraction from the inner corona of the central engine by the jet. The collimation results from the distribution of the total electromagnetic energy across the jet, as compared to the corresponding distribution of the thermal and gravitational energies. We use exact solutions of the fully relativistic magnetohydrodynamical (GRMHD) equations obtained by a nonlinear separation of the variables to study outflows from a Schwarzschild black hole corona. A strong correlation is found between the jet features and the energetic distribution of the plasma of the inner corona which may be related to the efficiency of the magnetic rotator. It is shown that observations of FRI and FRII jets may be partially constrained by our model for spine jets. The deceleration observed in FRI jets may be associated with a low magnetic efficiency of the central magnetic rotator and an important thermal confinement by the hot surrounding medium. Conversely, the strongly collimated and accelerated FRII outflows may be self collimated by their own magnetic field because of the high efficiency of the central magnetic rotator. △ Less

Submitted 10 June, 2010; originally announced June 2010.

Comments: Accepted for publication in the A&A

Jet simulations and Gamma-ray burst afterglow jet breaks

Authors: H. J. van Eerten, Z. Meliani, R. A. M. J. Wijers, R. Keppens

Abstract: The conventional derivation of the gamma-ray burst afterglow jet break time uses only the blast wave fluid Lorentz factor and therefore leads to an achromatic break. We show that in general gamma-ray burst afterglow jet breaks are chromatic across the self-absorption break. Depending on circumstances, the radio jet break may be postponed significantly. Using high-accuracy adaptive mesh fluid simul… ▽ More The conventional derivation of the gamma-ray burst afterglow jet break time uses only the blast wave fluid Lorentz factor and therefore leads to an achromatic break. We show that in general gamma-ray burst afterglow jet breaks are chromatic across the self-absorption break. Depending on circumstances, the radio jet break may be postponed significantly. Using high-accuracy adaptive mesh fluid simulations in one dimension, coupled to a detailed synchrotron radiation code, we demonstrate that this is true even for the standard fireball model and hard-edged jets. We confirm these effects with a simulation in two dimensions. The frequency dependence of the jet break is a result of the angle dependence of the emission, the changing optical depth in the self-absorbed regime and the shape of the synchrotron spectrum in general. In the optically thin case the conventional analysis systematically overestimates the jet break time, leading to inferred opening angles that are underestimated by a factor 1.32 and explosion energies that are underestimated by a factor 1.73, for explosions in a homogeneous environment. △ Less

Submitted 21 May, 2010; originally announced May 2010.

Comments: MNRAS submitted. 9 pages, 12 figures

Journal ref: MNRAS 410, 2016 (2010)

Gamma-ray burst afterglows from trans-relativistic blast waves

Authors: H. J. van Eerten, K. Leventis, Z. Meliani, R. A. M. J. Wijers

Abstract: We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysi… ▽ More We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics of shock acceleration, synchrotron self-absorption and electron cooling to a high-performance hydrodynamics simulation. △ Less

Submitted 2 December, 2009; originally announced December 2009.

Comments: small two-page contribution to conference proceedings "The Shocking Universe. Gamma-ray Bursts and High Energy Shock Phenomena". Venice, 14-18 September 2009

Decelerating relativistic two-component jets

Authors: Z. Meliani, R. Keppens

Abstract: Transverse stratification is a common intrinsic feature of astrophysical jets. There is growing evidence that jets in radio galaxies consist of a fast low density outflow at the jet axis, surrounded by a slower, denser, extended jet. The inner and outer jet components then have a different origin and launching mechanism, making their effective inertia, magnetization, associated energy flux and a… ▽ More Transverse stratification is a common intrinsic feature of astrophysical jets. There is growing evidence that jets in radio galaxies consist of a fast low density outflow at the jet axis, surrounded by a slower, denser, extended jet. The inner and outer jet components then have a different origin and launching mechanism, making their effective inertia, magnetization, associated energy flux and angular momentum content different as well. Their interface will develop differential rotation, where disruptions may occur. We here investigate the stability of rotating, two-component relativistic outflows typical for jets in radio galaxies. For this purpose, we parametrically explore the long term evolution of a transverse cross-section of radially stratified jets numerically With grid-adaptive relativistic MHD simulations, augmented with approximate linear stability analysis. We study the influence of dynamically important poloidal magnetic fields, with varying contributions of the inner component jet to the total kinetic energy flux of the jet, on their non-linear azimuthal stability. We demonstrate that two-component jets with high kinetic energy flux, and an inner jet effective inertia which is higher than the outer jet effective inertia are subject to the development of a relativistically enhanced, rotation-induced Rayleigh-Taylor type instability. This instability plays a major role in decelerating the inner jet and the overall jet decollimation. This novel deceleration scenario can partly explain the radio source dichotomy, relating it directly to the efficiency of the central engine in launching the inner jet component. The FRII/FRI transition could then occur when the relative kinetic energy flux of the inner to the outer jet grows beyond a certain treshold. △ Less

Submitted 2 October, 2009; originally announced October 2009.

Comments: 18 pages, ApJ accepted

Journal ref: Astrophys.J.705:1594-1606,2009

Gamma-ray burst afterglows from trans-relativistic blast wave simulations

Authors: H. J. van Eerten, K. Leventis, Z. Meliani, R. A. M. J. Wijers, R. Keppens

Abstract: We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics… ▽ More We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics of shock acceleration, synchrotron self-absorption and electron cooling to a high-performance hydrodynamics simulation. For the dynamics we find that the transition to the nonrelativistic regime generally occurs later than expected, that the Sedov-Taylor solution overpredicts the late time blast wave radius and that the analytical formula for the blast wave velocity from Huang (1999) overpredicts the late time velocity by a factor 4/3. For the radiation we find that the flux may differ up to an order of magnitude depending on the equation of state that is used for the fluid and that the counterjet leads to a clear rebrightening at late times for hard-edged jets. Simulating GRB030329 using predictions for its physical parameters from the literature leads to spectra and light curves that may differ significantly from the actual data, emphasizing the need for very accurate modelling. Predicted light curves at low radio frequencies for a hard-edged jet model of GRB030329 with opening angle 22 degrees show typically two distinct peaks, due to the combined effect of jet break, non relativistic break and counterjet. △ Less

Submitted 14 July, 2010; v1 submitted 13 September, 2009; originally announced September 2009.

Comments: 18 pages, 19 figures. Replaced with MNRAS accepted version, minor changes only

Journal ref: Mon.Not.Roy.Astron.Soc.403:300,2010

No visible optical variability from a relativistic blast wave encountering a wind-termination shock

Authors: H. J. van Eerten, Z. Meliani, R. A. M. J. Wijers, R. Keppens

Abstract: Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high perfor… ▽ More Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high performance adaptive mesh relativistic hydrodynamics code, AMRVAC, and the synchrotron emission is analyzed in detail with a separate radiation code. We find no bump in the resulting light curve, not even for very high density jumps. Furthermore, by analyzing the contributions from the different shock wave regions we are able to establish that it is essential to resolve the blast wave structure in order to make qualitatively correct predictions on the observed output and that the contribution from the reverse shock region will not stand out, even when the magnetic field is increased in this region by repeated shocks. This study resolves a controversy in recent literature. △ Less

Submitted 19 June, 2009; originally announced June 2009.

Comments: 4 figures, submitted to MNRAS letters

Journal ref: MNRAS (2009) 398, L63

Linear wave propagation in relativistic magnetohydrodynamics

Authors: R. Keppens, Z. Meliani

Abstract: The properties of linear Alfvén, slow, and fast magnetoacoustic waves for uniform plasmas in relativistic magnetohydrodynamics (MHD) are discussed, augmenting the well-known expressions for their phase speeds with knowledge on the group speed. A 3+1 formalism is purposely adopted to make direct comparison with the Newtonian MHD limits easier and to stress the graphical representation of their an… ▽ More The properties of linear Alfvén, slow, and fast magnetoacoustic waves for uniform plasmas in relativistic magnetohydrodynamics (MHD) are discussed, augmenting the well-known expressions for their phase speeds with knowledge on the group speed. A 3+1 formalism is purposely adopted to make direct comparison with the Newtonian MHD limits easier and to stress the graphical representation of their anisotropic linear wave properties using the phase and group speed diagrams. By drawing these for both the fluid rest frame and for a laboratory Lorentzian frame which sees the plasma move with a three-velocity having an arbitrary orientation with respect to the magnetic field, a graphical view of the relativistic aberration effects is obtained for all three MHD wave families. Moreover, it is confirmed that the classical Huygens construction relates the phase and group speed diagram in the usual way, even for the lab frame viewpoint. Since the group speed diagrams correspond to exact solutions for initial conditions corresponding to a localized point perturbation, their formulae and geometrical construction can serve to benchmark current high-resolution algorithms for numerical relativistic MHD. △ Less

Submitted 14 October, 2008; originally announced October 2008.

Comments: accepted in Physics of Plasmas

Faranoff-Riley type I jet deceleration at density discontinuities "Relativistic hydrodynamics with realistic equation of state"

Authors: Z. Meliani, R. Keppens, B. Giacomazzo

Abstract: The deceleration mechanisms for relativistic jets in active galactic nuclei remain an open question, and in this paper we propose a model which could explain sudden jet deceleration, invoking density discontinuities. This is particularly motivated by recent indications from HYMORS. Exploiting high resolution, numerical simulations, we demonstrate that for both high and low energy jets, always at… ▽ More The deceleration mechanisms for relativistic jets in active galactic nuclei remain an open question, and in this paper we propose a model which could explain sudden jet deceleration, invoking density discontinuities. This is particularly motivated by recent indications from HYMORS. Exploiting high resolution, numerical simulations, we demonstrate that for both high and low energy jets, always at high Lorentz factor, a transition to a higher density environment can cause a significant fraction of the directed jet energy to be lost on reflection. This can explain how one-sided jet deceleration and a transition to FR I type can occur in HYMORS, which start as FR II (and remain so on the other side). For that purpose, we implemented in the relativistic hydrodynamic grid-adaptive AMRVAC code, the Synge-type equation of state introduced in the general polytropic case by Meliani et al. (2004). We present results for 10 model computations, varying the inlet Lorentz factor from 10 to 20, including uniform or decreasing density profiles, and allowing for cylindrical versus conical jet models. As long as the jet propagates through uniform media, we find that the density contrast sets most of the propagation characteristics, fully consistent with previous modeling efforts. When the jet runs into a denser medium, we find a clear distinction in the decelaration of high energy jets depending on the encountered density jump. For fairly high density contrast, the jet becomes destabilised and compressed, decelerates strongly (up to subrelativistic speeds) and can form knots. We point out differences that are found between cylindrical and conical jet models, together with dynamical details like the Richtmyer-Meshkov instabilities developing at the original contact interface. △ Less

Submitted 18 August, 2008; originally announced August 2008.

Comments: accepted in A&A

A multidimensional grid-adaptive relativistic magnetofluid code

Authors: B. van der Holst, R. Keppens, Z. Meliani

Abstract: A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing Lax-Friedrichs scheme in a finite volume setting and is combined with a diffusive approach for controlling magnetic monopole errors. The consistency between the… ▽ More A robust second order, shock-capturing numerical scheme for multi-dimensional special relativistic magnetohydrodynamics on computational domains with adaptive mesh refinement is presented. The base solver is a total variation diminishing Lax-Friedrichs scheme in a finite volume setting and is combined with a diffusive approach for controlling magnetic monopole errors. The consistency between the primitive and conservative variables is ensured at all limited reconstructions and the spatial part of the four velocity is used as a primitive variable. Demonstrative relativistic examples are shown to validate the implementation. We recover known exact solutions to relativistic MHD Riemann problems, and simulate the shock-dominated long term evolution of Lorentz factor 7 vortical flows distorting magnetic island chains. △ Less

Submitted 4 July, 2008; originally announced July 2008.

Comments: accepted for publication in Computer Physics Communications

Extragalactic jets with helical magnetic fields: relativistic MHD simulations

Authors: R. Keppens, Z. Meliani, B. van der Holst, F. Casse

Abstract: Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind… ▽ More Extragalactic jets are inferred to harbor dynamically important, organized magnetic fields which presumably aid in the collimation of the relativistic jet flows. We here explore by means of grid-adaptive, high resolution numerical simulations the morphology of AGN jets pervaded by helical field and flow topologies. We concentrate on morphological features of the bow shock and the jet beam behind the Mach disk, for various jet Lorentz factors and magnetic field helicities. We investigate the influence of helical magnetic fields on jet beam propagation in overdense external medium. We use the AMRVAC code, employing a novel hybrid block-based AMR strategy, to compute ideal plasma dynamics in special relativity. The helicity of the beam magnetic field is effectively transported down the beam, with compression zones in between diagonal internal cross-shocks showing stronger toroidal field regions. In comparison with equivalent low-relativistic jets which get surrounded by cocoons with vortical backflows filled by mainly toroidal field, the high speed jets demonstrate only localized, strong toroidal field zones within the backflow vortical structures. We find evidence for a more poloidal, straight field layer, compressed between jet beam and backflows. This layer decreases the destabilizing influence of the backflow on the jet beam. In all cases, the jet beam contains rich cross-shock patterns, across which part of the kinetic energy gets transferred. For the high speed reference jet considered here, significant jet deceleration only occurs beyond distances exceeding ${\cal O}(100 R_j)$, as the axial flow can reaccelerate downstream to the internal cross-shocks. This reacceleration is magnetically aided, due to field compression across the internal shocks which pinch the flow. △ Less

Submitted 14 February, 2008; originally announced February 2008.

Comments: 16 pages, Astronomy and Astrophysics accepted for publication

Transverse stability of relativistic two-component jets

Authors: Z. Meliani, R. Keppens

Abstract: Context: Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius represents a second interface where disruptions may occur. Aims: We explore the stability of stratified, rotating, relativi… ▽ More Context: Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius represents a second interface where disruptions may occur. Aims: We explore the stability of stratified, rotating, relativistic two-component jets, in turn embedded in static interstellar medium. Methods: In a grid-adaptive relativistic hydrodynamic simulation with the AMRVAC code, the non-linear azimuthal stability of two-component relativistic jets is investigated. We simulate until multiple inner jet rotations have been completed. Results: We find evidence for the development of an extended shear flow layer between the two jet components, resulting from the growth of a body mode in the inner jet, Kelvin-Helmholtz surface modes at their original interface, and their nonlinear interaction. Both wave modes are excited by acoustic waves which are reflected between the symmetry axis and the interface of the two jet components. Their interaction induces the growth of near stationary, counterrotating vortices at the outer edge of the shear flow layer. The presence of a heavy external jet allows to slow down their further development, and maintain a collimated flow. At the outer jet boundary, small-scale Rayleigh-Taylor instabilities develop, without disrupting the jet configuration. Conclusion: We demonstrate that the cross-section of two-component relativistic jets, with a heavy, cold outer jet, is non-linearly stable. △ Less

Submitted 24 September, 2007; originally announced September 2007.

Comments: Accepted in A&A 24/09/2007

Two-flow magnetohydrodynamical jets around young stellar objects

Authors: Fabien Casse, Zakaria Meliani, C. Sauty

Abstract: We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. Our simulations show that the inner outflow is accelerated from the central object hot corona thanks… ▽ More We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. Our simulations show that the inner outflow is accelerated from the central object hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force. We also show that when a dense inner stellar wind occurs, the resulting disc-driven jet have a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass ejection rate in the disc-driven jet and a larger radial extension which is in better agreement with the observations besides being more consistent. △ Less

Submitted 22 June, 2007; originally announced June 2007.

Comments: Accepted for publication in Astrophysics & Space Science. Referred proceeding of the fifth Mont Stromlo Symposium Dec. 1-8 2006, Canberra, Australia. 5 pages, 3 figures. For high resolution version of the paper, please click here http://www.apc.univ-paris7.fr/~fcasse/publications.html

Journal ref: Astrophys.SpaceSci.311:57-61,2007

GRB blastwaves through wind-shaped circumburst media

Authors: Z. Meliani, R. Keppens

Abstract: Context:A significant fraction of progenitors for long gamma-ray bursts (GRBs) are believed to be massive stars. The investigation of long GRBs therefore requires modeling the propagation of ultra-relativistic blastwaves through the circumburst medium surrounding massive stars. We simulate the expansion of an isotropic, adiabatic relativistic fireball into the wind-shaped medium around a massive… ▽ More Context:A significant fraction of progenitors for long gamma-ray bursts (GRBs) are believed to be massive stars. The investigation of long GRBs therefore requires modeling the propagation of ultra-relativistic blastwaves through the circumburst medium surrounding massive stars. We simulate the expansion of an isotropic, adiabatic relativistic fireball into the wind-shaped medium around a massive GRB progenitor. The circumburst medium is composed of a realistically stratified stellar wind zone up to its termination shock, followed by a region of shocked wind characterized by a constant density. Aims: We followed the evolution of the blastwave through all its stages, including the extremely rapid acceleration up to a Lorentz factor 75 flow, its deceleration by interaction with stellar wind, its passage of the wind termination shock, until its propagation through shocked wind. Methods: We used the adaptive mesh refinement versatile advection code to follow the evolution of the fireball. Results: We show that the acceleration from purely thermal to ultra-relativistic kinetic regimes is abrupt and produces an internally structured blastwave. We resolved the structure of this ultra-relativistic shell in all stages, thanks to the adaptive mesh. We comment on the dynamical roles played by forward and reverse shock pairs in the phase of interaction with the free stellar wind and clearly identify the complex shock-dominated structure created when the shell crosses the terminal shock. Conclusion: We show that in our model where the terminal shock is taken relatively close to the massive star, the phase of self-similar deceleration of Blandford-McKee type can only be produced in the constant density, shocked wind zone. △ Less

Submitted 19 April, 2007; originally announced April 2007.

Comments: Accepted in A&A Letter (04/04/2007)

AMRVAC and Relativistic Hydrodynamic simulations for GRB afterglow phases

Authors: Z. Meliani, R. Keppens, F. Casse, D. Giannios

Abstract: We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the int… ▽ More We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $γ=100$ up to the almost Newtonian $γ\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart. △ Less

Submitted 15 January, 2007; originally announced January 2007.

Comments: 12 pages, accepted in MNRAS, 12/01/2007

Journal ref: Mon.Not.Roy.Astron.Soc.376:1189-1200,2007

Two-component magnetohydrodynamical outflows around young stellar objects Interplay between stellar magnetospheric winds and disc-driven jets

Authors: Z. Meliani, F. Casse, C. Sauty

Abstract: We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar magnetospheric winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. These innovative MHD simulations are devoted to the study of the interplay between a stellar wind (having different ejection mass rates) and an MHD disc-driven jet embedding the ste… ▽ More We present the first-ever simulations of non-ideal magnetohydrodynamical (MHD) stellar magnetospheric winds coupled with disc-driven jets where the resistive and viscous accretion disc is self-consistently described. These innovative MHD simulations are devoted to the study of the interplay between a stellar wind (having different ejection mass rates) and an MHD disc-driven jet embedding the stellar wind. The transmagnetosonic, collimated MHD outflows are investigated numerically using the VAC code. We first investigate the various angular momentum transports occurring in the magneto-viscous accretion disc. We then analyze the modifications induced by the interaction between the two components of the outflow. Our simulations show that the inner outflow is accelerated from the central object's hot corona thanks to both the thermal pressure and the Lorentz force. In our framework, the thermal acceleration is sustained by the heating produced by the dissipated magnetic energy due to the turbulence. Conversely, the outflow launched from the resistive accretion disc is mainly accelerated by the magneto-centrifugal force.}{The simulations show that the MHD disc-driven outflow extracts angular momentum more efficiently than do viscous effects in near-equipartition, thin-magnetized discs where turbulence is fully developed. We also show that, when a dense inner stellar wind occurs, the resulting disc-driven jet has a different structure, namely a magnetic structure where poloidal magnetic field lines are more inclined because of the pressure caused by the stellar wind. This modification leads to both an enhanced mass-ejection rate in the disc-driven jet and a larger radial extension that is in better agreement with the observations, besides being more consistent. △ Less

Submitted 31 August, 2006; v1 submitted 28 August, 2006; originally announced August 2006.

Comments: 16 pages, Accepted in A&A 04/08/2006