Гармонические сферические функции
Сферические функции
представляют собой угловую часть семейства ортогональных решений
уравнения Лапласа
, записанную в
сферических координатах
. Они широко используются для изучения физических
явлений в пространственных областях, ограниченных сферическими
поверхностями и при решении физических задач, обладающих
сферической симметрией.
Сферические функции имеют большое значение в теории
дифференциальных уравнений в частных производных
и
теоретической физике
, в частности в задачах расчёта
электронных орбиталей
в атоме,
гравитационного поля
геоида
,
магнитного поля
планет и интенсивности
реликтового излучения
.
Определение
Вещественные сферические функции Y
lm
,
l
=0…4 (сверху вниз),
m
=0…4 (слева направо). Функции отрицательного порядка Y
l-m
повёрнуты вокруг оси
Z
на 90/
m
градусов относительно функций положительного порядка.
Сферические функции
являются собственными функциями
оператора Лапласа
в
сферической системе координат
(обозначение
Y
l
m
(
θ
,
φ
)
{\displaystyle Y_{l}^{m}(\theta ,\varphi )}
). Они образуют
ортонормированную систему
в пространстве функций на сфере
S
2
{\displaystyle \mathbb {S} ^{2}}
в трёхмерном пространстве:
⟨
Y
l
m
;
Y
l
m
⟩
=
∬
|
Y
l
m
|
2
sin
θ
d
θ
d
φ
=
1
{\displaystyle \langle Y_{l}^{m};Y_{l}^{m}\rangle =\iint |Y_{l}^{m}|^{2}\sin {\theta }\,d\theta \,d\varphi =1}
⟨
Y
l
m
;
Y
l
′
m
′
⟩
=
∫
0
2
π
∫
0
π
Y
l
′
m
′
∗
Y
l
m
sin
θ
d
θ
d
φ
=
δ
l
l
′
δ
m
m
′
{\displaystyle \langle Y_{l}^{m};Y_{l'}^{m'}\rangle =\int \limits _{0}^{2\pi }\int \limits _{0}^{\pi }Y_{l'}^{m'*}Y_{l}^{m}\sin {\theta }\,d\theta \,d\varphi =\delta _{ll'}\delta _{mm'}}
,
где
*
обозначает
комплексное сопряжение
,
δ
l
l
′
{\displaystyle \delta _{ll'}}
—
символ Кронекера
.
Сферические функции имеют вид
Y
l
m
=
1
2
π
e
i
m
φ
Θ
l
m
(
θ
)
{\displaystyle Y_{l}^{m}={\frac {1}{\sqrt {2\pi }}}e^{im\varphi }\Theta _{lm}(\theta )}
,
где функции
Θ
l
m
(
θ
)
{\displaystyle \Theta _{lm}(\theta )}
являются решениями уравнения
1
sin
θ
d
d
θ
(
sin
θ
d
Θ
l
m
d
θ
)
−
m
2
sin
2
θ
Θ
l
m
+
l
(
l
+
1
)
Θ
l
m
=
0
{\displaystyle {\frac {1}{\sin {\theta }}}{\frac {d}{d\theta }}\left(\sin {\theta }{\frac {d\Theta _{lm}}{d\theta }}\right)-{\frac {m^{2}}{\sin ^{2}{\theta }}}\Theta _{lm}+l(l+1)\Theta _{l}^{m}=0}
и имеют вид
Θ
l
m
=
2
l
+
1
2
(
l
−
m
)
!
(
l
+
m
)
!
P
l
m
(
cos
θ
)
{\displaystyle \Theta _{l}^{m}={\sqrt {{\frac {2l+1}{2}}{\frac {(l-m)!}{(l+m)!}}}}P_{l}^{m}(\cos \theta )}
Здесь
P
l
m
(
cos
θ
)
{\displaystyle P_{l}^{m}(\cos \theta )}
— присоединённые
многочлены Лежандра
, а
m
!
{\displaystyle m!}
—
факториал
.
Присоединенные многочлены Лежандра с отрицательным
m
{\displaystyle m}
здесь вводятся как
P
ℓ
−
m
(
x
)
=
(
−
1
)
m
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
x
)
{\displaystyle P_{\ell }^{-m}(x)=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\ell }^{m}(x)}
Решение уравнения Лапласа в сферических координатах есть так называемая
, получаемая умножением сферической функции на решение радиального уравнения.
Вещественная форма
Вещественные сферические функции до шестого порядка
Для сферических функций форма зависимости от угла
φ
{\displaystyle \varphi }
— комплексная экспонента. Используя
формулу Эйлера
, можно ввести вещественные сферические функции. Иногда их удобнее использовать в связи с тем, что они могут быть наглядно показаны на иллюстрациях, в отличие от комплексных. Однако значимое удобство комплексных функций (утрачиваемое при переходе к вещественным) состоит в независимости квадрата их модуля
|
Y
l
m
|
2
{\displaystyle |Y_{l}^{m}|^{2}}
от угла
φ
{\displaystyle \varphi }
.
Y
ℓ
m
=
{
i
2
(
Y
ℓ
m
−
(
−
1
)
m
Y
ℓ
−
m
)
m
<
0
Y
ℓ
0
m
=
0
1
2
(
Y
ℓ
−
m
+
(
−
1
)
m
Y
ℓ
m
)
m
>
0.
=
{
2
(
−
1
)
m
Im
[
Y
ℓ
|
m
|
]
m
<
0
Y
ℓ
0
m
=
0
2
(
−
1
)
m
Re
[
Y
ℓ
m
]
m
>
0.
{\displaystyle {\begin{aligned}Y_{\ell m}&={\begin{cases}\displaystyle {i \over {\sqrt {2}}}\left(Y_{\ell }^{m}-(-1)^{m}\,Y_{\ell }^{-m}\right)&{\text{ }}\ m<0\\\displaystyle Y_{\ell }^{0}&{\text{ }}\ m=0\\\displaystyle {1 \over {\sqrt {2}}}\left(Y_{\ell }^{-m}+(-1)^{m}\,Y_{\ell }^{m}\right)&{\text{ }}\ m>0.\end{cases}}\\&={\begin{cases}\displaystyle {\sqrt {2}}\,(-1)^{m}\,\operatorname {Im} [{Y_{\ell }^{|m|}}]&{\text{ }}\ m<0\\\displaystyle Y_{\ell }^{0}&{\text{ }}\ m=0\\\displaystyle {\sqrt {2}}\,(-1)^{m}\,\operatorname {Re} [{Y_{\ell }^{m}}]&{\text{ }}\ m>0.\end{cases}}\end{aligned}}}
Обратное преобразование:
Y
ℓ
m
=
{
1
2
(
Y
ℓ
|
m
|
−
i
Y
ℓ
,
−
|
m
|
)
m
<
0
Y
ℓ
0
m
=
0
(
−
1
)
m
2
(
Y
ℓ
|
m
|
+
i
Y
ℓ
,
−
|
m
|
)
m
>
0.
{\displaystyle Y_{\ell }^{m}={\begin{cases}\displaystyle {1 \over {\sqrt {2}}}\left(Y_{\ell |m|}-iY_{\ell ,-|m|}\right)&{\text{ }}\ m<0\\\displaystyle Y_{\ell 0}&{\text{ }}\ m=0\\\displaystyle {(-1)^{m} \over {\sqrt {2}}}\left(Y_{\ell |m|}+iY_{\ell ,-|m|}\right)&{\text{ }}\ m>0.\end{cases}}}
Иногда вещественные сферические функции называют зональными, тессеральными и секториальными
.
Функции с
m
> 0 зависят от угла как косинус, а с
m
< 0 — как синус.
Y
ℓ
m
=
{
(
−
1
)
m
2
2
ℓ
+
1
4
π
(
ℓ
−
|
m
|
)
!
(
ℓ
+
|
m
|
)
!
P
ℓ
|
m
|
(
cos
θ
)
sin
(
|
m
|
φ
)
m
<
0
2
ℓ
+
1
4
π
P
ℓ
m
(
cos
θ
)
m
=
0
(
−
1
)
m
2
2
ℓ
+
1
4
π
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
cos
(
m
φ
)
m
>
0
.
{\displaystyle Y_{\ell m}={\begin{cases}\displaystyle (-1)^{m}{\sqrt {2}}{\sqrt {{2\ell +1 \over 4\pi }{(\ell -|m|)! \over (\ell +|m|)!}}}\ P_{\ell }^{|m|}(\cos \theta )\ \sin(|m|\varphi )&{\mbox{ }}m<0\\\displaystyle {\sqrt {2\ell +1 \over 4\pi }}\ P_{\ell }^{m}(\cos \theta )&{\mbox{ }}m=0\\\displaystyle (-1)^{m}{\sqrt {2}}{\sqrt {{2\ell +1 \over 4\pi }{(\ell -m)! \over (\ell +m)!}}}\ P_{\ell }^{m}(\cos \theta )\ \cos(m\varphi )&{\mbox{ }}m>0\,.\end{cases}}}
Повороты
Поворот вещественной сферической функции с m=0 и l=3. Коэффициенты не равны D-матрицам Вигнера, поскольку показаны вещественные функции, но могут быть получены при переразложении по комплексным функциям
Рассмотрим поворот системы координат
R
{\displaystyle {\mathcal {R}}}
, на
Углы Эйлера
α
,
β
,
γ
,
{\displaystyle \alpha ,\beta ,\gamma ,}
который преобрaзует единичный вектор
r
{\displaystyle \mathbf {r} }
в вектор
r
′
{\displaystyle {\mathbf {r} }'}
.
При этом углы
θ
′
,
φ
′
{\displaystyle \theta ',\varphi '}
вектора
r
′
{\displaystyle {\mathbf {r} }'}
в новой системе координат выражаются через углы в старой системе координат следующим образом
cos
θ
′
=
cos
θ
cos
β
+
sin
θ
sin
β
cos
(
φ
−
α
)
{\displaystyle \cos \theta ^{\prime }=\cos \theta \cos \beta +\sin \theta \sin \beta \cos(\varphi -\alpha )}
ctg
(
φ
′
+
γ
)
=
ctg
(
φ
−
α
)
cos
β
−
ctg
θ
sin
β
sin
(
φ
−
α
)
{\displaystyle \operatorname {ctg} \left(\varphi ^{\prime }+\gamma \right)=\operatorname {ctg} (\varphi -\alpha )\cos \beta -{\frac {\operatorname {ctg} \theta \sin \beta }{\sin(\varphi -\alpha )}}}
В новой системе координат сферическая функция с индексами
ℓ
{\displaystyle \ell }
и
m
{\displaystyle m}
будет представима в виде линейной комбинации всех функций с тем же номером
ℓ
{\displaystyle \ell }
и различными
m
{\displaystyle m}
. Коэффициентами в линейной комбинации являются комплексно- сопряженные
D-матрицы Вигнера
D
^
(
α
,
β
,
γ
)
Y
l
m
(
θ
,
φ
)
=
Y
ℓ
m
(
θ
′
,
φ
′
)
=
∑
m
′
=
−
ℓ
ℓ
[
D
m
m
′
(
ℓ
)
(
α
,
β
,
γ
)
]
∗
Y
ℓ
m
′
(
θ
,
φ
)
,
{\displaystyle {\hat {D}}(\alpha ,\beta ,\gamma )Y_{l}^{m}(\theta ,\varphi )=Y_{\ell }^{m}(\theta ',\varphi ')=\sum _{m'=-\ell }^{\ell }[D_{mm'}^{(\ell )}(\alpha ,\beta ,\gamma )]^{*}Y_{\ell }^{m'}(\theta ,\varphi ),}
Сферические функции с номером
ℓ
{\displaystyle \ell }
образуют базис неприводимого представления размерности
(
2
ℓ
+
1
)
{\displaystyle (2\ell +1)}
группы вращений SO(3).
Разложение плоской волны по сферическим функциям
Комплексная экспонента может быть представлена в виде разложения по сферическим функциям
e
i
k
⋅
r
=
4
π
∑
l
=
0
∞
i
l
j
l
(
k
r
)
∑
m
=
−
l
l
Y
l
m
∗
(
r
|
r
|
)
Y
l
m
(
k
|
k
|
)
{\displaystyle e^{i\mathbf {k} \cdot \mathbf {r} }=4\pi \sum _{l=0}^{\infty }i^{l}j_{l}(kr)\sum _{m=-l}^{l}Y_{l}^{m*}\left({\frac {\mathbf {r} }{|r|}}\right)Y_{l}^{m}\left({\frac {\mathbf {k} }{|k|}}\right)}
Здесь
j
n
(
x
)
=
π
2
x
J
n
+
1
2
(
x
)
{\displaystyle j_{n}(x)={\sqrt {\frac {\pi }{2x}}}J_{n+{\frac {1}{2}}}(x)}
— сферическая
функция Бесселя
Разложение произведений сферических функций
Разложения
Клебша-Гордана
для произведений двух сферических функций выглядят следующим образом
:
Y
ℓ
1
m
1
(
Ω
)
Y
ℓ
2
m
2
(
Ω
)
=
∑
L
,
M
(
2
ℓ
1
+
1
)
(
2
ℓ
2
+
1
)
4
π
(
2
L
+
1
)
⟨
ℓ
1
0
ℓ
2
0
|
L
0
⟩
⟨
ℓ
1
m
1
ℓ
2
m
2
|
L
M
⟩
Y
L
M
(
Ω
)
{\displaystyle Y_{\ell _{1}}^{m_{1}}(\Omega )Y_{\ell _{2}}^{m_{2}}(\Omega )=\sum _{L,M}{\sqrt {\frac {(2\ell _{1}+1)(2\ell _{2}+1)}{4\pi (2L+1)}}}\langle \ell _{1}\,0\,\ell _{2}\,0|L\,0\rangle \langle \ell _{1}\,m_{1}\,\ell _{2}\,m_{2}|L\,M\rangle Y_{L}^{M}(\Omega )}
См. также
Примечания
Литература
Приложения
Ссылки
Spherical harmonics applied to Acoustic Field analysis on
by
Stephen Wolfram
and
by Michael Trott,
Ссылки на внешние ресурсы
Словари и энциклопедии