计算化学公社

标题: 怎么获得相关的输入文件来计算相应分子的偶极矩 [打印本页]

作者
Author: feiwang3am 时间: 2018-7-19 17:11
标题: 怎么获得相关的输入文件来计算相应分子的偶极矩

Sobereva在
通过Multiwfn计算各个轨道的偶极矩一文中说明了怎么用Multiwfn来计算H2O的偶极矩。该输入文件是现成范例文件的。我的问题是如果计算其他分子的偶极矩，怎么获得这个输入文件？比如

C18H35NH2 （油胺），CH3(CH2)11NH2 (十二胺)。谢谢，在线等！

作者
Author: sobereva 时间: 2018-7-19 17:41
Gaussian里优化，末尾直接就输出了偶极矩
算偶极矩的基组至少用def2-SVPD，否则精度很烂

作者
Author: feiwang3am 时间: 2018-7-19 18:32

sobereva 发表于 2018-7-19 17:41
Gaussian里优化，末尾直接就输出了偶极矩
算偶极矩的基组至少用def2-SVPD，否则精度很烂

谢谢社长

作者
Author: Chenglong_li 时间: 2018-11-26 15:59
Sob老师好！我要计算有机分子的偶极矩，使用PBE0/def2-TZVPD来进行，输入文件如下，请问是否妥当？多谢

%chk=DAPCA-LYS_PBE0_Output.chk
#p opt freq pbe1pbe/gen
[空行]
Opt_DAPCA-LYS_PBE0
[空行]
0 1
C                -5.63074306 0.41933764 -0.47202404
N                -4.86553578 1.50659039 -0.26747676
C                -3.60358476 1.24690285 0.02991358
C                -3.05642979 -0.02263543 0.13554665
C                -3.95002204 -1.08858718 -0.17485804
N                -5.23247304 -0.84930595 -0.45872641
H                -2.97741424 2.12406199 0.16962406
N                -6.94284758 0.64836570 -0.71134478
H                -7.51943524 -0.11878318 -1.00775126
H                -7.23803291 1.59543572 -0.86955677
N                -3.55558850 -2.36947403 -0.18006429
H                -4.23174416 -3.06374046 -0.44795464
H                -2.57020126 -2.57607557 -0.09177886
C                -1.64484672 -0.28089587 0.49175093
O                -1.09157385 -1.35188443 0.25274881
N                -0.97102089 0.71926313 1.11518440
C                0.42961100 0.56421666 1.47404304
H                -1.46980804 1.53306062 1.42837180
C                1.36861741 0.78814779 0.29383573
C                2.82699450 0.57752011 0.68493294
H                1.22324599 1.80065074 -0.09515498
H                1.09191292 0.09088660 -0.50055465
H                2.95913157 -0.43909333 1.05847080
H                3.08696507 1.25505374 1.50718523
C                5.26424766 0.68220224 -0.16292334
H                3.63198384 1.85688570 -0.83898871
H                3.51992937 0.18282232 -1.32864668
N                6.06326519 1.44448145 -1.10775807
H                5.41909592 0.98558603 0.88654747
C                5.75271000 -0.75317543 -0.21483974
H                7.02341624 1.12062383 -1.11203774
O                6.81594971 -1.09598084 -0.65279196
O                4.88363618 -1.62582650 0.32328152
H                5.28901020 -2.50371158 0.27001718
[空行]
H    0
S 3 1.00
   34.061341000          0.60251978000E-02
   5.1235746000          0.45021094000E-01
   1.1646626000          0.20189726000
S 1 1.00
   0.32723041000       1.0000000
S 1 1.00
   0.10307241000       1.0000000
P 1 1.00
   0.80000000000       1.0000000
P 1 1.00
   0.95774129632E-01          1.0000000
****
C    0
S 6 1.00
  13575.3496820             0.22245814352E-03
2035.2333680             0.17232738252E-02
463.22562359          0.89255715314E-02
131.20019598          0.35727984502E-01
   42.853015891          0.11076259931
   15.584185766          0.24295627626
S 2 1.00
   6.2067138508          0.41440263448
   2.5764896527          0.23744968655
S 1 1.00
   0.57696339419       1.0000000
S 1 1.00
   0.22972831358       1.0000000
S 1 1.00
   0.95164440028E-01          1.0000000
S 1 1.00
   0.48475401370E-01          1.0000000
P 4 1.00
   34.697232244          0.53333657805E-02
   7.9582622826          0.35864109092E-01
   2.3780826883          0.14215873329
   0.81433208183       0.34270471845
P 1 1.00
   0.28887547253       1.0000000
P 1 1.00
   0.10056823671       1.0000000
D 1 1.00
   1.0970000000          1.0000000
D 1 1.00
   0.31800000000       1.0000000
D 1 1.00
   0.90985336424E-01          1.0000000
F 1 1.00
   0.76100000000       1.0000000
****
N    0
S 6 1.00
  19730.8006470             0.21887984991E-03
2957.8958745             0.16960708803E-02
673.22133595          0.87954603538E-02
190.68249494          0.35359382605E-01
   62.295441898          0.11095789217
   22.654161182          0.24982972552
S 2 1.00
   8.9791477428          0.40623896148
   3.6863002370          0.24338217176
S 1 1.00
   0.84660076805       1.0000000
S 1 1.00
   0.33647133771       1.0000000
S 1 1.00
   0.13647653675       1.0000000
S 1 1.00
   0.68441605847E-01          1.0000000
P 4 1.00
   49.200380510          0.55552416751E-02
   11.346790537          0.38052379723E-01
   3.4273972411          0.14953671029
   1.1785525134          0.34949305230
P 1 1.00
   0.41642204972       1.0000000
P 1 1.00
   0.14260826011       1.0000000
D 1 1.00
   1.6540000000          1.0000000
D 1 1.00
   0.46900000000       1.0000000
D 1 1.00
   0.12829642058       1.0000000
F 1 1.00
   1.0930000000          1.0000000
****
O    0
S 6 1.00
  27032.3826310             0.21726302465E-03
4052.3871392             0.16838662199E-02
922.32722710          0.87395616265E-02
261.24070989          0.35239968808E-01
   85.354641351          0.11153519115
   31.035035245          0.25588953961
S 2 1.00
   12.260860728          0.39768730901
   4.9987076005          0.24627849430
S 1 1.00
   1.1703108158          1.0000000
S 1 1.00
   0.46474740994       1.0000000
S 1 1.00
   0.18504536357       1.0000000
S 1 1.00
   0.70288026270E-01          1.0000000
P 4 1.00
   63.274954801          0.60685103418E-02
   14.627049379          0.41912575824E-01
   4.4501223456          0.16153841088
   1.5275799647          0.35706951311
P 1 1.00
   0.52935117943       1.0000000
P 1 1.00
   0.17478421270       1.0000000
P 1 1.00
   0.51112745706E-01          1.0000000
D 1 1.00
   2.3140000000          1.0000000
D 1 1.00
   0.64500000000       1.0000000
D 1 1.00
   0.14696477366       1.0000000
F 1 1.00
   1.4280000000          1.0000000
****
[空行]
[空行]

作者
Author: sobereva 时间: 2018-11-26 16:50

Chenglong_li 发表于 2018-11-26 15:59
Sob老师好！我要计算有机分子的偶极矩，使用PBE0/def2-TZVPD来进行，输入文件如下，请问是否妥当？多谢

...

输入文件没问题。优化和振动分析没必要带弥散，当然如果你不嫌耗时高也可以

作者
Author: Chenglong_li 时间: 2018-11-26 18:58

sobereva 发表于 2018-11-26 16:50
输入文件没问题。优化和振动分析没必要带弥散，当然如果你不嫌耗时高也可以

多谢sob老师。
我这个要计算偶极矩，因几何优化的结果中就有偶极矩，所以就直接带弥散进行优化了。
我也是看了《谈谈量子化学中基组的选择》博文，文中说“偶极矩、极化率、第一超极化率，算这三个必须带弥散函数，否则结果烂得根本没法用”，所以几何优化的时候带了弥散。

作者
Author: alonewolfyang 时间: 2018-11-26 20:53

Chenglong_li 发表于 2018-11-26 18:58
多谢sob老师。
我这个要计算偶极矩，因几何优化的结果中就有偶极矩，所以就直接带弥散进行优化了。
我 ...

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