Results are presented from a search for the Higgs boson decay H $\to$ Z$\gamma$, where Z $\to$$\ell^+\ell^-$ with $\ell$ = e or $\mu$. The search is performed using a sample of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV, recorded by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. Events are assigned to mutually exclusive categories, which exploit differences in both event topology and kinematics of distinct Higgs production mechanisms to enhance signal sensitivity. The signal strength $\mu$, defined as the product of the cross section and the branching fraction [$\sigma($pp $\to$ H$)\mathcal{B}($H $\to$ Z$\gamma)$] relative to the standard model prediction, is extracted from a simultaneous fit to the $\ell^+\ell^-\gamma$ invariant mass distributions in all categories and is found to be $\mu$ = 2.4 $\pm$ 0.9 for a Higgs boson mass of 125.38 GeV. The statistical significance of the observed excess of events is 2.7 standard deviations. This measurement corresponds to $\sigma($pp $\to$ H$)\mathcal{B}($H $\to$ Z$\gamma)$ = 0.21 $\pm$ 0.08 pb. The observed (expected) upper limit at 95% confidence level on $\mu$ is 4.1 (1.8). The ratio of branching fractions $\mathcal{B}($H $\to$ Z$\gamma) / \mathcal{B}($H $\to$ $\gamma\gamma)$ is measured to be 1.5 $^{+0.7}_{-0.6}$, which agrees with the standard model prediction of 0.69 $\pm$ 0.04 at the 1.5 standard deviation level.
Forum