identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】虛頻和軟模是一回事么？已有10人參與

一些位移型相變，比如BaTiO3的順電-鐵電相變，在相變點附近有軟模，這在第一性原理計算中怎么體現(xiàn)呢？是聲子譜中存在虛頻么？虛頻和軟模是一回事么？

[ Last edited by identation on 2010-7-1 at 20:02 ]

回復此樓

» 猜你喜歡

能源材料化學課題組招收碩士研究生8-10名已經(jīng)有6人回復
296求調(diào)劑已經(jīng)有5人回復
0703化學調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有8人回復
268求調(diào)劑已經(jīng)有8人回復
材料，紡織，生物（0856、0710），化學招生啦已經(jīng)有7人回復
301求調(diào)劑已經(jīng)有5人回復
281求調(diào)劑（0805）已經(jīng)有7人回復
材料與化工求調(diào)劑已經(jīng)有6人回復
299求調(diào)劑已經(jīng)有4人回復
277調(diào)劑已經(jīng)有6人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何選decoy set，做虛擬篩選？已經(jīng)有11人回復
【求助】請教大俠們：怎么將介電常數(shù)的實部和虛部在origin里面進行cole-cole 擬合？已經(jīng)有18人回復
【求助】dmol3 過渡態(tài)消虛頻！！已經(jīng)有4人回復
廣東省科技圖書館虛擬參考咨詢與原文傳遞平臺，讓你輕松找到想要的文獻已經(jīng)有4人回復
【求助】vasp優(yōu)化氣體分子，總是出現(xiàn)虛頻，該如何處理？已經(jīng)有5人回復
【轉(zhuǎn)帖】《PLoS綜合》：虛擬現(xiàn)實技術(shù)可使男人感到擁有女性身體已經(jīng)有3人回復
【分享】推薦100個好網(wǎng)站，從此上網(wǎng)不再虛度【已搜索無重復】已經(jīng)有228人回復

1樓 2010-04-25 21:20:34

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時
蟲號: 932458
注冊: 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
identation(金幣+5):謝謝參與

先用xcrysden找特殊K點,

再編譯gk.f, 網(wǎng)站上那個F77程式有錯誤
生成inp.kpt

[ Last edited by lono75 on 2010-4-25 at 22:06 ]

贊一下(1人)

回復此樓

由天命,成人事！

2樓2010-04-25 21:47:44

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gavinliu7390

木蟲 (著名寫手)

葉落鷹飛

1ST強帖: 3
應(yīng)助: 51 (初中生)
貴賓: 0.481
金幣: 1564.7
散金: 2975
紅花: 34
帖子: 1278
在線: 241.4小時
蟲號: 675997
注冊: 2008-12-17
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和

★
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點的坐標。
把下面這個程序編譯,執(zhí)行即可.
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 ---------------------------

網(wǎng)站上說的很清楚!
如果你想把兩個k點分成十份, 就把15換成10就可以了!

[ Last edited by gavinliu7390 on 2010-4-25 at 22:34 ]

贊一下(1人)

回復此樓

真理是一點點接近的！

3樓2010-04-25 22:22:42

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mudannanzi

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1959.1
帖子: 883
在線: 71.5小時
蟲號: 976416
注冊: 2010-03-20
性別: GG
專業(yè): 金屬材料的制備科學與跨學

★
identation(金幣+5):謝謝參與

我也是新手，學習中

贊一下(1人)

回復此樓

迷糊

4樓2010-04-25 22:25:28

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mxcrhr07jili

金蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 1532.9
散金: 77
紅花: 1
帖子: 543
在線: 510.9小時
蟲號: 868753
注冊: 2009-10-11
性別: GG
專業(yè): 固體力學

★
identation(金幣+5):謝謝參與

同參與，同學習

回復此樓

靜默以待，堅韌以行

5樓2010-04-26 07:49:07

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

tianai3385

銅蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 57.6
散金: 854
帖子: 523
在線: 30.5小時
蟲號: 926653
注冊: 2009-12-13
專業(yè): 光子與光電子器件

★
identation(金幣+5):謝謝參與

新手，學習中

回復此樓

6樓2010-04-26 08:14:35

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

qqpie

木蟲 (職業(yè)作家)

應(yīng)助: 0 (幼兒園)
金幣: 4657.9
帖子: 3245
在線: 227.6小時
蟲號: 243278
注冊: 2006-04-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

identation(金幣+5): 2010-04-26 11:58

似乎也能在materials studio 的cell文件中得到/

贊一下

回復此樓

7樓2010-04-26 09:11:14

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強帖: 16
應(yīng)助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

goldenfisher

金蟲 (著名寫手)

1ST強帖: 3
應(yīng)助: 55 (初中生)
貴賓: 0.153
金幣: 5779.1
散金: 2908
紅花: 21
帖子: 1552
在線: 351.5小時
蟲號: 117764
注冊: 2005-11-26
性別: GG
專業(yè): 原子和分子物理

這個問題貌似討論了很多次了哦。
最簡單的辦法就是一樓說的，用xcrysden直接產(chǎn)生，很傻瓜，也很強大

贊一下

回復此樓

9樓2010-04-26 22:12:12

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

enola

捐助貴賓 (正式寫手)

1ST強帖: 2
應(yīng)助: 13 (小學生)
金幣: 1485.9
散金: 1549
紅花: 14
帖子: 847
在線: 102.5小時
蟲號: 947564
注冊: 2010-01-24
專業(yè): 半導體光電子器件

identation(金幣+10):謝謝，郵箱站內(nèi)信發(fā)給你，嘿嘿 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
該網(wǎng)站中計算了能帶結(jié)構(gòu)，想問一下
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下的程序的，把郵箱給我吧，但是還是要結(jié)合特殊K點來的。

贊一下

回復此樓

timeflies..

10樓2010-04-27 11:11:24

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料專碩306英一數(shù)二 +7	z1z2z3879 2026-03-16	9/450	2026-03-17 17:31 by ccjequ
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +5	yangfz 2026-03-17	5/250	2026-03-17 17:01 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 308 085701 四六級已過求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 學碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[基金申請] 有必要更換申報口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 332求調(diào)劑 +3	zjy101327 2026-03-11	6/300	2026-03-13 22:48 by JourneyLucky
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 0856化學工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 0817化學工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao