My fellow scientists,

I have a question about getting coefficient of thermal expansion of a phase as function of temperature. I use the script below but the problem is even though system has only one phase (BCC in this case) at lower temperatures minimization process creates BCC_A2#1 and #2 which I don't need. Is there a way to prevent partitioning of all alloying elements?

Thanks in advance,