| 查看: 2940 | 回复: 17 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】关于QHA模块的使用已有4人参与
|
|||
|
以AlAs为例，想请教一下具体计算不同温度的热力学性质怎么实现。 以下是espresso-4.2中给出的用QHA计算alas例子中的Edit_Me文件 ############################################################################## # Optional parameters, any information specific for the system studied # SysInfo='AlAs' # Mandatory parameters # Specify SystemName and Force Constants matrix Sysname='AlAs' FC_file='alas444.fc' # # Specify lattice type (used to create ttrinp file). It should be the same as in scf.in file # Specify atoms in the unit cell as they specified in scf.in file # Specify atomic masses for these atoms in the same order as in scf.in # Specify the frequency step (delta_e) as well, but 0.75 is a good choice ibrav=2 atoms="Al As " mass="26.98 74.922 " delta_e=0.75 # Edit ONLY amass parameters # Please do not change flfrq='frequency' line # leave asr (acoustic sum rule) and flfrc lines cat >matdyn.init < amass(1)=26.98, amass(2)=74.922, asr='crystal', flfrc=$FC_file, flfrq='frequency' &end EOF # # In most cases there is no need to edit files listed below, but if you like ... # # Temperature range for thermodynamic properties # T_start, T_end, T_step for QHA calculations cat > Temperature < EOF # Debye Temperature calculations # Phonon DOS filename (total phonon DOS, not projected), leave it as PHDOS.out # accuracy (limited 1.d-5, more accuracy is not required ) # Low_Temp_start, Low_Temp_end, and Low_Temp_step for Low Temperature limit, up to 15-30K # Hihg temperature and T_step for HT limit cat >T_Debye.in < 0.0001 3 15 3 500 10 EOF 具体实施的时候，我的理解是这样的，不知道对不对，想请高人指教一下，感激不尽！ （1）先构建alas.scf.in文件，进行pw.x的自洽计算； （2）构建alas.ph.in，根据（1）的相关输出进行ph.x计算，得到力常数和动力学矩阵； （3）用QHA计算，但这里QHA计算的in文件的写法是怎样的呢？我的理解是如下的，但不知对否（里面的EOF和cat > Temperature < alas.qha.in文件 SysInfo='AlAs' Sysname='AlAs' FC_file='alas444.fc' ibrav=2 atoms="Al As " mass="26.98 74.922 " delta_e=0.75 cat >matdyn.init < amass(1)=26.98, amass(2)=74.922, asr='crystal', flfrc=$FC_file, flfrq='frequency' &end EOF cat > Temperature < EOF cat >T_Debye.in < 0.0001 3 15 3 500 10 EOF |
回复此楼» 猜你喜欢
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有8人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有4人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
COMSOL Multiphysics中文使用手册和化学反应工程模块教程
已经有233人回复- 关于用超胞方法计算声子的一些经验
已经有21人回复
- 如何通过声子计算求解热导率
已经有5人回复
- Castep模块CPU和内存使用率
已经有7人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
7楼2010-05-24 22:56:04
lono75
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6963.1
- 散金: 100
- 红花: 2
- 帖子: 2433
- 在线: 536.4小时
- 虫号: 932458
- 注册: 2009-12-23
- 专业: 凝聚态物性 II ：电子结构
2楼2010-05-23 20:40:55
lono75
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 6963.1
- 散金: 100
- 红花: 2
- 帖子: 2433
- 在线: 536.4小时
- 虫号: 932458
- 注册: 2009-12-23
- 专业: 凝聚态物性 II ：电子结构
3楼2010-05-23 20:43:24
4楼2010-05-24 10:17:21