各位师兄师姐好，我想用Amber做一个含铜蛋白的分子动力学模拟，这个蛋白是双铜中心，每个铜配位3个His上的N，就是请问铜的力场参数怎么确定？

http://ambermd.org/tutorials/advanced/tutorial1_orig/index.htm

我是想知道frcmod力场文件是怎么算的，就是铜和配位残基体系的参数。

应该是可以用Gaussian算的，不过读原子半径好像有点问题，也不知道是怎么解决的

楼上，请教一下，高斯难道可以动力学模拟？

高斯是算铜位点的力场参数的，然后在amber里建力场。

查了下，先优化，不算电荷，优化结果再单点能，可以读入原子半径了。。。
不过我的老变形，最后计算resp时老是出错， 估计这会是搞不定了，，，

好的，我研究下，谢谢