Estou fazendo simulações com duas moléculas distintas, mas acaba sempre caindo no mesmo problema, parece que o gaussian não equilibra as moléculas quando calculo o cation e o anion. Quando calulo a carga de ambos parece que tem preferencia por uma determinada molécula que não sei identificar qual é. Para o calculo neutro com uma e com duas moléculas os resultados são coerentes. Pensei em gerar uma superfície do potencial eletrostático mas isso não resolve meu problema.
Obs: estou desenhando as moléculas no Avogadro, não sei se isso interfere em algo, mas o mesmo gera um input para o gaussian
alguém tem alguma sugestão?