24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 890  |  回复: 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问题

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

弘德明志博学笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢总结,请标明已解决,谢谢合作 2010-12-03 17:53:38
引用回帖:
Originally posted by zhangguangping at 2010-11-03 04:14:47:
我使用的是siesta-3.0-rc2,做几何结构优化,在输入文件对SCF和MD做了如下设置:
# SCF options

DM.Tolerance          1.d-4         # Tolerance in maximum difference between input and output DM
MaxSC ...

问题得到解决。虽然手册上关于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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +26 111623 2026-04-04 29/1450 2026-04-08 00:10 by JourneyLucky
[考研] 266调剂 +5 daya sun 2026-04-07 6/300 2026-04-08 00:08 by daya sun
[考研] 材料工程322 +14 哈哈哈吼吼吼哈 2026-04-07 15/750 2026-04-07 22:51 by 来看流星雨10
[考研] 326求调剂 +5 9ahye 2026-04-02 6/300 2026-04-07 21:37 by lijunpoly
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 316求调剂 +7 yyx想调剂 2026-04-05 7/350 2026-04-07 14:31 by shdgaomin
[考研] 283求调剂 +18 A child 2026-04-04 18/900 2026-04-07 12:30 by 1018329917
[考研] 材料工程302分求调剂 +13 zyx上岸! 2026-04-04 13/650 2026-04-07 11:14 by 诗与自由
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 求调剂 +5 wos666 2026-04-03 5/250 2026-04-06 10:13 by 蓝云思雨
[考研] 一志愿华南师范361分,化学求调剂 +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 一志愿安徽大学计算机科学与技术学硕,331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
信息提示
请填处理意见