Search

Lattice calculation may play an important role to understand the cause of the long-standing puzzle among the determination of $|V_{cb}|$. The key element is the computation of the inclusive decay rate, for which the formalism is under development. The method has its own problem: the approximation of the phase-space factor and finite volume effect. I discuss the current status of such study and future prospects., Comment: 10 pages, 4 figures; a talk presented at EuroPLEx2023, Sep 2023, Berlin

Lattice calculation of $B\to D^{(*)}\ell\bar\nu$ decay form factors is now available from multiple lattice collaborations for both zero-recoil limit and non-zero recoil kinematics. Yet, there are a number of challenges we face to better control the systematic errors and to extend the application of lattice calculations to related quantities. The subjects are chosen from my own perspective and not meant to be conprehensive., Comment: 9 pages, 4 figures; a talk presented at EuroPLEx2023, Sep 2023, Berlin

With the development of lattice QCD in recent years, a determination of decay rates for inclusive semileptonic decays from lattice correlators now seems viable. We report on the calculation of the inclusive semileptonic decay of the $D_s$ meson on the lattice. We simulate the $D_s \rightarrow X_s\ell\nu_\ell$ process with M\"obius domain-wall charm and strange quarks, whose masses were approximately tuned to the physical values. We cover the whole kinematical region. We present a general overview on how the inclusive decays are treated on the lattice and discuss ongoing challenges, such as the estimation of systematic errors and future prospects of the project., Comment: 7 pages, 2 figures, Presented at the 12th Workshop on the CKM Unitarity Triangle, 18-22 September 2023, Santiago de Compostela

We study the $U(1)_A$ anomaly at high temperatures of $N_f=2+1$ lattice QCD with chiral fermions. Gauge ensembles are generated with M\"obius domain-wall (MDW) fermions, and the measurements are reweighted to those with overlap fermions. We report on the results for the Dirac spectra, the $U(1)_A$ susceptibility, and the topological susceptibility in the temperature range of $T=136$, $153$, $175$, and $204$ MeV, where the up and down quark masses are set to be near the physical point as well as at lighter or heavier masses., Comment: 10 pages, 4 figures, 1 table; Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermilab, Batavia, Illinois, USA

Based on simulations of 2+1 flavor lattice QCD with M\"obius domain wall fermions at high temperatures, we compute a series of spatial correlation functions to study the screening masses in mesonic states. We compare these masses with the symmetry relations for various quark masses and lattice sizes at temperatures above the critical point. Using these spatial correlation functions we examine the $SU(2)_L \times SU(2)_R$ symmetry as well as the anomalously broken axial $U(1)_A$ symmetry. Additionally we explore a possible and emergent chiral-spin symmetry $SU(2)_{CS}$., Comment: 10 pages, 5 figures, Contribution to the 40th International Symposium on Lattice Field Theory (Lattice 2023)

We present an update on the study of the QCD phase transition with 3 flavors of M\"obius domain wall fermions at zero baryon density. We performed simulations on lattices of size $36^3\times12\times16$ and $24^3\times12\times32$ with a variety of quark masses at a fixed lattice spacing $a=0.1361(20)$ fm, which correspond to a temperature 121(2) MeV. By analyzing the chiral condensate, chiral susceptibilitities and Binder cumulant on $36^3\times12\times16$ lattices together with the result obtained from our previous study on $24^3\times12\times16$ lattices, we identified a crossover occurring at quark mass around $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 3-4$ MeV for this temperature. Besides, we show the effects of residual chiral symmetry breaking on chiral condensate and chiral susceptibilities between $L_s=16$ and 32., Comment: 10 pages, 6 figures, Contribution to the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31 to August 4, 2023, Fermilab, USA

We report on our on-going study of the $B_s \to K\ell\nu$ decay in $N_f=2+1$ lattice QCD. We employ fully relativistic setup in which the M\"obius domain-wall action is used for all quark flavors. The lattice cutoff is $a^{-1} \sim 2.5$ GeV, where we take bottom quark masses up to $m_Q\!<\!0.7a^{-1}$ in order to control discretization errors. We present preliminary results for the relevant form factors extracted from correlator ratios by inspecting their ground state saturation., Comment: 8 pages, $ figures, Lattice 2023

We present a study on the nonperturbative calculation of observables for inclusive semileptonic decays of $B_{(s)}$ mesons using lattice QCD. We focus on the comparison of two different methods to analyse the lattice data of Euclidean correlation functions, specifically Chebyshev and Backus-Gilbert approaches. This type of computation may eventually provide new insight into the long-standing tension between the inclusive and exclusive determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{cb}|$ and $|V_{ub}|$. We report the results from a pilot lattice computation for the decay $B_s \rightarrow X_c \, l\nu_l$, where the valence quark masses are approximately tuned to their physical values using the relativistic-heavy quark action for the $b$ quark and the domain-wall formalism for the other valence quarks. We address the computation of the total decay rate as well as leptonic and hadronic moments, discussing similarities and differences between the two analysis techniques., Comment: 10 pages, 6 figures, 40th International Symposium on Lattice Field Theory - Lattice2023

We report on the calculation of the inclusive semileptonic decay of the $D_s$ meson on the lattice. We simulate the $D_s \rightarrow X_s\ell\nu_\ell$ process with M\"obius domain-wall charm and strange quarks, whose masses were approximately tuned to the physical values. We cover the whole kinematical region. The focus of this work is on the systematic error due to finite-volume effects. We construct a model of two-body final states to describe the data on a finite volume lattice of $L \simeq 0.055 \, \text{fm}$ to investigate the extrapolation to the infinite-volume limit.

In this contribution we describe a recent study focused on the lattice calculation of inclusive decay rates of heavy mesons. We show how the inclusive calculation can be achieved starting from four-point lattice correlation functions normalised appropriately. The correlators used in this project come from gauge ensembles provided by the JLQCD and ETM collaborations. An essential point of this method is the extraction of spectral densities from lattice correlators which is obtained using two of the most recent approaches in the literature. Our results are in remarkable agreement with analytical predictions from the operator-product expansion. This study represents the first step towards a full lattice QCD study of heavy mesons inclusive semileptonic decays., Comment: 6 pages, 2 figures, contribution to IFAE 2023, Catania (Italy), 12-14 April 2023

We address the nonperturbative calculation of the inclusive decay rate of semileptonic $B_{(s)}$-meson decays from lattice QCD. Precise Standard-Model predictions are key ingredients in searches for new physics, and this type of computation may eventually provide new insight into the long-standing tension between the inclusive and exclusive determinations of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements $|V_{cb}|$ and $|V_{ub}|$. We present results from a pilot lattice computation for $B_s \rightarrow X_c\, l \nu_l$, where the initial $b$ quark described by the relativistic-heavy-quark (RHQ) formalism on the lattice and the other valence quarks discretised with domain-wall fermions are simulated approximately at their physical quark masses. We compare two different methods for computing the decay rate from lattice data of Euclidean $n$-point functions, namely Chebyshev and Backus-Gilbert approaches. We further study how much the ground-state meson dominates the inclusive decay rate and indicate our strategy towards a computation with a more comprehensive systematic error budget., Comment: 43 pages, 18 figures. Add missing references, fix typos

We report on our finite temperature 2+1 flavor lattice QCD simulation to study the thermodynamic properties of QCD near the (pseudo) critical point employing $N_T=12$ and $16$. The simulation points are chosen along the lines of constant physics. The quark mass for M\"obius domain-wall fermion are tuned by taking into account the residual mass either by fine-tuning the input quark masses or by post-process using reweighting. In this talk, we focus on simulation details and present some preliminary results., Comment: 8 pages, contribution to the 39th International Symposium on Lattice Field theory (LATTICE2022), 8th-13th August, 2022, Bonn, Germany

We investigate the finite temperature QCD phase transition with three degenerate quark flavors using M\"{o}bius domain wall fermions. To explore the order of phase transition on the lower left corner of Columbia plot and if possible, to locate the critical endpoint we performed simulations at temperatures around 181 and 121 MeV with lattice spacing $a=0.1361(20)$ fm corresponding to temporal lattice extent $N_{\tau}=8,12$ with varying quark mass for two different volumes with aspect ratios $N_{\sigma}/N_{\tau}$ ranging from 2 to 3. By analyzing the volume and mass dependence of the chiral condensate, disconnected chiral susceptibility and Binder cumulant we find that there is a crossover at $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 44\, \mathrm{MeV}$ for $\mathrm{T_{pc}}\sim$ 181 MeV, At temperature 121 MeV, the binder cumulant suggests a crossover at $m_q^{\mathrm{\overline {MS}}}(2\, \mathrm{GeV}) \sim 3.7\, \mathrm{MeV}$, although a study of volume dependence would be important to confirm this., Comment: 9 pages, 5 figures, Contribution to the 39th International Symposium on Lattice Field Theory (Lattice 2022), 08-13 August 2022, Bonn, Germany

We perform a non-perturbative lattice calculation of the decay rates for inclusive semi-leptonic decays of charmed mesons. In view of the long-standing tension in the determination of the CKM matrix elements $|V_{ub}|$ and $|V_{cb}|$ from exclusive and inclusive processes, recently, the use of lattice QCD has been extended towards the description of inclusive decays. Since the determination of hadronic input parameters from QCD based methods require independent tests, we focus on the charm sector, since it not only offers experimental data, but also well determined CKM parameters. We carry out a pilot lattice simulation for the $D_s \rightarrow X_s \ell\nu$ and explore the improvement of existing techniques. Our simulation employs M\"obius domain-wall charm and strange quarks whose masses are tuned to be approximately physical and we cover the whole kinematical region. We report on our progress in analyzing different sources of systematic effects, such as the extrapolation of the kernel function chosen for the Chebsyhev approximation as well as the influence on the analysis in the region close to the kinematical limit., Comment: 9 pages, 6 figures, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

We address the non-perturbative calculation of the decay rate of inclusive semi-leptonic $B_{(s)}$ mesons decays from lattice QCD. Precise theoretical Standard Model predictions are key ingredients in searches for new physics. This type of computation may eventually provide new insight into the long-standing tension between the inclusive and exclusive determinations of the CKM matrix elements $|V_{cb}|$ and $|V_{ub}|$. We perform a pilot lattice computation for $B_s \rightarrow X_c \, l \nu_l$ and improve on existing techniques. The valence-quark masses in our simulations are approximately physical for the domain-wall strange and charm quarks as well as for the bottom quark, for which we use a relativistic heavy quark effective action. We report on our progress and discuss future plans towards a first study with fully controlled systematic effects., Comment: 10 pages, 6 figures, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

We present the lattice QCD calculation of inclusive semileptonic $B_s$-meson decays. We follow a recently proposed method, which is based on the extraction of smeared spectral densities from Euclidean correlation functions and on the numerical reconstruction of the integration kernel relevant for the inclusive decay rate calculation. We compute four-point Euclidean correlation functions using JLQCD and ETM gauge ensembles with unphysically light $b$-quark masses, and apply two different methods for the integration kernel reconstruction. Finally, we show that the lattice results obtained in this work are in good agreement with the analytic predictions of the operator-product-expansion. This opens the path for a future full lattice QCD calculation to be used as theoretical input for the determination of the magnitude of the CKM element $V_{cb}$., Comment: 10 pages, 4 figures, contribution to the 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Bonn, Germany. arXiv admin note: substantial text overlap with arXiv:2209.15494

We present a non-perturbative computation of inclusive rates of semileptonic decays of heavy mesons from lattice QCD simulations. The calculation is based on the extraction of smeared spectral functions obtained from four-point Euclidean correlation functions computed on configuration ensembles of the JLQCD and ETM collaborations. We compare our results for the inclusive decay rates with analytical predictions from the operator-product expansion, finding a good agreement for the calculation of the inclusive decay rate. This opens the path to the theoretical determination of the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element $V_{cb}$ to a level of precision competitive with the present experimental uncertainty., Comment: 7 pages, 1 figure, to appear in the proceedings of the 51st International Symposium on Multiparticle Dynamics (ISMD)

We study the $U(1)_A$ anomaly in the high-temperature phase of $N_f=2+1$ lattice QCD with chiral fermions. Gauge ensembles are generated with M\"obius domain-wall (MDW) fermions, and in the measurements the determinant is reweighted to that of overlap fermions. We report the results for the overlap Dirac spectrum, $U(1)_A$ susceptibility, and topological susceptibility at $T=204$ and $175$ MeV., Comment: 7 pages, 4 figures, 1 table; Proceedings of 38th International Symposium on Lattice Field Theory (LATTICE2021), 26-30 July 2021, Zoom/Gather@MIT, USA

We present an ab initio study of inclusive semileptonic decays of heavy mesons from lattice QCD. Our approach is based on a recently proposed method, that allows one to address the study of these decays from the analysis of smeared spectral functions extracted from four-point correlators on the lattice, where the smearing is defined in terms of the phase-space integration relevant to the inclusive decays. We present results obtained from gauge-field ensembles from the JLQCD and ETM collaborations, and discuss their relation with theoretical predictions from the operator-product expansion., Comment: 51 pages, 27 figures and 4 tables. Clarifications added, v2 to appear in JHEP

Using a fully relativistic lattice fermion action, we compute the form factors of the semileptonic decay $B\to\pi\ell\nu$, which is required for the determination of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{ub}|$. We employ the M\"{o}bius domain-wall fermion formalism for the generation of lattice ensembles with 2+1 sea quark flavours as well as for the valence heavy and light quarks. We compute the form factors at various values of the lattice spacing and multiple light and heavy quark masses, and extrapolate the results to the physical point. We combine our lattice results with the available experimental data to obtain $|V_{ub}| = (3.93\pm 0.41)\times 10^{-3}$., Comment: 34 pages, 20 figures, 12 tables. Plus supplemental material: 20 pages, 7 figures, 8 tables. Updated to reflect minor changes in published version, including additional references

Simulations for the thermodynamics of the 2+1 flavor QCD are performed employing chiral fermions. The use of M\"obius domain-wall fermions with stout-link smearing is more effective on the finer lattices where all the relevant chiral symmetries are realized more accurately. We report on the initial simulations near the (pseudo) critical point using the line of constant physics with an average $ud$ quark mass slightly heavier than physical at $a\lesssim 0.1$ fm., Comment: 9 pages, 8 figures, 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@MIT

Utilizing the approach recently proposed for the inclusive semi-leptonic decay rate on the lattice, we compute the differential decay rate of a $B_s$ meson for various kinematical channels. The results are compared with the contributions from the ground states (D and D^*) as well as from the orbitally excited states (D^**'s). The computation so far is carried out with an unphysically light bottom quark and strange spectator quark., Comment: 8 pages, 3 figures; contribution to The 38th International Symposium on Lattice Field Theory, LATTICE2021, 26th-30th July, 2021

We propose a new method to renormalize lattice operators. The method is based on the technique to compute the spectral sum appearing in the Shifman-Vainshtein-Zakharov QCD sum rule from lattice correlators. The application of this technique to the light quark system is useful for operator renormalization as well as for the test of perturbative QCD and OPE. We determine the renormalization constant of the vector current and discuss extensions to other current operators., Comment: 8 pages, 2 figures, 1 table, talk presented at The 38th International Symposium on Lattice Field Theory (LATTICE2021), 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Utilizing the approach recently proposed for the $\ell N$ inclusive scattering cross section on the lattice, we compute the differential scattering cross section for the charged current process $\ell p \rightarrow \nu n $ for various kinematical channels. The simulation is carried out on the 2+1 flavor $16^3 \times 32 $ ensemble with Iwasaki and Domain Wall Fermion action. The lattice results are compared with MINER$\nu$A result for the equivalent process., Comment: 9pages, LATTICE 2021

A recently proposed approach to compute the inclusive semi-leptonic decay rate on the lattice allows for the calculation of various quantities (differential distributions and moments) in different subchannels. We systematically compare the lattice QCD results for an unphysically light bottom quark with OPE-based predictions, including leading order perturbative and $\mathcal{O}\left(1/m^3 \right)$ power corrections, and explore possible strategies to decrease the systematic uncertainty., Comment: 8 pages, 4 figures, proceedings to the 38th International Symposium on Lattice Field Theory (LATTICE21), 26-30 July 2021

We perform a numerical computation of the anomalous magnetic moment ($g-2$) of the electron in QED by using the stochastic perturbation theory. Formulating QED on the lattice, we develop a method to calculate the coefficients of the perturbative series of $g-2$ without the use of the Feynman diagrams. We demonstrate the feasibility of the method by performing a computation up to the $\alpha^3$ order and compare with the known results. This program provides us with a totally independent check of the results obtained by the Feynman diagrams and will be useful for the estimations of not-yet-calculated higher order values. This work provides an example of the application of the numerical stochastic perturbation theory to physical quantities, for which the external states have to be taken on-shell., Comment: 35 pages, 5 figures, 2 tables, the version to appear in JHEP

We propose a method to use lattice QCD to compute the Borel transform of the vacuum polarization function appearing in the Shifman-Vainshtein-Zakharov QCD sum rule. We construct the spectral sum corresponding to the Borel transform from two-point functions computed on the Euclidean lattice. As a proof of principle, we compute the $s \bar{s}$ correlators at three lattice spacings and take the continuum limit. We confirm that the method yields results that are consistent with the operator product expansion in the large Borel mass region. The method provides a ground on which the OPE analyses can be directly compared with nonperturbative lattice computations., Comment: 27 pages, 13 figures, 1 table; data updated, 1 figure added; version to appear in PRD

We propose a fully non-perturbative method to compute inelastic lepton-nucleon ($\ell N$) scattering cross sections using lattice QCD. The method is applicable even at low energies, such as the energy region relevant for the recent and future neutrino-nucleon ($\nu N$) scattering experiments, for which perturbative analysis is invalidated. The basic building block is the forward Compton-scattering amplitude, or the hadronic tensor, computed on a Euclidean lattice. Total cross section is constructed from the hadronic tensor by multiplying a phase space factor and integrating over the energy and momentum of final hadronic states. The energy integral that induces a sum over all possible final states is performed implicitly by promoting the phase space factor to an operator written in terms of the transfer matrix on the lattice. The formalism is imported from that of the inclusive semileptonic B meson decay [P. Gambino, S. Hashimoto, arXiv:2005.13730]. It can be generalized to compute the $\ell N$ scattering cross sections and their moments, as well as the virtual correction to the nuclear $\beta$-decay. Necessary quark-line contractions for two current insertions corresponding to the Compton amplitude to be computed on the lattice are summarized., Comment: 44 pages, 16 figures; minor corrections, a few refs added

We develop a method to compute inclusive semi-leptonic decay rate of hadrons fully non-perturbatively using lattice QCD simulations. The sum over all possible final states is achieved by a calculation of the forward-scattering matrix elements on the lattice, and the phase-space integral is evaluated using their dependence on the time separation between two inserted currents. We perform a pilot lattice computation for the B_s -> X_c l nu decay with an unphysical bottom quark mass and compare the results with the corresponding OPE calculation. The method to treat the inclusive processes on the lattice can be applied to other processes, such as the lepton-nucleon inelastic scattering., Comment: 5 pages, 4 figures

We propose a method to reconstruct smeared spectral functions from two-point correlation functions measured on the Euclidean lattice. Arbitrary smearing function can be considered as far as it is smooth enough to allow an approximation using Chebyshev polynomials. We test the method with numerical lattice data of Charmonium correlators. The method provides a framework to compare lattice calculation with experimental data including excited state contributions without assuming quark-hadron duality., Comment: 23 pages

We report the current status of a study of charmonium contribution to $B \rightarrow K\ell^+\ell^-$ on the lattice. Our lattice calculation tests the factorization approximation for this contribution. In order to control the problem of the artificial divergence, we focus on the low $q^2$ region with a small b-quark mass. We also take into account the renormalization constants of relevant four-quark operators calculated through the temporal moments. Results suggest a violation of the factorization approximation., Comment: 7 pages, 5 figures

We propose a renormalization scheme that can be simply implemented on the lattice. It consists of the temporal moments of two-point and three-point functions calculated with finite valence quark mass. The scheme is confirmed to yield a consistent result with another renormalization scheme in the continuum limit for the bilinear operators. We apply a similar renormalization scheme for the non-perturbative renormalization of four-fermion operators appearing in the weak effective Hamiltonian., Comment: 7 pages, 3 figures, 2 tables; Proceedings of the 37th annual International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

We investigate the high-temperature phase of QCD using lattice QCD simulations with $N_f = 2$ dynamical M\"obius domain-wall fermions. On generated configurations, we study the axial $U(1)$ symmetry, overlap-Dirac spectra, screening masses from mesonic correlators, and topological susceptibility. We find that some of the observables are quite sensitive to lattice artifacts due to a small violation of the chiral symmetry. For those observables, we reweight the M\"obius domain-wall fermion determinant by that of the overlap fermion. We also check the volume dependence of observables. Our data near the chiral limit indicates a strong suppression of the axial $U(1)$ anomaly at temperatures $\geq$ 220 MeV., Comment: 7 pages, 4 figures, 1 table; Proceedings of the 37th International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

We report on a calculation of form factors for the semileptonic decay of $B$ meson to pion on $2+1$-flavour lattices with lattice spacings from 0.080 fm down to 0.044 fm. Using the M\"obius domain wall fermion action for both sea and valence quarks, we simulate pions with masses down to 225 MeV. By utilizing a range of heavy quark masses up to 2.44 times the mass of the charm quark we extrapolate to the physical $b$ quark mass. We discuss the dependence of the form factors on the pion mass, heavy quark mass, lattice spacing and the momentum-transfer. We extract the CKM matrix element $|V_{ub}|$ through a simultaneous fit with the $B \to \pi\ell\nu$ differential branching fractions provided by the Belle and BaBar collaborations after a chiral-continuum and physical $b$ quark extrapolations of our lattice data., Comment: 7 pages, Presented at the 37th International Symposium on Lattice Field Theory - Lattice2019, 16-22 June 2019, Wuhan, China

Using lattice QCD simulations with $N_f = 2$ dynamical fermions, we study the axial $U(1)$ symmetry, topological charge, and Dirac eigenvalue spectra in the high-temperature phase in which the chiral symmetry is restored. Our gauge ensembles are generated with M\"obius domain-wall fermions, but the measurements such as susceptibilities are reweighted to those for the overlap fermions by using overlap/domain-wall reweighting technique. We find that the $U(1)_A$ and topological susceptibilities are strongly suppressed in the small quark mass region, which is related to the reduction of chiral-zero and low-nonzero modes on the Dirac spectra. We also examine their volume dependence., Comment: 10 pages, 4 figures, 1 table; Proceedings of The 9th International Workshop on Chiral Dynamics (CD18), 17-21 September, 2018, Durham, NC, USA. arXiv admin note: text overlap with arXiv:1812.06621

Heavy flavor physics entered a new era when the Belle II experiment observed its first collision. There are several hints found so far by BaBar, Belle, and LHCb in particular, that suggest the physics beyond the Standard Model appearing in the loop processes at short distances. They will be further tested by higher precision experiments in the coming years, while the role of lattice QCD is to understand the long-distance physics quantitatively so that one can unambiguously isolate the short-distance physics from the experimental data. I summarize the status towards this goal and then look at the challenges we are facing., Comment: 20 pages, 4 figures, talk presented at the 36th Annual International Symposium on Lattice Field Theory - LATTICE2018, 22-28 July, 2018, Michigan State University, East Lansing, Michigan, USA

We report on a calculation of the charmonium contribution to the decay $B \rightarrow K\ell^+\ell^-$ using lattice simulations with 2+1 flavors of Mobius domain wall fermions. We focus on the region of $q^2$ below the $J/\psi$ resonance and test the factorization approximation to estimate the amplitude. We show a possible discrepancy between the lattice calculation and the factorization method for the non-factorizable contribution., Comment: 7pages, 5 figs. Proceedings of the 36th International Symposium on Lattice Field Theory (Lattice 2018), July 22-28, 2018, East Lancing, USA