| 查看: 956 | 回复: 2 | ||
lanlanwa9002铁虫 (小有名气)
|
[求助]
（求大神）SIESTA %block Zmatrix 疑问
|
|
我用SIESTA 计算表面的吸附性质，由于需要固定键长，所以原子坐标区块选择了内坐标形式 ：%block Zmatrix。在计算中却发现计算常常在并未到达受力收敛标准0.01 eV/Ang，的情况下就正常结束了。即便还在算，也总是达不到这个标准。想这是怎么回事呢？于是我放开对键长的限制，可还是达不到受力标准，可见并不像是约束导致的。然后我换了一个苯分子在某金属表面吸附的简单模型，发现依旧是受力不收敛。于是我想请教各位大牛：请问是不是使用了%block Zmatrix之后，在自洽场计算部分或结构优化部分还需要做什么特殊的设置啊？ 附（.fdf文件）： SystemName M SystemLabel 1 NumberOfAtoms 84 NumberOfSpecies 8 %block ChemicalSpeciesLabel 1 13 Al 2 14 Si 3 8 O 5 12 Mg 4 1 H 6 11 Na 7 8 O_surf 8 1 H_surf %endblock ChemicalSpeciesLabel PAO.BasisSize DZP PAO.EnergyShift 320 meV PAO.SplitNorm 0.23 LatticeConstant 1 Ang %block LatticeVectors 10.360 0.000 0.000 0.000 8.98 0.000 0.000 0.000 40.00 %endblock LatticeVectors %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 1 0.0 %endblock kgrid_Monkhorst_Pack XC.functional GGA # Exchange-correlation functional XC.authors PBE # Exchange-correlation version SpinPolarized F MeshCutoff 220.0 Ry MaxSCFIterations 300 DM.NumberPulay 3 DM.Tolerance 1.d-4 SolutionMethod diagon DM.MixingWeight 0.25 # New DM amount for next SCF cycle #ElectronicTemperature 25 meV # MD options MD.TypeOfRun CG MD.NumCGsteps 1000 ZM.UnitsLength Ang #the units of length used during Z-matrix input ZM.UnitsAngle deg #the units of angles used during Z-matrix input ZM.MaxDisplLength 0.1 Ang # controls the maximum change in a Zmatrix length during an optimisation step ZM.ForceTolLength 0.01 eV/Ang #controls the convergence with respect to forces on Z-matrix lengths ZM.ForceTolAngle 0.00356549 Ry/rad ZM.MaxDisplAngle 0.003 rad MD.VariableCell F AtomicCoordinatesFormat ScaledCartesian %block Zmatrix molecule 6 0 0 0 3.642743450 4.585977690 7.416304920 0 0 0 7 1 0 0 2.397 38.321 -122.002 0 1 1 8 1 2 0 3.135 13.167 58.378 1 1 1 8 1 2 3 2.651 21.66 115.885 1 1 1 cartesian 1 4.574144840 3.055682480 3.301943830 0 0 0 1 9.748525910 3.071737580 3.296439130 0 0 0 5 4.575629600 5.991863520 3.323686150 0 0 0 .................... 4 0.295759990 0.046124920 4.404757640 0 0 0 4 5.508497320 0.051660940 4.391599880 0 0 0 %endblock Zmatrix # Output options WriteCoorXmol T WriteCoorCerius True WriteCoorStep True #------------------------------------------------------------------------------- #Restart DM.UseSaveDM F MD.UseSaveXV F MD.UseSaveCG F |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
siesta基组优化问题
已经有7人回复- 请教如何加快siesta的收敛速度
已经有9人回复
- siesta中*.xyz和*.XV文件
已经有5人回复
- siesta运行出现问题，Error in Cholesky factorisation in cdiag，原子间距过小，求助
已经有5人回复
- siesta结构优化收敛标准
已经有5人回复
- 请大家帮忙看看siesta输入文件有哪些错误
已经有8人回复
- 求助用siesta计算COHP
已经有16人回复
- SIESTA计算能级不准怎么解决
已经有10人回复
- siesta计算出现forrtl: error (69): process interrupted (SIGINT)
已经有10人回复
- 请教，为什么会出现这样的情况？谢谢大家
已经有6人回复
- transiesta 计算错误
已经有15人回复
- 请问siesta计算BornCharge时为什么输出的张量不满足对称性要求呢？谢谢
已经有15人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- 【分享】SIESTA中Bader analysis
已经有43人回复
- 【求助】如何在siesta中设置ghost state原子的计算
已经有16人回复
- 【求助】H2 Binding Energy@SIESTA
已经有13人回复
- 【解决】讨论，siesta中作HOMO 与LUMO图的方法
已经有12人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】siesta输入文件中赝势文件的作用？
已经有8人回复
- 【求助】请教siesta3.0-beta如何计算PDOS、overlap population及coop
已经有12人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
lanlanwa9002
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2083.9
- 红花: 16
- 帖子: 253
- 在线: 80小时
- 虫号: 1754784
- 注册: 2012-04-14
- 专业: 理论和计算化学
|
下面是很让我无语的每一步的最大原子受力： Max 0.488463 constrained Max 0.288270 constrained Max 0.289708 constrained Max 0.279172 constrained Max 0.266928 constrained Max 0.252177 constrained Max 0.237016 constrained Max 0.221676 constrained Max 0.207109 constrained Max 0.191344 constrained Max 0.182107 constrained Max 0.182359 constrained Max 0.182625 constrained Max 0.182803 constrained Max 0.183053 constrained Max 0.183335 constrained Max 0.183581 constrained Max 0.183842 constrained Max 0.187210 constrained Max 0.197474 constrained Max 0.206927 constrained Max 0.216649 constrained Max 0.225810 constrained Max 0.235189 constrained Max 0.244180 constrained Max 0.252429 constrained Max 0.260475 constrained Max 0.268722 constrained Max 0.276296 constrained Max 0.283475 constrained Max 0.290631 constrained Max 0.298206 constrained Max 0.306024 constrained Max 0.313518 constrained Max 0.321264 constrained Max 0.329415 constrained Max 0.325507 constrained Max 0.227337 constrained Max 0.191721 constrained Max 0.229066 constrained Max 0.192127 constrained Max 0.337578 constrained Max 0.191278 constrained Max 0.216007 constrained Max 0.191841 constrained Max 0.191755 constrained Max 0.194218 constrained Max 0.197301 constrained Max 0.202302 constrained Max 0.199045 constrained Max 0.223329 constrained Max 0.190863 constrained Max 0.209688 constrained Max 0.256961 constrained Max 0.215010 constrained Max 0.216190 constrained Max 0.217720 constrained Max 0.220328 constrained Max 0.224456 constrained Max 0.222023 constrained Max 0.190453 constrained Max 0.240154 constrained Max 0.189719 constrained Max 0.295657 constrained Max 0.209734 constrained Max 0.255772 constrained Max 0.232474 constrained Max 0.221453 constrained Max 0.204881 constrained Max 0.190059 constrained Max 0.189692 constrained Max 0.234808 constrained Max 0.210606 constrained Max 0.252750 constrained Max 0.214895 constrained Max 0.212159 constrained Max 0.208620 constrained Max 0.212499 constrained Max 0.198380 constrained Max 0.194537 constrained Max 0.200486 constrained Max 0.228371 constrained Max 0.212400 constrained Max 0.222177 constrained Max 0.221101 constrained |
2楼2014-11-30 18:15:08
lanlanwa9002
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2083.9
- 红花: 16
- 帖子: 253
- 在线: 80小时
- 虫号: 1754784
- 注册: 2012-04-14
- 专业: 理论和计算化学
3楼2014-12-11 11:41:10