版塊導(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)劑公告

返回列表

hakuna

木蟲 (知名作家)

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

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

網(wǎng)上淘來的，向作者致敬！
不做能帶，不懂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ù)此樓

» 猜你喜歡

287求調(diào)劑已經(jīng)有7人回復(fù)
325求調(diào)劑已經(jīng)有5人回復(fù)
343求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑推薦材料 304 已經(jīng)有4人回復(fù)
316求調(diào)劑已經(jīng)有4人回復(fù)
351求調(diào)劑已經(jīng)有4人回復(fù)
349求調(diào)劑已經(jīng)有5人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有7人回復(fù)
342求調(diào)劑已經(jīng)有3人回復(fù)
調(diào)劑求收留已經(jīng)有7人回復(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)問題，請(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的回帖

sg18408926

至尊木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 10945.9
散金: 246
紅花: 1
帖子: 1075
在線: 431.4小時(shí)
蟲號(hào): 727956
注冊(cè): 2009-03-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

不錯(cuò)啊

回復(fù)此樓

2樓2010-12-26 20:43:52

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

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é)和

★
小木蟲(金幣+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的回帖

westmonster

金蟲 (正式寫手)

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

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

贊一下

回復(fù)此樓

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

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

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

akakcolin

金蟲 (著名寫手)

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

python 是個(gè)好語言啊

贊一下

回復(fù)此樓

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

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 hakuna 的主題更新

返回列表

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

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

[考研] 考研調(diào)劑 +8	呼呼？~+123456 2026-03-24	8/400	2026-03-26 23:04 by 不吃魚的貓
[考研] 305求調(diào)劑 +4	哇盧卡庫(kù) 2026-03-26	4/200	2026-03-26 20:58 by sanrepian
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 求調(diào)劑 +6	白QF 2026-03-21	6/300	2026-03-26 20:37 by fmesaito
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4	Higraduate 2026-03-22	8/400	2026-03-26 20:34 by Higraduate
[考研] 08開頭275求調(diào)劑 +3	拉誰不重要 2026-03-26	3/150	2026-03-26 20:22 by barlinike
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6	小小陳小小 2026-03-22	9/450	2026-03-26 15:58 by dick_runner
[考研] 0856調(diào)劑 +4	求求讓我有書讀?/a> 2026-03-26	5/250	2026-03-26 15:54 by dick_runner
[考研] 材料求調(diào)劑 +4	.m.. 2026-03-25	4/200	2026-03-25 21:30 by peike
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 0854電子信息求調(diào)劑 324 +4	Promise-jyl 2026-03-23	4/200	2026-03-25 11:36 by Sugarlight
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 070300化學(xué)求調(diào)劑 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 求調(diào)劑 +6	研研，接電話 2026-03-24	7/350	2026-03-24 17:01 by barlinike
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992

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

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

hakuna

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

» 猜你喜歡

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

sg18408926

gavinliu7390

westmonster

akakcolin

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

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

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