Two-dimensional particle-in-cell simulations for the interaction of a fs pulse with a low-density plasma were discussed. The temporal evolution and spatial distribution of terahertz electron current developed in a low-density plasma layer were also presented. It was found that laser and plasma parameter such as laser intensity, background plasma density, and pulse width were of key importance to the process. It was also observed that the radiation peaked at at 6.4 THz, with 900 fs duration and 9% bandwidth, can be generated in a plasma of density 5 × 1017 cm-3.