‘‘‘ ppGppSys1 = @reaction_network begin kin, ∅ --> G6P hill(G6P,v01,k01,1), G6P → F6P end odesys = convert(ODESystem, ppGppSys1; combinatoric_ratelaws=false, remove_conserved = true) latexify(odesys)‘‘‘ produces \begin{align} \frac{\mathrm{d} \mathtt{G6P}\left( t \ri…

I realise that the issue is not with “latexify” but rather the definition of Catalyst.mm and Catalyst.hill. If I define a “model” with few reactions using these rates and then use the following code to check the ode rate laws: rxs = reactions(lmodel) ν = oderatelaw.(rxs) I can see that the “mm” …

Can you please provide a full MWE that is properly formatted as code in your post so it can be copied and pasted to run?

Thanks - I figures another way, in which “mm” behaves as expected. Here is the code to illustrate. The first reaction (M0 to M1) is written, as I understand the “mm” usage from example in the docs and gives be the squared dependency on M0 concentration: mmModel = @reaction_network begin …

It seems like maybe there is a misunderstanding about how Catalyst arrows work. When you use --> Catalyst always generates a mass action term for the rate law defined by the substrate stoichiometry that multiplies whatever your “rate” term is. i.e. all the terms that multiply k here . In your case k …

So just to add a concrete example, mm(M0,v01,k01), M0 + 2A --> M1 should give the rate law of mm(M0,v01,k01) * M0 * (A^2/2) per that page I linked. If you only want the mm term you can use mm(M0,v01,k01), M0 + 2A => M1 which should just give a rate law of mm(M0,v01,k01)

Latexify/Catalyst introduces extra exponent - issue with mm and hill?

New to Julia

isaacsas August 25, 2025, 12:50pm 4

See PSA: how to quote code with backticks for info on formatting Discourse posts.

Topic		Replies	Views
Reaction rate laws of chemical reaction networks Modelling & Simulations question , diffeq	7	1233	March 3, 2021
Implement Matlab/Simbiology concepts (e.g. compartments, repeated assignments) in Julia/Catalyst General Usage question , modelingtoolkit , catalyst	9	403	May 15, 2024
Combining Catalyst DSL models and ModelingToolkit explicit ODE models Chemistry	18	1385	November 28, 2023
Programmatic Generation of Reaction (Catalyst.jl) with User-Defined Reaction Rate New to Julia	6	691	May 31, 2022
Anybody know why this error occurs when trying to define a Catalyst @reaction_network in Jupyter Notebooks? Modelling & Simulations catalyst	4	505	January 19, 2023

Latexify/Catalyst introduces extra exponent - issue with mm and hill?

Related topics