We present methods for the creation of semi-empirical quantum approaches and reactive force fields through force matching to quantum simulation data for materials under reactive conditions. Our methodologies overcome the extreme computational cost of standard Kohn–Sham Density Functional Theory (DFT) by mapping DFT computed simulation data onto functional forms with linear dependence on their parameters. This allows for quick parameterization of our models by avoiding the nonlinear fitting bottlenecks associated with most molecular dynamics model development. We illustrate our approach with two different systems: (i) determination of density functional tight binding models for aqueous glycine dimerization, and (ii) determination of the Chebyshev Interactional Model for Efficient Simulation (ChIMES) reactive force field for metallic liquid carbon. In each case, we observe that our approach is easy to parametrize and yields a model that is orders of magnitude faster than DFT while largely retaining its accuracy. Overall, our methods have potential use for studying complex long time and length scale chemical reactivity at extreme conditions, where there is a significant need for computationally efficient atomistic simulations methods to aid in the interpretation and design of experiments.