Abstract The effect of the shape of the applied voltage waveform on the energy efficiency of reactive species generation is investigated in an atmospheric pressure RF microplasma jet operated in a He/O2 mixture (99.5%/0.5%) based on a one-dimensional hybrid fluid-kinetic simulation method. Using a tailored waveform synthesized from four consecutive harmonics (with a base frequency of f b = 13.56 MHz and amplitudes of ( 160 / k ) V for the k-th harmonic), it is shown that by changing the identical phases of the even harmonics in the waveform, ϕ, the generation efficiencies of three specific reactive species (helium metastables, atomic and vibrationally excited oxygen), defined as the ratio of mean density and input plasma power, attain their maxima for different values of ϕ, due to changes in the Electron Energy Probability Function (EEPF). The phase control of the EEPF and its critical role in modulating generation energy efficiencies are explained in detail. The simulation results are verified by experimental (Phase Resolved Optical Emission Spectroscopy and Two Photon Absorption Laser Induced Fluorescence) data.