| 查看: 1290 | 回复: 16 | |||
padian木虫 (著名写手)
|
[求助]
请教一个最低能量的问题 已有1人参与
|
|
请问各位，我用QE计算不同a下Total Energy，由此来根据最低Total Energy来判断Lattice Parameter是多少，不知道这种思路是否正确？ QE计算的Total Energy包括很多部分，是看具体比较不同a下的某一个能量还是系统使用计算给出的Total Energy来进行比较呢？ ------------------------------- 我基本上是用example01进行修改的，对Cu在FCC结构下，使用一个Cu原子，变化a从3.5~3.8A（7~7.6），来进行scf计算。现在的计算结果是单调的，并没有存在能量最低点。 请教大家！ 谢谢 |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一个非常简单的拉式变换
已经有9人回复
- 请教一个simulink问题
已经有6人回复
- 请教一个电子能量的问题
已经有6人回复
- 求助：VASP单点能计算没有找到能量最低点
已经有6人回复
- 请教一个规律问题
已经有19人回复
- 请教一个关于单点能校正的问题
已经有11人回复
- 同样的一个样品，同样的热激励，有的地方比其他地方温度变化大说明什么？
已经有14人回复
- 固定一个键长的同时 扫描二面角？
已经有11人回复
- 小白问题：吉布斯自由能和活化能
已经有6人回复
- 请教NEB收敛判据中能量的设置
已经有6人回复
- 新手请教用NEB方法计算energy barrier.
已经有15人回复
- 【求助】请教一个固体物理的问题
已经有15人回复
- 计算形成能要用到水的能量，请教各位 怎么计算水的能量
已经有8人回复
- 【求助】关于玻尔兹曼公式
已经有29人回复
- 【求助】请教一个基层工作最低年限的问题
已经有5人回复
- 【求助】收敛性测试的问题
已经有40人回复
mywai520
铁杆木虫 (著名写手)
- 应助: 339 (大学生)
- 金币: 10425.6
- 红花: 45
- 帖子: 1664
- 在线: 296.9小时
- 虫号: 265142
- 注册: 2006-07-09
- 专业: 凝聚态物性 II ：电子结构
2楼2014-04-23 18:08:22
padian
木虫 (著名写手)
- 应助: 51 (初中生)
- 金币: 2399.2
- 散金: 841
- 红花: 7
- 帖子: 1141
- 在线: 229.3小时
- 虫号: 500995
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
3楼2014-04-23 18:35:35
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2014-04-23 21:10:49
mywai520
铁杆木虫 (著名写手)
- 应助: 339 (大学生)
- 金币: 10425.6
- 红花: 45
- 帖子: 1664
- 在线: 296.9小时
- 虫号: 265142
- 注册: 2006-07-09
- 专业: 凝聚态物性 II ：电子结构
5楼2014-04-23 21:40:40
mywai520
铁杆木虫 (著名写手)
- 应助: 339 (大学生)
- 金币: 10425.6
- 红花: 45
- 帖子: 1664
- 在线: 296.9小时
- 虫号: 265142
- 注册: 2006-07-09
- 专业: 凝聚态物性 II ：电子结构
6楼2014-04-23 21:41:38
padian
木虫 (著名写手)
- 应助: 51 (初中生)
- 金币: 2399.2
- 散金: 841
- 红花: 7
- 帖子: 1141
- 在线: 229.3小时
- 虫号: 500995
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
|
我又进一步扩充了a的范围，从3.0~4.0A（间隔为0.5A） 分别scf计算Total Energy 我只使用了1个Cu 原子 计算发现，能量最低点对应的晶格常数在3.2A附近，但是实际上Cu的晶格常数是在3.6A附近。 不知道为什么？我换了好几个赝势，都是差不多在3.2A能量最低。 我的文件如下： ----------------------------------------------------------------------- # self-consistent calculation a = 3.0 A cat > Cu300.scf.david.in << EOF &control calculation = 'scf' restart_mode='from_scratch', prefix='Copper300', tstress = .true. tprnfor = .true. pseudo_dir = '$PSEUDO_DIR/', outdir='$TMP_DIR/' / &system ibrav= 2, celldm(1) =6.00, nat= 1, ntyp= 1, ecutwfc =21.5, occupations='smearing', smearing='gaussian', degauss=0.02 / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES Cu 63.546 Cu.pz-d-rrkjus.UPF ATOMIC_POSITIONS Cu 0.00 0.00 0.00 K_POINTS (automatic) 12 12 12 0 0 0 |
7楼2014-04-23 22:00:38
padian
木虫 (著名写手)
- 应助: 51 (初中生)
- 金币: 2399.2
- 散金: 841
- 红花: 7
- 帖子: 1141
- 在线: 229.3小时
- 虫号: 500995
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
8楼2014-04-23 22:02:13
mywai520
铁杆木虫 (著名写手)
- 应助: 339 (大学生)
- 金币: 10425.6
- 红花: 45
- 帖子: 1664
- 在线: 296.9小时
- 虫号: 265142
- 注册: 2006-07-09
- 专业: 凝聚态物性 II ：电子结构
9楼2014-04-23 22:20:27
padian
木虫 (著名写手)
- 应助: 51 (初中生)
- 金币: 2399.2
- 散金: 841
- 红花: 7
- 帖子: 1141
- 在线: 229.3小时
- 虫号: 500995
- 注册: 2008-02-13
- 专业: 凝聚态物性 II ：电子结构
10楼2014-04-23 22:59:52