| 查看: 1868 | 回复: 3 | |||
[求助]
利用siesta做verlet MD时出现问题siesta: WARNING: Atoms
|
|
在利用siesta做verlet MD，用Si8做的时候顺利结束，但是换成LiBH4或者NH3BH3时，均出现siesta: WARNING: Atoms 1880 2366 too close: rij = 0.288938 Ang之类的错误警告，不知道是怎么回事？ 现将输入文件附上： SystemName libh4 SystemLabel libh4 NumberOfAtoms 24 NumberOfSpecies 3 %block ChemicalSpeciesLabel 1 1 H 2 3 Li 3 5 B %endblock ChemicalSpeciesLabel PAO.BasisSize SZ PAO.EnergyShift 300 meV LatticeConstant 6.4494 Ang %block LatticeVectors 1.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 0.822867243000000 0.000000000000000 0.000000000000000 0.000000000000000 0.820525940000000 %endblock LatticeVectors MeshCutoff 130.0 Ry %block kgrid_Monkhorst_Pack 6 0 0 0.0 0 6 0 0.0 0 0 6 0.0 %endblock kgrid_Monkhorst_Pack MaxSCFIterations 50 DM.MixingWeight 0.3 DM.NumberPulay 4 DM.Tolerance 1.d-4 DM.UseSaveDM SolutionMethod diagon ElectronicTemperature 25 meV MD.TypeOfRun Verlet MD.InitialTimeStep 1 MD.FinalTimeStep 200 MD.LengthTimeStep 3.0 fs MD.InitialTemperature 10.0 K WriteMDHistory T AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.0949999988079071 0.3339999914169312 0.0740000009536743 1 h 0.0949999988079071 0.8339999914169313 0.5740000009536742 1 h -0.0949999988079071 -0.3339999914169312 0.0740000009536743 1 h -0.0949999988079071 0.1660000085830689 0.5740000009536743 1 h 0.5949999988079071 -0.3339999914169312 0.0740000009536743 1 h 0.5949999988079071 0.1660000085830688 0.5740000009536743 1 h 0.4050000011920930 0.3339999914169312 0.0740000009536743 1 h 0.4050000011920930 0.8339999914169310 0.5740000009536743 1 h 0.2500000000000000 0.0496999993920326 0.1831000000238419 1 h 0.2500000000000000 0.3510000109672546 0.3802999854087830 1 h 0.2500000000000000 0.5496999993920326 0.6831000000238418 1 h 0.2500000000000000 0.8510000109672545 0.8802999854087830 1 h -0.2500000000000000 -0.0496999993920326 0.1831000000238419 1 h -0.2500000000000000 -0.3510000109672546 0.3802999854087830 1 h -0.2500000000000000 0.4503000006079674 0.6831000000238419 1 h -0.2500000000000000 0.1489999890327454 0.8802999854087829 1 h 0.0000000000000000 0.0000000000000000 0.0000000000000000 2 li 0.0000000000000000 0.5000000000000000 0.4999999999999999 2 li 0.5000000000000000 -0.0000000000000000 -0.0000000000000000 2 li 0.5000000000000000 0.4999999999999999 0.5000000000000000 2 li 0.2500000000000000 0.2712000012397766 0.1729999929666519 3 b 0.2500000000000000 0.7712000012397767 0.6729999929666518 3 b -0.2500000000000000 -0.2712000012397766 0.1729999929666519 3 b -0.2500000000000000 0.2287999987602234 0.6729999929666519 3 b %endblock AtomicCoordinatesAndAtomicSpecies 如需要输出文件，请留言告知 谢谢 |
回复此楼» 猜你喜欢
重要消息，中午系统在维护
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有101人回复
-
时间戳又变了8-15
已经有26人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- siesta 计算电子占据数的问题
已经有3人回复
- 用siesta画声子谱的问题
已经有12人回复
- 原来编译的siesta突然不能用，重新编译没错，运行出错
已经有13人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】siesta 做NOSE动力学出错，请高手指教
已经有11人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】siesta多个/种分子优化的问题
已经有7人回复
- 【求助】siesta3.0 安装，错误一箩筐
已经有11人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】关于siesta的几何结构优化中断后继续的问题【完结】
已经有6人回复
- 【求助】siesta使用B3lyp但是没有输出文件，所以找不出错误
已经有6人回复
hw45888792
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 5 (幼儿园)
- 金币: 848.5
- 散金: 45
- 红花: 3
- 帖子: 339
- 在线: 282.9小时
- 虫号: 711113
- 注册: 2009-02-28
- 性别: GG
- 专业: 理论和计算化学
2楼2011-06-28 10:11:12
3楼2011-06-28 10:20:33
hw45888792
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 5 (幼儿园)
- 金币: 848.5
- 散金: 45
- 红花: 3
- 帖子: 339
- 在线: 282.9小时
- 虫号: 711113
- 注册: 2009-02-28
- 性别: GG
- 专业: 理论和计算化学
4楼2011-06-28 10:27:50
