| 查看: 3614 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- COMSOL Multiphysics中文使用手册和化学反应工程模块教程
已经有233人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- 如何通过声子计算求解热导率
已经有5人回复
- Castep模块CPU和内存使用率
已经有7人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
6楼2010-05-24 22:31:42
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
40