Densidades de carga de polarización
De Laplace
Contenido |
1 Introducción
2 Potencial debido a una polarización
3 Densidades de carga de polarización
La expresión del potencial debido a un volumen polarizado puede transformarse en suma de dos integrales
![\phi(\mathbf{r}) =\frac{1}{4\pi\varepsilon_0}\int_\tau \mathbf{P}(\mathbf{r}'){\cdot}\frac{(\mathbf{r}-\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|^3}\,\mathrm{d}\tau' = \frac{1}{4\pi\varepsilon_0}\int_\tau\frac{\rho_p}{|\mathbf{r}-\mathbf{r}'|}\,\mathrm{d}\tau'+
\frac{1}{4\pi\varepsilon_0}\oint_{\partial\tau}\frac{\sigma_p}{|\mathbf{r}-\mathbf{r}'|}\,\mathrm{d}S'](/wiki/images/math/8/1/f/81fe70f4764ea86f10d178dc64eeedb3.png)
donde
![\rho_p(\mathbf{r}')=-\nabla'\cdot\mathbf{P}(\mathbf{r}')](/wiki/images/math/f/e/6/fe60bd8842dd3fda0c0457fa020956f3.png)
![\sigma_p(\mathbf{r}')=\mathbf{n}'\cdot\mathbf{P(\mathbf{r}')}](/wiki/images/math/8/6/f/86f24c8221d1f49f6d72a1ee59b70546.png)
son las llamadas densidades de carga de polarización. Aquí es la normal exterior a la superficie del material polarizado. Mantenemos la prima en el operador nabla y en el vector normal (aunque fuera de la integral no es imprescindible), para recordar que estamos hablando de una distribución de fuentes y que el volumen de integración es aquél en que hay polarización.
La ventaja de esta transformación en suma de dos integrales es que este tipo de integrales corresponden al potencial debido a distribuciones de carga eléctrica. De esta forma se reduce un problema nuevo (el potencial debido a una distribución continua de dipolos) a uno ya conocido (el potencial de distribuciones de cargas)
En términos físicos, esto equivale a decir que una polarización es equivalente a una superposición de dos distribuciones de carga, una volumétrica y otra superficial.
3.1 Volumétrica
La densidad volumétrica de carga de polarización tiene la expresión general
![\rho(\mathbf{r})=-\nabla\cdot\mathbf{P}](/wiki/images/math/7/c/7/7c79ffdb4093c3d87bc75b37a495d6a8.png)
Esta densidad de carga solo aparece en el interior de los dieléctricos, pues en el vacío y su divergencia es nula.
3.2 Superficial
Dado un bloque de material dieléctrico de volumen τ, la densidad superficial de carga de polarización en su frontera viene dada por la expresión
![\sigma_p = \mathbf{P}\cdot\mathbf{n}](/wiki/images/math/b/f/2/bf2d76f70ac31070fc31f4d55715a6cb.png)
siendo la normal exterior al volumen.
Esta definición se puede generalizar al caso de que en lugar de una frontera exterior con el vacío tengamos una interfaz entre dos dieléctricos con polarizaciones y
. En ese caso, la densidad superficial de carga será la suma de las correspondientes a cada dieléctrico
![\sigma_p = \sigma_{p1}+\sigma_{p2}=\mathbf{P}_1\cdot\mathbf{n}_1+\mathbf{P}_2\cdot\mathbf{n}](/wiki/images/math/8/1/0/8109695084fbd8e03a96a0389e7d1d0b.png)
Si, como es habitual a la hora de considerar condiciones de salto, llamamos a la normal a la interfaz que va del medio 1 al 2, tenemos que
![\mathbf{n}_1=\mathbf{n}](/wiki/images/math/b/a/3/ba34b15c86d11865f2254f4c092e48c5.png)
![\mathbf{n}_2=-\mathbf{n}](/wiki/images/math/8/4/1/841e3711059017d221e2eb2dd972745d.png)
![\Rightarrow](/wiki/images/math/d/f/0/df09aea884019cb88a2957126faba316.png)
![\sigma_p = \mathbf{n}\cdot\left(\mathbf{P}_1-\mathbf{P}_2\right) = -\mathbf{n}\cdot[\mathbf{P}]](/wiki/images/math/d/3/8/d3882a37820adf0c5452645bc7d55fb3.png)
siendo
![[\mathbf{P}]=\mathbf{P}_2-\mathbf{P}_1\,](/wiki/images/math/a/a/d/aada47df93e070fdb340e16fdf9c3733.png)
el salto en la polarización al pasar de un medio al otro. Tenemos entonces la fórmula general
![\sigma_p = -\mathbf{n}\cdot[\mathbf{P}]\,](/wiki/images/math/e/8/b/e8b3c1c177edbf64c9632e7be2224dee.png)
que tiene una mayor analogía con la expresión correspondiente a la densidad volumétrica de carga.
3.3 Demostración
Para demostrar la equivalencia entre ambas expresiones, simplemente usamos la relación vectorial
![\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3} = \nabla'\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)](/wiki/images/math/f/5/f/f5f6d6b6e1d2ba02664cc30a78fd2a43.png)
donde la prima sobre el operador nabla significa que las derivadas son respecto a las coordenadas con prima (las posiciones de las fuentes).
De esta forma el potencial debido a la polarización se puede escribir
![\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\int_\tau\mathbf{P}(\mathbf{r'})\cdot\nabla'\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right)\,\mathrm{d}\tau'](/wiki/images/math/0/b/7/0b7ace2a6ab0293105296c8a1c3941da.png)
Usando ahora el álgebra del operador nabla
![\nabla\cdot(\psi\mathbf{A}) = \nabla\psi\cdot\mathbf{A}+\psi\nabla\cdot\mathbf{A}](/wiki/images/math/9/1/a/91a296939a2d1bbb0a03fe6de402bb74.png)
podemos transformar la integral en
![\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\left(\int_\tau\nabla'\cdot\left(\frac{\mathbf{P}(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\right)\,\mathrm{d}\tau'-\int_\tau \frac{\nabla'\cdot\mathbf{P}}{|\mathbf{r}-\mathbf{r'}|}\mathrm{d}\tau'\right)](/wiki/images/math/0/a/2/0a2cf6790c8c1bc449778bb59bac74c4.png)
Aplicando el teorema de Gauss a la primera integral la transformamos en una de superficie
![\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\left(\oint_{\partial\tau}\frac{\mathbf{P}(\mathbf{r}')\cdot\mathrm{d}\mathbf{S}'}{|\mathbf{r}-\mathbf{r}'|}\,\mathrm{d}\tau'-\int_\tau \frac{\nabla'\cdot\mathbf{P}}{|\mathbf{r}-\mathbf{r}'|}\mathrm{d}\tau'\right)](/wiki/images/math/a/0/6/a061333eeb1df8386f22a51f508ef090.png)
Aplicando que obtenemos finalmente la expresión deseada
![\phi(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0}\oint_{\partial\tau}\frac{\sigma_p(\mathbf{r}')\,\mathrm{d}S'}{|\mathbf{r}-\mathbf{r}'|}\,\mathrm{d}\tau'+\frac{1}{4\pi\varepsilon_0}\int_\tau \frac{\rho_p(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}\mathrm{d}\tau'](/wiki/images/math/6/1/7/617330c6e04b176d85bb21f0d0389502.png)
4 Ventajas de usar σp y ρp
La introducción de las densidades de carga de polarización supone que transformamos una integral de volumen en dos nuevas integrales: una de volumen y una de superficie.
Aunque esto no parece un avance, la ventaja está en que las dos integrales son de un tipo conocido: potenciales debidos a densidades de carga eléctrica. Esto permite emplear una serie de métodos alternativos para el cálculo del potencial y del campo eléctrico, aparte de la integración directa. Así puede emplearse la ley de Gauss en situaciones de alta simetría (como en el caso de una esfera polarizada radialmente), o la solución de la ecuación de Poisson.
5 ¿Son reales las cargas de polarización?
Ésta es hasta cierto punto una pregunta metafísica, ya que en física no se dice si algo es real o no, sino si un modelo matemático se corresponde con al realidad o no. Y como todas las preguntas de este tipo, depende de qué entendamos por cargas reales.
Si consideramos cargas “reales” aquellas que corresponde a partículas cargadas (como electrones o protones) entonces no, las cargas de polarización no están asociadas a cargas individuales depositadas en partículas, sino a una cierta redistribución en un sistema neutro como un átomo.
Sin embargo, esa definición de carga es muy restrictiva, ya que en muchos casos ni siquiera cuando tenemos un cuerpo verdaderamente cargado (un conductor cargado, por ejemplo) podemos asociar esa carga con una partícula en concreto.
En su lugar, es siempre preferible dar una definición operacional, a partir de los efectos medibles que produce. Una carga eléctrica sería una fuente escalar que produce un campo eléctrico (medible a partir de las fuerzas que provoca). Con esta definición práctica, las densidades de carga de polarización son tan reales como cualquier otra, pues sí que producen campos eléctricos.
6 Densidad de carga libre
7 Carga y momento de ρp y σp
Como a toda densidad de carga, a la de polarización se le puede calcular sus momentos multipolares, con el fin de aproximar el potencial que produce una distribución de dipolos en puntos alejados de ella.
7.1 Carga neta
La carga neta (momento monopolar) de una distribución de carga de polarización es siempre nula
![Q_p=\int_\tau \rho_p\,\mathrm{d}\tau+\oint_{\partial\tau}\sigma_p\,\mathrm{d}S\equiv 0](/wiki/images/math/3/f/4/3f496f86be1c7dd3e7f22d716195004e.png)
Este resultado es una consecuencia de que la distribución de carga de polarización es equivalente a un conjunto de dipolos. Puesto que cada dipolo es eléctricamente neutro, la carga total del sistema es nula.
Podemos demostrar este resultado a partir de las definiciones de ρp y σp, por aplicación del teorema de Gauss. Por un lado tenemos que
![\oint_{\partial\tau}\sigma_p\,\mathrm{d}S=\oint_{\partial\tau}\mathbf{P}\cdot\mathbf{n}\,\mathrm{d}S=\oint_{\partial\tau}\mathbf{P}\cdot\mathrm{d}\mathbf{S}](/wiki/images/math/4/e/0/4e0344479fac36b39ea0ddd970ec274f.png)
y por otro
![\int_{\tau}\rho_p\,\mathrm{d}\tau=-\int_{\tau}\nabla\cdot\mathbf{P}\,\mathrm{d}\tau=-\oint_{\partial\tau}\mathbf{P}\cdot\mathrm{d}\mathbf{S}](/wiki/images/math/8/3/0/8305aa42a6302e33f9b9bb46bfe0f420.png)
Sumando los dos términos
![Q_p = \oint_{\partial\tau}\mathbf{P}\cdot\mathrm{d}\mathbf{S}-\oint_{\partial\tau}\mathbf{P}\cdot\mathrm{d}\mathbf{S} = 0](/wiki/images/math/9/6/8/968beaac227caa2f5eef081dc4dcc9f5.png)
7.2 Momento dipolar
El momento dipolar de la distribución de cargas de polarización equivale a la integral de la polarización
![\mathbf{p}= \int_\tau \rho_p\,\mathbf{r}\,\mathrm{d}\tau+\oint_{\partial\tau}\sigma_p\,\mathbf{r}\,\mathrm{d}S \equiv\int_\tau\mathbf{P}\,\mathrm{d}\tau](/wiki/images/math/c/9/5/c95c89d8a0c9550ecbeecfff0d685be0.png)
De nuevo este resultado es elemental: si la polarización es la densidad de momento dipolar, su integral nos dará el momento dipolar total.
La demostración es similar a la anterior pero, al tratarse de vectores, es un poco más complicada. Puede hacerse considerando por separado cada una de las componentes cartesianas del momento dipolar.