24小时热门版块排行榜

返回列表

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

秋雨子

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1578.8
散金: 13
帖子: 163
在线: 78.2小时
虫号: 1027042
注册: 2010-05-24
专业: 理论环境化学

[求助] homo和lumo批量提取的脚本

最近在弄一批分子，一万个左右，想提取其中的homo和lumo能量值，但是发现很难实现，因为才接触脚本半年，还是自学的，所以搞不定了，有哪位大侠给个小脚本？我只会一点shell，提取个Energy还行，homo和lumo不好提取，若可以分享一下不胜感激！

回复此楼

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

搜一搜，收一收

» 猜你喜欢

湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名已经有5人回复
东方理-中科大联合博士生招聘已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有156人回复
关于五水硝酸铋还原得到单质铋的问题已经有1人回复
阴离子交换膜电解二氧化碳还原已经有0人回复
山东省青年基金什么时候拨款已经有34人回复
CES投稿已经有2人回复
如何掌握一门外语？已经有0人回复
企业为什么一定要做“专精特新”认定？如何成为专精特新企业？已经有5人回复
求助ASIQwin新版本的软件已经有0人回复
期待科研合作，共同发表论文已经有197人回复

» 本主题相关商家推荐: (我也要在这里推广)

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

To the time to life, rather than to life in time

1楼 2011-12-05 22:47:40

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

asaka

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 392.2
红花: 2
帖子: 36
在线: 35.5小时
虫号: 1369381
注册: 2011-08-16
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★
zhou2009(金币+5): 2011-12-11 20:10:45
秋雨子(金币+20): 谢谢您的费心！只怪我学艺不精，只会shell！仍然非常感谢 2011-12-11 20:28:37

说一下我的想法啊：

首先用grep把有用的输出信息过滤一下：
一个是分割各个计算结果，免得把好几个计算的信息混在一起，
另外也是楼主要求的，只输出最后的结果，
另一个就是各个轨道能的数据。

然后把占据轨道的能级收集起来，取最大值，就是HOMO
以及把空轨道的能级收集起来，取最小值，就是LUMO
但是shell没有数学计算的功能，只好借用一下python命令

CODE:

#!/bin/sh

block='Population'

eigen='eigenvalues'

occst='occ'

virst='vir'

sepst='--'

for file in $*; do

  echo "File : $file"

  grep "\($block\|$eigen\)" $file > tmp

  while read line; do

    case $line in

      *$block*)

        eocc='' ;

        evir='' ;;

      *$occst*)

        data=`echo $line | sed "s/^.*$sepst\(.*\)$/\1/"` ;

        eocc=`echo $eocc $data` ;;

      *$virst*)

        data=`echo $line | sed "s/^.*$sepst\(.*\)$/\1/"` ;

        evir=`echo $evir $data` ;;

    esac

  done < tmp

  data=`echo $eocc | sed "s/ /,/g"`

  echo "  HOMO :" `python -c "print max($data)"`

  data=`echo $evir | sed "s/ /,/g"`

  echo "  LUMO :" `python -c "print min($data)"`

  rm tmp

done

以上程序需要bash和python就可以运行，测试LZ的文件，输出如下：
$ sh homolumo Row1.log
File : Row1.log
HOMO : -0.40657
LUMO : 0.13145
其中homolumo是上面shell的文件名。

当然楼主可能有其他类型的输出文件，
修改shell程序中部分的关键字和流程，应该可以应付其他情况。
block='Population' #不同计算间的分割行
eigen='eigenvalues' #包含轨道能级的输出行
occst='occ' #标记占据轨道的字符
virst='vir' #标记空轨道的字符
sepst='--' #分割信息和能级数据的字符