版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

hakuna

木蟲(chóng) (知名作家)

1ST強(qiáng)帖: 3
應(yīng)助: 82 (初中生)
金幣: 8735.9
散金: 158
紅花: 126
帖子: 5599
在線: 1367.2小時(shí)
蟲(chóng)號(hào): 167622
注冊(cè): 2006-01-15
專業(yè): 考古理論

[交流] 【分享】一個(gè)繪制band的python小程序已有4人參與

網(wǎng)上淘來(lái)的，向作者致敬！
不做能帶，不懂python，未經(jīng)測(cè)試

原文地址：http://people.vanderbilt.edu/~brandon.g.cook/vaspbandplot/vband.py

CODE:

!/usr/bin/env python

#make a simple plot of band structure from vasp calculation

# -Brandon Cook April 7, 2010

from numpy import *

import pylab

class OUTCAR():

    def __init__(self):

        f = open('OUTCAR', 'r')

        s = f.read()

        f.close()

        i = s.find('E-fermi')

        self.ef = float(s[i:i+19].split()[-1])

        i = s.find('NBANDS')

        self.nbands = int(s[i:i+15].split()[-1])

        i = s.find('NKPTS')

        self.nkpts = int(s[i:i+15].split()[-1])

class KPOINTS():

    def __init__(self):

        # kpts file should be in line mode

        # only take 1 line for now

        f = open('KPOINTS', 'r')

        s = f.readlines()

        f.close()

        l1 = s[-2].split()

        l2 = s[-1].split()

        self.p1 = array([ l1[0], l1[1], l1[2] ], float)

        self.p2 = array([ l2[0], l2[1], l2[2] ], float)

class BandStructure():

    # Just take all the points and make a scatter plot

    # todo: connect the lines

    def __init__(self):

        f = open('EIGENVAL', 'r')

        s = f.readlines()

        f.close()

        oc = OUTCAR()

        kp = KPOINTS()

        axis = kp.p2 - kp.p1

        axis /= linalg.norm(axis)

        nbands = oc.nbands

        nkpts  = oc.nkpts

        kpts = zeros( nkpts*nbands )

        en   = zeros( nkpts*nbands )

        istart = 7

        for k in range(0, nkpts):

            ik = istart + k*(nbands+2)

            kpt = s[ik].split()

            kpt = array([kpt[0], kpt[1], kpt[2]], float)

            x = dot(kpt, axis)

            i1 = k*nbands

            i2 = (k+1)*nbands

            kpts[i1:i2] = x

            for ib in range(0, nbands):

                j1 = ik + 1 + ib

                j2 = i1 + ib

                en[j2] = float(s[j1].split()[-1])

        pylab.scatter(kpts, en, s=1)

        # add a line for the fermi energy

        xmin = 0.0

        xmax = linalg.norm(kp.p2 - kp.p1)

        x = linspace(xmin, xmax, num=100, endpoint=True)

        y = oc.ef * ones(100)

        pylab.plot(x,y)

        pylab.show()

o = OUTCAR()

print 'fermi energy, nbands, nkpts'

print o.ef,o.nbands,o.nkpts

b = BandStructure()

回復(fù)此樓

» 猜你喜歡

262求調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑已經(jīng)有6人回復(fù)
張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 已經(jīng)有6人回復(fù)
考研生物與醫(yī)藥調(diào)劑已經(jīng)有3人回復(fù)
274求調(diào)劑已經(jīng)有5人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有26人回復(fù)
334求調(diào)劑已經(jīng)有3人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有5人回復(fù)
看看我！301求調(diào)劑 081700 已經(jīng)有3人回復(fù)
求調(diào)劑生物學(xué) 377分已經(jīng)有3人回復(fù)

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

推薦幾個(gè)vasp處理計(jì)算結(jié)果的python腳本已經(jīng)有66人回復(fù)
Python寫的arcgis中點(diǎn)坐標(biāo)txt轉(zhuǎn)面狀shp的小程序已經(jīng)有90人回復(fù)
求一個(gè)簡(jiǎn)單的讀寫文件的小程序。。。已經(jīng)有4人回復(fù)
用phonopy計(jì)算聲子譜已得到FORCE_SETS，之后我怎么得到聲子譜？已經(jīng)有7人回復(fù)
《python科學(xué)計(jì)算》已經(jīng)出版了已經(jīng)有36人回復(fù)
運(yùn)行一個(gè)python程序，出現(xiàn)問(wèn)題，請(qǐng)高手賜教已經(jīng)有9人回復(fù)
【原創(chuàng)】處理elk能帶的程序已經(jīng)有16人回復(fù)
【求助】求助安裝abinitbandstructuremaker.py? 已經(jīng)有3人回復(fù)

1樓 2010-12-26 17:32:54

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

gavinliu7390

木蟲(chóng) (著名寫手)

葉落鷹飛

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

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

執(zhí)行不了，錯(cuò)誤很多。

贊一下(1人)

回復(fù)此樓

真理是一點(diǎn)點(diǎn)接近的！

3樓2010-12-27 10:44:18

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

查看全部 5 個(gè)回答

westmonster

金蟲(chóng) (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 5637.3
帖子: 401
在線: 236.4小時(shí)
蟲(chóng)號(hào): 456224
注冊(cè): 2007-11-11
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

pylab是嘛？孤陋寡聞了，改改還是能用的，謝謝樓主分享。

贊一下

回復(fù)此樓

蟬，于其脫殼展翅，則蟄居地下，似無(wú)聲無(wú)息，實(shí)則以備有所為。

4樓2010-12-27 12:54:53

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

akakcolin

金蟲(chóng) (著名寫手)

應(yīng)助: 3 (幼兒園)
金幣: 685.9
散金: 51
紅花: 2
帖子: 1010
在線: 524.6小時(shí)
蟲(chóng)號(hào): 321213
注冊(cè): 2007-03-10
專業(yè): 力學(xué)

python 是個(gè)好語(yǔ)言啊

贊一下

回復(fù)此樓

為人應(yīng)當(dāng)誠(chéng)實(shí)正直不能心懷叵測(cè)心眼明亮才能迎來(lái)幸福避開(kāi)災(zāi)禍盲目者若有人指點(diǎn)才不會(huì)迷失道路若單獨(dú)上路就難以保證不誤入岐途

5樓2010-12-27 23:09:27

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

查看全部 5 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 313求調(diào)劑 +3	賣個(gè)關(guān)子吧 2026-03-31	3/150	2026-03-31 08:25 by yyyyry
[考研] 285求調(diào)劑 +6	AZMK 2026-03-30	9/450	2026-03-31 08:01 by jp9609
[考研] 材料求調(diào)劑一志愿哈工大總分298分，前三科223分 +10	dongfang59 2026-03-27	10/500	2026-03-30 23:42 by 果果媽咪
[考研] 調(diào)劑 +4	GK72 2026-03-30	4/200	2026-03-30 20:32 by dick_runner
[考研] 322求調(diào)劑 +10	宋明欣 2026-03-27	10/500	2026-03-30 18:47 by 544594351
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +9	XX小鄧 2026-03-29	9/450	2026-03-30 17:18 by limeifeng
[考研] 286求調(diào)劑 +5	Faune 2026-03-30	5/250	2026-03-30 17:15 by 無(wú)際的草原
[考研] 303求調(diào)劑 +7	DLkz1314. 2026-03-30	7/350	2026-03-30 16:05 by shuang5186
[考研] 085701求調(diào)劑初試286分 +5	secret0328 2026-03-28	5/250	2026-03-30 12:54 by fangnagu
[考研] 291求調(diào)劑 +5	Y-cap 2026-03-29	6/300	2026-03-29 13:18 by mumin1990
[考研] 352分-085602-一志愿985 +5	海納百川Ly 2026-03-29	5/250	2026-03-29 09:57 by Sjndkwm
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 壓國(guó)家一區(qū)線，求導(dǎo)師收留，有恩必謝！ +7	迷人的哈哈 2026-03-28	7/350	2026-03-28 16:47 by 催化大白
[考研] 085701環(huán)境工程，267求調(diào)劑 +16	minht 2026-03-26	16/800	2026-03-28 12:16 by zllcz
[考研] 085405 考的11408求各位老師帶走 +3	Qiu學(xué)ing 2026-03-28	3/150	2026-03-28 09:19 by 樂(lè)呵呵的追夢(mèng)人
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 275求調(diào)劑 +10	jjjjjjjjjjl 2026-03-27	10/500	2026-03-27 23:47 by barnett0632
[考研] 一志愿西北大學(xué) 總分282 英語(yǔ)一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

hakuna

[交流] 【分享】一個(gè)繪制band的python小程序 已有4人參與

» 猜你喜歡

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

gavinliu7390

westmonster

akakcolin

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 【分享】一個(gè)繪制band的python小程序已有4人參與

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