| 查看: 1145 | 回复: 1 | ||||
[交流]
【求助】siesta中PDOS文件的处理
|
| 请教高手，如题，在siesta中怎么处理PDOS文件！ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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/91
-
江苏师范大学化学及材料专业招收调剂研究生
+1/91
-
南京林业大学的土木工程是咋回事？
+1/91
-
085404 计算机技术 真诚求老师带走 杭电 336 六级已过
+1/88
-
坐标北京，诚征女友
+1/76
-
浙江理工大学蒋仲庆 物理/新能源/材料 招9名
+2/56
-
青岛农业大学海洋学院招调剂硕士研究生
+1/51
-
需要调剂化学专业的看过来
+1/40
-
复旦大学院士团队招聘药学博士后/科研助理
+1/34
-
北京理工大学（珠海）招收2026级考核制博士
+1/21
-
吉大生科博一博士联培自荐-已发表+在投一作SCI3篇-踏实肯干能加班
+1/21
-
华南师范大学环境学院大分子体系多尺度模拟与实验研究团队招收2026年考研调剂生
+1/20
-
河北大学分析化学招收多名调剂考生
+1/14
-
福建师范大学 2026硕士调剂 （化学、光电器件、物理、纳米材料方向）
+1/8
-
26年申博自荐-计算机视觉
+1/7
-
山东交通学院航运学院关于2026年硕士研究生招生考试调剂复试有关事项的公告
+1/5
-
北大口腔招人工智能、生物再造科研助理（优秀者可留校或推荐读博）
+1/5
-
【#上海调剂急录#】能接受985联合培养的速来！带你发一区文章！
+1/4
-
长江师范学院化学化工学院化学工程方向接收大量调剂考生
+1/4
-
【#上海调剂急录#】能接受985联合培养的速来！带你发一区文章！
+1/2
★ ★ ★ ★
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
