Tesis

L. Cancino
2023

Resumen
aaa

Índice
1. Introducción 2

2. Matriz S 2

3. Teorı́a de Regge y observables para el Pomeron/Odderon 5

3.1. Teorı́a de Regge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
3.2. Soft-Pomeron . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.3. Hard (BFKL) Pomeron . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.4. Odderon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

4. Grupo de Renormalización 14
4.1. Ecuación del Grupo de Renormalización Funcional . . . . . . . . . . . . . . . . . . . . . 14

5. Setup 17
5.1. Análisis Dimensional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
5.2. Función/Matriz de Green Libre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
5.3. Regulador y matriz de interacción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

6. Ecuación del Grupo de Renormalización Funcional 20

6.1. Integración en q . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
6.2. Integración en z . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
6.3. Solución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
6.4. Funciones β . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22

7. Dimensión Anómala 23
7.1. Derivada en q . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
7.2. Derivada en z . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
7.3. Integración en q y z . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26

8. Resultados Numéricos 27
8.1. Algoritmos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
8.2. Puntos Fijos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
8.3. Evolución en torno a los puntos fijos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
8.4. Diagonalización de matriz de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36

9. Conclusiones 41

1
A. Funciones β’s 42

1. Introducción
-ℏ=c=1
-Hablar un poco sobre el estudio de las interacciones fuertes:
- Necesidad de una interacción que explique la estabilidad del núcleo del átomo
- Yukawa teorizaba el pión que prevalece en la interacción fuerte (fue el primero en construir una
teorı́a cuantitativa de esta interacción) (1935)
-Antes de QCD se analizaba la matriz S para sacar información de estos procesos sin necesidad de
entender las partı́culas que las contituian. Pudiendo hacer predicciones
Para entender la importancia del Pomeron en la fı́sica, debemos entender el contexto donde fue
teorizada. Por lo que debemos hacer un repaso de conceptos de cálculos de la matriz de Scattering a
altas energı́as y sobre la teorı́a de Regge que la describı́a , esto antes de QCD.
-La teorı́a de Regge es una forma de ver las interacciones fuertes antes de que este fuera visto como
una teorı́a cuántica de campos.
-Hablar sobre el lı́mite de REGGE y los lı́mites IR Y UV.

2. Matriz S
Para entender las transiciones entre un estado de partı́culas (en un tiempo t = tA y 4-momentum
pA ) y otro (en t = tB y 4-momentum pB ), necesitamos un objeto que describa dicha evolución temporal,
la cual podemos representarla en el picture de Schrödinger como

|bf in ; tB ⟩ = exp(−iH(tB − tA ))|ain ; tA ⟩ = S|ain ; tA ⟩,

donde H es el Hamiltoniano y S es la matriz que cumple como mediadora de esta evolución temporal.
Si ambos estados coinciden, S debe ser igual a la matriz unidad, por lo que comúnmente se define
como

S = 1 + iT, donde ⟨bf in ; tB |T |ain ; tA ⟩ = (2π)4 δ 4 (pA − pB )A(a → b),

donde T es la matriz de transferencia que nos muestra las desviaciones del caso trivial y A es la
amplitud de probabilidad para el proceso.
Antes de que se aplicara una teorı́a cuántica de campos para las interacciones fuertes (QCD), los
fı́sicos trabajaban para sacar la mayor cantidad de propiedades mediante ciertas condiciones que debiese
cumplir la amplitud de probabilidad para el sistema, donde las singularidades en esta se interpretan
como partı́culas intercambiándose. La matriz S debe cumplir con ciertas propiedades que representan la
fı́sica del proceso que se esté estudiando, las cuales analizaremos a continuación y derivando ecuaciones
que nos serán útiles en las siguientes secciones.

S debe ser invariante de Lorentz

El que sea invariante de Lorentz implica que debe depender de productos escalares de los 4-
momentum y de las masas de las partı́culas de los estados inicial y final.Si pensamos en una colisión
donde tanto el estado inicial y final contienen 2 partı́culas (ver figura 1.a), podemos realizar un cambio
de variable propuesto por Mandelstan [MandelstamMandelstam1958], donde

s ≡ (p1 + p2 )2 = (p3 + p4 )2 ; t ≡ (p1 − p3 )2 = (p2 − p4 )2 ; u ≡ (p1 − p4 )2 = (p2 − p3 )2 , (1)

donde pi , con i = 1, 2, 3, 4, son los 4-momentum de las partı́culas de masa mi que participan en el
proceso mostrado en la figura 1.a.

2
Figura 1: Se muestra la colisión entre las partı́culas 1 y 2 con resultado de las partı́culas 3 y 4 en el
estado final. En a) el circulo al medio describe todas las posibles interacciones que pueden ocurrir entre
medio. En b) tenemos
√ el mismo proceso, pero hay un intercambio de una partı́cula en el canal t, esto
significa que q = t.

Si analizamos colisiones del tipo mostrado en la figura 1.b, en√ dicho esquema, s es el cuadrado
de la energı́a en el marco de referencia del centro de masa (i.e. s = Ec.m. ); t es el√ cuadrado del
momentum transferido entre partı́culas (i.e. el momentum del estado intermedio q = t) y u es una
variable dependiente que podemos definir a través de la ecuación
s + t + u = m21 + m22 + m23 + m24 ,
la cual se demuestra utilizando (1) y p2i = m2i . Al tener esta ecuación que relaciona 3 variables implica
que se tiene un grado de libertad menos. Ese grado de libertad menos nos permite pensar en la amplitud
como función solo de s y t (i.e. A(s, t)).

S debe ser unitaria

La unitariedad nos indica que se conserva la probabilidad. La conservación de la probabilidad en
una teorı́a cuántica implica que la probabilidad del estado inicial en el tiempo tA debe ser igual al del
estado final en el tiempo tB
⟨bf in ; tB |bf in ; tB ⟩ = ⟨ain ; tA |S † S|ain ; tA ⟩ = ⟨ain ; tA |ain ; tA ⟩,

dicho requerimiento nos permite comprobar que S † S = 1 lo cual es la definición de matriz unitaria.
La unitariedad de la matriz S también nos da unas relaciones entre las amplitudes a través de la
representación de T ,
1 = S † S = (1 − iT † )(1 + iT ) = 1 + iT − iT † + T † T ⇒ i(T † − T ) = T † T
⇒ ⟨bf in ; tB |i(T † − T )|ain ; tA ⟩ = ⟨bf in ; tB |T † T |ain ; tA ⟩
XZ Y
⇒ ⟨ain ; tA |iT † |bf in ; tB ⟩ − ⟨bf in ; tB |iT |ain ; tA ⟩ = d ⟨bf in ; tB |T † |X⟩⟨X|T |ain ; tA ⟩,
X X
P R Q
donde en el último paso utilizamos la relación de completitud 1 = X d X |X⟩⟨X|; X son todos los
Q Q d3 p 1
diferentes estados intermedios que pueden ocurrir y d X = jϵX (2π)j3 2E j
. Escribiendo esta relación
nuevamente, pero mostrando explı́citamente la amplitud, tenemos
i(2π)4 δ 4 (pin − pf in )(A∗ (f in → in) − A(f in → in)) =
XZ Y
d (2π)4 δ 4 (pin − pX )(2π)4 δ 4 (pf in − pX )A(in → X)A∗ (f in → X),
X X

3
la información que nos entregan las δ 4 (pin − pX )δ 4 (pf in − pX ) es la misma que nos entregarı́a δ 4 (pin −
pX )δ 4 (pin − pf in ), por lo que re-ordenamos la ecuación y cancelamos algunos término para tener
XZ Y
A∗ (f in → in) − A(f in → in) = i(2π)4 d δ 4 (pin − pX )A(in → X)A∗ (f in → X),
X X

el cual es el Teorema Óptico Generalizado. Este teorema es una gran consecuencia de la unitariedad
de la matriz S, ya que al cumplirse orden por orden nos relaciona por el lado derecho de la ecuación
términos cuadráticos de la amplitud de procesos intermedios con la parte imaginaria de loops (única
forma en que coincidan con los términos cuadráticos). Cabe recordar que este teorema va mas allá de
los casos perturbativos.
Un caso especial de este teorema es la de un proceso elástico (i.e. las partı́culas del estado inicial y
final son las mismas) y sin cambio de dirección (i.e θ = t = 0), el cual es
XZ Y
2Im[A(s, 0)] = −i(2π)4 d δ 4 (pin − pX )|A(in → X)| ≈ σt (2)
X X
P
donde σt = X σ(in → X) es la suma total de todas las secciones eficaces posibles para el estado
inicial.

S debe ser analı́tica

Que la matriz S sea analı́tica, significa que sus elementos son funciones continuas cuando asumimos
que sus variables cinemáticas son complejas. La excepción a esto lo dan singularidades dictadas por
la unitariedad (lo que permite que sigan siendo analı́ticas, a pesar de las singularidades). Dichas
discontinuidades se asocian a partı́culas intercambiándose en el canal donde exista esta singularidad.
Un ejemplo serı́a la colisión de dos partı́culas en el canal s, donde la partı́cula que intercambian
tiene propagador proporcional a
1
,
s − m + iε
por lo que en nuestra amplitud existirá un polo simple en s = m − iε. Dicho polo, junto al teorema del
residuo, nos dará una contribución a la amplitud en donde se intercambia esta partı́cula en el canal s.
Ademas, una propiedad bastante útil de las amplitudes, es su relación de simetrı́a al intercambiar
sus variables cinemáticas. Si en el canal s para el proceso a + b → c + d tenemos la amplitud A(s, t, u),
para el proceso en el canal t a + c̄ → b̄ + d tendremos la misma amplitud

Aa+c̄→b̄+d (s, t, u) = Aa+b→c+d (t, s, u), (3)

donde ā representa la anti-partı́cula de a.

En la próxima sección veremos la forma que contribuyen ciertas partı́culas compuestas (reggeones)
encontradas en polos del plano complejo usando la teorı́a de Regge y las propiedades que hemos dictado
en esta sección.

4
3. Teorı́a de Regge y observables para el Pomeron/Odderon
Usualmente, tenemos que el numero cuántico l que representa el espı́n (momentum angular) de las
partı́culas es un numero entero o semi-entero. En el caso de la colisión de hadrones (los cuales tienen
espı́n semi-entero), para que se conserven los números cuánticos deben intercambiar partı́culas con
espı́n entero (figura 1.b). En esta sección, veremos como contribuyen a la amplitud el intercambio de
algunas de dichas partı́culas.

3.1. Teorı́a de Regge

Regge [ReggeRegge1959] en 1959, para un modelo cuántico no-relativista, resolvió la ecuación de
Schrödinger bajo un potencial esférico. Para esto, expandió la amplitud en ondas parciales y realizó una
continuación analı́tica de l para valores no solamente no enteros, sino que complejos. Dichas amplitudes
contenı́an polos encontrados en el plano complejo de l (ademas de los polos en el eje real dados por
los valores fı́sicos de l).
Para el caso relativista, se sigue el mismo proceso pero con algunas consideraciones extra. Se
expande la amplitud en ondas parciales

∞  
X 2t
Aac̄→b̄d (s, t) ≈ (2l + 1)al (s)Pl 1 +
s
l=0

donde Pl son los polinomios de Legendre; al la amplitud de las ondas parciales y 1 + 2t

s ≡ θ es el angulo
generado por la colisión. Utilizando la relación dada por (3), tenemos que

∞  
X 2s
Aab→cd (s, t) ≈ (2l + 1)al (t)Pl 1 + .
t
l=0

Sommerfeld [SommerfeldSommerfeld1949], inspirado en Watson [WatsonWatson1918], re-escriben

estas series como integrales complejas con polos en los valores enteros del momentum angular, para
que pueda reproducir las series que tenemos anteriormente

I  
(2l + 1)a(l, t) 2s
A(s, t) ≈ dl P l, 1 + ,
sin(πl) t

donde a(l, t) = al (t) y P (l, 1 + 2s 2s

t ) = Pl 1 + t cuando l es un número entero. Al preguntarse si las
amplitudes a(l, t) son unicas, la manera para que esto sea ası́ es separar las contribuciones de l par e
impar. Por lo que introducimos η = ±1 y re-escribimos la amplitud como
I  
(2l + 1) X 2s
A(s, t) ≈ dl (η + e−iπl )aη (l, t)P l, 1 + .
sin(πl) η=±1 t

Al integrar según la figura 2, vamos a tener una contribución de los polos αiη (t) y otra de una sección
de la integral que no se anula (cuyos lı́mites mostraremos antes de eliminarla), pero que omitiremos a
primer orden.

X X (η + e−iπl )β η (t) 
2s
 Z −1/2+iε
A(s, t) ≈ i
P αiη (t), 1 + + dl...
η=±1 i
sin(παiη (t)) t −1/2−iε

donde βiη (t) son los residuos de los integrandos.

5
Figura 2: Plano complejo de l, donde se incluyen los polos αiη (t) asociados a los espı́n enteros y un
polo adicional complejo, notar que hay un branch point que no se incluye dentro del contorno y que
por ende no contribuye.

Los polinomios de Legendre en el lı́mite s ≫ |t| tienen la forma

2s Γ(2l + 1)  s l
Pl (1 + ) → 2 ,
t Γ (l + 1) 2t
donde Γ son las funciones Gamma. Por lo que finalmente tenemos una expresión para las amplitudes
de los reggeones que se intercambian para la colisión de hadrones en el lı́mite s ≫ |t| (solo considerando
un polo)
A(s, t) ≈ β(t)(η + e−iπα(t) )sα(t) ≈ sα(t) , (4)
donde α(t) es la trayectoria de Regge, esta se relaciona con la masa m y espı́n J de una familia de
partı́culas compuestas (resonancias) que se intercambian en el canal t, y cumple la relación α(m2 ) =
J para cada partı́cula que viva en dicha trayectoria. Cabe mencionar que la expresión (4), es la
contribución de una familia completa de partı́culas con los mismos números cuánticos, donde uno es
el estado basal y los demás son resonancias de el con diferente momentum angular. Las resonancias
que se intercambian en dicho canal son mesones (compuesto de un quark y un anti-quark) o partı́culas
compuestas de gluones.
¿La ecuación (4), cuando consideramos la suma sobre las posibles helicidades de la
colisión nos queda A(s, t) ≈ sα(t)−M , donde M es la helicidad (verificar) del Reggeon, si
M = 1 lo que viene a continuación adquiere sentido, si no es ası́ buscar otra.
Ya que hemos calculado como depende la amplitud de s y t, se puede relacionar con las secciones
eficaces. La sección eficaz diferencial es proporcional a dσ/dt ≈ |A(s, t)|2 , lo que implica que en el caso
de representar un Reggeon tiene la forma
dσ
≈ s2α(t)−2 ,
dt
por otro lado, para calcular la sección eficaz total utilizamos el teorema óptico (2), ası́ que depende de
la parte imaginaria de la amplitud elástica para t = 0 (i.e. σt ≈ Im[Ael (s, 0)]), por lo que
σt ≈ sα(0)−1 . (5)
Si expandimos la trayectoria de Regge en torno a t = 0, para valores pequeños de t tenemos
′
αR (t) = αR (0) + αR t,

6
 ′
donde αR (0) es el intercept y αR el slope del Reggeon. Como se ve en la figura 3, experimentalmente
tenemos que está acorde a los datos experimentales de dichos mesones. Cabe mencionar que en trabajos
de Godizov [GodizovGodizov2022] donde usa un modelo eikonal, utiliza otra representación de la
trayectoria de Regge, considerando ciertos lı́mites que le son útil.

Figura 3: Se muestra la posición de diferentes resonancias en una trayectoria de Regge[].

Donnachie y Landshoff, inspirados en datos experimentales de la época, muestran que la sección

eficaz diferencial de la colisión elastica entre pp y pp̄ a altas energı́as tiene dos regiones en t donde
procesos donde se intercambian cantidades distintas de gluones dominan en cada una.

Figura 4: En a) se muestra el intercambio del Pomeron entre el mar de quarks de cada hadrón que
domina en t pequeños []; en b) está el modelo de 3 gluones mostrado en [Donnachie LandshoffDonnachie
Landshoff1983] que domina para t grandes.

En 1979 [Donnachie LandshoffDonnachie Landshoff1979], calculan el proceso considerando que

cada quark de los hadrones intercambia un gluón entre ellos (como en la Fig. 4.b), mostrando el mismo
comportamiento que tenı́an los datos a valores de −t > 3,5 GeV 2 , lo cual indica que el intercambio de
3 gluones domina en esa región con la forma de la
dσ
≈ t−8 ,
dt
al ver que para la región de t pequeño necesitaban incluir mas gluones, decidieron representarlo por un
único objeto, el Pomeron. Inspirados en las técnicas que se muestrán al principio de la sub-seción del
soft-pomeron (ver Fig. 4.a), en 1983 [Donnachie LandshoffDonnachie Landshoff1983] logran mostrar
para la colisión elástica pp [] y pp̄ [] que dσ/dt tiene un comportamiento acorde a la teorı́a de Regge (4),
el cual se ajusta mejor a los datos para valores pequeños de t, como se muestra en la Fig. 5. Tambien

7
lograron mostrar que en el valle que se muestra en Fig. 5.b domina el intercambio de dos pomeron y
los tres gluones que habiamos mencionado que dominaban en la región de t grande.

Figura 5: Se muestra la sección eficaz diferencial para

√ la colisión pp; en a) se compara
√ el modelo [] que
intercambia un solo Pomeron a bajos t (en grande s = 52,8 GeV y en pequeño s = 550 GeV); en
b) se compara un modelo donde a bajo t domina el intercambio de un Pomeron, en el valle dominan el
intercambio de 2 Pomerones y 3 gluones, y para t grande domina el intercambio de 3 gluones []. Todos
estos modelos fueron comparados con los datos del IRS [] en 1982.

Con estos datos, modelados por el trabajo de Donnachie y Landshoff, podemos notar que en la teorı́a
de Regge pueden haber mas de un Reggeon compuesto de gluones. En las siguientes sub-secciones, nos
enfocaremos mas en el ajuste de los datos, desde los cuales se podrá calcular los intercepts de los
diferentes Reggeones.
Las colisiones de hadrones podemos separarlas en 2 regiones. En primer lugar, tenemos procesos
soft donde la escala de energı́a del momentum transferido t son de cientos de M eV (del orden del
radio del protón); las escalas de energı́a donde dominan los procesos soft los vuelve completamente
no perturbativos. Por otro lado, tenemos procesos hard el cual domina en escalas del momentum
transferido t > 1 GeV , en estos procesos tenemos una parte no perturbativa y otra que si lo es.
A continuación, mostraremos los procesos teóricos y análisis fenomenológicos de los reggeones que
participan en cada tipo de procesos.

3.2. Soft-Pomeron
En 1973, Jaroszkiewicz y Landshoff con la concepción de que los hadrones están compuestos de 3
quarks de valencia - los cuales le entregan sus propiedades (i.e. números cuánticos) - y un mar de pares
quarks-antiquarks proponen un modelo para la colisión protón-protón a altas energı́as y t pequeño, en
el cual el pomeron es intercambiado desde el mar de quarks, tal como se muestra en la Fig. 4.a . En
este trabajo, proponen que la sección eficaz diferencial debe ser
2αP (t)−2
(3βFp (t))4

dσ s
= ,
dt 4π sin2 (0,5παSP (t)) mp
donde mp es la masa del protón; β el acoplamiento entre Pomeron y quarks; y Fp (t) el factor de forma
del protón. En dicho trabajo, comparan su modelo con los datos experimentales de la época, los cuales
fueron tomados en el Intersecting Storage Rings (IRS) por el CERN [Barbiellini] a energı́as de s =
462, 949, 2013 y 2808 GeV 2 ; dicha comparación nos da un valor del slope del Pomeron de
′
αSP = 0,25 ± 0,02 GeV −2 ,

8
donde este valor del slope será utilizado y confirmado por los modelos y datos posteriores para el caso
del soft pomeron que domina en la región de t pequeño.
En 1992, el Particle Data Group propone modelar la sección eficaz de distintas colisiones entre
hadrones con la forma

σ = A + BP n + C(ln(P ))2 + Dln(P ),

En 1992, Donnachie y Landshoff [] mediante Flaminio obtienen los datos del Particle Data Group y
los compara con su propio modelo de Regge

σ = Xsε + Y s−η ,

donde ε = αSP (0) − 1 (como en (5) y η su análogo para los mesones. Ellos logran tener una incerteza
estadı́stica muy cercana a las del otro grupo, aunque con muchos grados de libertad menos y con una
relación fı́sica para cada término (el primero asociado a intercambio de Pomerones y el segundo al
intercambio de los mesones ρ, ω, f y a. Al saber que el Pomeron tiene números cuánticos nulos y tiene
un factor C = 1, se acoplan con la misma intensidad a partı́culas como a anti-partı́culas, teniendo
ası́ menos grados de libertad en este modelo (la mitad que el anterior, ya que solo Y cambia para las
colisiones con anti-partı́culas).
En la Figura 6 se pueden ver los 9 procesos y valores que calculan con los datos, obteniendo que

ε = αSP (0) − 1 = 0,0808,

también se debe notar, que el Kaón y el Pión, al tener un quark menos el pomeron se acopla con un
factor de 2/3 aproximadamente a diferencia de los hadrones.
Existen otros modelos que obtienen resultados diferentes, aunque no demasiado. Por ejemplo, en
1993 [] la colaboración CDF del FERMILAB calculó que ε = αP (0)−1 = 0,112±0,013 y en 1997 Cudell,
Kang y Kim [] filtrando datos (comparando resultados con datos dentro de 1 y 2 σ’s); preocupandose
de minimizar la incerteza estadistica y separando la contribución de los mesones en los que son par
e impar en carga, se obtiene el valor del intercept de 0,096+0,012
−0,009 . En 2022 [], Godizov en un modelo
eikonal que incluye al Hard pomeron obtiene un intercepto de 0,114.

3.3. Hard (BFKL) Pomeron

-
En el nacimiento de QCD, una tarea era encontrar partı́culas con los números cuánticos del vacı́o
(como modelo de Pomeron), ninguna partı́cula fundamental cumplı́a estos requisitos. Durante los 70
(década en la cual nació esta teorı́a) también se estudiaron las singularidades de la amplitud en teorı́as
de Yang-Mills y QED para colisiones a altas energı́as (inspirados en la teorı́a de Regge). En 1975,
Low [] presenta un bag model para el menor orden del Pomeron (lo que denomina bar pomeron) el
cual consiste de dos gluones intercambiándose entre quarks confinados. Esto para conseguir lo que
se observaba de las interacciones hadronicas en esos tiempos, que la sección eficaz total se mantiene
aproximadamente constante entre 10 y 300 GeV de energı́a. Ese mismo año, Nussinov [ ] desarrolla un
método pictórico para crear diagramas de quarks y gluones con color que aporten a la sección eficaz
en colisión de hadrones, esto incluye un apartado especifico para el intercambio de Pomeron. En 1976,
Nussinov [] resume las ideas planteadas por el y Low en un trabajo donde calcula los interceptos de
difrerentes representaciones pictoricas de procesos para la colisión de hadrones.
Siguiendo la metodologı́a desarrollada por Balitski, Fadin, y Lipatov (BFKL) [ - ] donde desarrollan
una ecuación integral para describir la evolución de una escalera de gluones a primer orden en una
aproximación logaritmica en ln(s), donde los gluones verticales están reggeizados y se conectan con los
gluones horizontales a traves del vertice efectivo no local que suma sobre contribuciones de diferentes
diagramas. Para el caso donde esa escalera de gluones fuera representada por un singlete, la ecuación
BFKL a primer orden dada por la figura 7 es

9
Figura 6: Se muestran las secciones eficaces dependiendo de la energı́a de la colisión del protón junto
con kaones, piones, neutrones, fotones y protones (y sus respectivas anti-partı́culas). Este modelo se
compara con datos de diferentes experimentos juntados por el Particle Data Group [] en 1992.

[ω − εG (−k12 ) − εG (−(k1 − q)2 )]f (ω, k1 , k2 , q)

q2 (k1 − q)2 k ′2
Z   
Nc αs 2 ′ 1
2
=δ (k1 − k2 ) − d k − ′ 1+ ′ f (ω, k ′ , k2 , q)
2π (k ′ − q)2 k12 (k − k1 )2 (k − q)2 k12
donde ¿... ? .
Para el caso donde no hay momentum transferido (i.e. t = q = 0), donde los términos ε los pasamos
al lado derecho de la ecuación y les realizamos un cambio de variable, se tiene que la expresión se
reduce a
ωf (ω, k1 , k2 , 0) = δ 2 (k1 − k2 ) + K ⊗ f (ω, k1 , k2 , 0), (6)
donde
d2 k ′ k12
Z  
Nc αs ′
K ⊗ f (ω, k1 , k2 , 0) = δ(k 1 − k ) − f (ω, k1 , k2 , 0). (7)
π2 (k1 − k ′ )2 k ′2 + (k1 − k ′ )2
Se encuentra la solución de la ecuación (6) calculando las auto-funciones de K
K ⊗ ϕa = ωa ϕa , (8)
′ 2 ′
P
que cumplen con la relación de completitud a ϕa (k)ϕa (k ) = δ (k − k ). Calculamos K ⊗ ϕa desde
(7) en coordenadas polares y ocupando el ansatz
1
ϕnν (|k|, θ) = √ (k 2 )−1/2+iν einθ
2π

10
Figura 7: En el lado izquierdo se muestra la contribución en el canal t, donde f corresponde a la suma
de escaleras de gluones. Este orden, representa la ecuación BFKL que se mostrará [].

obtenemos
k12 (k12 )−1/2+iν einθ
Z  
1 Nc αs ′2 1 ′2 −1/2+iν inϕ
K ⊗ ϕnν =√ dk dϕ 2 (k ) e − 2 ,
2π 2π 2 k1 − k ′2 − 2|k1 ||k ′ |cos(χ) k1 + 2k ′2 − 2|k1 ||k ′ |cos(χ)
donde χ ≡ ϕ − θ. Desde (8), despejamos ωa
" −1/2+iν #
k ′2 k12
Z
Nc αs ′2 1 inχ
ωn (ν) = dk dχ 2 e − 2 .
2π 2 k1 − k ′2 − 2|k1 ||k ′ |cos(χ) k12 k1 + 2k ′2 − 2|k1 ||k ′ |cos(χ)

Al resolver la integral angular, nos dan varios términos, si a algunos de ellos le realizamos el cambio
k2 ′2
de variable x ≡ k′21 y a otros x ≡ kk2 tendremos cancelaciones de ellos. Nos queda la solución
1

1
x(|n|+1)/2−iν − 1
Z 
2Nc αs
ωn (ν) = Re dx ,
π 0 1−x
el cual puede ser expresado en término dela función digamma
Z 1
d(ln(Γ(z))) xz−1
γ(z) ≡ = dx ,
dz 0 x−1
lo que nos deja la expresión final
2Nc αs
ωn (ν) = − Re[γ((|n| + 1)/2 − iν) − γ(1)].
π
La función tiene un maximo, dado por el primer orden (i.e. n = ν = 0). Este se relaciona con el
intercepto del Pomeron, ası́ que aproximadamente el intercepto es
2Nc αs 2Nc αs
ω0 (0) = − Re[γ(1/2) − γ(1)] = − (−ln(4)) = αs · 2,6476... ,
π π
donde concluyen que si αs ≈ 0,2 el intercepto debiese ser ω0 (0) ≈ 0,5.
Aunque, en 1986 Lipatov [] encuentra un limite inferior para el intercepto (a primer orden) del bare
pomeron. Dicho lı́mite es
h
αHP (0) − 1 ≳ 0,3,

11
lo cual está en concordancia con el resultado anterior, incluso si permitimos que αs tome otros valores
cercanos. Veamos ahora, que tal se ajustan con los datos experimentales.
Vimos anteriormente que en el lı́mite entre los procesos soft y hard, ademas de la contribución de los
tres gluones, contribuyen dos pomerones. En 1998 [], Donnachie y Landshoff comentan que los nuevos
datos muestran que debiese haber otro pomeron con intercept aproximadamente de 0.4; por lo que a
su parametrización para los procesos soft le agregan un nuevo factor asociado al hard pomeron. Allı́
utilizan dos parametrizaciones que le dan valores de 0,435 y 0,418 dependiendo del modelo. Tambien
argumentan que el slope debe ser mas pequeño para ajustarse a los datos, ellos utilizan en este y futuros
′
trabajos un valor de αHP = 0,1 GeV −2 . En 2001[] repiten el modelo con nuevos datos, obteniendo
valores del intercept de 0,4372 y 0,3936, todos acordes al lı́mite de Lipatov.
En 2013 [], Godizov analiza los datos del proceso γp y encuentra que el intercepto del hard pomeron
h
para ajustarse a la sección eficaz diferencial de αP (0) − 1 ≈ 0,317. En 2022 [], Godizov calcula teorica-
mente la sección eficaz diferencial para una colisión elastica difractiva de nucleones, en un formalismo
de Regge-eikonal. En este√modelo, incluyen a primer orden el soft-pomeron y el meson f reggeizado.
Para valores mayores de s ⩾ 7 T eV (Fig. 8) se ve que el modelo funciona para valores pequeños de
t, por lo que al incluir al hard pomeron con el intercepto calculado en el 2013 se ve en la figura 8 con
la linea punteada que el modelo se ajusta bastante bien con los ultimos datos de ATLAS y TOTEM.

Figura 8: Se ve la comparación entre los datos a diferentes energı́as, donde la linea solida es el modelo
sin hard pomeron y la linea cortada es ya incluyendolo. Datos de ATLAS y TOTEM a diferentes
energias enlos puntos negros.

3.4. Odderon
En 1975; Joynson, Leader, Nicolescu y Lopez [] analizan la discrepancia entre, los datos de colisiones
de hadrones donde existe un intercambio de carga entre ellos tal como π − p → π 0 n y los modelos ya
existentes donde se intercambian dos mesones ρ, donde uno de ellos está reggeizado. Analizan la
contribución de partı́culas con espı́n 1, ven los efectos de incluir una partı́cula que llaman Q y ademas
se incluyó un polo de Regge con intercepto α(0) = 1 el cual debiese estar compuesto de un número
impar de gluones, el cual llamaron odderon.
En 2021, la colaboración entre los experimentos DØ y TOTEM [] publicó el descubrimiento del
intercambio de una partı́cula compuesta de gluones, sin carga de color y C-impar (i.e. se acopla diferente

12
a la colisión pp que a la pp̄ y por ende contribuye negativamente a una de las amplitudes). Un grupo
extrapoló los datos de los experimentos de colisión pp a energı́as de 2.76, 7, 8 y 13 TeV detectados por
TOTEM hasta una energı́a de 1.96 TeV, esto se hace para compararlos (ver figura 9) con los datos de
la colisión pp̄ detectados por DØ que fueron a dicha energı́a. Obtuvieron un error en la extrapolación
de 3.4 σ’s, lo cual lo asocian a que el Odderon al acoplarse con signo diferente a ambos procesos
debe dominar en estas regiones y explica la diferencia de comportamientos entre estas colisiones.
Otro grupo [] analizo los datos a 13 TeV de T OT EM pero a ángulos muy pequeños (donde debiese
dominar el intercambio de reggeones). Encuentran que el radio entre la parte real e imaginaria de la
amplitud decae a altas energı́as a diferencia de los modelos usuales, el decrecimiento de este radio,
argumentan que se debe al intercambio de una partı́cula con C-impar. con una significancia de 3.4-4.6
σ. Combinando ambos datos de los grupos independientes tenemos una significancia total de al menos
5.2 σ, conviertiendose entonces en la primera evidencia experimental de la existencia del Odderon.

Figura 9: En a) se comparan los datos de la extrapolación hecha por TOTEM de la colisión pp (en
rojo) con los medidos en DØ de la colisión pp̄ (en azul). Se ve que en la colisión pp̄ se mantiene apro-
ximadamente constante cuando aumenta |t|, mientras pp tiene un punto silla mı́nimo y otro máximo
para luego decaer. En b) muestran la diferencia entre los puntos silla mı́nimo y máximo R.

Después de la publicación sobre el descubrimiento del Odderon, surgieron artı́culos cuestionando

la metodologı́a [] de esa declaración y contra-respuestas respaldándola[]. Ya que esas discusiones están
fuera de el objetivo de este documento, se le deja al lector referencias[].

13
4. Grupo de Renormalización
A finales de la década de los 30, los fı́sicos de partı́culas ya tenı́an una estructura bastante avanzada
en la teorı́a cuántica de campos, con varios resultados notables. Sin embargo, cuando trataban de hacer
cálculos mas allá de diagramas de árbol (i.e. incluyendo loops) aparecı́an infinitos en observables que
no tenı́an mucho sentido. Ya después de la segunda guerra mundial varios fı́sicos empezaron a trabajar
en encontrar una estructura en aquellos infinitos, desarrollando la idea de contra-términos en teorı́a
de perturbaciones y de como podemos medir diferentes observables efectivos a diferentes escalas de
energı́a [].
En ese mismo periodo de tiempo, las transiciones de orden de segundo grado eran un gran desafió
para la comunidad de materia condensada. Kadanoff [] en 1966 analiza el desarrollo de un modelo
de Ising cerca de la temperatura critica donde ocurre una transición de fase, introduce bloques de un
cierto tamaño para su modelo y muestra como los observables varı́an con el tamaño del bloque. En 1971
formaliza la idea de la dependencia de los observables de una escala, en la cual dejamos fuera modos
que puedan ocurrir a energı́as que no estén dentro de esa escala. La idea principal que ocuparemos en
este trabajo es la de funciones β’s, las cuales son las derivadas de los acoplamientos (λi ’s) con respecto
a dicha escala Λ
∂λ
βi = Λ .
∂Λ
Como vimos en esta pequeña introducción histórica, el grupo de renormalización une dos áreas de
la fı́sica como lo son la materia condensada y la fı́sica de partı́cula. Dichas áreas también comparten
herramientas como la mecánica estadı́stica, en la siguiente sub-sección utilizaremos propiedades de
ella superficialmente, para mas información dejamos algunos reviews [] donde se trata con mayor
profundidad estas técnicas.

4.1. Ecuación del Grupo de Renormalización Funcional

En 1992 Wetterich [], siguiendo trabajos previos, calcula una ecuación de evolución exacta para la
dependencia de la escala de una acción efectiva Γ, una de las caracterı́sticas mas importantes de esta
ecuación del grupo de renormalización es que aplica para regiones no perturbativas, que como vimos
anteriormente es una de las regiones que nos interesa para analizar el intercambio de reggeones. En
esta sub-sección calcularemos dicha ecuación, la que será clave en nuestro trabajo.
Empezamos con una teorı́a con N campos escalares χa , representada por una acción S[χa ] (para
campos escalares reales χa (−q) = χ∗a (q)). A la acción le añadimos un término cinético

dD q
Z
1
∆S = Rk (q)χ∗a (q)χa (q),
2 (2π)D

donde Rk (q) es el regulador, el cual se encarga de eliminar los modos que quedan fuera de la escala
(cutoff ) k. A la acción le agregamos un término de fuente
Z
Sk = S + ∆S − dD xJa (x)χa (x), (9)

por lo que acabamos de definir una acción dependiente de la escala k. Con esto podemos definir un
generador funcional k-dependiente
Z 
Wk [J] = ln(Zk [J]) = ln Dχ exp(−Sk [χ, J]) ,

donde ¿D ...? y Zk la función partición dependiente de k. En el lı́mite k → 0, Wk → W que genera

las funciones de correlación de n puntos para S. Calculamos la derivada funcional de Wk para obtener

14
dichas funciones de correlación
Z Z Z 
δWk 1 δZk [J] 1 δ(−Sk ) 1 δJb (z) b
= = Dχ exp(−Sk [χ, J]) = Dχ dz χ (z) exp(−Sk [χ, J])
δJa (x) Zk [J] δJa (x) Zk [J] δJa (x) Zk [J] δJa (x)
Z Z  Z
1 1
= Dχ dzδ(x − z)δba χb (z) exp(−Sk [χ, J]) = Dχ χa (x) exp(−Sk [χ, J]) ≡ ⟨χa (x)⟩
Zk [J] Zk [J]

encontrando que la primera derivada resulta ser el valor de expectación de los campos χ (notar que
cambiamos las etiquetas de (9), i.e. x, a → z, b). Definimos un campo macroscópico correspondiente al
valor de expectación de χ
δWk
φa (x) ≡ ⟨χa (x)⟩ = (10)
δJa (x)

donde φa (q) ≡ φa (q)[J] (definir la inversa y que también se quitará la dependencia en k de

φ y J). De la misma manera que antes, calculamos la segunda derivada funcional

δ(Zk−1 )
  Z Z  
δ δWk δZk 1 δ(−Sk )
= · b · Dχ χa (x) exp(−Sk [χ, J]) + · Dχ χa (x) exp(−Sk [χ, J])
δJ b (y) δJa (x) δZk δJ (y) Zk δJ b (y)
Z Z Z
1 1
= − 2 · Dχ χb (y) exp(−Sk [χ, J]) · Dχ χa (x) exp(−Sk [χ, J]) + · Dχ χa (x)χb (y) exp(−Sk [χ, J])
Zk Zk

obteniendo la función de correlación de 2 puntos

δ 2 Wk
Gk (x, y)ab = = ⟨χa (x)χb (y)⟩ − ⟨χa (x)⟩⟨χb (y)⟩.
δJa (x)δJ b (y)

Introducimos la acción efectiva Γ̃k [φk ] mediante una transformada de Legendre

Z
Γ̃k [φk ] = dD xJa (x)φak (x) − Wk [J],

tal que

δ Γ̃k
= Ja (x). (11)
δφa (x)

Diferenciamos (10) y (11) para buscar una relación entre ellos

δφa (x) δ 2 Wk
=
δJ b (y) δJ b (y)δJa (x)
δ 2 Γ̃k δJa (x)
=
δφb (y)δφa (x) δφb (y)
2
−1
δ 2 Wk

δ Γ̃k −1
⇒ = = (Gk (x, y)ab )
δφb (y)δφa (x) δJ b (y)δJa (x)
δ 2 Γ̃k
Z
⇒ dD Gk (x, y)ab · c
= δ(x − z)δca
δφ (z)δφb (y)

¿Para valores fijos de φk ?¿(la respuesta podrı́a estar en la relación k-dependiente entre
φk [J] y Jk [φ])?

∂t Γ̃k |φk = k∂k Γ̃k |φk = −∂t Wk |J = ⟨∂t ∆S⟩

15
 Calculamos la última expresión
Z Z
1 D D a 1
⟨∂t ∆S⟩ = ⟨ d xd y χa (x) (∂t Rk (x, y)) χ (y)⟩ = dD xdD y ∂t Rk (x, y)⟨χa (x)χa (y)⟩
2 2
Z Z
1 1
= dD xdD y ∂t Rk (x, y)[G(x, y) + ⟨χ(x)a ⟩⟨χ(y)a ⟩] = dD xdD y ∂t Rk (x, y)[G(x, y) + φa (x)φa (y)],
2 2
por lo tanto
Z Z
∂ ∂ 1 D D
Γ̃k |φk = Γ̃k [φ] = d xd y (∂t Rk (x, y)) G(x, y) + dD xdD y φ(x) (∂t Rk (x, y)) φ(y),
∂t ∂t 2

realizamos una re-definición Γ̃

Γk = Γ̃k − ∆S. (12)

Antes vimos que la función de Green correspondı́a a la inversa de las derivadas funcionales de Γ̃,
por lo que utilizamos (12) para ver su dependencia en Γk

δ2 (2)
G(x, y)−1 = (Γk + ∆S) = Γk + Rk .
δφ(x)δφ(y)

Si diferenciamos, tenemos
1
Z
1 1  −1
(2)
∂t Γk = dD xdD y (∂t Rk (x, y)) G(x, y) = T r[G ∂t Rk ] = T r[ Γk + Rk ∂t Rk ]
2 2 2

1  −1
(2)
∂t Γk = T r[ Γk + Rk ∂t Rk ]. (13)
2
La última expresión es la ecuación de Wetterich, la cual será la pieza principal en este trabajo para
encontrar las funciones β’s de la acción efectiva que proponemos en la siguiente sección.
La traza tiene una representación pictórica (ver figura 10) la cual se relaciona con los diagramas
one particle irreducible 1PI.

Figura 10: Representación pictórica de la traza de (13), donde la cruz es ∂t Rk y el circulo a la derecha
 −1
(2)
es el propagador dependiente de la escala Γk + Rk .

16
5. Setup
En esta sección, definiremos los elementos que se necesitan para poder resolver la ecuación de
Wetterich. Argumentar porque se utiliza esta acción y campos escalares (preguntar por
papers). Primer orden de la Derivative expansion para el término cinético.
La acción efectiva mas general con dos campos ψ y χ es:
Z
1 1
Γ[ψ † , ψ, χ† , χ] = dD xdτ [ZP ( ψ † ∂ˆτ ψ − αP′
ψ † ∇2 ψ) + ZO ( χ† ∂ˆτ χ − αO ′
χ† ∇2 χ) + m(ψχ† + ψ † χ)
2 2
− µp ψ † ψ − µo χ† χ − Vk [ψ † , ψ, χ† , χ; λj ] ], (14)
donde Zi se denomina la función de onda de la renormalización; αi′ es la pendiente de la trayectoria
del Reggeon correspondiente; λj los acoplamientos asociados al potencial y µi las masas de los campos
que se relaciona con el intercepto de la trayectoria de Regge como Z −1 µi = αi (0) − 1. Para este trabajo
consideramos dos campos interactuando (i.e. i = 1, 2) y 6 acoplamientos a tercer orden del potencial
(i.e. j = 1, ..., 6)
Siguiendo los requerimientos mostrados en [Bartels, Contreras VaccaBartels .2017], se expande
polinómicamente el potencial locamente (LPA’) sobre campos pequeños. En este caso consideraremos
las contribuciones cubicas de la expansión polinomial, los términos a orden superior se calculan de la
misma manera, por lo que esos ordenes pueden considerarse en trabajos futuros. Ası́ que el potencial
que utilizaremos es
V [ψ, ψ † , χ, χ† ] =iλ1 ψ † (ψ † + ψ)ψ + iλ2 χ† (χ† + χ)χ + iλ3 ψ † (χ† + χ)ψ (15)
† † † 2 2 † † 2 2 †
+ iλ4 χ (ψ + ψ)χ + iλ5 ((χ ) ψ + χ ψ ) + iλ6 ((ψ ) χ + ψ χ ). (16)
Nos gustarı́a trabajar con la contribución explicita en la acción de la escala k, por lo que realizaremos
un cambio a las variables adimensionales de todos los parámetros de la acción.

5.1. Análisis Dimensional

La dimensión de la acción debe ser nula, por lo que podemos re-definir los campos y las constantes
para trabajar con sus versiones adimensionales. Las variables espacio-temporales tienen dimensión
[x] = k −1 [τ ] = E −1
y ZX es adimensional, con X = P, O.
Como ejemplo, analizamos la dimensión de un solo término de la acción
Z
1
[ dD xdτ ZP ψ † ∂ˆτ ψ] =[xD ][τ ][ψ † ][∂ˆτ ][ψ] = k −D E −1 [ψ † ]E[ψ] = k −D [ψ † ][ψ]
2
⇒[ψ † ] = [ψ] = k D/2 ,
siguiendo la misma técnicas, las dimensiones de los otros parámetros de la acción son
[ψ † ] = [ψ] = [χ† ] = [χ] = k D/2 , [αp′ ] =[αo′ ] = Ek −2 , [m] = [µp ] = [µo ] = E
[λ1 ] = [λ2 ] = [λ3 ] = [λ4 ] =[λ5 ] = [λ6 ] = k −D/2 E
[V ] =k D E.
Con esta información, definimos los parámetros adimensionales como
p V
ψ̃ = Zp k −D/2 ψ , χ̃ = Zo k −D/2 χ , Ṽ = ′ D+2
p
αp k
m µp µo
m̃ = p , µ̃p = ′
, µ̃o =
′
Zp Zo αp k 2 Zp αp k 2 Zo αp′ k 2
λ1 k D/2 λ2 k D/2 λ3,6 k D/2 λ4,5 k D/2
λ̃1 = , λ̃2 = , λ̃3,6 = √ , λ̃ 4,5 =
3/2 3/2 Zp Zo αp′ k 2
p
Zp αp′ k 2 Zo αp′ k 2 Zo Zp αp′ k 2

17
los cuales se substituyen en la acción.
Mas adelante necesitaremos saber las derivadas en función de t = ln(k/k0 ) de los acoplamientos,
las cuales en su versión adimensional y utilizando que ∂t k D = Dk D son

ψ̃ † ψ̃
−(∂t µp )ψ † ψ = −(∂t (Zp k 2 αp′ µ̃p )) = k 2+D αp′ ((ηP + ξP − 2)µ̃p − ∂t µ̃p )ψ̃ † ψ̃ (17)
Zp k −D

−(∂t µo )χ† χ = k 2+D αp′ ((ηo + ξp − 2)µ̃o − ∂t µ̃o )χ̃† χ̃ (18)

ηp ηo
(∂t m)(ψχ† + ψ † χ) = k 2+D αp′ (− − − ξp + 2)m̃ + ∂t m̃(ψ̃ χ̃† + ψ̃ † χ̃) (19)
2 2
3 D
i(∂t λ1 )ψ † (ψ † + ψ)ψ = ik 2+D αp′ ((− ηp − ξp + 2 − )λ̃1 + ∂t λ̃1 )ψ̃ † (ψ̃ † + ψ̃)ψ̃ (20)
2 2
† † 2+D ′ 3 D
i(∂t λ2 )χ (χ + χ)χ = ik αp ((− ηo − ξp + 2 − )λ̃2 + ∂t λ̃2 )χ̃† (χ̃† + χ̃)χ̃ (21)
2 2
2 2+D ′ ηo D
i(∂t λ3,6 )ψ χ = ik αp ((−ηp − − ξp + 2 − )λ̃3,6 + ∂t λ̃3,6 )ψ̃ 2 χ̃ (22)
2 2

ηp D
i(∂t λ4,5 )χ2 ψ = ik 2+D αp′ ((−ηo − − ξp + 2 − )λ̃4,5 + ∂t λ̃4,5 )χ̃2 ψ̃. (23)
2 2

5.2. Función/Matriz de Green Libre

Para el caso libre (sin interacciones), realizamos una transformada de Fourier a la acción, resultando
ser
Z
Γ[ψ̄ † , ψ̄, χ̄† , χ̄] = dD qdω[ψ̄ † ZP (iω + αp′ q 2 )ψ̄ + χ̄† ZO (iω + αo′ q 2 )χ̄ + m(ψ̄ χ̄† + ψ̄ † χ̄) − µp ψ̄ † ψ̄ − µo χ̄† χ̄ ],

donde ψ̄ † , ψ̄, χ̄† , χ̄ son las transformadas de Fourier de sus respectivos campos y re-ordenamos los
(2)
términos de la acción para encontrar la segunda derivada funcional (Γ0 )kl ≡ δΨδΓ 0
k δΨ
T sobre los campos
l

Z
1 (2)
Γ0 = dD qdω ΨT Γ0 Ψ,
2
donde

ΨT = ψ̄ ψ̄ † χ̄ χ̄†

(−iω + αp′ q 2 − k2 αp′ µ̃p ) k2 αp′ m̃

 
0 0
′ 2 2 ′
(2) D (iω + αp q − k αp µ̃p ) 0 k2 αp′ m̃ 0
Γ0 =k  .

0 k2 αp′ m̃ 0 (−iω + αo′ q 2 − k2 αp′ µ̃o )
k2 αp′ m̃ 0 (iω + αo′ q 2 − k2 αp′ µ̃o ) 0

La función de Green libre (la cual utilizaremos mas adelante) resulta ser la inversa de la matriz
que hemos encontrado mas el regulador que viene de la ecuación de Wetterich
1 1 (2)
G0 = (2)
= (2)
Adj[Γ0 + R] (24)
Γ0 +R det[Γ0 + R]

5.3. Regulador y matriz de interacción

En este trabajo se utiliza el regulador de Litim [LitimLitim2001] el cual menciona que es el optimo
para campos bosonicos. Ademas, se ocupa este por coherencia al comparar nuestros resultados con
[ [Bartels, Contreras VaccaBartels .2016], [Bartels, Contreras VaccaBartels .2017]], aunque el uso de

18
otro regulador puede ser interesante para futuras investigaciones. El regulador que utilizamos permite
dejar fuera abruptamente los modos que están fuera de la escala k (existen otros reguladores que
suavizan esta función), lo cual ocurre al integrar la función Heaviside Θ(x). Como (al trabajar con
los campos adimensionales, se hace innecesario agregarles un ZX )

Rp = αp′ (k 2 − q 2 )Θ(k 2 − q 2 )
Ro = αo′ (k 2 − q 2 )Θ(k 2 − q 2 ) = rαp′ (k 2 − q 2 )Θ(k 2 − q 2 )

y sus derivadas

q2
   
1
R˙p = −αp′ (k 2 − q 2 )Θ(k 2 − q 2 ) (ζp + ηp ) + 2k 2 αp′ Θ(k 2 − q 2 ) = 2k 2 αp′ Θ k 2 − q 2


1− 1− (ζp + ηp )
2 k2
q2
   
1
Ṙo = 2rk 2 αp′ Θ k 2 − q 2 1 −


1 − 2 (ζo + ηo ) ,
2 k

las cuales forman las matrices reguladoras

 
0 RP 0 0
RP 0 0 0
R=  0
 (25)
0 0 Ro 
0 0 Ro 0
 
0 Ṙp 0 0
Ṙp 0 0 0
Ṙ = 
0
 (26)
0 0 Ṙo 
0 0 Ṙo 0

En el propagador, tendremos los siguientes términos

αi′ q 2 + Ri = αi′ q 2 + αi′ (k 2 − q 2 )Θ(k 2 − q 2 ) = αi′ k 2 ,

donde este ultimo resultado es valido cuando q 2 < k 2 , lo cual se justificará al integrar sobre q en la
sección 6.1.
El potencial dependiendo solo de parámetros adimensionales es

V [ψ̃, ψ̃ † , χ̃, χ̃† ] =αp′ k 2D [iλ̃1 ψ̃ † (ψ̃ † + ψ̃)ψ̃ + iλ̃2 χ̃† (χ̃† + χ̃)χ̃ + iλ̃3 ψ̃ † (χ̃† + χ̃)ψ̃
+ iλ̃4 χ̃† (ψ̃ † + ψ̃)χ̃ + iλ̃5 ((χ̃† )2 ψ̃ + χ̃2 ψ̃ † ) + iλ̃6 ((ψ̃ † )2 χ̃ + ψ̃ 2 χ̃† )]

con el cual construimos la matriz de interacción, la cual tiene de elementos las segundas derivadas del
potencial con respecto a los campos.

∂ψ̃ψ̃ V ∂ψ̃ψ̃† V ∂ψ̃χ̃ V ∂ψ̃χ̃† V

 
∂ψ̃† ψ̃ V ∂ψ̃† ψ̃† V ∂ψ̃† χ̃ V ∂ψ̃† χ̃† V 
V =
∂ V
.
χ̃ψ̃ ∂χ̃ψ̃† V ∂χ̃χ̃ V ∂χ̃χ̃† V 
∂χ̃† ψ̃ V ∂χ̃† ψ̃† V ∂χ̃† χ̃ V ∂χ̃† χ̃† V

Con todos los elementos necesarios definidos, se prosigue a resolver la ecuación de Wetterich para
la acción (14).

19
6. Ecuación del Grupo de Renormalización Funcional
Durante esta sección, cuando nos referimos a los parámetros y los campos, se estará mencionando a
sus versiones adimensionales, caso contrario a las secciones anterior y posterior. La ecuación del grupo
de renormalización funcional (Ecuación de Wetterich) es:
1
∂t Γ = T r[(Γ(2) + R)−1 Ṙ],
2
(2) R R dqD dω
donde Γ(2) = Γ0 −V y T r[ ... ] ≡ (2π)D (2π)
T r[ ... ]. Realizamos una expansión del tipo (A−B)−1 =
−1 −1 2 −2
−A (1 + BA + B A + ...) en torno a la matriz de interacción
(2)
(Γ0 + R − V )−1 = −G0 (1 + V G0 + V G0 V G0 + V G0 V G0 V G0 + ...).
iω
Realizamos un cambio de variable z ≡ α′p k2 y el propagador nos queda para q 2 < k 2
 
0 (z + r − µo ) 0 −m
(k 2+D αp )−1  0 0 0 0 
G0 =  
(z + 1 − µp ) (z + r − µo ) − m2  0 −m 0 (z + 1 − µp ) 
0 0 0 0
 
0 0 0 0
(k 2+D αp )−1  (−z + r − µo ) 0 −m 0 
+  
(−z + 1 − µp ) (−z + r − µo ) − m2  0 0 0 0 
−m 0 (−z + 1 − µp ) 0
 
1 Gp Gm
= 2+D + ,
k αp Det[+z] Det[−z]

siendo Det[±z] = (±z + 1 − µp ) (±z + r − µo ) − m2 .

6.1. Integración en q
Solo Ṙ contiene una dependencia de q en el caso q 2 < k 2 como se muestra en 25, y todos los términos
en
R Dintegral2 irán multiplicados por alguno de los elementos de R (i.e. solo tendremos integrales del tipo
la
dq (1 + q )...). Por lo que luego de integrar (¿Mostrar Calculo Explicito?) podemos re-definir
dicha matriz para que contenga los elementos que quedaron de la integración en q
 
0 ND AD (ζP , ηP ) 0 0
ND AD (ζP , ηP ) 0 0 0
Ṙ = k 2D αp′ 

, (27)
 0 0 0 rND AD (ζo , ηo )
0 0 rND AD (ζo , ηo ) 0

donde ND = (2√π)D1Γ(D/2) y AD (ζx , ηx ) = D 2

(1 − ζD+2
x ,ηx
).
Ya con la integral en q realizada, el RHS de la ecuación de Wetterich nos queda

dz ′
Z
1 (2) −1 D ′
T r[(Γ + R) Ṙ] = −ik αp T r[ṘG0 (1 + V G0 + V G0 V G0 + V G0 V G0 V G0 + ...)]
2 2π

6.2. Integración en z
Para el caso de la integración en z, la traza nos entregará potencias de z divididas por el determi-
nante de Γ(2) + R. Al re-escribir el determinante, podemos encontrar 4 diferentes polos, de los cuales
solo dos estarán dentro del contorno de integración (ver figura 11). Definimos el residuo de la integral,

20
por lo que para integrar habrı́a que re-emplazar en la traza los términos en z con sus respectivos
residuos.
Si nos damos cuenta, se vuelve un poco complicado ver los polos de Det[±z], por lo que expandire-
mos esa expresión y encontramos los polos utilizando la solución de una función cuadratica. Al final,
tenemos que los polos de Det[z] son

z = x ± y,
p
donde x ≡ (1 + r − µp − µo )/2 e y ≡ (1 − r − µo + µp )2 + (2m)2 /2. Los polos de Det[−z] son los
mismos pero negativo.
Escogemos que los dos polos que estarán dentro del contorno de integración serán los entregados
por Det[+z]. Las integrales a calcular son

dz ′ (z + r − µo )a (−z + r − µo )b (z + 1 − µP )c (−z + 1 − µP )d
Z X
−i r r s s
= Res(r, s, a, b, c, d), (28)
2π (−z + x + y) (−z + x − y) (z + x + y) (z + x − y)
donde la función Res(r,s,a,b,c,d) la definimos como

1 ds−1 (z + r − µo )a (−z + r − µo )b (z + 1 − µP )c (−z + 1 − µP )d

X  
Res(r, s, a, b, c, d) = lı́m
z→−(x+y) s! dz s−1 (−z + x + y)r (−z + x − y)r (z + x − y)s
s−1
(z + r − µo )a (−z + r − µo )b (z + 1 − µP )c (−z + 1 − µP )d
 
1 d
+ lı́m .
z→−(x−y) s! dz s−1 (−z + x + y)r (−z + x − y)r (z + x + y)s
Tomamos la traza que queremos analizar y realizamos las transformaciones que indicamos en esta
sección, agrupamos y encontramos las funciones β de nuestra teorı́a.

Figura 11: Se muestra el plano complejo de z, en rojo los polos de Res[z] y en azul de Res[−z].

6.3. Solución
El LHS de la ecuación de Wetterich para la acción (14) y el tercer orden de la expansión polinomial
local del potencial (15) es

∂t Γ = ṁ(ψχ† + ψ † χ) − µ̇p ψ † ψ − µo χ† χ + iλ̇1 ψ † (ψ † + ψ)ψ + iλ̇2 χ† (χ† + χ)χ + iλ̇3 ψ † (χ† + χ)ψ
+ iλ̇4 χ† (ψ † + ψ)χ + iλ̇5 ((χ† )2 ψ + χ2 ψ † ) + iλ̇6 ((ψ † )2 χ + ψ 2 χ† )],

21
donde cada término lo podemos remplazar por las ultimas ecuaciones (17-23) de la sección de análisis
dimensional. El RHS, hasta tercer orden de la matriz V , nos queda
1  
T r[(Γ(2) + R)−1 Ṙ] = αp′ k D T r[ṘGV GV G] + T r[ṘGV GV GV G] ,
2
donde realizamos el siguiente remplazo desarrollado en la sub-sección anterior

(z + r − µo )a (−z + r − µo )b (z + 1 − µP )c (−z + 1 − µP )d
→ ΣRes[r, s, a, b, c, d].
(−z + x + y)s (−z + x − y)s (z + x + y)r (z + x − y)r

Notamos que las contribuciones de k y αp′ de las matrices Ṙ, G y V se anulan entre si. Factorizando
el resultado final en orden de los campos e igualando los lados izquierdo y derecho de la ecuación en
cada orden, encontramos las funciones β’s de las constantes de acoplamiento.

6.4. Funciones β
Después de agrupar y filtrar nuestros resultados, obtenemos las funciones β’s que se estaban bus-
cando. Por el tamaño de dichas funciones, las presentaremos en una forma general, si forma completa
será mostrada en el apendice (¿?).
Las funciones beta’s tendrán dos partes; una relacionada a un término clásico dado por el análisis
dimensional de las ecuaciones (17-23) que depende de las dimensiones anomalas y el acoplamiento
relacionado a dicha función; en cambio, el otro término corresponderá a correciones radiativas que
dependen de todas las constantes de acoplamiento.
Siendo fk,ij ≡ fk,ij (µp , µo , r, ND , AD (ζp , ηp ), AD (ζo , ηo ))una cierta función, con polos en (1 − µp ),
(r − µo ), (1 + r − µo − µp ) y ((1 − µp )(r − µo ) − m2 ).

X
µ˙p =(−2 + ξp + ηp )µp + λi λj fp,ij
i,j
X
µ˙o =(−2 + ηo + ξp )µo + λi λj fo,ij
i,j
1 X
ṁ =(−2 + ξp + (ηp + ηo ))m + λi λj fm,ij
2 i,j

D 3 X
λ˙1 =(−2 + + ξP + ηp )λ1 + λi λj λk f1,ijk
2 2
i,j,k

D 3 X
λ˙2 =(−2 + + ξP + ηo )λ2 + λi λj λk f2,ijk
2 2
i,j,k

D 1 X
λ˙3 =(−2 + + ξP + ηp + ηo )λ3 + λi λj λk f3,ijk
2 2
i,j,k

D 1 X
λ˙4 =(−2 + + ξP + ηo + ηp )λ4 + λi λj λk f4,ijk
2 2
i,j,k

D 1 X
λ˙5 =(−2 + + ξP + ηo + ηp )λ5 + λi λj λk f5,ijk
2 2
i,j,k

D 1 X
λ˙6 =(−2 + + ξP + ηp + ηo )λ6 + λi λj λk f6,ijk .
2 2
i,j,k

Los términos radiativos tienen una representación pictórica, la cual explicaremos a continuación
para que se entienda cuales de ellos sobrevivirán o no. Como ejemplo, se utilizará la función β de µp , el
cual tiene un factor del potencial de µp ψψ † , los campos representan lineas externas, y sumamos sobre
todos los loops posibles que contengan dos vértices λi λj . Es claro que los vértices tienen que contener
al menos un campo ψ, lo cual descarta λ2 χ† (χ† + χ)χ como posible vértice.

22
7. Dimensión Anómala
En esta sección volvemos a la notación donde los parámetros y campos sin ningún distintivo son
dimensionales. Escogemos la siguiente forma de presentar la acción efectiva
Z
Γ[ψ̄ † , ψ̄, χ̄† , χ̄] = dD qdω[ψ̄ † ZP (iω + αp′ q 2 )ψ̄ + χ̄† ZO (iω + αo′ q 2 )χ̄ + m(ψ̄ χ̄† + ψ̄ † χ̄)

− µp ψ̄ † ψ̄ − µo χ̄† χ̄ + Vk [ψ̄ † , ψ̄, χ̄† , χ̄] ],

donde le realizamos una transformada de Fourier. Calculamos las derivadas funcionales para esta acción

δ2 Γ δ2 V (1,1)
†
= ZP (iω + αp′ q 2 ) − µp + ≡ ΓP (29)
δψδψ δψδψ †
δ2 Γ ′ 2 δ2 V (1,1)
= Z O (iω + αo q ) − µ o + ≡ ΓO .
δχδχ† δχδχ†

Para calcular las dimensiones anómalas, necesitamos las derivadas de ZP y αp′

(1,1) (1,1) (1,1)

∂ΓX ∂ Γ̇X 1 ∂ Γ̇X
lı́m = ZX αx′ ⇒ lı́m = ∂t (ZX αx′ ) ⇒ − ′
lı́m = η x + ξx
ω,q→0 ∂q 2 ω,q→0 ∂q 2 ZX αx ω,q→0 ∂q 2
(1,1) (1,1) (1,1)
∂ΓX ∂ Γ̇X i ∂ Γ̇X
lı́m = iZX ⇒ lı́m = i∂t (ZX ) ⇒ lı́m = ηx ,
ω,q→0 ∂ω ω,q→0 ∂ω ZX ω,q→0 ∂ω
(1,1)
por lo tanto necesitaremos conocer Γ̇X para calcular las dimensiones anómalas.
Para esto utilizamos la ecuación de Wetterich, a la cual le calculamos sus derivadas funcionales
(mostraremos explı́citamente los cálculos para los campos ψ’s, el caso para los χ’s es el mismo y con
igual resultado equivalente a intercambiar los campos)

(1,1) δ2 1 δ2 (2) −1 1 δ2
Γ̇X = ∂t Γ = T r[ Ṙ (Γ + R) ] = T r[ Ṙ (G0 + G0 V G0 + G0 V G0 V G0 + ...)]
δXδX † 2 δXδX
† 2 δXδX † 
(1,1) 1 δ δV δV δV
⇒ Γ̇X = T r[Ṙ G0 G 0 + G 0 G 0 V G 0 + G0 V G 0 G0 + ... ],
2 δX δX † δX † δX †

donde X en las derivadas funcionales representa ψ, χ. El potencial solo tiene términos cúbicos en los
campos, por lo que solo las primeras derivadas en el potencial contribuyen. Ademas, a tercer orden de
interacción hay términos que sobreviven, pero que están asociados a ordenes lineales en los campos, i.e.
el término con el potencial en (29) el cual no contribuye a la dimensión anómala, ya que se eliminará
al final después de derivar por q y ω. Por lo tanto, el resultado final que tendremos que derivar es
 
(1,1) 1 δV δV δV δV 3
Γ̇X = T r[Ṙ G0 G0 G0 + G0 G0 G0 + O(V ) ]. (30)
2 δX † δX δX δX †

Hasta el momento no hemos mostrado la dependencia explicita del propagador, ya que estamos
interesados en sus derivadas sobre q y ω, se introducirán para su análisis. Para hacerlo, hay que recordar
que las derivadas funcionales de la acción efectiva es la función de correlación de 1PI el cual su resultado
se asocia a sumas de diagramas de Feynmann. En nuestro caso, en vez de la acción efectiva tenemos su
derivada con respecto a la escala, aunque la visualización en diagramas de Feynmann seguirá siendo
útil. La Figura 12 nos muestra como cada linea interna (propagador) dependen de los ω’s y q’s; notar
que las variables sin apostrofe corresponden a estados on-shell y con apostrofe son variables of-shell.
Con el análisis anterior, solo nos queda hacer notar que las derivadas que realizaremos es en función
de variables on-shell, en cambio la traza es una suma sobre indices internos (i.e. la integral se realiza

23
Figura 12: Se muestra una representación pictórica de la ecuación (30), en donde el dibujo de abajo
δ2 V
representa un término ṘG0 δXδX † G0 el cual se elimina ya que al tener un potencial a 3 orden, la

segunda derivada funcional la anula (acordarse que la definición de la matriz V ya implica dos derivadas
funcionales).

sobre variables off-shell ). Finalmente, la dimensiones anómalas son

(1,1)  
i ∂ Γ̇X i 1 δV ∂ δV
ηx = lı́m = T r Ṙ G0 (ω ′ , q ′ ) †
(G0 (ω ′ + ω, q ′ − q)) G0 (ω ′ , q ′ )
ZX ω,q→0 ∂ω ZX 2 δX ∂ω δX

δV ∂ δV
+G0 (ω ′ , q ′ ) (G0 (ω ′ − ω, q ′ − q)) G 0 (ω ′ ′
, q )
δX ∂ω δX †
(1,1)  
1 ∂ Γ̇X 1 1 δV ∂ δV
η x + ξx = − lı́m =− T r Ṙ G0 (ω ′ , q ′ ) (G0 (ω ′ + ω, q ′ − q)) G0 (ω ′ , q ′ )
ZX αx′ ω,q→0 ∂q 2 ZX αx′ 2 δX † ∂q 2 δX

′ ′ δV ∂ ′ ′ δV ′ ′
+G0 (ω , q ) (G0 (ω − ω, q − q)) G0 (ω , q ) ,
δX ∂q 2 δX †

por lo que necesitamos conocer las derivadas del propagador para poder obtener las dimensiones
anómalas explı́citamente.
(2)
En G0 la dependencia en ω y q - como se vé en (24) - están dentro de Γ0 , por lo que en los proximos
(2)
calculos, las derivadas las realizaremos explicitamente en Γ0 sobre ω y q. Usaremos la siguiente forma
de la matriz
2
−z − µ̃p + h(q )
 
0 0 m̃
(2) z − µ̃p + h(q 2 ) 0 m̃ 0
Γ0 +R= k D+2 αp′  , (31)

0 m̃ 0 −z − µ̃o + rh(q 2 )
m̃ 0 z − µ̃o + rh(q 2 ) 0

1 1
donde h(q 2 ) ≡ k2 (q
2
+ (k 2 − q 2 )θ(k 2 − q 2 )) = 2
k2 (k θ(k
2
− q 2 ) + q 2 θ(q 2 − k 2 )).

7.1. Derivada en q
Al ver la derivada en q 2 nos damos cuenta que es respecto a un vector, por lo que por simplicidad,
escogemos un marco de referencia donde q = (qx , 0), por lo que el problema se reduce a derivadas de
escalares
∂ 1 ∂2
lı́m 2
= lı́m 2
ω,q→0 ∂qx ω,q→0 2 ∂qx

24
y por ende,
!
∂ ∂ 1 1 (2) 1
G0 = (2)
=− (2)
∂q (Γ0 + R) (2)
∂q ∂q Γ0 +R Γ0 +R Γ0 +R
2
∂ 1 (2) 1 (2) 1 1 (2) 1
G0 = 2 (2) ∂q (Γ0 + R) (2) ∂q (Γ0 + R) (2) − (2) ∂q2 (Γ0 + R) (2) ,
∂q 2 Γ0 + R Γ0 + R Γ0 + R Γ0 + R Γ0 + R

donde las derivadas del inverso del propagador son

 2 
0 ∂q h(q ) 0 0
2
(2) D+2 ′ ∂q h(q ) 0 0 0
∂q (Γ0 + R) = k αp  = k D+2 αp′ h′ O+

0 0 0 r∂q h(q 2 )
0 0 r∂q h(q 2 ) 0

∂q2 h(q 2 )
 
0 0 0
(2) ∂q2 h(q 2 ) 0 0 0
∂q2 (Γ0 + R) = K D+2 αp′  = k D+2 αp′ h′′ O+ ,
 
0 0 0 r∂q2 h(q 2 )
0 0 r∂q2 h(q 2 ) 0

donde
 
0 1 0 0
+
1 0 0 0
O ≡
0
.
0 0 r
0 0 r 0

Reordenando, las derivadas quedan

∂
lı́m G0 (ω + ω ′ , q + q ′ ) = − k D+2 αp′ h′ (q ′ )G0 (ω ′ , q ′ )O+ G0 (ω ′ , q ′ )
ω,q→0∂q
∂2 ′ ′
lı́m G0 (ω + ω ′ , q + q ′ ) =2k 2D+4 αp2 h 2 (q ′ )G0 (ω ′ , q ′ )O+ G0 (ω ′ , q ′ )O+ G0 (ω ′ , q ′ )
ω,q→0 ∂q 2

− k D+2 αp′ h′′ (q ′ )G0 (ω ′ , q ′ )O+ G0 (ω ′ , q ′ )

Discutir sobre h′ y h′′

∂ ′ ′ 1 ∂2 ′ ′
k D+2 αp′ ′′ ′
lı́m G0 (ω + ω , q + q ) = lı́m G 0 (ω + ω , q + q ) = − h (q )G0 (ω ′ , q ′ )O+ G0 (ω ′ , q ′ )
ω,q→0 ∂q 2 2 ω,q→0 ∂q 2 2

7.2. Derivada en z
iω
Realizamos un cambio de variable (z = α′p k2 ) utilizando la regla de la cadena

i
∂ω G0 = (∂z G0 )(∂ω z) = ∂z G0 .
αp′ k 2

Calculamos la derivada en z por separado

1 (2) 1
∂z G0 = − (2)
∂z (Γ0 + R) (2)
,
Γ0 +R Γ0 +R

25
donde
 
0 −1 0 0
(2) 1 0 0 0
∂z (Γ0 + R) = k D+2 αp′   = k D+2 αp′ O− ,

0 0 0 −1
0 0 1 0

por lo tanto,
∂ i ∂
lı́m G0 (ω + ω ′ , q + q ′ ) = ′ 2 lı́m G0 (z + z ′ , q + q ′ ) = −ik D G0 (z ′ , q ′ )O− G0 (z ′ , q ′ )
ω,q→0 ∂ω αp k z,q→0 ∂z

7.3. Integración en q y z
SI h’=0 y ambas trazas son iguales
Finalmente la dimensión anomala es
(1,1)  
i ∂ Γ̇P i δV ∂ δV
ηx = lı́m = T r ṘG0 (ω ′ , q ′ ) (G 0 (ω ′
+ ω, q ′
− q)) G0 (ω ′ ′
, q )
ZX ω,q→0 ∂ω ZX δX † ∂ω δX
(1,1)  
1 ∂ Γ̇X 1 ′ ′ δV ∂ ′ ′ δV ′ ′
η x + ξx = − ′
lı́m =− T r ṘG0 (ω , q ) (G0 (ω + ω, q − q)) G0 (ω , q ) ,
ZX αx ω,q→0 ∂q 2 ZX αx′ δX † ∂q 2 δX

reemplazando por los valores que obtuvimos en las sub-secciones anteriores y mostrando las integrales
explicitamente de la supertraza

−ik D+2 αp′ dq ′D dz

Z Z  
′ ′ δV ′ ′ − ′ ′ δV ′ ′
ηx = T r ṘG0 (ω , q ) † G0 (ω , q )O G0 (ω , q ) G0 (ω , q )
ZX (2π)D 2π δψ δψ
′
−ik D+4 αp2 dq ′D dz ′′ ′
Z Z  
′ ′ δV ′ ′ + ′ ′ δV ′ ′
ηx + ξx = h (q )T r ṘG0 (ω , q ) † G0 (ω , q )O G0 (ω , q ) G0 (ω , q ) .
2ZX αx′ (2π)D 2π δψ δψ

La integral en q ′ se calcula igual que en la sección 6.1, ya que la segunda derivada de h(q 2 ) debe
ser constante, ası́ que el cambio resulta en cambiar la matriz Ṙ por la dada en la ecuación (27). La
integración en z se realiza con las mismas técnicas vistas en la sección 6.2, aunque cambien algunas
matrices en la traza, la estructura seguirá siendo la dada por la ecuación (28).
Como último comentario en esta sección, hacemos notar que la traza nos arroja una constante igual
a (k D+2 αp′ )−1 , por lo que mostraremos la dimensión anómala, considerando que Ṙ está representada
por la matriz anteriormente dicha, y que ademas todas esas matrices de la traza están compuestas
solamente de las matrices (i.e. las constantes se sacarán de la integral) la dimensión anómala es
 
−i
Z
dz ′ ′ δV ′ ′ − ′ ′ δV ′ ′
ηx = T r ṘG0 (ω , q ) † G0 (ω , q )O G0 (ω , q ) G0 (ω , q )
ZX 2π δψ δψ
−ik 2 αp′ h′′
Z  
dz ′ ′ δV ′ ′ + ′ ′ δV ′ ′
η x + ξx = T r ṘG0 (ω , q ) † G0 (ω , q )O G0 (ω , q ) G0 (ω , q ) .
2ZX αx′ 2π δψ δψ

Por tamaño, la forma completa se mostrará en un apéndice.

26
8. Resultados Numéricos
Para encontrar puntos fijos uno de los métodos mas usados es el de Newton o su versión modificada
Newton-Raphson. Estos métodos necesitan estar cercanos al punto fijo para poder encontrarlo con
una gran precisión; al estar trabajando en un espacio 10-dimensional esto se vuelve muy complicado y
termina con muchas divergencias o encontrando el punto fijo Gaussiano (i.e. el origen del sistema de
coordenadas). Por ende, utilizaremos otro método con mejor convergencia, por lo que podremos estar
lejos de los puntos fijos, pero con poca precisión (orden de βi < (10−1 − 10−3 ), para luego utilizar
Newton-Raphson para mejorar la precisión.
Para cualquier método a utilizar, necesitamos condiciones iniciales desde donde nuestro algoritmo
empiece a iterar. Se generan valores aleatorios con la función random de Python para todos los aco-
plamientos, luego se divide en su modulo y obtenemos un punto en la superficie de una esfera unitaria
10-dimensional. Luego, generamos nuevos puntos permutando los signos de cada acoplamiento, lo que
nos genera 210 listas solo con un punto aleatorio; al imponer ligaduras (e.g. al querer masas positivas)
podemos reducir estas listas y concentrarnos en una región del espacio mas acotada.
En la sección 3 se mostraron valores que el intercepto puede tomar, lo cual para el caso del soft
pomeron es mas cercano a 0,1 y el hard pomeron al 0,35. Con esto en mente y revisando valores de
publicaciones anteriores [ [Bartels, Contreras VaccaBartels .2016], [Bartels, Contreras VaccaBartels
.2017]], a las masas las definimos en una esfera con radio 0.1 y los demas acoplamientos con radio 1.
Independiente de las condiciones iniciales, el método numérico que presentaremos a continuación es
muy bueno al encontrar puntos fijos, aún si la condición inicial no está cercana a ellos.

8.1. Algoritmos
El algoritmo con el que empezaremos, es desarrollado por Li y Zeng [Li ZengLi Zeng2008], el cual
muestra buena convergencia. Empezamos con uno de los puntos aleatorios generados anteriormente
(0)
Xi ,
donde i = 1, ..., 10 y representa cada acoplamiento. Luego, reemplazamos ese valor en la funciones β’s
y calculamos el error asociado (sabiendo que queremos que dichas funciones sean 0)
(0) (0) (0)
ei = −βi (X1 , X2 , ..., X10 ),
conociendo el error de cada función definimos la función objetivo
1 2
J= (e + e22 + ... + e210 ),
2 1
esta función objetivo es la que nos dirá en cada iteración si nos estamos acercando o alejando de nuestro
objetivo (i.e. J = 0), por lo tanto debemos elegir una tolerancia para esta función objetivo (T ol) que
nos permita parar la iteración (el cual representa el error acumulado al cuadrado), dicha función se
empieza con T ol = 10−1 y termina con T ol = 10−6 ).
Antes de mostrar la forma en la cual generaremos los nuevos puntos, este algoritmo nos pide definir
una función de optimización
∂β1 2 ∂β2 2 ∂β10 2
 
ηi = (1,0 − 1,6)/ + + ... + ,
Xi Xi Xi
donde el numero 1,0 − 1,6 se debe escoger para cada caso en particular, si su valor es muy grande J
empezará a oscilar y nunca llegará a su objetivo T ol, en cambio si es muy pequeño, su cambio será
muy pequeño y se demorará mucho en llegar al objetivo.
Ya con todas las funciones definidas, la forma en que el algoritmo funciona para entregarnos nuevos
puntos es
 
(k+1) (k) ∂β1 ∂β2 ∂β10
Xi = Xi + ηi ei · + + ... + ,
Xi Xi Xi

27
con estos nuevos puntos, volvemos a evaluar la función error ei y por ende la función objetivo J, si
J > T ol volvemos a iterar hasta que cumpla con nuestro objetivo o si supera una cantidad maxima
de iteraciones paramos el método y desechamos dicho punto para continuar con el otro. En cambio,
si J < T ol se termina la iteración y procedemos a utilizar los puntos encontrados con el método de
Newton-Raphson hasta alcanzar J < 10−14 .Incluir una mención a mathematica que nos entrega
precisiones mucho mayores cuando se está cerca.

8.2. Puntos Fijos

Volver a comparar los puntos fijos con los encontrados en Mathematica.
En un principio se verifica que se encuentran las soluciones mostradas en [Bartels, Contreras Vac-
caBartels .2016] y [Bartels, Contreras VaccaBartels .2017] en algunos puntos fijos. Los puntos fijos
mostrados en dichas publicaciones son calculados en el caso sin dimensión anómala, lo que en nuestro
caso no podemos hacer, ya que, se necesita de una ecuación diferencial asociada a la variable r para
poder ocupar los métodos de convergencia anteriores (e.g. el método de Newton-Raphson necesita de
una matriz invertible, lo que las no cuadradas no cumplen). Por lo que, se mostrarán los resultados
de puntos fijos para las funciones β’s mostradas en [Bartels, Contreras VaccaBartels .2016] en el caso
con dimensión anómala no nula, verificando que los puntos sean soluciones de dichas funciones β’s.
En la publicación [1512.07182], la notación de las constantes de acoplamiento se relacionan con este
trabajo según la siguiente transformación

(µ → µp,o , λ → λ1,2 ) = (0,1111, ±1,0504),

donde si se ocupa el indice izquierdo (derecho) en el intercepto, se debe ocupar el indice izquierdo
(derecho) en λ.
Ver si es necesario calcular el valor con dimensión anómala
En los puntos fijos encontrados tenemos los cuales reproducen estos puntos directamente

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.1333 0 0 -1.1676 0 0 0 0 0 0

como la reproducción de la solución aplicada para dos partı́culas donde ambas partı́culas no in-

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.1333 0.1333 0 -1.1676 -1.1676 0 0 0 0 1

teractuan entre ellas y reproducen por separado lo valores antes encontrados. Ademas, existen varios
puntos fijos muy parecidos, pero donde una de las partı́culas toma mas relevancia, modificando el valor
de r y también de sus acoplamientos. A continuación, mostramos dos ejemplos de ello donde al dividir

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.133333 0.36969 0 -1.16764 -3.23747 0 0 0 0 2.77267
0.133333 0.100925 0 -1.16764 -0.883826 0 0 0 0 0.756936

el intercepto de la partı́cula 2 con el de la partı́cula 1, tenemos el valor de r al menos en 5 decimales;

lo mismo ocurre con las constantes de acoplamiento. Ya que r en un principio se define solo para las
pendientes de la trayectorı́a de Regge, es interesante comprobar que para estos puntos fijos r también
es el radio entre ambas trayectorias. Aunque, se mostraron dos puntos, hubieron una cantidad bastante
considerables de puntos con estas caracteristicas, por lo que se abre una posibilidad de investigar si en

28
verdad es una colección de puntos o existe una trayectorı́a en el espacio de la teorı́a donde las derivadas
de los acoplamientos se anulan.
Para el caso ηx = 0, obtenemos los siguientes puntos, representando dos partı́culas sin interacción

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.105263 0.0123593 0 1.00006 -0.117421 0 0 0 0 0.117414
0.105263 0.0650316 0 1.00006 0.617838 0 0 0 0 0.6178

y una partı́cula sola

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.105263 0 0 1.00006 0 0 0 0 0 0

En el caso de [Bartels, Contreras VaccaBartels .2017], las constantes de acoplamiento en nuestra

notación serı́an Incluyendo la dimensión anómala, se obtuvo (verificar y disminuir el error) que

µp → µp µo → µo λ → λ1 λ2 → λ4 λ3 → λ5 r→r
0.1111 0.110753 1.05034 1.44665 0 0.92181

el punto fijo se deforma, quedando en Comenzando con las ecuaciones que se obtuvieron en la sección

µp → µp µo → µo λ → λ1 λ2 → λ4 λ3 → λ5 r→r
0.13299 0.13058 ± 1.1421 ± 1.55133 0 0.89455

(6.4), se obtienen los siguientes puntos que reproducen estos resultados

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.133333 0.131624 0 -1.16764 0 0 -1.59429 0 0 0.89934

Falta igualarles el error a ambos, son muy diferentes..

Para el caso de la dimensión anómala ηx = 0

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.124938 0.029937 0 ± 1.09179 0 0 ± 0.524002 ± 0.0107924 0 0.0168769
0.105263 0.0304597 0 1.00006 0 0 0.519648 0 0 0
0.113886 0.114224 0 0 1.0852 1.49513 0 0 0 1.08513

Flujos
Ya al haber comparado puntos encontrados por el algoritmo con otras investigaciones, se presentan
los puntos fijos asociados al problema de este trabajo. Aunque, antes queda comentar los filtros que
se usaron para desechar ciertos puntos que se encontraban en el proceso. En primer lugar se elimina
cualquier valor de los interceptos µp y µo negativos y mas pequeños que 10−2 , ya que buscamos valores
α′
que estén en acuerdo con los presentados en las seccion 2. Para r ≡ α′o , no estamos interesados en
p

29
valores de la pendiente αx′ negativas (acorde a la sección 3.), por lo que se filtran valores de r negativos,
ademas de valores mas grandes que 10 para controlar la proporción de las pendientes.
Por otro lado, nos interesa saber que tipo de direcciones tienen los flujos alrededor de los puntos
fijos, ya que esto nos entregará información si son atractores o repulsores en cierta región de energı́as.
Para realizar esto, se lineariza el sistema de ecuaciones diferenciales cerca de los puntos, calculando
la matriz de estabilidad

∂βi
Mij ≡
∂Λj

donde Λj - con i, j = 10 - son todos los parámetros de que dependen las funciones β’s. Los auto-
vectores de dicha matriz nos indican direcciones que definen el flujo en el espacio de la teorı́a, en
torno al punto fijo, el signo de los auto-valores definen si esa dirección es un atractor o repulsor. En el
caso IR (UV), si el auto-valor asociado a la dirección descrita por un auto-vector es positivo (negativo)
estaremos hablando de que esa dirección es atractora, en caso contrario será repulsora. En particular,
estamos interesados en puntos fijos que sean mayormente atractores en el IR, por lo que
se filtran los puntos, quedandose con los que tengan igual o mas cantidad de atractores
que de repulsores. También se encuentran puntos fijos con auto-valores complejos, los cuales para
este trabajo también se omitirán, al representar oscilaciones en torno al punto fijo.
Revisar y editar el filtro a usar, ya que en el Seminario cambió el objetivo y ahora
queremos relacionar puntos en el UV relacionados con los Hard Pomeron con puntos en
el IR relacionados con los Soft Pomeron, conectándolos a través del punto fijo relacionado
a la interacción entre el Soft Pomeron y el Hard Pomeron como estado intermedio.n O
ver si solo ponemos los filtros IR
Teniendo claro los filtros que se ocuparon, se muestra a continuación los puntos encontrados, sin
considerar los ya mostrados que siguen los lı́mites de otras publicaciones.

Atractor en IR

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.13842 0.13842 -0.00564237 -1.01175 -1.01175 -0.815009 -0.815009 0.196741 0.196741 1.
0.13842 0.13842 0.00564237 -1.01175 1.01175 0.815008 -0.815009 0.196741 -0.196741 1.
0.13842 0.13842 0.00564237 1.01175 -1.01175 -0.81501 0.815007 -0.196742 0.196741 1.
0.13842 0.13842 -0.00564237 -1.01175 -1.01175 -0.815008 -0.815009 0.196741 0.196741 1.
0.13842 0.13842 -0.00564237 1.01175 1.01175 0.815009 0.815008 -0.196741 -0.196741 1.

los cuales van asociados a los autovalores

2.49938 -2.22935 -2.08803 -2.03389 1.94173 1.40866 0.361209 -0.106694 0.0979625 0.0309539

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
1.57314 1.57314 1.57314 -1.74009 -0.00129862 -1.74138 -1.74138 -0.00129862 -1.74009 1.
0.552271 0.552271 0.552271 0.0324867 0.737446 -0.309563 -0.309563 0.737446 0.0324868 0.72014
0.766887 0.766887 -0.766887 -1.02402 0.0451113 -0.429861 0.429861 -0.0451113 1.02402 1.38861

con autovalores

30
 -16.4596 8.94694 3.18103 1.22764 0.795258 -0.772361 -0.772361 0.625268 0.625268 -0.535119
-22.1816 9.72513 3.82572 -3.80723 3.80618 -2. -2. 2. 0.802824 9.36178 ·10−8
-22.1816 9.72513 3.82572 -3.80723 3.80618 -2. -2. 2. 0.802824 1.13888 ·10−7

—
Para el caso ηi = 0

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.108359 0.108359 0 0.863122 0.863122 0.696732 0.696732 -0.166389 -0.166389 1.
0.108359 0.108359 0 0.863122 -0.863122 -0.696732 0.696732 -0.166389 0.166389 1.
0.108359 0.108359 0 -0.863122 0.863122 0.696732 -0.696732 0.166389 -0.166389 1.
0.108359 0.108359 0.00423506 0.863122 -0.863122 -0.696732 0.696732 -0.166389 0.166389 1.
0.108359 0.108359 0.00423506 -0.863122 0.863122 0.696732 -0.696732 0.166389 -0.166389 1.
0.108359 0.108359 -0.00423506 -0.863122 -0.863122 -0.696732 -0.696732 0.166389 0.166389 1.

con autovalores

4.25315 3.37136 1.95215 -1.8923 -1.73598 -1.69311 0.886658 0.623412 0.534125 -0.0252709

Atractor en UV
Los puntos asociados a resultados novedosos los dividiremos en sub grupos. Primero tendremos los
que µp = µo = |m|

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.189653 0.189653 0.189653 0.382113 0.382113 0.764225 0.764225 0.382113 0.382113 0.107708
0.2421 0.2421 0.2421 -0.1517 -0.5064 -0.6581 -0.6581 -0.5064 -0.1517 0.133
0.0695444 0.0695444 0.0695444 0.477846 0.477846 0.955692 0.955692 0.477846 0.477846 1.
0.0695444 0.0695444 -0.0695444 -0.477846 0.477846 0.955692 -0.955692 -0.477846 0.477846 1.
0.0695444 0.0695444 -0.0695444 0.477846 -0.477846 -0.955692 0.955692 0.477846 -0.477846 1.

con autovalores

-12.2448 6.41753 -2. -2. -1.83668 -0.73274 -0.73274 -0.73274 0.534519 -0.36637
-10.2679 8.18233 -2. -2. -1.52256 -0.943332 -0.943332 0.324997 -0.146443 0.113336
-2.09068 -2.09023 -2.09023 1.98022 -1.16291 -1.16291 -1.16291 -0.581455 -0.145364 -0.0902275
-2.09068 -2.09023 -2.09023 1.98022 -1.16291 -1.16291 -1.16291 -0.581455 -0.145364 -0.0902275
-2.09068 -2.09023 -2.09023 1.98022 -1.16291 -1.16291 -1.16291 -0.581455 -0.145364 -0.0902275

En segundo lugar tenemos los que solamente µp = µo

donde los dos primeros puntos parecen ser variaciones de [1608.08836]. Los dos últimos puntos son
un ejemplo de una colección de puntos fijos, los cuales tienen ambas masas igual a 0,1333, m = 0 y
r = 1, pero que sus acoplamientos varian un poco; λ1 y λ2 están estrechamente relacionados a dos
decimales y varian entre 1,08 y 1,16; λ3 y λ4 también lo están y fluctuan entre 0,04 y 0,64; finalmente
λ5 y λ6 varian entre 0,02 y 0,32. Tienen autovalores asociados de
y el resto

31
 µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.1333 0.1333 0 -1.1118 -1.1233 -0.4512 0.3186 0.1593 -0.2256 1
0.1333 0.1333 0 1.1563 -1.1574 -0.1963 -0.1677 -0.0839 -0.0982 1

2.45759 2.42011 -2.16667 -2.06395 -2.03451 -0.666667 -0.666667 -0.0429049 -0.0121239 0

2.45759 2.44941 -2.16667 -2.04093 -2.03451 -0.666667 -0.666667 -0.00954336 -0.00252623 0

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.1675 0.4249 0.2668 -0.1769 -0.7147 -0.5636 -0.8975 -0.4488 -0.2818 0.2739
0.1605 2.0132 0.1829 -1.0679 -0.2577 0.5181 2.2472 1.4385 0.1106 0.1582
0.2413 0.2724 0.2564 -0.1191 -0.5540 -0.6172 -0.6558 -0.5214 -0.1265 0.1530
1.1258 0.1437 -0.2387 -0.3262 -0.5141 0.8670 -0.4832 -0.6618 0.1979 1.5614
2.7554 0.09998 0.2811 -1.2174 -0.2203 0.7449 0.3994 -0.1513 -1.3123 1.3307

Caso η = 0
Por ahora colocaremos los puntos fijos con η = 0
Primero colocaremos los puntos fijos que no se le encuentra algún patron entre sus componentes; y
ademas, que fueron encontrados mas de una vez con una precisión de 5 decimales

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.124938 0.029937 0 ± 1.0918 0 0 ±0.524 ±0.0108 0 0.016877
0.262483 0.185318 0.220551 0.166613 0.438263 0.760747 0.639219 0.521586 0.139997 0.092263
0.010466 0.021695 0.01507 0.062694 0.187109 0.180529 0.259918 0.129959 0.09026 0.208445
0.011399 0.02366 0.016423 0.067912 0.203075 0.195681 0.281915 0.140957 0.09784 0.228453
0.016899 0.03462 0.024187 0.09901 0.290283 0.283408 0.405632 0.202816 0.141704 0.346425
0.01831 0.03716 -0.02608 -0.107197 0.310027 0.305457 -0.435201 -0.2176 0.152729 0.375794
0.02321 0.044944 0.0322986 0.136923 0.368921 0.381059 0.530247 0.265123 0.190529 0.472968
0.023557 0.04543 0.0327129 0.139094 0.372475 0.386308 0.536452 0.268226 0.193154 0.479469
1.46371 0.586534 0.926562 0.827321 0.20986 1.04742 0.663042 0.331521 0.523712 0.483232
0.03066 0.05352 -0.04051 -0.18597 0.428904 0.49142 -0.64926 -0.32463 0.245708 0.60243
0.07915 0.05485 0.06589 0.61003 0.35193 1.01567 0.84551 0.42276 0.50783 1.31277

El siguiente grupo de puntos corresponde a puntos donde µp = 0,105263, |λ1 | = 1,00006 y λ5 =

λ6 = 0

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.105263 0 0 1.00006 0 0 0 0 0 0
0.105263 0.0304597 0 1.00006 0 0 0.519648 0 0 0
0.105263 0.0123593 0 1.00006 -0.117421 0 0 0 0 0.117414
0.105263 0.0650316 0 1.00006 0.617838 0 0 0 0 0.6178

donde los últimos dos puntos son ejemplos de un grupo de puntos fijos que el modulo de λ2 coincide
con r.
El siguiente grupo de puntos corresponde a algunos ejemplos de una serie de puntos donde µp =
µo = 0,105263; m = 0 y r = 1.

32
 µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.105263 0.105263 0 -0.704548 -0.060613 -0.121226 -1.4091 0.704548 0.0606128 1
0.105263 0.105263 0 -0.0221863 0.706802 1.4136 -0.0443726 0.0221872 -0.706802 1
0.105263 0.105263 0 0.0429782 0.705843 1.41169 0.0859565 -0.0429783 -0.705843 1
0.105263 0.105263 0 0.809936 -0.640626 0.362456 1.05534 0.527672 0.181228 0.999999

Despues existe un gran grupo de puntos fijos que no aparecen mas de una vez; y que ademas, no
parece haber relación entre sus elementos. Solo se copiará acá un par que son de interes.

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.113886 0.114224 0 0 1.0852 1.49513 0 0 0 1.08513
0 0.016713 0.011648 0.049486 0.146181 0.14201 0.20376 0.10188 0.071004 0.15894
0 0.44553 0.06297 0 2.98537 0.11927 0.843874 0.42194 0.059635 0
0.066544 0 0.01939 -0.60599 -0.01499 -0.353153 -0.102904 -0.051452 -0.176576 0
0.036292 0.014972 0.02331 0.32303 0.085595 0.41496 0.26653 0.133264 0.20748 0
0.214874 0.065072 0.118247 1.46398 0.24398 1.61128 0.886703 0.44335 0.80564 0.031145
0.105928 0.05206 0.07426 0.94527 0.325683 1.32536 0.929135 0.464567 0.662678 9.56084

donde un par de elementos no cumplen con el criterio de que ambas masas deben ser distintas de
cero.

8.3. Evolución en torno a los puntos fijos

-Mostrar gráficos de evolución cerca de esos puntos fijos y terminar con un análisis de
ellos.
Escribir una motivación a analizar estos graficos. Argumentar que nos ayuda a analizar
los flujos y por ende entender la evolución y tratar de proponer conexiones con otros
puntos fijos a traves de alguna trayectoria.
Por ahora, se muestran los graficos para

µp µo m λ1 λ2 λ3 λ4 λ5 λ6 r
0.1384 0.1384 0.0056 -1.0118 1.0118 0.8150 -0.8150 0.1967 -0.1967 1
0,137 0,139 −0,0056424 1,0117501 1,0117497 0,8150081 0,815009 −0,19674 −0,19674 0,99999
0,1375 0,139 −0,0056424 1,0117501 1,0117497 0,8150081 0,815009 −0,19674 −0,19674 0,99999
0,138 0,139 −0,0056424 1,0117501 1,0117497 0,8150081 0,815009 −0,19674 −0,19674 0,99999

33
Running Couplings
Dimensional

34
Adimensional

3D

35
ScatterPlot con trayectoria

8.4. Diagonalización de matriz de masa

En esta sección el enfoque estará en analizar como la matriz de masa termina diagonalizandose (i.e.
ver la evolución con solución m ̸= 0 tornandose a soluciones con m = 0). Para esto nos enfocaremos
en la acción sin interacciones y por el momento sin tomar en cuenta los términos cineticos

Z
Γm [ψ̄ † , ψ̄, χ̄† , χ̄] = dD qdω[m(ψ̄ χ̄† + ψ̄ † χ̄) − µp ψ̄ † ψ̄ − µo χ̄† χ̄ ],

Z
1 (2)
Γm = dD qdω ΨT Γ0 Ψ,
2
donde

ΨT = ψ ψ† χ†


χ
 
0 µp 0 m
(2) µp 0 m 0
Γ0 = .
0 m 0 µo 
m 0 µo 0

Podemos diagonalizar esta matriz de 4x4, aunque primero consideraremos que esta matriz está
compuesta de matrices en bloque, ası́ que primero se diagonalizará la matriz

 
µp m
.
m µo

donde sus auto-valores son

 √ 
µp −µo − (µp −µo )2 +4m2
− 2m 1
 √ 
µp −µo + (µp −µo )2 +4m2
− 2m 1

a los cuales consideramos sus representaciones unitarias, por lo que cada vector lo dividimos por su
modulo

36
 v !2
u p
t1 + 1
u µp − µo ± (µp − µo )2 + 4m2
,
4 m

esto lo hacemos para diagonalizar la matriz, ya que se tiene la relación

O = P −1 M P

donde M es la matriz a diagonalizar; O es la matriz diagonalizada y P es la matriz que contiene a los

auto-vectores.
La matriz O resulta ser

 p 
1 µp + µo − (µp − µo )2 + 4m2
p 0 .
2 0 µp + µo + (µp − µo )2 + 4m2

Analizando el término de la acción dentro de la integral, se puede tambien definir los nuevos campos
despues de la diagonalización

ΨT M Ψ = ΨT P OP −1 Ψ = Ψ′T OΨ′ ,

donde

 
   
 ψ −µ + µ + p(µ − µ )2 + 4m2 + 2mχ 2mχ − ψ µ − µ + p(µ − µ )2 + 4m2 
p o p o p o p o
′
 
Ψ = , 
√ √
s s
  2 2 
 µp −µo − (µp −µo )2 +4m2 µp −µo + (µp −µo )2 +4m2 
m m +4 m m +4

Si realizamos un cambio de variables para ver la contribución de cada campo original en función
de las masas, podemos ver como evoluciona cada contribución al evolucionar en el espacio de la teorı́a.

ψ ′ = C11 ψ + C12 χ
χ′ = C21 ψ + C22 χ,

donde

37
  p 
−µp + µo + (µp − µo )2 + 4m2
C11 = s
√ 2
µp −µo − (µp −µo )2 +4m2
m m +4

2
C12 = s
√ 2
µp −µo − (µp −µo )2 +4m2
m +4
 p 
− µp − µo + (µp − µo )2 + 4m2
C21 = s
√ 2
µp −µo + (µp −µo )2 +4m2
m m +4

2
C22 = s .
√ 2
µp −µo + (µp −µo )2 +4m2
m +4

Se muestran
p los gráficos relacionados a estos parámetros
p adimensionales y dimensionales (con M 1 ≡
µp + µo − (µp − µo )2 + 4m2 y M 2 ≡ µp + µo + (µp − µo )2 + 4m2 ), donde como condición inicial
en t = 0 están los puntos definidos en la sub-sección anterior, y se muestra como evolucionan al IR (t
negativos) y al UV (t positivos).

Analizar el método de diagonalización y ver si dicha evolución muestra el cambio de

un punto fijo a otro.
Flujo de matriz de diagonalización con la condición que en kU V , C11 = 1; C12 = 0;
C21 = 0 y C22 = 1

Adimensional

38
Dimensional

39
¿Caso 4x4?

40
9. Conclusiones
—

41
A. Funciones β’s

µ˙p =λ21 fµ1 + λ1 λ3 fµ2 + λ1 λ5 fµ3 + λ23 fµ4 + λ3 λ5 fµ5 + λ25 fµ6
µ˙o =λ22 fµ7 + λ2 λ4 fµ8 + λ2 λ6 fµ9 + λ24 fµ10 + λ4 λ6 fµ11 + λ26 fµ12
ṁ =λ1 λ2 fm1 + λ1 λ4 fm2 + λ1 λ6 fm3 + λ2 λ3 fm4 + λ2 λ5 fm5 + λ3 λ4 fm6 + λ3 λ6 fm7 + λ4 λ5 fm8 + λ5 λ6 fm
λ˙1 =λ31 f1A + λ21 λ3 f1B + λ21 λ4 f1C + λ21 λ5 f1D + λ21 λ6 f1E + λ1 λ23 f1F + λ1 λ3 λ4 f1G + λ1 λ3 λ5 f1H
+ λ1 λ3 λ6 f1I + λ1 λ4 λ5 f1J + λ1 λ25 f1K + λ1 λ5 λ6 f1L + λ33 f1M + λ23 λ4 f1N + λ23 λ5 f1 + λ23 λ6 f1O
+ λ3 λ4 λ5 f1P + λ3 λ25 f1Q + λ3 λ5 λ6 f1R + λ4 λ25 f1S + λ25 λ6 f1T
λ˙2 =λ32 f2A + λ22 λ3 f2B + λ22 λ4 f2C + λ22 λ5 f2D + λ22 λ6 f2E + λ2 λ3 λ4 f2F + λ2 λ3 λ6 f2G + λ2 λ24 f2H + λ2 λ4 λ5 f2I
+ λ2 λ4 λ6 f2J + λ2 λ5 λ6 f2K + λ2 λ26 f2L + λ3 λ24 f2M + λ3 λ4 λ6 f2N + λ3 λ26 f2 + λ34 f2O + λ24 λ5 f2P + λ24 λ6 f2Q
+ λ4 λ5 λ6 f2R + λ4 λ26 f2S + λ5 λ26 f2T
λ˙3 =λ21 λ2 f3A + λ21 λ3 f3B + λ21 λ4 f3C + λ21 λ5 f3D + λ21 λ6 f3E + λ1 λ2 λ3 f3F + λ1 λ2 λ4 f3G + λ1 λ2 λ5 f3H + λ1 λ2 λ6 f3I
+ λ1 λ23 f3J + λ1 λ3 λ4 f3K + λ1 λ3 λ5 f3L + λ1 λ3 λ6 f3M + λ1 λ24 f3N + λ1 λ4 λ5 f3 + λ1 λ4 λ6 f3O + λ1 λ25 f3P + λ1 λ5 λ6 f3Q
+ λ1 λ26 f3R + λ2 λ23 f3S + λ2 λ3 λ4 f3T + λ2 λ3 λ5 f3U + λ2 λ3 λ6 f3V + λ2 λ4 λ5 f3W + λ2 λ25 f3zzz + λ2 λ5 λ6 f3X + λ33 f3Y + λ23 λ4 f3Z
+ λ23 λ5 f3AA + λ23 λ6 f3AB + λ3 λ24 f3AC + λ3 λ4 λ5 f3AD + λ3 λ4 λ6 f3AE + λ3 λ25 f3AF + λ3 λ5 λ6 f3AG + λ3 λ26 f3AH
+ λ24 λ5 f3AI + λ4 λ25 f3AI + λ4 λ5 λ6 f3AJ + λ35 f3yyy + λ5 λ26 f3AK
λ˙4 =λ1 λ22 f4A + λ1 λ2 λ3 f4B + λ1 λ2 λ4 f4C + λ1 λ2 λ5 f4D + λ1 λ2 λ6 f4E + λ1 λ3 λ4 f4F + λ1 λ3 λ6 f4G + λ1 λ24 f4H + λ1 λ4 λ5 f4I
+ λ1 λ4 λ6 f4J + λ1 λ5 λ6 f4K + λ1 λ26 f4L + λ22 λ3 f4M + λ22 λ4 f4XXX + λ22 λ5 f4Y Y Y + λ22 λ6 f4ZZZ + λ2 λ23 f4N + λ2 λ3 λ4 f4 + λ2 λ3 λ5 f4
+ λ2 λ4 λ5 f4R + λ2 λ4 λ6 f4S + λ2 λ25 f4ZZZ + λ2 λ5 λ6 f4T + λ2 λ26 f4U + λ23 λ4 f4V + λ23 λ6 f4W + λ3 λ24 f4X + λ3 λ4 λ5 f4Y + λ3 λ4 λ6 f4Z
+ λ3 λ5 λ6 f4AA + λ3 λ26 f4AB + λ34 f4AC + λ24 λ5 f4AD + λ24 λ6 f4AE + λ4 λ25 f4AF + λ4 λ5 λ6 f4AG + λ4 λ26 f4AH
+ λ25 λ6 f4AI + λ36 f4AJ
λ˙5 =λ1 λ22 f5A + λ1 λ2 λ3 f5B + λ1 λ2 λ4 f5C + λ1 λ2 λ5 f5D + λ1 λ2 λ6 f5E + λ1 λ3 λ4 f5F + λ1 λ3 λ6 f5G + λ1 λ24 f5H + λ1 λ4 λ5 f5I
+ λ1 λ4 λ6 f5J + λ1 λ5 λ6 f5K + λ22 λ3 f5L + λ22 λ5 f5M + λ2 λ23 f5 + λ2 λ3 λ4 f5O + λ2 λ3 λ5 f5P + λ2 λ3 λ6 f5Q
+ λ2 λ4 λ5 f5S + λ2 λ25 f5T + λ2 λ5 λ6 f5U + λ23 λ4 f5W + λ23 λ6 f5X + λ3 λ24 f5Y + λ3 λ4 λ5 f5Z
+ λ3 λ4 λ6 f5AA + λ3 λ5 λ6 f5AB + λ24 λ5 f5AC + λ4 λ25 f5AD + λ4 λ5 λ6 f5AE + λ25 λ6 f5AF
λ˙6 =λ21 λ2 f6A + λ21 λ4 f6B + λ21 λ6 f6C + λ1 λ2 λ3 f6D + λ1 λ2 λ4 f6E + λ1 λ2 λ5 f6F + λ1 λ2 λ6 f6G + λ1 λ3 λ4 f6H
+ λ1 λ3 λ6 f6J + λ1 λ24 f6K + λ1 λ4 λ5 f6L + λ1 λ4 λ6 f6M + λ1 λ5 λ6 f61 + λ1 λ26 f6O + λ2 λ23 f6P + λ2 λ3 λ4 f6Q
+ λ2 λ3 λ5 f6R + λ2 λ3 λ6 f6S + λ2 λ4 λ5 f6Y Y Y + λ2 λ5 λ6 f6XXX + λ23 λ4 f6U + λ23 λ6 f6W + λ3 λ24 f6X + λ3 λ4 λ5 f6Y + λ3 λ4 λ6 f6Z
+ λ3 λ5 λ6 f6AB + λ3 λ26 f6AC + λ24 λ5 f6AD + λ4 λ5 λ6 f6AF + λ5 λ26 f6AH

42
Referencias
[Bartels, Contreras VaccaBartels .2016] Bartels2015Bartels, J., Contreras, C. Vacca, GP. 2016. Could
reggeon field theory be an effective theory for QCD in the Regge limit? Could reggeon field theory
be an effective theory for QCD in the Regge limit? JHEP03201.

[Bartels, Contreras VaccaBartels .2017] Bartels2016Bartels, J., Contreras, C. Vacca, GP. 2017. Po-
meron - Odderon interactions in a reggeon field theory Pomeron - Odderon interactions in a reggeon
field theory. Phys. Rev. D951014013.
[Donnachie LandshoffDonnachie Landshoff1979] Donnachie1979Donnachie, A. Landshoff, PV. 1979.
Elastic Scattering at Large t Elastic Scattering at Large t. Z. Phys. C255. [Erratum: Z.Phys.C 2,
372 (1979)] 10.1007/BF01546237
[Donnachie LandshoffDonnachie Landshoff1983] Donnachie1983Donnachie, A. Landshoff, PV. 1983.
Multi - Gluon Exchange in pp Elastic Scattering Multi - Gluon Exchange in pp Elastic Scattering.
Phys. Lett. B123345–348. 10.1016/0370-2693(83)91215-7
[GodizovGodizov2022] Godizov2021Godizov, AA. 2022. High-energy elastic diffractive scattering of
nucleons in the framework of the two-Reggeon eikonal approximation (from U-70 to LHC) High-
energy elastic diffractive scattering of nucleons in the framework of the two-Reggeon eikonal appro-
ximation (from U-70 to LHC). Eur. Phys. J. C82156. 10.1140/epjc/s10052-022-10015-7
[Li ZengLi Zeng2008] Li2008Li, G. Zeng, Z. 2008. A neural-network algorithm for solving nonlinear
equation systems A neural-network algorithm for solving nonlinear equation systems. 2008 Inter-
national Conference on Computational Intelligence and Security 2008 international conference on
computational intelligence and security ( 1, 20–23).
[LitimLitim2001] Litim2001Litim, DF. 2001. Optimized renormalization group flows Optimized re-
normalization group flows. Phys. Rev. D64105007. 10.1103/PhysRevD.64.105007

[MandelstamMandelstam1958] Mandelstam1958Mandelstam, S. 1958. Determination of the pion -

nucleon scattering amplitude from dispersion relations and unitarity. General theory Determination
of the pion - nucleon scattering amplitude from dispersion relations and unitarity. General theory.
Phys. Rev.1121344–1360. 10.1103/PhysRev.112.1344
[ReggeRegge1959] Regge1959Regge, T. 1959. Introduction to complex orbital momenta Introduction
to complex orbital momenta. Nuovo Cim.14951. 10.1007/BF02728177
[SommerfeldSommerfeld1949] Sommerfeld1949Sommerfeld, A. 1949. Partial differential equations in
physics Partial differential equations in physics. Academic press.
[WatsonWatson1918] Watson1918Watson, GN. 1918. The diffraction of electric waves by the earth
The diffraction of electric waves by the earth. Proceedings of the Royal Society of London. Series
A, Containing Papers of a Mathematical and Physical Character9566683–99.

43