Possible Molecular States of $D^{*}_s\bar{D}^{*}_s$ System and Y(4140)

The interpretation of Y(4140) as a $D^{*}_s\bar{D}^{*}_s$ molecule is studied dynamically in the one boson exchange approach, where $\sigma$, $\eta$ and $\phi$ exchange are included. Ten allowed $D^{*}_s\bar{D}^{*}_s$ states with low spin parity are considered, we find that the $J^{PC}=0^{++}$, $1^{+-}$, $0^{-+}$, $2^{++}$ and $1^{--}$ $D^{*}_s\bar{D}^{*}_s$ configurations are most tightly bound. We suggest the most favorable quantum numbers are $J^{PC}=0^{++}$ for Y(4140) as a $D^{*}_s\bar{D}^{*}_s$ molecule, however, $J^{PC}=0^{-+}$ and $2^{++}$ can not be excluded. We propose to search for the $1^{+-}$ and $1^{--}$ partners in the $J/\psi\eta$ and $J/\psi\eta'$ final states, which is an important test of the molecular hypothesis of Y(4140) and the reasonability of our model. The $0^{++}$ $B^{*}_s\bar{B}^{*}_s$ molecule is deeply bound, experimental search in the $\Upsilon(1S)\phi$ channel at Tevatron and LHC is suggested.