计算化学公社

标题: 求助lammps中Ba0.67Sr0.33TiO3建模 [打印本页]

作者
Author: zhoulong 时间: 2021-8-30 21:00
标题: 求助lammps中Ba0.67Sr0.33TiO3建模
本帖最后由 zhoulong 于 2021-8-30 21:02 编辑

各位老师好，我想通过lammps命令建立Ba0.67Sr0.33TiO3模型。我的思路是这样的：首先建立BaTiO3模型(10x10x10),然后用Sr替换Ba原子，替换的概率是0.33；所以我的问题是：如何实现这样的操作？各位老师有什么好的办法或者其他思路吗？万分感谢。PS：本人邮箱95483431@qq.com

作者
Author: wzkchem5 时间: 2021-8-30 21:56
写个python脚本就能解决的问题。。。把BaTiO3的结构以任何纯文本格式保存到文件，然后依次读取每一行，如果这一行以Ba开头，那么产生一个0到1之间的随机数，如果这个随机数小于0.33，就把Ba换成Sr。替换过的坐标再用lammps relax一下就行了

作者
Author: zhoulong 时间: 2021-8-31 22:19

wzkchem5 发表于 2021-8-30 21:56
写个python脚本就能解决的问题。。。把BaTiO3的结构以任何纯文本格式保存到文件，然后依次读取每一行，如果 ...

难点在于不会编程语言

作者
Author: wzkchem5 时间: 2021-8-31 23:45

zhoulong 发表于 2021-8-31 15:19
难点在于不会编程语言

即使不会python，用bash也可以做，只不过会python的情况下python做起来可能方便一点。你但凡是在linux下跑的lammps，必然已经至少会一种shell语言了（比如bash）

作者
Author: zhoulong 时间: 2021-9-3 22:00

wzkchem5 发表于 2021-8-31 23:45
即使不会python，用bash也可以做，只不过会python的情况下python做起来可能方便一点。你但凡是在linux下 ...

好的，谢谢您的提示。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)