24小时热门版块排行榜    

查看: 1307  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

byin

铜虫 (正式写手)

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

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



这是文献中的
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 10 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 0703化学调剂 ,六级已过,有科研经历 +14 曦熙兮 2026-03-15 14/700 2026-03-21 19:12 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
信息提示
请填处理意见