| 查看: 5699 | 回复: 8 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 meizhixin 的 10 个金币 | |||
[交流]
ovito软件分析Voronoi
|
|||
|
使用ovito分析Voronoi，输入参数后，怎么得到结果？不明白？有会的吗？ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
兄弟们帮我看看这两个结构怎么合成
已经有5人回复
-
一审一个审稿人，小修，会怎么样呀？
已经有6人回复
-
院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度
已经有7人回复
-
没.赶.上.今.年.博.士.的.同.学.看.过.来.
已经有10人回复
-
美国教授开除中国研究生：我就不该录取你！
已经有6人回复
-
某系主任被同校教授举报包养情妇、长期嫖娼
已经有14人回复
-
特别资助审核状态
已经有6人回复
-
国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？
已经有5人回复
-
杭州国企和浙江高校如何选择？
已经有20人回复
-
发个文章，结果被某学报请的外审PUA了
已经有8人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25588.7
- 散金: 602
- 红花: 166
- 帖子: 3760
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
meizhixin: 金币+10 2017-02-28 13:32:37
meizhixin: 金币+10 2017-02-28 13:32:37
|
采用python script来计算吧： https://www.ovito.org/manual/pyt ... ute-voronoi-indices |
2楼2017-02-28 12:52:46
3楼2017-02-28 13:31:52
|
老师好，我把路径改了，在ovito图形界面，run scripts，但是ovito图形界面只是闪了下，没报错，也没结果，啥也没有，是脚本的问题吗？我用的是vasp的POSAR文件，下面是ovito计算voro的脚本，我直接复制的，不知道对不对 # Import OVITO modules. from ovito.io import * from ovito.modifiers import * # Import NumPy module. import numpy # Load a simulation snapshot of a Co-B metallic glass. node = import_file("E:/ovito/Co65B35/POSCAR"  # Set atomic radii (required for polydisperse Voronoi tessellation). atypes = node.source.particle_properties.particle_type.type_list atypes[0].radius = 1.25 # Co atomic radius (atom type 1 in input file) atypes[1].radius = 0.82 # B atomic radius (atom type 2 in input file) # Set up the Voronoi analysis modifier. voro = VoronoiAnalysisModifier( compute_indices = True, use_radii = True, edge_count = 6, # Length after which Voronoi index vectors are truncated edge_threshold = 0.1 ) node.modifiers.append(voro) # Let OVITO compute the results. node.compute() # Make sure we did not lose information due to truncated Voronoi index vectors. if voro.max_face_order > voro.edge_count: print("Warning: Maximum face order in Voronoi tessellation is {0}, " "but computed Voronoi indices are truncated after {1} entries. " "You should consider increasing the 'edge_count' parameter to {0}." .format(voro.max_face_order, voro.edge_count)) # Note that it would be possible to automatically increase the 'edge_count' # parameter to 'max_face_order' here and recompute the Voronoi tessellation: # voro.edge_count = voro.max_face_order # node.compute() # Access computed Voronoi indices as NumPy array. # This is an (N)x(edge_count) array. voro_indices = node.output.particle_properties['Voronoi Index'].array # 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))) |
4楼2017-03-09 14:30:58
tongaigao
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 108.5
- 帖子: 33
- 在线: 18.7小时
- 虫号: 4978420
- 注册: 2016-09-04
- 性别: MM
- 专业: 金属非晶态、准晶和纳米晶
5楼2017-05-18 15:49:20
gtssongchi
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4373.1
- 散金: 136
- 红花: 2
- 帖子: 157
- 在线: 99.5小时
- 虫号: 1137504
- 注册: 2010-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2017-11-15 08:31:05
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
7楼2019-03-13 08:16:21
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
8楼2019-04-17 19:13:01
cuso4852cuoh
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 223.4
- 散金: 20
- 帖子: 170
- 在线: 49.6小时
- 虫号: 1481571
- 注册: 2011-11-08
- 专业: 工程热物理相关交叉领域
9楼2019-04-17 19:13:52
10