本人是lammps菜鸟，最近在做固体推进剂粘结剂的分子动力学模拟，看了些文献，固化反应基本都是-N=C=O和-OH的反应，其中C原子与羟基的O原子相连，氢原子和N原子相连，C原子和氢原子我都能在生成的data文件里找到对应的编号，疑惑的是H原子的原子类型都一样，我该怎么做才能把-OH上的氢原子和其他碳原子上的H原子区分开呢？
另外我在文献中看到了交联过程结束后每一次循环结束都要让系统充满H原子，这一步的目的是什么呢？难道是把所有的H原子都设为一个组，然后让N、H成键吗？那这样会不会把其他C原子上的氢原子拿过来成键了呢？下面是文献的一部分：原文见Molecular dynamics simulation of bonding role of cyclic borate ester in improving mechanical properties of the HTPB propellant - ScienceDirect
The crosslinking flowchart is shown in Fig. 3. The cutoff distance increased gradually from 4 to 12 Å, with a 1 Å increment for each crosslinking loop. After each loop, the system was saturated with hydrogen atoms.
