24小时热门版块排行榜    

查看: 1513  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见