  The Q-state Potts model can be extended to noninteger and even complex Q in
the FK representation. In the FK representation the partition function,Z(Q,a),
is a polynomial in Q and v=a-1(a=e^-T) and the coefficients of this
polynomial,Phi(b,c), are the number of graphs on the lattice consisting of b
bonds and c connected clusters. We introduce the random-cluster transfer matrix
to compute Phi exactly on finite square lattices. Given the FK representation
of the partition function we begin by studying the critical Potts model
Z_{CP}=Z(Q,a_c), where a_c=1+sqrt{Q}. We find a set of zeros in the complex
w=sqrt{Q} plane that map to the Beraha numbers for real positive Q. We also
identify tilde{Q}_c(L), the value of Q for a lattice of width L above which the
locus of zeros in the complex p=v/sqrt{Q} plane lies on the unit circle. We
find that 1/tilde{Q}_c->0 as 1/L->0. We then study zeros of the AF Potts model
in the complex Q plane and determine Q_c(a), the largest value of Q for a fixed
value of a below which there is AF order. We find excellent agreement with
Q_c=(1-a)(a+3). We also investigate the locus of zeros of the FM Potts model in
the complex Q plane and confirm that Q_c=(a-1)^2. We show that the edge
singularity in the complex Q plane approaches Q_c as Q_c(L)~Q_c+AL^-y_q, and
determine the scaling exponent y_q. Finally, by finite size scaling of the
Fisher zeros near the AF critical point we determine the thermal exponent y_t
as a function of Q in the range 2<Q<3. We find that y_t is a smooth function of
Q and is well fit by y_t=(1+Au+Bu^2)/(C+Du) where u=u(Q). For Q=3 we find
y_t~0.6; however if we include lattices up to L=12 we find y_t~0.50.
