 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1451 | 回復: 3 | ||
[求助]
Openmx中Graphene的優(yōu)化問題！ 已有2人參與
|
|
小蝦最近在學習Openmx，想用它來優(yōu)化Graphene，元胞中使用了兩個C原子，但是一直得不到收斂的結果，程序一直運行到最大循環(huán)數(shù)然后被強制終止，請大俠幫忙看看問題出在哪了？謝謝！下面是我的輸入文件： # Graphene written by XXX System.CurrrentDirectory /home/wangdb/Workspaces/Openmx/work-wdb/ # default=./ System.Name graphene DATA.PATH /home/wangdb/Workspaces/Openmx/DFT_DATA06 # default=../DFT_DATA11 level.of.stdout 1 # default=1 (1-3) level.of.fileout 1 # default=1 (1-3) # Definition of Atomic Species Species.Number 1 <Definition.of.Atomic.Species C C6.0-s2p2d1 C_PBE Definition.of.Atomic.Species> # Atoms Atoms.Number 2 Atoms.SpeciesAndCoordinates.Unit FRAC # Ang|AU <Atoms.SpeciesAndCoordinates # Unit=Ang. 1 C 0.000000000000 0.000000000000 0.000000000000 2.0 2.0 2 C 0.333333333333 -0.333333333333 0.000000000000 2.0 2.0 Atoms.SpeciesAndCoordinates> Atoms.UnitVectors.Unit AU # Ang|AU <Atoms.UnitVectors # unit=Ang. 2.3045392835 3.99157912706038 0.000000000000 2.3045392835 -3.99157912706038 0.000000000000 0.000000000000 0.000000000000 20.9080342511016 Atoms.UnitVectors> # SCF or Electronic System scf.XcType LDA # LDA|LSDA-CA|LSDA-PW scf.SpinPolarization off # On|Off scf.ElectronicTemperature 300.0 # default=300 (K) scf.energycutoff 700.0 # default=150 (Ry) scf.maxIter 100 # default=40 scf.EigenvalueSolver band # Recursion|Cluster|Band scf.lapack.dste dstevx # dstegr|dstedc|dstevx, default=dstegr scf.Kgrid 24 24 1 # means 4x4x4 scf.Mixing.Type rmm-diisk # Simple|Rmm-Diis|Gr-Pulay scf.Init.Mixing.Weight 0.010 # default=0.30 scf.Min.Mixing.Weight 0.001 # default=0.001 scf.Max.Mixing.Weight 0.200 # default=0.40 scf.Mixing.History 15 # default=5 scf.Mixing.StartPulay 5 # default=6 scf.criterion 1.0e-6 # default=1.0e-6 (Hartree) # 1D FFT 1DFFT.NumGridK 900 # default=900 1DFFT.NumGridR 900 # default=900 1DFFT.EnergyCutoff 3600.0 # default=3600 (Ry) # SCF Order-N orderN.HoppingRanges 8.389 # default=5.0 (Ang) orderN.NumHoppings 2 # default=2 orderN.KrylovH.order 800 # default=400 # MD or Geometry Optimization MD.Type opt # Nomd|Opt|DIIS|NVE|NVT_VS|NVT_NH MD.Opt.DIIS.History 7 # default=7 MD.Opt.StartDIIS 5 # default=5 MD.maxIter 100 # default=1 MD.TimeStep 1.0 # default=0.5 (fs) MD.Opt.criterion 1.0e-7 # default=1.0e-4 (Hartree/bohr) |
回復此樓» 猜你喜歡
求調劑院校信息
已經有4人回復
-
085600材料與化工306
已經有4人回復
-
286求調劑
已經有10人回復
-
328求調劑，英語六級551，有科研經歷
已經有9人回復
-
一志愿北京化工大學070300 學碩336求調劑
已經有4人回復
-
286分人工智能專業(yè)請求調劑愿意跨考！
已經有8人回復
-
資源與環(huán)境 調劑申請(333分)
已經有5人回復
-
280求調劑
已經有12人回復
-
269專碩求調劑
已經有5人回復
-
材料學碩301分求調劑
已經有7人回復
【答案】應助回帖
★ ★
感謝參與，應助指數(shù) +1
fzx2008: 金幣+2, 謝謝指導 2014-03-17 10:32:06
感謝參與，應助指數(shù) +1
fzx2008: 金幣+2, 謝謝指導 2014-03-17 10:32:06
|
我覺得你直接去 OPENMX 官方發(fā)問會比較快 http://www.openmx-square.org/forum/patio.cgi 小木蟲上似乎很很很少人用 OPENMX OPENMX.... 不好上手 |
2樓2014-03-17 03:54:12
y1ding
鐵桿木蟲 (著名寫手)
- 1ST強帖: 1
- 應助: 61 (初中生)
- 貴賓: 0.33
- 金幣: 5959.3
- 散金: 1
- 紅花: 21
- 帖子: 1884
- 在線: 491.1小時
- 蟲號: 142265
- 注冊: 2005-12-21
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
wdaobin: 金幣+10, ★★★很有幫助, 謝謝！ 2014-03-19 21:06:53
感謝參與，應助指數(shù) +1
wdaobin: 金幣+10, ★★★很有幫助, 謝謝！ 2014-03-19 21:06:53
|
scf.energycutoff 200.0 # default=150 (Ry) <Definition.of.Atomic.Species C C6.0-s2p2d1 C_CA13 Definition.of.Atomic.Species> 改動2處，可以收斂 ******************************************************* Total Energy (Hartree) at MD = 4 ******************************************************* Uele = -3.916208267222 Ukin = 8.305197749056 UH0 = -115.945584489942 UH1 = 0.030973892655 Una = -5.898807821628 Unl = -1.050688172920 Uxc0 = -1.927500964985 Uxc1 = -1.927500964985 Ucore = 106.751941650432 Uhub = 0.000000000000 Ucs = 0.000000000000 Uzs = 0.000000000000 Uzo = 0.000000000000 Uef = 0.000000000000 UvdW = 0.000000000000 Utot = -11.661969122317 Note: Utot = Ukin+UH0+UH1+Una+Unl+Uxc0+Uxc1+Ucore+Uhub+Ucs+Uzs+Uzo+Uef+UvdW Uele: band energy Ukin: kinetic energy UH0: electric part of screened Coulomb energy UH1: difference electron-electron Coulomb energy Una: neutral atom potential energy Unl: non-local potential energy Uxc0: exchange-correlation energy for alpha spin Uxc1: exchange-correlation energy for beta spin Ucore: core-core Coulomb energy Uhub: LDA+U energy Ucs: constraint energy for spin orientation Uzs: Zeeman term for spin magnetic moment Uzo: Zeeman term for orbital magnetic moment Uef: electric energy by electric field UvdW: semi-empirical vdW energy (see also PRB 72, 045121(2005) for the energy contributions) ******************************************************* Computational times (s) at MD = 4 ******************************************************* DFT in total = 12.12758 Set_OLP_Kin = 1.44986 Set_Nonlocal = 2.16248 Set_ProExpn_VNA = 2.26113 Set_Hamiltonian = 1.70408 Poisson = 0.00246 diagonalization = 0.50167 Mixing_DM = 0.01076 Force = 1.85790 Total_Energy = 0.92433 Set_Aden_Grid = 0.02858 Set_Orbitals_Grid = 0.27365 Set_Density_Grid = 0.89989 RestartFileDFT = 0.02197 Mulliken_Charge = 0.00213 FFT(2D)_Density = 0.00316 ******************************************************* MD or geometry opt. at MD = 4 ******************************************************* <Steepest_Descent> SD_scaling= 0.944854433488 <Steepest_Descent> |Maximum force| (Hartree/Bohr) = 0.000000022448 <Steepest_Descent> Criterion (Hartree/Bohr) = 0.000000100000 atom= 1, XYZ(ang) Fxyz(a.u.)= -0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 atom= 2, XYZ(ang) Fxyz(a.u.)= 1.2195 -0.7041 -0.0000 0.0000 -0.0000 -0.0000 outputting data on grids to files... *********************************************************** *********************************************************** Computational Time (second) *********************************************************** *********************************************************** Min_ID Min_Time Max_ID Max_Time Total Computational Time = 4 67.461 0 68.100 readfile = 7 6.883 0 6.886 truncation = 5 1.936 1 2.428 MD_pac = 0 0.028 5 0.036 OutData = 1 1.956 0 2.593 DFT = 1 56.048 5 56.540 *** In DFT *** Set_OLP_Kin = 5 1.045 0 5.790 Set_Nonlocal = 0 8.654 5 13.400 Set_ProExpn_VNA = 5 6.023 0 8.983 Set_Hamiltonian = 1 11.370 6 11.370 Poisson = 7 0.014 1 0.015 Diagonalization = 0 3.395 1 3.395 Mixing_DM = 7 0.070 0 0.070 Force = 1 7.363 6 7.363 Total_Energy = 6 3.596 0 3.693 Set_Aden_Grid = 0 0.116 5 3.067 Set_Orbitals_Grid = 3 0.000 0 1.102 Set_Density_Grid = 4 5.358 0 5.364 RestartFileDFT = 5 0.066 7 0.895 Mulliken_Charge = 1 0.009 7 0.010 FFT(2D)_Density = 7 0.068 3 0.068 Others = 1 0.091 5 1.694 The calculation was normally finished. |
3樓2014-03-17 16:20:25
4樓2014-03-19 21:12:31
相關版塊跳轉
我要訂閱樓主 wdaobin 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學調劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 生物學調劑 +3 | Surekei 2026-03-21 | 3/150 |
|
|
|
[考研] 0703化學297求調劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
|
[考研] 085601調劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 一志愿西南交大，求調劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 296求調劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 22408 344分 求調劑 一志愿 華電計算機技術 +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南昌大學，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 321求調劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 求調劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 290求調劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 326求調劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|