计算化学公社

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

[Dalton] 关于dalton双光子吸收截面与双光子吸收波长的关系图要如何绘制

[复制链接 Copy URL]

68

帖子

0

威望

582

eV
积分
650

Level 4 (黑子)

各位老师好,我在使用dalton的时候,已经可以成功计算出分子结构的双光子最大吸收截面以及对应的吸收波长,但是我想请教一下各位老师要如何绘制像文献一样的双光子吸收截面与吸收波长之间的关系,即双光子吸收光谱图。我这边查阅官方的dalton手册里的双光子部分,也没有提到过双光子吸收光谱这一点,gaussian和multiwfn也不支持读取dalton的结果文件,同时双光子吸收光谱是双光子吸收截面与吸收波长之间的联系,与单光子吸收光谱不一样,而且从dalton最后的输入文件来看,它的Two-photon absorption summary 的第一列对应的就是最大双光子吸收截面以及它对应的吸收波长,再往后的话就是其他单重态的双光子吸收截面和吸收波长的相关信息,那这样的话dalton的双光子光谱要怎么绘制出来,从dalton输出文件的哪里获取光谱信息呢?

{7722708A-354E-4cb5-9AAF-B9101BFE575D}.png (93.96 KB, 下载次数 Times of downloads: 70)

文献里的TPA光谱图

文献里的TPA光谱图

{CC34BF9A-0440-4599-B1CC-D62775DA3D1C}.png (58.99 KB, 下载次数 Times of downloads: 62)

本人通过dalton计算得到的双光子结果

本人通过dalton计算得到的双光子结果

s0.out

343.31 KB, 下载次数 Times of downloads: 14

本人的dalton结果输出文件

22

帖子

0

威望

105

eV
积分
127

Level 2 能力者

14#
发表于 Post on 2024-12-11 20:57:13 | 只看该作者 Only view this author

老师,我根据您分享的内容,查看了muitiwfn手册的对应部分,里面有关于激发能、振子强度和摩尔吸收吸收系数的公式,里面我有两个问题,一个是反推的话,公式是加权求和并且高斯函数及其中的参数未知,第二个问题是我是通过利用激发能和双光子吸收截面在multiwfn中展宽得到光谱图,利用代入公式的是不是也是反推出激发能和双光子吸收截面?这个双光子吸收系数的反推还是不太清楚,能不能再辛苦老师指导一下

458

帖子

0

威望

1713

eV
积分
2171

Level 5 (御坂)

13#
发表于 Post on 2024-12-11 12:03:28 | 只看该作者 Only view this author
Yaoyl 发表于 2024-12-11 10:20
谢谢老师,老师我还想问一下我这边在multiwfn中利用高斯展宽得到了光谱图,因为高斯函数里面没摩尔吸收系 ...

参看multiwfn手册

multiwfn_mannul_p262.png (234.64 KB, 下载次数 Times of downloads: 55)

multiwfn_mannul_p262.png

22

帖子

0

威望

105

eV
积分
127

Level 2 能力者

12#
发表于 Post on 2024-12-11 10:20:00 | 只看该作者 Only view this author
cokie 发表于 2024-12-10 22:16
sob老师multiwfn绘制uv-vis光谱时得到的摩尔消光系数,是根据激发能+对应振子强度根据高斯展宽得到的。具 ...

谢谢老师,老师我还想问一下我这边在multiwfn中利用高斯展宽得到了光谱图,因为高斯函数里面没摩尔吸收系数,如何按照您说的利用现有的数据进行逆推啊?

高斯函数.png (36.76 KB, 下载次数 Times of downloads: 53)

高斯函数.png

spectrum_curve.txt

93.75 KB, 下载次数 Times of downloads: 3

spectrum_line.txt

2.26 KB, 下载次数 Times of downloads: 3

458

帖子

0

威望

1713

eV
积分
2171

Level 5 (御坂)

11#
发表于 Post on 2024-12-10 22:16:28 | 只看该作者 Only view this author
Yaoyl 发表于 2024-12-10 21:13
老师,您好,我想问一下怎么利用摩尔吸收系数转换成双光子吸收系数啊?这边有两张图片,以摩尔吸收系数为 ...

sob老师multiwfn绘制uv-vis光谱时得到的摩尔消光系数,是根据激发能+对应振子强度根据高斯展宽得到的。具体公式可以搜multiwfn手册里的关键词“Theory of CTS”。

你在用multiwfn绘制双光子吸收系数/吸收界面曲线时所用的.txt文件,实际上就是把绘制紫外时,所需要的激发能和对应的振子强度,换成了吸收系数和对应的吸收截面罢了。所以按uv-vis绘制出来光谱所得到的消光系数,实际上就是通过“吸收系数+对应吸收截面”进行高斯展宽得到的,所以根据公式逆推回去,自然就能把纵坐标推回相应的双光子吸收系数。

22

帖子

0

威望

105

eV
积分
127

Level 2 能力者

10#
发表于 Post on 2024-12-10 21:13:34 | 只看该作者 Only view this author
cokie 发表于 2024-7-10 16:03
半峰宽(FWHM)设窄一点就行。
我把你的log改了一下,定好了具体的半峰宽,直接加载,然后改一下横坐标 ...

老师,您好,我想问一下怎么利用摩尔吸收系数转换成双光子吸收系数啊?这边有两张图片,以摩尔吸收系数为纵坐标的是我根据dalton计算的out文件,提取数据写成txt格式,然后利用Mulitwfn处理出来的吸收光谱图,另一张是典型的双光子吸收系数光谱图,您看一下是我哪个地方不太对?下面附了相应的数据文件。

dislin.png (36.23 KB, 下载次数 Times of downloads: 46)

dislin.png

双光子.jpg (433.67 KB, 下载次数 Times of downloads: 59)

双光子.jpg

ta2.out

252.43 KB, 下载次数 Times of downloads: 4

ta2.txt

342 Bytes, 下载次数 Times of downloads: 2

22

帖子

0

威望

105

eV
积分
127

Level 2 能力者

9#
发表于 Post on 2024-12-10 11:20:13 | 只看该作者 Only view this author
您好,我想问一下您通过Multiwfn软件做出双光子光谱图后,它的纵坐标是摩尔吸收系数,怎么转换成双光子吸收系数啊?

68

帖子

0

威望

582

eV
积分
650

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-7-10 19:40:43 | 只看该作者 Only view this author
本帖最后由 zhouyulin 于 2024-7-10 19:48 编辑
cokie 发表于 2024-7-10 16:03
半峰宽(FWHM)设窄一点就行。
我把你的log改了一下,定好了具体的半峰宽,直接加载,然后改一下横坐标 ...

感谢老师的帮助,我这边按照您的指导尝试一下,但是遇到了一个问题:我在multiwfn上面调试好吸收光谱图后,通过2号选项“ Export X-Y data set of lines and curves to plain text file”,输出curve文件,但是我在用origin打开curve文件的时候,发现虽然吸收光谱图是一样的,但是它显示的部分好像是x轴与左边的y轴的函数图,并不是我想要的x轴的吸收波长与右侧y轴的双光子吸收截面的函数图,而且当我在origin把刻度和范围进行调整的时候,图像直接消失。我想请教一下有什么办法解决这个情况

{6167132D-7F4D-4fd9-A6DD-2474B6000B7C}.png (110.51 KB, 下载次数 Times of downloads: 71)

在multiwfn调整好的双光谱吸收图结果

在multiwfn调整好的双光谱吸收图结果

{D9134F50-8BBF-4cca-8176-261CA24BCC3A}.png (14.41 KB, 下载次数 Times of downloads: 66)

curve文件在origin中的显示状况

curve文件在origin中的显示状况

spectrum_curve.txt

96.68 KB, 下载次数 Times of downloads: 2

spectrum_line.txt

111 Bytes, 下载次数 Times of downloads: 2

458

帖子

0

威望

1713

eV
积分
2171

Level 5 (御坂)

7#
发表于 Post on 2024-7-10 16:03:46 | 只看该作者 Only view this author
本帖最后由 cokie 于 2024-7-10 16:05 编辑
zhouyulin 发表于 2024-7-10 15:05
老师,我昨天尝试了一下我得到的光谱,与文献数据还是差许多,是我自己的输入文件设置的有问题,还是图还 ...

半峰宽(FWHM)设窄一点就行。
我把你的log改了一下,定好了具体的半峰宽,直接加载,然后改一下横坐标范围就行了。另外其余9行数据既然没有数值,删了也可以。

EG.png (128.45 KB, 下载次数 Times of downloads: 76)

EG.png

test-new.log

47 Bytes, 下载次数 Times of downloads: 12

68

帖子

0

威望

582

eV
积分
650

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-7-10 15:05:42 | 只看该作者 Only view this author
本帖最后由 zhouyulin 于 2024-7-10 15:18 编辑
cokie 发表于 2024-7-9 09:42
你就用Multiwfn按紫外光谱的逻辑来绘制就行:

双光子吸收波长跟紫外吸收波长对应,都是第一列的内容( ...

老师,我昨天尝试了一下我得到的光谱,与文献数据还是差许多,是我自己的输入文件设置的有问题,还是图还需要在进行调整,我这边能用的数据绘制双光子吸收光谱的时候,需要的只有最重要的是那个最大的双光子吸收截面以及对应的吸收波长,也就是第一行的内容,其他的内容都不是需要的部分,想问一下是不是数据不够的原因

dislin.png (30.2 KB, 下载次数 Times of downloads: 80)

用multiwfn做出的图

用multiwfn做出的图

test.log

342 Bytes, 下载次数 Times of downloads: 9

multiwfn做光谱用的输入文件

DCM-OH.out

444.42 KB, 下载次数 Times of downloads: 9

文献里提供的双光子输出结果

68

帖子

0

威望

582

eV
积分
650

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-7-9 13:46:23 | 只看该作者 Only view this author
cokie 发表于 2024-7-9 09:42
你就用Multiwfn按紫外光谱的逻辑来绘制就行:

双光子吸收波长跟紫外吸收波长对应,都是第一列的内容( ...

我知道了,我再去试试,谢谢老师

458

帖子

0

威望

1713

eV
积分
2171

Level 5 (御坂)

4#
发表于 Post on 2024-7-9 09:42:43 | 只看该作者 Only view this author
本帖最后由 cokie 于 2024-7-9 09:47 编辑
zhouyulin 发表于 2024-7-8 23:58
但是老师我想几个问题问一下
(1)绘制光谱的输入文件对应的横纵坐标的怎么办,我这边的纵坐标和横坐标 ...

你就用Multiwfn按紫外光谱的逻辑来绘制就行:

双光子吸收波长跟紫外吸收波长对应,都是第一列的内容(填以eV为单位的能量,但eV不需要写);
纵坐标双光子吸收截面你就当振子强度填就行(把第二列填成双光子吸收截面的数值)
第三列的8.0是半峰宽的意思,根据你的实际需求设置(不过设置多少无所谓,你用Multiwfn载入之后,在选项8“ Input full width at half maximum (FWHM), current:   0.80000 eV”里还是可以改的)

按上述填好,用multiwfn按紫外波长读取,你就可以得到紫外的图,这时候横坐标是波长,纵坐标虽然你写的是双光子吸收截面,但是Multiwfn显示的Y轴是振子强度(数值是你自己填的双光子吸收截面的数值),这无所谓,你只需要选择2号选项“ Export X-Y data set of lines and curves to plain text file”,输出curve,自己在origin里作图,然后把纵坐标改成双光子吸收截面即可。

68

帖子

0

威望

582

eV
积分
650

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-7-8 23:58:28 | 只看该作者 Only view this author
本帖最后由 zhouyulin 于 2024-7-9 00:00 编辑
sobereva 发表于 2024-7-8 23:10
试图自行整理成Multiwfn绘制光谱通用的文本格式作为输入文件,在Multiwfn手册这一节说了

但是老师我想几个问题问一下
(1)绘制光谱的输入文件对应的横纵坐标的怎么办,我这边的纵坐标和横坐标分别是双光子吸收截面和吸收波长,正常的话吸收光谱对应的那部分是振子强度和吸收波长,这个是自己设置的吗,在哪里设置
(2)因为我这边绘制双光子吸收光谱的时候,需要的只有最重要的是那个最大的双光子吸收截面以及对应的吸收波长,也就是第一行的内容,其他的内容都不是需要的部分。这样的话,我的输入数据是只要输入第一行的双光子吸收截面以及吸收波长就可以了吗?还说还要输入别的内容
  (3)还有就这个multiwfn示例中这个8.0是什么意思,以及[color=rgba(0, 0, 0, 0.85)]UV-VIS的振子强度是无量纲的这句话是啥意思

{5D25B763-AA06-4db6-A161-E4614B27EA8A}.png (36.21 KB, 下载次数 Times of downloads: 70)

{5D25B763-AA06-4db6-A161-E4614B27EA8A}.png

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2024-7-8 23:10:52 | 只看该作者 Only view this author
试图自行整理成Multiwfn绘制光谱通用的文本格式作为输入文件,在Multiwfn手册这一节说了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

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

GMT+8, 2026-2-22 04:59 , Processed in 0.180299 second(s), 24 queries , Gzip On.

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