Context:Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion's optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the \texttt{Icarus} software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations.

We present the first measurement of the gravitational quadrupole moment of the companion star of a spider pulsar, namely the black widow PSR J2051$-$0827. To this end we have re-analysed radio timing data using a new model which is able to account for periastron precession caused by tidal and centrifugal deformations of the star as well as by general relativity. The model allows for a time-varying component of the quadrupole moment, thus self-consistently accounting for the ill-understood orbital period variations observed in these systems. Our analysis results in the first detection of orbital precession in a spider system at $\dot{\omega} = -68.6_{-0.5}^{+0.9}$ deg/yr and the most accurate determination of orbital eccentricity for PSR J2051$-$0827 with $e = (4.2 \pm 0.1) \times 10^{-5}$. We show that the variable quadrupole component is about 100 times smaller than the average quadrupole moment $\bar{Q} = -2.2_{-1}^{+0.6} \times 10^{41}$ kg.m$^2$. We discuss how accurate modelling of high precision optical light curves of the companion star will allow its apsidal motion constant to be derived from our results.