Program 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 initial 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. Type "make" in the main directory, where the makefile is also located, to create the executionable file.
STT.exe takes 5 inputs in-turn. 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.