版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

wdaobin

捐助贵宾 (正式写手)

应助: 1 (幼儿园)
金币: 143.5
红花: 1
帖子: 552
在线: 333.4小时
虫号: 538494
注册: 2008-04-02
专业: 传输与交换光子学

[求助] Openmx中Graphene的优化问题！已有2人参与

小虾最近在学习Openmx，想用它来优化Graphene，元胞中使用了两个C原子，但是一直得不到收敛的结果，程序一直运行到最大循环数然后被强制终止，请大侠帮忙看看问题出在哪了？谢谢！下面是我的输入文件：
# 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)

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有3人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2014-03-16 21:11:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢指导 2014-03-17 10:32:06

我覺得你直接去 OPENMX 官方發問會比較快
http://www.openmx-square.org/forum/patio.cgi
小木蟲上似乎很很很少人用 OPENMX
OPENMX.... 不好上手

赞一下(1人)

回复此楼

2楼2014-03-17 03:54:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

1ST强帖: 1
应助: 61 (初中生)
贵宾: 0.33
金币: 5959.3
散金: 1
红花: 21
帖子: 1884
在线: 491.1小时
虫号: 142265
注册: 2005-12-21
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wdaobin

捐助贵宾 (正式写手)

应助: 1 (幼儿园)
金币: 143.5
红花: 1
帖子: 552
在线: 333.4小时
虫号: 538494
注册: 2008-04-02
专业: 传输与交换光子学

谢谢y1ding，问题已解决，不过我发现在3.5上仍然不收敛，在3.7上运行到MD=13时计算正常结束了。以后还请多多指教！

赞一下

回复此楼

4楼2014-03-19 21:12:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 wdaobin 的主题更新

返回列表