Stellar Filaments in Self-Interacting Brans-Dicke Gravity

This paper is devoted to study cylindrically symmetric stellar filaments in self-interacting Brans-Dicke gravity. For this purpose, we construct polytropic filamentary models through generalized Lane-Emden equation in Newtonian regime. The resulting models depend upon the values of cosmological constant (due to scalar field) along with polytropic index and represent a generalization of the corresponding models in general relativity. We also investigate fragmentation of filaments by exploring the radial oscillations through stability analysis. This stability criteria depends only upon the adiabatic index.


Introduction
Filamentary structures have important implications in structure formation of the universe. These stellar configurations are omnipresent in the universe at various scales. On cosmological scales, the cosmic filament is associated with the cosmic web [1][2][3][4] where it behaves like a bridge connecting various dense regions (galaxies). On small scales, filaments are the features of the interstellar medium and instabilities within these filaments form a dense medium which turns into a star [5]. The stellar filamentary structures appear in a variety of astronomical contexts; their dynamical analysis can be considered as an effective tool to understand the behavior of galactic as well as interstellar structures in the universe.
It happens very often that the geometries of the proposed models are taken to be very simple or ideal, partly because of reasons like mathematical simplicity which may avoid complexities in the analysis, to obtain some results of physical interest, and to provide an ideal model which can act as a foundation and can be modified according to the required astrophysical as well as cosmological applications. The cylindrically symmetric configuration is one of the ideal a e-mail: msharif.math@pu.edu.pk b e-mail: rubab.manzoor@umt.edu.pk cases, which has widely been used to represent filamentary structures in the universe. In 1953, Chandrasekhar and Fermi [6] described the dynamics of a cylinder filled with a homogeneous and incompressible fluid. Stodolkiewicz [7] developed the magneto-hydrodynamics equilibrium of isothermal cylindrical filaments. Ostriker [8,9] studied the cylindrically symmetric filamentary structure with a polytropic equation of state and generalized the dynamical analysis of homogeneous compressible filaments. Afterward, many researchers explored cylindrical filamentary structures both analytically as well as numerically [10][11][12][13][14]. Recently, Breysse et al. [15] described polytropic filamentary structures through a stability analysis of the cylindrically symmetric self-gravitating fluid.
Modified theories of gravity constructed by modifying the Einstein-Hilbert action are considered to be candidates for the accelerating agent (dark energy) of the expanding universe [16][17][18][19]. There is a large body of literature [20][21][22][23][24][25] which has explored the dynamics of stellar structures in modified gravity to reveal the modification hidden in the structure formation of the universe. Brans-Dicke (BD) gravity, being a natural generalization of general relativity (GR) [26,27], is one of the fascinating examples of modified gravity, which is considered as a solution of various cosmic problems. One of the main features of this theory is that the gravitational force is described by a massless scalar field φ along with a curvature part (Ricci scalar). It also contains a coupling constant ω BD which serves as a tunable parameter to adjust the required results. This theory provides suitable solutions of many cosmic issues but remains unable to describe the "graceful exist" problem of old inflationary cosmology. The inflationary model of this theory is valid for a specific value of the coupling parameter ω BD ≤ 25 [28] which contradicts observational data [29]. Moreover, the observational limit of ω BD at small scale (weak field) [30,31] is inconsistent with those calculated at large scale [32]. These issues were resolved in self-interacting BD (SBD) gravity, which is developed by introducing a massive poten-tial function V (φ) in the Jordan framework of BD gravity [33]. The potential function represents a massive scalar field term which describes the potential of the scalar field (φ) or the potential of the energy density related to the scalar field in the system. This theory provides consistency between the weak field and the strong field regimes [34]. Recent literature indicates the dynamics of SBD gravity in many cosmic problems [35][36][37][38][39][40]. In recent papers [41][42][43], we have explored spherically as well as cylindrically symmetric selfgravitating fluids in SBD gravity and found some interesting results. We have also investigated hydrodynamics and oscillations of spherically symmetric gaseous distributions in the post-Newtonian approximations of this theory [44,45].
In this paper, we study filamentary structures of stellar systems in SBD theory. We construct cylindrical polytropic models of the filaments in the Newtonian (N) approximation and explore the fragmentation of filamentary structures through radial oscillations of self-gravitating fluids. The paper is organized as follows. Section 2 formulates SBD as well as dynamical equations in the N limit. In Sect. 3, we construct cylindrical polytropic filamentary models in N approximation. Section 4 is devoted to a discussion of radial oscillations of the filaments through a stability analysis. Finally, Sect. 5 summarizes the results.

Self-interacting Brans-Dicke gravity and dynamical equations
The action of SBD gravity [33] is given by where κ 2 = 8π G c 2 and L m shows the matter contribution. The variation of this action with respect to the scalar field (φ) and the metric tensor (g μν ) provides the SBD equations as follows: where T μν represents the energy-momentum tensor of matter, T = g μν T μν , and stands for the d'Alembertian operator. Equations (2) and (3) provide the field equations and evolution of the scalar field, respectively. We assume a matter contribution in the form of a perfect fluid, which can be compatible with the N regime, where ρ, p, u μ stand for density, pressure, and four velocity, respectively.

Newtonian approximation
The weak field approximated solutions of any relativistic theory describe the order of deviations or perturbation of the local system from its homogeneous and isotropic vacuum background. The N and parameterized post-Newtonian limits are widely used weak field approximated solutions that are derived by a Taylor expansion of the metric functions as follows: Here η μν indicates the Minkowski metric (describing a homogeneous and isotropic vacuum background of g μν ), h μν shows the deviation of g μν from its background values (η μν ), i, j = 1, 2, 3 and the superscripts (2), (3), and (4) represent the order of approximation (c −2 ), (c −3 ) as well as (c −4 ). The N limits require the information of g 00 ∼ η 00 + h (2) 00 , g i j ∼ η i j , whereas the parameterized post-Newtonian corrections use the approximations g 00 ∼ η 00 + h (2) 00 + h (4) 00 , g 0i ∼ h (3) oi , and g i j ∼ η i j + h (2) i j . Thus, the N limits of any system can be directly obtained from its known parameterized post-Newtonian approximations.
In order to discuss polytropic geometry in SBD gravity and check the compatibility of our results with the analysis of GR [15], we approximate the system in the N limits. For this purpose, we evaluate the N limits of SBD gravity from its known post-Newtonian approximate solutions. The parameterized post-Newtonian approximation of SBD (massive BD gravity) solutions has been evaluated by using the following Taylor expansion of the metric and the dynamical scalar field [46]: Here φ 0 describes the dynamical scalar field as a function of background cosmic time t 0 (which varies slowly with respect to cosmic time t 0 ), V 0 = V (φ 0 ) is the potential function of the scalar field at t 0 , and ϕ(t, x) represents local deviation of the scalar field from φ 0 . The lowest-order parameterized post-Newtonian corrections (O(c −2 )) of the SBD solutions are given by Here U = G eff M r (M is the Newtonian mass of the sun) shows the effective gravitational potential determined by the Poisson equation where The term γ BD represents the parameterized post-Newtonian parameter given by where m 0 represents mass of the massive scalar field with constraint (m 0 >> 1 r ) (wherer represents the scale of the experiment or observation testing the field). When the background value of this mass is very small (m 0 << 1 r ) the SBD system reduces to simple BD gravity (the massive scalar field reduces to a massless scalar field) having the post-Newtonian parameter That is why the BD theory (massless scalar field) is consistent with solar system constraints of the Cassini mission for ω BD > 40,000. However, for massive BD gravity (SBD gravity), the dynamics of the spatial part of φ is frozen on the solar system scales through the potential function and all values of ω BD are observationally acceptable [47][48][49]. The term V 0 r 2 6φ 0 c 2 = BD r 2 3c 2 describes the cosmological constant term which depends upon the potential of the scalar field. In order to be consistent with observational data ranging from the solar system to clusters of galaxies the contribution due to the scalar density should be very small and the following constraint must be satisfied: Here L indicates the length scale equal to or greater than the solar system. The N limits of any scalar-tensor theory require information of g 00 ∼ η 00 + h (2) 00 , g i j ∼ δ i j and φ ∼ φ 0 . Thus, from the parameterized post-Newtonian analysis of SBD theory the obtained N approximated SBD solutions are evaluated as From now on, we use h (2) 00 = h 00 for the sake of convenience. It may be noticed that, by applying the limits (m 0 << 1 r ) and V 0 φ 0 → 0, the above defined N approximations of SBD gravity can be shifted to N limits of BD solutions. Similarly, in the limits (m 0 << 1 r ) , V 0 φ 0 → 0 and ω BD → ∞, the obtained approximated system can be converted into the GR case.

Dynamical equations
The polytropic geometry of any configuration is based upon two dynamical equations, the Poisson equation and the equation of motion of the respective system [15]. Here, we calculate the generalized form of both dynamical equations for SBD gravity. Equation (2) can be rewritten as The temporal component of this equation in N approximation is given by where the contributions due toφ 0 andφ 0 are neglected because the term φ 0 behaves as almost a constant. The expansion of Ricci tensor in the N regime is approximated by Comparing these two values, we obtain the generalized Poisson equation as follows: where h 00 represents the gravitational potential due to matter as well as the massive scalar field. The generalized Euler equations can be obtained by using T μν ;ν = 0, whose time component yields the continuity equation, where v i (i = 1, 2, 3) stands for the components of the velocity. The spatial components give the generalized form of the equation of motion as follows:

Cylindrical polytropes
In order to discuss cylindrical polytropic filaments in SBD gravity, we consider the equilibrium configuration and the standard polytropic equation of state: where K is a constant and γ = n+1 n represents the polytropic exponent with n as a polytropic index [15,50]. In hydrostatic equilibrium, the cylindrical configurations of Eqs. (12) and (14) turn out to be 1 r Integration of Eq. (17) by using (15) provides Using a binomial expansion on g − 1 2(n+1) 00 and approximating up to O(c −2 ), we have Using Eqs. (19) into (18), we obtain the relation between the density and potential configurations as follows: Equations (16), (20), and the cosmological constant term provide a differential equation for a gravitational potential h 00 as follows: The above equation can be converted to a polytropic equation (Lane-Emden equation) by considering the following dimensionless variables: where h 00(c) and ρ c represent the potential as well as the density at the center of the cylinder (r = 0) [15,50]. With these assumptions, Eq. (21) turns out to be This is a modified form of the original Lane-Emden equation describing spherical polytropic model [50] and generalized form of the modified Lane-Emden equation that describes the cylindrically symmetric polytropic filament in GR [15]. Here β(s) is the Lane-Emden function describing the equilibrium potential and A is considered as a cosmological constant term due to the presence of the scalar field. This is a second order differential equation satisfying the following boundary conditions: at the center (r = 0) we have s = 0, h 00 ≈ h 00(c) , ρ ≈ ρ c , β ≈ 1 and dβ ds s=0 = 0.
The surface of the polytopic filament is represented by a specific value s = S for which the density ρ as well as the Lane-Emden function (β) becomes zero. Equation (23) can be solved analytically for different values of n. Some possible analytical solutions are given as follows: where J 0 (s) is a Bessel function of zeroth order. It is noted that the analytical results depend upon the values of s as well as A. We can also solve Eq. (23) 2 show polytropic models for n = 0 with A = 0, ±1/2, ±1, ±3. It is found that the polytropic filaments for A 1/2 have finite radii, i.e., after taking finite values of s, the polytropic function β becomes zero. For A > 1/2, the polytropic models have infinite radii. Figures  3 and 4 show that, for n = 1, we have finite radii models for A 0 and infinite radii models for A > 0. Figures 5 and 6 represent polytropic filaments for n = 3 where A 0 represents finite radii, while A > 0 gives infinite radii models. It can be noticed that under the constraint (m 0 << 1 r ) and A = 0, Eq. (23) (derived for SBD gravity) is converted into a non-linear homogeneous Lane-Emden type differential equation, which can be used to discuss polytropic filaments in BD gravity. Thus, numerical results obtained for A = 0 and n = 0, 1, 2, 3, 5 can indicate polytropic filaments in BD gravity. It is mentioned here that the polytropic models described by GR depend only upon the values of the polytropic index n. In the spherical case, the polytropic stars for n < 5 have finite radii, while n 5 provides infinite radii models [51]. In the cylindrical case, there exist models with finite radii for 0 n < ∞ [7]. Since we have derived filamentary polytopic models present in the universe from small scales to cosmic scales, the astrophysical length scale of our evaluated models are appropriate for both stellar size and cosmic large scales.

Radial oscillations of polytropic filaments
In this section, we discuss the stability as well as normal modes of radially oscillating filaments about their equilibrium configuration. For this purpose, we assume that initially the system is in complete hydrostatic equilibrium. Afterward, the stellar body starts oscillating (in the radial direction) and becomes perturbed [15,51].

Stability analysis
In order to explore radially oscillating filamentary structures for different matter distributions in SBD gravity, we assume Lagrangian coordinates (μ, t), where μ represents the mass per unit length and t stands for time. Let r = r (μ, t) be the distance of the mass from the center of the cylinder such that μ = πr 2 ρ satisfies the following relations [51]: .
To study the stability criteria for an oscillating cylindrical filament, we consider a thin shell of mass element dμ (per unit length). The shell faces a force F g per unit area from the gravitational field of where g = G (eff) μ r + BD 3 represents gravitational potential due to the matter and scalar field. The pressure gradient forces exert a force F P per unit area of the shell given by of eigenvalues (ω n ) satisfying the property ω 2 n+1 > ω 2 n . For each eigenvalue, there is a corresponding eigenfunctionr n , which represents an amplitude of oscillations with a number of n nodes in the range 0 < r 0 < R 0 . The lowest-order eigenfunctionr 0 shows no node and is known as the fundamental amplitude.
It is mentioned here that stability, instability, and marginal stability of any oscillating model depend upon the behavior of its frequency. If the frequency of the model is real, the perturbation is purely oscillatory with constant amplitude providing a dynamically stable equilibrium, while imaginary values of the frequency lead to periodic oscillations in which the amplitude increases exponentially in time and the system becomes dynamically unstable. If the frequency tends to zero, the system becomes marginally stable (neither stable nor unstable), i.e., the model will expand and contract with the homologous property [51].
In order to evaluate the stability criteria for a radial perturbation, we solve Eq. (40) for the fundamental mode with boundary condition p 0 (R 0 ) = 0. We insert the eigenfunctionr 0 in the SL equation and integrate over 0 < r 0 < R 0 ; it follows that Sincer 0 is finite everywhere, the first term in the above equation vanishes, yielding This equation provides the fundamental frequency of the filaments, which leads to the stability criterion of the system. If ω 2 0 > 0 then ω 2 n > ω 2 0 > 0 for all nodes n > 0 and the frequency (±ω n ) becomes real for all values of n, leading to dynamically stable filaments. Similarly, ω 2 0 < 0 ⇒ ω 2 n < 0 for a finite number of nodes and the frequency becomes imaginary, providing unstable configurations. Otherwise ω 2 0 = 0 gives a marginally stable configuration.
Since the fundamental amplitude represents no node, it has the same sign all over the cylindrical filament and hence the two integrals in Eq. (42) have the same sign. Thus, the stability criterion depends upon sign ω 2 0 = sign 2( (ad) − 1). If ω 2 0 > 0, and the corresponding model becomes stable. For ω 2 n = 0, and the system is marginally stable. When ω 2 0 < 0, and the system becomes dynamically unstable. Thus, the criterion of stability for the filamentary structures in the N limit of SBD gravity depends only upon the stiffness of the fluid ( (ad) ). The stability conditions are independent from the behaviors of dynamical variables related to the metric, the coupling constant (ω BD ), the scalar field, and the potential of the scalar field.
It is mentioned here that in GR, the stability of spherical as well as cylindrical polytropic models depends only upon the stiffness parameter. In the spherical case, (ad) < 4 3 leads to an unstable model, while cylindrically symmetric filaments remain stable for (ad) > 1 [15,50]. In the case of BD gravity, the spherically symmetric polytropic model remains unstable for (ad) > 4 3 [52], while cylindrically symmetric polytropic models are not discussed. Thus, our obtained results are consistent with GR in the N regime.

Normal modes of radial oscillations
The behavior of normal modes for polytropic filaments can be described by using polytropic quantities in the SL problem. We use Eqs. (15), (20), (22) in (40), and it follows that In order to solve this equation forr , the required boundary conditions are as follows: at the center of cylinder, we have dr ds s=0 = 0, at the boundary of polytrope s = S, we have p 0 ρ 0 <<< 1, and the values ofr remain finite. Thus, Eqs. (36) and (37) yield [ω 2 bS + g 0 (2 − 2 (ad) )]r (S) − g 0 (ad) dr ds s=S = 0.
Equation (46) along with boundary conditions provides different patterns of normal modes of oscillating polytropic filaments in SBD theory. It is mentioned here that the behavior of these normal modes depends upon the frequency, the central potential, the polytropic index, the adiabatic exponent, and the gravitational potential term due to the massive scalar field (g 0 ), i.e., for different values of these parameters, we have different modes of radial oscillations of polytropic filaments.

Final remarks
This paper investigates cylindrically symmetric filamentary structures in the N limit of SBD gravity. We have formulated a generalized form of the Lane-Emden equation in N approximations and obtained polytropic filament models analytically for n = 0, 1 as well as numerically for n = 0, 1, 3, 5.
We have found that the behavior of these models depend upon the values of polytropic index as well as the cosmological constant term A (due to the scalar field). For n = 0, the models have finite radii for A 1/2, otherwise they have infinite radii. For n = 1, 3, 5, we have approximated that A 0 represents finite radii polytropic filaments while A > 0 gives infinite radii models. We have also found that within m 0 << 1 r and A = 0, the models defined for A = 0 and n = 0, 1, 2, 3, 5 are same for BD gravity.
In order to study fragmentation of the filamentary structures, we have investigated the stability of radial oscillations of polytropic filaments. It is found that the stability criterion of cylindrical filaments in SBD gravity depends only upon the adiabatic index ( ad ). Generally in the weak field approximation of GR and BD gravity, the stability criterion depends upon the adiabatic index as well as on the dynamical variables related to matter and scalar field distributions. But in the polytropic case, the adiabatic index is responsible for the stability of the function. Finally, we have discussed possible normal modes of radial oscillations of polytropic filaments. It turns out that different values of parameters lead to different modes of radially oscillating filaments. It is found that the weak field approximation of SBD gravity is consistent with observations for all arbitrary values ω BD and hence our obtained results are valid for all arbitrary values of ω BD . It is interesting to mention here that our results provide a generalized form of cylindrical filament polytropic models of GR theory.