We study the protocol of entanglement harvesting when two local probes couple to the vacuum of a real scalar quantum field with arbitrary temporal profiles. We use a Hermite expansion to efficiently compute smeared field propagators in closed-form, recasting the negativity between the probes as a matrix product. We then optimize the protocol under different signalling conditions, enhancing entanglement harvesting by several orders of magnitude. This optimization would take current experimental proposals beyond the regime of second order perturbation theory.