Can calculate the average intensity just by |E|^2?

[beizhixing]

Hi, I want to reproduce the results in the reference below using MEEP.
(Slussarenko, S.,et al. Guiding light via geometric phases. Nature Photon 10, 571–575 (2016). https://doi.org/10.1038/nphoton.2016.138).

I use the code blow to obtain the ampulitude of E field, but I can not obtain the same picture as the reference. In my simulations, the waveguide eventually becomes multi-modal, as indicated by the appearance of breathing versus propagation.

Can you help me to solve it?

The code I used：

XY_data = [sim.get_array(center=mp.Vector3(0,0,0), size=mp.Vector3(d, gp, 0), component=c, cmplx=True) for c in [mp.Ex, mp.Ey, mp.Ez]]
XY_intensity = np.sqrt(np.abs(XY_data[0])**2 + np.abs(XY_data[1])**2 + np.abs(XY_data[2])**2)
plt.figure(dpi=1200)
plt.figure(figsize=(8, 4))
#plt.subplot(1,2,1)
plt.imshow(XY_intensity.transpose(), extent=[0, d,-gp/2, gp/2], cmap='jet')

Picture from the reference mentioned above.

Picture from my simulation.