Simon Boothroyd shared a script to analyze how the molecule populations distribute among (or move) the gradients for each of the parameters based on the first iteration, to find out whether we see multimodal distributions indicating a split of parameters.
Script is:
Here are the plots for fit7 for the TIG parameters, k1_bondorder1 and k1_bondorder2 for each parameter
TIG0:
TIG1c
TIG1d
TIG2
TIG3
TIG4
TIG5a
TIG5b
TIG6
TIG7
TIG8
TIG5a, TIG6, TIG7, TIG8 have gradients that vary on a large scale here.