Source author record

C. Palenzuela

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

ResearcherUnclaimed source record

gr-qc astro-ph.HE hep-th math-ph math.MP

Catalog footprint

What is connected

9works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Large Eddy Simulations of Magnetized Mergers of Neutron Stars with Neutrinos

Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is increasingly important in light of advancements in multi-messenger astronomy generally, and in gravitational wave astronomy in particular (such as the development of third-generation detectors). Several general relativistic codes have been incorporating some of these elements with different levels of realism. Here, we extend our code MHDuet, which can perform large eddy simulations of magnetohydrodynamics to help capture the magnetic field amplification during the merger, and to allow for realistic equations of state and neutrino cooling via a leakage scheme. We perform several tests involving isolated and binary neutron stars demonstrating the accuracy of the code.

preprint2022arXiv

Turbulent magnetic field amplification in binary neutron star mergers

Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.

preprint2020arXiv

Towards fidelity and scalability in non-vacuum mergers

We study the evolution of two fiducial configurations for binary neutron stars using two different general relativistic hydrodynamics(GRHD), distributed adaptive mesh codes. One code, HAD, has for many years been used to study mergers of compact object binaries, while a new code, MHDuet, has been recently developed with the experience gained with the older one as well as several novel features for scalability improvements. As such, we examine the performance of each, placing particular focus on future requirements for the extraction of gravitational wave signatures of non-vacuum binaries.

preprint2015arXiv

Effects of the microphysical Equation of State in the mergers of magnetized Neutron Stars With Neutrino Cooling

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics.

preprint2010arXiv

Action principle for Numerical Relativity evolution systems

A Lagrangian density is provided, that allows to recover the Z4 evolution system from an action principle. The resulting system is then strongly hyperbolic when supplemented by gauge conditions like '1+log' or 'freezing shift', suitable for numerical evolution. The physical constraint $Z_μ= 0$ can be imposed just on the initial data. The corresponding canonical equations are also provided. This opens the door to analogous results for other numerical-relativity formalisms, like BSSN, that can be derived from Z4 by a symmetry-breaking procedure. The harmonic formulation can be easily recovered by a slight modification of the procedure. This provides a mechanism for deriving both the field evolution equations and the gauge conditions from the action principle, with a view on using symplectic integrators for a constraint-preserving numerical evolution. The gauge sources corresponding to the 'puncture gauge' conditions are identified in this context.

preprint2004arXiv

Dynamical shift conditions for the Z4 and BSSN hyperbolic formalisms

A class of dynamical shift conditions is shown to lead to a strongly hyperbolic evolution system, both in the Z4 and in the BSSN Numerical Relativity formalisms. This class generalizes the harmonic shift condition, where light speed is the only non-trivial characteristic speed, and it is contained into the multi-parameter family of minimal distortion shift conditions recently proposed by Lindblom and Scheel. The relationship with the analogous 'dynamical freezing' shift conditions used in black hole simulations discussed.

preprint2003arXiv

General-covariant evolution formalism for Numerical Relativity

A general covariant extension of Einstein\'{}s field equations is considered with a view to Numerical Relativity applications. The basic variables are taken to be the metric tensor and an additional four-vector $Z_μ$. Einstein's solutions are recovered when the additional four-vector vanishes, so that the energy and momentum constraints amount to the covariant algebraic condition $Z_μ=0$. The extended field equations can be supplemented by suitable coordinate conditions in order to provide symmetric hyperbolic evolution systems: this is actually the case for either harmonic coordinates or normal coordinates with harmonic slicing.

preprint2002arXiv

A 3+1 covariant suite of Numerical Relativity Evolution Systems

A suite of three evolution systems is presented in the framework of the 3+1 formalism. The first one is of second order in space derivatives and has the same causal structure of the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) system for a suitable choice of parameters. The second one is the standard first order version of the first one and has the same causal structure of the Bona-Masso system for a given parameter choice. The third one is obtained from the second one by reducing the space of variables in such a way that the only modes that propagate with zero characteristic speed are the trivial ones. This last system has the same structure of the ones recently presented by Kidder, Scheel and Teukolski: the correspondence between both sets of parameters is explicitly given. The fact that the suite started with a system in which all the dynamical variables behave as tensors (contrary to what happens with BSSN system) allows one to keep the same parametrization when passing from one system to the next in the suite. The direct relationship between each parameter and a particular characteristic speed, which is quite evident in the second and the third systems, is a direct consequence of the manifest 3+1 covariance of the approach.