arXiv:2511.00006v2 Announce Type: replace-cross Abstract: Based on the multidimensional Leibniz integral rule, we propose a novel unbiased stochastic derivative estimation framework for sample performance functions that are discontinuous in the parameter of interest. In this framework, for general discontinuities, a change of variables followed by the Leibniz integral rule generally yields a multi-run estimator, whereas for indicator-type discontinuities, a suitable transformation of the sample performance leads directly to a single-run estimator. The new framework includes existing likelihood ratio (LR)-based approaches as special cases and applies more broadly. In particular, it encompasses generalized LR (GLR) methods and the new formulation permits us to obtain much milder and simpler sufficient regularity conditions for these estimators to be unbiased, compared with the conditions given in previous papers. We also provide guidance on how to construct efficient single-run estimators in settings where existing methods may require the simulation of multiple sample paths or may be inapplicable. We validate the effectiveness and robustness of the proposed framework through an inventory-price control problem.