24小时热门版块排行榜    

查看: 1493  |  回复: 4

伐木工

木虫 (初入文坛)


[交流] 【求助】怎么将siesta中的静电势ElectrostaticPotential的数据导出啊

如题,我刚刚学习siesta,想计算一下静电势,按照manu上将SaveElectrostaticPotential 设为True,可是计算成功后如何将有效势画成曲线呢,从哪里看静电势的数据啊
谢谢高手帮忙啦
回复此楼

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

Siesta

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
sunyang1988(金币+2): 谢谢指点 2011-04-02 22:44:38
伐木工(金币+10): 2011-04-03 09:38:59
引用回帖:
Originally posted by 伐木工 at 2011-04-02 13:18:53:
如题,我刚刚学习siesta,想计算一下静电势,按照manu上将SaveElectrostaticPotential 设为True,可是计算成功后如何将有效势画成曲线呢,从哪里看静电势的数据啊
谢谢高手帮忙啦

参照下面的两个帖子,估计会解决你的问题。
http://muchong.com/bbs/viewthread.php?tid=2482325
http://muchong.com/bbs/viewthread.php?tid=2386188
在那里,把电荷密度的文件替换为你现在的电势的文件.
2楼2011-04-02 21:56:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伐木工

木虫 (初入文坛)


引用回帖:
Originally posted by zhangguangping at 2011-04-02 21:56:33:
参照下面的两个帖子,估计会解决你的问题。
http://muchong.com/bbs/viewthread.php?tid=2482325
http://muchong.com/bbs/viewthread.php?tid=2386188
在那里,把电荷密度的文件替换为你现 ...

你好 谢谢你的回帖 这两个网址我仔细看过了 可是这里没有计算有效势的程序 我在siesta manu上看到说要用程序IORHO
Instructs to write the total electrostatic potential,
de ned as the sum of the hartree potential plus the local pseudopotential, at the mesh
used by DHSCF, in file SystemLabel.VH. This file can be read by routine IORHO, which may be used by an application program in later versions.
Use: File SystemLabel.VH is only written, not read, by siesta.
这是说现在的版本还不能处理有效势吗,可是在siesta原来的文献上不是计算了一个两个金电极 中间是真空层的有效势了吗 他们是怎么算的呢?有些迷茫了 请多多指教
3楼2011-04-29 19:31:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
gzqdyouxia(金币+2): 谢谢详细指导 2011-04-29 22:42:31
引用回帖:
Originally posted by 伐木工 at 2011-04-29 12:31:39:
你好 谢谢你的回帖 这两个网址我仔细看过了 可是这里没有计算有效势的程序 我在siesta manu上看到说要用程序IORHO
Instructs to write the total electrostatic potential,
de ned as the sum of the hartre ...

你不是已经有了.VH的文件了吗?如果你用第一个帖子的方法,可以将.VH文件中的数据转化成cube文件的格式,画等值面,等高线。
如果你用第二个帖子中我在7楼中附件中的程序,可以画出来某一条直线或者某一个平面,或者是某一个点上的电势的值。
当然,如果你只想沿着z方向画电势,但是对xy平面取平均,那么你需要利用siesta带的工具Macroave
我想你说的那个PRB,65,165401(2002)的文章中的图是用我说的第二种方法画的。
4楼2011-04-29 20:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伐木工

木虫 (初入文坛)


谢谢

[ Last edited by 伐木工 on 2011-5-2 at 09:07 ]
5楼2011-04-30 14:58:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 伐木工 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 337一志愿华南理工材料求调剂(有希望2吗?) +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 数二英二309分请求调剂 +3 dtdxzxx 2026-03-09 4/200 2026-03-09 19:56 by yuningshan
信息提示
请填处理意见