Researcher profile

Antonios Tsokaros

Antonios Tsokaros contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
2topics
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

8 published item(s)

preprint2022arXiv

Jet Launching from Binary Neutron Star Mergers: Incorporating Neutrino Transport and Magnetic Fields

We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of merging binary neutron stars incorporating neutrino transport and magnetic fields. Our new radiative transport module for neutrinos adopts a general relativistic, truncated-moment (M1) formalism. The binaries consist of two identical, irrotational stars modeled by the SLy nuclear equation of state (EOS). They are initially in quasicircular orbit and threaded with a poloidal magnetic field that extends from the stellar interior into the exterior, as in typical pulsars. We insert neutrino processes shortly after the merger and focus on the role of neutrinos in launching a jet following the collapse of the hypermassive neutron star (HMNS) remnant to a spinning black hole (BH). We treat two microphysical versions: one (a "warm-up") evolving a single neutrino species and considering only charged-current processes, and the other evolving three species $(ν_e, \barν_e, ν_{\rm x})$ and related processes. We trace the evolution until the system reaches a quasiequilibrium state after BH formation. We find that the BH + disk remnant eventually launches an incipient jet. The electromagnetic Poynting luminosity is $\sim 10^{53} \rm \, erg\, s^{-1}$, consistent with that of typical short gamma-ray bursts (sGRBs). The effect of neutrino cooling shortens the lifetime of the HMNS, and lowers the amplitude of the major peak of the gravitational wave (GW) power spectrum somewhat. After BH formation, neutrinos help clear out the matter near the BH poles, resulting in lower baryon-loaded surrounding debris. The neutrino luminosity resides in the range $\sim 10^{52-53} \rm \,erg\,s^{-1}$ once quasiequilibrium is achieved. Comparing with the neutrino-free models, we observe that the inclusion of neutrinos yields similar ejecta masses and is inefficient in carrying off additional angular momentum.

preprint2021arXiv

Black hole-neutron star coalescence: effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (sGRBs) and kilonovae, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin ($a_{\rm BH}/M_{\rm BH}=0$ and $0.75$), the NS spin ($a_{\rm NS}/M_{\rm NS}=-0.17,\,0,\,0.23$ and $0.33$), and the binary mass ratio ($q\equiv M_{\rm BH}:M_{\rm NS}=3:1$ and $5:1$). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time~$Δt\sim 3500-5500M\sim 88-138(M_{\rm NS}/1.4M_\odot)\,\rm ms$ after the peak gravitational wave amplitude, a magnetically--driven jet is launched only for $q=3:1$ regardless of the initial NS spin. The lifetime of the jets [$Δt\sim 0.5-0.8(M_{\rm NS}/1.4 M_\odot)\,\rm s$] and their outgoing Poynting luminosity [$L_{\rm Poyn}\sim 10^{51.5\pm 0.5}\,\rm erg/s$] are consistent with typical sGRBs luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic field collimation for the other systems we considered. The mass range of dynamically ejected matter is $10^{-4.5}-10^{-2}~(M_{\rm NS}/1.4M_\odot)M_\odot$, which can power kilonovae with peak bolometric luminosities $L_{\rm knova}\sim 10^{40}-10^{41.4}$ erg/s with rise times $\lesssim 6.5\,\rm h$ and potentially detectable by the LSST.

preprint2021arXiv

Magnetohydrodynamic simulations of self-consistent rotating neutron stars with mixed poloidal and toroidal magnetic fields

We perform the first magnetohydrodynamic simulations in full general relativity of self-consistent rotating neutron stars (NSs) with ultrastrong mixed poloidal and toroidal magnetic fields. The initial uniformly rotating NS models are computed assuming perfect conductivity, stationarity, and axisymmetry. Although the specific geometry of the mixed field configuration can delay or accelerate the development of various instabilities known from analytic perturbative studies, all our models finally succumb to them. Differential rotation is developed spontaneously in the cores of our magnetars which, after sufficient time, is converted back to uniform rotation. The rapidly rotating magnetars show a significant amount of ejecta, which can be responsible for transient kilonova signatures. However no highly collimated, helical magnetic fields or incipient jets, which are necessary for gamma-ray bursts, arise at the poles of these magnetars by the time our simulations are terminated.

preprint2021arXiv

Methods for relativistic self-gravitating fluids: From binary neutron stars to black hole-disks and magnetized rotating neutron stars

The cataclysmic observations of event GW170817, first as gravitational waves along the inspiral motion of two neutron stars, then as a short $γ$-ray burst, and later as a kilonova, launched the era of multimessenger astronomy, and played a pivotal role in furthering our understanding on a number of longstanding questions. Numerical modelling of such multimessenger sources is an important tool to understand the physics of compact objects and, more generally, the physics of matter under extreme conditions. In this review we present a unified view of various techniques used to obtain equilibrium and quasiequilibrium solutions for three astrophysically relevant relativistic, self-gravitating fluid systems: Binary neutron stars, black hole-disks, and magnetized rotating neutron stars. These solutions are necessary not only for modeling such compact objects, but equally important, for providing self-consistent initial data in numerical relativity simulations. Instead of presenting the full details of the formulations and numerical algorithms, we focus on painting the broadbrush picture of the methods developed to address these problems, and facilitate future work in the area.

preprint2020arXiv

Gravitational Waves from Disks Around Spinning Black Holes: Simulations in Full General Relativity

We present fully general-relativistic numerical evolutions of self-gravitating tori around spinning black holes with dimensionless spin $a/M = 0.7$ parallel or anti-parallel to the disk angular momentum. The initial disks are unstable to the hydrodynamic Papaloizou-Pringle Instability which causes them to grow persistent orbiting matter clumps. The effect of black hole spin on the growth and saturation of the instability is assessed. We find that the instability behaves similarly to prior simulations with non-spinning black holes, with a shift in frequency due to spin-induced changes in disk orbital period. Copious gravitational waves are generated by these systems, and we analyze their detectability by current and future gravitational wave observatories for large range of masses. We find that systems of $10 M_\odot$ - relevant for black hole-neutron star mergers - are detectable by Cosmic Explorer out to $\sim300$ Mpc, while DECIGO (LISA) will be able to detect systems of $1000 M_\odot$ ($10^5M_\odot$) - relevant for disks forming in collapsing supermassive stars - out to cosmological redshift of $z\sim5$ ($z\sim 1$). Computing the accretion rate of these systems we find that these systems may also be promising sources of coincident electromagnetic signals.

preprint2020arXiv

Great Impostors: Extremely Compact, Merging Binary Neutron Stars in the Mass Gap Posing as Binary Black Holes

Can one distinguish a binary black hole undergoing a merger from a binary neutron star if the individual compact companions have masses that fall inside the so-called mass gap of $3-5\ M_\odot$? For neutron stars, achieving such masses typically requires extreme compactness and in this work we present initial data and evolutions of binary neutron stars initially in quasiequilibrium circular orbits having a compactness $C=0.336$. These are the most compact, nonvacuum, quasiequilibrium binary objects that have been constructed and evolved to date, including boson stars. The compactness achieved is only slightly smaller than the maximum possible imposed by causality, $C_{\rm max}=0.355$, which requires the sound speed to be less than the speed of light. By comparing the emitted gravitational waveforms from the late inspiral to merger and postmerger phases between such a binary neutron star vs a binary black hole of the same total mass we identify concrete measurements that serve to distinguish them. With that level of compactness, the binary neutron stars exhibit no tidal disruption up until merger, whereupon a prompt collapse is initiated even before a common core forms. Within the accuracy of our simulations the black hole remnants from both binaries exhibit ringdown radiation that is not distinguishable from a perturbed Kerr spacetime. However, their inspiral leads to phase differences of the order of $\sim 5$ rad over an $\sim 81$ km separation (1.7 orbits) while typical neutron stars exhibit phase differences of $\geq 20$ rad. Although a difference of $\sim 5$ rad can be measured by current gravitational wave laser interferometers (e.g., aLIGO/Virgo), uncertainties in the individual masses and spins will likely prevent distinguishing such compact, massive neutron stars from black holes.

preprint2020arXiv

Locating ergostar models in parameter space

Recently, we have shown that dynamically stable ergostar solutions (equilibrium neutron stars that contain an ergoregion) with a compressible and causal equation of state exist [A. Tsokaros, M. Ruiz, L. Sun, S. L. Shapiro, and K. Uryū, Phys. Rev. Lett. 123, 231103 (2019)]. These stars are hypermassive, differentially rotating, and highly compact. In this work, we make a systematic study of equilibrium models in order to locate the position of ergostars in parameter space. We adopt four equations of state that differ in the matching density of a maximally stiff core. By constructing a large number of models both with uniform and differential rotation of different degrees, we identify the parameters for which ergostars appear. We find that the most favorable conditions for the appearance of dynamically stable ergostars are a significant finite density close to the surface of the star (i.e., similar to self-bound quark stars) and a small degree of differential rotation.

preprint2020arXiv

Magnetohydrodynamic Simulations of Binary Neutron Star Mergers in General Relativity: Effects of Magnetic Field Orientation on Jet Launching

Binary neutron star (NSNS) mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of spinning NSNSs initially on quasicircular orbits that merge and undergo delayed collapse to a black hole (BH). The NSNSs consist of two identical stars modeled as $Γ=2$ polytropes with spin $χ_{NS}= 0.36$ aligned along the direction of the total orbital angular momentum $L$. Each star is initially threaded by a dynamical unimportant interior dipole B-field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment $μ$ is either aligned or perpendicular to $L$ and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where $μ$ in both stars is aligned along $L$. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta are very sensitive to the B-field orientation. By contrast, the physical properties of the BH + disk remnant, such as the mass and spin of the BH, the accretion rate, and the electromagnetic luminosity, are roughly independent of the initial B-field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results show that an incipient jet emerges only when the seed B-field has a sufficiently large-scale poloidal component aligned to $L$. The lifetime [$Δt\gtrsim 140(M_{NS}/1.625M_\odot)\rm ms$] and Poynting luminosities [$L_{EM}\simeq 10^{52}$erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford--Znajek mechanism for launching jets.