24小时热门版块排行榜    

查看: 1674  |  回复: 5

maomao7910

金虫 (小有名气)

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

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

» 猜你喜欢

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

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

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 ]
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
4楼2011-07-12 18:22:31
已阅   回复此楼   关注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的回帖

iamikaruk

木虫 (著名写手)

【答案】应助回帖


maomao7910(金币+20): 2011-07-13 13:28:32
zzy870720z(金币+1): 鼓励交流 2011-07-13 16:18:32
问题的关键只在于表达出来的是直角坐标吗?如果是,那么菱方和六角的表达是一致的

[ 发自手机版 http://muchong.com/3g ]
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
6楼2011-07-13 12:47:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maomao7910 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 yaCh4X3wd045 2026-09-14 4/200 2026-09-17 13:17 by ObVzqQQh1rrE
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-17 05:50 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 04:14 by k06BKNblrGRH
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 03:52 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 02:04 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-17 02:02 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-16 07:12 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-16 02:50 by DgNGHc3h5tPl
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见