I am using Quantum Espresso to simulate a platinum surface with ligands, so far everything is working ok. I now need to introduce a solvated hydronium $H_3O^+(H_2O)_2$. Is there something I should be careful about when introducing charged species? I will add the keyword tot_charge = +1 but I see ...