A Wigner-Jordan (WJ) transformation is used to study the one-dimensional (1D) and two-dimensional (2D) quantum antiferromagnetic Heisenberg model. The advantage of using the Wigner-Jordan transformation is that it preserves all spin-commutation relations as well as the spin on-site exclusion principle. In the ID case a nearest-neighbor covalent-bonding state of the WJ fermions is found to have a ground-state energy (-0.4351J per site) comparable with that from the Bethe-ansatz solution (-0.4431J per site), and a linear energy spectrum at low energies with velocity 1.6366J, in close agreement with the velocity obtained by Haldane for the quantum antiferromagnetic Heisenberg model with 1/d2 interaction (1.5708J). The method used for studying the ID model is then applied to the 2D Heisenberg model in a square lattice. The resulting state at finite temperature is the in-phase flux state, i.e., a flux state of the Wigner-Jordan spinless fermions with an in-phase fermion orbital current circulating around each elementary plaquette. The single-particle excitation spectrum, i.e., the energy dispersion of reversing orientation of a spin in the system, of the in-phase flux state is overall similar to that of spin waves, with significant difference near the edge point k=(pi,0), at which the excitation energy of the WJ fermion is zero, whereas that of the spin-wave excitation is 2J. The specific heat of the in-phase flux state predicts a correct temperature dependence over the entire temperature range, namely a T2 dependence at low temperature, a peak near T/J=0.6, and a 1/T2 decreasing at high temperature. It also gives excellent agreement with the specific heat calculated from numerical methods. In contrast, the spin-wave theory only correctly predicts a T2 dependence at low temperature. The Raman spectrum of the in-phase flux state is calculated and shows significant improvement over that from spin-wave theories compared with the experimental spectrum of La2CuO4. The exchange parameter, J, obtained from the comparison is 1060 cm-1, in agreement with that obtained from analyzing the neutron-scattering data.