15 votos

Sutilezas de la solución exacta a la 1D cuántica XY modelo, en particular la transformación de Bogoliubov

Estoy escribiendo programas para la construcción de los espectros de los modelos conocidos con soluciones exactas, y pronto se percataron de algunas sutilezas que a menudo no se menciona en la mayoría de las referencias. Estas sutilezas no son importantes en el límite termodinámico, pero ya que estoy escribiendo recetas numérica de tamaño finito de los sistemas, que hacen de la materia. Ahora todavía hay un rompecabezas sobre Bogoliubov transformación (o eso creo) que no puedo resolver por mi cuenta.

Tomemos, por ejemplo, la 1D anisotrópico XY modelo en un ámbito transversal (incluyendo el caso del modelo de Ising), dado por la siguiente Hamiltoniano $$ H = - J \sum_{j=1}^N \left( g \sigma_j^z + \frac{1+\eta}{2} \sigma_j^x \sigma_{j+1}^x + \frac{1-\eta}{2} \sigma_j^y \sigma_{j+1}^y\right), $$ suponiendo periódica de la condición de límite $\vec{\sigma}_{N+1} = \vec{\sigma}_1$. Este modelo puede ser resuelto mediante la aplicación de Jordania-Wigner transformación, la transformada de Fourier, y Bogoliubov transformación sucesivamente. El final de la diagonal de la forma del Hamiltoniano, como se da en muchas referencias, es $$ H = \sum_k \varepsilon(k) (b^\dagger_k b_k - \frac{1}{2} ) $$ donde $ \varepsilon(k) = 2 J \sqrt{(g - \cos(ka))^2 + \eta^2 \sin^2(ka)}$.

Dos matices:

  1. El original espacio de Hilbert de periódico de tirada de la cadena se asigna a DOS sectores de Fock espacio: par/impar paridad sector ($n = \sum c^\dagger_j c_j$/impares) con periódico/anti-periódico de la condición de límite $c_1 = \pm c_{N+1}$ respectivamente.
  2. El impulso $k$ toma valores de entero/media-múltiplos enteros de $2\pi/Na$ en el par/impar paridad del sector, respectivamente.

Otro de los subtítulos que (casi) sigue inmediatamente (pero sin resolver) es que la relación de dispersión dada anteriormente en realidad sólo le da el valor absoluto de un solo fermión de la energía. La verdad solo fermión de energía debe ser dada por $$ \varepsilon(k) = \pm 2 J \sqrt{(g - \cos(ka))^2 + \eta^2 \sin^2(ka)}. $$ Este indeterminado signo que importa en un tamaño finito del sistema debido a que hay dos paridad sectores con diferentes opciones de $k$, de ahí la "presencia" y la "ausencia" de fermiones con "energías opuestas" no son equivalentes (en realidad, no puede tener exactamente energías opuestas en general). Incluso si este signo de la ambigüedad no afecta a la energía total $E_n$ sí, ciertamente aún perturba la relación de dispersión del espectro total $E_{n,k}$ (este es el caso cuando se $g=\eta=0$).

Yo creo que este signo se determina en el Bogoliubov transformación, pero no sé cómo. El Hamiltoniano de la derecha antes de la Bogoliubov transformación es $$ H = J \sum_k \left\{ 2(g- \cos(ka)) c^\dagger_k c_k + i \eta \sin(ka) (c^\dagger_{-k} c^\dagger_k + c_{-k} c_k ) \right\} + \text{const.} $$ Por ejemplo, en el caso isotrópico $\eta = 0$, mediante la comparación de los casos de prueba con los resultados exactos de diagonalización, puedo saber la correcta único fermión de energía está dado por $ \varepsilon(k) = 2J (g - \cos(ka)) $, no cualquier otras variantes. Sin embargo, al hacer la Bogoliubov transformación: $$ c_k = \cos \frac{\theta_k}{2} b_k + i \sin \frac{\theta_k}{2} b^\dagger_{-k}, $$ y suponiendo que $\cos \frac{\theta_k}{2} = \cos \frac{\theta_{-k}}{2}, \; \sin \frac{\theta_k}{2} = -\sin \frac{\theta_{-k}}{2} $, obtenemos $$ H_{\eta=0} = J \sum_k \left\{ (g-\cos(ka) ) \left[ 2 \cos \theta_k b^\dagger_k b_k + i \sin \theta_k ( b^\dagger_k b^\dagger_{-k} - b_{-k} b_k ) \right] \right\} + \text{const.} $$ La única restricción de que puedo ver es $\sin\theta_k = 0$. Si ese es el caso, entonces la $\cos\theta_k$ puede tomar $\pm 1$ sin embargo ella le gusta. Esto NO es lo que sucede! Sólo $\cos\theta_k = 1$ da la solución correcta!

Si alguien podría señalar lo que soy o puede ser que falte aquí, o me apunte a alguna referencia de la misma (la mayoría de las referencias no hable de tamaño finito de casos en esos detalles), eso sería genial. Gracias!

10voto

Andy Hume Puntos 15186

Después de días de reflexión, de búsqueda, los debates, y las pruebas, por fin puedo responder a mi propia pregunta ahora. La respuesta es mucho más complicado de lo que yo esperaba de un "simple" XY modelo (sólo para el modelo de Ising)!

Todas las "soluciones correctas/espectro", declaró a continuación se comparan con los resultados de exacta de diagonalización.


Simplemente, no hay nada malo con la elección de señas para el fermión de energía $\varepsilon(k)$, y no importa (ni) si nos atenemos a una sola rama de $\theta_k = \arctan \xi$. El aspecto crucial a tener en cuenta (y hacer las cosas bien) es que, después de hacer el Bogoliubov transformación, el $b$-fermión vacío $|0\rangle_b$ es diferente de la $c$-fermión vacío $|0\rangle_c$, y puede tener diferentes paridad $P_0$ e ímpetu $k_0$!

Para ilustrar mi afirmación anterior, en primer lugar quisiera aclarar y arreglar mi terminología y notaciones:

  • $c_j$ son el fermión operadores después de la Jordania-Wigner transformación, $c_k$ son el fermión operadores después de la transformada de Fourier, $b_k$ son el fermión operadores después de la Bogoliubov transformación.
  • $a$-fermión vacío de estado se refiere a la cero partículas estado $|0\rangle_a$ aniquilado por fermión operadores de $a$'s. Este NO es necesariamente el estado fundamental del sistema.
  • La paridad siempre se refiere a el operador $P \equiv \prod_j (-\sigma_j^z) = (-1)^{\sum_j c_j^\dagger c_j}$ o de sus correspondientes autovalores. Tenga en cuenta que se define en la vuelta de base, y tiene una simple representación de la $c$-fermión base, pero no necesariamente en el $b$-fermión.
  • $P_0$ $k_0$ se refieren a la paridad y el impulso de un vacío.

Como se indica en la pregunta, después de la Jordania-Wigner transformación, el original espacio de Hilbert de los periódicos de tirada de la cadena se asigna en dos sectores de Fock espacio con $P=\pm 1$ (par/impar), respectivamente. El $c$-fermión vacío estados $|0\rangle_c$ en ambos sectores han$P_0=1$$k_0=0$. La transformada de Fourier respeta la paridad y el número de partículas de conservación, y no cambia el estado de vacío. Si el Hamiltoniano es ya diagonal en este Fourier base (es decir, el caso isotrópico $\eta=0$), se puede construir el espectro completo (incluyendo relación de dispersión) el uso de $n$-partícula estado $c_{k_1}^\dagger c_{k_2}^\dagger \dots c_{k_n}^\dagger |0\rangle_c$ cuya paridad y el impulso dado por $P = (-1)^n$$k = k_1 + k_2 + \dots + k_n$. Con la elección adecuada de los $n$ (par/impar) y $k$ (enteros/media-enteros) en los dos sectores, no hay más matices.


Demostración con el $\eta=0, \theta_k = \pi$ caso

Las cosas se complican con la Bogoliubov transformación, debido a que $|0\rangle_b \neq |0\rangle_c$, y el último estado tiene diferentes $P_0$$k_0$, e incluso pueden no ser el mismo en los dos sectores! Para fines de demostración, vamos a tomar el $\eta=0$ isotrópica caso, y hacer un trivial Bogoliubov transformación con $\theta_k = \pi$: $$ c_k = i b_{-k}^\dagger, \quad c_k^\dagger = -i b_{-k}. $$ No hay nada malo con esta transformación. Sólo necesitamos ser cuidadosos acerca de la nueva vacío $|0\rangle_b$. En este caso, hay una relación simple entre los dos vacío estados: $$ |0\rangle_b = \prod_k c_k^\dagger |0\rangle_c. $$ Esto significa que, en el $b$-fermión base:

  • Si el número total de sitios de $N$ es incluso, $P_0=1$ en ambos sectores, $k_0 = 0$ o $\pi/2$ en la paridad par/impar sector, respectivamente.
  • Si el número total de sitios de $N$ es impar, $P_0=-1$ en ambos sectores, $k_0 = \pi/2$ o $0$ en la paridad par/impar sector, respectivamente.

La paridad y el impulso de un estado con $n$ $b$-fermiones está dado por $P = (-1)^n P_0$$k = k_0 + k_1 + \dots + k_n$. El uso de estas fórmulas y la elección del $P$ $k$'s en los dos sectores, todavía puedo construir el espectro correcto y la relación de dispersión con $\theta_k = \pi$.


La discusión sobre el caso general

El principio es el mismo que en el caso general ($\eta \neq 0$). Todo lo que tenemos que hacer es determinar $P_0$ $k_0$ $|0\rangle_b$ en los dos sectores, y el uso de $P = (-1)^n P_0$ $k = k_0 + k_1 + \dots + k_n$ a determinar la paridad y el impulso de un $n$ $b$-fermión estado. Sin embargo, la determinación del $P_0$ $k_0$ se convierte en algo muy trivial, porque no hay una relación simple entre el $|0\rangle_b$ $|0\rangle_c$ al $\eta \neq 0$ (sin embargo, ver referencia al final).

Como calcular correctamente en @cesaruliana la respuesta, tenemos $$ \tan\theta_k = \frac{ \eta \sin(ka)}{ g - \cos(ka)}, $$ todas las ramas de $\theta_k = \arctan\xi$ dar validez transformaciones, pero no es obvio que una(s) es la más útil. Nuestro objetivo es encontrar a $P_0$ $k_0$ de la tierra del estado. Una de las estrategias que podemos usar es la elección de $$ \cos \theta_k = \frac{ g - \cos(ka) }{ \sqrt{ (g-\cos(ka))^2 + \eta^2 \sin^2(ka)} }, \quad \sin \theta_k = \frac{ \eta \sin(ka) }{ \sqrt{ (g-\cos(ka))^2 + \eta^2 \sin^2(ka)} }, $$ por lo tanto, haciendo todo el single $b$-fermión energías positivas: $$ \varepsilon(k) = +2J \sqrt{ (g-\cos(ka))^2 + \eta^2 \sin^2(ka)}. $$ Ahora el vacío de estado $|0\rangle_b$ es el estado del suelo de las $b$-fermión sistema, por lo que el estado fundamental (de la fermión o el sistema de espín) mucho más fáciles de rastrear. A continuación, analizamos el estado de la vuelta de tuerca original modelo en diferentes límites.

Por ejemplo, en el Ising campo débil límite $ H = -J \sum \sigma_j^z \sigma_{j+1}^z, $ hay dos degenerados de tierra de los estados con relación a las paridades:$ |\pm\rangle = \prod_j \left| \leftarrow \right>_j \pm \prod_j \left| \rightarrow \right>_j$, $ P |\pm \rangle = \pm | \pm \rangle$ y ambos tienen el ímpetu $k_0=0$. Se nota que esta degeneración es la característica de una ruptura de la simetría de fase. Así que, hasta que una transición de fase se produce (en $g=1$ o $\eta=1$), el dos veces la degeneración debe estar siempre presente y exacta, y el número cuántico $P$ $k$ de los dos estados fundamentales debe ser invariable. Tener tal ("robusta") dos veces la degeneración de todos los valores de $0\leq g\leq 1$$0\leq \eta \leq 1$, los dos de tierra de los estados deben ser los dos al vacío estados $|0\rangle_b$ de los dos sectores. Esto significa que, para $0 \leq g\leq 1$, $0 \leq \eta \leq 1$, $P_0 = + 1$ en la paridad incluso del sector y $P_0 = -1$ en la paridad impar sector, y $k_0=0$ para ambos sectores. Por lo tanto, en ambos sectores, sólo los estados con un número de $b$-fermiones deben ser incluidos en la final del espectro. El espectro correcto es, de hecho, construido a partir de esta construcción.

$P_0$ $k_0$ puede cambiar después de una fase de transición. Similar análisis debe ser hecho de nuevo en otras fases.

El error que hacía era que me calcula la paridad y el impulso ingenuamente el uso de $b$-fermiones, siempre tomando en cuenta los estados con pares/impares número de $b$-fermiones en la paridad par/impar sector, esencialmente, siempre asumiendo $P_0=1, k_0=0$. Esto NO es cierto en general.


Agradecimientos y referencias

Me gustaría agradecer a mis compañeros de la estudiante Wen Wei muy útil para la discusión. Él descubrió que la solución correcta puede ser construido si tomamos el número de partículas en ambos sectores en la ruptura de la simetría de fase. Una muy útil la nota escrita por el Prof. Andreas Schadschneider y el Prof. Götz S. Uhrig que finalmente se resuelva en la confusión que se puede encontrar aquí:
Fuertemente Correlacionada de los Sistemas en la Física del Estado Sólido

Actualización: Después de que envió esta respuesta, me encontré con otro de referencia (tesis de master) por Erik Eriksson donde él rigurosamente calculado en detalle la relación entre el $|0\rangle_b$ ("Bogoliubov vacío") y $|0\rangle_c$ en el caso del modelo de Ising $\eta=1$. Creo que esta puede ser fácilmente generalizado para el caso de $\eta < 1$. Su tesis se puede encontrar aquí:
Quantum Transiciones de Fase en un Integrable unidimensional Spin Modelo

2voto

Julien N Puntos 1544

En el caso de $\eta=0$ el original de Hamilton ya es diagonal, por lo que al hacer la Bogoliubov transformación y, a continuación, tomar $\eta=0$ es poco definidas, y eso es lo que parece confuso.

A ver lo que está pasando supongamos $\eta\neq 0$. Al realizar la transformación e imponer que es diagonal usted debe conseguir (si no me han hecho ningún error)

$\tan \theta_k=\frac{\eta\sin{(ka)}}{g-\cos{(ka)}}$.

Por lo tanto, el ángulo se define por $\theta_k=\arctan \xi$ donde $\xi$ es la proporción de arriba. Aquí está el problema. Dado que las funciones trigonométricas no son bijective sus inversas de las funciones no están definidas en todo el eje real. Por lo que necesita para elegir una rama donde la función inversa es única valorada, la opción habitual es, por supuesto,$(-\pi/2,\pi/2)$, la denominada rama principal. Ahora, si definimos $\theta_k$ en este intervalo y tome $\eta\rightarrow 0$ sólo hay una solución, a saber,$\theta_k=0$, lo que implica $\cos \theta_k=1$.

Esto resuelve la señal de ambigüedades que usted menciona, ya que la energía también es determinada por otra función trigonométrica y en este intervalo de ellos tienen valores únicos, incluyendo signo.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X