24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

1ST强帖: 16
应助: 46 (小学生)
贵宾: 2.516
金币: 253.4
散金: 2603
红花: 44
帖子: 1283
在线: 500.5小时
虫号: 842155
注册: 2009-09-06
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】关于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 < &input
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 < 5 500 5
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 < PHDOS.out
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 < &input
amass(1)=26.98,
amass(2)=74.922,
asr='crystal',
flfrc=$FC_file,
flfrq='frequency'
&end
EOF

cat > Temperature < 5 500 5
EOF

cat >T_Debye.in < PHDOS.out
0.0001
3 15 3
500 10
EOF

回复此楼

» 猜你喜欢

上海工程技术大学张培磊教授团队招收博士生已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有7人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

COMSOL Multiphysics中文使用手册和化学反应工程模块教程已经有233人回复
关于用超胞方法计算声子的一些经验已经有21人回复
如何通过声子计算求解热导率已经有5人回复
Castep模块CPU和内存使用率已经有7人回复
【求助】siesta中transiesta模块的收敛标准是啥~~ 已经有20人回复
【求助】siesta中transiesta模块的收敛问题已经有13人回复

1楼 2010-05-23 15:49:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 6963.1
散金: 100
红花: 2
帖子: 2433
在线: 536.4小时
虫号: 932458
注册: 2009-12-23
专业: 凝聚态物性 II ：电子结构

identation(金币+5): 2010-05-24 22:56:14
identation(金币+18): 2010-05-25 15:35:17

# Parameters required for QHA calculations
# Total Phonon DOS file
# output file for C_V, S, Internal energy
# T_start, T_end, T_step for QHA calculations
#
cat >fqha.in < PHDOS.out
$Sysname.QHA.out
5 500 5
EOF

赞一下

回复此楼

由天命,成人事！

5楼2010-05-24 13:20:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 18 个回答

lono75

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 6963.1
散金: 100
红花: 2
帖子: 2433
在线: 536.4小时
虫号: 932458
注册: 2009-12-23
专业: 凝聚态物性 II ：电子结构

EOF (end of file) 作为标识

cat >T_Debye.in < PHDOS.out
0.0001
3 15 3
500 10
EOF

是将EOF前面的内容写入 T_Debye.in

qha.x < in.qha
qha.x 会读取in.qha 设定的参数和
matdyn.init Temperature T_Debye.in三个文件

[ Last edited by lono75 on 2010-5-23 at 20:59 ]

赞一下

回复此楼

由天命,成人事！

2楼2010-05-23 20:40:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖