Here I present the results on the study of surface properties of beta-HMX crystal utilizing molecular dynamics simulations. The surface polarity of three principal crystal surfaces, (011),(010),and(110), is investigated by measuring the water contact angles. The calculated contact angles are in excellent agreement with the values measured by experiment and show that the surface polarity of three crystal surfaces are different. The free energies and forces of detaching an Estane chain (with and without surrounding nitroplasticizer molecules) from the three principal crystal surfaces are also calculated using the umbrella sampling method. I find that the force for Estane detachment increases with the increasing HMX surface polarity. In addition, my results show that the nitroplasticizer also plays an important role in the adhesion between Estane and HMX surfaces.