phonopy计算声子谱时考虑LO-TO劈裂

2023-10-17 10:40

本文主要是介绍phonopy计算声子谱时考虑LO-TO劈裂,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

经验分享

在使用VASP+加phonopy计算材料声子色散结构时,可以考虑激活非分析项矫正,这时需要增加计算的工作量以获得BORN文件。以此计算获得的声子色散则会考虑LO-TO的劈裂效果。

在phonopy的手册中提到了可以通过自洽计算获得BORN,手册中推荐的INCAR为

  • PREC = Accurate

  • IBRION = -1

  • NELMIN = 5

  • ENCUT = 500

  • EDIFF = 1.000000e-08

  • ISMEAR = 0

  • SIGMA = 1.000000e-02

  • IALGO = 38

  • LREAL = .FALSE.

  • LWAVE = .FALSE.

  • LCHARG = .FALSE.

  • LEPSILON = .TRUE.

这是在使用有限位移法也就是冻声子法时,增加的额外的计算任务。当使用DFPT方法计算原子受力时,可直接在INCAR中添加下面这条tag,然后继续计算

  • LEPSILON = .TRUE.

在vasp的计算任务完成后,通过phonopy的代码可直接获得BORN文件,命令如下

  • outcar-born >BORN #旧版

  • phonopy-vasp-born>BORN #新版

在获得BORN文件后,在处理原子间力常数获得声子色散时的命令中加入--nac即可,如手册中的命令案例

  • phonopy -p --nac band.conf

在phonopy手册中给出考虑和不考虑LO-TO劈裂的NaCl的声子谱的对比,分别如下:

不考虑LO-TO劈裂

考虑LO-TO劈裂

愿有所成

IEchoQ

引喻失义 妄自菲薄

这篇关于phonopy计算声子谱时考虑LO-TO劈裂的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/224823

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

分布式系统的主要考虑

异构性:分布式系统由于基于不同的网路、操作系统、计算机硬件和编程语言来构造,必须要考虑一种通用的网络通讯协议来屏蔽异构系统之间的禅意。一般交由中间件来处理这些差异。缺乏全球时钟:在程序需要协作时,它们通过交换消息来协调它们的动作。紧密的协调经常依赖于对程序动作发生时间的共识,但是,实际上网络上计算机同步时钟的准确性受到极大的限制,即没有一个正确时间的全局概念。这是通过网络发送消息作为唯一的通信方式

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

GPU 计算 CMPS224 2021 学习笔记 02

并行类型 (1)任务并行 (2)数据并行 CPU & GPU CPU和GPU拥有相互独立的内存空间,需要在两者之间相互传输数据。 (1)分配GPU内存 (2)将CPU上的数据复制到GPU上 (3)在GPU上对数据进行计算操作 (4)将计算结果从GPU复制到CPU上 (5)释放GPU内存 CUDA内存管理API (1)分配内存 cudaErro

Java - BigDecimal 计算分位(百分位)

日常开发中,如果使用数据库来直接查询一组数据的分位数,就比较简单,直接使用对应的函数就可以了,例如:         PERCENT_RANK() OVER(PARTITION BY 分组列名 ORDER BY 目标列名) AS 目标列名_分位数         如果是需要在代码逻辑部分进行分位数的计算,就需要我们自己写一个工具类来支持计算了 import static ja

OpenStack离线Train版安装系列—2计算节点-环境准备

本系列文章包含从OpenStack离线源制作到完成OpenStack安装的全部过程。 在本系列教程中使用的OpenStack的安装版本为第20个版本Train(简称T版本),2020年5月13日,OpenStack社区发布了第21个版本Ussuri(简称U版本)。 OpenStack部署系列文章 OpenStack Victoria版 安装部署系列教程 OpenStack Ussuri版