It seemed to me that the net reaction rate only comes into play when the stiff chemistry solver is used, but that this part must be disabled when using DEFINE_VR_RATE?

I don't think I understand what you mean by "what do the net reaction rates of other components tel".
Below are the solved equations
