We consider Quantum Electrodynamics in $2{+}1$ dimensions with $N_f$ fermionic or bosonic flavors, allowing for interactions that respect the global symmetry $U(N_f/2)^2$. There are four bosonic and four fermionic fixed points, which we analyze using
the large $N_f$ expansion. We systematically compute, at order $O(1/N_f)$, the scaling dimensions of quadratic and quartic mesonic operators. We also consider Quantum Electrodynamics with minimal supersymmetry. In this case the large $N_f$ scaling dimensions, extrapolated at $N_f{=}2$, agree quite well with the scaling dimensions of a dual supersymmetric Gross-Neveu-Yukawa model. This provides a quantitative check of the conjectured duality.
We consider Quantum Electrodynamics with an even number $N_f$ of bosonic or fermionic flavors, allowing for interactions respecting at least $U(N_f/2)^2$ global symmetry. Both in the bosonic and in the fermionic case, we find four interacting fixed p
oints: two with $U(N_f/2)^2$ symmetry, two with $U(N_f)$ symmetry. Large $N_f$ arguments suggest that, lowering $N_f$, all these fixed points merge pairwise and become complex CFTs. In the bosonic QEDs the merging happens around $N_fsim 9{-}11$ and does not break the global symmetry. In the fermionic QEDs the merging happens around $N_fsim3{-}7$ and breaks $U(N_f)$ to $U(N_f/2)^2$. When $N_f=2$, we show that all four bosonic fixed points are one-to-one dual to the fermionic fixed points. The merging pattern suggested at large $N_f$ is consistent with the four $N_f=2$ boson $lra$ fermion dualities, providing support to the validity of the scenario.
