Re: [閒聊] 為什麼版上幾乎沒有討論非線性?

看板Physics作者 (如果狗狗飛上天)時間9年前 (2015/04/23 05:41), 9年前編輯推噓10(10036)
留言46則, 13人參與, 最新討論串3/5 (看更多)
※ 引述《peter308 (pete)》之銘言: : 我覺得非線性是非常有用的一門課 : 裡面的數學很簡單 大學生甚至高中生就會 : 可是可以應用的領域五花八門 : 從訊號 工程 生物 經濟 社交 化學 雷射 半導體 電漿 統計物理 相變 重整化群 : 等等都有 : 可是怎麼感覺版上很少人討論? 一點想法,其實我認為非線性絕對是物理中一個關鍵的部分.而且徹底的滲透進了 各個子學門中,只是非線性一詞意義太過廣泛,因此除非是專門研究數學的,不然一 般人很少會用非線性來描述自己的問題,即便他所做的事情就是非線性. 以凝態物理來說,類似的例子早就已經數不清了.例如透過自洽場求解密度泛函這 件事情本身就是在求解一個非線性方程.但是即便是當初這個工作的幾位主要貢獻 者也幾乎不會用非線性來描述自己的工作,縱然他們在做的事情本質上就是在求解 一個非線性方程. 舉一個典型的例子,不曉得版上有沒有人接觸過CPMD這種計算? 他的全名叫做 Car-Parrinello Molecular Dynamics,這是一個在凝態計算裡面很常用來求解 分子或是晶格最穩定結構的方法. 他主要處理的問題是,如果我們想知道求解 一個量子系統的結構隨時間的演化,最直覺的作法就是我先做電子結構的基態自 洽,求出了電子的基態之後,我們再去計算這個基態給原子的力.得出力之後,根據 這個力對原子給出位移,然後有了新的位移,我們再根據這個位移去重新做自洽求 出新的原子受力.所以整個流程就是: 自洽求基態->算受力->產生新位移->自洽求基態->算受力>產生新位移.... 這個做法想起來很簡單,但其實整個過程極其困難.原因是因為,你每把原子的位置 推動一小步,你都要重新求一次自洽.如果你把這整個過程所涉及到的方程式全部 寫下來,就會發現他不僅涉及了求解非線性方程,還涉及了自洽求解本徵值的問題. 整個問題運算量極大,而且每一次自洽都存在了可能基態求不對的風險使得誤差 不斷累積. 這在80年代,基本上是個不可解的問題. 直到Car跟Parinello兩個人的出現改變了這個局面.CP的想法是,他在原本的 Hamiltonian中,"毫無道理的"(是的,完全毫無道理,純粹是數學上瞎湊) 引進了一個跟電子密度有關的項(後來這個項稱作虛擬動能),但因為是胡亂 添加進來的項,所以它們刻意地在這個項前面乘上了一個很小的係數m,那麼我們 可以預期,只要m足夠小,小到可以忽略的程度.求解這個多出一項的Hamiltonian 問題,就會非常的近似於原本我們想求解的問題.(有興趣的人可以wiki一下) 這樣有什麼好處呢? 神奇的地方在於,這個"多了一項"的Hamiltonian,竟然可以把 電子的orbital的時間演化在"數學形式上"變成了一個古典的牛頓方程. 只要我的 m取的足夠小,整個問題就從一個需要不停地求自洽的非線性問題,變成了一個求解 一堆相互耦合的牛頓方程問題.再也不需要做自洽了!至此,運算被大幅減少,透過 數值運算量子系統的空間動力學問題才有了一個在具體上可解決的方法. 從此之後,這個方法被不斷的改進,加強.時至今日,你會在大部分的ab initio套裝 軟體裡面發現內建這個功能了.(當然,現在的計算能力,暴力法硬解也是很普通了, CP已經不是非靠它不可了,但許多量子結構動力學的計算還是基於CP). CP兩人也因為這個工作聞名世界,時至今日,當初他們提出的那篇文章已經被引 用了近萬次,而且每個月都還以數百次的數目再增加.可以說靠著一個工作,就奠定 了一代大師的地位.剩下的日子完全可以悠悠哉哉的看著citation破表. 我記得我當初在讀這個工作的時候,看到它們莫名其妙地引進一項沒意義的項進來 ,最後竟然可以把整個問題從量子的本徵值求自洽解的問題變成了一個牛頓力學的 問題時,幾乎感到不可置信,暗自拍案這簡直是神人的工作.甚至覺得他們從天上掉 一項下來竟然可以有這麼美妙的結果實在太沒道理了.是哪來的靈感? 但科學 常常就是這樣,努力固然很重要,但是那些最重要的工作,往往帶有著一點"天啟"的 成分. 話說從頭,我想說的是,其實我不覺得非線性在物理的領域裡面是一個冷課題,相反 的,很多重要的工作,往往都是建築在破解了一個非線性的問題上達成的.甚至你會 發現這些破解的方法後來往往造成了一個新的學門的產生(例如今日各種應用在DFT 上的方法都可以算是).只是這些問題往往不是以非線性的姿態出現,而是以一個工 具的姿態出現就是了. 美國的8,90年代,曾經有過一段物理人跑華爾街的風潮,我猜想很大部分就是因為 大多數的物理人都或多或少的被潛移默化許多非線性的觀念,即便這些觀念在物理 中多半不是以一門專門的學門來教育,但每個物理的領域或多或少都有一些屬於 他們自己的非線性問題可以移植到金融領域做分析. 不過之前看過一篇分析,認為物理人跑華爾街的時代已經過去了(以我身邊的觀察 來看,應該也確實如此,不過畢業跑去網路公司寫code,做data science的人倒是越 來越多了),很大的原因是因為那些曾經對金融業來說很新,很陌生的數學方法,早 已被玩金融的人學去了(畢竟人家也不比你笨),甚至有系統的發展成了屬於他們專 屬的計量學門,以至於其他領域的人毫無背景的跳進去,不僅學不來,也玩不贏了. -- ★人生中最溫暖的夏天是在紐約的冬天 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 128.200.93.30 ※ 文章網址: https://www.ptt.cc/bbs/Physics/M.1429738871.A.92C.html ※ 編輯: pipidog (128.200.93.30), 04/23/2015 05:50:53

04/23 06:55, , 1F
其實去華爾街的物理人還是有 只是華爾街要的是名校的
04/23 06:55, 1F

04/23 08:07, , 2F
聽起來好有趣 請問多的CP項後來有發現什麼物理意義嗎?
04/23 08:07, 2F

04/23 08:07, , 3F
可以將方程式簡化成這樣,還是就真的只是數學?
04/23 08:07, 3F

04/23 08:19, , 4F
純數學,完全沒有任何物理意義,但是讓問題變得可解了
04/23 08:19, 4F

04/23 08:25, , 5F
所以跑這類計算,你必須讓那個m非常小,結果才靠得住.
04/23 08:25, 5F

04/23 08:38, , 6F
@yarnball,其實沒錯,物理人跑華爾街的還是有,我有幾個朋
04/23 08:38, 6F

04/23 08:39, , 7F
有前陣子也剛面試完金融業的缺,但是他們現在要人,主要也
04/23 08:39, 7F

04/23 08:40, , 8F
是希望他的coding強,做高頻交易,對於物理則完全不在乎了
04/23 08:40, 8F

04/23 08:43, , 9F
所以物理的背景感覺加不了分,倒是CS背景走到哪都吃香
04/23 08:43, 9F

04/23 08:47, , 10F
每次上物理版最期待的就是 pipidog 大的文章
04/23 08:47, 10F

04/23 10:09, , 11F
你文章說的m取很小 是不是就是一種分離不同time scale的
04/23 10:09, 11F

04/23 10:10, , 12F
技巧?? strogaz 的 3.5 節 很像你在講的東西
04/23 10:10, 12F

04/23 10:11, , 13F
這個可以延伸到center manifold theory 或是slaving
04/23 10:11, 13F

04/23 10:11, , 14F
principle....
04/23 10:11, 14F

04/23 10:13, , 15F
自洽解我最近的感覺是很像 1-D maps.....
04/23 10:13, 15F

04/23 10:14, , 16F
1-D maps 有chaos 你的解就會有non-SCF convergene問題
04/23 10:14, 16F

04/23 12:53, , 17F
老實說,我根本不懂那是啥技巧,只覺得是天上掉下來的 @@~
04/23 12:53, 17F

04/23 23:20, , 18F
這位學長是作計算嗎@@
04/23 23:20, 18F

04/24 01:38, , 19F
04/24 01:38, 19F

04/24 04:29, , 20F
04/24 04:29, 20F

04/24 06:32, , 21F
推,我覺得這種靈感還是來自於一些數學上的背景知識
04/24 06:32, 21F

04/24 06:32, , 22F
如果數學有做過其他類似理論就會觸類旁通...
04/24 06:32, 22F

04/24 06:33, , 23F
當然也可能真的是蘋果砸到頭 XDDDD
04/24 06:33, 23F

04/24 10:29, , 24F
push~Everytime pipidog's posts have broadview
04/24 10:29, 24F

04/24 17:41, , 25F
我同意這整個點子很精彩,但是那一項本身沒有那麼神啦
04/24 17:41, 25F

04/24 17:41, , 26F
加進去的就是個標準的inertia term
04/24 17:41, 26F

04/24 17:43, , 27F
沒有inertia的系統有各種不連續,這不是什麼新聞。所以讓解
04/24 17:43, 27F

04/24 17:43, , 28F
析解簡單的這個性質,同時也可能讓數值解難搞到哭爹喊娘
04/24 17:43, 28F

04/24 17:45, , 29F
丟個inertia term進去做regularization是很直覺的作法
04/24 17:45, 29F

04/24 17:45, , 30F
這兩人了不起的地方在於,他們認識到你每一步都要重新求電
04/24 17:45, 30F

04/24 17:46, , 31F
子波函數的原因,正是Born-Oppenheimer approx這個
04/24 17:46, 31F

04/24 17:47, , 32F
inertia-less的假設引進的不連續性
04/24 17:47, 32F

04/24 22:17, , 33F
這裡說的inertia-less是指電子分布的取得跟原子核不同,
04/24 22:17, 33F

04/24 22:18, , 34F
並非透過求位移取得?
04/24 22:18, 34F

04/24 22:28, , 35F
Born-Oppenheimer講的是電子會對原子核位置瞬間作出反應
04/24 22:28, 35F

04/24 22:31, , 36F
像是一個只有電阻的電路,電流會瞬間隨著電壓改變
04/24 22:31, 36F

04/24 22:31, , 37F
這就是我所謂的沒有慣性
04/24 22:31, 37F

04/24 22:36, , 38F
啊,電流的例子不好,這個怎樣也沒有inertia。改成電容器和
04/24 22:36, 38F

04/24 22:36, , 39F
電荷好了。
04/24 22:36, 39F

04/24 22:38, , 40F
算 Q = CV 當然很簡單,但是如果你只准用數值積分算Q(t),
04/24 22:38, 40F

04/24 22:39, , 41F
那就會到處爆炸給你看
04/24 22:39, 41F

04/25 11:07, , 42F
這個「瞬間」是指求H本徵值解的時候沒有時間因素?
04/25 11:07, 42F

04/25 22:16, , 43F
不懂非線性的搞出LTCM破產 (創辦人還是經濟諾貝爾獎)
04/25 22:16, 43F

04/25 22:18, , 44F
James Simons數學家開的Renaissance Technologies還是很
04/25 22:18, 44F

04/25 22:18, , 45F
賺啊
04/25 22:18, 45F

05/01 14:34, , 46F
推pipidog 的討論和分享
05/01 14:34, 46F
文章代碼(AID): #1LE1Ltai (Physics)
討論串 (同標題文章)
文章代碼(AID): #1LE1Ltai (Physics)