Hi! I ran different models using the glmer function in the lme4 package and compared the performance of these models using the compare_performance function in the performance package. The model that best fit the data was:
I discuss some recommended best practices from guides on mixed models as well as my own suggestions in this answer on Cross Validated, with references that may be useful for you: