作者yuantsai (芋圓好吃)
看板Physics
標題[問題] WKB Reflection Coefficient計算
時間Tue Sep 1 21:42:36 2020
╭──────────── 提醒：版規三及版規四 ─────────────╮
│解題文需附上自己的解題過程或自己對題目的理解 │
│問題文需先自行查過網路資料並附上對資料的初步理解 │
│問題獲得回答後也勿任意刪除文章，並適度的向回答者表達感謝之意 │
╰──────────────────────────閱讀後可用ctrl+y刪除╯
【出處】(習題或問題的出處)
Michael Berry 半古典近似文章
【題目】(題目的文字敘述，如有圖片亦可提供圖片)
大家好，最近在看WKB相關論文對波函數 reflection coefficient 的計算，考慮一個粒子(或波)
從無窮遠處靠近一個 barrier，以potential V(x)表示。用WKB 近似的波函數可以得到
關於 reflection coefficient 的積分如下圖
https://imgur.com/BqvKIIo
其中 p 指得是粒子動量 p(x) = √2m[E-V(x)]，E是粒子的能量，m是粒子質量
這邊作者為了避免根號在複變會有多值的問題，直接令一新變數w，接著處理w-domain的
複變積分如下圖(2.30) 和 (2.31)式
https://imgur.com/R3PVibF
WKB 近似最大的問題就是在反曲點時 p=0，而這會造成上述被積分函數爆掉，所以作者
考慮複變積分，假設發生在Wj這個點時，考慮會環繞Wj的封閉積分路徑得到(2.32)
1.這邊的封閉積分路徑作者並沒有說明，我的認為是將原本從負無限大到正無限大的在實
軸上的積分做向上變形並形成無限大半徑的上半圓和包住奇異點Wj的小圓兩個積分，小
圓能得到(2.32)式，但他是怎麼確定 p'/2p 的積分在無線遠處一定會收斂到0的 ?
因為p'/2p 可寫成 (-1/4)[ V'(x)/(E-V)]，這樣來看在無窮遠處是否收斂至0決定在
分母的V(x)，偏偏分子是V'(x)，分母沒有高於分子兩次以上，所以在無窮遠處應該不
等於0，還是說我想的積分路徑有問題 ?
接下來我認為作者後面是轉回原本的x-domain來處理積分
2. 接著作者假設奇異的反曲點是 p^2 的簡單零點，我的認知是 p^2 (x) 可寫為
p^2(z) = (z-a)g(z) 其中a 是奇異點，g(z)是解析函數，所以 p'/2p 這個函數
可用 p^2 表示為 (p^2)'/(4p^2)，再將上式代入可得
(p^2)'/(4p^2) = (1/4){ [-a*g(z) + (z-a)g'(z)]/[(z-a)g(z)] }
上式乘上(z-a)並取極限 z=a 可得留數
R(a) = -a/4
所以得到(2.32)積分為 e^(2iWj)* 2πi*(-a/4) = e^(2iWj)*(-πai/2)，但我推出的
這個積分值和(2.34)式不同，不曉得作者是如何推出來的? 奇怪的是(2.34)式沒有包
含特定的奇異點，所以積分是independent of poles ? 這樣不就表示跟 potential
V(x)無關，這是我認為奇怪的點!
這兩個問題卡了好幾天了，還請各位前輩賜教! 謝謝!
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 110.50.159.229 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Physics/M.1598967758.A.670.html
※ 編輯: yuantsai (110.50.159.229 臺灣), 09/01/2020 21:49:30
1F：→ congeebone: 2. 你可能不小心搞錯了，(2.32) 的環積分是對 w 這個 09/04 11:54
2F：→ congeebone: 變數來做，而不是用你這邊的 z (也就是論文裡寫的 x) 09/04 11:54
3F：→ congeebone: 來做，所以 residue 也不會是你得到的那樣。然後 (2 09/04 11:54
4F：→ congeebone: .34) 怎麼會和 potential 無關？不然這些 complex tu 09/04 11:54
5F：→ congeebone: rning points w_j 和什麼有關？我有點不太確定你最後 09/04 11:54
6F：→ congeebone: 那幾句話的意思，因為最後不是就是把這些 pole 的貢 09/04 11:54
7F：→ congeebone: 獻 R_j 加起來嗎？ 09/04 11:54
8F：→ yuantsai: 我可能沒有說清楚，我的意思是2.34的積分值最後僅保留 09/04 12:41
9F：→ yuantsai: 到2.32的e^2iWj ，而後面被積分函數最後積出來得到iPI/3 09/04 12:41
10F：→ yuantsai: （這邊並不清楚作者是如何積的）表示後面的積分值和tur 09/04 12:41
11F：→ yuantsai: ning point 無關 不知道這樣認知是否有誤會 09/04 12:41
12F：→ congeebone: 我好像懂你意思了，但其實就是 (dp/dw)/2p 在 p^(z) 09/04 13:36
13F：→ congeebone: 有 simple zero 的點上有 simple pole 1/6w，dp/dw 09/04 13:36
14F：→ congeebone: 和 p相除所以沒有出現與 p(w)相關的係數也很正常。另 09/04 13:36
15F：→ congeebone: 外，我突然發現你 (z-a)g(z) 的微分算錯了，所以得到 09/04 13:36
16F：→ congeebone: 的答案 -a/4 因次上就不對了。 09/04 13:36
17F：→ congeebone: ...扯遠了，你還是先想看看怎麼把 p(x) 換成 p(w)，i 09/04 13:36
18F：→ congeebone: pi/3是怎麼出來的也就很簡單了。 09/04 13:36
19F：→ congeebone: 打錯...是 "p^2(z) 有 simple zero 的點上" 09/04 13:38
20F：→ yuantsai: 對 我的微分算錯了，residue 應該是1/4才對,但這樣也不 09/05 11:49
21F：→ yuantsai: 是1/6 雖然我還是在x domain 做計算,但應該不管是x 還是 09/05 11:49
22F：→ yuantsai: w 積出來的值是一樣的才對 ，不知道我哪裡miss掉了 謝 09/05 11:49
23F：→ yuantsai: 謝！ 09/05 11:49
24F：→ congeebone: 你要這樣算也可以，但你用x當變數，積分的路徑就不是 09/05 14:19
25F：→ congeebone: 完整的環積分了，以turning point為原點附近，△w～( 09/05 14:19
26F：→ congeebone: △x)^(3/2)，所以w積分繞turning point一圈，對應到x 09/05 14:19
27F：→ congeebone: 只有繞2/3圈，而因為這邊只是simple pole，你可以直 09/05 14:19
28F：→ congeebone: 接(2/3)*(1/4)=1/6。 09/05 14:19
29F：→ congeebone: 我原本不是這樣想，而是直接把w和x的關係泰勒展開， 09/05 14:19
30F：→ congeebone: 就可以得到 (dp/dw)/2p 的 Laurent series 的 1/w 項 09/05 14:19
31F：→ congeebone: 是 1/(6w)，所以積分的答案是 2*pi*i*(1/6)。 09/05 14:19
32F：推 congeebone: 更正：我上面很多寫1/w的地方我想講的是1/(w-w_j) 09/05 18:04