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

Antonios Tsokaros

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEhep-th

Catalog footprint

What is connected

18works
3topics
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

18 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Asymptotic vacua with higher derivatives

We study limits of vacuum, isotropic universes in the full, effective, four-dimensional theory with higher derivatives. We show that all flat vacua as well as general curved ones are globally attracted by the standard, square root scaling solution at early times. Open vacua asymptote to horizon-free, Milne states in both directions while closed universes exhibit more complex logarithmic singularities, starting from initial data sets of a possibly smaller dimension. We also discuss the relation of our results to the asymptotic stability of the passage through the singularity in ekpyrotic and cyclic cosmologies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Initial-data contribution to the error budget of gravitational waves from neutron-star binaries

As numerical calculations of inspiralling neutron-star binaries reach values of accuracy that are comparable with those of binary black holes, a fine budgeting of the various sources of error becomes increasingly important. Among such sources, the initial data is normally not accounted for, the rationale being that the error on the initial spacelike hypersurface is always far smaller than the one gained during the evolution. We here consider critically this assumption and perform a comparative analysis of the gravitational waveforms relative to essentially the same physical binary configuration when computed with two different initial-data codes, and then evolved with the same evolution code. More specifically, we consider the evolution of irrotational neutron-star binaries computed either with the pseudo-spectral code \lorene{}, or with the newly developed finite-difference code \cocal{}; both sets of initial data are subsequently evolved with the high-order evolution code \whiskythc{}. In this way we find that despite the initial data shows global (local) differences that are $\lesssim 0.02\%\ (1\%)$, the gravitational-wave phase at the merger time differs by $\sim 0.5$ radians after $\sim 3$ orbits, a surprisingly large value. Our results highlight the highly nonlinear impact that errors in the initial data can have on the subsequent evolution and the importance of using exactly the same initial data when comparative studies are done.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New code for quasiequilibrium initial data of binary neutron stars: Corotating, irrotational and slowly spinning systems

We present the extension of our \cocal~- Compact Object CALculator - code to compute general-relativistic initial data for binary compact-star systems. In particular, we construct quasiequilibrium initial data for equal-mass binaries with spins that are either aligned or antialigned with the orbital angular momentum. The Isenberg-Wilson-Mathews formalism is adopted and the constraint equations are solved using the representation formula with a suitable choice of a Green's function. We validate the new code with solutions for equal-mass binaries and explore its capabilities for a wide range of compactnesses, from a white dwarf binary with compactness $\sim 10^{-4}$, up to a highly relativistic neutron-star binary with compactness $\sim 0.22$. We also present a comparison with corotating and irrotational quasiequilibrium sequences from the spectral code \lorene [Taniguchi and Gourgoulhon, Phys. Rev. D {\bf 66}, 104019 (2002)] and with different compactness, showing that the results from the two codes agree to a precision of the order of $0.05\%$. Finally, we present equilibria for spinning configurations with a nuclear-physics equation of state in a piecewise polytropic representation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equilibrium solutions of relativistic rotating stars with mixed poloidal and toroidal magnetic fields

Stationary and axisymmetric solutions of relativistic rotating stars with strong mixed poloidal and toroidal magnetic fields are obtained numerically. Because of the mixed components of the magnetic field, the underlying stationary and axisymmetric spacetimes are no longer circular. These configurations are computed from the full set of the Einstein-Maxwell equations, Maxwell's equations and from first integrals and integrability conditions of the magnetohydrodynamic equilibrium equations. After a brief introduction of the formulation of the problem, we present the first results for highly deformed magnetized rotating compact stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asymptotic states of generalized universes with higher derivatives

We review ongoing research related to the asymptotic dynamics of isotropic universes in theories with higher derivatives, especially near the initial singularity. We treat two major cases, that is universes in vacuum, and also those filled with radiation using the method of asymptotic splittings of vector fields. Our solutions have the form of formal Fuchsian-type series for the basic unknowns and are valid near collapse or big rip regimes asymptotically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The initial state of generalized radiation universes

We use asymptotic methods to study the early time stability of isotropic and homogeneous solutions filled with radiation which are close initially to the exact, flat, radiation solution in quadratic lagrangian theories of gravity. For such models, we analyze all possible modes of approach to the initial singularity and prove the essential uniqueness and stability of the resulting asymptotic scheme in all cases except perhaps that of the conformally invariant Bach-Weyl gravity. We also provide a formal series representation valid near the initial singularity of the general solution of these models and show that this is dominated at early times by a form in which both curvature and radiation play a subdominant role. We also discuss the implications of these results for the generic initial state of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Binary black hole circular orbits computed with COCAL

In this work we present our first results of binary black hole circular orbits using {\sc cocal}, the Compact Object CALculator. Using the 3+1 decomposition five equations are being solved under the assumptions of conformal flatness and maximal slicing. Excision is used and the appropriate apparent horizon boundary conditions are applied. The orbital velocity is determined by imposing a Schwarzschild behaviour at infinity. A sequence of equal mass black holes is obtained and its main physical characteristics are calculated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

New code for equilibriums and quasiequilibrium initial data of compact objects. II. Convergence tests and comparisons of binary black hole initial data

COCAL is a code for computing equilibriums or quasiequilibrium initial data of single or binary compact objects based on finite difference methods. We present the results of supplementary convergence tests of COCAL code using time symmetric binary black hole data (Brill-Lindquist solution). Then, we compare the initial data of binary black holes on the conformally flat spatial slice obtained from COCAL and KADATH, where KADATH is a library for solving a wide class of problems in theoretical physics including relativistic compact objects with spectral methods. Data calculated from the two codes converge nicely towards each other, for close as well as largely separated circular orbits of binary black holes. Finally, as an example, a sequence of equal mass binary black hole initial data with corotating spins is calculated and compared with data in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new code for equilibriums and quasiequilibrium initial data of compact objects

We present a new code, named COCAL - Compact Object CALculator, for the computation of equilibriums and quasiequilibrium initial data sets of single or binary compact objects of all kinds. In the cocal code, those solutions are calculated on one or multiple spherical coordinate patches covering the initial hypersurface up to the asymptotic region. The numerical method used to solve field equations written in elliptic form is an adaptation of self-consistent field iterations in which Green's integral formula is computed using multipole expansions and standard finite difference schemes. We extended the method so that it can be used on a computational domain with excised regions for a black hole and a binary companion. Green's functions are constructed for various types of boundary conditions imposed at the surface of the excised regions for black holes. The numerical methods used in cocal are chosen to make the code simpler than any other recent initial data codes, accepting the second order accuracy for the finite difference schemes. We perform convergence tests for time symmetric single black hole data on a single coordinate patch, and binary black hole data on multiple patches. Then, we apply the code to obtain spatially conformally flat binary black hole initial data using boundary conditions including the one based on the existence of equilibrium apparent horizons.

0 citations0 reviewsNo code linkedOpen access