Magnetic dipole moments of the hidden-charm pentaquark states: Pc(4440)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_c(4440)$$\end{document}, Pc(4457)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_c(4457)$$\end{document} and Pcs(4459)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{cs}(4459)$$\end{document}

In this work, we employ the light-cone QCD sum rule to calculate the magnetic dipole moments of the Pc(4440)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_c(4440)$$\end{document}, Pc(4457)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_c(4457)$$\end{document} and Pcs(4459)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{cs}(4459)$$\end{document} pentaquark states by considering them as the diquark–diquark–antiquark and molecular pictures with quantum numbers JP=32-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P = \frac{3}{2}^-$$\end{document}, JP=12-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P = \frac{1}{2}^-$$\end{document} and JP=12-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^P = \frac{1}{2}^-$$\end{document}, respectively. In the analyses, we use the diquark–diquark–antiquark and molecular form of interpolating currents, and photon distribution amplitudes to obtain the magnetic dipole moment of pentaquark states. Theoretical examinations on magnetic dipole moments of the hidden-charm pentaquark states, are essential as their results can help us better figure out their substructure and the dynamics of the QCD as the theory of the strong interaction. As a by product, we extract the electric quadrupole and magnetic octupole moments of the Pc(4440)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_c(4440)$$\end{document} pentaquark. These values show a non-spherical charge distribution.

diquark-antiquark and molecular form of interpolating currents, and photon distribution amplitudes to obtain the magnetic dipole moment of pentaquark states. Theoretical examinations on magnetic dipole moments of the hidden-charm pentaquark states, are essential as their results can help us better figure out their substructure and the dynamics of the QCD as the theory of the strong interaction. As a by product, we extract the electric quadrupole and magnetic octupole moments of the P c (4440) pentaquark. These values show a non-spherical charge distribution.

Motivation
With the experimental discovery of the exotic X(3872) state, i.e., state that cannot be interpreted by the conventional meson or baryon picture, a new era began in high energy physics. Since the discovery of this state, numerous exotic states have been added to the family of particles that have been experimentally reported. Study of exotic particles that are composed of tetraquark, pentaquark or hybrid states is among the most attractive subjects of the hadron physics. Experimental data collected with various collaborations in recent years and theoretical developments obtained in different theoretical models form the field of rapidly growing exotic studies [1][2][3][4][5][6][7][8][9][10][11].
In 2019, LHCb collaboration reported three new narrow pentaquark states [12] as a e-mail: ulasozdem@aydin.edu.tr (corresponding author) But, spin and parity of the P cs (4459) state have not been determined yet. As regards to their decay products, one can easily conclude that these newly discovered four states consist of at least five quarks,ccuud orccuds, therefore they are perfect candidates of hidden-charm pentaquark states. After the experimental discovery several phenomenological models have been adopted to calculate the spectroscopic parameters, decays and production mechanisms of the pentaquarks, like the QCD sum rule, the meson-exchange model, the quark delocalization model, and so on . However, the substructure of these states are not determined yet.
In other words, in order to understand the internal structure and nature of these particles, different properties should be studied besides their decay channels and spectroscopic properties. For instance, investigating their electromagnetic form factors may provide important insights on this point. Electromagnetic form factors or multi-pole moments of hadrons are important parameters in study of their electromagnetic structure, and also they can ensure important knowledge about the dynamics of the QCD at low energy region. Electromagnetic multi-pole moments, especially magnetic dipole moment, are also a crucial part in the calculation of J/ψ photo-production cross sections, which can provide an independent analysis of the hidden-charm pentaquark states. Investigating electromagnetic features of exotic resonances is relatively new topic. However, in the literature there are a few studies where the electromagnetic properties of the hidden-charm pentaquark states are investigated [56][57][58][59][60]. In Ref. [56], the magnetic dipole moment of the hidden-charm pentaquark states have been extracted in the diquark-diquark-antiquark, diquark-triquark and molecular configuration with J P = 1 2 ± , 3 2 ± , 5 2 ± and 7 2 + quantum numbers in the different color-flavor structure. In Ref. [57], the magnetic dipole, electric quadrupole and magnetic octupole moments of the P c (4380) pentaquark with J P = 3 2 − quantum numbers have been obtained in the diquarkdiquark-antiquark and molecular pictures in the framework of the light-cone QCD sum rule (LCSR). In Ref. [58], they acquired the the ground state of hidden-charm pentaquarks with J P = 3 2 − quantum numbers and their associated magnetic dipole moments and electromagnetic couplings, of interest to pentaquark photoproduction experiments in the framework of the constituent quark model. In Ref. [59], the magnetic dipole moment of the P c (4312) pentaquark state have been extracted in the molecular and diquark-diquarkantiquark pictures via the LCSR with J P = 1 2 − quantum numbers. In Ref. [60], they achieved the magnetic dipole moment of the P c (4312) pentaquark state in the cD molecular picture by means of the QCD sum rule (QCDSR) in the external weak electromagnetic field.
In this study, the magnetic dipole moments of the pentaquark states P c (4440), P c (4457) and P cs (4459) (hereafter we will show these states as P c1 , P c2 and P cs , respectively) are obtained by employing the diquark-diquark-antiquark and molecular interpolating currents within the LCSR.
The LCSR is a powerful tool to obtain qualitative and quantitative information about hadron properties. In this method, hadrons are analytically studied by matching their behaviors at high and low energies. This matching provides characteristic properties of hadrons from QCD. The information obtained by LCSR has been used as an input to other theoretical and experimental approaches for a long time [61][62][63].
The rest of paper is structured as follows. In Sect. 2, the method used for the calculations are described. In Sect. 3, we carry out numerical analysis of the acquired LCSR for the magnetic dipole moment of pentaquark states. The explicit expressions of the magnetic moment of the P c1 pentaquark in the diquark-diquark-antiquark picture is presented in Appendix A.

Magnetic dipole moments of the pentaquark states via LCSR
In this section, we briefly introduce the LCSR method used to extract the magnetic dipole moments of the P c1 , P c2 and P cs pentaquark states. In order to obtain the magnetic dipole moment of the corresponding states with the QCD sum rule method, we begin by writing the correlation function suitable for the calculations. This correlation function is obtained in two representations which are called hadronic and QCD representations. The LCSR for the physical quantities are obtained from the matches of the coefficients of the same Lorentz structures achieved on both representations with the help of the quark-duality ansatz.
2.1 Formalism of the P c2 and P cs states The current correlation function for the P c2 and P cs pentaquark states is written as where J P c2(cs) (x) is the interpolating current of P c2 or P cs pentaquark state. In the diquark-diquark-antiquark and molecular pictures with quantum number J P = 1 2 − , they are written where C is the charge conjugation matrix; and a, b... are color indices.
As we mentioned at the beginning of this section, in LCSR studies we need to evaluate the correlation function at both hadron and quark-gluon degrees of freedom. At the hadron level, we embedding a complete set of intermediate pentaquark states into the correlation function to acquire the hadronic representation, and isolate the ground state pentaquark states, and get the results: The matrix element P c2(cs) ( p, s) | P c2(cs) ( p + q, s) γ entering Eq. (5) can be parameterized in terms of Lorentz invariant form factors as follows: where ε and q are the polarization vector and momentum of the photon, respectively. Substituting Eq. (6) in Eq. (5) for hadronic side we get with 0 | J P c2(cs) | P c2(cs) ( p, s) = λ P c2(cs) γ 5 u( p, s). The Lorentz invariant form factors f 1 (q 2 ) and f 2 (q 2 ) are related to the magnetic form factor by the relation, At q 2 = 0, the magnetic form factor is acquired associated with the functions f 1 (q 2 ) and f 2 (q 2 ) as: The magnetic dipole moment, μ P c2(cs) , is defined in the following way: We observe from Eq. (7) that the correlation function contains many structures, any of them can be selected in obtaining magnetic dipole moment of the P c2(cs) pentaquark, and thus we decide on the structure ε /q /. As a result, the correlation function can be expressed with regards to the magnetic dipole moment of P c2(cs) pentaquark as, At the QCD level, we contract the relevant quark fields in the correlation function with the help of the Wick's theorem, where with S q(c) (x) being the full light and charm quark propagators. The relevant propagators are given as [64] where K i are modified the second kind Bessel functions and G μν is the gluon field strength tensor. The QCD representation of the correlation function can be obtained with the help of photon distribution amplitudes (DAs) according to quark-gluon properties and after performing the Fourier transform to transfer the calculations to the momentum space.
As a final step, by applying the double Borel transform on the variables − p 2 and ( p +q) 2 and choosing the coefficients of the same Lorentz structures in both QCD and hadronic representations and matching them employing the quark-hadron duality approach, we obtain the desired LCSR for magnetic dipole moments of P c2 and P cs states: The functions are quite lengthy, therefore the explicit expression of this function are not presented here.

Formalism of the P c1 state
In this subsection we derive the LCSR for the magnetic dipole moment of the P c1 pentaquark state. For this purpose, we consider following correlation function, where J P c1 μ(ν) is the interpolating current of P c1 pentaquark with J P = 3 2 − quantum numbers. In the diquark-diquarkantiquark and molecular pictures, it is given as The hadronic side of the correlation function is written as, The matrix element of the interpolating current between the vacuum and the P c1 pentaquark is defined as where λ P c1 is the residue P c1 pentaquark and u μ ( p, s) is the Rarita-Schwinger spinor. The transition matrix element P c1 ( p) | P c1 ( p + q) γ entering Eq. (24) can be parameterized in terms of four Lorentz invariant form factors as follows [65][66][67][68]: In principle, we can obtain the final expression of the hadronic side of the correlation function using the above equations, but we encounter two difficulties: not all Lorentz structures are independent and the correlation function also includes contributions of spin-1/2 and these undesirable contributions must be eliminated. To remove undesirable contributions coming from the spin-1/2 particles and obtain only independent structures in the correlation function, we apply the ordering for Dirac matrices as γ μ p /ε /q /γ ν and remove terms with γ μ at the beginning, γ ν at the end and those proportional to p μ and p ν [69]. As a result, using Eqs. (22)- (26) the hadronic side take the form, The final form of the hadronic side in terms of the selected structures in momentum space is: Had where Had 1 and Had 2 are functions of the form factors F 1 (q 2 ) and F 2 (q 2 ), respectively; and other independent structures is represented by dots.
The magnetic, G M (q 2 ), form factor is defined in terms of the form factors F i (q 2 ) in the following way [65][66][67][68]: where τ = − q 2 4m 2 P c1 . At q 2 = 0, the magnetic dipole moment is obtained in terms of the functions F 1 (0) and F 2 (0) as: The magnetic dipole moment, (μ P c1 ), is defined in the following way: The next step is to calculate the correlation function in Eq. (22) in terms of quark-gluon parameters. When we apply the same procedures as in the previous subsection, we get the following result: As a result, the QCD side of the correlation function in terms of the selected structures is obtained as The following processes are applied as described in the previous subsection and magnetic dipole moment results are obtained in LCSR. The QCD and hadronic representations of the correlation function are then matched employing quarkduality assumption. By equating the coefficients of the structures g μν p /ε /q / and g μν ε /q /, respectively for the F 1 and F 2 we obtain sum rules for these two form factors. As a result, we get, The explicit expressions of the LCSR for the F Di 1 and F Di 2 are presented in the Appendix A. We are now ready to move on to numerical analysis.

Numerical analysis and discussions
This section is devoted to the numerical computations for the magnetic dipole moments of the P c1 , P c2 and P cs pentaquark states. We use m u = m d = 0, m s = 96 +8 −4 MeV, m c = 1.275±0.02 GeV [70] [48] and λ Mol P cs = 1.13 × 10 −3 GeV 6 . Another set of main input parameters are the photon wavefunctions of different twists, entering the DAs. These wavefunctions are given in Ref. [71].
Except the above mentioned input parameters, the estimations for the magnetic dipole moments of pentaquark states  iliary parameters. The continuum threshold is considered to be the point where the excited states and continuum begin to contribute to the correlation function. The upper and lower bound of the Borel parameter is decided by demanding that both the contributions of the higher states and continuum are adequately suppressed and the contributions coming from higher dimensional terms are small. Our numerical analysis leads to the conclusion that these requirements are fulfilled in the regions shown below for the considered pentaquark states.
By having the values of all input parameters, we can start carrying out numerical computations. In Fig. 1, as an example, we depict the dependencies of the magnetic dipole moments in diquark-diquark-antiquark picture on M 2 and s 0 . As is seen, the deviation of the results in connection with the s 0 is remarkable however there is much less dependence of the physical observables under consideration on the M 2 in its working window. Our final results for the magnetic dipole moments are presented in Table 1. Our results cover errors originated from the uncertainty of the determinations of auxiliary parameters (M 2 and s 0 ) and other input parameters used in the analyses. The numerical values of the magnetic dipole moments acquired by means of two configurations differ considerably from each other, which can be used to determine the fundamental structure of these pentaquark states. That is to say, as many theoretical models give compatible results on the decay channels and spectroscopic parameters with the experimental data preventing us assigning precise substructure for pentaquarks, the experimental measurement of the magnetic dipole moments of these pentaquarks indeed can help us exactly distinguish their substructures. Since there are no theoretical result and experimental data for the results we obtained for these pentaquark states, we cannot compare our results. However, we may compare the result of the P cs state with the P c (4312) state magnetic dipole moments. It can be seen that the quark configurations of the P cs and P c (4312) particles are similar and therefore the magnetic dipole moment results of these particles can be expected to be close to each other. In Ref. [59], the magnetic dipole moment of P c (4312) in the diquark-diquark-antiquark and molecular pictures are extracted as μ P c = 0.40 ± 0.15 μ N and μ P c = 1.98 ± 0.75 μ N , respectively. The magnetic dipole moment of the P c (4312) pantaquark state has also been extracted via the QCD sum rule and its extension in the weak electromagnetic field by employing a molecular type interpolating current in Ref. [60]. The numerical value was obtained as μ P c = 0.59 +10 −0.20 . As one can see from these estimations, the numerical values for the magnetic dipole moments of P c (4312) and P cs states acquired in the present study are close to each other and we see a reasonable SU(3) flavor violation which is roughly %15. The obtained results in both pictures are considerably different compared to the result of Ref. [60]. Comparing the results obtained using different theoretical models with our results can give an idea about the consistency of our predictions.
As a by product, we also obtain the electric quadrupole (Q P c1 ) and magnetic octupole (O P c1 ) moments of the P c1 pentaquark as Just like the magnetic dipole moment results, we can see that the electric quadrupole and magnetic octupole moments results obtained by using two different pictures are different from each other. The Q P c1 and O P c1 of the P c1 pentaquark state showing a non-spherical charge distribution. In summary, the observation of new pentaquark states such as P c (4312), P c (4440), P c (4457) and P cs (4459), ensures a new platform to investigate the exotic states in QCD. There are different interpretations of their inner structures and quantum numbers, and these should be shedded light on with further research. In the present work, stimulated by the observation of the hidden-charm pentaquark states we have achieved the magnetic dipole moments of the P c (4440), P c (4457) and P cs (4459) by considering them as diquark-diquarkantiquark and molecular pictures with quantum numbers J P = 3 2 − , J P = 1 2 − and J P = 1 2 − by means of the light-cone QCD sum rule, respectively. As a by product, the electric quadrupole and magnetic octupole moments of the P c (4440) pentaquark have also extracted. Our predictions for electromagnetic multi-pole moments of pentaquark states may be checked via different theoretical models such as Lattice QCD, chiral perturbation theory etc. With the latest developments in the experimental side, we hope that we will be able to measure the magnetic dipole moment of newly formed multi-quark states, especially pantaquark states in the future. Any experimental measurements of the electromagnetic multi-pole moments of the hidden-charm pentaquark states and comparison of the obtained results with the predictions of the this study may provide as helpful knowledge on the internal structure of the these states as well as the non-perturbative behaviors of the strong interaction at the low-energy region.
Acknowledgements We thank Hua-Xing Chen for providing us the numeric value of the residue of the Pcs(4459) state.

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 .

Appendix A: Explicit forms of the F Di 1 and F Di 2 functions
In this appendix, we give the explicit expressions for the F Di 1 and F Di 2 functions: We should also point out that in the above expressions, for simplicity we have only given the terms that give important contributions to the numerical values of the magnetic moments, and have not presented many higher dimensional contributions, although they have been considered in the numerical calculations.
The where A represents the corresponding photon distribution amplitudes.