Effective potential and dynamical symmetry breaking up to five loops
in a massless abelian Higgs model
A. G. Quinto
andresarturogomezquinto@mail.uniatlantico.edu.coR. Vega Monroy
ricardovega@mail.uniatlantico.edu.coFacultad de Ciencias Básicas, Universidad del Atlántico Km. 7, Via
a Pto. Colombia, Barranquilla, Colombia
A. F. Ferrari
alysson.ferrari@ufabc.edu.brCentro de Ciências Naturais e Humanas, Universidade Federal do ABC–
UFABC, Rua Santa Adélia, 166, 09210-170, Santo André, SP, Brazil
A. C. Lehum
lehum@ufpa.brFaculdade de Física, Universidade Federal do Pará, 66075-110, Belém,
Pará, Brazil
Abstract
In this paper, we investigate the application of the Renormalization
Group Equation (RGE) in the determination of the effective potential
and the study of Dynamical Symmetry Breaking (DSB) in a massless Abelian
Higgs (AH) model with an -component complex scalar field in
dimensional spacetime. The classical Lagrangian of this model has
scale invariance, which can be broken by radiative corrections to
the effective potential. It is possible to calculate the effective
potential using the RGE and the renormalization group functions that
are obtained directly from loop calculations of the model and, using
the leading logs approximation, information about higher loop orders
can be included in the effective potential thus obtained. To show
this, we use the renormalization group functions reported in the literature,
obtained with a four loop calculation, and obtain a five loop approximation
to the effective potential, in doing so, we have to properly take
into account the fact that the model has multiple scales, and convert
the functions that were originally calculated in the minimal subtraction
(MS) renormalization scheme to another scheme which is adequate for
the RGE method. This result is then used to study the DSB, and we
present evidence for a rich structure of classical vacua, depending
on the value of gauge coupling constant and number of scalar fields,
which are considered as free parameters.
I Introduction
The Abelian Higgs (AH) model is one of the most fundamental field
theories in both condensed matter and particle physics. As an example,
it is the prime textbook example for the superconducting transition
and the Anderson-Higgs mechanism (ZinnJustin1996, ; Peskin2018, ; Altland2010a, ; Herbut2007, ).
The AH model features a complex scalar field coupled to an
gauge field, and it displays two distinct phases separated by a sharp
transition: the symmetric phase and the phase with spontaneously broken
symmetry. In the context of superconductors, the symmetric phase is
related to the normal metallic state and the broken one to the superconducting
Meissner state. The transition between spontaneously broken and symmetric
phases is characterized by a dimensionful parameter that serves as
an order parameter.
Starting from a classical scale invariant Lagrangian, Coleman and
Weinberg (CW) demonstrated in (Coleman1973, ) how the order
parameter could be generated by radiative corrections, i.e., the spontaneous
symmetry breaking can occur as a dynamical mechanism, the radiative
corrections being entirely responsible for the appearance of the nontrivial
minima of the effective potential. This Dynamical Symmetry Breaking
(DSB) is a key concept that has many applications in particle physics (Elias2003, ; Chishtie2006, ; Chishtie2011, ; Steele2013, )
and condensate matter systems (Liu2013, ; Uchino2014, ; Burmistrov2020, ; Chodos1994, ; A.G.Quinto2021, ).
In order to study the CW mechanism, we need to calculate the effective
potential, a powerful tool to explore many aspects of the low-energy
sector of a quantum field theory. In many cases, the one-loop approximation
is good enough, but it can be improved by adding higher order contributions
to the loop expansion. A standard tool for improving a perturbative
calculation is the Renormalization Group Equation (RGE), which, together
with a reorganization of the perturbative series in terms of leading
logs, have been shown to be very effective in several instances (AHMADY2003221, ; PhysRevD.72.037902, ; Souza2020, ; Quinto2016, ; Dias2014, ; CHISHTIE2007, ; A.G.Quinto2021, ).
We refer the reader to section 3 in (Quinto2016, ) for a short
review of the method, and (Elias2003, ; Chishtie2006, ; Chishtie2011, ; Steele2013, )
for some of the interesting results that have been reported with the
use of the RG improvement, in the context of a scale-invariant approximation
of the Standard Model.
In this work we study the behavior of effective potential in a massless
AH model with -component complex scalar field in dimensional
space-time, which is scale invariant at the classical level. We observed
that the effective potential, computed up to five loops, leads to
an interesting phase structure arising from DSB. This result was achieved
by the use of RGE with the help of renormalization group functions,
and , which were calculated
up to four loops in the minimal subtraction (MS) scheme in (Ihrig2019, ).
From these renormalization group functions, we need to obtain the
corresponding functions in a different renormalization scheme, which
we call CW scheme, using the multi-scales techniques reported in (Chishtie2008, ).
This paper is organized as follows: in Sec. II,
we present our model, together with the renormalization group functions
found in the literature. In Section III
we obtain the corresponding functions in the CW scheme and we use
them in Section IV
for the calculation of the effective potential using the RGE approach.
This effective potential is used in Section V
to study different aspects of the DSB in our model. Section VI
presents our conclusions and perspectives.
II the massless abelian
higgs model in the MS scheme and its corresponding
and functions
We start with the -component massless Abelian Higgs (AH) model
defined in -dimensional Euclidean space-time by the Lagrangian
(1)
where describes the
-component complex scalar field with quartic self-interaction
. This scalar is minimally coupled to an gauge field
via covariant derivative ,
being the analogous to the “electric charge”. The field strength
tensor is defined as
and the gauge-fixing Lagrangian is ,
being the gauge-fixing parameter. For the case of a single
complex scalar field, , and in three spatial dimensions, this
model is used to describe transitions on superconductors (Ginzburg1950, )
and liquid crystals (Halperin1974, ).
The renormalized Lagrangian of the massless AH model is
(2)
where
and is a mass scale introduced by the (dimensional)
regularization scheme (Collins1984, ; Peskin2018, ; Altland2010, ).
The wave-function renormalization constants and
relate the bare and the renormalized fields in the Lagrangian through
and .
Also, we obtain the relations between bare and renormalized coupling
constants as
(3)
(4)
(5)
It is interesting to note that the gauge-fixing parameter is also
renormalized (Collins1984, ) in this case, meaning it will
have a corresponding function.
In the MS scheme, the beta functions of the model are defined by
(6)
These functions were computed in the Ref. (Ihrig2019, ) up
to four loops, from which we quote the expressions below. The contributions
to and ,
can be cast as, for the gauge coupling constant,
(7)
where
(8a)
(8b)
(8c)
(8d)
for the scalar self-interaction,
(9)
where
(10a)
(10b)
(10c)
(10d)
and finally, for the gauge parameter,
(11)
where
(12a)
(12b)
(12c)
(12d)
The anomalous dimensions are defined through the relation
(13)
and, up to four loops, the contributions to read
as
(14)
where
(15a)
(15b)
(15c)
(15d)
The superscript present in the previous expressions denotes the aggregate
power of coupling constants. So, for instance,
means the terms in which contain exactly
two powers of coupling constants.
In this section we will use the renormalization group function obtained
in section II to calculate
the and function in the CW scheme. We know the
effective potential will involve terms with logarithms, of the general
form
(16)
(17)
where is associated to some coupling constant present in the
model, and is the classical value of one of the components
of , the one which is shifted as
in order to study the symmetry breaking (see details in Section V).
We can obtain the relation between the renormalization group function
in the CW scheme from the knowledge of the corresponding function
in the MS scheme. This procedure is not straightforward because we
have multiple coupling constants and then we have to use the multi-scale
procedure described in (Chishtie2008, ). In order to do that,
we start with Eq. (17) applying to our model, i.e.,
, then we get
(18a)
(18b)
(18c)
where with are different scales.
Notice the gauge parameter appears here as if it were a coupling constant,
being dimensionless and having its own function as
and . If we compare with Eq. (16)
we obtain the following relations:
(19)
The renormalization group function in the CW scheme can be defined
as
(20)
(21)
(22)
(23)
Now, if we use the relations Eq. (19)
in the last set of equations with the condition ,
we get the final relation between the renormalization group functions
in the CW computed from MS scheme,
(24a)
(24b)
(24c)
(24d)
Notice that the minus sign come from the definition
of the renormalization group function in the MS scheme (see Eqs. (6)
and (13)).
We obtain the CW RG functions through an order by order comparison
of the previous expression. For example, for the lowest order, we
have
(25a)
(25b)
(25c)
(25d)
where these relations are expected because at this
order the renormalization group function in the CW and MS are easily
related (see for example, the section 3 of the Ref. (Quinto2016, )
for more details).
For the next order, i.e, and ,
the relation between the two schemes is more complex,
(26a)
(26b)
(26c)
(26d)
For the order and
we get
(27a)
(27b)
(27c)
(27d)
And finally, for the order and
,
(28a)
(28b)
(28c)
(28d)
Now, with the four-loop renormalization group
functions written in the CW scheme, in the next section we will compute
the effective potential and study the CW mechanism by using the RGE
in the leading log approximation up to five loops.
IV effective potential
in the leading logs approximation
The main object we shall be interested in studying is the effective
potential. In order to compute this object, we consider a shift in
the -th component of in (1),
(29)
where
with and being two real scalar fields and
is a constant expectation value of scalar field, called
background field. This scalar field has the same properties
of . If we substitute (29) into (1)
we can find the effective potential in the classical approximation
(30)
It is easy to see that is the minimum of ,
so there is no spontaneous symmetry breaking at classical level. Our
aim is to compute loop corrections to
in order to understand if these corrections are capable to induce
a spontaneous symmetry breaking and the corresponding generation of
mass as given below
(31)
(32)
(33)
Notice that the gauge dependence on the mass of is a
consequence of a -gauge, .
Actually, in the spontaneously symmetry broken phase, the
degree of freedom is absorbed by the photon, which becomes massive.
As discussed in (Quinto2016, ; A.G.Quinto2021, ), the knowledge
of is sufficient for investigating
the dynamical breaking of gauge symmetry. We will be able to calculate
it by using an ansatz for the RGE, motivated by dimensional analysis,
together with the renormalization group functions for the model found
in section III. Specifically,
we shall use for the ansatz
(34)
where
(35)
and and are defined as power series in the coupling
constants , and is defined in (16).
This ansatz follows from the conformal invariance at the tree-level,
leading to the fact that we can have only one type of logarithm appearing
in the quantum corrections. Comparison with (30)
show us that
(36)
Now, we need to calculate the dependent pieces of ,
involving . For this, we need to use the RGE, in order to
obtain this equation we start with
(37)
where is independent of mass scale .
By deriving Eq. (37) with respect to ,
finally, inserting Eq. (34)
into (39), we obtain an alternative
form for the RGE,
(40)
were we used the notation .
Inserting the ansatz (35) in (40),
and separating the resulting expression by orders of , we obtain
a series of equations,
(41)
(42)
(43)
As we can see in (41) the function is
only dependent of the coupling , see Eq. (36),
for this reason the others beta functions were dropped.
We now consider that all functions appearing in Eq. (41)
are defined as series in powers of the couplings,
(44)
where the numbers in the superscripts denote the power of global coupling
constant of each term. Since all terms of the previous equation start
at order , except the first, we conclude
that , and obtain the relation
(45)
This last equation fixes the coefficients of
in terms of (known) coefficients of
and , in the following
form,
(46)
If we repeat the same procedure done in order to obtain
we can find the others ´s with the helps of ´s
presents in (44), as a results
we obtain
(47)
(48)
(49)
The coefficients to appearing in this equation
are defined in the appendix A.
Now looking at Eq. (42) expanded in powers
of the couplings,
(50)
one may conclude that starts at
order , obtaining the relation,
(51)
from witch the coefficients of the form of
are calculated from known coefficients of the beta function, anomalous
dimension, and . The end result is as follows,
(52)
If we repeat the same procedure, we can find the others ´s
with the helps of ´s and s presents
in (50), as a results we get,
(53)
(54)
The coefficients to are presented in the appendix B.
Finally, looking at Eq. (43) expanded in powers
of couplings,
(55)
one may conclude that starts at
order , leading to the relation,
(56)
from witch the coefficients of the form of
are calculated from the beta function, anomalous dimension, and .
The end result is as follows,
(57)
Then, we can find with the helps of ´s
and s presents in (55), as a results
we get
(58)
The coefficients to are presented in the appendix C.
These results will be used, in the next section, to study the modification
introduced by the leading logs summation in the DSB in our model.
V Dynamical symmetric breaking
In this section we will study the DSB in our model, for this, we will
use the results obtained in the previous section for the effective
potential up to five loops which was calculated using the renormalization
group equation, in the following form,
(59)
where is a finite renormalization constant and
is the regularized effective potential up to five loops. The constant
is fixed using the CW normalization condition,
(60)
Requiring that has a minimum
at means imposing that
(61)
which can be used to determine the value of as a function
of free parameters and . Upon explicit
calculation, Eq. (61) turns out to be a
polynomial equation in , and among its solutions, we look
for real and positive values for , and correspond to a minimum
of the potential. i.e.,
(62)
Using a program created in MATHEMATICA it was possible to verify for
which values of the free parameters , and
we obtain a sensible value for , which means that the mechanism
of DSB is operational, and the symmetry is indeed broken by radiative
corrections.It is well-known that the effective potential can be gauge-dependent (Jackiw1974, ).
There is a sophisticated method to deal with this problem developed
by Nielsen (Nielsen1975, ), which for sake of simplicity, will
be properly addressed in a future work.
In order to suggest the rich structure of DSB in the model, we will
present some results for Feynman-t’Hooft gauge, i.e., . We
considered and as free parameters, varying in the ranges
and . Considering these values,
the parameter space in which the DSB occurs was analyzed, and the
summary of our findings is pictured in Fig. 1.
Figure 1: In the left hand side we show a region plot corresponding
to the result of scanning for DSB for different values of the free
parameters and , with . In our model we find
three regions: in the yellow region we have three possible solutions
for , in the brown one we have only one solution and in
the blue region we do not have solution for , meaning DSB
is not operational. In the right hand side, we show a set of plots
explicitly showing the behavior of the solutions for as
a function of the parameter, for specific values of
and .
We notice the existence of three regions of solutions for ,
where the yellow region is characterized by the existence of three
different solutions for , which means three different non-symmetric
vacua for each value of the parameters within this region. The brown
region corresponding to the existence of a single solution, and the
blue region with the absence of any solution which breaks symmetry,
i.e., for the case when DSB does not happen in our model.
We can analyze the minimum of the effective potential, Eq (59),
for example for values of , and , we
found the three values of ,
(63)
and when these are used together we can see that the minimum occurs
as show in figure 2.
Figure 2: This graph of
vs. shows the potential minimum for the gauge
parameter, . The effective potential was evaluated using ,
and the values of , as can be seen in Eq. (63).
The red rectangle in the botton-left graph is shown in a different
scale in the top-right one.
VI Conclusion
In this paper we have studied the behavior of effective potential
in a massless Abelian Higgs (AH) model with -component complex
scalar field in dimensional space-time, specifically concerning
its classical vacua structure, obtaining hints of a very rich structure
of DSB, depending on the free parameters of the model (the gauge coupling
constant and the number of scalar fields). These results were obtained
by calculating the effective potential using the RGE equation, and
the and functions already
reported in the literature. After adapting these renormalization group
functions, which were calculated in the minimal subtraction scheme,
to a renormalization scheme adequate for our purposes, we shown how
the RGE can be used to calculate, order by order in the logarithm
and the coupling constants,
the effective potential.
Our results points to some interesting prospects, that we intent to
approach in future publications. First, to investigate whether further
higher-order terms can be incorporated in the effective potential
by using a different summation, which could be implemented as a symbolic
package for MATHEMATICA, thus further improving our calculation (as
it was done in Quinto2016 ). Second, a full study of the
Nielsen identity in our model would be important to factor out the
possible gauge dependence of the results. Finally, we expect to apply
this formalism for other models with application in condensed matter
and particle physics in other to study their DBS properties.
Acknowledgements.
This work was supported by Fondo Nacional de Financiamiento
para la Ciencia, la Tecnología y la Innovación "Francisco
José de Caldas", Minciencias Grand No. 848-2019(AGQ), and Conselho Nacional de Desenvolvimento Científico e Tecnológico
(CNPq) Grant No. 305967/2020-7 (AFF).
Appendix A All coefficients of ´s
In this appendix we show the values of the coefficients associated
with the functions as a function
of the Riemann Zeta functions, . In our case we have
and presents in our results. So,
we can fix these with , known as Apéry’s constant,
and .
(64)
Appendix B All coefficients for ´s
In this appendix we show the values of the coefficients associated
with the functions ,
(65)
Appendix C All coefficients for ´s
In this appendix we show the values of the coefficients associated
with the function ,
(66)
References
[1]
J. Zinn-Justin.
Quantum field theory and critical phenomena.
Clarendon, Oxford, 1996.
[2]
M. E. Peskin.
An introduction to quantum field theory.
CRCPress, Boca Raton, FL, 2018.
[3]
A. Altland and B. D. Simons.
Condensed matter field theory.
Cambridge University Press, Cambridge, 2010.
[4]
I. Herbut.
A modern approach to critical phenomena.
Cambridge University Press, Cambridge, 2007.
[5]
Coleman Sidney and Weinberg Erick.
Radiative corrections as the origin of spontaneous symmetry breaking.
Phys. Rev. D, 7:1888–1910, Mar 1973.
[6]
V. Elias, R. B. Mann, D. G. C. McKeon, and T. G. Steele.
Radiative electroweak symmetry breaking revisited.
Phys. Rev. Lett., 91:251601, Dec 2003.
[7]
F. A. Chishtie, V. Elias, R. B. Mann, D. G. C. McKeon, and T. G. Steele.
Stability of subsequent-to-leading-logarithm corrections to the
effective potential for radiative electroweak symmetry breaking.
743:104–132.
[8]
F. A. Chishtie, T. Hanif, J. Jia, R. B. Mann, D. G. C. McKeon, T. N. Sherry,
and T. G. Steele.
Can the renormalization group improved effective potential be used to
estimate the higgs mass in the conformal limit of the standard model?
Phys. Rev. D, 83:105009, May 2011.
[9]
T. G. Steele and Zhi-Wei Wang.
Is radiative electroweak symmetry breaking consistent with a 125 gev
higgs mass?
Phys. Rev. Lett., 110:151601, Apr 2013.
[10]
Boyang Liu, Xiao-Lu Yu, and Wu-Ming Liu.
Renormalization-group analysis of -orbital bose-einstein
condensates in a square optical lattice.
Phys. Rev. A, 88:063605, Dec 2013.
[11]
Shun Uchino.
Coleman-weinberg mechanism in spinor bose-einstein condensates.
EPL (Europhysics Letters), 107(3):30004, aug 2014.
[12]
I. S. Burmistrov, Y. Gefen, D. S. Shapiro, and A. Shnirman.
Mesoscopic stoner instability in open quantum dots: Suppression of
coleman-weinberg mechanism by electron tunneling.
Phys. Rev. Lett., 124:196801, May 2020.
[13]
Alan Chodos and Hisakazu Minakata.
The gross-neveu model as an effective theory for polyacetylene.
Phys. Lett. A, 191:39, 1994.
[14]
R. Vega Monroy A. G. Quinto and A. F. Ferrari.
Renormalization group improvement of the effective potential in a
(1+1) dimensional gross-neveu model.
arXiv:2108.04079, 2021.
[15]
M.R. Ahmady, V. Elias, D.G.C. McKeon, A. Squires, and T.G. Steele.
Renormalization-group improvement of effective actions beyond
summation of leading logarithms.
Nuclear Physics B, 655(3):221–249, 2003.
[16]
V. Elias, R. B. Mann, D. G. C. McKeon, and T. G. Steele.
Higher order stability of a radiatively induced 220 gev higgs mass.
Phys. Rev. D, 72(3):037902, Aug 2005.
[17]
Huan Souza, L. Ibiapina Bevilaqua, and A. C. Lehum.
Renormalization group improvement of the effective potential in six
dimensions.
Phys. Rev. D, 102(4):045004, 2020.
[18]
A.G. Quinto, A.F. Ferrari, and A.C. Lehum.
Renormalization group improvement and dynamical breaking of symmetry
in a supersymmetric chern-simons-matter model.
Nuclear Physics B, 907:664–677, 2016.
[19]
A.G. Dias, J.D. Gomez, A.A. Natale, A.G. Quinto, and A.F. Ferrari.
Non-perturbative fixed points and renormalization group improved
effective potential.
Physics Letters B, 739:8–12, 2014.
[20]
F. A. CHISHTIE, V. ELIAS, R. B. MANN, D. G. C. MCKEON, and T. G. STEELE.
On the standard approach to renormalization group improvement.
International Journal of Modern Physics E, 16(06):1681–1685,
2007.
[21]
Bernhard Ihrig, Nikolai Zerf, Peter Marquard, Igor F. Herbut, and Michael M.
Scherer.
Abelian higgs model at four loops, fixed-point collision, and
deconfined criticality.
Phys. Rev. B, 100:134507, Oct 2019.
[22]
F. A. Chishtie, T. Hanif, D. G. C. McKeon, and T. G. Steele.
Unique determination of the effective potential in terms of
renormalization group functions.
Phys. Rev. D, 77:065007, Mar 2008.
[23]
V. L. Ginzburg and L. D. Landau.
On the theory of superconductivity.
Zh. Eksp. Teor. Fiz., 20(1064), 1950.
[24]
B. I. Halperin, T. C. Lubensky, and Shang-keng Ma.
First-order phase transitions in superconductors and smectic-
liquid crystals.
Phys. Rev. Lett., 32:292–295, Feb 1974.
[25]
John C. Collins.
Renormalization.
Cambridge University Press, 1984.
[26]
Alexander Altland and Ben Simons.
Condensed Matter Field Theory.
Cambridge, 2010.
[27]
A. G. Dias and A. F. Ferrari.
Renormalization group and conformal symmetry breaking in the
chern-simons theory coupled to matter.
Phys. Rev. D, 82:085006, Oct 2010.
[28]
R. Jackiw.
Functional evaluation of the effective potential.
Phys. Rev. D, 9:1686–1701, Mar 1974.
[29]
N. K. Nielsen.
On the gauge dependence of spontaneous symmetry breaking in gauge
theories.
Nuclear Physics B, 101(1):173–188, 1975.