Treating solar-wind ions as having a finite temperature instead of assuming they are cold changed the modeled electrical environment above the sunlit Moon, according to a new preprint. The warm-ion calculations produced a less strongly negative potential, a thinner sheath and a weaker overall electric field in the classical Debye-type case. The largest changes in the model’s more complex, non-monotonic potential structure appeared at higher lunar latitudes.
The study follows a near-surface layer containing photoelectrons, solar-wind electrons and ions. The authors suggest that changes in this structure could affect how dust becomes charged and moves above the lunar surface, although dust charging and motion were not calculated directly in the study.
A cold-ion shortcut under scrutiny
The study asked how ion temperature alters the lunar photoelectron sheath, particularly when the usual cold-ion approximation is replaced by a finite-temperature description. Its numerical model added a thermal-pressure term to the equation governing ion flow and solved Poisson’s equation—the equation linking charge in the plasma to electric potential—to calculate the sheath’s potential, electric field, thickness and particle-density profiles.
This was a deterministic modeling study, not an analysis of a participant or spacecraft sample. The calculations compared a cold-ion case, with Ti = 0, against warm-ion cases in which the ion temperature was finite. Under the stated nominal conditions, the model used an ion density of 5 cm−3, an electron temperature of 15 eV and an ion speed of 400 km s−1 at the sheath edge, while the ratio of ion temperature to electron temperature, Ti/Te, was varied from 0 to 10.
For the non-monotonic sheath solutions, the analysis required ions to reach the sheath edge at least at a specified modified Bohm speed, a threshold intended to allow that structure to persist. The calculations also varied latitude and the adiabatic index.
The modeled changes were strongest away from the lunar equator
In the classical Debye-type solutions, the finite-temperature treatment shifted several basic sheath properties together. Relative to cold ions, it gave a negative potential of smaller magnitude, reduced the sheath thickness and weakened the overall electric field. These are model outputs rather than direct measurements of conditions on the Moon.
The latitude dependence was clearest in the non-monotonic solutions. At higher latitudes, finite ion temperature was associated with a higher modeled surface potential and a less deep negative potential minimum. In other words, the shape of the potential profile changed most noticeably in those higher-latitude cases, rather than producing one uniform shift across all modeled locations.
The model also placed key features of the sheath at different heights as its parameters changed. Increasing Ti/Te and the adiabatic index reduced the calculated sheath thickness and shifted the locations of both the potential minimum and the point where the electric field became zero, with those shifts occurring on meter scales.
The particle populations changed alongside the potential. Warm-ion cases had higher modeled ion densities, while their weaker potential minima were linked in the calculation with less photoelectron trapping and less reflection of solar-wind electrons. Those results describe the behavior of the modeled populations within the assumed sheath solution.
One threshold highlighted by the analysis was the critical potential Vc, the value associated with obtaining real ion-density solutions in the model. Under the nominal electron temperature of 15 eV, Vc fell from 7.5 V for cold ions to 0.4 V for the case with ψ = 1 and 0.1 V for ψ = 3 when Ti/Te reached 10. The supplied analysis reports these as parameter-specific numerical results, without statistical comparisons or uncertainty intervals.
What the calculations can—and cannot—say
The authors interpret the combined changes in potential, field strength, thickness and particle populations as a material reshaping of the lunar photoelectron sheath, especially at high latitudes. They say the effect could be important for dust charging and dust dynamics above the lunar surface. They also connect the result with the detection of non-monotonic potentials.
The result does not show that warm ions produce the same changes in the actual lunar environment. The study did not directly measure the sheath, did not directly simulate dust charging or dust motion, and did not establish that its findings apply beyond the plasma and photoemission conditions used in the calculations.
Several assumptions narrow the result’s reach. The photoelectron and solar-wind electron populations were represented by Maxwellian distributions, and the sheath was treated with one-dimensional steady-state equations and zero net current. The calculations therefore do not test non-Maxwellian distributions, time-dependent behavior or multidimensional structure.
The paper calls for direct comparisons between the predicted high-latitude profiles and THEMIS or ARTEMIS measurements. It also identifies further modeling of dust charging and dynamics, along with tests using broader plasma descriptions, as steps needed to determine how robust the modeled pattern is.
A result awaiting observational checks
The work is an arXiv version-one preprint dated 20 August 2026. Its conclusions are conditional on the stated numerical setup rather than presented as a completed observational finding.
The work was supported by the Department of Space, Government of India, while NASA contract NAS5-02099 was acknowledged for use of THEMIS data. The paper reports that THEMIS and FISM2 data are available through the listed repositories and that other supporting data are contained within the article.
Paper data and sources
Original title: Effect of ion temperature on lunar photoelectron sheath
Authors: Trinesh Sana, Sanjay K. Mishra
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text