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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】虚频和软模是一回事么？
201/2返回列表12下一页
查看: 3496  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

identation

金虫 (著名写手)

[交流] 【求助】虚频和软模是一回事么？ 已有10人参与

一些位移型相变，比如BaTiO3的顺电-铁电相变，在相变点附近有软模，这在第一性原理计算中怎么体现呢？是声子谱中存在虚频么？虚频和软模是一回事么？

[ Last edited by identation on 2010-7-1 at 20:02 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-04-25 21:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)
★
identation(金币+5):谢谢参与
先用xcrysden找特殊K点,

再编译gk.f, 网站上那个F77程式有错误
生成inp.kpt

[ Last edited by lono75 on 2010-4-25 at 22:06 ]
赞一下(1人) 回复此楼
由天命,成人事！
2楼2010-04-25 21:47:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★
identation(金币+5):谢谢参与
identation(金币+40): 2010-04-25 23:36
这个一般来说都是自己写的, 很简单!
那个web site不是提供了程序.
它的输入文件为syml，输出文件为inp.kpt。(你会做能带吧,这个取点跟能带的取点宗旨是一样的.)
其中syml输入文件的格式如下：
8
15 15 15 15 15 15 15
G 0.0 0.0 0.0
K -0.33333333333 0.6666666666667 0.000000000
H -0.33333333333 0.6666666666667 0.500000000
A 0.0 0.0 0.5
G 0.0 0.0 0.0
M 0.0 0.5 0.0
L 0.0 0.5 0.5
A 0.0 0.0 0.5
第一行用来标记有多少个特殊k点，下面是这些特殊k点之间每个要分多少个k点，接着就是这些特殊k点的坐标。
把下面这个程序编译,执行即可.
c +---------------------------------------------------------
c For generating k-points along the high-symmetry lines in
c Brillouin zone and for calculate band-structures !
c +----------------------------------------------------------
C ---------'syml'---------
c 6 : nhighk
c 20 20 20 10 20 : ndiv(i)
c X 0.5 0.0 0.5 : labhk(1),phighk(1,1),........
c G 0.0 0.0 0.0
c L 0.5 0.5 0.5
c W 0.5 0.25 0.75
c K 0.375 0.375 0.75
c G 0.0 0.0 0.0
c direct & reciprocal lattice vectors over 'emin, emax' line
C -----------------------
c max k-points = 200
program gk
implicit real*8 (a-h,o-z)
character*2 labhk
dimension tkpt(200,3),pk(200,3),phighk(10,3)
dimension disk(200),dish(10),labhk(10)
dimension ndiv(10)
c
open(5,file='syml',status='old')
open(7,file='inp.kpt')
c
read(5,*) nhighk
read(5,*) (ndiv(i),i=1,nhighk-1)

do i=1,nhighk-1
ntkp=ntkp+ndiv(i)
enddo
ntotkpt=ntkp+1

if(nhighk>10)then
write(*,*)'Number of high-symmetry k points must < 10!'
STOP
endif
if(ntotkpt>200)then
write(*,*)'Total number of k points must <= 200!'
STOP
endif

do i=1, nhighk
read(5,*) labhk(i),(phighk(i,j),j=1,3)
enddo
write(*,*) (labhk(i),i=1,nhighk)
c
c----- generating k-points along high symmetric lines --------
c
c
pk(1,1)=phighk(1,1)
pk(1,2)=phighk(1,2)
pk(1,3)=phighk(1,3)
ii=1
do i = 2, nhighk
delx = (phighk(i,1) - phighk(i-1,1))/float(ndiv(i-1))
dely = (phighk(i,2) - phighk(i-1,2))/float(ndiv(i-1))
delz = (phighk(i,3) - phighk(i-1,3))/float(ndiv(i-1))
do j=1, ndiv(i-1)
ii = ii + 1
pk(ii,1) = pk(ii-1,1) + delx
pk(ii,2) = pk(ii-1,2) + dely
pk(ii,3) = pk(ii-1,3) + delz
enddo
enddo
c
10 format(A34)
weight=1.d0
do i=1,ntotkpt
write(7,200) pk(i,1),pk(i,2),pk(i,3),weight
enddo
200 format(3F10.6,F6.2)
stop
end
c----------------------- end ---------------------------

网站上说的很清楚!
如果你想把两个k点分成十份, 就把15换成10就可以了!

[ Last edited by gavinliu7390 on 2010-4-25 at 22:34 ]
赞一下(1人) 回复此楼
真理是一点点接近的！
3楼2010-04-25 22:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mudannanzi

木虫 (正式写手)
★
identation(金币+5):谢谢参与
我也是新手，学习中
赞一下(1人) 回复此楼
迷糊
4楼2010-04-25 22:25:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mxcrhr07jili

金虫 (正式写手)
★
identation(金币+5):谢谢参与
同参与，同学习
回复此楼
静默以待，坚韧以行
5楼2010-04-26 07:49:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianai3385

铜虫 (正式写手)
★
identation(金币+5):谢谢参与
新手，学习中
回复此楼
6楼2010-04-26 08:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqpie

木虫 (职业作家)
identation(金币+5): 2010-04-26 11:58
似乎也能在materials studio 的cell文件中得到/
赞一下 回复此楼
7楼2010-04-26 09:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by gavinliu7390 at 2010-04-25 22:22:42:
这个一般来说都是自己写的, 很简单!
那个web site不是提供了程序.
它的输入文件为syml，输出文件为inp.kpt。(你会做能带吧,这个取点跟能带的取点宗旨是一样的.)
其中syml输入文件的格式如下：
8
15 15 15 15 ...

您好，谢谢您帮忙，那您能提供一下编译好的exe文件么？fortran编译我不懂，而且没下到fortran的编译器。谢谢您
赞一下 回复此楼
8楼2010-04-26 21:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
这个问题貌似讨论了很多次了哦。
最简单的办法就是一楼说的，用xcrysden直接产生，很傻瓜，也很强大
赞一下 回复此楼
9楼2010-04-26 22:12:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)
identation(金币+10):谢谢，邮箱站内信发给你，嘿嘿 2010-04-27 12:26
引用回帖:
Originally posted by identation at 2010-04-25 21:20:34:
请问怎么可以获得K_POINTS的路径呢，例如，http://valenhou.blog.edu.cn/2006/133292.html
该网站中计算了能带结构，想问一下
K_POINTS
28
0.0 0.0 0.0 1.0
0.0 0.0 0.1 1.0
0.0 0.0 0.2 1.0
0 ...

我可以给你一个在win下的程序的，把邮箱给我吧，但是还是要结合特殊K点来的。
赞一下 回复此楼
timeflies..
10楼2010-04-27 11:11:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 材料专硕306英一数二 +4 z1z2z3879 2026-03-16 6/300 2026-03-16 19:38 by z1z2z3879
[考研] 本人考085602 化学工程 专硕 +12 不知道叫什么！ 2026-03-15 14/700 2026-03-16 16:45 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
信息提示
关闭
请填处理意见
关闭