| 查看: 2987 | 回复: 11 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
关于建模结构优化的一些问题
|
|||||
|
新手發問 请问一下我想做有特定形状的缺陷的2x2x2的 周期性 的结构，材料是使用硅，要计算态密度，声子谱，能带结构 有试着先做一个的3x3x3超晶胞10颗缺陷的结构，在CPU i5- 2400 内存为8G的电脑上面跑不过跑了一天的工作就失败了，看来电脑不够力 结构优化第一次点选ultralfine，180 eV的k点4x4x4所做出来的晶格常数为5.39 更改不同的截止能量后要看能量收敛的程度 看到文献都会用 截止能量 跟 能量 的关系 那么是要看 final energy 吗还是迭代之后的 total energy 呢 有试着用硅晶胞边常设为5.4309然后用原始晶胞去最佳化分别选用LDA跟GGA- PBE的截止能量各从100～1000 k点的1x1x1 赝势为ultrasoft 其他都是预设值，晶体越大所需k点数越少 看到一些教学提到硅材料截止能量不需要太高大概在200 eV的的截止能量就能得到准确的数值？ 而得出的总能量我是看最后一次出现的final energy 而总能量LDA的大约都在-102.xx eV的左右GGA大约在-102.90xx LDA的会低估晶格常数GGA会高估晶格常数    回归正题想问一下必须先用优化过后的结构去建立的2x2x2的超晶胞 再删除我想要的原子数，(還是說必须要有真空层吗，看到教学文件指出非周期性的的分子或晶体才需要使用真空层预防波函数重叠)，再进行一次结构最佳化计算我要的物性，選用的截止能量跟k-point都是沿用这样对吗？ 题外话，另外想问一下因为大概是一百颗原子，以现在的Core i7 3.4 GHz的四核8线的电脑大概需要多少时间有哪位高手知道吗 希望大家一起交流指教，谢谢 [ Last edited by yaoo_1 on 2011-9-9 at 01:01 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 MaterialStudio学习使用 | 第一性原理 | 杂七杂八 |
» 猜你喜欢
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- surflex-dock对接问题
已经有3人回复
- 计算Cl-在Fe表面吸附问题时体系加电荷方法
已经有27人回复
- 有关MS 的建模 优化及计算的一些问题
已经有6人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 关于ATK建模时原子掺杂的问题
已经有3人回复
- 【讨论】分子模拟还是同源建模
已经有16人回复
- 【讨论】如何利用gaussian来解释关于不同合成条件下得到的两种配合物（50金币）
已经有20人回复
★
uuv2010(金币+1): 欢迎参加交流 2011-09-11 10:21:11
uuv2010(金币+1): 欢迎参加交流 2011-09-11 10:21:11
|
看了帖子，进来学习，按照楼主说的： “超晶胞是对原胞的扩展，比如3*3，2*2等，扩展而成新的重复单元。超晶胞一般用于研究缺陷，比如空位浓度，或者掺杂原子等。 例如：如果你的晶胞有100个原子，空位浓度为0.01，那么就要把晶胞扩展10倍，然后，从supercell里面去除一个原子。这样的空位浓度才等于0.01” 我觉得按上面的做法来处理楼主的问题是行得通的。 另外，楼主要计算的体系包含100个左右的原子，8G的内存是必须的，可以看看下面这个帖子： http://muchong.com/bbs/viewthread.php?tid=2481753&fpage=1 帖子的作者测试了castep计算所需内存量与计算体系的关系，没有给出计算时间和体系的关系。我估计楼主的Core i7 3.4 GHz的四核8线电脑算包含100个左右原子、而且带缺陷的体系时间上会化的比较久。 |
8楼2011-09-10 20:12:26
★
uuv2010(金币+1): 欢迎参加交流 2011-09-17 09:08:56
uuv2010(金币+1): 欢迎参加交流 2011-09-17 09:08:56
|
不知道你这个问题解决没有，如果解决了，可以发帖和大家交流一下。 说说我的看法： 看了你发的文件，里面出现了2个warnings： Secondd_get_u0: Warning: Band structure did not converge for 1 out of 4 kpoints Consider increasing BS_MAX_ITER or BS_EIGENVALUE_TOL. Secondd_minimise_phonon_E2: Can not converge to requested precision in 100 iterations Best convergence achievable is 0.234633E-03eV/A**2 Consider increasing PHONON_MAX_CYCLES or convergence tolerance PHONON_ENERGY_TOL 两个“Consider”已经给出了解决方法。你给出的文件里前面部分的两个参数值如下： Info: BS_MAX_ITER modified from 60 to 250 Info: BS_EIGENVALUE_TOL modified from 0.100000E-05 eV to 0.100000E-08 eV 这个应该是你改动过的，可以再增大一点。 另外，你的计算中这两个参数的设置如下： phonon convergence tolerance : 0.5000E-05 eV/A**2 max. number of phonon cycles : 100 也可以适当改大些。 建议不要同时变动上面四个参数，先试试前两个，再改后两个。 PS：本人没有算过声子，上面也仅是看你给的文件而想到可能的解决方法。祝好运。 |
11楼2011-09-16 21:53:49