Searches for the exclusive decays of Higgs and $Z$ bosons into a vector quarkonium state and a photon are performed in the $\mu^+\mu^- \gamma$ final state with a proton$-$proton collision data sample corresponding to an integrated luminosity of $139$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN Large Hadron Collider. The observed data are compatible with the expected backgrounds. The 95% confidence-level upper limits on the branching fractions of the Higgs boson decays into $J/\psi \gamma$, $\psi(2S) \gamma$, and $\Upsilon(1S,2S,3S) \gamma$ are found to be $2.0\times10^{-4}$, $10.5\times10^{-4}$, and $(2.5,4.2,3.4)\times10^{-4}$, respectively, assuming Standard Model production of the Higgs boson. The corresponding 95% CL upper limits on the branching fractions of the $Z$ boson decays are $1.2\times10^{-6}$, $2.4\times10^{-6}$, and $(1.1,1.3,2.4)\times10^{-6}$. An observed 95% CL interval of $(-133,175)$ is obtained for the $\kappa_c/\kappa_\gamma$ ratio of Higgs boson coupling modifiers, and a 95% CL interval of $(-37,40)$ is obtained for $\kappa_b/\kappa_\gamma$.
Searches for the exclusive decays of the Higgs boson to an $\omega$ meson and a photon or a $K^{*}$ meson and a photon can probe flavour-conserving and flavour-violating Higgs boson couplings to light quarks, respectively. Searches for these decays, along with the analogous $Z$ boson decay to an $\omega$ meson and a photon, are performed with a $pp$ collision data sample corresponding to integrated luminosities of up to 134 fb$^{-1}$ collected at $\sqrt{s}=13$ TeV with the ATLAS detector at the CERN Large Hadron Collider. The obtained 95% confidence-level upper limits on the respective branching fractions are ${\cal B}(H\rightarrow\omega\gamma)< 5.5\times 10^{-4}$, ${\cal B}(H\rightarrow K^{*}\gamma)< 2.2\times10^{-4}$ and ${\cal B}(Z\rightarrow \omega\gamma)<3.9\times 10^{-6}$. The limits for $H\rightarrow \omega\gamma$ and $Z\rightarrow \omega\gamma$ are 370 times and 140 times the Standard Model expected values, respectively. The result for $Z\rightarrow \omega\gamma$ corresponds to a two-orders-of-magnitude improvement over the limit obtained by the DELPHI experiment at LEP.
Searches for the exclusive decays of the Higgs boson into $D^*\gamma$ and of the $Z$ boson into $D^0\gamma$ and $K^0_s\gamma$ can probe flavour-violating Higgs and $Z$ boson couplings to light quarks. Searches for these decays are performed with a $pp$ collision data sample corresponding to an integrated luminosity of $136.3$ fb$^{-1}$ collected at $\sqrt{s}=13$ TeV between 2016-2018 with the ATLAS detector at the CERN Large Hadron Collider. In the $D^*\gamma$ and $D^0\gamma$ channels, the observed (expected) 95$\%$ confidence-level upper limits on the respective branching fractions are ${\cal B}(H\rightarrow D^*\gamma)< 1.0 (1.2)\times 10^{-3}$, ${\cal B}(Z\rightarrow D^0\gamma)< 4.0 (3.4)\times 10^{-6}$, while the corresponding results in the $K^0_s\gamma$ channel are ${\cal B}(Z\rightarrow K^0_s\gamma)< 3.1 (3.0)\times 10^{-6}$.