Search

We present synthetic spectra corresponding to a 2.5D magnetohydrodynamical simulation of a rotating prominence in the Ca II 8542 \r{A}, H$\alpha$, Ca II K, Mg II k, Ly $\alpha$, and Ly $\beta$ lines. The prominence rotation resulted from angular momentum conservation within a flux rope where asymmetric heating imposed a net rotation prior to the thermal-instability driven condensation phase. The spectra were created using a library built on the Lightweaver framework called Promweaver, which provides boundary conditions for incorporating the limb-darkened irradiation of the solar disk on isolated structures such as prominences. Our spectra show distinctive rotational signatures for the Mg II k, Ly $\alpha$, and Ly $\beta$ lines, even in the presence of complex, turbulent solar atmospheric conditions. However, these signals are hardly detectable for the Ca II 8542 \r{A}, H$\alpha$, Ca II K spectral lines. Most notably we find only a very faint rotational signal in the H$\alpha$ line, thus reigniting the discussion on the existence of sustained rotation in prominences., Comment: Accepted in A&A Letters. Movies available at https://cloud.aip.de/index.php/s/cs9RKM4CwLrwR3N

Prominences and coronal rain are two forms of coronal condensations for which we still lack satisfactory details on the formation pathways and conditions under which the two come to exist. We compared prominences that formed via a steady versus stochastic type of heating. We performed 2.5D simulations using the open-source MPI-AMRVAC code. To further extend the work and allow for future direct comparison with observations, we used Lightweaver to form spectra of the filament view of our steady case prominence. With that, we analysed a reconnection event that shares certain characteristics with nanojets. We show how different forms of localised heating that induce thermal instability result in prominences with different properties. The steady form of heating results in prominence with a clear vertical structure stretching across the magnetic field lines. On the other hand, stochastic heating produces many threads that predominantly have a horizontal motion along the field lines. In the steady heating case, the prominence is relatively static; however, there is evidence of reconnection happening almost the entire time the prominence is present. In the case of stochastic heating, the threads are highly dynamic, with them also exhibiting a form of transverse oscillation (strongly resembling the decayless type). The fact that the threads in the stochastic heating case are constantly moving along the field lines suppresses any conditions for reconnection. It, therefore, appears that, to first order, the choice of heating prescription defines whether the prominence-internal dynamics are oriented vertically or horizontally. We closely inspected a sample reconnection event and computed the synthetic optically thick radiation using the open-source Lightweaver radiative transfer framework., Comment: 16 pages, 14 figures, accepted for publication

The aim of this paper is to present the implementation of two new CME models in the space weather forecasting tool, EUHFORIA. We introduce the two toroidal CME models analytically, along with their numerical implementation in EUHFORIA. One model is based on the modified Miller-Turner (mMT) solution, while the other is derived from the Soloviev equilibrium, a specific solution of the Grad-Shafranov equation. The magnetic field distribution in both models is provided in analytic formulae, enabling a swift numerical computation. After detailing the differences between the two models, we present a collection of thermodynamic and magnetic profiles obtained at Earth using these CME solutions in EUHFORIA with a realistic solar wind background. Subsequently, we explore the influence of their initial parameters on the time profiles at L1. In particular, we examine the impact of the initial density, magnetic field strength, velocity, and minor radius. In EUHFORIA, we obtained different thermodynamic and magnetic profiles depending on the CME model used. We found that changing the initial parameters affects both the amplitude and the trend of the time profiles. For example, using a high initial speed results in a fast evolving and compressed magnetic structure. The speed of the CME is also linked to the strength of the initial magnetic field due to the contribution of the Lorentz force on the CME expansion. However, increasing the initial magnetic field also increases the computation time. Finally, the expansion and integrity of the magnetic structure can be controlled via the initial density of the CME. Both toroidal CME models are successfully implemented in EUHFORIA and can be utilized to predict the geo-effectiveness of the impact of real CME events. Moreover, the current implementation could be easily modified to model other toroidal magnetic configurations.

The density decreases exponentially with height in the solar gravitationally stratified atmosphere, therefore the collisional coupling between the ionized plasma and the neutrals also decreases. Here, we investigate the role of collisions between ions and neutrals on the reconnection process occurring at various heights in the atmosphere. We perform simulations of magnetic reconnection induced by a localized resistivity in a gravitationally stratified atmosphere, where we vary the height of the initial reconnection X-point. We compare a magnetohydrodynamic (MHD) model and two two-fluid configurations: one where the collisional coupling is calculated from local plasma parameters and another where the coupling is decreased, so that collisional effects are enhanced. Simulations in a stratified atmosphere show similar structures in MHD and two-fluid simulations with strong coupling. However, when collisional effects are increased to attain representative parameter regimes, we find a nonlinear runaway instability, which separates the plasma-neutral densities across the current sheet (CS). With increased collisional effects, the initial decoupling in velocity heats the neutrals and this sets up a nonlinear feedback where neutrals migrate outside the CS, replacing charged particles which accumulate towards the center of the CS. The reconnection rate has a maximum value around 0.1, similar for both reconnection heights, and is consistent with the use of a localized enhanced resistivity used in all three models. The initial stages of plasmoid formation, observed near the end of our simulations, is influenced by the outflow from the primary reconnection point, rather than by collisions. We synthesize optically thin emission for both MHD and two-fluid models, which can show a very different evolution when the charged particle density is used instead of the total density.

With this study, we aim to understand the nature of prominences, governed by their formation process. We use a state-of-the-art threaded prominence model within a dipped magnetic arcade. The non-ideal magnetohydrodynamic (MHD) equations are solved using the open-source MPI-AMRVAC MHD toolkit. Unlike many previous 1D models, we study the full 2D dynamics in a fixed-shaped arcade. This allows for sideways field deformations and cross-field thermodynamic coupling. To achieve a realistic setup we consider field-aligned thermal conduction, radiative cooling and heating, wherein the latter combines a steady background and a localized stochastic component. The stochastic component simulates energy pulses localized in time and space at the footpoints of the magnetic arcade. We vary the height and amplitude of the localized heating and observe how it influences the prominence, its threads, and its overall dynamics. We show the importance of random localized heating in the evolution of prominences and their threaded structure. Random heating strongly influences the morphology of the prominence threaded structure, the area, the mass the threads reach, their minimum temperature and their average density. More importantly, the strength of the localized heating plays a role in maintaining the balance between condensation and draining, affecting the general prominence stability. Stronger sources form condensations faster and result in larger and more massive prominences. We show how the condensation rates scale with the amplitude of the heating inputs and quantify how these rates match with values from observations. We detail how stochastic sources determine counterstreaming flows and oscillations of prominence threads., Comment: 13 pages, 12 figures, accepted for publication

Aims. We here demonstrate how the recently developed Lightweaver framework makes non-LTE (NLTE) spectral synthesis feasible on a new 3D ab-initio magnetohydrodynamic (MHD) filament/prominence simulation, in a post-processing step. Methods. We clarify the need to introduce filament/prominent-specific Lightweaver boundary conditions that accurately model incident chromospheric radiation, and include a self-consistent and smoothly varying limb darkening function. Results. Progressing from isothermal/isobaric models to the self-consistently generated stratifications within a fully 3D MHD filament/prominence simulation, we find excellent agreement between our 1.5D non local thermodynamic equilibrium Lightweaver synthesis and a popular Hydrogen H{\alpha} proxy. We compute additional lines including Ca~\textsc{ii} 8542 alongside the more optically-thick Ca~\textsc{ii} H&K & Mg~\textsc{ii} h&k lines, for which no comparable proxy exists, and explore their formation properties within filament/prominence atmospheres. Conclusions. The versatility of the Lightweaver framework is demonstrated with this extension to 1.5D filament/prominence models, where each vertical column of the instantaneous 3D MHD state is spectrally analysed separately, without accounting for (important) multi-dimensional radiative effects. The general agreement found in the line core contrast of both observations and the Lightweaver-synthesised simulation further validates the current generation of solar filaments/prominences models constructed numerically with MPI-AMRVAC., Comment: 16+2 pages with appendix, 12+1 figures, accepted for publication in Astronomy & Astrophysics

Aims. We begin by exploring the influence of two classes of commonly adopted heating models on the formation behaviour of solar prominences. These models consider either an exponential variation dependent on height alone, or local density and magnetic field conditions. We highlight and address some of the limitations inherent to these early approximations by proposing a new, dynamic 2D flux rope heating model that qualitatively accounts for the 3D topology of the twisted flux rope field. Methods. We performed 2.5D grid-adaptive numerical simulations of prominence formation via the levitation-condensation mechanism. A linear force-free arcade is subjected to shearing and converging motions, leading to the formation of a flux rope containing material that may succumb to thermal instability. The eventual formation and subsequent evolution of prominence condensations was then quantified as a function of the specific background heating prescription adopted. For the simulations that consider the topology of the flux rope, reduced heating was considered within a dynamically evolving ellipse that traces the flux rope cross-section. This ellipse is centred on the flux rope axis and tracked during runtime using an approach based on the instantaneous magnetic field curvature. Results. We find that the nature of the heating model is clearly imprinted on the evolution and morphology of any resulting prominences: one large, low-altitude condensation is obtained for the heating model based on local parameters, while the exponential model leads to the additional formation of smaller blobs throughout the flux rope which then relocate as they tend towards achieving hydrostatic equilibrium. Finally, a study of the condensation process in phase space reveals a non-isobaric evolution with an eventual recovery of uniform pressure balance along flux surfaces., Comment: 19 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Magnetohydrodynamic turbulence is central to laboratory and astrophysical plasmas, and is invoked for interpreting many observed scalings. Verifying predicted scaling law behaviour requires extreme-resolution direct numerical simulations (DNS), with needed computing resources excluding systematic parameter surveys. We here present an analytic generator of realistically looking turbulent magnetic fields, that computes 3D ${\cal{O}}(1000^3)$ solenoidal vector fields in minutes to hours on desktops. Our model is inspired by recent developments in 3D incompressible fluid turbulence theory, where a Gaussian white noise vector subjected to a non-linear transformation results in an intermittent, multifractal random field. Our $B\times C$ model has only few parameters that have clear geometric interpretations. We directly compare a (costly) DNS with a swiftly $B\times C$-generated realization, in terms of its (i) characteristic sheet-like structures of current density, (ii) volume-filling aspects across current intensity, (iii) power-spectral behaviour, (iv) probability distribution functions of increments for magnetic field and current density, structure functions, spectra of exponents, and (v) partial variance of increments. The model even allows to mimic time-evolving magnetic and current density distributions and can be used for synthetic observations on 3D turbulent data cubes., Comment: 18 pages, 14 figures, Accepted for publication in Physical Review E

The chromosphere is a partially ionized layer of the solar atmosphere, the transition between the photosphere where the gas is almost neutral and the fully ionized corona. As the collisional coupling between neutral and charged particles decreases in the upper part of the chromosphere, the hydrodynamical timescales may become comparable to the collisional timescale, and a two-fluid model is needed. In this paper we describe the implementation and validation of a two-fluid model which simultaneously evolves charges and neutrals, coupled by collisions. The two-fluid equations are implemented in the fully open-source MPI-AMRVAC code. In the photosphere and the lower part of the solar atmosphere, where collisions between charged and neutral particles are very frequent, an explicit time-marching would be too restrictive, since for stability the timestep needs to be proportional to the inverse of the collision frequency. This is overcome by evaluating the collisional terms implicitly using an explicit-implicit (IMEX) scheme. The cases presented cover very different collisional regimes and our results are fully consistent with related literature findings. If collisional time and length scales are smaller than the hydrodynamical scales usually considered in the solar chromosphere, density structures seen in the neutral and charged fluids are similar, with the effect of elastic collisions between charges and neutrals being similar to diffusivity. Otherwise, density structures are different and the decoupling in velocity between the two species increases. The use of IMEX schemes efficiently avoids the small timestep constraints of fully explicit implementations in strongly collisional regimes. Adaptive Mesh Refinement (AMR) greatly decreases the computational cost, compared to uniform grid runs at the same effective resolution., Comment: accepted for publication in A&A

In this work, we study the causal relations between a localised energy release and a remote prominence oscillation, where the prominence has a realistic thread-like structure. We used an open source magnetohydrodynamic (MHD) code known as MPI-AMRVAC to create a multithreaded prominence body. We introduced an additional energy source from which a shock wave originates, thereby inducing prominence oscillation. We studied two cases with different source amplitudes to analyze its effect on the oscillations. Our results show that the frequently used pendulum model does not suffice to fully estimate the period of the prominence oscillation, in addition to showing that the influence of the source and the thread-like prominence structure needs to be taken into account. Repeated reflections and transmissions of the initial shock wave occur at the specific locations of multiple high-temperature and high-density gradients in the domain. This includes the left and right transition region (TR) located at the footpoints of the magnetic arcade, as well as the various transition regions between the prominence and the corona (PCTR). This results in numerous interferences of compressional waves. They contribute to the restoring forces of the oscillation, causing the period to deviate from the expected pendulum model, in addition to leading to differences in attributed damping or even growth in amplitude between the various threads. Along with the global longitudinal motion that result from the shock impact, small-scale transverse oscillations are also evident. Multiple high-frequency oscillations represent the propagation of magnetoacoustic waves. The damping we see is linked to the conversion of energy and its exchange with the surrounding corona. Our simulations demonstrate the exchange of energy between different threads and their different modes of oscillation., Comment: 13 pages, 11 figures, accepted for publication

In this paper, we revisit the governing equations for linear magnetohydrodynamic (MHD) waves and instabilities existing within a magnetized, plane-parallel, self-gravitating slab. Our approach allows for fully non-uniformly magnetized slabs, which deviate from isothermal conditions, such that the well-known Alfv\'en and slow continuous spectra enter the description. We generalize modern MHD textbook treatments, by showing how self-gravity enters the MHD wave equation, beyond the frequently adopted Cowling approximation. This clarifies how Jeans' instability generalizes from hydro to magnetohydrodynamic conditions without assuming the usual Jeans' swindle approach. Our main contribution lies in reformulating the completely general governing wave equations in a number of mathematically equivalent forms, ranging from a coupled Sturm-Liouville formulation, to a Hamiltonian formulation linked to coupled harmonic oscillators, up to a convenient matrix differential form. The latter allows us to derive analytically the eigenfunctions of a magnetized, self-gravitating thin slab. In addition, as an example we give the exact closed form dispersion relations for the hydrodynamical p- and Jeans-unstable modes, with the latter demonstrating how the Cowling approximation modifies due to a proper treatment of self-gravity. The various reformulations of the MHD wave equation open up new avenues for future MHD spectral studies of instabilities as relevant for cosmic filament formation, which can e.g. use modern formal solution strategies tailored to solve coupled Sturm-Liouville or harmonic oscillator problems., Comment: 20 pages, 1 figure, accepted for publication in MNRAS

Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavour. We develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. We build upon the finite volume code MPI-AMRVAC to solve the equations of hydrodynamics on multi-dimensional adaptive meshes and introduce a new module to handle the coupling with radiation. A non-equilibrium, flux-limiting diffusion approximation is used to close the radiation momentum and energy equations. The time-dependent radiation energy equation is then solved within a flexible framework, accounting fully for radiation forces and work terms and further allowing the user to adopt a variety of descriptions for the radiation-matter interaction terms (the 'opacities'). We validate the radiation module on a set of standard testcases for which different terms of the radiative energy equation predominate. As a preliminary application to a scientific case, we calculate spherically symmetric models of the radiation-driven and optically thick supersonic outflows from massive Wolf-Rayet stars. This also demonstrates our code's flexibility, as the illustrated simulation combines opacities typically used in static stellar structure models with a parametrised form for the enhanced line-opacity expected in supersonic flows. This new module provides a convenient and versatile tool to perform multi-dimensional and high resolution radiative-hydrodynamics simulations in optically thick environments with the MPI-AMRVAC code. The code is ready to be used for a variety of astrophysical applications, where a first target for us will be multi-dimensional simulations of stellar outflows from Wolf-Rayet stars., Comment: 14 pages, 10 figures Submitted to A&A

The class of Double-Double Radio Galaxies (DDRGs) relates to episodic jet outbursts. How various regions and components add to the total intensity in radio images is less well known. In this paper we synthesize synchrotron images for DDRGs based on special relativistic hydrodynamic simulations, making advanced approximations for the magnetic fields. We study the synchrotron images for: Three different radial jet profiles; Ordered, entangled or mixed magnetic fields; Spectral ageing from synchrotron cooling; The contribution from different jet components; The viewing angle and Doppler (de-)boosting; The various epochs of the evolution of the DDRG. To link our results to observational data, we adopt to J1835+6204 as a reference source. In all cases the synthesized synchrotron images show two clear pairs of hotspots, in the inner and outer lobes. The best resemblance is obtained for the piecewise isochoric jet model, for a viewing angle of approximately $\vartheta \sim -71^{\circ}$, i.e. inclined with the lower jet towards the observer, with predominantly entangled ($\gtrsim 70$ per cent of the magnetic pressure) in turbulent, rather than ordered fields. The effects of spectral ageing become significant when the ratio of observation frequencies and cut-off frequency $\nu_{\rm obs}/\nu_{\infty,0} \gtrsim 10^{-3}$, corresponding to $\sim 3 \cdot 10^2$ MHz. For viewing angles $\vartheta \lesssim -30^{\circ}$, a DDRG morphology can no longer be recognized. The second jets must be injected within $\lesssim$ 4 per cent of the lifetime of the first jets for a DDRG structure to emerge, which is relevant for Active Galactic Nuclei feedback constraints., Comment: 24 pages, 8 figures

The formation of shocks within the solar atmosphere remains one of the few observable signatures of energy dissipation arising from the plethora of magnetohydrodynamic waves generated close to the solar surface. Active region observations offer exceptional views of wave behavior and its impact on the surrounding atmosphere. The stratified plasma gradients present in the lower solar atmosphere allow for the potential formation of many theorized shock phenomena. In this study, using chromospheric Ca II 854.2nm spectropolarimetric data of a large sunspot, we examine fluctuations in the plasma parameters in the aftermath of powerful shock events that demonstrate polarimetric reversals during their evolution. Modern inversion techniques are employed to uncover perturbations in the temperatures, line-of-sight velocities, and vector magnetic fields occurring across a range of optical depths synonymous with the shock formation. Classification of these non-linear signatures is carried out by comparing the observationally-derived slow, fast, and Alfv\'en shock solutions to the theoretical Rankine-Hugoniot relations. Employing over 200,000 independent measurements, we reveal that the Alfv\'en (intermediate) shock solution provides the closest match between theory and observations at optical depths of log(tau) = -4, consistent with a geometric height at the boundary between the upper photosphere and lower chromosphere. This work uncovers first-time evidence of the manifestation of chromospheric intermediate shocks in sunspot umbrae, providing a new method for the potential thermalization of wave energy in a range of magnetic structures, including pores, magnetic flux ropes, and magnetic bright points., Comment: 17 pages, 9 figures, accepted by ApJ

The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars.

Ultra-luminous X-ray sources (ULX) have so high X-ray luminosities that they were long thought to be accreting intermediate mass black holes. Yet, some ULX have been shown to display periodic modulations and coherent pulsations, suggestive of a neutron star in orbit around a stellar companion and accreting at super-Eddington rates. In this letter, we propose that the mass transfer in ULX could be qualitatively the same as in Supergiant X-ray binaries (SgXB), with a wind from the donor star highly beamed towards the compact object. Since the star does not fill its Roche lobe, this mass transfer mechanism known as "wind Roche lobe overflow" can remain stable even for large mass ratios. Based on realistic acceleration profiles derived from spectral observations and modeling of the stellar wind, we compute the bulk motion of the wind to evaluate the fraction of the stellar mass outflow captured by the compact object. We identify the orbital and stellar conditions for a SgXB to transfer mass at rates matching the expectations for ULX and show that the transition from SgXB to ULX luminosity levels is progressive. These results indicate that a high stellar Roche lobe filling factor is not necessary to funnel large quantities of material into the Roche lobe of the accretor. Large stellar mass loss rates such as the ones from the Wolf-Rayet star in M101 ULX-1 or the late B9 Supergiant in NGC 7793 P13 are enough to lead to a highly beamed wind and a significantly enhanced mass transfer rate.

In Supergiant X-ray binaries (SgXB), a compact object captures a fraction of the wind of an O/B supergiant on a close orbit. Proxies exist to evaluate the efficiency of mass and angular momentum accretion but they depend so dramatically on the wind speed that given the current uncertainties, they only set loose constrains. Furthermore, they often bypass the impact of orbital and shock effects on the flow structure. We study the wind dynamics and the angular momentum gained as the flow is accreted. We identify the conditions for the formation of a disc-like structure around the accretor and the observational consequences for SgXB. We use recent results on the wind launching mechanism to compute 3D streamlines, accounting for the gravitational and X-ray ionizing influence of the compact companion on the wind. Once the flow enters the Roche lobe of the accretor, we solve the hydrodynamics equations with cooling. A shocked region forms around the accretor as the flow is beamed. For wind speeds of the order of the orbital speed, the shock is highly asymmetric compared to the axisymmetric bow shock obtained for a purely planar homogeneous flow. With net radiative cooling, the flow always circularizes for wind speeds low enough. Although the donor star does not fill its Roche lobe, the wind can be significantly beamed and bent by the orbital effects. The net angular momentum of the accreted flow is then sufficient to form a persistent disc-like structure. This mechanism could explain the proposed limited outer extension of the accretion disc in Cygnus X-1 and suggests the presence of a disc at the outer rim of the neutron star magnetosphere in Vela X-1, with dramatic consequences on the spinning up of the accretor.

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms., Comment: 38 pages, 12 figures, accepted by ApJ

Supergiant X-ray Binaries (SgXB) host a compact object, often a neutron star (NS), orbiting an evolved O/B star. Mass transfer proceeds through the intense line-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the NS. The subsequent accretion process onto the NS is responsible for the abundant X-ray emission from SgXB. They also display peak-to-peak variability of the X-ray flux by a factor of a few 10 to 100, along with changes in the hardness ratios possibly due to varying absorption along the line-of-sight. We use recent radiation-hydrodynamic simulations of inhomogeneities (aka clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. For this, we run 3D hydrodynamic simulations of the wind in the vicinity of the accretor to investigate the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the NS magnetosphere. In particular, we show that the impact of the wind clumps on the time-variability of the intrinsic mass accretion rate is severely tempered by the crossing of the shock, compared to the purely ballistic Bondi-Hoyle-Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line-of-sight and estimate the final effective variability of the column density and mass accretion rate for different orbital separations. Finally, we compare our results to the most recent analysis of the X-ray flux and the hardness ratio in Vela X-1.

We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensional and three-dimensional examples of spherically symmetric Bondi accretion, steady planar Bondi-Hoyle-Lyttleton flows, and wind accretion in Supergiant X-ray binaries. Another improvement is support for the generic splitting of any background magnetic field. We present several tests relevant for solar physics applications to demonstrate the advantages of field splitting on accuracy and robustness in extremely low plasma $\beta$ environments: a static magnetic flux rope, a magnetic null-point, and magnetic reconnection in a current sheet with either uniform or anomalous resistivity. Our implementation for treating anisotropic thermal conduction in multi-dimensional MHD applications is also described, which generalizes the original slope limited symmetric scheme from 2D to 3D. We perform ring diffusion tests that demonstrate its accuracy and robustness, and show that it prevents the unphysical thermal flux present in traditional schemes. The improved parallel scaling of the code is demonstrated with 3D AMR simulations of solar coronal rain, which show satisfactory strong scaling up to 2000 cores. Other framework improvements are also reported: the modernization and reorganization into a library, the handling of automatic regression tests, the use of inline/online Doxygen documentation, and a new future-proof data format for input/output, Comment: 33 pages, 18 figures

We analyze particle acceleration in explosive reconnection events in magnetically dominated proton-electron plasmas. Reconnection is driven by large-scale magnetic stresses in interacting current-carrying flux tubes. Our model relies on development of current-driven instabilities on macroscopic scales. These tilt-kink instabilities develop in an initially force-free equilibrium of repelling current channels. Using MHD methods we study a 3D model of repelling and interacting flux tubes in which we simultaneously evolve test particles, guided by electromagnetic fields obtained from MHD. We identify two stages of particle acceleration; Initially particles accelerate in the current channels, after which the flux ropes start tilting and kinking and particles accelerate due to reconnection processes in the plasma. The explosive stage of reconnection produces non-thermal energy distributions with slopes that depend on plasma resistivity and the initial particle velocity. We also discuss the influence of the length of the flux ropes on particle acceleration and energy distributions. This study extends previous 2.5D results to 3D setups, providing all ingredients needed to model realistic scenarios like solar flares, black hole flares and particle acceleration in pulsar wind nebulae: formation of strong resistive electric fields, explosive reconnection and non-thermal particle distributions. By assuming initial energy equipartition between electrons and protons, applying low resistivity in accordance with solar corona conditions and limiting the flux rope length to a fraction of a solar radius we obtain realistic energy distributions for solar flares with non-thermal power law tails and maximum electron energies up to 11 MeV and maximum proton energies up to 1 GeV., Comment: Matches published version in MNRAS, including Erratum

We investigate how night side cooling and surface friction impact surface temperatures and large scale circulation for tidally locked Earth-like planets. For each scenario, we vary the orbital period between $P_{rot}=1-100$~days and capture changes in climate states. We find drastic changes in climate states for different surface friction scenarios. For very efficient surface friction ($t_{s,fric}=$ 0.1 days), the simulations for short rotation periods ($P_{rot} \leq$ 10 days) show predominantly standing extra tropical Rossby waves. These waves lead to climate states with two high latitude westerly jets and unperturbed meridional direct circulation. In most other scenarios, simulations with short rotation periods exhibit instead dominance by standing tropical Rossby waves. Such climate states have a single equatorial westerly jet, which disrupts direct circulation. Experiments with weak surface friction ($t_{s,fric}=~10 -100$ days) show decoupling between surface temperatures and circulation, which leads to strong cooling of the night side. The experiment with $t_{s,fric}= 100$ days assumes climate states with easterly flow (retrograde rotation) for medium and slow planetary rotations $P_{rot}= 12 - 100$~days. We show that an increase of night side cooling efficiency by one order of magnitude compared to the nominal model leads to a cooling of the night side surface temperatures by 80-100~K. The day side surface temperatures only drop by 25~K at the same time. The increase in thermal forcing suppresses the formation of extra tropical Rossby waves on small planets ($R_P=1 R_{Earth}$) in the short rotation period regime ($P_{rot} \leq$ 10 days)., Comment: 25 pages, 21 figures, accepted by MNRAS

Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized \textit{Solar Dynamics Observatory}/Atmospheric Imaging Assembly (AIA) 131, 94~\AA~emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the \textit{Solar and Heliospheric Observatory}/Solar Ultraviolet Measurement of Emitted Radiation (SUMER) Fe XIX line emission, we confirm that these reflected slow mode waves are propagating waves., Comment: 10 pages, 5 figures

We extend our earlier multidimensional, magnetohydrodynamic simulations of coronal rain occurring in magnetic arcades with higher resolution, grid-adaptive computations covering a much longer ($>6$ hour) timespan. We quantify how in-situ forming blob-like condensations grow along and across field lines and show that rain showers can occur in limit cycles, here demonstrated for the first time in 2.5D setups. We discuss dynamical, multi-dimensional aspects of the rebound shocks generated by the siphon inflows and quantify the thermodynamics of a prominence-corona-transition-region like structure surrounding the blobs. We point out the correlation between condensation rates and the cross-sectional size of loop systems where catastrophic cooling takes place. We also study the variations of the typical number density, kinetic energy and temperature while blobs descend, impact and sink into the transition region. In addition, we explain the mechanisms leading to concurrent upflows while the blobs descend. As a result, there are plenty of shear flows generated with relative velocity difference around 80 km s$^{-1}$ in our simulations. These shear flows are siphon flows set up by multiple blob dynamics and they in turn affect the deformation of the falling blobs. In particular, we show how shear flows can break apart blobs into smaller fragments, within minutes.

We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips., Comment: 47 pages, 59 figures

In light of the recent detection of direct evidence for the formation of Kelvin-Helmholtz instabilities in the Orion nebula, we expand upon previous modelling efforts by numerically simulating the shear-flow driven gas and dust dynamics in locations where the H$_{II}$ region and the molecular cloud interact. We aim to directly confront the simulation results with the infrared observations. Methods: To numerically model the onset and full nonlinear development of the Kelvin-Helmholtz instability we take the setup proposed to interpret the observations, and adjust it to a full 3D hydrodynamical simulation that includes the dynamics of gas as well as dust. A dust grain distribution with sizes between 5-250 nm is used, exploiting the gas+dust module of the MPI-AMRVAC code, in which the dust species are represented by several pressureless dust fluids. The evolution of the model is followed well into the nonlinear phase. The output of these simulations is then used as input for the SKIRT dust radiative transfer code to obtain infrared images at several stages of the evolution, which can be compared to the observations. Results: We confirm that a 3D Kelvin-Helmholtz instability is able to develop in the proposed setup, and that the formation of the instability is not inhibited by the addition of dust. Kelvin-Helmholtz billows form at the end of the linear phase, and synthetic observations of the billows show striking similarities to the infrared observations. It is pointed out that the high density dust regions preferentially collect on the flanks of the billows. To get agreement with the observed Kelvin-Helmholtz ripples, the assumed geometry between the background radiation, the billows and the observer is seen to be of critical importance., Comment: 8 pages, 10 figures

Magnetohydrodynamic waves interact with structured plasmas and reveal the internal magnetic and thermal structures therein, thereby having seismological applications in the solar atmosphere. We investigate the evolution of fast magnetoacoustic pulses in randomly structured plasmas, in the context of large-scale propagating waves in the solar atmosphere. We perform one dimensional numerical simulations of fast wave pulses propagating perpendicular to a constant magnetic field in a low-$\beta$ plasma with a random density profile across the field. Both linear and nonlinear regimes are considered. We study how the evolution of the pulse amplitude and width depends on their initial values and the parameters of the random structuring. A randomly structured plasma acts as a dispersive medium for a fast magnetoacoustic pulse, causing amplitude attenuation and broadening of the pulse width. After the passage of the main pulse, secondary propagating and standing fast waves appear in the plasma. Width evolution of both linear and nonlinear pulses can be well approximated by linear functions; however, narrow pulses may have zero or negative broadening. This arises because a narrow pulse is prone to splitting, while a broad pulse usually deviates less from their initial Gaussian shape and form ripple structures on top of the main pulse. A linear pulse decays at almost a constant rate, while a nonlinear pulse decays exponentially. A pulse interacts most efficiently with a random medium which has a correlation length of about half of its initial pulse width. The development of a detailed model of a fast MHD pulse propagating in highly structured medium substantiates the interpretation of EIT waves as fast magnetoacoustic waves. Evolution of a fast pulse provides us with a novel method to diagnose the sub-resolution filamentation of the solar atmosphere., Comment: 25 pages, 8 figures

In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure., Comment: Accepted for publication in ApJS

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

We present the first multidimensional, magnetohydrodynamic simulations which capture the initial formation and the long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in-situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match with modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into $V$-shaped like features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views on blobs which evaporate in situ, or get siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings, to connect parametrized heating prescriptions with rain statistics, ultimately allowing to quantify the coronal heating input., Comment: 12 pages, 4 figures, accepted by ApJ Letters at 2013-05-14

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., Comment: 27 pages, 9 figures

We present in this paper both a linear study and numerical relativistic MHD simulations of the non-resonant streaming instability occurring in the precursor of relativistic shocks. In the shock front restframe, we perform a linear analysis of this instability in a likely configuration for ultra-relativistic shock precursors. This considers magneto-acoustic waves having a wave vector perpendicular to the shock front and the large scale magnetic field. Our linear analysis is achieved without any assumption on the shock velocity and is thus valid for all velocity regimes. In order to check our calculation, we also perform relativistic MHD simulations describing the propagation of the aforementioned magneto-acoustic waves through the shock precursor. The numerical calculations confirm our linear analysis, which predicts that the growth rate of the instability is maximal for ultra-relativistic shocks and exhibits a wavenumber dependence $\propto k_x^{1/2}$. Our numerical simulations also depict the saturation regime of the instability where we show that the magnetic amplification is moderate but nevertheless significant ($\delta B/B\leq 10$). This latter fact may explain the presence of strong turbulence in the vicinity of relativistic magnetized shocks. Our numerical approach also introduces a convenient means to handle isothermal (ultra-)relativistic MHD conditions., Comment: 14 pages, 6 figures, MNRAS (in press)

It is found that both microflare-sized impulsive heating at one leg of the loop and a suddenly imposed velocity perturbation can propel the prominence to oscillate along the magnetic dip. An extensive parameter survey results in a scaling law, showing that the period of the oscillation, which weakly depends on the length and height of the prominence, and the amplitude of the perturbations, scales with $\sqrt{R/g_\odot}$, where $R$ represents the curvature radius of the dip, and $g_\odot$ is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendulum, which implies that the field-aligned component of gravity is the main restoring force for the prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pressure gradient becomes non-negligible for short prominences. The oscillation damps with time in the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., $\tau\sim l^{1.63} D^{0.66}w^{-1.21}v_{0}^{-0.30}$, showing strong dependence on the prominence length $l$, the geometry of the magnetic dip (characterized by the depth $D$ and the width $w$), and the velocity perturbation amplitude $v_0$. The larger the amplitude, the faster the oscillation damps. It is also found that mass drainage significantly reduces the damping timescale when the perturbation is too strong., Comment: 17 PAGES, 8FIGURES

In the context of planet formation, anticyclonic vortices have recently received lots of attention for the role they can play in planetesimals formation. Radial migration of intermediate size solids toward the central star may prevent their growth to larger solid grains. On the other hand, vortices can trap the dust and accelerate this growth, counteracting fast radial transport. Multiple effects have been shown to affect this scenario, such as vortex migration or decay. The aim of this paper is to study the formation of vortices by the Rossby wave instability and their long term evolution in a full three dimensional protoplanetary disc. We use a robust numerical scheme combined with adaptive mesh refinement in cylindrical coordinates, allowing to affordably compute long term 3D evolutions. We consider a full disc stratified both radially and vertically that is prone to formation of vortices by the Rossby wave instability. We show that the 3D Rossby vortices grow and survive over hundreds of years without migration. The localized overdensity which initiated the instability and vortex formation survives the growth of the Rossby wave instability for very long times. When the vortices are no longer sustained by the Rossby wave instability, their shape changes toward more elliptical vortices. This allows them to survive shear-driven destruction, but they may be prone to elliptical instability and slow decay. When the conditions for growing Rossby wave-related instabilities are maintained in the disc, large-scale vortices can survive over very long timescales and may be able to concentrate solids., Comment: Accepted for publication in A&A

Starting from accurate MHD flux rope equilibria containing prominence condensations, we initiate a systematic survey of their linear eigenoscillations. To quantify the full spectrum of linear MHD eigenmodes, we require knowledge of all flux-surface localized modes, charting out the continuous parts of the MHD spectrum. We combine analytical and numerical findings for the continuous spectrum for realistic prominence configurations. The equations governing all eigenmodes for translationally symmetric, gravitating equilibria containing an axial shear flow, are analyzed, along with their flux-surface localized limit. The analysis is valid for general 2.5D equilibria, where either density, entropy, or temperature vary from one flux surface to another. We analyze the mode couplings caused by the poloidal variation in the flux rope equilibria, by performing a small gravity parameter expansion. We contrast the analytical results with continuous spectra obtained numerically. For equilibria where the density is a flux function, we show that continuum modes can be overstable, and we present the stability criterion for these convective continuum instabilities. Furthermore, for all equilibria, a four-mode coupling scheme between an Alfvenic mode of poloidal mode number m and three neighboring (m-1, m, m+1) slow modes is identified, occurring in the vicinity of rational flux surfaces. For realistically prominence equilibria, this coupling is shown to play an important role, from weak to stronger gravity parameter g values. The analytic predictions for small g are compared with numerical spectra, and progressive deviations for larger g are identified. The unstable continuum modes could be relevant for short-lived prominence configurations. The gaps created by poloidal mode coupling in the continuous spectrum need further analysis, as they form preferred frequency ranges for global eigenoscillations., Comment: Accepted by Astronmy & Astrophysics, 21 pages, 15 figures

Context. Prominence seismology exploits our knowledge of the linear eigenoscillations for representative magnetohydro- dynamic models of filaments. To date, highly idealized models for prominences have been used, especially with respect to the overall magnetic configurations. Aims. We initiate a more systematic survey of filament wave modes, where we consider full multi-dimensional models with twisted magnetic fields representative of the surrounding magnetic flux rope. This requires the ability to compute accurate 2.5 dimensional magnetohydrodynamic equilibria that balance Lorentz forces, gravity, and pressure gradients, while containing density enhancements (static or in motion). Methods. The governing extended Grad-Shafranov equation is discussed, along with an analytic prediction for circular flux ropes for the Shafranov shift of the central magnetic axis due to gravity. Numerical equilibria are computed with a finite element-based code, demonstrating fourth order accuracy on an explicitly known, non-trivial test case. Results. The code is then used to construct more realistic prominence equilibria, for all three possible choices of a free flux-function. We quantify the influence of gravity, and generate cool condensations in hot cavities, as well as multi- layered prominences. Conclusions. The internal flux rope equilibria computed here have the prerequisite numerical accuracy to allow a yet more advanced analysis of the complete spectrum of linear magnetohydrodynamic perturbations, as will be demonstrated in the companion paper., Comment: Accepted by Astronomy & Astrophysics, 15 pages, 15 figures

It has been established that cold plasma condensations can form in a magnetic loop subject to localized heating of the footpoints. In this paper, we use grid-adaptive numerical simulations of the radiative hydrodynamic equations to parametrically investigate the filament formation process in a pre-shaped loop with both steady and finite-time chromospheric heating. Compared to previous works, we consider low-lying loops with shallow dips, and use a more realistic description for the radiative losses. We demonstrate for the first time that the onset of thermal instability satisfies the linear instability criterion. The onset time of the condensation is roughly \sim 2 hr or more after the localized heating at the footpoint is effective, and the growth rate of the thread length varies from 800 km hr-1 to 4000 km hr-1, depending on the amplitude and the decay length scale characterizing this localized chromospheric heating. We show how single or multiple condensation segments may form in the coronal portion. In the asymmetric heating case, when two segments form, they approach and coalesce, and the coalesced condensation later drains down into the chromosphere. With a steady heating, this process repeats with a periodicity of several hours. While our parametric survey confirms and augments earlier findings, we also point out that steady heating is not necessary to sustain the condensation. Once the condensation is formed, it can keep growing also when the localized heating ceases. Finally, we show that the condensation can survive continuous buffeting by perturbations resulting from the photospheric p-mode waves., Comment: 43 pages, 18 figures

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., Comment: Accepted for publication in MNRAS

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., Comment: accepted for publication in A&A

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., Comment: MNRAS submitted. 9 pages, 12 figures

We simulate time-dependent particle acceleration in the blast wave of a young supernova remnant (SNR), using a Monte Carlo approach for the diffusion and acceleration of the particles, coupled to an MHD code. We calculate the distribution function of the cosmic rays concurrently with the hydrodynamic evolution of the SNR, and compare the results with those obtained using simple steady-state models. The surrounding medium into which the supernova remnant evolves turns out to be of great influence on the maximum energy to which particles are accelerated. In particular, a shock going through a $\rho \propto r^{-2}$ density profile causes acceleration to typically much higher energies than a shock going through a medium with a homogeneous density profile. We find systematic differences between steady-state analytical models and our time-dependent calculation in terms of spectral slope, maximum energy, and the shape of the cut-off of the particle spectrum at the highest energies. We also find that, provided that the magnetic field at the reverse shock is sufficiently strong to confine particles, cosmic rays can be easily re-accelerated at the reverse shock., Comment: 19 pages, 20 figures, accepted for publication in MNRAS

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., Comment: 18 pages, ApJ accepted

This work presents a new plasma cooling curve that is calculated using the SPEX package. We compare our cooling rates to those in previous works, and implement the new cooling function in the grid-adaptive framework `AMRVAC'. Contributions to the cooling rate by the individual elements are given, to allow for the creation of cooling curves tailored to specific abundance requirements. In some situations, it is important to be able to include radiative losses in the hydrodynamics. The enhanced compression ratio can trigger instabilities (such as the Vishniac thin-shell instability) that would otherwise be absent. For gas with temperatures below 10,000 K, the cooling time becomes very long and does not affect the gas on the timescales that are generally of interest for hydrodynamical simulations of circumstellar plasmas. However, above this temperature, a significant fraction of the elements is ionised, and the cooling rate increases by a factor 1000 relative to lower temperature plasmas., Comment: 11 pages, 6 figures. Typos fixed to match version on A&A 'forthcoming' website. Tables in text format online available at http://www.phys.uu.nl/~schure/cooling

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., Comment: 18 pages, 19 figures. Replaced with MNRAS accepted version, minor changes only

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., Comment: 4 figures, submitted to MNRAS letters

Observations show that the magnetic field in young supernova remnants (SNRs) is significantly stronger than can be expected from the compression of the circumstellar medium (CSM) by a factor of four expected for strong blast waves. Additionally, the polarization is mainly radial, which is also contrary to expectation from compression of the CSM magnetic field. Cosmic rays (CRs) may help to explain these two observed features. They can increase the compression ratio to factors well over those of regular strong shocks by adding a relativistic plasma component to the pressure, and by draining the shock of energy when CRs escape from the region. The higher compression ratio will also allow for the contact discontinuity, which is subject to the Rayleigh-Taylor (R-T) instability, to reach much further out to the forward shock. This could create a preferred radial polarization of the magnetic field. With an adaptive mesh refinement MHD code (AMRVAC), we simulate the evolution of SNRs with three different configurations of the initial CSM magnetic field, and look at two different equations of state in order to look at the possible influence of a CR plasma component. The spectrum of CRs can be simulated using test particles, of which we also show some preliminary results that agree well with available analytical solutions., Comment: 7 pages, 7 figures, to appear in Advances in Space Research