 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3522 | 回復(fù): 19 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
identation金蟲 (著名寫手)
|
[交流]
【求助】虛頻和軟模是一回事么？ 已有10人參與
|
||
|
一些位移型相變，比如BaTiO3的順電-鐵電相變，在相變點(diǎn)附近有軟模，這在第一性原理計(jì)算中怎么體現(xiàn)呢？是聲子譜中存在虛頻么？虛頻和軟模是一回事么？ [ Last edited by identation on 2010-7-1 at 20:02 ] |
回復(fù)此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 如何選decoy set，做虛擬篩選？
已經(jīng)有11人回復(fù)
- 【求助】請(qǐng)教大俠們：怎么將介電常數(shù)的實(shí)部和虛部在origin里面進(jìn)行cole-cole 擬合？
已經(jīng)有18人回復(fù)
- 【求助】dmol3 過(guò)渡態(tài)消虛頻！！
已經(jīng)有4人回復(fù)
- 廣東省科技圖書館虛擬參考咨詢與原文傳遞平臺(tái)，讓你輕松找到想要的文獻(xiàn)
已經(jīng)有4人回復(fù)
- 【求助】vasp優(yōu)化氣體分子，總是出現(xiàn)虛頻，該如何處理？
已經(jīng)有5人回復(fù)
- 【轉(zhuǎn)帖】《PLoS綜合》：虛擬現(xiàn)實(shí)技術(shù)可使男人感到擁有女性身體
已經(jīng)有3人回復(fù)
- 【分享】推薦100個(gè)好網(wǎng)站，從此上網(wǎng)不再虛度【已搜索無(wú)重復(fù)】
已經(jīng)有228人回復(fù)
goldenfisher
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 55 (初中生)
- 貴賓: 0.153
- 金幣: 5779.1
- 散金: 2908
- 紅花: 21
- 帖子: 1552
- 在線: 351.5小時(shí)
- 蟲號(hào): 117764
- 注冊(cè): 2005-11-26
- 性別: GG
- 專業(yè): 原子和分子物理
9樓2010-04-26 22:12:12
lono75
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6963.1
- 散金: 100
- 紅花: 2
- 帖子: 2433
- 在線: 536.4小時(shí)
- 蟲號(hào): 932458
- 注冊(cè): 2009-12-23
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-04-25 21:47:44
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★
identation(金幣+5):謝謝參與
identation(金幣+40): 2010-04-25 23:36
identation(金幣+5):謝謝參與
identation(金幣+40): 2010-04-25 23:36
|
這個(gè)一般來(lái)說(shuō)都是自己寫的, 很簡(jiǎn)單! 那個(gè)web site不是提供了程序. 它的輸入文件為syml，輸出文件為inp.kpt。(你會(huì)做能帶吧,這個(gè)取點(diǎn)跟能帶的取點(diǎn)宗旨是一樣的.) 其中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 第一行用來(lái)標(biāo)記有多少個(gè)特殊k點(diǎn)，下面是這些特殊k點(diǎn)之間每個(gè)要分多少個(gè)k點(diǎn)，接著就是這些特殊k點(diǎn)的坐標(biāo)。 把下面這個(gè)程序編譯,執(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)站上說(shuō)的很清楚! 如果你想把兩個(gè)k點(diǎn)分成十份, 就把15換成10就可以了! [ Last edited by gavinliu7390 on 2010-4-25 at 22:34 ] |
3樓2010-04-25 22:22:42
mudannanzi
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1959.1
- 帖子: 883
- 在線: 71.5小時(shí)
- 蟲號(hào): 976416
- 注冊(cè): 2010-03-20
- 性別: GG
- 專業(yè): 金屬材料的制備科學(xué)與跨學(xué)
4樓2010-04-25 22:25:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
