Función de Green (teoría de muchos cuerpos)

Correladores de operadores de campo.

En la teoría de muchos cuerpos , el término función de Green (o función de Green ) se utiliza a veces indistintamente con función de correlación , pero se refiere específicamente a correlacionadores de operadores de campo u operadores de creación y aniquilación .

El nombre proviene de las funciones de Green utilizadas para resolver ecuaciones diferenciales no homogéneas , con las que están vagamente relacionadas. (Específicamente, sólo las 'funciones de Green' de dos puntos en el caso de un sistema que no interactúa son funciones de Green en el sentido matemático; el operador lineal que invierten es el operador hamiltoniano , que en el caso que no interactúa es cuadrático en el caso de un sistema que no interactúa. campos.)

Caso espacialmente uniforme

Definiciones basicas

Consideramos una teoría de muchos cuerpos con operador de campo (operador de aniquilación escrito en base a la posición) . ${\displaystyle \psi (\mathbf {x} )}$

Los operadores de Heisenberg se pueden escribir en términos de operadores de Schrödinger como

$${\displaystyle \psi (\mathbf {x} ,t)=e^{iKt}\psi (\mathbf {x} )e^{-iKt},}$$

gran ${\displaystyle {\bar {\psi }}(\mathbf {x} ,t)=[\psi (\mathbf {x} ,t)]^{\dagger }}$ ${\displaystyle K=H-\mu N}$

De manera similar, para los operadores de tiempo imaginario ,

$${\displaystyle \psi (\mathbf {x} ,\tau )=e^{K\tau }\psi (\mathbf {x} )e^{-K\tau }}$$

$${\displaystyle {\bar {\psi }}(\mathbf {x} ,\tau )=e^{K\tau }\psi ^{\dagger }(\mathbf {x} )e^{-K\tau }.}$$

conjugado hermitiano ${\displaystyle {\bar {\psi }}(\mathbf {x} ,\tau )}$ ${\displaystyle \psi (\mathbf {x} ,\tau )}$

En tiempo real, la función -point Green está definida por ${\displaystyle 2n}$

$${\displaystyle G^{(n)}(1\ldots n\mid 1'\ldots n')=i^{n}\langle T\psi (1)\ldots \psi (n){\bar {\psi }}(n')\ldots {\bar {\psi }}(1')\rangle ,}$$

ordenamiento temporal ${\displaystyle j}$ ${\displaystyle (\mathbf {x} _{j},t_{j})}$ ${\displaystyle j'}$ ${\displaystyle (\mathbf {x} _{j}',t_{j}')}$ ${\displaystyle T}$

En tiempo imaginario, la definición correspondiente es

$${\displaystyle {\mathcal {G}}^{(n)}(1\ldots n\mid 1'\ldots n')=\langle T\psi (1)\ldots \psi (n){\bar {\psi }}(n')\ldots {\bar {\psi }}(1')\rangle ,}$$

${\displaystyle j}$ ${\displaystyle \mathbf {x} _{j},\tau _{j}}$ ${\displaystyle \tau _{j}}$ ${\displaystyle 0}$ ${\textstyle \beta ={\frac {1}{k_{\text{B}}T}}}$

Nota sobre los signos y la normalización utilizados en estas definiciones: Los signos de las funciones de Green se han elegido de modo que la transformada de Fourier de la función Green térmica de dos puntos ( ) para una partícula libre sea ${\displaystyle n=1}$

$${\displaystyle {\mathcal {G}}(\mathbf {k} ,\omega _{n})={\frac {1}{-i\omega _{n}+\xi _{\mathbf {k} }}},}$$

$${\displaystyle G^{\mathrm {R} }(\mathbf {k} ,\omega )={\frac {1}{-(\omega +i\eta )+\xi _{\mathbf {k} }}},}$$

$${\displaystyle \omega _{n}={\frac {[2n+\theta (-\zeta )]\pi }{\beta }}}$$

frecuencia de Matsubara

En todo momento, es para bosones y fermiones y denota un conmutador o anticonmutador , según corresponda. ${\displaystyle \zeta }$ ${\displaystyle +1}$ ${\displaystyle -1}$ ${\displaystyle [\ldots ,\ldots ]=[\ldots ,\ldots ]_{-\zeta }}$

(Consulte a continuación para obtener más detalles).

Funciones de dos puntos

La función de Green con un solo par de argumentos ( ) se conoce como función de dos puntos o propagador . En presencia de simetría traslacional tanto espacial como temporal, depende sólo de la diferencia de sus argumentos. Tomando la transformada de Fourier con respecto al espacio y al tiempo se obtiene ${\displaystyle n=1}$

$${\displaystyle {\mathcal {G}}(\mathbf {x} \tau \mid \mathbf {x} '\tau ')=\int _{\mathbf {k} }d\mathbf {k} {\frac {1}{\beta }}\sum _{\omega _{n}}{\mathcal {G}}(\mathbf {k} ,\omega _{n})e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')-i\omega _{n}(\tau -\tau ')},}$$

donde la suma está sobre las frecuencias ${\displaystyle (L/2\pi )^{d}}$

En tiempo real, indicaremos explícitamente la función ordenada en el tiempo con un superíndice T:

$${\displaystyle G^{\mathrm {T} }(\mathbf {x} t\mid \mathbf {x} 't')=\int _{\mathbf {k} }d\mathbf {k} \int {\frac {d\omega }{2\pi }}G^{\mathrm {T} }(\mathbf {k} ,\omega )e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')-i\omega (t-t')}.}$$

La función de Green de dos puntos en tiempo real se puede escribir en términos de funciones de Green "retardadas" y "avanzadas", que resultarán tener propiedades de analiticidad más simples. Las funciones Green retardada y avanzada están definidas por

$${\displaystyle G^{\mathrm {R} }(\mathbf {x} t\mid \mathbf {x} 't')=-i\langle [\psi (\mathbf {x} ,t),{\bar {\psi }}(\mathbf {x} ',t')]_{\zeta }\rangle \Theta (t-t')}$$

$${\displaystyle G^{\mathrm {A} }(\mathbf {x} t\mid \mathbf {x} 't')=i\langle [\psi (\mathbf {x} ,t),{\bar {\psi }}(\mathbf {x} ',t')]_{\zeta }\rangle \Theta (t'-t),}$$

Están relacionados con la función Verde ordenada en el tiempo por

$${\displaystyle G^{\mathrm {T} }(\mathbf {k} ,\omega )=[1+\zeta n(\omega )]G^{\mathrm {R} }(\mathbf {k} ,\omega )-\zeta n(\omega )G^{\mathrm {A} }(\mathbf {k} ,\omega ),}$$

$${\displaystyle n(\omega )={\frac {1}{e^{\beta \omega }-\zeta }}}$$

Bose-EinsteinFermi-Dirac

Ordenamiento en tiempo imaginario y periodicidad β

Las funciones térmicas de Green se definen sólo cuando ambos argumentos de tiempo imaginario están dentro del rango de . La función Green de dos puntos tiene las siguientes propiedades. (Los argumentos de posición o impulso se suprimen en esta sección). ${\displaystyle 0}$ ${\displaystyle \beta }$

En primer lugar, depende sólo de la diferencia de los tiempos imaginarios:

$${\displaystyle {\mathcal {G}}(\tau ,\tau ')={\mathcal {G}}(\tau -\tau ').}$$

${\displaystyle \tau -\tau '}$ ${\displaystyle -\beta }$ ${\displaystyle \beta }$

En segundo lugar, es (anti)periódico bajo cambios de . Debido al pequeño dominio dentro del cual se define la función, esto significa simplemente ${\displaystyle {\mathcal {G}}(\tau )}$ ${\displaystyle \beta }$

$${\displaystyle {\mathcal {G}}(\tau -\beta )=\zeta {\mathcal {G}}(\tau ),}$$

${\displaystyle 0<\tau <\beta }$

Estas dos propiedades permiten la representación de la transformada de Fourier y su inversa,

$${\displaystyle {\mathcal {G}}(\omega _{n})=\int _{0}^{\beta }d\tau \,{\mathcal {G}}(\tau )\,e^{i\omega _{n}\tau }.}$$

Finalmente, observe que tiene una discontinuidad en ; esto es consistente con un comportamiento de larga distancia de . ${\displaystyle {\mathcal {G}}(\tau )}$ ${\displaystyle \tau =0}$ ${\displaystyle {\mathcal {G}}(\omega _{n})\sim 1/|\omega _{n}|}$

Representación espectral

Los propagadores en tiempo real e imaginario pueden relacionarse con la densidad espectral (o peso espectral), dada por

$${\displaystyle \rho (\mathbf {k} ,\omega )={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}2\pi \delta (E_{\alpha }-E_{\alpha '}-\omega )|\langle \alpha \mid \psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle |^{2}\left(e^{-\beta E_{\alpha '}}-\zeta e^{-\beta E_{\alpha }}\right),}$$

| α ⟩H − μNE _α

El propagador de tiempo imaginario viene dado entonces por

$${\displaystyle {\mathcal {G}}(\mathbf {k} ,\omega _{n})=\int _{-\infty }^{\infty }{\frac {d\omega '}{2\pi }}{\frac {\rho (\mathbf {k} ,\omega ')}{-i\omega _{n}+\omega '}}~,}$$

propagador

$${\displaystyle G^{\mathrm {R} }(\mathbf {k} ,\omega )=\int _{-\infty }^{\infty }{\frac {d\omega '}{2\pi }}{\frac {\rho (\mathbf {k} ,\omega ')}{-(\omega +i\eta )+\omega '}},}$$

${\displaystyle \eta \to 0^{+}}$

El propagador avanzado viene dado por la misma expresión, pero con el denominador. ${\displaystyle -i\eta }$

La función ordenada en el tiempo se puede encontrar en términos de y . Como se afirmó anteriormente, y tienen propiedades de analiticidad simples: el primero (último) tiene todos sus polos y discontinuidades en el semiplano inferior (superior). ${\displaystyle G^{\mathrm {R} }}$ ${\displaystyle G^{\mathrm {A} }}$ ${\displaystyle G^{\mathrm {R} }(\omega )}$ ${\displaystyle G^{\mathrm {A} }(\omega )}$

El propagador térmico tiene todos sus polos y discontinuidades en el eje imaginario. ${\displaystyle {\mathcal {G}}(\omega _{n})}$ ${\displaystyle \omega _{n}}$

La densidad espectral se puede encontrar de forma muy sencilla a partir de , utilizando el teorema de Sokhatsky-Weierstrass. ${\displaystyle G^{\mathrm {R} }}$

$${\displaystyle \lim _{\eta \to 0^{+}}{\frac {1}{x\pm i\eta }}=P{\frac {1}{x}}\mp i\pi \delta (x),}$$

Pparte principal de Cauchy

$${\displaystyle \rho (\mathbf {k} ,\omega )=2\operatorname {Im} G^{\mathrm {R} }(\mathbf {k} ,\omega ).}$$

Esto implica además que obedece a la siguiente relación entre sus partes real e imaginaria: ${\displaystyle G^{\mathrm {R} }(\mathbf {k} ,\omega )}$

$${\displaystyle \operatorname {Re} G^{\mathrm {R} }(\mathbf {k} ,\omega )=-2P\int _{-\infty }^{\infty }{\frac {d\omega '}{2\pi }}{\frac {\operatorname {Im} G^{\mathrm {R} }(\mathbf {k} ,\omega ')}{\omega -\omega '}},}$$

${\displaystyle P}$

La densidad espectral obedece a una regla de suma,

$${\displaystyle \int _{-\infty }^{\infty }{\frac {d\omega }{2\pi }}\rho (\mathbf {k} ,\omega )=1,}$$

$${\displaystyle G^{\mathrm {R} }(\omega )\sim {\frac {1}{|\omega |}}}$$

${\displaystyle |\omega |\to \infty }$

Transformada de Hilbert

La similitud de las representaciones espectrales de las funciones de Green en tiempo imaginario y real nos permite definir la función

$${\displaystyle G(\mathbf {k} ,z)=\int _{-\infty }^{\infty }{\frac {dx}{2\pi }}{\frac {\rho (\mathbf {k} ,x)}{-z+x}},}$$

${\displaystyle {\mathcal {G}}}$ ${\displaystyle G^{\mathrm {R} }}$

$${\displaystyle {\mathcal {G}}(\mathbf {k} ,\omega _{n})=G(\mathbf {k} ,i\omega _{n})}$$

$${\displaystyle G^{\mathrm {R} }(\mathbf {k} ,\omega )=G(\mathbf {k} ,\omega +i\eta ).}$$

${\displaystyle G^{\mathrm {A} }}$

La relación entre y se conoce como transformada de Hilbert . ${\displaystyle G(\mathbf {k} ,z)}$ ${\displaystyle \rho (\mathbf {k} ,x)}$

Prueba de representación espectral

Demostramos la prueba de la representación espectral del propagador en el caso de la función térmica de Green, definida como

$${\displaystyle {\mathcal {G}}(\mathbf {x} ,\tau \mid \mathbf {x} ',\tau ')=\langle T\psi (\mathbf {x} ,\tau ){\bar {\psi }}(\mathbf {x} ',\tau ')\rangle .}$$

Debido a la simetría traslacional, sólo es necesario considerar para , dado por ${\displaystyle {\mathcal {G}}(\mathbf {x} ,\tau \mid \mathbf {0} ,0)}$ ${\displaystyle \tau >0}$

$${\displaystyle {\mathcal {G}}(\mathbf {x} ,\tau \mid \mathbf {0} ,0)={\frac {1}{\mathcal {Z}}}\sum _{\alpha '}e^{-\beta E_{\alpha '}}\langle \alpha '\mid \psi (\mathbf {x} ,\tau ){\bar {\psi }}(\mathbf {0} ,0)\mid \alpha '\rangle .}$$

$${\displaystyle {\mathcal {G}}(\mathbf {x} ,\tau \mid \mathbf {0} ,0)={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}e^{-\beta E_{\alpha '}}\langle \alpha '\mid \psi (\mathbf {x} ,\tau )\mid \alpha \rangle \langle \alpha \mid {\bar {\psi }}(\mathbf {0} ,0)\mid \alpha '\rangle .}$$

Dado que y son estados propios de , los operadores de Heisenberg se pueden reescribir en términos de operadores de Schrödinger, dando ${\displaystyle |\alpha \rangle }$ ${\displaystyle |\alpha '\rangle }$ ${\displaystyle H-\mu N}$

$${\displaystyle {\mathcal {G}}(\mathbf {x} ,\tau |\mathbf {0} ,0)={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}e^{-\beta E_{\alpha '}}e^{\tau (E_{\alpha '}-E_{\alpha })}\langle \alpha '\mid \psi (\mathbf {x} )\mid \alpha \rangle \langle \alpha \mid \psi ^{\dagger }(\mathbf {0} )\mid \alpha '\rangle .}$$

$${\displaystyle {\mathcal {G}}(\mathbf {k} ,\omega _{n})={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}e^{-\beta E_{\alpha '}}{\frac {1-\zeta e^{\beta (E_{\alpha '}-E_{\alpha })}}{-i\omega _{n}+E_{\alpha }-E_{\alpha '}}}\int _{\mathbf {k} '}d\mathbf {k} '\langle \alpha \mid \psi (\mathbf {k} )\mid \alpha '\rangle \langle \alpha '\mid \psi ^{\dagger }(\mathbf {k} ')\mid \alpha \rangle .}$$

La conservación del momento permite escribir el término final como (hasta posibles factores del volumen)

$${\displaystyle |\langle \alpha '\mid \psi ^{\dagger }(\mathbf {k} )\mid \alpha \rangle |^{2},}$$

La regla de la suma se puede probar considerando el valor esperado del conmutador,

$${\displaystyle 1={\frac {1}{\mathcal {Z}}}\sum _{\alpha }\langle \alpha \mid e^{-\beta (H-\mu N)}[\psi _{\mathbf {k} },\psi _{\mathbf {k} }^{\dagger }]_{-\zeta }\mid \alpha \rangle ,}$$

$${\displaystyle 1={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}e^{-\beta E_{\alpha }}\left(\langle \alpha \mid \psi _{\mathbf {k} }\mid \alpha '\rangle \langle \alpha '\mid \psi _{\mathbf {k} }^{\dagger }\mid \alpha \rangle -\zeta \langle \alpha \mid \psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle \langle \alpha '\mid \psi _{\mathbf {k} }\mid \alpha \rangle \right).}$$

Intercambiar las etiquetas en el primer término da

$${\displaystyle 1={\frac {1}{\mathcal {Z}}}\sum _{\alpha ,\alpha '}\left(e^{-\beta E_{\alpha '}}-\zeta e^{-\beta E_{\alpha }}\right)|\langle \alpha \mid \psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle |^{2}~,}$$

ρ

Caso que no interactúa

En el caso de no interacción, es un estado propio con energía (gran canónica) , donde se mide la relación de dispersión de una sola partícula con respecto al potencial químico. Por lo tanto, la densidad espectral se convierte en ${\displaystyle \psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle }$ ${\displaystyle E_{\alpha '}+\xi _{\mathbf {k} }}$ ${\displaystyle \xi _{\mathbf {k} }=\epsilon _{\mathbf {k} }-\mu }$

$${\displaystyle \rho _{0}(\mathbf {k} ,\omega )={\frac {1}{\mathcal {Z}}}\,2\pi \delta (\xi _{\mathbf {k} }-\omega )\sum _{\alpha '}\langle \alpha '\mid \psi _{\mathbf {k} }\psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle (1-\zeta e^{-\beta \xi _{\mathbf {k} }})e^{-\beta E_{\alpha '}}.}$$

De las relaciones de conmutación,

$${\displaystyle \langle \alpha '\mid \psi _{\mathbf {k} }\psi _{\mathbf {k} }^{\dagger }\mid \alpha '\rangle =\langle \alpha '\mid (1+\zeta \psi _{\mathbf {k} }^{\dagger }\psi _{\mathbf {k} })\mid \alpha '\rangle ,}$$

${\displaystyle [1+\zeta n(\xi _{\mathbf {k} })]{\mathcal {Z}}}$

$${\displaystyle \rho _{0}(\mathbf {k} ,\omega )=2\pi \delta (\xi _{\mathbf {k} }-\omega ).}$$

El propagador de tiempo imaginario es así

$${\displaystyle {\mathcal {G}}_{0}(\mathbf {k} ,\omega )={\frac {1}{-i\omega _{n}+\xi _{\mathbf {k} }}}}$$

$${\displaystyle G_{0}^{\mathrm {R} }(\mathbf {k} ,\omega )={\frac {1}{-(\omega +i\eta )+\xi _{\mathbf {k} }}}.}$$

Límite de temperatura cero

Como β → ∞ , la densidad espectral se vuelve

$${\displaystyle \rho (\mathbf {k} ,\omega )=2\pi \sum _{\alpha }\left[\delta (E_{\alpha }-E_{0}-\omega )\left|\left\langle \alpha \mid \psi _{\mathbf {k} }^{\dagger }\mid 0\right\rangle \right|^{2}-\zeta \delta (E_{0}-E_{\alpha }-\omega )\left|\left\langle 0\mid \psi _{\mathbf {k} }^{\dagger }\mid \alpha \right\rangle \right|^{2}\right]}$$

α = 0ω

Caso general

Definiciones basicas

Podemos usar 'operadores de campo' como arriba, u operadores de creación y aniquilación asociados con otros estados de una sola partícula, tal vez estados propios de la energía cinética (que no interactúa). Luego usamos

$${\displaystyle \psi (\mathbf {x} ,\tau )=\varphi _{\alpha }(\mathbf {x} )\psi _{\alpha }(\tau ),}$$

${\displaystyle \psi _{\alpha }}$ ${\displaystyle \alpha }$ ${\displaystyle \varphi _{\alpha }(\mathbf {x} )}$

$${\displaystyle {\mathcal {G}}_{\alpha _{1}\ldots \alpha _{n}|\beta _{1}\ldots \beta _{n}}^{(n)}(\tau _{1}\ldots \tau _{n}|\tau _{1}'\ldots \tau _{n}')=\langle T\psi _{\alpha _{1}}(\tau _{1})\ldots \psi _{\alpha _{n}}(\tau _{n}){\bar {\psi }}_{\beta _{n}}(\tau _{n}')\ldots {\bar {\psi }}_{\beta _{1}}(\tau _{1}')\rangle }$$

${\displaystyle G^{(n)}}$

Funciones de dos puntos

Éstos dependen sólo de la diferencia de sus argumentos temporales, de modo que

$${\displaystyle {\mathcal {G}}_{\alpha \beta }(\tau \mid \tau ')={\frac {1}{\beta }}\sum _{\omega _{n}}{\mathcal {G}}_{\alpha \beta }(\omega _{n})\,e^{-i\omega _{n}(\tau -\tau ')}}$$

$${\displaystyle G_{\alpha \beta }(t\mid t')=\int _{-\infty }^{\infty }{\frac {d\omega }{2\pi }}\,G_{\alpha \beta }(\omega )\,e^{-i\omega (t-t')}.}$$

Podemos nuevamente definir funciones retrasadas y avanzadas de la manera obvia; estos están relacionados con la función ordenada en el tiempo de la misma manera que antes.

Las mismas propiedades de periodicidad descritas anteriormente se aplican a . Específicamente, ${\displaystyle {\mathcal {G}}_{\alpha \beta }}$

$${\displaystyle {\mathcal {G}}_{\alpha \beta }(\tau \mid \tau ')={\mathcal {G}}_{\alpha \beta }(\tau -\tau ')}$$

$${\displaystyle {\mathcal {G}}_{\alpha \beta }(\tau )={\mathcal {G}}_{\alpha \beta }(\tau +\beta ),}$$

${\displaystyle \tau <0}$

Representación espectral

En este caso,

$${\displaystyle \rho _{\alpha \beta }(\omega )={\frac {1}{\mathcal {Z}}}\sum _{m,n}2\pi \delta (E_{n}-E_{m}-\omega )\;\langle m\mid \psi _{\alpha }\mid n\rangle \langle n\mid \psi _{\beta }^{\dagger }\mid m\rangle \left(e^{-\beta E_{m}}-\zeta e^{-\beta E_{n}}\right),}$$

${\displaystyle m}$ ${\displaystyle n}$

Las expresiones de las funciones de Green se modifican de la forma obvia:

$${\displaystyle {\mathcal {G}}_{\alpha \beta }(\omega _{n})=\int _{-\infty }^{\infty }{\frac {d\omega '}{2\pi }}{\frac {\rho _{\alpha \beta }(\omega ')}{-i\omega _{n}+\omega '}}}$$

$${\displaystyle G_{\alpha \beta }^{\mathrm {R} }(\omega )=\int _{-\infty }^{\infty }{\frac {d\omega '}{2\pi }}{\frac {\rho _{\alpha \beta }(\omega ')}{-(\omega +i\eta )+\omega '}}.}$$

Sus propiedades de analiticidad son idénticas y definidas en el caso traslacionalmente invariante. La demostración sigue exactamente los mismos pasos, excepto que los dos elementos de la matriz ya no son conjugados complejos. ${\displaystyle {\mathcal {G}}(\mathbf {k} ,\omega _{n})}$ ${\displaystyle G^{\mathrm {R} }(\mathbf {k} ,\omega )}$

Caso de no interacción

Si los estados particulares de una sola partícula que se eligen son 'estados propios de energía de una sola partícula', es decir

$${\displaystyle [H-\mu N,\psi _{\alpha }^{\dagger }]=\xi _{\alpha }\psi _{\alpha }^{\dagger },}$$

${\displaystyle |n\rangle }$

$${\displaystyle (H-\mu N)\mid n\rangle =E_{n}\mid n\rangle ,}$$

${\displaystyle \psi _{\alpha }\mid n\rangle }$

$${\displaystyle (H-\mu N)\psi _{\alpha }\mid n\rangle =(E_{n}-\xi _{\alpha })\psi _{\alpha }\mid n\rangle ,}$$

${\displaystyle \psi _{\alpha }^{\dagger }\mid n\rangle }$

$${\displaystyle (H-\mu N)\psi _{\alpha }^{\dagger }\mid n\rangle =(E_{n}+\xi _{\alpha })\psi _{\alpha }^{\dagger }\mid n\rangle .}$$

Por lo tanto tenemos

$${\displaystyle \langle m\mid \psi _{\alpha }\mid n\rangle \langle n\mid \psi _{\beta }^{\dagger }\mid m\rangle =\delta _{\xi _{\alpha },\xi _{\beta }}\delta _{E_{n},E_{m}+\xi _{\alpha }}\langle m\mid \psi _{\alpha }\mid n\rangle \langle n\mid \psi _{\beta }^{\dagger }\mid m\rangle .}$$

Luego reescribimos

$${\displaystyle \rho _{\alpha \beta }(\omega )={\frac {1}{\mathcal {Z}}}\sum _{m,n}2\pi \delta (\xi _{\alpha }-\omega )\delta _{\xi _{\alpha },\xi _{\beta }}\langle m\mid \psi _{\alpha }\mid n\rangle \langle n\mid \psi _{\beta }^{\dagger }\mid m\rangle e^{-\beta E_{m}}\left(1-\zeta e^{-\beta \xi _{\alpha }}\right),}$$

$${\displaystyle \rho _{\alpha \beta }(\omega )={\frac {1}{\mathcal {Z}}}\sum _{m}2\pi \delta (\xi _{\alpha }-\omega )\delta _{\xi _{\alpha },\xi _{\beta }}\langle m\mid \psi _{\alpha }\psi _{\beta }^{\dagger }e^{-\beta (H-\mu N)}\mid m\rangle \left(1-\zeta e^{-\beta \xi _{\alpha }}\right),}$$

$${\displaystyle \langle m\mid \psi _{\alpha }\psi _{\beta }^{\dagger }\mid m\rangle =\delta _{\alpha ,\beta }\langle m\mid \zeta \psi _{\alpha }^{\dagger }\psi _{\alpha }+1\mid m\rangle }$$

Finalmente, la densidad espectral se simplifica para dar

$${\displaystyle \rho _{\alpha \beta }=2\pi \delta (\xi _{\alpha }-\omega )\delta _{\alpha \beta },}$$

$${\displaystyle {\mathcal {G}}_{\alpha \beta }(\omega _{n})={\frac {\delta _{\alpha \beta }}{-i\omega _{n}+\xi _{\beta }}}}$$

$${\displaystyle G_{\alpha \beta }(\omega )={\frac {\delta _{\alpha \beta }}{-(\omega +i\eta )+\xi _{\beta }}}.}$$

Ver también

• Teorema de fluctuación
• Relaciones Verde-Kubo
• Función de respuesta lineal
• ecuación de Lindblad
• Propagador
• Función de correlación (teoría cuántica de campos)
• Continuación analítica numérica

Referencias

Libros

• Bonch-Bruevich VL, Tyablikov SV (1962): El método de la función verde en mecánica estadística. Compañía editorial de Holanda del Norte.
• Abrikosov, AA, Gorkov, LP y Dzyaloshinski, IE (1963): Métodos de la teoría cuántica de campos en física estadística Englewood Cliffs: Prentice-Hall.
• Negele, JW y Orland, H. (1988): Sistemas cuánticos de muchas partículas AddisonWesley.
• Zubarev DN , Morozov V., Ropke G. (1996): Mecánica estadística de procesos de desequilibrio: conceptos básicos, teoría cinética (Vol. 1). John Wiley e hijos. ISBN  3-05-501708-0 .
• Mattuck Richard D. (1992), Guía de diagramas de Feynman en el problema de muchos cuerpos , Publicaciones de Dover, ISBN 0-486-67047-3 . 

Documentos

• Bogolyubov NN , Tyablikov SV Funciones verdes retrasadas y avanzadas en física estadística, Física soviética Doklady, vol. 4 , pág. 589 (1959).
• Zubarev DN , Funciones verdes de doble tiempo en física estadística, Física soviética Uspekhi 3 (3), 320–345 (1960).

enlaces externos

• Funciones de respuesta lineal en Eva Pavarini, Erik Koch, Dieter Vollhardt y Alexander Lichtenstein (eds.): DMFT at 25: Infinite Dimensions, Verlag des Forschungszentrum Jülich, 2014 ISBN 978-3-89336-953-9