High-derivatives and massive electromagnetic models in the Lemaitre-Tolman-Bondi spacetime

The Maxwell electromagnetic theory embedded in an inhomogeneous Lema\^{\i}tre-Tolman-Bondi (LTB) spacetime background was described a few years back in the literature. However, terms concerning the mass or high-derivatives were no explored. In this work we studied the inhomogeneous spacetime effects on high-derivatives and massive electromagnetic models. We used the LTB metric and calculated the physical quantities of interest, namely the scale factor, density of the electromagnetic field and Hubble constant, for the Proca and higher-derivative Podolsky models. We found a new singularity in both models, and that the magnetic field must be zero in the Proca model.


I. INTRODUCTION
Homogeneity and isotropy, together with matter being treated as a specific gas, are the basic ingredients of the Friedmann-Lemaître-Robertson-Walker (FLRW) model. Since the pioneering work of E. Hubble, who suggested a correlation between the observed redshifts and distances of 24 galaxies, this is a very successfull model, being currently considered a very good measure to describe our cosmos.
In the last decades a great effort was dedicated to understanding the local inhomogeneities that occur in the Universe. Such efforts suggest that these inhomogeneities may be related to the expansion of the Universe. Inhomogeneities, as an alternative to dark energy, were first discussed in Ref. [1] as a means of explaning the observational results of the expansion of the Universe [2][3][4][5][6][7][8] without the need for postulating dark energy. Inhomogeneities inside astronomical objects can define different instability ranges, an effect that can describe distinct features of their evolution and structures formation.
The current standard model of cosmology, the Λ-Cold Dark Matter (ΛCDM) model, is a homogeneous solutions of the FLRW Einstein's field equations with only six free parameters that successfully accounts for most cosmological data, especially the characteristics of the Cosmic Microwave Background (CMB) and the structure formation on large scales considered through the theory of cosmological perturbations in homogeneous and isotropic background. However, in the last fifteen years "standard" inhomogeneous cosmological models that are generalizations of FLRW cosmologies have been the subject of growing interest in astrophysics community in order to investigate cosmological phenomena. Some authors have demonstrated that inhomogeneous models with spherical symmetry and dust source can be fitted to supernovae Ia (SNIa) data, as well as the position of the first peak of the CMB. These models show that the apparent accelerated expansion of the universe may not be a consequence of the repulsive gravity due to dark energy, but rather the result of inhomogeneities in the distribution of matter. In this context one needs to mention that an important cosmological model to describe the inhomogeneity universe is the Lemaître-Tolman-Bondi (LTB) spacetime [9][10][11][12][13][14][15][16][17], which is a spatially inhomogeneous description of a spherically symmetric distribution of dust matter in the Universe.
In the last few years several studies analyzing the scenario of electrodynamics embedded in the isotropic and homogeneous FLRW gravitational background models [18] have been produced. Astrophysical effects were explored together with an anisotropic expansion of the Universe. The interesting result of polarized electromagnetic radiation occurs when it travels through local anisotropic regions. Concerning the inhomogeneities, in Ref. [19] the authors investigated that the inhomogeneity with electromagnetic field caused a new scale factor. The propagation of photons was also affected, which is important phenomenum since most information obtained about the Universe is by means of photons. Ref. [20] studied the effects of Palatini f (R) gravity together with the so-called tilted observer on the dynamics of LTB spacetime embedded in an electromagnetic field.
In this work we studied the effect of a high-derivative electromagnetic field embedded in the inhomogeneous LTB geometry. We analyzed the Proca and high-derivatives Podolsky electromagnetic models. The electrodynamic contribution was inserted separately through the energy-momentum tensor (EMT) of the respective models. For the Proca model, besides the matter density there also are the electromagnetic contributions obtained from the Lagrangian of the Proca model in curved space-time. An analysis of the Proca model in curved space-time was carried out for a particular case of interest by Bekenstein [21]. Here, we have analyzed the electrodynamics effects in LTB cosmological model and calculate the scale factor in LTB universe. We also computed the luminosity distance in the presence of electromagnetic field.
Concerning the Podolsky electrodynamics effects in LTB cosmological model, the study of Podolsky model in curved space-time was made in Ref. [22] together with the analysis of the Bopp-Podolsky black holes. Therefore, we have obtained the line element and the equations that define the LTB model with Podolsky contributions. This paper is structured as follows. In Section 2 we review the main aspects of the LTB metric, where the Einstein tensor and the matter contribution for the EMT are obtained. In Section 3 we analyze the Proca model of electrodynamics in curved space-time, where the EMT and the Maxwell-Proca equations in curved space-time are derived. In Section 4 we solve the Einstein equations including the Proca contribution. We obtain the scale factor for this model and examine the inhomogeneities and luminosity distance. In section 5 we analyze the Podolsky electrodynamics contributions for LTB model. Section 6 discusses the results and presents some final considerations.

II. THE LTB COSMOLOGICAL MODEL: A BRIEF REVIEW
The LTB model depicts a self-gravitating spherically symmetric distribution of inhomogeneous nondissipative dust cloud where the EMT can be written as T µν = ǫ(τ, ρ)u µ u ν and u µ = u µ (τ, ρ) is the dust particle's four-velocity vector. The proper time is represented by τ and the several shells are labeled by ρ, which helps in the formation of the dust cloud. As mentioned above, the LTB metric is the most popular inhomogeneous cosmological model, since it is appropriate for both small and large scale inhomogeneities, and also is the simplest inhomogeneous solution of the Einstein equations.
In this section we will analyze the LTB model through the point of view of the field equations of general relativity, as well as the matter contribution concerning the EMT. In LTB cosmology, the line element is written as where A and R, unlike FLRW model, depend also on the r coordinate, namely, A≡A(r, t) and R≡R(r, t). Hence, the metric element is given by where The spatially homogeneous FLRW metric is a special case of the metric in Eq. (1), where and R(r, t) = a(t)r, where a(t) is the scale factor. The Einstein equations, where c = 1, are given by where G µν is the Einstein tensor G is the gravitational constant and T µν is the EMT. So, we can make use of the line element in Eq. (1) to obtain the non-zero components of the Einstein tensor, In LTB cosmology, the EMT is the zero pressure diagonal perfect fluid, (ρ, 0, 0, 0), where ρ is the mass density. So, the matter contribution for the EMT T M 00 is The electromagnetic contribution for the LTB cosmology will be obtained by adding on the right side of the Einstein equations. So, we have for the EMT, for both the Proca model and for the Podolsky model, that where T EM µν is the electromagnetic contribution for the EMT. Therefore, determining T EM µν , then we have the Einstein equation.

III. THE CURVED SPACETIME PROCA MODEL
The Proca model of electrodynamics is defined by the following Lagrangian where m is the mass of electromagnetic field, A µ is the four-potential and F µν is the electromagnetic tensor that, in curved space-time, is defined by where ∇ is the covariant derivative given by Moreover, the electromagnetic field and its dual can be obtained from the following expressions, and the Maxwell-Proca equations in curved spacetime are where J µ = (ρ, J) is the electromagnetic four-current density and * F µν is the dual electromagnetic tensor. From Eq. (10), the action for the Proca model is namely, and, from the EMT definition the Proca EMT in curved spacetime is In this work, we neglect the coupling term J αÃ α because most of the matter is electrically neutral. Thus, any fluctuations that may appear can be ignored. The trace of the energymomentum tensor is: i.e., unlike Maxwell's case, the tensor trace is not zero.
Using relations in Eq. (17), we find that the non-zero components of the EMT are The next step is introduce Eqs. (19) into the Einstein equations in Eq. (5). We will do this in the next section and we will analyze the effects of the Proca model on the LTB model.

IV. LTB SCENARIO PROCA ELECTRODYNAMICS
Now that we have the Proca's contribution in hand, we must analyze the terms of the Einstein tensor to obtain Einstein's equations. As we can see from Eq. (7), the only nonzero off-diagonal term in G µν is G 10 . The only way to satisfy this condition is to impose that and, consequently Taking into account the last conditions, the electromagnetic fields are as we can see, as consequence of conditions in Eq. (21), we have that the B µ field is zero. The only contribution of electromagnetic field is by E µ field. Now, imposing these same conditions in equations (19), the energy-momentum components are In order to write A in function of R, we impose that and the Einstein equations stay with matter e electromagnetic contribution in EMT are The last equation of Eqs. (25) allows us as to write A in terms of R ′ such as where k(r) is an arbitrary function. So, the line element is Moreover, we can rewrite the Einstein equations for the (00) and (11) components likė On the other hand, the non-zero components of the Maxwell-Proca equations in Eqs. (12) are So, we can write the electric field E(t, r) using the integration of Eq. (30), as where and ǫ(t) is the constant of integration. Now, from Eq. (31), using the fact that the charged matter belongs to the comoving matter, namely, the four-current appears as J α = (J 0 , 0) which leads us to Now, defining multiplying the equation (29) byṘ, we obtain thaṫ or yet, using the fact that ∂ 0 (RṘ) = 2AṘR, we have and using Eq. (28) we can write Eqs. (27), (40) and (41) define the LTB model with the Proca electrodynamics contributions. The singularities arising from R = 0, R ′ = 0 and k = 1. The R = 0 singularity we interpreted as the Big Bang singularity, and the R ′ = 0 singularity as the shell cross singularity. The last singularity, k = 1, that occurs if R ′ = 0, come from the Proca contribution. Now, defining where the Hubble constant is defined as We have imposed the boundary values at t 0 through A 0 (r)≡A(t 0 , r), H 0 (r)≡H(t 0 , r). Moreover, Ω k , Ω m and Ω σ are subjected to the constraint Ω k (r) + Ω σ (r) + Ω m (r) = 1 .
Therefore, we can write Eq. (40) aṡ and the Hubble constant is computed by and, considering the solution using thatṘ > 0, we have where M(t, r) = Ω 2 m (r) where R 0 (r) corresponds to the current shape of the scale factor, H 0 (r) is the current value of the Hubble constant in each point and Ω σ (r) is the density of the electromagnetic field, that, because of Proca's contributions, is purely electric.
Remember that the density of the electromagnetic field is given by and the Proca contribution for the density of the electromagnetic field is the mass term −2e −A m 2 A 2 0 . Now, let us examine the inhomogeneities and luminosity distance for the Proca electrodynamics perspective. Therefore, the geodesic requires that dθ = dφ = 0. Moreover, since light always travels along null geodesics, we have dS 2 = 0. So, using Eq. (27), we have Consider two light rays with solutions of Eq. (53) given by t 1 = t(u) and t 2 = t(u)+λ(u). Substituting these two into Eq. (53) we obtain Differentiating the definition of the redshift, z ≡ λ(0) − λ(u) /λ(u), we have, and, consequently, which determine the relation between the coordinates and the observable redshift. For the last equation, using the expression for k(r) we find dr dz = 1 + H 2 0 (r)(1 − Ω m (r) + Ωσ(r))R 2 0 (r) (1 + z)Ṙ ′ (t, r) .
Moreover, the relation between the redshift and the energy flux F, defined as d L ≡ L/(aπF ), where L is the total power distance radiated by the source, is and the angular distance diameter is given by which is a direct relation to the scale factor as a function of the redshift.

V. LTB GEOMETRY PODOLSKY ELECTRODYNAMICS
The Podolsky electrodynamics in curved space-time is given by or, where R σβ is the Ricci tensor. The Einstein-Podolsky action is which leads us, from the variation with respect to A µ , to the Einstein-Podolsky equation and The EMT for the Podolsky electrodynamics in curved space-time is given by and the trace of the EMT is Now, using the metric element Eq. (1), the non-null elements for the energy-momentum tensor are However, by Eq. (7), the only non-null off diagonal component of Einstein tensor is G 1 0 . Therefore, in order to have T 02 = T 03 = T 12 = T 13 = T 23 = 0 (79) we need to impose the conditions After imposing these conditions into Eqs. (69) − (78) we have that So, together with Eq. (7), we have Hence, in order to obtain a more concise notation, let us define Hence, that Eq. (88) can be written as where, together with the definition of covariant derivative in Eq. (11) we have that A(t, r) − E(t, r) , Integrating Eq. (90), we obtain where where α is the constant of integration and, as we can see, if the M is zero, we obtain the Maxwell case discussed in [19]. Therefore, defining χ ≡ 1 − k(t, r), we obtain and the line element in Eq. (1) can be write as The new result here is that the k-term depends on both r and t coordinate.

VI. CONCLUSION
The analysis of the electromagnetic field in an LTB background was explored by other authors in Ref. [19], however they studied only the Maxwell case, which means that the effects of a mass term or the complication of high-derivative electromagnetic term were not investigated until now. The objective of this paper is to fill this gap and to analyze the possible theoretical effects of these rich electrodynamic scenarios in a more realistic inhomogeneous cosmological background.
In this work, we have analyzed the Proca and the higher-derivative Podolsky models embedded in an inhomogeneous LTB background. In the case of Proca model we found a new singularity at k(r) = 0. Moreover, the magnetic field must be zero to satisfy the Einstein's equations. Considering the Proca model, we have determined the scale factor and provided an analysis of the luminosity distance.
In the Podolsky case, we analyzed the function k, which appears in the line element that defines the LTB model where the Podolsky electrodynamics is dependent of both t and r coordinate. This is different from Maxwell and Proca cases.
As a perspective, a more precise analysis of the Podolsky contributions for the cosmological LTB model can be carried out, as the analysis of black holes to generate the singularity at R ′ = 0, for example.