24小时热门版块排行榜    

查看: 567  |  回复: 7
当前主题已经存档。

ppshanshanqiu

铜虫 (小有名气)

[交流] 【求助】TD计算中峰强相差很大,怎么办

用TD-DFT方法计算了紫外光谱,直接拿算出来的振则强度对波长做了个简单的柱状图,结果表明,三个最大吸收峰位置相对实验值都蓝移了20nm,最大的问题是260nm以上的峰相对实验值处的峰而言,峰强增大了很多,这样的图是不是不合理的,但是又不知道问题出在哪里,怀疑是不是初始构型有问题。
by the way,我采用的计算方法是PCM,正庚烷体系,B3LYP 6-31 G+(d),实验体系是正己烷溶剂中,上传了图片(不知道有没有成功),求指点!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppshanshanqiu

铜虫 (小有名气)

我把整个计算思路列一下,
(1)基态,PCM,6-31G+d 优化
(2)TD ,PCM,同上基组,singlet,1-50
本来加上三重态一起算的,结果发现算出来的三重态都是振则强度为0,后来就直接全改成单重算TD了
2楼2008-11-02 11:23:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家


zzgyb(金币+1,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
是主族元素体系还是过度金属体系呢?开壳层体系还是闭壳层体系?
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2008-11-02 15:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppshanshanqiu

铜虫 (小有名气)

谢谢楼上的回帖,该分子中只有CHO和Br原子,是有机物,至于开壳和闭壳,我选的都是在单态下的计算,闭壳和开壳都有可能是单态吧,不知道这个在输入文件里应该怎么表示,比如我要做个闭壳体系,应该怎么表述呢?
4楼2008-11-02 15:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
TD-DFT计算的激发态就目前的水平为只能对上就很不错了,强度你不要管了,没有意义(对现在的计算水平)
蓝移了20nm也是很不错的结果了
5楼2008-11-02 20:26:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppshanshanqiu

铜虫 (小有名气)

谢谢楼上的回复,只是第三个(右数)峰强度相差的实在是有点离谱,而这个峰又尤其重要,担心这种结果将来投出去,会直接被拒掉,求人指点,或者给个合理的解释
6楼2008-11-02 21:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ppshanshanqiu

铜虫 (小有名气)

我看了很多JPCA上的关于TD-DFT计算的文献,人家的强度和峰位置都是符合的很好的,很迷茫,不知道我的这个结果到底出错在什么地方
7楼2008-11-02 21:50:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
6-31+g(d)对Br来说,没有考虑相对论效应。改用sdd基组可能会有改善。

如果仍不解决问题,还可能是所选的泛函不好。通常B3LYP是不错的选择,但也有很多例外的体系。TDDFT激发能受GGA/Hybrid-GGA的精确交换项HF影响很大。可以分别试试BLYP(no HF),BHandHLYP(0.5HF),加上B3LYP(0.2HF),可以看出第一个强峰的波长变化趋势,从而找到一个具有合适HF因子的泛函。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
8楼2008-11-02 22:21:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ppshanshanqiu 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-18 10:28 by yS4CkeWDhD6V
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 C79jjtjAKjEn 2026-09-15 4/200 2026-09-18 06:44 by manzzl
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-18 03:04 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +8 QUjhNVAcOSff 2026-09-13 9/450 2026-09-17 22:05 by cuPZTDXS3VOd
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 Aj1rhIDL5ixY 2026-09-14 6/300 2026-09-17 14:17 by ObVzqQQh1rrE
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 0pYnUiPDfdnk 2026-09-14 4/200 2026-09-17 14:15 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:38 by k06BKNblrGRH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:16 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 17:12 by i3KBRD19ERkG
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 QUjhNVAcOSff 2026-09-13 7/350 2026-09-16 14:36 by i3KBRD19ERkG
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
信息提示
请填处理意见