计算化学公社

 找回密码 Forget password
 注册 Register
Views: 505|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[其它程序] 使用rescu软件结构优化后,产生mat文件,不会转化为cif

[复制链接 Copy URL]

62

帖子

0

威望

511

eV
积分
573

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
# -*- coding: utf-8 -*-
"""
Extract lattice vectors and atom info from saved mat file
"""
from ase.atoms import Atom, Atoms
from ase.io import write
from pymatreader import read_mat

a = read_mat("H_scf.mat")
latvec = a["domain"]["latvec"]
poses = a["atom"]["xyz"]  # np.ndarray
ele_ind_list = a["atom"]["element"]  # [1,1,2,2]
ele_species_list = a["element"]["species"]  # ['H','Li']

atom_list = []
for i in range(len(ele_ind_list)):
    atom = Atom(ele_species_list[int(ele_ind_list[i]) - 1], poses[i])
    atom_list.append(atom)
ats = Atoms(atom_list, cell=latvec)
print(ats)
print("--> Writing to atoms.cif file...")

write("atoms.cif", ats)
以上为教程里的转化方法,可是我看不懂,有没有哪位大神可以给解释一下,感激。


本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 22:38 , Processed in 0.142694 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list