24小时热门版块排行榜    

查看: 6047  |  回复: 20

sunboyccc

铁杆木虫 (正式写手)


liliangfang(金币+1, 专家考核): 多谢指教 2012-01-11 20:24:32
引用回帖:
10楼: Originally posted by QuenZ at 2012-01-11 17:33:41:
VASP也不能直接得出偶极矩的数值?
这么说,你是只取一个BiO4X4结构单位,代表一个单胞BiOX进行计算偶极?具体的计算,是根据密立根电荷和键长分别计算出各个取向的偶极,再进行偶极的矢量合成,得到一个BiO4X4 ...

这有什么复杂,电荷Q是标量,令每个原子坐标为r=xi+yj+zk矢量形式,即使每两个近邻原子都是一对电偶极子(根据电偶极子定义,r必须很小,电四和电八又可忽略),体系总的电偶极矩不过是Q(r-r')求和,求和个数有限哈,记住用矢量形式更简洁哈
11楼2012-01-11 19:33:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

QuenZ

金虫 (小有名气)

引用回帖:
11楼: Originally posted by sunboyccc at 2012-01-11 19:33:03:
这有什么复杂,电荷Q是标量,令每个原子坐标为r=xi+yj+zk矢量形式,即使每两个近邻原子都是一对电偶极子(根据电偶极子定义,r必须很小,电四和电八又可忽略),体系总的电偶极矩不过是Q(r-r')求和,求和个数有限 ...

现在已经确定CASTEP不能直接算偶极矩了,大概只能自己进行矢量合成计算了,不过真要在晶体结构中去这样算,想死的心都有!因为对于很多多元素的复式格子,元素的化学环境和键的取向还是很复杂的,尽管它有对称性!
谢谢啊!
希望多多交流
12楼2012-01-12 21:37:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunboyccc

铁杆木虫 (正式写手)


liliangfang(金币+1): 谢谢交流 2012-01-13 10:15:43
引用回帖:
12楼: Originally posted by QuenZ at 2012-01-12 21:37:40:
现在已经确定CASTEP不能直接算偶极矩了,大概只能自己进行矢量合成计算了,不过真要在晶体结构中去这样算,想死的心都有!因为对于很多多元素的复式格子,元素的化学环境和键的取向还是很复杂的,尽管它有对称性 ...

我看了你说的复杂的格子,并不复杂,另外科研有时花笨功夫得到的,往往是你独有的东西,而商业软件计算得到的,通过相同的设置是大多数人都可以得到的,对于你的结构,所谓的矢量合成只是一些原子矢量坐标的加法而已,不用把每一对电偶极子的方向都写出来,所以与键的具体方向无关,因为已经内含在原子的矢量坐标中啦,矢量加法的结果已经包含两点之间的矢量方向,只考虑离子键共价键就行啦,考虑布居,也不过是每个原子坐标加两次而已,几十个原子小case啦,只是有些繁琐,而且对于晶体而言只算一个有限大的胞就可以啦,说的不对的地方请包涵。
13楼2012-01-13 08:23:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunboyccc

铁杆木虫 (正式写手)

谢谢你的金币,要是真的就好啦,哈哈
14楼2012-01-13 08:29:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingxiushi

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by navyzhang at 2012-01-10 19:51:31
你的方向跟我的方向很接近
我做的一篇工作
你可以参考一下
别忘了引用哦
Phys. Chem. Chem. Phys. 2012, 14, 1286-1292

我发的两片文章都引用了阁下的这篇文章,不过关于偶极矩的计算方法还不是太清楚,能不能麻烦你在这里多说两句话,解释一下。不胜感激
15楼2013-01-12 23:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

QuenZ

金虫 (小有名气)

引用回帖:
15楼: Originally posted by xingxiushi at 2013-01-12 23:39:39
我发的两片文章都引用了阁下的这篇文章,不过关于偶极矩的计算方法还不是太清楚,能不能麻烦你在这里多说两句话,解释一下。不胜感激...

你可以在小木虫里联系张海军,我后来问了他,他的意思是要算出正、负电荷中心(算法同数学里面质心的算法),你可以再问问
希望多多交流
16楼2013-01-15 19:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

navyzhang

木虫 (正式写手)

引用回帖:
15楼: Originally posted by xingxiushi at 2013-01-12 23:39:39
我发的两片文章都引用了阁下的这篇文章,不过关于偶极矩的计算方法还不是太清楚,能不能麻烦你在这里多说两句话,解释一下。不胜感激...

我的文章里有些东西没讲明白?
其实是算了各个原子上的密立根布居,根据电荷量q和原子间的几何距离r
然后利用简单的电偶极矩的定义,p=qr
当然BiO4X4里面不是简单的一个正负电荷中心
你可以找等效正负电荷中心,或者通过矢量叠加求
17楼2013-01-16 11:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

QuenZ

金虫 (小有名气)

引用回帖:
17楼: Originally posted by navyzhang at 2013-01-16 11:12:02
我的文章里有些东西没讲明白?
其实是算了各个原子上的密立根布居,根据电荷量q和原子间的几何距离r
然后利用简单的电偶极矩的定义,p=qr
当然BiO4X4里面不是简单的一个正负电荷中心
你可以找等效正负电荷中心 ...

截止p=qr,都明白。但始终得不到文献值,困惑。再细问两个问题:
1、等效电荷中心,即把4个O等效成一个负电中心,4个X等效成另一负电中心,一个Bi就是一个正电中心,Right?
【具体简单计算法:将BiO4X4单元放在以4个X为顶点的盒子里,则4个X可以等效到其所在面的中心,总贡献为0.5(1/8×4);4个O由于在棱上,也可以等效到面心,总贡献为1(1/4×4);Bi在体心,贡献为1。这样就简化为正上方一个O,中心一个Bi,正下方0.5个X的模型。事实上,即便是在一个单胞里,也可以取到一个O一个Bi和0.5个X的等效模型。最后将这些等效因子乘到电荷上去】
这样简化后我大约可以算出我的偶极矩的值,但用文献提供的数据却得不到文献的值。
2、矢量叠加法。取Bi为原点,分别算出每个Bi-O及每个Bi-Cl的偶极,得到8个p,矢量方向都取正电指向负电的方向,然后对这8个p求和。应该也是要考虑每个原子的贡献量的。

      以上分析是否正确?
希望多多交流
18楼2013-01-16 15:42:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

QuenZ

金虫 (小有名气)

引用回帖:
18楼: Originally posted by QuenZ at 2013-01-16 15:42:45
截止p=qr,都明白。但始终得不到文献值,困惑。再细问两个问题:
1、等效电荷中心,即把4个O等效成一个负电中心,4个X等效成另一负电中心,一个Bi就是一个正电中心,Right?
【具体简单计算法:将BiO4X4单元放在 ...

我感觉我那样等效有问题,因为X的贡献变成0.5了。话说为啥不以[BiO]作为一层,等效成正电中心,[X]作为一层,等效为负电中心。这样更符合BiOX的晶体结构诶!
再说这个BiOX这么反对称,激发后的电子空穴到底怎么迁移呢?按照偶极矩的计算模型,岂不是有一种载流子要堆积在X-X层间了?
希望多多交流
19楼2013-01-16 16:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

QuenZ

金虫 (小有名气)

希望多多交流
20楼2013-01-25 09:33:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 QuenZ 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 03:34 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 03:30 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 03:27 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 4/200 2026-08-04 03:25 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 02:10 by Ocx3OKOzGU4q
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 8/400 2026-08-03 16:02 by yanghuping
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见