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

[peter308 (pete)]

  因為一直以來 我都是用套裝軟體做計算

  但最近想說自己把dft的計算 用f77寫一個小程式 當作一個練習

  不知道有沒人有這方面的經驗  因為我雖然念過dft 知道k-s equation

  但要從零開始 自己寫出dft的code  忽然不知道麼下手

  我是想先計算最簡單的 氫原子系統就好

  就一個質子和一個電子

  能不能請各位給各方向  要怎麼下手??

  一些參考資料或建議也感謝不盡 ~~





--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 219.69.110.166
※ 編輯: peter308        來自: 219.69.110.166       (12/24 18:46)
→ waddler:http://tinyurl.com/37qplx5
推 uniserv1002:還以為是discrete fourier transform
推 robertshih:推樓上
→ MOONRAKER:看到FT只想得到fourier transform @@
→ Keelungman:你應該先弄清楚你要的exchange-correlation functional
→ peter308:基隆大 很清楚 LDA(Perdew-Zunger parameterisation )
→ peter308: 主要是怎麼弄成一個小程式 一些numerical 細節還在摸索
→ peter308:  以前用套裝軟體用習慣了 現在才發現根基很薄弱
→ peter308: density integrtion, 解K-Sequation 和 SCF算自洽解
→ peter308: 主要是這三部份.....我可能要先看看別人怎麼寫找靈感吧
→ peter308:但目前套裝軟體dft的subroutine都至少上百個 哈哈...
→ peter308: 可能要花些時間慢慢看吧....不過我想這過程應該無法逃避
→ peter308: 如果以後想做ab initio 方面的研究的話
推 Keelungman:先定好你的 basis 跟跑出 Hamiltonian 吧
→ Keelungman:我寫過的玩具, 主要是被處理晶體對稱的地方搞到頭大
→ Keelungman:如果你只跑單原子就比較不怕這問題
推 chendaolong:http://www.cpmd.org/manual.pdf Good!!!
→ pipidog:同意基隆男,我也覺得寫自洽不難,晶體的對稱性才頭大
→ pipidog:尤其是有磁性的系統,往往某些相的磁性周期性很複雜,簡直無
→ pipidog:從下手,完全不知道k該從哪邊走起
推 sunev:這不算小練習耶，如果是第一個程式不建議從這個開始。
→ pipidog:我覺得想玩自洽,其實從Hartree-Fock玩起就夠啦,最簡單的像
→ pipidog:純C組成的,例如graphene只要考慮一個orbital,高興的話還可
→ pipidog:以加個U進去看看,可以自洽做得很爽.
推 sunev:huckel --> hubbard --> PPP --> CNDO --> ZINDO
→ peter308: 感謝!!  我會試著寫看看
→ sneak: 以加個U進去看看,可以 https://noxiv.com
→ sneak:  以前用套裝軟體用習慣 https://daxiv.com
→ sneak: 如果你只跑單原子就比較 https://noxiv.com
→ sneak: 主要是怎麼弄成一個小程 https://muxiv.com
→ muxiv: 尤其是有磁性的系統,往 http://yofuk.com