24小时热门版块排行榜    

查看: 864  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见