Escuela Politécnica Nacional Mecánica Cuántica

CRISTALES DE TIEMPO
Evolución de los Cristales de tiempo: de la paradoja teórica a la realidad
cuántica

Kevin Alexander Muñoz Linares

Alexander Geovanny Chanatasig Vargas

20/07/2026

1. Introducción: Ruptura Espontánea de la Simetrı́a

de Traslación Temporal
En el marco formal de la mecánica cuántica y la teorı́a de campos, la conservación
de la energı́a es una consecuencia directa de la simetrı́a de traslación temporal continua
(T → T + ϵ). Un cristal espacial ordinario se define por la ruptura espontánea de la
simetrı́a de traslación espacial. En 2012, Shapere y Wilczek plantearon la posibilidad de
extender este fenómeno a la dimensión temporal [1, 2].
Clásicamente, las ecuaciones de Hamilton dictaminan que la energı́a E(p, q) = H(p, q)
se minimiza en un estado estático donde los gradientes se anulan (ṗj = q̇ j = 0). Para
evadir esta restricción y permitir movimiento en el estado fundamental, Shapere y Wilczek
demostraron que se requieren Lagrangianos no convencionales (ej. L = − κ2 ϕ̇2 + λ4 ϕ̇4 ), los
cuales generan Hamiltonianos multivaluados con ’cúspides’ donde las derivadas se anulan
en ϕ̇ ̸= 0 [2].
Llevando este concepto a la mecánica cuántica, Wilczek propuso un modelo basado en
una partı́cula de masa y carga unitarias confinada en un anillo unidimensional, el cual es
atravesado por un flujo magnético fraccionario de Aharonov-Bohm α. El Lagrangiano y
el Hamiltoniano del sistema están dados por:
1
L = ϕ̇2 + αϕ̇ (1)
2
1
H = (πϕ − α)2 (2)
2
∂
Donde el momento canónico es πϕ = −i ∂ϕ . Los autoestados del sistema son |l⟩ = eilϕ con
valores y números cuánticos enteros l. La energı́a del sistema es El = 12 (l − α)2 .
El estado fundamental ocurre para el entero l0 que minimice (l − α). Wilczek demostró
que, si el flujo α no es un número entero, el valor esperado de la velocidad angular en el
estado fundamental es distinto de cero [1]:
⟨l0 |ϕ̇|l0 ⟩ = l0 − α ̸= 0 (3)
Este resultado matemático fundamental sugiere la persistencia de una corriente en el
estado de mı́nima energı́a, lo que constituirı́a la base de un ’reloj cuántico’ espontáneo,
rompiendo ası́ la simetrı́a de traslación temporal.

1
Escuela Politécnica Nacional Mecánica Cuántica

2. Problema: El Teorema de Imposibilidad (No-Go

Theorem)
La propuesta original de Wilczek fue sometida a un riguroso escrutinio. En 2015,
Watanabe y Oshikawa demostraron que la definición ingenua de un cristal de tiempo
basada en el valor esperado de un observable ⟨Ô(t)⟩ es trivialmente constante en cualquier
ensamble de equilibrio (Gibbs) o en un autoestado Ĥ|n⟩ = En |n⟩, puesto que los factores
de evolución temporal e±iEn t se cancelan [3]:

Ô(t) = eiĤt Ô(0)e−iĤt (4)

Para ofrecer una definición rigurosa, formularon el orden cristalino temporal en fun-
ción del comportamiento asintótico de las funciones de correlación a larga distancia del
parámetro de orden Φ̂ [3]:
1 iĤt −iĤt
lı́m ⟨e Φ̂e Φ̂⟩ = f (t) (5)
V →∞ V 2

Si f (t) exhibe una oscilación periódica no trivial, el sistema es un cristal de tiempo.

Sin embargo, Watanabe y Oshikawa demostraron rigurosamente mediante la desigualdad
de Schwarz que, para un Hamiltoniano local Ĥ (con interacciones de corto alcance), la
variación de esta correlación en el estado fundamental está acotada por [3]:
1 −i(Ĥ−E0 )t t
⟨0|Âe B̂|0⟩ − ⟨0|ÂB̂|0⟩ ≤ C (6)
V2 V
Dado que en el lı́mite termodinámico (V → ∞) el lado derecho de la inecuación tiende
a cero para cualquier tiempo finito t, se concluye que f (t) debe ser una constante. Mediante
el uso de las cotas de Lieb-Robinson, extendieron esta prueba para temperaturas finitas
(T > 0). Este teorema (no-go theorem) clausuró matemáticamente la posibilidad de que
existan cristales de tiempo cuánticos en sistemas que se encuentren en equilibrio térmico.

3. Solución: Sistemas Floquet y Orden de Autoesta-

do
El teorema de Watanabe y Oshikawa dejó una única ventana matemática abierta: sus
restricciones no aplican a sistemas fuera del equilibrio. Esto motivó la transición hacia los
Sistemas Floquet, los cuales están descritos por Hamiltonianos dependientes del tiempo
con periodicidad H(t + T ) = H(t) [4].
La evolución temporal de estos sistemas está dictada por el operador de Floquet sobre
un periodo completo:
 Z T 
U (T ) = T exp −i H(t)dt ≡ e−iHF T (7)
0

En este régimen, la conservación de la energı́a se reemplaza por cuasi-energı́as ϵα defini-

das módulo 2π/T , y los estados estacionarios se sustituyen por autoestados de Floquet:
|ψα (t)⟩ = e−iϵα t |ϕα (t)⟩.
El principal obstáculo de los sistemas fuertemente impulsados es que obedecen la
Hipótesis de Termalización de Autoestados (ETH), absorbiendo energı́a hasta alcanzar

2
Escuela Politécnica Nacional Mecánica Cuántica

una ”temperatura infinita”(entropı́a máxima) donde cualquier correlación se destruye.

Como demostraron Khemani et al. (2016), la solución radica en introducir la Localización
de Muchos Cuerpos (MBL) mediante desorden en el sistema [4].
La fase MBL viola la ETH, impidiendo el transporte de energı́a y permitiendo que todo
el espectro de autoestados del espacio de Hilbert exhiba orden correlacionado (eigenstate
order ). Bajo estas condiciones, Khemani et al. identificaron una fase cuántica exótica
denominada Floquet π-Spin Glass (π-SG), donde los autoestados de Floquet se separan
por una cuasi-energı́a exacta de π/T . Esta separación espectral provoca que los observables
del sistema (como los espines matriciales) exhiban una respuesta sub-armónica, oscilando
con un periodo 2T . Ası́, el sistema rompe la simetrı́a temporal discreta del impulso externo
original, materializando el Cristal de Tiempo Discreto (DTC).

4. El Modelo Cuántico: Operador de Floquet

Para entender con rigor de mecánica cuántica por qué el perı́odo se duplica de T a
2T , desglosaremos cada paso algebraico utilizando únicamente álgebra de operadores de
Pauli e identidades exponenciales

4.1. 1. El Pulso que voltea los espines (Ûx )

π P
Ûx = e−i 2 j X̂j
(8)
Este operador aplica un campo magnético transversal muy fuerte por un instante. ¿Qué
hace fı́sicamente? Toma cada espı́n de la cadena y lo voltea perfectamente (si estaba en
|↑⟩, pasa a |↓⟩ y viceversa).

4.2. 2. La Interacción entre vecinos (ÛZZ )

P
ÛZZ = e−i j Jj Ẑj Ẑj+1
(9)
Este término representa la interacción clásica de Ising. Fı́sicamente, dice que los espi-
nes vecinos ”sienten”la presencia del otro. La constante Jj es aleatoria para cada par,
introduciendo el desorden necesario para evitar la termalización.

4.3. 3. El Campo magnético aleatorio (Ûz )

P
Ûz = e−i j hj Ẑj
(10)
Aquı́, cada espı́n siente un campo magnético local hj diferente y aleatorio. Esto añade
aún más desorden, asegurando que las partı́culas queden totalmente localizadas (MBL).

3
Escuela Politécnica Nacional Mecánica Cuántica

4.4. Paso 1: ¿Por qué el pulso Ûx actúa como un operador de

paridad?
Recordemos la definición del operador del pulso transversal:
L
PL
−i π2 π
Y
Ûx = e j=1 X̂j
= e−i 2 X̂j (11)
j=1

Para evaluar la exponencial de una matriz de Pauli, aplicamos la expansión en serie de

Taylor. Dado que X̂j2 = I:
∞
X (−iθX̂j )k
e−iθX̂j = = cos(θ)I − i sin(θ)X̂j (12)
k=0
k!

Sustituyendo el ángulo del pulso ideal θ = π2 :

π  π 
−i π2 X̂j
e = cos I − i sin X̂j = 0 − iX̂j = −iX̂j (13)
2 2
Al multiplicar sobre los L qubits de la cadena:
L
Y L
Y
L
Ûx = (−iX̂j ) = (−i) X̂j = (−i)L P̂X (14)
j=1 j=1

Donde P̂X = Lj=1 X̂j es el operador de paridad global de espı́n. Dado que (−i)L es
Q

simplemente una fase global escalar (eiϕ ), no afecta ningún valor de expectativa observable
ni las conjugaciones de operadores. Por lo tanto, operatoriamente Ûx ≡ P̂X .

4.5. Paso 2: Transformación de la matriz Pauli Ẑj bajo la pari-

dad P̂X
Evaluemos el producto de conjugación P̂X Ẑj P̂X† :
L
! L
!
Y Y
P̂X Ẑj P̂X† = X̂k Ẑj X̂m (15)
k=1 m=1

Dado que los operadores de Pauli en sitios distintos conmutan ([X̂k , Ẑj ] = 0 si k ̸= j),
todos los términos con k ̸= j se mueven libremente y se cancelan usando X̂k X̂k = I. Nos
queda únicamente el sitio j:
P̂X Ẑj P̂X† = X̂j Ẑj X̂j (16)
Utilizando la relación de anticonmutación fundamental {X̂j , Ẑj } = 0 =⇒ X̂j Ẑj =
−Ẑj X̂j :
X̂j Ẑj X̂j = (−Ẑj X̂j )X̂j = −Ẑj X̂j2 = −Ẑj (17)

4.6. Paso 3: Transformación de los operadores de evolución Ûz

y ÛZZ
†
Utilizaremos la identidad formal de conjugación de operadores AeB̂ A† = eAB̂A .

4
Escuela Politécnica Nacional Mecánica Cuántica

A. Para el campo magnético Ûz :

L
! L
!
X X
P̂X Ûz P̂X† = P̂X exp −i hj Ẑj P̂X† = exp −i hj (P̂X Ẑj P̂X† ) (18)
j=1 j=1

Sustituyendo P̂X Ẑj P̂X† = −Ẑj :

L
! L
!
X X
P̂X Ûz P̂X† = exp −i hj (−Ẑj ) = exp +i hj Ẑj = Ûz† (19)
j=1 j=1

B. Para la interacción de Ising ÛZZ : Primeramente, hallemos cómo transforma el

producto de dos espines adyacentes:
P̂X (Ẑj Ẑj+1 )P̂X† = (P̂X Ẑj P̂X† )(P̂X Ẑj+1 P̂X† ) = (−Ẑj )(−Ẑj+1 ) = +Ẑj Ẑj+1 (20)
Al no cambiar de signo el producto doble, la exponencial no se altera:
P̂X ÛZZ P̂X† = ÛZZ (21)

4.7. Paso 4: Cálculo explı́cito de la evolución tras dos perı́odos

(ÛF2 )
El operador de Floquet de un perı́odo es ÛF = Ûz ÛZZ P̂X . Al elevarlo al cuadrado:
ÛF2 = (Ûz ÛZZ P̂X )(Ûz ÛZZ P̂X ) (22)
Insertamos el operador identidad I = P̂X† P̂X justo después del primer grupo de operadores
para reordenar:  
†
ÛF = Ûz ÛZZ P̂X Ûz ÛZZ P̂X P̂X2
2
(23)

Dado que P̂X2 = j X̂j2 = I, y aplicando las conjugaciones del Paso 3:

Q

ÛF2 = Ûz ÛZZ (Ûz† ÛZZ ) (24)

Como tanto Ûz como ÛZZ son operadores construidos únicamente con matrices Ẑj , todos
son diagonales en la base computacional y conmutan perfectamente entre sı́ ([Ûz , ÛZZ ] =
0). Agrupando los términos:
ÛF2 = (Ûz Ûz† )(ÛZZ ÛZZ ) = I · ÛZZ
2 2
= ÛZZ (25)
1. En un perı́odo (t = T ), el operador contiene a P̂X , por lo que cualquier estado con
espines definidos en Z (como |↑↑ . . .⟩) se invierte completamente a |↓↓ . . .⟩.
2. En dos perı́odos (t = 2T ), la paridad se aplica dos veces (P̂X2 = I), y los campos
aleatorios localizadores hj se cancelan algebraicamente (Ûz Ûz† = I).
3. El sistema vuelve a la configuración inicial orientada hacia arriba, demostrando que
la frecuencia fı́sica observable es Ω/2, o equivalentemente, el perı́odo exacto es 2T .

5. Estado actual: Hitos Experimentales(MBL-DTC)

Demostración experimental de la fase MBL-DTC: Resultados del procesador Sycamore
(Google Quantum AI)

5
Escuela Politécnica Nacional Mecánica Cuántica

El Protocolo Experimental: Diseño del Circuito

El experimento se llevó a cabo utilizando una cadena lineal de 20 qubits superconduc-
tores tipo transmon. La dinámica del sistema se rige por la aplicación periódica (Floquet)
de un operador de evolución temporal unitario, denotado como ÛF . Este circuito consta
de tres partes fundamentales por cada ciclo:

1. Interacciones Ising desordenadas: Se aplican compuertas CPHASE ajustables

para generar interacciones Ẑi Ẑi+1 con fases aleatorias. Este desorden es el ingredien-
te crı́tico para inducir la Localización de Muchos Cuerpos (MBL) y evitar que el
sistema termalice.

2. Campos longitudinales aleatorios: Se aplican rotaciones Ẑ locales para romper

cualquier simetrı́a espacial remanente.

3. Pulso de rotación transversal (El “Kick”): Se aplica un pulso global que rota
todos los qubits un ángulo πg alrededor del eje X. Cuando g ≈ 1, el pulso invierte
casi perfectamente el estado de los espines en cada ciclo.

Resultados Experimentales

Figura 1: Resultados experimentales del MBL-DTC

Para confirmar la existencia de la fase MBL-DTC, el experimento debe demostrar

una respuesta subarmónica robusta (un perı́odo que es el doble del perı́odo de la fuerza
impulsora).

Análisis de la Polarización (Panel 2b)

El panel 2b de la Figura 2 muestra la medida de la polarización local del qubit, ⟨Ẑ(t)⟩,
en función del número de ciclos. Se contrastan dos regı́menes:

6
Escuela Politécnica Nacional Mecánica Cuántica

Fase Termal (g = 0,60): La polarización decae rápidamente hacia cero en menos de

10 ciclos. El sistema absorbe energı́a del pulso y termaliza rápidamente, perdiendo
la memoria de su estado inicial.

Fase MBL-DTC (g = 0,97): Se observa una oscilación persistente de ⟨Ẑ(t)⟩ que

alterna su signo en cada ciclo, manteniendo un perı́odo duplicado que sobrevive a
más de 100 ciclos.

Aislamiento de la Decoherencia (Paneles 2c y 2d)

En el panel 2b se observa que, aunque la oscilación persiste en g = 0,97, existe una en-
volvente de decaimiento. Para distinguir si este decaimiento proviene de una termalización
interna (fı́sica del sistema) o simplemente de decoherencia externa (ruido del hardware),
se implementó un circuito de inversión temporal o “eco”, ÛECHO .
Este circuito hace evolucionar el sistema hacia adelante y luego revierte la evolución.
Cualquier desviación de la identidad pura cuantifica el ruido externo del hardware. Al
normalizar el autocorrelador original (Ā) con la señal del eco (Ā0 ), el panel 2d muestra
que la señal corregida (Ā/Ā0 ) mantiene una amplitud casi máxima. Los resultados ob-
tenidos demuestran un Cristal de Tiempo Discreto porque cumplen con sus dos firmas
fundamentales simultáneamente:

El sistema rompe espontáneamente la simetrı́a de traslación temporal discreta del

pulso impulsor, oscilando con el doble de su perı́odo de forma indefinida.

La normalización con el circuito eco demuestra que la caı́da de la señal es puramen-

te ruido del hardware. Internamente, el desorden introducido por las compuertas
CPHASE crea una Localización de Muchos Cuerpos (MBL) que evade la termaliza-
ción, estabilizando el orden en todo el espectro de autoestados del sistema.

6. Conclusiones
El viaje de los cristales de tiempo, desde su prohibición teórica en el equilibrio hasta su
observación en procesadores cuánticos, representa un triunfo de la ingenierı́a de precisión.
Hemos aprendido que la MBL y el desorden no son imperfecciones, sino herramientas
para crear fases de la materia que solo existen dinámicamente. Este avance redefine nues-
tra capacidad tecnológica: estos sistemas actúan como memorias cuánticas robustas que
”protegen”su periodicidad frente al ruido. En última instancia, la realidad del cristal de
tiempo nos enseña que las leyes de lo posible son maleables cuando somos capaces de or-
questar la danza de los autoestados fuera de las fronteras del equilibrio térmico, abriendo
una nueva era en la fı́sica de la materia condensada no-equilibrio.

Referencias
[1] Wilczek F. Quantum time crystals. Phys Rev Lett. 2012;109(16):160401. doi:
10.1103/PhysRevLett.109.160401.

[2] Shapere A, Wilczek F. Classical time crystals. Phys Rev Lett. 2012;109(16):160402.
doi: 10.1103/PhysRevLett.109.160402.

7
Escuela Politécnica Nacional Mecánica Cuántica

[3] Watanabe H, Oshikawa M. Absence of quantum time crystals. Phys Rev Lett.
2015;114(25):251603. doi: 10.1103/PhysRevLett.114.251603.

[4] Khemani V, Lazarides A, Moessner R, Sondhi SL. Phase structure of driven quantum
systems. Phys Rev Lett. 2016;116(25):250401. doi: 10.1103/PhysRevLett.116.250401.

[5] Zhang J, Hess PW, Kyprianidis A, Becker P, Lee A, Smith J, et al. Observation of a
discrete time crystal. Nature. 2017;543(7644):217-20. doi: 10.1038/nature21413.

[6] Choi S, Choi J, Landig R, Kucsko G, Zhou H, Isoya J, et al. Observation of

discrete time-crystalline order in a disordered dipolar many-body system. Nature.
2017;543(7644):221-5. doi: 10.1038/nature21426.

[7] Mi X, Ippoliti M, Quintana C, Greene A, Chen Z, Gross J, et al. Time-

crystalline eigenstate order on a quantum processor. Nature. 2022;601(7894):531-6.
doi: 10.1038/s41586-021-04257-w.