gaussian加上电场书写格式问题

dashuo1

体系沿着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
空二行

sobereva

“必须用 Opt=Z-Matrix”这是手册里胡写，有过无数人问过这类问题。你的写法没问题

dashuo1

十分感谢sobereva老师的解答。

mamie

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

sobereva

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

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

mamie

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

谢谢老师！it helps a lot！

tangjinyu

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

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

sobereva

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

没法实现