版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

本帖产生 1 个 1ST强帖，点击这里进行查看

identation

金虫 (著名写手)

[交流] 【求助】声子与相变

这是PRL上的一篇文章，通过计算立方相BiMnO3的声子（已经预先优化），得到了虚频，然后作者写出了
The corresponding
eigenvectors are (Sr;Mn;Ojj;O?T ? e0:10; 0:31;0:52;0:56T for G-AFM and (0.03,0.41, 0:42,0:57) for FM, showing displacement of both the Sr and
Mn cations relative to a fairly rigid oxygen octahedral network.
和
The largest energy gain we
found was 6:2 meV=f:u:, for Rt4 ?001 rotation angle of 4.8, resulting in a I4=mcm structure
请问原子位移和稳定相的对称型（空间群）是怎么得到的？
有没有通过计算声子得到体系稳定相的资料供参考的？感激不尽！！

» 收录本帖的淘帖专辑推荐

第一原理资料汇编	科研技巧	第一性原理	量化,第一性,MD笔记
good resources	soft study	vasp

» 猜你喜欢

2025冷门绝学什么时候出结果已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

扩散与相变习题解答已经有3人回复
声子谱虚频求助已经有10人回复
一个相变会发生，而且发生的快慢不同，那么相变发生的本质是什么？已经有21人回复
第一性原理计算中钙钛矿结构相变的问题？已经有15人回复
特定分子振动的不可约表示已经有5人回复
求助一个热电材料一个晶格热导率与温度关系的问题已经有4人回复
声子是不是一种物质已经有26人回复
八面体畸变后的空间群已经有4人回复
高温相是否稳定？已经有6人回复
声子色散谱的分析已经有16人回复
非gamma点处的虚频意味着什么？已经有14人回复
！！！！测试Gama点的声子频率，怎么看测试结果；频率为0，还是两次之间差别很小..... 已经有3人回复
【求助】有关声子软化的处理已经有12人回复
【求助】-OH和苯环相连后对苯环电子云密度的影响已经有8人回复
【求助】什么是声子已经有14人回复
【求助成功】施加对称性的问题已经有13人回复
【求助】为什么优化结果总是有虚频已经有5人回复
【原创】DFT计算体会已经有20人回复
【求助】怎样由能带结构和态密度看计算的结构是不是可能存在超导现象已经有12人回复

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2010-12-03 21:40:56

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

回帖支持 ( 显示支持度最高的前 50 名 )

zhang668

木虫 (正式写手)

★ ★
zzy870720z(金币+2):谢谢指教 2010-12-10 12:27:32
identation(金币+15): 2010-12-10 15:45:29
identation(金币+5): 2010-12-10 17:34:29
bingmou(1ST强帖+1):详细解释，并且带有应用实例和理论，授予强贴 2010-12-10 22:05:30

有软模，不代表一定会发生相变，你可以计算下声子谱，会发现其它高对称K点也会有软模出现，一般数值最大的软模是最容易发生软化。比如你可以计算下钙钛矿的声子谱，发现除了Gamma点有软模，（对应铁电相变），R点或者M点也经常有软模，对应氧八面体的旋转。有时候氧八面体旋转模式软化后，会抑制铁电软模的软化，所以Gamma点有软模不一定发生铁电相变。
如果确定是某个软模软化造成的相变，根据软模可以确定空间群，比如Gamma点软模方向沿z方向，则从Pm-3m到P4mm，沿【111】方向，则到R3m。R点的氧八面体沿【111】软化时，则到R-3c。

赞一下(6人)

13楼2010-12-10 10:31:06

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

普通回帖

zhang668

木虫 (正式写手)

★ ★
identation(金币+8): 2010-12-04 13:50:50
sunyang1988(金币+2):谢谢交流 2010-12-04 14:59:25

原子位移就是计算声子得到的动力学矩阵元，稳定的对称性就是对立方相施加了软模的位移之后得到的结构

赞一下(1人)

2楼2010-12-04 13:49:16

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-04 13:49:16:
原子位移就是计算声子得到的动力学矩阵元，稳定的对称性就是对立方相施加了软模的位移之后得到的结构

谢谢哈
请问有这方面的资料么？
dielectric_hjzhao@126.com
能否发一份给我，或者你认为看看什么书能获得这方面的知识
谢谢哈

3楼2010-12-04 13:52:06

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

zhang668

木虫 (正式写手)

★ ★
identation(金币+15): 2010-12-04 14:12:38
sunyang1988(金币+2):谢谢交流 2010-12-04 14:59:40

介绍几篇文献吧：
PRB 60 836
PRB 42 6416
具体方法就是用软模方法来研究结构稳定性，认为高温立方相为原始结构，而基态相都是由各种结构畸变得到的，钙钛矿普遍存在的畸变有两种，一种是氧八面体的扭曲，属M或R点，不产生铁电性，另一种是A位或B位离子的移动，一般属Gamma点，方向可以沿4次轴或3次轴，甚至2次轴，一般产生铁电性。通过看计算声子得到的动力学矩阵元，可以看到A，B位位移大小，就知道在该模式中那个占主导！

赞一下(2人)

4楼2010-12-04 14:00:22

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-04 14:00:22:
介绍几篇文献吧：
PRB 60 836
PRB 42 6416
具体方法就是用软模方法来研究结构稳定性，认为高温立方相为原始结构，而基态相都是由各种结构畸变得到的，钙钛矿普遍存在的畸变有两种，一种是氧八面体的扭曲，属M或 ...

请问一下，黄昆的晶格动力学是不是就是介绍这些理论的教材？谢谢哈
想看看原理

5楼2010-12-04 14:16:51

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

zhang668

木虫 (正式写手)

★
identation(金币+15):谢谢哈；我看看 2010-12-04 14:33:47
sunyang1988(金币+1):谢谢交流 2010-12-04 14:59:52

黄昆的那本没介绍这方面，可以看下钟维烈的铁电体物理学

赞一下(1人)

6楼2010-12-04 14:20:50

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

sg18408926

至尊木虫 (著名写手)

应助: 0 (幼儿园)
金币: 10945.9
帖子: 1075
在线: 431.4小时
虫号: 727956

关注中

7楼2010-12-04 17:21:19

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

lzl8181

铁杆木虫 (著名写手)

1ST强帖: 3
应助: 10 (幼儿园)
金币: 6179.5
帖子: 1747
在线: 600.3小时
虫号: 777973

引用回帖:

Originally posted by zhang668 at 2010-12-04 14:20:50:
黄昆的那本没介绍这方面，可以看下钟维烈的铁电体物理学

搭车问一下这位兄弟，发现在声子计算中很多晶体即使用最高精度优化还是会出现软模，如何操作可以消除这些软模？

8楼2010-12-04 19:48:23

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

zhang668

木虫 (正式写手)

★
zzy870720z(金币+1):谢谢交流 2010-12-05 23:11:44

不知道计算得到的软模频率大吗？如果是很小的值基本就是计算问题了，如果达到几十或上百cm-1，要考虑是不是结构不稳定的软模了

赞一下(1人)

9楼2010-12-05 10:49:45

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-04 13:49:16:
原子位移就是计算声子得到的动力学矩阵元，稳定的对称性就是对立方相施加了软模的位移之后得到的结构

你好，请问动力学矩阵元.dyn文件里面哪些信息是原子位移，哪些又是晶格常数的变化呢？得到了动力学矩阵元，但是看不懂
以下的数据是什么含义呢？感激不尽！
  1  1
  0.12748015  0.00000000 0.00000000  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.12748015  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 0.12881256  0.00000000
  1  2
-0.02321693  0.00000000 0.00000000  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 -0.02321693  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.02301992  0.00000000
  1  3
-0.02333612  0.00000000 0.00065480  0.00000000 0.00000000  0.00000000
-0.00065480  0.00000000 -0.02333612  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.02606028  0.00000000
  1  4
-0.02088108  0.00000000 -0.00211841  0.00000000 0.00000000  0.00000000
  0.00211841  0.00000000 -0.02088108  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.01694545  0.00000000
  1  5
-0.00186596  0.00000000 0.00033727  0.00000000 0.00000000  0.00000000
-0.00007983  0.00000000 -0.00480048  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.01239771  0.00000000
  1  6
-0.00480048  0.00000000 0.00007983  0.00000000 0.00000000  0.00000000
-0.00033727  0.00000000 -0.00186596  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.01239771  0.00000000
  1  7
-0.00371702  0.00000000 0.00077140  0.00000000 0.00126688  0.00000000
-0.00086148  0.00000000 -0.00244921  0.00000000 -0.00126205  0.00000000
-0.00100066  0.00000000 0.00469964  0.00000000 -0.00204498  0.00000000
  1  8
-0.00371702  0.00000000 0.00077140  0.00000000 -0.00126688  0.00000000
-0.00086148  0.00000000 -0.00244921  0.00000000 0.00126205  0.00000000
  0.00100066  0.00000000 -0.00469964  0.00000000 -0.00204498  0.00000000
  1  9
-0.00244921  0.00000000 0.00086148  0.00000000 -0.00126205  0.00000000
-0.00077140  0.00000000 -0.00371702  0.00000000 -0.00126688  0.00000000
  0.00469964  0.00000000 0.00100066  0.00000000 -0.00204498  0.00000000
  1 10
-0.00244921  0.00000000 0.00086148  0.00000000 0.00126205  0.00000000
-0.00077140  0.00000000 -0.00371702  0.00000000 0.00126688  0.00000000
-0.00469964  0.00000000 -0.00100066  0.00000000 -0.00204498  0.00000000
  2  1
-0.02321693  0.00000000 0.00000000  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 -0.02321693  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.02301992  0.00000000
  2  2
  0.12748015  0.00000000 0.00000000  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.12748015  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 0.12881256  0.00000000
  2  3
-0.02333612  0.00000000 0.00065480  0.00000000 0.00000000  0.00000000
-0.00065480  0.00000000 -0.02333612  0.00000000 0.00000000  0.00000000
  0.00000000  0.00000000 0.00000000  0.00000000 -0.02606028  0.00000000
以下省略部分类似的内容
   Diagonalizing the dynamical matrix

   q = ( 0.000000000 0.000000000 0.000000000 )

**************************************************************************
   omega( 1) =    -0.977642 [THz] =    -32.610830 [cm-1]
(  0.000000  0.000000  0.000000  0.000000 -0.228436  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.228436  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.356796  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.344750  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.318843  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.318843  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.333974  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.333974  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.333974  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.333974  0.000000 )
   omega( 2) =    2.087454 [THz] =    69.630431 [cm-1]
(  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 )
(  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 )
(  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 )
(  0.000000  0.000000  0.000000  0.000000  0.000000  0.000000 )
(  0.000000  0.000000  0.000000  0.000000  0.005084  0.000000 )
(  0.000000  0.000000  0.000000  0.000000 -0.005084  0.000000 )
( -0.349920  0.000000  0.357131  0.000000  0.000000  0.000000 )
(  0.349920  0.000000 -0.357131  0.000000  0.000000  0.000000 )
( -0.357131  0.000000 -0.349920  0.000000  0.000000  0.000000 )
(  0.357131  0.000000  0.349920  0.000000  0.000000  0.000000 )
   omega( 3) =    2.330161 [THz] =    77.726308 [cm-1]
(  0.091196  0.000000 -0.148215  0.000000  0.000000  0.000000 )
(  0.091196  0.000000 -0.148215  0.000000  0.000000  0.000000 )
(  0.012317  0.000000 -0.417696  0.000000  0.000000  0.000000 )
(  0.092137  0.000000 -0.236277  0.000000  0.000000  0.000000 )
(  0.231703  0.000000  0.219166  0.000000  0.000000  0.000000 )
(  0.231703  0.000000  0.219166  0.000000  0.000000  0.000000 )
(  0.282183  0.000000  0.253460  0.000000  0.000000  0.000000 )
(  0.282183  0.000000  0.253460  0.000000  0.000000  0.000000 )
(  0.226328  0.000000  0.231156  0.000000  0.000000  0.000000 )
(  0.226328  0.000000  0.231156  0.000000  0.000000  0.000000 )
请问如何根据这些数据确定原子位移呢？同时比如说77.726308 [cm-1]下的振动模式是什么样的应该如何获得呢？比如说M-O的伸缩震动等等，如何标定呢？
谢谢

10楼2010-12-08 16:43:35

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

zhang668

木虫 (正式写手)

★ ★
sunyang1988(金币+2):谢谢回帖交流 2010-12-09 15:24:42
identation(金币+12): 2010-12-09 17:52:12
identation(金币+18): 2010-12-09 17:56:22

前面是动力矩阵，下面的才是动力学矩阵元，对应的就是原子位移。对于Gamma点来说，1，3，5列对应的就是每个原子的位移，2，4，6列总是0.知道原子位移画出就是振动模式。最好首先根据对称性分析看振动模式属于那个不可约表示，兼并度多少

赞一下(1人)

11楼2010-12-09 15:07:04

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-09 15:07:04:
前面是动力矩阵，下面的才是动力学矩阵元，对应的就是原子位移。对于Gamma点来说，1，3，5列对应的就是每个原子的位移，2，4，6列总是0.知道原子位移画出就是振动模式。最好首先根据对称性分析看振动模式属于那个 ...

那如果计算得到的结果中，几何优化精度已经很高了，仍然在gamma点出现虚频，那么可以断定会发生铁电相变吗？也就是说gamma点有虚频，而且这个虚频比较大，结果不稳定，能看出相变之后的空间群吗？

12楼2010-12-09 17:57:54

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

ustbwgs

木虫 (正式写手)

应助: 14 (小学生)
金币: 1449.6
帖子: 720
在线: 165.5小时
虫号: 1141993

★
identation(金币+1):谢谢参与

顶起来

14楼2010-12-10 12:44:11

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-10 10:31:06:
有软模，不代表一定会发生相变，你可以计算下声子谱，会发现其它高对称K点也会有软模出现，一般数值最大的软模是最容易发生软化。比如你可以计算下钙钛矿的声子谱，发现除了Gamma点有软模，（对应铁电相变），R点 ...

请问有没有系统介绍这些知识的教材呢？或者是教程什么的？
铁电体物理和相关的论文我看了一下，好像都是概述性地介绍，想看看系统的知识加强一下基础
谢谢！

16楼2010-12-10 16:09:30

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

zoujin129

禁虫 (知名作家)

本帖内容被屏蔽

17楼2010-12-10 16:43:40

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-09 15:07:04:
前面是动力矩阵，下面的才是动力学矩阵元，对应的就是原子位移。对于Gamma点来说，1，3，5列对应的就是每个原子的位移，2，4，6列总是0.知道原子位移画出就是振动模式。最好首先根据对称性分析看振动模式属于那个 ...

谢谢热心帮忙哈，再问一下哈，我用两种输入文件计算得到的Gamma点的声子相差很大，这是怎么回事呢？
输入文件分别如下：
phonon
&inputph
tr2_ph=1.0d-12,
prefix='PHONON',
fildvscf='PHONONdv',
amass(1)=*,
amass(2)=*,
amass(3)=*,
amass(4)=*,
outdir='./',
fildyn='PHONON.dyn',
elph=.false.,
trans=.true.,
ldisp=.true.,
nq1=1, nq2=1, nq3=1
/
另一种是
phonon
&inputph
tr2_ph=1.0d-14,
prefix='PHONON',
amass(1)=*,
amass(2)=*,
amass(3)=*,
amass(4)=*,
outdir='./',
fildyn='PHONON.dyn',
/
0.0 0.0 0.0

第一种方式得到的频率为
   omega( 1) =    -0.977642 [THz] =    -32.610830 [cm-1]
   omega( 2) =    2.087454 [THz] =    69.630431 [cm-1]
   omega( 3) =    2.330161 [THz] =    77.726308 [cm-1]
   omega( 4) =    2.330161 [THz] =    77.726308 [cm-1]
   omega( 5) =
   omega( 6) =
   omega( 7) =
   omega( 8) =
   omega( 9) =
   omega(10) =
   omega(11) =    4.517777 [THz] =    150.697833 [cm-1]
   omega(12) =    4.522795 [THz] =    150.865203 [cm-1]
   omega(13) =    4.715043 [THz] =    157.277958 [cm-1]
   omega(14) =    4.914551 [THz] =    163.932874 [cm-1]
   omega(15) =    4.914551 [THz] =    163.932874 [cm-1]
   omega(16) =    8.239041 [THz] =    274.826646 [cm-1]
   omega(17) =    8.239041 [THz] =    274.826646 [cm-1]
   omega(18) =    8.308145 [THz] =    277.131714 [cm-1]
   omega(19) =    11.747218 [THz] =    391.847617 [cm-1]
   omega(20) =    11.747218 [THz] =    391.847617 [cm-1]
   omega(21) =
   omega(22) =
   omega(23) =
   omega(24) =    13.043958 [THz] =    435.102505 [cm-1]
   omega(25) =    14.181002 [THz] =    473.030450 [cm-1]
   omega(26) =    14.240101 [THz] =    475.001786 [cm-1]
   omega(27) =    15.114911 [THz] =    504.182510 [cm-1]
   omega(28) =    15.589022 [THz] =    519.997271 [cm-1]
   omega(29) =    15.589022 [THz] =    519.997271 [cm-1]
   omega(30) =    23.166952 [THz] =    772.771485 [cm-1]
第二种方式得到的频率则是
   omega( 1) =    1.207097 [THz] =    40.264684 [cm-1]
   omega( 2) =    1.516525 [THz] =    50.586174 [cm-1]
   omega( 3) =    1.516525 [THz] =    50.586174 [cm-1]
   omega( 4) =    1.564170 [THz] =    52.175448 [cm-1]
   omega( 5) =
   omega( 6) =
   omega( 7) =
   omega( 8) =
   omega( 9) =
   omega(10) =
   omega(11) =    4.476623 [THz] =    149.325078 [cm-1]
   omega(12) =    4.506522 [THz] =    150.322389 [cm-1]
   omega(13) =    5.146817 [THz] =    171.680467 [cm-1]
   omega(14) =    5.338251 [THz] =    178.066077 [cm-1]
   omega(15) =    5.338251 [THz] =    178.066077 [cm-1]
   omega(16) =    8.249629 [THz] =    275.179818 [cm-1]
   omega(17) =    8.249629 [THz] =    275.179818 [cm-1]
   omega(18) =    8.269644 [THz] =    275.847472 [cm-1]
   omega(19) =    11.976885 [THz] =    399.508535 [cm-1]
   omega(20) =    12.010839 [THz] =    400.641144 [cm-1]
   omega(21) =
   omega(22) =
   omega(23) =
   omega(24) =    12.599531 [THz] =    420.277895 [cm-1]
   omega(25) =    14.205947 [THz] =    473.862535 [cm-1]
   omega(26) =    14.255378 [THz] =    475.511383 [cm-1]
   omega(27) =    15.127825 [THz] =    504.613293 [cm-1]
   omega(28) =    15.129552 [THz] =    504.670872 [cm-1]
   omega(29) =    15.129552 [THz] =    504.670872 [cm-1]
   omega(30) =    23.123039 [THz] =    771.306685 [cm-1]
相差很大，请问这是怎么回事呢？Gamma点的声子是属于上面的呢还是下面这种方式得到的呢？
   Mode symmetry, C_4h (4/m)  point group:

   omega(  1 -  1) =       40.3  [cm-1] --> ?
   omega(  1 -  1) =       40.3  [cm-1] --> 0E_g          R
   omega(  1 -  1) =       40.3  [cm-1] --> ?
   omega(  1 -  1) =       40.3  [cm-1] --> 0E_g*          R
   omega(  1 -  1) =       40.3  [cm-1] --> A_u          I
   omega(  1 -  1) =       40.3  [cm-1] --> ?
   omega(  1 -  1) =       40.3  [cm-1] --> 0E_u          I
   omega(  1 -  1) =       40.3  [cm-1] --> ?
   omega(  1 -  1) =       40.3  [cm-1] --> 0E_u*          I

   omega(  9 -  9) =       134.8  [cm-1] --> ?
   omega(  9 -  9) =       134.8  [cm-1] --> 0E_u*          I
   omega( 10 - 11) =       149.3  [cm-1] --> E_g          R
   omega( 10 - 11) =       149.3  [cm-1] --> E_g*          R
   omega( 12 - 12) =       150.3  [cm-1] --> B_g          R
   omega( 12 - 12) =       150.3  [cm-1] --> ?
   omega( 12 - 12) =       150.3  [cm-1] --> 0E_g          R
   omega( 12 - 12) =       150.3  [cm-1] --> ?
   omega( 12 - 12) =       150.3  [cm-1] --> 0E_g*          R
   omega( 12 - 12) =       150.3  [cm-1] --> ?
   omega( 12 - 12) =       150.3  [cm-1] --> 0E_u          I
   omega( 12 - 12) =       150.3  [cm-1] --> ?
   omega( 12 - 12) =       150.3  [cm-1] --> 0E_u*          I
   omega( 13 - 13) =       171.7  [cm-1] --> ?
   omega( 13 - 13) =       171.7  [cm-1] --> 0E_g          R
   omega( 13 - 13) =       171.7  [cm-1] --> ?
   omega( 13 - 13) =       171.7  [cm-1] --> 0E_g*          R
   omega( 13 - 13) =       171.7  [cm-1] --> A_u          I
   omega( 13 - 13) =       171.7  [cm-1] --> ?
   omega( 13 - 13) =       171.7  [cm-1] --> 0E_u          I
   o
   omega( 18 - 18) =       275.8  [cm-1] --> ?
   omega( 18 - 18) =       275.8  [cm-1] --> 0E_g*          R
   omega( 18 - 18) =       275.8  [cm-1] --> B_u
   omega( 18 - 18) =       275.8  [cm-1] --> ?
   omega( 18 - 18) =       275.8  [cm-1] --> 0E_u          I
   omega( 18 - 18) =       275.8  [cm-1] --> ?
   omega( 18 - 18) =       275.8  [cm-1] --> 0E_u*          I
   om
   omega( 19 - 19) =       399.5  [cm-1] --> 0E_u          I
   omega( 19 - 19) =       399.5  [cm-1] --> ?
   omega( 19 - 19) =       399.5  [cm-1] --> 0E_u*          I
   omega( 20 - 21) =       400.6  [cm-1] --> E_u          I
   omega( 20 - 21) =       400.6  [cm-1] --> E_u*          I

   omega( 23 - 24) =       420.3  [cm-1] --> E_g          R
   omega( 23 - 24) =       420.3  [cm-1] --> E_g*          R
   omega( 25 - 25) =       473.9  [cm-1] --> A_g          R
   omega( 25 - 25) =       473.9  [cm-1] --> ?
   omega( 25 - 25) =       473.9  [cm-1] --> 0E_g          R
   omega( 25 - 25) =       473.9  [cm-1] --> ?
   omega( 25 - 25) =       473.9  [cm-1] --> 0E_g*          R
   omega( 25 - 25) =       473.9  [cm-1] --> ?
   omega( 25 - 25) =       473.9  [cm-1] --> 0E_u          I
   omega( 25 - 25) =       473.9  [cm-1] --> ?
   omega( 25 - 25) =       473.9  [cm-1] --> 0E_u*          I
   omega( 26 - 26) =       475.5  [cm-1] --> B_g          R
   omega( 26 - 26) =       475.5  [cm-1] --> ?
   omega( 26 - 26) =       475.5  [cm-1] --> 0E_g          R
   omega( 26 - 26) =       475.5  [cm-1] --> ?

   omega( 27 - 27) =       504.6  [cm-1] --> ?
   omega( 27 - 27) =       504.6  [cm-1] --> 0E_u          I
   omega( 27 - 27) =       504.6  [cm-1] --> ?
   omega( 27 - 27) =       504.6  [cm-1] --> 0E_u*          I
   omega( 28 - 29) =       504.7  [cm-1] --> E_u          I
   omega( 28 - 29) =       504.7  [cm-1] --> E_u*          I
   omega( 30 - 30) =       771.3  [cm-1] --> A_g          R
   omega( 30 - 30
   omega( 30 - 30) =       771.3  [cm-1] --> ?
   omega( 30 - 30) =       771.3  [cm-1] --> 0E_u*          I
这列数据是什么意思呢？“？”处代表的是什么呢？
都很困惑，敬请帮忙，感激不尽！
是否有一些介绍这方面的书，可否推荐一下呵呵
感激不尽！

[ Last edited by identation on 2010-12-10 at 17:43 ]

18楼2010-12-10 17:40:17

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

zhang668

木虫 (正式写手)

★
gzqdyouxia(金币+1):谢谢提示 2010-12-10 19:22:28

引用回帖:

Originally posted by identation at 2010-12-10 16:09:30:

请问有没有系统介绍这些知识的教材呢？或者是教程什么的？
铁电体物理和相关的论文我看了一下，好像都是概述性地介绍，想看看系统的知识加强一下基础
谢谢！

最好看第一性原理计算铁电性的相关文献，都有具体的方法的，Ferroelectrics杂志有很多这方面的文章。关于铁电体的教材，我只看过钟伟烈的铁电体物理学，觉得还是在于看文献体会与积累吧。

赞一下(1人)

19楼2010-12-10 19:05:37

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

zhang668

木虫 (正式写手)

★
gzqdyouxia(金币+1):谢谢指教 2010-12-10 19:22:51

引用回帖:

Originally posted by identation at 2010-12-10 17:40:17:

谢谢热心帮忙哈，再问一下哈，我用两种输入文件计算得到的Gamma点的声子相差很大，这是怎么回事呢？
输入文件分别如下：
phonon
&inputph
tr2_ph=1.0d-12,
prefix='PHONON',
fildvscf='PHONONdv', ...

你的两种计算首先精度要求就不一样，第一种更多是在计算声子谱时候用到的，这个时候很少只取1个点的，平常计算单个K点最好用第二种。
两个都可以看。
omega( 1 - 1) = 40.3 [cm-1] --> 0E_g R
首先不知道你为什么一个模式书输出这么多遍。Eg中的E代表所属的不可约表示是二阶简并的，即会有两个振动模式频率一样。A，B是1阶，T是3阶。下表g、u分别代表中心反演对称和反对称等。对称的对应R即Raman激活，反对称u对应I，红外激活。这些你可以看下群论书，或者群论在物理中的应用

[ Last edited by zhang668 on 2010-12-10 at 19:27 ]

赞一下(1人)

20楼2010-12-10 19:17:24

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

zhang668

木虫 (正式写手)

★ ★
gzqdyouxia(金币+2):谢谢指教 2010-12-10 19:23:19

贴一个我算的立方晶系的Gamma点声子：

**************************************************************************
   omega( 1) =    -4.948777 [THz] = -165.074535 [cm-1]
   omega( 2) =    -4.948777 [THz] = -165.074535 [cm-1]
   omega( 3) =    -4.948777 [THz] = -165.074535 [cm-1]
   omega( 4) =    -0.218896 [THz] =    -7.301632 [cm-1]
   omega( 5) =    -0.218896 [THz] =    -7.301632 [cm-1]
   omega( 6) =    -0.218896 [THz] =    -7.301632 [cm-1]
   omega( 7) =    2.392257 [THz] =    79.797631 [cm-1]
   omega( 8) =    2.392257 [THz] =    79.797631 [cm-1]
   omega( 9) =    2.392257 [THz] =    79.797631 [cm-1]
   omega(10) =    6.340424 [THz] =    211.495176 [cm-1]
   omega(11) =    6.340424 [THz] =    211.495176 [cm-1]
   omega(12) =    6.340424 [THz] =    211.495176 [cm-1]
   omega(13) =    18.275622 [THz] =    609.613191 [cm-1]
   omega(14) =    18.275622 [THz] =    609.613191 [cm-1]
   omega(15) =    18.275622 [THz] =    609.613191 [cm-1]
**************************************************************************

   Mode symmetry, O_h (m-3m)  point group:

   omega(  1 -  3) =    -165.1  [cm-1] --> T_1u G_15  G_4- I
   omega(  4 -  6) =       -7.3  [cm-1] --> T_1u G_15  G_4- I
   omega(  7 -  9) =       79.8  [cm-1] --> T_2u G_25  G_5-
   omega( 10 - 12) =       211.5  [cm-1] --> T_1u G_15  G_4- I
   omega( 13 - 15) =       609.6  [cm-1] --> T_1u G_15  G_4- I
由于计算的体系是强关联体系，因此算出来是金属体系，没有LO-TO分裂，所以振动模式都是三度简并的。

赞一下(1人)

21楼2010-12-10 19:18:31

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

identation

金虫 (著名写手)

引用回帖:

Originally posted by zhang668 at 2010-12-10 19:17:24:

你的两种计算首先精度要求就不一样，第一种更多是在计算声子谱时候用到的，这个时候很少只取1个点的，平常计算单个K点最好用第二种。
两个都可以看。
omega( 1 - 1) = 40.3 [cm-1] --&g ...

对呀，为什么会输出那么多遍呢？这是个有自旋的体系，很困惑
然后还带有问号，是用PWSCF计算的

23楼2010-12-10 19:51:41

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

zhang668

木虫 (正式写手)

★
sunyang1988(金币+1):谢谢交流 2010-12-11 09:15:28

考虑换下并行参数试下，我觉得你出现的问号跟有些还没输出的模式，甚至还有o、 om这样的可能是程序输出结果错误或者什么造成的

赞一下(1人)

24楼2010-12-10 23:33:28

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

0710811132

荣誉版主 (文学泰斗)

应助: 0 (幼儿园)
贵宾: 1.944
金币: 20230
帖子: 80732
在线: 2802.2小时
虫号: 1047838

★
identation(金币+1):谢谢参与

不知道

25楼2011-01-10 14:24:42

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

我不是陈冠希

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 145.1
帖子: 111
在线: 4.7小时
虫号: 1183131

★
identation(金币+1):谢谢参与

不太清楚噢
本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2386813&fpage=1

赞一下(1人)

26楼2011-01-10 22:14:12

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

乐石||

木虫 (正式写手)

应助: 26 (小学生)
金币: 2416.4
帖子: 977
在线: 807.5小时
虫号: 1474843

★
identation(金币+1): 谢谢参与

引用回帖:

742807楼: Originally posted by zhang668 at 2010-12-09 15:07:04
前面是动力矩阵，下面的才是动力学矩阵元，对应的就是原子位移。对于Gamma点来说，1，3，5列对应的就是每个原子的位移，2，4，6列总是0.知道原子位移画出就是振动模式。最好首先根据对称性分析看振动模式属于那个不 ...

那么对于其他点呢？
除G点外，其他的点第偶数列是不为0的，其含义是什么？

27楼2012-05-28 11:40:13

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

迷途大书童

木虫 (正式写手)

应助: 9 (幼儿园)
金币: 1564
帖子: 974
在线: 282.1小时
虫号: 2393482

★
identation(金币+1): 谢谢参与

声子谱图，学过一段时间，还是不会分析！

28楼2013-12-09 09:16:17

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

简单回复

2010-12-10 12:56 回复

identation(金币+1):谢谢参与
zzy870720z:请不要再悬赏贴中纯表回复，相信你也不希望在求助时别人这么回复你，谢谢合作，下次就扣回所得了 2010-12-10 15:53:30

chdh88061822楼

2010-12-10 19:41 回复

identation(金币+1):谢谢参与
zzy870720z:请不要在悬赏贴中纯表回复，相信你也不希望在求助时得到这样的回复，谢谢合作 2010-12-10 21:35:54

奋力拼搏29楼

2015-01-28 16:16 回复

identation(金币+1): 谢谢参与

相关版块跳转我要订阅楼主 identation 的主题更新