Identify Hadronic Molecule States by Neural Network

Neural networks are trained to judge whether or not an exotic state is a hadronic molecule of a given channel according its line-shapes. This method performs well in both trainings and validation tests. As applications, it is applied to study $X(3872)$, $X(4260)$ and $Z_c(3900)$. The results show that $Z_c(3900)$ should be regarded as a $\bar{D}^* D$ molecular state but $X(3872)$ not. As for $X(4260)$, it can not be a molecular state of $\chi_{c0}\omega$. Some discussions on $X_1(2900)$ are also provided.


I. INTRODUCTION
Exotic hadron states, refer to those hadron states that do not appear to fit with the expectations for an ordinary qq or qqq hadron in the quark model. A situation is that an exotic state may have multiquark constituents, e.g., its valence quarks are qqqq or qqqqq. One explanation is that these states are compact or "elementary" states, where quarks and anti-quarks are building blocks, they form a compact core and interact with each other by exchanging gluons. However, since these exotic states usually appear near the two-hadron thresholds, a natural explanation is that they are hadronic molecules, which means the building blocks of these exotic states are usual hadrons and they interact with each other by exchanging color neutral forces. So, a central question of researching exotic states is to decide whether an exotic state is a molecular state or an "elementary" state (for reviews, see for example, Refs. [1,2]).
Neural network is a machine learning algorithm which works for regression and classification problems. It has been applied to many fields in physics, like nuclear physics [3], high energy physics experiments [4], as well as hydrodynamics [5]. Recently, it also been applied to high energy physics phenomenology to study, for example, intermediate states in πN scattering [6], N N scattering [7] as well as extracting scattering length and effective range of exotic hadron states [8]. Inspired by these applications, an idea is using neural networks to identify hadronic molecular states.
In this work, the neural networks are trained by the invariant mass spectra generated artificially with labels "0" for molecular states and "1" for elementary states. During training, the validation tests will be done at ev- * Electronic address: haochen0393@pku.edu.cn † Electronic address: 1701110076@pku.edu.cn ‡ Electronic address: zhenghq@scu.edu.cn ery training epoch to monitor the network performances and to avoid overfitting. After training, the invariant mass spectrum from real experimental data will be put into trained neural network and the output which is a number from 0 to 1 will describe the possibility to be an elementary state. That means, if the output is closer to 0, then the resonance is more like a molecular state. On the contrary, if the output is closer to 1, then the resonance is more like an elementary state.
In 2003, X(3872) was firstly detected by Belle [9] in the J/ψπ + π − final state. Then, it was also observed near theD * D threshold [10]. Similar situation happened on Z c (3900), which is both observed in the final state of J/ψπ [11] andD * D [12]. So, it is worthwhile to ask whether or not these two exotic states are hadronic molecules ofD * D. After X(3872) and Z c (3900) observed, many experiments are done to study their properties [13][14][15]. Furthermore, X(4260) has been detected in the J/ψππ [16] and χ c0 ω [17] final state. Since it is near the threshold of ωχ c0 , a valuable question is that wither it is a molecule of ωχ c0 . All these three states have data in at least two channels, especially in the final state whose threshold is near the resonance. So, in this work, the trained neural networks are applied on these three hadronic states. This paper is organized as follows, In Sec. II, the method of training data generation and labelling will be introduced. In Sec. III, the structure of neural network classifier is discussed. The results of training, validation tests, as well as applications are displayed in Sec. IV. Finally, some conclusions and outlooks are given.

II. TRAINING DATA GENERATION
A neural network should be trained by many data with labels, means that people should not only put lineshapes of resonances into neural network, but also know every lineshape in the training data standing for a molecular state or an elementary state. So, criteria for the nature of resonances are needed. The pole counting rule (PCR) [18] is a convenient choice. According to PCR, the nature of an S-wave resonance is connected to the number of poles near threshold. If there is only one pole near the threshold of the couple channel, the resonance is a molecular state and if there are two poles near the threshold, the resonance is elementary. This method has been used in many works to study the nature of exotic hadron states [19][20][21][22][23]. Then, it demands one to choose a parametrization for resonance amplitudes, in which the parameters can be adjusted easily and then the values of parameters will make difference in pole positions as well as the lineshapes. The Flatté-like parametrization is chosen in this work to generate training data. If two final states are considered, the parametrizations of invariant mass spectra are parametrized as, (1) where M is the bare mass for the resonance and bg is coherent background contributions which takes the formula of first order polynomial of √ s. φ 1,2 are coherent angles. ρ 1 (s), ρ 2 (s) represent two-body phase space factor for final state 1 (FS1) and final state 2 (FS2), respectively. The FS1 stands for lower couple channel whose threshold is much lower than the resonance mass while FS2 stands for the couple channel whose threshold is near the resonance. It means that if the resonance is a hadronic molecule, it should be the composite of the two hadrons in FS2. As an example, for Z c (3900), the FS1 stands for J/ψπ channel and FS2 stands forD * D channel.
In practice, the coupling constant g 1 , g 2 will be adjusted and then lineshapes as well as the pole positions in the complex s plane will change. Different sheets of complex s plane are defined by changing signs of phase space factors, as shown in Tab. I. Generally speaking, an elementary state means it hardly couples to FS1 nor FS2, so one should keep g 1 , g 2 small to ensure there are two poles near the threshold of FS2. A molecular state indicates that it strongly couples to the second channels, so one should make g 2 g 1 to ensure that only one pole is near the threshold of FS2. According to PCR, if the distance between two poles is larger than ∼ 200 MeV, one can say there is only one nearby pole [24]. Since the parameter are determined, the training data can be generated using Eq. (1). Firstly, one should choose a proper energy range and it can cover the widths of most exotic hadron states. In this paper, the data of FS1 are from (E th2 − 95) MeV to (E th2 + 95) MeV and the data of FS2 are from E th2 to (E th2 +100) MeV, where E thi , i = 1, 2 are the threshold energies of FS1 or FS2, respectively. It is stressed that the absolute value of E th is unimportant since only the lineshapes will be sent into neural network. The energy resolution is fixed at 1 MeV so that all the samples have the same size. In order to simulate experimental data, the effects of error bar are taken into consideration. The values at every energy point calculated from Eq. (1) will be taken as the average value of Gaussian distribution, and 10 percent of their values will be taken as the standard deviations. Then the Gaussian sampling with these average values and standard deviations gives the data in which the error bar effects are taken into account. At last, the normalization: makes points in every sample ranging from 0 to 1, where Ω training is the set of points in one sample.
Using the methods introduced above, 4 × 10 4 groups of labeled invariant mass spectra are generated in total, including 2 × 10 4 elementary state with label "1" and 2 × 10 4 molecular states with label "0". Each group consists of one invariant mass spectrum in FS1 and one invariant mass spectrum in FS2. In practice, 2.4 × 10 4 groups of invariant mass spectra are used to train the neural network and 1.6 × 10 4 groups are used to do validation test. In both training set and testing set, the number of elementary states and molecular states are the same.

III. CONSTRUCTION OF NEURAL NETWORK MODEL
The goal of the neural network is to establish a map between the input space of invariant mass spectra and the output elementary-molecular classification space. A typical neural network consists of an input layer, some hidden layers and an output layer. In this work, the units in the input layer are numerical values of invariant mass spectra in two final states. Once the inputs are given, they will be linearly transformed with some weights and biases. Then the transformed values will be sent into the activation function, which takes the form of ReLU (rectified linear unit), for units in hidden layers, and for units in output layer, it takes the form of sigmoid, The values passed through activation functions are fed to the units in next layer until meet the output layer, i.e., every unit value is obtained by combining all units in the previous layer linearly and passing through the activation function. The reason for choosing Eq. (4) as the activation function for output layer is that it can map arbitary value into a number ranging from 0 to 1, which represents the possibility of being an elementary state.
Once the output is given, which means the network gives the prediction of the nature of input sample, the discrepancy between the label and the prediction will be calculated using loss function. In this problem, the loss function takes the form of, where n is the batch size which means the number of samples in each input and p i is the output for ith sample, q i is the label of ith sample in the training batch. During the training, the weights and biases used to do the linear combination of unit values should be optimized to have more accurate predictions or smaller loss. Usually, the upgrading of weights and biases use gradient descent algorithm which has been integrated into many packages. The optimzer Adam [36] is used in this work. The process from the input to output called forward pass while the upgrading of weights and biases called backpropagation. With circulation of these two processes, the classifier will be trained accurate enough. All these algorithm are performed with PyTorch [37] and more details about neural network can be found in, e.g., Ref. [25].
Generally speaking, more hidden layers or hidden layers with more units means there are more parameters including weights and biases involved in the network, so that it can give more accurate predictions with respect to the training data. But on the other hand, if there are too many parameters, the problem of overfitting will arise which means that the network has poor predictive power when it is applied to new data. So, there exist a balance between accuracy and extrapolation. In this work, neural networks with different number of hidden layers and different number of units in these layers are constructed.
The process of a sample crossing the neural network is displayed in Fig. 1. The sample consists of two invariant mass spectrum and is represented by a vector. The weight matrix and the bias vector are donated as W, b, respectively. In Fig. 1, the subscription is the shape of this matrix where n i is the number of units in the i-th hidden-layer (i=1,2,...,m). The superscription (i) represents that this matrix belongs to the i-th hidden-layer and the superscription (o) means that this matrix belongs to the output layer. The output o is a number ranging from 0 to 1.

A. Training and validation tests
The structures of neural network models are listed in Tab. II. As mentioned in Sec. II, the performances during training are monitored by validation tests. The data in validation test are different from those in training, but they are generated by the same method. Fig. 2 shows the loss of training data as well as test data of different models. It is shown that after 500 training epochs, the losses in both training data as well as test data become small. The training loss usually becomes smaller with more training epochs, but it does not mean that the network will perform better and better with as many training epochs as one wants. If the test loss arrives a small value but become larger with more training epoch, then the neural network is thought to be well trained and more training epochs will cause overfitting, as shown in Fig. 2. It looks like converging to a local minimum of parameter space in conventional fitting. So, the network with the minimum test loss will be picked up. Further more, the output distributions of samples in test data are shown in Fig. 3.

B. Applications
After training and model selection, the trained models is applied to experimental data of X(3872), X(4260) as well as Z c (3900). For Z c (3900), the data observed in the J/ψπ andD * D final states from Ref. [11] and Ref. [12], respectively are used. As for X(3872), the data observed in the J/ψππ andD * D final states from Ref. [14] and Ref. [10] are used. For X(4260), the data observed in J/ψππ [16] and χ c0 ω [17] final states are used.
Before sending the experimental data into the neural network, it is noted that the energy resolution are different among experiments. So, these data should be supplied in order to have the same size as input of the neural network in which the energy resolution is fixed at 1 MeV. To be specific, if the energy resolution of experimental data is larger than 1 MeV (as shown in Fig. 4 (a), (b) and (d)), then the method of linear interpolation will be employed to supply extra points. For example, the energy resolution of Fig. 4 (a) is 10 MeV, so 9 points will be inserted using linear interpolation between two neighbouring experimental data. Of course, the effects of error bar are also taken into consideration, the final application data are obtained by Gaussian sampling in which the central values of experimental data are regarded as average values and the error bars are regarded as standard deviations. At last, 100 invariant mass spectra are obtained for every group of experimental data using Gaussian sampling and the energy resolution is 1 MeV after linear interpolation. After these procedures, the application data should be normalized into 0 to 1 as done in Eq. (2).
Then, the application data are substituted into trained models to obtain the predictions for the nature of Z c (3900), X(4260) and X(3872). The predictions for these three states are shown in Fig. 5. For Z c (3900), about all 100 samples in the application data have outputs near 0 in all four models, means that it strongly couples toD * D channel and should be regarded as a hadronic molecule ofD * D. This interpretation is consistent with Refs. [21,[26][27][28]. For X(3872), almost all outputs of 100 samples are near 1, that means the coupling between X(3872) andD * D channel is not strong enough to dominate its production. In other words, X(3872) behaves more like an elementary state. The neural network prediction for X(3872) is consistent with e.g., Refs. [19,[29][30][31]. The predictions for X(4260) are more coincident. The outputs from all four models are near 1, which means that X(4260) couples to ωχ c0 not strong enough and can not be regarded as a molecule of ωχ c0 . This prediction is consistent with e.g., Refs. [32,33]. Furthermore, the measurement of cross section e + e − → µ + µ − performed by BESIII [34] gives the muonic width of X(4260) to be from 1.09 to 1.53 KeV which strongly indicates that X(4260) has a charmonium nature [33].
Of course, many exotic hadron states have data only in one final state (usually only in FS1). At this time, one still can identify whether they are hadronic molecules by neural network. As an example, a network is trained by FS1 data in Eq. (1) alone. This network is used to identify the properties of X 1 (2900), which is regarded as aD 1 K molecule [23]. When the experimental data are sent into network, about a half samples are identified as elementary states and another are molecular states, which means the network is not enough yet to determine its nature. But this may caused by the big error bars of the experimental data (see Fig.6(a)). When generating samples by Gaussian sampling, their lineshapes can change a lot within error bars which are used as standard deviations. So, one can reduce the standard deviations in Gaussian sampling to avoid misleading. As shown in Fig. 6(b), with smaller standard deviations, the outputs are converged to a molecular state, i.e., scores of the most samples are < 0.5. That means, neural network can also identify hadronic molecules using only one invariant mass spectrum if the experimental error bars are small enough.

V. SUMMARY AND OUTLOOK
In this work, machine learning models based on neural network are developed and used to decide whether X(3872), X(4260) as well as Z c (3900) can be regarded as hadronic molecules of special channels. Resonance lineshape data observed in two final states are taken as input, and the output is a number which can be seen as the possibility of this resonance to be an elementary state. In other words, the resonance is more like a hadronic molecule of particles in FS2 whose threshold is near the resonance if the output is closer to 0. The well trained networks are picked up to do these tasks. The results show that Z c (3900) can be regarded a molecule ofD * D but X(3872) is not like a molecule ofD * D. Besides, if data of X(4260) in J/ψππ and ωχ c0 channels are taken as input, the predictions of neural networks suggest that it is not like a molecular state of ωχ c0 . These interpretations from neural networks are consistent with many previous phenomenological studies. So, the method employed in this work is practicable.
The philosophy of a neural network to do such a binary classification task is to represent the inputs as points in the parameter space and build a hypersurface to separate these points into two classes. The training process is nothing but to adjust the hypersurface to do the separation better. In other words, the output for a set of data reflects that how can the parameters take values to describe such an input. That sounds just like to do a conventional fit. In this work, the data in both two final states perform better compared to data in only one final state, it makes the classification like a joint fitting and the prediction are more reliable.
The method in this work provides a new way to identify the nature of hadron states. It is a data driven, model independent and general method. That means a trained model can be used to classify different resonances in different processes. In the future, the machine learning method can be developed better to understand the nature of exotic states better. For instance, Dalitz plots are more original data compared with invariant mass spectra and convolution neural network is powerful in image recognition. So, it can be used to classify elementary states and hadronic molecules based on Dalitz plots. It is belived that with the development of algorithms and computer technology, machine learning will become more and more popular in the field of particle physics.