Русская Википедия:Векторные сферические гармоники
Векторными сферическими гармониками являются векторные функции, преобразующиеся при вращениях системы координат так же, как скалярные сферические функции с теми же индексами, или определенные линейные комбинации таких функций.
Определения
1. Векторные сферические гармоники - векторные функции <math> \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)</math> , являющиеся собственными функциями операторов <math> \hat{\mathbf{J}}^{2}, \hat{J}_{z}, \hat{\mathbf{L}}^{2}, \hat{\mathbf{S}}^{2} </math>, где <math>\hat{\mathbf{L}}</math> - оператор орбитального углового момента, <math>\hat{\mathbf{S}}</math> - оператор спинового момента для спина 1, <math>\hat{\mathbf{J}}=\hat{\mathbf{L}}+\hat{\mathbf{S}}</math> - оператор полного углового момента. [1]
- <math> \begin{array}{l}{\hat{\mathbf{J}}^{2} \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)={J}({J}+1) \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)} \\ {\hat{J}_{z} \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)=M \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)} \\ {\hat{\mathbf{L}}^{2} \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)=L(L+1) \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)} \\ {\hat{\mathbf{S}}^{2} \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)=2 \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)}\end{array}
</math>
2. Часто (см., например, Рассеяние Ми) векторными гармониками называют фундаментальный набор решений векторного уравнения Гельмгольца в сферических координатах. [2] [3]
В этом случае векторные сферические гармоники порождаются скалярными функциями, являющимися решением уравнения Гельмгольца с волновым вектором <math> {\bf k}</math>.
- <math>
\begin{array}{l}{\psi_{e m n}=\cos m \varphi P_{n}^{m}(\cos \vartheta) z_{n}({k} r)} \\ {\psi_{o m n}=\sin m \varphi P_{n}^{m}(\cos \vartheta) z_{n}({k} r)}\end{array}
</math>
где <math>P_{n}^{m}(\cos \theta)</math> - присоединенные полиномы Лежандра, а <math> z_{n}({k} r) </math> - любая из сферических функций Бесселя.
Векторные гармоники выражаются как
- <math>
\mathbf{L}_{^e_o m n}=\mathbf{\nabla} \psi_{^e_o m n}
</math> - продольные гармоники
- <math>
\mathbf{M}_{^e_o m n}=\nabla \times\left(\mathbf{r} \psi_{^e_o m n}\right)
</math> - магнитные гармоники
- <math>
\mathbf{N}_{^e_o m n}=\frac{\nabla \times \mathbf{M}_{^e_o m n}}{\mathbf{k}}
</math> - электрические гармоники
Здесь вводятся производящие функции с вещественной угловой частью, но по аналогии можно ввести и комплексные гармоники.
3. Также часто вводятся шаровые векторы[4][5][6][7] , которые являются линейными комбинациями функций <math> \mathbf{Y}_{J M}^{L}(\vartheta, \varphi)</math>, но не являются собственными функциями квадрата орбитального углового момента, но определенным образом ориентированы относительно единичного орта <math> \mathbf{r} </math>. [1]. Определения и обозначения векторов этого типа в литературе широко варьируются, здесь приводится один из вариантов.
- <math>
\mathbf{X}_{JM}(\theta, \phi)=\frac{1}{\sqrt{l(l+1)}} \mathbf{L} Y_{JM} (\theta, \phi)=\frac{1}{i\sqrt{l(l+1)}} \mathbf{r}\times \nabla Y_{JM} (\theta, \phi)={\mathbf{Y}_{J M}^{(0)}(\vartheta, \varphi)=\mathbf{Y}_{J M}^{J}(\vartheta, \varphi)} </math> - векторы магнитного типа.
- <math>
i \mathbf{Y}_{J M}^{(0)}(\vartheta, \varphi) \times \mathbf{r}= {\mathbf{Y}_{J M}^{(1)}(\vartheta, \varphi)=\sqrt{\frac{J+1}{2 J+1}} \mathbf{Y}_{J M}^{J-1}(\vartheta, \varphi)+\sqrt{\frac{J}{2 J+1}} \mathbf{Y}_{J M}^{J+1}(\vartheta, \varphi)} </math> - векторы электрического типа
- <math>
\mathbf{r} Y_{J M}(\vartheta, \varphi)= {\mathbf{Y}_{J M}^{(-1)}(\vartheta, \varphi)=\sqrt{\frac{J}{2 J+1}} \mathbf{Y}_{J M}^{J-1}(\vartheta, \varphi)-\sqrt{ \fracШаблон:J+1{2 J+1}} \mathbf{Y}_{J M}^{J+1}(\vartheta, \varphi)}
</math> - продольный шаровой вектор
Для векторов этого типа производящими являются скалярные сферические функции <math>Y_{JM}(\vartheta, \varphi) </math> без радиальной части.
-
Электрические гармоники <math>\mathbf{N}_{^e_om1}</math>. <math>\mathbf{N}_{e01}</math> изображена дважды
-
Электрические гармоники <math>\mathbf{N}_{^e_om2}</math>. <math>\mathbf{N}_{e02}</math> изображена дважды
-
Электрические гармоники <math>\mathbf{N}_{^e_om3}</math>. <math>\mathbf{N}_{e03}</math> изображена дважды
-
Магнитные гармоники <math>\mathbf{M}_{^e_om1}</math>. <math>\mathbf{M}_{e01}</math> изображена дважды
-
Магнитные гармоники <math>\mathbf{M}_{^e_om2}</math>. <math>\mathbf{M}_{e02}</math> изображена дважды
-
Магнитные гармоники <math>\mathbf{M}_{^e_om3}</math>. <math>\mathbf{M}_{e03}</math> изображена дважды
Ортогональность
Решения векторного уравнения Гельмгольца подчиняются следующим отношениям ортогональности[3]:
- <math>
{\int_{0}^{2 \pi} \int_{0}^{\pi} \mathbf{L}_{^e_omn}\cdot\mathbf{L}_{^e_omn} \sin \vartheta d \vartheta d \varphi}{=(1+\delta_{m,0}) \frac{2 \pi}{(2 n+1)^{2}} \frac{(n+m) !}{(n-m) !} k^{2}\left\{n\left[z_{n-1}(k r)\right]^{2}+(n+1)\left[z_{n+1}(k r)\right]^{2}\right\}} </math>
- <math>
{\int_{0}^{2 \pi} \int_{0}^{\pi} \mathbf{M}_{^e_omn} \cdot \mathbf{M}_{^e_omn} \sin \vartheta d \vartheta d \varphi} {=(1+\delta_{m,0}) \frac{2 \pi}{2 n+1} \frac{(n+m) !}{(n-m) !} n(n+1)\left[z_{n}(k r)\right]^2} </math>
- <math>
\int_{0}^{2 \pi} \int_{0}^{\pi} \mathbf{N}_{^e_omn} \cdot \mathbf{N}_{^e_omn} \sin \vartheta d \vartheta d \varphi}{=(1+\delta_{m,0}) \frac{2 \pi}{(2 n+1)^{2}} \frac{(n+m) !}{(n-m) !} n(n+1)\left\{(n+1)\left[z_{n-1}(k r)\right]^{2}+n\left[z_{n+1}(k r)\right]^{2}\right\}
</math>
- <math>
{\int_{0}^{\pi} \int_{0}^{2 \pi} \mathbf{L}_{^e_omn} \cdot \mathbf{N}_{^e_omn} \sin \vartheta d \vartheta d \varphi} {=(1+\delta_{m,0}) \frac{2 \pi}{(2 n+1)^{2}} \frac{(n+m) !}{(n-m) !} n(n+1) k\left\{\left[z_{n-1}(k r)\right]^{2}-\left[z_{n+1}(k r)\right]^{2}\right\}} </math>
Все остальные интегралы по углам между различными функциями или функциями с различными индексами равны нулю.
Явный вид
Введем обозначение <math>\rho=kr</math>. Явный вид магнитных и электрических гармоник имеет следующую форму:
- <math>
\begin{aligned}{\mathbf{M}_{e m n}(k, \mathbf{r})={\frac{-m}{\sin (\theta)} \sin (m \varphi) P_{n}^{m}(\cos (\theta))} z_{n}(\rho)\mathbf{e}_{\theta}-} \\ {{-\cos (m \varphi) \frac{d P_{n}^{m}(\cos (\theta))}{d \theta}} }z_{n}(\rho) \mathbf{e}_{\varphi}\end{aligned} </math>
- <math>
\begin{aligned}{\mathbf{M}_{o m n}(k, \mathbf{r})={\frac{m}{\sin (\theta)} \cos (m \varphi) P_{n}^{m}(\cos (\theta)) }}z_{n}(\rho) \mathbf{e}_{\theta}- \\ {-\sin (m \varphi) \frac{d P_{n}^{m}(\cos (\theta))}{d \theta} z_{n}(\rho) \mathbf{e}_{\varphi}}\end{aligned} </math>
- <math>
\begin{aligned}{\mathbf{N}_{e m n}(k, \mathbf{r})=\frac{z_{n}(\rho)}{\rho} \cos (m \varphi) n(n+1) P_{n}^{m}(\cos (\theta)) \mathbf{e}_{\mathbf{r}}+} \\ {{{+\cos (m \varphi) \frac{d P_{n}^{m}(\cos (\theta))}{d \theta}}}} \frac{1}{\rho} \frac{d}{d \rho}\left[\rho z_{n}(\rho)\right] \mathbf{e}_{\theta}- \\ {{-m \sin (m \varphi) \frac{P_{n}^{m}(\cos (\theta))}{\sin (\theta)}}} \frac{1}{\rho} \frac{d}{d \rho}\left[\rho z_{n}(\rho)\right] \mathbf{e}_{\varphi}\end{aligned} </math>
- <math>
\begin{aligned} \mathbf{N}_{o m n} &(k, \mathbf{r})=\frac{z_{n}(\rho)}{\rho} \sin (m \varphi) n(n+1) P_{n}^{m}(\cos (\theta)) \mathbf{e}_{\mathbf{r}}+\\ &+\sin (m \varphi) \frac{d P_{n}^{m}(\cos (\theta))}{d \theta} \frac{1}{\rho} \frac{d}{d \rho}\left[\rho z_{n}(\rho)\right] \mathbf{e}_{\theta}+\\ &+{m \cos (m \varphi) \frac{P_{n}^{m}(\cos (\theta))}{\sin (\theta)}} \frac{1}{\rho} \frac{d}{d \rho}\left[\rho z_{n}(\rho)\right] \mathbf{e}_{\varphi} \end{aligned} </math> Можно видеть, что у магнитных гармоник отсутствует радиальная компонента. Для электрических гармоник радиальная компонента убывает быстрее, чем угловые, поэтому на больших <math>\rho</math> ей можно пренебречь. Кроме того, для электрических и магнитных гармоник с совпадающими индексами, угловые компоненты совпадают с точностью до перестановки полярного и азимутального единичных векторов, то есть при больших <math>\rho</math> векторы электрических и магнитных гармоник равны по модулю и перпендикулярны друг другу.
Явный вид продольных гармоник:
- <math>
\begin{aligned}\mathbf{L}_{^e_o{m n}}&(k, \mathbf{r})=\frac{\partial}{\partial r} z_{n}(k r) P_{n}^{m}(\cos \theta){^{\cos }_{\sin }} {m \varphi} \mathbf{e}_r +\\ &\frac{1}{r} z_{n}(k r) \frac{\partial}{\partial \theta} P_{n}^{m}(\cos \theta) {^{\cos }_{\sin}} m \varphi \mathbf{e}_{\theta}\mp \\ & \mp \frac{m}{r \sin \theta} z_{n}(k r) P_{n}^{m}(\cos \theta) {^{\sin }_{\cos}} m \varphi \mathbf{e}_{\varphi}\end{aligned}
</math>
Повороты и инверсия системы координат
При поворотах векторные сферические гармоники преобразуются друг через друга так же, как соответствующие скалярные сферические функции, которые являются производящими для конкретного типа векторнывх гармоник. Например, если производящими функциями являются обычные сферические функции, то векторные гармоники будут тоже преобразовываться с помощью D-матриц Вигнера[1][8][9]
- <math>
\hat{D}(\alpha, \beta, \gamma) \mathbf{Y}_{JM}^{(s)}(\theta, \varphi)= \sum_{m' = -\ell}^\ell [D^{(\ell)}_{MM'}(\alpha, \beta, \gamma)]^* \mathbf{Y}_{JM'}^{(s)}(\theta, \varphi), </math> Поведение при поворотах не отличается для электрических, магнитных и продольных гармоник.
При инверсии электрические и продольные сферические гармоники ведут себя так же, как скалярные сферические функции, то есть
- <math>
\hat{I}\mathbf{N}_{JM}(\theta, \varphi)= (-1)^J \mathbf{N}_{JM}(\theta, \varphi), </math> а магнитные обладают противоположной четностью:
- <math>
\hat{I}\mathbf{M}_{JM}(\theta, \varphi)= (-1)^{J+1} \mathbf{M}_{JM}(\theta, \varphi), </math>
Разложение плоской волны и интегральные соотношения
В этом параграфе будут использованы следующие обозначения
- <math>
\begin{aligned} Y_{e m n} &=\cos m \varphi P_{n}^{m}(\cos \theta) \\ Y_{o m n} &=\sin m \varphi P_{n}^{m}(\cos \theta) \end{aligned} </math>
- <math>
\mathbf{X}_{^e_o m n}\left(\frac{\mathbf{k}}{k}\right)=\nabla \times\left(\mathbf{k} Y_{^o_e m n}\left(\frac{\mathbf{k}}{k}\right)\right) </math>
- <math>
\mathbf{Z}_{^o_e m n}\left(\frac{\mathbf{k}}{k}\right)=i \frac{\mathbf{k}}{k} \times \mathbf{X}_{^e_o m n}\left(\frac{\mathbf{k}}{k}\right) </math> В случае, когда когда вместо <math>z_n</math> сферические функции Бесселя, с помощью формулы разложения комплексной экспоненты по сферическим функциям, можно получить следующие интегральные соотношения: [10]
- <math>
\mathbf {N}_{pmn}(k, \mathbf r) =\frac{i^{-n}}{4\pi} \int \mathbf Z_{pmn}\left(\frac{\mathbf{k}}{k}\right) e^{i \mathbf k \mathbf r} d\Omega_k </math>
- <math>
\mathbf {M}_{pmn}(k, \mathbf r) =\frac{i^{-n}}{4\pi} \int \mathbf X_{pmn}\left(\frac{\mathbf{k}}{k}\right) e^{i \mathbf k \mathbf r} d\Omega_k </math>
В случае, когда вместо <math>z_n</math> сферические функции Ханкеля, нужно использовать другие формулы разложения. [11] [10] Для векторных сферических гармоник получатся следующие соотношения:
- <math>
\mathbf{M}_{p m n}^{(3)}(k, \mathbf{r})=\frac{i^{-n}}{2 \pi k} \iint_{-\infty}^{\infty} d k_{ \|} \frac{e^{i\left(k_{x} x+k_{y} y \pm k_{z} z\right)}}{k_{z}}\left[\mathbf{X}_{p m n}\left(\frac{\mathbf{k}}{k}\right)\right] </math>
- <math>
\mathbf{N}_{p m n}^{(3)}(k, \mathbf{r})=\frac{i^{-n}}{2 \pi k} \iint_{-\infty}^{\infty} d k_{ \|} \frac{e^{i\left(k_{x} x+k_{y} y \pm k_{z} z\right)}}{k_{z}}\left[\mathbf{Z}_{p m n}\left(\frac{\mathbf{k}}{k}\right)\right] </math> где <math> k_{z}=\sqrt{k^{2}-k_{x}^{2}-k_{y}^{2}} </math>, а верхний индекс <math> (3) </math> означает, что используются сферические функции Ханкеля.
Ссылки
- ↑ 1,0 1,1 1,2 Варшалович Д. А., Москалёв А. Н., Херсонский В. К. Квантовая теория углового момента. Шаблон:Wayback — Л.: Наука, 1975.
- ↑ Шаблон:Книга
- ↑ 3,0 3,1 Шаблон:Книга
- ↑ Шаблон:Книга
- ↑ R.G. Barrera, G.A. Estévez and J. Giraldo, Vector spherical harmonics and their application to magnetostatics, Eur. J. Phys. 6 287-294 (1985)
- ↑ Шаблон:Книга
- ↑ R. Alaee, C. Rockstuhl, I. Fernandez‐Corbaton, Exact Multipolar Decompositions with Applications in Nanophotonics, Advanced Optical Materials 2019, 7, 1800783.
- ↑ H. Zhang, Yi. Han, Addition theorem for the spherical vector wave functions and its application to the beam shape coeffcients. J. Opt. Soc. Am. B, 25(2):255-260, Feb 2008.
- ↑ S. Stein, Addition theorems for spherical wave functions, Quarterly of Applied Mathematics, 19(1):15-24, 1961.
- ↑ 10,0 10,1 Шаблон:Cite web
- ↑ Шаблон:Cite web
