Statistical mechanics of a 1D multivalent Coulomb gas can be mapped onto non-Hermitian quantum mechanics. We use this example to develop the instanton calculus on Riemann surfaces. Borrowing from the formalism developed in the context of the Seiberg-Witten duality, we treat momentum and coordinate as complex variables. Constant-energy manifolds are given by Riemann surfaces of genus g ≥ 1. The actions along principal cycles on these surfaces obey the ordinary differential equation in the moduli space of the Riemann surface known as the Picard-Fuchs equation. We derive and solve the Picard-Fuchs equations for Coulomb gases of various charge content. Analysis of monodromies of these solutions around their singular points yields semiclassical spectra as well as instanton effects such as the Bloch bandwidth. Both are shown to be in perfect agreement with numerical simulations.