| 查看: 1128 | 回复: 1 | ||||
[交流]
【求助】siesta中PDOS文件的处理
|
| 请教高手，如题，在siesta中怎么处理PDOS文件！ |
回复此楼» 收录本帖的淘帖专辑推荐
 siesta问题 |
» 猜你喜欢
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有4人回复
-
今年春晚有几个节目很不错，点赞！
已经有9人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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/167
-
坐标广州，征女友
+2/130
-
天津科技大学海洋与环境学院殷焕顺团队招博士生1名---分析化学领域
+1/81
-
ChineseResearchLaTeX: 开源、免费的vibe coding辅助国自然写作
+1/80
-
一个陌生女人的来信
+1/59
-
上海理工大学2026年系统科学学科海外骨干教师招聘启事
+2/34
-
海法大学线上开放日
+1/33
-
清华大学深圳国际研究生院招聘-博士后（长期有效）
+1/27
-
海南大学海洋技术与装备学院-科研助理招聘（可读博）膜分离水处理方向
+1/27
-
英国布里斯托大学诚招博士生，博士后和联合培养生
+1/17
-
2026年天津科技大学“新能源催化与膜材料团队”研究生招生
+1/17
-
澳洲皇家墨尔本理工RMIT招收大模型安全/隐私PhD
+1/9
-
26年材料/化学相关专业-博士研究生招生-温州大学化学与材料工程学院(26年秋季入学)
+1/9
-
上海交通大学-宁波东方理工大学联合培养博士生
+1/8
-
湖南大学-分析检测技术和生物柔性传感器-招收1名博士研究生 (2026年，第二批)
+1/6
-
英国南安普顿大学禅铎课题组诚招气候动力方向博士后
+1/6
-
全奖博士 英国利物浦大学+台湾清华大学 双博士学位
+1/5
-
2月22日 以色列海法大学 线上开放日
+1/3
-
怎么发布了求助贴了， 一发就转到删除栏了
+1/3
-
南京大学能源与资源学院徐加陵课题组招聘：科研助理、硕士生、博士生
+1/1
★ ★ ★ ★
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