Reporting the results of a continuous moderating variable on a categorical independent variable I'm working with a dataset that was collected to test the following [adapted] research question: Does listening to music (categorical x; 2 conditions: pre_music - post_music) affect test scores (y; continuous variable), and is that relationship modified by IQ level (z; continuous variable)?
A reprex dataset:
dat <- data.table(participant_id = rep(1:30, each = 2), music_condition = rep(c("pre", "post"), 30), test_score = rnorm(60, 50, 15), iq = rep(rnorm(30, 100, 15), each = 2))

Note that iq was measured only once at the beginning of the experiment and is therefore repeated twice per each participant because the value does not change between conditions.
My first instinct was to make a model, m1 wherein the effects of music_condition on test_score were tested (or, equivalently, get the p-value from a paired t-test):
lmer(test_score ~ music_condition + (1 | participant_id), data = dat) %>% 
summary()

No significant effect.
Then, to answer the second question about IQ moderating the aforementioned relation, my thought was to fit a second model, m2, with music_condition*iq interaction:
lmer(test_score ~ music_condition*iq + (1 | participant_id), data = dat) %>% 
summary()

No moderating effect of IQ.
My question: is it preferred to report the p-value from m1 to address the first research question, and the p-value from m2 to address the second research question? Or should I use both p-values (music_conditionpre and music_conditionpre:iq) from m2? Also, in the context of my research question, what does the iq fixed effects represent in m2 and shouldn’t it be exactly 1 because it does not change at all between conditions?

Since you have paired pre-music & post-music test scores for each participant in the study, you can analyze the data by regressing the post-music test scores (the outcome) on the pre-music scores, the IQ scores and their interaction (the predictors):
$$
\begin{aligned}
\operatorname{E}\left\{ S_{\text{post}} \right\} = \beta_0 + \beta_1 S_{\text{pre}} + \beta_2 \text{IQ} + \beta_3 S_{\text{pre}} \times \text{IQ}
\end{aligned}
$$
where $S_{\text{pre}}$ and $S_{\text{post}}$ are the test score before and after music.
Aside: If some participants received the "treatment" (listened to music) and others didn't, we would include a treatment indicator variable. That experimental design would have been more convincing as you would have been able to compare the treated with the controls. However, it's not the design you chose, so we move forward.
How do we re-state the research question in terms of this model for $S_{\text{post}}$?

*

*Does listening to music affect (increase) test scores while adjusting for IQ: Is the pre-music score coefficient $\beta_1$ significantly different from (bigger than) 1?

*Does IQ moderate the relationship between pre- and post-music scores: Is the interaction term $\beta_3$ significant?

This analysis requires pivoting the data, so that the two scores for each participant are on the same row. There is no need for random participant effects. In R formula notation:
lm(score_post ~ score_pre * iq)

