Propagador Vestido del Gluón con una Ventana

Conforme en el Infrarojo y Confinamiento de

Quarks.

Marco Antonio Bedolla Hernández

Febrero de 2013
Agradecimientos

Antes que nada, quiero agradecer a mis padres, por todo el apoyo y cariño
que me han brindado siempre, por la confianza que han depositado en mi, los
cuales han resultado imprescindibles para poder llevar con éxito esta tarea.

A mi familia por el apoyo, paciencia y comprensión dados. A mi tía Gracia

por ver que nunca me faltara nada.

Al Dr. Alfredo Raya, por sus sugerencias y apoyo constante durante todo este
tiempo, además que ha sido una fuente importante de conocimientos y trabajo
académico, ha servido como una excelente guía no solo para la realización de
esta tesis, sino también para mi formación académica.

A los Drs. Mariana Kirchbach y Adnan Bashir, por su disposición para revi-
sar y aprobar la tesis a pesar de tener un corto margen de tiempo. Al Dr. Ulises
por su disponibilidad para ayudarnos a seguir con los trámites a pesar de las
inconveniencias que se presentaron.

A mis compañeros y amigos: Topo, Pancho, Nucita, Calambres, Dante, Abraham,

Khépani, Miguel, Mauricio, Saúl y Marisol. Mención especial para los últimos
dos, el primero por ayudarme con la portada de la tesis y con los trámites que
tan bien me caen. A Marisol por hacer que los días de mayor carga de trabajo
fueran mucho más llevaderos y además por sus consejos finales. A aquellos que
no mencioné, no se agüiten, ya saben los quiero.

Al Instituto de Física y Matemáticas y a la facultad de Cs. Físico-Matemáticas

de la UMSNH y a todo su cuerpo académico y administrativo por haber sido
de gran importancia para mi formación, así como las facilidades que me dieron
para reaizar este trabajo y por la hospitalidad que me han brindado.

Al CONACYT por el apoyo brindado a través de la beca que me otorgó para

poder estudiar esta Maestría en Ciencias.

Para finalizar me gustaría agradecer a Microsoft, a la comunidad de Ubun-

tu, a Wolfram, a Jimmy Wales y Larry Sanger, a Nintendo, a Epic Games, a
Valve y en especial a Bethesda por haber creado el mejor juego de la historia.

3
4
Índice general

1. Introducción 7

2. Modelos Potenciales 11
2.1. Potencial de Cornell . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2. Solución de la ecuación de Schrödinger con el potencial de Cornell 12
2.2.1. Perturbaciones a la Energía . . . . . . . . . . . . . . . . . 13
2.3. Potencial de Rosen-Morse Trigonométrico . . . . . . . . . . . . . 15
2.3.1. Solución de la ecuación de Schrödinger con el potencial
RMT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.4. Propagador Instantáneo . . . . . . . . . . . . . . . . . . . . . . . 21

3. Generación Dinámica de Masas y Confinamiento 27

3.1. Ecuaciones de Schwinger-Dyson . . . . . . . . . . . . . . . . . . . 27
3.1.1. Ecuación de GAP en QCD . . . . . . . . . . . . . . . . . 31
3.1.2. Otros Observables . . . . . . . . . . . . . . . . . . . . . . 33
3.1.3. Confinamiento en las Ecuaciones de Schwinger-Dyson . . 34
3.2. Modelo de Maris-Tandy para QCD . . . . . . . . . . . . . . . . . 34

4. Modelo Para el Propagador Vestido del Gluón 37

5. Conclusiones y Consideraciones Finales 47

A. Lagrangiano de QCD 49

B. Programa usado para resolver las ESD 51

5
6 ÍNDICE GENERAL
Capítulo 1

Introducción

Entender el espectro de los mesones y las resonancias bariónicas es crucial

para identificar el mecanismo del confinamiento y la dinámica de los quarks y
los gluones en la región no perturbativa de la cromodinámica cuántica (QCD,
por sus siglas en inglés, Quantum Chromodynamics). Para bariones, el explicar
los estados con espín alto es un verdadero reto en la espectroscopía de quarks
[1, 2]. En ese sentido, el número de resonancias observadas hasta el momento
[3] es considerablemente menor que aquellas predichas en modelos basados en
un espacio de Hilbert de seis grados de libertad de sabor y la simetría más baja
de las funciones de onda radiales [4]. Los candidatos obvios para reducir el nú-
mero de grados de libertad son los modelos quark-diquark (q-qq) [5] basados en
un diquark con momento angular acotado. El sistema q-qq ha sido promovido
en este contexto por varios autores [6, 7]. Se acostumbra imponer una simetría
en la función de onda mayor que SO(3)L , predominando la simetría conforme
como elección favorita. Esta simetría es –en buena medida– una simetría global
del Lagrangiano de QCD en el sector de sabores ligeros, razón por la que es
de esperarse que los datos de espectroscospía para el Nucleón y la ∆ puedan
ser examinados en el esquema de la correspondencia AdS/CF T . Los trabajos
pioneros en está conexión han construido exitosamente funciones de onda espa-
ciales para mesones y bariones ligeros [8]. Trabajos posteriores extendieron esta
idea [9, 10, 11]. En [12], una simetría conforme en el sector de bariones ligeros
dentro de la correspondencia AdS5 /CF T4 aparece en la degeneración entre pa-
rejas de paridad que es sistemáticamente observada en el espectro del Nucleón
y ∆. En ese trabajo, la simetría conforme se implementa en el Hamiltoniano de
quarks definido en un volumen finito. Existe una variedad de geometrías apro-
piadas para tomar en cuenta la simetría conforme en Hamiltonianos de este tipo,
siendo la hiperesfera S 3 una elección notable [13, 14, 15]. Se ha mostrado que
una geometría que contenga a S 3 , tal como S 1 ⊗ S 3 , puede obtenerse partiendo
directamente de AdS5 [14, 16, 17]. Este hecho fue usado en [18], obteniendo
el espectro de un sistema q-qq localizado directamente sobre el cono AdS5 en
su forma conformamente compactificada hacia un espacio-tiempo de Minkowski
compactificado conformemente, R1 ⊗ S 3 , sobre el cono AdS5 . Se observa que

7
8 CAPÍTULO 1. INTRODUCCIÓN

el movimiento cuántico del quark libre (por ahora considerado sin spin) sobre
S 3 y considerando al diquark ubicado en el polo Norte, se describe en términos
de la bien conocida ecuación escalar conforme, mientras que el potencial está-
tico dipolar entre q-qq se considera como el potencial “curvo” de Coulomb. El
potencial curvo de Coulomb fue estudiado por primera vez por E. Schrödinger
en [19]. Corresponde a la función cotangente (cot), que resulta ser armónica
en S 3 al resolver la ecuación de Laplace-Beltrami sobre S 3 . Debido a que el
grupo de isometrías de S 3 es SO(4), la función cot respeta esta simetría [20].
Además, el potencial cot encierra los rasgos principales de la dinámica en QCD,
como se muestra en [21], ya que se interpola entre el régimen perturbativo del
intercambio de un gluón (inverso con la distancia), pasando a través el régimen
no perturbativo de flujo de gluones (crecimiento lineal) y el crecimiento infi-
nito que confina a los quarks (pared infinita). El potencial cot es familiar en
estudios de de mecánica cuántica supersimétrica [22, 23, 24], donde es conocido
como potencial de Rosen-Morse, el cual, siendo de rango finito, ha capturado
el interés en modelar confinamiento espacial ya sea para electrones en puntos
cuánticos [25] o quarks en hadrones [18]. Conocer el potencial y las funciones
de onda en el espacio de momentos demanda el diseño de un espacio fase que
considere al potencial cot como una función del segmento de arco en S 3 . En [26],
la proyección sobre el disco ecuatorial del movimiento sobre el hemisferio Norte
del quark, perturbada por el potencial cot es considerada con el fin de efecutar
su transformada de Fourier hacia el espacio de momentos y así encontrar la for-
ma del potencial en este espacio. Por un lado, se ha encontrado que la función
Sinc (sin2 (q)/q 2 ) es la transformada de Fourier de la función cot, o la función
de Struve (H1 (q)/q 2 ) como el potencial efectivo de la transformada de Fourier
de la ecuación de movimiento completa. A estas funciones se les puede dar la
interpretación de propagadores instantáneos efectivos del gluón. Sabiendo que
el confinamiento de quark ligeros es un efecto de la teoría cuántica de campos
(TCC), en esta tesis se explora el impacto de un propagador vestido del gluón
que muestra confinamiento espacial, en el escenario de confinamiento de quarks
en un esquema de la TCC. Se ha seleccionado el marco de las Ecuaciones de
Schwinger-Dyson (ESD) para este propósito. Las ESD constituyen una herra-
mienta poderosa para la física de hadrones. Este esquema provee una explicación
satisfactoria de la relación entre el Rompimiento Dinámico de la Simetría Quiral
(RDSQ) y el Confinamiento, ambos fenómenos emergentes de QCD. Además,
el crecimiento infrarrojo de la función de masa de los quarks se entiende como
un efecto de la nube de gluones que viste al quark a bajos momentos [27]. Esos
hechos pueden ser explorados a través de las ESD. En esta tesis se explora el
RDSQ y Confinamiento con un modelo efectivo para la ecuación de GAP (hue-
co en el espectro). Bajo suposiciones minimalistas, es decir, usando el vértice
desnudo y un modelo para la función de revestimiento del gluón, se conduce
un análisis del propagador del quark en la norma de Landau. Para el gluón,
se emplea una función de revestimiento inspirada en el propagador instantáneo
del potencial de Coulomb curvo. Las predicciones son contrastadas con obser-
vables hadrónicos relevantes derivados de nuestro modelo con valores aceptados
fenomenológicamente, y se observa buena coincidencia dentro de la validez de
 9

nuestras suposiciones. La tesis está organizada de la siguiente forma: en el si-

guiente capítulo se resuelve la ecuación de Schrödinger para el potencial de
Cornell y el potencial Rosen-Morse trigonométrico; y luego se construye nuestra
propuesta para el propagador efectivo del gluón. En el capítulo 3 se da una breve
introducción a las ESD, se presentan los observables que pueden ser calculados;
la prueba confinamiento, el condensado quiral y la constante de decaimiento del
pión. Luego se presentan resultados para el modelo de Maris-Tandy (MT) [29],
que en los últimos quince años ha guiado las predicciones de las ESD para varios
observables hadrónicos. En el capítulo IV se pone a prueba nuestro propagador
efectivo del gluón y se compara con los resultados obtenidos en el modelo MT.
Finalmente, en el último capítulo se presentan las conclusiones generales sobre
este trabajo.
10 CAPÍTULO 1. INTRODUCCIÓN
Capítulo 2

Modelos Potenciales

La interacción fuerte entre quarks, los constituyentes fundamentales de los

hadrones, está regida por la cromodinámica cuántica (QCD), la cual es una
teoría de norma no Abeliana con los gluones como bosones de norma (en el
Apéndice A se presenta el Lagrangiano de QCD). Como consecuencia del ca-
rácter no Abeliano de QCD, las interacciones entre quarks van del intercambio
de uno a muchos gluones y autointeracciones del gluón, siendo responsable esto
último del confinamiento de quarks, que a bajas energías, permanecen confina-
dos, pero a altas energías y momentos se comportan casi como partículas casi
libres. Las ecuaciones de campo en QCD son no lineales y bastante complica-
das, debido a los procesos de autointeracción del gluón, y para poder resolverse,
se requiere el uso de técnicas altamente sofisticadas, como la discretización del
espacio-tiempo, Lattice QCD1 . Los cálculos en Lattice se han establecido como
una herramienta de confianza a la hora de hacer cálculos no perturbativos en
QCD. Uno de los resultados obtenidos mediante esta técnica es la forma cuali-
tativa para el potencial de confinamiento de los quarks, tomando como modelo
una interacción estática quark-diquark, que consiste en un potencial con un tér-
mino tipo Coulomb (inverso a la distancia) más un término lineal: Potencial de
Cornell [28]. Conociendo la forma del potencial, es posible resolver el problema
mecánico cuántico para obtener el espectro de masa y las funciones de onda para
la interacción entre quarks y diquarks dentro de los nucleones. A continuación,
se describen algunos ejemplos de dichos potenciales siendo resueltos usando la
ecuación de Schrödinger.

2.1. Potencial de Cornell

El potencial de Cornell tiene la forma

A
V (r) = − + Br + O(r2 ), (2.1)
r
1
Para un estudio detallado ver [30].

11
12 CAPÍTULO 2. MODELOS POTENCIALES

donde A y B son parámetros ajustables para lograr una descripción adecuada de

la dinámica de los quarks. El término inverso a r (término tipo Coulomb) puede
interpretarse como el intercambio de un gluón en el modelo de quarks dentro
del régimen perturbativo. Por otro lado, el término lineal en r corresponde a la
región no perturbativa, que se interpreta como un flujo de gluones coherentes
entre quarks, cuya intensidad va aumentando conforme los quarks se alejan,
similar al comportamiento de un cuerda que se estira.
El problema de este modelo consiste en que no existe una cota máxima para
la distancia entre quarks, lo que es un indicador de que sólo es válido para valores
pequeños de r. Esto sugiere que se debe encontrar un modelo más preciso para
la interacción de un quark-diquark dentro del nucleón. Sin embargo, se pueden
obtener resultados importantes usando esta primera aproximación.

2.2. Solución de la ecuación de Schrödinger con

el potencial de Cornell
La ecuación de Schrödinger con el potencial de Cornell ha sido ampliamente
estudiada en la mecánica cuántica. El método estándar –que será bosquejado
a continuación– consiste en tomar la parte lineal del potencial como una per-
turbación al potencial del tipo Coulomb y utilizar la teoría de perturbaciones
ordinaria. Considérese la ecuación de Schrödinger en tres dimensiones con un
potencial radial,
~2 2
 
− ∇ + V (r) Ψ(~r) = EΨ(~r), (2.2)
2µ
donde Ψ(~r) es la función de onda de la partícula, E son los eigenvalores de
energía y µ es la masa reducida del sistema quark-diquark. Para nuestro objetivo,
se necesita la solución de la parte radial de la ecuación (2.2) con el potencial
V (r) = −A/r en coordenadas esféricas [31]. La solución de este problema es
bien conocida, dado que la ecuación se reduce a la del problema del átomo de
hidrógeno [32]. La única diferencia en este problema es que la constante e2 /4π0
se cambia por A. La ecuación radial es
2µr2
   
d 2 dR A
r − 2 − − E R = l(l + 1)R. (2.3)
dr dr ~ r
Esta ecuación se simplifica al hacer el cambio de variable
u(r) = rR(r), (2.4)
obteniendo así
~2 d 2 u ~2 l(l + 1)
 
A
− + − + u = Eur. (2.5)
2µ dr2 r 2µ r2
√
Tomando k = −2µE/~, el cual es real para E < 0, se obtiene
1 d2 u
 
2µA 1 l(l + 1)
= 1 − + u, (2.6)
k 2 dr2 ~2 k kr (kr)2
2.2. SOLUCIÓN DE LA ECUACIÓN DE SCHRÖDINGER CON EL POTENCIAL DE CORNELL13

lo que sugiere los cambios de variables

2µA
ρ = kr, ρ0 = , (2.7)
~2 k
consiguiendo así una ecuación diferencial con dependencia en ρ,

d2 u
 
ρ0 l(l + 1)
= 1 − + u. (2.8)
dρ2 ρ ρ2

Nótese que la ecuación (2.8) ahora sí es idéntica a la ecuación de eigenvalores

para el átomo de hidrógeno [33]. Con este resultado, se puede establecer que el
número cuántico principal es ρ0 = 2n, donde n es entero positivo.
Por otro lado, las energías permitidas son de la forma
2
~2 k 2 ~2 2µA 2µA2

E=− =− = − . (2.9)
2µ 2µ ~2 ρ0 ~2 ρ20

Combinando los dos resultados previos, se obtiene que las energías permitidas
son  
µ 2 1 E1
En = − A = 2, n = 1, 2, 3, . . . , (2.10)
2~2 n2 n
donde
µ 2
E1 = A , (2.11)
2~2
es la energía del estado base. Finalmente, combinando los resultados obtenidos
para ρ0 , se tiene que

2µA µA 1
2n = 2
⇒k= = . (2.12)
~ k ~n an
Aquí,
~2
a= , (2.13)
µA
ha sido identificado como el Radio de Bohr. El parámetro A puede ajustarse
para calcular la energía del estado base no perturbado de modelos que se apro-
ximen al potencial de Cornell. Por ejemplo, para el Charmonio se tiene A ≈ 0.48
[34].

2.2.1. Perturbaciones a la Energía

Con el fin de obtener una solución problema para el potencial de Cornell, se
usa la teoría de perturbaciones usual para la mecánica cuántica, como se ve en
[35]. Considerando el Hamiltoniano

H = H0 + λH 0 , (2.14)
14 CAPÍTULO 2. MODELOS POTENCIALES

donde H 0 es la perturbación, la corrección perturbativa a primer orden para la

energía del estado base viene dada por
E01 = hψ0 |H 0 |ψ0 i, (2.15)
siendo ψ0 la función de onda del estado base. Recordando que dicha función y
el Hamiltoniano de perturbación están dados por [32]
1
ψ100 (~r) = √ e−r/a , (2.16)
πa3
r
H0 = B , (2.17)
λ
siendo B es un parámetro pequeño, se tiene que la corrección de la energía a
primer orden es
4B ∞
Z
E01 = hψ0 |Br|ψ0 i = 3 drr3 e2r/a
a 0
3
= Ba. (2.18)
2
Entonces, la energía corregida del estado base es
µ 2 3
E1 = − A + Ba, (2.19)
2~2 2
o, sustituyendo el valor de a,
µ 2 3 ~2
E1 = − A + B, (2.20)
2~2 2 µA
el cual coincide en la forma con el resultado obtenido mediante un procedimien-
to algebraico [36].

Es importante tomar en cuenta el hecho de que las soluciones obtenidas

anteriormente para el potencial de Cornell –usando teoría de perturbaciones
alrededor de un potencial de Coulomb– simplemente son una primera apro-
ximación para la interacción quark-diquark obtenida mediante Lattice. Como
se mencionó previamente, este potencial tiene rango infinito, lo cual es indicio
de que términos importantes en la expansión del potencial verdadero han sido
omitidos. Se ha visto que la forma del potencial requerida es del tipo − cot,
que puede aproximarse con el potencial de Cornell para valores pequeños de r,
pero además tiene la ventaja de tener un rango finito, el cual cumple con los
requisitos para un potencial entre quarks como se ve en [37]; el potencial cot
interpola el regimen perturbativo de intercambio de un gluón (término inverso
con la distancia), pasa a través del régimen no perturbativo de interacciones
(crecimiento lineal) y el crecimiento infinito que confina a los quarks (pared
infinita). El potencial cot es familiar dentro de los estudios de Mecánica Cuán-
tica Supersimétrica, donde se encuentra como un caso particular del potencial
Rosen-Morse Trigonométrico (RMT) [38, 39], el cual se analizará brevemente
en la siguiente sección.
2.3. POTENCIAL DE ROSEN-MORSE TRIGONOMÉTRICO 15

2.3. Potencial de Rosen-Morse Trigonométrico

El potencial de Rosen-Morse Trigonométrico (RMT) es

V (r) = −2b cot r + a(a + 1) csc2 r, (2.21)

el cual ha sido estudiado en la literatura y su solución ha sido reciente presenta-

da en [40] en términos de polinomios reales de Romanovski, en lugar de términos
de polinomios de Jacobi complejos. Para esta tesis, se resuelve nuevamente la
ecuación de Schödinger con dicho potencial, pero además se repite todo el pro-
cedimiento algebraico con el fin de corroborar que la solución presentada en [40].

El potencial RMT contiene el término − cot, cuya expansión en serie de

Taylor es de la forma

1 r r3
cot r ≈ − − + ··· , (2.22)
r 3 45
donde se puede ver que el potencial de Cornell está incluido en los primeros
dos términos de cot. La Figura 2.1 muestra la semejanza entre potenciales. Se

Figura 2.1: Comparación entre el potencial de Cornell y − cot.

resolverá la ecuación de Schrödinger para el potencial RMT usando el método

de Sturm-Liouville [41, 42] y enseguida se tomará el límite a → 0.

2.3.1. Solución de la ecuación de Schrödinger con el po-

tencial RMT
La parte radial de la ecuación de Schrödinger en tres dimensiones para el
potencial RMT es [32]

d2
 
l(l + 1) 2
− R(r) + − 2b cot r + a(a + 1) csc r −  R(r) = 0. (2.23)
dr2 r2
16 CAPÍTULO 2. MODELOS POTENCIALES

Nuestro enfoque será únicamente para el estado base, para lo cual se toma l = 0.
Así, la ecuación a resolver es
d2
− R(r) + (−2b cot r + a(a + 1) csc2 r − )R(r) = 0. (2.24)
dr2
Se propone una solución de la forma
R(r) = e−αr F (r), (2.25)
que al sustituir en la ecuación (2.24) y usando el hecho de que
R00 (r) = α2 e−αr F (r) − 2αe−αr F 0 (r) + e−αr F 00 (r), (2.26)
nos permite obtener
−α2 e−αr F (r) + 2αe−αr F 0 (r) − e−αr F 00 (r)
+ (−2b cot r + a(a + 1) csc2 −)e−αr F (r) = 0. (2.27)
Dividiendo por e−αr y factorizando los términos que incluyen derivadas de F ,
F 00 (r) − 2αF 0 (r) + (2b cot r − a(a + 1) csc2 r +  + α2 )F (r) = 0. (2.28)
Con el cambio de variable x = cot r, se obtienen los diferenciales
dx = − csc2 rdr = −(1 + x2 )dr,
d2
 
2 d 2 d
= −(1 + x ) − (1 + x )
dr2 dx dx
2
d d
= (1 + x2 )2 2 + 2x(1 + x2 ) ,
dx dx
que al sustituirlos en la ecuación (2.28), conducen a la ecuación diferencial final

(1 + x2 )2 f 00 (x) + 2(1 + x2 )(α + x)f 0 (x) + (2bx − a(a + 1)(1 + x2 ) +  + α2 )f (x) = 0,

(2.29)
donde se ha hecho el cambio f (x) = F (cot−1 r). Debido a la importancia para
conocer la solución de está ecuación, la mayoría de los pasos algebraicos son
mostrados en detalle. Primero, se propone una solución de la forma
1−β
f (x) = (1 + x2 )− 2 C(x), (2.30)
β es una constante que será calculada posteriormente. Las primeras derivadas
de la solución f (x) son:
(1 − β) 3−β 1−β
f 0 (x)
= −2x (1 + x2 )− 2 C(x) + (1 + x2 )− 2 C 0 (x) (2.31)
 2
00 (1 − β) (3 − β) 2 5−β
f (x) = 4 x (1 + x2 )− 2
2 2

(1 − β) 3−β
− −2 (1 + x2 )− 2 C(x)
2
3−β 1−β
− 2x(1 − β)(1 + x2 )− 2 C 0 (x) + (1 + x2 )− 2 C 00 (x), (2.32)
2.3. POTENCIAL DE ROSEN-MORSE TRIGONOMÉTRICO 17

que al sustituir en la ecuación (2.29), nos permite escribir


5−β
(1 + x ) (1 − β)(3 − β)x2 (1 + x2 )− 2 C(x)
2 2

3−β 3−β
− (1 − β)(1 + x2 )− 2C(x) − 2x(1 − β)(1 + x2 )− 2 C 0 (x)

1−β
+ (1 + x2 )− 2 C 00 (x)

3−β
+ 2(1 + x2 )(α + x) x(1 − β)(1 + x2 )− 2 C(x)

1−β
+ (1 + x2 )− 2 C 0 (x)
1−β
+ [2bx − a(a + 1)(1 + x2 ) +  + α2 ](1 + x2 )− 2 C(x) = 0. (2.33)

Con el fin de eliminar la dependencia de β en las potencias, la expresión anterior

es multiplicada por (1 + x2 )−(1+β)/2 , para así obtener
h
(1 + x2 )2 (1 − β)(3 − β)x2 (1 + x2 )−3 C(x) − (1 − β)(1 + x2 )−2 C(x)
i
− 2x(1 − β)(1 + x2 )−2 C 0 (x) + (1 + x2 )−1 C 00 (x)
h i
+ 2(1 + x2 )(α + x) − x(1 − β)(1 + x2 )−2 C(x) + (1 + x2 )−1 C 0 (x)
h i
+ 2bx − a(a + 1)(1 + x2 ) +  + α2 (1 + x2 )−1 C(x) = 0. (2.34)

Factorizando las derivadas de C(x), se ve que

0 = (1 + x2 )C 00 (x) + 2(α + βx)C 0 (x)

2x[−α(1 − β) + b] + α2 − (1 − β)2 + 
h 
− β(1 − β) + a(a + 1) − C(x).
1 + x2
(2.35)

Para obtener una ecuación más simple, se elimina la dependencia de x en el

término que multiplica a C(x). Se debe cumplir que

−α(1 − β) + b = 0, α2 − (1 − β)2 +  = 0, (2.36)

lo que conduce a una ecuación más simple,

(1 + x2 )C 00 (x) + 2(α + βx)C 0 (x) − [β(1 − β) + a(a + 1)]C(x) = 0. (2.37)

De acuerdo con la teoría de Sturm-Liouville [43], cualquier ecuación diferencial

puede escribirse en la forma

d2 Cm (x)
 
1 d(s(x)ω(x)) dCm (x)
S(x) + + λm C(x) = 0, (2.38)
dx2 ω(x) dx dx
18 CAPÍTULO 2. MODELOS POTENCIALES

donde
d2 s
 
dC1 (x) 1
λm = −m k1 + (m − 1) 2 , (2.39)
dx 2 dx
m
 
1 d
Cm (x) = ω(x)sm (x) . (2.40)
Km ω(x) dxm
Mediante un simple análisis a la ecuación (2.37), inmediatamente se identifican
los términos de la fórmula de Rodrigues
S(x) = (1 + x2 ), (2.41)
 
dC1
λm = −β(1 − β) − a(a + 1) = −m k + (m − 1) . (2.42)
dx
Por otro lado, usando el hecho de que
1 ds(x)ω(x)
C1 (x) = , (2.43)
kω(x) dx
se encuentra el valor de la constante C1 de la relación
kC1 (x) = 2(α + βx). (2.44)
De esta forma, se obtiene una ecuación más simple para β,
−β(1 − β) − a(a + 1) = −m(2β + m − 1). (2.45)
En resumen, las ecuaciones para calcular el valor de β son
b
β = 1− , (2.46)
α
 2
2 b
 = −α + , (2.47)
α
0 = β 2 + (2m − 1)β − a(a + 1) + m(m − 1). (2.48)
Sustituyendo el valor de β en la ecuación (2.48), se llega a la ecuación de segundo
orden
b2 b
− (2m + 1) + m(m + 1) − a(a + 1), (2.49)
α2 α
cuyas soluciones son
b
= m − a, (2.50)
α1
b
= m + 1 + a. (2.51)
α2
Tomando en cuenta la solución de la ecuación (2.51) y considerando m + 1 = n,
se obtiene
β = 1 − (n + a), (2.52)
2
b
 = (n + a)2 − . (2.53)
(n + a)2
2.3. POTENCIAL DE ROSEN-MORSE TRIGONOMÉTRICO 19

Finalmente, la función de peso ω(x) puede obtenerse a partir de

1 d(s(x)ω(x))
= 2(α + (1 − (n + a))x), (2.54)
ω(x) dx

de donde se obtiene

(1 + x2 )ω 0 (x) = 2(α − (n + a))xω(x), (2.55)

cuya solución es

−1
ωn (x) = (1 + x2 )−(n+a) e−2α cot x
. (2.56)

Los polinomios que resuelven la ecuación diferencial que se obtienen de la ecua-

ción (2.40) corresponden a los polinomios de Ruth-Romanovski [44, 45], los
cuales fueron redescubiertos en [40]. Si se reescribe dicha ecuación en términos
de n, ésta es

1 dn−1
Cn(a,b) (x) = (ωn (x)s(x)n−1 ), (2.57)
Kn ω(x) dxn−1

siendo Kn las constantes de normalización. Los primeros polinomios obtenidos

en esta forma son

(a,b) 1
C1 = , (2.58)
K1
 
(a,b) 2 b
C2 = − (1 + a)x , (2.59)
K2 2 + a
"
2b2

(a,b) 2 2(2a + 3)b
C3 = 2
− (1 + a) − x+
K3 (3 + a) 3+a
#
+ (1 + a)(2a + 3)x2 , (2.60)
20 CAPÍTULO 2. MODELOS POTENCIALES
"
(a,b) 4 2b3
C4 =
K4 (4 + a)3
2b2
  
(3a + 4)b
− − 3(2 + a) − (1 + a) x
4+a (4 + a)2
#
b 2 3
+3(a + 2)(2a + 3) x − (a + 1)(2a + 3)(2 + a)x , (2.61)
4+a
"
4b4 4b2

(a,b) 4
C5 = − (3a + 5) + 3(a + 2)(a + 1)
K5 (a + 5)4 (a + 5)2
2b3
 
b
−4(2a + 5) − (3a + 4) x
(a + 5)3 a+5
2b2
 
+6(a + 2)(2a + 5) − (a + 1) x2
(a + 5)2
b
−4(2a + 3)(a + 2)(2a + 5) x3
a+5
#
+(a + 1)(2a + 3)(a + 2)(2a + 5)x4 . (2.62)

La ortogonalidad de las funciones previas, se obtiene con la relación

Z ∞
dx (a,b)
(ωn (x))1/2 Cn(a,b) (ωn0 (x))1/2 Cn0 = δnn0 , (2.63)
−∞ s(x)
mediante la cual se pueden calcular las constantes de renormalización Kn 2 . Para
este caso, la ortogonalidad de las soluciones Rn (z) viene dada por
Z π
dzRn (z)(Rn0 (z))∗ = δnn0 . (2.64)
0

La ecuación (2.63) muestra que estas soluciones tienen propiedades de orto-

normalidad bien definidas sobre el eje real, lo cual las convierte en funciones
adecuadas para resolver problemas de mecánica cuántica.
Tomando el límite a → 0 se obtiene el espectro de energías para el potencial
v(r) = −2b cot r
b2
n = n 2 − 2 , (2.65)
n
y las funciones de onda son
(0,b)
R1 (r) = e−br sin rC1 (cot r), (2.66)
−br/2 2 (0,b)
R2 (r) = e sin rC2 (cot r), (2.67)
···
Rn (r) = e−br/n sinn rCn(0,b) (cot r). (2.68)
2
Un valor para la constante de normalización viene dado en [40], pero solamente es válido
para n = 1.
2.4. PROPAGADOR INSTANTÁNEO 21

En la Figura 2.2 se muestra el cuadrado de la función de onda para los primeros

niveles de energía. Dado que el máximo se encuentra muy cercano a uno, se
infiere que el diquark se encuentra altamente localizado. En la siguiente sección

(a) n = 1

(b) n = 2

Figura 2.2: Cuadrado de la norma de la función de onda para a = 0, b = 1,5.

se considera el potencial en el espacio de momentos. Esto es el punto de partida

para identificar un propagador instantantáneo del gluón.

2.4. Propagador Instantáneo

El contenido de esta sección tomado directamente de [12]; solamente se ha
detallado el cálculo de la transformada de Fourier del potencial.
Los potenciales en el espacio de los momentos obtenidos mediante una trans-
formación de Fourier de potenciales centrales son de gran interés dentro de una
22 CAPÍTULO 2. MODELOS POTENCIALES

amplia variedad de problemas físicos, desde materia condensada hasta física de

partículas. Pueden verse como propagadores instantáneos de los campos me-
diadores las respectivas interacciones y son especialmente importantes en los
cálculos de algunos cuerpos de Faddeev, que resultan ser más eficientes cuando
se realizan en el espacio de momentos. Desafortunadamente, los potenciales de
potencias de una amplia variedad de problemas como el oscilador linear y armó-
nico no tienen una integral de Fourier bien definida [46], siendo el inverso con
la distancia una entre pocas excepciones (1/r2 también tiene una transformada
de Fourier).

La interacción cot + csc2 fue introducida originalmente por Schrödinger en

[19] como una función angular en una superficie de tres dimensiones con una
curvatura constante y positiva, la cual es la hiperesfera S 3 . La hiperesfera S 3
se localiza dentro de un espacio Euclideano de 4 dimensiones E4 . Hasta una
constante aditiva, toma la forma

~2
VRM (χ) = −2B cot χ + l(l + 1) csc2 χ, (2.69)
2µd2

donde l es el valor del momento angular usual en 3D, d es una constante a ajustar
y χ es el segundo ángulo polar en E4 . Eligiendo la parametrización χ = r/R
para la variable angular, donde R es el radio fijo de S 3 , mientras r es la longitud
del arco en las superficie hiperesférica, VRM está usualmente dado en la forma
de un potencial en un espacio 3D curvo. El espacio Euclideano 3D E3 , embebido
en E4 , está descrito en términos del radio vector en valor absoluto |r|, definido
como |r| = R sin χ. Por lo tanto, la parametrización de χ correspondiente al
espacio 3D correctamente incrustado en E4 es χ = sin−1 |r|/R. La superficies
esférica se parametriza como x24 + r2 = R2 , con x4 = R cos χ y r = R sin χ.
La naturaleza del espacio, ya sea plano ó curvo, es de menor importancia para
el espectro de energías y las funciones de onda; todo se reduce a la interpretación
de R. En el espacio plano, R es visto como algún parámetro de longitud a ajustar,
mientras que en S 3 se considera como un parámetro que describe la curvatura.
Así, sin tomar en cuenta las transformaciones integrales como la transformada de
Fourier al espacio de momentos, la naturaleza del espacio adquiere importancia a
través de la definición de la integral de volumen. Si se intenta hacer la integral de
Fourier de ondas planas en un volumen en el espacio plano E3 , y se encuentra que
el potencial VRM (χ = r/R) es una función de arc tan r, el cual es inconsistente
y diverge. Por otro lado, considerando a VRM como una función del radio vector
del espacio plano E3 dentro de E4 , permite a la integral de Fourier ser tomada
de forma cerrada.
Tratar la interacción bajo discución como una función angular en S 3 es
posible debido a su simetría SO(4). La ventaja de tomar este punto de vista viene
de que la ecuación de Schrödinger con el potencial cot + csc2 está cercanamente
relacionado al problema de eigenvalores del momento angular 4D en S 3 . La
parte angular Π̂ del operador de Laplace-Beltrami es proporcional al operador
2.4. PROPAGADOR INSTANTÁNEO 23

de momento angular cuadrático 4D K2 ,

L2 (θ, ψ)
 
1 ∂ 2 ∂ 1
Π̂ = sin χ − = −κK2 , κ= . (2.70)
sin2 χ ∂χ ∂χ sin2 χ R2
Aquí, L2 (θ, ψ) es el operador de momento angular 3D usual en E3 , χ es el
segundo ángulo polar en E4 (χ ∈ [0, π]), mientras que κ es la constante de
curvatura. Como consecuencia, la ecuación de Schrödinger en S 3 es
 2 
~
κK2 − E(κ) ψ(χ, κ) = 0, (2.71)
2µ
donde µ puede tomarse como la masa reducida del diquark. El problema de
eigenvalores de K2 es [47]
 
2 K K
K |Klmi = K(K + 2)|Klmi, |Klmi ∈ , . (2.72)
2 2
Los niveles |Klmi pertenecen a la representación irreducible SO(4) del ti-
po (K/2, K/2), mientras que los números cuánticos K, l y m definen los ei-
genvalores de los operadores de momento angular en 4, 3 y 2 dimensiones,
respectivamente. Esos números cuánticos corresponden a la cadena de reduc-
ción SO(4) ⊃ SO(3) ⊃ SO(2) y satisfacen las reglas; l = 0, 1, 2 . . . , K y
m = −l, . . . , l. Multiplicando la ecuación (2.71) por (− sin2 χ) y cambiando la
variable a ψ(χ, κ) = sin χS(χ, κ), se obtiene la siguiente ecuación de Schrödinger

~2 d2
 
−κ + U l (χ, κ) = S(χ, κ) = E(κ)S(χ, κ),
2µ dχ2
~2
Ul (χ, κ) = κ l(l + 1) csc2 χ, (2.73)
2µ
con Ul (χ, κ) interpretado como el término de una barrera centrífuga en S 3 .
Una interpretación de las ecuaciones (2.71) y (2.73), es que el potencial csc2 ;
en representación de la barrera centrífuga en la hiperesfera 3D, tiene a so(4)
como álgebra potencial. Es importante observar que el álgebra potencial so(4)
permanece inalterada una vez que el término csc2 ha sido agregado a la fun-
ción armónica cot χ. Esto puede apreciarse del hecho de que la energía continúa
siendo una función que depende solo de los eigenvalores de K2 , K(K + 2), que
traslada al número cuántico principal n en n = K + 1. En efecto, la simetría
algebraica so(4) de la interacción cot + csc2 permite considerarla como una fun-
ción en S 3 , una circunstancia que facilitará sustancialmente su transformación al
espacio de momentos. Aquí se adopta la siguiente parametrización del potencial
RMT como una función del segundo ángulo polar χ, sobre S 3 y la curvatura:
√ ~2 l(l + 1)
V(χ) = −2G κ cot χ + κ . (2.74)
2µ sin2 χ
En coordenadas cartesianas, el término cot es igual a x4 /|r| y describe dos
potenciales difiriendo por un signo de S 3 , correspondientes a las interacciones
24 CAPÍTULO 2. MODELOS POTENCIALES

en los hemisferios norte y sur.

La transformada de Fourier al espacio de momentos es

√ 2µ ∞
Z Z 2π Z π
4πΠ(|q|) = −2G κ 2 d|x||x|3 δ(|x| − R) dψ dθ sin θ
~ 0 0 0
Z π/2,π sin χ
i|q|| √ | cos θ
× dχ sin2 χe κ~2 cot χ, (2.75)
0,π/2

donde la función δ(|x|−R) restringe E4 a S 3 . Aquí, el plano 4D ha sido evaluado

tomando el eje z como referencia del vector de momento y el vector de posición
de la partícula confinada con una proyección no cero en el eje extra de E4 :
|r|
i|q| sin χ
√ cos θ sin χ
eiqx/~ = ei|q| ~ cos θ
=e ~ κ , |r| = √ . (2.76)
κ

En S 3 , uno tiene que distinguir entre dos tipos de potenciales en el espacio de los
momentos. El primero, mostrado en la Figura 2.3, corresponde a x4 positivo que
va con χ ∈ [0, π/2] y describe un |r| que va incrementando. El siguiente paso es
obtener una expresión completa para el propagador instantáneo de la ecuación
(2.75) en el hemisferio norte. Para ello, nótese directamente que el valor de la
integral sobre |x| es R3 , mientras que el de la integral sobre ψ es 2π. Así, se
tiene
Z 1 Z π/2
√ µ i qu
√ sin χ
4πΠ(|q|) = −8πG κ 2 R3 duχ sin χ cos χe ~ κ , (2.77)
~ −1 0

donde se hizo el cambio cos θ →√u. Considérese la integral sobre χ en la ecuación

(2.77). Para ello, se toma qu/~ κ → α. Entonces,
Z π/2  2 Z iα
1
dχ sin χ cos χeiα sin χ = dyyey
0 iα 0
 iα
2  iα 
1  y
y

= ye  − e 
iα 0 0
 2  
1
= iαeiα − eiα + 1
iα
(1 − iα)eiα − 1
= . (2.78)
α2
Ahora solo falta hacer la integral sobre u. Antes de calcularla, nótese la expresión
que se tiene para el propagador

√ µ 3 1 ~2 κ
Z   
qu i ~qu
√
Π(|q|) = 2G κ 2 R du 2 2 i √ − 1 e κ +1 . (2.79)
~ −1 q u ~ κ

Al integrar el término que no está multiplicado por el factor exponencial, se

obtiene inmediatamente −2~2 κ/q 2 . Para integrar el término con el factor expo-
2.4. PROPAGADOR INSTANTÁNEO 25

Figura 2.3: Función de revestimiento en el espacio de los momentos sobre el he-

misferio norte. El hemisferio sur aparece como un reflejo con respecto
al plano horizontal.

nencial, primero se separa en la forma

1 1 √
~2 κ ~2 κ i ~qu
Z   Z  
qu i qu
√ ~ κ i ~qu
√ √
du 2 2 i √ − 1 e ~ κ = i e κ − 2 2e κ .
−1 q u ~ κ −1 qu q u
(2.80)

Integrando el primer sumando por partes,

1 √ 1 Z 1 2
~2 κ i ~qu
Z
~ κ i ~qu
√ √
 ~ κ i ~qu√
e κ = 2
e κ + 2 u2
e κ. (2.81)
−1 qu q u 
−1 −1 q

El segundo término de la ecuación (2.80) se cancela con el de la ecuación (2.81).

26 CAPÍTULO 2. MODELOS POTENCIALES

Así, se obtiene que

 2 1
√ µ 2~2 κ

~ κ i qu
√
2G κ 2 R3 2 e ~ κ  − 2

Π(|q|) =
~ q u −1 q
 2 
√ µ 3 ~ κ i √q 2~2 κ
 
q
−i ~√
= 2G κ 2 R e ~ κ +e κ − 2
~ q2 u q
√ µ 2
 
~ κ q
= 4G κ 2 R3 2 cos √ − 1
~ q ~ κ
√ µ 3 ~2 κ q
= 8G κ 2 R 2 sin2 √
~ q 2~ κ
2 q√
2Gµ sin 2~ κ
= 2
~2 κ

q√
2~ κ
q√
sin2 2~ κ 2Gµ
= c 2 , c= . (2.82)
q√
~2 κ
2~ κ

La expresión obtenida en la ecuación (2.82) corresponde al propagador del gluón

para el sistema quark-diquark con una interacción que tiene la simetría algebrai-
ca so(4) ⊂ so(2, 4), siendo so(2, 4) el álgebra conforme. Sin entrar en detalles,
es importante notar que el operador de Casimir del álgebra so(4), y en con-
secuencia de so(4, 2), no es el mismo que el operador de Casimir que genera
el de los grupos SO(4) ⊂ SO(2, 4). Se ha mostrado en [48] que al agregar el
potencial cot a K2 se obtiene un operador de Casimir de un álgebra so(4) en
una representación nueva y unitariamente no equivalente a la del grupo SO(4).
Esto debido a que la intensidad ~b2 /2µ del potencial cot introduce una escala
de masa, que se podría asociar como la masa del dilatón, la cual rompe la si-
metría del grupo conforme pero si preserva el álgebra. Esta última ya no genera
el grupo conforme, pero sigue definiendo la dinámica conforme en términos de
los eigenvalors del nuevo operador de Casimir. Es notorio que el álgebra es una
entidad matemática que, en una de las posibilidades, se podría relacionar con
grupos ortogonales o unitarios, pero en general esto no es posible. Los valores de
los parámetros G, µ y κ utilizados para ajustar el espectro de nucleones y ∆’s
son G = 0,0493GeV − f m, µ = 0,3213GeV y κ = 1,03826f m−2 [18]. En nues-
tro caso, la forma funcional de este propagador será usado en las ecuaciones de
Schwinger-Dyson, que además será comparado con un modelo bien establecido
conocido como modelo de Maris-Tandy [29].
En el siguiente capítulo se describen, en términos generales, las predicciones
de las ESD para la generación dinámica de masas y confinamiento en la física
hadrónica.
Capítulo 3

Generación Dinámica de
Masas y Confinamiento

La Generación Dinámica de Masa (GDM) –debida al rompimiento dinámico

de la simetría quiral– es un fenómeno no perturbativo que se estudia a través
de las Ecuaciones de Schwinger-Dyson (ESD) [49]. Por medio de las ESD, se
explica cerca del 99 % de la masa visible en el universo, que se produce debido
a fenómenos de auto-interacción en QCD, masa que no se produce mediante el
rompimiento espontáneo de la simetría quiral (mecanismo de Higgs [50]). Las
ESD se obtienen de manera formal en Teoría Cuántica de Campos [51], pero en
esta tesis se considera un enfoque más amigable, ya que para nuestro objetivo
únicamente se necesita entender el contenido físico de estas ecuaciones.

En este capítulo se estudian la generación dinámica de masas por medio

de las ESD y el fenómeno de confinamiento en QCD para el modelo de Maris-
Tandy (MT [29]). En primer lugar, se explica la manera en que surgen las ESD,
seguido de una breve presentación de los elementos que conforman las funciones
de Green que aparecen en las ESD. Posteriormente se presentan las ESD en QCD
en la aproximación arcoiris y se resuelven para el modelo MT. Para finalizar, se
estudia el confinamiento para dicho modelo.

3.1. Ecuaciones de Schwinger-Dyson

Los propagadores se modifican a través de autointeracciones. Para el caso del
quark, al propagador que representa la suma de todas las autointeracciones se
le conoce como como propagador completo, el cual posee todas las propiedades
de propagación de dicho fermión. Para el caso en que el propagador no emite ni
absorbe fotones se le nombra propagador desnudo. Existe una serie infinita de
correcciones radiativas al propagador que, mostrado en términos de diagramas,

27
28CAPÍTULO 3. GENERACIÓN DINÁMICA DE MASAS Y CONFINAMIENTO

1
es

= +

+ +

+ +

+ + + ...

Figura 3.1: Correcciones radiativas al propagador del fermión.

Estas correcciones están clasificadas en tres tipos: la del propagador del fer-
mión, del fotón y del vértice. Tal como aparecen en la Figura 3.1, éstas correccio-
nes no pueden sumarse. Dyson reunió todas esas contribuciones en una sola, a la
cual llamó autoenergía [53], que se representa por Σ(p) y esquemáticamente por

Figura 3.2: Autoenergía Σ(p).

La autoenergía involucra todas las correcciones a las funciones de Green in-

volucradas en el propagador, como es indicado por los puntos sobre las partes
correspondientes en la Figura 3.2. En términos de Σ(p), la expansión perturba-
tiva para el propagador del fermión se muestra en la Figura 3.3.
1
Se ejemplifica como obtener las ESD en la electrodinámica cuántica. Más tarde se dirá
como traducir estas ecuaciones para QCD.
3.1. ECUACIONES DE SCHWINGER-DYSON 29

= + + + ...

Figura 3.3: Propagador del fermión en teoría de perturbaciones en términos de

autoenergías.

Así, en términos de la expansión de Dyson, la Figura 3.3 se puede escribir

como

SF (p) = SF0 (p) + SF0 (p)Σ(p)SF0 (p) + · · ·

= SF0 (p)(1 + Σ(p)SF0 (p) + Σ(p)SF0 (p)Σ(p)SF0 (p))
" #
1
= SF0 , (3.1)
1 − Σ(p)SF0 (p)

(0)
donde SF (p) representa al propagador completo y SF (p) a su contraparte des-
nuda. Si se factoriza SF0 (p)Σ(p)SF0 (p) a partir del segundo término se obtiene

SF (p) = SF0 (p) + SF0 (p)Σ(p)SF0 (p)(1 + SF0 (p)Σ(p) + · · · )

( )
0 0 0
h 1 i
= SF + SF (p)Σ(p) SF (p) . (3.2)
1 − Σ(p)SF0 (p)

El término entre llaves se identifica nuevamente como el propagador completo

SF (p). Si se sustituye el término del propagador completo dentro de la ecuación
(3.2) se obtiene
SF (p) = SF0 (p) + SF0 (p)Σ(p)SF (p). (3.3)
La expresión anterior puede aún manipularse. Por ejemplo, si se multiplica
SF−1 (p) por la izquierda, se obtiene

I = SF0 (p)SF−1 (p) + SF0 (p)Σ(p), (3.4)

y luego se multiplica [SF0 (p)]−1 por la derecha

SF0 (p) = SF−1 (p) − Σ(p). (3.5)

Finalmente, los términos se rearreglan como

SF−1 (p) = (SF0 (p))−1 − Σ(p), (3.6)

la cual es la ESD que corresponde a la función de Green de dos puntos del

fermión. Esquemáticamente está representada por la Figura 3.4a
30CAPÍTULO 3. GENERACIÓN DINÁMICA DE MASAS Y CONFINAMIENTO

-1 -1
= -

(a)

-1 -1
= -

(b)

Figura 3.4: (a) ESD para el propagador del fermión. (b) ESD para el propagador
del fotón.

= +

Figura 3.5: ESD para el vértice fermión-fotón.

3.1. ECUACIONES DE SCHWINGER-DYSON 31

Cada una de las funciones de Green que aparecen en la autoenergía obede-

cen su propia ESD. Las correspondientes al propagador del fotón y al vértice se
representan en las Figuras 3.4b y 3.5.

Nótese que las funciones de Green de dos puntos están relacionadas con las
de tres puntos; la de tres puntos con la de dos y la de cuatro y así sucesivamente,
hasta formar una torre infinita de relaciones entre las funciones de Green. Para
poder extraer información física relevante de las ESD, se debe truncar esa torre
infinita de ecuaciones suponiendo una forma para las funciones de Green de
varios puntos. Para tener GDM usualmente se trunca dicha torre de ESD al
nivel de propagadores. En QCD es necesario conocer el propagador del quark
SF y el propagador del gluón Dµν en el espacio Euclideano. El propagador
completo del quark es
F (p2 ) F (p2 )(6p + M (p2 ))
S(p) = 2 = , (3.7)
−i6p + M (p ) p2 + M 2 (p2 )

donde 6p = γµ pµ 2 (siendo γµ las matrices de Dirac [51]). Por otro lado, F (p2 )
es un objeto que corrige a la función de onda llamado renormalización de la
función de onda del fermión. M (p2 ) es la función de masa.
Por otro lado, el propagador del gluón en la norma de Landau es
 
qµ qν
Dµν (q) = D(q 2 ) δµν − 2 , (3.8)
q
donde D(q 2 ) es un término de interacción efectivo.
En la siguiente sección se estudia la ESD para el propagador del quark. Basta,
para ello, reemplazar en los diagramas y expresiones anteriores al propagador
del fotón por el del gluón de QCD.

3.1.1. Ecuación de GAP en QCD

A continuación, se procede a estudiar la ecuación de gap3 . De acuerdo al
esquema estándar, el punto de partida consiste en expresar al propagador com-
pleto del quark en la forma general [54] en el espacio Euclideano
S −1 (p) = Z2 (i6p + mb ) + Z1 Σ(p), (3.9)
donde la autoenergía del quark es
d4 k 2 λa
Z
Σ(p) = g Dµν (k − p) γµ S(k)Γaν (k, p). (3.10)
(2π)4 2
La ecuación de gap de QCD se muestra esquemáticamente en la Figura 3.6.
Aquí, q = k − p, λa representa las matrices de Gell-Mann (ecuación (A.7)),
2
Se está trabajando en un espacio Euclideano, donde las matrices de Dirac satisfacen
{γµ , γν } = 2δµν .
3
A las ESD para el propagador del fermión también se le conoce como ecuación de gap o
“brecha”.
32CAPÍTULO 3. GENERACIÓN DINÁMICA DE MASAS Y CONFINAMIENTO

mb destaca como la masa de corriente de los quarks, g es la constante de aco-

plamiento del gluón, Dµν (q) es el propagador del gluón y Γaν (k, p) es el vértice
completo quark-gluón. Además, Z1 y Z2 son, respectivamente, las constantes de
renormalización de la función de onda del vértice y el quark. En esta sección se
exponen los resultados con la ecuación sin renormalizar. El tratamiento propia-
mente renormalizado es un trabajo en progreso y será presentado en [52].
El punto de partida es hacer un modelo para el producto g 2 Dµν (k − p)Γaν (k, p).

Figura 3.6: Ecuación de GAP en QCD.

Se considera la aproximación arcoiris, donde Γaν (k, p) = γν , entonces

a
λa
 
2 a 2λ 2 qµ qν
g Dµν (q)Γν (k, p) = g D(q)γν ∆µν = γν D(q ) δµν − 2 . (3.11)
2 2 q
Ahora, al sustituir el valor del vértice de la ecuación (3.11) en la ecuación (3.10)
se tiene
λa λa
Z  
qµ qν
S −1 (p) = S0−1 (p) − d4 k γµ S(k) γν D(q 2 ) δµν − 2 . (3.12)
2 2 q
El método usual para resolver la ecuación de gap en QCD consiste en separar
la ecuación (3.12) en dos funciones para M (p) y F (p). Primero, se escribe la
expresón completa para el propagador
a
λa
 
6p + M (p)
Z
4 λ F (k) 2 qµ qν
= 6p −mb − d k γµ 2 γν D(q ) δµν − 2 . (3.13)
F (p) 2 k + M 2 (k 2 ) 2 q
Para encontrar la función de masa M (p) se calcula la traza en la ecuación (3.13),
entonces
M (p) X λa λa Z d4 k G(q 2 ) F (k)

qµ qν

4 = 4mb + Tr[γ µ (6 k +M (k))γ ν ] δ µν − .
F (p) a
4 (2π)4 q 2 k 2 + M 2 (k 2 ) q2
P a a (3.14)
Usando el hecho de que a λ λ = (16/3)I y que la traza de productos de
matrices γ tiene las propiedades

 T r[#impar de γ 0 s] = 0,


T r[γα γβ ] = dδαβ , (3.15)

T r[γα γβ γµ γν ] = d[δαβ δµν − δαµ δβν + δαν δβµ ],


se cumple
X λa λa 16
Tr[γµ (6k + M (k))γν ] = M (k)δµν . (3.16)
a
4 3
3.1. ECUACIONES DE SCHWINGER-DYSON 33

Contrayendo los índices

 
qµ qν
δµν δµν − 2 = 3. (3.17)
q
Finalmente, la expresión para la función de masa
d4 k M (k)F (k)
Z
M (p) 4
= mb + D(q 2 ). (3.18)
F (p) 3 (2π)4 k 2 + M 2 (k 2 )
Para integrar expresiones en cuatro dimensiones se usan coordenadas hiperesfé-
ricas

d4 k = k 3 sin2 θ sin φdkdθdφdψ

q2 = (k − p)2 = k 2 + p2 − 2kp cos θ, (3.19)

donde θ el ángulo entre los momentos de los fermiones. Además

ψ ∈ [0, 2π], φ ∈ [0, π], θ ∈ [0, π], k 2 ∈ (0, Λ2 ), (3.20)

ya que se ha introducido una regularización de las integrales a la escala Λ.

Integrando sobre las coordenadas angulares independientes
Λ2 π
M (p2 ) 2
M (k 2 )F (k 2 )
Z Z
1 2k
= mb + 3 dk dθ sin2 θD(q 2 ), (3.21)
F (p2 ) 2π 0 k 2 + M 2 (k 2 ) 0

que es la ecuación para la función de masa.

Para calcular la ecuación para la función de renormalización de la función
de onda, la expresión dada por la ecuación (3.13) es multiplicada por 6p previo
al cálculo de la traza, obteniendo así
Z Λ2 Z π
1 1 k 2 F (k 2 )
= 1− 3 2 dk 2 2 dθ sin2 θD(q 2 ) ×
F (p2 ) 6π p 0 k + M 2 (k 2 ) 0
2(k 2 − k · p)(k · p − p2 )
 
× − k · p . (3.22)
q2
Nuestro primer objetivo será resolver las ecuaciones (3.21) y (3.22) el modelo
de Maris-Tandy.

3.1.2. Otros Observables

Algunos observables relevantes de los hadrones pueden obtenerse también
directamente del propagador del quark. Por ejemplo, el condensado quiral del
vacío, el cual se forma por las autointeracciones de las partículas y que por
lo tanto, proporciona un medio pegajoso en el cual los quarks de valencia se
propagan y adquieren masase calcula como
Z Λ 4
d k
−hψ̄ψi = lı́m NC TrD S(k), (3.23)
Λ→∞ (2π)4
34CAPÍTULO 3. GENERACIÓN DINÁMICA DE MASAS Y CONFINAMIENTO

donde NC = 3, además de que se calcula en el límite quital mb = 0. Un segundo

ejemplo es la constante de decaimiento del pión, que puede calcularse a través
de la ecuación de Pagel-Stokar [59].

NC Λ 2 k 2 F (k 2 ) k2 2
Z  
2 2 2 2 dM (k )
fπ = − 2 dk M (k ) − M (k ) . (3.24)
4π (k 2 + M (k 2 ))2 2 dk 2

3.1.3. Confinamiento en las Ecuaciones de Schwinger-Dyson

El fenómeno de confinamiento no se puede entender a partir de primeros
principios, debido a la imposibilidad de detectar estados de color, por lo que
solamente existe evidencia experimental de ello. Sin embargo, el confinamiento
del quark se encuentra decodificado en su propagador. Se puede observar a través
de la violación del axioma de reflexión de positividad de Osterwalder-Schrader
[55]. El propagador del quark puede escribirse en el espacio Euclideano como
S(p) = −i6pσν (p2 ) + σs (p2 ), (3.25)
donde las partes vectorial y escalar del propagador, σν (p2 ) y σs (p2 ), están ex-
presados en términos de la función de renormalización F (p2 ) y la función de
masa M (p2 ) como
F (p2 )
σν (p2 ) = , (3.26)
p2 + M 2 (p2 )
F (p2 )M (p2 )
σs (p2 ) = . (3.27)
p2 + M 2 (p2 )
Se define el promedio espacial de la parte escalar escalar del propagador [57, 58]

d4 p
Z Z
3
∆(t) = d x σs (p2 ). (3.28)
(2π)4
Si esta función de Schwinger es positiva definida, entonces describe una excita-
ción estable. Por ejemplo, para una partícula libre, ∆(t) = e−mt /2, donde m
es la masa de la partícula. En caso contrario, describe una partícula confinada
[55]. Un ejemplo típico es cuando la masa de la partícula desarrolla una parte
imaginaria, m = a + ib, de modo que ∆(t) = e−at cos(bt + δ). Si la partícula es
compuesta, b puede entenderse como un ancho de decaimiento, pero si se está
refiriendo a partículas fundamentales, entonces b 6= 0 indica que la partícula
está confinada.
Ahora ya se tienen las herramientas para calcular diferentes observables y
la tarea se reduce a proponer distintos modelos que se predigan o se ajusten a
resultados experimentales.

3.2. Modelo de Maris-Tandy para QCD

Con los elementos vistos en la sección anterior, ahora todo el trabajo se
reduce a calcular los observables mencionados para diferentes modelos. Para
3.2. MODELO DE MARIS-TANDY PARA QCD 35

esta sección, se pondrá a prueba el bien conocido modelo de Maris-Tandy [29].

Este modelo ha sido usado extensivamente en los estudios de fenomenología de
hadrones que involucran a las ESD. La función de revestimiento de este modelo
es
4π 2 Dq 2 −q2 /ω2 8π 2 γm Z(q 2 )
D(q) = e + , (3.29)
ω ln[τ + (1 + q 2 /ΛQCD )2 ]
donde
12
γm = , mz = 0.5GeV, NF = 4,
32 − 2NF
τ = e2 − 1, ΛQCD = 0.254GeV, ω = 0.4GeV,
s
−
D = 0.93GeV 2 , Z(s) = (1 − e 4m2z )/s.

En la Figura 3.29 se muestra la solución de la función de masa en el límite quiral

mb = 0. Se puede notar un crecimiento de la función de masa en la región infra-
rroja, que muestra el rompimiento dinámico de la simetría quiral. Se encuentra
que la masa generada dinámicamente aparece dentro del rango de 300-500MeV,
como es esperado de la fenomenología de quarks ligeros. Nótese también que los
efectos no perturbativos dominan en la región 0 < p2 < 1GeV2 . A altas energías,
la función de masa cae suavemente hacia la masa asintóticamente corriente de
los quarks libres, que es una predicción bien conocida desde hace años en los
estudios de las ESD.

Figura 3.7: Función de masa del modelo MT para el límite quiral.

La función de renormalización de la función de onda se muestra en la Figura 3.8.

En el siguiente sección se menciona su comportamiento dado que se compara
con la función de renormalización de otros modelos. Los valores obtenidos para
el condensado y fπ son

(−hψ̄ψi)3 = (0.261722GeV)3 , fπ = 0.0842529GeV. (3.30)

36CAPÍTULO 3. GENERACIÓN DINÁMICA DE MASAS Y CONFINAMIENTO

Figura 3.8: Función de renormalización para el modelo MT.

En la Figura 3.9 se muestra que nuestro modelo es confinante; el mínimo indica

que la función ∆(t) pasa por cero, esto indica que hay un cambio signo y por
tanto, no es positiva definida. El inverso de la posición donde se localiza el pico
ayuda a dar un parámetro de orden para el confinamiento; para este modelo co-
rresponde a una escala de unos cuantos MeV. En el siguiente capítulo se propone

Figura 3.9: Test de confinamiento para el Modelo MT.

un término de interacción efectivo basados en el propagador vestido del gluón

con una ventana comforme. Nuestra meta es reproducir lo más cercanamente
posible el escenario descrito en este capítulo.
Capítulo 4

Modelo Para el Propagador

Vestido del Gluón

Con los resultados de los capítulos previos, ahora se puede hacer un estudio
del propagador instantáneo del gluón obtenido en el Capítulo 2. El modelo que
será usado tratará de reproducir la altura de la función de masa y además se
dará una prueba de confinamiento como en el Capítulo 3 para el modelo MT.
Por otro lado, tres propagadores “hermanos” al presentado en el Capítulo 2 son
presentados y comparados con el modelo MT.
Motivados por la forma del propagador instantáneo del gluón (ecuación
(2.82)), se promueve la siguiente forma efectiva para la función de vestimiento
del gluón,
sin2 2~q√κ
2
D(q ) = c  2 , (4.1)
q√
2~ κ

más allá de la aproximación instantánea, ajustando el valor de c –que determina

la intensidad de la interacción– con el fin de obtener una altura de la función
de masa lo más cercana posible a la encontrada en el modelo de Maris-Tandy.
Además, se incluyen tres variantes para el revestimiento efectivo del gluón. Los
modelos son
H1 (q)
Modelo I : Struve : D(q 2 ) = c , (4.2)
q2
sin(q)2
Modelo II : Sinc : D(q 2 ) = c , (4.3)
q2
2/π
Modelo III : StruveLineal : D(q 2 ) = c 2 , (4.4)
q +3
q 2 + m2
Modelo IV : RGZ : D(q 2 ) = c 4 , (4.5)
q + q 2 m2 + m4
donde H1 (x) es la función de Struve. El Modelo II consite en una adaptación

37
38CAPÍTULO 4. MODELO PARA EL PROPAGADOR VESTIDO DEL GLUÓN

a la forma funcional de la función de revestimiento mostrada en la ecuación

(4.1); los cambios han sido absorbidos por la constante c. El Modelo III es una
linearización de la función H1 (x)/x2 y el Modelo IV corresponde al propagador
que captura la generación dinámica de masas para gluones en el escenario de
confinamiento de Gribov-Zwanzinger refinado [60, 61], con mb = 2.167GeV
(ecuación (4.6)).
q 2 + m2b
D(q 2 ) = 4 . (4.6)
q + mb q 2 + m4b
Las funciones de revestimiento se pueden ver en la Figura 4.1 con los valores c
mostrados en el Cuadro 4.1 que arrojan una masa dinámica para los quarks de
alrededor de 500M eV . De esa figura, se observa que el revestimiento instantáneo
es finito en el infrarrojo y cae de manera oscilatoria en la forma 1/x2 para los
Modelos I y II. Esas oscilaciones son el remanente de desdoblar un volumen finito
sobre uno infinito. Además, se espera que estos modelos sean válidos a una escala
de p2 ≤ 1GeV 2 [56]. Las ecuaciones (3.21) y (3.22) han sido resueltas usando el

Figura 4.1: Revestimiento efectivo del gluon para los modelos (4.3-4.5).

código del Apéndice B, sustituyendo las funciones de revestimiento (4.2-4.5). En

la Figura 4.2 se muestra la función de masa de los cuatro modelos junto con el
valor de la función de masa del modelo de MT. En dicha figura se puede apreciar
que los cinco modelos describen una masa dinámica de los quarks dentro de los
rangos fenomenológicos aceptados. La función de masa perteneciente al modelo
MT cae a cero de la forma dictada por Teoría de Perturbaciones antes que los
otros modelos, que van a cero de manera muy similar. Esto se traduce en un
valor mayor del condensado quiral como se muestra en el Cuadro 4.1.
Para la función de renormalización de la función de onda, se obtienen re-
sultados distintos para cada modelo, tal como puede verse en la Figura 4.3.
La función de renormalización para los modelos RGZ y Struve Lineal son muy
similares en el infrarojo, comenzando como una constante mayor que la uni-
 39

Modelo c M (0)(GeV ) (−hψ̄ψi(GeV ))3 fπ (GeV )

MT 1 0.49 (0.261722)3 0.0842529
Sinc 29.97 0.5 (0.375583)3 0.163578
Struve 21.72 0.5 (0.428904)3 0.2173
Struve Lineal 28.88 0.5 (0.443037)3 0.226299
RGZ 17.5 0.5 (0.446159)3 0.231806

Cuadro 4.1: Valor de los observables para los diferentes modelos al ajustar a la
altura de la función de masa de MT.

Figura 4.2: Función de masa para diferentes modelos.

40CAPÍTULO 4. MODELO PARA EL PROPAGADOR VESTIDO DEL GLUÓN

dad y disminuyendo suavemente hasta el límite ultravioleta. Mientras que, por

otro lado, los modelos Struve y Sinc tiene un comportamiento bastante dife-
rente a momentum pequeño; especialmente el Sinc, que muestra oscilaciones a
momentum intermedio. Algo a notar es que ninguna se asemeja a la función de
renormalización del modelo MT, que es una constante menor que 1 en el límite
infrarojo, pero converge a 1 conforme p crece. Esto indica que la renormalización
de la función de onda que se obtiene en nuestros modelos no será necesariamente
igual a la del modelo MT, excepto el límite ultravioleta, donde todos los modelos
covergen al mismo valor F (p2 ) = 1, como se espera en el límite perturbativo.

Figura 4.3: Renormalización de la función de masa para diferentes modelos al

ajustar a la altura de la función de masa de MT.

La diferencia cualitativa de la solución a la ecuación de gap con el propagador

gluónico en estudio se refleja en que no solo el condensado quiral, sino también
la constante de decaimiento leptónico del pión se sobreestiman, como se muestra
en la última columna del Cuadro 4.1. A groso modo, hψ̄ψi y fπ para nuestros
modelos son al menos el doble que en el modelo MT. Por otro lado, en la
Figura 4.4 se muestra la prueba de confinamiento hecha sobre los propagadores
obtenidos dinámicamente con los distintos modelos. En esta ocasión, se nota
que el modelo MT tiene una escala de confinamiento menor a la de los modelos
propuestos.
Debido a que en los observables hadrónicos no es la masa constituyente de
los quarks, sino el condensado quiral lo que se determina experimentalmente,
ahora se procede a ajustar la constante c para obtener un condensado quiral
similar al del modelo MT. Los resultados obtenidos se muestran en las Figuras
(4.5)-(4.8) y en el Cuadro 4.2.
En la Figura 4.5 se nota que cada modelo debe pesarse de forma distinta para
obtener valores similares del condensado quiral. La masa gluónica más grande
se obtiene con el Modelo IV. En la Figura 4.6 puede apreciar una dramática
 41

Figura 4.4: Test de confinamiento para los diferentes modelos al ajustar a la

altura de la función de masa de MT.

Modelo c M (0)(GeV ) (−hψ̄ψi(GeV ))3 fπ (GeV )

MT 1 0.49 (0.261722)3 0.0842529
Sinc 25.4 0.0916076 (0.270205)3 0.049967
Struve 20.23 0.0835128 (0.297258)3 0.0578297
Struve Lineal 27.5 0.0605863 (0.280298)3 0.0451153
RGZ 16.66 0.0530895 (0.273333)3 0.0412857

Cuadro 4.2: Resultados cuando se ajusta al condensado del modelo MT

Figura 4.5: Revestimiento efectivo del gluon ajustando al condensado de MT.

42CAPÍTULO 4. MODELO PARA EL PROPAGADOR VESTIDO DEL GLUÓN

Figura 4.6: Función de masa para diferentes modelos ajustando al condensado

de MT.

Figura 4.7: Renormalización de la función de masa para diferentes modelos ajus-

tando al condensado de MT.
 43

Figura 4.8: Test de confinamiento para los diferentes modelos ajustando al con-
densado de MT.

caída en la función de masa (menor de 0.1GeV para todos los modelos) y un

incremento bastante considerable en la escala de confinamiento (Figura 4.8).
Por otro lado, la constante de decaimiento del pión decrece aproximadamente a
la mitad del valor obtenido en el modelo MT (Cuadro 4.2).
Con el fin de tener una mejor aproximación en el infrarojo, se propone una
constante de acoplamiento que dependa del momento. Se va a tomar un caso en
que la constante c dependa del momento en la forma [62]

a0 + a1 q 2
 
πγm
c(q 2 ) = C + , (4.7)
1 + a2 q 2 + a3 q 4 + a4 q 6 ln(e + q 2 )

donde

a0 = 1.47 a1 = 0.881 a2 = 0.314

a3 = 0.00986 a4 = 0.00168, (4.8)

la constante C se ajusta para obtener la altura de la función de masa igual que

en el modelo de MT.
Los resultados obtenidos se muestran en las Figuras 4.9-4.12.
Se puede apreciar que en la función de masa (Figura 4.10) la caída de todos
los modelos se ha suavizado –especialmente en el Modelo II donde las oscilaciones
han desaparecido – y ocurre un poco antes que la mostrada en la Figura 4.2.
Por otro lado, la función de renormalización ha cambiado de manera bastante
notable para los Modelos I, III y IV y solamente tiene un ligero cambio en el
modelo II –las oscilaciones vuelven a desaparecer–.
Finalmente, en la prueba de confinamiento (Figura 4.12) no se nota cambio
significativo, pero se puede ver un ligero desplazamiento hacia la izquierda al
44CAPÍTULO 4. MODELO PARA EL PROPAGADOR VESTIDO DEL GLUÓN

Figura 4.9: Revestimiento efectivo del gluon para los modelos (4.3-4.5) con c =
c(q) y ajustando la altura de la función de masa a la de MT.

Modelo C M (0)(GeV ) (−hψ̄ψiGeV )3 fπ (GeV )

MT 1 0.49 (0.261722)3 0.0842529
Sinc 16.5 0.5 (0.338128)3 0.138744
Struve 23.58 0.5 (0.385949)3 0.184668
Struve Lineal 26.9 0.5 (0.384374)3 0.177719
RGZ 16.05 0.5 (0.39023)3 0.183999

Cuadro 4.3: Resultados obtenidos al proponer una constante que depende del
momento

comparar directamente el confinamiento para un modelos, dicha comparación

se muestra en la Figura 4.13 para el Modelo II. El Cuadro 4.3 muestra que
los valores observables obtenidos son más cercanos a aquellos obtenidos con
el Modelo MT que los mostrados en el Cuadro 4.1. Con esto se concluye la
presentación de nuestro trabajo.
 45

Figura 4.10: Función de masa para diferentes modelos con c = c(q 2 ).

Figura 4.11: Renormalización de la función de masa para diferentes modelos con

c = c(q 2 ).
46CAPÍTULO 4. MODELO PARA EL PROPAGADOR VESTIDO DEL GLUÓN

Figura 4.12: Test de confinamiento para los diferentes modelos con c = c(q 2 ).

Figura 4.13: Test de confinamiento para los diferentes modelos. La línea verde
representa al modelo con c = c(q 2 ).
Capítulo 5

Conclusiones y
Consideraciones Finales

En esta tesis, se han explorado dos diferentes caras del confinamiento: el pri-
mero en el sentido de una dinámica espacial efectiva para quarks y gluones sobre
un volumen finito delimitado por el potencial − cot; del cual se han obtenido
las primeras soluciones tomando l = 0. La otra cara es dentro de la descripción
teórica del confinamiento de quarks a través de las ESD. En cualquier sistema
dado, las restriciones espaciales siempre son consideradas en teorías clásicas sin
preocuparse mucho sobre la naturaleza de tales restricciones. Por otro lado, de
QCD se sabe que no hay estados de color libres, hecho confirmado por experi-
mentos. Sin embargo, la dinámica espacial de quarks y gluones es un problema
abierto. El potencial de Cornell puede considerarse como una primera aproxi-
mación al potencial − cot. La transformada de Fourier al espacio de momentos
puede verse como la función de revestimiento instantánea del gluón y toma la
forma vista en la ecuación (4.1).
Se ha adoptado tal función de revestimiento en una forma efectiva para
incluir confinamiento espacial en nuestro modelo, y entonces, se estudia el RDSQ
y confinamiento por medio de las ESD para el propagador del quark. Se observa
un crecimiento dinámico de la función de masa M (p2 ) en la región infrarroja
dirigida por el propagador efectivo del gluón, que también provee la caída 1/p2
esperada en el dominio perturbativo. Los efectos de desdoblar un volumen finito
dentro de un espacio infinito no modifica la forma asintótica de M (p2 ) como se
ve en la Figura 4.2. El propagador viola el axioma de reflexión de positividad,
como se observa en su promedio espacial (Figura 4.4), indicando que la función
de revestimiento efectivo es responsable para el confinamiento de quarks. Este
es el resultado principal que presenta esta tesis.
Por otro lado, cuando se trato de ajustar a la altura de la función de masa
del modelo MT se obtuvo un valor del condensado más grande y una constante
de decaimiento del pión un poco más del doble (Cuadro 4.1), además de que
da una escala de confinamiento relativamente mayor (Figura 4.4). Dado que,

47
48 CAPÍTULO 5. CONCLUSIONES Y CONSIDERACIONES FINALES

experimentalmente, es el condensado quiral y no la masa constituyente de los

quarks lo que afecta a las observables hadrónicas, cuando se trató de ajustar
el condensado (Cuadro 4.2), la función de masa cae dramáticamente (Figura
4.6), la escala de confinamiento crece demasiado pero el valor de la constante
del decaimiento del pión se vuelve aproximadamente la mitad a la del modelo
MT. La lección fundamental de este ejercicio es que las predicciones del modelo
solo deben considerarse dentro de su rango de validez. Dado que los modelos
potenciales discutidos en este trabajo son válidos a escalas menores que 1GeV ,
un acoplamiento que capture de manera adecuada los efectos infrarojos del mo-
delo y su transición al dominio de libertad asintótica es deseable para extender
el poder predictivo de estos modelos. Los resultados obtenidos en esta tesis
posteriormente serán mostrados en un articulo que se encuentra en preparación.
La motivación para trabajos futuros consiste en una refinación de la fun-
ción de revestimiento dada por la ecuación (4.1). Un primer intento se dió al
agregar una constante de acoplamiento con dependencia del momento dada por
la ecuación (4.7). Aunque la mejora en los resultados no es significativamente
importante, al menos indica que los modelos efectivos del gluón propuestos pue-
den mejorarse y además seguir mostrando confinamiento (Figura 4.12). Además,
se deben renormalizar apropiadamente las ecuaciones dinámicas para ajustar
apropiadamente los parámetros del modelo en el contexto de las ecuaciones de
Schwinger-Dyson. Como idea principal se tiene el resolver las ecuaciones (3.21)
y (3.22) a temperatura finita.
Apéndice A

Lagrangiano de QCD

La densidad Lagrangiana completa para QCD (Cromodinámica Cuántica,

más general para fermiones interactuando dentro del esquema de un grupo no
Abeliano, QCD corresponde cuando el grupo es G = SU (3)), está dado por
1 a µνa
LQCD = − Fµν F + iψ̄k γ µ (Dµ ψ)k − mψ̄k ψk , (A.1)
4
donde se puede poner m = 0 si los fermiones (quarks) son considerados sin
masa. Los elementos que conforman a la ecuación (A.1) son el tensor de fuerza
a
fuerte Fµν dado por
a
Fµν = ∂µ Aaν − ∂ν Aaµ − gf abc Abµ Acν = −Fνµ
a
, (A.2)

donde Aµ es el campo de norma, f abc son constantes reales completamente anti-

simétricas identificadas como las constantes de estructura del grupo de simetría;
Dµ es la derivada covariante que actúa sobre la función de onda como
a a a a
Dµ ψk = ∂µ ψk + igTkl Aµ ψl = (∂µ δkl + igTkl Aµ )ψl . (A.3)

Aquí, los T a ’s representan los generadores del grupo de simetría interno G y se

asumen Hermitianos. Los generadores satisfacen el álgebra de Lie del grupo de
Lie G que tiene la forma
[T a , T b ] = if abc T c , (A.4)
y la identidad de Jacobi asociada con esta álgebra es

[[T a , T b ], T c ] + [[T c , T a ], T b ] + [[T b , T c ], T a ] = 0. (A.5)

De acuerdo al modelo de quarks, cada sabor de quark viene en tres colores –

verde, azul y rojo. Por está razón los operadores Hermitianos son matrices 3 × 3,
las cuales pueden escribirse en términos de 9 números reales a1 , a2 , . . . a8 y θ

H = θ1 + ~λ · ~a, (A.6)

49
50 APÉNDICE A. LAGRANGIANO DE QCD

donde 1 es la matriz unitaria 3×3, θ representa un factor de fase y λ1 , λ2 , . . . , λ8

son las matrices de Gell-Mann denotadas por
     
0 1 0 0 −i 0 1 0 0
λ1 =  1 0 0  λ2 =  i 0 0  λ3 =  0 −1 0 
0 0 0 0 0 0 0 0 0
     
0 0 1 0 0 −i 0 0 0
λ4 =  0 0 0  λ5 =  0 0 0  λ6 =  0 0 1  (A.7)
1 0 0 i 0 0 0 1 0
   
0 0 0 1 0 0
1
λ7 =  0 0 −i  λ8 = √  0 1 0 
0 i 0 3 0 1 −2
y el producto punto denota la suma de 1 a 8
~λ · ~a = λ1 a1 + λ2 a2 + · · · + λ8 a8 . (A.8)
Apéndice B

Programa usado para

resolver las ESD

Para resolver numéricamente las ecuaciones (3.21) y (3.22) se ha usado el

siguiente código en Mathematica 7.0 es

m = 30;
pmin = 10^-5;
pmax = 500;
lambda = Log[10]/m;
bb = 0;
cc = \[Pi];
y = pmin;
nr = IntegerPart[(Log[10, pmax] - Log[10, pmin]) m + 1];
\[CapitalLambda] = 0.234;
w = 0.4;
Dd = 0.93;
mt = 0.5;
\[Tau] = E^2 - 1;
\[Gamma] = 12/(33 - 8);
\[Mu] = 19;
mq = 0;

Do[h = pmin Exp[(j - 1) lambda] (Exp[lambda] - 1); p[j] = N[y];

y = y + h, {j, 1, nr}]

Do[strip[j] = (p[j + 1] - p[j - 1])/2, {j, 2, nr - 1}]

strip[1] = (p[2] - p[1])/2;
strip[nr] = (p[nr] - p[nr - 1])/2;

MT[k_] := ((2 \[Pi])^2 Dd)/w^6 k*Exp[-k/w^2] +

2 (2 \[Pi])^2*\[Gamma]*(1 - Exp[-k/(4 mt^2)])/(k*

51
52 APÉNDICE B. PROGRAMA USADO PARA RESOLVER LAS ESD

Log[\[Tau] + (1 + k/\[CapitalLambda]^2)^2]);(*en está parte se cambia

el modelo que vaya a ser usado*)

ff[k_, p_, x_] := MT[k + p - 2 Sqrt[k p] Cos[x]]

For[i = 1, i <= nr, i++,

For[j = 1, j <= nr, j++,
angularM[i, j] =
NIntegrate[Sin[x]^2 ff[p[i], p[j], x], {x, 0, Pi}];
angularF[i, j] =
NIntegrate[
Sin[x]^2 ff[p[i], p[j],
x] (Sqrt[p[i] p[j]] Cos[x] +
2/(p[i] + p[j] -
2 Sqrt[p[i] p[j]] Cos[x]) (p[i] -
Sqrt[p[i] p[j]] Cos[x]) (Sqrt[p[i] p[j]] Cos[x] -
p[j])), {x, 0, Pi}]]];

Do[angularFren[i] =
NIntegrate[
Sin[x]^2 ff[p[i], \[Mu]^2,
x] (Sqrt[p[i] \[Mu]^2] Cos[x] +
2/(p[i] + \[Mu]^2 -
2 Sqrt[p[i] \[Mu]^2] Cos[x]) (p[i] -
Sqrt[p[i] \[Mu]^2] Cos[x]) (Sqrt[p[i] \[Mu]^2]
Cos[x] - \[Mu]^2)), {x, 0, Pi}];
angularMren[i] =
NIntegrate[Sin[x]^2 ff[p[i], \[Mu]^2, x], {x, 0, Pi}], {i, 1, nr}]

Do[{kerM[i, j] =
1/(2 (\[Pi]^3) ) (p[i]*M[i]*F[i]*F[j])/(
p[i] + M[i]^2) (angularM[i, j] - angularMren[i]),
kerF[i, j] =
1/(6 (\[Pi]^3) ) (p[i]*F[i] F[j])/(
p[i] + M[i]^2) (angularF[i, j]/p[j] -
angularFren[i]/\[Mu]^2)}, {i, 1, nr}, {j, 1, nr}];
Do[{kernintM[j] = mq F[j] + \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(nr\)]\(kerM[i, j]*
strip[i]\)\),
kernintF[j] = 1 - \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(nr\)]\(kerF[i, j]*
strip[i]\)\)}, {j, 1, nr}];

Minic = Interpolation[mini];
Finic = Interpolation[fini];
Ecs = Table[{F[i] == kernintF[i]}, {i, 1, nr}] \[Union]
 53

Table[M[i] == kernintM[i], {i, 1, nr}];

CondIni =
Table[{F[i], 0.6 Finic[p[i]], 1.1 Finic[p[i]]}, {i, 1,
nr}] \[Union]
Table[{M[i], 0.6 Minic[p[i]] , 1.1 Minic[p[i]] }, {i, 1, nr}];

Res = FindRoot[Ecs, CondIni, MaxIterations -> 1500]

ListLogLinearPlot[TablaF, PlotRange -> All]

ListLogLinearPlot[TablaM, PlotRange -> All]
ListLogLinearPlot[TablaZ, PlotRange -> All]

Do[\[Sigma]s[j] = (-M[j] F[j])/(p[j] + M[j]^2) /. Res, {j, nr}];

Tabla\[Sigma]srencb1750 = Table[{p[j], \[Sigma]s[j]}, {j, 1, nr}];
ListLogLinearPlot[Tabla\[Sigma]srencb1750, PlotRange -> All]

Condensado =
Re[(3/(8 \[Pi]) Sum[ \[Sigma]s[i] strip[i], {i, nr}])^(1/3)]

Promedio =
Table[{t,
Log[Abs[1/(2 \[Pi])
Sum[1/Sqrt[p[i]] Cos[t Sqrt[p[i]]] \[Sigma]s[i] strip[i], {i,
nr}]]]}, {t, 1, 20, 0.01}];
ListLogLinearPlot[Promedio,
AxesLabel -> {t, "ln|\[CapitalDelta][t]|"}, PlotRange -> All]

MM = Interpolation[Table[{p[i], -M[i] /. Res}, {i, 1, nr}]]

FF = Interpolation[Table[{p[i], F[i] /. Res}, {i, 1, nr}]]

fpi = (3/(4 \[Pi]^2)

NIntegrate[(
k FF[k] MM[k])/(k + MM[k]^2)^2 (MM[k] - k/2 MM’[k]), {k, p[1],
p[nr]}])^(1/2)
54 APÉNDICE B. PROGRAMA USADO PARA RESOLVER LAS ESD
Bibliografía

[1] V. D. Burkert and T. S. H. Lee , Int. J. Mod. Phys. E 13, 1035 (2004).

[2] S. S. Afonin, Int. J. Mod. Phys. A 22, 4537 (2007).

[3] S. G. Amsler et al., Phys. Lett. B 667, 1 (2008).

[4] S. Capstick and N. Isgur, Phys. Rev. D 34, 2809 (1986).

[5] R. L. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003); R. L. Jaffe,
Phys. Rep. 409, 1 (2005).

[6] E. Santopinto, Phys. Rev. C 72, 022201(R) (2005).

[7] E. Klempt and J. M. Richard, Rev. Mod. Phys. 82, 1095 (2010).

[8] Guy F. de Téramond and Stanley J. Brodsky, Phys. Rev. Lett. 94, 201601
(2005); Nucl. Phys. B, Proc. Suppl. 199, 89 (2010); Stanley J. Brodsky and
Guy F. de Teramond, Phys. Rev. Lett. 96, 201601 (2006); International
School of Subnuclear Physics, 2007 (unpublished), arXiv:0802.0514.

[9] E. Klempt, Eur. Phys. J. A 38, 187 (2008).

[10] Jaemo Park and Piljin Yi, J. High Energy Phys. 06, 011 (2008).

[11] K. Nawa, H. Suganama, and T. Kojo, Phys. Rev. D 79, 026005 (2009).

[12] C.B. Compean and M. Kirchbach, J. Phys, A: Math. Theor 42, 365301
(5pp), (2009).

[13] M. Lüscher and G. Mack, Commun. Math. Phys. 41, 203 (1975).

[14] G.W. Gibbons and A. R. Steit, Phys. Lett. B 346, 255 (1995).

[15] J. L. F. Barbon and C. A. Fuertes, JHEP, 09, 030 (2008).

[16] H. Nicolai, E. Sezgin, and Y. Tanii, Nucl. Phys. B305, 483 (1988).

[17] C.W. Gibbons and P. K. Townsend, Phys. Rev. Lett. 71, 3754 (1993).

[18] M. Kirchbach and C. B. Compean, Phys. Rev. D 82, 034008 (2010).

55
56 BIBLIOGRAFÍA

[19] E. Schrödinger, Proc. R. Ir. Acad. A 46 9-16 (1940).

[20] P. W. Higgs, J. Phys. A: Math. Gen. 12, 309-23 (1979).
[21] C.B. Compean and M. Kirchbach, Eur. Phys. J., A 33 1-4 (2007).
[22] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251, 267-385 (1995).
[23] G. Levai, J. Phys. A: Math. Gen. 22, 689-702 (1989).
[24] R. De, R. Dutt and U. Sukhatme, J. Phys. A: Math. Gen. 25, L843-50
(1992).
[25] V. V. Gritsev and A. Kurochkin Yu, Phys. Rev. B 64, 035308 (2001).
[26] C. B. Compean and M. Kirchbach, J. Phys. A 44, 015304 (2011).
[27] M. S. Bhagwat, I. C. Cloet and C. D. Roberts, (arXiv:0710.2059 [nucl-th]),
in Proceedings of the Workshop on Exclusive Reactions at High Momentum
Transfer, Newport News, Virginia, 21-24 May 2007, Eds. A. Radyushkin
and P. Stoler (World Scientific, Singapore, 2007).
[28] E. Eichten, K. Gottfried, T. Konoshita, K.D. Lane, and T.-M. Yan, Phys.
Rev. D17, 3090 (1978); D21, 203 (1980).
[29] M. Maris, P.C. Tandy, Phys. Rev. C60 055214, (1999).
[30] G.S. Bali, Phys. Rev. 343, 1 (2001).
[31] Arfken, George, Mathematical Methods For Physycists, Third Edition, Aca-
demic Press Inc., 1985, San Diego California 92101, pp 102.
[32] Griffiths, George J. Introduction to Quantum Mechanics, Prentice Hall,
Inc., 1995, Upper Saddle River, NJ 07458, pp 121-140.
[33] Cohen-Tannoudji C., Diu B., Laloe F. Quantum Mechanics vol. 1, pp. 799.
[34] Andreev Oleg, Zakharov Valentin I. Heavy Quark Potentials and AdS/QCD
arXiv:hep-ph/0604204v3 26 26 Jul 2006.
[35] Cohen-Tannoudji C., Diu B., Laloe F. Quantum Mechanics vol. 2, pp. 1096-
1104.
[36] M. Kirchbach, M. Moshinsky, Yu. F. Smirnov, Phys. Rev. D 64, 114005
(2001).
[37] C.B. Compean and M. Kirchbach, J. Phys., A: Math. Gen. 39, 547-57.
(2009).
[38] Ed. por: R. Bijker, O. Castaños, D. Fernández, H. Morales-Técotl, L. Urru-
tia, C. Villareal, Latin American School of Physicis XXXV ELAF, Supersy-
metries in Physics and its Applications, AIP Conference Proceedings 744,
New York, (2005).
BIBLIOGRAFÍA 57

[39] De R., R. Dutt, U. Sukhatme, J. Phys. A: Math, Gen, L843-L850 (1992).

F. Cooper, A. Khare, U.P. Sukhatme, Super Simmetry in Quantum Me-
chanics (World Scientific, Singapore, 2001).

[40] C.B. Compean, M. Kirchbach, J. Phys. A. Math. Gen 38(2005).

[41] P. Dennery and A. Krzywicki, Mathematics for Physicist New York: Dover
(1996).

[42] Abramowitz M. and Stegun I.A., Handbook of Mathematical Functions with

Formulas, Graphs and Mathematical Tables 2nd Ed., New York: Dover
(1972).

[43] [Link]

[44] E. J. Routh, Proc. London Math. Soc. 16, 245 (1884).

[45] V. I. Romanovski, Compt. Rend. Ac. Sci. Paris, 188, 1023 (1929).

[46] H. Garcilazo, Phys. Rev., C 67, 055203, (2003).

[47] Y.S. Kim and M.E. Noz, Theory and Application of the Poincaré Group
(Dordrecht: Reidel), (1986).

[48] A. Pallares-Rivera, M. Kirchbach, J. Phys. A 44, 445302 (2011).

[49] D. Craig Roberts and Anthony G. Williams, [Link].33:477-

575,(1994).

[50] P. W. Higgs, Phys. Lett. 12, 132 (1964).

P. W. Higgs, Phys. Rev. Lett. 13, 508 (1964).
P. W. Higgs, Phys. Rev. 145, 1156 (1966).
T. W. B. Kibble, Phys. Rev. 155, 1554 (1967).

[51] Das Ashok, Lectures of Quantum Field Theory World Scientific Publishing
Company, 2008., pp. 690.

[52] A. Bashir, M.A. Bedolla, C. Compean-Jaso y A. Raya, Trabajo en progreso.

[53] Dyson F.J., Phys. Rev. 75, 1736 (1949).

[54] C.D. Roberts and A.G. Williams, Prog. Part. Nucl. Phys. 33, 477 (1994).

[55] K. Osterwalder and R. Schrader, Commun. Math. Phys. 31, 83 (1973); ibid
Commun. Math. Phys. 42, 281 (1975).

[56] M. Kirchbach, Comunicación privada.

[57] L. C. L. Hollenberg, C.D. Roberts and B.H.J. McKellar, Phys Rev. C 46,
2057 (1992); F.T. Hawes, C.D. Roberts and A.G. Williams, Phys. Rev. D
49, 4683 (1994).
58 BIBLIOGRAFÍA

[58] C.P. Hofmann, A. Raya and S. Sanchez Madrigal, Phys. Rev. D 82, 090611
(2010).
[59] J.C. Taylor, Nucl. Phys. B33, 436 (1971).
[60] V. N. Gribov, Nucl. Phys. B 139 (1978).

[61] D. Zwanziger, Nucl. Phys. B 323 (1989).

[62] P. Maris, A. Raya, C.D. Roberts and S.M. Schmidt, Eur. Phys. J. A 18,
231-235 (2003) y referencias dentro.