Treating electron and ion kinetics on the same physics footing or in a symmetric framework, dispersion relations (ωr, k) for ion and electron modes in a 1D Vlasov-Poisson (VP) plasma were obtained in the limit of weak flattening of the electron and ion distributions. Using this information of (ωr, k), a nonlinear VP solver with a low amplitude, harmonic external forcing has been solved. A possibility of simultaneous excitation of all the electrostatic normal modes was demonstrated in Trivedi et al. [Phys. Plasmas 25, 112102 (2018)]. The main focus of the paper is on the excitation of normal modes by applying a small amplitude external, monochromatic, electric field drive. In order to drive as well as to identify Phase Space Vortices (i.e., Bernstein-Greene-Kruskal/Cnoidal modes, etc.) formed out of the driven nonlinear VP system, assuming local flattening of electron and ion distributions, the real part of the plasma dispersion function for arbitrary ratios of ion to electron masses and temperatures is solved, which is known to be accurate in the low amplitude, harmonic limit, as pointed out by Schamel [Phys. Plasmas 19, 020501 (2012)].