Source author record

Maximiliano Ujevic

Maximiliano Ujevic 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

Elliptica: a new pseudo-spectral code for the construction of initial data

Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, Elliptica deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of Elliptica and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude $\sim 0.8$.

preprint2022arXiv

High-accuracy high-mass ratio simulations for binary neutron stars and their comparison to existing waveform models

The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical-relativity simulations for four different physical setups with mass ratios $q$ = $1.25$, $1.50$, $1.75$, $2.00$, and total gravitational mass $M = 2.7M_\odot$ . Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately 2nd order converging results for the dominant $(2,2)$, but also subdominant $(2,1)$, $(3,3)$, $(4,4)$ modes, while, generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model to describe tidal effects is a valid description for high-mass ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs.

preprint2022arXiv

Incorporating a radiative hydrodynamics scheme in the numerical-relativity code BAM

To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, e.g., through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory based Equations of State and a simple description of neutrino interactions through a Neutrino Leakage Scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB$Λϕ$. Overall, our test simulations have good agreement with those reported in the literature.

preprint2021arXiv

Investigating the mass-ratio dependence of the prompt-collapse threshold with numerical-relativity simulations

The next observing runs of advanced gravitational-wave detectors will lead to a variety of binary neutron star detections and numerous possibilities for multi-messenger observations of binary neutron star systems. In this context a clear understanding of the merger process and the possibility of prompt black hole formation after merger is important, as the amount of ejected material strongly depends on the merger dynamics. These dynamics are primarily affected by the total mass of the binary, however, the mass ratio also influences the postmerger evolution. To determine the effect of the mass ratio, we investigate the parameter space around the prompt-collapse threshold with a new set of fully relativistic simulations. The simulations cover three equations of state and seven mass ratios in the range of $1.0 \leq q \leq 1.75$, with five to seven simulations of binary systems of different total mass in each case. The threshold mass is determined through an empirical relation based on the collapse-time, which allows us to investigate effects of the mass-ratio on the threshold mass and also on the properties of the remnant system. Furthermore, we model effects of mass ratio and equation of state on tidal parameters of threshold configurations.

preprint2020arXiv

Gravitational waves and mass ejecta from binary neutron star mergers: Effect of the spin orientation

We continue our study of the binary neutron star parameter space by investigating the effect of the spin orientation on the dynamics, gravitational wave emission, and mass ejection during the binary neutron star coalescence. We simulate seven different configurations using multiple resolutions to allow a reasonable error assessment. Due to the particular choice of the setups, five configurations show precession effects, from which two show a precession ("wobbling") of the orbital plane, while three show a "bobbing" motion, i.e., the orbital angular momentum does not precess, while the orbital plane moves along the orbital angular momentum axis. Considering the ejection of mass, we find that precessing systems can have an anisotropic mass ejection, which could lead to a final remnant kick of $\sim 40 \rm km/s$ for the studied systems. Furthermore, for the chosen configurations, antialigned spins lead to larger mass ejecta than aligned spins, so that brighter electromagnetic counterparts could be expected for these configurations. Finally, we compare our simulations with the precessing, tidal waveform approximant IMRPhenomPv2_NRTidalv2 and find good agreement between the approximant and our numerical relativity waveforms with phase differences below 1.2 rad accumulated over the last $\sim$ 16 gravitational wave cycles.

preprint2015arXiv

Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

We study equal and unequal-mass neutron star mergers by means of new numerical relativity simulations in which the general relativistic hydrodynamics solver employs an algorithm that guarantees mass conservation across the refinement levels of the computational mesh. We consider eight binary configurations with total mass $M=2.7\,M_\odot$, mass-ratios $q=1$ and $q=1.16$, and four different equation of states (EOSs), and one configuration with a stiff EOS, $M=2.5M_\odot$ and $q=1.5$. We focus on the post-merger dynamics and study the merger remnant, dynamical ejecta and the postmerger gravitational wave spectrum. Although most of the merger remnants form a hypermassive neutron star collapsing to a black hole+disk system on dynamical timescales, stiff EOSs can eventually produce a stable massive neutron star. Ejecta are mostly emitted around the orbital plane; favored by large mass ratios and softer EOS. The postmerger wave spectrum is mainly characterized by non-axisymmetric oscillations of the remnant. The stiff EOS configuration consisting of a $1.5M_\odot$ and a $1.0M_\odot$ neutron star shows a rather peculiar dynamics. During merger the companion star is very deformed; about~$\sim0.03M_\odot$ of rest-mass becomes unbound from the tidal tail due torque; and the merger remnant forms stable neutron star surrounded by a massive accretion disk $\sim0.3M_\odot$. Similar configurations might be particularly interesting for electromagnetic counterparts. Comparing results obtained with and without the conservative mesh refinement algorithm, we find that post-merger simulations can be affected by systematic errors if mass conservation is not enforced in the mesh refinement strategy. However, mass conservation also depends on grid details and on the artificial atmosphere setup. [abridged]

preprint2014arXiv

Quasiuniversal properties of neutron star mergers

Binary neutron star mergers are studied using nonlinear 3+1 numerical relativity simulations and the analytical effective-one-body (EOB) model. The EOB model predicts quasiuniversal relations between the mass-rescaled gravitational wave frequency and the binding energy at the moment of merger, and certain dimensionless binary tidal coupling constants depending on the stars Love numbers, compactnesses and the binary mass ratio. These relations are quasiuniversal in the sense that, for a given value of the tidal coupling constant, they depend significantly neither on the equation of state nor on the mass ratio, though they do depend on stars spins. The spin dependence is approximately linear for small spins aligned with the orbital angular momentum. The quasiuniversality is a property of the conservative dynamics; nontrivial relations emerge as the binary interaction becomes tidally dominated. This analytical prediction is qualitatively consistent with new, multi-orbit numerical relativity results for the relevant case of equal-mass irrotational binaries. Universal relations are thus expected to characterize neutron star mergers dynamics. In the context of gravitational wave astronomy, these universal relations may be used to constrain the neutron star equation of state using waveforms that model the merger accurately.

preprint2009arXiv

Potential flows in the Reissner-Nordström-(anti) de Sitter metric: numerical results

Numerical solutions for the integral curves of the velocity field (streamlines), the density contours, and the accretion rate of a steady-state flow of an ideal fluid with $p=K n^γ$ equation of state are presented. The streamlines and velocity fields associated with a black hole and a rigid sphere in a Reissner-Nordström-(anti) de Sitter spacetimes are studied in some detail. For each case the fluid density is studied using contour lines. For $γ\neq 2$, we found that the studied properties of the fluid are more sensitive to variations of the electric charge and the cosmological constant. Also, the accretion rate was found to increase or decrease when the cosmological constant increases or decreases respectively.

preprint2002arXiv

Potential flows in a core-dipole-shell system: numerical results

Numerical solutions for: the integral curves of the velocity field (streamlines), the density contours, and the accretion rate of a steady-state flow of an ideal fluid with p=K n^(gamma) equation of state orbiting in a core-dipole-shell system are presented. For 1 < gamma < 2, we found that the non-linear contribution appearing in the partial differential equation for the velocity potential has little effect in the form of the streamlines and density contour lines, but can be noticed in the density values. The study of several cases indicates that this appears to be the general situation. The accretion rate was found to increase when the constant gamma decreases.