关于高斯优化有两个问题：
1.给定一个分子，想固定其中某个二面角，优化分子的其它部分，应该用什么命令实现？
2.进一步的一个问题：假如想扫描某个二面角的势能面，其它部分优化，就是scan一个relaxed PES，应该如何实现？
目前我只会硬扫二面角势能面，举例如下：
# B3LYP/3-21G opt=modredundant
molecule
0 1
A x,y,z
B x,y,z
C x,y,z
D x,y,z
D 1 2 3 4 += 0.000000 S 7 30.000000
这样就能实现分子中二面角ABCD 从0到180度，步长为30度的势能面。除了二面角，其余部分不优化。
现在求软扫的方法？请高人指点！ |
回复此楼
