计算化学公社

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

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

[复制链接 Copy URL]

40

帖子

0

威望

221

eV
积分
261

Level 3 能力者

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

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

文献里的TPA光谱图

文献里的TPA光谱图

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

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

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

s0.out

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

本人的dalton结果输出文件

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

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

40

帖子

0

威望

221

eV
积分
261

Level 3 能力者

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: 23)

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

338

帖子

0

威望

1107

eV
积分
1445

Level 4 (黑子)

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里作图,然后把纵坐标改成双光子吸收截面即可。

40

帖子

0

威望

221

eV
积分
261

Level 3 能力者

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

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

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

40

帖子

0

威望

221

eV
积分
261

Level 3 能力者

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: 20)

用multiwfn做出的图

用multiwfn做出的图

test.log

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

multiwfn做光谱用的输入文件

DCM-OH.out

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

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

338

帖子

0

威望

1107

eV
积分
1445

Level 4 (黑子)

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: 21)

EG.png

test-new.log

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

40

帖子

0

威望

221

eV
积分
261

Level 3 能力者

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: 22)

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

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

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

curve文件在origin中的显示状况

curve文件在origin中的显示状况

spectrum_curve.txt

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

spectrum_line.txt

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

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 05:21 , Processed in 0.175100 second(s), 24 queries , Gzip On.

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