[Users] hybrid star setup in the TOVSolver thorn

[Wolfgang Kastaun]

Hi,

there is another TOV solver available in the "reprimand" thorn. It only
solves the NS equations, however, but does not set up any initial data.
It provides ordinary C++ functions for EOS handling and TOV solving,
but there is no cactus interface. To set up initial data, you would
need to write a small thorn or replace the TOV solver part in the 
available TOVSolver thorn with calls to reprimand.
The thorn is a based on the standalone C++ library also named reprimand
(which also has a python interface that can be pip-installed).
https://github.com/wokast/RePrimAnd
It is documented here
https://wokast.github.io/RePrimAnd/index.html

It might be well suited for your special use case since it is designed
for dealing also with discontinuous EOS, as described in this article:
https://arxiv.org/abs/2404.11346


Cheers,
Wolfgang.

On Thu, 2024-10-31 at 20:58 +0000, CJ Osakwe wrote:
> 
> Hello,
> 
> 
> I am trying to set up a stable hybrid star (with a quark matter core
> and hadronic crust, each with their own equation of state) in the
> Einstein Toolkit. My initial approach was to adapt the TOVSolver
> thorn to achieve this, but the code is removing the density
> discontinuity at the hadronic-quark matter interface when it
> interpolates from 1D to 3D. As well the discontinuity doesn't seem to
> be reflected in the pressure anyway. I'm wondering if there is a
> thorn that is better suited for setting up a stable hybrid star?
> 
> 
> Cheers,
> CJ Osakwe
> PhD candidate, Department of Physics and Astronomy
> University of Calgary
> 
> 
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users