| 查看: 2943 | 回复: 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 |
回复此楼» 猜你喜欢
AI论文写作工具：是科研加速器还是学术作弊器？
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
COMSOL Multiphysics中文使用手册和化学反应工程模块教程
已经有233人回复- 关于用超胞方法计算声子的一些经验
已经有21人回复
- 如何通过声子计算求解热导率
已经有5人回复
- Castep模块CPU和内存使用率
已经有7人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
又有问题了
|
图2示出了XXXX直流电导率与温度倒数（180─300K）的关系点。XXXX在180K的直流电导率在10-4S/m左右，在该温度XXXX呈现电阻特性。XXXX的电导率随温度上升而提高，并在室温时达到10-2S/m。计算得到的激活能为Ea=0.237eV（±0.001 eV）（用Arrhenius关系σ=σ0exp(-Ea/kBT)拟合得到，σ0是指前因子、Ea是激活能、kB是Boltzmann常数、T是温度），这是由自由电荷载流子（如电子、特定种类的缺陷、自由离子）的长程迁移所贡献的。如此低的激活能（<0.3 eV）通常表明，一般与钙钛矿型化合物高温烧结有关的氧空位不参与长程电导。带电载流子，特别是局域于Mn4+离子的两个自由电子，其长程跳跃可能对STMO的电导率产生显著的贡献。 请教如下问题 1、这里的激活能Ea=0.237eV（±0.001 eV）是带隙么？ 2、PWSCF的GGA+U计算只指定ldau=.true.和Hubbard_U(i)是不是就足够了？至于starting_ns_eigenvalue(m,ispin,I)需要指定么？ [ Last edited by identation on 2010-6-1 at 08:45 ] |
10楼2010-06-01 08:42:20
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