Simulating phase transition (structure change) with temperature

Thanks Steve and Carlos !

Just as you said, the "iso" "couple xyz" so they have to change proportionally. I got the correct answer after I changed to "aniso", which "couple none" so x, y and z are independent.

Cheers,
Xing