24小时热门版块排行榜

>论坛更新日志 (4303)
>文献求助 (507)
>导师招生 (435)
>虫友互识 (394)
>考博 (185)
>基金申请 (159)
>招聘信息布告栏 (150)
>硕博家园 (146)
>论文投稿 (132)
>考研 (117)
>休闲灌水 (111)
>博后之家 (78)
>教师之家 (59)
>找工作 (39)
>物理 (35)
>公派出国 (34)

返回列表

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

秋雨子

金虫 (小有名气)

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

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

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

回复此楼

» 猜你喜欢

中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事已经有0人回复
九江学院2026年最新高层次人才招聘公告已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有134人回复
深圳信息职业技术学院2026年教师招聘公告（最新）已经有0人回复
研究生招生已经有0人回复
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生已经有3人回复
一志愿南师大0703化学 275求调剂已经有6人回复
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向已经有0人回复
欢迎调剂长江师范学院汤老师课题组的硕士研究生已经有0人回复
中国电机工程学报的投稿系统已经有0人回复

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

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