24小时热门版块排行榜    

查看: 2641  |  回复: 10
【悬赏金币】回答本帖问题,作者lirunhao将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lirunhao

铜虫 (初入文坛)

[求助] 求大神帮忙看一下 怎么解决这样高斯计算问题已有2人参与

简单说我是想计算 2个有机分子的堆积结构优化  主要是想得到2个分子间的HOMO LUMO   分子结构模型我是用Hg软件打开单晶文件然后根据弱相互作用力 点出的第二个分子   

我用的是混合基组  
# opt b3lyp/gen pseudo=read  

title

0 1
Au                 4.42900000    0.54800000    7.17180000
C                  5.60410000   -0.70960000    8.24600000
C                  5.79210000   -0.55160000    9.59080000
C                  6.70200000   -1.30070000   10.32030000
C                  7.41090000   -2.28340000    9.68540000
C                  7.20380000   -2.48930000    8.34850000
C                  6.33390000   -1.72640000    7.65450000
C                  3.34470000    1.85120000    6.14010000
C                  1.96080000    3.62460000    5.05570000
C                  1.22410000    4.39390000    5.94690000
H                  1.27000000    4.25630000    6.86580000
C                  0.40240000    5.39050000    5.39480000
H                 -0.11190000    5.92510000    5.95870000
C                  0.35460000    5.58040000    4.03820000
C                  1.09840000    4.81860000    3.18640000
H                  1.06100000    4.97330000    2.26760000
C                  1.91820000    3.80390000    3.68720000
H                  2.41970000    3.26610000    3.11940000
C                 -1.23550000    7.35800000    4.27480000
H                 -0.70460000    7.90650000    4.87030000
H                 -1.82370000    6.79680000    4.80330000
C                 -2.03550000    8.22340000    3.32050000
H                 -2.44810000    7.65470000    2.65400000
H                 -2.74860000    8.65870000    3.81740000
C                 -0.97000000   10.31680000    3.49400000
H                 -0.53570000   10.02130000    4.31030000
H                 -1.78350000   10.78630000    3.73060000
C                 -0.06340000   11.21310000    2.72900000
H                  0.29250000   11.89710000    3.31650000
H                  0.68180000   10.70090000    2.37800000
C                  0.12330000   12.43480000    0.73350000
H                  0.70450000   11.76470000    0.34310000
H                  0.67810000   13.08890000    1.18700000
C                 -0.66700000   13.08890000   -0.31550000
H                 -1.33120000   12.46580000   -0.64670000
H                 -0.08420000   13.32580000   -1.05290000
Cl                -1.49870000   14.57310000    0.27760000
F                  5.07030000    0.35640000   10.25920000
F                  6.86960000   -1.11390000   11.61300000
F                  8.27780000   -3.00360000   10.39130000
F                  7.91670000   -3.44210000    7.72550000
F                  6.21790000   -1.93340000    6.33650000
N                  2.75500000    2.60560000    5.59990000
O                 -0.38360000    6.55030000    3.43090000
O                 -1.26630000    9.21030000    2.67140000
O                 -0.76250000   11.82240000    1.66420000
Au                -0.22380000    0.54800000    7.17180000
C                  0.95130000   -0.70960000    8.24600000
C                  1.13930000   -0.55160000    9.59080000
C                  2.04920000   -1.30070000   10.32030000
C                  2.75810000   -2.28340000    9.68540000
C                  2.55100000   -2.48930000    8.34850000
C                  1.68110000   -1.72640000    7.65450000
C                 -1.30810000    1.85120000    6.14010000
C                 -2.69200000    3.62460000    5.05570000
C                 -3.42870000    4.39390000    5.94690000
H                 -3.38280000    4.25630000    6.86580000
C                 -4.25040000    5.39050000    5.39480000
H                 -4.76470000    5.92510000    5.95870000
C                 -4.29820000    5.58040000    4.03820000
C                 -3.55440000    4.81860000    3.18640000
H                 -3.59180000    4.97330000    2.26760000
C                 -2.73460000    3.80390000    3.68720000
H                 -2.23310000    3.26610000    3.11940000
C                 -5.88830000    7.35800000    4.27480000
H                 -5.35740000    7.90650000    4.87030000
H                 -6.47650000    6.79680000    4.80330000
C                 -6.68830000    8.22340000    3.32050000
H                 -7.10090000    7.65470000    2.65400000
H                 -7.40140000    8.65870000    3.81740000
C                 -5.62280000   10.31680000    3.49400000
H                 -5.18850000   10.02130000    4.31030000
H                 -6.43630000   10.78630000    3.73060000
C                 -4.71620000   11.21310000    2.72900000
H                 -4.36030000   11.89710000    3.31650000
H                 -3.97100000   10.70090000    2.37800000
C                 -4.52950000   12.43480000    0.73350000
H                 -3.94830000   11.76470000    0.34310000
H                 -3.97470000   13.08890000    1.18700000
C                 -5.31980000   13.08890000   -0.31550000
H                 -5.98400000   12.46580000   -0.64670000
H                 -4.73700000   13.32580000   -1.05290000
Cl                -6.15150000   14.57310000    0.27760000
F                  0.41750000    0.35640000   10.25920000
F                  2.21680000   -1.11390000   11.61300000
F                  3.62500000   -3.00360000   10.39130000
F                  3.26390000   -3.44210000    7.72550000
F                  1.56510000   -1.93340000    6.33650000
N                 -1.89780000    2.60560000    5.59990000
O                 -5.03640000    6.55030000    3.43090000
O                 -5.91910000    9.21030000    2.67140000
O                 -5.41530000   11.82240000    1.66420000

C H N O F Cl 0
6-31G*
****
Au 0
LANL2DZ
****

Au 0
LANL2DZ

最后报错L502 如图。。。   
  请各位专家帮帮忙  应该如何解决

求大神帮忙看一下 怎么解决这样高斯计算问题
1.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

引用回帖:
10楼: Originally posted by lirunhao at 2020-05-17 21:26:43
老师  那我现在的情况下  您觉得 是应该用wb97xd方法吗   关键词是这样吗 # opt wb97xd/genecp scf(vshift=500,xqc)...

我觉得可以试试看    或者你用 b3lyp 不加d3 也可以算算看
11楼2020-05-18 10:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

天天进步啊

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
1、按楼主所说,是优化2个分子的结构,那么得到的HOMO,LUMO只可能是每一个分子的,不存在2个分子之间的HOMO,LUMO这一说法。
2、如果这2个分子之间有化学键,那么则视为一个分子,如果没有,可以视为2个片段,要在输入文件中标明每一个原子是属于哪个片段。
3、楼主不应该用手机拍照输出文件,而应该直接上传输出文件,原因有2个:一是看不清,二是可能真正出错的原因在前面
分子筛、分子动力学模拟,有问题可咨询我的知乎主页 https://www.zhihu.com/people/rao777
2楼2020-05-17 21:10:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nilman

金虫 (著名写手)

引用回帖:
2楼: Originally posted by 天天进步啊 at 2020-05-17 21:10:10
1、按楼主所说,是优化2个分子的结构,那么得到的HOMO,LUMO只可能是每一个分子的,不存在2个分子之间的HOMO,LUMO这一说法。
2、如果这2个分子之间有化学键,那么则视为一个分子,如果没有,可以视为2个片段,要在输 ...

工整合理,逻辑清楚,赞!

发自小木虫Android客户端
3楼2020-05-17 21:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
从你的上传的截图中,勉强辨认出这个是SCF不收敛错误 。 下次确实应该直接上传log文件好一些。或者用截图功能。

这个要算两个分子的dimmer complex这没问题。但是要做好准备这个结构直接摆上就算。 一开始SCF是肯定不好收敛的。需要加关键词帮助收敛才行。这个分子不但原子众多 而且包含重金属。

一些常用的操作见我下面截图中的内容。
求大神帮忙看一下 怎么解决这样高斯计算问题-1
QQ截图20200517102603.png

4楼2020-05-17 22:27:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 论大家对6070后普通教授导师的看法 +7 SNaiL1995 2024-05-28 11/550 2024-06-02 13:16 by 小虫子咔咔
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +33 kaoyan250 2024-05-27 44/2200 2024-06-02 13:06 by 苏东坡二世
[基金申请] 讨论贴:好本子一定会被打招呼的淘汰吗 +40 ioe3937 2024-05-27 50/2500 2024-06-02 12:12 by ECOblance
[基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 5+3 宜兰 2024-05-31 9/450 2024-06-02 09:42 by Xiaolin81
[基金申请] 化学口B0109(高分子合成),拿青年基金一般需要怎样的文章水平? +13 salmon95 2024-05-30 22/1100 2024-06-01 21:22 by salmon95
[考研] 研0二导师分到新来的博士后靠谱吗 +7 sone9 2024-05-31 7/350 2024-06-01 19:10 by 梦燕园
[硕博家园] 各位同学能否分享一下实验室的学生劳务发放标准呀? +18 ma3252788 2024-05-30 18/900 2024-06-01 17:02 by 落雨风霜
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[论文投稿] 纠结选哪一个期刊,电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[高分子] MMA预聚体光固化发雾问题求助 +3 惠亚金总 2024-05-29 10/500 2024-05-31 14:59 by 惠亚金总
[有机交流] 液相纯度高,但产品析不太出来 10+4 cui19236 2024-05-27 9/450 2024-05-30 07:45 by yuanjijoy
[论文投稿] 《控制理论与应用》期刊收版面费吗? +4 ygj2015 2024-05-29 4/200 2024-05-30 06:59 by nono2009
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[博后之家] 2024公派博后申请 +4 326lhpqk 2024-05-27 5/250 2024-05-29 20:03 by @古月胡
[基金申请] 信息学部函评结束了吗? +6 ducan21 2024-05-28 7/350 2024-05-29 12:10 by WORLD0256
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[有机交流] 机理求助 120+4 15147165026 2024-05-26 10/500 2024-05-28 14:42 by 江东闲人
[基金申请] 面上基金会评专家,有回避机制吗? +4 huang1991js 2024-05-27 4/200 2024-05-27 19:08 by 星火12
[硕博家园] 我是很理想化一人 +6 hahamyid 2024-05-26 6/300 2024-05-27 18:13 by 大飞鱼鱼鱼
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
信息提示
请填处理意见