#### Abstract

We study vortices in generalized Maxwell-Higgs models, with the inclusion of a quadratic kinetic term with the covariant derivative of the scalar field in the Lagrangian density. We discuss the stressless condition and show that the presence of analytical solutions helps us to define the model compatible with the existence of first order equations. A method to decouple the first order equations and to construct the model is then introduced and, as a bonus, we get the energy depending exclusively on a function of the fields calculated from the boundary conditions. We investigate some specific possibilities and find, in particular, a compact vortex configuration in which the energy density is all concentrated in a unit circle.

#### 1. Introduction

Vortices are localized structures that appear in two spatial dimensions. They are present in many areas of nonlinear science and were firstly investigated in the context of fluid mechanics [1, 2]. These objects also appear in type II superconductors [3] when one deals with the Ginzburg-Landau theory of superconductivity [4] and may also be present as magnetic domains in magnetic materials and in many other applications in condensed matter [5, 6].

In high energy physics, in particular, vortices firstly appeared in the Nielsen-Olesen work [7], which is perhaps the simplest relativistic model that supports these structures. The model consists of a Maxwell gauge field minimally coupled to a complex scalar field under the Abelian symmetry in the Minkowski spacetime. An interesting feature of the Nielsen-Olesen vortices is that they are electrically neutral and engender quantized magnetic flux. Their equations of motion are of second order and present couplings between the fields. To simplify the problem, first order equations that solve the equations of motion were found in [8, 9]. In this case, the first and second order equations are only compatible if the potential is of the Higgs type, a potential that engenders spontaneous symmetry breaking. It is worth mentioning that, even with the Bogomol’nyi procedure, the analytical solutions that describe the vortices remain unknown.

Vortices have also been investigated in generalized models with distinct motivations in several works; see, e.g., [10–27]. In particular, k-vortices, which are vortices in models with generalized kinematics, similar to the models studied before in [28–31], were investigated in [12, 13], without the presence of a first order formalism and analytical solutions, but with the search for new effects. Another motivation relies on the possibility of specifying the form of potential, imposed by the first order formalism. For instance, in [23], modifications in the magnetic permeability allowed to develop a route to make the vortex compact. Also, in [27], we have developed a method to obtain vortices and to construct a class of models that supports analytical solutions. Recently, in [32], we have found vortices with internal structure, which arise in generalized models with the magnetic permeability controlled by the addition of a neutral field, enlarging the symmetry to become .

Motivated by several works that appeared with generalized dynamics, we have developed a first order formalism for these models in [26]. This investigation focused on the search for the conditions that could lead to first order equations in a case similar to the one considered before in [12], with the inclusion of a quadratic kinetic term that involves the covariant derivative of the scalar field in the Lagrangian density. In the current work we further explore the subject, extending the previous results of [26, 27] to this much harder class of models. The main results show how the presence of analytical solutions can be used to construct the model, if one imposes that its equations of motion are solved by solutions of first order differential equations compatible with the stressless condition.

Although we are working in the dimensional spacetime with the Minkowski metric, we think that the results of the current work are also of interest to General Relativity (RG), in particular to the case of the so-called Ricci-based theories of gravity (RBG) formulated within the metric-affine approach. For instance, in the recent work [33], the authors unveiled an interesting correspondence between the space of solutions of RBG and RG, under certain circumstances. The results show that it is sometimes possible to map complicated nonlinear models into simpler ones, and we think that the models to be explored in the current work can provide novel possibilities of current interest to the scenario explored in [33, 34].

To study the subject, the work is organized in a way such that in Section 2 we present the model and the procedure, showing the requirements to make it work in the presence of first order equations. In Section 3, we illustrate our findings with some new models that support analytical solutions. In particular, we also calculate the magnetic field, energy density, and total energy of the vortex analytically and investigate the possibility of building compact solutions. Finally, in Section 4 we end the work with some conclusions and an outlook for future investigations.

#### 2. Model and Procedure

We consider the generalized action for a complex scalar field coupled to a gauge field under the local symmetry in a three-dimensional Minkowski spacetime with metric tensor . The Lagrangian density to be investigated has the formIn the above expression, , , and are nonnegative functions that modify the dynamics of the model and is the potential. The minus sign in the term is to keep the vortex energy nonnegative. Also, and define the kinetic terms of the scalar and gauge fields, respectively, aswhere , , and the overline stands for the complex conjugation. The equations of motion for this model arewhere is the conserved current, given by the expression . Also, we are using the notation , etc.

The energy-momentum tensor for the generalized model (1) is We then consider static configurations; take and work with the usual ansatz for vorticesin which and are the polar coordinates and is the vorticity. The boundary conditions for and arewhere is the symmetry breaking parameter which is supposed to be present in the model under investigation. The ansatz (5a) and (5b) makes and be written aswhere the prime denotes the derivative with respect to . The magnetic field is given by . This can be used to show that the magnetic flux is quantized; that is,The ansatz (5a) and (5b) can be plugged in the equations of motion (3a) and (3b), which take the formwhere , etc. The components of the energy-momentum tensor areUp to this point, the scenario is quite similar to the one investigated before in [12]. Here, however, we want to go further and search for a first order framework that helps us to find analytical solutions. We then follow [26] and take the stressless conditions, , which ensure stability of the solution under radial rescaling. This requires the solutions to obey the following first order equationsThey allow us to write . The above equations, however, must be compatible with the equations of motion (10a) and (10b). Similarly to the case that was shown in [26], for and constant, this requirement leads to a constraint that depends on , , and . Therefore, it is hard to obtain a constraint in terms of and reconstruct the model by finding the explicit form of the potential in terms of , , and , as in the case that was carefully investigated in [27]. The main issue appears because does not depend exclusively on , but also on and ; see (8). Nevertheless, if the analytical solutions, as well as their inverses, are known, we may write exclusively in terms of , which we call . By substituting (12a) and (12b) into (10b), the following constraint arisesOne may wonder if the compatibility of (12a) and (12b) with (10a) does not imply another constraint. Nonetheless, as it was demonstrated in [26], once the above constraint is satisfied and the solutions solve (12a), (12b), (10a) becomes an identity. In our model, the choice of the functions , , and must be done in a way that allows the symmetry breaking of the potential to match with the boundary conditions in (6).

The energy density is given by (11a). By using the first order equations (12a) and (12b), it can be written asHere, we follow the procedure developed in [26] and introduce an additional function , defined by where and . By combining the first order equations (12a) and (12b) and the constraint (13), one can show thatIn this case, we can write the energy density as which can be integrated all over the plane to provide the energy

Now, we follow the route suggested in [27] and develop a procedure to build analytical solutions. This can be achieved by decoupling the first order equations (12a) and (12b), as we describe below. For simplicity, we consider dimensionless fields and take ; also, we work with unity vorticity, setting , which means to consider only the upper signs in (12a) and (12b).

In order to decouple the first order equations, we introduce the generating function such thatTherefore, for a given , we can solve the above equation and obtain obeying the boundary conditions (6). By using this into (12a) and (12b) we obtain

We also introduce another function, , which is defined by . By using this and the constraint in (13), we getOne can show that is obtained in terms of the given function from (12b):where is the inverse of . This procedure is valid if is written only as a function of . Using the definition in (20), we findWe can also take advantage of the function to write the magnetic field asand (17) as , which leads to the total energyThis procedure decouples the first order equations in a manner that the solutions depend only on the generating function . As depends only on and , we see from (22a) and (22b) that we have two equations that constrain the functions , , , and . This means that there are several models that support the same analytical solutions defined by (20). Therefore, to find the explicit form of the models, we need to suggest two of the aforementioned functions. Even though these functions lead to the same solutions and magnetic field, they modify the energy density in (14). Thus, one must choose functions that lead to a well defined energy.

We also highlight here that the above procedure to construct the model, described by (22a), (22b), (23), and (24), is only valid in the interval , which is the one where the solution exists, according to the boundary conditions (6). Nonetheless, it is important to suggest nonnegative functions and a potential that supports a minimum at , in order to include spontaneous symmetry breaking and avoid instabilities and negative energies.

#### 3. Specific Examples

Let us now illustrate our procedure with some examples. We firstly suggest an that leads to analytical solutions and then apply the method in (22a), (22b), (23), and (24) to construct the model.

##### 3.1. First Example

The first example arises from the generating functionThis function was previously considered in [27], but with a model in which , which kills the term in the Lagrangian density. By substituting the above expression in (20) and (21) we get the solutionswhich satisfy the boundary conditions (6). The inverse function of the solution in (28), combined with (23) and (24), allows us to writeNotice that these equations and the solutions in (28) are exclusively determined by the function given in (27). This also occurs with the magnetic field, given by (25), which leads toIn Figure 1, we display the solutions (28) and the magnetic field given above. Notice that their behavior is similar to the one for the Nielsen-Olesen case [7, 9].

In order to construct a model that supports the solutions in (28), we use (22a) and (22b). Firstly, though, we need to suggest an explicit form for two of the functions among , , , and . We consider the potentialwhere is a real number. It presents a set of minima at and a local maximum at as illustrated in Figure 2. The other function that we suggest iswhere is a real, nonnegative parameter. The case investigated in [27] is obtained for . By substituting the above and the potential (31) in (22a) and (22b), we obtainIn order to avoid negative coefficients in the above functions, we impose the condition . The functions in (31), (32), (33a), and (33b) determine model (1). We want to emphasize here that this model can only be obtained explicitly because we know the analytical solutions before its construction.

The energy density can be calculated from (14), which leads us toBy a direct integration, one can show that the energy is , which matches with the result obtained by (26). Notice that only the parameter modifies the energy. The above energy density can be seen in Figure 3.

Another model can be generated straightforwardly from the same choice of in (27), which presents well defined , , , and for all .

##### 3.2. Second Example

Here, we consider a generalization of the previous example by considering the generating function to bewhere is a nonnegative real parameter. This function was also investigated in [27], but with . From (20) and (21), we get the analytical solutionswhich satisfy the boundary conditions (6). From the inverse of the solution , combined with (23) and (24), we obtainAs in the previous model, these equations and the solutions in (36) are solely determined by the in (27). The same is valid for the magnetic field in (25), which leads toOne can show that, as increases, the solutions in (36) tend to compactifyand the same happens for the magnetic field in (38), which for very large tends towhere is the Dirac delta function. In Figure 4, we depict the solutions (36) and the magnetic field given above for several values of , including the compact limit in (39a) and (39b).

Again, to find the functions , , , and we must suggest two of them and use (22a) and (22b). We take the potential in the formwhere is a real number. This potential presents minima at for any . The point is a maximum for and a minimum for . This behavior is shown in Figure 5. Together with the potential in (41), we keep the same lines of the previous example and suggest the term in the Lagrangian density to be modified by where is a real parameter. Substituting and in (22a) and (22b), we obtain To avoid the presence of negative coefficients in the above expressions, we impose that .

The energy density is calculated from (14), which leads toOne may integrate it to get the total energy , which matches with the value obtained by (26). Again, only the parameter modifies the energy of the vortices, meaning that the term in the Lagrangian density (1) plays a significant role in the model. Following a similar procedure that was done in [27], one can show that the energy density tends to compactify into a ringlike region of unit radius in the plane, described byIn Figure 6, we display the energy density for several values of , including the compact limit given above. Its behavior, even with the presence of the parameter , is qualitatively similar to the one found in [35] for the compactification of vortices in a generalized Chern-Simons-Higgs model.

#### 4. Comments and Conclusions

In this work, we have developed a procedure that allows us to construct k-vortex models that support a first order framework. As we discussed above, the method is important because the constraint that dictates the form of the potential cannot be solved in general in the presence of the squared kinetic term of the scalar field, , in the Lagrangian density. Thus, it seems to be very hard to start from a model with this term and find the potential that leads to the first order equations compatible with the stressless condition, vital to the stability of the system.

Nevertheless, we got inspiration from the recent works [26, 27] and noticed that, if an analytical solution is known, we can construct a model that satisfies the stressless condition and find the energy depending exclusively on a function of the fields calculated from the boundary conditions. In order to achieve this, we have introduced the generating function that decouples the first order equations. It is interesting feature of this procedure that it shows there is a class of models that leads to the same analytical, stressless solutions and their respective magnetic fields, which only depend on the generating function. However, the energy density as well as the total energy depends on the model to be chosen, so we have to properly define the model, to make it behave adequately.

It is worth commenting that a similar method can be developed for the more general Lagrangian density . Thus, among the myriad of possibilities, one may develop a construction method for the kinetic term of the scalar field being of the Born-Infeld type, for instance. Other perspectives should include the possibility of considering the case in which the dynamics of the gauge field is driven by the Chern-Simons term, which cannot be multiplied by if one wants to keep gauge invariance. Since the magnetic permeability of the model is generalized, one may also investigate the presence of vortices in metamaterials; see, e.g., [36–38]. Furthermore, as the model supports the in (17), one may seek for supersymmetric extensions, to investigate how the supersymmetry works in this scenario to lead us with first order differential equations. One may also try to extend these results to other topological structures, such as monopoles and skyrmions. We hope to report on some of the above issues in the near future.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

We would like to acknowledge the Brazilian agency CNPq for partial financial support. D. Bazeia appreciates the support from grant 306614/2014-6, L. Losano appreciates the support from grant 303824/2017-4, M. A. Marques appreciates the support from grant 140735/2015-1, and R. Menezes appreciates the support from grant 306826/2015-1.