#### Abstract

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, namely, and . 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 consistency of the respective generalized models. 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.

#### 1. 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, 3], 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 [4]. Regarding higher-dimensional scenarios, the vortex stands for a planar configuration solving some radially symmetric Abelian-Higgs models [5, 6], whilst the magnetic monopole is a three-dimensional spherically symmetric object appearing in the non-Abelian-Higgs case [7, 8]. 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 has 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, [9–17]. 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 [18–35]. Many authors have also studied interesting applications of these new solutions within several different scenarios, specially involving the accelerated inflationary phase of the universe [36] via the so-called -essence models [37], strong gravitational waves [38], tachyon matter [39], dark matter [40], and others [41–43].

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, self-dual analytical monopoles were achieved in [44] in the context of some generalized Yang-Mill models [45]. 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 paper aims at investigating the existence of analytical BPS vortex solutions within the framework of the generalized Maxwell-Higgs model [46].

The letter is organized as follows. In Section 2, we review some important details regarding the generalized Maxwell-Higgs scenario. Section 3 is devoted to describing 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 to the ANO profiles, allowing the identification of the main properties acquired by them. Finally, in Section 4, we present our final remarks and conclusions.

#### 2. The Nonstandard Model

We begin by reviewing the (1 + 2)-dimensional generalized Maxwell-Higgs model introduced in [46], whose Lagrangian density is where is the usual field strength tensor and stands for the covariant derivative. Moreover, and are positive functions which change the dynamics of the overall model, being called dielectric functions because they mimic some effective electrodynamics in continuous media, as already mentioned in the literature. So far, the possibility of such generalizing functions to provide exactly solvable models for vortex configurations was not explored, being the main motivation of this paper. Here, for simplicity, all fields, coordinates, and parameters are supposed to be dimensionless.

The corresponding Euler-Lagrange equation for the gauge field is where is the generalized current vector, which is also conserved (). The stationary Gauss law then reads being trivially verified by , revealing that the static configurations of the generalized model (1) generate no electric field.

Stationary Ampère's law can be written as (already using ) whilst the equation controlling the Higgs field is Here, 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; that is, which can also be written in the form The energy is minimized by imposing which are the generalized self-dual BPS equations. Considering (8), the BPS energy is then reduced to and static Ampère's law is rewritten as With it, the BPS energy becomes 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, (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 stands for the vorticity of the resulting configuration, and the magnetic field is The profile functions and are constrained to behave according to the standard boundary conditions giving rise to regular solutions possessing finite energy, as desired.

Now, we come back to (11) defining it in terms of the energy density related to the BPS solutions as where with the auxiliary function being given by This function is finite at origin, , and fulfills . Observing these boundary conditions for , the resulting total energy (16) is One also remarks that is proportional to , the winding number characterizing the vortex solution.

In terms of and , the BPS equations read which solve the Euler-Lagrange equations of motion. In order to perform such verification explicitly, we first write Ampère's law (4) in its axially symmetric form which becomes when (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 self-dual equations with Ampère's law. In addition, one can express (5) for the Higgs field in terms of and ; that is, It simply provides when saturated by the self-dual equations, which coincides with (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 and , 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 , , and satisfying (26), regular solutions and 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 . The proportionality constant is finite and related to the value of near the origin. Moreover, it is worthwhile to point out that is the BPS energy density (17), which becomes positive whenever (26) ensures a positive (for a given pair of functions and 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 , , , and (17) are depicted and compared with the usual (numerical) ANO solution. Furthermore, the main features of the new analytical vortices are highlighted.

#### 3. 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 , 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 (20), (21), and (17) only. Also, for simplicity, we set . Our prescription to find analytical self-dual vortices can be described as follows. Firstly, we choose the potential supporting the spontaneous symmetry breaking of the local gauge symmetry inherent to model (1). In the sequel, we choose an analytical function satisfying the boundary conditions (14). Then, we use to solve (20), which allows obtaining the corresponding profile for fulfilling the boundary conditions (15). In the end, we use (21) to evaluate the expression for the generalizing function , writing it as a function of the radial variable (i.e., regarding the analytical models as effective ones).

A general observation about the kinetic functions and is that they are presented as functions of the radial variable , not of the field variable , that is, . Expressing and in terms of 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 (19). It is also worthwhile to remember that the corresponding 's and 's are positive, as required. These new solutions are shown in Figures 1, 2, 3, and 4, from which we highlight their main features.

##### 3.1. -Models

We first investigate some models defined by the usual fourth-order Higgs potential where the coupling constant for the scalar-matter self-interaction was supposed to be equal to the unity, for simplicity.

The -models here presented possess a generalizing function finite at the boundaries (i.e., for and asymptotically). In this case, the BPS equations (20) and (21) reduce to In particular, one clearly sees that leads us back to the standard case which yields the well-known Abrikosov-Nielsen-Olesen numerical vortices [5, 6]. Moreover, families containing nontrivial numerical solutions of the same kind were also studied in [46].

We now proceed looking for analytical solutions of (29) and (30). In this sense, by following our prescription, the first model we introduce is defined by the BPS solution which trivially obeys (14). Replacing it in (29), one achieves the gauge field profile satisfying the boundary conditions (15), with . Now, by using (32) and (33) in (30), one gets the fact that the corresponding function reads as being a smooth and positive function with values and . Finally, we obtain the auxiliary function by substituting (28), (33), and (34) into (18). The resulting expression is a smooth function for which the values at the boundaries read as and , as desired. Finally, the magnetic field associated with (33) is given by whose profile is a lump centered at the origin; see the dash-dotted red line in Figure 3. Therefore, the solutions (32) and (33) represent analytical BPS Maxwell-Higgs vortices possessing total energy equal to , according to (19).

The second model, inherent to the -potential, is defined by the BPS solution which also fulfills the conditions (14). From (29), we attain the gauge field profile satisfying the boundary conditions (15) (also with ). Then, by combining (37) and (38) in (30), the following generalizing function is achieved: This is a positive and finite function whose values at the boundaries are and . Moreover, the corresponding is given by providing and , as required. The magnetic field, also is a lump centered at the origin; see the dashed green line in Figure 3. This way, the profiles equations (37) and (38) describe analytical self-dual vortices whose BPS total energy is . In particular, it means that these analytical solutions saturate the very same Bogomol'nyi bound fulfilled by the usual (numerical) ANO vortex.

##### 3.2. -Models

As was shown in [46], the generalized model allows finding BPS vortices even in the presence of a higher-order potential describing the scalar-matter self-interaction. 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 -model is represented by a dot surrounded by a circle, with the dot standing for a symmetric vacuum. As a consequence, under a suitable choice of the boundary conditions to be satisfied by the profile functions and , 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 [47]. It is worthwhile to point out that nontopological vortices do not occur in the -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 -potential (42) ensures the self-duality of the usual Chern-Simons-Higgs (CSH) model, whose topological vortices possess both electric and magnetic fields. Despite the fact that our generalized -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 -Maxwell-Higgs models were already investigated in [46]. Here, for completeness, we consider only the analytical solutions for it.

Returning to our prescription, under the -potential (42), the BPS equations (20) and (21) can be written as In addition, as we demonstrate below, the resulting generalizing function is not necessarily finite at the boundaries.

The first analytical model has the profile defined by from which one gets, according to (43), the corresponding gauge solution whilst (44) gives the following generalizing function: It has the following behavior at the boundaries: . 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 and arising from (46) and (45). The magnetic field presents a ring-like profile, a typical magnetic behavior related to the -vortices; see the dotted blue line in Figure 3. The auxiliary function reads as which, at the boundaries, assumes the values and , as desired. As a result, the analytical profiles (45) and (46) give rise to a BPS vortex whose total energy is .

The second -model is a little bit more sophisticated than the previous ones, with the Higgs profile being given by with the function standing for the exponential integral Nevertheless, the gauge field has a simpler solution, whereas the corresponding functions and are defined by from which we get the fact that and , whilst and , being Euler's constant (). We see that 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 . The magnetic field of the resulting self-dual configuration, also presents the aforecited typical ring-like behavior; see the long-dashed gold line in Figure 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 are depicted in Figure 1. The profile (32) is plotted with the dash-dotted red line, whilst (37) is represented by the dashed green line; these solutions correspond to the generalized -models. On the other hand, the dotted blue line stands for (45) and the long-dashed gold line represents (50), both belonging to the noncanonical -models. The standard (numerical) 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 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 Figure 1.

The gauge profiles are plotted in Figure 2. There, we see that the profiles (45) and (50), related to the -potential, have developed a plateau close to the origin, with such structure being a common feature presented in the self-dual -scenarios. On the other hand, the profiles (33) and (38), corresponding to the -models, are lumps centered at the origin. The amplitude in corresponds to the winding number , as already commented. For large radius, all profiles vanish monotonically.

We show the profiles we have found for the magnetic field in Figure 3. The solutions regarding the fourth-order 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 ring-like behavior, as expected when considering the magnetic fields belonging to the self-dual -vortices. We also observe differences both on the amplitudes and on the characteristic lengths of the magnetic profiles.

At last, the BPS energy densities are depicted in Figure 4. We see that all the profiles are lumps centered at origin, as expected for vortices (including the -ones). 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.

#### 4. Ending Comments

We have investigated the existence of analytical BPS vortices within the nonstandard Maxwell-Higgs scenario proposed in [46], with such model being generalized by two positive functions, namely, and , which change the overall dynamics of the original theory; see (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 and the functions and .

The existence of a well-defined Bogomol'nyi bound (19) is ruled by an auxiliary function obeying appropriated boundary conditions which guarantee a finite total energy. So, the well-behaved selects the function defining the generalized model. With this prescription, we have obtained analytical profiles to some self-dual vortex configurations within both - and -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 Figures 1, 2, and 3, respectively, whilst the BPS energy density is plotted in Figure 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 ring-like magnetic field related to the self-dual -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 . Such procedure is clarified by the following example related to the -models. We propose with , , as a generalization for the Higgs profile given by (45). This leads to compatible with the boundary conditions for a vortex possessing any integer winding number (i.e., ). The resulting magnetic field is whilst the corresponding and are respectively, with . However, only the case provides positive and , nonsingular magnetic field and finite energy. Obviously, the choice or reduces this general proposal to the case (45). 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 (42); see [47]. In parallel, some of us are also working in a general formulation of the deformation method [48–50] applicable to field models possessing generalized dynamics [51]. These two fronts are now under investigation, and we expect interesting results for a future report.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The Brazilian authors thank CAPES, CNPq, and FAPEMA for partial financial support.