记得我小学和中学的时候，由于是在天主教学校就读，每天早上都要念经唱歌赞美主。我当时就想，这个主是什么玩意，怎样跟历史书中所形容的昏君一样，这么喜欢听这种阿谀奉承之词呢。

也不知道是有意还是无意，学校的生物课安排得非常不认真，找了个教化学的老师，用二十年前出版，又贵又重的课本来上生物课。到最后，由于这样种种原因（还有我太懒），生物课基本上一点没有学会，突别是进化论。Dobzhansky讲过，“Nothing in biology makes sense except in the light of evolution”，没有学会进化论，就无法真正学会生物学。

不过由于我这个人一直特别讨厌拍马屁和装模作样的东西，所以学校那些赞美主的仪式，反而让我不可能去信天主教。你不是每天强迫我赞美主吗？我就天天咒骂它，于是每天我都过得很高兴，把怨气都发泄到主身上了。

然而，这种反对宗教的立场是不牢固的。缺少了科学思想的洗礼，就无法看透这些愚昧迷信的玩意的本质，不能做到举一反三。所以在当时，我尽管反对天主教，但对中国那些传统的迷信活动，如风水，气功和中医，还有其他什么武功，人体特异功能，UFO，世界XX个未解之迷等，都抱有一种“可能有这么一回事”的想法。

不过在科学面前，任何迷信活动都不堪一击。当年宗教就是人生的全部，宗教是社会规范，是道德标准，也解释世间万物。后来科学出现了，宗教就节节败退，到最后穷途末路了，就声称宗教和科学是两个不同的领域，企图把宗教和现实分离出去，宗教是所谓的“信”，是faith，是另外一个领域。可笑的是，宗教却处处企图“回到现实”，不断曲解科学发现，污蔑进化论和宇宙大爆炸理论，别人笑他地球只有6000年的理论，就气到跳起来。可以想象如果有一天科学证实了地球是宇宙的中心（尽管已经证明不是了……），宗教界人士一定欢恩鼓舞，说：“你看，科学不是证明了我们是正确的了吗？”完全忘了当时所谓的“科学和宗教是不同领域”这种推诿之词了。

其他迷信活动也是类似的。科学太强，又太无情与残酷了，一次又一次戮破了人们的美梦。如果迷信活动单单是自己占了自己的老家，宣布独立，成立XX共和国，自己做大总统，那么就随便你。但那些迷信活动要生存下去（可能也是进化论在驱使他们这么做），总免不了要骗人。只要是骗人的东西我就要反对，暨大的小朋友们被我骂了，只能说句不好意思了。

好了，就写这么多，有时间再写写这方面的感想。

CORRESPONDENCE:
SENT    RECEIVED    DESCRIPTION
01May08         Editorially approved for publication
10Apr08 30Apr08 Review request to referee; report received
30Apr08         Reminder to referee [others sent (not shown) at 1-2 week intervals]
07Apr08         Correspondence (miscellaneous) sent to author
20Mar08 07Apr08 Ed. decision and/or ref. comments to author; response rcvd
31Jan08 11Mar08 Review request to referee; report received
31Jan08 26Feb08 Review request to referee; report received
21Feb08         Reminder to referee [others sent (not shown) at 1-2 week intervals]
21Feb08         Reminder to referee [others sent (not shown) at 1-2 week intervals]
30Jan08         Acknowledgment sent to author
29Jan08         Correspondence (miscellaneous) sent to author

整個過程非常順利,審稿人第一次的意見就沒提出甚麼批評.審稿人A婉轉地提出要引用他的幾篇文獻,審稿人B則提出
了幾點針對文中presentation的更改意見.審稿人A見我引了他的文獻,於是就二話不說,放過我了~
雖然我以前都發過一篇文章在PRB,但那是comment,不是正式的文章.另外,我還有一篇第二作者(導師第一)的文章
投到PRB的brief report去了,現在搞了N久,審稿人又說最好寫成regular paper,不知這可憐的文章還要搞多久才
出能出來(搞了足足一年了)...

6月2日新增:文章出来了，可看Phys. Rev. B 77, 174305 (2008).

近来试用Maple 11的Physics Package，主要是关于量子力学计算的部份。不知道是我的水平低，还是Physics Package的功能局限，我觉得算一个么正变换，感觉用手算比用Maple 11算要方便……Maple 11那个Physics Package好像有许多古古怪怪的问题，我也搞不清楚。以下计算在压缩态|s> = SD|0>下算符Y1 = a exp(-i * phi)) + a^+ exp(i * phi)的不确定度sqrt( - ^2)，其中S是压缩态算符，D为位移算符，|0>为真空态，a和a^+分别为湮灭和产生算符，phi为实数。我的计算过程非常繁琐，希望高人多多指点，给我一些简单的解决方法。相关计算过程请参考D.F.Walls, G.J.Milburn, Quantum Optics, (Springer-Verlag, Berlin, 1994).
> restart;
> with(Physics):
设S和D为量子算符：
> Setup(quantumoperator = {S, D});
定义湮灭和产生算符：
> am := Annihilation(n, 1); ap := Creation(n, 1);
                         [quantumoperator = {D, S}]
                                     a-
                                     a+
定义几个么正变换的结果：
> DamD := Dagger(D) . am . D = am + alpha;
> DapD := Dagger(D) . ap . D = ap + conjugate(alpha);
> SamS := Dagger(S) . am . S = am * cosh(r) - ap * exp(2 * I * phi) * sinh(r);
> SapS := Dagger(S) . ap . S = ap * cosh(r) - am * exp(-2 * I * phi) * sinh(r);
还要设定S的么正性：
> SUI := Dagger(S) . S = 1;
                     ((Dagger(D)) . a-) . D = a- + alpha
                     ((Dagger(D)) . a+) . D = a+ + alpha
        ((Dagger(S)) . a-) . S = a- cosh(r) - a+ exp(2 I phi) sinh(r)
       ((Dagger(S)) . a+) . S = a+ cosh(r) - a- exp(-2 I phi) sinh(r)
                             (Dagger(S)) . S = 1
定义真空态：
> VS := Ket(n, 0);
                                  Ket(n, 0)
定义Y1：
> Y1 := am * exp(-I * phi) + ap * exp(I * phi);
                       a- exp(-I phi) + a+ exp(I phi)
准备计算S^+ D^+ Y1 D S：
> tmp := Dagger(S) . (Dagger(D) . Y1 . D) . S;
        exp(-I phi) ((Dagger(S)) . (((Dagger(D)) . a-) . D)) . S + exp(I
phi) ((Dagger(S)) . (((Dagger(D)) . a+) . D)) . S
把相关么正变换的结果代入去：
> tmp := subs({DamD, DapD}, tmp);
              exp(-I phi) ((Dagger(S)) . (a- + alpha)) . S
                        /              /          \\   
           + exp(I phi) \(Dagger(S)) . \a+ + alpha// . S
> tmp := subs({SamS, SapS, SUI}, tmp);
       exp(-I phi) (a- cosh(r) - a+ exp(2 I phi) sinh(r) + alpha)
                       /                                             \
          + exp(I phi) \a+ cosh(r) - a- exp(-2 I phi) sinh(r) + alpha/
便计算完了：
> ry1 := Dagger(VS) . tmp . VS;
                    exp(-I phi) alpha + exp(I phi) alpha
的计算过程也是类似的，也要定义一大堆么正变换的结果，如要令D^+ a a D = D^+ a D D^+ a D，非常麻烦，计算结果太长的我就不作显示了。
> tmp := Dagger(S) . (Dagger(D) . (Y1 . Y1) . D) . S;
> DamapD := Dagger(D) . (am @ ap) . D = (Dagger(D) . am . D) . (Dagger(D) . ap . D);
> DapamD := Dagger(D) . (ap @ am) . D = (Dagger(D) . ap . D) . (Dagger(D) . am . D);
> DamamD := Dagger(D) . (am . am) . D = (Dagger(D) . am . D) . (Dagger(D) . am . D);
> DapapD := Dagger(D) . (ap . ap) . D = (Dagger(D) . ap . D) . (Dagger(D) . ap . D);
> tmp := subs({DamapD, DapamD, DamamD, DapapD}, tmp);
> tmp := subs({DamD, DapD}, tmp);
> SamapS := Dagger(S) . (am @ ap) . S = (Dagger(S) . am . S) . (Dagger(S) . ap . S);
> SapamS := Dagger(S) . (ap @ am) . S = (Dagger(S) . ap . S) . (Dagger(S) . am . S);
> SamamS := Dagger(S) . (am . am) . S = (Dagger(S) . am . S) . (Dagger(S) . am . S);
> SapapS := Dagger(S) . (ap . ap) . S = (Dagger(S) . ap . S) . (Dagger(S) . ap . S);
> tmp := subs({SamapS, SapamS, SamamS, SapapS}, tmp);
> tmp := subs({SamS, SapS, SUI}, tmp);
把算出来了：
> ry2 := Dagger(VS) . tmp . VS;
不确定度sqrt( - ^2)的计算就是好几个代数化简：
> expand(ry2 - ry1^2);
                                           2          2
                  -2 sinh(r) cosh(r) + sinh(r)  + cosh(r)
> sqrt(factor(%));
                                               (1/2)
                         /                  2 \    
                         \(cosh(r) - sinh(r)) /    
> convert(%, exp);
                                          (1/2)
                              /         2\    
                              \(exp(-r)) /    
> simplify(%) assuming r >= 0;
                                   exp(-r)
于是得出了和Walls书中相同的结果。可见，这种计算过程非常麻烦，主要原因有几个：１、不知道怎么定义算符的么正性，使得如S^+ S这种本该等于1的结果也要自己去做代入；２、不知道如何定义么正变换（或其他什么变换）的结果。如果能用Setup函数去定义D^+ a D的结果，并且能自动根据D的么正性得出D^+ a a D = D^+ a D D^+ a D，免得又要自己设定一个变量去做代入。这种人为代入很容易发生错误，而且设定这种代入式子也很麻烦；３、如果直接算，而不是先算出S^+ D^+ Y1 D S再代入<0|..|0>内，会发生代入不了那些么正变换结果的情况。具体把式子代入S^+ D^+ Y1 D S时，也会发生一些奇怪的代入不成功的情况，使得代入过程要分两步去做。