24小时热门版块排行榜

>论坛更新日志 (4238)
>文献求助 (437)
>导师招生 (381)
>虫友互识 (381)
>休闲灌水 (259)
>考博 (172)
>硕博家园 (127)
>招聘信息布告栏 (120)
>博后之家 (98)
>论文投稿 (93)
>公派出国 (89)
>基金申请 (73)
>考研 (68)
>找工作 (60)
>论文道贺祈福 (58)
>教师之家 (41)

返回列表

maomao-666

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 4915.6
帖子: 74
在线: 122.7小时
虫号: 920616

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

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

回复此楼

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

siesta问题

» 猜你喜欢

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

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人回复

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

查看全部散金贴

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