I’ve been trying to implement a model with multiple generations of a real scalar “phi”, which leads to a parameter with four indices if the quartic self-interaction
It's likely that this is not supported. I haven't ever worked with a model having that feature or have one at hand to test it.
The simple solution is to define independent fields. The only drawback is that the Lagrangian will be a bit lengthy.