Dimensión de hamiltoniano y diagonalizabilidad

A menudo, en la literatura de física de la materia condensada, uno se encuentra con un hamiltoniano que dice algo así como:

H = i = 1 norte j i   S i z S i + 1 z ,
dónde j i son las constantes de acoplamiento, S i z representando al operador de espín S z actuando en i el sitio y norte siendo el número total de sitios (quizás con una identificación en forma de anillo). Ahora, S z es un 2 × 2 matriz con el espín en cualquier sitio en particular siendo un 2 × 1 matriz o vector columna.

Si uno supiera la función de onda como | ψ = | i | i + 1 , entonces sería fácil simplemente aplicar el hamiltoniano y hacer lo que uno quiera. Pero supongamos que no tengo eso | ψ y necesita diagonalizar directamente H para calcular sus valores propios y vectores propios (por ejemplo, al simular un modelo de este tipo en una red en la computadora), ¿cómo lo haría?

¿Cómo se aseguraría uno de que S i z es diferente de S i + 1 z ? La única forma de hacerlo es operarlos en algún giro supuesto (aleatorio) existente en sus sitios. Pero si uno hace eso, entonces el hamiltoniano es una suma del producto de dos 2 × 1 matrices ( [ 2 × 2 ]     [ 2 × 1 ] [ 2 × 1 ] matriz) que no tiene sentido!

¿Alguien podría explicar dónde me estoy equivocando? Básicamente, quiero diagonalizar un hamiltoniano de la forma anterior en una computadora (usando, digamos, Mathematica, que tiene funciones incorporadas para dar los valores propios y los vectores propios) sin conocer las funciones propias.

Editar: después de pensarlo más, al decir que el giro en cualquier sitio en particular es un 2 × 1 matriz, ¿estoy asumiendo implícitamente que es un estado no tensorial? Incluso si eso es cierto, ¿cómo ayuda a dar la dimensión correcta del hamiltoniano para diagonalizarlo?

Respuestas (1)

Dimensión del espacio de Hilbert

El espacio de Hilbert de una cadena de espín de longitud norte viene dada por el producto tensorial de norte copias del espacio de Hilbert de un solo giro. Suponiendo que estamos hablando de spin- 1 2 , la dimensión del espacio es 2 norte .

En este contexto, un "operador que actúa en un solo giro" significa un 2 norte -operador dimensional obtenido como el producto tensorial de un operador 2d con norte 1 copias de la identidad 2d, por lo que la acción en los otros giros es trivial.

Por ejemplo, dada la ( 2 × 2 ) Matriz de Pauli σ z ,

S i z = I 2 I 2 i 1 σ z I 2 I 2 norte i .

De manera similar, un operador de "dos cuerpos" S i z S i + 1 z se puede escribir explícitamente

S i z S i + 1 z = I 2 I 2 i 1 σ z σ z I 2 I 2 norte i 1 ,
etcétera.

Diagonalización numérica

Diagonalizar tal hamiltoniano numéricamente es un problema difícil. En principio, con las expresiones anteriores puede construir explícitamente la matriz en su lenguaje de programación favorito (por ejemplo, en Mathematica usando KroneckerProducty PauliMatrix). Sin embargo, como se señaló anteriormente, la dimensión crece exponencialmente con el número de giros, lo que significa que ni siquiera podrá llegar a norte 20 .

Por lo tanto, uno necesita encontrar formas más inteligentes de hacer las cosas; por ejemplo, un punto de partida sería representar H como una matriz dispersa. Este es un campo activo de investigación, e incluso con algoritmos de última generación y grandes computadoras, las personas pueden resolver el problema con unas pocas decenas de giros como máximo.