¿La ecuación de Poisson 1-D tiene potenciales monótonos si ρ=ρ(ϕ(z))ρ=ρ(ϕ(z))\rho=\rho(\phi(z))?

Estoy resolviendo la ecuación de Poisson 1-D:

d 2 ϕ d z 2 = 4 π ρ ( ϕ )

con el requisito adicional de que ρ ( ϕ ( z = 0 ) ) = 0 . Si comienzo multiplicando cada lado por d ϕ d z e integramos de 0 a ϕ yo obtengo

1 2 ( d ϕ d z ) 2 = 4 π 0 ϕ ρ ( ϕ ) d ϕ d z d z

así que tengo

1 2 ( d ϕ d z ) 2 = 4 π 0 ϕ ρ ( ϕ ) d ϕ

Debido al cuadrado del lado izquierdo, solo obtendré soluciones cuando el lado derecho sea positivo. Encuentro este requisito extraño porque siempre obtendrás que la derivada del potencial nunca cambia de signo y siempre es monótona cuando la densidad de carga se escribe en la forma ρ ( ϕ ) . ¿Estoy haciendo algo mal aquí o he malinterpretado el resultado?

Más específicamente, si proporciono una densidad de carga que cambia de signo con ϕ ¿Seré incapaz de encontrar una solución?

También para cualquiera que esté familiarizado con esto, esta primera serie de pasos se usa en la derivación de la ley de Child-Langmuir para el espesor de una cubierta de plasma, pero en el problema de la cubierta de iones ρ no cambia de signo (Ver el Criterio de Bohm si está interesado).

Sugerencia para la pregunta (v1): si asumimos la primera ecuación, entonces falta una 'constante' de integración en la segunda ecuación. Expresado de manera equivalente: falta el límite inferior de la antiderivada.
Si tengo una integral definida, ¿por qué hay una constante de integración?
Ah, probablemente olvidé mencionar eso. ρ ( 0 ) = 0 lo cambiaré
No hay necesidad de asumir ρ ( 0 ) = 0 . El problema está en otra parte.
¿Podrías explicar un poco más? La última vez que hice esto en clase, multiplicamos por d ϕ / d z e integrado por partes por lo que después de que el LHS debe ser 1 / 2 [ ( d ϕ / d z ) ϕ = ϕ 2 ( d ϕ / d z ) ϕ = 0 2 ]
Para ρ ( ϕ ) := 0 la solucion es ϕ ( z ) = a · z + b y la ecuacion dice 0 = 0 . Dices que te "integras" y terminas con 1 2 a 2 = 0 . Magia.
Esta no es la ecuación de Poisson. En la ecuación de Poisson ρ es una función dada y no una función de ϕ .
@NikolajK El requisito no era ρ ( ϕ ) = 0 fue ρ ( ϕ = 0 ) = 0
@Anode: ¿Entonces? ρ ( ϕ ) := 0 cumple ρ ( ϕ = 0 ) = 0 . Entonces estoy señalando que para el caso trivial, la relación concluida ya está rota. Esa es una herramienta general para investigar matemáticas: si es escéptico de un resultado obtenido con parámetros libres, regrese y elija ejemplos particulares para esos parámetros para ver en qué punto de la derivación el sistema falla.
@V.Moretti ϕ es una función libre de z que estoy resolviendo dada una densidad de carga que depende del potencial.
@NikolajK ¿Qué relación concluida se rompe? En el caso trivial \rho no cambia de signo. no es de extrañar que ϕ ( z ) = 0 en tu ejemplo
Permítanme ampliar el comentario de Qmechanic :
0 z 0 d d z ( d ϕ d z ) 2 d z ( d ϕ d z ) 2
En cambio
0 z 0 d d z ( d ϕ d z ) 2 d z = ( d ϕ d z ) 2 | 0 z 0
La expresión "incorrecta" parece "manifiestamente positiva". El correcto no.
@WillieWong ¿No es esto lo que escribí en un comentario anterior?
@Ánodo: no. Usted evaluó en ϕ = 0 , que es una tontería. Integrando una derivada en z debes evaluar en z = z 0 y z = 0 . Por cierto, también hace su elección de límite superior de integración " ϕ " se ve muy extraño. (Además, en este nivel es un problema de cálculo , y tiene relativamente poco que ver con la física ).
@WillieWong Ah, veo el problema. debería haber escrito ρ ( ϕ ( z = 0 ) ) = 0 .
@Anode: "No sorprende que ϕ ( z ) = 0 en tu ejemplo". ¿Qué quieres decir con eso? La solución es ϕ ( z ) = a · z + b y estos son validos para todos a , b . La conclusión que debes sacar no es que ϕ ( z ) = 0 sigue.

Respuestas (1)

El sistema resuelve las ecuaciones de Euler-Lagrange del Lagrangiano:

L ( ϕ , d ϕ d z ) = 1 2 ( d ϕ d z ) 2 4 π 0 ϕ ( z ) ρ ( s ) d s
ya que las ecuaciones de Lagrange correspondientes son:
d 2 ϕ d z 2 = 4 π ρ ( ϕ ) .
Como el lagrangiano no depende explícitamente del "tiempo" z , la "energía" se conserva a lo largo de las soluciones (variando z ):
1 2 ( d ϕ d z ) 2 + 4 π 0 ϕ ( z ) ρ ( s ) d s = mi .
El constante mi se puede determinar a partir de las condiciones iniciales, por ejemplo ϕ ( 0 ) y d ϕ / d z | z = 0 . Una vez mi ha sido fijada, la ecuación se convierte en:
d ϕ d z = ± 2 mi 8 π 0 ϕ ( z ) ρ ( s ) d s .
El signo, de nuevo puede ser fijo desde la condición inicial ya que es la de d ϕ / d z | z = 0 .

Me parece que, al menos localmente, esa ecuación siempre se puede resolver si ρ es continuo, a menos que d ϕ / d z | z = 0 = 0 , situación que requiere mucho cuidado ya que el RHS (mientras siga siendo continuo) no es localmente Lipschitz en torno a la correspondiente condición inicial ϕ ( 0 ) .

¿Se supone que el signo en la primera línea es un +?
Tal vez me equivoque con los signos, pero no lo creo: d / d z ( L / ( d ϕ / d z ) ) = L / ϕ produce tu ecuación...