留言板
- [8]dengfengli
- 很有可能，我来试一下，我把自洽计算时的K点设为32 32 1，这样的话声子计算需不需要重新计算。声子ph.x计算时需要wcf和save文件，这些文件是自洽计算时的还是电子计算时的，或者说两个都需要，谢谢
- [7]dengfengli
- 你好，请教你一下关于QE算电声耦合方面的，我算的是个二维片，我自洽计算时，K点设为24 24 1，电子计算K点设为16 16 1,声子计算时设的q点为8 8 1，计算过程没有问题，不过算力常数时最后报错， 输入文件是： zasr='simple', fildyn='al.dyn', flfrc='Al444.fc', la2F=.true.，报错放在下面的，请给点指导，谢谢
q-space grid ok, #points = 64
fft-check warning: sum of imaginary terms = .2414214E-07
Preparing gamma for a2F
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Error in routine init (1):
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Error in routine init (1):
nc already filled: wrong q grid or wrong nr - 我的回复(2013-12-7 00:22)：不大清楚，感觉k、q点设置，可能都需要成倍数关系
- [6]dengfengli
- 你好，请教你一下关于QE算电声耦合方面的，我算的是个二维片，我自洽计算时，K点设为24 24 1，电子计算K点设为16
- [5]黄林科
- 你好，由Quasiharmonic Debye Approximation（QDA）方法计算固体相关的热力学参数。
在这一模型中，先得到helmholtz自由能F
恒温下，F对V求偏导得到P。
最后，G=F+PV，这里的压强P什么含义呢？
计算固体Gibbs自由能的时候，P指的是什么压强？
如果固体膨胀，产生负压P（P<0）,计算G的时候，P如何定义？ - 我的回复(2013-8-14 22:24)：没做过相关计算，不过觉得，非高压条件下，对于固体，PV项应该可以忽略。关注此问题的最终答案。
- [4]杨创华
- 你好！我看了你的关于DFPT声子的两个帖子，http://blog.sciencenet.cn/blog-567091-675179.html http://blog.sciencenet.cn/blog-567091-675178.html 一个帖子说“即使用包含一个原子的Primitive cell，计算得到的Dynamical Matrix也是很准确的” ，一个帖子说“采用Phonopy软件，扩胞，一般应包含接近100个原子，超胞的晶格尺寸应达到10A”，感觉有点矛盾，如果单包的精度就很高了，没有必要增加计算量。我现在就是用DFPT算声子的时候出现了虚频，所以很矛盾，不知道有么有必要增加超胞的大小，如果增加的话内存又不够。
- 我的回复(2013-6-25 03:43)：你好，我也在学习中。你提到的两句话，分别对应不同的计算方法，一个是有限位移法，一个是线性响应的方法。出现虚频多半是由于结构不合理引起的，可以考虑对结构进行多次优化，直到完全relaxed.
- [3]叶小球
- 这个似乎有点复杂：
1） 网上有相关程序，但我没用过，你可试下；
2） 我自己用的是组里有人修改过的VASP源程序，怎么改的不大清楚，但这显然不是好方法。
- [2]xcid
- 请问如何用VASP画费米面？
- [1]李曼静
- 您好！
我是华中师范大学的研二学生，这是本人读研以来为第一篇论文设计的问卷。是以导师的项目为主题，研究学术虚拟社区用户的持续使用行为，需要大量的问卷调查数据为研究基础。此项研究有助于了解在线学术交流的现状、趋势和影响在线学术交流的因素，丰富与深化互联网环境中科研知识交流与共享机制研究的理论体系，从而为科研性的网络服务在中国的发展提供理论指导。本调查完全采用匿名方式进行，您的个人回答只供学术研究使用，绝对不会对您造成任何不良影响。对您的参与我将不胜感激，祝您学习、工作顺利！问卷地址: http://www.sojump.com/jq/2307285.aspx
