lammps .eam势文件 metal单位制转real单位制的python程序

最新推荐文章于 2025-09-12 12:38:11 发布

原创

最新推荐文章于 2025-09-12 12:38:11 发布 · 1.7k 阅读

收录于

当前文章被以下社区和专栏收录：

本文介绍了一个用于将原子模拟中使用的势能文件从一种能量单位转换到另一种能量单位的Python程序。该程序读取原始文件，进行必要的单位换算，并将转换后的数据写入新文件。涉及的换算包括从电子伏特每摩尔到千卡每摩尔。

1.程序

测试内容有限，不保证适用于所有文件

# 1 eV/mol = 1.6021766208×10^-19 J
# 1 J*6.02214076×10^23 = J/mol
# 1KJ/mol = 0.2389029576186153 Kcal/mol
C= 1.6021766208*6.02214076*10*0.2389029576186153

file_name_read = "Cu_u3.eam"
file_name_write = "/home/zxy/Documents/lammps/eam_test/Cu_u3_real.eam"


import re
import os
import numpy as np
import math 
import time
os.system("clear")
f_r= open(file=file_name_read,mode= "r")
f_w= open(file=file_name_write,mode= "w")
print("Start read file...")

# ------------------------- 读取原子信息 ------------------------
line= f_r.readline()
f_w.write(time.asctime( time.localtime(time.time()) )+'\n')
line= f_r.readline()
f_w.write(line)
patten= r'([0-9\.]+)' 
match= re.findall(patten, line)
atom_num= int(match[0])
atom_mass= float(match[1])
atom_lattric_c= float(match[2])

patten= r'([a-zA-Z]+)' 
match= re.findall(patten, line)
atom_lattric_t= match[0]
print("atom_number= %d; atom_mass= %f; lattric_constant=

标签