24小时热门版块排行榜    

查看: 2756  |  回复: 15

guohui_Hao

银虫 (小有名气)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
sunyang1988(金币+2): 谢谢交流 2011-04-02 00:34:42
引用回帖:
Originally posted by yasistar at 2011-04-01 15:05:59:
你说的,先计算偶极距导数矩阵。

偶极距本来就是个 数阵,你怎么求导数呢?

一忙写错了。改了以后是这样

d(Dx)/dx1 d(Dy)/dy1 d(Dz)/dz1
d(Dx)/dx2 d(Dy)/dy2 d(Dz)/dz2
...
d(Dx)/dxn d(Dy)/dyn d(Dz)/dzn

1,2,...,n是原子的编号x,y,或z方向。Dx,Dy,Dz是偶极矩矢量在直角坐标系下的三个分量

[ Last edited by beefly on 2011-4-2 at 11:48 ]
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
14楼2011-04-01 22:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ym23

铁杆木虫 (著名写手)


ellsaking(金币+1):鼓励回帖交流~ 2010-12-20 23:23:13
看看help中有没有解释
博观约取、厚积薄发---壁立千仞、无欲则刚
2楼2010-12-20 23:12:52
已阅   回复此楼   关注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的回帖

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的回帖

guohui_Hao

银虫 (小有名气)

zzy870720z:如果问题解决,请帮忙吧标题【求助】修改为【求助成功】,谢谢合作 2010-12-22 14:32:29
非常感谢!! 帮我解决了很大的难题!!
6楼2010-12-22 13:16:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohui_Hao

银虫 (小有名气)

你好,可否Q聊下?
7楼2010-12-23 11:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohui_Hao

银虫 (小有名气)

254281974
8楼2010-12-23 11:06:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohui_Hao

银虫 (小有名气)

beefly 可否Q聊下
9楼2010-12-23 12:40:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guohui_Hao

银虫 (小有名气)

偶极矩大小也是3 x 3N吗  另外M的含义不是很明白? 麻烦你再解析下
10楼2010-12-23 13:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 guohui_Hao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +8 Ldrop2023 2026-08-13 8/400 2026-08-16 11:28 by jiege890620
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见