跳至主要内容

OVITO技巧+MS 旋转 redefine 格式化元素序号与力场一致,vasp文件格式元素替换

 问题:ovito读取数据后默认将最小的边长设置成align X轴(因为Restricted triclinic simulation cell format),但是碰到一些情况我们需要将其余两个晶格align X/Y/Z 目的为了沿着这个晶格矢量做变形,这时候需要用到affine transformation这个模块中的transform to target cell,这个矩阵的信息可以由MS很方便获得,MS可以随意旋转或者格式化标准取向 然后在properties里的Lattice 3D可以读取VectorA/B/C即可 很方便,输入到ovtio后输出data文件的时候 不要勾选Restricted triclinic simulation cell format,不然的话又回到标准格式了。

而且一切结构都可以正交化 但是会牺牲掉另外两个不align坐标轴晶格方向的周期性最后用P S S边界条件,在ovito里面通过slice切出矩形的柱就行了,通过G(r)观察最终结构矩形柱的腰围略大一点点避免原子靠近。

问题:自己构建的或者研究的结构元素序号与力场不一致,需要手动修改结构文件里的元素序号么?

文件:需要POSCAR(.vasp)文件里的元素名称,所以整个过程需要POSCAR格式文件(vesta可以保留元素名称)。

步骤:

1.将需要格式化的文件转成POSCAR

2.复制一个副本POSCAR,并且直接文件编辑里面的元素名称一行,到准确的顺序,这个只是作为模板用,不用考虑结构合理性问题。

3.将这个副本文件拖入OVITO,记得是完全从新导入,可以看到particle types 内容是正确的顺序。

4.最后将需要格式化的文件全部转换成POSCAR,然后拖入OVITO选择保留现在的设置(即覆盖原结构文件即可)。元素就自动格式化成想要的顺序了.

快速实现元素替换和修改,采用OVITO转存为vasp格式,直接修改元素即可,后续再保存为data格式


ovito解决位移随时间变化曲线,难点是时间函数,采用Generate trajectory lines+Scatter plot 可以解决.

OVITO-Displacement vectors 无法识别unwarp结果,导致位移MSD统计有误,建议用trajectory lines给出绝对坐标重新计算或者后续用 lammps输出xu yu zu等坐标值

Data Table Text File(需要上面的技巧Scatter plot画出所需要的数据 再后续python脚本简单处理即可): 最后Data Table部分特指画图部分的数据可以单独输出每一帧。

Table of value(最常见直接):中间Global Attributes里的值与时间的关系

Data Plot File: 图片格式pdf或者png


评论

此博客中的热门博文

lammps 压痕划痕模拟设置参考

  Molecular dynamics study on the effect of electric current on electrically-assisted scratching for crystal copper - IOPscience 原因深入分析如下: ✅ 切削 / 摩擦 / 划痕:局部剧烈变形 → 热量集中 这类过程模拟的是工具与材料 接触区域的强烈局部非平衡过程 ; 如果对整个系统控温,会 严重抹平局部发热、滑移带的应变能耗散 等重要现象; 所以 只在边界区域(如底部、侧边)设 thermostat,起到“热沉”作用 ; 文献经典设置就是: 底部固定 ; 边缘 slab 控温 ; 接触区完全不控温,自由演化 。 ✅ 拉伸 / 压缩 / 剪切:全局加载 → 热传导充分 是材料整体在受力,不存在特别“集中”的能量输入区域; 局部发热相对温和,且在 bulk 系统中可以通过自身结构进行导热 ; 实验中常常是等温加载(准静态过程); 所以 很多文献就直接用整体 fix nvt 控温 ,保持恒温环境,简化模拟; 注意有些更精细的研究会改为: 只在两端 slab 控温,中间 Newtonian 自由演化 。 📚 二、典型模拟场景下的控温策略总览 场景类别 控温方式 控温范围 控温方法 控温目的 注意事项 ✅ 平衡态热力学性质 (如热容、扩散、应力) 整体控温 全体系 fix nvt , fix npt 模拟室温等温状态 标准EMD方法 ✅ 热导率(Green-Kubo) 整体控温 全体系 fix nvt (前期平衡), 后期 nve 采集能流自相关函数 采样期不能控温 ✅ 热导率(NEMD) 区域控温 热源/热沉 fix langevin , fix heat 人为施加温差形成热流 中间区自由演化 ✅ 拉伸 / 压缩 / 剪切 整体控温(常用)或局部控温 全体系或上下 slab fix nvt 或 langevin slab 保持恒温,避免非真实升温 全控温可抹平热应变 ✅ 应力松弛 / 加热冷却过程 整体控温 全体系 fix nvt 或温度渐变 等温退火、升温或冷却 控温方式决定退火速率 ✅ 位错运动 / 缺陷扩散 局部控温 边界或部分 slab ...

lammps Pdamp,Tdamp的设置经验

 一张小抄(固体/位错/二维材料都适用) fix             11    all npt temp 0.1 0.1 0.5 tri 0.0 0.0 5   drag 2 tchain 3 pchain 3 保持 Pdamp ≫ Tdamp(通常 10× 左右)。 固体/低温:Tdamp 取 0.5–1 ps,Pdamp 取 5–15 ps;需要更稳就再加大 Pdamp。 所以推荐如下: 0.1K 用1 10 300K 固体 用0.5 5  高温用0.2 2 液体用0.1 1 液体/高温:Tdamp 0.2–0.5 ps,Pdamp 2–5 ps 往往够。 2D(石墨烯等,只控平面 x/y):Pdamp 常用 10–20 ps 起步,z 固定或 z NULL。 观察到体积/压力振荡大:增加 Pdamp 或加 drag 2–3,必要时把 dt 临时降到 0.5 fs。 drag 2、tchain/pchain 3 保留没坏处,确实能再抑制一点振荡;不是硬性必须,但在固体+低温+(可能还有 tri 或剪切)的组合里,“更稳”>“更快”,我一般会开着。

dpgen训练经验

最新的dpgen参考PtCuP /work/qan/kluo/PtCuP 0.1K的第0代采样很重要 可以多重复几次确保99以上的准确率,它是后续高温的基础 单点能计算  ISYM = 2 nohup dpgen run param.json machine.json 1>log 2>err& nohup dpgen init_bulk param.json machine.json  1>log 2>err& 初始数据集产生,只能一个POSCAR的计算 但是可以同时提多个任务,但是每个任务都需要 POTCAR POSCAR 一 一对应才行。  elements 和POSCAR POTCAR保持一致。 " type_map " : [ "Ti" , "C" , "V" , "Cr" , "Nb" , "Mo" ],都写全才行 POSCAR 不需要 改成特定顺序 程序最终生成数据集的时候会根据 type_map自动统一匹配 usage: dpgen [-h] {init_surf,init_bulk,auto_gen_param,init_reaction,run,run/report,collect,simplify,autotest,db,gui} ... dpgen is a convenient script that uses DeepGenerator to prepare initial data, drive DeepMDkit and analyze results. This script works based on several sub-commands with their own options. To see the options for the sub-commands, type "dpgen sub-command -h". positional arguments:   {init_surf,init_bulk,auto_gen_param,init_reaction,run,run/report,collect,simpli...