text:
# Self-Dual Cosmic Strings and Gravitating Vortices in Gauged Sigma Models

## 1 Introduction

Topological defects are generally believed to have been formed during a series of phase transitions in the early universe. In particular, many field-theoretical models suggest the formation of cosmic strings (for a review, see for instance Vilenkin and Shellard ).

The BOOMERANG results , especially concerning the location of the acoustic peaks , essentially rule out GUT scale ($`10^{16}`$GeV) topological defects being the sole source of the cosmic microwave background (CMB) anisotropy, which reflects the early universe density fluctuations, eventually leading to galaxy formation. However, recent analysis  shows that a mixture of cosmic strings and inflation is consistent with current CMB data.

Independently of a possible role in the large scale structure formation, and even if their mass scale is lower than the GUT scale, cosmic strings may have also a number of significant observable astrophysical effects. Just to mention a few: Cosmic strings could be sources of double images , of gravitational waves  and of ultra high energy cosmic rays . Recently it was even suggested  that they can serve as gamma ray burst engines. Thus, cosmic strings are still of considerable interest in cosmology and astrophysics.

The most popular field-theoretical system, used to model cosmic string generation , is the Abelian Higgs model. There exist, however, other systems which may be used for the same purpose and are not more complicated. As long as we model cosmic strings as static cylindrically symmetric sources coupled to gravity, which is essentially a two dimensional system, we may borrow circular solutions of other well-known two dimensional systems, namely non-linear sigma models.

Non-linear sigma models  are very much used as effective theories describing various systems such as low energy effective QCD, and are also used as ”toy models” for gauge theories. A generalization of these, where the global symmetry group or its subgroup is promoted to be local, are gauged sigma models which appear naturally in supersymmetric field theories.

All these models have been extensively studied for decades. Most of the results were obtained in flat (Minkowski) spacetime, but more recently gravitating solutions were studied too. However, the gravitating sigma model solutions are usually taken to be spherically symmetric as in the cases of gravitating Skyrmions and textures (see e.g. Volkov and Gal’tsov ), although some authors gave attention to cylindrically symmetric solutions of the non-gauged  as well as the gauged models with either Maxwell  or Chern-Simons terms .

Here we will analyze the gravitating gauged sigma models with a general two-dimensional target space defined in terms of two field dependent ”dielectric functions”, which may be viewed as generalized Abelian Higgs models: Maxwell type and Chern-Simons type, according to the gauge field term in the Lagrangian. We will find the general conditions for self-dual cosmic string solutions and get several kinds of interesting solutions.

We will start by considering, in the next section, static and translationally invariant solutions of a generic gauged gravitating non-linear sigma model of a Maxwell type, and find simple conditions for the existence of self-dual solutions. In sections 3 and 4, we concentrate in the O(3) model and add rotational symmetry to obtain self-dual cosmic string solutions. In section 5, we move to the Chern-Simons type of model, and consider stationary (not static) solutions. First we obtain the self-duality conditions for this kind of solutions, and then we specialize to spinning cosmic vortices. In section 6, we discuss the relation between the two types of theories and solutions, and we end with a discussion of the topological charges in section 7.

## 2 The Generalized Abelian Higgs Model

We start by considering a generalization of the Abelian Higgs model defined by the action:

$`S={\displaystyle \mathrm{d}^4x\sqrt{|g|}\left(\frac{1}{2}_1(|\mathrm{\Phi }|)(D_\mu \mathrm{\Phi })^{}(D^\mu \mathrm{\Phi })U(|\mathrm{\Phi }|)\frac{1}{4}_2(|\mathrm{\Phi }|)F_{\mu \nu }F^{\mu \nu }+\frac{1}{16\pi G}R\right)}`$ (2.1)
where $`_1(|\mathrm{\Phi }|)`$ and $`_2(|\mathrm{\Phi }|)`$ are non-negative dimensionless functions, which may be interpreted as Weyl factors of a conformally-flat target space metric of a non-linear sigma model with local U(1) symmetry . The function $`_2(|\mathrm{\Phi }|)`$ plays further the role of a dielectric function . In what follows, we use the collective notation $`_a,(a=1,2`$), and refer to both as ”dielectric functions”. The action (2.1) also generalizes the ones considered by Lohe . Furthermore, the action (2.1) arises in various extended supergravity theories . In the general case, we may also add the term $`_3(|\mathrm{\Phi }|)ϵ^{\mu \nu \rho \sigma }F_{\mu \nu }F_{\rho \sigma }`$, but this term vanishes identically for the configurations considered in the present paper and has no effect on the field equations.

The field equations derived from the action (2.1) are:

$$_1(|\mathrm{\Phi }|)D_\mu D^\mu \mathrm{\Phi }+\frac{\mathrm{\Phi }^{}}{2|\mathrm{\Phi }|}\frac{d_1}{d|\mathrm{\Phi }|}D_\mu \mathrm{\Phi }D^\mu \mathrm{\Phi }+\frac{\mathrm{\Phi }}{|\mathrm{\Phi }|}\frac{dU}{d|\mathrm{\Phi }|}+\frac{\mathrm{\Phi }}{4|\mathrm{\Phi }|}\frac{d_2}{d|\mathrm{\Phi }|}F_{\mu \nu }F^{\mu \nu }=0$$
(2.2)
$$_\mu \left(_2(|\mathrm{\Phi }|)F^{\mu \nu }\right)=j^\nu =\frac{i}{2}e_1(|\mathrm{\Phi }|)(\mathrm{\Phi }^{}(D^\nu \mathrm{\Phi })\mathrm{\Phi }(D^\nu \mathrm{\Phi })^{})$$
(2.3)
$`{\displaystyle \frac{1}{8\pi G}}R_{\mu \nu }+{\displaystyle \frac{1}{2}}_1(|\mathrm{\Phi }|)\left((D_\mu \mathrm{\Phi })^{}(D_\nu \mathrm{\Phi })+(D_\nu \mathrm{\Phi })^{}(D_\mu \mathrm{\Phi })\right)U(|\mathrm{\Phi }|)g_{\mu \nu }+`$
$`_2(|\mathrm{\Phi }|)\left(g^{\kappa \lambda }F_{\mu \kappa }F_{\lambda \nu }+{\displaystyle \frac{1}{4}}F_{\kappa \lambda }F^{\kappa \lambda }g_{\mu \nu }\right)=0`$ (2.4)
Conventions: $`D_\mu =_\mu ieA_\mu `$, signature $`(+,,,)`$ and $`R_{\lambda \mu \nu }^\kappa =_\nu \mathrm{\Gamma }_{\lambda \mu }^\kappa _\mu \mathrm{\Gamma }_{\lambda \nu }^\kappa +..`$.

Usually (but not always) we will be interested in potentials, which ensure spontaneous symmetry breaking and which lead to a massive gauge field. That is to say, potentials with a circle of degenerate minima at (say) $`|\mathrm{\Phi }|=v`$ and with positive second derivative at $`v`$. We further normalize the potential such that it vanishes in the vacuum. Thus, altogether:

$$U=0,U^{}=\frac{dU}{d|\mathrm{\Phi }|}=0;\text{for}|\mathrm{\Phi }|=v$$
(2.5)
With these normalizations, the Higgs mass and gauge boson mass are, respectively:

$$m_H^2=\frac{U^{\prime \prime }(v)}{_1(v)},m_A^2=\frac{e^2v^2_1(v)}{_2(v)}$$
(2.6)
In order to study cosmic string solutions, we assume the metric and matter fields to be static and symmetric under $`z`$-translations. Thus we assume that $`A_\mu `$ and $`\mathrm{\Phi }`$ depend only on the two transverse coordinates $`x^k`$, and for the metric we use the following form:

$$ds^2=N^2(x{}_{}{}^{k})dt^2\gamma _{ij}(x{}_{}{}^{k})dx{}_{}{}^{i}dx{}_{}{}^{j}K^2(x{}_{}{}^{k})dz^2$$
(2.7)
We also require the presence of a magnetic field only, i.e. that the gauge potential will have the form:

$$A_\mu dx^\mu =A_i(x{}_{}{}^{k})dx^i$$
(2.8)
such that the Maxwell tensor will contain a single magnetic component $`B`$:

$$F_{\mu \nu }dx^\mu dx^\nu =B\sqrt{|\gamma |}ϵ_{ij}dx^idx^j$$
(2.9)
where $`|\gamma |=|\text{det}(\gamma _{ij})|`$. In order to get the field equations for static solutions, we compute also the components of the Ricci tensor:

$`R_{00}`$ $`=`$ $`{\displaystyle \frac{N}{K}}_i(K^iN)`$
$`R_{33}`$ $`=`$ $`{\displaystyle \frac{K}{N}}_i(N^iK)`$
$`R_{ij}`$ $`=`$ $`R_{ij}(\gamma )+{\displaystyle \frac{1}{N}}_i_jN+{\displaystyle \frac{1}{K}}_i_jK`$
$`R_{i0}`$ $`=`$ $`R_{i3}=0,R_{03}=0`$ (2.10)
where $`_i`$ is the covariant derivative with respect to the two-dimensional metric $`\gamma _{ij}`$ and $`R_{ij}(\gamma )`$ the corresponding Ricci tensor.

A significant simplification of this system is obtained if self-duality conditions are satisfied, i.e. if the system admits a Bogomolnyi limit . It is well-known  that in the usual Higgs model, the flat space considerations can be carried over to curved background if $`N(x^i)`$ and $`K(x^i)`$ are constants, say, 1. We will see now that the present generalized Higgs model has also a Bogomolnyi limit if we keep $`N(x^i)=K(x^i)=1`$. If we use these conditions, we find that the $`(00)`$ and $`(33)`$ components of Einstein equations will be satisfied only if:

$$U(|\mathrm{\Phi }|)=\frac{1}{2}_2(|\mathrm{\Phi }|)B^2$$
(2.11)
Now we turn to the $`(ij)`$ components of the Einstein equations, or even better to $`G_{ij}`$ which vanish identically. Consequently, $`T_{ij}=0`$ as well and we get:

$`{\displaystyle \frac{1}{2}}_1(|\mathrm{\Phi }|)\left((D_i\mathrm{\Phi })^{}(D_j\mathrm{\Phi })+(D_j\mathrm{\Phi })^{}(D_i\mathrm{\Phi })\right)`$
$`\left({\displaystyle \frac{1}{2}}_1(|\mathrm{\Phi }|)\gamma ^{kl}(D_k\mathrm{\Phi })^{}(D_l\mathrm{\Phi })+U(|\mathrm{\Phi }|){\displaystyle \frac{1}{2}}_2(|\mathrm{\Phi }|)B^2\right)\gamma _{ij}=0`$ (2.12)
Using (2.11), this condition simplifies further and it follows that it is equivalent to the curved spacetime version of the self-duality condition:

$$D_i\mathrm{\Phi }=i\eta \sqrt{|\gamma |}ϵ_{ij}\gamma ^{jk}D_k\mathrm{\Phi }$$
(2.13)
where $`\eta =\pm 1`$ corresponds to self-dual or anti self-dual solutions. This is a first order equation for the Higgs field. Analogously, Eq.(2.11) is a first order equation for the gauge potential.

We also find the following expression for the two dimensional Ricci scalar:

$$R(\gamma )=8\pi G\left(_1(|\mathrm{\Phi }|)\gamma ^{ij}(D_i\mathrm{\Phi })^{}(D_j\mathrm{\Phi })+4U(|\mathrm{\Phi }|)\right)$$
(2.14)
which serves as the Einstein equation for the two metric $`\gamma _{ij}`$.

Equation (2.2) for the Higgs field becomes a consistency condition, which constrains the form of the potential $`U(|\mathrm{\Phi }|)`$:

$$\eta e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)B+\frac{dU}{d|\mathrm{\Phi }|}+\frac{B^2}{2}\frac{d_2}{d|\mathrm{\Phi }|}=0$$
(2.15)
Maxwell equations (2.3) give:

$$_j(_2(|\mathrm{\Phi }|)B)=\frac{\eta e}{2}_1(|\mathrm{\Phi }|)_j|\mathrm{\Phi }|^2$$
(2.16)
which is a second order equation for the gauge potential, but it is not an independent one, as is easily shown: If we think of $`B`$ as a function of $`|\mathrm{\Phi }|`$, we may use (2.11) and get from (2.15) the following relation, which we will refer to as the ”Bogomolnyi constraint”:

$$\frac{d}{d|\mathrm{\Phi }|}(_2(|\mathrm{\Phi }|)B)+\eta e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)=0$$
(2.17)
This is actually again the Maxwell equation (2.16) in disguise, and can be used together with (2.11) to get the function $`U(|\mathrm{\Phi }|)`$ for any set of given ”dielectric functions” $`_a,a=1,2`$. In order to do it directly, we need to express the magnetic field $`B`$ in terms of the potential $`U`$. If $`B`$ has a definite sign, we infer from the Bogomolnyi constraint that $`_2(|\mathrm{\Phi }|)B(|\mathrm{\Phi }|)0`$ for $`\eta =+1`$, and similarly for negative values, so we may take a square root of (2.11):

$$B=\eta \sqrt{2U/_2}$$
(2.18)
and get a simple relation between the potential and the two dielectric functions, which is a necessary condition for the action (2.1) to have self-dual solutions, i.e. a Bogomolnyi limit:

$$\frac{d}{d|\mathrm{\Phi }|}\sqrt{2_2(|\mathrm{\Phi }|)U(|\mathrm{\Phi }|)}+e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)=0$$
(2.19)
The simplest case where we can apply (2.19) is of course $`_a=1`$, which reproduces immediately the standard Higgs potential:

$$U(|\mathrm{\Phi }|)=\frac{a}{4}(v^2|\mathrm{\Phi }|^2)^2$$
(2.20)
where the Bogomolnyi relation between the coupling constants holds:

$$\alpha e^2/a=2$$
(2.21)
Note that here, the vacuum expectation value $`v`$ enters into the potential as an integration constant, and has also a role of the value of $`|\mathrm{\Phi }|`$ for which the magnetic field vanishes.

The next case is the non-linear O(3) sigma model, which is obtained by taking $`_1`$ to be the usual conformal factor for $`\text{S}^2`$ with a radius $`\mu /2`$ :

$$_1(|\mathrm{\Phi }|)=1/(1+|\mathrm{\Phi }|^2/\mu ^2)^2$$
(2.22)
The parameter $`\mu `$ sets a second energy scale in the system and renders $`_1`$ dimensionless. It is very simple to integrate the Bogomolnyi constraint (2.17), or equivalently Eq.(2.19), also in this case and get the following form of $`B(|\mathrm{\Phi }|)`$ and potential, which allows a Bogomolnyi limit in the O(3) generalized Abelian Higgs model:

$$B(|\mathrm{\Phi }|)=\frac{\eta e}{2(1+v^2/\mu ^2)}\frac{1}{_2(|\mathrm{\Phi }|)}\frac{v^2|\mathrm{\Phi }|^2}{1+|\mathrm{\Phi }|^2/\mu ^2}$$
(2.23)
$$U(|\mathrm{\Phi }|)=\frac{e^2}{8(1+v^2/\mu ^2)^2}\frac{1}{_2(|\mathrm{\Phi }|)}\frac{(v^2|\mathrm{\Phi }|^2)^2}{(1+|\mathrm{\Phi }|^2/\mu ^2)^2}$$
(2.24)
Here again $`v`$ is an integration constant, which parametrizes the minimum of the potential and the field value for which the magnetic field vanishes. Generally there is no relation between this scale and the scale $`\mu `$. We therefore have a one-parameter family of potentials in this O(3) generalized Abelian Higgs model for any given $`_2(|\mathrm{\Phi }|)`$, which we may take at our will. Some special cases with $`_2(|\mathrm{\Phi }|)=1`$ were already analyzed by several authors mainly in flat space , but also with coupling to gravity .

A more ”symmetric” picture may be obtained by using the angular variable on $`\text{S}^2`$ defined by:

$$|\mathrm{\Phi }|=\mu \mathrm{tan}(\mathrm{\Theta }/2)$$
(2.25)
Generically, the potential exhibits a U(1)-symmetry breaking minimum at $`|\mathrm{\Phi }|=v`$, but there are two special limits $`v=0,\mathrm{}`$, where the ground state is only a point in target space (the north or south poles of $`\text{S}^2`$); namely, no symmetry breaking occurs. This pattern of symmetry breaking is reflected by the Higgs and gauge boson masses, which turn out to have the following equal (due to the self-duality) values:

$$m_H=m_A=\frac{ev}{1+v^2/\mu ^2}$$
(2.26)
The two cases of symmetric vacuum correspond to vanishing Higgs and gauge masses.

Just as a final check we note that the Higgs potential ((2.20) with (2.21)) is obtained in the limit $`\mu \mathrm{}`$ from (2.24) with $`_2(|\mathrm{\Phi }|)=1`$. The masses which are obtained in this limit are the usual ones.

Now we return to the Einstein equation (2.14) and notice that for self-dual solutions, its right hand side is actually a two-dimensional divergence (it must be as we will see in section 7), and may be simplified to:

$$R(\gamma )=\frac{16\pi G}{\eta e}_i\left(\gamma ^{ij}_2(|\mathrm{\Phi }|)B_j\mathrm{log}|\mathrm{\Phi }|\right)$$
(2.27)
Without loss of generality, we may take the two-dimensional transverse part of the metric tensor to be conformally flat, i.e. $`\gamma _{ij}=H^2(x^k)\delta _{ij}`$ so the Einstein equation (2.14) reduces further to:

$$\delta ^{ij}_i\left(_j\mathrm{log}H\frac{8\pi G}{\eta e}_2(|\mathrm{\Phi }|)B_j\mathrm{log}|\mathrm{\Phi }|\right)=0$$
(2.28)
This equation is not of first order, but it may get the form of a two dimensional Laplace equation, if we introduce a ”super dielectric function” $`𝒜(|\mathrm{\Phi }|)`$, which solves the following equation:

$$|\mathrm{\Phi }|\frac{d𝒜}{d|\mathrm{\Phi }|}=\eta _2(|\mathrm{\Phi }|)B(|\mathrm{\Phi }|)=\sqrt{2_2(|\mathrm{\Phi }|)U(|\mathrm{\Phi }|)}$$
(2.29)
or equivalently the second order equation:

$$\frac{1}{|\mathrm{\Phi }|}\frac{d}{d|\mathrm{\Phi }|}\left(|\mathrm{\Phi }|\frac{d𝒜}{d|\mathrm{\Phi }|}\right)+e_1(|\mathrm{\Phi }|)=0$$
(2.30)
Now we can use the function $`𝒜`$, in order to give Eq. (2.28) the form of a Laplace equation:

$$\delta ^{ij}_i_j\left(\mathrm{log}H\frac{8\pi G}{e}𝒜\right)=0$$
(2.31)
Special solutions will be discussed in section 4.

## 3 Flat Space Solutions

First we discuss flat space solutions within this general framework. As mentioned above, some of the solutions are already known but new ones can be easily obtained.

In order to study a single cosmic string solution, we take the usual cylindrically symmetric Nielsen-Olesen ansatz for $`n`$ flux units:

$$\mathrm{\Phi }=\mu f(r)e^{in\phi },A_\mu dx^\mu =A(r)d\phi $$
(3.1)
where $`\mu `$ is a second energy scale, which is generally independent of $`v`$. We will also assume that the dielectric functions depend on $`|\mathrm{\Phi }|`$ through the dimensionless ratio $`|\mathrm{\Phi }|/\mu =f`$ only. For further use we define here $`v/\mu \beta `$.

As for the boundary conditions, we will see that they will not be the same for all systems, but will rather have to be adapted to the specific system. However, since we are interested in solutions with finite energy per unit length and finite flux, the boundary conditions at infinity should ensure asymptotically vanishing energy density. The usual Nielsen-Olesen conditions should be generalized such that all three contributions (scalar gauged kinetic term, potential and Maxwell term) will vanish asymptotically. That is:

$$\underset{r\mathrm{}}{lim}U(f(r))=0,\underset{r\mathrm{}}{lim}_2(f(r))B^2(r)=0,\underset{r\mathrm{}}{lim}f^2(r)_1(f(r))(eA(r)n)^2=0$$
(3.2)
In order to proceed, we will concentrate in the O(3) model. The field equations (with $`\eta =+1`$ which we will mostly use from now on) are the following simple first order set:

$$\frac{rf^{}}{f}=eA(r)n$$
(3.3)
$$\frac{A^{}}{r}=\frac{e\mu ^2}{2(1+\beta ^2)}\frac{\beta ^2f^2}{_2(f)(1+f^2)}$$
(3.4)
where a prime denotes differentiation with respect to $`r`$.

We still have a freedom in the function $`_2`$, but we choose until further notice $`_2(f)=1`$. The simplest case, which yields the closest to the Abelian Higgs flux tube, is the case $`\beta =1`$ which corresponds to a potential with a minimum along the equator of $`\text{S}^2`$. We thus impose the following additional boundary conditions:

$$f(0)=0,A(0)=0$$
(3.5)
and the two last conditions at (3.2) may be replaced by:

$$\underset{r\mathrm{}}{lim}A(r)=n/e$$
(3.6)
Note that the boundary conditions we took enforce $`n`$ to be negative and $`A(r)`$ to be a non-positive decreasing function. The magnetic field $`B(r)`$ is non-negative and the flux (in units of $`2\pi /e`$) is $`eA(\mathrm{})=n`$, which is positive.

There is no known analytical solution of this system, but it is very easy to get numerical solutions. We will comment a little more about it at the end of this section. Fig.1 contains the field variables ($`\mathrm{\Theta }`$, $`A`$ and $`B`$) in this $`\beta =1`$ case. This is actually the solution discussed already by Mukherjee . There is another related solution, which is obtained by reflecting this solution in the $`\text{S}^2`$ target space with respect to the equator, i.e. starting at $`f=\mathrm{}`$ on the $`z`$-axis ($`r=0`$) and decreasing with $`r`$ to $`f=1`$:

$$\underset{r0}{lim}f(r)=\mathrm{},\underset{r\mathrm{}}{lim}f(r)=1$$
(3.7)
The reflection property is much more transparent in terms of the angular field variable $`\mathrm{\Theta }`$, Eq.(2.25), in terms of which the boundary conditions on the scalar field are:

$$\mathrm{\Theta }(0)=\pi ,\underset{r\mathrm{}}{lim}\mathrm{\Theta }(r)=\pi /2$$
(3.8)
This reflected solution has obviously negative magnetic flux.

Another solution, which was discussed already by Schroers , corresponds to the limiting case of potential with a minimum which does not break the U(1) symmetry. This corresponds, in our terminology, to $`\beta =0`$ while $`\mu `$ stays finite so $`v=0`$. In this case, the minimum of the potential is on the north pole of $`\text{S}^2`$. Finite energy (per unit length) solutions still exist, but they have now different boundary conditions namely:

$$\mathrm{\Theta }(0)=\pi ,A(0)=0$$
(3.9)
$$\underset{r\mathrm{}}{lim}\mathrm{\Theta }(r)=0,\underset{r\mathrm{}}{lim}B(r)=0$$
(3.10)
The last boundary condition in (3.2) is automatically satisfied as $`f0`$ asymptotically. This gives a solution, which maps the $`z`$-axis ($`r=0`$) to the south pole of $`\text{S}^2`$ and the circle $`r\mathrm{}`$ to the north pole. Since there is no symmetry breaking at spatial infinity, the magnetic flux is not quantized but rather takes continuous values.

In order to compare these solutions with the $`\beta =1`$ ones, it is more instructive to present the ”anti-Schroers” solutions, which are solutions to the same problem with the reflected potential ($`\beta \mathrm{}`$). The boundary conditions should be therefore also reflected:

$$\mathrm{\Theta }(0)=0,A(0)=0$$
(3.11)
$$\underset{r\mathrm{}}{lim}\mathrm{\Theta }(r)=\pi ,\underset{r\mathrm{}}{lim}B(r)=0$$
(3.12)
where again the last boundary condition in (3.2) is automatically satisfied, as now $`f^2_10`$ asymptotically. There are two unusual features of this family of solutions: First, as mentioned above, there is no flux quantization, and second, there does not exist a solution with $`|n|=1`$ but with $`|n|2`$ . Consequently, the energy density is maximal not on the symmetry axis, but rather on a cylindrical surface whose radius is of the order of $`1/e\mu `$. Similar behavior is seen in the analogous $`|n|2`$ Nielsen-Olesen flux tubes .

Another kind of interesting potential is a triple well potential, whose minima are at the north and south pole and at the equator of the target space. In term of the angular field $`\mathrm{\Theta }`$, it has the following simple form:

$$U(\mathrm{\Theta })=\frac{e^2\mu ^4}{128}\mathrm{sin}^2(2\mathrm{\Theta })$$
(3.13)
This potential is obtained by using $`_2(|\mathrm{\Phi }|)=\frac{1}{4}(|\mathrm{\Phi }|/\mu +\mu /|\mathrm{\Phi }|)^2`$ and $`v=\mu `$ in Eq.(2.24). Actually, $`_2`$ may be multiplied by an arbitrary constant, thus having the height of the potential free. Since the ”polar” minima (at $`\mathrm{\Theta }=0,\pi `$) do not break the U(1) symmetry, we do not expect flux quantization for solutions that approach asymptotically one of these minima. On the other hand, solutions which have $`\mathrm{\Theta }=\pi /2`$ as their asymptotic value, will have of course quantized magnetic flux.

There are thus several kinds of finite energy (per unit length) solutions in this system. For $`\mathrm{\Theta }(0)=0`$, all three kinds with asymptotic values $`\mathrm{\Theta }=0,\pi /2,\pi `$ are realized. The fields $`\mathrm{\Theta }`$, $`\overline{A}`$ and $`\overline{B}`$ in figures 3,5,7 which actually correspond to the Chern-Simons type system analyzed in sec. 5 depict also the solutions of the present system with $`\mathrm{\Theta }(\mathrm{})=\pi /2,\pi ,0`$ respectively. We will say more about it in sec. 6.

Similar solutions exist with $`\mathrm{\Theta }(0)=\pi `$. A different kind of solution starts with a ”non-polar” value (i.e. $`\mathrm{\Theta }(0)0,\pi `$), but for continuity must have $`n=0`$. The asymptotic values correspond to each of the three possible vacua.

The numerical procedure we have used is essentially the same as in previous works , and we will not elaborate about it more than the following remark regarding the treatment of boundary conditions. A common practice in numerical analysis is to replace the boundary conditions at infinity with the same boundary conditions at some finite point, $`r_{\mathrm{}}`$. This works well on the system (3.3)-(3.4), provided that $`\beta `$ is not $`0`$ or $`\mathrm{}`$. If we try to shoot from $`r_{\mathrm{}}`$ to $`0`$ in these two cases, it is easy to see (from the equations in terms of the angular field, $`\mathrm{\Theta }`$), that the solutions become constant. We are therefore forced to replace the boundary condition on $`\mathrm{\Theta }`$ with $`\mathrm{\Theta }(\mathrm{})=\beta \pm \delta `$ where $`\delta `$ is a (positive) parameter, which is to be determined numerically. It turns out, that the continuum of solutions is found by giving $`\delta `$ slightly different values.

Another remark concerns the classification of solutions according to the possible choices of boundary conditions. As is already clear from the flat space solutions, the O(3) system has a reflection symmetry with respect to the $`\text{S}^2`$ equatorial plane (i.e. $`\mathrm{\Theta }\pi \mathrm{\Theta }`$ or $`f1/f`$) together with $`\beta 1/\beta `$. We may therefore limit ourselves to solutions with $`f(0)=0`$ (which for $`\eta =+1`$ are also $`n<0`$ solutions) and classify them according to the asymptotic values $`f(\mathrm{})`$, which are the various minima of the potential function. The $`n>0`$ solutions for the same $`\eta `$ and $`\beta `$ (and hence the same potential), which start at the ”south pole” of target space, are easily obtained from the $`n<0`$ solutions with $`\beta 1/\beta `$ and further obvious changes as $`nn`$ and $`AA`$.

## 4 Gravitating Solutions

In this section we consider gravitating O(3) cosmic strings. It is in fact very easy to get the gravitating self-dual cylindrical solutions (i.e. single cosmic string solutions) in the Bogomolnyi limit, since the field equations for the matter (scalar and gauge) fields are very similar to those in flat spacetime. More precisely, they are given by the following first order equations:

$$\frac{rf^{}}{f}=\eta \left(eA(r)n\right)$$
(4.1)
which is identical with the Minkowskian one (for $`\eta =+1`$ compare (3.3)), and:

$$\frac{A^{}}{H^2r}=\frac{\eta e\mu ^2}{2(1+\beta ^2)}\frac{\beta ^2f^2}{_2(f)(1+f^2)}$$
(4.2)
Einstein equations for the metric field (2.28) get in the cylindrical case the following form:

$$\left[r\left((\mathrm{log}H)^{}\frac{8\pi G}{\eta e}B_2(f)(\mathrm{log}f)^{}\right)\right]^{}=0$$
(4.3)
Since we are interested here in cases where the matter fields tend asymptotically to their vacuum values, the geometry of space will evidently be conic with a deficit angle given by the usual relation:

$$\frac{\delta \phi }{2\pi }=\underset{r\mathrm{}}{lim}(r(\mathrm{log}H)^{})$$
(4.4)
We may obtain more information about the deficit angle by integrating equation (4.3) once and using the result at both ends, $`r=0`$ and $`r\mathrm{}`$:

$$\frac{\delta \phi }{2\pi }=\frac{8\pi G}{e}|nB(0)|_2(f(0))$$
(4.5)
where we used the boundary conditions (3.5) and (3.6) supplemented by $`H(0)=1,H^{}(0)=0`$, and the fact that $`nB(0)`$ is always (i.e. for $`\eta =\pm 1`$) positive. The factor $`B(0)_2(f(0))`$ in the last equation can be expressed in terms of the dielectric function $`_1(f)`$ using (2.17) and assuming a vanishing magnetic field as $`f\beta `$:

$$\frac{\delta \phi }{2\pi }=8\pi G\mu ^2|n|_0^\beta f_1(f)𝑑f$$
(4.6)
For the O(3) sigma model $`_1(f)`$ under consideration here, this integral can be explicitly calculated, or we may rather use the already given expression of the magnetic field - Eq.(4.2). Both ways we get:

$$\frac{\delta \phi }{2\pi }=\frac{4\pi Gv^2|n|}{1+\beta ^2}$$
(4.7)
Note that these two last equations are valid for the $`f(0)=0`$ solutions - see final note of previous section. The well-known usual Higgs result  is obtained for $`\beta =0`$, while holding $`v`$ fixed. On the other hand, we may take another limit to obtain the deficit angle for the ”anti-Schroers” solutions by $`\beta \mathrm{}`$, while holding $`\mu `$ fixed. This gives $`(\delta \phi /2\pi )_S=4\pi G\mu ^2|n|`$. This result applies also to the Schroers solutions.

In Fig.2, we show the fields ($`\mathrm{\Theta },A,B,H`$) for the case $`\beta =1`$, $`8\pi G\mu ^2=0.5`$; compare with Fig.1.

## 5 Chern-Simons Type of the Generalized Higgs Model

The possibility of a second field dependent ”dielectric function” $`_2`$, used at the end of section 3, is much more interesting from another aspect, which is a connection with $`D=3`$ Chern-Simons theory. It turns out that the self-dual solutions for a $`D=3`$ system of a gauged sigma model coupled to pure Chern-Simons theory, are related to those of our generalized Higgs system. In other words, we may replace the Maxwell term with a Chern-Simons term in the $`D=3`$ version of the action (2.1) to get a generalized Chern-Simons-Higgs (GCSH) model. We will find interesting relations between the self-dual solutions of this GCSH system and the previous one, which from now on we call generalized Maxwell-Higgs (GMH) model.

There is of course a physical difference between time-independent solutions in both models, which is the presence of electric charge and electric field in the GCSH model. However, due to the linear relation between the magnetic field and the electric charge density, a $`|\mathrm{\Phi }|`$-dependent magnetic energy density appears, which in flat space mimics exactly the Maxwellian energy density in the GMH system, as will be shown below. This sheds new light on the well-known self-dual solutions of the Chern-Simons gauged O(3) sigma model . The self-dual solutions of the Chern-Simons type of the Abelian Higgs system , fit also to this framework by the obvious choice $`_1(|\mathrm{\Phi }|)=1`$. In order to demonstrate these relations we write the $`D=3`$ action:

$`S={\displaystyle \mathrm{d}^3x\sqrt{|g|}\left(\frac{1}{2}_1(|\mathrm{\Phi }|)(D_\mu \mathrm{\Phi })^{}(D^\mu \mathrm{\Phi })U(|\mathrm{\Phi }|)+\frac{\kappa }{4}\frac{ϵ^{\lambda \mu \nu }}{\sqrt{|g|}}A_\lambda F_{\mu \nu }+\frac{1}{16\pi G_3}R\right)}`$ (5.1)
where all the geometrical terms here are 3-dimensional, and $`\kappa `$ is a positive parameter (it cannot be a field-dependent function, as in the GMH case, because it would break gauge invariance) with dimensions of 1/length. Actually, the fields here have different dimensionalities than in $`D=4`$, but the self-dual solutions of the GMH system are identical in $`D=3`$ and $`D=4`$, so the comparison becomes trivial.

The field equations of this system are:

$$_1(|\mathrm{\Phi }|)D_\mu D^\mu \mathrm{\Phi }+\frac{\mathrm{\Phi }^{}}{2|\mathrm{\Phi }|}\frac{d_1}{d|\mathrm{\Phi }|}D_\mu \mathrm{\Phi }D^\mu \mathrm{\Phi }+\frac{\mathrm{\Phi }}{|\mathrm{\Phi }|}\frac{dU}{d|\mathrm{\Phi }|}=0$$
(5.2)
$$\frac{\kappa }{2}\frac{ϵ^{\lambda \mu \nu }}{\sqrt{|g|}}F_{\lambda \mu }=j^\nu =\frac{i}{2}e_1(|\mathrm{\Phi }|)(\mathrm{\Phi }^{}(D^\nu \mathrm{\Phi })\mathrm{\Phi }(D^\nu \mathrm{\Phi })^{})$$
(5.3)
$$\frac{1}{8\pi G_3}R_{\mu \nu }+\frac{1}{2}_1(|\mathrm{\Phi }|)\left((D_\mu \mathrm{\Phi })^{}(D_\nu \mathrm{\Phi })+(D_\nu \mathrm{\Phi })^{}(D_\mu \mathrm{\Phi })\right)2U(|\mathrm{\Phi }|)g_{\mu \nu }=0$$
(5.4)
Spontaneous symmetry breaking occurs if the potential has a circle of degenerate minima , just as for the GMH system. The mass of the Higgs field is given by the same expression as for GMH - Eq.(2.6), but the mass (squared) of the gauge field is now <sup>1</sup><sup>1</sup>1Compare Deser and Yang .:

$$m_A^2=\left(\frac{e^2v^2_1(v)}{\kappa }\right)^2$$
(5.5)
In order to find static self-dual solutions in flat space, one may proceed by dropping the curvature terms and look for conditions for Eq.(2.13) to be satisfied. However, we prefer to take the more general approach of dealing with the gravitating system from the beginning. We therefore try to repeat what we did in the previous sections and require $`g_{00}=1`$ in a $`D=3`$ version of the static metric (2.7). It turns out that this cannot be done ”naively”, but some modifications are required. The reason is the well-known property of Chern-Simons vortices, which should carry also electric charge so the Maxwell tensor must contain in this case also an electric field, and the vortex carries angular momentum as well. This addition has important consequences if we treat the gravitational field as dynamical, since the angular momentum of the Chern-Simons field forces the gravitational field to be stationary and not simply static. We therefore parametrize the $`D=3`$ metric by:

$$ds^2=N^2(x{}_{}{}^{k})(dt+L_i(x{}_{}{}^{l})dx{}_{}{}^{i})^2\gamma _{ij}(x{}_{}{}^{k})dx{}_{}{}^{i}dx^j$$
(5.6)
and write for the gauge field:

$$A_\mu dx^\mu =A_0(x{}_{}{}^{k})dt+A_i(x{}_{}{}^{k})dx^i$$
(5.7)
In order to get the field equations for stationary solutions, we introduce also the following notation:

$`\overline{A}_i`$ $`=`$ $`A_iA_0L_i`$
$`\overline{D}_i`$ $`=`$ $`_iie\overline{A}_i`$
$`L_{ij}`$ $`=`$ $`_iL_j_jL_i=\sqrt{|\gamma |}ϵ_{ij}\mathrm{}`$
$`\overline{F}_{ij}`$ $`=`$ $`_i\overline{A}_j_j\overline{A}_i=\sqrt{|\gamma |}ϵ_{ij}\overline{B}`$ (5.8)
and compute the components of the Ricci tensor:

$`R_{00}=N_i^iN{\displaystyle \frac{N^4}{2}}\mathrm{}^2`$
$`R_0^i={\displaystyle \frac{1}{2N\sqrt{|\gamma |}}}ϵ^{ij}_j(N^3\mathrm{})`$ (5.9)
$`R^{ij}={\displaystyle \frac{1}{N}}^i^jN+{\displaystyle \frac{1}{2}}(R(\gamma )N^2\mathrm{}^2)\gamma ^{ij}`$
where $`_i`$ is the covariant derivative with respect to the two-dimensional metric $`\gamma _{ij}`$ and $`R(\gamma )`$ the corresponding Ricci scalar.

Now we impose $`N=1`$ and find that the $`(00)`$ component of Einstein equations (5.4) will be satisfied only with the following condition on the potential:

$$U(|\mathrm{\Phi }|)=\frac{e^2}{2}|\mathrm{\Phi }|^2_1(|\mathrm{\Phi }|)(A_0)^2\frac{\mathrm{}^2}{32\pi G_3}$$
(5.10)
From the $`(ij)`$ Einstein equations (or better from the $`G^{ij}`$ equations), we find the form of the self-duality equation in this case:

$$\overline{D}_i\mathrm{\Phi }=i\eta \sqrt{|\gamma |}ϵ_{ij}\gamma ^{jk}\overline{D}_k\mathrm{\Phi }$$
(5.11)
as well as the following expression for the two dimensional Ricci scalar:

$$R(\gamma )=\mathrm{}^28\pi G_3\left(_1(|\mathrm{\Phi }|)\gamma ^{ij}(\overline{D}_i\mathrm{\Phi })^{}(\overline{D}_j\mathrm{\Phi })+4U(|\mathrm{\Phi }|)\right)$$
(5.12)
Due to the self-duality, the other field equations simplify as follows. The non-diagonal Einstein equations give:

$$_i\mathrm{}=16\pi G_3\eta e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)A_0_i|\mathrm{\Phi }|$$
(5.13)
The spatial Chern-Simons equations give:

$$_iA_0=\frac{\eta e}{\kappa }|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)_i|\mathrm{\Phi }|$$
(5.14)
which is used in the time component of the Chern-Simons equations to obtain an expression for $`\overline{B}`$:

$$\overline{B}=(\mathrm{}+\frac{e^2}{\kappa }|\mathrm{\Phi }|^2_1(|\mathrm{\Phi }|))A_0$$
(5.15)
Finally, the equation for the Higgs field gives another functional relation between the potential and the other quantities:

$$\eta e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)\overline{B}+\frac{dU}{d|\mathrm{\Phi }|}\frac{e^2}{2}(A_0)^2\frac{d}{d|\mathrm{\Phi }|}(|\mathrm{\Phi }|^2_1)=0$$
(5.16)
Actually it is not an independent equation, since by substitution of (5.10) in (5.16), one may get back (5.15).

The two differential equations for $`\mathrm{}(x^k)`$ and $`A_0(x^k)`$ may be also converted into equations for $`\mathrm{}(|\mathrm{\Phi }|)`$ and $`A_0(|\mathrm{\Phi }|)`$:

$$\frac{d\mathrm{}}{d|\mathrm{\Phi }|}=16\pi G_3\eta e|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)A_0,$$
(5.17)
$$\frac{dA_0}{d|\mathrm{\Phi }|}=\frac{\eta e}{\kappa }|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)$$
(5.18)
and we find a simple expression of $`\mathrm{}`$ in terms of $`A_0`$:

$$\mathrm{}=8\pi G_3\kappa (cA_0^2).$$
(5.19)
where $`c`$ is an integration constant, which should be non-negative for solutions with finite angular momentum (where $`\mathrm{}`$ vanishes asymptotically).

As in the gravitating GMH model, the right hand side of (5.12) is a two-dimensional divergence and we find:

$$R(\gamma )=16\pi G_3\kappa _i\left(\frac{1}{\eta e}\gamma ^{ij}A_0_j\mathrm{log}|\mathrm{\Phi }|+\frac{cϵ^{ij}}{\sqrt{|\gamma |}}L_j\right)$$
(5.20)
If we use a conformally flat metric, (5.20) can be rewritten as:

$$\delta ^{ij}_i\left(_j\mathrm{log}H\frac{8\pi G_3\kappa }{\eta e}A_0_j\mathrm{log}|\mathrm{\Phi }|\right)=8\pi G_3\kappa cH^2\mathrm{}$$
(5.21)
Now we consider some special cases that allow self-dual solutions. We will then concentrate in rotationally symmetric solutions, i.e. gravitating vortices.

### 5.1 $`_1(\mathrm{\Phi }|)=1`$

This is the Chern-Simons version of the usual Higgs model, which we may consider either in flat space or coupled to gravity. Integration of Eq.(5.18) is trivial and gives:

$$A_0=\frac{\eta e}{2\kappa }(v^2|\mathrm{\Phi }|^2)$$
(5.22)
where $`v`$ is an integration constant. In flat space $`\mathrm{}=0`$ and Eq.(5.10) gives immediately the well-known  sixth order potential:

$$U(|\mathrm{\Phi }|)=\frac{e^4}{8\kappa ^2}|\mathrm{\Phi }|^2(v^2|\mathrm{\Phi }|^2)^2$$
(5.23)
In the case of a gravitating system, we solve for $`\mathrm{}(|\mathrm{\Phi }|)`$ and find:

$$\mathrm{}(|\mathrm{\Phi }|)=\frac{2\pi G_3e^2}{\kappa }\left(\sigma (v^2|\mathrm{\Phi }|^2)^2\right)$$
(5.24)
where $`\sigma =4\kappa ^2c/e^2`$. The potential is easily found to be:

$$U(|\mathrm{\Phi }|)=\frac{e^4}{8\kappa ^2}\left[|\mathrm{\Phi }|^2(v^2|\mathrm{\Phi }|^2)^2\pi G_3\left(\sigma (v^2|\mathrm{\Phi }|^2)^2\right)^2\right]$$
(5.25)
which has an unbounded additional term, thus rendering the whole system possibly unstable. Notwithstanding the possible instability, this potential was discussed by several authors  and vortex solutions were also obtained. There is a simple and natural way to cure the ill-behaved potential, which is generalizing this system to the O(3) sigma model.

### 5.2 $`_1(\mathrm{\Phi }|)=1/(1+|\mathrm{\Phi }|^2/\mu ^2)^2`$

This is the Chern-Simons gauged O(3) sigma model considered in flat background by some authors . Within the present framework the analysis is straightforward. First we integrate Eq.(5.18) and find:

$$A_0=\frac{\eta e\mu ^2}{2\kappa (1+\beta ^2)}\frac{\beta ^2|\mathrm{\Phi }|^2/\mu ^2}{1+|\mathrm{\Phi }|^2/\mu ^2}$$
(5.26)
In flat space $`\mathrm{}=0`$ and Eq.(5.10) gives immediately the following potential:

$$U(|\mathrm{\Phi }|)=\frac{e^4\mu ^4}{8\kappa ^2(1+\beta ^2)^2}\frac{|\mathrm{\Phi }|^2\left(\beta ^2|\mathrm{\Phi }|^2/\mu ^2\right)^2}{(1+|\mathrm{\Phi }|^2/\mu ^2)^4}$$
(5.27)
Some authors have already discussed the self-dual solutions in flat space for special cases of the potential, which from this point of view are just special values of $`\beta `$: $`\beta =0`$  and $`\beta =1`$ . This general form of the potential appeared already (in a different parametrization) in Ref., but only the case $`\beta =0`$ was discussed there. Kimm et al.  have studied the general case in flat space and found three kinds of flux tube solutions characterized by the asymptotic value of $`f(r)`$, which (for $`f(0)=0`$) may be either $`0`$, $`\beta `$ or $`\mathrm{}`$.

If this Chern-Simons gauged O(3) model is coupled to gravity, we find that the ill-behaved potential from the Higgs system becomes now bounded from below. First we integrate (5.17) either directly or using (5.19) and (5.26) to get:

$$\mathrm{}(|\mathrm{\Phi }|)=\frac{2\pi G_3e^2\mu ^4}{\kappa (1+\beta ^2)^2}\left(\lambda \frac{(\beta ^2|\mathrm{\Phi }|^2/\mu ^2)^2}{(1+|\mathrm{\Phi }|^2/\mu ^2)^2}\right)$$
(5.28)
where the role of the integration constant is played by $`\lambda `$ defined by $`\lambda =(1+\beta ^2)^2\sigma /\mu ^4`$. Then we find from (5.10) a new potential with $`\lambda `$ as an additional (non-negative) free parameter:

$$U(|\mathrm{\Phi }|)=\frac{e^4\mu ^4}{8\kappa ^2(1+\beta ^2)^2}\left[\frac{|\mathrm{\Phi }|^2\left(\beta ^2|\mathrm{\Phi }|^2/\mu ^2\right)^2}{(1+|\mathrm{\Phi }|^2/\mu ^2)^4}\frac{\pi G_3\mu ^4}{(1+\beta ^2)^2}\left(\lambda \frac{(\beta ^2|\mathrm{\Phi }|^2/\mu ^2)^2}{(1+|\mathrm{\Phi }|^2/\mu ^2)^2}\right)^2\right]$$
(5.29)
This potential is clearly bounded from below, thus restoring the stability. Moreover, it has always a local minimum at $`|\mathrm{\Phi }|=v=\beta \mu `$, which breaks the U(1) local symmetry. The extremal points at $`|\mathrm{\Phi }|=0`$ and $`|\mathrm{\Phi }|\mathrm{}`$, which are U(1)-symmetric, may be either minima or maxima depending on the following conditions:

$`\beta ^2(1+\beta ^2)+4\pi G_3\mu ^2(\beta ^4\lambda )>0\text{minimum at}|\mathrm{\Phi }|=0`$
$`1+\beta ^2+4\pi G_3\mu ^2(1\lambda )>0\text{minimum at}|\mathrm{\Phi }|\mathrm{}`$ (5.30)
The potential of the gravitating Chern-Simons type of the usual Higgs system, is obtained in the limit $`\mu \mathrm{}`$ (with $`\lambda \mu ^4`$ and $`\beta \mu `$ kept finite).

Actually a special case of this potential for $`\lambda =0`$ was recently obtained by Abou-Zeid and Samtleben  from a different direction of three dimensional $`N=2`$ supergravity theory. We will see that $`\lambda 0`$ yields interesting solutions as well.

Now we concentrate in rotationally symmetric solutions, i.e. we use the Nielsen-Olesen ansatz including the time component of the gauge potential:

$$\mathrm{\Phi }=\mu f(r)e^{in\phi },A_\mu dx^\mu =A_0(r)dt+A(r)d\phi $$
(5.31)
with the additional requirement that all the metric components depend on $`r`$ only.

The field equations for self-dual rotationally symmetric solutions are quite similar to the ones of the GMH model with two additional equations: For $`A_0`$ and for the metric component $`L_\phi (r)`$ ($`L_r`$ now vanishes). The dependence of $`A_0`$ on $`f`$ is obvious from Eq.(5.26). For $`L_\phi (r)`$ we have:

$$\frac{L_\phi ^{}}{H^2r}=\mathrm{}$$
(5.32)
where the dependence of $`\mathrm{}`$ on $`f`$ is obtained from (5.28). The self-duality equation for the scalar field is now ($`\eta =+1`$ as usual):

$$\frac{rf^{}}{f}=e\overline{A}(r)n$$
(5.33)
where for brevity we denote $`\overline{A}=AA_0L_\phi `$. The equation for $`\overline{A}(r)`$ may be easily obtained from (5.15), expressing $`\mathrm{}`$ and $`A_0`$ in terms of $`f`$:

$$\frac{\overline{A}^{}}{H^2r}=\frac{e^3\mu ^4}{2\kappa ^2(1+\beta ^2)}\frac{\beta ^2f^2}{1+f^2}\left[\frac{2\pi G_3\mu ^2}{(1+\beta ^2)^2}\left(\lambda \frac{(\beta ^2f^2)^2}{(1+f^2)^2}\right)+\frac{f^2}{(1+f^2)^2}\right]$$
(5.34)
Einstein equations for the metric field (5.21) reduce in this case to the following:

$$\left[r\left((\mathrm{log}H)^{}\frac{8\pi G_3\kappa }{e}A_0(\mathrm{log}f)^{}\right)8\pi G_3\kappa cL_\phi \right]^{}=0$$
(5.35)
The reflection symmetry observed in the GMH system exists in the GCSH system as well, provided $`\lambda `$ is also rescaled according to its $`\beta `$ dependence. We may therefore limit ourselves to solutions with $`f(0)=0`$ as before. Clearly there are three kinds of solutions classified by the three possible values which $`f(\mathrm{})`$ may take, according to the different minima of the potential. However, these three values $`U(0)`$, $`U(v)`$ and $`U(\mathrm{})`$ are generally different from each other, and each of them vanishes for a different value of $`\lambda `$. Since vanishing value of the potential minimum is an essential property of our localized solutions, we encounter here a situation which differs from that in flat space , where all three possible boundary conditions are realized for the same potential. In the gravitating case, each kind is realized for a different value of $`\lambda `$. The flux tube solutions with quantized flux and the usual boundary condition $`f(\mathrm{})=\beta `$ exist for $`\lambda =0`$ only. Representative solutions are depicted in Figs.3-4. The two other solution types with $`f(\mathrm{})=0,\mathrm{}`$ exist for $`\lambda =\beta ^4,1`$, respectively and are depicted in Figs.5-8. We may therefore summarize the situation by the following relations:

$`\lambda =\beta ^4`$ $``$ $`\underset{r\mathrm{}}{lim}f(r)=0`$
$`\lambda =0`$ $``$ $`\underset{r\mathrm{}}{lim}f(r)=\beta `$ (5.36)
$`\lambda =1`$ $``$ $`\underset{r\mathrm{}}{lim}f(r)=\mathrm{}`$
A special case is $`\lambda =\beta =1`$ which allows both boundary conditions ($`f(\mathrm{})=0,\mathrm{}`$) for the same potential.

Note that all three kinds of boundary conditions ensure also a vanishing asymptotic value of the angular momentum density $`\mathrm{}`$, which is a necessary condition for a finite angular momentum. The gauge potential $`A_0`$ need not always vanish asymptotically. It does only in the $`f(\mathrm{})=\beta `$ case.

For all these three kinds of boundary conditions, the matter fields tend asymptotically to their vacuum values. Thus asymptotically, the geometry of space-time will be rotating conic metric. The angular deficit may be easily obtained from (5.35) to be:

$$\frac{\delta \phi }{2\pi }=\frac{8\pi G_3\kappa }{e}\left[nA_0(0)+A_0(\mathrm{})(eA(\mathrm{})n)\right]$$
(5.37)
For vanishing $`c`$ (or $`A_0(\mathrm{})`$ or $`\lambda `$, i.e. quantized flux) we find:

$$\frac{\delta \phi }{2\pi }=\frac{8\pi G_3\kappa }{e}nA_0(0)=8\pi G_3\mu ^2|n|_0^\beta f_1(f)𝑑f$$
(5.38)
that is the same as the result (4.6) for the GMH model. For the O(3) case under consideration here, we find the same expression as we had in the case of O(3)-GMH:

$`{\displaystyle \frac{\delta \phi }{2\pi }}`$ $`=`$ $`{\displaystyle \frac{4\pi G_3v^2|n|}{1+\beta ^2}},\underset{r\mathrm{}}{lim}f(r)=\beta `$ (5.39)
The angular deficit for the other solutions without flux quantization (i.e. non-vanishing $`c`$ or $`\lambda `$) can be easily calculated in the O(3) model from (5.37):

$`{\displaystyle \frac{\delta \phi }{2\pi }}`$ $`=`$ $`{\displaystyle \frac{4\pi G_3\mu ^2\beta ^2e|A(\mathrm{})|}{1+\beta ^2}},\underset{r\mathrm{}}{lim}f(r)=0`$ (5.40)
$`{\displaystyle \frac{\delta \phi }{2\pi }}`$ $`=`$ $`4\pi G_3\mu ^2|{\displaystyle \frac{eA(\mathrm{})}{1+\beta ^2}}n|,\underset{r\mathrm{}}{lim}f(r)=\mathrm{}`$ (5.41)
Note the dependence on the unquantized magnetic flux represented here by $`A(\mathrm{})`$. Thus we have in these cases a continuum of values for the angular deficit for a given potential and $`n`$ and $`f(\mathrm{})`$ values.

We end this section by computing the angular momentum $`J`$ of the rotationally symmetric solutions, which we get by integrating $`T_\phi ^0`$ over all two-space. The angular momentum also determines the asymptotic value of the metric component $`L_\phi (r)`$ such that:

$$(1\frac{\delta \phi }{2\pi })L_\phi (\mathrm{})=4GJ$$
(5.42)
This relation is most easily obtained in a Gaussian normal coordinate system, in which the line element is written as: $`ds^2=d\rho ^2+h_{\alpha \beta }dx^\alpha dx^\beta `$ ($`x^\alpha =t,\phi `$). Since $`h_{\alpha \beta }`$ depends only on $`\rho `$, the components of the Ricci tensor have the following simple form in terms of the extrinsic curvature $`K_\beta ^\alpha `$:

$$R_\beta ^\alpha =\frac{1}{\sqrt{|h|}}\frac{}{\rho }(\sqrt{|h|}K_\beta ^\alpha ),K_\beta ^\alpha =\frac{1}{2}h^{\alpha \gamma }\frac{h_{\gamma \beta }}{\rho }$$
(5.43)
It turns out that we can compute the angular momentum directly without even the self-duality assumption, due to the simple identity which holds for solutions of the flux tube form (5.31):

$$eT_\phi ^0=(neA(r))j^0$$
(5.44)
Now we use the time component of the Chern-Simons equation (5.3), to write $`j^0`$ in terms of $`B`$ and get:

$$J=\frac{2\pi \kappa }{e}_0^{\mathrm{}}(neA(r))A^{}𝑑r=\frac{\pi \kappa }{e}(2neA(\mathrm{}))A(\mathrm{})$$
(5.45)
Solutions which break the U(1) symmetry asymptotically have quantized angular momentum of:

$$J=\frac{\pi \kappa }{e^2}n^2$$
(5.46)
while otherwise, there is no flux quantization and we cannot say about the angular momentum more than the right hand side of (5.45). Note that the angular momentum in (5.46) is a very general result independent of all details of the theory, except the existence of symmetry breaking vacuum in the potential. We note further that all solutions have the same sign of angular momentum, which is a manifestation of parity violation.

## 6 Correspondence between Maxwell and Chern-Simons Type of the Generalized Higgs Model

We note that in flat space, the electric energy density of the Chern-Simons field plays an equivalent role to that of the magnetic energy density in the GMH model. Moreover, due to (5.15), there is a relation between the magnetic and electric terms, thus enabling us to eliminate the scalar (electric) potential from the equations and to give them a form identical to those of the GMH model.

Comparison of the Bogomolny constraints in both cases, equations (2.17) and (5.18), gives a relation between the magnetic field in the GMH model and the time component of the gauge potential in the GCSH model:

$$_2B=\kappa A_0$$
(6.1)
Further comparison of the magnetic field appearing in this equation with the one in (5.15) (with $`\mathrm{}=0`$) or the potentials ((2.11) and (5.10) (with $`\mathrm{}=0`$)), yields the following characteristic relation for the generalized Maxwell-Higgs model:

$$\frac{\kappa ^2}{_2(|\mathrm{\Phi }|)}=e^2|\mathrm{\Phi }|^2_1(|\mathrm{\Phi }|)$$
(6.2)
Thus, all the flat space self-dual solutions to the $`D=3`$ system, studied in section 5, are also self-dual solutions to the $`D=4`$ GMH model of section 3, provided we use $`_2`$ consistent with (6.2). Consequently, the curves in figures 3,5,7 represent also the corresponding fields of the GMH model with the triple well potential (3.13) and the appropriate dielectric function mentioned together with that potential.

However, the presence of an electric field in the $`D=3`$ Chern-Simons theory, in addition to the magnetic field, is responsible for the existence of angular momentum, which calls for corrections to this simple correspondence when gravity is considered dynamical. It is straightforward (although tedious) to show that there does not exist spinning gravitating self-dual flux tubes in the GMH model. Therefore, we have only static versus stationary correspondence, and all we can hope for is equivalence of the two-dimensional spatial metrics in addition to the scalar and vector correspondence. This turns out to be the case, provided we also impose $`c=\sigma =\lambda =0`$.

Equation (6.1) is evidently still valid in the gravitating case and should be accompanied by $`B=\overline{B}`$ where $`B`$ is the GMH magnetic field defined by (2.9) and $`\overline{B}`$ \- in (5.8). Thus we get the relation instead of (6.2):

$$\frac{\kappa ^2}{_2(|\mathrm{\Phi }|)}=e^2|\mathrm{\Phi }|^2_1(|\mathrm{\Phi }|)+\kappa \mathrm{}$$
(6.3)
Note that now the potentials are not equal in both cases, but there is an $`\mathrm{}`$-dependent difference:

$$U_M=U_{CS}\frac{\mathrm{}^2}{32\pi G_3}$$
(6.4)
As for the metric tensor, it is easy to see that both expressions for the Ricci scalar, equations (2.27) and (5.20), become identical due to equation (6.1), thus resulting in equal 2-metrics in the two cases.

Actually, those $`\mathrm{}`$-dependent modifications should be expressed in terms of functions that appear in the Lagrangian, and should be considered as functions of $`|\mathrm{\Phi }|`$. We will not do this explicitly, since it is straightforward and does not add any more insight into the physical picture.

## 7 Topological Charge

The quantization of magnetic flux is intimately connected with the fact that the minimum of the potential breaks the local U(1) symmetry. Consequently, the magnetic flux number is the index, or winding number, of the map defined by the scalar fields from infinite distance from the flux tube to the vacuum manifold.

If however, the target space is compact, there exists a further possibility of homotopy classification of the maps defined by the fields from ”real” space to target space. We can use, as a ground of the present discussion, the well-known notions from sigma models defined in three spacetime dimensions . The GCSH model is naturally three dimensional, but the GMH model should be thought of here as three dimensional as well. In this case, the scalar fields map the two dimensional physical space into target space. We have therefore an identically conserved topological current defined by:

$$K^\lambda =\frac{ϵ^{\lambda \mu \nu }}{2\mathrm{\Omega }_T\sqrt{|g|}}\left(i_1(|\mathrm{\Phi }|)(D_\mu \mathrm{\Phi })^{}(D_\nu \mathrm{\Phi })+e(|\mathrm{\Phi }|)F_{\mu \nu }\right)$$
(7.1)
where $`\mathrm{\Omega }_T`$ is the target space volume, and $`(|\mathrm{\Phi }|)`$ is defined by:

$$\frac{d}{d|\mathrm{\Phi }|}=|\mathrm{\Phi }|_1(|\mathrm{\Phi }|)$$
(7.2)
This is a gauge invariant generalization of the standard non-gauged sigma model current, whose time component integral over all space is the winding number, or index, of the map defined by the sigma model fields. Note that $`(|\mathrm{\Phi }|)`$ is connected to $`𝒜(|\mathrm{\Phi }|)`$ (defined in Eq.(2.30)) by $`e(|\mathrm{\Phi }|)=|\mathrm{\Phi }|d𝒜/d|\mathrm{\Phi }|`$, up to a possible additive constant. By the field equations of the GCSH model, it is also related to $`A_0`$ by $`\kappa A_0(|\mathrm{\Phi }|)=\eta e(|\mathrm{\Phi }|)`$, up to an additive constant. Unlike previous authors , we fix the additive constant of $``$ such that $``$ vanishes as $`|\mathrm{\Phi }|\mathrm{}`$, irrespective of the boundary conditions imposed on the solutions. This way, we can treat uniformly all sigma models differing only by the potential term. For the special O(3) ”dielectric function”, given in (2.22), we may write the topological charge (and current) in terms of $`_1(|\mathrm{\Phi }|)`$ only, since $`(|\mathrm{\Phi }|)=\mu ^2\sqrt{_1(|\mathrm{\Phi }|)}/2`$. In this case, we have also $`\mathrm{\Omega }_T=\pi \mu ^2`$ and the apparent dependence on the scale $`\mu `$ disappears.

The fact that this gauge invariant current has still a topological meaning is clearly seen by the fact that it is actually the non-gauged sigma model current with an addition, which is a divergence of an anti-symmetric tensor:

$$K^\lambda =\frac{ϵ^{\lambda \mu \nu }}{2\mathrm{\Omega }_T\sqrt{|g|}}\left(i_1(|\mathrm{\Phi }|)(_\mu \mathrm{\Phi })^{}(_\nu \mathrm{\Phi })+2e_\mu ((|\mathrm{\Phi }|)A_\nu )\right)$$
(7.3)
Indeed, the first term of (7.3) can also be written as a divergence of an anti-symmetric tensor, so we have the alternative form:

$$K^\nu =_\mu K^{\mu \nu },K^{\mu \nu }=\frac{(|\mathrm{\Phi }|)ϵ^{\mu \nu \lambda }}{2\mathrm{\Omega }_T\sqrt{|g|}}\left(\frac{i}{|\mathrm{\Phi }|^2}(\mathrm{\Phi }^{}_\lambda \mathrm{\Phi }\mathrm{\Phi }_\lambda \mathrm{\Phi }^{})+2eA_\lambda \right)$$
(7.4)
Integration of the time component of the topological current provides us with the topological charge $`T`$, which characterizes the solutions from the homotopy point of view. For our stationary solutions, we write the two equivalent expressions:

$`T={\displaystyle \frac{1}{2\mathrm{\Omega }_T}}{\displaystyle \mathrm{d}^2x\left(iϵ^{ij}_1(|\mathrm{\Phi }|)(D_i\mathrm{\Phi })^{}(D_j\mathrm{\Phi })2e\sqrt{|\gamma |}B(|\mathrm{\Phi }|)\right)}`$
$`={\displaystyle \frac{i}{2\mathrm{\Omega }_T}}{\displaystyle \mathrm{d}^2xϵ^{ij}_1(|\mathrm{\Phi }|)(_i\mathrm{\Phi })^{}(_j\mathrm{\Phi })}{\displaystyle \frac{e}{\mathrm{\Omega }_T}}(w){\displaystyle _r\mathrm{}}A_i𝑑x^i`$ (7.5)
where $`w`$ is the asymptotic value of $`|\mathrm{\Phi }|`$. We stress the difference between $`w`$, which corresponds to a property of a solution, and $`v`$ which is a parameter in the potential. $`w`$ need not be equal to $`v`$ (although it may be), as happens for solutions of the GCSH model discussed in section 5, or solutions to the GMH model with the potential (3.13), which may have all three asymptotic values $`w=0,v,\mathrm{}`$.

Note that $`T`$ cannot be expressed as a surface integral (actually line integral in $`D=2+1`$) since the winding number, i.e. the (non-gauged) sigma model contribution of $`K^{i0}`$, does not satisfy Stokes theorem. The second term is however a boundary term, and we note that it vanishes for solutions for which $`|\mathrm{\Phi }|\mathrm{}`$ as $`r\mathrm{}`$. In the O(3) model, these are solutions which do not break the U(1) symmetry asymptotically. Another possibility is solutions which do break the U(1) symmetry asymptotically, and thus do not cover all target space. This kind of solutions have non-integer winding number, but the topological charge has also a magnetic flux contribution, which may compensate for this fact.

This is exactly what happens for the self-dual cylindrically symmetric solutions of the O(3) GMH model, where we find that in all cases $`T=n`$. If the solution tends to a symmetry breaking vacuum, the boundary term of (7.5) does not vanish and evaluates to $`n/(1+\omega ^2)`$, where we define $`\omega =w/\mu `$. The winding number, which is the first term, is easily found to be $`n\omega ^2/(1+\omega ^2)`$. They both add up to give $`T=n`$. In the case $`\omega \mathrm{}`$, where the solution does not break asymptotically the U(1) symmetry, the flux contribution vanishes but the first is an integer and we get $`T=n`$ again. The third possibility, $`\omega =0`$ is however an exception: In this case, the winding number vanishes but we have a flux contribution to get $`T=eA(\mathrm{})`$.

Another property, which we can note here, is the relation with the Euler number of the 2-surface generated by the self-dual GMH solutions. The Euler number is given by:

$$\chi =\frac{1}{4\pi }d^2x\sqrt{|\gamma |}R(\gamma )$$
(7.6)
while for (say) self-dual solutions ($`\eta =1`$), we may write the topological charge as:

$`T={\displaystyle \frac{1}{2\mathrm{\Omega }_T}}{\displaystyle \mathrm{d}^2x\sqrt{|\gamma |}\left(_1(|\mathrm{\Phi }|)\gamma ^{ij}(D_i\mathrm{\Phi })^{}(D_j\mathrm{\Phi })+2_2(|\mathrm{\Phi }|)B^2\right)}{\displaystyle \frac{e}{\mathrm{\Omega }_T}}(v){\displaystyle _r\mathrm{}}A_i𝑑x^i`$ (7.7)
By the field equations (2.11) and (2.14), we see that the first term is just proportional to the Euler number, so we have a simple relation between the topological charge, the Euler number and the magnetic flux, which we denote here by $`\mathrm{\Psi }`$:

$$T+\frac{\chi }{4G\mathrm{\Omega }_T}\frac{e(v)\mathrm{\Psi }}{\mathrm{\Omega }_T}=0$$
(7.8)
For the O(3) system, (7.8) yields the following relation:

$$2\pi T+\frac{\chi }{2G\mu ^2}+\frac{e\mathrm{\Psi }}{1+\beta ^2}=0$$
(7.9)
For asymptotically conic space, the Euler number is related to the deficit angle by $`\chi =\delta \phi /2\pi `$, so Eq.(7.8) is easily verified by using (4.7).

Next we turn to the GCSH model, where things are quite similar. The topological current and charge are still given by Eqs.(7.1)-(7.5). However, since the self-duality condition is now modified by replacing $`A_i`$ by $`\overline{A}_i`$, it is useful to write the topological charge also in a way that is ready for direct use of the modified self-duality. We therefore write the topological charge density as:

$$K^0=\frac{ϵ^{ij}}{2\mathrm{\Omega }_T\sqrt{|g|}}\left(i_1(|\mathrm{\Phi }|)(\overline{D}_i\mathrm{\Phi })^{}(\overline{D}_j\mathrm{\Phi })+2e_i((|\mathrm{\Phi }|)A_0L_j)\right)+\frac{e(|\mathrm{\Phi }|)\overline{B}}{\mathrm{\Omega }_T}$$
(7.10)
Note however that the surface term vanishes only for quantized flux solutions ($`c=\lambda =0`$, $`A_0(w)=0`$). If we proceed along similar lines as for the GMH case, we easily find the same relation between the topological charge, Euler number and the magnetic flux, i.e. Eq.(7.8). It simplifies again for the O(3) model to (7.9).

Finally, we use this relation in order to calculate the topological charge for the solutions with quantized as well as non-quantized flux, which is obtained directly from Eq.(7.5):

$`T`$ $`=`$ $`\{\begin{array}{cc}eA(\mathrm{})\hfill & ,w=0\hfill \\ n\hfill & ,w=v\hfill \\ n\hfill & ,w\mathrm{}\hfill \end{array}`$ (7.14)
A nice check of our results is that we can easily reproduce by substitution in (7.9) the expressions (5.39)-(5.41) for the deficit angles.

## 8 Outlook

We have shown that the inclusion of dielectric functions gives the possibility of dramatic changes in the long range behavior of scalar, gauge and metric fields around a cosmic string or vortex. This was obtained even for very simple and symmetric models, like O(3) with the potential (3.13). We believe that the astrophysical and cosmological consequences deserve further studies.

Another continuation and generalization of our work is to consider theories with a more extended field content. One natural possibility is to gauge a U(1)$`\times `$U(1) subgroup of the global symmetry group of a CP(2) non-linear sigma model noticing that CP(2) can be parametrized by two complex coordinates. One can therefore expect that in analogy with the way the O(3) model (with $`\text{S}^2`$ CP(1) as target space) generalizes the ordinary cosmic string, the gauged CP(2) model generalizes the superconducting cosmic string .