Source author record

Federico Carrasco

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

Gamma-radiation sky maps from compact binaries

We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.

preprint2022arXiv

Skylight: a new code for general-relativistic ray-tracing and radiative transfer in arbitrary spacetimes

To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from astrophysical models of compact objects as seen by distant observers. We incorporate two different schemes, namely Monte Carlo radiative transfer, integrating geodesics from the astrophysical region to distant observers, and camera techniques with backwards integration from the observer to the emission region. The code is validated by successfully passing several test cases, among them: thin accretion disks and neutron star hot spot emission.

preprint2020arXiv

Gradient sub-grid-scale model for relativistic MHD Large Eddy Simulations

MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations.

preprint2020arXiv

GRMHD large eddy simulations with gradient subgrid-scale model

The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations with a self-consistent subgrid-scale gradient model, that we generalized to the special relativistic case in a previous work and now extend to the general relativistic case. Adapted from nonrelativistic scenarios, the so-called gradient model allows to capture part of the effects of the hidden dynamics on the resolved scales, by means of a physically-agnostic, mathematically-based Taylor expansion of the nonlinear terms in the conservative evolution equations' fluxes. We assess the validity of this approach in bounding-box simulations of the magnetic Kelvin-Helmholtz instability. Several resolutions and a broad range of scenarios are considered in order to carefully test the performance of the model under three crucial aspects: (i) highly curved backgrounds, (ii) jumps on the fluid density profiles and (iii) strong shocks. The results suggest our extension of the gradient subgrid-scale model to general relativistic magnetohydrodynamics is a promising approach for studying binary neutron stars mergers, and potentially to other relevant astrophysical scenarios.

preprint2020arXiv

Magnetosphere of an orbiting neutron star

We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed within its surrounding magnetosphere. A detailed study of the resulting magnetospheric configurations in our setting is presented, incorporating as well the effects due to neutron star spin and the misalignment of the magnetic and orbital axes. We find many features resembling those of pulsar magnetospheres for the orbiting neutron star --even in the absence of spin--, being of particular interest the development of a spiral current sheet that extends beyond the light cylinder. Then, we use recent advances in pulsar theory to estimate electromagnetic emissions produced at the reconnection regions of such current sheets.

preprint2016arXiv

Covariant Hyperbolization of Force-free Electrodynamics

Force-Free Flectrodynamics (FFE) is a non-linear system of equations modeling the evolution of the electromagnetic field, in the presence of a magnetically dominated relativistic plasma. This configuration arises on several astrophysical scenarios, which represent exciting laboratories to understand physics in extreme regimes. We show that this system, when restricted to the correct constraint submanifold, is symmetric hyperbolic. In numerical applications is not feasible to keep the system in that submanifold, and so, it is necessary to analyze its structure first in the tangent space of that submanifold and then in a whole neighborhood of it. As already shown by Pfeiffer, a direct (or naive) formulation of this system (in the whole tangent space) results in a weakly hyperbolic system of evolution equations for which well-possednes for the initial value formulation does not follows. Using the generalized symmetric hyperbolic formalism due to Geroch, we introduce here a covariant hyperbolization for the FFE system. In fact, in analogy to the usual Maxwell case, a complete family of hyperbolizers is found, both for the restricted system on the constraint submanifold as well as for a suitably extended system defined in a whole neighborhood of it. A particular symmetrizer among the family is then used to write down the pertaining evolution equations, in a generic ($3+1$)-decomposition on a background spacetime. Interestingly, it turns out that for a particular choice of the lapse and shift functions of the foliation, our symmetrized system reduces to the one found by Pfeiffer. Finally we analyze the characteristic structure of the resulting evolution system.

preprint2015arXiv

Nonlinear electrodynamics as a symmetric hyperbolic system

Nonlinear theories generalizing Maxwell's electromagnetism and arising from a Lagrangian formalism have dispersion relations in which propagation planes factor into null planes corresponding to two effective metrics which depend on the point-wise values of the electromagnetic field. These effective Lorentzian metrics share the null (generically two) directions of the electromagnetic field. We show that, the theory is symmetric hyperbolic if and only if the cones these metrics give rise to have a non-empty intersection. Namely that there exist families of symmetrizers in the sense of Geroch which are positive definite for all covectors in the interior of the cones intersection. Thus, for these theories, the initial value problem is well-posed. We illustrate the power of this approach with several nonlinear models of physical interest such as Born-Infeld, Gauss-Bonnet and Euler-Heisenberg.

preprint2013arXiv

A Holographic Path to the Turbulent Side of Gravity

We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.

preprint2012arXiv

Turbulent flows for relativistic conformal fluids in 2+1 dimensions

We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbulent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. Our motivation for this study is to understand the implications for gravitational solutions through the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the corresponding black brane/black hole spacetimes (for sufficiently large scales/temperatures) will display a similar cascade towards longer wavelengths.