Estaba tratando de resolver la función de onda de un electrón en un átomo de hidrógeno confinado a solo dos dimensiones (en tierra plana). En primer lugar, descubrí la energía potencial electrostática en Flatland. Esto resultó ser
La ecuación de Schrödinger en dos dimensiones es
La ecuación angular es
Puedo resolver y cuantificar fácilmente la ecuación angular, pero he estado atascado en la ecuación radial durante días. No sé cómo resolver la ecuación diferencial radial (no soy muy bueno para resolver ED de segundo orden). Tampoco puedo averiguar cómo cuantificar la energía total del electrón.
Preguntas:
¿Cómo resolver la ecuación de la función de onda radial y cómo se cuantifica la energía en dicho átomo?
Cuando busqué fuentes de ayuda en línea, descubrí que todos consideraban que la energía potencial era de la forma en lugar de . Según la ley de Gauss, el potencial electrostático será logarítmico en dos dimensiones. ¿Por qué consideran que el potencial es de la forma ?
Estoy de acuerdo en que esta ecuación probablemente no admita una solución en términos de funciones elementales. Sin embargo, investigando un poco encontré este artículo: Atabek et. Alabama. física Rev. A 1974 que cubre sus preguntas (usan un potencial 2D logarítmico y discuten el espectro de valores propios de energía).
Una cosa útil que hacen es sustituir eliminar el término proporcional a en la ecuación y generar una ecuación de la forma: que es más agradable para trabajar.
Con respecto a una solución numérica:
Lo que tiene aquí es un "problema de valor límite de 2 puntos (BVP)" (ya que especifica el valor de , o equivalente , a y quieres que desaparezca como ) que a menudo se tratan con los llamados "Códigos de disparo".
Si está familiarizado con métodos numéricos más simples (de 1 paso, explícitos) como RK4, una cosa simple que puede hacer es iterar sobre los valores de con una conjetura para y verifique el valor de la solución ( ) en algún valor final grande de . Dependiendo de sus unidades, incluso algo como parece suficiente. Luego puede ajustar manualmente su cuadrícula de búsqueda de es hasta este valor final, , es aproximadamente cero.
Codifiqué este esquema iterativo de RK4 solo para verificar la forma cualitativa de las funciones de onda radial (con ) informado en el artículo vinculado de 1974 y parecen estar de acuerdo (cf. Fig. 4 del artículo).
Una advertencia : no preste atención a los valores numéricos de o en la siguiente figura, las unidades probablemente estén un poco confundidas. Además, normalicé los picos de a 1 para la comparación entre 's.
Sin embargo, para la inspección cualitativa, aquí hay tres funciones propias de la ecuación, una para , y dos para . Estos se encontraron utilizando el método iterativo RK4.
Dudo que esta ecuación tenga solución, aunque sería prudente consultar un libro sobre funciones especiales, como Abramovitz & Stegun o Gradshtein & Ryzhik.
Su ecuación parece ser correcta para una llanura . Sin embargo, muchos problemas se relacionan con el mundo real en 3D, donde el movimiento se limita a dos o incluso a una dimensión. En este caso el potencial sigue siendo 3D: . El caso unidimensional es notable, ya que la energía de enlace diverge, lo que durante mucho tiempo planteó un problema para analizar excitones en nanotubos de carbono.
G. Smith
Gert
Manvendra Somvanshi