We investigated Einstein–Gauss–Bonnet (EGB) 4D massive gravity coupled to nonlinear electrodynamics (NED) in an anti-de Sitter (AdS) background and found an exact magnetically charged black hole solution. The metric function was analyzed for different values of massive gravity parameters. The first law of black hole thermodynamics and generalized Smarr formula were verified, where we treated the cosmological constant as thermodynamic pressure. We defined vacuum polarization as the conjugate to NED parameter. To analyze the local stability of the black hole, we computed specific heat. We investigated the van der Waals-like/re-entrant phase transition of the black holes and estimated the critical points. We observed small black hole (SBH)/large black hole (LBH) and SBH/intermediate black hole (IBH)/LBH phase transitions. The Joule–Thomson coefficient, inversion, and isenthalpic curves were discussed. Finally, the minimum inversion temperature and the corresponding event horizon radius were obtained using numerical techniques.
