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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 计算石墨烯不收敛的问题
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 2446  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fubao-411

木虫 (小有名气)

[求助] 计算石墨烯不收敛的问题

Sample Text    我用ATK 计算石墨烯的透射谱，电极选用一个单胞时，计算是收敛的，后来我又算了下电极选用2个单胞时的情况，计算不收敛了，运行了4天后提示“exceeded the maximum number of the iteration”。理论上，电极是越长越好，但这里为什么会不收敛呢
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
轻松走一回！
1楼 2011-08-01 10:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

【答案】应助回帖

★ ★
franch(金币+2): 谢谢回帖交流 2011-08-01 21:28:46
fubao-411(金币+3): 谢谢，今天补上 2011-08-03 16:16:37
是不是自洽步数设置的太小了?增大一点试一试。
赞一下(1人) 回复此楼
2楼2011-08-01 12:12:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fubao-411

木虫 (小有名气)
我也想过可能是这个问题，自洽默认是100步，建模时我都是精确到了小数点后4位，并且是在ATK内优化完成的，1单胞电极没出现这个问题，2单胞电极就出现这个问题，
赞一下 回复此楼
轻松走一回！
3楼2011-08-02 10:41:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fubao-411

木虫 (小有名气)
按下“回帖评分”没反应，草的......
赞一下 回复此楼
轻松走一回！
4楼2011-08-02 10:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

太白击剑

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
fubao-411(金币+7): 谢谢，我再试试 2011-08-03 16:17:19
gzqdyouxia(金币+2): 谢谢指导 2011-08-04 20:20:19
还有可能是k点取少了，可以增大一点试试；另外一种方法是增加温度，先在在较高的温度下计算一个.nc文件，然后再将温度调到你所需要的值（一般是300k），再以这个文件为初始电荷分布计算体系的.nc文件
赞一下(1人) 回复此楼
5楼2011-08-02 17:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohuazhon

木虫 (正式写手)
★ ★
gzqdyouxia(金币+2): 谢谢指导 2011-08-04 20:19:57
为了节省计算时间，加快收敛速度，可以尝试以下设置：
• Basis set parameters: Type = SingleZetaPolarized (saves time and memory, without compromising the accuracy substantially for gold)
• Brillouin zone integration parameters: Number of k-points = 5 x 5 x 100
• Electron density parameters: Mesh cut-off = 200 Ry
• Eigenstate occupation parameters: Electron temperature = 1000 K
• Iteration mixing parameters: Diagonal mixing parameter =0.05
楼上说的很有道理。
赞一下(1人) 回复此楼
6楼2011-08-04 19:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohuazhon

木虫 (正式写手)
★ ★
gzqdyouxia(金币+2): 谢谢提示 2011-08-04 20:19:35
你还可以用下面的脚本来检验K点设置多少合适，计算时N取值到多少看你自己了，这里只取到10
from ATK.TwoProbe import *
from ATK.MPI import processIsMaster
checkpoint_file = "×.nc"
scf = restoreSelfConsistentCalculation(checkpoint_file)
for N in range(1,10):
kpoints = brillouinZoneIntegrationParameters((N, N))
transmission_spectrum = calculateTransmissionSpectrum(
self_consistent_calculation = scf,
energies = (0.0,)*electronVolt,
brillouin_zone_integration_parameters = kpoints, )
if processIsMaster():
print N, transmission_spectrum.coefficients()[0]
赞一下(1人) 回复此楼
7楼2011-08-04 20:00:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fubao-411

木虫 (小有名气)
多谢楼上几位  
回复此楼
轻松走一回！
8楼2011-08-06 11:05:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaypu

铁杆木虫 (小有名气)
★ ★
franch(金币+2): 谢谢回帖交流 2011-08-30 11:05:06
引用回帖:
1楼: Originally posted by fubao-411 at 2011-08-01 10:17:00:
Sample Text    我用ATK 计算石墨烯的透射谱，电极选用一个单胞时，计算是收敛的，后来我又算了下电极选用2个单胞时的情况，计算不收敛了，运行了4天后提示“exceeded the maximum number of t ...

交换关联 可以用LDA 我记得有篇文章提到过这个问题LDA和GGA计算石墨烯的时候差别不大，不知道有没有记错，但是 LDA比较容易收敛，你可以试一下
赞一下(1人) 回复此楼
9楼2011-08-28 10:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fubao-411

木虫 (小有名气)
引用回帖:
9楼: Originally posted by kaypu at 2011-08-28 10:25:25:
交换关联 可以用LDA 我记得有篇文章提到过这个问题LDA和GGA计算石墨烯的时候差别不大，不知道有没有记错，但是 LDA比较容易收敛，你可以试一下

多谢
回复此楼
轻松走一回！
10楼2011-08-30 08:38:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fubao-411 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 安徽大学专硕生物与医药专业(086000)324分，英语已过四六级，六级521，求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 0856调剂 +5 求求让我有书读� 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 一志愿武汉理工，总分321，英一数二，求老师收留。 +5 nnnnnnn5 2026-03-25 5/250 2026-03-27 04:42 by wxiongid
[考研] 317求调剂 +7 蛋黄咸肉粽 2026-03-26 7/350 2026-03-27 02:29 by fmesaito
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 求调剂 +5 芦lty 2026-03-25 6/300 2026-03-26 20:49 by 不吃魚的貓
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 286求调剂 +13 Faune 2026-03-21 13/650 2026-03-26 19:52 by peike
[考研] 生物学学硕，一志愿湖南大学，初试成绩338 +4 YYYYYNNNNN 2026-03-26 4/200 2026-03-26 19:00 by macy2011
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
信息提示
关闭
请填处理意见
关闭