 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 2396 | 回復(fù): 3 | |||
| 【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 cnzs 的 7 個(gè)金幣 ，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì) | |||
[交流]
VASP數(shù)據(jù)使用ovito官網(wǎng)腳本分析voronio
|
|||
|
1.能不能用ovito處理AIMD的XDATCAR文件？ 2.如果可以為什么我用官網(wǎng)提供的Python腳本只能得出第一個(gè)離子步的數(shù)據(jù)? 3.有沒有大佬可以幫忙改進(jìn)一下Python腳本？ # Import OVITO modules. from ovito.io import * from ovito.modifiers import * from ovito.pipeline import * # Import NumPy module. import numpy # Load a simulation snapshot of a Cu-Zr metallic glass. pipeline = import_file("D:/ovito/400/XDATCAR",multiple_frames = True) # Set atomic radii (required for polydisperse Voronoi tessellation). atom_types = pipeline.source.data.particles['Particle Type'].types atom_types[0].radius = 1.28 # atomic radius (atom type 1 in input file) atom_types[1].radius = 1.55 # atomic radius (atom type 2 in input file) # Set up the Voronoi analysis modifier. voro = VoronoiAnalysisModifier( compute_indices = True, use_radii = True, edge_threshold = 0.1 ) pipeline.modifiers.append(voro) # Let OVITO compute the results. data = pipeline.compute() # Access computed Voronoi indices. # This is an (N) x (M) array, where M is the maximum face order. voro_indices = data.particles['Voronoi Index'] # This helper function takes a two-dimensional array and computes a frequency # histogram of the data rows using some NumPy magic. # It returns two arrays (of equal length): # 1. The list of unique data rows from the input array # 2. The number of occurences of each unique row # Both arrays are sorted in descending order such that the most frequent rows # are listed first. def row_histogram(a): ca = numpy.ascontiguousarray(a).view([('', a.dtype)] * a.shape[1]) unique, indices, inverse = numpy.unique(ca, return_index=True, return_inverse=True) counts = numpy.bincount(inverse) sort_indices = numpy.argsort(counts)[::-1] return (a[indices[sort_indices]], counts[sort_indices]) # Compute frequency histogram. unique_indices, counts = row_histogram(voro_indices) # Print the ten most frequent histogram entries. for i in range(10): print("%s\t%i\t(%.1f %%)" % (tuple(unique_indices), counts, 100.0*float(counts)/len(voro_indices))) result： (0, 0, 0, 2, 8, 4, 0, 0) 11 (5.5 %) (0, 0, 0, 3, 6, 4, 0, 0) 7 (3.5 %) (0, 0, 0, 0, 12, 0, 0, 0) 6 (3.0 %) (0, 0, 0, 1, 10, 3, 0, 0) 5 (2.5 %) (0, 0, 0, 2, 8, 1, 0, 0) 5 (2.5 %) (0, 0, 0, 4, 6, 4, 0, 0) 5 (2.5 %) (0, 0, 0, 2, 8, 5, 0, 0) 5 (2.5 %) (0, 0, 0, 1, 10, 2, 0, 0) 5 (2.5 %) (0, 0, 0, 2, 8, 0, 0, 0) 4 (2.0 %) (0, 0, 0, 4, 4, 6, 0, 0) 3 (1.5 %) |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有116人回復(fù)
-
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
上海優(yōu)秀男生誠征結(jié)婚女友符合微信聯(lián)系
+3/784
-
海南大學(xué)肖永昊老師團(tuán)隊(duì)招收2026年博士研究生（第二批）
+5/555
-
面對核三廢，威格屏蔽手套箱如何筑牢第一道安全防線？
+1/90
-
招收調(diào)劑生，動(dòng)物醫(yī)學(xué)、水產(chǎn)養(yǎng)殖、水產(chǎn)動(dòng)物醫(yī)學(xué)等專業(yè)
+1/89
-
Journal of Environmental Management投稿
+1/83
-
江西南部某一本高校招收環(huán)境、生物、計(jì)算機(jī)、農(nóng)業(yè)等專業(yè)調(diào)劑生
+1/81
-
中國科學(xué)院-材料計(jì)算小組尋求2026級(jí)研究生（包括調(diào)劑生/聯(lián)合培養(yǎng)等）
+5/60
-
中國石油大學(xué)（北京）國家級(jí)大人才團(tuán)隊(duì)博士招生3名：化學(xué)、材料、石油工程：油田化學(xué)
+1/30
-
坐標(biāo)北京，86年女
+1/25
-
申請取消論壇的歸檔操作
+1/20
-
中國科學(xué)院蘇州納米所院士團(tuán)隊(duì)博士后崗位招聘
+1/20
-
接收2026調(diào)劑研究生
+1/18
-
青島科技大學(xué)0860 一志愿./調(diào)劑 招生
+1/11
-
魯東大學(xué)水產(chǎn)養(yǎng)殖0908 有調(diào)劑名額，需要調(diào)劑聯(lián)系
+1/9
-
26 fall & 27 spring考博求助，控制/機(jī)器人/人工智能方向
+1/5
-
福建師范大學(xué)化學(xué)與材料學(xué)院杜克釗團(tuán)隊(duì)招生
+1/5
-
2026年東北石油大學(xué)“頁巖油氣鉆采高效井眼清潔”創(chuàng)新團(tuán)隊(duì)招碩士生
+1/5
-
復(fù)試可能過不了想調(diào)劑
+1/4
-
2026年東北石油大學(xué)“頁巖油氣鉆采高效井眼清潔”創(chuàng)新團(tuán)隊(duì)招生
+1/3
-
095400 森林培育 303求調(diào)劑
+1/2
4樓2022-06-22 14:07:45
簡單回復(fù)
2020-09-16 14:37
回復(fù)
cnzs(金幣+1): 謝謝參與
1 發(fā)自小木蟲IOS客戶端
psylhh3樓
2020-09-17 17:32
回復(fù)
cnzs(金幣+1): 謝謝參與
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cnzs 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 282分材料專業(yè)求調(diào)劑院校 +10 | 楓橋ZL 2026-03-09 | 13/650 |
|
|---|---|---|---|---|
|
[碩博家園] 2026級(jí)碩士研究生招生/調(diào)劑 +3 | 知足常樂的樂 2026-03-06 | 5/250 |
|
|
|
[考研] 材料調(diào)劑，307分 +8 | 張泳銘1 2026-03-09 | 8/400 |
|
|
|
[考研] 調(diào)劑 +4 | 調(diào)劑的考研學(xué)生 2026-03-09 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級(jí)已過 +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|
|
[考研] 求調(diào)劑，數(shù)一英一274分 +4 | 小菲會(huì)努力 2026-03-08 | 4/200 |
|
|
|
[考研] 337求調(diào)劑 +3 | 睡醒，。 2026-03-09 | 3/150 |
|
|
|
[考研] 材料工程330分求調(diào)劑，一志愿985 +3 | 小材化本科 2026-03-07 | 3/150 |
|
|
|
[考研] 22408 275分求調(diào)劑 +3 | 宇智波比 2026-03-03 | 4/200 |
|
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
|
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3 | 0703楊悅305分 2026-03-05 | 3/150 |
|
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
|
[考研] 282求調(diào)劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|