| 查看: 1031 | 回复: 1 | ||||
zhangguangping木虫 (著名写手)
|
[交流]
【已解决】siesta-3.0-rc2中几何结构优化时DM的extrapolation问题 已有1人参与
|
|
我使用的是siesta-3.0-rc2,做几何结构优化，在输入文件对SCF和MD做了如下设置： # SCF options DM.Tolerance 1.d-4 # Tolerance in maximum difference between input and output DM MaxSCFIterations 300 # Maximum number of SCF iter DM.UseSaveDM T # to use continuation files DM.MixingWeight 0.05 # New DM amount for next SCF cycle DM.NumberPulay 5 DM.MixSCF1 F DM.PulayOnFile F # Store in memory ('F') or in files ('T') # NeglNonOverlapInt T # 'F'=do not neglect SolutionMethod Diagon # OrderN or Diagon ElectronicTemperature 300 K # Temp. for Fermi smearing # MD options MD.TypeOfRun CG # Type of dynamics: MD.VariableCell F MD.NumCGsteps 300 # Number of CG steps for coordinate optimization 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.ForceTolLength 0.02 eV/Ang #controls the convergence with respect to forces on Z-matrix lengths ZM.ForceTolAngle 0.00356549 Ry/rad #controls the convergence with respect to forces on Z-matrix angles ZM.MaxDisplLength 0.1 Ang # controls the maximum change in a Zmatrix length during an optimisation step ZM.MaxDisplAngle 0.003 rad # controls the maximum change in a Z-matrix angle during an optimisation step ZM.CalcAllForces T # Default value 但是在计算的时候，为什么每一个 CG move的时候关于新的DM出现如下的输出： New_DM. Step: 5 Re-using DM from previous geometry... Re-using DM without extrapolation (not allowed) Density Matrix sparsity pattern changed. 为什么出现： Re-using DM without extrapolation (not allowed) 我没有限制进行外推啊？ 怎么设置能够进行外推呢？ 我使用3.0-b的时候正常的外推，同样的设置 New_DM. Step: 85 Re-using DM from previous geometry... Extrapolating Density Matrix... Density Matrix sparsity pattern changed. 请教高手！ [ Last edited by zzy870720z on 2010-12-3 at 20:28 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
中青基了要发朋友圈吗？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
国社科申报系统有变化
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 谁做过有机硅树脂与SiO2溶胶的共混实验？两种物质混合能否制成互穿网络结构？
已经有6人回复
- 关于蓝宝石晶体生长过程中的散射中心以及镶嵌结构的消除
已经有20人回复
- 利用概率密度函数生成随机变量过程中出现超几何分布函数问题
已经有4人回复
- 请教一下那款免费软件的结构优化功能比较好用？
已经有7人回复
- 各位大侠，我可以优化结构，但只要一计算频率就出现下面的问题？
已经有4人回复
- 几何优化
已经有7人回复
- 今天计算发现了一些问题，几何成功就一定收敛了吗？
已经有6人回复
- 表面分子吸附的几何优化
已经有5人回复
- 气态结构的溶剂化，还是溶剂模型中优化？
已经有3人回复
- SIESTA几何优化不收敛是不是就说明了这个结构是不稳定的?
已经有4人回复
- 【求助】结构优化之后的静态自洽计算
已经有12人回复
- 【求助】利用侯博的方法优化结构出现的问题（求验证）
已经有11人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★
zzy870720z(金币+2):谢谢总结，请标明已解决，谢谢合作 2010-12-03 17:53:38
zzy870720z(金币+2):谢谢总结，请标明已解决，谢谢合作 2010-12-03 17:53:38
|
问题得到解决。虽然手册上关于DM.AllowReuse,DM.AllowExtrapolation 的设置默认都是true，但是在new_dm.F中写道： ! For "CG" calculations, the default is not to extrapolate the ! DM (unless requested by setting 'DM.AllowExtrapolation' to ! "true"  . The previous step's DM is reused.也就说当你的做CG方法的几何优化的时候，除非你明确写道 DM.AllowExtrapolation T 否则，对该方法这个选项为：F |
2楼2010-12-03 17:47:30
投票:
30