Swarms of small robots are an attractive way to survey lunar caves and lava tubes. Missions such as NASA's CADRE will fly cooperating rovers for distributed mapping, but only in daylight. The 354 h lunar night drops the surface towards 95 K, and the radioisotope heaters that let single rovers survive it cannot be supplied at swarm scale. In this study we investigate whether a swarm can instead use its own members as thermal shielding. We train one policy, shared by all twenty robots, by reinforcement learning with an evolution-strategies estimator. The policy receives a single reward, once, at the end of the night: the mean temperature of the five coldest robots. Nothing in that reward refers to position, spacing, exposure or motion. Two behaviors nonetheless emerge. The swarm aggregates, and robots circulate between the sheltered core and the exposed rim, the same pair of behaviors with which emperor penguins survive the Antarctic winter. We measure the survival threshold, the smallest stored energy per robot at which every robot is still alive at sunrise. Across ten runs it is 155.6 +/- 1.3 pcm-h, where one pcm-h is 7.2 kJ of stored heat. That is 37% below a motionless huddle and within 13% of a physics-based lower bound, and at a budget where the motionless huddle loses six robots, every robot survives. The behavior needs no shared spatial reference and no robot-robot contact, and transfers from twenty robots to forty without retraining.
Figures & tables
Fig. 1: Twenty simulated robots during the 354h lunar night, from one rollout of the learned policy. Ring color is each robot’s sky exposure χ on the scale of Fig. 3 . Top: the scattered start at hour 0. Every robot is red: with no neighbor close enough to block it, each one radiates to the full sky. Bottom: hour 73, after the swarm has gathered. The interior is blue, shielded by the bodies around it, and only the rim stays warm-colored; all twenty are alive.
Parameter
Value
Parameter
Value
N
20
Tset
273.15K
R
0.085m
Tenv
95K
M
2.0kg
Tdeath
233.15K
tnight
354h
Cth
2000JK−1
Δt
54.1s
Gsky
0.0108WK−1
vcmd
0.015ms−1
Qmax
2.0W
TABLE I: Model parameters.
Fig. 2: Overview of the learning system. Left: fourteen local observation channels enter a 64-unit LSTM whose state persists across thermal steps; a skip connection passes the input directly to the output layer, which produces a velocity command. ES perturbs the shared parameters θ and updates them from a single scalar returned after the full night. Middle: the simulated swarm, colored by sky exposure χ . Right: the control loop deployed on every robot, using the trained θ unchanged.
Training budget
Training budget
Resulting threshold
Static
( × bound)
(pcm-h)
(pcm-h)
survivors
1.10
151.9
155.8
14/20
1.18
163.0
155.7
14/20
1.40
193.3
176.4
19/20
TABLE II: Sensitivity of the survival threshold to the training budget. The budget is an input, the threshold an outcome; both are in pcm-h. One policy per budget, same seed.
Fig. 3: Emergence from a scattered start; all rows share one initial configuration. Color is each robot’s time-averaged exposure ⟨χ⟩ up to that frame; crosses are dead robots. Under each frame: robots alive and the swarm’s mean ⟨χ⟩ . The thresholds of Table III are a separate quantity, bisected on the budget rather than read off these runs. Top: the learned policy gathers within five hours and circulates, every average converges, 20/20 survive. Middle: the gather-only rule, which aggregates but never circulates, forms the same huddle, but a core–rim gradient persists all night and nine robots are lost. Bottom: under zero command the swarm stays scattered and all are lost by 225 h.
Condition
n
Threshold
× bound
Dispersed, no command
—
328.5
2.38
Motionless huddle
—
246.0
1.78
Gather-only rule
—
250.3
1.81
Learned, packed, shared ref.
3
155.7±1.3
1.13
Learned, packed, no shared ref.
2
154.3±1.6
1.12
Learned, scattered, shared ref.
2
155.4±0.6
1.13
TABLE III: Survival thresholds in pcm-h. Learned entries pool converged runs.
Fig. 4: Transfer to larger swarms. Change in survival requirement relative to N=20 , the size at which all five policies were trained; they were evaluated at N=30 and N=40 without retraining. Squares are the mean across policies with bars spanning best and worst; faint points are the individual policies. The grey region is the change available from the equalization bound.
Channel zeroed
s0
s1
s2
None
157.1
155.8
155.4
Neighbor stored energy
161.4
157.5
155.8
Neighbor temperature
179.8
154.9
159.6
Own sky exposure
157.1
159.2
242.6
Both neighbor channels
175.5
154.9
268.4
TABLE IV: Observation ablation results, reported per policy.