| 查看: 2401 | 回复: 3 | |||
| 【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 cnzs 的 7 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会 | |||
[交流]
VASP数据使用ovito官网脚本分析voronio
|
|||
|
1.能不能用ovito处理AIMD的XDATCAR文件? 2.如果可以为什么我用官网提供的Python脚本只能得出第一个离子步的数据? 3.有没有大佬可以帮忙改进一下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 %) |
» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
计算三阶力常数用thirdorder扩胞时参数如何选择?临近原子数怎么取值
已经有7人回复
物理学I论文润色/翻译怎么收费?
已经有227人回复
2026年循环经济功能材料国际会议(ICFMCE 2026)
已经有0人回复
2026年第五届电气、电子与信息工程国际会议(ISEEIE 2026)
已经有0人回复
哈尔滨理工大学物理系招收物理学考研调剂
已经有3人回复
0702一志愿吉大B区求调剂
已经有5人回复
求调剂
已经有0人回复
0702一志愿吉大B区求调剂有论文
已经有0人回复
请问还有没有用Latex写文章的小伙伴们?
已经有0人回复
» 抢金币啦!回帖就可以得到:
冷冻行星球磨机、低温高温球磨机
+1/87
招聘有机合成人员,新成立部门,上升空间大
+1/79
福建师范大学化学与材料学院杜克钊团队招生
+1/40
湖北省属重点大学招收2026年材料、化学、物理背景硕士
+2/28
北京理工大学-化学与化工学院-招收2026级博士生 [申请-考核制]
+1/21
学术/专业硕士调剂——化学、资源与环境
+1/16
招博士,申请考核制,环境化学、理论计算、环境毒理/健康方向
+1/15
湖南大学材料科学与工程学院招收博士研究生
+1/9
海南大学邱天宇老师团队招收博士生与招聘博士后
+5/5
【神经生物学考博求助】
+1/5
请问为什么我的贴子自动到“已删除的主题”里而没有被发布?
+1/5
大湾区大学-哈尔滨工业大学(深圳)联培博士招生
+1/3
欢迎报考浙江工业大学埃克塞特联合学院硕士研究生
+1/3
温州大学化材学院王娟课题组招生
+1/2
333求调剂
+1/2
东北农业大学(211)水利工程招生
+1/2
重庆交通大学山区桥梁及隧道工程国家重点实验室——Smart Concrete Lab招生
+1/2
南方科技大学基础免疫与微生物学系招聘科研助理1-2名,从事微生物与免疫学方向研究
+1/2
招收环工、环科、化工、应化、化学、生工、生物、材化等相关专业研究生,要求过国家线
+1/2
澳大利亚科学院院士Prof.Gordon Wallace, University of Wollongong 全奖博士机会
+1/1
4楼2022-06-22 14:07:45
简单回复
2020-09-16 14:37
回复
cnzs(金币+1): 谢谢参与
1 发自小木虫IOS客户端
psylhh3楼
2020-09-17 17:32
回复
cnzs(金币+1): 谢谢参与













回复此楼