This paper introduces a novel parallel shift operator finite-difference time-domain (SO-FDTD) method for plasma in the dispersive media. We calculate the interaction between the electromagnetic wave of various frequencies and non-magnetized plasma by using the parallel SO-FDTD method. Then, we compare the results，which are calculated with serial and parallel SO-FDTD executions to obtain the speedup ratio and validate the parallel execution. We conclude that the parallel SO method has almost the same precision as the serial SO method, while the parallel approach expands the scope of memory and reduces the CPU time.