Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 904 https://internationalpubls.com Anisotropic Charged Spheres in Higher Dimensional Rosen's Bimetric Relativity Theory R. T. Katre1*, R. K. Agrawal2, S. D. Tade3 1Department of Mathematics, Bajaj College of Science, Wardha, India 2Department of Mathematics, Hislop College, Nagpur, India 3Ex-Department of Mathematics, Jawaharlal Nehru College, Wadi, Nagpur (India) *E-mail: katrert9@gmail.com Article History: Received: 1-11-2024 Revised: 17-12-2024 Accepted: 25-12-2024 Abstract: In the present study, we presented the field equations for anisotropic charged spheres in higher dimensional spherically symmetric space-time in the framework of Rosen’s bimetric relativity theory (Rosen 1979, 1980). To obtain the determinate solutions of the higher dimensional field equations, the strategy we adopted is to assume that the space-time admits one parameter group of conformal motions generated by the vector field 𝝃𝝁 if Lπœ‰π‘”πœ‡πœˆ = πœ“β€„π‘”πœ‡πœˆ . Keywords: Bimetric theory, spherically symmetric model 1. Introduction Rosen [1, 2, 3] proposed a new theory of gravitation called the bimetric general relativity theory. In this theory, in addition to the metric tensor 𝑔𝑖𝑗 (the gravitational or physical metric), a second metric tensor 𝛾𝑖𝑗 (the background metric) has been introduced describing the geometry of space-time, which one would have if there were no matter in the universe. In this respect the theory is similar to an earlier version of bimetric theory formulated by Rosen [4, 5, 6]. The field equations of the present theory of gravitation are taken to be the same as in general relativity, except for the fact that ordinary derivatives of the physical metric are replaced by covariant derivatives with respect to the background metric. The background metric 𝛾𝑖𝑗 is taken to describe a space-time of constant curvature. It was found by Rosen [2] that these field equations could be written in the form of Einstein field equations, but with an additional term on the right-hand side, 𝐺𝑖𝑗 = 𝑆𝑖𝑗 βˆ’  8πœ‹β€„π‘‡π‘–π‘—, (1) where the Einstein tensor 𝐺𝑖𝑗 = 𝑅𝑖𝑗 βˆ’ 1 2 𝑔𝑖𝑗𝑅 and 𝑇𝑖𝑗 is the energy-momentum tensor of the matter distribution and 𝑆𝑖𝑗 = 3 π‘Ž2 (𝛾𝑖𝑗 βˆ’ 1 2 π‘”π‘–π‘—π‘”π‘™π‘˜π›Ύπ‘™π‘˜) , (2) where, π‘Ž is a constant scale parameter, of the order of the size of the universe. mailto:katrert9@gmail.com Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 905 https://internationalpubls.com From this one sees that, in general, if one is dealing with a physical system that is small compared to the size of the universe, e.g., the solar system, then the term 𝑆𝑖𝑗 is negligible in the field equations (1) and the present theory gives agreement with Einstein general relativity theory. In recent years, the study of higher dimensional space-time has gained momentum because of the fact that the cosmos at eats early stage of evolution might have had a higher dimensional era. The extra space reduced to volume with passage of time, which is beyond the ability of experimental observations at the moment. Also, the possibility that the world may have more than 4-dimensions is due to Kaluza [7] and Klein [8], who used one extra dimension to unify gravity and electromagnetism in a theory that was essentially 5-dimensional in general relativity. Sabbatta [9], Lee [10], Appelquist and Chodos [11], and Collins et al. [12] accepted this idea and constructed cosmological models in higher dimensions by using various phenomena of particle physics and cosmology. Overduin and Wesson [13] presented an excellent review of Kaluza Klein theory and higher-dimensional unified theories, in which the cosmological and astrophysical implications of extra dimensions were studied. Mohanty and Mahanta [14, 15], Mohanty and Samanta [16] are some of the authors who have studied higher-dimensional cosmology in various theories of gravitation. There is now extensive literature dealing with different aspects of higher-dimensional cosmology. Particularly, Karade and khadekar [17] have shown non-existence of higher dimensional axially symmetric field in Rosen's theory of gravitation. Reddy and Venkateswara [18] have shown non- existence of higher dimensional anisotropic cosmological model in bimetric theory of gravitation. Katore et al. [19] have shown the non-existence of n-dimensional static plane symmetric solutions in bimetric relativity theory. Khadekar and Tade [20] studied string cosmological models in five dimensional bimetric theory of gravitation. Herrera et al. [21] studied the consequences of inclusion of one parameter group of conformal motion of anisotropic matter in the Einstein’s general relativity and obtained analytical solutions of field equations for static and spherically symmetric distribution of isotropic and anisotropic matter. The study of static anisotropic spheres is important in relativistic astrophysics [22]. Sah and Prakash Chandra [23] have studied class of charged fluid balls in general relativity. Hasmani and Pandya [24] have investigated cosmological models for the static spherically symmetric spacetime with charged anisotropic fluid distribution in (n+2)-dimensions in context of Rosen's Bimetric General Relativity. Mukesh Kumar Madhukar [25] presented the study of some charged fluid spheres in General Relativity solutions and conditions for charged fluid spheres in General Relativity. Francesco Bajardi and et. al. [26] studied the higher dimensional static and spherically symmetric solutions in Extended Gauss–Bonnet Gravity. Recently, Zahraa and Mardanb [27] have presented five-dimensional analysis of electromagnetism with heat flow in the post-quasi-static approximation. With the above motivation, in this paper, we have presented the field equations in higher dimensional static spherically symmetric space-time in the framework of bimetric theory of relativity proposed by Rosen [2, 3] for anisotropic distribution of matter and also obtain the solutions for the physical metric with the assumption that the space-time admits a one-parameter group of conformal motions. This work is a generalization of the results of Hererra et al. [21] obtained in general relativity to the Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 906 https://internationalpubls.com higher dimensional bimetric theory of relativity proposed by Rosen. The solution gives agreement with Einstein’s general relativity for physical system known in the universe, such as the solar system. The work done, in this paper, is organized as follows: In section-2, we obtain the field equations in the presence of an anisotropic charged spherically symmetric matter in the frame work of higher dimensional bimetric theory of relativity proposed by Rosen [2, 3]. Section-3 deals with the solutions of the field equations under the assumption that the space time admits a one-parameter group of conformal motions [Herrera et al. (1985)] generated by the vector field πœ‰πœ‡ if Lπœ‰π‘”πœ‡πœˆ = πœ“β€„π‘”πœ‡πœˆ , (1) where L.H.S. of (1) is the Lie derivative operator and πœ“ is conformal factor. Section-4 contains some conclusions. 2. Metric and Field Equations We consider the higher dimensional static distribution of matter represented by charged spherically symmetric fluid which may be anisotropic. In terms of Schwarzschild coordinates the physical metric takes the form 𝑑𝑠2 = π‘’πœˆπ‘‘π‘‘2 βˆ’ π‘’πœ†π‘‘π‘Ÿ2 βˆ’ π‘Ÿ2𝑑𝛺2, (2) where 𝑑𝛺2 = π‘‘πœƒ1 2 + 𝑠𝑖𝑛2 πœƒ1 π‘‘πœƒ2 2 + 𝑠𝑖𝑛2 πœƒ1  𝑠𝑖𝑛2 πœƒ3 π‘‘πœƒ3 2 + β‹― + [ ∏ 𝑖=1 π‘›βˆ’1 𝑠𝑖𝑛2 πœƒπ‘–] β€„π‘‘πœƒπ‘› 2 and πœ†, 𝜈 are the functions of π‘Ÿ only. With this choice of coordinates the total energy-momentum tensor 𝑇𝑖𝑗 is assumed to be the sum of two parts, 𝑇𝑖𝑗 𝐴 and 𝑇𝑖𝑗 𝐸 , i.e. 𝑇𝑖𝑗 = 𝑇𝑖𝑗 𝐴 + 𝑇𝑖𝑗 𝐸 . (3) Here π‘»π’Šπ’‹ 𝑨 is an anisotropic spherically symmetric matter distribution and is given by 𝑇𝑖𝑗 𝐴 = (πœŒπ‘š + 𝑃βŠ₯)β€„π‘ˆπ‘–β€„π‘ˆπ‘— βˆ’ 𝑃βŠ₯𝑔𝑖𝑗 + (π‘ƒπ‘Ÿ βˆ’ 𝑃βŠ₯)β€„πœ’π‘–β€„πœ’π‘—, (4) where π‘ˆπ‘– is the time like (𝑛 + 2) -velocity π‘ˆπ‘– = 𝛿(𝑛+2) 𝑖 π‘’βˆ’πœˆ/2, πœ’π‘– a unit space-like vector in the radial direction πœ’π‘– = 𝛿1 π‘–π‘’βˆ’πœ†/2, πœŒπ‘š the energy density, π‘ƒπ‘Ÿ the radial pressure in the direction of πœ’π‘– and 𝑃βŠ₯ the pressure orthogonal to πœ’π‘–. Similarly, the electromagnetic contribution can be written as π‘»π’Šπ’‹ 𝑬 = βˆ’ 𝟏 πŸ’π… (π‘­π’Š πœΆπ‘­π’‹πœΆ βˆ’ 𝟏 πŸ’ π’ˆπ’Šπ’‹β€„π‘­π’π’Œπ‘­π’π’Œ), (5) where π‘­π’Šπ’‹ is the electromagnetic field tensor defined in terms of the (𝒏 + 𝟐)-potential π‘¨π’Š as π‘­π’Šπ’‹ = π‘¨π’Š, 𝒋 βˆ’ 𝑨𝒋,π’Š . (6) Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 907 https://internationalpubls.com As the frame is considered to be the rest frame, we adopt the gauge π‘¨π’Š = π‘¨π’Š(𝟎, β€„πŸŽ, β€„πŸŽ, ⋯ (𝒏 + 𝟏)β€„π’•π’Šπ’Žπ’†π’”β€„, 𝝋(𝒓)). Also, π‘­π’Šπ’‹ satisfy the Einstein-Maxwell field equations 𝐹𝑖𝑗;π‘˜ + πΉπ‘—π‘˜;𝑖 + πΉπ‘˜π‘–;𝑗 = 0, (7) 𝐹;𝑗 𝑖𝑗 = βˆ’4πœ‹β€„π½π‘–, (8) where𝐽𝑖 being the charge-current density that becomes ie i uJ βˆ’ = ( e βˆ’  is the proper charge density). In the commoving system, we choose π‘ˆπ‘– = ( 0,  0,  0, ⋯ , (𝑛 + 1)0, π‘ˆπ‘›+2) πœ’π‘– = (πœ’1,  0,  0, ⋯ ,0(𝑛 + 1)) From π‘ˆπ‘–β€„π‘ˆπ‘– = βˆ’πœ’π‘–β€„πœ’π‘– = 1, we obtain π‘ˆ(𝑛+2) = π‘’βˆ’πœˆ/2,β€„β€„πœ’1 = π‘’βˆ’πœ†/2. (10) The non-vanishing components of energy momentum tensors 𝑇𝑗 𝑖 are =βˆ’====βˆ’= + +βŠ₯ + + 2n 2n AA 1n 1n A 3 3 A 2 2r A 1 1 T,PTTT,PT  and 21n 1n E 4 4 E 3 3 E 2 2 E 22n 2n E 1 1 E ETTTT,ETT βˆ’======= + + + +  . Now let us tern to the gravitational field equations. In bimetric theory of reletivity they have the form [Rosen, (1980)] οƒ· οƒΈ οƒΆ    +ο°βˆ’=βˆ’βˆ’ο‚Ί E ij A ijijijijij TT8SRg 2 1 RG , (12) where the additional term 𝑆𝑖𝑗 is given by 𝑆𝑖𝑗 = 3 π‘Ž2 (𝛾𝑖𝑗 βˆ’ 1 2 𝑔𝑖𝑗𝑔𝛼𝛽𝛾𝛼𝛽). (13) If we take the background metric 𝛾𝑖𝑗 in a higher dimensional static de-Sitter form, the line element in terms of polar coordinates (π‘₯1, π‘₯2 β‹― , π‘₯𝑛+1,  π‘₯𝑛+2) = (π‘Ÿ, πœƒ1, β€„πœƒ2,  ⋯ , πœƒπ‘›,  𝑑) is given by π‘‘πœŽ2 = (1 βˆ’ π‘Ÿπ‘› π‘Ž2)  𝑑𝑑2 βˆ’ (1 βˆ’ π‘Ÿπ‘› π‘Ž2) βˆ’1 π‘‘π‘Ÿ2 βˆ’β€„ π‘Ÿ2𝑑𝛺2. (14) For a region very small as compared to π‘Ž (i.e. for π‘Ÿ << π‘Ž), this line element has the flat space form π‘‘πœŽ2 = 𝑑𝑑2 βˆ’ π‘‘π‘Ÿ2 βˆ’ π‘Ÿ2𝑑𝛺2. (15) (9) (11) Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 908 https://internationalpubls.com Let us consider a region with π‘Ÿ << π‘Ž and neglect the terms, which are small throughout this region, we can write the non-vanishing components of 𝑆𝑖 𝑗 in higher dimensions (Falik and Rosen, 1980) βˆ’π‘†1 1 = βˆ’π‘†2 2   =   βˆ’π‘†3 3 = β‹― = βˆ’π‘†π‘›+1 𝑛+1 = 𝑆𝑛+2 𝑛+2 = 3 2π‘Ž π‘’βˆ’πœˆ. (16) Rosen’s (1980) bimetric field equations given by equation (12) for the for the line-element (2) (for π‘Ÿ << π‘Ž) can be written as βˆ’π‘’βˆ’πœ† [ 𝑛(π‘›βˆ’1) 2π‘Ÿ2 βˆ’ π‘›πœ†β€² 2π‘Ÿ ] + 𝑛(π‘›βˆ’1) 2π‘Ÿ2 = βˆ’ 3 2π‘Ž2 π‘’βˆ’πœˆ + .8πœ‹πœŒ + 𝐸2 (17) βˆ’π‘’βˆ’πœ† [ 𝑛(π‘›βˆ’1) 2π‘Ÿ2 + π‘›πœˆβ€² 2π‘Ÿ ] + 𝑛(π‘›βˆ’1) 2π‘Ÿ2 = 3 2π‘Ž2 π‘’βˆ’πœˆ βˆ’ 8πœ‹π‘ƒπ‘Ÿ + 𝐸2 (18) π‘’βˆ’πœ† [ πœˆβ€³ 2 + πœˆβ€²2 4 βˆ’ (𝑛 βˆ’ 1)(πœ†β€² βˆ’ πœˆβ€²) 2π‘Ÿ βˆ’ πœ†β€²πœˆβ€² 4 + (𝑛 βˆ’ 1)(𝑛 βˆ’ 2) 2π‘Ÿ2 ] βˆ’ (𝑛 βˆ’ 1)(𝑛 βˆ’ 2) 2π‘Ÿ2         = βˆ’ 3 2π‘Ž2 π‘’βˆ’πœˆ + 8πœ‹π‘ƒβŠ₯ + 𝐸2, (19) where a prime denotes a derivative with respect to π‘Ÿ. We now define the effective density and pressures πœŒπ‘šπ‘’  ,β€„β€„π‘ƒπ‘Ÿπ‘’ and 𝑃βŠ₯𝑒 (Harpaz and Rosen, 1985) as πœŒπ‘šπ‘’ = πœŒπ‘š βˆ’ 3 16πœ‹π‘Ž2 π‘’βˆ’πœˆ , β€„β€„β€„π‘ƒπ‘Ÿπ‘’ = π‘ƒπ‘Ÿ βˆ’ 3 16πœ‹π‘Ž2 π‘’βˆ’πœˆ and 𝑃βŠ₯𝑒 = 𝑃βŠ₯ βˆ’ 3 16πœ‹π‘Ž2 π‘’βˆ’πœˆ, so that the field equations (17) – (19) take the form βˆ’π‘’βˆ’πœ† [ 𝑛(π‘›βˆ’1) 2π‘Ÿ2 βˆ’ π‘›πœ†β€² 2π‘Ÿ ] + 𝑛(π‘›βˆ’1) 2π‘Ÿ2 = 8πœ‹πœŒπ‘šπ‘’ + 𝐸2 (21) βˆ’π‘’βˆ’πœ† [ 𝑛(π‘›βˆ’1) 2π‘Ÿ2 + π‘›πœˆβ€² 2π‘Ÿ ] + 𝑛(π‘›βˆ’1) 2π‘Ÿ2 = βˆ’8πœ‹π‘ƒπ‘Ÿπ‘’ + 𝐸2 (22) π‘’βˆ’πœ† [ πœˆβ€³ 2 + πœˆβ€²2 4 βˆ’ (𝑛 βˆ’ 1)(πœ†β€² βˆ’ πœˆβ€²) 2π‘Ÿ βˆ’ πœ†β€²πœˆβ€² 4 + (𝑛 βˆ’ 1)(𝑛 βˆ’ 2) 2π‘Ÿ2 ] βˆ’ (𝑛 βˆ’ 1)(𝑛 βˆ’ 2) 2π‘Ÿ2         = 8πœ‹π‘ƒβŠ₯𝑒 + 𝐸2. (23) These equations look like the usual Einstein equations with πœŒπ‘š, π‘π‘Ÿ and 𝑝βŠ₯ replaced by πœŒπ‘šπ‘’ , β€„β€„π‘π‘Ÿπ‘’ and 𝑝βŠ₯𝑒 respectively. It is very difficult to obtain the solution of equations (21)-(23) due to the non-linearity of the field equations and therefore one has to make certain simplifying assumption to derive the useful results. 3. Conformal motions and solutions of field equations In this section, we assume that the physical metric (2) admits a one-parameter group of conformal motions and obtained solution of the field equations for static and spherically symmetric fluid distributions of anisotropic matter. The field equations (21)-(23) resulting in the new form (20) Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 909 https://internationalpubls.com 𝑛(π‘›βˆ’1) 2π‘Ÿ2 [1 βˆ’ πœ“2 𝐢3 2] βˆ’ π‘›πœ“πœ“β€² 𝐢3 2π‘Ÿ = 8πœ‹πœŒπ‘’ + 𝐸2, (24) βˆ’ π‘›πœ“2 𝐢3 2π‘Ÿ2 [1 βˆ’ (π‘›βˆ’1) 2 ] + 𝑛(π‘›βˆ’1) 2π‘Ÿ2 = 8πœ‹π‘ƒπ‘Ÿπ‘’ + 𝐸2, (25) π‘›πœ“πœ“β€² π‘ŸπΆ3 2 + 𝑛(π‘›βˆ’1)πœ“2 2π‘Ÿ2𝐢3 2 βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 2π‘Ÿ2 = 8πœ‹π‘ƒβŠ₯𝑒 + 𝐸2, (26) We now define a function π›₯, as a measure of anisotropy, such that π›₯ = 4πœ‹(𝑃βŠ₯𝑒 βˆ’ π‘ƒπ‘Ÿπ‘’ ). (27) Also, if we take 𝑋 = πœ“2 𝐢3 2, then in terms of the function defined in (27), the equations (24) –(26) reduces to : 𝑛(π‘›βˆ’1) 4π‘Ÿ2 βˆ’ 𝑛𝑋 2π‘Ÿ2 + 𝑛𝑋′ 4π‘Ÿ βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 βˆ’ π›₯ = 𝐸2 (28) 𝑛(π‘›βˆ’1) 4π‘Ÿ2 βˆ’ 3𝑛𝑋′ 4π‘Ÿ βˆ’ 𝑛(π‘›βˆ’2)𝑋 4π‘Ÿ2 βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 + π›₯ = 8πœ‹πœŒπ‘šπ‘’ (29) 𝑛𝑋 2π‘Ÿ2 βˆ’ 𝑛(π‘›βˆ’1) 4π‘Ÿ2 + 𝑛𝑋′ 4π‘Ÿ βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 βˆ’ π›₯ = 8πœ‹π‘ƒπ‘Ÿπ‘’ . (30) From which we obtain 𝑛(π‘›βˆ’1) 4π‘Ÿ2 βˆ’ 𝑛𝑋 2π‘Ÿ2 + 𝑛𝑋′ 4π‘Ÿ βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 βˆ’ π›₯ = 𝐸2 (31) 𝑛(π‘›βˆ’1) 4π‘Ÿ2 βˆ’ 3𝑛𝑋′ 4π‘Ÿ βˆ’ 𝑛(π‘›βˆ’2)𝑋 4π‘Ÿ2 βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 + π›₯ = βˆ’3 2π‘Ž2 π‘’πœˆ + 8πœ‹πœŒπ‘šπ‘’ (32) 𝑛𝑋 2π‘Ÿ2 βˆ’ 𝑛(π‘›βˆ’1) 4π‘Ÿ2 + 𝑛𝑋′ 4π‘Ÿ βˆ’ (π‘›βˆ’1)(π‘›βˆ’2) 4π‘Ÿ2 βˆ’ π›₯ = 3 2π‘Ž2 π‘’πœˆ βˆ’ 8πœ‹π‘ƒπ‘Ÿπ‘’ . (33) Therefore, the line element (2), in terms of π‘’πœˆ = 𝐢2 2π‘Ÿ2, πœ“ = 𝐢3 2π‘’βˆ’πœ†/2 and 𝑋 = πœ“2 𝐢3 2 now becomes 𝑑𝑠2 = 𝐢2 2π‘Ÿ2𝑑𝑑2 βˆ’ ( 1 𝑋 ) π‘‘π‘Ÿ2 βˆ’ π‘Ÿ2𝑑𝛺2. (34) 4. Conclusions In this paper, we have presented the field equations for anisotropic charged spheres in higher dimensional Rosen’s bimetric theory of relativity. To get a determinate solution of the field equations, we have assumed that the space-time admits one parameter group of conformal motion. It is observed that the solutions obtained in bimetric relativity in higher dimensions will reduce to the corresponding solutions of Einstein’s field equations obtained by Herrera et al. [21] of a physical system, which is very small compared to the size of the universe, when the term 3 2π‘Ž2 π‘’βˆ’πœˆ that appears in the field equations of bimetric relativity is negligible. Communications on Applied Nonlinear Analysis ISSN: 1074-133X Vol 31 No. 8s (2024) 910 https://internationalpubls.com References [1] Rosen, N.: Lett. Nuovo Cim. 25, 266 (1979) [2] Rosen, N.: Found. Phys. 10, 673 (1980a) [3] Rosen, N.: Gen. Rel. Grav. 12, 493 (1980b) [4] Rosen, N.: Ann. Phys. 84, 455 (1974) [5] Rosen, N.: Astrophys. J. 221, 284 (1978b) [6] Rosen, N.: Gen. Rel. Grav. 9, 339 (1978a) [7] Kaluza, T. Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys.), 1, 966-972 (1921) [8] Klein, O. Zeits. Phys. 37, 895-906 (1926) [9] De Sabbatta, V.; Schmutzer, E. Unified Field Theories of More than Four Dimensions Including Exact Solutions; World Scientific: Singapore (1983) [10] Lee, H. C. An Introduction to Kaluza-Klein Theories; World Scientific: Singapore (1984) [11] Appelquist, T.; Chodos, A.; Freund, P. G. O. Modern Kaluza-Klein Theories; Addison-Wesley: Boston, MA, USA (1987) [12] Collins, P. D. B.; Martin, A. D.; Squires, E. J. Particle Physics and Cosmology; Wiley: New York, NY, USA (1989) [13] Overduin, J. M.; Wesson, P. S. Phys. Rep. 283, 303-378 (1997) [14] Mohanty, G.; Mahanta, K. L. Turk. J. Phys. 31, 299-306 (2007) [15] Mohanty, G.; Mahanta, K. L. Turk. J. Phys., 32, 299-303 (2008) [16] Mohanty, G.; Samanta, G. C. Turk. J. Phys. 32, 251-260 (2008) [17] Karade, T. M.; G. S. Khadekar: Czechoslovak J. of Physics, 39, 9, 962-967 (1989) [18] Reddy D.R.K., Rao Venkateswara: Astrophys. Space Sci., 271, 2, 165-169 (2000) [19] Katore S. D., Adhav K.S., Thakare S.V.: Astrophys. Space Sci., 299, 3, 233 (2005) [20] Khadekar G. S., Tade S. D.: Astrophys. Space Sci., 310, 1, 47-51 (2007) [21] Herrera, L. and Ponce de Leon, J. (1985) J. Math. Phys. 26, 2018. [22] Bowars, R. and Liang, E. : Astrophys. J. 188, 657 (1974) [23] Sah, A.; Prakash Chandra: International Journal of astronomy and Astrophysics, 17, 4, 494-511 (2016) [24] Hasmani, A.H., Pandya, D.N.: Journal of Dynamical System and Geometric Theories, 17(1), 13- 21 (2019) [25] Mukesh Kumar Madhukar: Journal of Advances and Scholarly Researches in Allied Education, 17, 1, 606 - 608 (3) (2020) [26] Francesco Bajardi, Konstantinos F., Dialektopoulos, Capozziello, S.: Symmetry, 12 (3), 372 (2020) [27] Zahraa, A.; Mardanb, S.A.: The European Physical Journal C, 83:231 (2023)