Scientific machine learning (SciML) through differentiable simulation is a method that is faster and more robust than techniques like physics-informed neural networks (PINNs), neural operators, and other techniques. However, doing differentiable simulation correctly requires a deep understanding of automatic differentiation and the numerical properties of simulation. In this talk, we dive deep into the numerical stability of derivatives of simulation processes, showing how naive applications of differentiable programming can give surprisingly incorrect results, and how one may need to modify simulations in order to perform robust automated model discovery and calibration.Originally part of the JuliaHEP 2023 Workshop at the ECAP (Erlangen Ce...