24小时热门版块排行榜    

查看: 806  |  回复: 6

Gongmingxing

铜虫 (小有名气)

[求助] 温度效应的输入文件怎么编写

我现在想要编写一个三聚体在不同温度下的计算输入文件,可是我不知道要怎么编,有单体的例子,看不太明白,请各位高手指点,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsihang

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gmy1990: 金币+1 2013-11-01 18:56:06
Gongmingxing: 金币+1, 有帮助 2013-11-04 10:50:25
可以直接在命令行输入temperature=N。
2楼2013-11-01 16:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
Gongmingxing: 金币+1, 有帮助 2013-11-04 10:50:43
# b3lyp/6-31G    Temmperature=X (高斯默认单位是K不需要加单位)             opt
思想有多远,我们就能走多远~
3楼2013-11-04 09:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gongmingxing

铜虫 (小有名气)

我是要计算一系列不同的温度,不可能每一个文件都只是换一个温度再重新算一遍吧,有没有可能就是将所有的温度放在一个文件里,用link隔开呢?请指教!
4楼2013-11-04 10:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gongmingxing

铜虫 (小有名气)

我试了一下两个温度下的计算,还是算不过去,显示的是:
Base         20480
End          66048
End1         66048
Wr Pntr      20480
Rd Pntr      20480
Length       45568
Error termination in NtrErr:
NtrErr Called from FileIO.
输入文件如下:

%Chk=Azone
#MP2/6-31G Opt Freq=ReadIso

Azone

0 1
N                     3.97565  -3.69012  -0.7192
C                     2.90239  -3.06002  -0.20121
C                     3.34438  -2.37186   1.08169
C                     4.88323  -2.34593   0.97577
C                     5.21063  -3.56031   0.06286
H                     3.90109  -4.23759  -1.54899
O                     1.77868  -3.08191  -0.6834
H                     2.88465  -1.40058   1.13605
H                     3.01414  -2.97696   1.91977
H                     5.19523  -1.43467   0.4832
H                     5.37536  -2.41444   1.93392
H                     6.05786  -3.37343  -0.58001
H                     5.39675  -4.45402   0.64633

298.15 1.0 0.9804
14
12
12
12
12
1
16
1
1
1
1
1
1

--Link1--
%Chk=Azone
%NoSave
# MP2/6-31G(d,p) SP Guess=Read Geom=AllCheck Test

300.00 1.0 0.9804
14
12
12
12
12
1
16
1
1
1
1
1
1

--Link1--
%Chk=Azone
%NoSave
# MP2/6-31G(d,p) SP Guess=Read Geom=AllCheck Test
请各位高手指点!谢谢啦!
5楼2013-11-07 17:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
Gongmingxing: 金币+3, 有帮助 2013-11-11 19:24:15
引用回帖:
5楼: Originally posted by Gongmingxing at 2013-11-07 17:15:12
我试了一下两个温度下的计算,还是算不过去,显示的是:
Base         20480
End          66048
End1         66048
Wr Pntr      20480
Rd Pntr      20480
Length       45568
Error termination ...

输入文件看上去貌似没有错,但我一般不用readisotopes,而是像2楼所说用temperature=,scale=,等等,可以试试。
另外你的第二个计算命令行并没有指定用readisotopes,我不知道是不是可以这么用,为保险起见一般我都会再写一遍。

» 本帖已获得的红花(最新10朵)

6楼2013-11-08 07:11:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何飘同学

新虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by 枪下游魂 at 2013-11-08 07:11:10
输入文件看上去貌似没有错,但我一般不用readisotopes,而是像2楼所说用temperature=,scale=,等等,可以试试。
另外你的第二个计算命令行并没有指定用readisotopes,我不知道是不是可以这么用,为保险起见一般我 ...

用temperature=,得到的结果温度貌似被改变一点点,比如我设置的是:temperature=353.150,pressure=1.000
在热力学计算结构部分,显示的温度是:353.149, 当然得到的能量也上有所差别,请教是何原因呢?
7楼2013-11-12 09:38:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gongmingxing 的主题更新
信息提示
请填处理意见