Polarization and quantum entanglement effects in $B^\pm_c\to J/\psi+\pi^\pm +\pi^0$ process

Motivated by the very recent observation of the $B^+_c\to J/\psi+\pi^+ +\pi^0$ decay using proton-proton collision data by the LHCb collaboration, we study the four-body angular distributions and the quantum entanglement effects in the $B^+_c\to J/\psi+\pi^+ +\pi^0$ associated with $J/\psi\to \mu^++\mu^-$. The helicity angular distributions are given in the QCD effective theory and the von Neumann entropy is obtained in $B^\pm_c\to J/\psi(\to \mu^+\mu^-)+\rho^\pm(\to \pi^\pm \pi^0)$ decay process.


INTRODUCTION
Beauty-charm mesons composed of two different heavy flavor quarks are yet to be fully explored compared with charmonia and bottomonia which are composed of two identical heavy flavors.Up to now, only three members of beauty-charm meson family, i.e., B c , B c (2S) and B * c (2S) are discovered in particle physics experiments [1][2][3][4].Other beauty-charm states are still missing and waiting for discovery in experiments [5][6][7][8][9].
The beauty-charm mesons below the BD threshold major decay into the ground B c state via strong or electromagnetical interactions, while the ground B c state has to weak decay.Therein the transition of the ground B c state into J/ψ is very important because the J/ψ → µ + + µ − decay has clear signal at hadron collider experiments and it sheds light on the rare B c decays and leads to the discoveries of tens of rare B c decay channels at LHCb.
Very recently, the observation of B c decays into J/ψ associated with two Pion mesons has been reported for the first time at the LHCb experiment using protonproton collision data at centre-of-mass energies of 7, 8, and 13TeV corresponding to an integrated luminosity of 9f b −1 .Though the reconstruction efficiency for this channel is low, over one thousand events are yielded using LHCb run-I and run-II data.The measured ratio of the branching fractions of the In history, the first exclusive decay B + c → J/ψ + π + was discovered until 2005 by the CDF collaboration [11], while the multi-body exclusive decays B + c → J/ψ + π + + π + + π − and B + c → J/ψ + 3π + + 2π − were discovered by the CMS and LHCb collaborations in 2014 [12,13], respectively.Then the B c → J/ψ + nπ decay with Pion meson numbers n = 1, 2, 3, 5 are all observed in particle physic experiments [14].The discovery of B c decays into J/ψ associated with two Pion mesons with relative large ratio shall indicate the experimental feasibility of B c decays into J/ψ associated with four Pion mesons.
On the other hand, the phenomenon of quantum entanglement and the violation of Bell inequality are widely studied in quantum mechanics via electromagnetic forces.However, there are few studies of quantum entanglement and the violation of Bell inequality via strong or/and weak interactions which would tell us the essential behavior of the Quantum Chromodynamics (QCD) field theory and the standard electroweak unified theory of particle physics.The quantum entanglement effects in hadron decays are studied in literatures [15][16][17][18][19][20][21].
In this paper, we study the polarization and quantum entanglement effects in B ± c → J/ψ + π ± + π 0 process.First, we analyze the invariant mass distribution of two Pion mesons in the B + c → J/ψ + π + + π 0 decay and find that the mediate vector rho meson dominates the decay of B + c → J/ψ + π + + π 0 .Next we investigate the helicity angle distribution in B ± c → J/ψ(→ µ + µ − ) + ρ ± (→ π ± π 0 ) decay based on the Jacob-Wick theory.Using the decay amplitudes in nonrelativistic Quantum Chromodynamics (NRQCD) effective theory, the helicity angle distributions are plotted.These angle distributions can be employed to test the QCD effective theory and the fundamental properties of the initial beauty-charm mesons.Then we study the quantum spin entanglement of two vector mesons J/ψ and ρ in B ± c → J/ψ + π ± + π 0 process.The von Neumann entropy to describe the degree of entanglement is derived in the end.

ANGULAR DISTRIBUTIONS
The invariant mass distribution of two Pion mesons in B + c → J/ψ + π + + π 0 have been studied previously in literatures [22][23][24] before the experimental discovery of this channel [10].Using the method in Ref. [24] but with different parametrization form, the spectral function ρ 2π T (m 2π ) can de written as ( Refitting the LHCb data, the parameters can be determined as a = 8.25 × 10 4 GeV , b = 0.37GeV −2 , c = 0.158GeV and m 1 = 0.7518GeV with χ 2 /dof = 0.96.From the curve, the ρ meson plays the dominant role in the B + c → J/ψ + π + + π 0 decay.The mass and decay width of ρ can be obtained through the parameter c and m 1 .The fitting curve in comparison with the LHCb data is plotted in Fig. 1.FIG. 1: Invariant mass distribution of π + π 0 in the B + c → J/ψ + π + + π 0 decay, where the LHCb data is from Ref. [10].
Thus we can study the dominant two-body B + c → J/ψ + ρ + decay with ρ + → π + + π 0 in place of B + c → J/ψ + π + + π 0 decay.The pure three-body decays without mediate resonances in B + c → J/ψ +π + +π 0 are small and ignored.The decay amplitude of B c → J/ψ + ρ can be written as where the helicities satisfy with λ ∈ (±1, 0) due to the conservation of momentum and angular momentum in this process.m i , p i and ϵ i are the mass, momentum and polarization vector for ρ and J/ψ, respectively.a, b and c can be calculated in certain theoretical frameworks.The weak effective Hamiltonian H ef f governing B c → J/ψ + ρ can be written as where V ud and V cb are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements.Q 1,2 (µ) are the effective fourfermion operators while C 1,2 (µ) are the corresponding Wilson coefficients [25].
After the decay B c → J/ψ + ρ, the ρ further decays into two Pion mesons almost 100% while J/ψ can further decay into µ + µ − with the branching ratio 5.96% [14].All the decay chains are depicted in Fig. 2, where the helicity angles θ 1 , θ 2 and ϕ are defined.According to the Jacob-Wick theory [26], the angular distributions for B c → J/ψ(→ µ + µ − ) + ρ(→ ππ) can be written as where H λλ ′ = h λ h * λ ′ is the density matrix elements for the B c → J/ψ + ρ decay.M is the mass of initial B c meson.The parameter p m is written as FIG. 2: Helicity angles in the sequential decay process B ± c → J/ψ +ρ ± → (µ + µ − )(π ± π 0 ), where θ1 is the Pion meson angle in the rest of ρ while θ2 is the muon angle in the rest of J/ψ.ϕ is the angle between the two subprocess decay planes.
The relationship between the parameters a,b,c and the helicity amplitude can be expressed as where In naive factorization based on the NRQCD effective theory, the amplitude parameter a, b, c can be expressed by form factors [25,[27][28][29][30] Then the values of the three helicity amplitudes can be obtained as Using the helicity amplitudes, we plot the normalized angular distributions in Fig. 3, Fig. 4 and Fig. 5. From these results, the distribution curves for helicity angles θ 1 and θ 2 are different because the two vector mesons decay to final states with different spin structures.For the secondary ρ + → π + + π 0 decay in B + c → J/ψ + ρ + process, the maximum probability of helicity angle θ 1 is around π/4, while the secondary J/ψ → µ + + µ − decay in B + c → J/ψ + ρ + process, the maximum probability of helicity angle θ 2 is around π/2.The secondary decays of two vector mesons refer to occur in the identical plane.

QUANTUM ENTANGLEMENT AND VON NEUMANN ENTROPY
If the polarization state of the two vector meson system for J/ψ and ρ from B c is pure, the wave function in the rest frame of B c can be written as [31] where the +1, − 1, 0 in the bracket of the vector mesons denote the helicity values.The normalized constant H is defined as Then it is easy to write down the polarization density matrix ϱ = |Ψ⟩⟨Ψ| in terms of the polarization amplitude These helicity amplitudes can be calculated according to certain theoretical frames presented in this paper.
In order to study the entanglement of two polarized vector particles with mass, the von Neumann entropy is a good choice to quantify the entanglement.The von Neumann entropy is defined as [32] ϱ A and ϱ B are reduced density matrices of two subsystems, which can be obtained by tracing one of the subsystems by the density operator.In this article, the two subsystems are J/ψ and ρ.For a two-qutrit system, the von Neumann entropy satisfies 0 ≤ ε ≤ √ 3.When the wave function of the two systems can be written as a direct product state, the minimum value is taken.When the entanglement of the two systems is maximum, the maximum value is taken.It is easy to figure out that the reduced density matrices of two sub systems satisfy Combined with Eqs.(9)(10)(11), the von Neumann entropy of subsystems J/ψ and ρ can be obtained The result of von Neumann entropy proved that there is quantum spin entanglement between two daughter vector particles in the B c decay.

CONCLUSION
In this paper, the angular distributions of B ± c → J/ψ(→ µ + µ − )+ρ ± (→ π ± π 0 ) process are given in the helicity form of Jacob-Wick theory and the quantum spin entanglement effects among vector mesons are studied.To describe the degree of quantum spin entanglement in this process, the von Neumann entropy is calculated.Through the measurement of the invariant mass distribution of two Pion mesons, the secondary decay angular distributions in ρ ± → π ± + π 0 and J/ψ → µ + + µ − , and the distribution with respect to the angle between the two secondary decay planes, the nature of the initial beauty-charm meson can be determined.On the other hand, these measurements of angular distributions shall provide new views to hunting for the missing members of beauty-charm meson family.The study of the quantum spin entanglement effects in hadron cascade decays will be useful for us to understand the quantum correlations and non-locality in standard field theory of particle physics.