We simulate four flavor noncompact lattice QED using the Hybrid Monte Carlo
algorithm on 104 and 164 lattices. Measurements of the monopole
susceptibility and the percolation order parameter indicate a transition at
β=1/e2=.205(5) with critical behavior in the universality class of
four dimensional percolation. We present accurate chiral condensate
measurements and monitor finite size effects carefully. The chiral condensate
data supports the existence of a power-law transition at β=.205 in the
same universality class as the chiral transition in the two flavor model. The
resulting equation of state predicts the mass ratio mπ2/mσ2 in
good agreement with spectrum calculations while the hypothesis of a
logarithmically improved mean field theory fails qualitatively.Comment: 17 pages, 10 figure