Interactions of $B_{c}$ Meson in Relativistic HeavyIon Collisions
Abstract
We calculate the dissociation crosssections of mesons by and mesons including anomalous processes using an effective hadronic Lagrangian. The enhancement of production is expected due to QGP formation in heavyion experiments. However it is also expected that the production rate of meson can be affected due to the interaction with comovers. These processes are relevant for the experiments at RHIC. Thermal average crosssections of are evaluated with a form factor when a cut off parameter in it is 1 and 2 GeV. Using these thermal average crosssections in the kinetic equation we study the time evolution of mesons due to dissociation in the hadronic matter formed at RHIC.
Keywords: Relativistic heavy ion collisions, MesonMeson interaction, QGP.
PACS number(s): 25.75.q, 13.75.Lb, 14.40.Nd
1 Introduction
In 1986 Matsui and Satz [1] hypothesized that in a deconfined medium color screening would have dissociated the , resulting in a suppressed yield of . This deconfined state is called QuarkGluon Plasma (QGP). Thus for the existence of QGP, suppression of could be considered as a probe. Anomalously large suppression of events was observed by NA50 experiment at CERN [2] with moderate to large transfer energy from the Pb+Pb collision at GeV/c. However, this observed suppression may also occur due to absorption by comoving hadrons mainly and , especially if the dissociation cross section is at least few mb [3, 4, 5, 6, 7, 8]. To calculate these cross sections, quark potential models, perturbative QCD [9], QCD sumrule approach [10, 11] and flavor symmetric effective Lagrangian [12, 13, 14, 15] has been used.
Analogous to charmonium, suppression of bottomonium states is also predicted during the formation of QGP [1]. Recently it was observed by CMS in Pb+Pb collisions that excited states of bottomonium are strongly suppressed [16]. To have unambiguous interpretation of the the observed signal, the information of dissociation cross section is also needed [13, 17]. It was suggested that the production rate of heavy mixed flavor hadrons would be affected in the presence of QGP [18, 19]. For calculating the rate of production of these hadrons comprehensive information is required to distinguish QGP affected hadron production and suppression due to dissociation by comovers. It is expected that production could be enhanced in the presence of QGP [19, 20]. QGP contains many unpaired and quarks due to color Debye screening. These unpaired and quarks upon encounter could form or mesons and due to relatively large binding energy, mesons probably survive in QGP [20]. However, observed production rate would also depend upon the dissociation cross section by hadronic comovers.
In Ref. [20] absorption by nucleons was examined with the mesonbaryon exchange model. The calculated cross sections were in the range of a few millibarn. Recently in Ref. [21], using the same couplings and hadronic Lagrangian within meson exchange model the dissociation of meson by meson were examined. The range of the resulting cross sections involving the form factors were mb and mb for the processes and , respectively. In Ref. [22], the dissociation of meson by mesons were examined. For the processes and the resultant cross sections with the form factor were in the range of and mb, respectively .
In this paper we investigate the dissociation by and mesons including anomalous couplings like PVV, PPPV and PVVV which were ignored in the previous studies. Inclusion of these couplings results in opening of new dissociation channels and addition of new processes and extra diagrams. The contribution of anomalous couplings is found to be significant for calculating cross sections of charmonium dissociation with and meson in Ref. [23], mesons in Ref. [24] and dissociation of meson by nucleons in Ref. [25]. We also calculate the thermal average cross sections and study the time evolution of meson at RHIC using a schematic expanding fireball model with an initial abundance determined by the statistical model.
The paper is organized as follows. In Sec. 2, the interaction Lagrangian terms which are relevant for the description of the dissociation of by and mesons including anomalous processes are given and also analytical expressions of the amplitudes for the dissociation of meson are reported. In Sec. 3, we calculate the cross sections with and without form factor and thermal average cross sections. In Sec. 4, we study time evolution of the meson abundance at RHIC in a schematic model. In Sec. 5, we present the summary and discussion.
2 Interaction Lagrangian and Amplitudes of meson dissociation
2.1 Interaction Lagrangian
We consider the following reactions using an effective Hadronic Lagrangian.

(1) 
The processes in the first and second column, and also in the third and fourth column have same cross sections as being charge conjugation of each other. The generic form for the 1st reaction is given as

(2) 
For calculating the cross sections of the above reactions, relevant interaction Lagrangian terms are required. The required interaction Lagrangian for normal processes (for which the relevant couplings are dimensionless) are obtained using the method described in Refs. [21, 22] and are given as follows.
(3a)  
(3b)  
(3c)  
(3d)  
(3e)  
(3f)  
(3g)  
(3j)  
(3k) 
In addition to the above normal terms there are anomalous terms as well which are required to give a complete description of the hadronic processes. The required interaction Lagrangian for the anomalous processes (for which the relevant couplings are not dimensionless) are obtained using the method described in Ref. [23] and are given as follows.
(4a)  
(4b)  
(4c)  
(4d)  
(4e)  
(4f)  
(4g)  
(4h)  
In Eqs. (3) and (4) represents Pauli spin matrices, and and represent isospin triplets,
while vector and pseudoscalar charm and bottom meson doublets are given as

2.2 Amplitudes for meson dissociation
For calculating the cross section for meson dissociation by and mesons, we use the effective Lagrangian given in Eqs. (3) and (4). In this paper we are only reporting the scattering amplitudes of anomalous processes and of additional diagrams which are dependent on the anomalous couplings. Absorption amplitudes of other diagrams which depend only on normal couplings are given in Refs. [21, 22]. Diagrams of the process are shown in Fig. 1 (2a to 2c) and the amplitudes of the diagrams are
(5c)  
And the full amplitude is written as  
(5d) 
Diagrams of the process are shown in Fig. 1 (3a to 3c) and the amplitudes of the diagrams are
(6c)  
And the full amplitude is written as  
(6d) 
Diagrams of the process are shown in Fig. 1 (4a to 4e). The amplitudes of diagram 4d and 4e which depend on anomalous couplings are
And the full amplitude is written as  
(7c) 
Now we report the absorption amplitudes of the anomalous processes of by . Diagrams of the process are shown in Fig. 2 (5a to 5c). The amplitudes of these diagrams are
(8c)  
And the full amplitude is written as  
(8d) 
Diagrams of the process are shown in Fig. 2 (6a to 6d). The amplitudes of the anomalous diagram 6d is given as
And the full amplitude is written as  
(9b) 
Diagrams of the process are shown in Fig. 2 (7a to 7d). The amplitude of the anomalous diagram 7d is given as
And the full amplitude is written as  
(10b) 
Diagrams of the process are shown in Fig. 2 (8a to 8e), the amplitudes of these diagrams are
(11b)  
(11e)  
And the full amplitude is written as  
(11f) 
We define the fourmomenta of the incoming particles as and and those of the final particles as and , which then defines and . Here , , and represent the , , and mesons masses, respectively. The polarization vector of a vector meson with momentum is represented by . After averaging (summing) over initial (final) spins and including isospin factor, we calculate the cross sections by using the total amplitudes specified in above equations. The isospin factor for calculating these cross section is 2 for all the processes.
3 Dissociation CrossSections of Meson
3.1 Numerical values of input parameters
Coupling constant  Value  Method of Derivation 

9.08  Heavy quark symmetries  
2.34  Heavy quark symmetries  
6.134  Heavy quark symmetries  
21.56  symmetry  
21.56  symmetry  
2.82  VMD  
2.58  Heavy quark symmetries  
21.56  symmetry  
67  symmetry 
Numerical values of all the meson masses are taken from Particle Data Group [26]. Estimation of the coupling constants of effective Lagrangian is required for calculating the cross sections. To fix the couplings for the normal processes, we follow the methods of Refs. [12, 27]; we refer to Ref. [12] for details. In a similar way we have determined the couplings for the anomalous interactions, which are reported in this paper whereas normal couplings are given in Refs. [21, 22]. The coupling which has a dimension of is fixed by applying the heavy quark spin symmetry. We follow Ref. [23] in which this coupling is given as
(12) 
where represents the average mass of and .
For couplings, we can apply the VMD (Vector Meson Dominance) model [12] to the radiative decays of into , i.e., . We use the same method as in ref. [23]; this leads to