A candidate giant planet companion to the massive, young white dwarf GALEX J071816.4+373139 informs the occurrence of giant planets orbiting B stars

It has been suggested that giant planet occurrence peaks for stars with $M_{\ast}~\approx~3~M_{\odot}$ at a value a factor of 4 higher than observed for solar-mass stars. This population of giant planets predicted to frequently orbit main-sequence B stars at $a$ ~ 10au is difficult to characterize during the few hundred million years while fusion persists in their host stars. By the time those stars become massive, young white dwarfs, any giant planets present would still be luminous as a consequence of their recent formation. From an initial sample of 2195 Gaia-identified massive, young white dwarfs, we use homogeneous Spitzer Infrared Array Camera (IRAC) photometry to search for evidence of unresolved giant planets. For 30 systems, these IRAC data provide sensitivity to objects with $M~\lesssim$ 10 $M_{\text{Jup}}$, and we identify one candidate with ~ 4 $M_{\text{Jup}}$ orbiting the white dwarf GALEX J071816.4+373139. Correcting for the possibility that some of the white dwarfs in our sample result from mergers, we find a giant planet occurrence $\eta_{\text{GP}}~=~0.11_{-0.07}^{+0.13}$ for stars with initial masses $M_{\ast}~\gtrsim~3~M_{\odot}$. Our occurrence inference is consistent with both the Doppler-inferred occurrence of giant planets orbiting 2 $M_{\odot}$ giant stars and the theoretically predicted factor of 4 enhancement in the occurrence of giant planets orbiting 3 $M_{\odot}$ stars relative to solar-mass stars. Future James Webb Space Telescope NIRCam observations of our sample would provide sensitivity to Saturn-mass planets and thereby a definitive estimate of the occurrence of giant planets orbiting stars with $M_{\ast}~\gtrsim~3~M_{\odot}$.