24小时热门版块排行榜    

查看: 2274  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

guohui_Hao

银虫 (小有名气)

[交流] 【求助】Dmol3 晶体计算 已有5人参与

Dmol3 对于晶体光谱计算, 强度计算需要difference-diople caculations  不太明白怎么操作,急求各位帮帮忙!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢理论指导 2010-12-21 08:16:35
红外强度的计算需要用导偶极矩导数矩阵(对N个原子在x,y,z三个不同方向求偶极矩的导数),再乘以不同的简正模式(来自频率计算),就是该振动模式的红外强度。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2010-12-20 23:52:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

ym23

铁杆木虫 (著名写手)


ellsaking(金币+1):鼓励回帖交流~ 2010-12-20 23:23:13
看看help中有没有解释
博观约取、厚积薄发---壁立千仞、无欲则刚
2楼2010-12-20 23:12:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohui_Hao

银虫 (小有名气)

不是很明白 能否再详细一点 计算的properity 里除了frequency 还需要什么
4楼2010-12-21 13:33:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhang668(金币+2):多谢指导! 2010-12-22 12:58:42
先计算偶极矩导数矩阵D,矩阵大小是3 x 3N。第一个3表示偶极矩及其导数有x,y,z三个分量;3N表示对N个原子分别在x,y,z方向做位移并求导数。如果是高对称体系,很多矩阵元为零。

Q是某个振动频率的简正模式,是一个大小为3N的数组。

M是大小为3N的数组,每个元素对应原子量的平方根。

有了D,Q,M以后,计算红外强度的Matlab命令如下:

c=Q./M;
f=31.22329;
d=D*c;
d=d*f;
IR=sum(d.*d)
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2010-12-22 12:14:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见