24小时热门版块排行榜    

查看: 597  |  回复: 3
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

swallow319

木虫 (小有名气)

[交流] 求助ins文件damp命令的含义

求助ins文件damp命令的含义
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xi2004

至尊木虫 (职业作家)

DAMP  damp[0.7]  limse[15]

The DAMP parameters take different meanings for L.S. and CGLS refinements.  For L.S., damp is usually left at the default value unless there is severe correlation, e.g. when trying to refine a pseudo-centrosymmetric structure, or refining with few data per parameter (e.g. from powder data). A value in the range 1-10000 might then be appropriate.  The diagonal elements of the least-squares matrix are multiplied by (1+damp/1000) before inversion; this is a version of the Marquardt (1963) algorithm.  A side-effect of damping is that the standard deviations of poorly determined parameters will be artificially reduced; it is recommended that a final least-squares cycle be performed with little or no damping in order to improve these estimated standard deviations.  Theoretically, damping only serves to improve the convergence properties of the refinement, and can be gradually reduced as the refinement converges; it should not influence the final parameter values.  However in practice damping also deals effectively with rounding error problems in the (single-precision) least-squares matrix algebra, which can present problems when the number of parameters is large and/or restraints are used (especially when the latter have small esd's), and so it may not prove possible to lift the damping entirely even for a well converged refinement.

Note the use of 'DAMP 0 0' to estimate esds but not apply shifts, e.g. when a final L.S. 1 job is performed after CGLS refinement.

For CGLS refinements, damp is the multiplicative shift factor applied in the first cycle.  In subsequent CGLS cycles it is modified based on the experience in the previous cycles.  If a refinement proves unstable in the first cycle, damp should be reduced from its default value of 0.7.

If the maximum shift/esd for a L.S. refinement (excluding the overall scale factor) is greater than limse, all the shifts are scaled down by the same numerical factor so that the maximum is equal to limse.  If the maximum shift/esd is smaller than limse no action is taken.  This helps to prevent excessive shifts in the early stages of refinement.  limse is ignored in CGLS refinements.
xi2004-emuch@163.com
4楼2008-11-01 17:11:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

xi2004

至尊木虫 (职业作家)

解释在shelx的说明书中有.

就是有时会发散, 可以用damp 来限制一下, 叫它修得慢一点. 大致如此.
xi2004-emuch@163.com
2楼2008-11-01 16:45:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swallow319

木虫 (小有名气)

那么damp 0 0.001表示什么意思?
3楼2008-11-01 16:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见