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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】虚频和软模是一回事么？
201/2返回列表12下一页
查看: 3610  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +3 汪！？！ 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 269求调剂 +6 啊啊我我 2026-04-07 6/300 2026-04-08 20:04 by 我减肥1
[考研] 土木水利专硕276分求调剂 +6 我想上学！！6 2026-04-05 9/450 2026-04-08 17:45 by 宋小宝HQ
[考研] 材料科学与工程320求调剂，080500 +12 黄瓜味薯片 2026-04-06 12/600 2026-04-08 16:26 by luoyongfeng
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +5 爱读书的小鳄鱼 2026-04-02 5/250 2026-04-08 14:55 by CAS小花工
[考研] 085404，285分求调剂 +4 薇薇考研 2026-04-07 4/200 2026-04-08 14:32 by zhq0425
[考研] 266求调剂 +18 阳阳哇塞 2026-04-07 18/900 2026-04-08 12:59 by grayjzr
[考博] 申博 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 一志愿郑州大学材料与化工085600，求调剂 +34 吃的不少 2026-04-02 34/1700 2026-04-07 20:01 by lrll?l
[考研] 作栽330调剂 +3 我要上好学 2026-04-02 4/200 2026-04-07 19:54 by biomichael
[考研] 0703调剂 +16 拾玖壹 2026-04-04 18/900 2026-04-07 12:49 by flydream1314
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸！ 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 26考研调剂0710 0860 +9 补补不补 2026-04-03 14/700 2026-04-04 23:32 by 果冻大王
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
信息提示
关闭
请填处理意见
关闭