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

登录注册

应《网络安全法》要求，自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

» 猜你喜欢

写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有3人回复
需要合成515-64-0，50g，能接单的留言已经有3人回复
最近几年招的学生写论文不引自己组发的文章已经有10人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复
临港实验室与上科大联培博士招生1名已经有8人回复
26申博自荐已经有7人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有4人回复
带资进组求博导收留已经有9人回复
求助大佬们，伤口沾上了乙腈已经有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的回帖

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的回帖

查看全部 29 个回答

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的回帖

查看全部 29 个回答

24小时热门版块排行榜

identation

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

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

» 猜你喜欢

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

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

identation

zhang668

identation

zhang668

» 抢金币啦！回帖就可以得到:

查看全部散金贴