Abstract

In this manuscript, we study the relativistic quantum mechanics of an electron in external fields in the spinning cosmic string spacetime. We obtain the Dirac equation and write the first- and second-order equations from it, and then, we solve these equations for bound states. We show that there are bound state solutions for the first-order equation Dirac. For the second-order equation, we obtain the corresponding wave functions, which depend on the Kummer functions. Then, we determine the energies of the particle. We examine the behavior of the energies as a function of the physical parameters of the model, such as rotation, curvature, magnetic field, Aharonov-Bohm flux, and quantum numbers. We find that, depending on the values of these parameters, there are energy nonpermissible levels.

1. Introduction

Symmetry is key ingredient in the description of natural phenomena. The notion of symmetry is an essential feature in several areas of physics. In this context, the well-known Noether’s theorem [1] establishes a connection between symmetry and conservation laws of relevant physical quantities. In quantum mechanics, we often use symmetry to obtain crucial results concerning angular momenta operators [2]. Likewise, symmetry is relevant in the topic of quantum information [3]. Also, symmetry is essential in the framework of relativity [4], and for this reason, it is indispensable in research areas such as particle physics [5] and cosmology [6].

A pertinent question in the research areas cited above refers to think about the consequences of symmetry-breaking in a given physical system. It can happen in several contexts, like in phase transitions, for example [79]. Another example of symmetry-breaking occurs in solids: we can create topological defects like disclinations and dislocations [10, 11] from the Volterra process [12].

Topological defects can emerge in a large number of physical systems covering themes such as liquid crystals [13], graphene physics [14, 15], magnetism [16], and cosmology [17]. Recent studies also have reported the importance of topological defects in Life Sciences [18, 19]. In cosmology, defects in the spacetime topology can be viewed as a possible consequence of the evolution of the early universe, which has suffered phase transitions due to the temperature decreasing and the process of expansion [20, 21].

In this contribution, we are involved in studying the topological defect known as a cosmic string. A cosmic string is a linear defect, similar to a flux tube in type-II superconductors [20]. The spacetime around such defect has a conical symmetry, being identical to the case of a disclination [22]. The concept of a cosmic string was introduced in the literature by Kibble [23]. An intriguing facet in this subject refers to the quantum mechanical description of a particle in a region of the spacetime containing this defect. There are several works dealing with this issue, both in the nonrelativistic and relativistic scenarios. Reference [24], for instance, analyzes the hydrogen atom in a spacetime of a cosmic string. Reference [25] considers the problem of a relativistic electron in the presence of both Coulomb and scalar potentials in the cosmic string spacetime. Results about vacuum polarization in a cosmic string spacetime were reported in Ref. [26]. Again, the cosmic string spacetime was considered as a background to examine relativistic oscillators [2730], quantum phases [31], and fermionic currents [32].

A relevant issue in this context consists of taking into consideration the influence of electromagnetic fields in the quantum particle motion. Landau levels [33] and the Aharonov-Bohm effect [34, 35], for instance, are essential ingredients in the investigation of quantum systems even in a flat spacetime. Landau levels are a quantum analog of classical cyclotron motion, while the Aharonov-Bohm effect reveals the significance of the vector potential in the quantum world. Then, studying the contribution of these magnetic interactions to the quantum mechanical description of a system in spacetime having a topological defect is a natural development. In Ref. [36], for example, it was considered the interaction of a point charge with a magnetic field in a spacetime with a distortion. References [3744] are examples of studies dealing with Landau levels and the Aharonov-Bohm effect in the presence of topological defects. Reference [45] deals with a quantum ring in graphene with a topological defect and a magnetic field. Reference [46] investigates an analog of the Aharonov-Bohm effect for bound states in a spacetime containing a spacelike dislocation. In particular, the inclusion of electromagnetic interactions in the case of a cosmic string background also has been considered. For instance, Ref. [47] analyzes the quantum dynamics of a charged particle in the presence of a magnetic field and scalar potential. Reference [48] examines different configurations of confined magnetic fields and investigates the existence of induced vacuum fermionic currents.

On the other hand, we can be interested in analyzing the behavior of rotating systems. It is a pertinent question, since large structures in the Universe, such as black holes and galaxies, have angular momentum due to rotation [49]. Also, it is possible to conceive cosmologic models with rotation [50]. The study of quantum systems in spinning spacetimes also has been reported in the literature [5153]. Noninertial effects in quantum systems can provide novel theoretical predictions and feasible experimental developments. For instance, quantum phases in rotating systems can emerge, in analogy to the Aharonov-Bohm effect [54, 55]. It is also possible to establish a relation between the Hall effect and the inertial forces [56]. Besides, if a given system is rotating, physical properties like spin transport [57, 58] and electronic structure [59] are modified. A system can present magnetization due to rotation, like in the Barnett effect [60]. While a magnetic field produces a spin-field coupling, resulting in the anomalous Zeeman effect [61], rotation produces an analog effect, due to the spin-rotation coupling [62]. Thus, rotation can contribute similarly to a magnetic field in the dynamics of a quantum system. More, noninertial effects are an interesting issue when the spacetime contains topological defects. In this case, we can include noninertial effects and the topological defects in the quantum mechanical description by employing the same tool. More specifically, we can use a metric tensor to a spinning spacetime with a topological defect [63]. The spacetime of a spinning cosmic string has been considered as background for several problems involving quantum systems [64]. For instance, Ref. [65] deals with the Schrödinger equation in that spacetime. Reference [66] analyzed the bound states for neutral particles in a rotating frame of a cosmic string. Likewise, Ref. [67] investigates rotating effects on a Landau-Aharonov-Casher System in the spacetime of a cosmic string.

As we already have mentioned, rotation can present similarities to the electromagnetic fields. This way, it is also an attractive question examining how the electromagnetic interactions affect the particle quantum motion of a rotating system in the presence of a topological defect. Reference [68], for instance, is a recent work dealing with both topological and noninertial effects in the presence of an Aharonov-Bohm potential. Reference [69] addresses the problem of a spinless relativistic particle subjected to a uniform magnetic field in the spinning cosmic string spacetime. The Dirac oscillator in the spacetime of a cosmic string in the presence of the Aharonov-Casher effect was analyzed in Ref. [70]. The presence of noninertial effects also was considered in this context [71]. Reference [72] treats the problem of a charged half-spin particle depicted by the Dirac equation in the presence of a uniform magnetic field in the rotating cosmic string spacetime. A meaningful aspect in this context consists of analyzing how different configurations of magnetic fields affect the quantum particle motion. In particular, in this contribution, we choose to deal with a superposition of two different configurations, due to their fundamental aspects: a uniform magnetic field, which is related to the raising of Landau levels in quantum systems, and an Aharonov-Bohm flux, since it provides a new significance to the role of the electromagnetic interactions in the quantum theory, which can be manifested even for bound states [73]. Besides, several analogs of the AB effect can emerge in spacetimes with topological defects [7476].

In this paper, we study the relativistic quantum mechanics of an electron in the presence of both a uniform magnetic field and Aharonov-Bohm potential in the spinning cosmic string spacetime. In other words, we solve the Dirac equation in this scenario and investigate how the rotation, curvature, and external magnetic fields affect the wave functions and energies of the electron. Although the existence of a previous study deals with the Dirac equation in the presence of noninertial effects and uniform magnetic field, in Section 3 of Ref. [72], here, we introduce some new ingredients, as well as a different framework: more explicitly, we have the following contributions: (i) the tetrad basis we adopt leads to a zero affine connection in the flat space, which does not occur in the previous study; (ii) the presence of the Aharonov-Bohm flux tube. We show explicitly that the presence of this effect modifies the energy levels of the particle; (iii) we also study first-order solutions for the Dirac equation and find particular bound state solutions for it; (iv) we employ an ansatz for the solutions whose spinor is an eigenstate of the momentum angular operator [47].

The manuscript is organized as follows. In Section 2, we present some algebraic elements necessary to construct the field equations in curved spacetime and write the Dirac equation describing the quantum motion of the electron in the presence of external magnetic fields in the spinning cosmic string background. In Section 3, we deal with first-order solutions and study the existence of isolated solutions for the particular case of a particle at rest. In Section 4, we take our attention to the case when the energy of the particle is different from its rest energy. We map the Dirac equation problem in curved space with minimal coupling into a Sturm-Liouville problem for the upper component of the Dirac spinor, and using an appropriate ansatz, we derive the radial equation. We solve the radial equation and find the wave functions and energies of the particle. We make a detailed discussion of the results and also comparisons with other studies in the literature. In Section 5, we present our conclusions. In our work, we use natural units, .

2. Dirac Equation in the Spinning Cosmic String Spacetime

In this section, we briefly present the tools needed to construct the Dirac equation in the conical spacetime in the presence of noninertial effects. The first step consists of taking a look at the metric tensor characterizing this geometry. Next, we will choose an appropriate tetrad basis and implement the field configuration involved through the performing of a minimal substitution. The spacetime induced by a rotating cosmic string is described by the metric [72, 77, 78]. where , , and . The parameter is related to the linear mass density of the cosmic string through the relation , and it runs in the interval . The quantity is the rotation parameter, with representing the angular momentum of the spinning cosmic string. A remarkable feature of the spinning cosmic string spacetime is related to the appearance of closed timelike curves (CTCs) [7981]. Such structures can take place depending on the relation between the radial coordinate, the rotation parameter, and the curvature parameter. Writing the metric (1) in the form such a structure appears if . CTCs are one of the most controversial aspects of modern physics [82]. It is possible to find in the literature works dealing with the concepts of causality [83], chronology protection [84], time travel [85], and flow information [86] in the framework of closed timelike curves. In our study, we consider the quantum motion only in regions such that . More explicitly, we will consider the following restriction on the radial coordinate:

This type of restriction on the radial coordinate also occurs, for example, in the case of a spacetime containing a time dislocation, as pointed in Ref. [51].

The relativistic quantum dynamics of a spin-1 particle interacting with external magnetic fields in the rotating cosmic string spacetime is governed by the Dirac equation where is the mass of the particle and is the covariant derivative for fermion fields defined as and are the Dirac matrices in the rotating cosmic string spacetime, which are defined in terms of the tetrad fields and Dirac matrices in the flat space in the following way: where are the standard Dirac matrices and are the usual Pauli matrices. The matrices of (6) satisfy the following relation:

Also, in Equation (5), is the spin affine connection given by where are the Christoffel symbols of the second kind and is the tetrad field. The tetrad basis satisfies the relations

Here, the Greek letters are used for tensor indices while the Latin letters are denoting Minkowski indices. We use the tetrad basis and its inverse defined as [77]

For this choice, it can be shown that the nonvanishing affine connection is given by where

By using the tetrad basis (14), the matrices of (6) can be written explicitly as with being the Pauli matrices in the curved spacetime.

Since we are first interested in studying the solutions of the Dirac equation in its present form (Equation (4)), we need to write the corresponding system of first-order coupled differential equations. Let us assume the time-dependence of the wave functions together with the decomposition of the fermion field in the form with

The system we will analyze takes into account the particle is immersed in a region where there is a uniform magnetic field and also the potential due to a thin long solenoid along the -axis. Having this field configuration in mind, we study the physical implications due to noninertial effects and the Aharonov-Bohm potential on the relativistic Landau quantization. We also take into account the translational invariance of the system along the -direction, which allows us to eliminate the third direction (), and, consequently, we can consider only the planar motion [8790]. Here, the four-potential of the electromagnetic field in Equation (4) has only the spatial component, , and the potential vector in the Coulomb gauge is specified by with where is the magnetic field magnitude, , is the magnetic flux, and is the quantum of magnetic flux along the solenoid. This configuration provides an superposition of magnetic fields in the -direction with

Note that the particle only interacts with the magnetic field due to the potential vector . Here, we are focused on studying the electron motion only in the region, so that we can neglect the point interaction and, consequently, consider only regular wave functions.

Using the results above, the Dirac equation (4) can be written as

At this point, we are ready to solve Equations (28) and (29) by considering two distinct circumstances: (i)Take our attention to isolated solutions of the first-order Dirac equation by imposing the condition (ii)By imposing the condition , we search the solutions of the second-order Dirac equation

We will show in the next two sections that there are bound state solutions for both cases and discuss their main physical properties. To distinguish each case in (i), in the next section, we use the superscripts () to label the quantities corresponding to .

3. Solution of the Equation of Motion to

To study the existence of isolated solutions of the Dirac equation (4), we must set in Equations (28) and (29). In literature, such solutions are known to be excluded from the Sturm-Liouville problem. The search for isolated solutions of the Dirac equation has been performed in different physical contexts [9195]. Since condition (3) establishes a lower limit to the radial coordinate, given by , the bound state solution must satisfy the normalization condition where , which using the ansatz (20) and the matrices (18) and (17), takes the form

Equations (30) and (31) must be finite to ensure the existence of bound state solutions. By making in Equations (28) and (29) and using Equations (21) and (22), we get which using Equation (25) can be written as with

The solution of the coupled linear differential equation systems (36)–(39) is given by with where are upper incomplete Gamma functions [96], and , , , and are constants. Analyzing the solutions (42) and (44), we note that dominates over for any value of , in such way both solutions converge when and . This will not occur for the function in the solutions (43) and (45). Moreover, since the incomplete Gamma functions and always diverge, then the function will only converge when if while the function will always diverge when and, therefore, will not be a square-integrable function. Thus, the only allowed solution for the equation systems (36)–(39) results with , and the requirement . Solution (49) satisfies Equation (31) and constitutes a bound state solution for the case , i.e., an isolated solution to the Dirac equation (4) in the metric spacetime (1). Note there is another particular solution for this case when we reverse the signal of the magnetic field or the charge, i.e., . In this case, the only nonzero component of the normalizable solution is (Equation (43) with ).

Proceeding in an analogous way, now we make in Equations (28) and (29). We find the system of equations with

The solution of the coupled linear ordinary differential equation systems (51)–(54) is given by with where

By making the same analysis of the solutions as we have made for the case , i.e., analyzing the behavior of the functions for , we find that the only solution that admits bound states is the (59) one. Thus, the appropriate solution for the case satisfying the normalization condition (31) is given by with and . Note that the solutions (49) and (63) are affected by rotation through Equations (40) and (55), respectively. Similarly to the previous case, if , the only nonzero component of the normalizable solution is (Equation (60) with )

Moreover, from the discussion of the solutions for each case above, we immediately see that only the first line of the normalization condition (31) is relevant in the current study.

4. Solution of the Equation of Motion to

In this section, we solve the second-order equation to that we find from Equations (28) and (29). The solution of this equation is different from that one calculated in the previous section and allows us to obtain an expression for the particle energies. By isolating in Equation (29) and replacing in Equation (28), we are able to write the second-order differential equation for as

Using the decomposition of the fermion field (21) together with Equations (23)–(25), we obtain the radial equation for where

Equation (66) is the confluent hypergeometric equation, and its solution is well known. Thus, it can be shown that the solution to is where denotes the confluent hypergeometric function of the first kind or Kummer’s function and the normalization constant. It can be shown that the hypergeometric function has a divergent behavior for large values of . Because of this, bound state solutions for Equation (69) are only possible if we impose that this function becomes a polynomial of degree . For this to be accomplished, we require that , where , with denoting the set of the nonnegative integers. However, as we can see in Equation (68), the absolute value of the effective angular moment is defined in terms of the energy . In this way, to obtain the energy eigenvalues from the above condition, we must consider and , respectively, and then solve them for . By making this, we get with the following requirement: where

In Equations (70) and (71), the superscripts () refer to the energies calculated for and , respectively. A similar equation for can be immediately obtained by making and in Equation (66). We find where

The function is found by solving Equation (74). The result is

The energies obtained from (74) are with the following requirement:

where

Note that two other radial equations can be derived from the equation for from Equations (28) and (29). However, there is no need to solve them here because they have both energy and wave functions similar to those obtained above. Let us study the energies (70) and (71). First, we see that the energy depends only on the quantum number and the magnetic field . If we make in Equation (70), we get . As mentioned above, since we are solving the problem with the requirement that , then these energies are excluded in the present case.

Now let us investigate the energy levels we find. In all the energy profiles we address here, we use , . We also use solid lines to represent the energy of the particle and dashed lines for the energy of the antiparticle. In Figure 1, we show the profile of for this particular case, where we see clearly that increases when the magnetic field is increased. The energies (70) and (71) denote the relativistic Landau levels in the present context. These energies can be directly compared with those obtained for the relativistic oscillator (Dirac oscillator) addressed in Ref. [77]. Although that scenario is different from the one we are exploring here, there are similarities between the profiles of the energy levels in both models. For example, the energy (47) of the Ref. [77] depends only on the frequency of the oscillator and the quantum number . In our case, by defining the cyclotron frequency , Equation (70) results which makes such a similarity clear. We can see that the energies (71) and (78) are the only ones that depend on all the physical parameters involved in the current problem. However, it is possible to verify that such spectra are similar. Because of this, we analyze in detail only the energy (71). One of the objectives in our analysis is to show that the combined effect of the external magnetic field, the Aharonov-Bohm flux, the rotation, and the curvature modifies the energy spectrum of the particle. As mentioned in Section 1, one of the new ingredients in the present model is the inclusion of the Aharonov-Bohm flux tube. Then, let us focus on the physical implications of this. We can access varied forms of energy profiles by attributing different values to the parameters of the model. Figure 2 shows the profile of as a function of by considering four different parameter choices. In these four configurations, we consider and . For , and , and the energies are permissible only for values of flux larger than (blue line), (red line), and (green line) (Figure 2(a)). Each value of flux corresponds to the energies , , and , respectively. The values of are chosen in such a way the separation between the energies becomes more notable. When this configuration is changed to , and (Figure 2(b)), the energy levels are more spaced out and a gap between the particle and antiparticle energies emerges. By increasing the magnetic field to and maintaining the other parameters, we notice that the antiparticle energy changes faster than the energy of the particle and the spacing between the energy levels decreases (Figure 2(c)). For increasing values, these effects tend to increase. In the fourth configuration, we use , and (Figure 2(d)). When approaches , we can see the energy is only permissible for values of larger than while and are defined for any value of . In the four settings displayed in Figure 2, we note that increases when and are increased and the other parameters are kept fixed.

Another interesting configuration is the behavior of as a function of for , , , , and some values of (Figure 3). This configuration shows that when is close to , the energy is defined to all the values of in the considered range. However, for small values of , a range of where the energy is nonpermissible arises. In this process, the most pronounced variations occur in the antiparticle energy.

We can also observe interesting effects in the profile of as a function of . For this case, we consider , , , and and investigate it for two sets of values of as illustrated in Figure 4. By carefully studying this configuration, we can observe several effects for values of smaller than (Figure 4(a)). For (blue line) and (red line), the energy tends to decrease when is increased while for (green line), the energy decreases when is increased to . At this value, an overlap between particle and antiparticle energy occurs. The energy increases for values of higher than . For (brown line), the energy of the particle tends to decrease when is increased to and then the behavior remains approximately stable. When we look at the energy of the antiparticle, we observe the opposite process. By increasing to , the energy remains practically stable and then tends to increase as is increased. When is changed to (orange line), both the particle and antiparticle energies increase when is increased. However, the variation on the antiparticle energy is larger. On the other hand, when we consider different values of higher than , the profiles are similar (Figure 4(b)). In this case, increases when both and are increased.

In Figure 5, we investigate the profile of as a function of for , , , and different values of by considering the cases where there is no Aharonov-Bohm flux, and also, considering it is present. In Figure 5(a), we consider . We can see that the presence of Aharonov-Bohm flux reveals that the energies are well defined for the considered values. In the absence of Aharonov-Bohm flux, there are intervals of associated with nonpermissible energies (Figure 5(b)). For values of close to , the changes in the energies are not too evident, while for smaller values of , the modifications in the energy are more evident, being the antiparticle energy the most affected. In both cases, the energy increases when increases, and decreases.

It is also important to examine the energy as a function of , for some values of . Figure 6 illustrates this situation. For , , , and , there are intervals of leading to forbidden energies. This effect becomes more notable to larger values of . Also, the energy increases when increases. It is also interesting to explore the energy profiles by performing energy sketches as a function of two quantities of the model. For a given set of fixed parameters, for example, , , , and , we have the profile of the energy levels as a function of and (Figure 7). We can clearly see that increases with and . The green solid bars denote the discrete energy values for a given and . On the other hand, when we investigate the behavior of (71) as a function of and for specific values of the other parameters, we see that the negative spectrum changes more rapidly when compared with the positive one (Figure 8). In the positive spectrum, both rotation and curvature lead to a linear change, except in the region with and arbitrary . In the negative spectrum, we see that the curvature effects are more predominant in the region where has values smaller than . In this region, any variation in the rotation parameter implies in an abrupt change in the energy spectrum. Modifications in the energies for small values of are an expected manifestation in our analysis. Its physical implication is inherent in metric (1) and is an immediate consequence of the topological cone, since the curvature becomes larger for smaller values of .

When we study the profile of as a function of and , besides the characteristics displayed in Figures 2 and 3, we also see modifications due to noninertial effects (Figure 9). For this particular case, we see that for small values of and , some energy levels are forbidden. This characteristic is observed for , , and (Figure 9(a)). By changing to , this effect disappears, and the energies are set to all values of and in the considered range (Figure 9(b)).

To complete our analysis, we investigate the profile of the energy (71) as a function of magnetic field and the magnetic flux through the solenoid, (Figure 10). Similarly to Figure 8, by fixing the other parameters, we see that the energy of the antiparticle varies more rapidly when compared to the energy of the particle (Figure 10). Clearly, we observe that the energy of the particle varies very slowly throughout the region of flux and magnetic field. As a final commentary, we clarify that the cases discussed in Figures 7, 8, and 10 can be investigated for other fixed parameter values. In this way, it can be shown that there are forbidden energies, depending on the values of the parameters considered. In general, this occurs when both the parameter and the rotation parameter are smaller than and the other parameters assuming higher values than those we use here.

5. Conclusions

In the present manuscript, we have addressed the problem of the relativistic quantum motion of an electron in the spinning cosmic string background considering the presence of a uniform magnetic field and the Aharonov-Bohm potential. We have shown that this combination of potentials allows bound state configurations in the scenario of first-order solutions as well as in the case of second-order solutions of the Dirac equation. It is worth noting the role played by the two different terms in the vector potential. As already known in the literature, we have shown that the uniform field is responsible for a behavior analog to a harmonic oscillator, which leads to the relativistic Landau quantization. Also, we saw the Aharonov-Bohm flux contributes to the angular momentum and energy spectrum of the particle. In the case of first-order solutions, which were obtained by solving Equations (49) and (63) for and , respectively, the oscillator-like behavior provided by the uniform magnetic field guarantees the convergent first-order solutions and, consequently, the existence of bound states. The isolated solutions obtained (Equations (49) and (63)) are particular solutions of the Dirac equation (4).

We have also studied the more general problem by solving the second-order equation implied by Equations (28) and (29) for the upper component of the Dirac spinor for . Using appropriate solutions (Equation (20)), we have derived the radial equation and shown that its solution depends on the Kummer functions. From these functions, we have extracted the expression for the energy levels of the particle (Equations (70) and (71)). For the field configuration considered, we have found that the effective angular momentum of the electron depends on its energy and the Aharonov-Bohm flux tube. Besides, the potential vector corresponding to the uniform field leads to a charged oscillator. Thus, such field superposition provides distinct effects on the motion of the particle. Additionally, in some cases, the rotation produces a combined effect with both the uniform magnetic field and the curvature (see Equation (71)). We have shown that the energy levels of the particle and antiparticle depend on the values of the physical parameters involved. In the case of energy (71), its validity is conditioned to Equation (72). Depending on the choice we make for the parameters, we can obtain forbidden energies. Sketches in Figures 7, 8, and 10 illustrate the profiles of the particle and antiparticle energies and show that they belong to the same spectrum. The effects of curvature and rotation are more evident when , being the antiparticle energy the most affected. Summarizing our results, we can state that the rotation breaks the symmetry in the energies (71) and (78) about and , respectively, while the energies (70) and (77) are symmetrical. When or in (71) and (78), this symmetry tends to be recovered. In the case of a strong magnetic field, the separation between the particle and antiparticle energies increases in comparison with the case of a weak field. When is large and , the energies increase rapidly. We saw also that the inclusion of the Aharonov-Bohm flux can produce notable modifications in the energy profile. For instance, it can affect even the behavior of the energy as a function of the uniform magnetic field. As a final comment, we would like to emphasize that the model studied in this article generalizes others found in the literature, such as those of Refs. [68, 72] for the case including a superposition of external magnetic fields and the investigation of isolated solutions of the Dirac equation. Furthermore, we present a detailed discussion on the energy levels of the particle which, in general, is not found in the literature.

Data Availability

No data were used to support this study.

Disclosure

A previous version of this manuscript was presented in arXiv:2006.06511.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

We would like to thank E.R.B. Mello (Universidade Federal da Paraíba, PB, Brazil) for his remarks and comments. This work was partially supported by the Brazilian agencies CAPES, CNPq, and FAPEMA. EOS acknowledges CNPq Grant 307203/2019-0 and FAPEMA Grant 01852/14. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nvel Superior-Brasil (CAPES)-Finance Code 001. MMC acknowledges CAPES Grant 88887.358036/2019-00.