umxGxE with quadratic mean effects

Is it possible to modify the umxGxE function to include quadratic mean effects in the case where the twins do not share the moderator? If not, can anyone refer me to an R script that estimates Purcell’s original ACE-XYZ-M model with quadratic mean effects?

Best regards,
Emilie