计算非中性离子的福井函数求助

weiyiran

老师好！想计算某种染料（比如罗丹明B，CAS81-88-9）的福井函数来预测和自由基的反应位点，因为他的结构中是一个有机阳离子和一个阴离子，有以下两个问题不理解：
1）计算ｐN（体系中N个电子的电子密度）采用的分子结构，应该是阳离子，还是包含阳离子和阴离子呢？查找手册和文献，计算的都是中性分子，没有看到计算离子的。
2）计算ｐN+1，即电荷+1，是整个分子设置电荷加1，还是直接去掉阴离子呢？电荷-1也是一样的问题。从原理上直接去掉离子感觉是不对的，其实和第一个问题是密切联系的
辛苦老师指点一下！

wzkchem5

1）这取决于你是否认为那个阴离子是罗丹明B的一部分。换句话说，在你们的小领域内，如果你发现罗丹明B的阴离子反应了，你写文章能否说罗丹明B反应了？如果在你们小领域内确实能这么说，就加阴离子，否则不加。
2）是否保留阴离子，必须和pN体系一致，只改变电荷和自旋多重度。Fukui函数是在原子核的数目、种类以及坐标不变的情况下，只改变电子数所导致的电子密度变化，前提条件很重要。并不是用任意手段改变电荷都行

weiyiran

感谢您的耐心解答！这个阴离子就是一个氯离子，染料降解一般是指有机部分的反应，对罗丹明B而言就是阳离子。如果不加阴离子，只计算阳离子部分，那么他的电荷我是应该定为0还是+1呢？按照福井函数的含义，似乎必须是中性分子，如果只计算阳离子部分，写成电荷为0不符合实际情况啊。但如果阴阳离子都填上，因为阴阳离子之间是离子健，优化会比较困难。所以有些两难，还请老师再指点一下。