24小时热门版块排行榜    

查看: 1245  |  回复: 4

0502114073

木虫 (正式写手)

[求助] 高斯优化中,能否固定一部分,而只优化另外一部分呢?

比如石墨烯吸附CO,我怎么固定石墨烯,只优化CO分子,请教各位高手,谢谢!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

★ ★
zhou2009: 金币+2 2016-09-26 08:41:47
固定的部分可以写-1(在坐标和元素符号之间)
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
2楼2016-09-26 08:32:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

0502114073

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 小范范1989 at 2016-09-26 08:32:05
固定的部分可以写-1(在坐标和元素符号之间)

您的意思 直接在坐标和元素符号之间写-1,比如下面例子Fe(NH2)4  我要固定一个NH2就写成如下格式:

%chk=0007_H8N4Fe_neutral.chk
%mem=16GB
%nprocshared=16
# opt=gdiis freq b3lyp/6-311+g* geom=connectivity scf=(maxcycle=200,xqc)

Title Card Required

0 5
Fe                -0.91754600   -0.46831100    0.27157400
H                 -1.45244500   -0.76586600    1.83722800
N                  0.16000600    1.48253700    0.21099400
H                  1.14777100    1.24446300    0.16141600
H                 -3.16465200    1.08078500   -0.70341000
N                  4.18178500   -0.54864800   -0.29920800
H                  0.02178700    2.02434800    1.05949300
H                 -0.06401700   -1.31732100   -0.90479000
N                  3.37090400    0.17099000   -0.13735000
H                 -0.03924300    2.08668100   -0.58176900
H       -1        -3.54974800   -0.38967500   -0.09381000
N       -1       -2.87137700    0.10755200   -0.66468000
H       -1        -2.93247200   -0.27433300   -1.60356000

1
2
3 4 1.0 7 1.0 10 1.0
4
5 12 1.0
6 9 3.0
7
8
9
10
11 12 1.0
12 13 1.0
13
3楼2016-09-26 08:56:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
3楼: Originally posted by 0502114073 at 2016-09-26 08:56:59
您的意思 直接在坐标和元素符号之间写-1,比如下面例子Fe(NH2)4  我要固定一个NH2就写成如下格式:

%chk=0007_H8N4Fe_neutral.chk
%mem=16GB
%nprocshared=16
# opt=gdiis freq b3lyp/6-311+g* geom=connect ...

对,不过你的任务是不是想优化+freq?我感觉,计算的时候,表头写的越简单越好,而且你这个maxcycle=200压根就没有用。
我感觉就写opt freq就行。乱七八糟的命令不要写。除非出现的问题,调节收敛等等。个人见解

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

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
4楼2016-09-26 09:25:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

0502114073

木虫 (正式写手)

送红花一朵
引用回帖:
4楼: Originally posted by 小范范1989 at 2016-09-26 09:25:35
对,不过你的任务是不是想优化+freq?我感觉,计算的时候,表头写的越简单越好,而且你这个maxcycle=200压根就没有用。
我感觉就写opt freq就行。乱七八糟的命令不要写。除非出现的问题,调节收敛等等。个人见解...

好的,谢谢,我优化过程中总是出现不收敛,这个是出现了震荡 所以写成这样形式,我随便拿了一个当例子的,受教了。非常感谢!!!怎么给你金币?
5楼2016-09-26 09:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 0502114073 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 北京林业大学硕导招生广告 +3 kongweilin 2026-03-26 3/150 2026-03-26 00:53 by 15723253818
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 308求调剂 +5 墨墨漠 2026-03-25 5/250 2026-03-25 22:19 by 544594351
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 321求调剂 +3 璞玉~~ 2026-03-25 3/150 2026-03-25 19:07 by Zhanglab-TJU
[考研] 329求调剂 +3 钮恩雪 2026-03-25 3/150 2026-03-25 14:43 by 糖加冰
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 081700 调剂 267分 +10 迷人的哈哈 2026-03-23 10/500 2026-03-25 12:12 by userper
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
信息提示
请填处理意见