New algorithms to obtain analytical solutions of Einstein's equations in isotropic coordinates

TL;DR

This paper introduces two novel methods inspired by the MGD approach to derive new anisotropic solutions of Einstein's equations in isotropic coordinates, demonstrating their effectiveness through four analytical solutions.

Contribution

It presents two inequivalent algorithms for solving Einstein's equations with anisotropy in isotropic coordinates, expanding the toolkit for finding exact solutions.

Findings

Four analytical solutions obtained, two are physically acceptable.

Methods are inspired by the MGD approach in Schwarzschild coordinates.

Demonstrates the applicability of new algorithms to generate solutions.

Abstract

The main objective of this work, is to show two inequivalent methods to obtain new spherical symmetric solutions of Einstein's Equations with anisotropy in the pressures in isotropic coordinates. This was done inspired by the MGD method, which is known to be valid for line elements in Schwarzschild coordinates. As example, we obtained four analytical solutions using Gold III as seed solution. Two solutions, out of four, (one for each algorithm), satisfy the physical acceptability conditions.

Figures

Equations

$$d{s}^2=e^{\nu (r)}d{t}^2-\frac{1}{\mu (r)}d{r}^2-r^2\left(d{\theta }^2+{\sin }^2\theta d{\varphi }^2\right).$$

$$R_{\mu \nu }-\frac{1}{2}R{g}_{\mu \nu }=-8\pi T_{\mu \nu }^{.},$$

$$8\pi T_0^0$$

$$-8\pi T_1^1$$

$$-8\pi T_2^2$$

$${\nabla }_{\mu }{T}^{\mu \nu }=0,$$

$$8\pi (T_1^1{)}^{\prime }+\frac{16\pi }{r}(T_1^1-T_2^2)+4\pi {\nu }^{\prime }(T_1^1-T_0^0)=0$$

$$T_{\mu \nu }=T_{\mu \nu }^{(\mathrm{PF})}+\alpha {\theta }_{\mu \nu },$$

$$T_{\mu \nu }^{(\mathrm{PF})}=(\rho +p)u_{\mu }{u}_{\nu }-p{g}_{\mu \nu },$$

$$\mu (r)\equiv 1-\frac{8\pi }{r}{\int }_0^r{x}^2\rho dx=1-\frac{2m(r)}{r},$$

$$\nu$$

$$\mu$$

$$8\pi \rho$$

$$8\pi p$$

$$p^{\prime }+\frac{{\nu }^{\prime }}{2}\left(\rho +p\right)=0,$$

$$8\pi {\theta }_0^0$$

$$8\pi {\theta }_1^1$$

$$8\pi {\theta }_2^2$$

$${\left({\theta }_1^1\right)}^{\prime }-\frac{{\nu }^{\prime }}{2}\left({\theta }_0^0-{\theta }_1^1\right)-\frac{2}{r}\left({\theta }_2^2-{\theta }_1^1\right)=0.$$

$$d{s}^2=e^{\nu (r)}d{t}^2-\frac{1}{\omega }(d{r}^2+r^{2}d{\Omega }^2),$$

$$8\pi T_0^0$$

$$-8\pi T_1^1$$

$$-8\pi T_2^2$$

$$r=A\exp \left\{\int {\mu }^{-1/2}\frac{d{r}_1}{r_1}\right\},$$

$$8\pi T_0^0$$

$$-8\pi T_1^1$$

$$-8\pi T_2^2$$

$$d{s}^2=\left(1-\frac{2M}{r_1}\right)d{t}^2-{\left(1-\frac{2M}{r_1}\right)}^{-1}d{r}_1^2+r_1^{2}d{\Omega }^2,$$

$$e^{{\nu }_{\Sigma }}$$

$$\frac{r_{\Sigma }}{2}{\left(\frac{2}{r_{\Sigma }}-\frac{w^{\prime }}{w}\right)}_{\Sigma }$$

$$P_r(r_{\Sigma })$$

$$M=\frac{r_{\Sigma }}{2w^{1/2}}{\left[1-\frac{r_{\Sigma }^2}{4}{\left(\frac{2}{r_{\Sigma }}-\frac{w^{\prime }}{w}\right)}^2\right]}_{\Sigma }$$

$$8\pi \tilde{\rho }(r)$$

$$8\pi {\tilde{p}}_r(r)$$

$$8\pi {\tilde{p}}_t(r)$$

$$\tilde{\rho }\equiv \rho +\alpha {\theta }_0^0,{\tilde{p}}_r\equiv p_r-\alpha {\theta }_0^0,{\tilde{p}}_t\equiv p_t-\alpha {\theta }_2^2.$$

$$\omega \mapsto \tilde{\omega }=\omega +\alpha f,$$

Full text

Anisotropic Fluids, Isotropic Coordinates, Interior Solutions.

I Introduction

An interesting problem for both, mathematicians and physicists, is to find new analytical solutions of Einstein’s equations. This is not an easy task because, in general, it corresponds to a system of coupled differential equations. Moreover, in order to be a realistic solution, it must satisfy certain conditions for physical acceptability. For these reasons, most of the known solutions to Einstein’s equations in the literature, are related to specific restrictions about the geometry of the space-time or particular distributions of matter (like the isotropic perfect fluid) Stephani ; lake2 ; Delgaty .

Recently, a method was discovered in the context of Randall-Sundrum Brane World Randall1 ; Randall2 named Minimal Geometric Deformation (MGD) Ovalle1 ; Ovalle2 . It has many different applicationsOvalle4 ; Ovalle5 ; Ovalle6 ; Ovalle7 ; Ovalle8 ; Ovalle9 ; Ovalle10 ; Ovalle11 ; Casadio ; Darocha1 ; Darocha2 ; Darocha3 ; Contreras ; Contreras2 ; Contreras3 ; Contreras4 ; Contreras5 ; Contreras6 and it allows us to find new analytical solutions to Einstein’s equations by decoupling gravitational sources in GR Ovalle . This solution could describe the inner and outer regions of compact objects, black holes and the coupling of the Einstein’s equation with a scalar field, among many other possibilities. It is important to mention that these solutions are in general anisotropic in the pressures. This is really important when dealing with an internal solution. There are many theoretical reasons to consider local anisotropy of the pressures Herrera2 ; Herrera3 ; Herrera4 ; Herrera5 in order to describe more realistic gravitational objects. In lhl we found 2 new families of anisotropic solutions, and we obtained an interval of values where the free parameters shall dwell in order to reproduce a physically acceptable solution. There are a good number of references that shows the efficiency of this method to obtain useful solutions to describe realistic systems that considers anisotropy in the pressures Gabbanelli ; Tello1 ; Morales1 ; Morales2 ; Estrada1 ; Estrada2 . However, MGD requires spherical symmetry, and considers solutions of Einstein’s equations whose line element, is written in Schwarzschild’s coordinates (also known by the name of canonical coordinates).

Nevertheless, there are other known exact solutions of Einstein’s equations written in other coordinates systems. In some of these cases, the transformation to move to the Schwarzschild coordinates is not always well defined. This becomes evident when we work in isotropic coordinates. It has been shown (see for example Nariai ) that every line element given in Schwarzschild-like coordinates can be transformed into the isotropic one, but the reverse process is not always possible. This seems to indicate that the isotropic coordinates are more general that canonical coordinates. The main feature of the isotropic line elements is that it treats the 3 spatial dimensions in the same way. Also, this system of coordinates has been adopted by physicist and astronomers, over the Schwarzschild coordinates in order to describe the solar system Green .

In this work we present two inequivalent algorithms to obtain new analytical and anisotropic solutions of Einstein’s equations, starting with known solutions in isotropic coordinates. These two algorithms are based in the same ansatz of the MGD method. However, there is no decoupling of Einstein’s equation in the first one, while in the second one, it is necessary to assume additional conditions that are not included in the MGD aproach. In order to check the validity of the methods, we have chosen a known seed solution in isotropic coordinates, as Gold III Gold . We have founded four different analytical solutions. A few different algorithms that also allow us to obtain new solutions from a line element in isotropic coordinates are explained in Maharaj ; Maharaj2 ; Maharaj3 . However, our method is simpler and it will lead us, generally, to solutions with anisotropy in the pressures, which is not the case in Maharaj .

This paper is organized as follows: in section 2 we briefly review Einstein’s equations in standard coordinates. In section 3 we describe the solution-generating technique known as MGD, which is useful with line elements in Schwarzschild-like coordinates. We briefly discuss Einstein’s equations on isotropic coordinates in section 4. The acceptability conditions for the solutions of Einstein’s equations are presented in the section 5. In section 6 we proposed two different, MGD inspired, algorithms that allow us to obtain new analytical solutions in isotropic coordinates, with anisotropy in the pressures. In order to verify the usefulness of the methods, we extend the Gold III solution to the anisotropical domain with both methods and using two different conditions. The physical acceptability conditions of the four new and exact anisotropic solutions are also discussed in Section 7.

II Einstein’s equations in standard coordinates

The line element in standard coordinates, also known as Schwarzschild like coordinates, is given by

[TABLE]

This metric must be a solution of Einstein’s equations

[TABLE]

which take the following form

[TABLE]

where the prime indicates derivatives respect to variable $r$.

Finally the conservation equation, derived from the latter system of equations is

[TABLE]

whose radial component lead to the equilibrium equation of the matter distribution

[TABLE]

III The MGD Approach

The starting point of this method is to assume that the energy-momentum tensor has the specific form

[TABLE]

where $\alpha$ is a coupling constant and $T^{\rm(PF)}$ is the matter-energy content associated to a perfect fluid

[TABLE]

with the fluid 4-velocity given by $u^{\mu}=e^{-\nu/2}\,\delta_{0}^{\mu}$.

Now we can consider a perfect fluid solution of Einstein’s equations ($\alpha=0$) with a line element written in standard coordinates as (1), and define

[TABLE]

which is the standard expression for the mass function in GR. The next step is to take account of the anisotropy introduced by the gravitational source $\theta_{\mu\nu}$ in our system. This will be done by assuming that the contribution of the $\alpha$ parameter, in the perfect fluid solution (1), is encoded in the deformations $g$ and $f$ of the temporal and radial metric component respectively.

[TABLE]

The Minimal Geometric Deformation will be associated to the specific case where $g=0$. In this case ($g=0$) is easy to check that, using (8) and (12), Einstein’s equations (3)-(5) splits in two systems. The first one coincides with Einstein’s equations system for a perfect fluid

[TABLE]

with the correspondent conservation equation

[TABLE]

that turns out to be equation (6) with the energy momentum tensor associated to a perfect fluid (9)(eq (8) with $\alpha=0$).

The second system of equations, which is not an Einstein’s equations system due to the lackness of a $\frac{1}{r}$ on the right side of the first two equations, reads

[TABLE]

and the conservation equation associated with the source is

[TABLE]

In order to find the solution of Einstein’s equations for an energy-momentum tensor of the form (8) we have to solve the equations systems (13)-(15) and (17)-(19). In the case when we start with a known perfect fluid solution, then is only necessary to solve the second system. Now, in both cases there are more unknown functions that equations so additional information is required in order to solve the systems. This information can be given in the form of equation of state or other expressions that relates the physical variables of the system under study.

Finally it can be seen that redefining an energy momentum tensor $\tilde{\theta}_{\mu\nu}$, the system of equations (17)-(19) associated to the source, is equivalent to a Einstein’s equations system for an anisotropic fluid. Also, it can be verified that the conservation equations (16) and (20), implies that the interaction of the perfect fluid with the source is purely gravitational.

IV Einstein’s equations in isotropic coordinates

In this paper we shall study the interior of static and spherical symmetric matter distributions, but using the line element in isotropic coordinates given by

[TABLE]

in which Einstein’s equations take the following form

[TABLE]

This coordinates seems to be more general than the Schwarzschild like ones (1), due to the fact that there is always possible to transform the line element from the standard form to the isotropic one by

[TABLE]

where $\mu^{-1}$ and $r_{1}$ are the radial component and the coordinate associated to the line element in standard coordinates (1). However, is not always possible to perform the reverse process. Therefore, there is a chance to obtain solutions to Einstein’s equations that can not be found using the Schwarzschild like coordinates.

If we consider that $\omega=(\tilde{A}(r))^{2}$, the eqs (22)-(24) reads

[TABLE]

At this point is easy to see that the system (22)-(24) (or (26)-(28)) will not be decoupled if we choose an energy momentum tensor of the form (8) and consider the particular ansatz for the metric (12) of the MGD method . Nevertheless, as we will show in next sections, this conditions can be used to obtain new internal analytical solutions (physically acceptable) of Einstein’s equations in isotropic coordinates.

Finally, in order to avoid the appearance of singular behaviour of the physical variables on the surface of our distribution, we must impose the well known matching conditions between the interior and the exterior space-time geometries. The inner region is defined by the metric (21) and we will consider that the outer region is described by the vacuum Schwarzschild solution

[TABLE]

where M denote the total mass of the distributions.

Now, using (21) and (29) the matching condition takes the following form

[TABLE]

were the subscript $\Sigma$ indicates that the quantity is evaluated at the boundary of the distribution. Then its possible to obtain an expression form the total mass of the distribution given by

[TABLE]

V Physical acceptability conditions

Now, solving Einstein’s equations does not ensure that the solution will describe any physical system. Indeed, among all the known solutions of Einstein’s equations, only a part of them fulfill the physically acceptable conditions (see for example Delgaty ).

Then, in order to ensure that the solutions of Einstein’s equations are physically acceptable, they must satisfy the following conditions

• •

$P_{r}$, $P_{t}$ and $\rho$ are positive and finite inside the distribution.

• •

$\frac{dP_{r}}{dr}$, $\frac{dP_{t}}{dr}$ and $\frac{d\rho}{dr}$ are monotonically decreasing.

• •

Dominant energy condition: $\frac{P_{r}}{\rho}\leq 1$  ,   $\frac{P_{t}}{\rho}$ $\leq 1$.

• •

Causality condition: $0<\frac{dP_{r}}{d\rho}<1$ ,  $0<\frac{dP_{t}}{d\rho}<1$.

• •

The local anisotropy of the distribution should be zero at the center and increasing towards the surface.

VI The algorithms for obtain new solutions

In this section we present two possibles and inequivalent algorithms to find new analytical anisotropic solutions of Einstein’s equations in isotropic coordinates based in the conditions (8) and (12) of the MGD method.

VI.1 The first algorithm

In order to solve Einstein’s equations we will consider an energy-momentum tensor of the form (8). In this case, as the geometry remains untouched, Einstein’s equations reads as before

[TABLE]

but with an effective radial pressure $\tilde{p}_{r}$, effective tangential pressure $\tilde{p}_{t}$ and effective density $\tilde{\rho}$ that are given by

[TABLE]

At this point it is easy to see that we have an anisotropic system with three independent equations and seven unknown functions, the same 4 as the isotropic perfect fluid case (13)-(15) plus $\theta_{0}^{0}$ , $\theta_{1}^{1}$ , $\theta_{2}^{2}$.

Now, as we mention before, if we try to consider the deformation

[TABLE]

is not possible to decouple the system of equations (33)-(35) due to the presence of $\alpha^{2}$ terms in $\omega^{\prime 2}$ and $\omega^{\prime\prime}$.

However, if we consider a specific combination of equations (33)-(35) given by

[TABLE]

and we introduce the minimal geometric deformation

[TABLE]

we obtained that this equation splits in two: one equation associated to a perfect fluid and another one related to the gravitational source

[TABLE]

respectively.

If we now impose a linear constraint of the form

[TABLE]

we have from (41) that

[TABLE]

where $C$ is a constant. Although it is not possible to decouple the system of equations (33)-(35) after we introduce the deformation, is remarkable that knowing $f$, $f^{\prime}$ and $f^{\prime\prime}$, we can obtain $\theta_{0}^{0}$,$\theta_{1}^{1}$ and $\theta_{2}^{2}$ for each value of $\alpha$, from equations (33), (34) and (35) respectively. In fact

[TABLE]

where we have used eqs (22)-(24), and

[TABLE]

with

[TABLE]

Then, knowing $\theta_{0}^{0}$, $\theta_{1}^{1}$ and $\theta_{2}^{2}$, we can write the expressions for the effective pressures and energy density using (36). This represent an analytical and (in general) anisotropic solution of Einstein’s Equations in isotropic coordinates, different from the original one.

VI.2 The second algorithm

Starting with the line element in isotropic coordinate given by

[TABLE]

and considering that the energy momentum tensor has the form

[TABLE]

where $T^{PF}_{\mu\nu}$ is the energy momentum tensor of a perfect fluid. Then Einstein’s equations leads to the system

[TABLE]

where

[TABLE]

Now assuming that the system with $\alpha\not=0$ is characterized by

[TABLE]

and supposing that the sets of functions $\{\nu,A,\rho,P\}$ and $\{\nu,f,H^{0}_{0},H^{1}_{1},H^{2}_{2}\}$ are solutions of Einstein’s equations we have that the above equation system is decomposed in the following systems: the first

[TABLE]

the second

[TABLE]

which are both systems of Einstein’s equations, and a third one

[TABLE]

which is not.

Now in order to solve the systems and find the deformation function $f$ we could propose any relation of the components of the sources. In the MGD method the solutions were found imposing some simple mimic constrains. In this case, a simple possibility to solve the system is to assume

[TABLE]

which leads to

[TABLE]

whose solution is

[TABLE]

where $C$ is an integration constant. Then, if we choose an known solution of Einstein’s equations in isotropic coordinates is possible to obtain new analytical solutions with local anisotropy in pressures.

VII New exact anisotropic solutions of Einstein’s equations in isotropic coordinates

VII.1 Using the first method

In order to verify the first algorithm, let us choose a known perfect fluid solution in isotropic coordinates as Gold III Gold solution

[TABLE]

where $D$, $B$, $a$ and $b$ are constants. This solution will be regular at the origin only if

[TABLE]

nevertheless this equation will not be necessary satisfied in the new solutions obtained from this one.

Solution 1

If we consider $F(r)=0$, for simplicity, it can be check that the deformation function takes the form

[TABLE]

with $C$ an integration constant. Therefore

[TABLE]

where

[TABLE]

and

[TABLE]

We can write now the effective energy density and pressures as (36).

Now, is easy to see that for $a=2$, $b=1$ and $\tilde{\alpha}=D\alpha=0.05$, we found that $r_{\Sigma}=0.885$, $r_{1\Sigma}=0.7853$, $D=0.881$ and $M=0.126$. The pressures and the energy density are plotted in figures (1)-(2) while the acceptability condition are plotted in figures (3)-(6).

Solution 2

The latter solution can be generalized to the case where $F(r)\neq 0$. In particular, if we choose

[TABLE]

with $G$ a proportionality constant, then

[TABLE]

and it can be seen that

[TABLE]

where $\tilde{C}=2K+C$ and $K$ is an integration constant. Is easy to see that it contains the solution 1 ($K=G=0$).

It can be check that

[TABLE]

Now, from the matching conditions, is easy to see that for$a=2$, $b=1$ and $\tilde{\alpha}=D\alpha=0.1$, we found that $r_{\Sigma}=0.5955$, $r_{1\Sigma}=0.6527$ , $M=0.07$ and $D=1.1496$. The pressures, energy density and acceptability conditions are plotted in figures (7)-(12).

VII.2 Using the second method

In order to check the method, we choose as seed solution the Gold III solution again, which is given by eqs (74)-(78) where is easy to see that

[TABLE]

Now, using these expressions and the differential equation for $f$ (Eq. (72)) we can rewrite the systems (65)-(67) and (68)-(69) as

[TABLE]

and

[TABLE]

where

[TABLE]

Solution 1

Let us now study the most simple case were $\tilde{F}=0$ in which

[TABLE]

where we can take, without losing generality, $C=1$.

Now, introducing (107) in the systems (98)-(100) and (101)-(103) , we can compute the effective energy density and pressures

[TABLE]

Finally in order to give an example of a physical acceptable solution we can take the following values for the constants, $\alpha D=0.020$, $a=2$ and $b=0.4$. Then from the matching conditions we obtain that $r_{\Sigma}=1.270$, $r_{1\Sigma}=1.142$, $M=0.175$, $D=0.922$ and $\alpha=0.022$.The pressures can be found in the figure (13) while the energy density is in the figure (14). All the acceptability condition are plotted in the figures (15)-(18).

Solution 2

In order to see if we can improve the behavior of the previous solution let us now take a more complicated case in which $\tilde{F}\not=0$. In particular we will suppose that

[TABLE]

where $G$ is a constant that allow us to connect with the previous case ($\tilde{F}=0$). Then we can obtain

[TABLE]

Now, in this case the regularity in the center of the distribution implies

[TABLE]

As in the previous case, introducing (111)-(113) in the systems (98)-(100) and (101)-(103) we can obtain the effective pressures and energy density from

[TABLE]

Finally, to give an example we choose the values $G=D$, $D\alpha=0.304$, $a=3.3$, $b=1.02$ and $C=0$. Then with the matching conditions it is easy to obtain $r_{\Sigma}=0.952$, $r_{1\Sigma}=0.913$, $M=0.047$, $D=0.953$ and $\alpha=0.319$.In this case the pressures and the energy density are in the figures (19) and (20), while the acceptability conditions are in the figures (21)-(24).

VIII Conclusions

The Minimal Geometric Deformation is an efficient method to obtain new analytical internal solutions to Einstein’s equations. However, the main advantage of MGD is not present if one considers a line element in isotropic coordinates, because it is not possible to decouple the system of equations related to the gravitational sources.

In this work, we present two inequivalent methods to obtain analytical and anisotropical solutions of Einstein’s equations in isotropic coordinates, starting with a known (isotropic or anisotropic in pressures) solution. Although this is not exactly the known Minimal Geometric Deformation, both algorithms uses the same ansatz for the metric deformations. The differences with MGD is the following: in the first method there is no decoupling of Einstein’s equations and we only use convenient combinations of the equations to solve the complete system. In the second method, there is a decoupling of the system of equations. Nevertheless, to achieve this, is necessary to include an additional source ($H_{\mu\nu}$) that also have to satisfy Einstein’s equations.

In order to test the two methods we select Gold III as a seed solution and with each algorithm we extend it to the anisotropical domain using two different conditions over the new sources. Indeed, we have found four analytical solutions in isotropic coordinates, completely different from the initial ones, and with local anisotropic pressures. It is important to mention that the inequivalence of both methods is evident. This can be easily verified by selecting the same seed solution and the same values for the free parameters. This way, it can be verified that the resulting solutions are different. In fact, selecting the values for the parameters used in the second method in order to obtain a physical acceptable solution, the first method throws a solution that does not satisfy the conditions for physical acceptability.

As we can observe from figures (1) and (13), the solution obtained from Gold III with $F=0$ in the first method and $\tilde{F}=0$ in the second one exhibit regions, in the internal part of the distributions, in which the radial pressure is greater than the tangential pressure. This behavior could suggest that this models are unstable under perturbations. In fact the change of sign in the force contribution coming from the local anisotropy indicates the possible appearance of cracking (see cn ; cn1 ; cn2 ; cn3 ; cn4 ; cn5 ; cn6 ; cn7 ) under perturbations. It had been widely discussed in the literature that the presence of cracking can affect the evolution of compact objects in a drastically way when these depart from the equilibrium. Indeed cracking phenomenon can be associated to different scenarios that can be related with observations (for example 33 ; 34 ; 36 ; 37 ; 38 ; 39 ; 40 ; 41 ; 42 ).

However, in both methods we show that this behavior can be corrected when $F\not=0$ in the first one and $\tilde{F}\not=0$ in the second one. Indeed, we show that the solutions displayed in the figures (9)-(12) and (19)-(24) satisfied all the physical acceptability conditions. Another possibility to avoid the odd behavior of the solutions obtained here with $F=0$ ( $\tilde{F}=0$ in the second method), which we do not discuss, consist in doing the matching with an exterior solution different from the Schwarzschild’s vacuum.

Finally, we what to mention that the complicate form that solutions have, suggest that this may be new solutions of the Einstein’s equations not reported in the literature, at least compared with other solutions in isotropic coordinates. Nevertheless, to probe that these are indeed new solutions of the field equations (including the rewriting in any other coordinates) it is necessary to perform a deeper analysis (computing the scalar invariants of our solutions and comparing them with those obtained from all the others known solutions of the Einstein’s equation’s) which is beyond the purpose of this work. Our main goal in this work is not to study any particular model. We only center our attention in presenting two different simple algorithms to find analytical solutions to Einstein’s equations in isotropic coordinates. In order to describe any particular model, it is necessary to consider $F(r)$ (or $\widetilde{F}(r)$) related with the system under study, in order to solve the resulting equations.

IX Acknowledgements

We want to say thanks for the financial help received by the Projects ANT1855 and ANT1856 of the Universidad de Antofagasta. P.L wants to say thanks for the financial support received by the CONICYT PFCHA / DOCTORADO BECAS CHILE/2019 - 21190517 and also to Professor M. Asorey of the Physics Department of the Universidad de Zaragoza where part of this work was done. C.L.H wants to say thanks for the financial support received by CONICYT PFCHA / DOCTORADO BECAS CHILE/2019 - 21190263, and to Professor B. Fiol from the Physics Department of the Universidad de Barcelona for his hospitality during the realization of part of this work. P.L and C.L.H are also grateful with Project Fondecyt Regular 1161192 and Semillero de Investigación SEM 18-02 from Universidad de Antofagasta.