In this work, we show a proof of concept that optimization of plasma physics simulations using gradient-based methods can work. We also show that by applying the correct loss function, one can tease out previously unknown physics. This is an example of AI-based physics discovery.