| 查看: 2228 | 回复: 10 | ||
fubao-411木虫 (小有名气)
|
[求助]
计算石墨烯不收敛的问题
|
| Sample Text 我用ATK 计算石墨烯的透射谱，电极选用一个单胞时，计算是收敛的，后来我又算了下电极选用2个单胞时的情况，计算不收敛了，运行了4天后提示“exceeded the maximum number of the iteration”。理论上，电极是越长越好，但这里为什么会不收敛呢 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 氧化石墨烯的问题
已经有39人回复
- 为何用晶格常数的实验值计算声子？
已经有18人回复
- 请教关于石墨烯XRD的问题
已经有26人回复
- 石墨烯干燥问题
已经有35人回复
- 石墨烯转变温度的问题
已经有13人回复
- 求助有关石墨烯的相关问题
已经有5人回复
- 【求助】石墨烯吸附气体分子的优化问题
已经有23人回复
- 【求助】准确的建石墨烯结构的方法
已经有26人回复
- 【求助成功】石墨烯能带计算的问题
已经有4人回复
- 【请教】关于石墨烯合成的几个问题
已经有22人回复
- 【讨论】关于石墨烯的问题
已经有25人回复
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3679.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 915.2小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
2楼2011-08-01 12:12:08
fubao-411
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5841.5
- 帖子: 228
- 在线: 82.7小时
- 虫号: 733126
- 注册: 2009-03-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-08-02 10:41:16
fubao-411
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5841.5
- 帖子: 228
- 在线: 82.7小时
- 虫号: 733126
- 注册: 2009-03-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2011-08-02 10:47:25
太白击剑
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 6482.2
- 散金: 10
- 红花: 1
- 帖子: 372
- 在线: 349.9小时
- 虫号: 970413
- 注册: 2010-03-14
- 专业: 凝聚态物性 II ：电子结构
5楼2011-08-02 17:26:51
guohuazhon
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 12479.2
- 散金: 10
- 红花: 2
- 帖子: 931
- 在线: 92.6小时
- 虫号: 807030
- 注册: 2009-07-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
gzqdyouxia(金币+2): 谢谢指导 2011-08-04 20:19:57
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 楼上说的很有道理。 |
6楼2011-08-04 19:56:52
guohuazhon
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 12479.2
- 散金: 10
- 红花: 2
- 帖子: 931
- 在线: 92.6小时
- 虫号: 807030
- 注册: 2009-07-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
gzqdyouxia(金币+2): 谢谢提示 2011-08-04 20:19:35
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] |
7楼2011-08-04 20:00:56
fubao-411
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5841.5
- 帖子: 228
- 在线: 82.7小时
- 虫号: 733126
- 注册: 2009-03-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2011-08-06 11:05:25
kaypu 
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 399
- 散金: 130
- 帖子: 292
- 在线: 113.4小时
- 虫号: 864391
- 注册: 2009-10-07
- 专业: 凝聚态物性 II ：电子结构
9楼2011-08-28 10:25:25
fubao-411
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 5841.5
- 帖子: 228
- 在线: 82.7小时
- 虫号: 733126
- 注册: 2009-03-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2011-08-30 08:38:34