12 votos

¿Cuáles son las ventajas de la expansión multipolar de los potenciales?

Cuando veo las ecuaciones de las expansiones multipolares me "parecen" más duras que las expresiones originales. Por ejemplo:

Expansión multipolar - forma esférica, en Wikipedia

Apuesto a que no es así, pero quiero saber por qué. Por ejemplo, cuando se usan como integradores, ¿son estas expresiones ventajosas para la integración analítica o para la numérica?

Mi principal preocupación es sobre su uso en la mecánica cuántica (específicamente sobre las expansiones magnéticas multipolares). He encontrado muchas referencias sobre la matemática involucrada pero no puedo darme cuenta de por qué es útil. He leído la referencia:

Ideas de la química cuántica , Lucjan Piela (Elsevier, 2014). Anexo X ; Enlace a Google Books

pero no me queda claro "por qué hay una pérdida de la precisión" (como se dice en la primera página) al usar la expresión original y cómo el uso de la expansión multipolar evita eso.

Me parece que la ventaja está en un intercambio de variables pero no me queda claro.

Gracias anticipadas a cualquier ayuda.

6voto

Nathan Feger Puntos 7675

Para mí, las principales ventajas de las expansiones multipolares son típicamente conceptuales, en el sentido de que ayudan a encapsular y separar los efectos de la forma de la distribución de la carga y la posición de la carga de prueba en la interacción electrostática entre ambas.

Consideremos, por ejemplo, la interacción entre una carga de prueba en la posición $\mathbf r$ y una distribución de carga $\rho(\mathbf r)$ que está localizada dentro de una esfera de radio $R$ alrededor del origen. El potencial electrostático para tal situación viene dado, en general, por la integral del núcleo de Coulomb sobre la nube de carga, $$ V(\mathbf r)=\int_\Omega\frac{\rho(\mathbf r')\mathrm d\mathbf r'}{|\mathbf r-\mathbf r'|}. \tag 1 $$ Este es un objeto un poco grande y desordenado, ya que requiere llevar a cabo una larga y complicada integral numérica tridimensional cada vez que se quiere saber el valor de $V(\mathbf r)$ .

Entonces, ¿cómo ayuda una expansión multipolar? Voy a omitir todos los detalles de cómo se obtiene, pero es un resultado bastante estándar que cuando $r>R$ puede expresar el potencial en $(1)$ como $$ V(\mathbf r) =\int_\Omega\frac{\rho(\mathbf r')\mathrm d\mathbf r'}{|\mathbf r-\mathbf r'|} =\sum_{l=0}^\infty\sum_{m=-l}^l Q_{lm}\frac{Y_{lm}(\theta,\varphi)}{r^l} \tag 2 $$ donde $$ Q_{lm}=\frac{4\pi}{2l+1}\int Y_{lm}^*(\theta',\varphi')r'^l\rho(\mathbf r')\mathrm d\mathbf r' \tag 3 $$ son los momentos multipolares del sistema. A primera vista, esto parece tan feo como antes, y sólo he cambiado la integral desordenada por una serie desordenada, que todavía tendré que sumar numéricamente cada vez que quiera calcular $V(\mathbf r)$ .

El punto aquí es que en series multipolares como $(2)$ En el caso de la hipótesis de que la contribución de los términos de la ecuación sea menor que la de la ecuación de la ecuación, normalmente son sólo unos pocos términos los que contribuyen de forma significativa, lo que significa que sólo tengo que calcular unos pocos momentos multipolares. Esto es importante. Significa que puedo precalcular sólo algunas de las integrales numéricas tridimensionales completas en $(3)$ y todavía tienen un excelente dominio del potencial.

El punto de precomputación es el importante. He conseguido aislar las pocas cantidades de la distribución de carga - los momentos multipolares $Q_{lm}$ - que controlan el peso de los diferentes componentes multipolares del campo, y después de eso todo lo que necesito hacer es calcular unos cuantos potenciales (que son numérica y conceptualmente simples) y superponerlos para obtener la interacción completa.


Además de necesitar calcular menos integrales, estas integrales tienden a ser también más sencillas y benignas numéricamente, y esto enlaza con la preocupación expresada por Piela en el libro que enlazas. En particular, las integrales en $(1)$ a menudo implican muchos términos que casi se cancelan, lo que significa que se necesita una gran precisión en ambos $a$ y $b$ para obtener sólo una precisión mediocre en $a-b$ .

Para fundamentar esto en un ejemplo, considere dos cargos $\pm q$ en el $z$ eje en $\pm \delta=1\:\mathrm{mm}$ y su interacción con una partícula de prueba en el mismo eje a $z=1\:\mathrm{km}$ . Si quieres calcular el potencial electrostático allí, el equivalente de la integral en $(1)$ es sumar directamente las dos contribuciones: \begin {align} V &= \frac {q}{z+ \delta }- \frac {q}{z- \delta } = \frac {q}{1000.001\: \mathrm m}- \frac {q}{999.999\: \mathrm m} \\ &= \frac {q}{ \mathrm {m}} \left (0.00099999900000999 \cdots -0.00100000100000100 \cdots\right ) \\ & \approx -2 \times 10^{-9} \frac {q}{ \mathrm {m}} \end {align} Esto significa, en particular, que necesito los nueve decimales de cada uno de los potenciales -cinco cifras significativas- para obtener una sola cifra significativa de precisión en el resultado. Si esta es la única manera de hacer las cosas, entonces es una situación de fruncir el ceño y soportarlo, pero si quiero calcular el potencial en un montón de posiciones diferentes, entonces habrá que fruncir el ceño, sobre todo si en lugar de simples sumas tengo que estar haciendo integrales numéricas de alta precisión todo el tiempo, para obtener sólo un par de cifras significativas en mi potencial.

En cambio, si primero calculo el momento dipolar del sistema y uso el campo dipolar, sé que $$Q_{q,z}=q\delta-q(-\delta)=2q\delta=2q\times 1\:\mathrm{mm},$$ por lo que obtengo la misma cifra significativa de precisión para $V$ mientras se utiliza una sola cifra significativa de precisión para $z$ y $\delta$ . Este tipo de titulación numérica, en la que el nivel de precisión deseado en el resultado se corresponde con el nivel de precisión requerido en el cálculo, es el sello de un esquema numérico deseable. Además, sé que el siguiente término de la serie multipolar es del orden de $\delta/z\approx 10^{-6}$ con respecto al término dipolar, así que si quieres un resultado más exacto, entonces para las primeras cinco cifras significativas más o menos sabes que debes concentrarte en el término dipolar, y es sólo después de eso que necesitas incluir los términos de orden superior.


Esto me lleva a otra ventaja conceptual de las expansiones multipolares, y es la separación de escalas que está presente en $(2)$ Cada término de la escalera se escala como una potencia mayor de $1/r$ que su predecesor; además, la parte inferior $l$ Los términos tienden a tener una estructura angular muy "difusa", y sólo en los $l$ que el detalle angular entra en juego. Esta es la forma de precisar la observación de que, desde lejos, un sistema cargado es prácticamente indistinguible de una sola bola de carga, y que los detalles de la distribución de la carga sólo adquieren importancia al acercarse.

Esto establece una jerarquía de aproximaciones que es muy importante en términos de cómo conceptualizamos y aplicamos la electrodinámica en una amplia gama de contextos - desde la electrodinámica en medios continuos hasta la interacción de los átomos y las moléculas con la radiación, pasando por la descripción de las fuerzas intermoleculares, la óptica de enfoque de los microscopios y aceleradores de electrones y una serie de otros escenarios. En todos estos entornos, la separación de las interacciones electromagnéticas en diferentes órdenes multipolares aporta claridad conceptual a la importancia relativa de las diferentes partes de la interacción, y ese es su principal valor.

3voto

John R Ramsden Puntos 143

Como ejemplo, utilizaré principalmente la electrostática con su ecuación de Poisson $\Delta \phi=\rho$ (además, se omiten las constantes). Esto es fácilmente escalable al magnetismo estacionario ( $\Delta \vec{A}=\vec{j}$ ) o la gravedad ( $\Delta \Phi = 4 \pi \rho_{\rm matter}$ ).


Los potenciales multipolares suelen surgir "sin ningún esfuerzo" al tratar de resolver las ecuaciones de campo del vacío (sin fuente) por la método de separación de variables . Como tales, representan una familia infinita de " armónicos " que puede caracterizar completamente cualquier campo de vacío.

Sin embargo, un detalle importante es que estos armónicos no carecen completamente de fuente, sino que suelen tener singularidades en las que se localizan densidades infinitas de carga (o un tipo análogo de fuente).

Por ejemplo, los armónicos esféricos, que se obtienen por separación de la ecuación de Laplace en coordenadas radiales, tienen el primer "armónico monopolar" $\phi_{\rm m}\sim Q/r$ . Este armónico tiene $\Delta \phi_{\rm m}=0$ en todas partes excepto $r=0$ donde no sabemos realmente lo que está pasando porque el potencial diverge. Un análisis más sofisticado nos da que $\Delta \phi_{\rm m}\sim Q\delta$ , donde $\delta$ es el pico infinito Función delta de Dirac .

El dipolo $\sim 1/r^2$ correspondería entonces a dos picos en dirección opuesta uno al lado del otro (o dos cargas puntuales de signo opuesto una al lado de la otra) lo que se caracteriza por la "derivada de Dirac-delta" $\delta'$ . Es posible encontrar una caracterización similar en términos de fuentes "finas" o "puntuales" para cada potencial multipolar. Es decir, los multipolos corresponden de hecho a fuentes idealizadas infinitamente empaquetadas, en el caso de los multipolos esféricos alrededor de $r=0$ .


Pero ahora considera lo siguiente: Desde cualquier El campo del vacío se caracteriza por los armónicos, al entender los armónicos, se entiende automáticamente el comportamiento típico del campo. Si situamos el origen en algún lugar del centro de la nube de carga y sólo estudiamos el campo fuera de ella, sabremos que el campo exterior está representado por una serie de funciones proporcionales a $1/r^n$ . Si estamos muy lejos de la nube ( $r$ muy grande), un $\sim 1/r^{10}$ probablemente sea bastante insignificante con respecto a la $\sim 1/r$ término. Esto significa que podemos aproximar bastante bien el campo exterior incluyendo sólo algunos primeros multipolos.


Pero, ¿cómo calculamos los coeficientes de los multipolos? Con el aumento del orden de los multipolos, esto puede resultar cada vez más difícil, así que consideremos por ahora los dos primeros.

El coeficiente de orden principal correspondiente al monopolo $\sim 1/r$ es sólo la cantidad total de carga en la nube. Si, por ejemplo, la cantidad total de carga en la nube es cero (como es el caso de un átomo o molécula neutra), se siempre saber que el campo lejos de la nube caerá al menos tan $1/r^2$ .

Para calcular el dipolar $\sim 1/r^2$ primero hay que calcular el vector dipolo $$\vec{d} = \sum q_i \vec{r}_i$$ donde $q_i$ son cargos de fuentes individuales en la nube y $\vec{r}_i$ sus posiciones. El coeficiente del dipolo es entonces proporcional a $\vec{d}\cdot \vec{r}$ , donde $\vec{r}$ es el punto en el que estamos evaluando el potencial. Es decir, si se tiene una nube de 1000 cargas con carga total cero, basta con ejecutar esta única suma y se obtiene una aproximación fiel del potencial electrostático lejos de la nube. Esto puede resultar muy útil.


En cuanto a la citada utilidad en el Ideas de la química cuántica La cuestión es que si representamos los números en un ordenador con una cantidad finita de dígitos, la sustracción de números cercanos reduce el número de dígitos significativos. Consideremos, por ejemplo, una representación de 5 dígitos. Ahí, se podría decir que por ejemplo $1.0000$ representa de hecho todo el número de $0.99995$ a $1.000049999...$ es decir, cuando el ordenador le da $1.0000$ es de hecho $1.0000 \pm 0.00005$ . Consideremos ahora la sustracción $$1.0001-1.0000 = 0.0001$$ deberíamos añadir barras de error $$1.0001\pm 0.00005-1.0000\pm 0.00005 = 0.0001\pm 0.0001$$ Es decir, ¡¡¡el error relativo del resultado es del 100%!!!

Por ello, a menudo puede ser útil evitar la evaluación directa de las sustracciones de números cercanos utilizando varias aproximaciones. Una de ellas es utilizar un potencial dipolar para el potencial de cargas muy cercanas e iguales.

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