能带反折叠计算

最新推荐文章于 2024-05-14 15:32:36 发布
最新推荐文章于 2024-05-14 15:32:36 发布
阅读量8.2k 收藏 15
点赞数 3
分类专栏： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wwxy1995/article/details/103078186
版权

做一下投影看看轨道成分

第一步：准备以下文件：

  1. POSCAR文件，必须是原胞结构（primitive cell）。以下是MoS2单层结构文件

    MoS2-H.POSCAR
  1.00000000000000000
3.190406376    0.000000    0.0000000
-1.59520318   2.7629729    0.0000000
0.000000000    0.000000   17.6585355
  Mo   S
    1     2
Direct
0.33333334 0.66666668 0.50000000
0.66666668 0.33333334 0.41170027
0.66666668 0.33333334 0.58829969

     

  2. 包含原胞能带路径的KPATH.in文件，可通过vaspkit 302或303产生，也可手动调整能带路径。这里我们选择的MoS2能带路径是从M-K-Gamma。

    KPATH for MoS2
20
Line mode
Rec
0.0000000 0.5000000 0.0000000 # M
0.3333333 0.3333333 0.0000000 # K

0.3333333 0.3333333 0.0000000 # K
0.0000000 0.0000000 0.0000000 # GAMMA

  3. 用于静态计算的INCAR文件，并设置LWAVE =.TRUE.，为了避免读取WAVECAR时出现'ERROR! Computed NPLANE= ** != input NPLANE= **'错误，可选择较低的截断能或令PREC=N。INCAR的典型参数为

    SYSTEM = MoS2
ISTART = 1
LREAL = F
PREC = N
LWAVE = .TRUE.
LCHARG = F
ISMEAR = 0
SIGMA = 0.05
NELM = 60
NELMIN = 6
EDIFF = 1E-08

     

  4. 用于产生超胞的变化矩阵TRANSMAT.in文件（可选文件）。我们先看看TRANSMAT.in的格式

    Read transformation matrix from the TRANSMAT.in file if it exists. 第一行是注释行
2 0 0   # must be three integers
0 3 0   # must be three integers
0 0 3   # must be three integers

第二行至第四行的整数构成3*3变换矩阵M，超胞晶格基矢A和原胞晶格基矢a满足以下的变换关系：

如果M中非对角元素都为0，则产生的超胞晶格基矢与原胞的平行，这种超胞也是最常见的，也可通过vaspkit-401命令产生。

 

MoS2原胞属于六角结构，这里我们想产生一个矩形超胞，我们可以通过vaspkit-400命令实现。注意如果TRANSMAT.in不存在时，运行vaspkit-400后变换矩阵M每个元素通过对话框输入，如下图命令所示。否则，直接从TRANSMAT.in文件中读取换矩阵M，这样更方便快捷。

vaspkit -task 400
+---------------------------------------------------------------+
|           VASPKIT Version: 1.00.RC (17 Jun. 2019)             |
|       A Pre- and Post-Processing Program for VASP Code       |
|           Running VASPKIT Under Command-Line Mode           |
+---------------------------------------------------------------+
-->> (01) Reading Structural Parameters from POSCAR File...
Enter the new lattice verctor a in terms of old:
(MUST be three integers, e.g., 1 2 3)
  4   0   0
Enter the new lattice verctor b in terms of old:
  2   4   0
Enter the new lattice verctor c in terms of old:
  0   0   1
+-------------------------- Summary ----------------------------+
The Transformation Matrix P is:
  4   0   0
  2   4   0
  0   0   1
Lattice Constants in New Cell: 12.762 11.052 17.659
Lattice Angles in New Cell:   90.00   90.00   90.00
Total Atoms in New Cell: 48
Volume of New Cell is 16 times of the Old Cell
+---------------------------------------------------------------+
-->> (02) Written SUPERCELL.vasp File

第三步：cp SUPERCELL.vasp POSCAR和cp TRANSMAT TRANSMAT.in（当TRANSMAT.in文件不存在时）, 并运行vaspkit-281生成KPOINTS文件；

------------>>
28
+-------------------------- Warm Tips --------------------------+
      See some examples in vaspkit/examples/band_unfolding.
        Only Support KPOINTS file Generated by VASPKIT.
          Please Set LWAVE= .TRUE. in the INCAR file.
+---------------------------------------------------------------+
================== Band-Unfolding ==============================
281) Generate KPOINTS File for Band-Unfolding Calculation
282) Calculate Effective Band Structure

0)   Quit
9)   Back
------------>>
281

-->> (01) Reading Structural Parameters from POSCAR File...
+---------------------------------------------------------------+
|               Selective Dynamics is Activated!               |
+---------------------------------------------------------------+
======================= K-Mesh Scheme ==========================
1) Monkhorst-Pack Scheme
2) Gamma Scheme

0)   Quit
9)   Back
------------->>
2
+-------------------------- Warm Tips --------------------------+
Input Resolution Value to Determine K-Mesh for SCF Calculation:
(Typical Value: 0.02-0.03 is Generally Precise Enough)
------------>>
0.03
Input Resolution Value along K-Path for Band Calculation:
(Typical Value: 0.02-0.04 for DFT and 0.04-0.06 for hybrid DFT)
------------>>
0.03
+---------------------------------------------------------------+
-->> (02) Readin Transformation Matrix from TRANSMAT.in file...
-->> (03) Reading K-Paths From KPATH.in File...
+-------------------------- Summary ----------------------------+
K-Mesh for SCF Calculation:   3   3   1   # 注意这部分根据超胞的倒格矢长度产生
The Number of K-Point along K-Path 1: 21     # 注意这部分根据原胞的倒格矢长度产生
The Number of K-Point along K-Path 2: 43     # 注意这部分根据原胞的倒格矢长度产生
+---------------------------------------------------------------+

第四步：提交vasp作业；

第五步：运行vaspkit-282；

------------>>
282
+-------------------------- Warm Tips --------------------------+
  Current Version Only Support the Stardard Version of VASP code.
+---------------------------------------------------------------+
-->> (01) Reading the header information in WAVECAR file...
+--------------------- WAVECAR Header --------------------------+
SPIN =   1
NKPTS =   68
NBANDS = 252
ENCUT =   280.00
Coefficients Precision: Complex*8
Maximum number of G values: GX =   18, GY =   16, GZ =   25
Estimated maximum number of plane-waves:   30159
+---------------------------------------------------------------+
-->> (02) Start to read WAVECAR file, take your time ^.^
          Percentage complete: 25.0%
          Percentage complete: 50.0%
          Percentage complete: 75.0%
          Percentage complete: 100.0%
-->> (03) Readin Transformation Matrix from TRANSMAT.in file...
-->> (04) Reading Fermi-Energy from DOSCAR File...
ooooooooo The Fermi Energy will be set to zero eV ooooooooooooooo
-->> (05) Reading KPOINTS_MAPPING_TABLE.in file...
-->> (06) Reading K-Paths From KPATH.in File...
-->> (07) Start to Calculate Effective Band Structure...
          Percentage complete: 25.0%
          Percentage complete: 50.0%
          Percentage complete: 75.0%
          Percentage complete: 100.0%
-->> (08) Written EBS.dat File!
-->> (09) Written KLABELS File!

第六步：最后，可以利用vaspkit/examples/band_unfolding/eps_plot.py脚本画图

 

注意，这个过程只需进行一次静态自洽计算，不需要非自洽取点

wwxy261
关注
  • 3
    点赞
  • 踩
  • 15
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
vaspkit使用_VASPKIT中轨道投影能带反折叠功能在二维异质结中的应用案例
weixin_42679995的博客
01-02 3257
二维异质结，尤其是旋转二维异质结是近年来十分热门的研究方向，接下来我们演示如何应用VASPKIT 1.20及之后版本中的轨道投影能带反折叠功能计算二维异质结的电子结构性质。为了便于比较，我们重复了Matsushita等人关于MoS2/graphene二维异质结能带计算的计算结果[Physical Review Materials, 2(1), 010801]，该模型由4×4 MoS2和...
能带反折叠计算能带
wwxy1995的博客
12-12 2991
没有畸变的晶格，能带反折叠计算
“掺杂三个原子后能带不连续了？” 何解？
icehoqion的博客
06-08 910
最近有小伙伴求助笔者，他在算能带反折叠的过程中遇到了一些困难：计算纯相模型的时候所得到的能带是连续的，但是扩胞后并掺杂了三个原子后，能带出现了不连续的情况。如图所示：这里没有标注费米能级以及高对称点位置。我们可以在横坐标x=2的位置左边看到有很明显的一处能带处于不连续的地方。 根据后续了解，这里描述了一下这位伙伴计算的一些细节。✓计算细节原胞模型：金刚石结构掺杂原子数：3计算类型：能带反折叠能带路径（KPATH）：与原胞相同使用软件：VASP、VASPKIT 看过计算细节其实差不多可以知道问题所
vaspkit生成vasp的KPOINT文件
最新发布
liouver的博客
05-14 1628
总结：如果要让vasp程序以Monkhorst-Pack和Gamma方法自动生成k点，则可以自己按编写上述5行内容的KPOINTS文件。选择vaspkit的第三种方法 3）Irreducible K-Points with Gamma Scheme生成k点，精度0.02。输入1，选择Monkhorst-Pack Scheme （最常用的方法，对于六方晶系则选 2 Gamma Scheme）第三行：Monkhorst-Pack表示以此方法生成k点，可以改成Gamma表示以Gamma方法生成k点。
独家解析：Comsol能带折叠演示全面剖析与应用实践
DEPBapO的博客
04-13 433
通过合理设置展示参数和添加适当的注释，折叠演示功能可以帮助我们更好地展示和传达自己的工作成果，提升工作效率和沟通效果。在得到仿真结果后，我们可以通过折叠演示功能将其展示出来。在展示模型时，我们可以通过折叠演示功能对关键的物理场进行高亮显示，突出展示分析的重点和结果。折叠演示是Comsol中一项非常实用的功能，它能够将复杂的仿真过程简化为一张图或一段动画，方便工程师、设计师等针对不同场景进行快速展示和交流。通过设置不同的时间步长和观察角度，我们可以观察到仿真结果随时间变化的过程，了解物理场的动态演化规律。
VSe2 CDW态能带计算
wwxy1995的博客
11-27 1755
第一步，得到单胞的POSCAR和POTCAR 第二部，推导CDW超胞与单胞的晶格矢量关系，写TRANSMAT.in文件，如下 Read transformation matrix from the TRANSMAT.in file if it exists. 第一行是注释行 2 0 0 # must be three integers 0 3 0 # must be three in...
vaspkit自动绘图设置教程
icehoqion的博客
08-14 1384
vaspkit自1.2.5版本后即更新了自动绘图功能，绘图的范围包括布里渊区示意图、能带图、态密度图、3d能带等等，均可在vasp计算完成后处理时自动绘制能带，在最新版本更新了可将能带和态密度绘制在一张图片上的功能。在每次绘图后，在文件夹中会存有一个名为PLOT.In 的文件，如对图片绘制有要求，可以自行更改PLOT.In文件然后重新执行命令，vaspkit会按该文件里的设置重新绘图。部分都是画图的设置，部分设置为vaspkit pro的设置内容，即便修改，对普通版vaspkit也不会生效。
【能带计算】基于平面波展开法实现二维声子晶体能带计算附Matlab代码
m0_57702748的博客
12-23 949
基于平面波展开法实现二维声子晶体能带计算是固体物理学中一个重要的研究课题。声子晶体是一种具有周期性结构的材料，其能带结构对材料的热学、光学和电子性质有着重要影响。因此，准确地计算声子晶体的能带结构对于理解和预测材料的性质具有重要意义。平面波展开法是计算固体电子结构和声子结构的常用方法之一。它通过将波函数或位移场展开为平面波的线性组合，将固体的晶格周期性结构考虑在内，从而可以较为准确地描述固体的能带结构。在二维声子晶体的能带计算中，平面波展开法可以很好地适用于描述声子的振动特性和能带结构。
ILspy 反编译软件
04-23
5. **颜色编码和代码折叠**：为了增强代码的可读性，ILSpy提供了代码高亮和折叠功能，使得复杂的代码结构更易于理解和管理。 6. **插件支持**：通过扩展插件，ILSpy可以增加更多功能，例如与其他版本控制系统集成，...
JAVA反编译工具jd-gui&luyten
01-30
luyten不仅拥有jd-gui的所有特性，还增加了代码折叠、颜色高亮、书签、代码对比、自定义热键等高级功能。此外，luyten支持插件扩展，可以集成到其他开发环境中，如IntelliJ IDEA或Eclipse，提升开发者的效率。 使用...
免费的.NET反编译工具ILSpy
07-11
6. **颜色编码和折叠代码**：源代码视图提供了颜色编码的语法高亮和代码折叠功能，便于阅读和理解复杂的代码结构。 7. **可视化依赖关系**：通过类图和调用图，ILSpy可以帮助用户直观地看到类之间的依赖关系和函数...
VASPKIT中文手册1
08-04
1. VASPKIT简介 4. 用户自定义界面（测试功能） 1. VASPKIT简介 4.生成晶体的能带路径（包括杂化泛函），并处理能带
反汇编工具
12-03
- ILSpy：免费且开源的.NET反编译器，提供代码高亮和折叠功能，支持C#和VB.NET。 - dotPeek：由JetBrains开发的强大工具，提供快速的反编译和导航功能。 - Reflector：早期的.NET反编译器，现在由Red Gate公司...
2021-2022计算机二级等级考试试题及答案No.5071.docx
11-02
17. 加密是通过加密算法将明文转化为密文，解密则是反过来。 18. WWW服务基于HTTP（超文本传输协议）。 19. SQL中的DELETE命令标记要删除的记录，但并不立即删除，需要COMMIT或ROLLBACK事务来确认或撤销。 20. ...
计算能带自动产生k_VASPKIT 1.2 支持全晶系KPath自动生成
weixin_39815286的博客
01-15 6383
相关往期回顾：VASPKIT 1.2正式版发布深入分析能带结构(一)深入分析能带结构(二)-VASPKIT能带图计算深入分析能带结构(三)-Origin画能带图深入分析能带结构(四)-高通量生成能带图pymatgen深入分析能带结构(五)-魔角石墨烯的构建与其能带结构浅析深入分析能带结构(六)-对赝原子做投影，电子化合物(下)深入分析能带结构(七)-VASPKIT计算有效质量深入分析能带...
科研进展7/3
wwxy1995的博客
07-03 274
总结： 1. 高温相结构如图1没问题，可以和APRES实验对上。 2. 从原始费米面nesting可以解释sqrt(3)xsqrt(7)的低温CDW 3. 可以用第一性的方法研究，不同CDW序，在做一个2*1的结果 4. 尝试用应力解释2*1 5.用能带反折叠的方法拿去和实验APRES对比。 ...
CDW态的能带计算
wwxy1995的博客
11-26 553
基本一致，Garma点简并性需要仔细考虑 布里渊区均匀变化，不需要反折叠技术。 能带反折叠技术 、
html折叠菜单栏附带css样式
05-31
实现HTML折叠菜单栏可以使用HTML、CSS和JavaScript。下面是一个简单的示例代码： HTML： ```html <nav class="nav"> <ul> <li> <a href="#">菜单项1</a> <ul> <li><a href="#">子菜单项1</a></li> <li><a href="#">子菜单项2</a></li> <li><a href="#">子菜单项3</a></li> </ul> </li> <li><a href="#">菜单项2</a></li> <li><a href="#">菜单项3</a></li> </ul> </nav> ``` CSS： ```css .nav ul { margin: 0; padding: 0; list-style: none; } .nav li { position: relative; margin: 0; padding: 0; } .nav a { display: block; padding: 10px; background-color: #f7f7f7; color: #333; text-decoration: none; } .nav a:hover { background-color: #ddd; } .nav ul ul { display: none; position: absolute; top: 100%; left: 0; background-color: #f7f7f7; z-index: 999; } .nav ul ul li { display: block; width: 100%; } .nav ul ul a { padding: 10px; color: #333; } .nav ul ul a:hover { background-color: #ddd; } ``` JavaScript： ```javascript var nav = document.querySelector('.nav'); var navItems = nav.querySelectorAll('li'); for (var i = 0; i < navItems.length; i++) { var navItem = navItems[i]; var subNav = navItem.querySelector('ul'); if (subNav) { navItem.classList.add('has-subnav'); navItem.addEventListener('click', function(event) { event.preventDefault(); if (this.classList.contains('active')) { this.classList.remove('active'); } else { var activeNavItem = nav.querySelector('.active'); if (activeNavItem) { activeNavItem.classList.remove('active'); } this.classList.add('active'); } }); } } ``` 以上代码实现了一个具有折叠效果的菜单栏，点击菜单项时可以展开或折叠其子菜单。需要注意的是，CSS中的绝对定位和JavaScript中的类名操作都是用来控制菜单项的展开和折叠的。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值