Authors’ abstract: For a given polynomial \(P\) with simple zeros, and a given semiclassical weight \(w\), we present a construction that yields a linear second-order differential equation (ODE), and in consequence, an electrostatic model for zeros of \(P\). The coefficients of this ODE are written in terms of a dual polynomial that we call the electrostatic partner of \(P\). This construction is absolutely general and can be carried out for any polynomial with simple zeros and any semiclassical weight on the complex plane. An additional assumption of quasi-orthogonality of \(P\) with respect to \(w\) allows us to give more precise bounds on the degree of the electrostatic partner. In the case of orthogonal and quasi-orthogonal polynomials, we recover some of the known results and generalize others. Additionally, for the Hermite-Padé or multiple orthogonal polynomials of type II, this approach yields a system of linear second-order differential equations, from which we derive an electrostatic interpretation of their zeros in terms of a vector equilibrium. More detailed results are obtained in the special cases of Angelesco, Nikishin, and generalized Nikishin systems. We also discuss the discrete-to-continuous transition of these models in the asymptotic regime, as the number of zeros tends to infinity, into the known vector equilibrium problems. Finally, we discuss how the system of obtained second-order ODEs yields a third-order differential equation for these polynomials, well described in the literature. We finish the paper by presenting several illustrative examples.