应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » DFT Scan之后使用MP2计算能量，如何建立批量输入文件？
71/1返回列表
查看: 2403  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ChemiAndy

木虫 (正式写手)

[求助] DFT Scan之后使用MP2计算能量，如何建立批量输入文件？

如题，使用B3LYP Scan优化了很多构型，然后需要使用MP2依次计算所有构型的结合能。由于点很多，一个一个构型建立输入文件并运行太麻烦。请问有没有简单的办法建立一个输入文件实现上述功能？谢谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-12-29 04:32:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZDBWHZ

金虫 (正式写手)
去学学shell吧
回复此楼
2楼2011-12-30 00:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
谢谢，正在搞Shell
赞一下 回复此楼
3楼2011-12-31 01:28:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeozhou

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
ChemiAndy(金币+20): ★★★很有帮助 好主意，我试试 2012-01-02 14:04:16
可以尝试在输入中使用
----link1---
将构型隔开。然后，可以把所有构型都放到一个Input中
注意mp2通常chk文件比较大，如果担心输出文件太大的话，可以均分为几个输入文件完成

不知道这个方法适不适合你？
赞一下 回复此楼
4楼2011-12-31 17:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★ ★ ★ ★ ★
小红豆(金币+5): 感谢分享 2012-01-16 21:43:23
写了一个bash脚本（只能在linux, bash类型的shell下运行），供参考：

下面的脚本把高斯MP2单点Scan(注意不是Scan优化)生成的输出文件中的结构全部读出，生成DFTB所需的输入文件。你可以修改它实现其它目的。
引用回帖:

#!/bin/bash

#

# 1. Extract 21 configurations from scan-calculation of gaussian

# 2. Create 21 DFTB 'gen' files

# 3. Create 21 tinker 'xyz' files (TBA)

#

# Xijun Wang, Jan. 5, 2012

# xijun(at)cermm.concordia.ca

# Ref:1.Cu G. Phung, Gaustool, http://www.ccl.net/cca/software/UNIX/gaussian-utilities/README



MOL='g-w'   // here you name your system

NATOM=13    // here you control how many atoms in the system

NMOL=21    //  here you control how many configurations will be read from the gaussian file

echo "MOL name is $MOL, NMOL is $MOL, NATOM is $NATOM"



rm $MOL*   // delete residue files of previous running



# Extract configurations from gaussian output file;

# Only Z-matrix orientation coordinates will be extracted;

sed '/Z-Matrix orientation/,/Distance/!d' $1 >> temp1

echo "Gaussian structure has been extracted ..."



# Delete un-useful lines

grep -v 'Num' temp1|grep -v 'Dis'|grep -v '\-\-\-'|grep -v 'Z'|grep -v 'Center' >> temp2

echo "Some lines have been deleted ..."



# Replace element number with their types: 6-C; 7-N; 1-H; 8-O;

# Examine the gaussian file for details of the replace below:

sed -i -e 's/6           0 /C /g' -e 's/7           0 /N /g' \

       -e 's/1           0 /H /g' -e 's/8           0 /O /g'  temp2

echo "Element number has been replaced by element name ..."



# change element type: C-1; N-2; H-3; O-4;

sed -i -e 's/C/1/g' -e 's/N/2/g' -e 's/H/3/g' -e 's/O/4/g' temp2

echo "Element name have been replaced by number in order of C, N, H, O ..."



# Separate temp2 into 21 files with name g-w.xx.gen

for ((IMOL=0;IMOL<$NMOL;IMOL++)); do

  let  m=13*$IMOL+1

  let  n=$m+$NATOM-1

  let  fn=$IMOL+1

# Add the first two lines to the gen file

  echo "  $NATOM C" > $MOL.$fn.gen   

  echo "  C N H O" >> $MOL.$fn.gen

  echo "Print from lines $m to line $n into file $MOL.$fn.gen..."

  sed -n''"$m"','"$n"'p' temp2 >> $MOL.$fn.gen

done

rm temp*

exit



This script only works for gausian output of single-point scan calculations. PAY ATTENTION that it does not work for optimization calc. since all structures during the optimization will be output.

Use of the script: Save it with name like g2gen.sh

$ chmod +x g2gen.sh           // make it executable

$ g2gen.sh gaussianfile.out      // replace "gaussianfile.out" with your gaussian output file

Send me email if you have problem.
[/FONT]

http://hi.baidu.com/xijunw/blog/ ... 6efb5dfbf2c01f.html
赞一下(1人) 回复此楼
5楼2012-01-16 10:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★ ★ ★ ★ ★
小红豆(金币+5): 感谢参与应助 2012-01-16 21:43:36
另外一个更加通用点的脚本

该脚本能够读取Gaussian03/09输出文件中的最后一个结构，生成xyz坐标文件，然后利用openBabel转换成Gaussian输入文件。因此使用前请确保安装OpenBabel先。

1. 生成的gaussian输入文件为Cartessian坐标。欲得到Z-Matrix坐标的输入文件，可修改Babel行为“babel -ixyz $name.xyz -ogzmat $name.gjf”

2. 由于需要替换原子序数为原子元素名称，而这里只转换C,H,N,O,Cl等5个元素，因此，你的体系如果包含其它元素，请仿照原格式修改其中的“# Change the numbers into elements names”部分，其实就是在那一句中添加gsub(元素序号,"元素名称",$2);即可（顺序无关），别忘了;号。其实元素序号/名称转换也可以从输出文件中提取。只是暂没时间完成这一步。
引用回帖:


#!/bin/bash
#
# Xijun wang, 2012.01.14
# xijun(at)gmail.com

# get the prefix of the gaussian output file
name=`echo $1|sed "s/\..*$//"`

# Find charge and multiplicity, get the first value in case multiple results found
charge=`awk '/Charge =/{print $3}' $1|awk 'NR==1'`
multip=`awk '/Charge =/{print $6}' $1|awk 'NR==1'`

# Find the block of the coordination, delete extra lines
sed '/Input orientation/,/Distance/!d' $1| grep -v -E 'I|C|N|D|\-\-\-' > $name.temp.01

# Find how many atoms by reading the first string of the last line
NN=`awk 'END{print $1}' $name.temp.01`

# Output the last structure from all structures found in the output file
tail -$NN $name.temp.01 > $name.temp.02

# Change the numbers into elements names, to be inproved
awk '{gsub(1,"H",$2); gsub(6,"C",$2); gsub(7,"N",$2); gsub(8,"O",$2); gsub(17,"Cl",$2); \
      printf "%5s %10.6f %10.6f %10.6f\n", $2, $4, $5, $6}' $name.temp.02 > $name.temp.03

# Delete previous resulting file
rm $name.xyz $name.gjf

# Create xyz file
echo $NN >> $name.xyz
echo $name >> $name.xyz
cat $name.temp.03 >> $name.xyz

# Transform xyz file into gaussian input file using openBabel
babel -ixyz $name.xyz -ogjf $name.gjf

# Replace some content of the input files
sed '1 i %chk='"$name.chk"'\n%mem=1024MB' $name.gjf > $name.temp.gjf
sed -i 's/#Put Keywords Here, check Charge and Multiplicity./#MP2\/6-311++g** SCRF=(PCM, solvent=water, DoVacuum, ExternalIteration, SMD)/' $name.temp.gjf
sed -i "s/0  1/$charge $multip/" $name.temp.gjf

# Rename to get the final file and delete temp files
mv $name.temp.gjf $name.gjf
rm $name.temp.*

Use of the script:

save it as g2g.sh; add excutable: chmod +x g2g.sh

then: ./g2g gauss.out // replace gauss.out with your output file name.
赞一下(1人) 回复此楼
6楼2012-01-16 10:24:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
上面第二个脚本生成的gaussian输入文件所带的预设关键词是用PCM模型计算溶剂化自由能的。计算后直接读G(Sol)那一行的数值。
赞一下 回复此楼
7楼2012-01-16 10:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ChemiAndy 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 19:34 by lXgU0RvNTDC9
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:34 by lXgU0RvNTDC9
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:22 by lXgU0RvNTDC9
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 18:55 by lXgU0RvNTDC9
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 10/500 2026-08-23 07:17 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭