应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【已完结】如何计算macroscopic average
161/1返回列表
查看: 4132  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

huangll99

木虫 (职业作家)

[交流] 【已完结】如何计算macroscopic average

Siesta中有一个很好的tool-->macroave,可以用来计算一些与表面和界面相关的物理量。主要就是计算planar average 和 macroscopic average。
由于我对Siesta完全不懂，想要请教一下各位，Siesta中保存电荷密度和静电势的文件的格式是怎样的，如何把VASP中的CHGCAR 和 LOCPOT 的格式转换成Siesta中的格式，然后借助macroave这个工具来做分析

弄好了

[ Last edited by huangll99 on 2012-5-14 at 16:56 ]
回复此楼

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

@个人收集@ first principle pwscf Siesta

» 猜你喜欢

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

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

查看全部散金贴
1楼 2012-05-13 09:19:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

huangll99

木虫 (职业作家)
上传一个自己找到的资料
赞一下(1人) 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : macroave.f
    • 2012-05-13 09:21:19, 21.24 K
  • 附件 2 : macroave.pdf
    • 2012-05-13 09:21:19, 68.42 K
2楼2012-05-13 09:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

huangll99

木虫 (职业作家)
没人理，自己顶一顶

把LOCPOT 转换成PWscf中可以读的格式，然后借助average.f90这个小工具来处理有没有可行性呢，

请goldenfisher专家来帮帮忙  @goldenfisher

[ Last edited by huangll99 on 2012-5-14 at 08:27 ]
赞一下 回复此楼
3楼2012-05-14 08:15:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+1, 专家考核 2012-05-14 12:27:05
引用回帖:
3楼: Originally posted by huangll99 at 2012-05-14 08:15:20:
没人理，自己顶一顶

把LOCPOT 转换成PWscf中可以读的格式，然后借助average.f90这个小工具来处理有没有可行性呢，

请goldenfisher专家来帮帮忙  goldenfisher

siesta我不太了解，没用过。 pwscf里面average.f90貌似是可以计算你自己说的系综平均的，你可以试试
赞一下 回复此楼
4楼2012-05-14 11:56:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
引用回帖:
4楼: Originally posted by goldenfisher at 2012-05-14 11:56:28:
siesta我不太了解，没用过。 pwscf里面average.f90貌似是可以计算你自己说的系综平均的，你可以试试

谢谢指教。我读了average.f90的代码，是可以做的。主要是我不会用PWscf,看了一下PWscf的论坛，应该是用average.x来处理pp.x产生的文件。现在主要是不知道PWscf里面保存电荷密度和静电势的文件是什么，格式是什么样的。我再学习一下PWscf吧。
PS:我用average.f90主要是想处理静电势文件，计算异质节的能带带阶
赞一下 回复此楼
5楼2012-05-14 12:26:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
哇哈哈，PWscf中的average.f90果然可以，就是不知道PWscf中的静电势文件里的这些内容是什么意思
                                                                           
      15      15     144      15      15     144       4       1
     6  5.42350901  0.00000000  6.36396103  0.00000000  0.00000000  0.00000000
       44.7046033528        4.0000000000       15.0000000000    11
   1   Al    3.00
   1       0.000000000    0.000000000    4.242640687    1
   2       0.500000000    0.500000000    3.535533906    1
   3       0.000000000    0.000000000    2.828427125    1
   4       0.500000000    0.500000000    2.121320344    1

传一张效果图


俺一个人自言自语，版主可不要给偶删帖啊



[ Last edited by huangll99 on 2012-5-14 at 13:09 ]
赞一下 回复此楼
6楼2012-05-14 13:03:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
huangll99: 金币+10, 非常感谢^-^ 2012-05-14 14:41:24
franch: 专家考核 2012-05-14 15:48:16
我不太明白你说的静电势在这里为何意？
存储电荷密度的文件就是pp.x计算后产生的.charge文件，就是pp.x计算的输入文件中命名的那个文件。
上面的意思就是什么格点数目啊，晶格常数啊，等等，这些基本的意思。
赞一下 回复此楼
7楼2012-05-14 14:25:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
引用回帖:
7楼: Originally posted by goldenfisher at 2012-05-14 14:25:22:
我不太明白你说的静电势在这里为何意？
存储电荷密度的文件就是pp.x计算后产生的.charge文件，就是pp.x计算的输入文件中命名的那个文件。
上面的意思就是什么格点数目啊，晶格常数啊，等等，这些基本的意思。

average.f90这个程序既可以处理电荷密度的文件，也可以处理potential (我说的静电势就是把total potential 里的交换关联势去掉后的势),这两个文件的格式类似。
这几行的意思俺还没有搞懂，还需要学习，相信肯定能把VASP的文件转成这种格式
多谢 goldenfisher 专家指点了
赞一下 回复此楼
8楼2012-05-14 14:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
huangll99: 金币+10, 非常感谢^-^ 2012-05-15 17:11:59
看来LZ不喜欢多文件对比哦～其实查文件内容的规律是很有意思的一件事的。
检查Al.scf.in，可以看到，ibrav = 6，并且
CODE:
  celldm(1)   = 5.4235090117D0,
  celldm(3)   = 6.3639610306789276D0,

事实上，PWscf的celldm(i)有6个，不需要标明的是根据ibrav计算的，所以把ibrav和celldm(i)全列出来，就是：
CODE:
6  5.42350901  0.00000000  6.36396103  0.00000000  0.00000000  0.00000000

至于第一列，注意Al.scf.out中：
CODE:
FFT grid: ( 15, 15, 144)

以及Al.scf.in中：
CODE:
  ecutwfc     = 15.D0,
  nr3         = 144,

所以推测前面六个数字15      15     144      15      15     144分别是：nr1  nr2  nr3  nr1s  nr2s  nr3s（对FFT不熟悉，所以不保证推测正确）
第一列最后的4  1应该是四个原子，一种元素的意思。
赞一下 回复此楼
9楼2012-05-15 17:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)
引用回帖:
9楼: Originally posted by souledge at 2012-05-15 17:01:35:
看来LZ不喜欢多文件对比哦～其实查文件内容的规律是很有意思的一件事的。
检查Al.scf.in，可以看到，ibrav = 6，并且  celldm(1)   = 5.4235090117D0,
  celldm(3)   = 6.3639610306789276D0,
事实上，PWscf的 ...



说的好清楚！谢谢souledge专家
赞一下 回复此楼
10楼2012-05-15 17:13:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoboy185

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，请问您的问题搞定了吗  能否分享一下
赞一下 回复此楼
11楼2012-10-30 12:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leonidasn

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
同问楼主能否分享下呢？
赞一下 回复此楼
12楼2012-11-07 16:52:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

光尘如水

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
偶也想做vasp的macroscopic average ，楼主指点一下呀
赞一下 回复此楼
13楼2013-07-01 16:54:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyuewen

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
我不知道楼主最终是以何种方式解决这个问题的，但是我个人觉得不能用average.x来处理。因为用这个的时候得需要用到pwscf里面的上一步pp.x处理得到的一些文件。没错。PWSCF中的CHARGE和VASP得出的LOCPOT简直是一样的，除了以上侯老师指出的部分，但是pwscf里面需要一个输入的in文件，这样可能就不易实现了。
赞一下 回复此楼
14楼2013-10-16 22:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vinbert

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
LZ，能传我一个计算macroscopic average potential的程序么？ 邮箱：vinberg.lee@mail.dlut.edu.cn
谢谢先～
赞一下 回复此楼
15楼2014-04-24 18:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luomin_zyj

捐助贵宾 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主，macroscopic average怎么处理的，我也是VASP算LOCPO

怎么将LOCPOT转成Siesta的格式 ?x谢谢楼主
赞一下 回复此楼
16楼2016-05-30 14:48:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huangll99 的主题更新
161/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 18:46 by lXgU0RvNTDC9
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +3 process2012 2026-08-23 4/200 2026-08-23 12:54 by huixian257
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭