Role of vanishing complexity factor in generating spherically symmetric gravitationally decoupled solution for self-gravitating compact object

In this work, we study the role of the vanishing complexity factor in generating self-gravitating compact objects under gravitational decoupling technique in f(Q)-gravity theory. To tackle the problem, the gravitationally decoupled action for modified f(Q) gravity has been adopted in the form S=SQ+Sθ∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr {S}}={{\mathscr {S}}_{Q}}+{{\mathscr {S}}^{*}_{\theta }}$$\end{document}, where SQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr {S}}_Q$$\end{document} denotes the Lagrangian density of the fields which appears in the f(Q) theory while Sθ∗(=αSθ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathscr {S}}^{*}_{\theta } (=\alpha {\mathscr {S}}_{\theta }$$\end{document}, where α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} is just a coupling parameter which controls the deformation) describes the Lagrangian density for a new kind of gravitational sector which has not been included in f(Q) gravity. After that, we developed an important relation between gravitational potentials via a systematic approach (Contreras and Stuchlik in Eur Phys J C 82:706, 2022) using the vanishing complexity factor condition in the context of f(Q) theory. We have used the Buchdahl model along with the mimic-to-density constraints approach for generating the complexity-free anisotropic solution. The qualitative physical analysis has been done along with the mass-radius relation for different compact objects via M-R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M-R$$\end{document} curves to validate our solution. It is noticed that the coupling constant β1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _1$$\end{document} has a definite impact on constraining the mass and radii of the object that are shown in M-R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M-R$$\end{document} curves. The obtained results show that the compactness of the objects can be controlled by the coupling parameters.


Introduction
Einstein's General Relativity (GR), has been working as the basic platform to understand the mechanism of cosmological as well as astrophysical phenomena. However, recently it is under a few crucial issues, like the observational evidences of late-time accelerated phase of the Universe as evident from CMBR [1,2], high redshift supernovae [3], baryon acoustic oscillations [4], Planck data [5] and supernovae of type Ia [6,7]. To explain this observed dynamics of cosmology, the simplest prescription was to modify the field equations of GR by introducing the cosmological constant (Λ) which was originally suggested by Einstein himself to construct a static model of the universe [8,9].
A Λ-cosmological model of GR provides the best fit to the observed data, however, needs further assumption of Dark Matter [10,11] as well as Dark Energy. The latter one has a huge discrepancy in its value with the theoretical value of Λ [12,13]. Thus, there was an emergency need to modify GR via different extension schemes of the Einstein-Hilbert action which eventually gave birth to several alternative theories of gravity which are evidently quite different from the traditional dark energy depiction. Till date there are dozens of gravity theories available with numerous modification of GR which present a realistic perspective of the quantum mechanics and gravitational field at the high energy physical arena [14,15]. Therefore, in the initial phase modifying the Lagrangian density by a simplistic function f (R) in the Einstein-Hilbert action, an extended GR has been provided in f (R) gravity [16][17][18][19]. Later on, several other alternative gravity theories have been proposed by other scientists, such as the f (G) gravity [20,21], f (R, G) gravity [22], f (T ) gravity [23][24][25][26], Brans-Dicke gravity [27][28][29][30][31][32] and so on, where the symbols R, G and T are, respectively, the Ricci scalar, Gauss-Bonnet scalar and the torsion scalar.
On the other hand, another modified theory of gravity that has a deep-rooted historical background needs here to mention and discuss now. In the year 1928 Einstein [69] proposed a new theory in relation to his fundamental or classical GR. This is originally known as Teleparallelism which nowadays comes in the picture as teleparallel gravity, where an attempt was done by Einstein in connection to unification between the forces like gravity and electromagnetism. Basically, in the mathematical structure a spacetime was characterized by a curvature-free linear connection in conjunction with a metric tensor field and has received much attention in recent research on astrophysics as well as cosmology [70,71].
In the last few years, there has been greater emphasis on teleparallel theories as one of the alternate gravity theories. Although these theories most often work in a torsion-based environment, however the enthusiasm has shifted to symmetric teleparallelism, where gravity is effectively defined by nonmetricity rather than curvature or torsion. It is worth noting that it is the affine connection, rather than the physical manifold, which takes major role as the mechanism for mediating gravity. This is realized by the fact that, in GR, the curvature is a property of the connection and not of the manifold itself, and hence it can also be characterized by other connection properties like non-metricity. Symmetric teleparallel gravity is formulated in terms of a spacetime metric and Levi-Civita connection or its teleparallel formulation and offers an interesting geometric interpretation of gravitation. It exploits a connection that is torsion-free and curvature-free, but not metric-compatible to describe gravity. By using the so-called 'coincident gauge', [72,73] this connection can be minimized to a partial derivative. From a geometrical point of view, this means that the vectors on a manifold stay parallel over large distances [74]. The rest of the gravitational information is contained in the non-metricity contributions [72,[75][76][77] by requiring the connection to be torsion-free and the curvature to disappear.
The original GR being associated with spacetime curvature, a plausible explanation of the accelerated dynamics of the Universe was to entreat an exotic agent, known as dark energy (DE) [8,9] and for that, it is necessary to give relaxation to the torsion and nonmetricity. Therefore, as twin baby from the mother GR, one may get two equivalent prescriptions for an alternative gravity: (i) teleparallel equivalent of GR (TEGR) where torsion T has definite role [78][79][80] whereas (ii) symmetric teleparallel general relativity (STGR) where nonmetricity Q has a vital role [72,[81][82][83].
However, in the present work, we are employing STGR theory where the nonmetricity describes the gravity instead of the joint effect of the curvature and the torsion. Recently, f (Q) have gone through several successful applications in cosmology [84][85][86][87][88][89][90][91][92][93] as well as in astrophysics [94][95][96][97][98][99][100][101]. We would like to mention here an important work on f (Q) theory of gravity [102] where the authors considered several cosmological relevance by using the evolution equations and constraining specific functional forms in the Lagrangian L . Another work that can claim special mention is that Wang et al. [103] particularly have shown that there does not exist an exact Schwarzschild solution for the nontrivial f (Q) function. Therefore, out of the above-mentioned total scenario we are treating compact stars (specifically neutron and/or strange stars) under the teleparallel f (Q)-gravity as an alternative theory.
In the context of the compact stars, it is now evident that these are endpoints of normal stars due to tremendous gravi-tational inward pulling under the general relativistic effect of Einstein. They are usually clubbed as a white dwarf, neutron star and black hole where the following internal constitution of the matter content either of the formations into endpoint can take place. From the basic astrophysical point of view, to form a white dwarf the Chandrasekhar limit [142] should be maintained. However, all the ions of a normal star can be inverted into neutral particles to form a neutron star or beyond that the ultra-dense compact objects made up out of u, d and s type quark particles to form strange stars [143,144]. Eventually one may get a black hole from which nothing can be escaped due to its tremendous potential and hence the escape velocity becomes unattainable for a photon. This obviously makes the compact star invisible, i.e. black to an observer though the Hawking radiation [145] may work from the event horizon because of the relativistic quantum effects.
The outline of the present study is as follows. In Sect. 2 we have provided field equations (i.e. in 2.1 set of field equations in pure f (Q)-gravity whereas in 2.2 set of equations for new source θ ν ). In Sect. 3 qualitative physical analysis has been performed, e.g. physical behavior of density and pressures (Sect. 3.1), energy conditions (Sect. 3.2), causality and sound velocity (Sect. 3.3), adiabatic stability criterion (Sect. 3.4). The impact of constant β 1 on constraining the mass via equi-mass diagram is discussed in Sect. 3.5 while the impact of deformation through constant α as well as β 1 on mass-radii relations via M − R curves are presented in Sect. 3.6. However, the impact of β 1 on physical parameters such as density and pressures is present in Sect. 3.7. We have devoted Sect. 4 for the overall summary and a few final concluding remarks on the investigation. Therefore, the constraints on torsionless and curvatureless provide the affine-connection as Furthermore, one can opt for a special coordinate system, e.g. coincident gauge, under which Γ λ ν = 0. The nonmetricity Eq. (4) then gets the following form The above form does help immensely to perform the calculation as the fundamental variable being here only the metric function. However, it is to note that the action in this case no longer remains a diffeomorphism invariant unlike the standard GR [93]. Therefore, as an easy path one may undergo the covariant formulation of f (Q) gravity. The affine connection as expressed in Eq. (12) being purely inertial, we can perform the covariant formulation by determining the affine connection when gravity is absent [146].
In the ongoing study, our interest is to find the gravitationally decoupled solutions under the theory of f (Q) gravity for self-bound compact stellar objects. For this we have considered a Schwarzschild-type static spherically symmetric metric as where Φ(r ) and μ(r ), respectively, describe the metric potentials which are mainly functions of r only.
Furthermore, the formula for Q in connection to the spacetime (14), is determined as On the other hand, we also consider that the interior of the compact star under f (Q) gravity is composed of anisotropic matter distributions. In this situation, one can provide the original source T ν in the context of f (Q) theory as (16) where energy density (ρ), radial pressure ( p r ) and tangential pressure ( p t ) are the main matter variables in the pure f (Q)-gravity wheres η ν and χ ν satisfies further definition as: u 0 u 0 = −1 and χ 1 χ 1 = 1.
The components for effective energy-momentum tensor are defined as As a result, the θ -sector components can be provided as Using above we get and effective anisotropy is defined as In the above expression it is clear that the effective anisotropy is the combination of two anisotropies corresponding to T ν and θ ν . Here the anisotropy Δ θ has been invoked due to gravitational decoupling which may enhance the effective anisotropy.
Due to the anisotropic fluid (17), the independent components of the equations of motion (8) under f (Q) gravity can be given as where f Q (Q) = ∂ f (Q) ∂ Q . Furthermore, according to Wang et al. [149] the exact Schwarzschild anti-de Sitter solution exists for the constraint f Q Q = 0. On the other hand, functional form derived from f Q Q = 0 is also compatible with the line element (14) and the coincident gauge [146]. Therefore, we must derive the functional form of f (Q) by considering only f Q Q to be zero in order to obtain the solution for self-gravitating compact star models. Hence where β 1 and β 2 are constants. By inserting of Eqs. (14) and (25), the Eqs. (21)-(23) give the following set of equations of motion The linear combination of equations (26)- (28) gives, Interestingly, the above-mentioned Eq. (29) is called as the modified TOV equation [147,148] under f (Q)-gravity theory [149,150] satisfying the relation (25) which is equivalent to the conservation equation in classical general relativity. Now, our motivation in the present work is to obtain an exact solution of equations of motion (26)-(28) under the vanishing complexity formalism in f (Q)-gravity context. For this purpose, we obtain the condition of vanishing complexity in f (Q)-gravity using the definition of complexity proposed by Herrera [122][123][124][125]. Based on his definition, the complexity factor (Y Q T F ) is defined for the system (26)- (28) as On inserting of P eff r , P eff t and ρ eff in above equation, we Then Since β 1 can not be zero therefore the second factor must be zero, which gives The simplified form of the above differential equation can be written after performing the integration as Lastly, we derive the bridge equation between Φ and μ as where A and B are two constants of integration. It is very interesting to note that the above relation (35) in f (Q)gravity theory is similar to the relation derived by Contreras and Stuchlik [151] in the context of Einstein's GR. As mentioned above, our primary objective was to get an exact solution to the field equations (26)-(28) describing a compact star beyond the f (Q)-gravity theory. Therefore, we use the well-known technique of gravitational decoupling through the minimal geometric deformation (MGD) approach via a specific transformation where η(r ) denotes the deformation function along the temporal component of spacetime while the deformation function Ψ (r ) is along the radial metric component.
In the MGD case Ψ (r ) = 0 while defromation function η(r ) = 0 and which we have shown in Fig. 2. Moreover, the MGD technique basically divides the decoupled system (22)- (24) into two subsystems under f (Q)-gravity, where equations of motion for T ν is denoted as first set while the second set is for the new matter source θ ν as, will be given in the next two subsections.

Systemic field equations in f (Q)-gravity
We have the following set of equations and according to Eq. (28) one can get which is nothing but a TOV equation for the system (31)- (34) and its solution can be provided by the spacetime The linear combination of equations of motion (43)-(45) yields From the above equation, one important point may be noted that it is not required any energy exchange among the sources for decoupling.
Using Eq. (37), the condition (35) reduce to Now we move on to the solution of field equations for both systems. We have nine unknowns and six independent equations along with vanishing complexity condition (36). Therefore, we need two auxiliary conditions to solve the system completely. For this purpose, we consider well-behaved ansatz for the radial metric function W (r ) in pure f (Q)gravity [152] connected to the spacetime (42) and mimicking of density constraint approach (ρ = ρ θ ) as where L and N are constants with dimension length −2 and length −4 , respectively (Fig. 1). Now the differential Eq. (49) along with Eq. (48) leads the solution for deformation function Ψ (r ) as Now the new form of the gravitational potential e μ(r ) is given by e μ(r ) = 6β 1 Nr 2 + 1 6β 1 1 + (α + 1)Lr 2 − α N r 2 + αβ 2 r 2 Nr 2 + 1 .
Before proceeding further, it is highlighted that the most suitable exterior spacetime under the functional form Eq. (25) in the f (Q)-gravity theory can be given by exact Schwarzschild (Anti-) de Sitter spacetime. Therefore, the solution's constant must be calculated by joining of the deformed interior spacetime with the Schwarzschild (Anti-) de Sitter exterior spacetime at the boundary surface r = R. The Schwarzschild (Anti-) de Sitter spacetime and the suitable boundary conditions are: Now use the Darmois-Israel boundary conditions [153,154] in order to express the first and second fundamentals mathematically as where M is the total mass and Λ the cosmological constant. It is noted that cosmological constant Λ can be written in the form of constant β 1 and β 2 as Λ = β 2 /2β 1 . As we know that the value of the cosmological constant Λ in the present universe according to the current observational evidence is approx 10 −46 /km 2 . Therefore, there is no effect of Λ on the current stellar models and then it can be considered zero. Due to this assumption, we get β 2 = 0. By substituting of Eq. (51) into the condition (36) and integrate, we get where

Qualitative physical analysis
In this Section, we present the analysis of the gravitationally decoupled anisotropic solution to test the regularity conditions of the obtained models in f (Q)-gravity theory. Note that we have checked the model behavior for α = 0.16, 0.18, 0.20, and 0.22. It is obvious that when α = 0, then the pure f (Q)-gravity model is recovered. It is to be mentioned that the solution is not physically valid when α < 0.16 (due to increasing tangential pressure within the objects) as well as α > 0.22 (due to not preserving the causality condition, especially radial velocity will be more than the speed of light). Therefore, the following range of the decoupling constant 0.16 ≤ α ≤ 0.22 yields a physically viable model in f (Q)-gravity theory.

Physical behavior of density and pressures
In Figs. 2, 3 and 4, the variations of effective energy density, effective radial and tangential pressures with respect to radial distance r for different α are shown. It is clear that energy density and pressures are monotonically decreasing for increasing r for each α as expected for a physically viable solution. On the other hand, the central values of density and pressures are increasing when α increases which implies that a more dense compact object can be obtained in the context of gravitational decoupling. This can be also observed from the increasing value of mass-radius ratios with increasing α given in Table 1.

Energy conditions
For investigating the structural properties of our solution, we emphasize the energy conditions (EC) which are really predictable that all of the standard models must satisfy the energy conditions in all kinds of gravity. The necessary energy conditions are groups of inequalities subjected to the energy-momentum tensor of matter. The most viable conditions are Weak energy condition (WEC) (ρ eff + P eff r ≥ 0 and ρ eff + P eff t ≥ 0), Strong energy condition (SEC) (ρ eff + P eff r + P eff t ≥ 0), and Dominant energy condition (DEC) (ρ eff − P eff r ≥ 0 and ρ eff − P eff t ≥ 0). To verify all the above conditions, we first look at Figs. 2, 3, 4, the pressure and density are positive throughout the stellar object. Then WEC and SEC are already satisfying. Now we need to check only DEC. For this purpose, we plotted Fig. 6 to show the variation of quantities ρ eff − P eff r and ρ eff − P eff t within the stellar object which clearly show that both quantities are positive at each point of the object. This implies that our solution also satisfies the dominant energy condition (Fig. 5).   for all decoupling constant α. This implies that the causality conditions preserves throughout the model.

Adiabatic stability criterion
The adiabatic index (γ ) plays an important role in investigating the dynamical stability of a compact object. This adiabatic index Γ for anisotropic matter distribution for the decoupled system is given by where the derivation is accomplished at constant entropy S while d P eff r dρ eff S denote the velocity of sound. This shows that the sound velocity is an important measure connected to the adiabatic index, e.g. the adiabatic index for the Schwarzschild  [155] show that the central value of the adiabatic index should be more than 4/3 for a stable polytropic star which will obviously depend on the amount of central velocity as well. Furthermore, based on current investigations, it is argued that the range of Γ r should be 2 to 4 in massive neutron stars [156]. In Fig. 9 we plot adiabatic index Γ r as a function of radial distance r for different α. The figure shows that the adiabatic index Γ r > 4/3 when α ≥ 0.19. In this connection, Moustakidis [157] proposed a stronger condition on the stability of a selfgravitating system which is Γ > Γ crit , where Γ crit denotes  Table 1 contains the values of Γ r for different values of α, which shows that Γ > Γ crit satisfies when α ≈ 0.22. Finally, we conclude that the gravitational decoupling technique may help in stabilizing the stellar configuration. Fig. 11 The M − R curve shows the mass-radius relation for different β 1 with fixed α = 0.14

Fig. 12
The M − R curve shows the mass-radius relation for different α with fixed β 1 = 1.2 km 2 3.5 Impact of coupling parameter β 1 in constraining the mass via equi-mass diagrams We discuss here on the impact of coupling parameters β 1 on constraining the mass of the object via equi-mass diagram. In Fig. 10 we have shown the equi-mass diagram on R − β 1 plane. One can notice from this figure that the pattern mass distribution on R −β 1 plane shows that the mass is increasing in this situation when β 1 increases with fixed R or R increases with fixed β 1 . But the effect of β 1 can be noticed clearly near the boundary. Table 2 The predicted radii of compact stars LMC X-4, PSR J 1614+2230, PSR J 0740+6620, and GW190814 for the different β 1 with fixed α = 0.14 and different α with fixed β 1 = 1. In this Section, we focus on determining the masses of four neutron star candidates with masses ranging from 1.29 +0.05 −0.05 to 2.14 +0.2 −0.17 . Furthermore, we also discuss the estimation of massive neutron stars with masses in the range of 2.5-2.67 M came from the GW190814 event. Observations of the masses and their radii for these four neutron stars can be determined from Figs. 11 and 12 for the values of different β 1 with fixed α and different α with fixed β 1 , respectively.
It is observed from Fig. 11 that our complexity-free solution expects the presence of low and high-mass neutron stars observed by authors [158][159][160][161]. As β 1 increases with fixed α, our solution possesses the higher mass neutron star with larger radii.  Table 2 includes the predicted radii of different compact objects for different β 1 and α. Finally, we can say that higher mass objects can be expected in f (Q)gravity theory.
It has been noted recently by Rincón et al. [162,163] that the M − R profile can be obtained in the context of the vanishing complexity factor and it was also confirmed by the group that the solutions thus obtained exhibit more compact objects than GR. In this regard, we have extended our study to testify self-gravitating compact objects via M − R profile in f (Q)-gravity theory under complexity formalism. Based on these arguments and the above analysis, we conclude that we get more compact objects when β 1 is decreasing while higher β 1 gives more massive objects. It is obvious that if β 1 = 1, then f (Q)-gravity theory will be equivalent to GR. The effective energy density distribution on the r − β 1 plane via the contour diagram is given in Fig. 13. It clearly shows that the density is decreasing monotonically when we move toward the boundary for any fixed value of β 1 ∈ [0.8, 1.4] and attains its minimum value at the boundary. On the other hand, when we increase β 1 , then the magnitude of density increase as well, and its maximum value is obtained at the core of the star. Figures 14 and 15 show the effective radial and tangential pressures ( p eff r and p eff t ) variations on the r − β 1 plane. For fixing β 1 ∈ [0.8, 1.4], both pressures are decreasing mono-tonically with increasing radial distance r while magnitude of p eff r and p eff t increases when β 1 increases. The maximum pressure is observed at the center for β 1 ∈ [0.8, 1.2] as expected in our analysis.
Therefore, we may assert that the constant β 1 play important roles in constraining the energy density and pressure in f (Q)-gravity theory.

Concluding remarks
In this article, we have successfully explored the role of the vanishing complexity factor in generating self-gravitating compact objects in f (Q)-gravity theory. Here, we have mainly adopted the gravitationally decoupled action for modified f (Q) gravity in the form S = S Q + αS θ , which is equivalent to S = S Q + S * θ , where S Q designates the Lagrangian density of the fields occurring in the f (Q) gravity theory while S * θ (= αS θ ) represents the Lagrangian density of a new gravitational sector, not described by f (Q) gravity. Since the new sector αθ ν is equivalent to θ * ν because α is just a control parameter. What's more, using a systematic approach suggested by Contreras and Stuchlik [151] and the vanishing complexity factor condition stated by Herrera [122], we found a significant relationship between gravitational potentials. The complexity-free anisotropic solution is then generated by combining the Buchdahl model with the mimic-to-density constraints approach.
The eventual model was carefully examined to assess the acceptance of the gravitationally decoupled anisotropic solution via satisfying stringent regularity and stability constraints in f (Q)-gravity in the presence of deformation function via decoupling parameter α. It is interesting to note that the pure f (Q)-gravity theory can be achieved in the absence of deformation, whereas we get a physically viable solution beyond the f (Q)-gravity theory in the presence of deformation function. In the presence of deformation function, the physical behavior of the effective state variables (ρ eff , P eff r , P eff t ), effective anisotropy (Δ eff ) and the effective energy conditions namely, WEC, SEC, and DEC are analyzed within the stellar object. It is revealed that the limit necessary for the physical viability of the models beyond the f (Q)-gravity theory is satisfied by our anisotropic solution. We have also tested the stability of the resulting solution using the cracking technique and the adiabatic index, and we have discovered that the solution upholds the model-wide stability criterion.
We have then analyzed and made a few notable observations about the impact of coupling parameter β 1 in constraining the mass via equi-mass diagrams. We have also shown the effect of the presence of new sector and constant β 1 in constraining the Radii for different objects via M − R curves along with the impact of constant β 1 on density and pressures. We specifically observed the impact of β 1 on constraining the mass of the object via equi-mass diagrams shown in Fig. 10. We noticed from the mass distribution on R − β 1 plane that if β 1 increases with fixed R or R increases with fixed β 1 , the mass magnitude increases. Furthermore, it is evident from the pattern of the mass distribution on R−β 1 plane that the increment in mass can be observed when R ≥ 6 and high mass can only be achieved for high values of β 1 at R = 10 km.
With an emphasis on the masses of four neutron star candidates with masses ranging from 1.29 +0.05 −0.05 to 2.14 +0.2 −0.17 and the estimation of massive neutron stars with masses in the range of 2.5 to 2.67 M , came from the GW190814 event, we have discussed the impact of coupling parameters in constraining the Radii via M − R curves for different values of α with a fixed β 1 and vise-versa. On the one hand, we notice that our solution has the higher mass neutron star with larger radii when β 1 increases. For β 1 = 1.4, our model predicts the massive neutron star with masses in the range of 2.5 to 2.67 with a radius of 14.23 +0. 21 −0.02 km, which is a secondary component of the GW190814 event. However, the masses and radii for other neutron stars are highlighted that secondary component of the GW190814 is excluded for low values of β 1 = 0.8. On the other hand, we observe that the mass and radii both are increasing when α increases with fixed β 1 . It should be highlighted that in the absence of deformation, the secondary component of the GW190814 event with masses 2.5 − 2.67 is ruled out. However, in presence of deformation with α ≥ 0.14, we forecast the existence of a secondary component of the GW190814 event with mass 2.5. Moreover, the constant β 1 has more impact on the radii than the constant α (See the predicted radii in Table 2 for more details). Interestingly, we have clearly shown the impact of the constant β 1 on the effective energy density distribution and the effective radial and tangential pressure variations on the r − β 1 plane via the contour diagram. We notice that all the physical variables are decreasing monotonically with increasing radial distance r while their magnitude increases when β 1 increases and their maximum values are obtained at the core of the star. Therefore, we can also state that the constant β 1 play important roles in constraining the energy density and pressure in f (Q)-gravity theory.
Finally, it can be claimed that the role of the vanishing complexity factor in generating self-gravitating compact objects under f (Q)-gravity theory is viable for constructing astrophysical models that are congruent with experimental events.
its Visiting Research Associateship Programme as well as the facilities under ICARD, Pune at CCASS, GLA University, Mathura.

Data Availability Statement
No new data were generated or analyzed in support of this research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.