10 votos

Ayudar a que el cálculo de una integral de la función de Green de un Laplaciano discreto en un cuadrado de celosía

Necesito calcular la siguiente integral: $$ \int_0^1 \int_0^1 \frac{1-\cos(2 \pi k_1 x) \cos(2 \pi k_2 y)}{4 \sin(\pi k_1)^2 + 4 \sin( \pi k_2)^2} dk_1 dk_2 $$ He intentado utilizar algunos de contorno integral de las técnicas, pero sin mucho éxito. Cualquier ayuda será apreciada.

2voto

Evan Anderson Puntos 118832

De acuerdo a estos dos artículos: Celosía de funciones de Green en todas las dimensiones, y Discretos dispersión de la teoría: la función de Green para una plaza de celosía, la entrada del Laplaciano de la matriz de las dos dimensiones de la plaza de celosía: $$ G_{xy} = \frac{1}{2}\int_0^1 \int_0^1 \frac{1-\cos(2 \pi k_1 x) \cos(2 \pi k_2 y)}{2 - \cos(2 \pi k_1 ) - \cos(2 \pi k_2 )} dk_1 dk_2\\ = \frac{1}{4}\int_0^1 \int_0^1 \frac{1-\cos(2 \pi k_1 x) \cos(2 \pi k_2 y)}{ \sin^2(\pi k_1) + \sin^2( \pi k_2) } dk_1 dk_2. $$ OP de la fórmula parece un poco engañoso porque la plaza de seno es en el paréntesis.

Al $x = y = n$, $$ G_{xy} =\frac{1}{4}\int_0^1 \int_0^1 \frac{1-\cos(2 \pi k_1 n) \cos(2 \pi k_2 n)}{ \sin^2(\pi k_1) + \sin^2( \pi k_2) } dk_1 dk_2 \\ = \frac{1}{4}\int_0^1 \int_0^1 \frac{2- \cos\big(2 n\pi (k_1-k_2)\big)- \cos\big(2 n\pi (k_1+k_2) \big)}{ 2 - \cos(2 \pi k_1 ) - \cos(2 \pi k_2 ) } dk_1 dk_2 \\ \frac{1}{4}\int_0^1 \int_0^1 \frac{2- \cos\big(2 n\pi (k_1-k_2)\big)- \cos\big(2 n\pi (k_1+k_2) \big)}{ 2 - 2 \cos\big(\pi (k_1+k_2) \big) \cos\big( \pi (k_1-k_2) \big) } dk_1 dk_2. $$ Vamos $u = k_1 + k_2\in (0,2)$, $v = k_1- k_2\in (-1,1)$, $|\det \partial(k_1,k_2)/\partial(u,v)| = 1/2 $, por lo tanto $$ \begin{aligned} G_{n,n} &= \frac{1}{16} \int^2_0 \left(\int^{1}_{-1}\frac{2- \cos(2 n\pi v)- \cos(2 n\pi u)}{ 1 - \cos (\pi u) \cos( \pi v ) } dv\right)du \\ &= \frac{1}{16} \int^{1}_{-1} \left(\int^{1}_{-1}\frac{2- \cos(2 n\pi v)- \cos(2 n\pi u)}{ 1 + \cos (\pi u) \cos( \pi v ) } dv\right)du \\ &=\frac{1}{16} \int^{1}_{-1} \left(\int^{1}_{-1}\frac{2- 2\cos(2 n\pi u)}{ 1 + \cos (\pi u) \cos( \pi v ) } dv\right)du && \text{(By symmetry)} \\ &=\frac{1}{8} \int^{1}_{-1}\frac{1- \cos(2 n\pi u)}{ \sin(\pi u)} du \\ &=\frac{1}{4\pi} \int^{1}_{-1}\frac{ \sin^2( n\pi w)}{ \sin^2(\pi w)} d(\cos\pi w ) && \text{(Let %#%#%)} \\ &= \frac{1}{4\pi} \int^{1}_{-1} \left(\frac{\sin\big(n \arccos x \big)}{\sqrt{1-x^2}} \right)^2 dx. \end{aligned} $$ Y este polinomio de Chebyshev de la segunda clase $w= -u$. Por $$ \begin{aligned} &U^2_n(x) - U^2_{n-1}(x) = \big(U_n(x) + U_{n-1}(x)\big)\big(U_n(x) - U_{n-1}(x)\big) \\ =&\frac{\sin\big((n+1)w\big) + \sin\big(nw\big)}{\sin w}\cdot \frac{\sin\big((n+1)w\big) - \sin\big(nw\big)}{\sin w} \\ =&\frac{\sin\big((2n+1)w/2\big)\cos (w/2 )}{\cos (w/2 )\sin (w/2 )}\cdot \frac{\sin (w/2 )\cos\big((2n+1)w/2\big)}{\cos (w/2 )\sin (w/2 )} \\ =& \frac{\sin\big((2n+1)w/2\big)\cos\big((2n+1)w/2\big)}{\cos (w/2 )\sin (w/2 )} \\ =& \frac{\sin\big((2n+1)w/2\big)}{\sin w} = U_{2n}(x). \end{aligned} $$ Por lo tanto $$ G_{nn} = \frac{1}{2\pi} \int^{1}_{0} \sum^{n}_{k=1} U_{2k-2}(x)dx = \frac{1}{2\pi} \sum^{n}_{k=1} \frac{1}{2k-1}. $$ Ahora por la relación de recurrencia se puede pasar a calcular $U_{n-1}(x)$ de los casos. Por favor, consulte la referencia de que el segundo documento en el que se enumeran los puntos a: Funciones de Green en la Física Cuántica. Mi biblioteca me permite tener una vista previa gratuita (o Springer sólo la tiene), por favor refiérase a la fórmula (5.39) en el enlace de arriba, el autor utiliza una integral elíptica para representar a $x = n,y =n+1$ $G_{n,n+1}$ llamó a $G_{n,n}$$G_{1,z}$, así que aquí yo, básicamente, se trató de derivar una más explícita fórmula para $y-x=1$.

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