计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7362|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] g-c3n4的自旋密度分布不均问题

[复制链接 Copy URL]

173

帖子

0

威望

1036

eV
积分
1209

Level 4 (黑子)

各位老师好,关于加一个正电荷后的g-c3n4自旋密度分布图我有一个问题不太明白,我的优化过程(b3lyp/6-31g* em=gd3bj scrf=(smd,solvent=water))是这样的,考虑边界问题,固定顶端和边界的共6个N原子,并补上H原子自行优化,优化成功后,把电荷数和自旋电子数改为 1 2,在继续优化。优化成功后读取fchk文件,画出以下自旋密度图。图中电荷主要只分布在里面三角形的一侧N原子上,这让我感到很困惑,这种对称结构不应该是三个边上的N原子上都有分布吗?

spindensity.png (49.31 KB, 下载次数 Times of downloads: 32)

spindensity.png

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

2#
发表于 Post on 2021-8-1 19:16:42 | 只看该作者 Only view this author
说明你的体系不是严格对称的。你量一下键长键角,就会发现不一样
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

173

帖子

0

威望

1036

eV
积分
1209

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-8-1 20:35:31 | 只看该作者 Only view this author
wzkchem5 发表于 2021-8-1 19:16
说明你的体系不是严格对称的。你量一下键长键角,就会发现不一样

多谢wzkchem5老师,我猜也是这样了,体系最后的结构是自行优化出来的,是不是证明我的优化有问题,需要加上空间群去优化吗

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

4#
发表于 Post on 2021-8-1 21:01:59 | 只看该作者 Only view this author
din5g 发表于 2021-8-1 13:35
多谢wzkchem5老师,我猜也是这样了,体系最后的结构是自行优化出来的,是不是证明我的优化有问题,需要加 ...

不对称不一定是错的,可能是因为Jahn-Teller效应导致对称性自发破缺的
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

4290

帖子

4

威望

9558

eV
积分
13928

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2021-8-2 01:16:06 | 只看该作者 Only view this author
din5g 发表于 2021-8-1 20:35
多谢wzkchem5老师,我猜也是这样了,体系最后的结构是自行优化出来的,是不是证明我的优化有问题,需要加 ...

你把这个体系截成了团簇用高斯算是吧?那哪里有地方给你考虑空间群。。。周期性都没了哪来空间群
自动做多参考态计算的程序MOKIT

255

帖子

1

威望

2682

eV
积分
2957

Level 5 (御坂)

6#
发表于 Post on 2021-8-2 12:38:10 | 只看该作者 Only view this author
如果做频率肯定会发现三个环共面的此结构有虚频。部分文献是没有固定原子,优化至没有虚频的结构。这样就有个疑问,是固定原子不考虑虚频,还是放开原子消除虚频,哪种会更合理一些呢
我唯一知道的就是我一无所知,但我是化学小迷弟

1万

帖子

0

威望

9896

eV
积分
22146

Level 6 (一方通行)

7#
发表于 Post on 2021-8-2 14:51:56 | 只看该作者 Only view this author
ChemG 发表于 2021-8-2 05:38
如果做频率肯定会发现三个环共面的此结构有虚频。部分文献是没有固定原子,优化至没有虚频的结构。这样就有 ...

理论上应该把被固定原子对Hessian的贡献投影掉,再做频率分析。不过我不知道高斯是不是这么做的,或者说是否支持这么做
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

173

帖子

0

威望

1036

eV
积分
1209

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-8-2 19:45:09 | 只看该作者 Only view this author
zjxitcc 发表于 2021-8-2 01:16
你把这个体系截成了团簇用高斯算是吧?那哪里有地方给你考虑空间群。。。周期性都没了哪来空间群

多谢zjxitcc老师,我的意思是在结构优化前,调整结构使其满足对称性再优化

33

帖子

0

威望

181

eV
积分
214

Level 3 能力者

9#
发表于 Post on 2021-8-2 21:19:57 | 只看该作者 Only view this author
请问你固定N原子后优化出来是平面吗?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-27 04:59 , Processed in 0.196570 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list