Thanks for the explanation Gianluca! Plotting it out, it seems like the solution is in the 4th quadrant. In my model, ix and iy can't be negative numbers, I modified the code you had to include the constraints for ix and iy, and the plot turned out empty. I'm not sure if this is because I didn't put the constraint in the right area of the code, or does a positive solution just not exist for ix and iy? (and how would I be able to tell for the latter?) Thank you so much!
f[iy_] = ix /.
First@Solve[
Simplify[(-r +
9/40 b ((7 Sqrt[c] Sqrt[r] +
Sqrt[q + b iy q + 49 c r])^2/(1 +
b ix)^2 + (2 q^2)/(Sqrt[
q + b ix q + 49 c r] (-14 Sqrt[c] Sqrt[r] +
Sqrt[q + b ix q + 49 c r] +
Sqrt[q + b iy q + 49 c r]))) == 0 && ix >= 0 &&
iy >= 0) /. {q -> 50, r -> 50, c -> 1/100, b -> 1/10}], ix,
Reals];
Plot[f[iy], {iy, -15, 20}]