24小时热门版块排行榜

返回列表

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

秋雨子

金虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

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

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

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

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

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

mchen10

金虫 (正式写手)

应助: 13 (小学生)
金币: 1508.2
散金: 41
红花: 6
帖子: 425
在线: 159.5小时
虫号: 695973
注册: 2009-02-04
性别: GG
专业: 理论和计算化学

【答案】应助回帖

引用回帖:

18楼: Originally posted by asaka at 2011-12-09 20:33:17:
以上作者只考虑闭壳层计算的情况，你有开壳层的计算吗？那时候你需要总的HOMO,LUMO还是分别alpha以及beta的HOMO和LUMO?
他们的程序都只有一点小问题，相信LZ可以搞定。
mchen10 的程序少了一个 |，多了一个奇怪 ...

CODE:

#!/bin/bash

# print homo lumo energy of gaussian output

# usage: print_homo_lumo filename(s)

for file in "$@"

do

echo $file

tac $file | tr "\n" "@" | egrep -o "@ Alpha virt\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues.{0,75}@ Alpha  occ\. eigenvalues" | tr "@" "\n" | head -n3 | awk '{if ($2 == "virt.") print $1" LUMO "$5; if($2 == "occ."){i=NF; print $1" HOMO "$i} }'

done

抱歉，tio2.rev是测试产物去掉就好了（留着也不影响）
至于开壳闭壳的问题，我们给个思路楼主如果什么都让别人来搞定就算程序可以用真正的收益也就那么一点点