BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. A. Aloy

M. A. Aloy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-phastro-ph.IMphysics.comp-phastro-ph.COastro-ph.SRmath.NAmath-phmath.MP

Catalog footprint

What is connected

36works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2021arXiv

Computational General Relativistic Force-Free Electrodynamics: I. Multi-Coordinate Implementation and Testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole shadows from the galactic center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as black holes and neutron stars. By this, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. We provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the Einstein Toolkit. The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of black hole magnetospheres.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Computational General Relativistic Force-Free Electrodynamics: II. Characterization of Numerical Diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe's most extreme flows of energy. The discretization of Maxwell's equations - needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling - introduces numerical diffusion. It acts as a source of dissipation independent of the system's physical constituents. Understanding the numerical diffusion of scientific codes is the key to classify their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics, by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field ($\mathbf{j}_\parallel$), in combination with the break-down of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g. the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of $\mathbf{j}_\parallel$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational-wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO-Virgo detector noise and the correctness of our analyses as applied to the resulting data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive black hole in such binaries is well approximated by models with no more than 1% of black holes more massive than $45\,M_\odot$, and a power law index of $α= {1.3}^{+1.4}_{-1.7}$ (90% credibility). We also show that BBHs are unlikely to be composed of black holes with large spins aligned to the orbital angular momentum. Modelling the evolution of the BBH merger rate with redshift, we show that it is at or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of $R = {53.2}^{+55.8}_{-28.2}$ Gpc$^{-3}$ yr$^{-1}$ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of black holes via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of black holes across cosmic time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentric orbits, rendering the search less efficient to eccentric binary systems. Here we present results of a search for binary black hole mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. The search uses minimal assumptions on the morphology of the transient gravitational waveform. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $χ_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Striped Blandford/Znajek jets from advection of small scale magnetic field

Black hole - accretion disc systems are the central engines of relativistic jets from stellar to galactic scales. We numerically quantify the unsteady outgoing Poynting flux through the horizon of a rapidly spinning black hole endowed with a rotating accretion disc. The disc supports small-scale, concentric, flux tubes with zero net magnetic flux. Our General Relativistic force-free electrodynamics simulations follow the accretion onto the black hole over several hundred dynamical timescales in 3D. For the case of counter-rotating accretion discs, the average process efficiency reaches up to $\left\langleε\right\rangle\approx 0.43$, compared to a stationary energy extraction by the Blandford/Znajek process. The process efficiency depends on the cross-sectional area of the loops, i.e. on the product $l\times h$, where $l$ is the radial loop thickness and $h$ its vertical scale height. We identify a strong correlation between efficient electromagnetic energy extraction and the quasi-stationary setting of ideal conditions for the operation of the Blandford/Znajek process (e.g. optimal field line angular velocity and fulfillment of the so-called Znajek condition). Remarkably, the energy extraction operates intermittently (alternating episodes of high and low efficiency) without imposing any large-scale magnetic field embedding the central object. Scaling our results to supermassive black holes, we estimate that the typical variability timescale of the system is of the order of days to months. Such timescales may account for the longest variability scales of TeV emission observed, e.g. in M87.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynamics of the progenitor and the properties of the outflow. We searched for associated emission of gravitational waves and high-energy neutrinos from astrophysical transients with minimal assumptions using data from Advanced LIGO from its first observing run O1, and data from the ANTARES and IceCube neutrino observatories from the same time period. We focused on candidate events whose astrophysical origin could not be determined from a single messenger. We found no significant coincident candidate, which we used to constrain the rate density of astrophysical sources dependent on their gravitational wave and neutrino emission processes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical simulations of the jetted tidal disruption event Swift J1644+57

In this work we focus on the technical details of the numerical simulations of the non- thermal transient Swift J1644+57, whose emission is probably produced by a two-component jet powered by a tidal disruption event. In this context we provide details of the coupling between the relativistic hydrodynamic simulations and the radiative transfer code. First, we consider the technical demands of one-dimensional simulations of a fast relativistic jet, and show to what extent (for the same physical parameters of the model) do the computed light curves depend on the numerical parameters of the different codes employed. In the second part we explain the difficulties of computing light curves from axisymmetric two dimensonal simulations and discuss a procedure that yields an acceptable tradeoff between the computational cost and the quality of the results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method

The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations ($Au=b$) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of $P$ relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of $M$ consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method. Our method to estimate the weights has the advantage that the explicit computation of the maximum and minimum eigenvalues of the matrix $A$ is replaced by the (much easier) calculation of the maximum and minimum frequencies derived from a von Neumann analysis. This set of weights is also optimal for the general problem, resulting in the fastest convergence of all possible SRJ schemes for a given grid structure. We also show that with the set of weights computed for the optimal SRJ scheme for a fixed cycle size it is possible to estimate numerically the optimal value of the parameter $ω$ in the Successive Overtaxation (SOR) method in some cases. Finally, we demonstrate with practical examples that our method also works very well for Poisson-like problems in which a high-order discretization of the Laplacian operator is employed. This is of interest since the former discretizations do not yield consistently ordered $A$ matrices. Furthermore, the optimal SRJ schemes here deduced, are advantageous over existing SOR implementations for high-order discretizations of the Laplacian operator in as much as they do not need to resort to multi-coloring schemes for their parallel implementation. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The influence of circumnuclear environment on the radio emission from TDE jets

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Scheduled Relaxation Jacobi method: improvements and applications

Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficiency in the reduction of the residual increases with the number of levels employed in the algorithm. Applying the original methodology to compute the algorithm parameters with more than 5 levels notably hinders obtaining optimal SRJ schemes, as the mixed (non-linear) algebraic-differential equations from which they result become notably stiff. Here we present a new methodology for obtaining the parameters of SRJ schemes that overcomes the limitations of the original algorithm and provide parameters for SRJ schemes with up to 15 levels and resolutions of up to $2^{15}$ points per dimension, allowing for acceleration factors larger than several hundreds with respect to the Jacobi method for typical resolutions and, in some high resolution cases, close to 1000. Furthermore, we extend the original algorithm to apply it to certain systems of non-linear ePDEs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The radio afterglow of Swift J1644+57 reveals a powerful jet with fast core and slow sheath

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accompanying synchrotron radio emission is produced by the shock interaction between the jet and the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, however, is the origin of the late-time radio rebrightening, which occurred well after the peak of the jetted X-ray emission. Here, we systematically explore several proposed explanations for this behavior by means of multi-dimensional hydrodynamic simulations coupled to a self-consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one-dimensional top-hat angular structure. However, a more complex angular structure comprised of an ultra-relativistic core (Lorentz factor $Γ\sim 10$) surrounded by a slower ($Γ\sim $ 2) sheath provides a reasonable fit to the data. Such a geometry could result from the radial structure of the super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of the ejecta that we infer of $\sim$ few $10^{53}\,$erg requires a highly efficient jet launching mechanism. Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is used to predict the radio light curves for off-axis viewing angles. Implications for the presence of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects for detecting orphan TDE afterglows with upcoming wide-field radio surveys and resolving the jet structure with long baseline interferometry, are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A method for computing synchrotron and inverse-Compton emission from hydrodynamic simulations of supernova remnants

The observational signature of supernova remnants (SNRs) is very complex, in terms of both their geometrical shape and their spectral properties, dominated by non-thermal synchrotron and inverse-Compton scattering. We propose a post-processing method to analyse the broad-band emission of SNRs based on three-dimensional hydrodynamical simulations. From the hydrodynamical data, we estimate the distribution of non-thermal electrons accelerated at the shock wave and follow the subsequent evolution as they lose or gain energy by adiabatic expansion or compression and emit energy by radiation. As a first test case, we use a simulation of a bipolar supernova expanding into a cloudy medium. We find that our method qualitatively reproduces the main observational features of typical SNRs and produces fluxes that agree with observations to within a factor of a few. allowing for further use in more extended sets of models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

General Relativistic Simulations of the Collapsar Scenario

We are exploring the viability of the collapsar model for long-soft gamma-ray bursts. For this we perform state-of-the-art general relativistic hydrodynamic simulations in a dynamically evolving space-time with the CoCoNuT code. We start from massive low metallicity stellar models evolved up to core gravitational instability, and then follow the subsequent evolution until the system collapses forming a compact remnant. A preliminary study of the collapse outcome is performed by varying the typical parameters of the scenario, such as the initial stellar mass, metallicity, and rotational profile of the stellar progenitor. 1D models (without rotation) have been used to test our newly developed neutrino leakage scheme. This is a fundamental piece of our approach as it allows the central remnant (in all cases considered, a metastable high-mass neutron star) to cool down, eventually collapsing to a black hole. In two dimensions, we show that sufficiently fast rotating cores lead to the formation of Kerr black holes, due to the fall-back of matter surrounding the compact remnant, which has not been successfully unbounded by a precedent supernova shock.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The influence of the magnetic field on the spectral properties of blazars

We explore the signature imprinted by dynamically relevant magnetic fields on the spectral energy distribution (SED) of blazars. It is assumed that the emission from these sources originates from the collision of cold plasma shells, whose magnetohydrodynamic evolution we compute by numerically solving Riemann problems. We compute the SEDs including the most relevant radiative processes and scan a broad parameter space that encompasses a significant fraction of the commonly accepted values of not directly measurable physical properties. We reproduce the standard double hump SED found in blazar observations for unmagnetized shells, but show that the prototype double hump structure of blazars can also be reproduced if the dynamical source of the radiation field is very ultrarelativistic both, in a kinematically sense (namely, if it has Lorentz factors \gtrsim 50) and regarding its magnetization (e.g., with flow magnetizations σ\simeq 0.1). A fair fraction of the {\em blazar sequence} could be explained as a consequence of shell magnetization: negligible magnetization in FSRQs, and moderate or large (and uniform) magnetization in BL Lacs. The predicted photon spectral indices (\gph) in the γ-ray band are above the observed values (Γ_{\rm ph, obs} \lesssim 2.6 for sources with redshifts 0.4\le z \le 0.6) if the magnetization of the sources is moderate (σ\simeq 10^{-2}).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative signature of magnetic fields in internal shocks

Common models of blazars and gamma-ray bursts assume that the plasma underlying the ob- served phenomenology is magnetized to some extent. Within this context, radiative signatures of dissipation of kinetic and conversion of magnetic energy in internal shocks of relativistic magnetized outflows are studied. We model internal shocks as being caused by collisions of homogeneous plasma shells. We compute the flow state after the shell interaction by solving Riemann problems at the contact surface between the colliding shells, and then compute the emission from the resulting shocks. Under the assumption of a constant flow luminosity we find that there is a clear difference between the models where both shells are weakly magne- tized (σ<\sim0.01) and those where, at least, one shell has a σ>\sim0.01. We obtain that the radiative efficiency is largest for models in which, regardless of the ordering, one shell is weakly and the other strongly magnetized. Substantial differences between weakly and strongly magne- tized shell collisions are observed in the inverse-Compton part of the spectrum, as well as in the optical, X-ray and 1GeV light curves. We propose a way to distinguish observationally between weakly magnetized from magnetized internal shocks by comparing the maximum frequency of the inverse-Compton and synchrotron part of the spectrum to the ratio of the inverse-Compton and synchrotron fluence. Finally, our results suggest that LBL blazars may correspond to barely magnetized flows, while HBL blazars could correspond to moderately magnetized ones. Indeed, by comparing with actual blazar observations we conclude that the magnetization of typical blazars is σ <\sim 0.01 for the internal shock model to be valid in these sources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An unusual stellar death on Christmas Day

Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- sociated with supernovae (Woosley et al. 2006). They release ultra-relativistic jets producing non-thermal emission through synchrotron radiation as they interact with the surrounding medium (Zhang et al. 2004). Here we report observations of the peculiar GRB 101225A (the "Christmas burst"). Its gamma-ray emission was exceptionally long and followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously ejected material and thus creating the observed black-body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy (Campana et al. 2011).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Efficiency of internal shocks in magnetized relativistic jets

We study the dynamic and radiative efficiency of conversion of kinetic-to-thermal/magnetic energy by internal shocks in relativistic magnetized outflows. A parameter study of a large number of collisions of cylindrical shells is performed. We explore how, while keeping the total flow luminosity constant, the variable fluid magnetization influences the efficiency and find that the interaction of shells in a mildly magnetized jet yields higher dynamic, but lower radiative efficiency than in a non-magnetized flow. A multi-wavelength radiative signature of different shell magnetization is computed assuming that relativistic particles are accelerated at internal shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Variable Ly alpha sheds light on the environment surrounding GRB 090426

Long duration gamma-ray bursts are commonly associated with the deaths of massive stars. Spectroscopic studies using the afterglow as a light source provide a unique opportunity to unveil the medium surrounding it, probing the densest region of their galaxies. This material is usually in a low ionisation state and at large distances from the burst site, hence representing the normal interstellar medium in the galaxy. Here we present the case of GRB 090426 at z=2.609, whose optical spectrum indicates an almost fully ionised medium together with a low column density of neutral hydrogen. For the first time, we also observe variations in the Ly alpha absorption line. Photoionisation modeling shows that we are probing material from the vicinity of the burst (~80 pc). The host galaxy is a complex of two luminous interacting galaxies, which might suggest that this burst could have occurred in an isolated star-forming region outside its host galaxy created in the interaction of the two galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Roe-type Riemann solver based on the spectral decomposition of the equations of Relativistic Magnetohydrodynamics

In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we have developed, and some note on the (non-)convex character of the relativistic MHD equations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical efficiency of collisionless magnetized shocks in relativistic jets

The so-called internal shock model aims to explain the light-curves and spectra produced by non-thermal processes originated in the flow of blazars and gamma-ray bursts. A long standing question is whether the tenuous collisionless shocks, driven inside a relativistic flow, are efficient enough to explain the amount of energy observed as compared with the expected kinetic power of the outflow. In this work we study the dynamic efficiency of conversion of kinetic-to- thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We find that the collision between shells with a non-zero relative velocity can yield either two oppositely moving shocks (in the frame where the contact surface is at rest), or a reverse shock and a forward rarefaction. For moderately magnetized shocks (magnetization σ ~ 0.1), the dynamic efficiency in a single two-shell interaction can be as large as 40%. Hence, the dynamic efficiency of moderately magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. We find that the efficiency is only weakly dependent on the Lorentz factor of the shells and, thus internal shocks in the magnetized flow of blazars and gamma-ray bursts are approximately equally efficient.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Local simulations of the magnetized Kelvin-Helmholtz instability in neutron-star mergers

Context. Global MHD simulations show Kelvin-Helmholtz (KH) instabilities at the contact surface of two merging neutron stars. That region has been identified as the site of efficient amplification of magnetic fields. However, these global simulations, due to numerical limitations, were unable to determine the saturation level of the field strength, and thus the possible back-reaction of the magnetic field onto the flow. Aims. We investigate the amplification of initially weak fields in KH unstable shear flows, and the back-reaction of the field onto the flow. Methods. We use a high-resolution ideal MHD code to perform 2D and 3D local simulations of shear flows. Results. In 2D, the magnetic field is amplified in less than 0.01ms until it reaches locally equipartition with the kinetic energy. Subsequently, it saturates due to resistive instabilities that disrupt the KH vortex and decelerate the shear flow on a secular time scale. We determine scaling laws of the field amplification with the initial field strength and the grid resolution. In 3D, this hydromagnetic mechanism may be dominated by purely hydrodynamic instabilities limiting the amplification. We find maximum magnetic fields of 10^16 G locally, and r.m.s. maxima within the box of 10^15 G. However, such strong fields exist only for a short period. In the saturated state, the magnetic field is mainly oriented parallel to the shear flow for strong initial fields, while weaker initial fields tend to lead to a more balanced distribution of the field energy. In all models the flow shows small-scale features. The magnetic field is at most in equipartition with the decaying shear flow. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiwavelength afterglow light curves from magnetized GRB flows

We use high-resolution relativistic MHD simulations coupled with a radiative transfer code to compute multiwavelength afterglow light curves of magnetized ejecta of gamma-ray bursts interacting with a uniform circumburst medium. The aim of our study is to determine how the magnetization of the ejecta at large distance from the central engine influences the afterglow emission, and to assess whether observations can be reliably used to infer the strength of the magnetic field. We find that, for typical parameters of the ejecta, the emission from the reverse shock peaks for magnetization $σ_0 \sim 0.01 - 0.1$ of the flow, and that it is greatly suppressed for higher $σ_0$. The emission from the forward shock shows an achromatic break shortly after the end of the burst marking the onset of the self-similar evolution of the blast wave. Fitting the early afterglow of GRB 990123 and 090102 with our numerical models we infer respective magnetizations of $σ_0 \sim 0.01$ and $σ_0 \sim 0.1$ for these bursts. We argue that the lack of observed reverse shock emission from the majority of the bursts can be understood if $σ_0 \simmore 0.1$, since we obtain that the luminosity of the reverse shock decreases significantly for $σ_0 \sim 1$. For ejecta with $σ_0 \simmore 0.1$ our models predict that there is sufficient energy left in the magnetic field, at least during an interval of ~10 times the burst duration, to produce a substantial emission if the magnetic energy can be dissipated (for instance, due to resistive effects) and radiated away.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Searching for differences in Swift's intermediate GRBs

Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Detection of the high z GRB 080913 and its implications on progenitors and energy extraction mechanisms

Aims: We present multiwavelength observations of one of the most distant gamma-ray bursts detected so far, GRB 080913. Based on these observations, we consider whether it could be classified as a short-duration GRB and discuss the implications for the progenitor nature and energy extraction mechanisms. Methods: Multiwavelength X-ray, near IR and millimetre observations were made between 20.7 hours and 16.8 days after the event. Results: Whereas a very faint afterglow was seen at the 3.5m CAHA telescope in the nIR, the X-ray afterglow was clearly detected in both Swift and XMM-Newton observations. An upper limit is reported in the mm range. We have modeled the data assuming a collimated $θ_0$ $\gtrsim$ 3$^\circ$ blast wave with an energy injection at 0.5 days carrying $5\sim 10^{52}$ erg or approximately 12 times the initial energy of the blast wave. We find that GRB 080913 shares many of the gamma-ray diagnostics with the more recent burst GRB 090423 for being classified as short had they ocurred at low redshift. If the progenitor were a compact binary merger, it is likely composed by a NS and BH. The Blandford-Znajek (BZ) mechanism is the preferred one to extract energy from the central, maximally-rotating BH. Both the magnetic field close to the event horizon (B) and the BH mass ($M_{bh}$) are restricted within a relatively narrow range, such that $(B / 3\times 10^{16} \rm{G}) (M_{bh} / 7 M_\odot) \sim 1$. Similar constraints on the central BH hold for collapsar-like progenitor systems if the BZ-mechanism works for the system at hand.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the dynamic efficiency of internal shocks in magnetized relativistic outflows

We study the dynamic efficiency of conversion of kinetic-to-thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We model internal shocks as being caused by collisions of shells of plasma with the same energy flux and a non-zero relative velocity. The contact surface, where the interaction between the shells takes place, can break up either into two oppositely moving shocks (in the frame where the contact surface is at rest), or into a reverse shock and a forward rarefaction. We find that for moderately magnetized shocks (magnetization $σ\simeq 0.1$), the dynamic efficiency in a single two-shell interaction can be as large as 40%. Thus, the dynamic efficiency of moderately magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. If the slower shell propagates with a sufficiently large velocity, the efficiency is only weakly dependent on its Lorentz factor. Consequently, the dynamic efficiency of shell interactions in the magnetized flow of blazars and gamma-ray bursts is effectively the same. These results are quantitatively rather independent on the equation of state of the plasma. The radiative efficiency of the process is expected to be a fraction $f_r<1$ of the estimated dynamic one, the exact value of $f_r$ depending on the particularities of the emission processes which radiate away the thermal or magnetic energy of the shocked states.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver

We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Spectral evolution of superluminal components in parsec-scale jets

(Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectral evolution, with observational consequences only above radio frequencies. Larger values of the magnetic field (\ab \sim 6\times 10^2), such that synchrotron losses are moderately important at radio frequencies, present a larger ratio of shocked-to-unshocked regions brightness than the models without radiative losses, despite the fact that they correspond to the same underlying hydrodynamic structure. We also show that jets with a positive photon spectral index result if the lower limit γ_min of the non-thermal particle energy distribution is large enough. A temporary increase of the Lorentz factor at the jet inlet produces a traveling perturbation that appears in the synthetic maps as a superluminal component. We show that trailing components can be originated not only in pressure matched jets, but also in over-pressured ones, where the existence of recollimation shocks does not allow for a direct identification of such features as Kelvin-Helmholtz modes, and its observational imprint depends on the observing frequency. If the magnetic field is large (\ab \sim 6\times 10^2), the spectral index in the rarefaction trailing the traveling perturbation does not change much with respect to the same model without any hydrodynamic perturbation. If the synchrotron losses are considered the spectral index displays a smaller value than in the corresponding region of the quiescent jet model.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Relativistic outflows from remnants of compact object mergers and their viability for short gamma-ray bursts

We present the first general relativistic hydrodynamic models of the launch and evolution of relativistic jets and winds, driven by thermal energy deposition, possibly due to neutrino-antineutrino annihilation, in the close vicinity of black hole-accretion torus systems. The latter are considered to be the remnants of compact object mergers. Our two-dimensional simulations establish the link between such mergers and future observations of short gamma-ray bursts (GRBs) by the SWIFT satellite. They show that ultrarelativistic outflow with maximum terminal Lorentz factors (Gamma) around 1000 develops for polar energy deposition rates above some 1e48 erg/s per steradian, provided the merger environment has a sufficiently low baryon density. Due to the collimation by the dense accretion torus the typical semi-opening angles of the Gamma > 100 cone are 5-10 degrees, corresponding to about 0.4-1.5% of the hemisphere and apparent isotropized energies (kinetic plus internal) up to ~1e51 erg. 10-30% of the deposited energy are transferred to the outflow with Gamma > 100. Our models confirm the viability of post-merger BH-torus systems as engines of short, hard GRBs and can explain the durations of all observed short GRBs, because different propagation velocities of the front and rear ends lead to a radial stretching of the ultrarelativistic fireball before transparency is reached. The ultrarelativistic flow reveals a highly non-uniform structure with Lorentz factor variations up to factors of a few, caused by the action of Kelvin-Helmholtz instabilities that originate at the fireball-torus interface (abbreviated).

0 citations0 reviewsNo code linkedOpen access