根据纯物质状态的质量密度、分子质量，换算成加入的有机分子和水分子数

绝对不要用constraints         = all-bonds，公社论坛之前我回复过很多次

dt                  =  0.0002莫名其妙，用这么小步长干嘛，一般都用0.002

若无特殊情况，rlist也设成1.0

体系都还没平衡，甭用PR压浴

甭用gen_vel             =  yes。老老实实从0K线性升温

periodic_molecules  = yes根本就不是用在当前场合的


你的mdp问题太多，不懂的不要瞎写