计算化学公社

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

[综合讨论] 关于Co-crystal的计算

[复制链接 Copy URL]

156

帖子

0

威望

3571

eV
积分
3727

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 superrice 于 2015-6-5 23:46 编辑

1. 请问有办法计算co-crystal的吸收光谱吗?现在组里用的Crystal14,不过手册里好像没有光谱的计算,请问其他软件有吗?
2. 这个co-crystal里有氢键存在,在结构优化的时候氢原子会从acceptor移动到donor,请问这样的话是用实验数据还是用优化过之后的数据呢?
3. 一般如果体系有氢键存在,会计算什么性质?体系里两种分子,怎么区分哪个是酸或者哪个是碱?有方法可以计算吗?
4. 比如valance band,能知道这个band是由分子的哪个轨道,比如HOMO, HOMO-1 之类组成的吗?

11

帖子

0

威望

133

eV
积分
144

Level 2 能力者

16#
发表于 Post on 2022-1-10 23:40:54 | 只看该作者 Only view this author
我计算了金刚石的,但是没有IR活性。事实上在3-7微米的时候,确实是有吸收,这个不知道能不能计算?我现在一直在摸索,没有计算出来?还有,请问您有crystal的安装包吗?想要个盗版,嘿嘿,谢谢您能回复我。

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

15#
发表于 Post on 2022-1-10 23:34:41 | 只看该作者 Only view this author
lhd 发表于 2022-1-10 21:52
用castep可以计算晶体的红外吸收吗?实验测得的金刚石晶体的红外透过谱,能用CASTEP计算吗?谢谢

红外吸收谱应该可以算,你可以看下性质计算有IR、Raman的计算选项,不过需要注意的是可能只支持使用模守恒赝势。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

11

帖子

0

威望

133

eV
积分
144

Level 2 能力者

14#
发表于 Post on 2022-1-10 21:53:17 | 只看该作者 Only view this author
楼主现在还在吗?这个crystal的安装是需要购买使用还是有别的渠道呢?谢谢

11

帖子

0

威望

133

eV
积分
144

Level 2 能力者

13#
发表于 Post on 2022-1-10 21:52:20 | 只看该作者 Only view this author
卡开发发 发表于 2015-6-16 20:02
手里的是7.0的,对于8.0的情况不明。
前面看了下8.0的help文档,上面还是这么写着的。也去查了一下CASTE ...

用castep可以计算晶体的红外吸收吗?实验测得的金刚石晶体的红外透过谱,能用CASTEP计算吗?谢谢

464

帖子

11

威望

3954

eV
积分
4638

Level 6 (一方通行)

12#
发表于 Post on 2021-9-21 22:39:08 | 只看该作者 Only view this author
CP2K也支持TDDFT计算,楼上的老师们都用过这个没呀?其优势和短板在哪里呢?
自由发挥,野蛮生长

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

11#
发表于 Post on 2015-6-17 10:03:06 | 只看该作者 Only view this author

可靠,不过需要设置比较大的nstate

418

帖子

4

威望

2912

eV
积分
3410

Level 5 (御坂)

10#
发表于 Post on 2015-6-16 21:25:21 | 只看该作者 Only view this author
jiewei 发表于 2015-6-16 20:00
层主有所不知,8.0 可用td-dft计算周期性的吸收光谱

方法可靠吗?

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

9#
发表于 Post on 2015-6-16 20:02:11 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2015-6-16 20:26 编辑
jiewei 发表于 2015-6-16 19:54
你用的是哪个版本?

手里的是7.0的,对于8.0的情况不明。
前面看了下8.0的help文档,上面还是这么写着的。也去查了一下CASTEP开发组的信息,对这个问题也没有明确说明。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

8#
发表于 Post on 2015-6-16 20:00:51 | 只看该作者 Only view this author
万里云 发表于 2015-6-16 14:27
现在TDDFT只适用于非周期性结构。周期性结构需要非常特殊的交换关联泛函,只有少数程序支持,并且指明这项 ...

层主有所不知,8.0 可用td-dft计算周期性的吸收光谱

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

7#
发表于 Post on 2015-6-16 19:54:41 | 只看该作者 Only view this author
卡开发发 发表于 2015-6-16 11:39
CASTEP没办法使用TDDFT做,手册给出计算限制:
Only a Gamma point calculation is allowed, which shou ...

你用的是哪个版本?

418

帖子

4

威望

2912

eV
积分
3410

Level 5 (御坂)

6#
发表于 Post on 2015-6-16 14:27:40 | 只看该作者 Only view this author
现在TDDFT只适用于非周期性结构。周期性结构需要非常特殊的交换关联泛函,只有少数程序支持,并且指明这项功能是"EXPERIMENTAL"。可以尝试exciting / elk这两个程序。不过这两个程序是全电子程序,估计计算量很大。

现在算周期性结构介电函数最常用的是DFT+RPA,并且视材料不同,可以用杂化泛函或者DFT+U。VASP和CASTEP都支持此种计算。

更为精确的计算方法是解Bethe-Salpeter方程,不过这个计算量惊人。

对于第四点:固体中波函数用n和k两个量子数标识,一条能带囊括了所有n为定制的波函数,因此“一条能带由哪些原子波函数贡献”的说法其实不准确。
VASP可以通过PROCAR得到此种信息,pwscf也可以。

评分 Rate

参与人数
Participants 2
eV +9 收起 理由
Reason
sobereva + 4
卡开发发 + 5 我很赞同

查看全部评分 View all ratings

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

5#
发表于 Post on 2015-6-16 11:39:32 | 只看该作者 Only view this author
jiewei 发表于 2015-6-16 10:15
co-crystal的吸收光谱吗?

用castep的TD-DFT方法可算。

CASTEP没办法使用TDDFT做,手册给出计算限制:
Only a Gamma point calculation is allowed, which should be sufficiently accurate for a "molecule in a box" geometry。

但是CASTEP可以通过计算Optical properties来获得介电常数换算到吸收谱,与VASP的过程没啥区别,只是这个换算过程程序能够直接实现而已。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

207

帖子

0

威望

2324

eV
积分
2531

Level 5 (御坂)

4#
发表于 Post on 2015-6-16 10:15:21 | 只看该作者 Only view this author
co-crystal的吸收光谱吗?

用castep的TD-DFT方法可算。

3754

帖子

3

威望

1万

eV
积分
19678

Level 6 (一方通行)

围观吃瓜群众

3#
发表于 Post on 2015-6-6 20:06:19 | 只看该作者 Only view this author
1、VASP可以计算,但需要计算得到介电常数的实部和虚部进行换算,网上有小程序可以进行换算,但使用并不方便,自己用Fortran写一下也并不复杂,公式黄昆先生的书上有,如果采取这样的方式搞不定我们可以再进行交流;
2、赞同楼上的,H一般没法测,一般看计算结果;
3、这个我不是很熟悉;
4、有些程序应该可以做分子轨道投影态密度,但哪些程序能够实现我并不清楚(印象中GPAW是可以做的)这样能够得到大致的分布,对于单个的能带如果写程序原则上应该也能实现,但我没做过相关的处理。

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
chenghao + 3 精品内容
sobereva + 3

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

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

GMT+8, 2025-8-16 19:43 , Processed in 0.157512 second(s), 22 queries , Gzip On.

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