[問題] 有人自己寫過DFT的subroutine嗎??

看板Physics作者 (pete)時間15年前 (2010/12/24 18:46), 編輯推噓6(6推 0噓 25→)
留言31則, 11人參與, 6年前最新討論串1/1
因為一直以來 我都是用套裝軟體做計算 但最近想說自己把dft的計算 用f77寫一個小程式 當作一個練習 不知道有沒人有這方面的經驗 因為我雖然念過dft 知道k-s equation 但要從零開始 自己寫出dft的code 忽然不知道麼下手 我是想先計算最簡單的 氫原子系統就好 就一個質子和一個電子 能不能請各位給各方向 要怎麼下手?? 一些參考資料或建議也感謝不盡 ~~ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 219.69.110.166 ※ 編輯: peter308 來自: 219.69.110.166 (12/24 18:46)
→
12/24 18:49, , 1F
12/24 18:49, 1F
推
12/24 19:25, , 2F
還以為是discrete fourier transform
12/24 19:25, 2F
推
12/24 21:42, , 3F
推樓上
12/24 21:42, 3F
→
12/24 22:07, , 4F
看到FT只想得到fourier transform @@
12/24 22:07, 4F
→
12/24 22:11, , 5F
你應該先弄清楚你要的exchange-correlation functional
12/24 22:11, 5F
→
12/24 22:20, , 6F
基隆大 很清楚 LDA(Perdew-Zunger parameterisation )
12/24 22:20, 6F
→
12/24 22:22, , 7F
主要是怎麼弄成一個小程式 一些numerical 細節還在摸索
12/24 22:22, 7F
→
12/24 22:23, , 8F
以前用套裝軟體用習慣了 現在才發現根基很薄弱
12/24 22:23, 8F
→
12/24 22:28, , 9F
density integrtion, 解K-Sequation 和 SCF算自洽解
12/24 22:28, 9F
→
12/24 22:29, , 10F
主要是這三部份.....我可能要先看看別人怎麼寫找靈感吧
12/24 22:29, 10F
→
12/24 22:29, , 11F
但目前套裝軟體dft的subroutine都至少上百個 哈哈...
12/24 22:29, 11F
→
12/24 22:30, , 12F
可能要花些時間慢慢看吧....不過我想這過程應該無法逃避
12/24 22:30, 12F
→
12/24 22:31, , 13F
如果以後想做ab initio 方面的研究的話
12/24 22:31, 13F
推
12/25 02:03, , 14F
先定好你的 basis 跟跑出 Hamiltonian 吧
12/25 02:03, 14F
→
12/25 02:05, , 15F
我寫過的玩具, 主要是被處理晶體對稱的地方搞到頭大
12/25 02:05, 15F
→
12/25 02:06, , 16F
如果你只跑單原子就比較不怕這問題
12/25 02:06, 16F
推
12/25 10:00, , 17F
12/25 10:00, 17F
→
12/25 15:21, , 18F
同意基隆男,我也覺得寫自洽不難,晶體的對稱性才頭大
12/25 15:21, 18F
→
12/25 15:22, , 19F
尤其是有磁性的系統,往往某些相的磁性周期性很複雜,簡直無
12/25 15:22, 19F
→
12/25 15:22, , 20F
從下手,完全不知道k該從哪邊走起
12/25 15:22, 20F
推
12/25 16:54, , 21F
這不算小練習耶，如果是第一個程式不建議從這個開始。
12/25 16:54, 21F
→
12/25 22:07, , 22F
我覺得想玩自洽,其實從Hartree-Fock玩起就夠啦,最簡單的像
12/25 22:07, 22F
→
12/25 22:08, , 23F
純C組成的,例如graphene只要考慮一個orbital,高興的話還可
12/25 22:08, 23F
→
12/25 22:08, , 24F
以加個U進去看看,可以自洽做得很爽.
12/25 22:08, 24F
推
12/25 23:17, , 25F
huckel --> hubbard --> PPP --> CNDO --> ZINDO
12/25 23:17, 25F
→
01/19 10:08, , 26F
感謝!! 我會試著寫看看
01/19 10:08, 26F
→
08/13 16:00, , 27F
以加個U進去看看,可以 https://noxiv.com
08/13 16:00, 27F
→
09/17 14:00, , 28F
以前用套裝軟體用習慣 https://daxiv.com
09/17 14:00, 28F
→
11/09 12:12, , 29F
如果你只跑單原子就比較 https://noxiv.com
11/09 12:12, 29F
→
01/02 14:38, 7年前 , 30F
主要是怎麼弄成一個小程 https://muxiv.com
01/02 14:38, 30F
→
07/06 22:47, 6年前 , 31F
尤其是有磁性的系統,往 http://yofuk.com
07/06 22:47, 31F
‣ 更多 peter308 的文章
文章代碼(AID): #1D57dlx9 (Physics)