24小时热门版块排行榜    

查看: 1183  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见