利用LAMMPS中自带的命令,很难生成比较复杂的模型,此时需要借助外部模拟生成data文件,如使用Materials Studio(MS)
新建项目
File→NEW... →3D ATOMISTICS
![57c37de864ab161ed6a92c39a55f97a9.png](https://i-blog.csdnimg.cn/blog_migrate/a5629fa6449dc4b0110d586e874a654f.jpeg)
绘制环状结构
![b0637d37536ce3d5e450336d65b94904.png](https://i-blog.csdnimg.cn/blog_migrate/165e1e1f039c611b2f9410957c381203.png)
![47952de6d9200fc5582789fa84caee72.png](https://i-blog.csdnimg.cn/blog_migrate/c4d3b422161fed3e05b8e156230d6ac3.png)
右键→Display style 选择Ball and stick(球棍模型) ,其中 Stickradius 为棍的半径,Ball radius 为球的半径
![a2c0903760740f1e407d517bf9d58ce7.png](https://i-blog.csdnimg.cn/blog_migrate/645d74939deea7353043cbaf79772829.png)
![d59ec0257bcea969b7885c9a990d32a5.png](https://i-blog.csdnimg.cn/blog_migrate/1d035a72f7ef818a8c0213281836616b.png)
选择工具栏上的指针可以选中原子或键,按ctrl可以多选
![9f49502f5fdae58d3bf82303a68f4c5f.png](https://i-blog.csdnimg.cn/blog_migrate/25ba4a55afeef492fc09a3bdc47adaad.png)
![c7c9f87951eb5e9057b8a30a09b684f1.png](https://i-blog.csdnimg.cn/blog_migrate/c2cc1c475123041ee330cadb392b0989.png)
选择Modify Bond Type 可以改变原子之间成键的类型,将上述三键改为Double Bond
![00237c0d7457c9d5f3d4b0bce565f4dc.png](https://i-blog.csdnimg.cn/blog_migrate/cffeeb6222d72cdab07d986c1e641480.png)
![37964f8fc4aaf616cf19a060c1ac153e.png](https://i-blog.csdnimg.cn/blog_migrate/234d97df8c79dea8f269093cc14a838d.png)
选择Sketch Atom 绘制单个原子,选择绘制碳原子(按Esc可以在绘制结束后结束绘制)
![fd19c9c47b190fa0f892d3fc3d59b47e.png](https://i-blog.csdnimg.cn/blog_migrate/2084e4ccdbc648ca2454071b3a8ce3ba.png)
再在刚刚绘制的碳原子上绘制一个氧原子,并将碳氧键改为双键
![a69054b32b1bdc95e8014194ad988f7e.png](https://i-blog.csdnimg.cn/blog_migrate/2b2bfd7b1a0d03b6a277efa38dee09eb.png)
![562b16c1cd9ae841ae3f04e40e0d2848.png](https://i-blog.csdnimg.cn/blog_migrate/e49098bd8a2b6917b73a90d3d1e7a58d.png)
![03a7a98f97411265eaa62a60952a8926.png](https://i-blog.csdnimg.cn/blog_migrate/d4d081e23b2855bbfd914a92a789685a.png)
利用Clean工具可以将原子调整为合适的、合理的形状,利用Adjust Hydrogen工具可以给当前原子添加氢原子
![6ab89de2ee5c9b23ce3698ed6c0e2b9f.png](https://i-blog.csdnimg.cn/blog_migrate/aa76c96dfacd999530c1e22a7dc56efa.jpeg)
绘制氧气
用同样的方法可以绘制出氧气分子
![6ad1913450999a004f7ffdc9cbe9968f.png](https://i-blog.csdnimg.cn/blog_migrate/01b0a6d59d3b3e23239ccf0eb0ee5fac.png)
进行混合
利用Amorphous Cell Tools可以将生成的模型按一定比例混合
![c424cfa5ec9167e5d93dc17dfaa0b405.png](https://i-blog.csdnimg.cn/blog_migrate/710eaecfa5ed36dfeca0e7b7ea3d3437.png)
![fc3c606b5707bab4506dd73755ce1bd2.png](https://i-blog.csdnimg.cn/blog_migrate/8ee7ec3d0eed9005383824852310512b.jpeg)
设置好各种分子的数量,和体系的密度,MS会自动生成一个盒子,盒子大小在窗口最后一行可以看见
![0d25c734a26f7f213bf7a550e4f492f2.png](https://i-blog.csdnimg.cn/blog_migrate/a0b511a49397a0e7e3be67adb3e6adbe.jpeg)
此外,在导出数据之前需要添加势函数,选为cvff,charge选为 charge using qeq
![ab37402630a5c217817636cba9d67a46.png](https://i-blog.csdnimg.cn/blog_migrate/fa730dc9900d7899637ccb34e08d9479.jpeg)
导出文件
File → Export
保存类型选择为car 和 cor
![1ed726a31ee2883913e5679acbc60fb0.png](https://i-blog.csdnimg.cn/blog_migrate/1aae262561bff936f9c63779a63b087a.png)
![0015f0202c9d8d76b45bd57122ef6496.png](https://i-blog.csdnimg.cn/blog_migrate/f02258e6cb228ac16dc6d4a8a8b34f8d.png)
在lammps中安装msi2lmp工具
工具所在目录如下:
![01869564c355c1f1692ea4d98f45d9e2.png](https://i-blog.csdnimg.cn/blog_migrate/578014fe73bac1a9cb9cdc8ab6ce6cf9.png)
进入src文件夹进行安装:
cd src
make
安装成功之后如下:
![1f2b67f2997d009326f8fcf43f2866bc.png](https://i-blog.csdnimg.cn/blog_migrate/d7c7db6288059afce6c2f351424bd6e3.png)
进入frc_files文件夹,将MS导出的两个文件夹放入其中
执行命令:
./msi2lmp C7 -i -class I -frc cvff > data.C7
将C7换成你所生成的car文件的名字即可。执行上述命令之后,会生成两个文件:C7.data和data.C7,其中C7.data为所需的模型文件。(具体使用需要看手册)