计算化学公社

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

[VASP] 关于电荷密度、差分电荷密度以及电荷转移的问题

[复制链接 Copy URL]

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

本帖最后由 茶味果冻 于 2014-10-24 10:45 编辑

今天不经意翻了一篇文献,文件中建立的模型如下图1模型
这里文献是采用VASP计算的分子动力学的relax,之后有计算了电子密度(我理解的是)结果如下图图2计算结果

这里我有几个问题这三个图计算的是整个体系的电子密度还是差分电荷密度?中间图和最右边的图中的电子密度计算一般怎么实现?比如中间图的电子密度计算时候是把上面原子固定不计算,直接计算下面的原子的电荷密度?左边图是固定下面原子不计算,直接计算上边原子的电荷密度?是这样实现的么??最后,从这三个图中怎么观察电子转移情况?这个模型是金团簇吸附半导体二氧化钛表面,想分析光致激发,表面电荷转移的情况。本人新手刚学第一性原理计算不久,求高人解答啊 谢谢了!


ja5001592.pdf

4.01 MB, 阅读权限: 10, 下载次数 Times of downloads: 320

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
ca0yan9 + 4 谢谢

查看全部评分 View all ratings

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2014-11-20 19:05:55 | 只看该作者 Only view this author
damao4361556 发表于 2014-11-20 13:11
请问楼主用pwscf吗, 能帮我 解答下 这个问题吗 http://emuch.net/bbs/viewthread.php?tid=8185305&fpage=1 ...

ibrav=0的话,celldm(1)这个是必须有的pwscf中默认单位是bohr,所以a需要进行单位换算的,Bohr=0.52917706 angstrom;你换算就行了,后边cell_parameter指定晶胞的基矢,MS中建立好模型就能获得,注意单位

21

帖子

0

威望

108

eV
积分
129

Level 2 能力者

17#
发表于 Post on 2014-11-20 13:11:38 | 只看该作者 Only view this author
本帖最后由 damao4361556 于 2014-11-20 13:14 编辑

请问楼主用pwscf吗, 能帮我 解答下 这个问题吗 http://emuch.net/bbs/viewthread.php?tid=8185305&fpage=1
但是 本人还有几点不理解,
采用 ibrav=0 时, celldm(1) 究竟怎么取值呢?这个值是怎么得来的。 设置完以后 cell_parameter的 基矢该如何设置呢?
采用ibrav=0时, celldm(1)也可以不取值吗? 不取值的话 cell_parameter 的 基矢又该怎么设置呢?


谢谢

30

帖子

0

威望

151

eV
积分
182

Level 3 能力者

16#
发表于 Post on 2014-10-28 23:54:47 | 只看该作者 Only view this author
jiewei 发表于 2014-10-28 11:04
今天抽空看了一下原文。下图所有的都是含时电子密度。 Fig4a 指激发态的电子密度,Fig4b指激发态经过时间弛 ...

原来电子密度还分这么多种,学习了,看来还要好好看看

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2014-10-28 18:53:09 | 只看该作者 Only view this author
jiewei 发表于 2014-10-28 11:04
今天抽空看了一下原文。下图所有的都是含时电子密度。 Fig4a 指激发态的电子密度,Fig4b指激发态经过时间弛 ...

好的,学习了。如果要实现这样的计算,VASP自带的功能可以实现嘛?不过我现在学的是pwscf,里面用到TDDFT的好像是计算吸收光谱的,感觉文献里做的计算好深啊,我找到作者Oleg V. Prezhdo做的好几遍文章,都是用的一样的理论,做的计算。

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

14#
发表于 Post on 2014-10-28 11:04:34 | 只看该作者 Only view this author
今天抽空看了一下原文。下图所有的都是含时电子密度。 Fig4a 指激发态的电子密度,Fig4b指激发态经过时间弛豫后的电子密度,Fig4c指激发态经过时间弛豫变成基态的电子密度(HOMO,也就是最高能带)。a图需要进行激发态的优化计算,然后得到电子密度,右边两图在优化结构上进行MD模拟。

122

帖子

0

威望

3397

eV
积分
3519

Level 5 (御坂)

13#
发表于 Post on 2014-10-24 12:17:48 | 只看该作者 Only view this author

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2014-10-24 10:49:10 | 只看该作者 Only view this author
jiewei 发表于 2014-10-23 20:47
另外,我怎么看都觉得你给的图是说的能带电子密度图吧? 要不你上传文献,我去仔细看看。

文献已经上传,附件

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2014-10-24 10:48:38 | 只看该作者 Only view this author
jiewei 发表于 2014-10-23 20:44
为什么要固定原子呢?  差分电荷密度 你可以去看看sob老师的 mutiwfn。
至于最后一个问题,表面你阅读文 ...

她写的我看了,但是看不懂。。我是学暖通的,突然就被老师抓来做这方面的东西了,基础太差,时间又紧现在基本都是照葫芦画瓢的做。。。

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2014-10-24 10:46:21 | 只看该作者 Only view this author
panger 发表于 2014-10-21 20:17
能分享下文献吗?

附件

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2014-10-24 10:46:13 | 只看该作者 Only view this author
aqhuangry 发表于 2014-10-21 20:54
能分享下文件吗?。谢谢!

附件

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

8#
发表于 Post on 2014-10-23 20:47:45 | 只看该作者 Only view this author
jiewei 发表于 2014-10-23 20:44
为什么要固定原子呢?  差分电荷密度 你可以去看看sob老师的 mutiwfn。
至于最后一个问题,表面你阅读文 ...

另外,我怎么看都觉得你给的图是说的能带电子密度图吧? 要不你上传文献,我去仔细看看。

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

7#
发表于 Post on 2014-10-23 20:44:23 | 只看该作者 Only view this author
茶味果冻 发表于 2014-10-22 10:08
还有就是,图二中两边的图电子密度计算是我说的那样子嘛,就是要分别固定原子,计算其它的电子密度??从 ...

为什么要固定原子呢?  差分电荷密度 你可以去看看sob老师的 mutiwfn。
至于最后一个问题,表面你阅读文献量不够。所以你应该多看看有关电荷转移的文献就懂了。 如果我给你讲解可能你就懒了,不利于你获取知识,这也是sob老师没有回答你问题的原因(我估计)

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2014-10-23 08:45:02 | 只看该作者 Only view this author
:):):)能帮助解答下问题嘛,谢谢了@sobereva

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2014-10-22 10:08:08 | 只看该作者 Only view this author
jiewei 发表于 2014-10-21 22:01
我仔细看了看,图2和图1的截图方向似乎不一样,图二是从图1的左边看过去的。
另外无论咋个看图二都是计算 ...

还有就是,图二中两边的图电子密度计算是我说的那样子嘛,就是要分别固定原子,计算其它的电子密度??从图上怎么观察电子转移情况?这个您能帮解答吗谢谢!:)

111

帖子

0

威望

459

eV
积分
570

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2014-10-22 10:01:59 | 只看该作者 Only view this author
jiewei 发表于 2014-10-21 22:01
我仔细看了看,图2和图1的截图方向似乎不一样,图二是从图1的左边看过去的。
另外无论咋个看图二都是计算 ...

好的,谢谢您的解答:),我再仔细看看

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

GMT+8, 2025-8-15 07:18 , Processed in 0.185501 second(s), 25 queries , Gzip On.

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