The reaction kinetic parameters were fitted with a single-and two-step reaction network (Figure ). The developed Julia software is able to fit reaction orders, the Arrhenius pre-exponential factor (A) and activation energy (Ea). For better comparison between each amidation reaction a two-step reaction network was used and the reaction orders were fixed to 1 for ester, amine, TBD and intermediate. In first reaction the ester is activated by TBD to form an active N-acyl intermediate. This intermediate then reacts with the amine to corresponding the amide product, regenerating TBD. In total 10 different amidation reactions were investigated, including a comparison of different benzoic esters: methyl vs ethyl vs isopropyl ester (Figure ). The activation energy (Ea1) and pre-exponential factor (A1) for the rate determining step (activation of ester) is displayed in Figure .