A numerical algorithm for modelling of boson-fermion stars in dilatonic gravity (Q1612391)

For the summary: We investigate numerically a class of models of the static spherically symmetric boson-fermion stars in the scalar-tensor theory of gravity with massive dilaton field. The proper mathematical model of such stars is interpreted as a nonlinear two-parametric eigenvalue problem. The first of the parameters is the unknown internal boundary (the radius of the fermionic part of the star) \(R_s\), and the second one represents the frequency \(\Omega\) of the time oscillations of the bosonic field. To solve this problem, the whole space \([0,\infty)\) is split into two domains: internal \([0,R_s]\) (inside the star) and external \([R_s,\infty)\) (outside the star). In each domain the physical model leads to two nonlinear boundary value problems with respect to metric functions, the functions describing the fermionic and bosonic matter, and the dilaton field. These boundary value problems have different dimensions inside and outside the star, respectively. The solutions in these regions are obtained separately and matched using the necessary algebraic continuity conditions including \(R_s\) and \(\Omega\). The continuous analogue of the Newton method for solving both the nonlinear differential and algebraic problems is used.