[交流] 【求助】siesta中PDOS文件的处理

请教高手，如题，在siesta中怎么处理PDOS文件！

回复此楼

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-10 16:41:37

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

xiaochun37

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 206.7
帖子: 19
在线: 1.6小时
虫号: 1121966

★ ★ ★ ★
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人)

回复此楼

2楼2010-11-10 17:30:45

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

[文学芳草园] 厉娜 +3	张丹枫1991 2026-11-09	5/250	2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5	kvewwGtKktqB 2026-11-04	12/600	2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5	ttH3f1oA6ace 2026-11-04	15/750	2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5	5sJFDKp2kwuI 2026-11-06	16/800	2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5	5sJFDKp2kwuI 2026-11-06	15/750	2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4	5sJFDKp2kwuI 2026-11-06	10/500	2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3	5sJFDKp2kwuI 2026-11-06	12/600	2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4	5sJFDKp2kwuI 2026-11-06	8/400	2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	5sJFDKp2kwuI 2026-11-06	8/400	2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5	5sJFDKp2kwuI 2026-11-06	9/450	2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3	yexuqing 2026-11-04	5/250	2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5	ttH3f1oA6ace 2026-11-04	10/500	2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6	ttH3f1oA6ace 2026-11-04	12/600	2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4	ttH3f1oA6ace 2026-11-04	11/550	2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5	kvewwGtKktqB 2026-11-04	14/700	2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6	kvewwGtKktqB 2026-11-04	16/800	2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4	ttH3f1oA6ace 2026-11-04	7/350	2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4	ttH3f1oA6ace 2026-11-04	10/500	2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	ttH3f1oA6ace 2026-11-04	9/450	2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4	hNpp4Y2CGvSC 2026-11-04	12/600	2026-11-06 01:53 by qL1cU46YIdK0

24小时热门版块排行榜

maomao-666

[交流] 【求助】siesta中PDOS文件的处理

» 收录本帖的淘帖专辑推荐

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到: 查看全部散金贴

xiaochun37

» 抢金币啦！回帖就可以得到:

查看全部散金贴