| 查看: 1737 | 回复: 22 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
Abinit 中计算 PDOS 和PJPDOS
|
|||
|
前面有人问abinit中如何计算声子谱，通过这两天的学习，现在把计算声子谱的方法贴出来，供大家使用（以abinit自带的AlAs教程为例） 1. 原来的老办法 在l计算热力学性质的Log文件中提取DOS信息，然后作图 2. Abinit 5.6.3中的方法，我是根据今天在abinit forum里看到程序作者的建议来操作的。大家可以去forum mailing list 看看 重新编译了anaddb后，计算出了声子态密度和分态密度。 下面给出我两个方法计算出来的态密度图 PS：发现最新的abinit版本中很多功能是我们在手册上是没有的，需要我们多多交流，多多研究 |
回复此楼» 猜你喜欢
酰胺脱乙酰基
已经有9人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有3人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
21楼2008-12-19 20:58:11
★ ★
fegg7502(金币+2,VIP+0):thanks
fegg7502(金币+2,VIP+0):thanks
|
mailing list中的邮件原文，大家根据邮件来改anaddb就可以了 Digging into my files I found an anaddb input that can be used to calculate the PHDOS. The flags that have to be added in the input file are: ifcflag 1 #IFC's are needed #DOS related variables prtdos 1 # 1 for gaussian, 2 for tetrahedrons ng2qpt 35 35 7 # Dense q-mesh for Fourier interpolation q2shft 0.0 0.0 0.0 # Shift for the dense mesh dosdeltae 9.112670E-7 # 0.1 cm-1 frequency step for PHDOS dosmear 1.366900E-5 # 3 cm-1 broadening (only for Gaussian method) At the end of the calculation the results are stored in the PHDOS file There are several columns: # omega PHDOS IPHDOS PJDOS[1] IPJDOS[1] .. PHDOS refers to the total phonon DOS, IPHDOS is the integrated DOS, then we have the contribution due to all the atoms of the first type as well as the corresponding integrated dos. 3 Comments: 1) To call mkphdos you have to modify a bit 17ddb/invars9 commenting out the call to leave_new in the following piece of code if(anaddb_dtset%prtdos/=0 .and. anaddb_dtset%ifcflag/=1) then write(message, '(5a)' )& & ' invars9 : ERROR -',ch10,& & ' ifcflag must be 1 when the calculation of the phonon DOS is required',ch10,& & ' Action : correct ifcflag in your input file.' call wrtout(6,message,'COLL') call leave_new('COLL') end if 2) If I remember well in old versions of abinit there's a small bug (just a missing factor) in the calculation of the integrated type-projected PHDOS. So don't rely on the columns reporting the type-projected IPHDOS. If you need, you can easily integrate the curve using a plotting software such as xmgrace. 3) Be careful when using the tetrahedron method. The CPU time doesn't scale well wrt to the number of Fourier interpolated q-points. This is a general problem I found also in case of electronic DOS calculations. One should do some serious profiling of the tetrahedron routines and the fortran loops to individuate the bottleneck |
2楼2008-12-04 14:13:07
3楼2008-12-04 14:17:20
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
5楼2008-12-04 19:37:52