求助：gromacs能不能实现Drude振子模型中的SCF方法

木贞

使用Drude振子模型进行模拟主要有两种方法。第一个是SCF方法，即一旦核构型发生变化，电子自由度就会立即放松到它们的基态。在采用SCF方法的MD模拟中，在每个积分步长都通过能量最小化来求解Drude振子的位置，使得Drude振子到达能量最低处。目前采用更多的是第二种方法，扩展拉格朗日法，即赋予Drude振子一定质量，并使用双恒温器法进行模拟，通过Gromacs可以实现。但是因为想得到每个核构型下的最低能量，所以想用SCF方法，不知道Gromacs有没有相关功能可以实现。