版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

秋雨子

金蟲(chóng) (小有名氣)

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

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

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

回復(fù)此樓

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

搜一搜，收一收

» 猜你喜歡

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

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

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

重金求助HOMO/LUMO問(wèn)題已經(jīng)有5人回復(fù)
【解決】討論，siesta中作HOMO 與LUMO圖的方法已經(jīng)有12人回復(fù)
【求助】Gaussian計(jì)算HOMO-LUMO問(wèn)題已經(jīng)有17人回復(fù)
【求助】關(guān)于HOMO和LUMO軌道的確定問(wèn)題已經(jīng)有8人回復(fù)
black dye 的HOMO與LUMO 已經(jīng)有6人回復(fù)
【求助】計(jì)算homo和lumo軌道出錯(cuò) 已經(jīng)有14人回復(fù)
【求助】怎么看HOMO LUMO 已經(jīng)有11人回復(fù)
【求助】關(guān)于HOMO和LUMO的問(wèn)題已經(jīng)有5人回復(fù)
【求助】請(qǐng)問(wèn)，如何在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的回帖

asaka

銀蟲(chóng) (初入文壇)

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

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

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

說(shuō)一下我的想法�。�

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

然后把占據(jù)軌道的能級(jí)收集起來(lái)，取最大值，就是HOMO
以及把空軌道的能級(jí)收集起來(lái)，取最小值，就是LUMO
但是shell沒(méi)有數(shù)學(xué)計(jì)算的功能，只好借用一下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就可以運(yùn)行，測(cè)試LZ的文件，輸出如下：
$ sh homolumo Row1.log
File : Row1.log
HOMO : -0.40657
LUMO : 0.13145
其中homolumo是上面shell的文件名。

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

贊一下(2人)

回復(fù)此樓

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

秋雨子

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

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

查看全部 25 個(gè)回答

秋雨子

金蟲(chóng) (小有名氣)

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

沒(méi)有人回復(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的回帖

秋雨子

金蟲(chóng) (小有名氣)

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

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

abbott

金蟲(chóng) (著名寫(xiě)手)

不要用QQ問(wèn)我東西

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

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

查看全部 25 個(gè)回答

普通表情龍兔虎貓

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

[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[考研] 303求調(diào)劑 +3	睿08 2026-03-17	5/250	2026-03-17 22:01 by 睿08
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 材料與化工304求B區(qū)調(diào)劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學(xué)習(xí)賺彩禮
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty