關於角動量算符們
這是一個Mathematica程式碼備忘錄,量子力學基礎略過
l
=
5
;
Jz
=
DiagonalMatrix
[
Reverse
[
Table
[
x
,
{
x
,
-
l
,
l
}]]];
Data
=
Table
[
Sqrt
[
l
(
l
+
1
)
-
x
^
2
+
Abs
[
x
]],
{
x
,
-
l
,
l
}];
Mat
=
DiagonalMatrix
[
Delete
[
Data
,
FirstPosition
[
Data
,
Max
[
Data
]]]];
Jm
=
Reverse
[
Map
[
Append
[
#
,
0
]
&
,
Append
[
Reverse
[
Mat
],
Table
[
0
,
{
x
,
2
l
}]]]];
Jp
=
Transpose
[
Jm
];
Jx
=
(
Jp
+
Jm
)
/
2
;
Jy
=
(
Jp
-
Jm
)
/
(
2
I
);
Jp
//
MatrixForm
Jm
//
MatrixForm
Jz
//
MatrixForm
Jy
//
MatrixForm
Jx
//
MatrixForm
為角量子數,
為約化普朗克常數
時,
即電子自旋角動量算符
時,即各種教材上熟知的:
以此類推,
時:
時:
時:(在知乎Markdown公式編輯器的效能邊緣試探)
不同的角量子數取值可以導致不同的算符分量之間的特殊關係
唯一兩條與角量子數取值無關的分量關係為:
(對易關係)
(本徵方程)
(我覺得自己早應該先從矩陣的角度理解角動量算符)
還有上面那個龐大的矩陣主對角線兩側的數值分佈規律
在角量子數取值很大的情況下,散點圖即為:
l=100
(其實基本上就是個圓了)
挖坑存程式碼,方便你我他
2019年2月28日續更
如果像自旋中對於Pauli矩陣的定義,重新定義算符:
由此定義得到
,
,
對於任意一個三維空間向量
記此向量模長
定義一個新的矩陣
則該矩陣的本徵值數量為
,且均勻地分佈在此向量模長兩倍的線段上
例如,
時,
的本徵值只有
時,本徵值有
時,本徵值有
時,本徵值有
以此類推