Geometrı́a de la Información Cuántica

El Espacio de Parámetros y el Tensor Geométrico Cuántico a partir de

integrales de Feynman

Miguel Ángel Flores Morales

Cómputo Cuántico y Álgebras de Clifford: Proyecto Final

Hamiltonianos y Funciones de Onda

Para un sistema cuántico, la forma más inmediata para caracterizarlo es a través de su Lagrangiano o su Hamiltoniano.
Entonces, supongamos que tenemos un Hamiltoniano con un potencial que depende de los parámetros λi ; i = 1, ..., n:

p2
H(p, q, t; {λi }) = + V (q, t; {λi }) (1)
2m
Dichos parámetros son cantidades caracterı́sticas del sistema, y podemos encontrar la dinámica del mismo a partir de
la ecuación de Schödringer:
∂
iℏ Ψ(q, t; {λi }) = Ĥ Ψ(q, t; {λi }) (2)
∂t
Entonces, notemos que la dependencia en los parámetros del Hamiltoniano se heredan a la función de onda (que se
encuentra parametrizada por los mismos valores).
Entonces, para ahorrar notación, omitimos las dependencias de q, t de los estados cuánticos del sistema y nos enfocamos
principalmente en la dependencia de los parámetros. Eso nos resulta interesante porque un cambio en los parámetros nos
resulta en un Hamiltoniano de un sistema del mismo tipo pero que presenta diferencias en las formas del potencial. Esto
quiere decir que podemos interpretar a los estados Ψ({λi }) , Ψ({λi + δλi }) como dos estados que pertenecen a sistemas
semejantes pero que no son el mismo.
Además, ası́ como los estados cuánticos definen un estado de Hilbert y obtenemos la amplitud de probabilidad a partir
de la función de onda que podemos expresar como Ψ(q) = ⟨q|Ψ⟩, y para esta forma, los vectores de estado forman una
variedad de amplitudes de probabilidad; nuestra nueva perspectiva de parámetros resulta el equivalente a cambiar el atlas
coordenado bajo la cual describimos dicha variedad. A esta descripción de la variedad de amplitudes de probabilidad de
estados cuánticos le llamamos el espacio de parámetros.
Ası́ pues, definimos la fidelidad cuántica como la amplitud de transición entre dos estados |Ψ⟩ , |Ψ′ ⟩ y la calculamos de
tal forma que no nos afecten las transformaciones de fase de los estados[1]:

F (Ψ, Ψ′ ) = |⟨Ψ|Ψ′ ⟩| (3)

Aunque ahora trabajaremos principalmente con estados puros, también existe una forma de definir la fidelidad cuántica
para estados mixtos a partir de sus operadores de densidad. Si tenemos dos estados mixtos con operadores de densidad
ρ, ρ′ , entonces la fidelidad cuántica se define como:
np o
F (ρ, ρ′ ) = Tr ρ1/2 ρ′ ρ1/2 (4)

Si tomamos los operadores de densidad de estados puros y aplicamos esta definición, recuperamos la expresión de (3).
Y podemos ver a la fidelidad como la magnitud de un producto interno entre los elementos de un espacio de Hilbert que,
además, resulta ser una medida de la cantidad de información transferida durante un proceso adiabático.
Entonces, podemos obtener mucha información sobre un sistema si aplicamos la fidelidad cuántica en nuestro espacio
de parámetros.

1
La Fidelidad Cuántica en el Espacio de Parámetros
Ahora, en búsqueda para caracterizar la geometrı́a del espacio de parámetros para un sistema, supongamos un sistema
descrito por el hamiltoniano (1) con n = 1; es decir, que el sistema depende únicamente de un parámetro λ, y obtenemos
otro estado introduciendo una ligera variación en el parámetro. Con lo que tenemos los estados: |Ψ(λ)⟩ , |Ψ(λ + δλ)⟩
En nuestro camino por encontrar la fidelidad cuántica, calculamos el traslape entre ambos estados para saber cómo
cambia la información del sistema al variar los parámetros. Si hacemos una expansión en serie de Taylor sobre λ:

∞
∂kΨ
(δλ)k
X 

⟨Ψ(λ)|Ψ(λ + δλ)⟩ = 1 + Ψ
k! ∂λk
k=1
(δλ)2 ∂2Ψ
   
∂Ψ
= 1 + δλ Ψ + Ψ + ...
∂λ 2 ∂λ2
(δλ)2 ∂2Ψ
   
∂Ψ
⇒ ⟨Ψ(λ)|Ψ(λ + δλ)⟩ ≈ 1 + δλ Ψ + Ψ (5)
∂λ 2 ∂λ2

Esto implica que:

2 †
|⟨Ψ(λ)|Ψ(λ + δλ)⟩| = ⟨Ψ(λ)|Ψ(λ + δλ)⟩ ⟨Ψ(λ)|Ψ(λ + δλ)⟩
= ⟨Ψ(λ + δλ)|Ψ(λ)⟩ ⟨Ψ(λ)|Ψ(λ + δλ)⟩
(δλ)2 ∂ 2 Ψ (δλ)2 ∂2Ψ
         
∂Ψ ∂Ψ
≈ 1 + δλ Ψ + Ψ 1 + δλ Ψ + Ψ
∂λ 2 ∂λ2 ∂λ 2 ∂λ2
        2   2 
∂Ψ ∂Ψ 2 ∂Ψ ∂Ψ 1 ∂ Ψ ∂ Ψ
= 1 + δλ Ψ + Ψ + (δλ ) Ψ Ψ + Ψ + Ψ
∂λ ∂λ ∂λ ∂λ 2 ∂λ2 ∂λ2
    2   2 
∂Ψ ∂Ψ 1 ∂ Ψ ∂ Ψ
= 1 + (δλ2 ) Ψ Ψ + Ψ + Ψ
∂λ ∂λ 2 ∂λ2 ∂λ2
    
∂Ψ ∂Ψ ∂Ψ ∂Ψ
= 1 − (δλ)2 − Ψ Ψ (6)
∂λ ∂λ ∂λ ∂λ

Dado a que este es el cálculo de la fidelidad al cuadrado, factorizamos la raı́z cuadrada y despreciando términos
superiores al orden cuadrático en δλ, con lo que nos queda:

(δλ)2
    
∂Ψ ∂Ψ ∂Ψ ∂Ψ
F (Ψ(λ), Ψ(λ + δλ)) = 1 − − Ψ Ψ (7)
2 ∂λ ∂λ ∂λ ∂λ

Aquı́ hicimos una aproximación a segundo término, por lo que si ahora tuviéramos nuevamente un sistema que depende
de un conjunto de n parámetros {λi }, usando la convención de la suma de Einstein, la ecuación (7) resulta:
1
F (Ψ(λ), Ψ(λ + δλ)) = 1 − gij δλi δλj (8)
2
Donde:       
1 ∂Ψ ∂Ψ ∂Ψ ∂Ψ ∂Ψ ∂Ψ
gij = + − Ψ Ψ (9)
2 ∂λi ∂λj ∂λj ∂λi ∂λi ∂λj

Entonces, llamamos a gij la métrica de la información cuántica, pues es el tensor métrico que nos va a inducir una
forma de medir distancias en el espacio de parámetros.

2
El Tensor Geométrico Cuántico
Aunque ya tenemos una forma de calcular que tan ”[Link] una función de onda de otra en el espacio de parámetros.
A pesar de esto, tenemos un problema único en las variedades de funciones y nuestra forma de medir distancias entre
ellas: las fases relativas.
Supongamos que tenemos dos estados |ψ⟩ , |ϕ⟩ los cuales son linealmente dependientes y entre ellos existe una fase
relativa eiα ; es decir:

|ψ⟩ = eiα |ϕ⟩ (10)

En tal caso, notemos que los estados viven en el mismo espacio de Hilbert y, escencialmente, son el mismo estado pues
dado un operador hermitiano Â:

⟨ψ|ψ⟩ = ⟨ϕ|e−iα eiα |ϕ⟩ = ⟨ϕ|ϕ⟩ = 1 (11)

−iα iα
⇒ ⟨ψ|Â|ψ⟩ = ⟨ϕ|e Âe |ϕ⟩ = ⟨ϕ|Â|ϕ⟩ (12)

A pesar de que, basándonos en las ecuaciones previas, podemos decir que tenemos los mismos estados, realmente hay
una diferencia crucial entre ellos: la interferencia cuántica.
Si tenemos un sistema compuesto por la superposición de estos dos estados, entonces el término de interferencia i es
tal que

I ∝ 4 cos2 (α/2)
Entonces ahora sı́ tenemos un problema porque tenemos estados que sı́ pueden diferenciarse en el sentido fı́sico pero
que nuestra métrica de la información cuántica va a tratar de forma indistinta. Para facilitar el desarrollo, trabajemos
únicamente con el estado base del sistema |0(t)⟩ tal que Ĥ(γ(t)) |0(γ(t))⟩ = E(γ(t)) |0(γ(t))⟩,(aunque es posible extender
la deducción para estados excitados).
Para esto tomemos una curva cerrada dependiente del tiempo en la N-variedad definida por los parámetros (no confundir
con el espacio de parámetros) γ(t); γ(0) = γ(T ) y ahora denotamos:

Ĥ(t) := Ĥ(γ(t)) (13)

E(t) := E(γ(t)) (14)
|ψ(t)⟩ = |0 ◦ γ(t)⟩ := |0(γt )⟩ (15)

Ahora bien, suponiendo que el Hamiltoniano es dependiente del tiempo, que nos encontramos en un proceso adiabático,
entonces [2]:
i
Rt
|ψ(t)⟩ = e− ℏ 0
E(τ )dτ iφ(t)
e |0(γt )⟩ (16)
A partir de donde se deduce en el apéndice A.2 que:

φ̇ = i 0(γt ) 0̇(γt ) (17)

Este ’ajuste’ que toma en cuenta la fase relativa puede verse como una 1-forma A en el espacio de parámetros. Si ahora
denotamos al estado base como |Ψ⟩, entonces:
 
∂Ψ
A = i ⟨Ψ|dΨ⟩ ; A = i Ψ dλk = Ak dλk (18)
∂λk
Y a partir de esto, nos damos cuenta que iAk , de hecho, es una conexión en el espacio de parámetros. Es por eso que
definimos la 2-forma:
   
∂Ψ ∂Ψ ∂Ψ ∂Ψ
F = dA = i ⟨dΨ|∧|dΨ⟩ ; F = i − dλi ∧ dλj = Fij dλi ∧ dλj (19)
∂λi ∂λj ∂λj ∂λi
Debido a que esta 2-forma nos caracteriza la curvatura propia del espacio de parámetros a partir de las tansformaciones
de norma es que llamamos a Ak la fase de Berry y a Fij la curvatura de Berry.

3
 Ahora que ya tenemos todo para caracterizar las distancias en el espacio de parámetros (y parece que no se nos
escapa nada) podemos definir un tensor que describe toda la geometrı́a del espacio de parámetros asociados a un sistema.
Naturalmente, llamamos a este objeto el Tensor Geométrico Cuántico y sus componentes se definen como sigue:
    
∂Ψ ∂Ψ ∂Ψ ∂Ψ
Gij = − Ψ Ψ (20)
∂λi ∂λj ∂λi ∂λj
Este es un objeto crucial en la Geometrı́a de Información Cuántica, pues lleva codificada toda la geometrı́a del espacio
de parámetros. Resulta ser que la métrica de la información cuántica y la curvatura de Berry son proporcionales a la parte
simétrica y antisimétrica del tensor respectivamente:

gij = Re{Gij } (21)

Fij = −2 Im{Gij } (22)

Por lo que también podemos expresar las componentes del tensor geométrico cuántico como:
1
Gij = gij + Fij (23)
2i

Tensor Geométrico Cuántico con integrales de Feynman

Presentaremos una forma de incorporar un método perturbativo que nos facilitará la forma de calcular las componentes
del Tensor Geométrico Cuántico a costa de tener que usar métodos aproximados y renunciar a obtener soluciones exactas.
Sea un sistema para el cual se tiene un Lagrangiano L0 que depende de N parámetros {λi }N i=1 y en un tiempo t = 0 se
modifica adiabáticamnete el sistema por lo que a partir de ese tiempo el Lagrangiano del sistema es L′ = L0 + Oi (t)δλi .
Y sean H0 , H ′ las transformadas de Legendre de L, L′ ; es decir, sus Hamiltonianos asociados, y tomamos los eigenestados
de H, H ′ como |n⟩ , |n′ ⟩ respectivamente.
Ahora supongamos que sabemos que la partı́cula se encuentra en un estado inicial |q⟩ a un tiempo t0 < 0 y en un
estado final |q ′ ⟩ en el momento t′ > 0. Buscamos saber cuál es la amplitud de transición de ese estado al otro, y notemos
que ya no estamos trabajando con las funciones de onda de los sistemas, pero ya tomamos en cuenta el cambio en la
descripción de ambos sistemas a partir de los parámetros. Entonces podemos expresar dicha amplitud como:

i ′ ′ i
⟨q ′ , t′ |q0 , t0 ⟩ = ⟨q ′ | e ℏ H t e− ℏ H0 t0 |q0 ⟩
X i ′ ′ i
= ⟨q ′ | e− ℏ H t |n′ ⟩ ⟨n′ | e ℏ H0 t0 |q0 ⟩
n′
X i ′ ′ i
= ⟨q ′ | e− ℏ H t
|n′ ⟩ ⟨n′ |n⟩ ⟨n| e ℏ H0 t0 |q0 ⟩
n,n′
X i ′ i
= ⟨q ′ | e− ℏ En′ t |n′ ⟩ ⟨n′ |n⟩ ⟨n| e ℏ En t0 |q0 ⟩ (24)
n,n′

Ahora vamos a realizar una rotación de Wick tal que t → −iτ y nos pasaremos a unidades naturales (ℏ = 1) para
facilitar el cálculo. Entonces:
X ′
⟨q ′ , τ ′ |q0 , τ0 ⟩ = ⟨q ′ | e−En′ τ |n′ ⟩ ⟨n′ |n⟩ ⟨n| eEn τ0 |q0 ⟩ (25)
n,n′

Y notemos que si nos situamos mucho antes de la ’interacción’ (τ0 → −∞) y mucho después de ésta (τ ′ → ∞) entonces
las exponenciales de la energı́a se harán menores, haciendo que predomine por mucho la energı́a más pequeña; es decir,
prácticamente nos quedaremos únicamente con el estado base, cuyas energı́as denotamos E, E ′ . Entonces:

⟨q ′ , ∞|q0 , −∞⟩ = ⟨q ′ , ∞|Ψ′ ⟩ ⟨Ψ′ |Ψ⟩ ⟨Ψ|q0 , −∞⟩ (26)

4
 Aquı́ volvimos a aplicar el operador de evolución temporal en los estados iniciales y finales para hacer tender τ ′ , τ0 a
infinito y menos infinito respectivamente. Si despejamos de la ecuación anterior:

⟨q ′ , ∞|q0 , −∞⟩
⟨Ψ′ |Ψ⟩ = (27)
⟨q , ∞|Ψ′ ⟩ ⟨Ψ|q0 , −∞⟩
′

Y hacemos uso de la integral de Feynman para expresar el numerador y denominador de la ecuación (27). Recordemos
que:
Z qN  Z tN 
⟨qN , tN |q0 , t0 ⟩ = Dq exp i dtL[q] (28)
q0 t0

Entonces, aplicando esto tenemos:

Z q′  Z ∞ Z 0 
⟨q ′ , ∞|q0 , −∞⟩ = Dq exp − dτ L′ − dτ L0
q0 0 −∞
Z  Z ∞ Z 0 
i
 
= Dq exp − dτ L0 + Oi (τ )δλ − dτ L0
0 −∞
Z  Z ∞ Z ∞ 
= Dq exp − dτ L − dτ Oi (τ )δλi (29)
−∞ 0

Por otro lado, tenemos que podemos definir a partir de la integral de Feynman una función de partición temporal para
el sistema que va como:
Z  Z ∞ 
Z[L] = Dq exp − dτ L[q] (30)
−∞

Entonces, abusando un poco de la notación, podemos escribir lo siguiente:

Z  Z ∞ 
p
⟨q ′ , ∞|Ψ′ ⟩ = Dq exp − dτ L′ = Z[L′ ] (31)
0
Z  Z
p 0 
⟨Ψ|q0 , −∞⟩ = Dq exp − dτ L0 = Z[L0 ] (32)
−∞
′ ′
p
⇒ ⟨q , ∞|Ψ ⟩ ⟨Ψ|q0 , −∞⟩ = Z[L′ ]Z[L0 ] (33)

Y, a su vez, el valor esperado de un operador para nuestro sistema puede definirse a partir de las integrales de camino
como:
Z  Z ∞ 
1
⟨X⟩ = Dq exp − dτ L0 X(q) (34)
Z[L0 ] −∞

Con esto podemos expresar la ecuación (29) como:

 Z  ∞ 
′ i
⟨q , ∞|q0 , −∞⟩ = Z[L0 ] exp − dτ Oi (τ )δλ (35)
0

Con esto podemos escribir el traslape en la ecuación (27) de la siguiente manera:

 R∞  R∞
′ exp − 0 dτ Oi (τ )δλi exp − 0 dτ Oi (τ )δλi
⟨Ψ |Ψ⟩ = q ′
= rR
Z[L ] R∞
Dq exp{− −∞ dτ L′ }
Z[L0 ]
Z[L0 ]
 R∞  R∞
exp − 0 dτ Oi (τ )δλi exp − 0 dτ Oi (τ )δλi
= rR R∞
=q  R∞ (36)
Dq exp{− −∞ dτ L′ } exp − 0 dτ Oi (τ )δλi
Z[L0 ]

5
 2
Ahora, en un procedimiento análogo al cálculo de la fidelidad cuántica en (6) calculamos ⟨Ψ′ |Ψ⟩ , donde los cálculos
con más detalle se realizan en el apéndice A.4:

 R∞ D n R oE
i 0 j
2
exp − 0
dτ O i (τ )δλ exp − −∞
dτ O j (τ )δλ
⟨Ψ′ |Ψ⟩ = D n R
∞
oE
exp − −∞ dτ Ok (τ )δλk
R∞ R∞ D R0 R0 E
1 − 0 dτ Oi (τ )δλi + 1/2 0 dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj 1 − −∞ dτ Oj (τ )δλj + 1/2 −∞ dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
≈ D R∞ R∞ E
1 − −∞ dτ Ok (τ )δλk + 1/2 −∞ dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
1 ∞
 Z 0 Z
1
=1+ dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩
2 −∞ 2 0
Z Z ∞Z 0 
1
− dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ + dτ1 dτ2 ⟨Oi (τ1 )⟩ ⟨Oj (τ2 )⟩ δλi δλj
2 0 −∞
Z ∞Z 0
=1− dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ − ⟨Oi (τ1 )⟩ ⟨Oj (τ2 )⟩ δλi δλj (37)
0 −∞

Como aquı́ tomamos en cuenta la fase al no tomar el módulo del traslape al cuadrado, por lo que no estamos deduciendo
únicamente la métrica. Más aún:
1
⟨Ψ′ |Ψ⟩ = 1 − Gij δλi δλj (38)
2
Por lo que al juntar las ecuaciones (37) y (38) encontramos una expresión para calcular las componentes del Tensor
Geométrico Cuántico deducida a partir de integrales de Feynman:
Z 0 Z ∞
Gij = dτ1 dτ2 [⟨Oi (τ1 )Oj (τ2 )⟩ − ⟨Oi (τ1 )⟩ ⟨Oj (τ2 )⟩] (39)
−∞ 0

Notemos que, efectivamente, en ninguna parte de la deducción requerimos de conocer la forma explı́cita de la función de
onda del sistema; y, aunque los valores de expectación sı́ se encuentran determinados por las funciones de onda, podemos
encontrar una proporcionalidad entre dichos valores de expectación y las funciones de Green de interacción del sistema.
Sobre todo, cuando los coeficientes se pueden poner en monomios de q(τ ) podemos asegurar que:

⟨q m (τ1 )q n (τ2 )⟩ = G′m+n (τ1m , τ2n ) (40)

Donde G′n+m (τ1m , τ2n ) es la función de Green de interacción del sistema de grado n + m. Entonces, podemos facilitar
aún más el cálculo haciendo teorı́a de perturbaciones sobre el sistema para poner las funciones de Green de interacción
en términos de perturbaciones de la función de Green de un sistema con un Lagrangiano base que puede ser tan simple
como queramos (por ejemplo, una partı́cula libre o un oscilador armónico).

Conclusión
El Tensor Geométrico Cuántico es una herramienta muy útil para caracterizar la información de un sistema y para
poder encontrar transiciones de fase del mismo; e incorporar la Geometrı́a Diferencial y el Cálculo de Variedades de una
manera más formal a la Teorı́a de Información Cuántica nos permite obtener una nueva perspectiva de las propiedades
cuánticas de nuestros objetos de interés.

Además, las integrales de camino de Feynman nos aportan una muy poderosa herramienta para realizar cálculos
aproximados y facilitar la obtención de incógnitas claves para describir sistemas de información cuántica.

6
1. Referencias

[1] Álvarez Jimenez, J. (2019). Elementos de la Geometrı́a del Espacio de Parámetros Cuántico. [Tesis de Doctorado],
UNAM

[2] Bustos Juárez, S. (2019). Estensión del Tensor Geométrico Cuántico bajo el Formalismo de Integrales de Trayectoria.
[Tesis de Licenciatura], UNAM

[3] Zweibach, B. (2022). Mastering Quantum Mechanics: essentials, theory and applications. Imprenta del MIT, Cam-
bridge.

[4] MacKenzie, R. (2000). Path Integral Methods and Applications. arXiv: Quantum Physics.
Disponible en [Link]

[5] Austrich-Olivares, J. A., & Vergara, J. D. (2022). The Quantum Geometric Tensor in a Parameter-Dependent
Curved Space. Entropy, 24(9), 1236. MDPI AG.
Disponible en [Link]

7
A. Apéndice: Cálculos importantes para el desarrollo
A.1. Derivadas del traslape
Recordemos que por las condiciones de normalización de las funciones de onda:
   
∂ ∂Ψ ∂Ψ
⟨Ψ|Ψ⟩ = 1 ⇒ ⟨Ψ|Ψ⟩ = Ψ + Ψ =0 (41)
∂λ ∂λ ∂λ

Ası́ pues, si calculamos la segunda derivada respecto al parámetro:

∂2
 
∂ ∂
⟨Ψ|Ψ⟩ = ⟨Ψ|Ψ⟩
∂λ2 ∂λ ∂λ
   
∂ ∂Ψ ∂Ψ
= Ψ + Ψ
∂λ ∂λ ∂λ
 2
∂2Ψ
    
∂ Ψ ∂Ψ ∂Ψ
= Ψ +2 + Ψ =0 (42)
∂λ2 ∂λ ∂λ ∂λ2

Y la relación (42) implica directamente que:

 2
∂2Ψ
    
1 ∂ Ψ ∂Ψ ∂Ψ
Ψ + Ψ =− (43)
2 ∂λ2 ∂λ2 ∂λ ∂λ

A.2. Cálculo de la fase temporal de Berry

Derivando la ecuación (16) tenemos que:
 
∂ − ℏi 0t E(τ )dτ iφ(t)
R i − ℏi 0t E(τ )dτ iφ(t)
R i
Rt
|ψ⟩ = iφ̇e e − E(t)e e |0(γt )⟩ + e− ℏ 0 E(τ )dτ eiφ(t) 0̇(γt ) (44)
∂t ℏ
Y a partir de la ecuación de Schrödringer tenemos que:

∂ i
|ψ⟩ = − Ĥ(t) |ψ⟩
∂t ℏ
i i Rt
= − e ℏ 0 E(τ )dτ eiφ(t) Ĥ(t) |0(γt )⟩
ℏ
i i
Rt
= − E(t)e ℏ 0 E(τ )dτ eiφ(t) |0(γt )⟩ (45)
ℏ
Usando este desarrollo en (44), nos queda:
i
Rt i
Rt
iφ̇e− ℏ 0
E(τ )dτ iφ(t)
e |0(γt )⟩ + e− ℏ 0
E(τ )dτ iφ(t)
e 0̇(γt ) = 0 (46)
i
Rt
E(τ )dτ −iφ(t)
Y si multiplicamos toda la ecuación por ⟨ψ| = e ℏ 0 e ⟨0(γt )|, tenemos finalmente que:

iφ̇ = − 0(γt ) 0̇(γt )

∴ φ̇ = i 0(γt ) 0̇(γt ) (47)

8
A.3. Deducción de la curvatura de Berry
Partimos de la 1-forma definida en (18) con lo que, por definición y usando que el producto cuña es antisimétrico,
tenemos que:
 
∂Ai i j ∂Aj j i ∂Ai ∂Aj
dA = dλ ∧ dλ + dλ ∧ dλ = − dλi ∧ dλj (48)
∂λj ∂λi ∂λj ∂λi
      
∂Ai ∂Aj ∂ ∂Ψ ∂ ∂Ψ
− = i j Ψ −i i Ψ
∂λj ∂λi ∂λ ∂λi ∂λ ∂λj
2
∂2Ψ
       
∂Ψ ∂Ψ ∂ Ψ ∂Ψ ∂Ψ
=i + Ψ − − Ψ (49)
∂λj ∂λi ∂λj ∂λi ∂λi ∂λj ∂λi ∂λj
Cancelando términos y recuperando la ecuación (48) tenemos que:
   
∂Ψ ∂Ψ ∂Ψ ∂Ψ
F = dA = i − dλi ∧ dλj (50)
∂λj ∂λi ∂λi ∂λj

A.4. Cálculo del cuadrado del traslape con integrales de Feynmann

Recuperamos a partir de la segunda deducción de la ecuación (37):

R∞ R∞ D R0 R0 E
2
1− 0
dτ Oi (τ )δλi + 1/2 0
dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj 1 − −∞ dτ Oj (τ )δλj + 1/2 −∞ dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
⟨Ψ′ |Ψ⟩ ≈ D R∞ R∞ E
1 − −∞ dτ Ok (τ )δλk + 1/2 −∞ dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
(51)

DEntonces, suponiendo δλk ≪ 1 podemos aplicar la linealidad

E del valor de expectación y hacer una serie de Taylor para
R∞ k
R∞ i j
− −∞ dτ Ok (τ )δλ − /2 −∞ dτ1 dτ2 Oi (τ1 )δλ Oj (τ2 )δλ . Entonces:
1

1
D R∞ R∞ E
1− −∞
dτ Ok (τ )δλk + 1/2 −∞
dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
1
= DR
∞ R∞ E
1 − −∞ dτ Ok − (τ )δλk
−∞
dτ1 dτ2 Oi
1/2
(τ1 )δλi O (τ )δλj
j 2
Z ∞ Z ∞ 
=1+ dτ Ok (τ )δλk − 1/2 dτ1 dτ2 Oi (τ1 )δλi Oj (τ2 )δλj
−∞ −∞
Z ∞ Z ∞ 2
k i j
+ dτ Ok (τ )δλ −
dτ1 dτ2 Oi (τ1 )δλ Oj (τ2 )δλ
1/2 + ...
−∞ −∞
∞ Z Z ∞   
k i j
≈1+ dτ Ok (τ )δλ − /2 1 dτ1 dτ2 Oi (τ1 )δλ Oj (τ2 )δλ
−∞ −∞
Z ∞
1 ∞
Z
=1+ dτ ⟨Ok (τ )⟩ δλk − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ δλi δλj
−∞ 2 −∞
Ahora bien, haciendo el producto de la ecuación (51) tenemos que:

∞
1 ∞
 Z Z 
′ 2 i i j
⟨Ψ |Ψ⟩ = 1 − dτ ⟨Oi (τ )⟩ δλ + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ δλ δλ ·
0 2 0
Z 0
1 0
 Z 
1− dτ ⟨Oj (τ )⟩ δλj + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ δλi δλj ·
−∞ 2 −∞
Z ∞
1 ∞
 Z 
1+ dτ ⟨Ok (τ )⟩ δλk − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ δλi δλj (52)
−∞ 2 −∞

9
 Ahora bien, notemos que todos los términos lineales se van a anular porque tenemos:

Z ∞ Z ∞ Z 0 Z ∞ Z ∞
dτ ⟨Ok (τ )⟩ δλk − dτ ⟨Oj (τ )⟩ δλj − dτ ⟨Oi (τ )⟩ δλi = dτ ⟨Ok (τ )⟩ δλk − dτ ⟨Ok (τ )⟩ δλk = 0
−∞ 0 −∞ −∞ −∞

Y si truncamos el producto hasta términos cuadráticos nos resulta que:

1 ∞
 Z 0 Z
′ 2 1
⟨Ψ |Ψ⟩ = 1 + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩
2 −∞ 2 0
Z Z ∞Z 0 
1
− dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ + dτ1 dτ2 ⟨Oi (τ1 )⟩ ⟨Oj (τ2 )⟩ δλi δλj (53)
2 0 −∞

Y usamos una identidad del cálculo multivariable:

Sea una función ξ(x, y) integrable en todo el plano. Entonces:
Z ∞ Z ∞ Z ∞ Z ∞ Z 0 Z 0 Z ∞ Z 0 Z 0 Z ∞
dxdy ξ = dxdy ξ + dxdy ξ + dxdy ξ + dxdy ξ
−∞ −∞ 0 0 −∞ −∞ 0 −∞ −∞ 0

Entonces, usando que trabajamos con una teorı́a con simetrı́a bajo inversión temporal y tomando los primeros tres
términos entre corchetes de (53) :
Z 0 Z ∞ Z
1 1 1
dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ + dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ =
2 −∞ 2 0 2
Z ∞Z 0 Z 0 Z ∞ Z 0 Z ∞
1 1
− dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ = − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩
2 0 −∞ 2 −∞ 0 −∞ 0

Con lo que concluimos que la ecuación (53) resulta en:

Z ∞ Z 0
′ 2
⟨Ψ |Ψ⟩ = 1 − dτ1 dτ2 ⟨Oi (τ1 )Oj (τ2 )⟩ − ⟨Oi (τ1 )⟩ ⟨Oj (τ2 )⟩ δλi δλj (54)
0 −∞

10