I'm trying to simulate the sulfate anion structure in Gaussian, but my double bonds disappeared after the calculations. I had the same problem when I was trying to simulate the acetate anion.
For the sulfate, I put two double bounds in my structure and I decreased their length compared to the two remaining single length.
The following are my file details.