24小时热门版块排行榜    

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

vegalew

金虫 (小有名气)

[交流] 【求助】 CASTEP计算单原子能量不收敛问题(帖出了para,cell文件)

很多时候比如算原子在表面吸附能,或者内聚能的时候要计算单原子的能量

我计算单个Pt原子在真空的能量的时候把Pt原子放在边长为1nm的正方盒子中心计算的时候SCF算了1000多步还是不能收敛,Pt原子的能量在-700多eV左右震荡

求助为什么会SCF不收敛啊 一般30步之内都会收敛的。。。是不是模型有问题?

cell文件如下:
REMARK   Materials Studio PDB file
REMARK   Created:  Thu Jan 10 10:08:59 China Standard Time 2008
CRYST1   10.000   10.000   10.000  90.00  90.00  90.00 P1
ATOM      1 Pt   MOL     2       5.000   5.000   5.000  1.00  0.00          Pt   
TER


para文件如下:
comment : CASTEP calculation from Materials Studio
task : SinglePoint
xc_functional : PW91
spin_polarized : false
opt_strategy : Default
page_wvfns :        0
cut_off_energy :      280.000000000000000
grid_scale :        2.300000000000000
finite_basis_corr :        0
elec_energy_tol :   5.000000000000000e-007
max_scf_cycles :     9999
fix_occupancy : false
metals_method : dm
mixing_scheme : Pulay
mix_charge_amp :        0.500000000000000
mix_charge_gmax :        1.500000000000000
mix_history_length :       20
nextra_bands :        4
smearing_width :        0.010000000000000
calculate_stress : false
popn_calculate : false
calculate_densdiff : false
pdos_calculate_weights : false
num_dump_cycles : 0

体系是一个Pt原子在10A的立方盒子中心

[ Last edited by csfn on 2008-1-11 at 08:41 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csfn

荣誉版主 (知名作家)

优秀超版优秀版主优秀区长

引用回帖:
Originally posted by yhjiangshiyu at 2008-1-10 16:49:
我在1.5nm的格子里用dmol计算的时候基本上是交上去就完了啊,你的参数是不是有问题。但是我觉得计算出来的结果好像不对。

dmol算单个原确实不对~
何况你还加了个1.5nm的格子
7楼2008-01-10 16:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

vegalew

金虫 (小有名气)

还想问一下真空部分对计算复杂程度还有影响啊
2楼2008-01-10 16:25:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csfn

荣誉版主 (知名作家)

优秀超版优秀版主优秀区长

★ ★ ★ ★
小红豆(金币+4,VIP+0):感谢回答!
应该没设置对吧,否则不可能会出现这个情况
把你的parm文件和cell文件黏贴上来给大家看下吧~
3楼2008-01-10 16:25:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csfn

荣誉版主 (知名作家)

优秀超版优秀版主优秀区长

★ ★
spur(金币+2,VIP+0):感谢帮助!!
引用回帖:
Originally posted by vegalew at 2008-1-10 16:25:
还想问一下真空部分对计算复杂程度还有影响啊

有影响,不过对复杂程度不大~至于多少合适,一般来说文献都会说10A,大点点可以,小点就不大好。具体取多少,如果你不放心,可以做一个测试,当真空层大到一定的层次了对总能量影响应该是很小,也就说他会随着真空层增加出现一个总能收敛的现象
4楼2008-01-10 16:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 4/200 2026-08-18 18:31 by Niubile666
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +50 医学老男孩 2026-08-13 109/5450 2026-08-18 16:56 by 医学老男孩
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见