The newly observed Zcs(3985)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{cs}(3985)^- $$\end{document} state: in vacuum and a dense medium

We report on some properties of the newly observed charged hidden-charmed open strange Zcs(3985)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{cs}(3985)^- $$\end{document} state by BESIII Collaboration. Assigning the quantum numbers JP=1+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ J^{P} = 1^{+}$$\end{document} and the quark composition cc¯su¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c \bar{c} s\bar{u} $$\end{document} and considering it as the strange partner of the famous Zc(3900)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_c(3900) $$\end{document} state, we estimate the mass of the Zcs(3985)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{cs}(3985)^- $$\end{document} resonance in vacuum and compare it with the experimental data. We also investigate its mass, current coupling and vector-self energy in a medium with finite density. Our result on the mass in vacuum agrees well with the experimental data. We estimate the mass and current coupling of the b-partner of this state, Zbs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{bs}$$\end{document}, in the vacuum as well. For its mass we get mZbs=10732-46+97MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ m_{Z_{bs}}= 10732^{+97}_{-46}~\hbox {MeV}$$\end{document}, which may be checked via other nonperturbative approaches as well as future experiments. We present the dependence of the spectroscopic parameters of Zcs(3985)-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Z_{cs}(3985)^- $$\end{document} state on density and observe that these parameters are linearly changed with increasing in the density.


Introduction
The recently discovered charged resonances are good candidates of tetraquark systems as they can not be put in the spectrum of the usual cc mesons by any means. The observation of the ground state Z c (3900) was reported by the BESIII [1] and Belle [2] Collaborations in 2013, simultaneously. After these observations many theoretical studies were performed on different properties of these states (as examples see Refs. [3][4][5][6][7] and references therein). The next state in this group is the Z (4430) which discovered by the Belle and LHCb Collaborations [8][9][10][11][12] and fully studied in Refs. [4,[13][14][15][16][17][18] and references therein. The mass difference between Z (4430) and Z c (3900) is roughly equal to the mass a e-mail: kazem.azizi@ut.ac.ir (corresponding author) b e-mail: nuray@ibu.edu.tr difference between ψ(2S) and J/ψ, that respectively appear in their decay products, was the reason for consideration of Z (4430) to be the first radial excitation of the Z c (3900) state [4,18].
Very recently, as the first candidate, the BESIII Collaboration reported the observation of a charged hiddencharm tetraquark with strangeness, decaying into D − s D * 0 and D * − s D 0 in the process e + e − → K + (D − s D * 0 + D * − s D 0 ) [19]. This new member is called Z cs (3985) − and its quark composition is most likely ccsū. The measured mass is just above the thresholds of D sD * and D * sD . This value together with the value of the width reported by the experiment are: The charmonium-like open strange four quark systems were already in agenda of different theoretical studies [20][21][22][23][24]. In Ref. [20], assuming SU(3) flavor symmetry, the existence of strange partners to the Z c states denoted by Z cs was predicted. In this study, as partners of Z c (4100) and Z c (4200) hadrocharmonium states, the hypothetical strange hadrocharmonium resonances with quantum numbers J P = 0 + and 1 + were referred as Z cs (4250) and Z cs (4350), respectively. In Ref. [21], the mass of the the possible D s D * molecular state, as the strange partner of the observed Z ± c (3900) states, called Z + cs was predicted. In Ref. [22], using the QCD sum rule approach, and considering a state near to the D sD * /D * sD threshold, its four different decay channels were studied and by calculating the related strong coupling constants and decay widths, the total width of the proposed state was predicted. In Ref. [23], while studying the hidden-charm di-kaon decays of higher charmonia and charmonium-like state, the authors predicted the existence of hidden-charm and openstrange channels near the DD * s /D * D s and D * D * s /D * D * s thresholds. In Ref. [24], the authors calculated the spectrum of strange hidden-charm and bottom tetraquarks in the hadroquarkonium and the compact tetraquark models. However, there appeared various studies after the BESIII Collaboration observation report. Different assumptions and assignments have been made to clarify the nature and structure of the observed state [25][26][27][28][29][30][31][32][33]. In Ref. [25], by performing a combined analysis for the line shape of the recoil mass distribution of K + at five energy points, the authors find that the Z cs (3985) − can be explained as a reflection structure of charmed-strange meson D * s2 (2573), which is produced from the open-charm decay of Y (4660) with a D * s meson. In Ref. [26], the authors introduced the G U/V parity to construct the flavor wave functions of the Z cs states and predicted their masses. They also predicted the masses of the tetraquark states Z bs and Z bs , which are the U/V -spin partner states of Z b (10610) ± and Z b (10650) ± , respectively. In Ref. [27], using the one-boson-exchange model, it was concluded that the Z cs (3985) could not be a pure hadronic molucular state ofD s D * andD * s D and could consist of large components of compact nature. In Ref. [28], the quantum numbers I (J P ) = 1 2 (1 + ) were assigned to the new resonance Z cs by the mass analyses and the comparision of the results with the experimental data. After producing the phase shifts for all the considered channels, the results of Ref. [29], exclude the In Ref. [30], the states Z c (3900) and Zc(4020), as well as their strange partners Z cs (3985) and Z cs (4120) were studied. With an overall short-ranged contact potential in SU(3) flavor symmetry the authors showed that the pole structures of these states could be extracted from the experimental data. In Ref. [31], the authors studied the photoproduction of exotic states Z cs (3985) and Z cs and explored the possibility of searching for them in future electron-ion colliders (EIC). In Ref [32] the authors studied the structure of Z − cs (3985) state in the picture of D ( * )− s D ( * )0 molecular state and got a consistent result with the experiment. In Ref. [33], the authors calculated the mass spectra for theD ( * ) s D ( * ) molecular states and scqc tetraquark states with J P = 0 + , 1 + , 2 + ) and got consistent results with the experiment in both scenarios. In both the molecular and diquark-antidiquark pictures, their results suggested that there may exist two almost degenerate states, as the strange partners of the X (3872) and Zc(3900).
In the present study, we investigate the Z cs (3985) − (hereafter Z cs ) state both in the vacuum and a medium with finite density. Based on the provided information by the experiment that the new state is composed of ccsū and the measured central value for its mass is sightly higher than the thresholds D sD * and D * sD , the Z cs state can be considered as a compact tetraquark rather than the loosely bound molecular states D sD * and D * sD . With this information and the results of theoretical studies before the experimental report, Refs. [4,6,[20][21][22][23][24], we treat the Z cs as a compact tetraquark and assign the quantum numbers J P = 1 + and quark composition ccsū to this state and consider it as the strange partner of the famous Z c (3900) state. We calculate its mass and current coupling in vacuum. We also investigate this particle in a dense medium to estimate its vector self-energy, mass and current coupling at nuclear matter saturation density and study the changes of its parameters with respect to the density of the medium in order to more clarify its nature and structure. We use the technique of the QCD sum rules [34,35] to calculate the spectroscopic parameters of this state. Our calculations let us also estimate some parameters of the bpartner of Z cs state, namely Z bs , which may be in agenda of future experiments.
In in-medium sum rules, the medium effects are included via various operators in the nonperturbative parts of the operator product expansion (OPE) (see for example [36,37]). The new progresses recorded in the heavy ion collision experiments like the Relativistic Heavy Ion Collider (RHIC) at the Bookhaven National Laboratory (BNL) and the Large Hadron Collider (LHC) and some in-medium experiments likeP AN D A at Facility for Antiproton and Ion Research (FAIR) at GSI and the Nuclotron-based Ion Collider Facility at JINR (NICA) at Dubna, it will be possible to study not only the ground state particles but also the excited exotic states [38][39][40]. Hence, these experiments may provide good opportunities to investigate Z c (3900), Z (4430) and Z cs in order to clarify their nature and substructure and fix their quantum numbers.
In next section we provide some details of the calculations to obtain sum rules for the physical quantities under study. Sect. 3 encompasses the numerical analyses of the obtained expressions. The last section is devoted to the summary of the results and conclusions.

Sum rules for in-medium parameters of Z cs
The information provided by the BESIII collaboration lead us to assign the quantum numbers J P = 1 + to Z cs state and consider it as the strange partner of the famous Z c (3900) particle with quark composition ccsū. The following current can interpolate the Z cs state with the assumed nature: where abc and dec are anti-symmetric Levi-Civita symbols with a, b, c, d and e being color indices. Here, T stands for the transpose and C represents the charge conjugation operator.
The basic function that we need to evaluate is the following in-medium correlation function (CF): where p is the four-momentum of the Z cs state, T is the time ordering operator, and |ψ 0 is the ground state of nuclear medium. We will calculate this two-point CF once in terms of in-medium hadronic parameters and the second in terms of QCD fundamental parameters and the in-medium quark and gluon condensates. By matching the obtained results in two different pictures we will able to calculate the inmedium hadronic parameters in terms of the QCD degrees of freedom. We should note that although the QCD sum rule method [34,35] was originally formulated for calculations of different parameters of the standard mesons and baryons, it has been successfully applied for the tetraquarks and other multiquark systems giving consistent results with the experiments on the observed states. This approach has helped physicists to describe the results of experiments on the XYZ tetraquark systems provided by LHCb, BESIII, Belle and other related experiments and fix the nature, quark-gluon organization and quantum numbers of these states (for some of these studies see for instance [3][4][5][6][41][42][43][44][45][46] and references therein). Despite the successes of the QCD sum rules in applications to the tetraquark systems, there are some doubts on the applied approaches to the tetraquark molecular states. In Refs. [47,48], it is stated that the contributions at the orders O(1) and O(α s ) in OPE are exactly canceled out by the meson-meson scattering states at the hadronic side and the tetraquark molecular states begin to receive contributions at the order O(α 2 s ). However, in Ref. [49], by choosing an axialvector and a tensor currents, it is shown that the mesonmeson scattering states do not give contribution to the CF, while the tetraquark molecular states can saturate it. It is also stated that the Landau equation is of no use to study the Feynman diagrams in the QCD sum rules for the tetraquark molecular states and the OPE of these states receive contributions at leading orders.
The hadronic representation of CF is obtained by saturating it with a full set of hadronic state. By separating the ground state contribution, we get the following expression after performing the four-intergral over coordinate space: where p * is the in-medium momentum, m * Z cs is the inmedium mass and dots stand for the contributions of the higher states and continuum. The above expression can be further simplified by introducing the in-medium decay constant or current coupling f * Z cs defined by where ε μ is the polarization four-vector of Z cs . Using this definition leads to the summation over the polarization vector. Hence, Had is obtained for the hadronic side. To proceed, we need to define two kinds of self-energies. First one is the scalar selfenergy defined by the shift in the hadronic mass due to the nuclear medium, s = m * Z cs − m Z cs , which is a Lorentz scalar and the second one is the vector self-energy υ that transforms as a Lorentz four-vector and denotes the shift in the four-momentum of the particle due to the dense medium. The vector self energy enters to the calculations as p * μ = p μ − υ u μ (for more details about the scalar and vector self energies see [37]). Here, u μ is the velocity four-vector of the nuclear matter. We perform the calculations in the rest frame of the medium, where u μ = (1, 0). We should remark that we also choose to work at zero momentum limit of the particle's three momentum, for which the longitudinal and transverse polarizations become degenerate, and there is no way to distinguish these two polarizations from each other. For the finite and non-zero three momentum in a dense medium, the longitudinal and transverse polarizations are distinguishable and they should be calculated, separately. In terms of the new quantities, the hadronic part takes the form where μ 2 = m * 2 Z cs − 2 υ + 2 p 0 υ , with p 0 = p.u being the energy of the quasi-particle in the dense medium. We apply the Borel transformation to suppress the contributions of the higher states and continuum based on the standard prescriptions of the method. As a result, we get where M 2 is the Borel mass parameter. The QCD side of the CF is calculated in deep Euclidean region by the help of OPE. In this representation, the shortdistance (perturbative) effects are separated from the longdistance (nonperturbative) effects. To this end, we insert the explicit expression of the interpolating current J μ (x) into the CF and perform all the possible contractions among the quark fields using the Wick's theorem. As a result, we get an expression in terms of the in-medium light and heavy quark propagators: whereS q(c) = C S T q(c) C, and q stands for u, or s quark. We make use of the in-medium light and heavy quark propagators in coordinate space. Against our previous studies on the properties of Z c (3900) and Z (4430) states [6,18], we keep the light quark mass in the calculations in order to take into account especially the strange quark mass effects: and where χ i q andχ j q are the Grassmann background quark fields and we use the short-hand notation with F A μν being the classical background gluon fields. Here, t i j,A = λ i j,A 2 and λ i j,A are the Gell-Mann matrices. Now, we insert the explicit forms of the light and heavy quark propagators into the CF (9) and perform all the fourintegrals over x to go to the momentum space. We then perform the Borel transformation in Euclidean region and apply the continuum subtraction procedure. Here, we do not present the details of calculations, for which one may refer to Refs. [6,36,50]. After these standard but lengthy procedures, one gets the QCD side of the CF in terms of different Lorentz structures as where the functions ϒ QC D i (i = 1, ..., 5) in Eq. (13) are given by with s * 0 being the in-medium continuum threshold. The spectral densities ψ QC D i (s, ρ) are very lengthy functions having various components: perturbative ( pert) as well as twoquark (qq), two-gluon (gg) and mixed quark-gluon (qgq) condensates, i.e, The explicit forms of the above spectral densities, as an example, corresponding to the structure g μν are presented in the "Appendix".
The next step is to determine the working regions of the auxiliary parameters M 2 and s * 0 that the sum rules depend on. To fix them, the standard criteria of the method such as, mild dependence of the results on the auxiliary parameters, convergence of the OPE and pole dominance (for details see for instance Refs. [6,18]) are used. These requirements lead to the intervals 3 GeV 2 ≤ M 2 ≤ 5 GeV 2 and 18.3 GeV 2 ≤ s * 0 ≤ 20.1 GeV 2 for the Borel parameter and continuum threshold in the charmed channel, respectively. Making use of all the input parameters and the working windows for the auxiliary parameters, we plot the vacuum mass and current coupling of Z cs state (at ρ → 0 limit) with respect to M 2 at different fixed values of the continuum threshold in figures 1 and 2. From these figures, we see Comparing the value of mass obtained in the vacuum with the experimental value of the BESIII Collaboration, 3982.5 +1.8 −2.6 ± 2.1 MeV, we see a nice agreement between the two values. The value provided for the current coupling constant in above equation can be served as one of the main inputs in investigation of decay products of the Z cs state. Note that as it is clear from the above values for the mass and current coupling, the uncertainty in the value of current coupling is much more compared to that of the mass. This can be attributed to the fact that the mass is found from the ratio of two sum rules in Eq. (16) that kill the uncertainties of each other and lead to a more accurate result, while the current coupling is obtained using only one sum rule that leads to a relatively large uncertainty. However, all uncertainties remain within the limits allowed by the method.
At saturation density of the nuclear matter, the mass reduces considerably. At this point, we get showing that the particle gains a large vector self-energy, referring to the considerable vector repulsion of this state by the medium at saturated nuclear matter density. These results also indicate that the particle receives a scalar self energy of −143 MeV. This negative shift in the mass shows the scalar attraction of the Z cs state by the dense medium. We show the dependence of m * Z cs and υ on ρ/ρ sat at mean values of the continuum threshold and Borel parameter in figure 3. We should note that our sum rules give reliable results up to roughly ρ/ρ sat = 1.5, after which the results are not decent and we use some fit functions to extend the results up to ρ = 5ρ sat , which corresponds to the density of the cores of neutron stars (see [18] for more details). The fit functions for the in-medium mass, current coupling and vector self-energy are obtained as where x = ρ/ρ sat . These functions and also the related figure, 3, show that the mass, current coupling and vector selfenergy linearly depend on the density: the mass and current coupling decrease with increasing in the density, while the vector self-energy increases with the density. Finally, we extend our analyses to predict the mass and decay constant of the b-partner of the observed Z cs (3985) − state. Our calculations show that the working windows 9 GeV 2 ≤ M 2 ≤ 12 GeV 2 and 126 GeV 2 ≤ s 0 ≤ 130 GeV 2 satisfy all the requirements of the method. Using the input values in b-channel we get in vacuum. The future experiments may search for the charged hidden-bottom open strange state, Z − bs , with quantum numbers J P = 1 + and quark composition bbsū around the 10730 MeV mass. The Belle Collaboration has previously reported the observation of two narrow struc-tures in the mass spectra of ϒ(nS)π + π − , n = 1, 2, 3 and h b (m P)π + π − , m = 1, 2 pairs called Z b (10610) and Z b (10650) in the decays ϒ(5S) → ϒ(nS)π + π − and ϒ(5S) → h b (m P)π + π − [66,67] and assigned the spinpartiy J P = 1 + to these states. The measured masses for these states are The Belle Collaboration may search for the strange Z bs states, as well.

Summary and conclusions
Inspired by the recent observation of the charged hiddencharmed open strange Z cs (3985) − state by the BESIII Collaboration, we performed a QCD sum rule analysis for the spectroscopic parameters of this state both in vacuum and a medium with finite density. For its mass in vacuum we obtained m Z cs = 3989 +36 −17 MeV, which is in accord with the experimental value, m = 3982.5 +1.8 −2.6 ± 2.1 MeV. From this result, we conclude that the assignments J C = 1 + and ccsū for the quantum numbers and quark structure of this state works well. We calculated the current coupling constant of this state in vacuum as well, which can be used in investigation of the electromagnetic, weak or strong decays of Z cs state as an important input parameter. We also estimated the mass and current coupling of the b-partner of this state, Z bs , in vacuum and in the same picture. The obtained value for the mass of this hypothetical resonance, m Z bs = 10732 +97 −46 MeV, may be checked via different nonperturbative approaches as well as future experiments.
We also reported the values of mass, current coupling constant and vector self-energy of Z cs state at the saturation density of a nuclear medium. We investigated the behavior of the spectroscopic parameters of this state in terms of the density and observed that they are linearly changed by increasing in the value of the density: the mass slightly decreases while the vector self-energy considerably increases with respect to the density. Our results may be helpful for future heavy ion collision experiments and those aiming to study the behavior of the exotic states at finite densities.
Acknowledgements The authors thank TUBITAK for the partial support provided under the Grant No. 119F094.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All the numerical and mathematical data have been included in the paper and we have no other data regarding this paper.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .