24小时热门版块排行榜    

查看: 1554  |  回复: 9

byin

铜虫 (正式写手)

[交流] 【求助】请问有人算过铁电体的软模吗? 已有3人参与

计算了钛酸钡的声子散射谱,可是怎么也算不出软模。不知道怎么回事?还请高手们多多指教。



这是文献中的
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢回复 2011-03-23 08:42:09
我用VASP和pwscf算过可以的,你可以把输入文件贴出来看看

[ Last edited by zhang668 on 2011-3-22 at 23:39 ]
2楼2011-03-22 23:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)

引用回帖:
Originally posted by zhang668 at 2011-03-22 23:37:29:
我用VASP和pwscf算过可以的,你可以把输入文件贴出来看看

[ Last edited by zhang668 on 2011-3-22 at 23:39 ]

SYSTEM=BaTiO3_phonon
PREC=high
ISTART=0
ICHARG=2
ENCUT=800
EDIFF=0.1E-06
ISMEAR=-5
SIGMA=0.20
NSW=100
IBRION=5
ISIF=3
NFREE=2
POTIM=0.01
LWAVE=.FALSE.
LCHARG=.FALSE.
**********************************************
Automatic mesh
0
Monkhorst-pack
7 7 7
0 0 0
**********************************************
BaTiO3_phonon
4.04
1.000000000000000       0.000000000000000       0.000000000000000
0.000000000000000       1.000000000000000       0.000000000000000
0.000000000000000       0.000000000000000       1.000000000000000
3 1 1
Selective
Direct
0.0000000000000000   0.5000000000000000   0.5000000000000000        T        T        T
0.5000000000000000   0.0000000000000000   0.5000000000000000        T        T        T
0.5000000000000000   0.5000000000000000   0.0000000000000000        T        T        T
0.5000000000000000   0.5000000000000000   0.5000000000000000        T        T        T
0.0000000000000000   0.0000000000000000   0.0000000000000000        T        T        T
**********************************************


还请多多指点
3楼2011-03-23 08:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢提示 2011-03-23 14:01:37
这种有限差分法只能计算Gamma点的声子,不知道你怎么得来整个声子谱的。
4楼2011-03-23 09:38:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)

引用回帖:
Originally posted by zhang668 at 2011-03-23 09:38:33:
这种有限差分法只能计算Gamma点的声子,不知道你怎么得来整个声子谱的。

是的,只有GAMMA点。那个图是我用castep算的。

您不是提到vasp了吗,我用vasp也算过,同样也没有虚频。请版主帮我看看vasp的输入文件是否正确。

多谢多谢
5楼2011-03-23 09:48:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢提示 2011-03-23 14:01:51
虚频在OUTCAR里对应  f/i= 这样的形式。你在OUTCAR找下。
6楼2011-03-23 10:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)

引用回帖:
Originally posted by zhang668 at 2011-03-23 10:21:30:
虚频在OUTCAR里对应  f/i= 这样的形式。你在OUTCAR找下。

原来那就是虚频啊,那就有了,而且很多。
但是有个问题,如果我只考虑z方向的运动,那与Ti在一个x-y平面内的两个O应该振幅一致才是。可是我的计算结果中没有一个是这样的。

如果想要重复Cohen,1990,PRB上面的

需要如何设置vasp参数,请您多多指教。
7楼2011-03-23 11:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang668

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-24 08:15:46
有时候是计算精度问题,你把dx,dy,dz里很小的值看成0就行。你的ENCUT太大了,收敛精度也够了。POSCAR没必要设置成分子动力学,一般的就可以。
8楼2011-03-23 20:28:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)

引用回帖:
Originally posted by zhang668 at 2011-03-23 20:28:56:
有时候是计算精度问题,你把dx,dy,dz里很小的值看成0就行。你的ENCUT太大了,收敛精度也够了。POSCAR没必要设置成分子动力学,一般的就可以。

按照您的提示,我计算了一下。输入文件和结果:
***************************************
SYSTEM=BaTiO3_phonon
PREC=ACCURATE
ISTART=0
ICHARG=2
ENCUT=500
EDIFF=1E-07
ISMEAR=-5
SIGMA=0.20
NSW=100
IBRION=5
ISIF=3
NFREE=2
POTIM=0.01
LWAVE=.FALSE.
LCHARG=.FALSE.
***************************************
Automatic mesh
0
Monkhorst-pack
7 7 7
0 0 0
***************************************
BaTiO3_phonon
4.04
1.00       0.00       0.00
0.00       1.00       0.00
0.00       0.00       1.00
3 1 1
Selective
Direct
0.00   0.50   0.50        F        F        T
0.50   0.00   0.50        F        F        T
0.50   0.50   0.00        F        F        T
0.50   0.50   0.50        F        F        T
0.00   0.00   0.00        F        F        T
***************************************
***************************************
Eigenvectors and eigenvalues of the dynamical matrix
----------------------------------------------------


   1 f  =   16.339658 THz   102.665101 2PiTHz  545.032319 cm-1    67.575419 meV
             X         Y         Z           dx          dy          dz
      0.000000  2.020000  2.020000            0           0   -0.810971  
      2.020000  0.000000  2.020000            0           0   -0.255736  
      2.020000  2.020000  0.000000            0           0    0.412193  
      2.020000  2.020000  2.020000            0           0    0.326851  
      0.000000  0.000000  0.000000            0           0    0.013817  

   2 f  =   11.622032 THz    73.023380 2PiTHz  387.669245 cm-1    48.064878 meV
             X         Y         Z           dx          dy          dz
      0.000000  2.020000  2.020000            0           0    0.462840  
      2.020000  0.000000  2.020000            0           0   -0.823027  
      2.020000  2.020000  0.000000            0           0    0.306192  
      2.020000  2.020000  2.020000            0           0    0.116969  
      0.000000  0.000000  0.000000            0           0    0.031148  

   3 f  =    5.897317 THz    37.053934 2PiTHz  196.713306 cm-1    24.389350 meV
             X         Y         Z           dx          dy          dz
      0.000000  2.020000  2.020000            0           0   -0.272114  
      2.020000  0.000000  2.020000            0           0   -0.333826  
      2.020000  2.020000  0.000000            0           0   -0.354800  
      2.020000  2.020000  2.020000            0           0   -0.516373  
      0.000000  0.000000  0.000000            0           0    0.649608  

   4 f/i=    0.753578 THz     4.734869 2PiTHz   25.136649 cm-1     3.116548 meV
             X         Y         Z           dx          dy          dz
      0.000000  2.020000  2.020000            0           0    0.225961  
      2.020000  0.000000  2.020000            0           0    0.312941  
      2.020000  2.020000  0.000000            0           0    0.245035  
      2.020000  2.020000  2.020000            0           0    0.464423  
      0.000000  0.000000  0.000000            0           0    0.758471  

   5 f/i=   11.562640 THz    72.650207 2PiTHz  385.688129 cm-1    47.819251 meV
             X         Y         Z           dx          dy          dz
      0.000000  2.020000  2.020000            0           0   -0.054783  
      2.020000  0.000000  2.020000            0           0   -0.218754  
      2.020000  2.020000  0.000000            0           0   -0.741902  
      2.020000  2.020000  2.020000            0           0    0.630204  
      0.000000  0.000000  0.000000            0           0   -0.039623  
***********************************
第四个看做0吗?
最后一个是虚频,但是还是那个问题,最上面两个O的唯一不一样啊~~
怎么解决呢?

[ Last edited by byin on 2011-3-24 at 23:54 ]
9楼2011-03-24 23:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjttalento

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问下楼主这个问题解决了么…最近我也在计算这个,但是,Ba原子的dz的值却特别的小……

发自小木虫Android客户端
10楼2017-05-01 01:16:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 byin 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +7 wangmingqi 2026-08-28 8/400 2026-08-28 18:42 by zhuzg0628
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见