The levitation of charged dust, which may cause serious hazards to astronauts and lunar rovers, has been one of the most significant challenges in lunar exploration. Here we simulate lunar sheath potentials in different solar wind conditions and solar zenith angles (SZAs) on the lunar surface by the particle-in-cell method. The simulated potentials exhibit two types of distributions as a function of height, depending on the SZAs. For SZA ∼ 0°–70°, the nonmonotonic distribution with positive surface potential dominates in the photoelectron sheath. For SZA >∼81°, the monotonic distribution with negative surface potential is observed in the plasma sheath. With the calculated potentials and the assumption that the dust radius distribution exponentially decreases, we further investigate spatial distributions of the dust levitated above the surface. It is found that number density of the levitating lunar dust is enhanced at the terminator (SZA ∼ 81°) in the plasma sheath. In the photoelectron sheath it gradually decreases as the SZA increases from 0° to 70°. Further calculations of the potential and the derived electrostatic field suggest that the dust spatial distributions can be influenced by the bulk velocity, number density, and temperature of the solar wind. Those findings deepen our understanding of lunar surface charging and the mechanism of lunar dust levitation, which can provide technical support for lunar explorations.