24小时热门版块排行榜    

查看: 678  |  回复: 4
当前主题已经存档。

oldpig

铁虫 (小有名气)

[交流] 【求助】在做TDDFT计算时如何得到激发态的总电荷密度

一般得到的是基态的电荷密度,如要想要得到处于激发态是的总电子密度(如电子从HOMO->LUMO或LUMO+1跃迁后的总电荷密度),应该怎样设置输入或者如何从输出文件的提取。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 12-22 13:35
想要得到处于激发态是的总电子密度(如电子从HOMO->LUMO或LUMO+1跃迁后的总电荷密度),设置输入时建议用关键词cube=(Orbitals,cards),指定你要算的轨道LUMO或LUMO+1的轨道序号。这样就得到了它们的MO的cube数据文件。将这cube数据文件输入到gsgrid程序,对这MO的cube数据进行平方运算,就得到跃迁后的LUMO或LUMO+1的电子密度。
当你的MO可以在一个平面上展示的时候,可以用GsGrid程序对它取一个最有代表性的截面,作ρ的等值线图。当我们考察的ρ,不能在一个平面上展示的时候,可以用GaussView程序对它作等值面图。

因为G03不能作一个MO的density,作了Orbitals又不能进入cubman作平方运算作成ρ,因为cubman缺少平方运算。所以需要进入GsGrid程序对它们作运算。

最近网友lihb734求助帖,上传了一段文献上作分子电子激发前后的Δρ(CDD)的例子。
文献详细地说明了这种Δρ的计算和作图的操作步骤,当然也包含你要的激发态是的电子密度,值得学习:
http://faculty.ycp.edu/~jforesma ... l/DiffGaussView.htm

激发态,在弛豫前,激发时构型来不及改变,激发态的构型就可以取基态的构型来计算。
2楼2009-12-22 08:49:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oldpig

铁虫 (小有名气)

非常感谢
3楼2009-12-23 15:47:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主


erylingjet(金币+1,VIP+0):感谢交流 12-24 15:22
引用回帖:
Originally posted by zhou2009 at 2009-12-22 08:49:
因为G03不能作一个MO的density,作了Orbitals又不能进入cubman作平方运算作成ρ,因为cubman缺少平方运算。所以需要进入GsGrid程序对它们作运算

周先生,从Gaussian03的B.01版开始,cubman中有square选项,可以做平方运算
4楼2009-12-23 17:13:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
谢谢yjcmwgk 的指正。

我用的版本还没有。

[ Last edited by zhou2009 on 2009-12-24 at 17:14 ]
5楼2009-12-23 17:36:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oldpig 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见