Thank you so much, Steve.
I know that some atoms are too close in the boundaries and it is ok to get rid of them.
I overlooked the sequence of thermo_style and thermo_modify. The program is now running.
Thank you and all others again for your valuable time!
Yujie