antechamber sqm Error

报错内容

antechamber -i XXXXXX.mol2  -fi mol2 -o ligand.mol2 -fo mol2 \
	-c bcc -s 2 -at gaff2 -nc 0 -dr n -j 4 -rn LIG -pf y

Info: Total number of electrons: 101; net charge: 0
Info: The number of electrons is odd (101).

错误原因分析

分子的电子总数是奇数，而设置的净电荷为 0(中性的分子应该拥有偶数个电子)，导致电子总数与净电荷和自旋多重性设置不匹配。

确认分子总电荷和自旋多重性

1. 使用pymol确认分子总电荷

# 加载分子文件
load XXXXXX.mol2

# 初始化变量
total_charge = 0.0

# 迭代所有原子并累加电荷
iterate all, total_charge += partial_charge

# 打印总电荷
print("Total charge:", total_charge)

如果分子是中性的(净电荷为 0)，但电子总数是奇数，说明可能需要调整净电荷。例如：

• 如果分子应该是带一个正电荷(阳离子)，设置 "-nc 1"。
• 如果分子应该是带一个负电荷(阴离子)，设置" -nc -1"。

2. 确认自旋多重性

自旋多重性 -m 的设置应该与电子总数一致。通常，对于奇数个电子，自旋多重性应为 2(即单电子)；对于偶数个电子，自旋多重性应为 1(即成对电子)。

3. 解决方法

# 尝试设置正电荷
antechamber -i XXXXXXXXX.mol2 -fi mol2 -o ligand.mol2 -fo mol2 -c bcc -s 2 -at gaff2 -nc 1 -m 2 -dr n -j 4 -rn LIG -pf y
# 或 尝试设置负电荷
antechamber -i XXXXXXXXX.mol2 -fi mol2 -o ligand.mol2 -fo mol2 -c bcc -s 2 -at gaff2 -nc -1 -m 2 -dr n -j 4 -rn LIG -pf y