We investigate the decay mechanism in the B − → φφK − decay with the φφ invariant mass below the charm threshold and in the neighborhood of the η c invariant mass region. Our approach is based on the use of factorization model and the knowledge of matrix elements of the weak currents. For the B meson weak transition we apply form factor formalism, while for the light mesons we use effective weak and strong Lagrangians. We find that the dominant contributions to the branching ratio come from the η, η ′ and η(1490)) pole terms of the penguin operators in the decay chains B − → η(η ′ , η(1490))K − → φφK − . Our prediction for the branching ratio is in agreement with the Belle's result.
I. INTRODUCTIONIt is a very fruitful era in B meson physics. A lot of experimental data on B meson decays is coming from the B meson factories. Many of their results are still not explained. Recently, Belle collaboration has announced the observation of the BR(B ± → φφK ± ) = (2.6for a φφ invariant mass below 2.85 GeV. This is the first of the three-body B decays with two vector mesons and one pseudoscalar meson in the final state that has been observed. The B meson decays into three pseudoscalar mesons have been studied [2,3] within heavy quark symmetry accompanied by chiral symmetry. One might explain the observed rates using heavy quark symmetry for the strong vertices, while for the weak transition we rely on the existing knowledge of the form factors [2]. The three-body decay with two vector meson states and one pseudoscalar is much more difficult to approach.The additional insight on the decay mechanism might come from the analysis of the B meson two-body decays. Particularly interesting are the decaysThey have been extensively studied using different existing techniques: the naive factorization [4,5,6], the QCD factorization [7] and the SU(3) symmetry [8]. Each of these decay modes is rather difficult to explain theoretically. The decays B ± → φK ± and B ± → K * ± φ might have significant annihilation contribution [4, 6, 9], but it is not simple to have a consistent treatment of it. There is an interesting proposal [6] in which the angular distributions of the final outgoing particles can be used to estimate the magnitude of annihilation contribution to the amplitude. However, we have to wait for the new experimental data to extract the size of the annihilation contribution. The B ± → η(η ′ )K ± decay rate has not been easy to explain. It accounts for the wellknown problem of the η − η ′ mixing [10,11] as can be seen from a variety of approaches used for 4,12,13,14]. In the B ± → η(η ′ )K ± decay mode, it seems that the annihilation contribution is not very significant [4,13].One has to expect that the above described difficulties in these decay modes might appear in the three-body decay we discuss. Based on the current knowledge of twobody transitions, we build a simple model which might clarify the role of the non-charm contributions in the BR(B ± → φφK ± ) decay. In our study of the B ± → φφK ± decay mechanism, we f...