In many astrophysical applications, the cost of solving a chemical network represented by system ordinary differential equations (ODEs) grows significantly with size network, and can often represent significant computational bottleneck, particularly in coupled chemo-dynamical models. Although standard numerical techniques complex solutions tailored to thermochemistry somewhat reduce cost, more ...