当前位置:网站首页>On the problems of born charge and non analytical correction in phonon and heat transport calculations
On the problems of born charge and non analytical correction in phonon and heat transport calculations
2022-07-06 22:35:00 【odin_ linux】
When calculating the phonon structure or heat transport properties of polar materials , In order to describe accurately in dispersion gamma point LO-TO The classification of , In the process of solving the dynamic matrix , Add the correction term of polar charge action to the dynamic matrix , As in many programs involving phonon calculation phononpy, alamode as well as shengbte And so on non-analytic correct , To achieve this correction , It is necessary to give the BORN Effective charge , That is, the tensor form of the equivalent charge in motion , Each atom is 3*3 Tensor , And dielectric constant ε tensor , The specific dynamic matrix is modified as follows :
among ε by BORN Dielectric constant in the document ,Z Is the effective charge tensor of the corresponding atom . This BORN The document must pass DFPT Perturbation method calculation , Can pass QE or VASP Realization .
But there is usually a problem here , It involves the conversion of units in the implementation of the above formula , such as phononpy in , In the production BORN file , The first line is a unit Converted values , Yes VASP yes 14.399652 ( namely Ha*bohrbohr To eV*ang), What is actually converted here is the unit of dielectric tensor ,born The charge is in e In units of , The international unit of dielectric constant is F/m, Atomic units are usually used in material calculations , namely 1/(eV*ang), In the correction of the above expression is 1/ε Dimension is eV*ang, So when phonon calculation defaults eV and ang In units , Dielectric function ( speculation VASP Use in the result Ha and Bohr) Put it into the above formula , If the above formula is used by default (ev*ang), At this point, you need to do unit conversion , Equivalent to the above formula, you need to multiply by a conversion coefficient (2*13.6*0.529177), just as phononpt What is needed in .
This also shows that VASP The calculated dielectric unit is (1/(Ha*Bohr).
But yes. QE,phononpy Targeted at QE The force constant unit of is also adjusted , So the final conversion coefficient is only 2, Equivalent to Ha To Ry. This also shows that QE The unit of output dielectric constant and VASP identical , however QE The internal unit is usually Ry and Bohr, Aimed at GaN, Use them separately VASP and QE Calculate dielectric and BORN Effective charge , The result is :
VASP: Dielectric xx The value is 5.5429220 QE: 5.82133
VASP: Ga Atomic effective charge xx value 2.5749225 QE: 2.570464
The comparison shows that , Considering the difference between lattice point and pseudopotential in calculation , The difference is very small , It is speculated that the unit of output dielectric of the two is the same , All for (1/(Ha*Bohr)).
边栏推荐
- Mise en place d'un environnement de développement OP - tee basé sur qemuv8
- [IELTS speaking] Anna's oral learning record part1
- Use ECs to set up an agent
- The difference between enumeration and define macro
- UE4蓝图学习篇(四)--流程控制ForLoop和WhileLoop
- Improving Multimodal Accuracy Through Modality Pre-training and Attention
- (十八)LCD1602实验
- Aardio - does not declare the method of directly passing float values
- POJ 1258 Agri-Net
- config:invalid signature 解决办法和问题排查详解
猜你喜欢
如何用程序确认当前系统的存储模式?
[linear algebra] determinant of order 1.3 n
How to confirm the storage mode of the current system by program?
Netxpert xg2 helps you solve the problem of "Cabling installation and maintenance"
Aardio - construct a multi button component with customplus library +plus
Aardio - does not declare the method of directly passing float values
CocosCreator+TypeScripts自己写一个对象池
NPDP认证|产品经理如何跨职能/跨团队沟通?
Export MySQL table data in pure mode
Clip +json parsing converts the sound in the video into text
随机推荐
MATLAB小技巧(27)灰色预测
视图(view)
剪映+json解析将视频中的声音转换成文本
手写ABA遇到的坑
如何实现文字动画效果
Daily question 1: force deduction: 225: realize stack with queue
The SQL response is slow. What are your troubleshooting ideas?
How to use flexible arrays?
枚举与#define 宏的区别
Advantages of link local address in IPv6
Crawler obtains real estate data
Build op-tee development environment based on qemuv8
three.js绚烂的气泡效果
memcached
How do I write Flask's excellent debug log message to a file in production?
0 basic learning C language - digital tube
The difference between enumeration and define macro
树的先序中序后序遍历
Aardio - 通过变量名将变量值整合到一串文本中
Attack and defense world miscall