Here, numerical data described in the paper "Two boson stars in equilibrium", arXiv:2210.01833, is made available for public use.
The attachment "Data_files.zip", contains the data files:
- configuration-I.dat
- configuration-II.dat
- configuration-III.dat
These files correspond to the three representative dipolar BS (DBS) solutions among those described in the paper above. These solutions are:
Configuration I :
a DBS belonging to the main branch before the critical solution with maximal mass;
- the input parameter: \omega=0.9;
- ADM mass=0.9613; Noether charge Q=0.9884
========================================================
Configuration II :
a DBS belonging to the main branch after the critical solution with maximal mass;
- the input parameter: \omega=0.8
- ADM mass= 1.0156 ; Noether charge Q= 1.0484
========================================================
Configuration III :
a DBS belonging to the secondary branch (after the minimal frequency);
- the input parameter: \omega=0.75
- ADM mass= 0.6578 ; Noether charge Q= 0.5715
========================================================
The data is presented in the following order, in the files (F_i, \phi, are the metric and scalar functions used in the paper):
--------------------------------
x_1 theta_1 F_1 F_2 F_0 \phi
x_2 theta_1 F_1 F_2 F_0 \phi
...
x_251 theta_1 F_1 F_2 F_0 \phi
x_1 theta_2 F_1 F_2 F_0 \phi
x_2 theta_2 F_1 F_2 F_0 \phi
...
x_251 theta_2 F_1 F_2 F_0 \phi
x_1 theta_50 F_1 F_2 F_0 \phi
x_2 theta_50 F_1 F_2 F_0 \phi
...
x_251 theta_50 F_1 F_2 F_0 \phi
--------------------------------
where the grid points are:
x_k=(k-1)/250 (k=1,..,251)
and
theta_k=pi/2*(k-1)/49 (k=1,..,50)
The corresponding values for pi/2< \theta \leq pi result from the reflection symmetry of the solutions along the equatorial plane.
x=r/(1+r) is a compactified radial coordinate, 0 \leq x \leq 1.