MS计算MOFs晶体对于C2H2的吸附点位时碰到的一些问题

YTBW

大佬们好，本人计算小白一个，最近在用自己合成的一个MOFs晶体的单晶数据尝试进行计算想确定一下这个晶体对于C2H2气体分子的吸附点位，碰到了一些问题。
1.我用DMol3模块对晶体结构进行优化的时候一直失败，查了一下结果文件显示原因是我使用的QEq电荷由于“This periodic system bears a net charge (40.000000e) and the
QEq calculation cannot be performed”这个问题导致无法使用QEq电荷来进行结构优化，导致后续无法查看原子电荷情况。我想问一下如何解决这种情况？
2.使用sorption模块中的fix pressure进行模拟时如果使用QEq电荷则无法计算，如果选择use current则因为之前的结构优化失败（从CIF导入的结构中原子电荷都显示为0）默认结构中所有原子的电荷为0，感觉算出来的结果不可靠。
3.本来想问模拟出最低能量结构之后如何以此结构出发进行单点能计算得出对应的主客体结合能，现在前两步无法进行也就做不下去了
补充一下，我有尝试使用forcite进行结构优化，但是同样会因为QEq无法使用从而导致计算失败，现在就有点迷糊了