Xe gas bubble superlattice formation is observed in irradiated uranium-10 wt% molybdenum (U10Mo) fuels. However, the thermodynamic properties of the bubbles(the relationship among bubble size, equilibrium Xe concentration, and bubble pressure)and the mechanisms of bubble superlattice formation are not well known.In this work, the molecular dynamics (MD) method is used to study these properties and mechanisms. The results provide important inputs for quantitative mesoscale models of gas bubble evolution and fuel performance. In the MD simulations, the embedded-atom method (EAM) potential of U10Mo-Xe[1] is employed. Initial gas bubbles with a low Xeconcentration(underpressured) are generated in a body-centered cubic(bcc)U10Mo single crystal.Then Xeatoms aresequentiallyadded into the bubbles one by one, and the evolution of pressure and dislocation emission around the bubbles is analyzed. The relationship between pressure, equilibrium Xe concentration, and radiusofthe bubblesis established.It was found that an overpressuredgas bubble emits partial dislocations witha Burgers vector along the <111> direction and a slip plane of (11-2). Meanwhile, dislocation loop punchout was not observed.The overpressured bubble also induces an anisotropic stress field. A tensile stress was found along <110> directions around the bubble, favoring the nucleation and formation of aface-centered cubic bubble superlattice in bccU10Mo fuels.