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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】虚频和软模是一回事么？
201/2返回列表12下一页
查看: 3479  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
关闭
请填处理意见
关闭