Researcher profile

Elias R. Most

Elias R. Most contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 published item(s)

preprint2023arXiv

Causal, stable first-order viscous relativistic hydrodynamics with ideal gas microphysics

We present the first numerical analysis of causal, stable first-order relativistic hydrodynamics with ideal gas microphysics, based in the formalism developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK theory). The BDNK approach provides definitions for the conserved stress-energy tensor and baryon current, and rigorously proves causality, local well-posedness, strong hyperbolicity, and linear stability (about equilibrium) for the equations of motion, subject to a set of coupled nonlinear inequalities involving the undetermined model coefficients (the choice for which defines the "hydrodynamic frame"). We present a class of hydrodynamic frames derived from the relativistic ideal gas "gamma-law" equation of state which satisfy the BDNK constraints, and explore the properties of the resulting model for a series of (0+1)D and (1+1)D tests in 4D Minkowski spacetime. These tests include a comparison of the dissipation mechanisms in Eckart, BDNK, and Muller-Israel-Stewart theories, as well as investigations of the impact of hydrodynamic frame on the causality and stability properties of Bjorken flow, planar shockwave, and heat flow solutions.

preprint2023arXiv

Impact of a mean field dynamo on neutron star mergers leading to magnetar remnants

We investigate the impact of a mean field model for the $αΩ$-dynamo potentially active in the post-merger phase of a binary neutron star coalescence. We do so by deriving equations for ideal general relativistic magnetohydrodynamics (GRMHD) with an additional $α-$term, which closely resemble their Newtonian counterpart, but remain compatible with standard numerical relativity simulations. We propose a heuristic dynamo closure relation for the magnetorotational instability-driven turbulent dynamo in the outer layers of a differentially rotating magnetar remnant and its accretion disk. As a first demonstration, we apply this framework to the early stages of post-merger evolution ($\lesssim 50\, \rm ms$). We demonstrate that depending on the efficacy of the dynamo action, magnetically driven outflows can be present with their amount of baryon loading correlating with the magnetic field amplification. These outflows can also contain precursor flaring episodes before settling into a quasi-steady state. For the dynamo parameters explored in this work, we observe electromagnetic energy fluxes of up to $10^{50}\, \rm erg/s$, although larger amplification parameters will likely lead to stronger fluxes. Our results are consistent with the expectation that substantial dynamo amplification (either during or after the merger) may be necessary for neutron-star remnants to power short gamma-ray bursts or precursors thereof.

preprint2022arXiv

Characterizing the breakdown of quasi-universality in the post-merger gravitational waves from binary neutron star mergers

The post-merger gravitational wave (GW) emission from a binary neutron star merger is expected to provide exciting new constraints on the dense-matter equation of state (EoS). Such constraints rely, by and large, on the existence of quasi-universal relations, which relate the peak frequencies of the post-merger GW spectrum to properties of the neutron star structure in a model-independent way. In this work, we report on violations of existing quasi-universal relations between the peak spectral frequency, f_2, and the stellar radius, for EoSs models with backwards-bending slopes in their mass-radius relations (such that the radius increases at high masses). The violations are extreme, with variations in f_2 of up to ~600 Hz between EoSs that predict the same radius for a 1.4 Msun neutron star, but that have significantly different radii at higher masses. Quasi-universality can be restored by adding in a second parameter to the fitting formulae that depends on the slope of the mass-radius curve. We further find strong evidence that quasi-universality is better maintained for the radii of very massive stars (with masses 2 Msun). Both statements imply that f_2 is mainly sensitive to the high-density EoS. Combined with observations of the binary neutron star inspiral, these generalized quasi-universal relations can be used to simultaneously infer the characteristic radius and slope of the neutron star mass-radius relation.

preprint2022arXiv

Conservative finite volume scheme for first-order viscous relativistic hydrodynamics

We present the first conservative finite volume numerical scheme for the causal, stable relativistic Navier-Stokes equations developed by Bemfica, Disconzi, Noronha, and Kovtun (BDNK). BDNK theory has arisen very recently as a promising means of incorporating entropy-generating effects (viscosity, heat conduction) into relativistic fluid models, appearing as a possible alternative to the so-called Müller-Israel-Stewart (MIS) theory successfully used to model quark-gluon plasma. The major difference between the two lies in the structure of the system of PDEs: BDNK theory only has a set of conservation laws, whereas MIS also includes a set of evolution equations for its dissipative degrees of freedom. The simpler structure of the BDNK PDEs in this respect allows for rigorous proofs of stability, causality, and hyperbolicity in full generality which have as yet been impossible for MIS. To capitalize on these advantages, we present the first fully conservative multi-dimensional fluid solver for the BDNK equations suitable for physical applications. The scheme includes a flux-conservative discretization, non-oscillatory reconstruction, and a central-upwind numerical flux, and is designed to smoothly transition to a high-resolution shock-capturing perfect fluid solver in the inviscid limit. We assess the robustness of our new method in a series of flat-spacetime tests for a conformal fluid, and provide a detailed comparison with previous approaches of Pandya & Pretorius (2021).

preprint2022arXiv

Electromagnetic precursor flares from the late inspiral of neutron star binaries

The coalescence of two neutron stars is accompanied by the emission of gravitational waves, and can also feature electromagnetic counterparts powered by mass ejecta and the formation of a relativistic jet after the merger. Since neutron stars can feature strong magnetic fields, the non-trivial interaction of the neutron star magnetospheres might fuel potentially powerful electromagnetic transients prior to merger. A key process powering those precursor transients is relativistic reconnection in strong current sheets formed between the two stars. In this work, we provide a detailed analysis of how the twisting of the common magnetosphere of the binary leads to an emission of electromagnetic flares, akin to those produced in the solar corona. By means of relativistic force-free electrodynamics simulations, we clarify the role of different magnetic field topologies in the process. We conclude that flaring will always occur for suitable magnetic field alignments, unless one of the neutron stars has a magnetic field significantly weaker than the other.

preprint2022arXiv

How Do Spherical Black Holes Grow Monopole Hair?

Black holes in certain modified gravity theories that contain a scalar field coupled to curvature invariants are known to possess (monopole) scalar hair while non-black-hole spacetimes (like neutron stars) do not. Therefore, as a neutron star collapses to a black hole, scalar hair must grow until it settles to the stationary black hole solution with (monopole) hair. In this paper, we study this process in detail and show that the growth of scalar hair is tied to the appearance and growth of the event horizon (before an apparent horizon forms), which forces scalar modes that would otherwise (in the future) become divergent to be radiated away. We prove this result rigorously in general first for a large class of modified theories, and then we exemplify the results by studying the temporal evolution of the scalar field in scalar Gauss-Bonnet gravity in two backgrounds: (i) a collapsing Oppenheimer-Snyder background, and (ii) a collapsing neutron star background. In case (i), we find an exact scalar field solution analytically, while in case (ii) we solve for the temporal evolution of the scalar field numerically, with both cases supporting the conclusion presented above. Our results suggest that the emission of a burst of scalar field radiation is a necessary condition for black hole formation in a large class of modified theories of gravity.

preprint2022arXiv

Impact of extreme spins and mass ratios on the post-merger observables of high-mass binary neutron stars

The gravitational-wave events GW170817 and GW190425 have led to a number of important insights on the equation of state of dense matter and the properties of neutron stars, such as their radii and the maximum mass. Some of these conclusions have been drawn on the basis of numerical-relativity simulations of binary neutron-star mergers with vanishing initial spins. While this may be a reasonable assumption in equal-mass systems, it may be violated in the presence of large mass asymmetries accompanied by the presence of high spins. To quantify the impact of high spins on multi-messenger gravitational-wave events, we have carried out a series of high-mass binary neutron-star mergers with a highly spinning primary star and large mass asymmetries that have been modelled self-consistently using two temperature-dependent equations of state. We show that, when compared with equal-mass, irrotational binaries, these systems can lead to significant differences in the remnant lifetime, in the dynamical ejecta, in the remnant disc masses, in the secular ejecta, and on the bulk kilonova properties. These differences could be exploited to remove the degeneracy between low- and high-spin priors in the detection of gravitational waves from binary neutron-star mergers.

preprint2022arXiv

Magnetar bursts due to Alfvén wave nonlinear breakout

The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from $10^{36}$ to $10^{43}\ {\rm erg}\,{\rm s}^{-1}$. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, $\lesssim 1$ sr if the quake launching region is small, $\lesssim 0.01$ sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154.

preprint2022arXiv

Modeling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modeling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron-ion or electron-positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann-Vlasov equation. Crucially, this new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. In particular, we keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion disks. This new formulation will be instrumental in the construction of a large class of next-generation simulations of relativistic transient phenomena produced around black holes and neutron stars.

preprint2021arXiv

GW170817 and GW190814: tension on the maximum mass

The detection of the binary events GW170817 and GW190814 has provided invaluable constraints on the maximum mass of nonrotating configurations of neutron stars, $M_{_{\rm TOV}}$. However, the large differences in the neutron-star masses measured in GW170817 and GW190814 has also lead to a significant tension between the predictions for such maximum masses, with GW170817 suggesting that $M_{_{\rm TOV}} \lesssim 2.3\,M_{\odot}$, and GW190814 requiring $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$ if the secondary was a (non- or slowly rotating) neutron star at merger. Using a genetic algorithm, we sample the multidimensional space of parameters spanned by gravitational-wave and astronomical observations associated with GW170817. Consistent with previous estimates, we find that all of the physical quantities are in agreement with the observations if the maximum mass is in the range $M_{_{\rm TOV}} = 2.210^{+0.116}_{-0.123} \,M_{\odot}$ within a $2\textrm{-}σ$ confidence level. By contrast, maximum masses with $M_{_{\rm TOV}} \gtrsim 2.5\,M_{\odot}$, not only require efficiencies in the gravitational-wave emission that are well above the numerical-relativity estimates, but they also lead to a significant under-production of the ejected mass. Hence, the tension can be released by assuming that the secondary in GW190814 was a black hole at merger, although it could have been a rotating neutron star before.

preprint2020arXiv

Beyond second-order convergence in simulations of magnetised binary neutron stars with realistic microphysics

We investigate the impact of using high-order numerical methods to study the merger of magnetised neutron stars with finite-temperature microphysics and neutrino cooling in full general relativity. By implementing a fourth-order accurate conservative finite-difference scheme we model the inspiral together with the early post-merger and highlight the differences to traditional second-order approaches at the various stages of the simulation. We find that even for finite-temperature equations of state, convergence orders higher than second order can be achieved in the inspiral and post-merger for the gravitational-wave phase. We further demonstrate that the second-order scheme overestimates the amount of proton-rich shock-heated ejecta, which can have an impact on the modelling of the dynamical part of the kilonova emission. Finally, we show that already at low resolution the growth rate of the magnetic energy is consistently resolved by using a fourth-order scheme.

preprint2020arXiv

Electromagnetic precursors to gravitational wave events: Numerical simulations of flaring in pre-merger binary neutron star magnetospheres

The detection of gravitational waves from neutron star merger events has opened up a new field of multi-messenger astronomy linking gravitational waves events to short-gamma ray bursts and kilonova afterglows. A further - yet to be discovered - electromagnetic counterpart is precursor emission produced by the non-trivial interaction of the magnetospheres of the two neutron stars prior to merger. By performing special-relativistic force-free simulations of orbiting neutron stars we discuss the effect of different magnetic field orientations and show how the emission can be significantly enhanced by differential motion present in the binary, either due to stellar spins or misaligned stellar magnetospheres. We find that the built-up of twist in the magnetic flux tube connecting the two stars can lead to the repeated emission of powerful flares for a variety of orbital configurations. We also discuss potential coherent radio emission mechanisms in the flaring process.

preprint2020arXiv

The heavier the better: how to constrain mass ratios and spins of high-mass neutron-star mergers

The first binary neutron-star merger event, GW170817, and its bright electromagnetic counterpart have provided a remarkable amount of information. By contrast, the second event, GW190425, with $M_{\rm tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic counterpart, has hardly improved our understanding of neutron-star physics. While GW190425 is compatible with a scenario in which the merger has lead to a prompt collapse to a black hole and little ejected matter to power a counterpart, determining the mass ratio and the effective spin $\tildeχ$ of the binary remains difficult. This is because gravitational waveforms cannot yet well constrain the component spins of the binary. However, since the mass of GW190425 is significantly larger than the maximum mass for nonrotating neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the heavier star would not be gravitationally stable. Making use of universal relations and a large number of equations of state, we provide limits in the $(\tildeχ,q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tildeχ_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\, M_\odot$. Finally, we show how future observations of high-mass binaries can provide a lower bound on $M_{_{\rm TOV}}$.

preprint2019arXiv

On the Deconfinement Phase Transition in Neutron-Star Mergers

We study in detail the nuclear aspects of a neutron-star merger in which deconfinement to quark matter takes place. For this purpose, we make use of the Chiral Mean Field (CMF) model, an effective relativistic model that includes self-consistent chiral symmetry restoration and deconfinement to quark matter and, for this reason, predicts the existence of different degrees of freedom depending on the local density/chemical potential and temperature. We then use the out-of-chemical-equilibrium finite-temperature CMF equation of state in full general-relativistic simulations to analyze which regions of different QCD phase diagrams are probed and which conditions, such as strangeness and entropy, are generated when a strong first-order phase transition appears. We also investigate the amount of electrons present in different stages of the merger and discuss how far from chemical equilibrium they can be and, finally, draw some comparisons with matter created in supernova explosions and heavy-ion collisions.