Hi everyone,
I am trying to set up a reaction network using Catalyst where one species remains at a constant value. I’m modelling a photo reaction and want to track the concentration of photons.
In a simplified form, it looks like this, where Ph (Photon) should be handled as a constant species:
rn = @reaction_network begin
k1, A + Ph --> B
k2, C + Ph --> D
end
The idea is that k2 >> k1, so only when C is used up will there be photons left for A to react. Essentially, C acts as a quencher for the reaction of A.
I’m not sure if this is the best way to approach it or if it’s even possible using Catalyst, but I was hoping someone could point me in the right direction.
Thank you in advance!