Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index

The role of the equation of state for a perfectly conducting, relativistic magnetized fluid is the main subject of this work. The ideal constant $\Gamma$-law equation of state, commonly adopted in a wide range of astrophysical applications, is compared with a more realistic equation of state that better approximates the single-specie relativistic gas. The paper focus on three different topics. First, the influence of a more realistic equation of state on the propagation of fast magneto-sonic shocks is investigated. This calls into question the validity of the constant $\Gamma$-law equation of state in problems where the temperature of the gas substantially changes across hydromagnetic waves. Second, we present a new inversion scheme to recover primitive variables (such as rest-mass density and pressure) from conservative ones that allows for a general equation of state and avoids catastrophic numerical cancellations in the non-relativistic and ultrarelativistic limits. Finally, selected numerical tests of astrophysical relevance (including magnetized accretion flows around Kerr black holes) are compared using different equations of state. Our main conclusion is that the choice of a realistic equation of state can considerably bear upon the solution when transitions from cold to hot gas (or viceversa) are present. Under these circumstances, a polytropic equation of state can significantly endanger the solution.Comment: 14 pages, 14 figure

Two-component jet simulations: II. Combining analytical disk and stellar MHD outflow solutions

Theoretical arguments along with observational data of YSO jets suggest the presence of two steady components: a disk wind type outflow needed to explain the observed high mass loss rates and a stellar wind type outflow probably accounting for the observed stellar spin down. Each component's contribution depends on the intrinsic physical properties of the YSO-disk system and its evolutionary stage. The main goal of this paper is to understand some of the basic features of the evolution, interaction and co-existence of the two jet components over a parameter space and when time variability is enforced. Having studied separately the numerical evolution of each type of the complementary disk and stellar analytical wind solutions in Paper I of this series, we proceed here to mix together the two models inside the computational box. The evolution in time is performed with the PLUTO code, investigating the dynamics of the two-component jets, the modifications each solution undergoes and the potential steady state reached.Comment: accepted for publication in A&

The Athena Astrophysical MHD Code in Cylindrical Geometry

A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athena's original developers and has been designed to alter the existing Cartesian-coordinates code as minimally and transparently as possible. The numerical equations in cylindrical coordinates are formulated to maintain consistency with constrained transport, a central feature of the Athena algorithm, while making use of previously implemented code modules such as the Riemann solvers. Angular-momentum transport, which is critical in astrophysical disk systems dominated by rotation, is treated carefully. We describe modifications for cylindrical coordinates of the higher-order spatial reconstruction and characteristic evolution steps as well as the finite-volume and constrained transport updates. Finally, we present a test suite of standard and novel problems in one-, two-, and three-dimensions designed to validate our algorithms and implementation and to be of use to other code developers. The code is suitable for use in a wide variety of astrophysical applications and is freely available for download on the web

TPCI: The PLUTO-CLOUDY Interface

We present an interface between the (magneto-) hydrodynamics code PLUTO and the plasma simulation and spectral synthesis code CLOUDY. By combining these codes, we constructed a new photoionization hydrodynamics solver: The PLUTO-CLOUDY Interface (TPCI), which is well suited to simulate photoevaporative flows under strong irradiation. The code includes the electromagnetic spectrum from X-rays to the radio range and solves the photoionization and chemical network of the 30 lightest elements. TPCI follows an iterative numerical scheme: First, the equilibrium state of the medium is solved for a given radiation field by CLOUDY, resulting in a net radiative heating or cooling. In the second step, the latter influences the (magneto-) hydrodynamic evolution calculated by PLUTO. Here, we validated the one-dimensional version of the code on the basis of four test problems: Photoevaporation of a cool hydrogen cloud, cooling of coronal plasma, formation of a Stroemgren sphere, and the evaporating atmosphere of a hot Jupiter. This combination of an equilibrium photoionization solver with a general MHD code provides an advanced simulation tool applicable to a variety of astrophysical problems.Comment: 13 pages, 10 figures, accepted for publication in A&

Young stellar object jet models: From theory to synthetic observations

Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.Comment: accepted for publication A&A, 20 pages, 11 figure

TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.Comment: ApJS, 197, 1

The Dynamics of Radiative Shock Waves: Linear and Nonlinear Evolution

The stability properties of one-dimensional radiative shocks with a power-law cooling function of the form $\Lambda \propto \rho^2T^\alpha$ are the main subject of this work. The linear analysis originally presented by Chevalier & Imamura, is thoroughfully reviewed for several values of the cooling index $\alpha$ and higher overtone modes. Consistently with previous results, it is shown that the spectrum of the linear operator consists in a series of modes with increasing oscillation frequency. For each mode a critical value of the cooling index, $\alpha_\textrm{c}$, can be defined so that modes with $\alpha < \alpha_\textrm{c}$ are unstable, while modes with $\alpha > \alpha_\textrm{c}$ are stable. The perturbative analysis is complemented by several numerical simulations to follow the time-dependent evolution of the system for different values of $\alpha$. Particular attention is given to the comparison between numerical and analytical results (during the early phases of the evolution) and to the role played by different boundary conditions. It is shown that an appropriate treatment of the lower boundary yields results that closely follow the predicted linear behavior. During the nonlinear regime, the shock oscillations saturate at a finite amplitude and tend to a quasi-periodic cycle. The modes of oscillations during this phase do not necessarily coincide with those predicted by linear theory, but may be accounted for by mode-mode coupling.Comment: 33 pages, 12 figures, accepted for publication on the Astrophysical Journa

Model-independent determination of the cosmic expansion rate. I. Application to type-Ia supernovae

Aims: In view of the substantial uncertainties regarding the possible dynamics of the dark energy, we aim at constraining the expansion rate of the universe without reference to a specific Friedmann model and its parameters. Methods: We show that cosmological observables integrating over the cosmic expansion rate can be converted into a Volterra integral equation which is known to have a unique solution in terms of a Neumann series. Expanding observables such as the luminosity distances to type-Ia supernovae into a series of orthonormal functions, the integral equation can be solved and the cosmic expansion rate recovered within the limits allowed by the accuracy of the data. Results: We demonstrate the performance of the method applying it to synthetic data sets of increasing complexity, and to the first-year SNLS data. In particular, we show that the method is capable of reproducing a hypothetical expansion function containing a sudden transition.Comment: 9 pages, 8 figures; accepted by A&A; subsection 3.6 added, new references and minor change

A model of Mira's cometary head/tail entering the Local Bubble

We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.Comment: 12 pages, 8 figures (4 in color). Accepted for publication in The Astrophysical Journa

Large scale magnetic fields in viscous resistive accretion disks. I. Ejection from weakly magnetized disks

Cold steady-state disk wind theory from near Keplerian accretion disks requires a large scale magnetic field at near equipartition strength. However the minimum magnetization has never been tested. We investigate the time evolution of an accretion disk threaded by a weak vertical magnetic field. The strength of the field is such that the disk magnetization falls off rapidly with radius. Four 2.5D numerical simulations of viscous resistive accretion disk are performed using the magnetohydrodynamic code PLUTO. In these simulations, a mean field approach is used and turbulence is assumed to give rise to anomalous transport coefficients (alpha prescription). The large scale magnetic field introduces only a small perturbation to the disk structure, with accretion driven by the dominant viscous torque. A super fast magnetosonic jet is observed to be launched from the innermost regions and remains stationary over more than 953 Keplerian orbits. The self-confined jet is launched from a finite radial zone in the disk which remains constant over time. Ejection is made possible because the magnetization reaches unity at the disk surface, due to the steep density decrease. However, no ejection is reported when the midplane magnetization becomes too small. The asymptotic jet velocity remains nevertheless too low to explain observed jets due to the negligible power carried away by the jet. Astrophysical disks with superheated surface layers could drive analogous outflows even if their midplane magnetization is low. Sufficient angular momentum would be extracted by the turbulent viscosity to allow the accretion process to continue. The magnetized outflows would be no more than byproducts, rather than a fundamental driver of accretion. However, if the midplane magnetization increases towards the center, a natural transition to an inner jet dominated disk could be achieved.Comment: Accepted by Astronomy and Astrophysic