| 查看: 2001 | 回复: 22 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
Abinit 中计算 PDOS 和PJPDOS
|
|||
|
前面有人问abinit中如何计算声子谱，通过这两天的学习，现在把计算声子谱的方法贴出来，供大家使用（以abinit自带的AlAs教程为例） 1. 原来的老办法 在l计算热力学性质的Log文件中提取DOS信息，然后作图 2. Abinit 5.6.3中的方法，我是根据今天在abinit forum里看到程序作者的建议来操作的。大家可以去forum mailing list 看看 重新编译了anaddb后，计算出了声子态密度和分态密度。 下面给出我两个方法计算出来的态密度图 PS：发现最新的abinit版本中很多功能是我们在手册上是没有的，需要我们多多交流，多多研究 |
回复此楼» 猜你喜欢
浙江师范大学是怎么坑我的
已经有10人回复
-
浙江师范大学是怎么坑我的
已经有3人回复
-
求教各位大神：2026教育部人文社科青年基金项目何时公示呀？
已经有6人回复
-
面上有专家说实验设备不是我们单位的
已经有8人回复
-
科研求助
已经有4人回复
-
博士去二本高校当辅导员，如何调整心态？
已经有11人回复
-
江苏省自然基金 什么时候出结果
已经有6人回复
-
闲聊
已经有3人回复
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
8楼2008-12-10 17:25:00
★ ★
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