| 查看: 1245 | 回复: 1 | ||||
[交流]
【求助】siesta中PDOS文件的处理
|
| 请教高手，如题，在siesta中怎么处理PDOS文件！ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- siesta-3.1编译成功，但有几个问题求助
已经有7人回复
- 求助Siesta计算熟手，想学习此软件
已经有8人回复
- SIESTA中能带文件*.bands的处理
已经有49人回复
- origin 做DOS 图
已经有7人回复
- 【求助】H2 Binding Energy@SIESTA
已经有13人回复
- 【求助】请教在ATK中怎样使用自己的或者Siesta的赝势
已经有4人回复
- 【求助】siesta输入文件中赝势文件的作用？
已经有8人回复
- 【求助】siesta2.0.2 安装有些警告
已经有4人回复
- 【求助】如何用origin做pdos图!
已经有9人回复
- 【求助】siesta3.0 安装，错误一箩筐
已经有11人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】请教siesta3.0-beta如何计算PDOS、overlap population及coop
已经有12人回复
★ ★ ★ ★
cenwanglai(金币+2):谢谢回复！ 2010-11-10 17:33:18
sunyang1988(金币+2):呵呵，辛苦了 2010-11-10 19:00:47
maomao-666(金币+5): 2010-11-10 19:03:48
cenwanglai(金币+2):谢谢回复！ 2010-11-10 17:33:18
sunyang1988(金币+2):呵呵，辛苦了 2010-11-10 19:00:47
maomao-666(金币+5): 2010-11-10 19:03:48
|
百度来的哈~~ 一： pdosxm 我主要参考的http://fisica.ehu.es/ag/siesta-extra/issues.html 使得siesta输出的PDOS文件中没有＊。 另外对于nspin=2的情况，修改了pdosxml的代码，使得其能用。 （1） 在arch.make中加入 FPPFLAGS= -DWXML_INIT_FIX $make clean , $make 这样siesta输出的PDOS文件中就没有＊。 （2）推荐使用/Util/pdosxml，来读取PDOS文件，并转化成gnuplot能处理的形式 使用参见/Util/pdosxml/README 主要步骤如下: 1 在m_orbital_chooser.f90 中选好你要的轨道 2 make 3 ./pdos h2o_dos.PDOS >out 采用原来的代码，计算nspin=2时会出错，http://fisica.ehu.es/ag/siesta-extra/issues.html给出的 m_pdos.f90文件 其实和原来的一样。我们要自己改一下代码 在 m_pdos.f90 大概150多行的地方： if (spin_polarized) then dos1 = dos1 + data(1:ndata:2) dos2 = dos2 + data(2:ndata:2) else dos1 = dos1 + data endif ｜｜ 替换成 ｜｜ ｜｜ \ / if (spin_polarized) then !shanghui add here to make nspin=2's pdos do i=1,ndata/2 dos1(i) = dos1(i) + data(2*i-1) dos2(i) = dos2(i) - data(2*i) enddo !shanghui end add here else dos1 = dos1 + data endif 这样你就能做出正负分开的PDOS了。 上面修改好的siesta放在 cluster4 /public/soft/siesta/siesta_Util_pdosxml 二 ： readpdos 我写的一个小代码，放在 210.45.121.200 ~/share/Seminar_backup/siesta_tutorial/soft_siesta/readpdos 可以对有＊号的PDOS文件都能处理，但是每次体系，都要重新修改代码 |
2楼2010-11-10 17:30:45
