Y(4626) as a molecular state from interaction Ds∗D¯s1(2536)-DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)-{D}_s{\bar{D}}_{s1}(2536)$$\end{document}

Recently, a new structure Y(4626) was reported by the Belle Collaboration in the process e+e-→Ds+Ds1(2536)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow D_s^+D_{s1}(2536)^-$$\end{document}. In this work, we propose an assignment of the Y(4626) as a Ds∗D¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)$$\end{document} molecular state, which decays into the Ds+Ds1(2536)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_s^+D_{s1}(2536)^-$$\end{document} channel through a coupling between Ds∗D¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)$$\end{document} and DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_s{\bar{D}}_{s1}(2536)$$\end{document} channels. With the help of the heavy quark symmetry, the potential of the interaction Ds∗D¯s1(2536)-DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)-{D}_s{\bar{D}}_{s1}(2536)$$\end{document} is constructed within the one-boson-exchange model, and inserted into the quasipotential Bethe–Salpeter equation. The pole of obtained scattering amplitude is searched for in the complex-energy plane, which corresponds to a molecular state from the interaction Ds∗D¯s1(2536)-DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)-{D}_s{\bar{D}}_{s1}(2536)$$\end{document}. The results suggest that a pole is produced near the Ds∗D¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)$$\end{document} threshold, which exhibits as a peak in the invariant mass spectrum of the DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}_s{\bar{D}}_{s1}(2536)$$\end{document} channel at about 4626 MeV. It obviously favors the Y(4626) as a Ds∗D¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)$$\end{document} molecular state. In the same model, other molecular states from the interaction Ds∗D¯s1(2536)-DsD¯s1(2536)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D}^*_s{\bar{D}}_{s1}(2536)-{D}_s{\bar{D}}_{s1}(2536)$$\end{document} are also predicted, which can be checked in future experiments.