秋雨子

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時
蟲號: 1027042
注冊: 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

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

最近在弄一批分子，一萬個左右，想提取其中的homo和lumo能量值，但是發(fā)現(xiàn)很難實現(xiàn)，因為才接觸腳本半年，還是自學(xué)的，所以搞不定了，有哪位大俠給個小腳本？我只會一點shell，提取個Energy還行，homo和lumo不好提取，若可以分享一下不勝感激！

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

搜一搜，收一收

» 猜你喜歡

求標(biāo)準(zhǔn)卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340 已經(jīng)有2人回復(fù)
比利時魯汶大學(xué)與國家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費? 已經(jīng)有132人回復(fù)
量子點電致發(fā)光器件（QLED）：藍光量子點（QDs）和納米氧化鋅電子傳輸材料已經(jīng)有0人回復(fù)
需要一個圓偏振光誘導(dǎo)化學(xué)反應(yīng)，請問使用圓偏振光作為光源，應(yīng)該需要哪些設(shè)備搭建已經(jīng)有0人回復(fù)
半透性容器進行高濕或者低濕研究已經(jīng)有0人回復(fù)
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復(fù)
中國科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)
求ENDNOTE破解軟件安裝包已經(jīng)有4人回復(fù)
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實驗室丁勇團隊歡迎有志青年！已經(jīng)有3人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

重金求助HOMO/LUMO問題已經(jīng)有5人回復(fù)
【解決】討論，siesta中作HOMO 與LUMO圖的方法已經(jīng)有12人回復(fù)
【求助】Gaussian計算HOMO-LUMO問題已經(jīng)有17人回復(fù)
【求助】關(guān)于HOMO和LUMO軌道的確定問題已經(jīng)有8人回復(fù)
black dye 的HOMO與LUMO 已經(jīng)有6人回復(fù)
【求助】計算homo和lumo軌道出錯已經(jīng)有14人回復(fù)
【求助】怎么看HOMO LUMO 已經(jīng)有11人回復(fù)
【求助】關(guān)于HOMO和LUMO的問題已經(jīng)有5人回復(fù)
【求助】請問，如何在MS中分析HOMO與LUMO，謝謝您的幫助已經(jīng)有12人回復(fù)

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

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

asaka

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 392.2
紅花: 2
帖子: 36
在線: 35.5小時
蟲號: 1369381
注冊: 2011-08-16
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
zhou2009(金幣+5): 2011-12-11 20:10:45
秋雨子(金幣+20): 謝謝您的費心！只怪我學(xué)藝不精，只會shell！仍然非常感謝 2011-12-11 20:28:37

說一下我的想法�。�

首先用grep把有用的輸出信息過濾一下：
一個是分割各個計算結(jié)果，免得把好幾個計算的信息混在一起，
另外也是樓主要求的，只輸出最后的結(jié)果，
另一個就是各個軌道能的數(shù)據(jù)。

然后把占據(jù)軌道的能級收集起來，取最大值，就是HOMO
以及把空軌道的能級收集起來，取最小值，就是LUMO
但是shell沒有數(shù)學(xué)計算的功能，只好借用一下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的文件名。

當(dāng)然樓主可能有其他類型的輸出文件，
修改shell程序中部分的關(guān)鍵字和流程，應(yīng)該可以應(yīng)付其他情況。
block='Population' #不同計算間的分割行
eigen='eigenvalues' #包含軌道能級的輸出行
occst='occ' #標(biāo)記占據(jù)軌道的字符
virst='vir' #標(biāo)記空軌道的字符
sepst='--' #分割信息和能級數(shù)據(jù)的字符

贊一下(2人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

秋雨子

9樓2011-12-09 04:13:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

asaka

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 392.2
紅花: 2
帖子: 36
在線: 35.5小時
蟲號: 1369381
注冊: 2011-08-16
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
zhou2009(金幣+5): 2011-12-11 20:11:12

上面的帖子不能修改。稍微修改了一下程序，
速度快一點，另外避免了數(shù)據(jù)多于255個時python出錯的情況。

CODE:

#!/bin/sh

block='Population'

eigen='eigenvalues'

occst='occ'

virst='vir'

sepst='--'

for file in $*; do

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

  while read line; do

    case $line in

      *$block*)

        eocc='' ;

        evir='' ;;

      *$occst*)

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

        eocc="$eocc $data" ;;

      *$virst*)

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

        evir="$evir $data" ;;

    esac

  done < tmp

  rm tmp

  echo "File : $file"

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

  [ ! -z "$data" ] && echo "  HOMO :" `python -c "print max(($data))"`

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

  [ ! -z "$data" ] && echo "  LUMO :" `python -c "print min(($data))"`

done

贊一下(2人)

回復(fù)此樓

10樓2011-12-09 06:02:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

秋雨子

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時
蟲號: 1027042
注冊: 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

沒有人回復(fù)呢？！

回復(fù)此樓

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

2樓2011-12-06 14:11:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋雨子

金蟲 (小有名氣)

3樓2011-12-07 11:23:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

abbott

金蟲 (著名寫手)

不要用QQ問我東西

應(yīng)助: 16 (小學(xué)生)
金幣: 1071.6
散金: 2787
紅花: 10
帖子: 1015
在線: 105.2小時
蟲號: 452267
注冊: 2007-11-05
性別: GG
專業(yè): 計算機硬件技術(shù)

★ ★
小紅豆(金幣+2): 感謝交流。 2011-12-14 18:07:53

你先貼一個樣板文件,,,,
否則誰知道你要從什么文檔中提取你要的東西啊.

贊一下(1人)

回復(fù)此樓

Chemistry[]==[]Chem[]is[]try!!!

4樓2011-12-08 10:57:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ZDBWHZ

金蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 1882
紅花: 7
帖子: 736
在線: 228.7小時
蟲號: 655151
注冊: 2008-11-15
專業(yè): 理論和計算化學(xué)

最近不小心寫了個shell腳本，也是糾結(jié)了好久才完成。

贊一下

回復(fù)此樓

5樓2011-12-08 22:14:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋雨子

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時
蟲號: 1027042
注冊: 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

引用回帖:

4樓: Originally posted by abbott at 2011-12-08 10:57:14:
你先貼一個樣板文件,,,,
否則誰知道你要從什么文檔中提取你要的東西啊.

謝謝提醒~，貼個輸出文件，就是最簡單的優(yōu)化和頻率計算，因為優(yōu)化會計算好幾圈，想提取最后一次優(yōu)化的homo和lumo，但是找不到唯一區(qū)別的關(guān)鍵詞，不知道說清楚沒有，我是入門新手希望交流！

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : Row1.log

2011-12-08 22:46:50, 2.38 M

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

6樓2011-12-08 22:46:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秋雨子

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1578.8
散金: 13
帖子: 163
在線: 78.2小時
蟲號: 1027042
注冊: 2010-05-24
專業(yè): 理論環(huán)境化學(xué)

引用回帖:

5樓: Originally posted by ZDBWHZ at 2011-12-08 22:14:39:
最近不小心寫了個shell腳本，也是糾結(jié)了好久才完成。

是提取homo，lumo的嗎？可以分享下嗎？

我的shell幾乎是自學(xué)的，寫著很蹩腳，簡單還行，復(fù)雜就不行了

贊一下

回復(fù)此樓

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

7樓2011-12-08 22:49:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ZDBWHZ

金蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 1882
紅花: 7
帖子: 736
在線: 228.7小時
蟲號: 655151
注冊: 2008-11-15
專業(yè): 理論和計算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★
小紅豆(金幣+3): 感謝交流。 2011-12-14 18:08:11

fortran and shell，各有一個

贊一下(1人)

回復(fù)此樓

8樓2011-12-08 23:14:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 268求調(diào)劑 +7	好運連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 289求調(diào)劑 +6	步川酷紫123 2026-03-11	6/300	2026-03-17 10:23 by Sammy2
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655