We address a method of limiting neutron–mirror neutron mixing ($$\epsilon _{nn'}$$
ϵ
n
n
′
) by analyzing its effect on neutron star (NS) heating. This method employs observational bounds on the surface temperature of NSs to constrain $$\epsilon _{nn'}$$
ϵ
n
n
′
. It has been suggested that the bound obtained this way is so stringent that it would exclude any discovery of $$n-n'$$
n
-
n
′
oscillation in the currently planned terrestrial experiments at various laboratories. This conclusion motivated us to critically analyze this suggestion in more detail. In this note, we point out a very interesting new effect present in nearly exact mirror models, which can significantly affect this bound. The new element is that in nearly exact mirror models there is the mirror analog of $$\beta $$
β
decay, i.e. $$n' \rightarrow p' + e' + {\bar{\nu }}'_e$$
n
′
→
p
′
+
e
′
+
ν
¯
e
′
, which creates a cloud of mirror particles $$n'$$
n
′
, $$p'$$
p
′
, $$e'$$
e
′
, $$D'$$
D
′
and He$$'$$
′
inside the NS. The resulting $$e'$$
e
′
can “rob” the energy generated by the $$n \rightarrow n'$$
n
→
n
′
transition from the NS, via $$e-e'$$
e
-
e
′
scattering enabled by the presence of a (minute) millicharge in mirror particles. Such a tiny millicharge on mirror particles is highly likely in these models. This results in energy being emitted as unobserved mirror photons via fast mirror bremsstrahlung. whose effect is to relax the stringent bounds on $$\epsilon _{nn'}$$
ϵ
n
n
′
.