Suppose the market demand is $P(Q) = a-bQ$ and the manufacturing cost of an item is $c_1q_1$ for firm $1$ and $c_2q_2$ for firm $2$.
In the Cournot model, we maximize the profit functions and plug in the best response function of one firm into the other's best response function.
Why can't we do this instead: If $\pi_1 = (a-q_1-q_2)q_2 - c_1q_1$ is the profit of firm $1$, why don't we plug in $q_2^* = \text{BR}_2(q_1)$ and solve for firm $1$'s BR function. That is, we differentiate $\pi_i = (a-q_i-\text{BR}_{-i}(q_i))q_i - c_iq_i$ wrt $q_i$ for both $i=1,2$.
I tried this and it didn't work. This is also how the Stackelberg model is solved. So my question is, why does plugging in the BR function before and after create a difference when in the end both firms are choosing their best response functions? I don't get the math behind how plugging in the BR later creates the "first-mover advantage" while it doesn't if we do it the usual Cournot way.
Notation: If $i=1$, BR$_{-i}(q_1)$ corresponds to BR$_2(q_1)$. Similarly for $i=2$.