24小时热门版块排行榜    

查看: 882  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +11 生如夏花… 2026-03-11 13/650 2026-03-17 09:02 by 雾散后相遇lc
[考研] 考研化学学硕调剂,一志愿985 +3 张vvvv 2026-03-15 5/250 2026-03-16 20:25 by 张vvvv
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
请填处理意见