Analytical BPS Maxwell-Higgs vortices

We have established a prescription for the calculation of analytical vortex solutions in the context of generalized Maxwell-Higgs models whose overall dynamics is controlled by two positive functions of the scalar field. We have also determined a natural constraint between these functions and the Higgs potential allowing the existence of axially symmetric Bogomol'nyi-Prasad-Sommerfield (BPS) solutions possessing finite energy. Furthermore, when the generalizing functions are chosen suitably, the nonstandard BPS equations can be solved exactly. We have studied some examples, comparing them with the usual Abrikosov-Nielsen-Olesen (ANO) solution. The overall conclusion is that the analytical self-dual vortices are well-behaved in all relevant sectors, strongly supporting the generalized models they belong themselves. In particular, our results mimic well-known properties of the usual (numerical) configurations, as localized energy density, while contributing to the understanding of topological solitons and their description by means of analytical methods.


I. INTRODUCTION
In the context of classical field theories, structures possessing topologically nontrivial profiles are usually described as the static solutions of the Euler-Lagrange equations in the presence of finite energy boundary conditions [1]. In some special cases, by requiring the minimization of the corresponding energy functional, such structures can also be described via a set of first-order differential equations also known as BPS equations [2], which provide genuine solutions of the Euler-Lagrange ones.
The kink is a one-dimensional topological object arising within the simplest field model containing only a single real scalar field [3]. Regarding higher dimensional scenarios, the vortex stands for a planar configuration solving some radially symmetric Abelian-Higgs models [4], whilst the magnetic monopole is a threedimensional spherically symmetric object appearing in the non-Abelian-Higgs case [5]. All these solutions possess the minimum energy possible, being stable against decaying into their respective mesons. Moreover, it is well-known that, in order to give rise to topological fields, the corresponding model must allow for the spontaneous symmetry breaking mechanism, with its potential term presenting at least two asymmetric vacua, since topological defects are known to be formed during symmetry breaking phase transitions.
During the last years, a new kind of topologically nontrivial objects have been intensively studied in connection with field models endowed with noncanonical kinetic terms which change the dynamics of the overall system in a nonusual way. It is worthwhile to point out that the motivation regarding exotic dynamics arises in a rather natural way in the context of the string theories. In particular, given some special constraints, field models possessing nonusual dynamics also support minimum energy solutions; see, for instance, Ref. [6]. Moreover, except for their nontrivial nature, these solutions behave in the same general way as their standard counterparts do. On the other hand, exotic kinetic terms also induce slightly variations on the shape of the corresponding field solutions, changing their amplitudes and/or characteristic lengths. Nonstandard field models are defined by introducing generalizing functions on usual field models. Detailed investigations regarding topological defects in the context of these models are found in Ref. [7]. Many authors have also studied interesting applications of these new solutions within several different scenarios, specially involving the accelerated inflationary phase of the universe [8] via the so-called k-essence models [9], strong gravitational waves [10], tachyon matter [11], dark matter [12], and others [13].
Besides the variations on the defect amplitudes and characteristic lengths, these generalizing functions also provide new features for some models, as for example, self-dual analytical solutions which certainly enriches our understanding about integrable systems. Recently, selfdual analytical monopoles were achieved in Ref. [14] in the context of some generalized Yang-Mill models [15]. These new analytical solutions, unattainable in the absence of the modifying functions, were divided into two different classes according to their capability of recovering (or not) the standard 't Hooft-Polyakov result. By following the purpose of achieving analytical solutions for topological defects, the present manuscript aims at investigating the existence of analytical BPS vortex solutions within the framework of the generalized Maxwell-Higgs model [16].
The letter is organized as follows. In Sec. II, we review some important details regarding the generalized Maxwell-Higgs scenario. Sec. III is devoted to describe the prescription implemented to find the analytical BPS solutions of the generalized model. The consistence of our approach is verified by investigating some explicit examples. In the sequel, the new solutions are compared ANO profiles, allowing the identification of the main properties acquired by them. Finally, in Sec. IV, we present our final remarks and conclusions.

II. THE NONSTANDARD MODEL
We begin by reviewing the (1+2)-dimensional generalized Maxwell-Higgs model introduced in Ref. [16], whose Lagrangian density is where J µ = iew φD µ φ − φD µ φ is the generalized current vector, which is also conserved (∂ µ J µ = 0). The stationary Gauss law then reads being trivially verified by A 0 = 0, revealing that the static configurations of the generalized model (1) generate no electric field. The stationary Ampère's law can be written as (already using A 0 = 0) whilst the equation controlling the Higgs field is Here, B = ǫ jk ∂ j A k represents the magnetic field. In order to obtain the first order self-dual equations of the model (1), we start from the expression for the generalized total energy, i.e., which can also be written in the form The energy is minimized by imposing which are the generalized self-dual BPS equations. Considering Eqs. (8), the BPS energy is then reduced to and the static Ampère's law is rewritten as With it, the BPS energy becomes (11) The point to be clarified here is that the integrand in (11) can be reduced to a total derivative only when considering axially symmetric configurations. In this context, Eq. (10) stands for the key condition for attaining self-duality in this generalized Maxwell-Higgs theory.
Hence, from now on, the fields are supposed to be described by the usual axially symmetric vortex Ansatz where n = ±1, ±2, ±3... stands for the vorticity of the resulting configuration, and the magnetic field is The profile functions g (r) and a (r) are constrained to behave according to the standard boundary conditions g (0) = 0 and g (∞) = 1, (14) a (0) = n and a (∞) = 0, giving rise to regular solutions possessing finite energy, as desired. Now, we come back to Eq. (11) defining it in terms of the energy density ε bps related to the BPS solutions as with the auxiliary function H (r) being given by This function is finite at origin, H (0) = H 0 , and fulfills H (∞) = 0. Observing these boundary conditions for H (r) , the resulting total energy (16) is One also remarks that H 0 is proportional to n, the winding number characterizing the vortex solution.
In terms of g (r) and a (r), the BPS equations read which solve the Euler-Lagrange equations of motion. In order to perform such verification explicitly, we first write the Ampère's law (4) in its axially symmetric form which becomes when Eqs. (20) and (21) are used. The form (23) recovers the very same condition (10) that assures the self-duality of the overall model, revealing the consistence of the selfdual equations with the Ampère's law. In addition, one can express Eq. (5) for the Higgs field in terms of g (r) and a (r), that is, It simply provides when saturated by the self-dual equations, which coincides with Eq. (21). In this way, we have explicitly shown that the self-dual equations solve the stationary Euler-Lagrange equations of motion. Due to the arbitrariness of f (g) and w (g), the search for solutions to the axially symmetric Euler-Lagrange equations can be a quite hard task, even in the presence of the suitable boundary conditions (14) and (15). A way to circumvent this point is focusing the attention on the self-dual equations (20) and (21). However, it is worthwhile to reinforce that such equations only hold when the model is constrained by condition (10), also expressed as We can summarize in the following way: given a set of functions f , w and U satisfying (26), regular solutions g (r) and a (r) can be found by solving (20) and (21) using the boundary conditions (14) and (15). The resulting configurations stand for topological vortices possessing finite energy given by (19), which remains proportional to the magnetic flux Φ B = 2πn/e. The proportionality constant is finite and related to the value of f √ 2U near the origin. Moreover, it is worthwhile to point out that is the BPS energy density (17), which becomes positive whenever Eq. (26) ensures a positive w (for a given pair of functions U and f conveniently chosen). The next Section introduces some effective Maxwell-Higgs models for which the BPS equations (20) and (21) can be solved analytically (instead of numerically, as usually done). The analytical profiles representing g (r), a (r), B (r) and ε bps (17) are depicted and compared with the usual (numerical) ANO solution. Furthermore, the main features of the new analytical vortices are highlighted.

III. ANALYTICAL BPS VORTICES
Now, we present the main goal of this work by introducing generalized Maxwell-Higgs models for which the BPS equations (20) and (21) can be solved analytically according the finite energy boundary conditions (14) and (15). Here, for simplicity, we only consider those configurations for which the winding number is equal to the unity (n = 1), although it is also possible to find solutions with higher vorticity, as it will be explained below.
Along this section, we work with the upper signs in Eqs. (20), (21) and (17) only. Also, for simplicity, we set e = v = 1. Our prescription to find analytical self-dual vortices can be described as follows. Firstly, we choose the potential U (g) supporting the spontaneous symmetry breaking of the U (1) local gauge symmetry inherent to the model (1). In the sequel, we choose an analytical function g (r) satisfying the boundary conditions (14). Then, we use g (r) to solve Eq. (20), which allows to obtain the corresponding profile for a (r) fulfilling the boundary conditions (15). In the end, we use Eq. (21) to evaluate the expression for the generalizing function f , writing it as a function of the radial variable r (i.e., regarding the analytical models as effective ones).
A general observation about the kinetic functions f 2 (|φ|) and w (|φ|) is that they are presented as functions of the radial variable r, not of the field variable g, i.e., |φ|. Expressing f and w in terms of g gives very long expressions, when possible. In general, it becomes a very difficult task.
The analytical profiles here obtained provide a set of self-dual vortices possessing finite total energy given by Eq. (19). It is also worthwhile to remember that the corresponding f 's and w's are positive, as required. These new solutions are shown in figs. 1, 2, 3 and 4, from which we highlight their main features.

A. |φ| 4 -models
We first investigate some models defined by the usual fourth-order Higgs potential where the constant the scalar-matter self-interaction was supposed equal to the unity, for simplicity. The |φ| 4 -models here presented possess a generalizing function f (r) finite at the boundaries (i.e., for r = 0 and asymptotically). In this case, the BPS equations (20) and (21) reduce to In particular, one clearly sees that f = 1 leads us back to the standard case dg dr = ag r , which yields the well-known Abrikosov-Nielsen-Olesen numerical vortices [4]. Moreover, families containing nontrivial numerical solutions of the same kind were also studied in Ref. [16]. We now proceed looking for analytical solutions of Eqs. (29) and (30). In this sense, by following our prescription, the first model we introduce is defined by the BPS solution g (r) = tanh (r) , which trivially obeys (14). Replacing it in Eq. (29), one achieves the gauge field profile a (r) = 2r sinh (2r) , whose profile is a lump centered at the origin; see the dash-dotted red line in Fig. 3. Therefore, the solutions (33) and (34) represent analytical BPS Maxwell-Higgs vortices possessing total energy equal to E bps = 3π/2, according to the Eq. (19). The second model, inherent to the |φ| 4 -potential, is defined by the BPS solution which also fulfills the conditions (14). From Eq. (29), we attain the gauge field profile a (r) = r 2 e r 2 − 1 , satisfying the boundary conditions (15) (also with n = 1). Then, by combining Eqs. (38) and (39) in Eq. (30), the following generalizing function is achieved: This is a positive and finite function whose values at the boundaries are f (r = 0) = 1 and f (r = ∞) = 0. Moreover, the corresponding H (r) is given by providing H (0) = −1 and H (∞) = 0, as required. The magnetic field, also is a lump centered at the origin; see the dashed green line in Fig. 3. This way, the profiles Eqs. (38) and (39) describe analytical self-dual vortices whose BPS total energy is E bps = 2π. In particular, it means that these analytical solutions saturate the very same Bogomol'nyi bound fulfilled by the usual (numerical) n = 1 ANO vortex.
B. |φ| 6 -models As it was shown in Ref. [16], the generalized model allows to find BPS vortices even in the presence of a higher order potential describing the scalar-matter selfinteraction. In this sense, we now go further by introducing analytical self-dual vortices arising in the presence of a sixth-order potential, defined by The vacuum manifold of the corresponding |φ| 6 -model is represented by a dot surrounded by a circle, the dot standing for a symmetric vacuum. As a consequence, under a suitable choose of the boundary conditions to be satisfied by the profile functions g (r) and a (r), the model also supports nontopological self-dual structures possessing finite energy. Indeed, some of us have already obtained such objects, with results being forthcoming reported [17]. It is worthwhile to point out that nontopological vortices do not occur in the |φ| 4 -model (28) because it has no a symmetric vacuum (i.e., its vacuum manifold is a circle). On the other hand, it is well-known that the |φ| 6potential (43) ensures the self-duality of the usual Chern-Simons-Higgs (CSH) model, whose topological vortices possess both electric and magnetic fields. Despite our generalized |φ| 6 -Maxwell-Higgs model supports only noncharged self-dual solutions, they are expected to behave in the same general way the CSH ones do. Some numerical self-dual |φ| 6 -Maxwell-Higgs models were already investigated in Ref. [16]. Here, for completeness, we consider only the analytical solutions for it.
Returning to our prescription, under the |φ| 6 -potential (43), the BPS equations (20) and (21) can be written as In addition, as we demonstrate below, the resulting generalizing function f (r) is not necessarily finite at the boundaries. The first analytical model has the profile g (r) defined by from which one gets, according to Eq. (44), the corresponding gauge solution, whilst Eq. (45) gives the following generalizing function It has the following behavior at the boundaries: f (0) = f (∞) = ∞. However, even in this case, one still achieves well-behaved solutions possessing finite total energy. To clarify the way it happens, we calculate the fields B (r) and H (r) arising from Eqs. (47) and (46). The magnetic field, presents a ringlike profile, a typical magnetic behavior related to the |φ| 6 -vortices; see the dotted blue line in Fig. 3. The auxiliary function H (r) reads as which, at the boundaries, assumes the values H (0) = −1/4 and H (∞) = 0, as desired. As a result, the analytical profiles (46) and (47) give rise to a BPS vortex whose total energy is E bps = π/2. The second |φ| 6 -model is a little bit more sophisticated than the previous ones, the Higgs profile being given by with the function Ei (1, r) standing for the exponential integral Nevertheless, the gauge field has a simpler solution, whereas the corresponding functions f (r) and H (r) are defined by is divergent at the origin, whilst vanishing asymptotically. Notwithstanding, as in the previous case, the Bogomol'nyi bound for the BPS total energy saturates at E bps = π √ e γ . The magnetic field of the resulting selfdual configuration, also presents the aforecited typical ringlike behavior; see the long-dashed gold line in Fig. 3.
In the sequel, we depict all the analytical solutions together with the standard ANO profile, from which we highlight their main features and also the differences of the generalized solutions in comparison with the usual Maxwell-Higgs ones.
The analytical solutions defining the Higgs profiles g (r) are depicted in Fig. 1. The profile (33) is plotted with the dash-dotted red line, whilst Eq. (38) is represented by the dashed green line; these solutions corresponding to the generalized |φ| 4 -models. On the other hand, the dotted blue line stands for (46) and the long-dashed gold line represents (51); both belonging to the noncanonical |φ| 6 -models. The standard (numerical) n = 1 ANO Higgs profile is drawn with the solid black line. The overall conclusion is that the analytical solutions behave in the same general way the standard one does. However, the new profiles saturate the asymptotic value g (r = ∞) = 1 faster, so that the new Higgs profiles are more localized, and the corresponding bosons are more massive than the ANO ones. From now on, we follow the same line/color definitions established in Fig.  1. The gauge profiles a (r) are plotted in Fig. 2. There, we see that the profiles (46) and (51), related to the |φ| 6 -potential, have developed a plateau close to the origin, such structure being a common feature presented in the self-dual |φ| 6 -scenarios. On the other hand, the profiles (34) and (39), corresponding to the |φ| 4 -models, are lumps centered at the origin. The amplitude in r = 0 corresponds to the winding number n = 1, as already commented. For large radius, all profiles vanish monotonically.
We show the profiles we have found for the magnetic field B (r) in Fig. 3. The solutions regarding the fourthorder potential are lumps centered at the origin, just as the ANO magnetic field. On the other hand, the solutions related to the sixth-order potential present a ringlike behavior, as expected when considering the magnetic fields belonging to the self-dual |φ| 6 -vortices. We also observe differences both on the amplitudes and on the characteristic lengths of the magnetic profiles.
By last, the BPS energy densities are depicted in Fig.  4. We see that all the profiles are lumps centered at origin, as expected for n = 1 vortices (including the |φ| 6ones). Besides the different characteristic lengths, it is also interesting to note that the profiles coming from the fourth-order potential have achieved greater amplitudes than those ones related to the sixth-order potential.
In the next Section, we present our final considerations and the perspectives regarding future works.

IV. ENDING COMMENTS
We have investigated the existence of analytical BPS vortices within the nonstandard Maxwell-Higgs scenario proposed in Ref. [16], such model being generalized by two positive functions, namely f (|φ|) and w (|φ|), which change the overall dynamics of the original theory; see Eq.(1). By imposing the radially symmetric ansatz (12) for the gauge and scalar fields, we have reviewed the BPS framework and achieved the first-order equations whose corresponding solutions have finite energy. The self-duality arises when the condition (26) is satisfied by the potential U and the functions f and w.
The existence of a well-defined Bogomol'nyi bound (19) is ruled by an auxiliary function H (r) obeying appropriated boundary conditions which guarantee a finite total energy. So, the well-behaved H (r) selects the function f (r) defining the generalized model. With this prescription, we have obtained analytical profiles to some n = 1 self-dual vortex configurations within both |φ| 4and |φ| 6 -models saturating different energy bounds. In order to compare the analytical self-dual solutions with the ANO ones, we have depicted the profiles for the scalar, gauge and magnetic fields, in figs. 1, 2 and 3, respectively, whilst the BPS energy density is plotted in Fig. 4. The overall conclusion is that the new profiles are well-behaved in all relevant sectors, assuring the consistence of the generalized models here proposed. Furthermore, all the solutions have provided localized energy densities and magnetic fields, as expected. In general, they mimic the behavior of the well-known numerical configurations, as the ringlike magnetic field related to the self-dual |φ| 6 -vortices. Concerning the possibility of obtaining analytical solutions with higher winding numbers, no obvious route for it seems to be available within the models we have studied in this work. Nevertheless, it is possible to construct analytical BPS vortices possessing higher vorticity but different generalizing functions f (r). Such procedure is clarified by the following example related to the |φ| 6 -models. We propose g (r) = r n (1 + r m ) n/m , with m, n > 0, as a generalization for the Higgs profile given by Eq. (46). This leads to a(r) = n 1 + r m , compatible with the boundary conditions for a vortex possessing any integer winding number (i.e., n = +1, +2, +3...). The resulting magnetic field is whilst the corresponding f (r) and H (r) are f (r) = r n−m+2 nm R 2−3n/m R 2n/m − r 2n 2 , H (r) = r 2n−m+2 m R 1−6n/m R 2n/m − r 2n 2 , respectively, with R = 1 + r m . However, only the case m = 2n + 2 provides positive f and w, nonsingular magnetic field and finite energy. Obviously, the choice n = 1 and m = 4 reduces this general proposal to the case (46). The generalization of the vortex configurations presenting higher winding numbers within other field scenarios can follow this general idea. Regarding future works, interesting issues include the search for the nontopological self-dual vortices arising in the generalized Maxwell-Higgs scenario (1) when endowed by the sixth-order potential (43); see Ref. [17].
In parallel, some of us are also working in a general formulation of the deformation method [18] applicable to field models possessing generalized dynamics [19]. These two fronts are now under investigation, and we expect interesting results for a future report.
The Brazilian authors thank CAPES, CNPq and FAPEMA for partial financial support.