@article{UjevicRashtiGiegetal.2022,
  author    = {Ujevic, Maximiliano and Rashti, Alireza and Gieg, Henrique Leonhard and Tichy, Wolfgang and Dietrich, Tim},
  title     = {High-accuracy high-mass-ratio simulations for binary neutron stars and their comparison to existing waveform models},
  series = {Physical review : D, Particles, fields, gravitation, and cosmology},
  volume    = {106},
  journal   = {Physical review : D, Particles, fields, gravitation, and cosmology},
  number    = {2},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {2470-0010},
  doi       = {10.1103/PhysRevD.106.023029},
  pages     = {10},
  year      = {2022},
  abstract  = {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.7 M???. Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately second-order converging results for the dominant (2,2) mode, but also the subdominant (2,1), (3,3), and (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 used 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.},
  language  = {en}
}
@article{RashtiFabbriBruegmannetal.2022,
  author    = {Rashti, Alireza and Fabbri, Francesco Maria and Br{\"u}gmann, Bernd and Chaurasia, Swami Vivekanandji and Dietrich, Tim and Ujevic, Maximiliano and Tichy, Wolfgang},
  title     = {New pseudospectral code for the construction of initial data},
  series = {Physical review D},
  volume    = {105},
  journal   = {Physical review D},
  number    = {10},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {2470-0010},
  doi       = {10.1103/PhysRevD.105.104027},
  pages     = {18},
  year      = {2022},
  abstract  = {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 pseudospectral 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 ∼0.8.},
  language  = {en}
}
@article{DudiDietrichRashtietal.2022,
  author    = {Dudi, Reetika and Dietrich, Tim and Rashti, Alireza and Br{\"u}gmann, Bernd and Steinhoff, Jan and Tichy, Wolfgang},
  title     = {High-accuracy simulations of highly spinning binary neutron star systems},
  series = {Physical review : D, Particles, fields, gravitation, and cosmology},
  volume    = {105},
  journal   = {Physical review : D, Particles, fields, gravitation, and cosmology},
  number    = {6},
  publisher = {American Physical Society},
  address   = {College Park},
  issn      = {2470-0010},
  doi       = {10.1103/PhysRevD.105.064050},
  pages     = {13},
  year      = {2022},
  abstract  = {With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and antialigned spins within a range of dimensionless spins of chi similar to [-0.28, 0.58]. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 arc not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.},
  language  = {en}
}