- 博客(93)
- 收藏
- 关注
RSS订阅原创 脚本分享：快速作图对比wannier拟合能带python脚本
本脚本通过Python实现电子能带结构数据的快速作图，能够从两个不同的数据文件（BAND.dat 和 wannier90_band.dat）中提取有效数据，并在同一坐标系下绘制对比图。准备工作：使用VASPKIT处理获得能带数据BAND.dat，通过wannier90.x拟合获得紧束缚模型的电子能带 数据wannier90_band.dat。注：BAND.dat中能量修改为未减去费米能级的绝对能量，如文件夹中保存有KLINES.dat文件，则用于绘制能带高对称点。
2025-04-25 21:49:25
115
原创 论文精读：NC CoNb3S6拓扑能斯特效应 附Peer Review
作为偏态密度中上自旋态和下自旋态之间的分离能，平均自由程lmfp由费米速度vF和弛豫时间τ的乘积计算，前者由角分辨光电能谱（ARPES）估算，其中，对于CoNb3S6 布里渊区边缘的Co衍生带，带色散的线性斜率定义为vF =（1/ℏ）ΔE/Δk，其中vF = 2.2·105 m/s。相关的2cWyckoff位置包括点r1 =（1/3，2/3，1/4）和r2 =（2/3，1/3，3/4），单位为原始晶格向量a1，a2，a3，定义为。在低温下，Sxx（T）中的拐点与在磁化测量中观察到的TN一致。
2025-04-06 20:04:35
920
原创 oneAPI编译VASP 中 icx icpx和ifx设置 适用于icc和ifort Command not found
近期VASP更新出6.5版本，在6.5.0的安装包中，makefile.include已更新出适用于新版oneapi的版本makefile.include.oneapi，同时还包括openmpi版本的makefile.include.oneapi_omp。最新版的安装包将不再支持旧版命令，比如笔者在以前帖子中已经提到过的icc，icpc的问题。在此makefile.include中，调整了对于ifx和icx icpx的设置。附makefile.include.oneapi。
2025-03-21 19:32:40
411
原创 libxc+VASP编译（补充）和杂化泛函使用
按照前文提到的编译方法，在更改泛函后进行VASP计算时会出现保存提示 Libxc was not compiled with the option --disable-fhc。查询VASP wiki后可知，VASP官网推荐的使用libxc.5.2.0以后的版本，并在编译过程中使用 --disable-fhc选项。但在实际测试中仍存在一些问题，本篇基于前期内容进行修正并附加测试。同时必须使用VASP.6.3.0或更高版本的VASP程序。在VASP的编译文件makefile.include中加入。
2025-03-21 18:45:37
421
原创 libxc+VASP编译
Libxc是密度泛函理论的交换相关泛函库。目的是提供一套可移植的、经过良好测试的、可靠的交换和相关泛函，供各种程序使用。本文使用libxc 6.1.0版本，下载后进行解压和生成configure文件。本文使用前提：已用于完整Intel编译器并可正常编译VASP程序。在VASP的编译文件makefile.include中加入。随后生成makefile文件，指定好安装目录。check无问题后进行安装。
2025-03-10 19:36:43
206
原创 VASPKIT 3D 能带计算及动态3D能带绘图脚本
how_to_visual.m对应于MATLAB代码， how_to_visual.py 对应于python代码，均可进行对BAND_HOMO.grd BAND_LUMO.grd KX.grd KY.grd 数据文件的读取和绘图。所生成的KP点包含1加权KP点和0加权KP点(注意，为了获得平滑的3D能带，用于生成KP点文件的分辨率值应该设置在0.01左右)3D能带即将常见绘制在布里渊区高对称点路径上的电子能带绘制在倒空间内，通过倒空间轴的立体坐标，将能带立体化。
2025-03-04 20:53:59
630
原创 论文精读： JACS 、 GdAlSi石墨化低维磁性
GdAlSi多晶型的自旋极化基态计算是通过维也纳从头计算模拟包（VASP）中周期性边界条件的投影仪增强波方法（PAW），采用广义梯度近似的Perdew−Burke−Ernzerhof（PBE）交换关联泛函，并进行DFT+U，参数U和J分别设置为6.7和0.7 eV，以改进Gd原子中4f态。根据（000n）XRD反射的位置计算出的三角GdAlSi的参数c为4.148（5）Å。制造蜂窝材料的一种选择是采用自发石墨化的概念：据预测，块体中已知的许多立方晶格更喜欢低于某个临界度的层状石墨结构，大约为十几ML。
2025-02-23 22:32:36
1380
原创 四声子计算软件FourPhono安装教程
FourPhonon的使用步骤和流程都与ShengBTE保持一致：用户只需要准备一个参数设定文件CONTROL，以及力常数文件（二阶力常数，三阶力常数，四阶力常数）。与ShengBTE类似，通过修改 arch.make进行编译，需要指定spglib的库文件libsymspg.a。与Thirdorder类似，参考shengbte中编译方法，在完成spglib编译后，修改setup.py。论文链接：https://doi.org/10.1016/j.cpc.2021.108179。获得ShengBTE文件。
2025-02-23 15:47:07
353
原创 Debye-Einstein-模型拟合比热容Python脚本
德拜模型将晶体中的原子振动视为连续弹性介质中传播的弹性波。固体的热容主要由低频的声学支声子贡献，存在截止频率，并未考虑光学支声子的贡献。在低温区与实验结果高度一致。爱因斯坦模型基于三种假设：1.晶格中的每一个原子都是三维量子谐振子；3.所有的原子都以相同的频率振动（与德拜模型不同）。固体比热模型中的德拜模型和爱因斯坦模型是固体物理学中用于估算固体热容的两种重要原子振动模型。在高温下，爱因斯坦模型与实验结果一致，特别是与杜隆-珀替定律相符。HC模型混合（将D和E模型填入）磁熵或相变熵值计算和统计。
2025-01-13 23:45:12
305
原创 BoltzGnu Boltztrap数据绘图脚本
需要Gnuplot版本大于5.0 ，依据File.trace & Doping.data 文件绘图。BoltzGnu包含四个Gnuplot脚本，允许绘制BoltzTraP输出数据。通过修改intrans文件并给出掺杂(空穴或电子)的浓度，案例如下。然后执行boltztrap程序可得到具有处女座效果的数据文件文件。
2025-01-08 22:06:58
195
原创 ubuntu安装vasp时gcc版本设置
当前超算服务器和个体服务器系统使用Ubuntu系统越来越多，尤其是在当前使用较多的系统CentOS7停止更新后，如果使用windows系统自带子系统进行计算练习时，ubuntu的安装与获取更加方便简单。在ubuntu初始环境下安装vasp需要配置足够合适的编译环境，对于gcc和g++编译器也是有着对于基于系统版本不同的使用方法。注：该例子使用ubuntu最高gcc仅为gcc-11,，24.0版本的ubuntu已经默认到了gcc-13，请视情况修改代码。注意请准确写入.h文件路径。同时更改默认版本gcc。
2025-01-05 03:08:31
494
原创 安装Intel-oneapi后 ifort : command not found，附2023版链接
笔者近期使用oneAPI2025版时发现fortran编译器执行文件ifort也被更迭，https://www.intel.com/content/www/us/en/docs/fortran-compiler/get-started-guide/2025-0/overview.html。链接: https://pan.baidu.com/s/1Mrd-hdaXDAaw3V3Lt9CWww?为了便于VASP等软件安装，建议使用oneAPI2023版本，即未对icc，icpc和ifort进行删改的版本。
2025-01-05 01:17:56
1703
3
原创 Phono3py hdf5文件数据读取与处理
在phonopy和phono3py的计算过程中，会将过程文件和结果文件写入到hdf5的格式文件中，以此来方便程序读入和写入，phonopy有关的诸多输出量可以通过额外设置获得文本文件，并用于数据处理，phono3py则不能自主输出文本数据。Phono3py是一个主要用python写的声子-声子相互作用相关性质的模拟包，可以基于有限位移算法实现三阶力常数和晶格热导率的计算过程，同时输出包括声速，格林奈森常数，声子寿命和累积晶格热导率等参量。本文则以此来介绍读取phono3py计算完晶格热导率后的数据读出。
2024-12-22 20:57:06
440
原创 Phono3py 计算累积晶格热导率及绘图
Phono3py是一个主要用python写的声子-声子相互作用相关性质的模拟包，可以基于有限位移算法实现三阶力常数和晶格热导率的计算过程，同时输出包括声速，格林奈森常数，声子寿命和累积晶格热导率等参量。本文以最新版phono3py（版本3.10）的example为基础，介绍如何基于力常数计算晶格热导率后，获得累积晶格热导率。有关计算结果会被写入到 kappa-m111111.hdf5中，通过读取该文件并生成累积晶格热导率。相关介绍和安装请参考往期推荐。以硅计算的文件为例。使用gnuplot绘制。
2024-12-22 20:55:56
569
原创 Anaconda3 pypi 清华大学TUNA镜像源使用帮助
Anaconda 是一个用于科学计算的 Python 发行版，支持 Linux, Mac, Windows, 包含了众多流行的科学计算、数据分析的 Python 包。Anaconda 安装包可以到以下链接下载。使用方法镜像站提供了 Anaconda 仓库与第三方源（conda-forge、msys2、pytorch 等，各镜像站镜像的第三方源并不相同，可以参考下方「第三方镜像源」一节）的镜像，各系统都可以通过修改用户目录下的.condarc文件来使用镜像站。
2024-12-20 21:10:26
786
1
原创 论文精读： NC Kagome ScV6Sn6竞争CDW起源
Nat Commun 15, 10428 (2024).https://doi.org/10.1038/s41467-024-54702-3摘要节选双层kagome金属ScV6Sn6中的竞争电荷密度波（CDW）序的起源仍未被充分理解。实验结果表明，冷却过程中，波矢量q2 =（ 1/3, 1/3, 1/2）的短程序，被长程q3 = （1/3, 1/3, 1/3）的 CDW序的在较低温度下抑制。然而，理论上，q2 CDW被预测为基态，这中现象难以被解释。本文利用非简谐
2024-12-20 01:07:04
913
原创 VASPKIT 功能及教程贴汇总
本文意在依据VASPKIT当前版本已有功能列表，按照功能模块编号顺序，整理相关已有教程，限于精力有限，只涵盖部分第一性原理计算相关公众号发表过的文章，更多详情请自行查阅和鉴别。您可以在接下来的几天内通过https://vaspkit.com/features.html获取该列表PDF版本。VASPKIT下载网址：https://sourceforge.net/projects/vaspkit/files/Binaries。VASPKIT官方论坛（中文）：http://bbs.vaspkit.cn。
2024-12-05 08:43:12
3239
原创 VASP结构优化专题贴汇总
本文意在归类汇总微信公众平台上已有较为细致的有关VASP结构优化相关内容的文章贴，部分取截图来突出，排序稍欠逻辑。以供 初学者 更为广泛和深入学习，请多加比较与甄别，取其精华。如有需要，可直接点开。
2024-12-03 09:21:04
592
原创 论文精读：NC VSe2/NbSe2异质结近藤晶格实验与理论计算
计算得到的2×√3和√7×√3 vdW结构的高度分别为0.61和0.62 nm，而插入的√3×√3 vdW结构的高度为0.58 nm，与实验测量结果完全吻合。近藤晶格起源于局部磁矩和传导电子的周期性阵列之间的杂化，导致局部磁矩的限制和伴随的低能电子有效质量的强烈增强，这构成了重费米子材料的基本理论模型。比较了有V原子插入和没有V原子插入的单层VSe2的形貌、CDW模式和光谱，未插层的VSe2表现为2×√3和√7 ×√3周期的CDW模式共存，与文献中报道的在HOPG底物上生长的相似。
2024-11-28 17:59:44
1415
原创 论文精读： PRB LiVSe2 Zigzag链序实验与理论计算
相反，zigzag链态在高压下稳定.。此外，虽然在环境压力下，V-V原子间距离表现出强烈的温度依赖性，但在压力下，8和300 K之间的位移大小没有差异，表明zigzag链在压力下迅速稳定。在30k下的LiVS2数据中，在~ 2 Å处出现了一个V-S距离对应的峰，在2.5 ~ 3 Å处出现了两个V-V距离对应的峰。在LiVO2和LiVS2的低温相中持续出现的钒三聚化在LiVSe2中消失，并且在LiVSe2中实现了规则的三角形晶格，但本研究表明，在LiVSe2中出现了有限相关长度的 Zigzag链局部畸变。
2024-11-20 23:54:32
642
原创 FCO-C6的光照可控Weyl节点和费米弧
本文基于第一性原理计算、基于最大局域化万尼尔函数的紧密结合模型和Floquet定理，研究了在线偏振光照射下，碳同素异形体C6结晶面心正交结构（FCO-C6）（具有两对Weyl节点的理想Weyl半金属）的光控演化过程。结果表明，光照射fco-C6的对称性可以被选择性地保持，这保证了光操纵的Weyl节点只能在动量空间的高对称性平面上运动。通过求解式(8)的零能模式，得到Weyl节点的位置为（+kx0, 0, kz1±）和（-kx0, 0, kz2±），其中。fco-C6光照射下（010）表面态和对应的费米面。
2024-11-16 17:24:34
188
原创 ATAT-mcsqs生成准随机结构(SQS)更新
而由于高熵合金内部多元、混合等特性，ATAT工具包内的mcsqs模块所涉及的特殊准随机结构Special quasi-random structure（SQS）于此有了非常大的应用需求。将vaspkit生成的原胞文件PRIMCELL.vasp复制为POSCAR，然后通过vaspkit 的414功能。先使用ATAT获得in文件的对称性信息，-2=x需要根据结构实际情况修改，过小会导致ATAT出现段错误。例如这里准备了一个面心立方的结构，组分单一，面心和顶角元素相同，这里设置为Si。
2024-11-12 21:30:21
1929
原创 论文精读：NC kagome FeGe 自旋声子耦合驱动CDW 实验与理论计算
基于QL = (0,0.5, 2.5)， Q∥M =(0,0.5, 3)和Q⊥A=(0,0,2.5)和（0,0,4.5）的温度相关超晶格峰，在TCDW （90 K）以下和（116 K）以上的θ - 2θ扫描以及穿过TCDW与温度相关的峰强度和峰宽度。A点和M点与温度有关的声子宽度。基于动量空间L=2.5平面上的CDW矢量Q =(0.5, -0.5, 2.5)，(-0.5,0.5,2.5)，0.5, 0,2.5)，（-0.5,0,2.5），(0, -0.5, 2.5)和（0,0.5,2.5）的HKL扫描。
2024-11-09 23:19:43
1159
原创 Kane-Mele X4Y2Z6材料自旋电子和谷电子理论研究
QVHS和SVHS分别表示量子谷霍尔态和自旋谷霍尔态。Ni4Hg2(TeSe)3 、 Pd4Hg2(TeSe)3 、Ni4Hg2(TeS)3 、 Pd4Hg2(TeS)3 、Ni4Hg2(SSe)3 和 Pd4Hg2(SSe)3 单层膜自旋分量S不同方向贡献。Pt4Hg2Se6, Pt4Hg2Se3(S/Te)3, Pt4Hg(Cd/Zn)Se6 和 Pt2(Pd/Ni)2Hg2Se6的晶体结构和电荷密度分布图。不同拉伸和压缩量下Pt4Hg2(SeTe)3和Pt4Hg2(SSe)3能谱的演变。
2024-11-05 21:53:43
643
原创 12种新型Full-Heusler化合物X2YTe （X = Na， K, Rb, Cs）热电第一性原理研究
值得注意的是，除Na2ZnTe和Na2HgTe外，其他材料均超过了长期存在的小于3的优值记录，其中Rb2ZnTe在700 K，n型掺杂浓度为2 × 10 19时，，优值达到8.11。12种Full-Heusler化合物在0-700 K温度范围内，计算了HA+3ph、HA+3,4ph、SCP+3ph和SCP+3,4ph四种模式下的温度依赖的𝜅L。12种Full-Heusler化合物3ph， 4ph和3ph+4ph的SRs在SCP模型下，在300 K下随Y原子的变化而变化。在T = 0和700k之间。
2024-11-04 21:35:48
924
原创 论文精读：PRL RuO2中不存在交错磁自旋分裂
压缩应变（-3.1%）和拉伸应变（+2.3%）分别沿着[001]和[110]方向产生,费米能级有轻微的变化，这与薄膜和单晶之间的ARPES结果中看到的整体能量转移相对应，突出了它们的电子结构的微小调整。三种不同动量下的S001极化与结合能（c）和自旋分辨能量分布曲线（EDCs）（e），外加S001极化在EF以下约0.1 eV处作为动量的函数（d）表明，体带是高度自旋极化的，其面内偏振矢量S001在Γ-Z-M'高对称平面上是反对称的。因此，这两种方法都为在体带α中观测到的自旋极化现象提供了一致和确证的解释。
2024-11-02 02:15:26
1232
原创 DFT新手教程：VASP中ISYM(对称性设置)取值设置
分析得到以下信息：晶胞的Bravais晶格类型，晶胞的点群对称和空间群(静态和动态)，如果晶胞为超胞则生成原胞，电荷对称化所需的所有周期性平移，使用自动k网格生成时，则是k点的对称不可约集合，标记对称等价离子。取而代之的是，电荷密度是通过对轨道应用相关的对称运算来构造的k-布里渊区不可约部分中的点。程序自动对称的识别标准：确定的空间群的总电荷密度，确定的空间群对离子施加的力，确定的空间群的应力张量。当打开对称性时，VASP确定POSCAR中的空间群和点群对称性，同时考虑在INCAR中通过。
2024-10-30 20:41:08
1657
原创 论文精读：PRL 交变磁MnTe中的手性分裂磁振子
通过加权最小二乘法拟合NPs中的谱，得到最佳解为J1= 3.99(3) meV, J2= −0.120(2) meV, J3= 0.472(3) meV，J10+J11= 0.0931(8) meV，D= 0.0482(5) meV。交替传播的磁振子色散也是交替的手性分裂。先前在MnTe上的INS实验受到相对宽松的能量分辨率的限制，只探测了第一、第二和第三近邻交换相互作用（记为J1、J2和J3），而忽略了磁振子色散的分裂特征。因此，即使在没有相对自旋轨道耦合的情况下，交变磁体中的电子带也是自旋分裂的。
2024-10-20 02:52:29
1162
原创 论文精读：拓扑超导体PdBi2Te4和PdBi2Te5计算
另一方面，也可以构建一个将SOC从非SOC第一原理计算中排除的slab Hamiltonian H_nsoc_slab，并手动添加H_SOC和H_surf来模拟拓扑电子态和TSC中的可变SOC和表面效应，其中H_surf将包括表面电位差、上下表面杂化等。本文与前人对Class D TSC4的研究平行，开发了一个程序，从从头计算的角度模拟二维平板系统的超导性质，并表征其拓扑不变性，其中包括基于SCTM实现手性TSC的必要成分，如体带结构、超导间隙、Zeeman分裂、Rashba自旋轨道耦合和化学势。
2024-09-30 02:18:31
1136
原创 VESTA 可视化声子谱声子振动模式
这里我们介绍一个简易python脚本，用方便处理和生成可用于VESTA打开的文件，导入到VESTA中便可以显示对应q点相应频率的声子模式下的个原子振动图像。phonopy在处理原子间力常数和声子色散后可以导出部分位置的声子传输矢量和模型中各原子的振动模式，但文本内容复杂并不能很好查看。将扩包前的POSCAR导入到VESTA中，左上角点击file—save，保存为vesta格式文件，用于脚本的输入文件。根据声子频率的序号或者频率值来选择需要绘图的文件。在声子谱的设置文件中，加入。
2024-09-13 18:25:12
1785
原创 晶体晶格热导率的快速可解公式
a，格内森参数γ，热容量CV，晶格常数a，原子半径ra, rb，原子质量ma, mb，电子亲和度Eea, Eeb，每原子基态能Ea0, Eb0，价电子数V Ea, V Eb，电负性χa， χb，原子序数Za, Zb。在这项工作中，将Slack公式简化为更精确和紧凑的形式，并将公式的计算结果与AFLOW数据库、Slack公式和密度泛函理论(DFT)进行了比较，以验证公式的准确性。在这项研究中，引入了一个公式，可以实现高精度(平均相对误差=8.97%)，并提供快速预测，可在广泛的无机二元和三元材料中获得κL。
2024-09-11 08:58:28
1067
原创 Mn基kagome平带和铁磁波动实验和理论验证
对150 K以上的χ(T)−χ0进一步进行居里-魏斯(CW)分析，得到有效磁矩μab eff= 0.86 μB/Mn，μc eff=0.87μB/Mn，CW温度θab CW = 421 K，θc CW = 368 K。这种行为与观察到的ρ(T)在25 K以下的上升和χ(T)的幂律增长是一致的。从M(H, T)中减去线性项，并绘制磁化曲线Mdiff = M(H, T)−M(H, T = 19 K)的结果差。随着温度的降低，χ(T)急剧增加，没有饱和或异常的迹象，从而排除了远程磁有序的发生。
2024-09-09 15:31:06
1117
原创 论文精读：Dirac半金属反常能斯特效应设计
为了在费米能级上获得增强的ANC值，需要AHC曲线是一个奇函数，在狭窄的能量区间内，两个相反的峰位于费米能级的上方和下方。在这项工作中，从理论上证明了一对Dirac节点在Zeeman场下表现出奇数分布，双峰异常霍尔电导率曲线，其化学势和补偿载流子特征，导致与具有两个Weyl节点的简单Weyl半金属相比，ANC增强。在单峰Weyl模型中，ANC曲线在费米能级为零处接近零，而在双峰Dirac模型中，ANC曲线恰好在费米能级处显示出最大值，几乎是Weyl模型最大值的300%。
2024-09-03 22:59:00
1266
原创 论文精读：JACS —— Sb2Si2Te6与Sc2Si2Te6热电性能
在简并导带最小值处(Σ/ Σ ' /Σ″)，Sc dxy、dyz、dxz和dx 2 - y 2(但没有dz 2)轨道的贡献是相等的，还有一些Si pz的贡献，形成了杂化轨道状态。与Sb2Si2Te6相比，Sc2Si2Te6表现出更高的能带简并度(n型和p型分别为6和2，而Sb2Si2Te6中的两种载流子类型分别为1)和载流子有效质量，不同的颜色表示载波群速度的大小。对应于Sc2Si2Te6中Σ(绿色)，Σ'(红色)和Σ''(蓝色)k点的实空间方向(即6倍简并导带谷的位置)，沿着层(a)和垂直于层(b)。
2024-08-31 14:12:06
1010
原创 Doped code 介绍
再现性、制表、自动化兼容性/健全性检查、应变/位移分析、浅层缺陷/特征值分析、高通量兼容性、Wyckoff 分析…doped是一款Python软件，用于缺陷超单元计算的生成、前/后处理和分析，以高效、可重复、用户友好、功能强大且完全可定制的方式实施缺陷模拟工作流程。缺陷分析:自动解析计算输出，以计算缺陷形成能量、有限尺寸修正(FNV & eFNV)、对称性、简并度、跃迁能级等。绘图:生成缺陷形成能量、化学势极限、缺陷/载流子浓度、费米能级、电荷校正等出版物质量的图表。教程页面提供了演示代码功能和用法，
2024-08-31 14:10:24
567
原创 ThermoParser 介绍
ThermoParser是一个工具包，用于简化专业材料科学代码产生的数据分析，以热电学为中心，但也适用于任何与电子和/或声子传输有关的内容。ThermoParser是一个Python库，它包含数据检索、操作和绘图的函数，只需几行代码就可以轻松地使用它，生成各种高质量的绘图。ThermoParser还包含一套命令行工具，可以在一个命令中检索特定数据、保存派生属性和绘制图形。
2024-08-31 14:07:34
235
原创 vaspup2.0介绍
(光学/高频介电常数)相对于n带宽，使用furthüller等人的方法计算(LOPTICS = True)收敛测试。(离子对静态介电常数的贡献)关于ENCUT和k点密度，用密度泛函微扰理论(DFPT)计算收敛测试。最初的vaspup是由Alex Ganose开发，用于基态能量收敛测试和POTCAR生成。实时软件库:https://github.com/kavanase/vaspup2.0。基态能量相对于ENCUT和k点密度的收敛测试。原子力相对于ENCUT和k点密度的收敛测试。
2024-08-31 14:05:03
343
原创 BoltzTraP2稳定安装教程
BoltzTraP2是形成原始和广泛使用的BoltzTraP代码的基础的电子波段平滑傅立叶插值算法的现代实现。然而，许多其他的特征是可用的，包括基于重建带的费米曲面的三维图。建议使用conda建立一个新的安装环境，以此来稳定安装和使用BoltzTraP2。但是实际上因为安装环境的错综复杂，以及软件依赖较多，会出现各种各样的问题。BoltzTraP2官方安装安装教程较为简单，通过pip可以一键安装，激活btp2环境后，通过pip安装boltztrap2的依赖包。然后可再通过pip安装boltztrap2。
2024-07-16 09:53:10
1118
原创 VASP-DFPT计算报错SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG
笔者以往写过一则SETUP_DEG_CLUSTERS: NB_TOT exceeds NMAX_DEG报错提示的解决办法，即弃用密度泛函微扰，改用有限位移方法。使用密度泛函微扰方法计算材料力常数时，由于使用超胞结构和高精度的设置等等，计算时会产生诸多报错，其中。修改第一处NMAX_DEG的赋值，默认数值为48，我这里修改为480。这次分享一则通过修改和重新编译vasp的方法，来解决这一问题。在已经编译好的vasp安装 目录中，找到src文件夹下的。文件相关，不会影响vasp的正常使用。
2024-06-04 10:23:12
741
原创 VASP_AIMD+VASPKIT计算含温力学性质
材料的力学性质一直是DFT计算的重要方向，笔者在以往已有针对于静态结构的力学性质诸如弹性常数的相关计算，同时可通过VASPKIT借助相关方程导出力学性能。手动生成VASPKIT计算所需要的设置文件INPUT.in，内容如下，参数详细说明附在后面，可根据需要自行更改。本文依据VASPKIT的EXAMPLE，简单介绍使用AIMD方法和VASPKIT软件计算含温力学性能。在计算AIMD时首先进行扩胞，这里使用2×2×2的超胞结构进行下一步计算。准备好VASP的输入文件后可运行VASPKIT来生成计算文件。
2024-05-09 11:15:59
1844
3
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人