计算化学公社

标题: (已解决)聚合物限制性剪切建模问题 [打印本页]

作者
Author: 15939474177 时间: 2026-3-9 16:33
标题: (已解决)聚合物限制性剪切建模问题
本帖最后由 15939474177 于 2026-3-21 09:59 编辑

1.建立环氧树脂和固化剂AC packing盒子，然后交联之前加30真空层，交联完成以后再删除真空层（第一次交联后做界面建模有原子断开，所以调整了方法）
2.对盒子采用结构优化，收敛，然后进行ANNEAL NVT系综（五次循环，300-500k，600ps退火），抽出结果中能量最低结构，继续进行NPT动力学弛豫，298k，标准大气压，500ps过程。
3.抽出最低能量构型，跟铁分子建立LAYER三层做限制性剪切，结果如下，中间层整个跟随上层移动是怎么回事？

作者
Author: yuzc 时间: 2026-3-10 15:39
shear的原理在于空间切片+动量交换，这里不涉及键的事情。
所以我这里有几种可能供你参考：
1. 你不小心手滑点击了calculate bonds；
2. 交联脚本聚合度没做上去，里面有寡聚物或者小分子，给你造成了断键的错觉。
3. 你错误地把build工具里建立用于shear的构型这个功能，当成了做了一次剪切。实际上剪切应该单独在Forcite里做。建模时，可能会由于区域判断问题，导致建出来的模型对于聚合物链的索引不完全，造成链断裂，这看上去是一种“断键”。
我对这一点给的建议是，右键，寻找主链所有原子，复制，复制到新的xsd文件里，然后再build。
------------------------------------------------------------------------------------------------------------------------------------------------------
此外你这套工作流后边的部分和公开的MS shear计算教程相同，但是前面交联到去构象这块感觉是邪修。至少对于图示的盒子做NPT平衡我认为是毫无道理的。此外聚合物松弛的标准EQ21工作流也没体现在你的perl script里。

作者
Author: 15939474177 时间: 2026-3-10 20:58

yuzc 发表于 2026-3-10 15:39
shear的原理在于空间切片+动量交换，这里不涉及键的事情。
所以我这里有几种可能供你参考：
1. 你不小心 ...

首先谢谢老师能在百忙之中答复学生，万分感激
1、确实没点击了calculate bonds
2、交联度都在80以上，建立两层界面做吸附的时候一开始没采用加真空层的构象时，确实出现了断键，加了以后再去构建两层界面时就没再出现聚合物原子断键了，所以我理所应当想着建立三层用这个模型应该也没断键吧，结果确实如上出现了断键、
3、build建模以后(也是出现断键聚合物)，想试试我把这个模型拿去forcite进行限制性剪切了，如上图，但是目前是失败了，可以指正一下吗？您的意思是把未交联的分子都删除吗？
最后，做NPT是因为要平衡交联以后的密度，您说的聚合物松弛的标准EQ21工作流是什么意思呢？学生比较愚钝，刚接触MS不久，因为之前一直在做GROMACS
4、再次感谢老师

作者
Author: yuzc 时间: 2026-3-11 16:59

15939474177 发表于 2026-3-10 20:58
首先谢谢老师能在百忙之中答复学生，万分感激
1、确实没点击了calculate bonds
2、交联度都在80以上， ...

没贴报错不能帮你分析原因。
我又仔细看了一下你的问题，建模这块的问题我认为可以关闭了。你的回复也没答到点子上，这个靠你后期自己排查就好了。
此外跑剪切这块的问题描述非常不明确，参数，设置文件，后处理，报错，z方向上动量统计信息都没贴。
最后就是聚合物松弛的标准工作流，如果你不清楚的话，可以看下相关的文献。可能你是接触高分子模拟不久，或者以前不是高分子物理背景出身的。
高分子松弛的参考文章：
https://link.springer.com/article/10.1007/s00214-013-1334-z 只看支持信息里对应的lammps.in

作者
Author: 15939474177 时间: 2026-3-11 21:13

yuzc 发表于 2026-3-11 16:59
没贴报错不能帮你分析原因。
我又仔细看了一下你的问题，建模这块的问题我认为可以关闭了。你的回复也没 ...

谢谢老师指证，建模这块我又摸索了一遍，已经建出来没有断键的模型了，是我这边自己理解建模方法有问题，应该拿交联后出来的模型，直接做限制性建模
剪切模块就是forcite confined sheer 剪切速度是0.05，时间500ps 温度353K，并没有报错信息，剪切结果就是上图那样，Z方向的压力分布以及聚合物相对原子浓度分析我贴上了老师，
这块我还在尝试新建立的无断键模型的剪切，先做了几何优化，再进行ANNEAL NVT系综（五次循环，300-500k，600ps退火），然后再进行forcite confined sheer剪切

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