Articulo de referencia

función de onda de Coulomb

Función de onda de Coulomb irregular G representada gráficamente de 0 a 20 con interacciones repulsivas y atractivas en Mathematica 13.1. Se ha añadido una imagen de un gráfico ...

Función de onda de Coulomb irregular G representada gráficamente de 0 a 20 con interacciones repulsivas y atractivas en Mathematica 13.1.
Función de onda de Coulomb irregular G representada gráficamente de 0 a 20 con interacciones repulsivas y atractivas en Mathematica 13.1.
Se ha añadido una imagen de un gráfico complejo de la función de onda de Coulomb regular.

En matemáticas , una función de onda de Coulomb es una solución de la ecuación de onda de Coulomb , que recibe su nombre de Charles-Augustin de Coulomb . Se utilizan para describir el comportamiento de partículas cargadas en un potencial de Coulomb y pueden expresarse en términos de funciones hipergeométricas confluentes o funciones de Whittaker de argumento imaginario.

ecuación de onda de Coulomb

La ecuación de onda de Coulomb para una sola partícula cargada de masametro{\displaystyle m}es la ecuación de Schrödinger con potencial de Coulomb [ 1 ]

(222metro+Zdoαr)ψk(r)=2k22metroψk(r),{\displaystyle \left(-\hbar ^{2}{\frac {\nabla ^{2}}{2m}}+{\frac {Z\hbar c\alpha }{r}}\right)\psi _{\vec {k}}({\vec {r}})={\frac {\hbar ^{2}k^{2}}{2m}}\psi _{\vec {k}}({\vec {r}})\,,}

dóndeZ=Z1Z2{\displaystyle Z=Z_{1}Z_{2}}es el producto de las cargas de la partícula y de la fuente del campo (en unidades de la carga elemental ,Z=1{\displaystyle Z=-1}para el átomo de hidrógeno),α{\displaystyle \alpha }es la constante de estructura fina y2k2/(2metro){\displaystyle \hbar ^{2}k^{2}/(2m)}es la energía de la partícula. La solución, que es la función de onda de Coulomb, se puede encontrar resolviendo esta ecuación en coordenadas parabólicas.

ξ=r+rk^,ζ=rrk^(k^=k/k).{\displaystyle \xi =r+{\vec {r}}\cdot {\hat {k}},\quad \zeta =r-{\vec {r}}\cdot {\hat {k}}\qquad ({\hat {k}}={\vec {k}}/k)\,.}

Dependiendo de las condiciones de contorno elegidas, la solución tiene diferentes formas. Dos de las soluciones son [ 2 ] [ 3 ]

ψk(±)(r)=Γ(1±iη)miπη/2miikrMETRO(iη,1,±ikrikr),{\displaystyle \psi _{\vec {k}}^{(\pm )}({\vec {r}})=\Gamma (1\pm i\eta )e^{-\pi \eta /2}e^{i{\vec {k}}\cdot {\vec {r}}}M(\mp i\eta ,1,\pm ikr-i{\vec {k}}\cdot {\vec {r}})\,,}

dóndeMETRO(a,b,z)1F1(a;b;z){\displaystyle M(a,b,z)\equiv {}_{1}\!F_{1}(a;b;z)}es la función hipergeométrica confluente ,η=Zmetrodoα/(k){\displaystyle \eta =Zmc\alpha /(\hbar k)}yΓ(z){\displaystyle \Gamma (z)}es la función gamma . Las dos condiciones de contorno utilizadas aquí son:

ψk(±)(r)miikr(kr±),{\displaystyle \psi _{\vec {k}}^{(\pm )}({\vec {r}})\rightarrow e^{i{\vec {k}}\cdot {\vec {r}}}\qquad ({\vec {k}}\cdot {\vec {r}}\rightarrow \pm \infty )\,,}

que corresponden ak{\displaystyle {\vec {k}}}Estados asintóticos de onda plana orientados antes o después de su aproximación a la fuente de campo en el origen, respectivamente. Las funcionesψk(±){\displaystyle \psi _{\vec {k}}^{(\pm )}}están relacionados entre sí por la fórmula

ψk(+)=ψk().{\displaystyle \psi _{\vec {k}}^{(+)}=\psi _{-{\vec {k}}}^{(-)*}\,.}

Expansión de ondas parciales

La función de ondaψk(r){\displaystyle \psi _{\vec {k}}({\vec {r}})}se pueden expandir en ondas parciales (es decir, con respecto a la base angular) para obtener funciones radiales independientes del ángulo.w(η,ρ){\displaystyle w_{\ell}(\eta,\rho)}. Aquíρ=kr{\displaystyle \rho =kr}.

ψk(r)=4πr=0metro=iw(η,ρ)Ymetro(r^)Ymetro(k^).{\displaystyle \psi _{\vec {k}}({\vec {r}})={\frac {4\pi }{r}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }i^{\ell }w_{\ell }(\eta ,\rho )Y_{\ell }^{m}({\hat {r}})Y_{\ell }^{m\ast }({\hat {k}})\,.}

Un único término de la expansión puede aislarse mediante el producto escalar con un armónico esférico específico.

ψkmetro(r)=ψk(r)Ymetro(k^)dk^=Rk(r)Ymetro(r^),Rk(r)=4πiw(η,ρ)/r.{\displaystyle \psi _{k\ell m}({\vec {r}})=\int \psi _{\vec {k}}({\vec {r}})Y_{\ell }^{m}({\hat {k}})d{\hat {k}}=R_{k\ell }(r)Y_{\ell }^{m}({\hat {r}}),\qquad R_{k\ell }(r)=4\pi i^{\ell }w_{\ell }(\eta ,\rho )/r.}

La ecuación para una sola onda parcialw(η,ρ){\displaystyle w_{\ell}(\eta,\rho)}se puede obtener reescribiendo el laplaciano en la ecuación de onda de Coulomb en coordenadas esféricas y proyectando la ecuación sobre un armónico esférico específico.Ymetro(r^){\displaystyle Y_{\ell }^{m}({\hat {r}})}

d2wdρ2+(12ηρ(+1)ρ2)w=0.{\displaystyle {\frac {d^{2}w_{\ell }}{d\rho ^{2}}}+\left(1-{\frac {2\eta }{\rho }}-{\frac {\ell (\ell +1)}{\rho ^{2}}}\right)w_{\ell }=0\,.}

Las soluciones también se denominan funciones de onda de Coulomb (parciales) o funciones de Coulomb esféricas.z=2iρ{\displaystyle z=-2i\rho }transforma la ecuación de onda de Coulomb en la ecuación de Whittaker , de modo que las funciones de onda de Coulomb se pueden expresar en términos de funciones de Whittaker con argumentos imaginarios.METROiη,+1/2(2iρ){\displaystyle M_{-i\eta,\ell +1/2}(-2i\rho)}yWiη,+1/2(2iρ){\displaystyle W_{-i\eta,\ell +1/2}(-2i\rho)}Esto último puede expresarse en términos de las funciones hipergeométricas confluentes .METRO{\displaystyle M}yU{\displaystyle U}. ParaZ{\displaystyle \ell \in \mathbb {Z} }, se definen las soluciones especiales [ 4 ]

H(±)(η,ρ)=2i(2)miπη/2mi±iσρ+1mi±iρU(+1±iη,2+2,2iρ),{\displaystyle H_{\ell }^{(\pm )}(\eta ,\rho )=\mp 2i(-2)^{\ell }e^{\pi \eta /2}e^{\pm i\sigma _{\ell }}\rho ^{\ell +1}e^{\pm i\rho }U(\ell +1\pm i\eta ,2\ell +2,\mp 2i\rho )\,,}

dónde

σ=argΓ(+1+iη){\displaystyle \sigma _{\ell }=\arg \Gamma (\ell +1+i\eta )}

Se denomina desplazamiento de fase de Coulomb. También se definen las funciones reales.

F(η,ρ)=12i(H(+)(η,ρ)H()(η,ρ)),{\displaystyle F_{\ell }(\eta ,\rho )={\frac {1}{2i}}\left(H_{\ell }^{(+)}(\eta ,\rho )-H_{\ell }^{(-)}(\eta ,\rho )\right)\,,}
Función de onda de Coulomb regular F representada gráficamente de 0 a 20 con interacciones repulsivas y atractivas en Mathematica 13.1.
Función de onda de Coulomb regular F representada gráficamente de 0 a 20 con interacciones repulsivas y atractivas en Mathematica 13.1.
GRAMO(η,ρ)=12(H(+)(η,ρ)+H()(η,ρ)).{\displaystyle G_{\ell }(\eta ,\rho )={\frac {1}{2}}\left(H_{\ell }^{(+)}(\eta ,\rho )+H_{\ell }^{(-)}(\eta ,\rho )\right)\,.}

En particular uno tiene

F(η,ρ)=2miπη/2|Γ(+1+iη)|(2+1)¡ρ+1miiρMETRO(+1+iη,2+2,2iρ).{\displaystyle F_{\ell }(\eta ,\rho )={\frac {2^{\ell }e^{-\pi \eta /2}|\Gamma (\ell +1+i\eta )|}{(2\ell +1)!}}\rho ^{\ell +1}e^{i\rho }M(\ell +1+i\eta ,2\ell +2,-2i\rho )\,.}

El comportamiento asintótico de las funciones de Coulomb esféricasH(±)(η,ρ){\displaystyle H_{\ell }^{(\pm )}(\eta ,\rho )},F(η,ρ){\displaystyle F_{\ell }(\eta ,\rho )}, yGRAMO(η,ρ){\displaystyle G_{\ell }(\eta ,\rho )}en libertadρ{\displaystyle \rho }es

H(±)(η,ρ)mi±iθ(ρ),{\displaystyle H_{\ell }^{(\pm )}(\eta ,\rho )\sim e^{\pm i\theta _{\ell }(\rho )}\,,}
F(η,ρ)pecadoθ(ρ),{\displaystyle F_{\ell }(\eta ,\rho )\sim \sin \theta _{\ell }(\rho )\,,}
GRAMO(η,ρ)porqueθ(ρ),{\displaystyle G_{\ell }(\eta ,\rho )\sim \cos \theta _{\ell }(\rho )\,,}

dónde

θ(ρ)=ρηregistro(2ρ)12π+σ.{\displaystyle \theta _{\ell }(\rho )=\rho -\eta \log(2\rho )-{\frac {1}{2}}\ell \pi +\sigma _{\ell }\,.}

Las solucionesH(±)(η,ρ){\displaystyle H_{\ell }^{(\pm )}(\eta ,\rho )}corresponden a ondas esféricas entrantes y salientes. Las solucionesF(η,ρ){\displaystyle F_{\ell }(\eta ,\rho )}yGRAMO(η,ρ){\displaystyle G_{\ell }(\eta ,\rho )}son reales y se denominan funciones de onda de Coulomb regulares e irregulares. En particular, se tiene la siguiente expansión de onda parcial para la función de onda.ψk(+)(r){\displaystyle \psi _{\vec {k}}^{(+)}({\vec {r}})}[ 5 ]

ψk(+)(r)=4πρ=0metro=imiiσF(η,ρ)Ymetro(r^)Ymetro(k^),{\displaystyle \psi _{\vec {k}}^{(+)}({\vec {r}})={\frac {4\pi }{\rho }}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }i^{\ell }e^{i\sigma _{\ell }}F_{\ell }(\eta ,\rho )Y_{\ell }^{m}({\hat {r}})Y_{\ell }^{m\ast }({\hat {k}})\,,}

En el límiteη0{\displaystyle \eta \to 0}funciones de onda de Coulomb regulares/irregularesF(η,ρ){\displaystyle F_{\ell }(\eta ,\rho )},GRAMO(η,ρ){\displaystyle G_{\ell }(\eta ,\rho )}son proporcionales a las funciones de Bessel esféricas y a las funciones de Coulomb esféricas.H(±)(η,ρ){\displaystyle H_{\ell }^{(\pm )}(\eta ,\rho )}son proporcionales a las funciones de Hankel esféricas

F(0,ρ)/ρ=j(ρ){\displaystyle F_{\ell }(0,\rho )/\rho =j_{\ell }(\rho )}
GRAMO(0,ρ)/ρ=y(ρ){\displaystyle G_{\ell }(0,\rho )/\rho =-y_{\ell }(\rho )}
H(+)(0,ρ)/ρ=ih(1)(ρ){\displaystyle H_{\ell }^{(+)}(0,\rho )/\rho =i\,h_{\ell }^{(1)}(\rho )}
H()(0,ρ)/ρ=ih(2)(ρ){\displaystyle H_{\ell }^{(-)}(0,\rho )/\rho =-i\,h_{\ell }^{(2)}(\rho )}

y se normalizan igual que las funciones de Bessel esféricas.

0jl(kr)jl(kr)r2dr=0F(±1a0k,kr)krF(±1a0k,kr)krr2dr=π2k2δ(kk){\displaystyle \int \limits _{0}^{\infty }j_{l}(k\,r)j_{l}(k'r)\,r^{2}dr=\int \limits _{0}^{\infty }{\frac {F_{\ell }\left(\pm {\frac {1}{a_{0}k}},k\,r\right)}{k\,r}}{\frac {F_{\ell }\left(\pm {\frac {1}{a_{0}k'}},k'r\right)}{k'r}}\,r^{2}dr={\frac {\pi }{2k^{2}}}\delta (k-k')}

y similar para los otros 3.

Propiedades de la función de Coulomb

Las partes radiales para un momento angular dado son ortonormales. Cuando se normalizan en la escala del número de onda ( escala k ), las funciones de onda radiales del continuo satisfacen [ 6 ] [ 7 ].

0Rk(r)Rk(r)r2dr=δ(kk){\displaystyle \int _{0}^{\infty }R_{k\ell }^{\ast }(r)R_{k'\ell }(r)r^{2}dr=\delta (k-k')}

Otras normalizaciones comunes de las funciones de onda continuas se encuentran en la escala de número de onda reducido (k/2π{\displaystyle k/2\pi }-escala),

0Rk(r)Rk(r)r2dr=2πδ(kk),{\displaystyle \int _{0}^{\infty }R_{k\ell }^{\ast }(r)R_{k'\ell }(r)r^{2}dr=2\pi \delta (k-k')\,,}

y en la escala energética

0Rmi(r)Rmi(r)r2dr=δ(mimi).{\displaystyle \int _{0}^{\infty }R_{E\ell }^{\ast }(r)R_{E'\ell }(r)r^{2}dr=\delta (E-E')\,.}

Las funciones de onda radial definidas en la sección anterior están normalizadas a

0Rk(r)Rk(r)r2dr=(2π)3k2δ(kk){\displaystyle \int _{0}^{\infty }R_{k\ell }^{\ast }(r)R_{k'\ell }(r)r^{2}dr={\frac {(2\pi )^{3}}{k^{2}}}\delta (k-k')}

como consecuencia de la normalización

ψk(r)ψk(r)d3r=(2π)3δ(kk).{\displaystyle \int \psi _{\vec {k}}^{\ast }({\vec {r}})\psi _{{\vec {k}}'}({\vec {r}})d^{3}r=(2\pi )^{3}\delta ({\vec {k}}-{\vec {k}}')\,.}

Las funciones de onda de Coulomb del continuo (o de dispersión) también son ortogonales a todos los estados ligados de Coulomb [ 8 ].

0Rk(r)Rnorte(r)r2dr=0{\displaystyle \int _{0}^{\infty }R_{k\ell }^{\ast }(r)R_{n\ell }(r)r^{2}dr=0}

debido a que son autoestados del mismo operador hermitiano (el hamiltoniano ) con diferentes autovalores.

Lecturas adicionales

  • Bateman, Harry (1953), Funciones trascendentales superiores (PDF) , vol.  1, McGraw-Hill, archivado del original (PDF) el 11 de agosto de 2011 , consultado el 30 de julio de 2011..
  • Jaeger, JC; Hulme, HR (1935), "La conversión interna de rayos γ con la producción de electrones y positrones", Actas de la Royal Society de Londres. Serie A, Ciencias Matemáticas y Físicas , 148 (865): 708–728 , Bibcode : 1935RSPSA.148..708J , doi : 10.1098/rspa.1935.0043 , ISSN 0080-4630 , JSTOR 96298  
  • Slater, Lucy Joan (1960), Funciones hipergeométricas confluentes , Cambridge University Press , MR 0107026 .

Referencias

  1. Hill, Robert N. (2006), Drake, Gordon (ed.), Handbook of atomic, molecular and optical physics , Springer New York, pp. 153– 155, doi : 10.1007/978-0-387-26308-3 , ISBN  978-0-387-20802-2
  2. Landau, LD; Lifshitz, EM (1977), Curso de física teórica III: Mecánica cuántica, Teoría no relativista (3.ª ed.), Pergamon Press, pág. 569  
  3. Messiah, Albert (1961), Mecánica cuántica , North Holland Publ. Co., pág. 485 
  4. Gaspard, David (2018), "Fórmulas de conexión entre funciones de onda de Coulomb", J. Math. Phys. , 59 (11): 112104, arXiv : 1804.10976 , doi : 10.1063/1.5054368
  5. Messiah, Albert (1961), Mecánica cuántica , North Holland Publ. Co., pág. 426 
  6. ^ Formánek, Jiří ( 2004), Introducción a la teoría cuántica I (en checo) (2ª ed.), Praga: Academia, págs. 128-130  
  7. Landau, LD; Lifshitz, EM (1977), Curso de física teórica III: Mecánica cuántica, Teoría no relativista (3.ª ed.), Pergamon Press, pág. 121  
  8. Landau, LD; Lifshitz, EM ( 1977), Curso de física teórica III: Mecánica cuántica, Teoría no relativista (3.ª ed.), Pergamon Press, págs. 668–669