根据分子动力学模拟推测TBHQ在三油酸甘油酯中的溶解度参数和均方位移

LIUXJ

求助！！！参照文章“Effect of caffeic acid esters on antioxidant activity and oxidative stability of sunflower oil: molecular simulation and experiments”，“天然酚酸抗氧化剂在酯类润滑剂中的抗氧化作用的分子模拟与实验研究”和“基于分子动力学的橡胶沥青相容机理研究”的方法及参数进行分子动力学模拟，但结果与实际情况相差甚远（模拟的具体参数见附件）。主要体现在，在NPT系综下进行分子动力学模拟后，①发现盒子体积增大（长度有40增加至50）；②密度降低，原AC盒子密度设置为0.915，其模拟结果却不到0.5；③其溶解度参数也与文章报道相差甚远。
推测其原因可能是前期分子处理的问题，但尚不知该如何解决。