应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp用菱方结构得到的BORN有效电荷矩阵和光学谱的方向问题
51/1返回列表
查看: 1718  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

maomao7910

金虫 (小有名气)

[求助] vasp用菱方结构得到的BORN有效电荷矩阵和光学谱的方向问题

各位，我要用vasp算一个R3c的结构，可以有六角胞表示，也可以用菱方胞来表示。我们一般计算时习惯用菱方结构来计算，这样原子数少，计算较快。不过，在用dfpt法得到波恩有效电荷3介矩阵和光学的实部和虚部时，是往那个方向投影的呢？是菱方的a，b，c轴方向，还是六角胞中的x，y，z方向呢？如果是前者，那么计算这些性质的时候还得转化成六角胞再算一次呢，原子数增加3倍，计算量不小。如果是后者，我急需有这方面经验的朋友给我个确切的答案。多谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-07-12 13:47:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao7910

金虫 (小有名气)
ion    1
    1     6.22415    -0.62964     0.00088
    2     0.63312     6.23359     0.00232
    3     0.00000     0.00000     4.75161
ion    2
    1     2.56974    -0.04776    -0.00010
    2     0.04966     2.57012    -0.00190
    3     0.00000     0.00000     2.35818
从我算的born电荷结果来看。在12和21分量上有小小不对称，不过这两个量相对较小，影响应该不大吧。我的问题是：比如离子1的这个矩阵我是用菱方结构计算得到的，和对应的六角的是一样的么？
  我根据MS里面六角和菱方的基矢矩阵仔细体会了一下。他的六角相如果保持c沿z轴，b在yz平面会得到一个相应的菱方的基矢。如果调整六角相在直角系的投影方向，会得到另一套基矢。所以我目前得到的结论是，如果你用其中一套菱方的基矢来算波恩电荷和光学介电性质。得到的是和这个菱方的基矢相对应的六角相的结果一致，不知道我的想法是否正确？请多指点，真的想弄明白这个问题。
另：iamikaruk请用回应求助，我可以把金币加给你。
赞一下 回复此楼
5楼2011-07-13 11:18:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

iamikaruk

木虫 (著名写手)
vasp用的是笛卡尔坐标吧？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
2楼2011-07-12 17:00:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao7910

金虫 (小有名气)
是对，vasp是把菱方的a，b，c基矢分解到了笛卡尔坐标系内。那是否意味着求出的结果也是沿笛卡尔坐标系分解的呢？
赞一下 回复此楼
3楼2011-07-12 18:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-07-13 09:33:33
玻恩有效电荷和介电常数张量应该满足菱方的对称性，你可以用这个验证一下。不过我记得VASP的有效电荷和介电常数数值上不是很严格的对称性

[ 发自手机版 http://muchong.com/3g ]
赞一下(1人) 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
4楼2011-07-12 18:22:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +6 Isak0320 2026-10-12 9/450 2026-10-12 21:13 by Isak0320
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么？ +3 好运来～9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭