Code for obtaining neutron star models for tabulated equations of state1, in the context of Damour-Esposito-Farese (DEF) scalar tensor theory model and R^2 gravity model. The system of ODEs is solved by implementing a C++11 version of LSODA2 and a C++ version of the Nelder-Mead method3, which is used for the optimization of the boundary conditions. Both of these libraries are slightly modified to fit the STT class in stt.h. Additionaly, some modules from the BOOST C++ library4 are used to interpolate the tabulated EoS. The code prints the model's parameters e.g. mass, radius, scalar charge, scalar at center and boundary minimization accuracy, for the corresponding input parameters.
The C++ code is wrapped using SWIG5, into a python library which can be imported and used in any way, as shown in swig/test.py.
Compile, using the provided makefile.
STT.exe takes 5 inputs in total. The parameters are specified using the following flags:
- -m theory (DEF or R2)
- -f eos name (EoS file name).
- -c coupling (model's coupling constant).
- -e central density (central energy density in CGS/10^15).
For example,
./STT -m DEF -f ppsly4.cold.rns1.1.txt -c -5.0 -e 1.5
./STT -m R2 -f ppsly4.cold.rns1.1.txt -c 20.0 -e 1.5
typical input parameter space:
- central density ~ 0.5-4.0 (10^15 gr/cm^3)
- coupling ~ (-5.0,-4.5) (DEF), (5,10^3) (R^2)
GR solutions may also be obtained simply by choosing "DEF" as theory and setting the coupling to zero.