Next Article in Journal
The Optogalvanic Spectrum of Neutral Lanthanum between 5610 and 6110 Å
Previous Article in Journal
A Laboratory Astrophysics Problem: The Lifetime of Very Long-Lived Levels in Low-Charge Ions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Branching Ratio and CP Violation in B s 0 ¯ ϕϕ Decay in the Framework of QCD Factorization

1
Department of Physics and Astronomy, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
2
Physics Discipline, Science, Engineering and Technology School, Khulna University, Khulna 9208, Bangladesh
*
Author to whom correspondence should be addressed.
Atoms 2020, 8(2), 22; https://doi.org/10.3390/atoms8020022
Submission received: 16 March 2020 / Revised: 29 April 2020 / Accepted: 14 May 2020 / Published: 16 May 2020

Abstract

:
We present an analytic calculation of Branching Ratio (BR) and Charge-Parity (CP) violating asymmetries of the B s 0 ¯ meson decay into the two light vectors ϕ ϕ . In doing this we calculate the helicity amplitude of the present decay in the framework of QCD factorization approach. We find the BR of B s 0 ¯ ϕ ϕ = ( 1.56 ± 0.23 ) × 10 5 . We also calculate the direct CP violation, CP violation in mixing and CP violation due to interference which are A C P dir = 0.00355 ± 0.00152 , A C P mix = 0.00629 ± 0.03119 and A C P Δ Γ = 0.99997 ± 0.00019 , respectively. Our results are in agreement with the recent theoretical predictions and experimental measurements.

1. Introduction

The Standard Model (SM) of Particle Physics describes the fundamental building blocks of matter and their interactions. One of the deficiencies in the SM is that it does not accommodate the matter-antimatter asymmetry of the universe. One of three proposed conditions to explain this asymmetry is the violation of Charge Conjugation-Parity (CP) symmetry [1]. In 1973, Kobayashi and Maskawa (KM) proposed an explanation of CP violation in the B meson and in doing so, predicted the existence of the third generation of quarks [2]. In recent years, sizable CP violation has been observed in the decay of B d meson [3,4,5]. The charmless two-body non-leptonic B s decay is another important choice in exploring CP violation. At the B-factories [6,7,8,9,10] ( Belle, BaBar, CDF and LHCb) both B d - B ¯ d and B s - B ¯ s systems are produced. Both systems exhibit the particle-antiparticle mixing phenomenon. Because of the B s - B ¯ s system has higher mass difference ( Δ m s = 17.69 ± 0.08 ps 1 ) than that of B d - B ¯ d system ( Δ m d = 0.510 ± 0.004 ps 1 ), the B s meson has faster oscillations than B d meson [11]. The study of the CP violation in the B s system with the benefit of faster oscillation offers an excellent opportunity to detect the possible deviations from SM predictions and may lead to a new physics beyond the SM. Many authors have studied the decays of B s P P , P V , V V [3,5,12,13,14,15,16,17,18,19,20,21,22,23,24], with V being a light vector meson and P being a pseudoscalar meson. The decay of B s V V reveals more dynamics than B s P V or B s P P .
In order to get a clear idea of CP violation, one needs to know the exact BR of the decay modes which motivates us to make an analytic calculation of the B R ( B s 0 ¯ ϕ ϕ ) . The analytic calculation of the BR of B s V V decays is achieved by many approaches as Quantum Chromodynamics Factorization (QCDF) [5,18], the Perturbative QCD (PQCD) [20,21], the Soft-Collinear Effective Theory (SCET) [22,23] and Factorization-Assisted Topological amplitude (FAT) [24].
This paper focuses on the calculation of the BR and CP violation of the B s 0 ¯ ϕ ϕ decay in the framework of QCDF approach. In 2003, X. Li et al. [12] studied B s 0 ¯ ϕ ϕ decay and predicted B R ( B s 0 ¯ ϕ ϕ ) as 1.79 × 10 5 and 3.68 × 10 5 within the naive factorization (NF) and QCDF approaches, respectively. The present decay mode can be used as normalization for the studies of some other channels of charmless B s 0 meson decays [10]. The first observation of the B R ( B s 0 ¯ ϕ ϕ ) = 1.4 0.5 + 0.6 ( stat . ) ± 0.6 ( syst . ) × 10 5 was performed by CDF in 2005 [8]. Later on, they updated their result to 2.32 ± 0.18 ( stat . ) ± 0.82 ( syst . ) × 10 5 in 2011 [9]. Recently, LHCb collaboration reported that B R ( B s 0 ¯ ϕ ϕ ) = 1.84 ± 0.05 ( stat . ) ± 0.07 ( syst . ) ± 0.11 ( f s / f d ) × 10 5 [10]. A recent theoretical calculation was performed by Yan et al. [20] and they found that B R ( B s 0 ¯ ϕ ϕ ) = 1.88 0.38 + 0.49 × 10 5 with a A C P dir ( % ) = 0.7 ± 0.2 .
In this paper, we report the calculation of the BR and CP violation of B R ( B s 0 ¯ ϕ ϕ ) within QCDF using Mathematica packages.
Since B R ( B s 0 ¯ ϕ ϕ ) is a vector decay, we need to find the helicity amplitude to calculate the BR. We formulate the helicity amplitude with neglecting both the annihilation and chiral contributions since they are negligibly small. The calculated B R ( B s 0 ¯ ϕ ϕ ) is ( 1.56 ± 0.23 ) × 10 5 . Furthermore, we calculated the CP violation in the SM for B s 0 ¯ ϕ ϕ decay via pure penguin diagram within QCDF. We find that direct CP violation ( A C P dir ), CP violation in mixing ( A C P mix ) and CP violation due to interference ( A C P Δ Γ ) are 0.00355 ± 0.00152 , 0.00629 ± 0.03119 and 0.99997 ± 0.00019 , respectively. We find that the reported results are consistent with other available predictions and experimental observation [5,8,9,10,18,20,21,23,24,25].

2. Theoretical Framework

2.1. The Effective Hamiltonian

The effective Hamiltonian ( H eff ) describing the transition amplitude of an initial state to a final state follows the Fermi’s Golden Rule. In terms of the effective Hamiltonian, the BR of B s 0 ¯ decays to two vector mesons can be written as [12,16]:
B R ( B s 0 ¯ V 1 V 2 ) = τ B s 0 p c s 8 π m B s 0 2 V 1 ( h 1 ) V 2 ( h 2 ) | H eff | B s 0 ¯ 2 ,
where h 1 , h 2 are the helicities of the final-state vector mesons V 1 and V 2 with four-momentum p 1 and p 2 , respectively, the m B s 0 and τ B s 0 are the mass and lifetime of B s 0 ¯ meson, the statistical factor s = 1 2 , 1 for two identical and different meson final states, respectively. In the rest frame of B s 0 ¯ system, since B s 0 ¯ meson has spin 0, we have h 1 = h 2 = h and p c = | p 1 | = | p 2 | is the momentum of either of the two outgoing vector mesons.
Only the penguin operators can contribute to the B s 0 ¯ ϕ ϕ decay channel which is b s transition, so the relevant H eff can be written as [12]
H eff = G F 2 λ t ( s ) i = 3 10 C i ( μ ) O i ( μ ) ,
where G F is the Fermi constant and λ t ( s ) = V t b V t s * is the CKM factor. The O i ( μ ) are the local four-fermion operators (Wilson operators). The C i ( μ ) are the effective Wilson coefficients which have been reliably evaluated to the next-to-leading logarithmic order (NLL) with μ m b being the renormalization scale. The O 3 , , 6 and O 7 , , 10 are the QCD and electroweak penguin operators, respectively and can be expressed as follows [12,26]:
O 3 = ( s ¯ b ) V A q ( q ¯ q ) V A , O 4 = ( s ¯ i b j ) V A q ( q ¯ j q i ) V A , O 5 = ( s ¯ b ) V A q ( q ¯ q ) V + A , O 6 = ( s ¯ i b j ) V A q ( q ¯ j q i ) V + A , O 7 = 3 2 ( s ¯ b ) V A q e q ( q ¯ q ) V + A , O 8 = 3 2 ( s ¯ i b j ) V A q e q ( q ¯ j q i ) V + A , O 9 = 3 2 ( s ¯ b ) V A q e q ( q ¯ q ) V A , O 10 = 3 2 ( s ¯ i b j ) V A q e q ( q ¯ j q i ) V A ,
where i and j are the color indices, q denotes all the active quarks at the scale μ m b , i.e., q u , d , s , c , b and e q are the corresponding quark charges. Also, the ( q ¯ 1 q 2 ) V ± A = q ¯ 1 γ μ ( 1 ± γ 5 ) q 2 are the right and left-handed vector-axial currents, respectively.

2.2. The Factorizable Amplitude for B s 0 ¯ V 1 V 2

The BR for B s 0 ¯ V 1 V 2 decays can be calculated by inserting Equation (2) into Equation (1). The calculation of the resulting hadronic matrix elements of the local four fermion operators i.e., V 1 ( h ) V 2 ( h ) | O i | B s 0 ¯ represents a theoretical challenge. In order to solve this problem, naive factorization (NF) in which the hadronic matrix elements is replaced by the product of the matrix elements of two currents is carried out as follows [12,16]:
V 1 V 2 | ( q ¯ 2 q 3 ) V ± A ( q ¯ 1 b ) V A | B s 0 ¯ V 2 | ( q ¯ 2 q 3 ) V ± A | 0 V 1 | ( q ¯ 1 b ) V A | B s 0 ¯ ,
with the vector meson V 2 being factored out. The quark content of the two vectors in the above equation are ( q 1 q 4 ¯ ) V 1 and ( q 2 q 3 ¯ ) V 2 . The Equation (4) represents factorizable amplitude, χ ( B s 0 ¯ V 1 , V 2 ) which is expressed by
χ ( B s 0 ¯ V 1 , V 2 ) = V 2 | ( q ¯ 2 q 3 ) V ± A | 0 V 1 | ( q ¯ 1 b ) V A | B s 0 ¯ .
The massive vector of spin 1 has three z-components s z = ± 1 , 0 corresponding to three possible helicity states h = ± , 0 . So, it can exist in three possible orthogonal polarization states ε * μ ( ± ) and ε * μ ( 0 ) . These states represent two transverse polarization modes ε * μ ( ± ) corresponding to s z = ± 1 and a longitudinal one ε * μ ( 0 ) corresponding to s z = 0 [27]. For a vector meson V, let p μ , ε * μ , f V and m V be its momentum, polarization vector, decay constant and mass, respectively. Let p B s 0 and m B s 0 be the momentum and mass of the B s 0 ¯ meson. If we choose, the coordinate systems in the Jackson convention; that is, in the B s 0 ¯ rest frame, one of the vector mesons is moving along the + z -axis of the coordinate system and the other along the z -axis, while the x-axes of both daughter particles are parallel [17]:
ε 1 * μ ( 0 ) = p c , 0 , 0 , E 1 / m V 1 , ε 2 * μ ( 0 ) = p c , 0 , 0 , E 2 / m V 2 , ε 1 * μ ( ± ) = 1 2 0 , 1 , i , 0 , ε 2 * μ ( ± ) = 1 2 0 , 1 , + i , 0 , p 1 μ = ( E 1 , 0 , 0 , p c ) , p B s 0 μ = ( m B s 0 , 0 , 0 , 0 ) ,
with
p c = m B s 0 2 ( m V 1 m V 2 ) 2 m B s 0 2 ( m V 1 + m V 2 ) 2 2 m B s 0 , E = p c 2 + m 2 .
By combining Equations (6) and (7), we obtain
ε 1 * . p 1 = 0 , and ε * . p B s 0 = m B s 0 p c m V ,
ε 1 * ( 0 ) . ε 2 * ( 0 ) = m B s 0 2 m V 1 2 m V 2 2 2 m V 1 m V 2 and ε 1 * ( ± ) . ε 2 * ( ± ) = 1 ,
( ε 1 * ( 0 ) . p B s 0 ) ( ε 2 * ( 0 ) . p B s 0 ) = m B s 0 2 p c 2 m V 1 m V 2 and ( ε 1 * ( ± ) . p B s 0 ) ( ε 2 * ( ± ) . p B s 0 ) = 0 .
The first part of Equation (5) is given by [12]
V 2 ( p 2 , ε 2 * ) | q 2 γ μ q ¯ 3 | 0 = i f V 2 ε 2 μ * m V 2 ,
where V 2 | q 2 ¯ γ μ γ 5 q 3 | 0 = 0 . Whereas, the second part of Equation (5) can be written as [12,17,19]
V 1 ( p 1 , ε 1 * ) ( q ¯ 1 b ) V A B s 0 ¯ ( p B s 0 ) = ε 1 μ * ( m B s 0 + m V 1 ) A 1 B s 0 ¯ V 1 ( q 2 ) + ( p B s 0 + p 1 ) μ ( ε 1 * . p B s 0 ) A 2 B s 0 ¯ V 1 ( q 2 ) m B s 0 + m V 1 + q μ ( ε 1 * . p B s 0 ) 2 m V 1 q 2 [ A 3 B s 0 ¯ V 1 ( q 2 ) A 0 B s 0 ¯ V 1 ( q 2 ) ] i ϵ μ ν α β ε 1 * ν p B s 0 α p 1 β V B s 0 ¯ V 1 ( q 2 ) m B s 0 + m V 1 ,
where A 0 B s 0 ¯ V 1 , A 1 B s 0 ¯ V 1 , A 2 B s 0 ¯ V 1 and A 3 B s 0 ¯ V 1 are the transition form factors of the B s 0 ¯ V 1 V 2 decay via the axial current while V B s 0 ¯ V 1 is the transition form factor via the vector one. To cancel the poles at q 2 = 0 in Equation (12) we have the relation
A 3 B s 0 ¯ V 1 ( 0 ) = A 0 B s 0 ¯ V 1 ( 0 ) .
Using Equation (13) in Equation (12) we obtain
V 1 ( p 1 , ε 1 * ) ( q ¯ 1 b ) V A B s 0 ¯ ( p B s 0 ) = ε 1 μ * ( m B s 0 + m V 1 ) A 1 B s 0 ¯ V 1 ( q 2 ) + ( p B s 0 + p 1 ) μ ( ε 1 * . p B s 0 ) 2 A 2 B s 0 ¯ V 1 ( q 2 ) m B s 0 + m V 1 i ϵ μ ν α β ε 1 * ν p B s 0 α p 1 β 2 V B s 0 ¯ V 1 ( q 2 ) m B s 0 + m V 1 .
Thus, the factorizable amplitude χ ( B s 0 ¯ V 1 , V 2 ) is the product of Equations (11) and (14) which can be written as
χ h ( B s 0 ¯ V 1 , V 2 ) = i f V 2 m V 2 [ ε 1 * ( h ) . ε 2 * ( h ) ( m B s 0 + m V 1 ) A 1 B s 0 ¯ V 1 ( m V 2 2 ) ( ε 1 * ( h ) . p B s 0 ) ( ε 2 * ( h ) . p B s 0 ) 2 A 2 B s 0 ¯ V 1 ( m V 2 2 ) m B s 0 + m V 1 + i ϵ μ ν α β ε 2 * μ ( h ) ε 1 * ν ( h ) p B s 0 α p 1 β 2 V B s 0 ¯ V 1 ( m V 2 2 ) m B s 0 + m V 1 ] ,
where the momentum transfer is q = p B s 0 p 1 and the totally antisymmetric Levi-Civita tensor is normalized by ϵ 0123 = 1 . By the sum over the non-zero 24 components of the Levi-Civita tensor with the definitions in Equation (6), and for h = 0 we find
i ϵ μ ν α β ε 2 * μ ( 0 ) ε 1 * ν ( 0 ) p B α p 1 β = 0
and for h = ± the only survived terms are the two terms of { μ , ν , α , β } = { 1 ( 2 ) , 2 ( 1 ) , 0 , 3 } . Thus we can easily get
i ϵ μ ν α β ε 2 * μ ( ± ) ε 1 * ν ( ± ) p B α p 1 β = ± m B s 0 p c .
Substituting from Equations (8)–(10), (16) and (17) into Equation (15), one can get
χ 0 ( B s 0 ¯ V 1 , V 2 ) = i f V 2 b 2 m V 1 a A 1 B s 0 ¯ V 1 ( m V 2 2 ) c 2 A 2 B s 0 ¯ V 1 ( m V 2 2 ) , χ ± ( B s 0 ¯ V 1 , V 2 ) = i f V 2 m V 2 b A 1 B s 0 ¯ V 1 ( m V 2 2 ) ± c V B s 0 ¯ V 1 ( m V 2 2 ) ,
with
a = m B s 0 2 m V 1 2 m V 2 2 , b = m B s 0 + m V 1 , c = 2 m B s 0 p c b .

2.3. The Helicity Amplitude of B s 0 ¯ ϕ ϕ

In general, the B s 0 ¯ V 1 V 2 amplitude can be decomposed into three independent helicity amplitudes A 0 , A + and A corresponding to h = 0 , + and −, respectively. We use the notation
A h = V 1 ( h ) V 2 ( h ) H eff B s 0 ¯ .
Then, A h ( B s 0 ¯ ϕ ϕ ) can be written as
A h = G F 2 λ t ( s ) ϕ ϕ | i = 3 10 C i ( μ ) O i ( μ ) | B s 0 ¯ .
The dynamical details of the decay process are coded in the so-called effective parameters a i ( μ ) which are related to the Wilson coefficients C i ( μ ) through the relation
a i ( μ ) = C i ( μ ) + C i ± 1 ( μ ) N c ,
where N c = 3 is the number of colors and the upper (lower) sign apply when i is odd (even). The helicity amplitude can be written as a linear combination of the a i ( μ ) parameters as follows
A h = G F 2 λ t ( s ) i = 3 10 a i h χ h ( B s 0 ¯ ϕ , ϕ ) .
To calculate the exact formula of A h ( B s 0 ¯ ϕ ϕ ) , we need to write down the corresponding factorizable amplitude. Moreover, the even and odd terms of a i that contribute to Equation (23) must be calculated. From Equation (5), the factorizable amplitude of B s 0 ¯ ( b s ¯ ) ϕ ( s s ¯ ) ϕ ( s s ¯ ) process with V 1 = V 2 = ϕ can be written as:
χ h ( B s 0 ¯ ϕ , ϕ ) = ϕ | ( s ¯ s ) V ± A | 0 ϕ | ( s ¯ b ) V A | B s 0 ¯ .
Since there are some of the quark flavors of the operators in Equation (3) do not match the quark flavor in Equation (24). So, they must be Fierz-transformed to contribute to the decay using the following transformations [28]
( m ¯ n ) V A ( k ¯ l ) V A ( k ¯ n ) V A ( m ¯ l ) V A , ( m ¯ n ) V A ( k ¯ l ) V + A 2 ( k ¯ n ) S P ( m ¯ l ) S + P ,
where ( m ¯ n ) S ± P = m ¯ ( 1 ± γ 5 ) n . Also, the color singlet-singlet term of the operators can be obtained by [29]
( m ¯ i n j ) V A ( k ¯ j l i ) V ± A = 1 N c ( m ¯ n ) V A ( k ¯ l ) V ± A .
In the calculation of the A h ( B s 0 ¯ ϕ ϕ ) , following up H. Y. Cheng et al. [16,18], the contribution of the odd and even terms of the effective operators a i ( μ ) in Equation (23) can be derived as following:
For the odd terms of a i , one can directly use the penguin operators from Equation (3) to get
O 3 = O 5 = χ h ( B s 0 ¯ ϕ , ϕ ) , O 4 = O 6 = 1 N c χ h ( B s 0 ¯ ϕ , ϕ ) , O 7 = O 9 = 1 2 χ h ( B s 0 ¯ ϕ , ϕ ) , O 8 = O 10 = 1 2 N c χ h ( B s 0 ¯ ϕ , ϕ ) ,
where the short-hand notation O i stands for ϕ ϕ O i B s 0 ¯ . Then, substitution from Equation (27) into (21) and using Equation (22) yields
a 3 h + a 5 h 1 2 a 7 h 1 2 a 9 h χ h ( B s 0 ¯ ϕ , ϕ ) .
For the even terms of a i , we can rewrite the operators O 3 , O 4 , O 9 and O 10 in Equation (3) in the form
O 3 = q ( s ¯ i b j ) V A ( q ¯ j q i ) V A , O 4 = q ( q ¯ b ) V A ( s ¯ q ) V A , O 9 = q 3 2 e q ( s ¯ i b j ) V A ( q ¯ j q i ) V A , O 10 = q 3 2 e q ( q ¯ b ) V A ( s ¯ q ) V A .
According to Equation (24), only the quark flavors with q = s contribute to the present decay. A straightforward calculation with using Equations (25) and (26) yields
O 3 = 2 O 9 = 1 N c χ h ( B s 0 ¯ ϕ , ϕ ) , O 4 = 2 O 10 = χ h ( B s 0 ¯ ϕ , ϕ ) ,
Again, substitution from Equation (30) into (21) and using Equation (22) yields
a 4 h 1 2 a 10 h χ h ( B s 0 ¯ ϕ , ϕ ) .
Finally, substituting from Equations (28) and (31) into Equation (23), one can get
A h = G F 2 λ t ( s ) a 3 h + a 4 h + a 5 h 1 2 ( a 7 h + a 9 h + a 10 h ) χ h ( B s 0 ¯ ϕ , ϕ ) .
Now using Equations (1) and (32) we can write the BR formula of the present decay mode with expanding the helicity amplitude as [12]:
B R ( B s 0 ¯ ϕ ϕ ) = τ B s 0 p c s 8 π m B s 0 2 A 0 2 + A + 2 + A 2 .
For the calculation of the B R ( B s 0 ¯ ϕ ϕ ) in the framework of QCDF, we define the effective parameters a i h ( μ ) with i = 3 , 4 , 5 , 7 , 9 , 10 in the next subsection.

2.4. QCD Factorization for B s 0 ¯ ϕ ϕ Process

In the framework of QCDF approach the general form of the effective parameters a i h ( μ ) in the naive dimensional regularization (NDR) scheme at the next-to-leading order (NLO) is given by [12]
a i h ( μ ) = C i ( μ ) + C i ± 1 ( μ ) N c + N i C i ± 1 ( μ ) N c α s C F 4 π f 1 h + f 2 h + P i = 4 h ,
with
N i = 1 for i = 1 , 2 , 3 , 4 , 9 , 10 1 for i = 5 , 7 0 for i = 6 , 8 ,
where the quantities f 1 h account for one loop vertex corrections, f 2 h for hard spectator interactions and P i = 4 h for penguin contribution which has been calculated only for i = 4 . We now write down the explicit expressions for a i h ( μ ) with i = 3 , 4 , 5 , 7 , 9 , 10 which are given by [12]
a 3 h ( μ ) = C 3 + C 4 N c + α s C F 4 π N c C 4 ( f 1 h + f 2 h ) , a 4 h ( μ ) = C 4 + C 3 N c + α s C F 4 π N c C 3 ( f 1 h + f 2 h ) + α s C F 4 π N c C 3 C 9 2 Q h ( β s ) + Q h ( β b ) 4 3 2 3 α s C F 4 π N c C 1 λ u λ t Q h ( β u ) + λ c λ t Q h ( β c ) + 2 3 1 3 + α s C F 4 π N c ( C 4 + C 6 ) + 3 2 C 8 + C 10 Q h ( β s ) + G 8 g Q g h , a 5 h ( μ ) = C 5 + C 6 N c α s C F 4 π N c C 6 ( f 1 h + f 2 h ) , a 7 h ( μ ) = C 7 + C 8 N c α s C F 4 π N c C 8 ( f 1 h + f 2 h ) , a 9 h ( μ ) = C 9 + C 10 N c + α s C F 4 π N c C 10 ( f 1 h + f 2 h ) , a 10 h ( μ ) = C 10 + C 9 N c + α s C F 4 π N c C 9 ( f 1 h + f 2 h ) ,
where C F = 4 3 and β q = m q 2 / m b 2 . In the expression a 4 h ( μ ) , the upper (lower) value in parenthesis corresponds to h = 0 ( ± ) state.
In Equation (36), the contribution from the vertex corrections f 1 h are given by
f 1 0 = 12 log μ m b 18 + 0 1 Φ V 2 ( u ) ϱ ( u ) d u , f 1 ± = 12 log μ m b 16 + 0 1 d u ξ V 2 ± ( u ) ϱ ( u ) + 2 0 1 κ ± ( u , x , y ) d x d y ,
where
ϱ ( u ) = 3 1 2 u 1 u log ( u ) i π ,
ξ V ± ( z ) = g ( υ ) V ( z ) η g ( a ) V ( z ) 4 ,
κ ± ( u , x , y ) = 1 x y x y u x u + y ( 1 x ) u y ( x u + y ) ,
with z { u , v } is the light-cone momentum fraction of the quark in the vector. For V = V 2 in Equation (39), the η = + or − corresponding to ( V A ) ( V A ) or ( V A ) ( V + A ) current, respectively.
For hard spectator interactions f 2 h are given by [12]
f 2 0 = 4 π 2 N c i f B s 0 f V 1 f V 2 χ 0 ( B s 0 V 1 , V 2 ) m B s 0 λ B s 0 0 1 d υ d u Φ V 1 ( υ ) Φ V 2 ( u ) u υ ¯ , f 2 ± = 4 π 2 N c i f B s 0 f V 1 f V 2 m V 2 λ B s 0 χ ± ( B s 0 V 1 , V 2 ) 2 ( 1 ± 1 ) 0 1 d υ d u Φ V 2 ( υ ) ξ V 2 ± ( u ) υ ¯ 2 4 π 2 N c i f B s 0 f V 1 f V 2 m V 1 m V 2 m B s 0 λ B s 0 χ ± ( B s 0 V 1 , V 2 ) 0 1 d υ d u ξ V 1 ± ( υ ) ξ V 2 ± ( u ) u + υ ¯ u υ ¯ 2 ,
where z ¯ = 1 z and the quantity λ B s 0 is the parametrization parameter of the distribution amplitude of the B s 0 meson. Also, the functions Φ V ( z ) , Φ V ( z ) , g ( υ ) V ( z ) and g ( a ) V ( z ) are the light-cone distribution amplitudes (LCDAs) of the vector meson and we adopt them in the following asymptotic form [12]
Φ V ( z ) = Φ V ( z ) = g ( a ) V ( z ) = 6 z z ¯ , g ( υ ) V ( z ) = 3 4 1 + ( 2 z 1 ) 2 .
Also
g ( a ) V ( z ) = d g ( a ) V ( z ) d z .
The non-factorizable corrections induced by local four-quark operators O i can be described by the function Q h ( β q ) which is given by [12,19]
Q 0 ( β q ) = 2 3 + 4 3 log μ m b 4 0 1 d u Φ V 2 ( u ) g ( u , β q ) , Q ± ( β q ) = 2 3 + 2 3 log μ m b 2 0 1 d u ξ V 2 ± ( u ) g ( u , β q ) ,
with the function
g ( u , β q ) = 0 1 d x x x ¯ log [ β q x x ¯ u ¯ i ϵ ] .
In Equation (36), we also take into account the contributions of the dipole operator O 8 g which will give a tree-level contribution described by the function Q g h defined as [12,19]
Q g 0 = 0 1 d u Φ V 2 ( u ) u ¯ , Q g + = 0 1 d u ξ V 2 ( u ) , and Q g = 0 1 d u ξ V 2 ( u ) u ¯ .
Finally, one can calculate the effective parameters a 3 h ( μ ) , a 4 h ( μ ) , a 5 h ( μ ) , a 7 h ( μ ) , a 9 h ( μ ) and a 10 h ( μ ) by substitution the equations from (37) to (46) into Equation (36). Consequently, the the three helicity amplitudes A 0 , A + and A can be calculated from Equation (32) which leads to the calculation of the BR by Equation (33).

2.5. CP Violation

We derive the equations for time dependence of the CP asymmetry ( A C P ( t ) ), direct CP violation ( A C P dir ), CP violation due to mixing ( A C P mix ) and CP violation due to interference ( A C P Δ Γ ) as [4]:
A C P ( t ) = A C P dir cos Δ m s t + A C P mix sin Δ m s t cosh Δ Γ s t / 2 A C P Δ Γ sinh Δ Γ s t / 2 ,
where Δ Γ s = ( Γ L ( s ) Γ H ( s ) ) / 2 = τ B s 0 1 is the difference decay width of the B s system with Γ H s and Γ L s being the decay widths of the “heavy” and “light” mass eigenstates of the B s system, respectively. The Δ m s = m s L m s H is the mass difference of B s system. The three CP observables in Equation (47) are defined by [4]
A C P dir ( B s f ) 1 ζ f ( s ) 2 1 + ζ f ( s ) 2 ,
A C P mix ( B s f ) 2 Im ( ζ f ( s ) ) 1 + ζ f ( s ) 2 ,
A C P Δ Γ ( B s f ) 2 Re ( ζ f ( s ) ) 1 + ζ f ( s ) 2 ,
where
ζ f ( s ) = e i ϕ s A ¯ ( B s f ) A ( B s f ) .
In the SM, for the amplitude of the penguin B d decays, we can write [4]
A ( B d f ) 1 + λ 2 b e i θ e i γ
Using Equation (52) and in analogy to the penguin B d system [4], one can rewrite Equation (51) as follows
ζ f ( s ) = e i ϕ s 1 + λ 2 b e i θ e i γ 1 + λ 2 b e i θ e + i γ ,
where ϕ s is the mixing phase of B s system, b e i θ is a penguin parameter, γ is the angle of the unitarity triangle of the CKM matrix and λ ( = V u s ) is the Wolfenstein parameter with V u s is the CKM matrix element. The values of ϕ s , γ and V u s are given in Table 1. According to the Ref. [4], one can define the parameter b e i θ by
b e i θ R b 1 λ 2 A h u A h c .
Here R b = 1 λ 2 / 2 V u b / V c b / λ , where V u b and V c b are the C K M matrix elements which are given in Table 1. Substitution from Equation (32) into Equation (54) yields
b e i θ R b 1 λ 2 a 3 u , h + a 4 u , h + a 5 u , h 1 2 ( a 7 u , h + a 9 u , h + a 10 u , h ) a 3 c , h + a 4 c , h + a 5 c , h 1 2 ( a 7 c , h + a 9 c , h + a 10 c , h ) .
From Equation (36), only the definition of the a 4 p , h ( μ ) parameter can be rewritten as follows
a 4 p , h ( μ ) = C 4 + C 3 N c + α s C F 4 π N c C 3 ( f 1 h + f 2 h ) + α s C F 4 π N c C 3 C 9 2 Q h ( β s ) + Q h ( β b ) 4 3 2 3 α s C F 4 π N c C 1 λ p λ t Q h ( β p ) + 2 3 1 3 + α s C F 4 π N c ( C 4 + C 6 ) + 3 2 C 8 + C 10 Q h ( β d ) + G 8 g Q g h ,
where p = u , c . Moreover, we use the same definitions for the other parameters in Equation (36) so that
a i p , h ( μ ) = a i h ( μ ) for i = 3 , 5 , 7 , 9 , 10 .
Hence, the penguin parameter b e i θ can easily be calculated by combining Equations (36) and (55)–(57).

3. Numerical Results and Discussions

3.1. Numerical Results of the Branching Ratio for B s 0 ¯ ϕ ϕ

To calculate the effective parameters a i h with i = 3 , 4 , 5 , 7 , 9 , 10 of Equation (36), we use the Mathematica packages. The used input parameters are given in Table 1 and Table 2 where the renormalization scale μ = m b is used in the calculations. Since there are a logarithmic infrared divergence integrals in the f 2 ± expression in Equation (41), so we use the following approximations [5,12]
0 1 d z z = log m b Λ h and 0 1 d z z 2 = m b Λ h ,
with Λ h = 0.5 GeV. The results of a i h with i = 3 , 4 , 5 , 7 , 9 , 10 corresponding to the three helicities h = 0 , ± are listed in Table 3.
Table 4 shows the calculated values of the helicity amplitudes A 0 , A + and A . One can calculate them by computing the corresponding factorizable amplitudes χ 0 ( B s 0 ¯ ϕ , ϕ ) , χ + ( B s 0 ¯ ϕ , ϕ ) and χ ( B s 0 ¯ ϕ , ϕ ) from Equation (18) with using the results of Table 3 into Equation (32).
Finally, we get the B R ( B s 0 ¯ ϕ ϕ ) by using Equation (33) and Table 4 to be
B R ( B s 0 ¯ ϕ ϕ ) = ( 1.56 ± 0.23 ) × 10 5 .
The main sources of uncertainty of the calculated B R come from the uncertainties of both the decay constants and the form factors.
Table 5 shows a comparison between the present calculation of B R ( B s 0 ¯ ϕ ϕ ) in Equation (59) and the available experimental and other theoretical values. This table shows that the present predicted results are very much consistent with the theoretical and experimental ones.

3.2. Numerical Results for the CP Violation for B s 0 ¯ ϕ ϕ

To calculate the observables of the CP violation, we have estimated the effective parameters using the Equation (56). For h = ± , we found that the results of these observables are negligibly small, so we report only the calculations for h = 0 . In this case the estimated effective parameters using Equation (56) are
a 4 u , 0 = 0.0373402 0.000321103 i , a 4 c , 0 = 0.0537056 0.00735794 i ,
and according to Equation (57), the parameters a i p , 0 with i = 3 , 5 , 7 , 9 , 10 are listed in Table 3. Substitution from Equation (60) into Equation (55), one can get the corresponding value of the penguin parameter to be
b e i θ = 0.282421 0.0369209 i .
By imposing the value of Equation (61) into Equation (53) we obtain
ζ ϕ ϕ ( s ) = 0.996437 0.00626556 i .
Using Equation (62) with Equations (48)–(50) one can get
A C P dir ( B s 0 ¯ ϕ ϕ ) = 0.00355 ± 0.00152 ,
A C P mix ( B s 0 ¯ ϕ ϕ ) = 0.00629 ± 0.03119 ,
A C P Δ Γ ( B s 0 ¯ ϕ ϕ ) = 0.99997 ± 0.00019 .
The main sources of uncertainty of the C P asymmetries in Equations (63)–(65) come from the uncertainties of the mixing parameters (i.e., Δ m s , Δ Γ s , ϕ s , and γ ).
Table 6 summarizes the results for A C P dir of the present calculation and for comparison we list the available theoretical values. This table shows that the present calculated value is consistent with the other theoretical ones.
The three C P asymmetries A C P dir , A C P mix and A C P Δ Γ are related observables by the relation [4]
A C P dir 2 + A C P mix 2 + A C P Δ Γ 2 = 1 .
Because the B s 0 ¯ ϕ ϕ mode is induced only by penguin operators, its direct C P asymmetry ( A C P dir ) is naturally zero at leading order (LO) contributions [18,20]. After the inclusion of the NLO contributions, its direct C P asymmetry is nonzero but still very small. Also, in the standard model mixing-induced C P asymmetry of pure penguin decays is predicted to be small [18]. Then, the value of A C P Δ Γ is naturally large as shown in Equation (65).
The asymmetry as a function of time is shown in Figure 1. To show the asymmetry we use our calculated values in Equations (63)–(65). From this figure one can see a clear CP asymmetry with time for the present decay process.
Even though some previous theorists employed QCDF approach [5,18] but the present work is different in several aspects from them. In the present work, the structure of the effective parameters ( a i h ) is different from the one of Refs. [5,18]. We also have reported the CP violation which has not been done before for the present decay channel. Finally we include the error analysis of the BR and CP violation calculations. Therefore we think that the reported results in this paper will be helpful for future experiment as well as theoretical studies.

4. Conclusions

In this work, we studied the B s 0 ¯ ϕ ϕ decay mode in the framework of the QCDF approach. In doing this we employed QCDF approach in the NDR scheme at NLO contributions. In this study, we calculated the branching ratio and the C P violation. After numerical evaluation, we found that the BR of B s 0 ¯ ϕ ϕ decay is ( 1.56 ± 0.23 ) × 10 5 . This value for B R is consistent with the experimental values as well as with the other theoretically predicted ones. We also calculated the CP violation asymmetries and we found that the present decay mode is governed by the longitudinal amplitude (i.e., for h = 0 ). The calculated values for A C P dir , A C P mix and A C P Δ Γ are 0.00355 ± 0.00152 , 0.00629 ± 0.03119 and 0.99997 ± 0.00019 , respectively. These values are consistent with the available theoretical ones.

Author Contributions

Conceptualization, methodology, investigation, formal analysis, and writing of the original draft were carried out by E.H.R. and H.R.K. Also, E.H.R. was responsible for theoretical calculations, data curation, and plotting as well as tabulating them. H.R.K. has revised the original draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research (DSR) of King Saud University for funding and supporting this research through the initiative of DSR Graduate Students Research Support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sakarov, A.D. Violation of CP invariance, C asymmetry, and baryon asymmetry of the universe. JETP Lett. 1967, 5, 24–27. [Google Scholar]
  2. Kobayashi, M.; Maskawa, T. CP-violation in the renormalizable theory of weak interaction. Prog. Theor. Phys. 1973, 49, 652–657. [Google Scholar] [CrossRef] [Green Version]
  3. Du, D.; Yang, D.; Zhu, G. Analysis of the decays Bππ and πk with QCD factorization in the heavy quark limit. Phys. Lett. B 2000, 488, 46–54. [Google Scholar] [CrossRef] [Green Version]
  4. Fleischer, R. Flavour Physics and CP Violation. arXiv 2006, arXiv:hep-ph/0608010. [Google Scholar]
  5. Beneke, M.; Rohrer, J.; Yang, D. Branching fractions, polarisation and asymmetries of BVV decays. Nucl. Phys. B 2007, 774, 64–101. [Google Scholar] [CrossRef] [Green Version]
  6. Abe, K.; Abe, R.; Adachi, I.; Ahn, B.S.; Aihara, H.; Akatsu, M.; Alimonti, G.; Asai, K.; Asai, M.; Asano, Y.; et al. Observation of large CP violation in the neutral B meson system. Phys. Rev. Lett. 2001, 87, 091802. [Google Scholar] [CrossRef] [Green Version]
  7. Aubert, B.; Boutigny, D.; De Bonis, I.; Gaillard, J.M.; Jeremie, A.; Karyotakis, Y.; Lees, J.P.; Robbe, P.; Tisserand, V.; Palano, A.; et al. Measurement of CP-violating asymmetries in B0 decays to CP eigenstates. Phys. Rev. Lett. 2001, 86, 2515. [Google Scholar] [CrossRef] [Green Version]
  8. Acosta, D.; Adelman, J.; Affolder, T.; Akimoto, T.; Albrow, M.G.; Ambrose, D.; Amerio, S.; Amidei, D.; Anastassov, A.; Anikeev, K.; et al. Evidence for B s 0 ϕ ϕ Decay and Measurements of Branching Ratio and ACP for B+ϕK+. Phys. Rev. Lett. 2005, 95, 031801. [Google Scholar] [CrossRef] [Green Version]
  9. Aaltonen, T.; González, B.Á.; Amerio, S.; Amidei, D.; Anastassov, A.; Annovi, A.; Antos, J.; Apollinari, G.; Appel, J.A.; Apresyan, A.; et al. Measurement of Polarization and Search for CP Violation in B s 0 ϕ ϕ Decays. Phys. Rev. Lett. 2011, 107, 261802. [Google Scholar] [CrossRef] [Green Version]
  10. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. (LHCb Collaboration) Measurement of the B s 0 ϕ ϕ branching fraction and search for the decay B 0 ϕ ϕ . JHEP 2015, 10, 053. [Google Scholar] [CrossRef] [Green Version]
  11. Olsen, S.L. Flavor Physics in the LHC Era. arXiv 2014, arXiv:1409.0273. [Google Scholar]
  12. Li, X.; Lu, G.; Yang, Y.D. Charmless B s 0 V V decays in QCD factorization. Phys. Rev. D 2003, 68, 114015, Erratum in 2005, 71, 019902. [Google Scholar] [CrossRef] [Green Version]
  13. Beneke, M.; Buchalla, G.; Neubert, M.; Sachrajda, C.T. QCD factorization in BπK,ππ decays and extraction of Wolfenstein parameters. Nucl. Phys. B 2001, 606, 245–321. [Google Scholar] [CrossRef] [Green Version]
  14. Beneke, M.; Neubert, M. Flavor-singlet B-decay amplitudes in QCD factorization. Nucl. Phys. B 2003, 651, 225–248. [Google Scholar] [CrossRef] [Green Version]
  15. Beneke, M.; Neubert, M. QCD factorization for BPP and BPV decays. Nucl. Phys. B 2003, 675, 333–415. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, Y.H.; Cheng, H.Y.; Tseng, B.; Yang, K.C. Charmless hadronic two-body decays of Bu and Bd mesons. Phys. Rev. D 1999, 60, 094014. [Google Scholar] [CrossRef] [Green Version]
  17. Cheng, H.Y.; Yang, K.C. Branching ratios and polarization in BVV, VA, AA decays. Phys. Rev. D 2008, 78, 094001. [Google Scholar] [CrossRef]
  18. Cheng, H.Y.; Chua, C.K. QCD factorization for charmless hadronic Bs decays revisited. Phys. Rev. D 2009, 80, 114026. [Google Scholar] [CrossRef] [Green Version]
  19. Kamal, A.; Khan, H.R.; Alhendi, H.A. Analytic Calculations of the Branching Ratio and CP Violation in B s 0 K * 0 K ¯ * 0 Decay. Phys. Atom. Nuclei 2019, 82, 299–308, Erratum in 2019, 82, 549. [Google Scholar] [CrossRef]
  20. Yan, D.C.; Liu, X.; Xiao, Z.J. Anatomy of BsVV decays and effects of next-to-leading order contributions in the perturbative QCD factorization approach. Nucl. Phys. B 2008, 935, 17–39. [Google Scholar] [CrossRef]
  21. Zou, Z.T.; Ali, A.; Lü, C.D.; Liu, X.; Li, Y. Improved estimates of the BsVV decays in perturbative QCD approach. Phys. Rev. D 2015, 91, 054033. [Google Scholar] [CrossRef] [Green Version]
  22. Bauer, C.W.; Pirjol, D.; Rothstein, I.Z.; Stewart, I.W. BM1M2: Factorization, charming penguins, strong phases, and polarization. Phys. Rev. D 2004, 70, 054015. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, C.; Zhou, S.H.; Li, Y.; Lü, C.D. Global analysis of charmless B decays into two vector mesons in soft-collinear effective theory. Phys. Rev. D 2017, 96, 073004. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, C.; Zhang, Q.A.; Li, Y.; Lü, C.D. Charmless B(s)VV decays in factorization-assisted topological-amplitude approach. Eur. Phys. J. C 2017, 77, 333. [Google Scholar] [CrossRef] [Green Version]
  25. Patrignani, C.P.D.G.; Weinberg, D.H.; Woody, C.L.; Chivukula, R.S.; Buchmueller, O.; Kuyanov, Y.V.; Blucher, E.; Willocq, S.; Höcker, A.; Lippmann, C.; et al. Review of particle physics. Chin. Phys. C 2016, 40, 100001. [Google Scholar]
  26. Buchalla, G.; Buras, A.J.; Lautenbacher, M.E. Weak decays beyond leading logarithms. Rev. Mod. Phys. 1996, 68, 1125. [Google Scholar] [CrossRef] [Green Version]
  27. Thomson, M. Modern Particle Physics, 1st ed.; Cambridge University Press: Cambridge, UK, 2013; pp. 528–530. ISBN 978-1-107-03426-6. [Google Scholar]
  28. Datta, A.; London, D. Triple-product correlations in BVV decays and new physics. Int. J. Mod. Phys. A 2004, 19, 2505–2544. [Google Scholar] [CrossRef] [Green Version]
  29. Ali, A.; Greub, C. Analysis of two-body nonleptonic B decays involving light mesons in the standard model. Phys. Rev. D 1998, 57, 2996. [Google Scholar] [CrossRef] [Green Version]
  30. Tanabashi, M.; Hagiwara, K.; Hikasa, K.; Nakamura, K.; Sumino, Y.; Takahashi, F.; Tanaka, J.; Agashe, K.; Aielli, G.; Amsler, C.; et al. Review of particle physics. Phys. Rev. D 2018, 98, 030001. [Google Scholar] [CrossRef] [Green Version]
  31. Amhis, Y.; Banerjee, S.; Ben-Haim, E.; Bona, M.; Bozek, A.; Bozzi, C.; Brodzicka, J.; Chrzaszcz, M.; Dingfelder, J.; Egede, U.; et al. (HFLAV) Averages of b-hadron, c-hadron and τ-lepton properties as of 2018. arXiv 2019, arXiv:1909.12524. [Google Scholar]
Figure 1. Time dependent CP asymmetry A C P (t).
Figure 1. Time dependent CP asymmetry A C P (t).
Atoms 08 00022 g001
Table 1. Inputs Parameters.
Table 1. Inputs Parameters.
ParameterThe ValueRef.
Mass of B s 0 meson ( m B s 0 ) 5366.88 ± 0.17 MeV[30]
Mean life time of B s 0 ( τ B s 0 ) 1.510 ± 0.004 ps[30]
Mass of ϕ meson ( m ϕ ) 1019.461 ± 0.016 MeV[30]
Fermi coupling constant ( G F ) 1.1663787 × 10 5 GeV 2 [30]
The fine-structure coupling ( α ) 1 / 128 [30]
The strong coupling constant ( α s ) 0.01181 ± 0.0011 [30]
The CKM matrix element V u b 0.00394 ± 0.00036 [30]
The CKM matrix element V u s 0.2243 ± 0.0005 [30]
The CKM matrix element V c b 0.0422 ± 0.0008 [30]
The CKM matrix element V c s 0.997 ± 0.017 [30]
The CKM matrix element V t b 1.019 ± 0.025 [30]
The CKM matrix element V t s 0.0394 ± 0.0023 [30]
The mass difference of B s 0 system ( Δ m s ) 17.757 ± 0.021 ps 1 [30]
The decay width difference of B s 0 system ( Δ Γ s ) 0.090 ± 0.005 ps 1 [30]
Decay constant of B s 0 meson ( f B s 0 ) 240 ± 30 MeV[21]
Longitudinal decay constant of ϕ meson ( f ϕ ) 215 ± 5 MeV[21]
Transverse decay constant of ϕ meson ( f ϕ ) 186 ± 9 MeV[21]
The form factor A 1 B s ϕ 260 ± 10 [18]
The form factor A 2 B s ϕ 230 ± 10 [18]
The form factor V B s ϕ 300 ± 10 [18]
The parametrization parameter λ B s 0 350 ± 150 MeV[15]
Mass of b quark ( m b )4660 MeV[12]
Mass of c quark ( m c )1470 MeV[12]
Mass of u quark ( m u )0[12]
Mass of d quark ( m d )0[12]
Mass of s quark ( m s )0[12]
The phase mixing ϕ s 0.021 ± 0.031 rad[31]
The C K M angle γ ( 71 . 1 5.3 + 4.6 ) o [31]
Table 2. Wilson coefficients in the NDR scheme at NLO with μ = m b [12].
Table 2. Wilson coefficients in the NDR scheme at NLO with μ = m b [12].
Wilson CoefficientThe Value
C 1 1.078
C 2 0.176
C 3 0.014
C 4 0.034
C 5 0.008
C 6 0.039
C 7 / α 0.011
C 8 / α 0.055
C 9 / α 1.341
C 10 / α 0.264
C 8 g 0.146
Table 3. The effective parameters in QCDF approach at NLO.
Table 3. The effective parameters in QCDF approach at NLO.
a i h h = 0 h = + h =
a 3 h 0.00529251 + 0.00133847 i 0.354452 + 0.00133847 i 0.00587474 + 0.00133847 i
a 4 h 0.0536064 0.00756553 i 0.0957453 0.00233928 i 0.0406903 0.0188806 i
a 5 h 0.008012 0.0015353 i 11.765 0.0015353 i 0.00654938 0.0015353 i
a 7 h 0.0000904768 + 0.0000169154 i 0.12962 + 0.0000169154 i 0.0000699551 + 0.0000169154 i
a 9 h 0.00994835 0.0000811938 i 0.0118744 0.0000811938 i 0.00927093 0.0000811937 i
a 10 h 0.000620574 + 0.000412427 i 0.11147 + 0.000412427 i 0.00406159 + 0.000412427 i
Table 4. Helicity amplitudes in QCDF approach at NLO.
Table 4. Helicity amplitudes in QCDF approach at NLO.
A h The Value
A 0 0.00521217 + 0.0335507 i
A + 0.000572798 2.45607 i
A 0.000571254 + 0.000988547 i
Table 5. Theoretical and experimental B R ( B s 0 ¯ ϕ ϕ ) in units of 10 5 .
Table 5. Theoretical and experimental B R ( B s 0 ¯ ϕ ϕ ) in units of 10 5 .
Theory/ExperimentThe Value
Present Work 1.56 ± 0.23
QCDF [5] 2.18 1.71 + 3.04
QCDF [18] 1.67 0.91 + 1.19
PQCD [20] 1.88 0.38 + 0.49
PQCD [21] 1.67 0.38 + 0.49
SCET [23] 1.90 ± 0.65
FAT [24] 2.64 ± 0.76
CDF Experiment [8] 1.4 0.5 + 0.6 ( stat . ) ± 0.6 ( syst . )
CDF Experiment [9] 2.32 ± 0.18 ( stat . ) ± 0.82 ( syst . )
LHCb Experiment [10] 1.84 ± 0.05 ( stat . ) ± 0.07 ( syst . )
Particle Data Group [25] 1.87 ± 0.15
Table 6. The predicted values of the direct CP asymmetry A C P dir (%).
Table 6. The predicted values of the direct CP asymmetry A C P dir (%).
Approach A C P dir (%)
Present Work 0.35 ± 0.15
QCDF [18] 0.2 0.3 + 0.6
PQCD [20] 0.7 ± 0.2
SCET [23] 0.39 ± 0.44
FAT [24] 0.83 ± 0.28

Share and Cite

MDPI and ACS Style

Raslan, E.H.; Khan, H.R. Branching Ratio and CP Violation in B s 0 ¯ ϕϕ Decay in the Framework of QCD Factorization. Atoms 2020, 8, 22. https://doi.org/10.3390/atoms8020022

AMA Style

Raslan EH, Khan HR. Branching Ratio and CP Violation in B s 0 ¯ ϕϕ Decay in the Framework of QCD Factorization. Atoms. 2020; 8(2):22. https://doi.org/10.3390/atoms8020022

Chicago/Turabian Style

Raslan, E. H., and H. R. Khan. 2020. "Branching Ratio and CP Violation in B s 0 ¯ ϕϕ Decay in the Framework of QCD Factorization" Atoms 8, no. 2: 22. https://doi.org/10.3390/atoms8020022

APA Style

Raslan, E. H., & Khan, H. R. (2020). Branching Ratio and CP Violation in B s 0 ¯ ϕϕ Decay in the Framework of QCD Factorization. Atoms, 8(2), 22. https://doi.org/10.3390/atoms8020022

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop