计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6619|回复 Reply: 1

[综合讨论] 利用photoshop计算晶体表面孔洞面积

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89866

管理员

公社社长+计算化学玩家

发表于 Post on 2019-8-31 00:30:15 | 显示全部楼层 Show all |阅读模式 Reading model
利用photoshop计算晶体表面孔洞面积

文/Sobereva@北京科音   2019-Aug-31

今天思想家公社QQ群里有人发了一张图片,问文献中标注的MoS2的孔洞面积是怎么算的
1.png


其实用photoshop就可以实现计算。拿上图右下角那个体系为例,我们实际上要统计的就是图中孔洞对应的白色区域的面积。

把图片载入ps,用魔棒工具选择那块白色区域,再点击“窗口”-“直方图”,就看到了选区的像素数为1630。
2.png


然后按照下图,用矩形选择工具,从图中最上面一层的最左边的硫的中心位置拉到最右边的硫的中心位置,总共横跨12个S-S键。
3.png


然后按ctrl+Q把选区截取出来,再按ctrl+alt+I,可见横向像素是209。
4.png


笔者量了一下M$里自带的MoS2的结构中相邻的S之间距离,是3.166埃。因此12个S-S键总长度是3.166*12=37.992埃,因此上图中每个像素对应37.992/209=0.18178埃。因此1630像素面积就折合1630*0.18178*0.18178=53.86 埃^2。

还有一种计算方式是根据原子半径算。我们将S-S键长3.166的一半作为S原子球半径,图中孔洞相当于缺了7个硫,因此面积是7*pi*(3.166/2)^2=55.11 埃^2,和我们用photoshop那种方式算的基本吻合。

至于他人传到群里的那张图里标注的孔洞面积是37.345 埃^2,笔者不知怎么来的,我没有那篇文章,可能文中的数据不合理,反正本文的计算方式无疑是很合理的。值得提醒的是,用photoshop量面积的方式计算的时候,应当使用正交视角显示晶胞,否则由于透视视角的近大远小会导致算出来的面积没有意义。另外,在产生晶体图像的时候最好把像素分辨率设得高一些,这样通过上文的做法统计误差才会较小。

评分 Rate

参与人数
Participants 3
eV +14 收起 理由
Reason
灰天皮匠 + 5 牛!
rpestana94 + 4 <font style="vertical-align: inh
zsu007 + 5 不明觉厉

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口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!

877

帖子

36

威望

4800

eV
积分
6397

Level 6 (一方通行)

发表于 Post on 2019-9-1 22:25:56 | 显示全部楼层 Show all

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 00:38 , Processed in 0.197589 second(s), 25 queries .

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