where L is the length of the system, but I am not sure if this is correct.
Will the order of the coefficients change depending on the region of the phase diagram that I am studying?

They will only depend on the charge conservation you're using.

I guess you're aware of this, but just to be clear: when you do this, you write down all 2^L states, so you could just as well use exact diagonalization....