autodataming的个人博客分享 http://blog.sciencenet.cn/u/autodataming

博文

如何批量提取分子的HOMO lumo的值

已有 7780 次阅读 2014-9-12 10:04 |个人分类:DrugDesign|系统分类:科研笔记

如果需要显示HOMO LUMO 轨道, 电脑中需要同时安装Gaussian 和GView.







背景:

课题中被要求提取HOMO lumo的值

有人说从fchk(chk)中提取,又有人说可以从log中提取。

他们的做法是用gauss view打开fchk文件,或者log文件。

然后右击--》edit--》MOs,没必要选中任何原子,因为homo lumo是针对分子的而不是原子的。

我在chk文件中没有找到HOMO lumo的值,在log文件中找到了HOMO LUMO的值

方法

从log文件中自动提取HOMO lumo的值,主要是用正则表达式。


extractHOMOLUMO.pl


最新版

extractHOMOLUMO.pl





#!/usr/bin/perl -w

#usage:perl extractHOMOLUMO.pl   xxx.log
#author: dddc chen zhaoqiang
#blog: http://blog.sciencenet.cn/home.php?mod=space&uid=950202




use strict;

$/=undef;
open FH,$ARGV[0];
#open FH,"1.log";
my $text=<FH>;

# Alpha  occ. eigenvalues --   -0.29762  -0.26789
#Alpha virt. eigenvalues --   -0.02175   0.00021   0.00256   0.00998   0.01628
 
if($text=~/.*Alpha  occ. eigenvalues .+?\s(\S+)\n\sAlpha virt. eigenvalues --   (\S+)/ms)
{
   
   print "homo ",$1," lumo ",$2;
}


###########批量提取log文件中的HOMO LUMO##################################

我建议把刚才的脚本放到环境变量PATH所在的目录中;

这样就可以把extractHOMOLUMO.pl当成linux中的一个命令使用,

for i in $(ls *log);do echo -n "$i ";extractHOMOLUMO.pl $i;echo "";done

##########################################################################



https://wap.sciencenet.cn/blog-950202-827091.html

上一篇:multiwfn 进行AIM分析
下一篇:提取高斯log文件中的晗 熵 吉布斯自由能
收藏 IP: 202.127.19.*| 热度|

1 刘东坡

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 09:11

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部