24小时热门版块排行榜    

查看: 590  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wkfg

新虫 (初入文坛)


[交流] Pd-DPPF螯合物Gaussian结构优化

Gaussian小白求助
最近想优化DPPF和Pd(OAc)2的螯合物,查了文献发现要用混合基组,但是这一部分没有接触过。查了一下小木虫其他人做的写出了以下输入文件,不知道对不对。
因为混合机组有人写的是b3lyp/gen,还有人写b3lyp/genecp或者b3lyp/gen pseudo=read,目前写成第一种可以计算下去,但是还没有算完。希望大神指点一下这些写法有什么区别,同时看一下我这么写输入文件对不对。本来想放一下gauss view 画的结构,但是点击上传图片没有反应,实在是对不住了。


%chk=D:/Cai/gaussian/Furan Artical/Pd-DPPF.chk
%mem=1GB
%nproc=1

# opt freq b3lyp/gen

title

0 1

C
C                  1              B1
C                  2              B2    1              A1
C                  1              B3    2              A2    3              D1
C                  4              B4    1              A3    2              D2
C                  5              B5    4              A4    1              D3
C                  6              B6    5              A5    4              D4
C                  6              B7    5              A6    4              D5
C                  7              B8    6              A7    5              D6
C                  8              B9    6              A8    5              D7
Fe                 9             B10    7              A9    6              D8
Pd                 1             B11    2             A10    3              D9
O                 12             B12    1             A11    2             D10
O                 12             B13    1             A12    2             D11
C                 13             B14   12             A13    1             D12
O                 15             B15   13             A14   12             D13
C                 14             B16   12             A15    1             D14
O                 17             B17   14             A16   12             D15
C                 15             B18   13             A17   12             D16
H                 19             B19   15             A18   13             D17
H                 19             B20   15             A19   13             D18
H                 19             B21   15             A20   13             D19
C                 19             B22   15             A21   13             D20
H                 23             B23   19             A22   15             D21
H                 19             B24   15             A23   13             D22
H                 19             B25   15             A24   13             D23

   B1             1.48193120
   B2             1.38379624
   B3             1.48146474
   B4             1.38372071
   B5             3.43192569
   B6             1.38224965
   B7             1.48650858
   B8             1.47915197
   B9             1.38292968
   B10            2.15774562
   B11            2.17519582
   B12            1.94000000
   B13            1.94000000
   B14            1.43000000
   B15            1.25840000
   B16            1.43000000
   B17            1.25840000
   B18            1.54000000
   B19            1.07000000
   B20            1.07000000
   B21            1.07000000
   B22            3.33290692
   B23            1.07000000
   B24            2.60277666
   B25            3.95106972
   A1           117.95612818
   A2            93.09498549
   A3           118.27987077
   A4            80.51882551
   A5            80.22491113
   A6            90.19610491
   A7           117.99822155
   A8           103.93841516
   A9            59.23246337
   A10          115.47097779
   A11          113.87300177
   A12          116.48312337
   A13          109.47122063
   A14          120.00000000
   A15          109.47122063
   A16          120.00000000
   A17          120.00000000
   A18          109.47122063
   A19          109.47122063
   A20          109.47122063
   A21           70.52936823
   A22          129.58891871
   A23           82.21695687
   A24           80.42015143
   D1            16.44059389
   D2           -16.50729926
   D3            99.12701850
   D4             0.70075937
   D5          -103.45690189
   D6           -99.49752318
   D7            80.04534697
   D8            67.52764485
   D9           136.28990203
   D10         -174.74631147
   D11           65.49379626
   D12         -130.18186358
   D13           90.00000000
   D14           89.56736649
   D15          -30.00000000
   D16          -90.00000000
   D17           90.00000000
   D18         -150.00000000
   D19          -30.00000000
   D20          131.37680951
   D21           67.05203180
   D22          141.49577983
   D23          121.38559077

C 0
6-31g(d,p)
****
H 0
6-31g(d,p)
****
O 0
6-31g(d,p)
****
Fe 0
6-31g(d,p)
****

Pd 0
lanl2dz
****
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wkfg

新虫 (初入文坛)


送红花一朵
引用回帖:
4楼: Originally posted by paramecium86 at 2022-05-05 22:23:54
混合基组写法有误。
1. 首先你的体系中含有重金属。那么一定要带有赝势的,所以要么是genecp 要么是 gen pseudo=read

2. 赝势基组包括两部分 一个是基组部分 一个赝势部分都要书写

3. 基组部分 写在前面  每 ...

懂了,谢谢版主,刚刚查了一下才知道gen是混合基组 genecp和gen pseudo=read才是混合+赝势
6楼2022-05-06 09:49:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wkfg(金币+2): 谢谢参与
wkfg: 金币+10 2022-05-06 09:50:41
混合基组写法有误。
1. 首先你的体系中含有重金属。那么一定要带有赝势的,所以要么是genecp 要么是 gen pseudo=read

2. 赝势基组包括两部分 一个是基组部分 一个赝势部分都要书写

3. 基组部分 写在前面  每一个基组后要用**** 即使是最后一个基组也要写****, 不同基组之间不要空行。且使用相同基组的原子可以写在一起。

4. 赝势和基组部分之间要有空行

5. 赝势部分不要写****



所以综上所述你的输入文件应该写成
%chk=d:/cai/gaussian/furan artical/pd-dppf.chk
%mem=1gb
%nproc=1

# opt freq b3lyp/genecp

title

0 1

c
c                  1              b1
c                  2              b2    1              a1
c                  1              b3    2              a2    3              d1
c                  4              b4    1              a3    2              d2
c                  5              b5    4              a4    1              d3
c                  6              b6    5              a5    4              d4
c                  6              b7    5              a6    4              d5
下面坐标我省略了

pd 0
lanl2dz
****
c h o fe 0
6-31g(d,p)
****

pd 0
lanl2dz

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

4楼2022-05-05 22:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 刚刚收到科研之友邮件 +9 olivermiaoer 2024-06-19 11/550 2024-06-20 06:49 by kobe0107
[找工作] 药学硕士找不到工作,打算去做科研助理了 +9 pom戴墨镜 2024-06-14 22/1100 2024-06-20 05:47 by 投必得科研顾问
[考博] 这个博士要读吗 +14 Sea Breeze 2024-06-16 23/1150 2024-06-20 05:41 by 投必得科研顾问
[基金申请] 面青地会评时间??? +5 Axvdvbfs 2024-06-19 5/250 2024-06-19 21:11 by 人民海军
[基金申请] 2024国社科通讯评审 +21 qsd10086 2024-06-13 43/2150 2024-06-19 20:57 by 05111401022
[基金申请] 75批博后基金 +11 kyukitu 2024-06-13 14/700 2024-06-19 18:44 by Siuhong Yiu
[基金申请] 面上博士后 +3 jsqy 2024-06-19 3/150 2024-06-19 17:58 by bxwisdom_886
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +8 babu2015 2024-06-14 8/400 2024-06-19 15:25 by 风今25
[论文投稿] 求机械类四区sci推荐 5+4 迷茫小旷 2024-06-14 5/250 2024-06-19 14:08 by tangjie12345
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+4 BruceChum 2024-06-15 22/1100 2024-06-19 08:00 by kanyechris
[公派出国] CSC德国博后每个月资助多少呀?够用吗 +4 326lhpqk 2024-06-16 7/350 2024-06-19 02:03 by PLHOU
[教师之家] 请问事业编制和年薪制冲突吗? +12 ZHONGWU_U 2024-06-14 12/600 2024-06-18 19:31 by fangyl2005
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 6/300 2024-06-18 14:15 by 茕茕恭煮
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[硕博家园] 博士毕业高校和就业的相关问题 +7 SCITOPPP 2024-06-14 11/550 2024-06-18 07:51 by yinxing1995
[论文投稿] 投稿被一个审稿人恶意评审了怎么样? +5 1chen 2024-06-14 7/350 2024-06-15 23:15 by xy66xy
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
信息提示
请填处理意见