-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathsourcecode.tex
8294 lines (7506 loc) · 612 KB
/
sourcecode.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
% 导言区
% !TEX program = xelatex
\documentclass[a4paper,zihao=-4,linespread=1]{ctexrep}
\usepackage[]{mathrsfs}
\usepackage[]{animate}
\usepackage{color}
\usepackage{booktabs}
\usepackage{titlesec}
\usepackage{titletoc}
\usepackage[titletoc]{appendix}
\usepackage[normalem]{ulem} % 关于下划线的一个宏包
\usepackage[]{indentfirst} % 强制段首缩进
\usepackage[]{bm}
\usepackage[]{amsmath}
\allowdisplaybreaks[0] % 允许公式换页 0表示可以换页但是尽量不换页
\usepackage[]{esint}
\usepackage[]{amssymb}
\usepackage[]{extarrows}
\usepackage{graphicx}
\usepackage[]{subfigure}
\usepackage[colorlinks,linkcolor=blue,bookmarksopen=true,bookmarksnumbered=true,citecolor=blue]{hyperref} % 目录点击跳转
\usepackage[dvipsnames,svgnames]{xcolor}
\usepackage{tcolorbox}
\tcbuselibrary{breakable} % 分页
\tcbuselibrary{skins} % 箱子分页
\usepackage[]{tikz}
\usepackage[]{pgfplots}
\pgfplotsset{compat=1.18}
\usepackage{newtxtext}
\usepackage[strict]{changepage} % 提供一个 adjustwidth 环境
\usepackage{framed} % 实现方框效果
\usetikzlibrary{positioning} % tikz绘图相对位置
\usetikzlibrary{decorations.pathreplacing} % tikz画大括号
\usetikzlibrary{arrows.meta} % tikz画各种箭头
\usetikzlibrary{patterns} % 条纹
\usetikzlibrary {angles,quotes}
\tikzset{elegant/.style={smooth,thick,samples=50,cyan}} % 新建一种绘图样式
\usepackage[]{newtxmath}
\usepackage[]{bbm}
\usepackage[]{relsize}
\usepackage[]{braket}
\usepackage{siunitx}
\DeclareSIUnit\angstrom{\text {Å}}
\usepackage{diagbox}
\usepackage[final]{pdfpages}
\newcommand{\DOI}[1]{doi: \href{https://doi.org/#1}{#1}} % 引用文献超链接
\usepackage{mhchem} % 化学宏包
\usepackage{cancel} % 打叉叉的宏包
\usepackage{ntheorem}
\usepackage{arydshln} % 分块矩阵横竖线绘制
\usepackage{epigraph} % 章节引用名人名言
\usepackage{floatrow}
\usetikzlibrary{avoidpath}% Github上面的一个绘制复变函数积分围道的tikz宏包
\usepackage{xpinyin} % 汉语拼音注音
\usepackage{ytableau} % 绘制杨图和杨表
\usepackage{autobreak} %公式自动断行
\floatsetup{heightadjust=all, floatrowsep=columnsep}
\newfloatcommand{figurebox}{figure}[\nocapbeside][\dimexpr(\textwidth-\columnsep)/2\relax]
\newfloatcommand{tablebox}{table}[\nocapbeside][\dimexpr(\textwidth-\columnsep)/2\relax] % 三行实现图表并排
% 证明环境
\newtheorem*{proof}{证明}
% 证明结束标志符号
\newcommand\qed[0]{\hfill$\boxed{\text{搞定}}$}
% 定义公式环境
\newenvironment{lequation}{\large\begin{equation}}{\end{equation}}
\newenvironment{sequation}{\small\begin{equation}}{\end{equation}}
\newenvironment{tequation}{\tiny\begin{equation}}{\end{equation}}
% 向下的长箭头
\newcommand\longdownarrow[1]{\mathrel{\scalebox{1}[#1]{$\downarrow$}}}
% 定义方框
% 陈述性语句
\newtcolorbox[]{proposition}[2][]{title = \textbf{#2}, colback=Salmon!20, colframe=Salmon!90!Black,breakable, enhanced jigsaw}
% 历史注记
\newtcolorbox[]{history}[2][]{title = \textbf{#2}, colback=Emerald!10, colframe=cyan!40!black,breakable, enhanced jigsaw}
% 定理
\newtcolorbox[]{theorem}[2][]{title=\textbf{#2}, colback=red!5,colframe=red!75!black,breakable, enhanced jigsaw}
% 定义
\newtcolorbox[]{define}[2][]{title=\textbf{#2},colback=SeaGreen!10!CornflowerBlue!10,colframe=RoyalPurple!55!Aquamarine!100!,breakable, enhanced jigsaw}
% 例子
\newtcolorbox[]{example}[2][]{title=\textbf{例:#2}, colback=green!5,colframe=green!75!black,breakable, enhanced jigsaw}
% 思考注记
\definecolor{formalshade}{rgb}{0.99,0.97,0.93} % 文本框颜色
\newenvironment{thinknote}{%
\def\FrameCommand{%
\hspace{1pt}%
{\color{BurlyWood}\vrule width 2pt}%
{\color{formalshade}\vrule width 4pt}%
\colorbox{formalshade}%
}%
\MakeFramed{\advance\hsize-\width\FrameRestore}%
\noindent\hspace{-4.55pt}% disable indenting first paragraph
\begin{adjustwidth}{}{7pt}%
\vspace{2pt}\vspace{2pt}%
}
{%
\vspace{2pt}\end{adjustwidth}\endMakeFramed%
}
\setlength{\parindent}{4em}
\title{Notes of Griffiths' Introduction to QM} % 添加标题
\author{郑卜凡} % 添加作者
\date{\today} %最后更新日期
\begin{document}
% \maketitle %制作封面
\includepdf[height=29.7cm,width=21.0cm]{titlepage/fourthtitlepage.pdf} % 插入封面 height参数将封面拉伸为a4尺寸
\tableofcontents %制作目录
\setcounter{page}{0}
\thispagestyle{empty}
\chapter{波函数}
\section{薛定谔方程}
\begin{lequation}\label{S-1-D}
\boxed{
i\hbar\frac{\partial \Psi\left(x,t\right)}{\partial t}=-\frac{\hbar^{2}}{2m}\frac{\partial^{2}\Psi\left(x,t\right)}{\partial x^{2}}+
V\left(x,t\right)\Psi\left(x,t\right)
}
\end{lequation}
\begin{lequation}\label{S-3-D}
\boxed{
i\hbar\frac{\partial \Psi\left(\bm{r},t\right)}{\partial t}=\left[-\frac{\hbar^{2}}{2m}\nabla^{2}+
V\left(\bm{r},t\right)\right]\Psi\left(\bm{r},t\right)
}
\end{lequation}
上面给出的公式中第一个是一维形式, 第二个是三维一般形式。对于某些公式推导上, 使用薛定谔方程时, 常常是对方程两边进行共轭操作\footnote{物理实质可以理解
为时间反演对称性}(以一维形式为例)也即:
\begin{lequation}
-i\hbar\frac{\partial \Psi^{*}\left(x,t\right)}{\partial t}=-\frac{\hbar^{2}}{2m}\frac{\partial^{2}\Psi^{*}\left(x,t\right)}{\partial x^{2}}+
V\left(x,t\right)\Psi^{*}\left(x,t\right)
\end{lequation}
注意到上式的导出我们假定势能函数是实变函数, 这是有道理的, 但是书后面的习题\footnote{详见第三版Problem1.17}也给出了一个例子, 那就是在不稳定的系统中, 找到粒子的概率不是守恒的, 也就是说$P$依赖于时间。
这时如果引入含有虚部项的势能就可以很好地解释这一点。在量子力学中Schr{\"o}dinger方程的地位和牛顿第二定律一样, 现在只是描述粒子位置的函数变成了波函数。
Born后面给波函数一个统计上的解释, 这也就说明了在量子力学中的不确定性, 我们无法再像牛顿运动定律一样精确的预言一个粒子之后的运动, 我们只能给出它之后在某处的\uwave{概率}是多少
\begin{proposition}{波函数的统计诠释}
以一维情形为例, Born在统计上给出了对于波函数的解释, 他认为当一个微观粒子处于状态$\Psi\left(\bm{r},t\right)$时, 表示在$t$时刻在$x$处发现粒子的\uwave{概率}, 更准确的说\footnote{一个常用的代换是$\left|\Psi\left(x,t\right)\right|^{2}=
\Psi^{*}\left(\left(x,t\right)\right) \Psi\left(x,t\right) $}
\begin{lequation}
\int_{a}^{b}\left|\Psi\left(x,t\right)\right|^{2}dx=\left\{\text{$t$时刻在$\left[a,b\right]$内\\发现粒子的概率}\right\}
\end{lequation}
\end{proposition}
自然的, 我们会问, 测量时我们会发现粒子处于某点(C点), 那么测量之前粒子在哪?历史上有三种观点
\begin{history}{测量前粒子在哪?}
\textbf{1.现实主义学派:粒子还是在C点}, 这种观点完全否定了量子理论的不确定性, 也是爱因斯坦一直坚信的观点;\\
\textbf{2.正统学派:粒子哪也不在}, 这种观点认为正是我们的测量\uwave{迫使}粒子在C点, 这个观点被广泛接受, 但到底什么是测量还有待讨论;\\
\textbf{3.不可知论学派:拒绝回答}, 这种观点认为\uwave{测量前}本身就是难以定义的, 去讨论测量前粒子的位置也是没有意义的。
\end{history}
现代量子理论在实验上说明了正统学派的正确性\footnote{John Bell在1964年派排除了不可知论}, 有一点需要注意, 测量会导致波函数的坍塌(\ref{fig-1.1}), 坍塌成了一个类似于狄拉克delta函数的图像, 在波函数还没有按照薛定谔方程重新弥散开来的时候继续测量, 我们会发现测量结果不变, 也
就是说\uwave{测量完全改变了波函数}, 导致连续的测量得到的结果是一样的。
既然波函数在统计上可以解释为概率密度分布函数, 那么一定要满足\textbf{归一化条件}
\begin{proposition}{波函数归一化条件}
\begin{lequation}\label{normalized}
\int_{-\infty}^{\infty}\left|\Psi\left(x,t\right)\right|^{2}dx=1
\end{lequation}
\end{proposition}
而且从薛定谔方程的线性性可以看出, 如果$\Psi$是方程的解, 那么$A\Psi$也一定是方程的解, 这里的A就类似于微分方程通解里面的系数, 你需要使用归一化条件去确定它, 求解出来的$A$是不需要考虑相位问题的, 不会在物理上产生任何影响, 你只需要确定
它的模长就可以了, 下面一个关于归一化的定理让我们能更简单的对波函数进行归一化。
\begin{theorem}{$A$是一个与时间无关的常数}
\begin{lequation}
\label{normalized-independent-time}
\frac{d}{dt}\int_{-\infty}^{\infty}\left|\Psi\left(x,t\right)\right|^{2}dx=0
\end{lequation}
\end{theorem}
由于这个定理的正确性, 我们找到波函数的一个可能解后, 只需要任意代入一个$t$的值, 然后将波函数乘上一个常数因子$A$对波函数进行全空间积分解出$A$的大小即得到了波函数的真正有物理意义的解, \textbf{任何无法进行归一化的解(比如$\Psi=0$)都要舍去}。
\section{力学量的期望值和标准差}
\begin{define}{数学上的定义: 平均值和标准差}
\begin{center}
\begin{math}
\displaystyle
\left \langle x \right \rangle \overset{\text{def}}{=}\int_{-\infty}^{\infty}x\rho(x)dx \qquad\qquad\qquad
\sigma _{x}\overset{\text{def}}{=}\int_{-\infty }^{\infty } (x-\left \langle x \right \rangle )^2\rho(x)dx
\end{math}
\end{center}
其中$\rho(x)$是概率密度函数, 把上面的$x$换成$f(x)$就可以得到某个一般量的平均值和标准差
\end{define}
根据上面的定义我们可以得到一个更加常用的计算标准差的公式: $\sigma(x)=\sqrt{\left \langle x^2 \right \rangle-\left \langle x \right \rangle^2}$。
在进一步说明力学量的平均值的时候要先明确平均值的意义, 就比如说发现粒子所处位置的平均值, 你不能将其理解成连续测量一个系统很多次之后计算得到的平均值, 因为前面就说过
波函数会由于测量而坍缩, 连续多次对一个系统的测量得到的结果是一致的! 这里对平均值的定义是对于\uwave{系综}的, 也就是你需要对大量相同状态下的系统进行测量来求平均值, 或
者简单一点, 对一个系统测量很多次, 但每次测量要隔一段时间要等待波函数重新回到测量前未坍缩的样子。
\begin{proposition}{位置和动量的算符}
\begin{center}
\begin{math}
\displaystyle
\hat{x}=[x] \qquad\qquad\qquad\qquad \hat{p}=\left[-i\hbar\frac{\partial}{\partial x}\right]
\end{math}
\end{center}
值得一提的是, 上面的动量算符是利用$\left \langle p \right \rangle=m\frac{d\left \langle x \right \rangle}{dt}$得出来的
\footnote{这里的$\left \langle x \right \rangle$必须要事先写成关于$t$的函数, 另见Problem1.16(c)}
\end{proposition}
\begin{theorem}{任意一个力学量的统计量}
\begin{lequation}
\left \langle Q\left(x,p\right) \right \rangle=\int \Psi^{*}\left[Q\left(x,-i\hbar\frac{\partial}{\partial x}\right)\right]\Psi dx
\end{lequation}
\begin{equation}
\sigma_Q = \sqrt{\left \langle Q^2 \right \rangle-\left \langle Q \right \rangle^2}
\end{equation}
\end{theorem}
上面的定理说明了在量子力学中\uwave{算符}显得尤为重要, 你要计算一个力学量$Q$的平均值只要把这个力学量的算子\uwave{夹在}波函数中间再积分即可。确定一个力学量的算子时, 先将
这个力学量由经典力学的公式表示成关于位置$x$和动量$p$的函数然后将$x$和$p$全部换成对应的算符即可。\footnote{我们提倡使用算符这个新的工具去计算, 但有时候你会发现, 得到了$\left \langle x \right \rangle(t) $后直接使用
$\left \langle p \right \rangle=\frac{d\left \langle x \right \rangle}{dx}$更快。有时候直接使用定义(波函数是实质是测量到粒子在位置$x$的概率密度分布函数), 用$\int x\left| \Psi(x,t)\right|dx$直接计算也能达到事半功倍
的效果}
\section{海森堡不确定性原理}
这个原理在历史上曾经被称作\uwave{测不准原理}, 有很大的误导性, 实际上这个原理与测量误差毫无关系, 这是一个量子力学本身决定的原理。它表明了你对系统位置了解的越多, 比如说你把系统限制在某个
确定的轨道狭槽内, 那么你对系统动量的了解程度一定越低, 测出来的动量分布肯定越是分散。其中动量与波函数之间的关系最早由德布罗意(de Broglie)给出:
\begin{center}
\begin{math}
\displaystyle
\boxed{
p=\frac{h}{\lambda}=\frac{2\pi\hbar}{\lambda}
}
\end{math}
\end{center}
\begin{theorem}{Heisenberg Uncertainty Principle}
\begin{lequation}
\sigma_x\sigma_p \geq \frac{\hbar}{2}
\end{lequation}
\end{theorem}
\newpage
\begin{figure}[htbp]
\label{fig-1.1}
\centering
\begin{tikzpicture}[scale=0.4] %缩放,还可以设置xscale, yscale
\draw[->](4,-2)--(4,8) node[above]{$\Psi (x,0^-)$};
\draw[->](0,2)--(12,2) node[below]{$x$}; % 画坐标轴
\draw (4,2) node[below right]{$O$};
\draw[elegant,domain=0:12] plot(\x,{5.5*exp(-(\x-6)^2/4.56)+2});
\draw[dashed,red] (8,2)--(8,8);
\draw (8,2) node[below]{$C$};
\end{tikzpicture}
\hspace{4em}
\begin{tikzpicture}[scale=0.4] %缩放,还可以设置xscale, yscale
\draw[->](4,-2)--(4,8) node[above]{$\Psi (x,0^+)$};
\draw[->](0,2)--(12,2) node[below]{$x$}; % 画坐标轴
\draw (4,2) node[below right]{$O$};
\draw[dashed,red] (8,2)--(8,8);
\draw (8,2) node[below]{$C$};
\draw[->,red] (8,2)--(8,6);
\draw (8,6) node[right]{$\infty$};
\end{tikzpicture}
\caption{在$t=0$时刻测量, 波函数将在$C$处坍塌为狄拉克$\delta$函数}
\end{figure}
\begin{history}{关于sch\"{o}dinger方程}
薛定谔方程在现在的物理意义上来看应该认为是量子力学的一条基本假设, 因为这个公式并不是推出来的, 而是在受到德布罗意物质波的启发下\textbf{猜}出来的。
薛定谔的博士导师试图让它根据德布罗意的观点找出一个对应的波函数描述, 最开始想从相对论来构建方程, 后来走不通, 有的人说薛定谔是受到经典力学的波动方程
启发, 类比构造薛定谔方程的; 还有一种观点也比较有根据, 认为薛定谔是受到了它导师的导师玻尔兹曼的启发, 利用热力学第二定律中的玻尔兹曼熵$S=k\ln W$,一步步
构建出来的薛定谔方程, 具体做法就是将$k$换成$\hbar$, 然后根据能量量纲将$S$换为作用量$H$, 由于是波动方程, 所以还要在指数项上加一个虚数单位, 方程就
变成了这样:$$W=e^{\frac{iH}{\hbar}}$$然后经典力学里面刚好在正则变换中有一个Hamilton-Jacobi方程(这里的$S$不是熵)$$\frac{\partial S}{\partial t}+H=0$$
然后带入后将$W$换成更加\textbf{量子力学}的$\Psi$, 将$H$写成算符$\hat{H}$。
\setlength\parindent{2em}对的, 整个过程就是这么的没有道理, 所以我们说这个方程是猜出来的, 后来海森堡证明了波动力学与其创立的矩阵力学是等价的。薛定谔自己其实也没搞清楚他搞出来
的波函数有什么意义, 也是后来波恩给出了个统计解释, 不禁感叹新兴学科的发展总是迂回曲折!
\end{history}
\section{概率流密度}
用$\Psi^*$左乘\ref{S-3-D}得到:
\begin{equation}
i\hbar\Psi^*\frac{\partial}{\partial t}\Psi=-\frac{\hbar^2}{2m}\Psi^*\nabla^2\Psi+\Psi^*V\Psi
\end{equation}
然后再对上式取复共轭
\begin{equation}
-i\hbar\Psi\frac{\partial}{\partial t}\Psi^*=-\frac{\hbar^2}{2m}\Psi\nabla^2\Psi^*+\Psi V\Psi^*
\end{equation}
再将上面两式相减:
\begin{equation*}
\begin{aligned}
&i\hbar \left(\Psi^*\frac{\partial}{\partial t}\Psi+\Psi\frac{\partial}{\partial t}\Psi^*\right)=-\frac{\hbar^2}{2m}\left(\Psi^*\nabla^2\Psi-\Psi\nabla^2\Psi^*\right)\\
\Leftrightarrow&\frac{\partial}{\partial t}|\Psi|^2+\frac{i\hbar}{2m}\left(\Psi\nabla^2\Psi^*-\Psi^*\nabla^2\Psi\right)=0\\
\Leftrightarrow&\frac{\partial\rho}{\partial t}+\nabla\cdot\left[\frac{-i\hbar}{2m}\left(\Psi^*\nabla\Psi-\Psi\nabla\Psi^*\right)\right]=0
\end{aligned}
\end{equation*}
最终可以写成下面的连续性方程的形式:
\begin{equation}
\label{eq:1.12}
\boxed{
\frac{\partial\rho}{\partial t}+\nabla\cdot\mathbf{j}=0}
\end{equation}
其中$\rho\equiv |\Psi|^2$意义是概率密度分布
\begin{equation}
\boxed{
\mathbf{j}\equiv \frac{-i\hbar}{2m}\left(\Psi^*\nabla\Psi-\Psi\nabla\Psi^*\right)=\operatorname{Re}\left(\Psi^*\hat {\mathbf{v}}\Psi\right)
}
\end{equation}
上式中$\hat {\mathbf{v}}\equiv \hat {\mathbf{p}}/m$定义为速度算符。类似于流体力学中或是电动力学中的连续性方程,这里\ref{eq:1.12}也可以看成是某种“概率流”的连续性方程。$\mathbf{j}$就是对应的\textbf{概率流密度}。
根据Stokes公式我们任意选取一个区域$\Omega$,连续性方程可以写成下面的积分形式:
\[\frac{\partial }{\partial t}\int_{\Omega}\rho dV+\oiint_{\partial\Omega}\mathbf{j}\cdot d\mathbf{S}=0\]
第一项是在闭区域内发现粒子的概念,第二项可以看作是“概率流”流入这个区域的量。对于下一章我们要将的定态,$\partial \rho/\partial t=0$,那这就意味着概率不会“堆积”在某处,和其它领域内的连续性方程意义是一致的,只是现在不是那些具象的电荷或是水流了,而是抽象的概率。可以把$\mathbf{j}\cdot d\mathbf{\sigma}$解释为单位时间内粒子通过面元$ d\mathbf{\sigma}$的概率。
之后我们在讨论量子力学的散射问题时会用到这个概念。
\newpage
\chapter{定态Schr\"{o}dinger方程}
\section{定态和分离变量法}
列出Sch\"{o}dinger方程后的下一步就是解出系统对应的波函数, 首先我们应该从最最初等的分离变量法来解方程, 毫无疑问我们只能通过它来\uwave{猜出}解空间
的一小部分子集, 但是这种方法得出来的解却具有非常重要的物理意义, 后面的一系列讨论都默认势能函数不随时间变化(保守场)。
分离变量法的基本思想就是将$\Psi$拆分为两个函数, 一个只关于$x$, 一个只关于$t$, 即$\Psi(x,t)=\psi(x)\phi(t)$, 带入到Sch\"{o}dinger方程(\ref{S-1-D})我们可以得到
下面的式子\footnote{$\phi(t)$附带的常数项我们合并到$\psi(x)$里面了,反正最后是对$\Psi$进行归一化}。
\begin{lequation}
\label{wiggle-function}
\boxed{
\phi(t)=e^{-\frac{iEt}{\hbar}}
}
\end{lequation}
\begin{lequation}
\label{time-independent-equation}
\boxed{
-\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2}+V(x)\psi(x)=E\psi(x)
}
\end{lequation}
式\ref{wiggle-function}你可以称之为\uwave{wiggle-function}, 它只和时间相关而且你如果使用Euler公式\footnote{$e^{i\theta}=\cos\theta+i\sin\theta$}, 你
会发现它是按照正弦规律进行\uwave{振动}的, 暂且就先把它当成是波函数里面的振荡项。后面的式子(\ref{time-independent-equation})及其重要, 它就是我们要
谈的\textbf{定态薛定谔方程}, 不含时间, 其中的$E$实际上代表着系统的能量(哈密顿量\footnote{哈密顿量就是动能加势能, 与系统的拉格朗日量对应, 后者是动能减势能})。
\begin{define}{哈密顿算子}
\begin{center}
\begin{math}
\displaystyle
\hat{H}=\frac{\hat{p}^2}{2m}+[V]=-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}+V
\end{math}
\end{center}
\end{define}
使用哈密顿算子可以简化方程的书写为$\hat{H}\psi=E\psi$。下面列出来这些定态解重要的性质。
\begin{theorem}{定态解$\Psi(x,t)=\psi(x)e^{iEt/\hbar}$}
定态薛定谔方程的每个解都对应某个能量为$E$的\uwave{定态系统}:\\
1.每个定态解对应系统能量(Hamiltonian)无论何时测量都是$E$,很容易通过证明它的方差和平均值都不随时变来证明它, 其它力学量的平均值也不随时变;\\
2.定态解的时间项对概率密度没有任何贡献, i.e.用$\psi$替代$\Psi$去计算物理量的平均值不会产生任何差异, 因此有时候也刻意不区分两者;\\
3.由于定态解能量的选取是任意的\footnote{实际上不完全任意}, 所以存在无穷多个定态解, 一般的解就是这些由分离变量得出来的定态解的\uwave{线性组合}
\begin{lequation}
\label{find_c}
\boxed{
\Psi(x,t)=\sum_{n=1}^{\infty}c_n\psi_n e^{-iE_nt/\hbar}
}
\end{lequation}
其中常数$c$\footnote{可以是复数}需要使用初始波函数去确定。\\
特别注意, 现在由定态解合成的一般解\textbf{不一定}是定态的, 而且比较遗憾的是一般你不能把它写成一个闭公式, 只能用无穷级数表示。
\end{theorem}
由于波函数需要满足归一化条件(eq.\ref{normalized}), 所以所有不满足归一化条件的解我们都应舍去, 他们是没有物理意义的, 所以这一步可以帮我们排除很多定态解
的可能性。
\begin{thinknote}
这里实际上定态并不是我们真正要求的解, 我们要求的解是这些解的线性组合, 但是我们还是要将他们($\psi$)进行归一化处理, 但没有关系, 反正归一化只是乘上
一个常数, 但在这里我们进行归一化后可以使得解的物理意义更加明朗, 更重要的是它会简化我们后续的数学对解的性质上的讨论。叠加后, 定态解前面的系数也更加
具有物理意义。
\end{thinknote}
\begin{proposition}{定态解可归一化必要条件}
1.能量$E$的虚部为零, 是一个实数;\\
2.能量$E$要大于$V(x)$的最小值, 也就是说的能量不能一直小于势能\footnote{这个蛮好理解的, 因为动能项始终非负}。
\end{proposition}
\begin{proposition}{$c_n$的物理意义}
一般解前面对于定态解$\Psi_n(x,t)=\psi(x)e^{-iE_nt/\hbar}$的权重$\left|c_n\right|^2$表示\textbf{测量到系统能量值为$E_n$的概率}, 这也意味着, 你
对一个系统进行测量, 只可能得到它所包含的定态解分量对应的能量值, 一定是一组离散的值。
\end{proposition}
自然的, $c_n$也要满足归一化条件(eq.\ref{normalized-c})。注意, 我们上面提到是说在某次测量时发现系统能量为$E_n$, 而不是说系统在测量时处于能量为$E_n$的
一个定态, 系统的态始终是没有发生改变的。而且我们也注意到系统的能量并不是一成不变的, 不同的测量显现出来的能量也不同, 但是系统的能量的平均值是随时间或者测量改变
的(eq.\ref{energy-conservation}), 这一个性质就是在微观量子效应下的\textbf{能量守恒定律}。
\begin{lequation}
\label{normalized-c}
\sum_{n=1}^{\infty}\left|c_n\right|^2=1
\end{lequation}
\begin{lequation}
\label{energy-conservation}
\left \langle H \right \rangle=\sum_{n=1}^{\infty}\left|c_n\right|^2E_n
\Rightarrow \frac{\partial H}{\partial t}=0\footnote{利用这个性质, 你可以进一步简化平均能量的计算, 即计算$\left \langle H \right \rangle=\int \Psi(x,t)^*\hat{H}\Psi(x,t) dx$, 你可以直接代入$t=0$进一步
简化运算, 因为这个性质保证了平均能量与时间无关。}
\end{lequation}
\section{一维无限深方势阱}
\begin{define}{The Infinite Square Well}
粒子如果处于无限深势阱中, 那么它的势能可以写成下面的这种形式:
\begin{center}
\begin{math}
\displaystyle
V(x) = \begin{cases}
0, &0\leq x \leq a\\
\infty, & \text{otherwise}
\end{cases}
\end{math}
\end{center}
就像是在$x=0$和$x=a$处由两堵墙, 粒子被困在两堵墙之间, 然而在两堵墙之间它可以\uwave{自由移动}。
\end{define}
首先我们来看一下边界条件, 显然, 在势阱之外$\Psi(x,t)\equiv 0$, 同样也应该有$\Psi(0,t)=\Psi(a,t)=0$。根据定态薛定谔方程(eq.\ref{time-independent-equation}), 我
们要去求解的是一个分段函数形式给出的微分方程, 不妨分段去考虑它, $(0,a)$之外波函数恒等于$0$, 对于$(0,a)$那一部分, 势能为$0$, 定态薛定谔方程写成\footnote{注意到$E\leq0$时的解是平凡的或者不能归一化的}:
\begin{lequation}
\boxed{
-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2}=k^2\psi,\quad k=\frac{\sqrt{2mE}}{\hbar}
}
\end{lequation}
\begin{thinknote}
实际上这里的$k$从经典力学上看可以理解为\textbf{角波数}, $k=2\pi/\lambda\xlongequal[]{de-Broglie-eq.}p/\hbar$, 使用经典力学的观点来看,粒
子的能量为$E$, 那么动量为$\sqrt{2mE}$。
\end{thinknote}
这就是一个自由振子的数学模型, 代入边界条件\footnote{$\psi(0)=\psi(a)=0$, 而且这里悄然已经利用了$\psi(x)$的连续性, 后面会再次提到}并且归一化可以得到:
\begin{lequation}
\boxed{
\psi_n(x)=\sqrt{\frac{2}{a}}\sin\left(\frac{n \pi}{a}x\right),\quad k=\frac{n \pi}{a}
}
\end{lequation}
\begin{lequation}
\boxed{
E_n=\frac{n^2\pi^2\hbar^2 }{2ma^2}
}
\end{lequation}
其中\textbf{$n$是正整数}。计算的时候你会发现, 两个边界条件实际上不能完全确定微分方程解的两个常数项, 但是它却限制了$ka=n\pi(n=1,2,3,\ldots)$, 而$k$
的值和能量是相关的, 所以这也就说明了\textbf{量子力学系统中的定态解的能量这个时候只能取一系列离散数值}\footnote{对于束缚态是离散的, 对散射态不是}, 又考虑到系统的真实解是定态解的叠加, 而且每次测量系统的能量
时都会返回某个定态解的能量, 这也就体现了量子力学中\textbf{能量的不连续性}, 真正确定另一个常数的条件是归一化条件\footnote{回想一下对$\Psi$归一化等价于对$\psi$归一化}。
完整写下定态的波函数应该为:
\begin{lequation}
\Psi(x,t)=\sqrt{\frac{2}{a}}\sin\left(\frac{n \pi}{a}x\right)e^{in^2\pi^2\hbar t/2ma^2}
\end{lequation}
$n=1$时我们称作\textbf{基态}, 其它的我们称为\textbf{激发态}。观察不同解在$t=0$时的波函数(fg.\ref{fig-2.1})我们可以发现, 这些定态本质上就是一系列\uwave{驻波}, 势阱的
两个边界是它的两个波节。
\begin{figure}[htbp]
\centering
\subfigure[$n=1$]{
\begin{tikzpicture}[scale=0.6] %缩放,还可以设置xscale, yscale
\draw[->](0,-3)--(0,3) node[above]{$\Psi (x,0)$};
\draw[->](-2,0)--(8,0) node[below]{$x$}; % 画坐标轴
\draw (0,0) node[below right]{$O$};
\draw (6,0) node[below]{$a$};
\draw[elegant,domain=0:6] plot(\x,{3*sin(\x*pi/6 r)});
\draw[-,red] (-2,0)--(0,0);
\draw[-,red] (6,0)--(8,0);
\end{tikzpicture}
}
\subfigure[$n=2$]{
\begin{tikzpicture}[scale=0.6] %缩放,还可以设置xscale, yscale
\draw[->](0,-3)--(0,3) node[above]{$\Psi (x,0)$};
\draw[->](-2,0)--(8,0) node[below]{$x$}; % 画坐标轴
\draw (0,0) node[below right]{$O$};
\draw (6,0) node[below]{$a$};
\draw[elegant,domain=0:6] plot(\x,{3*sin(2*\x*pi/6 r)});
\draw[-,red] (-2,0)--(0,0);
\draw[-,red] (6,0)--(8,0);
\end{tikzpicture}
}
\subfigure[$n=3$]{
\begin{tikzpicture}[scale=0.6] %缩放,还可以设置xscale, yscale
\draw[->](0,-3)--(0,3) node[above]{$\Psi (x,0)$};
\draw[->](-2,0)--(8,0) node[below]{$x$}; % 画坐标轴
\draw (0,0) node[below right]{$O$};
\draw (6,0) node[below]{$a$};
\draw[elegant,domain=0:6] plot(\x,{3*sin(3*\x*pi/6 r)});
\draw[-,red] (-2,0)--(0,0);
\draw[-,red] (6,0)--(8,0);
\end{tikzpicture}
}
\caption{不同的$n$对应不同的驻波波形}
\label{fig-2.1}
\end{figure}
下面要谈到的性质虽然是根据无限深势阱总结来的, 但是他们却是普遍适用的\footnote{$\delta$符号的相关定义见附录A}。
\begin{theorem}{定态解的重要性质定理}
% 第一点和第二点可以删掉
1.不管势能本身是否具有对称性, 定态波函数总是关于势阱的中心成奇函数或偶函数(即关于中心轴对称或中心对称), 且随着$n$的变化交替出现;\\
2.不管势能本身形状如何, 波函数的波节(零点)个数总是随着能量的增大而增大, 且公差为$1$;\\
$\bigstar$3.\textbf{定态解之间相互正交}\footnote{回忆一下我们已经对$\psi$进行归一化}\footnote{积分为全空间, 对于无限深势阱积分限为$(0,a)$。
实际上这里使用术语\uwave{本征函数}相互正交更好, 因为我们描述的是$\psi$而非$\Psi$之间的关系。}
\begin{lequation}
\label{orthogonal}
\int \psi_m(x)^*\psi_n(x) dx=\delta_{mn}
\end{lequation}
$\bigstar$4.\textbf{定态解集合是完备的}:这个意思就是说你可以使用定态解的线性组合来表示\textbf{任何连续函数}, 对于这里的无限深势阱, 实质就是
\uwave{傅里叶级数}。
\end{theorem}
上面的性质第四点说明了我们始终可以找到一组合适的$c_n$去满足初始波函数, 而如何去求这些$c_n$又是基于定态解的正交性(eq.\ref{orthogonal}), 使用傅里叶方法我们可以
很容易的得到$c_n$\footnote{$\Psi(x,0)$是初始波函数, 如果初始波函数不是在$t=0$时刻给定的, 你可能要额外考虑一下wiggle-function项(eq.\ref{wiggle-function})(更改时间原点也是一个不错的选择, 复习一下波动学里面的操作)}:
\begin{lequation}
\boxed{
c_n=\int \psi_n(x)^*\Psi(x,0) dx
}
\end{lequation}
使用正交性你还可以去证明\footnote{分别基于初始波函数的归一化和定态薛定谔方程的算子表示法}eq.\ref{normalized-c}和eq.\ref{energy-conservation}。
\begin{thinknote}
上面的两个证明书上都有, 实际上你也可以证明按照上面的方法解出来的系统的波函数, $\Psi$已经自动归一化。但这是多此一举的, 因为无论是初始波函数, 还是确定$c_n$后的完整波函数, 波函数都是满足薛定谔方程的, 在
前面的章节我们就说明了如果$\Psi$在某一时刻是归一化的, 那么之后任一时刻它也是归一化的(eq.\ref{normalized-independent-time})
\end{thinknote}
在计算能量平均值的级数时, 经常会涉及黎曼函数, 这里不做深入展开, 仅仅列出几个常用的和式。
\begin{define}
{黎曼$\zeta$函数}
\begin{center}
\begin{math}
\displaystyle
\zeta(s)\overset{def}{=}\frac{1}{1^s}+\frac{1}{2^s}+\frac{1}{3^s}+\frac{1}{4^s}+\cdots
\end{math}
\end{center}
\end{define}
可以很容易的发现下面的等式成立:
\begin{lequation}
\boxed{
\frac{1}{1^s}+\frac{1}{3^s}+\frac{1}{5^s}+\frac{1}{7^s}+\cdots = (2^s-1)\zeta(s)
}
\end{lequation}
$s$为偶数时, 可以求出$\zeta$的精确值, 但对于奇数情况却异常复杂, 下面列出具体表达式和几个值供参考:
\begin{theorem}{$s$为偶数时的解}
\begin{center}
\begin{math}
\displaystyle
\zeta(2n)=\eta_n \pi^{2n},\quad
\eta_1=\frac{1}{ 6},\quad\eta_n=\sum_{k=1}^{n-1}(-1)^{n-1}\frac{\eta_{n-k}}{(2k+1)!} + (-1)^{n+1}\frac{n}{(2n+1)!}
\end{math}
\end{center}
\end{theorem}
\begin{table}[htbp]
\centering
\resizebox{\textwidth}{10mm}{
\begin{tabular}{cccccccc}
\hline
$s$ & $2$ & $4$ & $6$ & $8$ & $10$ & $12$ & $14$ \\
\hline
value & $\frac{\pi^2}{6}$ & $\frac{\pi^4}{90}$ & $\frac{\pi^6}{945}$ & $\frac{\pi^8}{9450}$ & $\frac{\pi^{10}}{93555}$ & $\frac{\pi^{12}}{638512875}$ & $\frac{2\pi^{14}}{18243225}$ \\
\hline
\end{tabular}}
\end{table}
\section{简谐振子}
通过上一节我们已经大致知道了如何去求解特定势能下的波函数, 大致来说就是找到所有的定态解, 每一个解对应一个常量(能量), 这些常量的选取是离散的, 然后我们再
对求出来的定态解进行叠加, 使用傅里叶方法定下权重即可。现在我们要碰到的势能函数形式是一个二次式, 但方程的求解却困难许多, 但这个工作是很有意义的, 因为任何
势能函数的驻点附近, 使用Taylor展开, 你都可以将它处理成一个简谐振子的模型。
\begin{define}{The Harmonic Oscillator}
势能形式是关于$x$的二次式, 我们使用类似离心势能的形式写出:
\begin{center}
\begin{math}
\displaystyle
V(x)=\frac{1}{2}m\omega^2x^2
\end{math}
\end{center}
\end{define}
定态薛定谔方程相应的写成$$-\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2}+\frac{1}{2}m\omega^2x^2\psi(x)=E\psi(x)$$
这个方程可以使用幂级数解法去解决, 我们先介绍一种比较物理的方法去求解, 即\uwave{升降阶算符法}。
\subsection{代数方法}
定态薛定谔方程的算子写法为$\hat{H}\psi=E\psi$, 代数方法的基本思路是去分解哈密顿算符:
$$\hat{H}=\frac{1}{2m}\left[\hat{p}^2+(m\omega x)^2\right]$$
如果是普通的复数$u^2+v^2$我们可以分解为$(u+iv)(u-iv)$, 这也就启发了我们像下面一样去定义升降阶算符, 前面多出来的$\frac{1}{\omega\hbar}$因子可以让
后面的形式更加美观。
\begin{define}{产生/湮灭算子}
\begin{lequation}
\hat{a}_\pm\overset{def}{=}\frac{1}{\sqrt{2m\omega\hbar}}(\mp i\hat{p}+m\omega x)
\end{lequation}
\end{define}
\begin{thinknote}
在计算算符时, 要格外小心, 算符一般情况下并不满足交换律, 你需要先将整个算符作用于一个函数上, 然后按照运算顺序逐个计算, 最后再进行化简。比如
$x$和$\hat{p}$就不满足交换律\footnote{量子力学很多与经典力学相冲突的地方就是因为这两个的不可交换性, 有的人将它当作公理去推导其它定理}, $
x\hat{p}f(x)=-ix\hbar\frac{df}{dx}$但是$\hat{p}xf(x)=-i\hbar\frac{d}{dx}\left[xf(x)\right]=-i\hbar[f(x)+x\frac{df(x)}{dx}]$。
\end{thinknote}
\begin{define}{对易子}
\begin{lequation}
\left[\hat{A},\hat{B}\right]\overset{def}{=}\hat{A}\hat{B}-\hat{B}\hat{A}
\end{lequation}
\end{define}
显然对易子是反对称的, 即$\left[\hat{A},\hat{B}\right]=-\left[\hat{B},\hat{A}\right]$, 计算可以得到下面很有用的关系式:
\begin{lequation}
\boxed{
\left[x,\hat{p}\right]=i\hbar\qquad
\left[\hat{a}_-,\hat{a}_+\right]=1
}
\end{lequation}
第一个式子也常称作\textbf{正则对易关系}, 我们继续使用算子重写定态薛定谔方程:
\begin{lequation}
\boxed{
\hat{H}=(\hat{a}_-\hat{a}_+-\frac{1}{2})\omega\hbar=(\hat{a}_+\hat{a}_-+\frac{1}{2})\omega\hbar
}
\end{lequation}
\begin{lequation}
\label{eq:2.17}
\boxed{
\omega\hbar(\hat{a}_\pm\hat{a}_\mp\pm\frac{1}{2})\psi(x)=E\psi(x)
}
\end{lequation}
\begin{theorem}{谐振子的定态解}
如果$\psi(x)$是谐振子在能量为$E$时的定态解, 那么$A\hat{a}_+\psi(x)$就是在能量为$E+\frac{1}{2}\omega\hbar$时的定态解\footnote{前面乘上
常数是归一化需要, 我们前面说过, 你求出来的解都需要进行归一化}, 同理, $A\hat{a}_-\psi(x)$是在能量为$E-\frac{1}{2}\omega\hbar$时的定态解。
根据这个性质我们便可知道谐振子是离散谱, 我们用$\ket{\psi_n}$标记这些定态。
\end{theorem}
直接根据湮灭和产生算子的定义可以很快地证明上述定理。现在很自然的就可以发现, 因为能量是不能一直递减的, 它至少要大于势能的最小值, 所以我们或许可以找到
一个具有最低能量的定态解, 不能再使用湮灭算子产生新的解, 这个解我们称为$\psi_0(x)$, 是基态解。构造这个解的思路就是它的下一级是没有物理意义, 不能归一化的。
也就是说有条件$$\hat{a}_-\psi_0(x)=0$$使用这个条件就可以得到基态解(不要忘了归一化)以及激发态解:
\begin{lequation}
\boxed{
\psi_0(x)=\left(\frac{m\omega}{\pi\hbar}\right)^{\frac{1}{4}}e^{-\frac{m \omega}{2\hbar}x^2}
}
\end{lequation}
\begin{lequation}
\label{eq:2.19}
\boxed{
\psi_n(x)=\frac{1}{\sqrt{n!}}\left(\hat{a}_+\right)^n\psi_0(x),E_n=\left(n+\frac{1}{2}\right)\omega \hbar
}
\end{lequation}
前面由归一化条件所决定的系数可以使用递推关系来确定, 再次强调, $\hat{a}_+\psi_n(x)=c_n\psi_{n+1}(x)$, $c$要根据归一化去确定, 并不是说使用产生算符可以
直接得到高一个能级的解, 前面还有一个归一化条件确定的待定系数。产生湮灭算符之间还有如下的非常有用的关系式, 也很容易证明。\footnote{在下一章讲了Dirac符号之后再来回过头来看这些}
\begin{theorem}{$\hat{a}_-$与$\hat{a}_+$是厄密共轭的(相互为伴随算子)}
\begin{lequation}
\int f^*(\hat{a}_\pm g)dx=\int (\hat{a}_\mp f)^* gdx
\end{lequation}
\end{theorem}
根据\ref{eq:2.17}, 以及\ref{eq:2.19}中给出的能级表达式我们可以得知:
\[\hat{a}_+\hat{a}_-\psi_n=n\psi_n,\quad\hat{a}_-\hat{a}_+\psi_n=(n+1)\psi_n\]
再通过直接计算$\left \| \hat{a}_+\psi_n(x) \right \|$和$\left \| \hat{a}_-\psi_n(x) \right \|$, 并利用上面的等式便
可以确定递推关系前面的系数(老规矩, 归一化只能确定模长, 但是我们只取最简单的那个值):
\begin{lequation}
\boxed{
\label{eq:2.21}
\hat{a}_+\psi_n=\sqrt[]{n+1}\psi_{n+1},\quad\hat{a}_-\psi_n=\sqrt[]{n}\psi_{n-1}
}
\end{lequation}
使用算符去表达可以使得计算更加简便, 你可以很容易的验证这些定态解之间相互正交, 所以可以使用傅里叶方法(你要求哪个常数, 你就在初态前面乘上对应的定态解, 然后
全空间内积分, $t\neq0$时要考虑一下wiggle-function项)去定系数。计算力学量平均值的时候也可以使用产生湮灭算符去重新描述$x$和$\hat{p}$简化运算:
\begin{lequation}
\boxed{
x =\sqrt[]{\frac{\hbar}{2m\omega}}\left(\hat{a}_++\hat{a}_-\right), \quad
p =i\,\sqrt[]{\frac{\hbar m \omega}{2}}\left(\hat{a}_+-\hat{a}_-\right)
}
\end{lequation}
根据上面的等式可以比较容易的求出下面两个用的比较多的公式:
\begin{align*}
\left\langle n|x| n^{\prime}\right\rangle & =\sqrt{\frac{\hbar}{2 m \omega}}\left\langle n\left|\left(a_{+}+a_{-}\right)\right| n^{\prime}\right\rangle=\sqrt{\frac{\hbar}{2 m \omega}}\left[\sqrt{n^{\prime}+1}\left\langle n \mid n^{\prime}+1\right\rangle+\sqrt{n^{\prime}}\left\langle n \mid n^{\prime}-1\right\rangle\right] \\
& =\sqrt{\frac{\hbar}{2 m \omega}}\left(\sqrt{n^{\prime}+1} \delta_{n, n^{\prime}+1}+\sqrt{n^{\prime}} \delta_{n, n^{\prime}-1}\right)=\boxed{\sqrt{\frac{\hbar}{2 m \omega}}\left(\sqrt{n} \delta_{n^{\prime}, n-1}+\sqrt{n^{\prime}} \delta_{n, n^{\prime}-1}\right) }\\
\left\langle n|p| n^{\prime}\right\rangle & =\boxed{i \sqrt{\frac{m \hbar \omega}{2}}\left(\sqrt{n} \delta_{n^{\prime}, n-1}-\sqrt{n^{\prime}} \delta_{n, n^{\prime}-1}\right)}
\end{align*}
下面\ref{fig:2.2}画出前三个能级对应波函数的能量, 不难看出图像具有的规律性和一维无限深势阱相同。
\begin{figure}[htbp]
\label{fig:2.2}
\centering
\subfigure[$n=1$]{
\includegraphics[scale=0.7]{fig/2-3-a.eps}
}
\subfigure[$n=2$]{
\includegraphics[scale=0.7]{fig/2-3-b.eps}
}
\subfigure[$n=3$]{
\includegraphics[scale=0.7]{fig/2-3-c.eps}
}
\caption{可以看到, 图像之间的递推联系还是满足前面的规则的。}
\end{figure}
\subsection{解析方法}
\begin{define}{两个无量纲数}
\begin{center}
\begin{math}
\displaystyle
\xi = \sqrt{\frac{m\omega}{\hbar}}x \qquad\qquad K=\frac{2E}{\omega\hbar}
\end{math}
\end{center}
\end{define}
重写方程为:
\begin{lequation}
\label{rewrite}
\frac{d^2\psi}{d\xi^2}=(\xi^2-K)\psi
\end{lequation}
\begin{thinknote}
你可能对这个形式的得到有所疑问, 实际上, 这就是一个求导的链式法则的问题, 注意, 换元前的方程$\frac{d\psi}{dx}$表示把$\psi$写成$f(x)$后再求导, 而
换元后的方程中$\frac{d\psi}{d\xi}$表示将$\psi$表示为$\phi(\xi)$后再进行求导。而$\xi=\varphi(x)$, 使用链式法则便可以理解这里的换元了。
\end{thinknote}
使用幂级数方法求解的第一步就是去求其\uwave{渐近解}, 也就是说去观察$\xi \to \pm\infty$的时候方程的行为。这个方程在$\xi\gg K$时
$$\frac{d^2\psi}{d\xi^2}\approx\xi^2\psi$$这个方程的通解形式是$$\psi\approx Ae^{-\xi^2/2}+Be^{\xi^2/2}$$显然$B\neq 0$时波函数不能进行归一化, 所
以, 方程\ref{rewrite}的解的形式应该是$\psi = h(\xi)e^{-\xi^2/2}$, 代入后有$$\frac{d^2h}{d\xi^2}-2\xi\frac{dh}{d\xi}+(K-1)h=0$$这是一个二阶方程
但是是变系数, 所以也很难求解, 只能考虑使用幂级数解法, 两边进行Taylor展开后解得$$h(\xi)=\sum_{j=0}^\infty a_j\xi^j,\qquad a_{j+2}=\frac{2j+1-K}{(j+1)(j+2)}a_j$$
我们只要确定$a_0$和$a_1$也就可以确定解, 二阶方程刚好两个待定系数, 但是我们只有波函数归一化这一个方程似乎无法去完整的确定两个待定系数。这里的原因是能量
的量子化取值, 导致了最后方程只会含有一个待定系数。
我们可以证明, $h(\xi)$只能在某些特定的情况下在$\xi \to \pm \infty$时收敛, 这也就决定了\textbf{$K$的取值是量子化的}。事实上, 收敛的充要条件是, 上面的数列$\{a_j\}$
会终止于$0$。
\begin{proposition}{$a_j$会终止于$0$}
\begin{center}
\begin{math}
\displaystyle
K=2n+1,\quad (n=0,1,2,\ldots)
\end{math}
\end{center}
如果$n$是奇数, 那么$a_{2n+1}$会终止于$0$, 但是$a_0=0$也就是说必须有$a_{2n}\equiv 0$, $n$为偶数时类似
\end{proposition}
这样方程最终得到的解就只会有一个待定系数了, 而且由递推公式得知是正比于$a_0$或$a_1$的, 再利用归一化条件便可以得到解
\begin{define}{厄米多项式}
把$h(\xi)$的$a_0$或者$a_1$因子去掉并且乘上一个数将最高次项前面的因子化为$2^n$, 得到的多项式为\textbf{厄米多项式}, 记作$H_n(\xi)$\\
例如:$h_2(\xi)=a_0(1-2\xi^2)=-\frac{1}{2}a_0(4\xi^2-2)$则$H_2(\xi)=4\xi^2-2$
\end{define}
我们不加证明地列出波函数归一化后的解与厄米多项式之间的关系\footnote{前面因子的相位按照惯例按最简单的取}:
\begin{lequation}
\boxed{
\psi_n(x)=\left(\frac{m\omega}{\pi\hbar}\right)^{\frac{1}{4}}\frac{1}{\sqrt{2^nn!}}H_n(\xi)e^{-\xi^2/2}
}
\end{lequation}
\begin{theorem}{$H_n(\xi)$的诸多性质}
\begin{itemize}
\item
\begin{math}
\displaystyle
H_n(\xi)=(-1)^ne^{\xi^2}\left(\frac{d}{d\xi}\right)^ne^{-\xi^2}
\end{math}
\item
\begin{math}
\displaystyle
H_{n+1}(\xi)=2\xi H_n(\xi)-2nH_{n-1}(\xi)
\end{math}
\item
\begin{math}
\displaystyle
\frac{d H_n}{d\xi}=2n H_{n-1}(\xi)
\end{math}
\item
\begin{math}
\displaystyle
e^{-z^2+2z\xi}=\sum_{n=0}^{\infty}\frac{z^n}{n!}H_n(\xi)\quad\text{(generate function)}
\end{math}
\end{itemize}
\end{theorem}
量子效应下, 谐振子的行为和经典力学非常不同, 不再有\uwave{振幅}这一概念, 粒子可以在无穷远处被发现, 不违背能量守恒定律正是因为量子力学中, 我们只谈系综的
力学量的平均值, 不再对于某个粒子有诸如动能这些的定义了, 我们只讲系综的平均效应, 只谈概率, 不谈确定性。
\subsection{*谐振子的相干态}
下面介绍一个量子光学中比较常用的概念——相干态。下面讨论中我们采用与量子场论教材接轨的记号,$\hat a_+\to A^\dagger,\hat a_-\to A$。
首先我们在占有(居位)数表象\footnote{也就是粒子数算符$N=A^\dagger A$的本征矢$\ket{n}(\psi_n)$构成的表象,后面讲二次量子化时会再次提到}下求湮灭算符的本征态,注意$A$不是厄米算符,所以其本征值可能是任意复数。
\begin{equation}
\begin{aligned}
A\ket{z}=z\ket{z}&\Rightarrow \Braket{n|A|z}=z\braket{n|z}\\
&\Rightarrow \sqrt{n+1}\braket{n+1|z}=z\braket{n|z}\\
&\Rightarrow \ket{z}=C\sum_{n=0}^{\infty}\frac{z^n}{\sqrt{n!}}\ket{n}
\end{aligned}
\end{equation}
根据归一化条件$\braket{z|z}=1$得到相干态的完整表达式:
\begin{equation}
\boxed{
\ket{z}=\mathrm{e}^{-\frac{1}{2}|z|^2}\sum_{n=0}^{\infty}\frac{z^n}{\sqrt{n!}}\ket{n}
}
\end{equation}
注意$\ket{n}$是定态,随时间演变不变(只差一个指数因子),但是叠加起来的$\ket{z}$没有这一性质:
\begin{equation}
\begin{aligned}
\ket{z(t)}&=\mathrm{e}^{-\frac{1}{2}|z|^2}\sum_{n=0}^{\infty}\frac{z^n}{\sqrt{n!}}\mathrm{e}^{\frac{-iE_n}{\hbar}t}\ket{n}\\
&=\mathrm{e}^{-\frac{1}{2}i\omega t}\mathrm{e}^{-\frac{1}{2}|z\mathrm{e}^{-i\omega t}|^2}\sum_{n=0}^{\infty}\frac{\left(z\mathrm{e}^{-i\omega t}\right)^n}{\sqrt{n!}}\ket{n}\\
&=\mathrm{e}^{-\frac{1}{2}i\omega t}\ket{z\mathrm{e}^{-i\omega t}}\sim\ket{z\mathrm{e}^{-i\omega t}}
\end{aligned}
\end{equation}相干态在演变时我们只要作替换$z\to z\ket{z\mathrm{e}^{-i\omega t}}$就好。
\subsubsection{相干态的性质}
我们计算下相干态下的$\Braket{X}$和$\Braket{H}$,我们只需要先计算$t=0$时刻的值,最后再作替换$z\to z\ket{z\mathrm{e}^{-i\omega t}}$就好。
\begin{equation}
\begin{aligned}
\Braket{z|X|z}&=\sqrt{\frac{\hbar}{2m\omega}}\Braket{z|A^\dagger+A|z}\\
&=\sqrt{\frac{\hbar}{2m\omega}}\mathrm{e}^{-|z|^2}\sum_{m=0}^\infty\sum_{n=0}^\infty\Braket{m|A^\dagger+A|n}\frac{\left(z^*\right)^mz^n}{\sqrt{m!n!}}\\
&=\sqrt{\frac{\hbar}{2m\omega}}\mathrm{e}^{-|z|^2}\sum_{m=0}^\infty\sum_{n=0}^\infty\left(\sqrt{n+1}\delta_{m,n+1}+\sqrt{n}\delta_{m,n-1}\right)\frac{\left(z^*\right)^mz^n}{\sqrt{m!n!}}\\
&=\sqrt{\frac{\hbar}{2m\omega}}(z+z^*)\\
\Rightarrow&\Braket{z(t)|X|z(t)}=\lambda\left(z\mathrm{e}^{-i\omega t}+z^*\mathrm{e}^{i\omega t}\right)
\end{aligned}
\end{equation}
\begin{equation}
\begin{aligned}
\Braket{z|H|z}&=\omega \hbar\Braket{z|A^\dagger A|z}+\frac{1}{2}\omega\hbar\\
&=\omega \hbar \mathrm{e}^{-|z|^2}\sum_{m=0}^\infty\sum_{n=0}^\infty\Braket{m|A^\dagger A|n}\frac{\left(z^*\right)^mz^n}{\sqrt{m!n!}}+\frac{1}{2}\omega\hbar\\
&= \omega\hbar |z|^2+\frac{1}{2}\omega\hbar
\end{aligned}
\end{equation}
这里$\frac{1}{2}\omega\hbar$是量子力学特有的,用量子场论的方法叫\textbf{真空零点能}。在量子力学里面谐振子能力最低也不为0,但是我们可以把这一项丢掉,因为我们最终能测量的都是能级差,用量子场论来说就是取正规编序后的哈密顿量$:H:=\omega\hbar A^\dagger A$为体系的哈密顿量。我们最终得到:
\begin{equation}
\Braket{z(t)|:H:|z(t)}=2m\omega^2\lambda^2|z|^2
\end{equation}
现在取$\lambda z=\frac{1}{2}x_0\mathrm{e}^{-i\varphi}$得到:
\begin{equation}
\braket{X}=\frac{1}{2}x_0\left(\mathrm{e}^{-i(\omega t+i\varphi )}+\mathrm{e}^{i(\omega t+i\varphi )}\right)=x_0\cos(\omega t+\varphi)
\end{equation}
\begin{equation}
\braket{H}=\frac{1}{2}m\omega^2 x_0^2
\end{equation}
不难发现在相干态下$\braket{X}$和$H$的演化和经典谐振子十分相似,这是相干态的重要性质,历史上也作为提出相干态概念的来源。
$\bullet$虽然相干态是归一的,但是\textbf{任意两个不同的态不正交}:
\begin{proof}
\begin{align*}
\braket{z|z^\prime}&=\mathrm{e}^{-\frac{1}{2}\left(|z|^2+|z^\prime|^2\right)}\sum_{m=0}^\infty\sum_{n=0}^\infty\braket{m|z}\frac{\left(z^*\right)^m\left(z^\prime\right)^n}{\sqrt{m!n!}}\\
&=\mathrm{e}^{-\frac{1}{2}\left(|z|^2+|z^\prime|^2\right)+z^*z^\prime}
\end{align*}
\qed
\end{proof}
$\bullet$\textbf{全部相干态是线性相关的,但是有限(可数)个相干态构成的集合线性无关的}:
\begin{proof}
在复平面上积分我们取极坐标$\mathrm{d}^2z=r\mathrm{d}r\mathrm{d}\theta$,我们为每个相干态加上权重$z^n$后积分发现:
\begin{align*}
\int z^m\ket{m}\mathrm{d}^2z=\sum_{n=0}^\infty\frac{1}{\sqrt{n!}}\ket{n}\int_0^\infty r^{m+n+1}\mathrm{e}^{-\frac{1}{2}r^2}\mathrm{d}r{\color{red}\int_0^{2\pi}\mathrm{e}^{i(m+n)\theta}\mathrm{d}\theta}=0
\end{align*}
可数个相干态集合线性无关可以从作为$z$的函数$\braket{z|z_k}$之间线性无关看出。\qed
\end{proof}
$\bullet$\textbf{全体相干态具有完全性}:
\begin{equation}
\boxed{\frac{1}{\pi}\int\ket{z}\bra{z}\mathrm{d}^2z}
\end{equation}
\begin{proof}
只用证明$\frac{1}{\pi}\int\braket{\psi|z}\bra{z|\varphi}\mathrm{d}^2z=\braket{\psi|\varphi}$即可。
tips:利用$\ket{n}$的完全性以及积分公式$\int_0^\infty r^{k}\mathrm{e}^{-r}\mathrm{d}r=k!$\qed
\end{proof}
$\bullet$\textbf{相干态使得海森堡不确定性关系取等号}:
\[\Delta x\Delta p=\frac{\hbar}{2}\]
\begin{proof}
直接代入相干态波函数计算即可。\qed
\end{proof}
\subsubsection{相干态表象}
既然相干态是归一完备的(但不正交),那么就可以作为Hilbert空间的一组基底将态矢量展开:
\[\ket\psi=\frac{1}{\pi}\int \ket z\braket{z|\psi}\mathrm{d}^2z\equiv\frac{1}{\pi}\int \ket z\psi(z)\mathrm{d}^2z\]
下面具体来求在相干态表象下,产生湮灭算符作用在态函数上的具体形式。
\begin{equation}
\begin{aligned}
\braket{z|A|z^\prime}&=z^\prime\braket{z|z^\prime}=z^\prime \mathrm{e}^{-\frac{1}{2}\left(|z|^2+|z^\prime|^2\right)+z^*z^\prime}\\&=\left(\frac{z}{2}+\frac{\partial}{\partial z^*}\right)\braket{z|z^\prime}
\end{aligned}
\end{equation}
再插入完全性关系:\footnote{这里我们给$A$加上了$\hat{ }$,强调这是函数空间中的算符,前面都是作用在态矢量空间的算符,或者说是矩阵形式讨论的,不过很多情况下我们不会区分这两者,比如前面我们讨论$\hat{p},\hat {x}$都没有刻意强调是态函数空间算符(详见$\S 3$和附录$\S B$)。}
\begin{equation}
\begin{aligned}
\hat {A}\psi(z)&=\braket{z|A|\psi}=\frac{1}{\pi}\int\braket{z|A|z^\prime}\braket{z^\prime|\psi}\mathrm{d}^2z^\prime\\
&=\frac{1}{\pi}\int\left(\frac{z}{2}+\frac{\partial}{\partial z^*}\right)\braket{z|z^\prime}\braket{z^\prime|\psi}\mathrm{d}^2z^\prime\\
&=\left(\frac{z}{2}+\frac{\partial}{\partial z^*}\right)\frac{1}{\pi}\int\braket{z|z^\prime}\braket{z^\prime|\psi}\mathrm{d}^2z^\prime\\
&=\left(\frac{z}{2}+\frac{\partial}{\partial z^*}\right)\psi(z)
\end{aligned}
\end{equation}
对于产生算符,计算相对容易:
\begin{equation}
\begin{aligned}
&A\ket{z}=z\ket{z}\Rightarrow\bra{z}A^\dagger=z^*\bra{z}\\
\Rightarrow&\hat{A}^\dagger\psi(z)=\braket{z|A^\dagger|\psi}=z^*\psi(z)
\end{aligned}
\end{equation}
最终我们得到相干态表象下的产生湮灭算符形式:
\begin{equation}
\boxed{\hat{A}=\left(\frac{z}{2}+\frac{\partial}{\partial z^*}\right)\qquad \hat{A}^\dagger=z^*}
\end{equation}
\section{自由粒子}
自由粒子情况下即$V \equiv 0$, 不难发现定态薛定谔方程和无限深势阱的形式是一样的, 即$$\frac{d^2\psi}{dx^2}=-k^2\psi$$解的形式为$$Ae^{ikx}+Be^{-ikx}$$
注意, \uwave{我们不再有边界条件去直接确定$A$和$B$的值}, \textbf{这个时候定态解的能量的取值并不是离散的!可以取到任何大于$0$的值!}其实本身量子力学就是不排斥
连续性的, 离散可以有很多种, 不一定就表明某个量一定是离散取值。
考虑wiggle-function后写成下面的形式\footnote{这里$\psi_k$无法进行归一化, 我们这里添上系数${\raise0.5ex\hbox{$\scriptstyle 1$}\kern-0.1em/\kern-0.15em
\lower0.25ex\hbox{$\scriptstyle {\sqrt {2\pi } }$}}$是狄拉克意义下的正交归一化的考虑。}:
\begin{lequation}
\boxed{
\Psi_k(x,t)=Ae^{i\left(kx-\frac{k^2\hbar t}{2m}\right)}, \psi_k(x)=\frac{1}{\sqrt{2\pi}}e^{ikx}
}
\end{lequation}
\begin{thinknote}
实际上我们更常见的写法是用:
\[\psi_p(x)=\frac{1}{\sqrt{2\pi\hbar}}e^{i\frac{p}{\hbar}x}\]
作为自由粒子体系的能量表象, 这和动量表象表达式相同。详见下一章讨论。
\end{thinknote}
考虑到最后反正要对解进行线性叠加, 而两项仅仅只在$e$指数上面差了一个负号, 所以我们将负号纳入$k$后得到:
$$k\overset{def}{=}\pm\frac{\sqrt{2mE}}{\hbar}$$观察到每一项就代表一个\uwave{行波}, $k>0$的时候正向传播(向右), 反之负向传播, 不再是前面势阱模型里面的
驻波。
\begin{thinknote}
对于某个确定的振幅对应了一个$x$, 也就是说概率波上面的某一个确定的点, 其$x$和$t$之间满足$x\pm vt=\text{const}$的演化关系, 也即$x=\mp v +\text{const}$。所以可以认为概率
波上的每个点都以$v$的速度在\uwave{平移运动}, 也很容易根据速度前的正负号确定是左行波还是右行波了。顺带一提, 这种形式下的行波, 模长平方(概率密度)与位置
无关(不同位置找到粒子概率都一样), 都等于指数项前面振幅的平方, 这也是后面$\delta$势阱散射态求透射系数的基础。
\end{thinknote}
与经典力学中的波动方程相对应, 前面讲过$k$可以理解为角波数, 那么波速$v=\frac{|k|\hbar}{2m}=\sqrt{\frac{E}{2m}}$, 角频率$\omega=\frac{k^2\hbar}{2m}$, 过会我们会再回到波速的问题上来, 目前
来看似乎有自由粒子速度$v_p=\sqrt{\frac{2E}{m}}=2v_{wave}$。
自由粒子的波函数的定态解与前面两个模型最大的不同就是\textbf{它是无法归一化的}。所以现在我们求出来的定态解完全只是数学上的一个过程\footnote{虽然不存在这样一个定态, 但是我们前面介绍的解薛定谔方程的一般手段仍旧不变, 只是这些“定态”不可归一化且是连续谱}, 没有实际的物理意义
不存在一个状态, 其中自由粒子属于能量始终不变的定态(无论如何测量都是一个值)。但是没事, 虽然不存在定态, 但是我们还是可以使用定态解的线性组合来构造符合初值的
解\footnote{这些合成的解沿用波动学的观点, 称为\uwave{波包}, 是可归一化的}, 这些“定态”解也是归一化和完备的。
我们仍旧对定态解进行线性叠加, 注意是对定态解叠加, 也就是对定态解的波函数进行叠加不是$\psi$而是$\Psi$, 所以不要忘记了每个$\psi$后面的关于时间的指数项。
\begin{lequation}
\label{2.26}
\boxed{
\Psi=\int_{-\infty}^{\infty}\phi(k)\psi_ke^{-\frac{ik^2\hbar}{2m}t}dk=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}\phi(k)e^{i\left(kx-\frac{k^2\hbar}{2m}t\right)}dk
}
\end{lequation}
和\ref{find_c}对比一下就会发现, $\phi(k)$取代了$c_n$, 这一点很好理解, 因为$\phi_k(x)$变成了一个关于$k$的连续函数, 而数列我们也通常称为\uwave{整标函数}, 所
以, 前面变成连续函数来加权, 求和也变成了积分。
初始值还是在$t=0$处给定\footnote{若不是, 你对$t$进行一个换元平移一下计时零点即可}, 那么可以得到:
\begin{lequation}
\Psi(x,0)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}\phi(k)e^{ikx}dk
\end{lequation}
数学上$\Psi$就是对$\phi$的\uwave{傅里叶逆变换}, 下面的定理可以很方便的求出$\phi(k)$。
\begin{theorem}{Plancherel's theorem}
\begin{center}
\begin{math}
\displaystyle
f(x)=\frac{1}{\sqrt{2\pi}} \int\limits_{-\infty}^{\infty}F(k)e^{ikx}dk\Longleftrightarrow F(k)=\frac{1}{\sqrt{2\pi}} \int\limits_{-\infty}^{\infty}f(x)e^{-ikx}dx
\end{math}
\end{center}
\end{theorem}
其中$F(x)$表示$f(x)$的傅里叶变换, 反之即为逆变换。
\begin{lequation}
\boxed{
\phi(k)=\frac{1}{\sqrt{2\pi}} \int _{-\infty}^{\infty}\Psi(x,0)e^{-ikx}dx
}
\end{lequation}
遗憾的是, 很多函数都只能写出积分后进行数值模拟, 无法用基本初等函数表示。
\begin{define}{群速度和相速度}
\begin{itemize}
\item \textbf{群速度}:波包的移动速度
\item \textbf{相速度}:波包里面有很多小峰, 这些小峰的移动速度便是相速度, 可以理解为一个在波上的质点跟随波的平移速度
\end{itemize}
\end{define}
\animategraphics[scale=0.75,autoplay=True,loop]{24}{videos/wave_group/wave_group-}{0}{250}
\begin{center}
上面的动画中, {\color{red}红色}表示相速度, {\color{green}绿色}表示群速度
\end{center}
关于群速度和相速度的讨论应该是波动学的内容, 这里我们不想讨论过多, 只是先说明两个速度的计算公式, 在波包很明显也就是有很好的群速度的定义的时候, 单单从
自由粒子的波动方程可以推出下面这一点。
\begin{theorem}{群速度和相速度的计算}
\begin{lequation}
\boxed{
v_g=\frac{d\omega}{dk}\qquad\qquad v_p=\frac{\omega}{k}
}
\end{lequation}
\end{theorem}
其中, 波包是由一系列定态解叠加而成的, 具有下面的形式(对比\ref{2.26}):
$$\Psi=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}\phi(k)e^{i\left(kx-\omega t\right)}dk,\quad \omega=\frac{k^2\hbar}{2m}$$
正是由于多列行波的叠加才形成了波包, 而每个分量(行波)中的频率又是和$k$相关的, 这样就会导致群速度不等于相速度, 这时我们称之为\textbf{色散}, 自由粒子的
波函数恰好就是多列行波的叠加, $\omega$是关于$k$的二次式, 那么, 代表粒子的群速度是代表行波前进的相速度的两倍这一论断就不难看出了。
\section{\texorpdfstring{$\delta$}.函数势阱}
\subsection{束缚态和散射态}
最好引入束缚态和散射态的方法是使用经典类比, 我们前面提到的谐振子和无限深势阱都是束缚态, 他们每个定态解是可归一化的, 而且是分立谱, 粒子的行为像驻波一
般, 会因为某些特定的边界条件使得波函数被限制在一个局域内。但是后面讲的自由粒子是连续谱, 不存在定态, 定态解不可归一化, 它的行为更像是一种行波, 可以自
由延伸到无穷远处, 但是由定态解的线性组合构成的波包是可以归一化的, 这也就意味着, 一旦我们讨论的是自由粒子这种散射态, 他一定是很多种
定态能量的叠加态, 不可能是测量后只会出现一种能量的某个定态, 一定是很多很多定态的组合, 因为对于散射态定态解是没有物理意义的, 这一点我们后面会再次提到。
从经典力学来看, 由于物体的动能始终是一个定值, 所以在势能图上, 粒子永远不可能达到势能大于粒子总能量的地方(能量守恒), 从图上来看就是粒子被限制在两个点
之间来回折返(\ref{turning_point}), 除非物体的能量大于势能的最大值, 这样粒子就会越过“山峰”, 不断前进。
\begin{figure}[htbp]
\centering
\label{turning_point}
\includegraphics[scale=0.5]{fig/2-5-1.eps}
\caption{经典力学中, 粒子会被“困在”$x=a$和$x=b$之间}
\end{figure}
还是以上面的图像为例, 在经典力学看来(如果初始位置在$ab$之间的话), 这个图像会产生一个束缚态, 在量子力学中, 马上你将会看到, 由于所谓的\textbf{量子隧穿效应}
, 任何有限高的势能都不能阻挡粒子, 粒子完全可以越过最高点大于其总能量的“山峰”。所以上面的图像对于量子理论来说也是一个散射态, 只有像前面所说的无限深势阱和
谐振子, 才会导致粒子的束缚态的产生。
\begin{theorem}{束缚态和散射态的判定准则}
\begin{center}
\begin{math}
\displaystyle
\begin{cases}
E>V(+\infty) \vee E>V(-\infty) \Rightarrow \text{Bound State}\\
E<V(+\infty) \wedge E<V(-\infty) \Rightarrow \text{Scattering State}
\end{cases}
\end{math}
\end{center}
一般情况下, $V(x)$会在$x\to\pm\infty$时趋近于$0$(谐振子和无限深势阱就是特例), 这个时候可以再次简化判定准则为:
\begin{center}
\begin{math}
\displaystyle
\begin{cases}
E>0 \Rightarrow \text{Bound State}\\
E<0 \Rightarrow \text{Scattering State}
\end{cases}
\end{math}
\end{center}
\end{theorem}
注意, 我们前面在求解薛定谔方程时, 都说过$E>V_{min}$, 很遗憾, 这个条件在散射态时, 不起作用, 因为散射态我们并不严格要求定态解一定是可归一化的, 所以这
个条件只能对于束缚态使用, 前面我们对谐振子和无限深势阱使用这个条件是完全正确的, 因为无论$E$是大于$0$还是小于$0$, 这两个势阱产生的都是束缚态, 那么束缚
态的定态解必须是可归一化的, 这也就间接要求了$E$的非负性。事实上, 你如果在$E<0$的假定下求解这两种模型的解, 你会发现无论是定态还是线性组合后的解都是不可归一化的
, 没有实在的物理意义。 在这个时候我们就说$E>0$的解已经构成了一个\uwave{完备集}, 初始的波函数一定是可以归一化的, 也一定是可以由对应的定态解组合而成, 然而
$E<0$的解都不可以归一化, 所以这两种势能模型中仅当粒子$E>0$时有物理意义, 这一点实际上根据经典力学我们可以很容易类比得出。
\begin{history}{依赖于直觉的推理}
下面我们用一种偏向于物理直觉的方法来理解束缚态和散射态。前面提到了束缚态导致离散谱, 散射态导致连续谱, 实际上束缚态也可以\uwave{大致}上理解为波函数局限在某个有限范围内,
就像高斯波包一样。那么如果束缚态的谱线是连续谱, 定态解用连续的参数$\alpha$标记:$\ket{\psi(\alpha)}$。而且考虑参数微小的连续变化, 态也会产生一个微小的连续变化:
\[\ket{\psi(\alpha+\delta\alpha)}\approx\ket{\psi(\alpha)}+\ket{\delta\psi}\]
由于哈密顿算符是厄米算符, 那么$\ket{\psi(\alpha+\delta\alpha)}$和$\ket{\psi(\alpha)}$应该相互正交。\footnote{这个是属于Dirac意义下的正交归一化, 问题不大, “正交”和离散情况还是一样的, “归一”改了改。}
\[
\braket{\psi(\alpha+\delta\alpha)|\psi(\alpha)}=\braket{\psi(\alpha)|\psi(\alpha)}+\braket{\delta\psi|\psi(\alpha)}
\]
上式中第一项是正的\footnote{这里我并不想直接说它是“归一”的, 虽然这里只是一种直觉上的推理, 但是我还是想更一般的去讨论它。}。而后面的一项是一个任意小的值, 从波函数空间内的
内积定义出发将它写成积分形式可以更好的看出这一点, 被积函数是任意小的, 而积分区间又因为是束缚态被限制了(也即波函数集中在一个有限区间上), 所以总的来说如果束缚态也可以有连续谱的话, 是违背了
“正交性”的, 所以束缚态都是离散谱。而散射态积分区域就是全空间, 就没有这个问题了。当然, 体系也可以是部分连续谱部分散射谱, 下面要讲的$\delta$势阱就是一类, 还有比如氢原子的电离和非电离情况也是一类。
\end{history}
\subsection{\texorpdfstring{$\delta$}d函数势阱}
注意这里我们讨论的是\uwave{势阱}, 主要是为了和后面的\uwave{势垒}区分, 你可以将势阱想象成一口井, 但势垒是一座山。
\subsubsection*{狄拉克$\delta$函数}
这个函数的定义你可以立即为将克罗内克符号扩充为一个连续函数, 也可以理解为是以个函数序列的极限。
\begin{define}{$\delta$ function}
\begin{center}
\begin{math}
\displaystyle