LAMMPS计算 LiS体膨胀系数 | Bohrium-玻尔科研空间站

空间站广场

论文

Notebooks

比赛

课程

Apps

我的主页

我的Notebooks

我的论文库

我的足迹

我的工作空间

任务

节点

文件

数据集

镜像

项目

数据库

公开

LAMMPS计算 LiS体膨胀系数

《计算材料学》组队共读

1310819006@qq.com

发布于 2024-03-21

推荐镜像 :Basic Image:ubuntu:22.04-py3.10-cuda12.1

推荐机型 :c2_m4_cpu

双击即可修改

代码

文本

[ ]

%%writefile in.LiS

#能源材料模拟利分子动力学模拟电池材料

#清洁能源，电池，太阳能、氢能、风能

#LiS电池体积膨胀的模拟：使用LAMMPS模拟S电极随Li嵌入而产生的体积膨胀。

#尝试撰写 Notebook，描述你的LAMMPS脚本及结果分析方法。

#2019年诺贝尔化学奖

#s电极随着Li的嵌入产生的体积膨胀，

#建立晶胞后插入不同数量的Li离子

#升温融化和快速退火

#在NPT下进行计算体积膨胀率

#寻找Li-S势函数

#openkim 模型包含电荷信息，随机插入Li

units real

boundary p p p

atom_style charge

read_data ****.data

pair_style reax/c NULL checkqeq no #检查电荷收敛性，NO 为不检查

pair_coeff * * ffiled.reax.Li S S Li

neighbor 3.0 bin

min_style cg

minimize 1e -4 1e-6 100 1000

velocity all create 300 31231

thermo 1000

thermo_style custom step temp press

density vol

compute myRDF all rdf 50 1 1 1 2 2 1 2 2 #S Li之间的RDF径向分布函数厚度为dr的区间的原子密度与总密度的比值

#表征有序无序，中子散射可以测量

fix 2 all ave/time 100 10 10000 c_myRDF[*] file initial.rdf mode vector

fix 3 all nvt temp 300 300 10

run 10000

unfix 2

unifx 3

fix 4 all npt temp 300 1800 100 iso 1.0 1.0 1000 drag 1.0

run 50000

unifx 4

fix 5 npt all temp 1800 1800 100 iso 1.0 1.0 1000 drag 1.0

fix 6 all ave/time 100 10 5000 c_myRDF[*]file melted.rdf mode vector

run 5000

unfix 5

unfix 6

fix 7 all npt temp 1800 300 1 iso 1.0 1.0 1000 drag 0.1

thermo 10

run 100

unifx 7

fix 8 all ave/time 100 10 5000 c_myRDF[*] file amorphous.rdf mode vector

fix 9 all npt temp 300 300 100 iso 1.0 1.0 1000 drag 1.0

thermo 1000

run 50000

代码

文本

《计算材料学》组队共读

点个赞吧