GULP采用经典势能模型输出shengBTE计算用力常数文件

墨衍会员 · AI 创作全网分发 墨衍智能分发 本文由作者通过墨衍一键同步至各平台 1分发平台 400+累计阅读 472平均单平台已同步平台优快云微信公众号微博知乎掘金百家号博客园抖音小红书 你的文章也可以这样分发 墨衍会员 ¥399起/年，写一次发全网我也要 → 码龄8年

关注

原创最新推荐文章于 2026-05-19 11:23:42 发布 472 阅读 AI 写同款· GEO 优化 ›

Fourthorder:与shengBTE中的thirdorder.py类似，开发了Fourthorder.py以计算四阶原子间力常数（4th-IFC）。请在FourPhonon存储库的“讨论”部分中发布您的问题四阶 Zherui韩枭阳吴莉天利锋秀林阮四阶脚本源自三阶脚本，用于帮助[包含四声子模块的ShengBTE]用户以高效便捷的方式创建FORCE_CONSTANTS_4TH文件。更具体地说，它执行两项任务：它解析出一组不可约的原子位移，从中可以计算出完整的非谐四阶原子间力常数（IFC）矩阵。置换后的超级单元将保存到输入文件中，该文件可输入到第一性原理DFT代码中，以计算原子位移产生的力。当前支持的DFT代码是VASP（Fourthorder_vasp.py）和Quantum ESPRESSO（Fourthorder_espresso.py）。从DFT代码创建的输出文件中，Fourthorder重构完整的四阶IFC矩阵，并将其以正确的格式写入FORCE_CONSTANTS_4TH。汇编四阶是一组Python脚本。它是使用Python 2.7.3开发的，但应该可以使用稍旧的立即下载

仅需在option参数中使用supercell控制2阶力常数的扩胞倍数，以及output shengBTE:

但是注意，gulp这里存在一个bug，输出的2阶力常数文件中超胞中原子次序与phonopy以及shengBTE用不一致：shengBTE中FORCE_CONSTANT中原子次序为：ucatom*uz*uy*ux；但是gulp中为ucatom*ux*uy*uz，对于各向同性系统，则计算没有区别，但是对于各向异性系统，会产生计算错误，需要转换2阶力常数的中的原子次序。

收录于

当前文章被以下社区和专栏收录：

标签

#gulp #node.js #前端

sousaw-ShengBTE-Kit.zip ShengBTE 立即下载

相关推荐

【AMD显卡在WIndow10/11部署带GPU支持的深度学习环境(Pytorch-DirectML篇)】

这篇是之前Tensorflow-DirectML的姊妹篇，如果对TF感兴趣的可以搜索我的博文。首先我们丢一个官网地址，但是官网那个会误导新人如果你按那个操作的话。 https://docs.microsoft.com/en-us/windows/ai/directml/gpu-pytorch-windows 这里用我的方法操作：首先声明一下，因为我的显卡是AMD RX6600，不带XT(性能和RTX3060 12G版本半斤八两)，所以目前没有最新的WSL2驱动，如果有了的话欢迎告知我，所以我这里就.

weixin_44029053的博客 2万+

ShengBTE简单介绍、安装与使用

在ShengBTE和thirdorder安装文件中均带有其“使用说明书”，且都命名为README.md，里面有一些基本的安装、使用介绍，尤其是一些参数的介绍，初学者在使用之前要认真阅读并理解，对以后遇到问题有很大帮助，这也是对于软件的最基础的了解。这里主要是补充了我们make编译出的libsymspg.a文件的路径，以及更改了编译器为mpiifort，在下方也调整了MKL的配置，大家在安装的时候可以复制上方的内容，手动更改libsymspg.a文件的路径，然后可以尝试进行安装。

DFT_M的博客 7456

R语言vegan包实战：用varpart函数搞定生态学VPA分析（附完整代码与数据集）

本文详细介绍了如何使用R语言vegan包中的varpart函数进行生态学方差分解分析（VPA），包括数据准备、模型构建、结果解读和可视化呈现。通过完整的代码示例和数据集，帮助研究者理解不同环境因子对生物群落变化的相对贡献，并掌握显著性检验等高级技巧。

weixin_42525005的博客 447

shengBTE计算热导率数值不对的可能原因

如使用cp2k计算力常数，通常使用QE等平面波软件，对于超胞体系，随着扩胞，在较大的截断半径下，所需计算量非常大，使用CP2K中的QS模块也能进行第一性原理计算，具有相对较高的计算效率，但是CP2K不同于平面波的QE等，计算中仅有gamma点计算，需要采用实空间较大超胞，如果使用和QE等相同的超胞，则会导致原子受力精度不高，一般在~0.0003Rh/bohr^2量级的误差，也会对某些体系，特别是二维体系的热导率有很大的影响。2.shengBTE的control文件中晶格参数单位默认nm，有转换系数的参数；

odin_linux的博客 3059

山东大学软件学院计算机组成原理课程设计实验二

文章目录一、实验时间二、实验内容1.设计异步模16加1计数器2.设计4位移位器3.超前进位16位ALU的设计一、实验时间 2021年4月1日星期四，第五周二、实验内容 1.设计异步模16加1计数器基本功能：使用D触发器（dff元件）和异或门构建从0开始计数的模16加1计数器对应的电路图（这里已经实现了预置初值的功能）拓展功能：（1）实现模256加1计数器（2）实现计数器预置初值的功能这里复用了上面的模16加1计数器这里异步触发器预置初值的功能是这样实现的：（1）首先将CLR端调为0，

叶卡捷琳堡的博客 5250

Win10 LTSC 2021 CPU占用率飙升？KB5017308补丁修复全攻略

本文针对Windows 10 LTSC 2021系统中因wsappx进程导致的CPU占用率异常飙升问题，提供了详细的解决方案。核心在于安装微软官方发布的KB5017308补丁，该补丁能从根本上修复系统组件缺陷，使CPU占用恢复正常。文章详解了补丁的下载、安装步骤、注意事项及安装后的验证方法，是解决此问题的权威指南。

python9snake的博客 1262

基于VASP+phonopy+shengbte计算声子相干的热学性质

声子相干的热学性质在计算材料领域有着非常重要的地位，我们可以通过第一性原理计算得到任一种晶体的热力学性质，但是很多小伙伴尤其是初学者还比较陌生，本文介绍基于DFT计算与声子相关的性质，是用VASP+phonopy+shengbte计算一种材料的声子谱，声子态密度，热力学性质（包括准简谐近似下的亥姆霍子自由能，等压热容Cp，热膨胀系数，格林爱森参数，和考虑高阶声子的声子散射，群速度，声子自由程，热导率等等）共分为如下几个部分： (1). 生成4个输入文件: POSCAR POTCAR INCAR..

320科技工作室 1万+

YOLOv11改进 - C3k2融合 | C3k2融合CBSA 收缩 - 广播自注意力:轻量级设计实现高效特征压缩，优化处理效率 | NeurIPS 2025

本文介绍收缩 - 广播自注意力（CBSA）（（Contract-and-Broadcast Self-Attention））机制，并将其集成到YOLOv11中。传统注意力机制存在黑盒难理解、计算复杂度高的问题，CBSA通过算法展开推导出本质上可解释且高效的注意力机制。它先从输入数据中选出少量代表性tokens，接着对代表进行收缩计算，再将结果广播给所有原始数据。该机制计算量线性增长，具有明确数学解释，还能统一多种注意力机制。我们将CBSA代码集成到YOLOv11的C3k2模块中。实验证明，改进后的YOLOv

魔改工程师的博客 852

FORCE_CONSTANTS中3阶力常数大小与原子间距的分析脚本

shengBTE所需的3阶力常数文件格式为FORCE_CONSTANTS_3RD, 想要得知所有力常数数值与原子间距的变化关系，或者输出更小截断半径下新的FORCE_CONSTANTS_3RD以验证热导率随阶段半径的变化，需要分析FORCE_CONSTANTS_3RD内数据，这里提供一个PYTHON脚本，能够读取FORCE_CONSTANTS_3RD文件，并绘制力常数强度和原子间距的散点图，以及输出新的阶段半径下的力常数文件，供大家参考，注意python的缩进。 #####################

odin_linux的博客 2838

大模型部署指南：SGLang 框架

本文介绍使用 SGLang 框架部署大模型

举世誉之而不加劝，举世非之而不加沮，定乎内外之分，辩乎荣辱之境，斯已矣。 1611

【Conda虚拟环境】新建 + 激活 + 安装 + 查看 + 退出 + 删除 + 复制 + 导出 + 导入 + 更改

（1）创建虚拟环境（2）激活虚拟环境（4）查看虚拟环境（3）退出虚拟环境（5）删除虚拟环境（6）复制虚拟环境（7）导出虚拟环境（8）导入虚拟环境（9）更改虚拟环境的安装路径（10）在指定路径下，配置虚拟环境

📷 图像处理践行者 | 传统算法 · 深度学习 · 多模态 | 分享实战经验与技术洞见，探索从像素到智能的无限可能。 4万+

热输运物性的理论计算方法

在固体物理中，除了大家广泛关注的电子结构和光电性质等物理性质外，晶格中热流输运和调控也是至关重要的一个方面，这不仅影响微纳器件的散热问题，同时在热电等领域的应用中，较低的材料热导率是影响热电转换效率的重要因素。这里简单汇总目前能够开展热输运计算的程序平台、原理和使用方法： 1.shengBTE 基于第一性原理的声子玻尔兹曼输运方程法，本质是基于玻尔兹曼输运方程，这个方程其实是空间热流守恒的微分方程，其中微分点处的热流会聚等于该处声子散射等造成的热流消减。求解这个方程，需要知道声子散射等的寿命，基于三声子

odin_linux的博客 2867

MuJoCo 入门教程（四）建模

在 MJCF 中，一个体可以包含多个关节，因此无需引入虚拟体来创建复合关节。相反，只需在同一主体中定义构成所需复合关节的所有原始关节即可。例如，可以使用两个滑块和一个铰链来模拟在平面内运动的肢体。其他 MJCF 元素可以在嵌套的主体元素所创建的树中定义，特别是关节、几何体、场地、摄像机、灯光。当一个元素被定义在一个肢体内时，它就固定在该肢体的局部坐标系上，并始终随其移动。涉及多个机构或根本不涉及机构的元素，将在运动树之外的单独部分中定义。

机器人控制与仿真应用 4162

GULP不同版本导出ShengBTE用力常数文件，存在一些问题

2. 但是，GULP-6版本中导出的2阶力常数文件正确，3阶力常数文件存在错误，表现在三阶力常数文件中对超胞内相邻胞的位置x和y存在颠倒，反倒GULP-5版本中的是正确的。3. 可以自己通过thirdorder.py来计算，这时候，gulp输出中力的单位为eV, 对正交晶格，需要有一个转换系数，-1*输出力（eV）/相应方向的晶格模量。1. GULP-5版本，导出的2阶力常数文件，存在原子扩胞后超胞内原子排序和phononpy等的不一致，从而导致热导率计算问题，这个在GULP-6版本中已经纠正；

odin_linux的博客 735

visio安装包

Microsoft Office Visio 是一款专业的制图软件。利用经过更新的形状、协作工具和数据链接图表，创建专业级图表以简化复杂的信息。

华为的C语言编程规范.pdf

华为技术有限公司内部技术规范。本规范制定了编写C语言程序的基本原则、规则和建议。从代码的清晰、简洁、可测试、安全、程序效率、可移植各个方面对C语言编程作出了具体指导。

ShengBTE安装processes. f90报错的解决办法

而我早起安装教程是依据旧版，在2022年1月份完成，所以依据我去年的教程安装最新版的软件会出现一些原本不该发生的问题，但笔者当时已经说明了具体修改的思路，如果按照此思路来的话也是可以避免出错的（甩锅）。在arch.make.example文件中软件作者更新了一些参数，尤为重要的是在FFLAGS后添加了一个-qonemp，而参照我的教程的小伙伴不少将其去掉，软件作者在今年有对软件进行更新。在软件安装的过程中可在arch.make保留此参数并进行下一步，即可使软件正常安装。引喻失义妄自菲薄。

DFT_M的博客 631

Devexpress16.2.6破解工具

Devexpress16.2.6破解工具Devexpress16.2.6破解工具Devexpress16.2.6破解工具Devexpress16.2.6破解工具

bash+vasp+ShengBTE自动计算材料热运输性质脚本

之前我们已经详细介绍了ShengBTE软件与其安装教程，详细请查看ShengBTE简单介绍、安装与使用由于在计算三阶力学常数时，根据扩胞比例以及近邻原子受力情况的不同，生成的3RD.POSCAR数量也尽不相同，笔者最近一次的计算总共计算了700+个任务，如果在此过程中手动布置任务则费时费力，效率低下。而在后续的ShengBTE计算时，通过bash脚本进行任务布置也可以实现多任务顺序计算，节省精力与时间。本文针对材料的热输运性质计算主要参考天玑算所提供的教程，详细请见https://phadwiki.com/

DFT_M的博客 3670

精雕软件(JDpaint) 5.5 绿色版

精雕软件(JDpaint) 5.5 绿色版精雕软件(JDpaint) 5.5 绿色版精雕软件(JDpaint) 5.5 绿色版

VASP结合vaspkit+ShengBTE计算热电优值（一）

使用 vaspkit-681命令生成高密度的KPOINTS，然后进行静态计算（注意只有使用这项功能生成KPOINTS计算的结果才能继续使用vaspkit命令计算下一步，使用M-S方法自动生成K点的计算结果无法进行下一步）。这一步一般需要非常长的时间去进行计算，因此扩胞的大小可以稍微比二阶力常数矩阵的计算时的大小小一点。其中Phonopy用于计算声子谱及二阶力常数，Thirdorder用于计算三阶力常数，ShengBTE用于结合前面两者的结果计算晶格热导率。准备好计算的材料对应的POSCAR。

320科技工作室 3003

VASP结合vaspkit+ShengBTE计算热电优值（二）

注意BTE.kappa_tensor文件中，每列对应的方向顺序分别为XX,XY,XZ,YX,YY,YZ,ZX,ZY,ZZ.且最后一行为有效的收敛值，使用时只需要取最后一行。这两种方法的计算结果都是可用的。原子层厚度=POSCAR 真空层方向相距最远的两个原子的范德华半径之和+POSCAR 真空层方向相距最远的两个原子的距离。目前有两种方法计算电子热导率，一种是使用BoltzTraP计算得到的电子热导减去TσS2。因此要计算ZT值，我们需要分别计算电导率σ，塞贝克系数S，电子热导率κe，晶格热导率κc。

320科技工作室 2459

vasp手册_VASP计算笔记-ShengBTE计算热导率流程(1)

注：对计算结果的准确性不做任何保证。对于参数的设置没怎么验证，可能会有些错误。主要使用ShengBTE、thirdorder(两者在ShengBTE网站下载) 和phonopy软件。ShengBTE计算热导率前所需要准备的文件FORCE_CONSTANTS_2ND #见二阶力常数计算 FORCE_CONSTANTS_3RD #见三阶力常数计算注：基于玻尔兹曼输运方程计算晶格热导率的软件一般都需要计...

weixin_39948439的博客 6196

ShengBTE自动生成CONTROL文件和数据处理脚本

将自动读取CONTROL文件中设置的温度，并扫描当前文件夹中是否有未设置或者其他计算的文件结果。每个温度的处理结果将放置在result文件夹中。数据处理方法，直接python执行。准备文件：力常数文件和结构文件。产生的CONTROL文件。

DFT_M的博客 437

使用GULP(General Utility Lattice Program)计算声子结构和晶格动力学下热导率

GULP(General Utility Lattice Program)基于经典的势能模型，能够以此计算材料的一系列物性，包括晶格优化，原子位置优化，经典势能模型参数拟合，声子计算，缺陷计算等。这里主要介绍

odin_linux的博客 4111

GULP案例3：如何利用经验势函数计算弹性常数？

弹性常数描述了晶体对外加应变的响应的刚度。在材料的微小线性变形范围内，体系的应力与应变满足胡克定律。对于足够小的变形，应力与应变成正比，三维材料的弹性刚度常数矩阵是6×6：以六角晶系为例进行说明，其独立的弹性常数有C11，C12，C13，C33和C44：复现内容： ZnO-WZ结构的弹性常数和模量期刊名字： Journal of Physical Chemistry C 影响因子（IF）：4.126 中科院分区：2区文献题目： New ab initio based

nmrwxw的博客 1717

VASP+phonopy计算声子群速与数据导出

通常在声子计算中声子群速等也具有一定的分析意义，笔者在翻阅一些相关内容的经验贴上并结合phonopymanual内容总结了一下有关声子群速的计算，并借助vaspkit导出有关声速的数据。而且具体使用过程中，关于执行命令还有执行文件的tag内容还有一些可替换的关系,--gv,--group_velocity可在执行计算命令时替代执行文件中的GROUP_VELOCITY=.TRUE.内容，--gv_delta_q可在执行计算命令时替代执行文件中的(GV_DELTA_Q)内容。...

DFT_M的博客 4192

lapack安装_Linux安装ShengBTE(以本人超算为例）

第一步：去此网址下载ShengBTE(https://bitbucket.org/sousaw/shengbte)上传至安装目录:/lustre/ISSP2/xm/software;安装好spglib，比如我的目录：/lustre/ISSP2/xm/software/sousaw-thirdorder-b6962b2088ec/spglib第二步：根据arch.make.example编辑arch...

weixin_39625747的博客 1254

上一篇：关于声子和热输运计算中BORN电荷和non-analytic修正的问题

下一篇： shengBTE计算热导率数值不对的可能原因

【AMD显卡在WIndow10/11部署带GPU支持的深度学习环境(Pytorch-DirectML篇)】

weixin_44029053的博客

11-28

2万+