24小时热门版块排行榜    

CyRhmU.jpeg
查看: 6498  |  回复: 48
【奖励】 本帖被评价22次,作者ilpisces增加金币 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ilpisces

新虫 (初入文坛)


[资源] 用Gnuplot画出分子轨道

理论化学专业的学生经常需要解析分子轨道(MO),而常用的量子化学程序(比如Gaussian,Gamess,ADF 等),通常配有相应的应用程序将分子轨道可视化。其实,借用绘图软件(比如 Gnuplot)和一些脚本代码,自己动手也可画出分子轨道。而这一过程对于理解量子化学一些基本概念,比如原子轨道(AO),基组,分子轨道等,不失为一个好的训练方法。本文中例在Gnuplot 4.6版本运行。

[ Last edited by ilpisces on 2017-5-25 at 10:19 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 用Gnuplot画出分子轨道.pdf
  • 2017-05-24 08:45:28, 1.38 M

» 收录本帖的淘帖专辑推荐

粒子宇宙生命 第一性原理杂 量子化学和计算化学 Photonic Crystal

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

ilpisces

新虫 (初入文坛)


引用回帖:
2楼: Originally posted by scfslyzkf at 2017-05-24 09:53:41
这是你写的论文么?

是我写的,但不是论文。
身边一个学生问我什么是基组和分子轨道等,我为他写了这些东西。
如果其他学生有这方面疑惑,可以参考一下。
3楼2017-05-24 10:02:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
4楼: Originally posted by scfslyzkf at 2017-05-24 10:13:53
最近在学习自洽场方法,也涉及到基函数和轨道。有几个问题想请教您。
请问氢原子有精确解,基组与氢原子的精确解是什么关系呢?

谢谢回复!
我很乐意回答您的问题,但在论坛我只能定性地解释,很难回答清楚,反而可能会误导您,对不起!
因为在我看来,计算化学中很多概念,如果脱离程序和数值计算,仅用文字是很难清晰描述的。
您的问题在我上传的资料里可以回答一部分,建议您先看一下能否解决。
5楼2017-05-24 10:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
4楼: Originally posted by scfslyzkf at 2017-05-24 10:13:53
最近在学习自洽场方法,也涉及到基函数和轨道。有几个问题想请教您。
请问氢原子有精确解,基组与氢原子的精确解是什么关系呢?

不过谢谢您的提问!
我会考虑将针对它的回答抽时间写在新的文件中。
6楼2017-05-24 11:06:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
4楼: Originally posted by scfslyzkf at 2017-05-24 10:13:53
最近在学习自洽场方法,也涉及到基函数和轨道。有几个问题想请教您。
请问氢原子有精确解,基组与氢原子的精确解是什么关系呢?

氢原子有精确解,是因为它只有一个电子。
用类氢原子模型,可以解出s, p, d等单电子轨道,表示成函数形式,是Slater型函数,zeta为参数。
但实际原子大多含有多个电子,所处的原子轨道与类氢原子模型解出的轨道不会完全一样,但依然是Slater型函数,有s, p, d之分,不同的是zeta应该根据具体原子被优化。
描述一个原子所有原子轨道的函数集合,以及zeta值,都是基组应包含的信息。
以上是针对于Slater型基组。

如果是Gaussian型基组,则完全是为了计算上的便利,用Gaussian型函数去拟合Slater型函数,参数由zeta变为alpha。
这样,描述一个原子所有原子轨道的Gaussian型函数集合,以及alpha值,将包含在Gaussian型基组内。
12楼2017-05-26 15:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
13楼: Originally posted by scfslyzkf at 2017-05-26 16:33:51
谢谢!slater型波函数以及在SCF中得到的轨道与类氢离子波函数只是有效核电荷数Z不同,可以这样理解吗?用变分法确定那个Z。昨天把Gauss型分子积分又看了一遍,用几个Gauss型函数拟合一个slater函数,在分子自洽场计 ...

---------------------------------------------------------
slater型波函数以及在SCF中得到的轨道与类氢离子波函数只是有效核电荷数Z不同,可以这样理解吗?
---------------------------------------------------------
请问,“SCF中得到的轨道”是什么意思?
SCF之后得到的轨道是分子轨道,而Slater型函数是用来描述原子轨道。

---------------------------------------------------------
用几个Gauss型函数拟合一个slater函数,在分子自洽场计算过程中,组合系数和gauss函数的轨道系数都要优化吗?还是在做分子自洽场计算前已经优化好了?
---------------------------------------------------------
SCF之前已经优化好了,不同的contraction coefficients(d)和exponents(alpha)正是不同基组的区别所在。
14楼2017-05-26 17:21:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
13楼: Originally posted by scfslyzkf at 2017-05-26 16:33:51
谢谢!slater型波函数以及在SCF中得到的轨道与类氢离子波函数只是有效核电荷数Z不同,可以这样理解吗?用变分法确定那个Z。昨天把Gauss型分子积分又看了一遍,用几个Gauss型函数拟合一个slater函数,在分子自洽场计 ...

不好意思我想问一下,我上传的这份文件您是否已经看了呢?
因为您的这些问题,我以为在文件里是可以找到答案的。
也许可能是有不好理解的地方,如果这样的话,您可以告诉我。
15楼2017-05-26 17:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
16楼: Originally posted by scfslyzkf at 2017-05-26 20:04:46
SCF的解不是单电子函数吗?单电子函数又称为分子轨道。我想问这个单电子函数是类氢离子波函数么?
...

不是。
类氢原子解出的波函数是单电子波函,是原子轨道。
SCF解出的也是单电子波函,但是分子轨道。
18楼2017-05-26 20:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
17楼: Originally posted by scfslyzkf at 2017-05-26 20:08:30
看了,您写的内容很丰富,我还要仔细的读一读。不同的人看同一本书体会也不一样。可能是我刚入门不久,很多东西理解不准确,还请老师多多指导。(毕竟这是量子化学,专业性很强的东西,不像小说一遍就能明白)
...

你可以把里面的脚本在gnuplot下运行一遍。
19楼2017-05-26 20:11:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
20楼: Originally posted by scfslyzkf at 2017-05-26 21:22:47
独立粒子近似下,原子自洽场和分子自洽场的解是类似的吧?只是原子中单粒子对应的平均势场是中心对称的,而分子中单电子对应的平均势场符合分子点群的对称性。分子轨道,分子的电子组态按其所属点群的不可约表 ...

我的文件里没有原子自洽场,hfscf是用来解分子轨道的。如果你对原子自洽场感兴趣,这个文件里没有讨论。

[ 发自手机版 http://muchong.com/3g ]
21楼2017-05-26 22:27:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


在文件里用Gaussian函数去拟合slater函数,只是提了一句用最小二乘法,使一个积分最大。如果在这里会令人产生疑惑的话,我会再修改一下,因为有的基组中收缩系数和轨道指数会用原子自洽场去求得。但在这个文件中,hfscf是没有讨论的,就连分子轨道系数也是直接给出的。如果要写hfscf的话,那又要增加很大篇幅了。

[ 发自手机版 http://muchong.com/3g ]
22楼2017-05-26 22:49:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
23楼: Originally posted by scfslyzkf at 2017-05-26 23:42:38
量化书上说,分子自洽场继承了原子自洽场,对原子自洽场讨论得到的结论,在分子自洽场中同样适用。这里面分子关于无扰态,谱项,谱项波函数与原子中相关概念(除了少许变动)几乎一样。
还有个问题,分子轨道是分子 ...

--------------------------------------------
这里面分子关于无扰态,谱项,谱项波函数与原子中相关概念(除了少许变动)几乎一样。
--------------------------------------------
我觉得这句话还是值得斟酌的,比如原子力场中有角动量守恒,分子力场未必。
但我不想在这里讨论它太多,单靠文字描述,文字越多越容易引起误解。

-------------------------------------------
点群不可约表示最多是三维,所以分子轨道最多三重简并,不存在更高的简并度,这样的判断对吗?
-------------------------------------------
不对。有不可约表示大于三维的点群。
26楼2017-05-27 09:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
24楼: Originally posted by scfslyzkf at 2017-05-27 00:23:02
您的文章说的是如何根据基组和分子轨道系数画出分子轨道图。我可以这么理解吗?
主要是我最近在学自洽场,而且分子自洽场这部分书上写的很简略,很难看懂;如果您详细讲解,并配上几个不同点群分子的例子,我觉得非 ...

-----------------------------------------------------------
您的文章说的是如何根据基组和分子轨道系数画出分子轨道图。我可以这么理解吗?
-----------------------------------------------------------
可以。所以在这个帖子中,我希望尽量只回答与这份文件有关的问题。

-------------------------------------------------------------------
主要是我最近在学自洽场,而且分子自洽场这部分书上写的很简略,很难看懂;如果您详细讲解,并配上几个不同点群分子的例子,我觉得非常有意义,可以知道来龙去脉,毕竟化学生经常把软件当黑箱。
-------------------------------------------------------------------
HFSCF本身并不一定要讨论点群。

--------------------------------------------------------------------
不过那样的话,估计工作量会很大。即便不可,也希望您在有空的时候能够给予指导
--------------------------------------------------------------------
我相信您手里有不少量子化学的教科书,如果单是靠公式和文字来解释概念的话,任何一本书都比我写的详细。
但依然有学生无法将理论与程序连接起来,我如果继续写下去,肯定不会在理论部分做更多解释,而是侧重于理论是如何变为程序的。
虽然这份文件里的内容还没有进入到量化计算(比如积分的求解,矩阵对角化等),只是外围的解析而已,但如果您对程序感兴趣,而且没有写代码的经历,我建议您先把文件里的例子重复一遍(比如您可以将重复完成后的.eps文件上传上来,然后再问其它问题,我觉得这样更合适一点,您觉得呢?)。

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

27楼2017-05-27 09:33:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilpisces

新虫 (初入文坛)


引用回帖:
29楼: Originally posted by scfslyzkf at 2017-05-27 11:31:05
是的,从理论到程序有很大的鸿沟,感谢您的贡献和答疑。
如果我想学习Gnuplot, 看懂您的代码,应该如何入手呢?
主要编程这方面我完全空白,但大致意思我能明白,就是知道了函数的解析式,然后取值,在坐标 ...

---------------------------------------------------
如果我想学习Gnuplot, 看懂您的代码,应该如何入手呢?
---------------------------------------------------
1. 先安装gnuplot。
2. 将文件中脚本运行一遍。

至于另外两个问题,在您解决第一个问题后再讨论比较好。
30楼2017-05-27 11:43:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ilpisces 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见