40 votos

Condiciones límite periódicas vs. abiertas

En la materia condensada, se suelen utilizar condiciones de contorno periódicas para realizar cálculos sobre las propiedades del grueso de un material. En general, se argumenta que en la $N\rightarrow\infty$ limitar las condiciones de contorno no afectan a las propiedades del grueso, por lo que se pueden utilizar condiciones de contorno periódicas para calcular las propiedades del grueso de sistemas con fronteras abiertas.

¿Existe alguna prueba matemática formal de este hecho? Estoy pensando en afirmaciones como

  • Como $N\rightarrow\infty$ cualquier observable que sólo mire al grueso no se ve afectado por las condiciones de contorno.
  • Como $N\rightarrow\infty$ el solapamiento del estado básico del sistema periódico y el estado básico del sistema abierto llega a 1.
  • Como $N\rightarrow\infty$ La matriz de densidad reducida en la masa no depende de las condiciones de contorno.

O algo similar. No busco argumentos intuitivos, sino pruebas en la literatura si es que existen.


Un pensamiento que tuve después de pensar en esto por un tiempo: Creo que debe ser cierto que el hamiltoniano con condiciones de contorno periódicas y el hamiltoniano con condiciones de contorno abiertas deben estar conectados adiabáticamente.

Un ejemplo sencillo es el efecto de los campos magnéticos límite en el modelo de Ising. Si modifico el hamiltoniano del modelo de Ising en la frontera, puedo cambiar la física en la masa. Consideremos

$$ H=-\sum_{n=1}^{N-1} J \sigma^z_n\sigma^z_{n+1} +h\sigma_1 $$ Si $h$ es positivo, el estado básico es todo espín abajo; si $h$ es negativo, es todo giro. El cambio del Hamiltoniano en la frontera cambió la física en la masa. ¿Por qué ocurrió esto? Es porque hay un cruce de niveles cuando se sintoniza $h$ de positivo a negativo; un estado excitado se convierte en el estado de tierra, y viceversa. En general, si el cambio de las condiciones de contorno induce un cruce de niveles, deberíamos esperar que la física de la masa cambie.

Por lo tanto, creo que cualquier teorema que diga que la física de la masa entre los dos sistemas es idéntica sólo debe ser cierto si los sistemas abierto y periódico están conectados adiabáticamente.

13voto

BigManning Puntos 1

La respuesta a esta pregunta varía en dificultad de un sistema a otro. Sin embargo, toda esta historia es rigurosamente conocida y probada en el ejemplo más simple no trivial: el modelo Ising clásico 2D (el argumento también funciona en el caso más simple del modelo Ising clásico 1D, pero entonces el fenómeno no describirá el efecto de frontera de ningún modelo cuántico de interés):

\begin{align*} Z&=\sum_{\{s_{ij}\}}e^{-\beta S}\\ &~~~~S= \sum_{i,j\,\in \,L}\,(J_x\,s_{ij}\,s_{ij+1}+J_y\,s_{ij}\,s_{i+1\,j})\\ \end{align*}

El entramado $L$ es semi-infinito (es decir, tiene un límite a la izquierda) y se visualiza a continuación:

enter image description here

Podemos entonces imaginarnos tomando un observable local $\mathcal O$ (cuyo soporte se visualiza a continuación:) y formando la traslación $T_j(\mathcal O)$ del observable $j$ unidades en la masa:

enter image description here

Si imaginamos tomar el valor de la expectativa de tal observable en el límite que $j$ tiende a infinito, recuperamos el valor de la expectativa de ese mismo observable, pero evaluado en el plano completo. Esta es la valor límite termodinámico del observable, y se conoce la tasa de aproximación: \begin{align*} \langle T_j(\mathcal O)\rangle-\langle T_\infty(\mathcal O)\rangle\sim\,\,\, O(e^{-j/\zeta})\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\tag{1}\\ \end{align*} donde $\zeta$ es la longitud de correlación en el modelo de masa. Esta es la primera afirmación que querías demostrar:

(1) Como $N\to \infty$ cualquier observable que sólo mire al grueso no se ve afectado por las condiciones de contorno.

$$$$

$$$$

Esquema de la demostración matemática de (1)

Daré un esbozo de la prueba de (1). Comenzamos reescribiendo la función de partición en términos de la matriz de transferencia, que actúa sobre el espacio de Hilbert de las configuraciones de una sola columna en la red: \begin{align*} Z&=\lim_{N\to\infty}\sum_{\{s_0\},\, \{s_N\}}\, \langle \{s_0\}|T^N|\{s_N\}\rangle,\,\,\,\,\,\,\,\,\\ &\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,T:=e^{-\beta \sum_jJ_x\sigma^x_j\sigma^x_{j+1}}\cdot e^{-\beta \sum_jJ_y\sigma^z_j}\\ \end{align*} Aquí, $\{s_0\},\{s_N\}$ denotan las configuraciones de los espines en la columna del límite y en la columna del bulto, respectivamente. Si tenemos un observable local en el plano semi-infinito, entonces su valor puede expresarse en términos de la matriz de transferencia $T$ de la siguiente manera: enter image description here En las ecuaciones: \begin{align*} \langle T_j(\mathcal O)\rangle&=\lim_{N\to\infty}\frac{\sum_{\{s_0,s_N\}}\langle\{s_0\}|T^j\,\widehat O T^{N-j} |s_N\rangle}{\sum_{\{s_0\}}\langle\{s_0\}|T^N|\{s_N\}\rangle}\\ \end{align*} Inmediatamente, visualmente podemos ver por qué el límite es irrelevante: en la fase de alta temperatura del modelo, $T$ tiene un hueco y un único vector propio máximo. Realizando el límite termodinámico medio-infinito se sustituye $T$ con proyección sobre su vector propio máximo, que denotamos por $|GS\rangle$ : \begin{align*} T^j=|GS\rangle\langle GS|+O(e^{-j/\zeta}) \end{align*} donde $1/\zeta$ sea la diferencia entre el mayor y el segundo mayor valor propio de $T$ . Por lo tanto, sustituyendo esto en el observable \begin{align*} \langle T_j(\mathcal O)\rangle&=\frac{\sum_{\{s_0,s_N\}}\langle\{s_0\}|GS\rangle \langle GS|\,\widehat O|GS\rangle\langle GS |s_N\rangle}{\sum_{\{s_0\}}\langle\{s_0\}|GS\rangle \langle GS|\{s_N\}\rangle} + O(e^{-j/\zeta})=\langle GS|\widehat{\mathcal O}|GS\rangle+O(e^{-j/\zeta}) \end{align*} $$$$ Let's compare this with the result we would have obtained in the bulk: \begin{align*} \langle \mathcal O\rangle_\text{bulk}&:=\lim_{j\to \infty}\langle T_j( \mathcal O)\rangle = \langle GS|\widehat{\mathcal O}|GS\rangle. \end{align*} By subtracting this from the finite-$ j $ result, this finishes the proof of (1) in the simple case $ T>T_c$, mostrando que el rango característico de los efectos de frontera en los observables locales es igual a la longitud de correlación del bulto.

$$$$

Del mismo modo, utilizando este argumento de convergencia de la matriz de transferencia, se puede establecer:

(2) Como $N\to \infty$ el solapamiento del estado básico del sistema periódico y el estado básico del sistema abierto llega a 1.

(3) Como $N\to \infty$ La matriz de densidad reducida en la masa no depende de las condiciones de contorno.

De nuevo, como en la prueba de (1), sustituimos $T^j \sim |GS\rangle \langle GS|$ para $j$ suficientemente grande. De nuevo, el ingrediente clave aquí es la convergencia de las potencias de la matriz de transferencia $T^j$ que pierde cualquier "memoria" de efectos/ordenación de tamaño finito como $j\to \infty$ .

$$$$

$$$$

Extensión a un sistema cuántico a temperatura cero

Después de haber esbozado una prueba rápida de la intuición "los efectos de frontera no importan" en un entorno clásico, en realidad es trivial generalizar esto para cuantificar el efecto de frontera en un estado básico cuántico. En realidad, no necesitamos rehacer la prueba; en su lugar, la trasladamos al entorno cuántico utilizando la correspondencia cuántica-clásica que establece, a grandes rasgos, que

$$ \{\text{$ d+1 $-dim'l Classical system at $ T>0 $}\}\leftrightarrow \{\text{$ d $-dim'l Quantum system at $ T=0 $}\}$$

enter image description here

Por lo tanto, para demostrar (1-3) para un estado básico cuántico, basta con utilizar el mapeo cuántico-clásico. Por ejemplo, para el caso del modelo de Ising de campo transversal 1+1D (TFIM), aplicamos la transformación de Suzuki-Trotter con paso de tiempo $\Delta \tau>0$ a la función de partición cuántica $$Z=\lim_{\beta\to \infty}\text{tr}(e^{-\beta H_{TFIM}}),$$ que produce una familia de acciones eficaces $\{S[\Delta\tau]\}_{\Delta\tau}$ que describe las correlaciones estadísticas de un campo de Ising discreto en una red cilíndrica de espaciotiempo. La familia de acciones efectivas es: $$$$ \begin{align*} S[\Delta\tau]\underset{\,\,\Delta\tau\to \,0\,\,}{\sim} \sum_{j\tau}(J[\Delta\tau]\,s_{j\tau}s_{j+1\tau}+J_\perp[\Delta\tau]\, s_{j\tau}s_{j\tau+\delta\tau})\,\,\,\,\,\,\,\,\,\,\,\,\,.\\ \end{align*} Por supuesto, el efecto límite en el modelo cuántico es entonces precisamente el efecto límite en el modelo clásico 2D, que acotamos en los argumentos anteriores (usando la matriz de transferencia).

$$$$

$$$$

Más allá de este ejemplo ilustrativo

Espero que a estas alturas esto esté claro: utilizando la poderosa combinación de la correspondencia cuántico-clásica con el método de la matriz de transferencia (todos ellos métodos estándar en la caja de herramientas de la materia condensada), se pueden cuantificar los efectos de frontera en una amplia clase de estados básicos cuánticos, mucho más allá del ejemplo aquí con el TFIM 1+1D. En cualquier caso, espero que esta respuesta explique por qué los físicos tienen la intuición que tienen, y también demuestre que una prueba matemática del efecto de frontera, al menos en un entorno bastante general de sistemas de espín con hueco "trotterizable" con interacciones de vecino más cercano, está bien establecida (o al menos debe ser) en la literatura.

2voto

Vadim Puntos 377

Matemáticamente, una función periódica que satisface ciertas condiciones generales puede expandirse en Serie de Fourier . Cuando se trata de una función no periódica, definida en un intervalo, todavía se puede continuar periódicamente más allá del intervalo de interés, y expandir en series de Fourier, que todavía representaría fielmente esta función en el intervalo. Se puede entonces proceder de dos maneras:

  • Enfoque matemático es tomar el límite de un intervalo infinitamente amplio, en cuyo caso pasamos de la serie de Fourier a Transformación de Fourier . La transición no está exenta de dificultades, pero ha sido tratada con rigor en muchos libros de texto de matemáticas, especialmente en los de análisis complejo.
  • Enfoque físico es considerar un intervalo lo suficientemente grande como para que la discreción y otros efectos de las condiciones de contorno no importen -de hecho, en física siempre tenemos escalas de energía, tiempo, espacio y otras que limitan nuestra precisión-, ya sea debido a algunas propiedades del sistema (interacciones o efectos ignorados en el modelo, aka la temperatura cero está por encima de la temperatura Kondo ) o debido a las limitaciones inherentes a las mediciones. Aunque esto pueda parecer menos riguroso que el razonamiento matemático, está destinado a darnos resultados (físicamente) correctos.

Observaciones: el razonamiento anterior es fácilmente aplicable a los sistemas de celosía, donde la transformada de Fourier se convierte en transformada discreta de Fourier .

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