版块导航: 正在加载中...

客户端APP下载

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wtscrystal

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 214.2
散金: 250
红花: 3
帖子: 155
在线: 69.7小时
虫号: 252525
注册: 2006-05-20
专业: 无机合成和制备化学

[交流] 【讨论】ADF问题请教！

我用ADF优化结构
#! /bin/sh

# ==============================
# The Molecule
# ==============================

"$ADFBIN/adf" < ATOMS
1 xx    -1.770410000000    -3.066430000000    -4.964340000000 b=M51Vfinal
2 xx    -3.078970000000    -5.332940000000    -2.728710000000 b=M51Vfinal

378 281 285 1.0
END

CHARGE -21

UNRESTRICTED

BASIS
type DZ
core Large
createoutput None
END

XC
LDA SCF VWN
END

GEOMETRY
optim Delocalized
END

SAVE  TAPE21 TAPE13

SCF
iterations 500
END

NOPRINT LOGFILE

eor

＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃
最后出现以下结果：

<13:11:30> cart. step rms 0.00370285 0.00666667 T
Coordinates in Geometry Cycle 30
..................................

<13:23:04> 2 ErrMat 0.02700647 MaxEl 0.00139724

<13:23:32> 3 ErrMat 0.02097968 MaxEl 0.00081368

<13:24:01> 4 ErrMat 0.08667922 MaxEl -0.00441872

<13:24:27> 5 ErrMat 0.01938159 MaxEl -0.00076675

<13:24:50> 6 ErrMat 0.01986987 MaxEl 0.00155846

<13:25:11> 7 ErrMat 0.01017430 MaxEl -0.00099106

<13:25:27> 8 ErrMat 0.00705066 MaxEl 0.00070761

<13:25:42> 9 ErrMat 0.00354366 MaxEl -0.00021697

<13:25:58> 10 ErrMat 0.00172884 MaxEl -0.00006357

<13:26:13> 11 ErrMat 0.00078621 MaxEl -0.00003336

<13:26:15> SCF converged

<13:26:29> 12 ErrMat 0.00046677 MaxEl 0.00002342

<13:26:39> >>>> TOTEN

<13:27:26> >>>> POPAN

<13:27:34> >>>> DEBYE

<13:27:48> >>>> AMETS

<13:29:17> >>>> ENGRAD

<13:47:11> current energy -74.30022571 Hartree

<13:47:11> abs of energy change 0.00187059 0.00100000 F

<13:47:11> constrained gradient max 0.02527662 0.00100000 F

<13:47:11> constrained gradient rms 0.00217207 0.00066667 F

<13:47:11> gradient max 0.02527662

<13:47:11> gradient rms 0.00217207

<13:47:11> cart. step max 0.03342936 0.01000000 F

<13:47:11> cart. step rms 0.00842956 0.00666667 F

<13:47:11> ERROR: GEOMETRY DID NOT CONVERGE

<13:47:13> NORMAL TERMINATION

<13:47:23> END
Job v9 has finished

请问: SCF 收敛，但报错说： ERROR: GEOMETRY DID NOT CONVERGE
这是为何，是轮数不够，我用的30轮。
而发现结构中的一个二甲基胺发生了变形，如下图所示：

有何高招没有，各位仁兄！

» 猜你喜欢

1楼2010-03-28 10:06:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wtscrystal

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 214.2
散金: 250
红花: 3
帖子: 155
在线: 69.7小时
虫号: 252525
注册: 2006-05-20
专业: 无机合成和制备化学

引用回帖:

Originally posted by dukelr at 2010-04-01 13:59:48:
您的体系至少有285个原子 SCF需要的圈数多应该算是正常了吧

此君为正解，确实不少于285个原子，但问题是它总不收敛，是不是我的基组用的不对：
我里面有Mo, V W, C,N,O,h等，我用的是DZ，有人讲基组不对，现改为:TZ2P, spin polarization=2, 不知行不？

15楼2010-04-01 16:28:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 15 个回答

qphll

金虫 (正式写手)

应助: 18 (小学生)
金币: 2708.7
散金: 2294
红花: 8
帖子: 554
在线: 227.1小时
虫号: 89654
注册: 2005-08-29
性别: GG
专业: 碳素材料与超硬材料

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1):谢谢提示 2010-03-31 09:17

显然是在你的30 Geometry Cycle后,没有满足收敛条件. 或者增多步数, 或者减弱收敛要求.

赞一下(2人)

Life, Love, Laugh.

2楼2010-03-29 08:46:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wtscrystal

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 214.2
散金: 250
红花: 3
帖子: 155
在线: 69.7小时
虫号: 252525
注册: 2006-05-20
专业: 无机合成和制备化学

也就是说只是时间不够，或步数不够并不是因为别的什么原因？？？

3楼2010-03-29 12:00:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

random789

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 454.6
散金: 115
帖子: 523
在线: 20.8小时
虫号: 945039
注册: 2010-01-19
专业: 理论和计算化学

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+1):谢谢提示 2010-03-31 09:17

你用最后一圈不收敛的集合结构代入输入文件中重算，顺便把几何优化的圈数扩大

赞一下(2人)

好好学习，天天向上！

4楼2010-03-29 14:49:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 15 个回答