24小时热门版块排行榜

返回列表

当前主题已经存档。

oldpig

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 40.6
帖子: 67
在线: 8.8小时
虫号: 125345
注册: 2005-12-06
性别: GG
专业: 理论和计算化学

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

一般得到的是基态的电荷密度，如要想要得到处于激发态是的总电子密度（如电子从HOMO->LUMO或LUMO+1跃迁后的总电荷密度），应该怎样设置输入或者如何从输出文件的提取。

回复此楼

» 猜你喜欢

1楼 2009-12-21 20:23:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhou2009

版主 (著名写手)

QC强帖: 15
应助: 337 (大学生)
贵宾: 4.053
金币: 42654.6
红花: 153
帖子: 1729
在线: 1274.7小时
虫号: 676800
注册: 2008-12-18
性别: GG
专业: 理论和计算化学
管辖: 量子化学

★ ★ ★ ★
小木虫(金币+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

激发态，在弛豫前，激发时构型来不及改变，激发态的构型就可以取基态的构型来计算。