# Examples in STTs

Here we present a few examples to show how to work with the code and the related performances. All
cases have been run on a laptop.

## Non-rotating, unmagnetised NS in STT with the APR EoS

This is a model of an unmagnetised NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the APR EoS. It has a J-frame central density $\rho _\mathrm{c}=1.4\times 10^{-3}$ in code units (corresponding to $8.684\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.086$ M$ _\odot$. The circumferential radius in the J-frame is $10.963$ km, and the E-frame scalar charge is $0.532$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 10, RMAXSTR = 100,
RMAX = 100, REQMAX = 11.50, RHOINI = 1.40E-3, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .TRUE., FILEEOS = 'APR_resampled.dat'
```

| ![alt-text-1](./images/stt/apr/starplot_rho_b-6_unmag.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/apr/starplot_chi_b-6_unmag.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/apr/rho_b-6_unmag.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/p_b-6_unmag.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/chi_b-6_unmag.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/alpha_b-6_unmag.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/psi_b-6_unmag.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

## Non-rotating NS in STT with the APR EoS and a purely poloidal magnetic field

This is a model of an NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the APR EoS, endowed with a purely poloidal field. It has a J-frame central density $\rho _\mathrm{c}=1.4\times 10^{-3}$ in code units (corresponding to $8.684\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.098$ M$ _\odot$. The circumferential radius in the J-frame is $11.072$ km, and the E-frame scalar charge is $0.551$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 10, RMAXSTR = 100,
RMAX = 100, REQMAX = 11.50, RHOINI = 1.40E-3, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .TRUE., FILEEOS = 'APR_resampled.dat',
IMAG = .TRUE., IPOL = .TRUE., KBPOL = 0.2, NPOL = 0.0, CSI = 0.0
```

| ![alt-text-1](./images/stt/apr/starplot_rho_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/apr/starplot_chi_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/apr/starplot_b_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the poloidal magnetic field strength and its field lines.* |

| ![alt-text-1](./images/stt/apr/rho_b-6_polo.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/p_b-6_polo.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/chi_b-6_polo.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/alpha_b-6_polo.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/psi_b-6_polo.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/bpol_b-6_polo.png) |
|:--:|
| *Radial profile of the poloidal magnetic field strength at the pole and at the equator.* |

## Non-rotating NS in STT with the APR EoS and a purely toroidal magnetic field

This is a model of an NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the APR EoS, endowed with a purely toroidal field. It has a J-frame central density $\rho _\mathrm{c}=1.4\times 10^{-3}$ in code units (corresponding to $8.684\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.078$ M$ _\odot$. The circumferential radius in the J-frame is $11.060$ km, and the E-frame scalar charge is $0.515$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 10, RMAXSTR = 100,
RMAX = 100, REQMAX = 11.50, RHOINI = 1.40E-3, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .TRUE., FILEEOS = 'APR_resampled.dat',
IMAG = .TRUE., ITOR = .TRUE., BCOEF = 0.7, NPOL = 0.0, CSI = 0.0
```

| ![alt-text-1](./images/stt/apr/starplot_rho_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/apr/starplot_chi_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/apr/starplot_b_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the toroidal magnetic field strength.* |

| ![alt-text-1](./images/stt/apr/rho_b-6_toro.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/p_b-6_toro.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/chi_b-6_toro.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/alpha_b-6_toro.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/psi_b-6_toro.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

| ![alt-text-1](./images/stt/apr/btor_b-6_toro.png) |
|:--:|
| *Radial profile of the toroidal magnetic field strength at the pole and at the equator.* |

## Non-rotating, unmagnetised NS in STT with the POL2 EoS

This is a model of an unmagnetised NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the analytic POL2 EoS. It has a J-frame central density $\rho _\mathrm{c}=1.36\times 10^{-3}$ in code units (corresponding to $8.44\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.296$ M$ _\odot$. The circumferential radius in the J-frame is $14.156$ km, and the E-frame scalar charge is $0.645$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 10, RMAXSTR = 100,
RMAX = 100, REQMAX = 15.0, RHOINI = 1.36E-3, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .FALSE., K1 = 110.0, GAMMA = 2.0
```

| ![alt-text-1](./images/stt/pol2/starplot_rho_b-6_unmag.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/pol2/starplot_chi_b-6_unmag.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/pol2/rho_b-6_unmag.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/p_b-6_unmag.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/chi_b-6_unmag.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/alpha_b-6_unmag.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/psi_b-6_unmag.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

## Non-rotating NS in STT with the POL2 EoS and a purely poloidal magnetic field

This is a model of an NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the analytic POL2 EoS, endowed with a purely poloidal field. It has a J-frame central density $\rho _\mathrm{c}=8.30\times 10^{-4}$ in code units (corresponding to $5.15\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.356$ M$ _\odot$. The circumferential radius in the J-frame is $16.718$ km, and the E-frame scalar charge is $0.564$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 20, RMAXSTR = 100,
RMAX = 100, REQMAX = 25.0, RHOINI = 8.30E-4, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .FALSE., K1 = 110.0, GAMMA = 2.0,
IMAG = .TRUE., IPOL = .TRUE., KBPOL = 0.44, NPOL = 0.0, CSI = 0.0
```

| ![alt-text-1](./images/stt/pol2/starplot_rho_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/pol2/starplot_chi_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/pol2/starplot_b_b-6_polo.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the poloidal magnetic field strength and its field lines.* |

| ![alt-text-1](./images/stt/pol2/rho_b-6_polo.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/p_b-6_polo.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/chi_b-6_polo.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/alpha_b-6_polo.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/psi_b-6_polo.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/bpol_b-6_polo.png) |
|:--:|
| *Radial profile of the poloidal magnetic field strength at the pole and at the equator.* |

## Non-rotating NS in STT with the POL2 EoS and a purely toroidal magnetic field

This is a model of an NS in STTs with $\alpha _0 = -2.0\times 10^{-4}$ and $\beta _0 = -6$ described by the analytic POL2 EoS, endowed with a purely toroidal field. It has a J-frame central density $\rho _\mathrm{c}=1.36\times 10^{-3}$ in code units (corresponding to $8.44\times 10^{14}$ gcm$^{-3}$) and a Komar mass in the E-frame of $1.467$ M$ _\odot$. The circumferential radius in the J-frame is $20.916$ km, and the E-frame scalar charge is $0.473$ M$ _\odot$.

Here are the main parameters used to compute this model. The detailed log with all the parameters can be found in the corresponding LogFile.dat.

```
NR = 900, NTH = 100, NRREG = 600, MLS = 20, NGQ = 50, RREG = 20, RMAXSTR = 100,
RMAX = 100, REQMAX = 25.0, RHOINI = 1.36E-3, ALPHA0 = -2.0E-4, BETA0 = -6, GR = .FALSE.,
STRETCH = .TRUE., ANALYTIC = .TRUE., CONVHELP = .FALSE., QFACTOR = 0.85, QFACTORCHI = 0.45,
QFACTORMETRIC = 0.35, QRELAX = 0.30, QAPHI = 0.50, EOSINT = .FALSE., K1 = 110.0, GAMMA = 2.0,
IMAG = .TRUE., ITOR = .TRUE., BCOEF = 4.0, NPOL = 0.0, CSI = 0.0
```

| ![alt-text-1](./images/stt/pol2/starplot_rho_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the density.* |

| ![alt-text-1](./images/stt/pol2/starplot_chi_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the scalar field.* |

| ![alt-text-1](./images/stt/pol2/starplot_b_b-6_toro.png) |
|:--:|
| *Section of the NS in the $x-z$ plane along with the contours the toroidal magnetic field strength.* |

| ![alt-text-1](./images/stt/pol2/rho_b-6_toro.png) |
|:--:|
| *Radial profile of the NS density at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/p_b-6_toro.png) |
|:--:|
| *Radial profile of the NS pressure at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/chi_b-6_toro.png) |
|:--:|
| *Radial profile of the scalar field at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/alpha_b-6_toro.png) |
|:--:|
| *Radial profile of the lapse function at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/psi_b-6_toro.png) |
|:--:|
| *Radial profile of the conformal factor at the pole and at the equator.* |

| ![alt-text-1](./images/stt/pol2/btor_b-6_toro.png) |
|:--:|
| *Radial profile of the toroidal magnetic field strength at the pole and at the equator.* |