Establishing the heavy quark spin and light flavor molecular multiplets of the $X(3872)$, $Z_c(3900)$ and $X(3960)$

Recently, the LHCb Collaboration reported a near-threshold enhancement, $X(3960)$, in the $D_s^+D_s^-$ invariant mass distribution. We show that the data can be well described by either a bound or a virtual state below the $D_s^+D_s^-$ threshold. The mass given by the pole position is $(3928\pm3)$ MeV. Using this mass and the existing information on the $X(3872)$ and $Z_c(3900)$ resonances, a complete spectrum of the $S$-wave hadronic molecules formed by a pair of ground state charmed and anticharmed mesons is established. Thus, pole positions of the partners of the $X(3872)$, $Z_c(3900)$ and the newly observed $D_s^+D_s^-$ state are predicted. Calculations have been carried out at the leading order of nonrelativistic effective field theory and considering both heavy quark spin and light flavor SU(3) symmetries, though conservative errors from the breaking of these symmetries are provided.