24小时热门版块排行榜    

查看: 1923  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 编辑是什么意思 15+3 s090604054 2024-06-15 3/150 2024-06-16 10:29 by bobvan
[有机交流] 车间生产,真空度很高,温度很高,但减压蒸馏速度很慢。 10+12 召唤鬼泣lL 2024-06-13 36/1800 2024-06-16 09:20 by ddc805
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中! +6 恐龙爸爸 2024-06-14 6/300 2024-06-16 08:41 by shl2112501
[考博] 34岁读博士晚吗 +26 emitdne 2024-06-13 26/1300 2024-06-16 07:16 by liyeqik
[基金申请] 博后面上今天有bug可以看到是否资助? +20 lyfbangong 2024-06-12 31/1550 2024-06-15 21:18 by since—2010
[基金申请] Nature 11日发文,中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[论文投稿] 审稿人含糊拒稿,还需要回复吗?如何回复? 20+3 BruceChum 2024-06-15 17/850 2024-06-15 20:19 by arthas_007
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[论文投稿] 投了一篇4区的SCI,审稿人一个拒稿,一个小修,编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[论文投稿] 审稿问题:为什么荧光激发波长和紫外吸收波长差的大? 10+4 sdawege 2024-06-14 8/400 2024-06-14 22:39 by 东北读书人
[基金申请] 博士后基金需要结题吗? +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[论文投稿] ACS Nano投稿后分配到副编辑手里12天了,能确定送审了吗? +5 潇洒怡惜 2024-06-12 10/500 2024-06-14 09:56 by 潇洒怡惜
[论文投稿] 最近写了一篇控制优化领域的文章,可以投哪里啊?有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[硕博家园] 机械研究生如何拿到年薪40+w +13 阿巴阿巴哦哦 2024-06-11 15/750 2024-06-13 15:40 by 113745685
[基金申请] 博后特助这周出结果吗?往年都是啥时候啊? +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
信息提示
请填处理意见