QE个从编译到应用（不定期更新）

万里云 · 发表于 Post on 2015-9-22 14:11:27

本帖最后由万里云于 2015-9-22 14:12 编辑

茶味果冻发表于 2015-9-22 10:44
您好！按照说的我调了，但是运算有出现错误了，附件是我的band计算的输入文件和scf、bands计算的输出文 ...

从input_tmpband.in来看，CELL_PARAMETERS和K_POINTS后面都缺少参数，导致程序认错了跟在后面的坐标类型。

这很危险。建议看手册把缺少的参数补上。

茶味果冻 · 发表于 Post on 2015-9-22 15:07:33

万里云发表于 2015-9-22 14:11
从input_tmpband.in来看，CELL_PARAMETERS和K_POINTS后面都缺少参数，导致程序认错了跟在后面的坐标类型 ...

OK~

茶味果冻 · 发表于 Post on 2015-12-21 11:38:47

楼主好久不更新了啊，我能催更嘛

，相看看后续的具体计算

冰山 · 发表于 Post on 2015-12-21 13:42:05

楼主，麻烦问两个问题：
1.
http://theossrv1.epfl.ch/Main/Pseudopotentials 里面的相对论norm conserving pseudopotentials ，好用么？
需要特别设置什么么？
2.
计算soc的时候，如果必须用模守恒的赝势，有什么好的赝势推荐么？有什么注意事项？

万里云 · 发表于 Post on 2015-12-21 14:26:08

冰山发表于 2015-12-21 13:42
楼主，麻烦问两个问题：
1.
http://theossrv1.epfl.ch/Main/Pseudopotentials 里面的相对论norm conser ...

1. 没用过这里的赝势，无法给出参考；

2. QE赝势选择界面上可以看到每种赝势的信息，其中标有"Classification controlled by xxxxx“的比较可靠。

冰山 · 发表于 Post on 2015-12-23 10:12:10

万里云发表于 2015-12-21 14:26
1. 没用过这里的赝势，无法给出参考；

2. QE赝势选择界面上可以看到每种赝势的信息，其中标有"Classif ...

多谢楼主啦。
我用的元素需要用到的摸守恒赝势，没有完全相对论修正的。
QE网站上还推荐了一些赝势来源，赝势页面左侧 More about pseudopotentials
这里面有些赝势库还不错。可惜也没有我要用到的，

打算转openmx试试了。

万里云 · 发表于 Post on 2015-12-23 19:35:47

冰山发表于 2015-12-23 10:12
多谢楼主啦。
我用的元素需要用到的摸守恒赝势，没有完全相对论修正的。
QE网站上还推荐了一些赝势来源 ...

哦。那祝好。

didi_dudu · 发表于 Post on 2016-4-12 16:28:58

大神求组我最近在算AlN的phonon frequency，计算时候发现设置不同的k点，在ph那一部分输出的对称性不一样。比如如果用自动取点，ph输出的对称性就是C1，然后我用了这种取点方式
==============
K_POINTS
  20
0.1250000  0.1250000  0.1250000 1.00
0.1250000  0.1250000  0.3750000 1.00
0.3750000  0.1250000  0.1250000 2.00
0.1250000  0.1250000  0.6250000 1.00
0.6250000  0.1250000  0.1250000 2.00
0.1250000  0.1250000  0.8750000 1.00
0.8750000  0.1250000  0.1250000 2.00
0.1250000  0.3750000  0.3750000 2.00
0.3750000  0.3750000  0.1250000 1.00
0.1250000  0.3750000  0.6250000 2.00
0.6250000  0.3750000  0.1250000 2.00
0.1250000  0.6250000  0.3750000 2.00
0.1250000  0.3750000  0.8750000 2.00
0.3750000  0.8750000  0.1250000 2.00
0.1250000  0.8750000  0.3750000 2.00
0.1250000  0.6250000  0.6250000 2.00
0.6250000  0.6250000  0.1250000 1.00
0.3750000  0.3750000  0.3750000 1.00
0.3750000  0.6250000  0.3750000 2.00
0.3750000  0.3750000  0.6250000 3.00
=============================
ph部分输出的对称性就是c3v。
然而文献中aln的对称性是c6v。
想请问一下k点的选取和对称性有关系么？如果有的话怎么设置k点可以得到c6v的对称性呢。或者qe里面是否有可以指定对称性的功能。我查了相关手册等，都没有找到讲解。
谢谢！

万里云 · 发表于 Post on 2016-4-14 09:29:38

didi_dudu 发表于 2016-4-12 16:28
大神求组我最近在算AlN的phonon frequency，计算时候发现设置不同的k点，在ph那一部分输出的对称性不一样 ...

我群论学得很烂的。

1. AlN似乎是六方晶系的，为了同时取到M点和K点，MP网格应该是6的倍数，最小是6*6*2这样的，从你的数据看不是六等分的；
2. 和中心反演、镜面对称等相关的对称性和体系几何中心相关，把体系几何中心放到（0，0，0）和(0.5,0.5,0.5)处可能会提高对称性；

didi_dudu · 发表于 Post on 2016-4-15 15:52:06

万里云发表于 2016-4-14 09:29
我群论学得很烂的。

1. AlN似乎是六方晶系的，为了同时取到M点和K点，MP网格应该是6的倍数，最小是 ...

意思是用自动取点的方式取6-6-2这样？我尝试了自动取点6-6-2，输出中还是不能得到C6V对称性。大神能不能举个简单例子？另外，体系几何中心放到（0,0,0）处是怎么体现的呢？
------------------------------------------
如能回答万分感谢！
不能回答也万分感谢！

万里云 · 发表于 Post on 2016-4-16 09:58:15

didi_dudu 发表于 2016-4-15 15:52
意思是用自动取点的方式取6-6-2这样？我尝试了自动取点6-6-2，输出中还是不能得到C6V对称性。大神能不能 ...

就是求原子坐标的平均值，如果平均值不是(0,0,0)，就把所有原子坐标都平移平均值和(0,0,0)的差值。

shaofeng · 发表于 Post on 2016-12-20 22:47:47

万里云发表于 2015-6-24 21:17
一、QE源码下载、编译及测试

1.1 源码下载

非常好的编译方法，但是我用mkl库编译完以后，采用lda+U计算含有铁和砷的晶体时，结果中会出现很多Intel MKL ERROR: Parameter 10 was incorrect on entry to DGEMM . 的错误，但是最终优化的结构是合理的。请教下是什么原因。采用的是模守恒赝势，但是采用超软赝势时没有问题。由于模守恒赝势更接近实验值，希望能采用模守恒赝势计算。

obaica · 发表于 Post on 2017-6-9 16:50:34

楼主求更新啊，拯救我等小白！谢谢啊！

ChemG · 发表于 Post on 2021-5-21 22:36:42

真实用好方法，虽然并没有进行更多比对，使用i7-9750H的Y7000，仅通过两个原子si原胞vc_relax测试，相比gnu编译的44s，此方法编译的能提升至仅需13s。后面会在工作站上安装再测试

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Quantum ESPRESSO] QE个从编译到应用（不定期更新）

评分 Rate

评分 Rate

评分 Rate