计算化学公社

标题: 能带计算必须在原胞上吗 [打印本页]

作者
Author: Theor_Comp 时间: 2020-9-12 22:21
标题: 能带计算必须在原胞上吗
最近在学vasp能带计算，有些疑惑：1. 看到有些说计算能带时，由于只能在原胞上选择高对称K点，所以计算能带时要使用原胞的结构。那convetional cell或者自己构建的超胞啥的，就不能计算能带了吗？
2. PBE泛函计算能带使用的K点是line-mode形式，而HSE这种计算能带的话，K点使用的是IBZKPT里面的吗？

作者
Author: Theor_Comp 时间: 2020-9-26 20:18
@卡开发发请问您可以帮忙看下吗？非常感谢。

作者
Author: kyuu 时间: 2020-9-28 00:06
那就去问有些啊

作者
Author: 卡开发发 时间: 2020-9-28 08:54

Theor_Comp 发表于 2020-9-26 20:18
@卡开发发请问您可以帮忙看下吗？非常感谢。

1、超晶格得到的能带是折叠的，固体物理的书对这个问题都有论述。
2、可以把产生的k拷贝过来用。

作者
Author: Theor_Comp 时间: 2020-9-28 22:28

卡开发发发表于 2020-9-28 08:54
1、超晶格得到的能带是折叠的，固体物理的书对这个问题都有论述。
2、可以把产生的k拷贝过来用。

谢谢您。
1、意思是超胞也可以计算能带，只不过它是折叠的，对吗？
2、PBE泛函和HSE这两种，计算能带时都可以把IBZKPT里产生的k直接拷贝过来用吗？还是说只适用于HSE计算能带？

作者
Author: 卡开发发 时间: 2020-9-28 23:37

Theor_Comp 发表于 2020-9-28 22:28
谢谢您。
1、意思是超胞也可以计算能带，只不过它是折叠的，对吗？
2、PBE泛函和HSE这两种，计算能带时 ...

1、对。
2、用MP网格加权重，高对称点+路径用0权重的。

作者
Author: Theor_Comp 时间: 2020-9-29 18:55

卡开发发发表于 2020-9-28 23:37
1、对。
2、用MP网格加权重，高对称点+路径用0权重的。

1. “用MP网格加权重，高对称点+路径用0权重”，这个K点，是专门针对使用HSE杂化泛函的K点吧？
2. “MP网格加权重”，这个可以直接从IBZKPT中得到，那“高对称点+路径”这个怎么得到呢？
非常感谢您！

作者
Author: 卡开发发 时间: 2020-9-30 22:20

Theor_Comp 发表于 2020-9-29 18:55
1. “用MP网格加权重，高对称点+路径用0权重”，这个K点，是专门针对使用HSE杂化泛函的K点吧？
2. “MP ...

1、对。
2、其他程序生成，或者知道高对称点用plain DFT算完能带应该在IBZKPT得到吧。

作者
Author: Theor_Comp 时间: 2020-10-1 08:24

卡开发发发表于 2020-9-30 22:20
1、对。
2、其他程序生成，或者知道高对称点用plain DFT算完能带应该在IBZKPT得到吧。

再次感谢。
但现在生成高对称点路径的程序，像seek-path，好像都是针对原胞而找出的高对称点。那对于超胞的高对称点，该如何生成呢？

作者
Author: 卡开发发 时间: 2020-10-1 09:27

Theor_Comp 发表于 2020-10-1 08:24
再次感谢。
但现在生成高对称点路径的程序，像seek-path，好像都是针对原胞而找出的高对称点。 ...

按照文献http://bbs.keinsci.com/thread-37-1-1.html和晶型去对照着设置。

作者
Author: Theor_Comp 时间: 2020-10-5 10:12

卡开发发发表于 2020-10-1 09:27
按照文献http://bbs.keinsci.com/thread-37-1-1.html和晶型去对照着设置。

谢谢您。
这个文献是只针对超胞吗，适用于原胞/单胞吗？

作者
Author: 卡开发发 时间: 2020-10-7 21:29

Theor_Comp 发表于 2020-10-5 10:12
谢谢您。
这个文献是只针对超胞吗，适用于原胞/单胞吗？

原则上是可以，只要对称性满足要求就行，但超晶格的情况如上面所说会出现能带折叠。

作者
Author: 珊珊来迟 时间: 2020-10-21 18:31

卡开发发发表于 2020-10-7 21:29
原则上是可以，只要对称性满足要求就行，但超晶格的情况如上面所说会出现能带折叠。

老师我看了这篇文献，还有使用seek-path网站设置高对称点，都是要把输入的结构先转换到standard primitive cell，然后如果要使用他所提供的k-path，就要基于这个standard primitive cell来计算能带。我现在的结构是实验来的单晶，我要计算这个单晶的电输运性质，这其中我要分析能带结构，还要基于能带结构计算有效质量。但用standard primitive cell和我的实验晶体计算出来的能带结构是不一样的，我该分析哪一个呢？
如果用standard primitive cell能带的结果，就不能与我实际晶体的堆积结构特点对应起来分析（因为在转化为primitive cell的时候基矢的顺序都改变了），但如果用实验的晶体结构，我不知道如何得到他的高对称点和k-path（实验晶体不是标准的原胞，有一个晶体还不是最小重复单元，无法与这篇文献中给出的晶型去对照设置）？

现在我是用的我的实验晶体结构，路径选择的是seek-path网站给出的路径计算的能带，还基于这个能带计算了有效质量，拿这个能带和有效质量的结果去进行讨论会有什么问题吗？可能我基于这个能带得到的vbm和cbm并不是真正的vbm/cbm?所以计算出来的有效质量也不一定正确？
问题困扰学生多时，期待您的回复！

作者
Author: 卡开发发 时间: 2020-10-22 23:01
本帖最后由卡开发发于 2020-10-22 23:05 编辑

k点采样足够的情况，超晶格和单胞的VBM、CBM的位置应该是不会变的。只是从能带图看，超晶格会有能带折叠现象，但不会改变价带顶或导带底附近的色散行为。不过因为对称性变化了，这样高对称点也会发生变化，稳妥起见做3维的有效质量会好些。@珊珊来迟

作者
Author: sciwizard 时间: 2023-12-5 17:13

卡开发发发表于 2020-10-22 23:01
k点采样足够的情况，超晶格和单胞的VBM、CBM的位置应该是不会变的。只是从能带图看，超晶格会有能带折叠现 ...

老师您好，我在对掺杂体系进行能带计算时，结构优化后，对称性改变了，然后高对称点使用的是原包的高对称点。将静态的Incar改成能带计算的Incar，提交任务之后，不开始电子步的迭代，请问这种情况应该如何解决？

作者
Author: 卡开发发 时间: 2023-12-5 19:20

sciwizard 发表于 2023-12-5 17:13
老师您好，我在对掺杂体系进行能带计算时，结构优化后，对称性改变了，然后高对称点使用的是原包的高对称 ...

不开始电子步迭代的原因很多，这里的信息不足以判断。

作者
Author: sciwizard 时间: 2023-12-5 19:31
本帖最后由 sciwizard 于 2023-12-5 20:40 编辑

卡开发发发表于 2023-12-5 19:20
不开始电子步迭代的原因很多，这里的信息不足以判断。

您好，这是我的一些计算信息，刚开始接触计算，不大清楚是哪里的问题，麻烦您帮忙看一下。

LDA part: xc-table for Pade appr. of Perdew
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
WAVECAR not read
charge-density read from file:
WARNING: random wavefunctions but no delay for mixing, default for NELMDL

===================================================================================
= BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
= RANK 0 PID 1837375 RUNNING AT c4111
= KILLED BY SIGNAL: 9 (Killed)

作者
Author: 卡开发发 时间: 2023-12-5 21:08

sciwizard 发表于 2023-12-5 19:31
您好，这是我的一些计算信息，刚开始接触计算，不大清楚是哪里的问题，麻烦您帮忙看一下。

LDA part: ...

那是直接挂了，但原因不清楚，你先看看OUTCAR有啥提示。

作者
Author: sciwizard 时间: 2023-12-6 10:14

卡开发发发表于 2023-12-5 21:08
那是直接挂了，但原因不清楚，你先看看OUTCAR有啥提示。

刚开始学这个东西，不大能搞懂哪里出错了。能麻烦您帮忙看一下嘛？OUTCAR写的这个对称性和我给的高对称点的路径不能对上，我给出的是原包的高对称点。

vasp.5.4.4.18Apr17-6-g9f103f2a35 (build May 27 2019 14:10:10) complex

executed on          LinuxIFC date 2023.12.05  18:45:19
running on 64 total cores
distrk:  each k-point on 64 cores, 1 groups
distr:  one band on NCORES_PER_BAND= 8 cores, 8 groups

Analysis of symmetry for initial positions (statically):
=====================================================================
Subroutine PRICEL returns:
Original cell was already a primitive cell.

Routine SETGRP: Setting up the symmetry group for a
simple monoclinic supercell.

Subroutine GETGRP returns: Found  1 space group operations
(whereof  1 operations were pure point group operations)
out of a pool of  4 trial point group operations.

The static configuration has the point symmetry C_1 .

Analysis of symmetry for dynamics (positions and initial velocities):
=====================================================================
Subroutine PRICEL returns:
Original cell was already a primitive cell.

Routine SETGRP: Setting up the symmetry group for a
simple monoclinic supercell.

Subroutine GETGRP returns: Found  1 space group operations
(whereof  1 operations were pure point group operations)
out of a pool of  4 trial point group operations.

The dynamic configuration has the point symmetry C_1 .

Subroutine INISYM returns: Found  1 space group operations
(whereof  1 operations are pure point group operations),
and found    1 'primitive' translations

Static calculation
charge density and potential remain constant during run
non-spin polarized calculation
Variant of blocked Davidson
Davidson routine will perform the subspace rotation
perform sub-space diagonalisation
after iterative eigenvector-optimisation
no mixing
using additional bands          46
reciprocal scheme for non local part
use partial core corrections
calculate Harris-corrections to forces
(improved forces if not selfconsistent)
use gradient corrections
WARNING: stress and forces are not correct
  (second derivative of E(xc) not defined)
use of overlap-Matrix (Vanderbilt PP)
Gauss-broadening in eV    SIGMA  = 0.03

maximum and minimum number of plane-waves per node :    7353    7244

maximum number of plane-waves:    58580
maximum index in each direction:
IXMAX= 24 IYMAX= 24 IZMAX= 26
IXMIN=  -25 IYMIN=  -25 IZMIN=  -26

parallel 3D FFT for wavefunctions:
minimum data exchange during FFTs selected (reduces bandwidth)
parallel 3D FFT for charge:
minimum data exchange during FFTs selected (reduces bandwidth)

total amount of memory used by VASP MPI-rank0  2626912. kBytes
=======================================================================

base    :    30000. kBytes
nonl-proj : 1595139. kBytes
fftplans  :    9690. kBytes
grid    :    18168. kBytes
one-center:       360. kBytes
wavefun :    973555. kBytes

   INWAV:  cpu time 0.0000: real time 0.0000
initial charge density was supplied:
number of electron    386.9999998 magnetization
keeping initial charge density in first step

--------------------------------------------------------------------------------------------------------

Maximum index for augmentation-charges       1376 (set IRDMAX)

作者
Author: 卡开发发 时间: 2023-12-6 13:16

sciwizard 发表于 2023-12-6 10:14
刚开始学这个东西，不大能搞懂哪里出错了。能麻烦您帮忙看一下嘛？OUTCAR写的这个对称性和我给的高对称点 ...

这点我暂时不太能确定，也许你可以试一下在执行vasp之前执行ulimit -s unlimited。
对称性的问题你也许可以先确认一下之前算出来的结构对称性的情况，以及适当调整symprec试试。

作者
Author: sciwizard 时间: 2023-12-6 13:59

卡开发发发表于 2023-12-6 13:16
这点我暂时不太能确定，也许你可以试一下在执行vasp之前执行ulimit -s unlimited。
对称性的问题你也许 ...

好的，多谢，麻烦您啦

作者
Author: walxf1987 时间: 2024-9-29 09:32

sciwizard 发表于 2023-12-5 17:13
老师您好，我在对掺杂体系进行能带计算时，结构优化后，对称性改变了，然后高对称点使用的是原包的高对称 ...

请问下老师我也遇到过这种情况，要怎么解决呢

作者
Author: walxf1987 时间: 2024-10-11 08:15

sciwizard 发表于 2023-12-5 17:13
老师您好，我在对掺杂体系进行能带计算时，结构优化后，对称性改变了，然后高对称点使用的是原包的高对称 ...

你好，请问你这个问题解决了吗

作者
Author: weilu 时间: 2024-10-31 19:06

sciwizard 发表于 2023-12-5 17:13
老师您好，我在对掺杂体系进行能带计算时，结构优化后，对称性改变了，然后高对称点使用的是原包的高对称 ...

您好请问问题解决了吗？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)