Ни физически, ни математически задача не поставлена, по крайней мере, до конца.
У тебя же решением является чисто сферическая расходящаяся волна. Я вижу сферические волны, выглядывающие из-за камышей, но камыши-то - по ту сторону реки. Никогда это решение не будет общим решением, а от того, что будет за добавка, может очень многое зависеть. Вообще, где гран. условия? В одномерном случае понятно как они изменятся по отношению к стандартным, а в трехмерном? Задумалась и над физическим смыслом такого потенциала. 1Д - это, например, идеальный бесконечный одномерный "кристалл". Впринципе, любой
очень сильный и
очень быстро убывающий потенциал можно аппроксимировать дельта-функцией, но так можно и на кладбище угодить. Причем "потенциал"
![$U=a \delta (\vec r)$ $U=a \delta (\vec r)$](https://dxdy-02.korotkov.co.uk/f/d/8/2/d82006cdc6bbf69cca170134c698e8aa82.png)
в смысле возмущения в первом порядке стоит отличать от одномерного дельта-потенциала и потенциала трехмерной дельта-сферы при одинаковых обозначениях. Любое же физически "плохое" решение можно убрать, сказав, что оно противоестественно. В любом случае в пространстве появляется выделенная точка
![$r=0$ $r=0$](https://dxdy-02.korotkov.co.uk/f/1/b/0/1b0129678603dff06573a74c88866a5a82.png)
, сферическая симметрия налицо. Уравнение Шредингера через оператор углового момента будет:
![$\left[\frac{1}{r}\frac{d^2}{dr^2}r + \frac{2m}{{\hbar}^2}\left(E_{n_r l}-\frac{{\hbar}^2l(l+1)}{2mr^2}-U(r)\right)\right]R_{n_r l}(r)=0$ $\left[\frac{1}{r}\frac{d^2}{dr^2}r + \frac{2m}{{\hbar}^2}\left(E_{n_r l}-\frac{{\hbar}^2l(l+1)}{2mr^2}-U(r)\right)\right]R_{n_r l}(r)=0$](https://dxdy-04.korotkov.co.uk/f/f/8/5/f85dfbaa0745fef0b99f9f80317559ff82.png)
. То, что тебя интересует, - трехмерный аналог одномерного дельта-потенциала, - называется
потенциалом нулевого радиуса, который действует лишь на частицу с моментом
![$l=0$ $l=0$](https://dxdy-02.korotkov.co.uk/f/d/8/d/d8d907e87692bcdea2713804f65a134082.png)
и свойства состояний частицы в котором слабо зависят от конкретного вида
![$U(r)$ $U(r)$](https://dxdy-04.korotkov.co.uk/f/3/1/5/315ed515a7438fe959a31cef786171d282.png)
. На волновую функцию налагается граничное условие вида
![$\frac{(r\psi(r))'}{r\psi(r)}\to - a$ $\frac{(r\psi(r))'}{r\psi(r)}\to - a$](https://dxdy-04.korotkov.co.uk/f/7/f/0/7f0951e69f8f479b095e16b8a38f3faa82.png)
при
![$r\to 0$ $r\to 0$](https://dxdy-03.korotkov.co.uk/f/e/b/0/eb0ee41bf821c1d0f89bf72eba8cb1fd82.png)
, т.о.
![$\psi=\left(-\frac{1}{ar}+1+...)$ $\psi=\left(-\frac{1}{ar}+1+...)$](https://dxdy-02.korotkov.co.uk/f/d/f/f/dff85903eb15e0ff55df7eac4088d27a82.png)
. Для связанных состояний
![$E<0$ $E<0$](https://dxdy-04.korotkov.co.uk/f/7/0/a/70acff151a1c7329f30c0b7d59cf8df282.png)
при
![$r\ne0$ $r\ne0$](https://dxdy-01.korotkov.co.uk/f/0/8/3/083ff7c4f5e2ac6f266196ffca7554d782.png)
решение
![$\psi(r)=\frac{Ae^{-kr}}{\sqrt{4\pi}r}$ $\psi(r)=\frac{Ae^{-kr}}{\sqrt{4\pi}r}$](https://dxdy-04.korotkov.co.uk/f/3/6/d/36d87e80ddbe3e56729a25db33f5a33382.png)
, где
![$k=\sqrt{-\frac{2mE}{{\hbar}^2}}>0$ $k=\sqrt{-\frac{2mE}{{\hbar}^2}}>0$](https://dxdy-02.korotkov.co.uk/f/d/e/2/de2ec948915bb546b101d1bf5bd2fc1682.png)
,
![$l=0$ $l=0$](https://dxdy-02.korotkov.co.uk/f/d/8/d/d8d907e87692bcdea2713804f65a134082.png)
, нормировка
![$A=\sqrt{2a}$ $A=\sqrt{2a}$](https://dxdy-02.korotkov.co.uk/f/5/0/f/50fa9aaac515724672df2e7f4bd7c62e82.png)
. При
![$\psi=\frac{A}{\sqrt{4\pi}}\left(\frac{1}{r}-k+...\right)$ $\psi=\frac{A}{\sqrt{4\pi}}\left(\frac{1}{r}-k+...\right)$](https://dxdy-04.korotkov.co.uk/f/7/4/c/74c154f1d630c79508442008bf4f525c82.png)
, поэтому
![$k=a$ $k=a$](https://dxdy-02.korotkov.co.uk/f/9/2/e/92e721e293592a72d57d82f5500d3e1282.png)
. Cвязанные состояния будут только в случае
![$a>0$ $a>0$](https://dxdy-04.korotkov.co.uk/f/f/b/6/fb619bc416586dec067ec756bf572e6982.png)
, причем, как и в одномерном пространстве, будет только один уровень энергии
![$E=-\frac{{\hbar}^2a^2}{2m}$ $E=-\frac{{\hbar}^2a^2}{2m}$](https://dxdy-04.korotkov.co.uk/f/f/d/1/fd177b7e8c62022e1d80f7136e82b86782.png)
. Из
![$\triangle\psi(\vec r)=-4\pi\delta (\vec r)$ $\triangle\psi(\vec r)=-4\pi\delta (\vec r)$](https://dxdy-02.korotkov.co.uk/f/1/b/e/1becae9594327af6e91095fca23c3f8a82.png)
при
![$r\to0$ $r\to0$](https://dxdy-03.korotkov.co.uk/f/e/b/7/eb7e9f381e7ae52698050501f474066a82.png)
(
![$\triangle \left(\frac{1}{r}\right)=-4\pi\delta (\vec r)$ $\triangle \left(\frac{1}{r}\right)=-4\pi\delta (\vec r)$](https://dxdy-02.korotkov.co.uk/f/d/1/5/d159cf9d2b0f37111b3215582de9a77782.png)
)
![$\bar T=\infty$ $\bar T=\infty$](https://dxdy-01.korotkov.co.uk/f/0/a/0/0a02654f0fce1e664dc7c4681da7af7682.png)
и соответственно
![$\bar U=-\infty$ $\bar U=-\infty$](https://dxdy-01.korotkov.co.uk/f/0/f/2/0f2c552c1bd751ecabeba0b24f5c0ee482.png)
(
![$\bar T+\bar U=E$ $\bar T+\bar U=E$](https://dxdy-01.korotkov.co.uk/f/4/3/3/4338f4c6468ad105371cafb3ed19202c82.png)
). Очень много замечаний - и что потенциал нулевого радиуса всегда является потенциалом притяжения и что условие, которому должна удовлетворять волновая функция, является общим и получается из условия самосопряженного расширения эрмитова оператора. Интересно, читай монографии Островского, Демкова "Метод потенциалов нулевого радиуса в атомной физике" и Смирнова "Асимптотический методы в теории атомных столкновений". В "Обобщенных функциях" Гельфанда есть какие-то случаи решения уравнений, в которых обобщенные функции являются не свободными членами, а коэффициентами. Вообще-то уравнениями в частных производных с такими сложностями рассматривал еще Соболев, но мне, пока что, это не понадобилось, несмотря на то, что в ГП в эффективном гамильтониане контактное взаимодействие
![Razz :P](./images/smilies/icon_razz.gif)
.
---------------------------------------------------------------------------------------------------------------
Сегодня решила для потенциала
![$U=-\lambda\delta(r-a)$ $U=-\lambda\delta(r-a)$](https://dxdy-01.korotkov.co.uk/f/4/c/7/4c72787e3b7f9fc55678d0bb94f69b2a82.png)
для s-состояний (волновая функция не зависит от углов). А знаешь как? Ну не разделяются переменные, так перейдем от координатного представления к импульсному в уравнении Шредингера. Оператор кинетической энергии - оператор умножения, потенциальной - интегральный оператор с ядром. Ну интегрировать надо с дельта-функцией в сферических координатах, ну интеграл, считаемый с помощью теории вычетов, появляется - не смертельно все это. Потом напишу, если надо. Думаю, задачи эти - "стары как мир" - нам просто не попадаются.
----------------------------------------------------------------------------------------------------------------
Да, указанные монографии я не читала и даже не видела. Смотрела по обобщенным функциям литературу только. И...сидела и думала, как же решить такое простое, на первый взгляд, уравнение.