I investigate a class of models with scalar and pseudoscalar solutions to the g−2 anomaly for both the muon and the electron over the mass range of perturbativity (mϕ≲50 GeV), with Yukawa couplings proportional to the lepton’s mass. In particular, I investigate the constraints from BaBar, beam dump experiments, Z decay measured quantities, LEP mono-γ searches, ee⟶ττγ searches, and solar and horizontal branch (HB) star bounds. For a pseudoscalar, I find that no region in the parameter space can simultaneously provide a solution for both the electron and the muon anomalies while maintaining the required form of the couplings, and therefore, the pseudoscalar solution is disfavored. On the other hand, I find for the scalar case that there is an...