Finding new solutions of physical interest to the Einstein field equations is not an easy task due to the highly nonlinear behavior of its equations. In this regard, in (2017) was proposed the Gravitational Decoupling Method [1]}, which represents an easy algorithm to decouple gravitational sources in General Relativity. This algorithm involves a Minimal Geometric Deformation (MGD) to the metric tensor together with a decoupling of sources. The method was described in reference [2]} as follows: “given two gravitational sources: a source A and an extra source B, standard Einstein’s equations are first solved for A, and then a simpler set of quasi-Einstein equations are solved for B. Finally, the two solutions can be combined in order to derive the complete solution for the total system". The source A represents to the seed sector and the source B represents to the extra sector. Since its appearance, this method has served to find several new solution of physical interest, as for example well behaved solutions that could represent stellar distributions [2]}, [4]}, [5]}, [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, black hole solutions [12]}, [13]}, [14]}, \(f({\mathcal {G}})\)  gravity [15]}, Einstein Klein Gordon System [16]}, Pure Lovelock gravity [17]}, \(f(R,T)\)  gravity [18]}, \(f(R)\)  gravity [19]}. See other applications [20]}, [21]}, [22]}, [23]}, [24]}, [25]}, [26]}, [27]}, [28]}, [29]}, [30]}, [31]}, [32]}, [33]}
