计算化学公社

标题: gaussian加上电场书写格式问题 [打印本页]

作者
Author: dashuo1 时间: 2019-5-9 21:32
标题: gaussian加上电场书写格式问题
体系沿着X轴加上2*10^-8a.u.电场，根据gaussian使用说明，电场下优化，必须用 Opt=Z-Matrix NoSymm关键词，并用传统的Z矩阵坐标或带符号的笛卡尔坐标定义输入结构.我这里只用笛卡尔坐标，没用带符号的形式，请问各位高手，如下面这样编程序，可以吗。
%nprocshared=20
%mem=44gb
%chk=grae.chk
#p opt b3lyp/6-31g(d) nosymm field=read

grae

0 1
C                4.20961549 2.55178712 -4.02429112
C                3.00698664 1.83591275 -3.99799319
C                1.79115294 2.52724515 -3.95755066
C                3.01881940 0.43247437 -4.01201484
C                1.80705726 -0.27510072 -3.98733746
C                0.58396839 1.82092970 -3.93234583
C                0.58865799 0.42028330 -3.94796432
C                -0.62513265 -0.28389583 -3.92429653
C                -1.84154844 0.41484716 -3.88445758
C                -1.83912139 1.81835078 -3.86762437
C                -0.62789275 2.51845031 -3.89153699
C                -3.05930703 -0.28598376 -3.86148696
C                -4.26918793 0.43510940 -3.82133812
C                -4.26017116 1.83207812 -3.80422560
C                -3.04653031 2.52535966 -3.82698653
C                -5.48965476 -0.23801502 -3.79815951
C                -6.69622588 0.47521703 -3.75806305
C                -6.68743736 1.83269061 -3.74181168
C                -5.47202660 2.53309030 -3.76421179
C                6.63494288 2.57753280 -4.08965838
C                5.42791557 1.86751367 -4.06320200
C                5.44637432 0.47062772 -4.07570866
C                4.24117091 -0.25937238 -4.05066045
C                -3.05474781 -1.70012131 -3.87887065
C                -4.26002671 -2.43002565 -3.85721544
C                -5.48476545 -1.76576073 -3.81737966
C                -0.62055703 -1.68456237 -3.94046454
C                -1.83247988 -2.39212141 -3.91829595
C                1.81170113 -1.67570431 -4.00199212
C                0.59788188 -2.37997114 -3.97894777
C                4.24598064 -1.67347649 -4.06381017
C                3.02823978 -2.37440025 -4.03959323
C                7.85520866 1.88609152 -4.12739728
C                7.87320671 0.52858514 -4.13884172
C                6.67138276 -0.19347805 -4.11339766
C                6.67659462 -1.72125015 -4.12633610
C                5.45615281 -2.39447255 -4.10097205
C                -6.68667443 -2.48765429 -3.79547799
C                -4.24194336 -3.82690995 -3.87537683
C                -5.44917810 -4.53669654 -3.85370839
C                -6.66913558 -3.84497335 -3.81356487
C                -1.82090947 -3.79552350 -3.93557058
C                -3.02374252 -4.51132971 -3.91487734
C                0.60240643 -3.78068638 -3.99444055
C                -0.60510048 -4.48696671 -3.97361257
C                3.02587915 -3.77788917 -4.05297316
C                1.81448528 -4.47819812 -4.03077203
C                5.44733157 -3.79150700 -4.11281868
C                4.23350070 -4.48485830 -4.08866306
C                7.88332421 -2.43437137 -4.16302741
C                7.87486234 -3.79188869 -4.17402581
C                6.65948192 -4.49253407 -4.14914883
H                8.79931373 -4.32998904 -4.20182589
H                -7.58986951 -4.38970321 -3.79706112
H                -7.61177619 2.37094097 -3.71175307
H                -7.62704584 -0.05221608 -3.74025928
H                -7.62086784 -1.96696624 -3.76439262
H                8.81403458 -1.90677614 -4.18235887
H                8.80764477 0.00796370 -4.16741834
H                8.77582316 2.43119132 -4.14720724
H                -5.47089484 3.60300186 -3.75071318
H                -3.04054274 3.59517112 -3.81342618
H                -0.62854768 3.58833099 -3.87876156
H                1.78471820 3.59715600 -3.94588710
H                4.19644751 3.62161779 -4.01455130
H                6.62641556 3.64745175 -4.08106900
H                -5.44109125 -5.60655595 -3.86778067
H                -3.01057307 -5.58116545 -3.92921637
H                -0.59890341 -5.55690717 -3.98665213
H                1.81522892 -5.54815184 -4.04155619
H                4.22755680 -5.55477519 -4.09754215
H                6.65859877 -5.56249267 -4.15790037
空一行
0.0002 0.0000 0.0000
空二行

作者
Author: sobereva 时间: 2019-5-9 21:34
“必须用 Opt=Z-Matrix”这是手册里胡写，有过无数人问过这类问题。你的写法没问题

作者
Author: dashuo1 时间: 2019-5-9 21:36
十分感谢sobereva老师的解答。

作者
Author: mamie 时间: 2019-10-10 15:16
您好！想请教一下如果按照您的写法，是添加了“0.0002 0.0000 0.0000”的矢量方向，即电场（由正指向负）沿着这个矢量方向，大小为0.0002*10^-4au的电场是嘛？

作者
Author: sobereva 时间: 2019-10-11 21:22

mamie 发表于 2019-10-10 15:16
您好！想请教一下如果按照您的写法，是添加了“0.0002 0.0000 0.0000”的矢量方向，即电场（由正指向负）沿 ...

错。field=read时读取的矢量没有10^-4
方向没错

作者
Author: mamie 时间: 2019-10-12 01:48

sobereva 发表于 2019-10-11 21:22
错。field=read时读取的矢量没有10^-4
方向没错

谢谢老师！it helps a lot！

作者
Author: tangjinyu 时间: 2023-8-7 22:24

sobereva 发表于 2019-10-11 21:22
错。field=read时读取的矢量没有10^-4
方向没错

老师，想问您一下，如果我想在空间坐标系中某个点比如（2,2,2）施加+0.001a.u.的电场，那应该怎么写啊？

作者
Author: sobereva 时间: 2023-8-8 01:58

tangjinyu 发表于 2023-8-7 22:24
老师，想问您一下，如果我想在空间坐标系中某个点比如（2,2,2）施加+0.001a.u.的电场，那应该怎么写啊？

没法实现

欢迎光临计算化学公社 (http://bbs.keinsci.com/)