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:
- 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.
- 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!