24小时热门版块排行榜

返回列表

【奖励】本帖被评价37次，作者neweroica增加金币 28.3 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

neweroica

木虫 (著名写手)

应助: 53 (初中生)
金币: 5007.5
帖子: 2473
在线: 216.6小时
虫号: 36775

[资源] 【neweroica个人文集】操作VASP的CHGCAR文件的程序：chgcarman（v0.0）

在前面的帖子中，我发布将CHGCAR转换为cube文件的程序：chgcar2cube。见：http://muchong.com/bbs/viewthread.php?tid=1809216

有时候，我们需要对多个CHGCAR文件进行加减乘除等操作（比如计算电荷密度的变化）。Gaussian提供了cubeman程序用于操作cube文件，但我发现有时会出现bug。所以，我干脆编写了这个chgcarman程序，用于直接对CHGCAR文件进行进行操作。

目前chgcarman的主要功能有：

1. 输出CHGCAR文件的属性，包括：最大值，最小值，平均值，RMS值，正负值的RMS值。这样在绘制等值面图时，有助于选取合适的等值面的数值。
2. 将两个CHGCAR文件相加。
3. 将两个CHGCAR文件相减。
4. 将CHGCAR文件乘以一个缩放因子。
5. 将CHGCAR文件平方。
6. 将CHGCAR文件取常用对数（以10为底）。

建议在Linux下进行编译和使用，则可用标准C编译器进行编译，如：
gcc -o chgcarman -lm chgcarman.c

如果你非要用Windows的话，请先设置源文件中的MS_Window值为1，然后再用C编译器进行编译。这很重要！！！否则程序虽然编译通过，但转换结果会发生严重错误。

注意：目前该程序只在Linux下测试通过，Windows下我还没有测试过，也欢迎大家测试后向我反馈结果。

程序使用帮助：
      chgcarman [option] [ file ...]
      Options:
         -a
               Add two CHGCAR files to produce a new one.
         -h
               Show information and usage of this program and exit.
         -l
               Compute the logarithm of a CHGCAR file to produce a new one.
         -p
               Get properties of a CHGCAR file  without producing a new one.
         -q
               Compute the square of a CHGCAR file to produce a new one.
         -s
               Substract two CHGCAR files to produce a new one.
         -x factor
               Scale a CHGCAR file by a scaling factor to produce a new one.
      Input files:
         VASP's CHGCAR or LOCPOT files

[ Last edited by yjcmwgk on 2010-6-16 at 21:41 ]

回复此楼

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

材料计算模拟实用技巧	第一性原理相关文档	vasp计算资料	Photochemistry
量化经验	淘淘

» 猜你喜欢

求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已经有2人回复
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有123人回复
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料已经有0人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
英国利物浦大学屠昕和雍雪课题组招收CSC博士生已经有26人回复
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）已经有0人回复
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建已经有0人回复
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人已经有45人回复
东方理-中科大联合博士生招聘已经有0人回复

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

CHGCAR的文件含义已经有7人回复
请教一个运行vasp的问题已经有7人回复
求gaussian98 windows的源文件已经有5人回复
[font=宋体]修改.cpp文件怎么操作[/font] 已经有3人回复
量子化学课题组vasp经典教程已经有64人回复
关于VASP计算得到CHGCAR文件的数据处理以及数据单位的已经有16人回复
vasp算出来的文件基本都是0字节的已经有24人回复
vasp中CHGCAR是PS电荷还是AE电荷已经有6人回复
对CHGCAR进行加减已经有3人回复
为什么没有产生 CHG 文件？（MedeA-VASP5.2）已经有5人回复
【求助】CHGCAR用什么来读取已经有10人回复
【求助】帮忙VASP写个小程序已经有21人回复
【求助】帮忙编个用于VASP小程序已经有22人回复
【neweroica个人文集】转换VASP的CHGCAR为cube文件的程序chgcar2cube（v0.0）已经有45人回复

1楼2010-01-24 10:45:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

stefsui

木虫 (小有名气)

应助: 11 (小学生)
金币: 2252.3
帖子: 155
在线: 472.7小时
虫号: 872985

引用回帖:

5楼: Originally posted by bingmou at 2010-01-26 15:41:50:
谢了，用linux的shell在这一点还是不行，awk只能进行单个文件的处理（或者我不懂），所以这个软件太棒了！一定要学习一下，以后自己可以写。希望看不懂问lz，lz可以回复。谢谢

awk可以处理，单纯相减的话
awk很强大.....