Showing 1–50 of 81 results for author: Ishikawa, T

$K^-d\rightarrowπぱいΛらむだN$ reaction with in-flight kaons for studying the $ΛらむだN$ interaction

Authors: Shunsuke Yasunaga, Daisuke Jido, Takatsugu Ishikawa

Abstract: The $ΛらむだN$ invariant mass spectra for the reactions $K^-d\rightarrowπぱい^-Λらむだp$ and $K^-d\rightarrowπぱい^0Λらむだn$ are calculated for experimental study of isospin symmetry breaking in the $ΛらむだN$ scattering at low energies, the difference in the scattering lengths and effective ranges of $Λらむだp$ and $Λらむだn$ systems. The calculations are performed for in-flight kaons with a momentum of 1000 MeV/c with employing partial w… ▽ More The $ΛらむだN$ invariant mass spectra for the reactions $K^-d\rightarrowπぱい^-Λらむだp$ and $K^-d\rightarrowπぱい^0Λらむだn$ are calculated for experimental study of isospin symmetry breaking in the $ΛらむだN$ scattering at low energies, the difference in the scattering lengths and effective ranges of $Λらむだp$ and $Λらむだn$ systems. The calculations are performed for in-flight kaons with a momentum of 1000 MeV/c with employing partial wave analysis up to the p-wave for meson-baryon amplitudes and the spin-flip term for baryon-baryon amplitudes. Kinematic selection is utilized to suppress the background processes by selecting forward-emitting pions and higher momentum nucleons. It is worth noting that isospin symmetry breaking in the $ΛらむだN$ system can be extracted from the difference of the $ΛらむだN$ invariant mass spectra between the $K^-d\rightarrowπぱい^-Λらむだp$ and $K^-d\rightarrowπぱい^0Λらむだn$ reactions. △ Less

Submitted 28 June, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Analytic and numerical approaches for depictive 3-loop integrals using sector decomposition

Authors: Elise de Doncker, Tadashi Ishikawa, Kiyoshi Kato, Fukuko Yuasa

Abstract: Four 3-loop two-point functions are studied analytically and numerically using a simplified sector decomposition method. The coefficients of the ultraviolet divergent part are determined analytically, and those of the finite part are computed numerically. The energy dependence of the integrals is shown explicitly, and a discussion of its behavior is presented. Four 3-loop two-point functions are studied analytically and numerically using a simplified sector decomposition method. The coefficients of the ultraviolet divergent part are determined analytically, and those of the finite part are computed numerically. The energy dependence of the integrals is shown explicitly, and a discussion of its behavior is presented. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 14 pages

$K^- d\rightarrow πぱいΛらむだN$ reaction for studying charge symmetry breaking in the $ΛらむだN$ interaction

Authors: Yutaro Iizawa, Daisuke Jido, Takatsugu Ishikawa

Abstract: We discuss charge symmetry breaking in the $ΛらむだN$ interaction. In order to investigate the $Λらむだp$ and $Λらむだn$ interactions at low energies, we propose to utilize the $K^-d\to πぱい^-Λらむだp$ and $K^-d\to πぱい^0Λらむだn$ reactions. They are symmetric under the exchange of both the pion and nucleon isospin partners in the final states. This advantage allows us to study charge symmetry breaking in the $ΛらむだN$ interaction. We ca… ▽ More We discuss charge symmetry breaking in the $ΛらむだN$ interaction. In order to investigate the $Λらむだp$ and $Λらむだn$ interactions at low energies, we propose to utilize the $K^-d\to πぱい^-Λらむだp$ and $K^-d\to πぱい^0Λらむだn$ reactions. They are symmetric under the exchange of both the pion and nucleon isospin partners in the final states. This advantage allows us to study charge symmetry breaking in the $ΛらむだN$ interaction. We calculate the differential cross sections of these reactions with stopped kaons theoretically and discuss the possibility to extract the scattering length and effective range of the $ΛらむだN$ scattering. With stopped kaons, the $ΛらむだN$ interaction takes place dominantly in the spin-triplet state thanks to the deuteron spin and $s$-wave dominance of the scattering amplitudes at the low energy. We find that the ratio of the differential cross sections as a function of the $ΛらむだN$ invariant mass between the two reactions is useful for revealing how charge symmetry breaking, or isospin symmetry breaking appears in the low-energy $ΛらむだN$ scattering. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 13 pages, 19 figures, 1table

Journal ref: Phys.Rev.C 106, 045201 (2022)

Spectral sum from Euclidean lattice correlators and determination of renormalization constants

Authors: Tsutomu Ishikawa, Shoji Hashimoto, Takashi Kaneko

Abstract: 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… ▽ More 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. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 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

Report number: KEK-CP-0388

Kondo effect with Wilson fermions

Authors: Tsutomu Ishikawa, Katsumasa Nakayama, Kei Suzuki

Abstract: We investigate the Kondo effect with Wilson fermions. This is based on a mean-field approach for the chiral Gross-Neveu model including four-point interactions between a light Wilson fermion and a heavy fermion. For massless Wilson fermions, we demonstrate the appearance of the Kondo effect. We point out that there is a coexistence phase with both the light-fermion scalar condensate and Kondo cond… ▽ More We investigate the Kondo effect with Wilson fermions. This is based on a mean-field approach for the chiral Gross-Neveu model including four-point interactions between a light Wilson fermion and a heavy fermion. For massless Wilson fermions, we demonstrate the appearance of the Kondo effect. We point out that there is a coexistence phase with both the light-fermion scalar condensate and Kondo condensate, and the critical chemical potentials of the scalar condensate are shifted by the Kondo effect. For negative-mass Wilson fermions, we find that the Kondo effect is favored near the parameter region realizing the Aoki phase. Our findings will be useful for understanding the roles of heavy impurities in Dirac semimetals, topological insulators, and lattice simulations. △ Less

Submitted 29 December, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 11 pages, 8 figures; published version

Report number: DESY 21-110, KEK-TH-2338

Journal ref: Phys. Rev. D 104, 094515 (2021)

Spectral sum of current correlators from lattice QCD

Authors: Tsutomu Ishikawa, Shoji Hashimoto

Abstract: 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 conti… ▽ More 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. △ Less

Submitted 12 November, 2021; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: 27 pages, 13 figures, 1 table; data updated, 1 figure added; version to appear in PRD

Report number: KEK-CP-0381

Journal ref: Phys.Rev.D 104, 074521 (2021)

Neutron-neutron scattering length from $πぱい^+$ photoproduction on the deuteron

Authors: Satoshi X. Nakamura, Takatsugu Ishikawa, Toru Sato

Abstract: We discuss the possibility of extracting the neutron-neutron scattering length $a_{nn}$ and effective range $r_{nn}$ from cross section data ($d^2σしぐま/dM_{nn}/dΩおめが_πぱい$), as a function of the $nn$ invariant mass $M_{nn}$, for $πぱい^+$ photoproduction on the deuteron ($γがんまd\to πぱい^+nn$). The analysis is based on a $γがんまd\to πぱい^+nn$ reaction model in which realistic elementary amplitudes for $γがんまp\to πぱい^+n$, $NN\to NN$,… ▽ More We discuss the possibility of extracting the neutron-neutron scattering length $a_{nn}$ and effective range $r_{nn}$ from cross section data ($d^2σしぐま/dM_{nn}/dΩおめが_πぱい$), as a function of the $nn$ invariant mass $M_{nn}$, for $πぱい^+$ photoproduction on the deuteron ($γがんまd\to πぱい^+nn$). The analysis is based on a $γがんまd\to πぱい^+nn$ reaction model in which realistic elementary amplitudes for $γがんまp\to πぱい^+n$, $NN\to NN$, and $πぱいN\to πぱいN$ are built in. We show that $M_{nn}$ dependence (lineshape) of a ratio $R_{\rm th}$, $d^2σしぐま/dM_{nn}/dΩおめが_πぱい$ normalized by $dσしぐま/dΩおめが_πぱい$ for $γがんまp\toπぱい^+ n$ and the nucleon momentum distribution inside the deuteron, at the kinematics with $θしーた_πぱい=0^\circ$ and $E_γがんま\sim 250$ MeV is particularly useful for extracting $a_{nn}$ and $r_{nn}$ from the corresponding data $R_{\rm exp}$. It is found that $R_{\rm exp}$ with 2% error, resolved into the $M_{nn}$ bin width of 0.04 MeV (corresponding to the $p_πぱい$ bin width of 0.05 MeV$/c$), can determine $a_{nn}$ and $r_{nn}$ with uncertainties of $\pm 0.21$ fm and $\pm 0.06$ fm, respectively, for the case of $a_{nn}=-18.9$ fm and $r_{nn}=2.75$ fm. The requirement of such narrow bin widths indicates that the momenta of the incident photon and the emitted $πぱい^+$ have to be measured with high resolutions. This can be achieved by utilizing virtual photons of very small $Q^2$ from electron scattering at Mainz MAMI facility. The proposed method for determining $a_{nn}$ and $r_{nn}$ from $γがんまd\to πぱい^+ nn$ has a great experimental advantage over the previous one utilizing $πぱい^- d\toγがんまnn$ for being free from the formidable task of controlling the neutron detection efficiency and its uncertainty. △ Less

Submitted 16 March, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 22 pages, 15 figures, several refinements (not major change)

Reconstruction of smeared spectral function from Euclidean correlation functions

Authors: Gabriela Bailas, Shoji Hashimoto, Tsutomu Ishikawa

Abstract: 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 calculatio… ▽ More 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. △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: 23 pages

Report number: KEK-CP-374

Journal ref: Prog Theor Exp Phys (2020)

Charmonium contribution to $B \rightarrow K\ell^+\ell^-: testing the factorization approximation on the lattice

Authors: Katsumasa Nakayama, Tsutomu Ishikawa, Shoji Hashimoto

Abstract: 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-quar… ▽ More 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. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 7 pages, 5 figures

$D$ mesons as a probe of Casimir effect for chiral symmetry breaking

Authors: Tsutomu Ishikawa, Katsumasa Nakayama, Daiki Suenaga, Kei Suzuki

Abstract: We propose $D$ mesons as probes to investigate finite-volume effects for chiral symmetry breaking at zero and finite temperature. By using the $2+1$-flavor linear-sigma model with constituent light quarks, we analyze the Casimir effects for the $σしぐま$ mean fields: The chiral symmetry is rapidly restored by the antiperiodic boundary for light quarks, and the chiral symmetry breaking is catalyzed by th… ▽ More We propose $D$ mesons as probes to investigate finite-volume effects for chiral symmetry breaking at zero and finite temperature. By using the $2+1$-flavor linear-sigma model with constituent light quarks, we analyze the Casimir effects for the $σしぐま$ mean fields: The chiral symmetry is rapidly restored by the antiperiodic boundary for light quarks, and the chiral symmetry breaking is catalyzed by the periodic boundary. We also show the phase diagram of the $σしぐま$ mean fields on the volume and temperature plane. For $D$ mesons, we employ an effective model based on the chiral-partner structure, where the volume dependence of $D$ mesons is induced by the $σしぐま$ mean fields. We find that $D_s$ mesons are less sensitive to finite volume than $D$ mesons, which is caused by the insensitivity of $σしぐま_s$ mean fields. An anomalous mass shift of $D$ mesons at high temperature with the periodic boundary will be useful in examinations with lattice QCD simulations. The dependence on the number of compactified spatial dimensions is also studied. △ Less

Submitted 26 August, 2019; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 15 pages, 13 figures, 3 tables; published version

Report number: KEK-TH-2134, OU-HET-1019

Journal ref: Phys. Rev. D 100, 034016 (2019)

Casimir effect for nucleon parity doublets

Authors: Tsutomu Ishikawa, Katsumasa Nakayama, Kei Suzuki

Abstract: Finite-volume effects for the nucleon chiral partners are studied within the framework of the parity-doublet model. Our model includes the vacuum energy shift for nucleons, which is the Casimir effect. We find that for the antiperiodic boundary the finite-volume effect leads to chiral symmetry restoration, and the masses of the nucleon parity doublets degenerate. For the periodic boundary, the chi… ▽ More Finite-volume effects for the nucleon chiral partners are studied within the framework of the parity-doublet model. Our model includes the vacuum energy shift for nucleons, which is the Casimir effect. We find that for the antiperiodic boundary the finite-volume effect leads to chiral symmetry restoration, and the masses of the nucleon parity doublets degenerate. For the periodic boundary, the chiral symmetry breaking is enhanced, and the masses of the nucleons also increase. We also discuss the finite-temperature effect and the dependence on the number of compactified spatial dimensions. △ Less

Submitted 18 March, 2019; v1 submitted 28 December, 2018; originally announced December 2018.

Comments: 14 pages, 11 figures, 3 tables; published version

Report number: KEK-TH-2097

Journal ref: Phys. Rev. D 99, 054010 (2019)

$ηいーた$ photoproduction off the deuteron and low-energy $ηいーた$-nucleon interaction

Authors: S. X. Nakamura, H. Kamano, T. Ishikawa

Abstract: We study $ηいーた$ photoproduction off the deuteron ($γがんまd\toηいーたpn$) at a special kinematics: $\sim 0.94$ GeV of the photon beam energy and $\sim 0^\circ$ of the scattering angle of the proton. This kinematics is ideal to extract the low-energy $ηいーた$-nucleon scattering parameters such as $a_{ηいーたN}$ (scattering length) and $r_{ηいーたN}$ (effective range) because the $ηいーた$-nucleon elastic scattering is significantly enh… ▽ More We study $ηいーた$ photoproduction off the deuteron ($γがんまd\toηいーたpn$) at a special kinematics: $\sim 0.94$ GeV of the photon beam energy and $\sim 0^\circ$ of the scattering angle of the proton. This kinematics is ideal to extract the low-energy $ηいーた$-nucleon scattering parameters such as $a_{ηいーたN}$ (scattering length) and $r_{ηいーたN}$ (effective range) because the $ηいーた$-nucleon elastic scattering is significantly enhanced. We show that if a ratio $R$, the $γがんまd\toηいーたpn$ cross section divided by the $γがんまp\toηいーたp$ cross section convoluted with the proton momentum distribution in the deuteron, is measured with 5% error, ${\rm Re}[a_{ηいーたN}]$ (${\rm Re}[r_{ηいーたN}]$) can be determined at the precision of $\sim\pm$0.1 fm ($\sim\pm$0.5 fm), significantly narrowing down the currently estimated range of the parameters. The measurement is ongoing at the Research Center for Electron Photon Science (ELPH), Tohoku University. △ Less

Submitted 1 December, 2018; originally announced December 2018.

Comments: 4 pages, 5 figures, Contribution to the Proceedings for 8th International Conference on Quarks and Nuclear Physics (QNP2018), November 13-17, 2018, Tsukuba, Japan

Report number: LFTC-18-16/37

Numerical calculation of the full two-loop electroweak corrections to muon (g-2)

Authors: Tadashi Ishikawa, Nobuya Nakazawa, Yoshiaki Yasui

Abstract: Numerical calculation of two-loop electroweak corrections to the muon anomalous magnetic moment ($g$-2) is done based on, on shell renormalization scheme (OS) and free quark model (FQM). The GRACE-FORM system is used to generate Feynman diagrams and corresponding amplitudes. Total 1780 two-loop diagrams and 70 one-loop diagrams composed of counter terms are calculated to get the renormalized quant… ▽ More Numerical calculation of two-loop electroweak corrections to the muon anomalous magnetic moment ($g$-2) is done based on, on shell renormalization scheme (OS) and free quark model (FQM). The GRACE-FORM system is used to generate Feynman diagrams and corresponding amplitudes. Total 1780 two-loop diagrams and 70 one-loop diagrams composed of counter terms are calculated to get the renormalized quantity. As for the numerical calculation, we adopt trapezoidal rule with Double Exponential method (DE). Linear extrapolation method (LE) is introduced to regularize UV- and IR-divergences and to get finite values. The reliability of our result is guaranteed by several conditions. The sum of one and two loop electroweak corrections in this renormalization scheme becomes $a_μみゅー^{EW:OS}[1{\rm+}2{\rm -loop}]= 151.2 (\pm 1.0)\times 10^{-11}$, where the error is due to the numerical integration and the uncertainty of input mass parameters and of the hadronic corrections to electroweak loops. By taking the hadronic corrections into account, we get $a_μみゅー^{EW}[1{\rm+}2 {\rm -loop}]= 152.9 (\pm 1.0)\times 10^{-11}$. It is in agreement with the previous works given in PDG within errors. △ Less

Submitted 24 March, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 12 pages, 6 figures. arXiv admin note: substantial text overlap with arXiv:1709.03284

Journal ref: Phys. Rev. D 99, 073004 (2019)

One-loop effects of MSSM particles in e^-e^+ \to Zh and e^-e^+ \to νにゅー\barνにゅーh at the ILC

Authors: Yusaku Kouda, Tadashi Kon, Yoshimasa Kurihara, Tadashi Ishikawa, Masato Jimbo, Kiyoshi Kato, Masaaki Kuroda

Abstract: The 1-loop effects of the MSSM at the ILC are investigated through numerical analysis. We studied the higgs production processes $e^-e^+\rightarrow Zh$ and $e^-e^+\rightarrow νにゅー\barνにゅーh$ at the ILC. It is found that the magnitude of the MSSM contribution through the 1-loop effects is sizable enough to be detected. In the study, three sets of the MSSM parameters are proposed, which are consistent with… ▽ More The 1-loop effects of the MSSM at the ILC are investigated through numerical analysis. We studied the higgs production processes $e^-e^+\rightarrow Zh$ and $e^-e^+\rightarrow νにゅー\barνにゅーh$ at the ILC. It is found that the magnitude of the MSSM contribution through the 1-loop effects is sizable enough to be detected. In the study, three sets of the MSSM parameters are proposed, which are consistent with the observed higgs mass, the muon $g$-$2$, the dark matter abundance and the decay branching ratios of $B$ mesons. In the $e^-e^+\rightarrow Zh$ process, the 1-loop effects of the MSSM are visible and the distinction of the parameter sets is partially possible. For the study of $e^-e^+\rightarrow νにゅー\barνにゅーh$, we used the equivalent $\it W$-boson approximation in the evaluation of the 1-loop cross section. While the 1-loop effect of the MSSM is visible, the distinction of the parameter sets might not be possible in this process under the value of realistic luminosity at the ILC. △ Less

Submitted 6 July, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Journal ref: Prog Theor Exp Phys (2018)

A development of an accelerator board dedicated for multi-precision arithmetic operations and its application to Feynman loop integrals II

Authors: H Daisaka, N Nakasato, T Ishikawa, F Yuasa, K Nitadori

Abstract: Evaluation of a wide variety of Feynman diagrams with multi-loop integrals and physical parameters and its comparison with high energy experiments are expected to investigate new physics beyond the Standard Model. We have been developing a direct computation method of multi-loop integrals of Feynman diagrams. One of features of our method is that we adopt the double exponential rule for numerical… ▽ More Evaluation of a wide variety of Feynman diagrams with multi-loop integrals and physical parameters and its comparison with high energy experiments are expected to investigate new physics beyond the Standard Model. We have been developing a direct computation method of multi-loop integrals of Feynman diagrams. One of features of our method is that we adopt the double exponential rule for numerical integrations which enables us to evaluate loop integrals with boundary singularities. Another feature is that in order to accelerate the numerical integrations with multi-precision calculations, we develop an accelerator system with Field Programmable Gate Array boards on which processing elements with dedicated logic for quadruple/hexuple/octuple precision arithmetic operations are implemented. In addition, we also develop a programming interface designed for easy use of the system. The development is continued for practical use of the system. We present the current development of our system, and the numerical results of higher-loop diagrams performed using our system. △ Less

Submitted 19 March, 2018; originally announced March 2018.

Comments: 7 pages, 4 figures, in proceedings of the 18th International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT 2017)

Direct numerical computation and its application to the higher-order radiative corrections

Authors: K Kato, E de Doncker, T Ishikawa, F Yuasa

Abstract: The direct computation method(DCM) is developed to calculate the multi-loop amplitude for general masses and external momenta. The ultraviolet divergence is under control in dimensional regularization. In this paper we report on the progress of DCM to several scalar multi-loop integrals after the presentation in ACAT2016. Also the discussion is given on the application of DCM to physical 2-loop pr… ▽ More The direct computation method(DCM) is developed to calculate the multi-loop amplitude for general masses and external momenta. The ultraviolet divergence is under control in dimensional regularization. In this paper we report on the progress of DCM to several scalar multi-loop integrals after the presentation in ACAT2016. Also the discussion is given on the application of DCM to physical 2-loop processes including numerator functions. △ Less

Submitted 14 March, 2018; originally announced March 2018.

Comments: 5 pages, in proceedings of the 18th International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT 2017)

Perturbative matching of continuum and lattice quasi-distributions

Authors: Tomomi Ishikawa

Abstract: Matching of the quasi parton distribution functions between continuum and lattice is addressed using lattice perturbation theory specifically with Wilson-type fermions. The matching is done for nonlocal quark bilinear operators with a straight Wilson line in a spatial direction. We also investigate operator mixing in the renormalization and possible O(a) operators for the nonlocal operators based… ▽ More Matching of the quasi parton distribution functions between continuum and lattice is addressed using lattice perturbation theory specifically with Wilson-type fermions. The matching is done for nonlocal quark bilinear operators with a straight Wilson line in a spatial direction. We also investigate operator mixing in the renormalization and possible O(a) operators for the nonlocal operators based on a symmetry argument on lattice. △ Less

Submitted 7 February, 2018; originally announced February 2018.

Comments: 8 pages, 2 figures, talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

Possible $ηいーた' d$ bound state and its $s$-channel formation in the $γがんまd \to ηいーたd$ reaction

Authors: Takayasu Sekihara, Hiroyuki Fujioka, Takatsugu Ishikawa

Abstract: We theoretically investigate a possibility of an $ηいーた^{\prime} d$ bound state and its formation in the $γがんまd \to ηいーたd$ reaction. First, in the fixed center approximation to the Faddeev equations we obtain an $ηいーた^{\prime} d$ bound state with a binding energy of 25 MeV and width of 19 MeV, where we take the $ηいーた^{\prime} N$ interaction with a coupling to the $ηいーたN$ channel from the linear $σしぐま$ model. Then, in o… ▽ More We theoretically investigate a possibility of an $ηいーた^{\prime} d$ bound state and its formation in the $γがんまd \to ηいーたd$ reaction. First, in the fixed center approximation to the Faddeev equations we obtain an $ηいーた^{\prime} d$ bound state with a binding energy of 25 MeV and width of 19 MeV, where we take the $ηいーた^{\prime} N$ interaction with a coupling to the $ηいーたN$ channel from the linear $σしぐま$ model. Then, in order to investigate the feasibility from an experimental point of view, we calculate the cross section of the $γがんまd \to ηいーたd$ reaction at the photon energy in the laboratory frame around 1.2 GeV. As a result, we find a clear peak structure with the strength $\sim$ 0.2 nb/sr, corresponding to a signal of the $ηいーた^{\prime} d$ bound state in case of backward $ηいーた$ emission. This structure will be prominent because a background contribution coming from single-step $ηいーた$ emission off a bound nucleon is highly suppressed. In addition, the signal can be seen even in case of forward $ηいーた$ emission as a bump or dip, depending on the relative phase between the bound-state formation and the single-step background. △ Less

Submitted 21 April, 2018; v1 submitted 18 December, 2017; originally announced December 2017.

Comments: 9 pages, 13 eps figures, version accepted for publication in PRC, description of the discussion part is improved

Journal ref: Phys. Rev. C 97, 045202 (2018)

Parton distributions and lattice QCD calculations: a community white paper

Authors: Huey-Wen Lin, Emanuele R. Nocera, Fred Olness, Kostas Orginos, Juan Rojo, Alberto Accardi, Constantia Alexandrou, Alessandro Bacchetta, Giuseppe Bozzi, Jiunn-Wei Chen, Sara Collins, Amanda Cooper-Sarkar, Martha Constantinou, Luigi Del Debbio, Michael Engelhardt, Jeremy Green, Rajan Gupta, Lucian A. Harland-Lang, Tomomi Ishikawa, Aleksander Kusina, Keh-Fei Liu, Simonetta Liuti, Christopher Monahan, Pavel Nadolsky, Jian-Wei Qiu , et al. (7 additional authors not shown)

Abstract: In the framework of quantum chromodynamics (QCD), parton distribution functions (PDFs) quantify how the momentum and spin of a hadron are divided among its quark and gluon constituents. Two main approaches exist to determine PDFs. The first approach, based on QCD factorization theorems, realizes a QCD analysis of a suitable set of hard-scattering measurements, often using a variety of hadronic obs… ▽ More In the framework of quantum chromodynamics (QCD), parton distribution functions (PDFs) quantify how the momentum and spin of a hadron are divided among its quark and gluon constituents. Two main approaches exist to determine PDFs. The first approach, based on QCD factorization theorems, realizes a QCD analysis of a suitable set of hard-scattering measurements, often using a variety of hadronic observables. The second approach, based on first-principle operator definitions of PDFs, uses lattice QCD to compute directly some PDF-related quantities, such as their moments. Motivated by recent progress in both approaches, in this document we present an overview of lattice-QCD and global-analysis techniques used to determine unpolarized and polarized proton PDFs and their moments. We provide benchmark numbers to validate present and future lattice-QCD calculations and we illustrate how they could be used to reduce the PDF uncertainties in current unpolarized and polarized global analyses. This document represents a first step towards establishing a common language between the two communities, to foster dialogue and to further improve our knowledge of PDFs. △ Less

Submitted 8 February, 2018; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 80 pages, 15 figures, 28 tables, minor typos corrected, published in Progress in Particle and Nuclear Physics

Report number: DESY 17-185, IFJPAN-IV-2017-19, INT-PUB-17-042, MSUHEP-17-017, Nikhef-2017-047, OUTP-17-15P, SMU-HEP-17-08

Gaussian-weighted Parton Quasi-distribution

Authors: Tomomi Ishikawa, Luchang Jin, Huey-Wen Lin, Andreas Schäfer, Yi-Bo Yang, Jian-Hui Zhang, Yong Zhao

Abstract: We propose a revised definition of quasi-distributions within the framework of large-momentum effective theory (LaMET) that improves convergence towards the large-momentum limit. Since the definition of quasi-distributions is not unique, each choice goes along with a specific matching function, we can use this freedom to optimize convergence towards the large-momentum limit. As an illustration, we… ▽ More We propose a revised definition of quasi-distributions within the framework of large-momentum effective theory (LaMET) that improves convergence towards the large-momentum limit. Since the definition of quasi-distributions is not unique, each choice goes along with a specific matching function, we can use this freedom to optimize convergence towards the large-momentum limit. As an illustration, we study quasi-distributions with a Gaussian weighting factor that naturally suppresses long-range correlations, which are plagued by artifacts. This choice has the advantage that the matching functions can be trivially obtained from the known ones. We apply the Gaussian weighting to the previously published results for the nonperturbatively renormalized unpolarized quark distribution, and find that the unphysical oscillatory behavior is significantly reduced. △ Less

Submitted 10 May, 2019; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 9 pages, 3 figures, version published in SCPMA

Journal ref: Sci.China Phys.Mech.Astron. 62 (2019) no.9, 991021

Symmetry Properties of Nonlocal Quark Bilinear Operators on a Lattice

Authors: Jiunn-Wei Chen, Tomomi Ishikawa, Luchang Jin, Huey-Wen Lin, Jian-Hui Zhang, Yong Zhao

Abstract: Using symmetry properties, we determine the mixing pattern of a class of nonlocal quark bilinear operators containing a straight Wilson line along a spatial direction. We confirm the previous study that mixing among the lowest dimensional operators, which have mass dimension equals three, can occur if chiral symmetry is broken in the lattice action. For higher dimensional operators, we find that t… ▽ More Using symmetry properties, we determine the mixing pattern of a class of nonlocal quark bilinear operators containing a straight Wilson line along a spatial direction. We confirm the previous study that mixing among the lowest dimensional operators, which have mass dimension equals three, can occur if chiral symmetry is broken in the lattice action. For higher dimensional operators, we find that the dimension three operators will always mix with dimension four operators even if chiral symmetry is preserved. Also, the number of dimension four operators involved in the mixing is large hence it is impractical to remove the mixing by the improvement procedure. Our result is important to determining the Bjorken-$x$ dependence parton distribution functions using the quasi-distribution method on a Euclidean lattice. The requirement of using large hadron momentum in this approach makes the control of errors from dimension four operators even more important. △ Less

Submitted 19 July, 2019; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 10 pages, 1 figure; version to be published in Chinese Physics C

Report number: MSUHEP-17-015, MIT-CTP/4942

Automatic calculation of two-loop ELWK corrections to the muon ($g$-2)

Authors: Tadashi Ishikawa, Nobuya Nakazawa, Yoshiaki Yasui

Abstract: Two-loop electroweak corrections to the muon anomalous magnetic moment are automatically calculated by using GRACE-FORM system, as a trial to extend our system for two-loop calculation. We adopt the non-linear gauge (NLG) to check the reliability of our calculation. In total 1780 two-loop diagrams consisting of 14 different topological types and 70 one-loop diagrams composed of counter terms are c… ▽ More Two-loop electroweak corrections to the muon anomalous magnetic moment are automatically calculated by using GRACE-FORM system, as a trial to extend our system for two-loop calculation. We adopt the non-linear gauge (NLG) to check the reliability of our calculation. In total 1780 two-loop diagrams consisting of 14 different topological types and 70 one-loop diagrams composed of counter terms are calculated. We check UV- and IR-divergences cancellation and the independence of the results from NLG parameters. As for the numerical calculation, we adopt trapezoidal rule with Double Exponential method (DE). Linear extrapolation method (LE) is introduced to regularize UV- and IR- divergence and to get finite values. △ Less

Submitted 23 September, 2017; v1 submitted 11 September, 2017; originally announced September 2017.

Improved Parton Distribution Functions at Physical Pion Mass

Authors: Huey-Wen Lin, Jiunn-Wei Chen, Tomomi Ishikawa, Jian-Hui Zhang

Abstract: We present the first lattice results on isovector unpolarized and longitudinally polarized parton distribution functions (PDFs) at physical pion mass. The PDFs are obtained using the large-momentum effective field theory (LaMET) framework where the full Bjorken-$x$ dependence of finite-momentum PDFs, called "quasi-PDFs", can be calculated on the lattice. The quasi-PDF nucleon matrix elements are r… ▽ More We present the first lattice results on isovector unpolarized and longitudinally polarized parton distribution functions (PDFs) at physical pion mass. The PDFs are obtained using the large-momentum effective field theory (LaMET) framework where the full Bjorken-$x$ dependence of finite-momentum PDFs, called "quasi-PDFs", can be calculated on the lattice. The quasi-PDF nucleon matrix elements are renormalized nonperturbatively in RI/MOM-scheme. However, the recent renormalized quasi-PDFs suffer from unphysical oscillations that alter the shape of the true distribution as a function of Bjorken-$x$. In this paper, we propose two possible solutions to overcome this problem, and demonstrate the efficacy of the methods on the 2+1+1-flavor lattice data at physical pion mass with lattice spacing 0.09~fm and volume $(5.76\mbox{ fm})^3$. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Report number: MSUHEP-17-014

Journal ref: Phys. Rev. D 98, 054504 (2018)

On the Renormalizability of Quasi Parton Distribution Functions

Authors: Tomomi Ishikawa, Yan-Qing Ma, Jian-Wei Qiu, Shinsuke Yoshida

Abstract: Quasi-parton distribution functions have received a lot of attentions in both perturbative QCD and lattice QCD communities in recent years because they not only carry good information on the parton distribution functions, but also could be evaluated by lattice QCD simulations. However, unlike the parton distribution functions, the quasi-parton distribution functions have perturbative ultraviolet p… ▽ More Quasi-parton distribution functions have received a lot of attentions in both perturbative QCD and lattice QCD communities in recent years because they not only carry good information on the parton distribution functions, but also could be evaluated by lattice QCD simulations. However, unlike the parton distribution functions, the quasi-parton distribution functions have perturbative ultraviolet power divergences because they are not defined by twist-2 operators. In this paper, we identify all sources of ultraviolet divergences for the quasi-parton distribution functions in coordinate-space, and demonstrate that power divergences, as well as all logarithmic divergences can be renormalized multiplicatively to all orders in QCD perturbation theory. △ Less

Submitted 2 November, 2017; v1 submitted 10 July, 2017; originally announced July 2017.

Comments: 11 pages, 9 figures, version accepted by PRD

Journal ref: Phys. Rev. D 96, 094019 (2017)

Parton Distribution Function with Non-perturbative Renormalization from Lattice QCD

Authors: Jiunn-Wei Chen, Tomomi Ishikawa, Luchang Jin, Huey-Wen Lin, Yi-Bo Yang, Jian-Hui Zhang, Yong Zhao

Abstract: We present lattice results for the isovector unpolarized parton distribution with nonperturbative RI/MOM-scheme renormalization on the lattice. In the framework of large-momentum effective field theory (LaMET), the full Bjorken-$x$ dependence of a momentum-dependent quasi-distribution is calculated on the lattice and matched to the ordinary lightcone parton distribution at one-loop order, with pow… ▽ More We present lattice results for the isovector unpolarized parton distribution with nonperturbative RI/MOM-scheme renormalization on the lattice. In the framework of large-momentum effective field theory (LaMET), the full Bjorken-$x$ dependence of a momentum-dependent quasi-distribution is calculated on the lattice and matched to the ordinary lightcone parton distribution at one-loop order, with power corrections included. The important step of RI/MOM renormalization that connects the lattice and continuum matrix elements is detailed in this paper. A few consequences of the results are also addressed here. △ Less

Submitted 23 January, 2018; v1 submitted 5 June, 2017; originally announced June 2017.

Comments: 11 pages, 6 figures, version accepted by Phys. Rev. D

Report number: MSUHEP-17-007, MIT-CTP/4912

Journal ref: Phys. Rev. D 97, 014505 (2018)

Low-energy $ηいーた$-nucleon interaction studied with $ηいーた$ photoproduction off the deuteron

Authors: S. X. Nakamura, H. Kamano, T. Ishikawa

Abstract: We develop a reaction model for $ηいーた$ photoproduction off the deuteron ($γがんまd\toηいーたpn$), and study the reaction at a special kinematics, where the photon beam energy is $\sim 0.94$ GeV and the scattered proton is detected at $\sim 0^\circ$, for the purpose of determining the $ηいーた$-nucleon scattering length ($a_{ηいーたN}$) and effective range ($r_{ηいーたN}$). In this kinematics, the $ηいーた$-nucleon elastic rescattering… ▽ More We develop a reaction model for $ηいーた$ photoproduction off the deuteron ($γがんまd\toηいーたpn$), and study the reaction at a special kinematics, where the photon beam energy is $\sim 0.94$ GeV and the scattered proton is detected at $\sim 0^\circ$, for the purpose of determining the $ηいーた$-nucleon scattering length ($a_{ηいーたN}$) and effective range ($r_{ηいーたN}$). In this kinematics, the $ηいーた$-nucleon elastic rescattering is significantly enhanced while other background mechanisms being suppressed. We show that a ratio $R$, the $γがんまd\toηいーたpn$ cross section divided by the $γがんまp\toηいーたp$ cross section convoluted with the proton momentum distribution in the deuteron, has a very good resolving power of $a_{ηいーたN}$ and $r_{ηいーたN}$. We conclude that the $R$ data with 5% error, binned in 1 MeV width of the $ηいーた$-neutron invariant mass, can determine ${\rm Re}[a_{ηいーたN}]$ (${\rm Re}[r_{ηいーたN}]$) at the precision of $\sim\pm$0.1 fm ($\sim\pm$0.5 fm), significantly narrowing down the previously estimated ranges of the parameters. To arrive at the conclusion, it is essential to use the $γがんまd\toηいーたpn$ reaction model equipped with elementary amplitudes that are well constrained by $πぱいN$ and $γがんまN$ reaction data through a sophisticated coupled-channels analysis. This result strongly motivates the Research Center for Electron Photon Science (ELPH) at Tohoku University to measure $R$. △ Less

Submitted 26 October, 2017; v1 submitted 23 April, 2017; originally announced April 2017.

Comments: 6 pages, 6 figures, published version, 1 figure added, several references added, discussion on the previous works elaborated

Report number: KEK-TH-1975, J-PARC-TH-0094

Journal ref: Phys. Rev. C 96, 042201 (2017)

Matching issue in quasi parton distribution approach

Authors: Tomomi Ishikawa, Yan-Qing Ma, Jian-Wei Qiu, Shinsuke Yoshida

Abstract: In recent years, the quasi parton distribution has been introduced for extracting the parton distribution functions from lattice QCD simulations. The quasi and standard distribution share the same perturbative collinear singularity and the renormalized quasi distribution can be factorized into the standard distribution with a perturbative matching factor. The quasi parton distribution is known to… ▽ More In recent years, the quasi parton distribution has been introduced for extracting the parton distribution functions from lattice QCD simulations. The quasi and standard distribution share the same perturbative collinear singularity and the renormalized quasi distribution can be factorized into the standard distribution with a perturbative matching factor. The quasi parton distribution is known to have power-law UV divergences, which do not exist in the standard distribution. We discuss in this talk the nonperturbative renormalization scheme for the power divergence. We also demonstrate the perturbative matching of the quasi quark distribution between continuum and lattice at the one-loop. △ Less

Submitted 25 March, 2017; originally announced March 2017.

Comments: 7 pages, 2 figures, talk presented at the 34th annual International Symposium on Lattice Field Theory (Lattice 2016), 24-30 July 2016, University of Southampton, UK

Report number: RBRC-1225

Journal ref: PoS(LATTICE2016)177

One loop effects of natural SUSY in indirect searches for SUSY particles at the ILC

Authors: M. Jimbo, T. Kon, Y. Kouda, M. Ichikawa, Y. Kurihara, T. Ishikawa, K. Kato, M. Kuroda

Abstract: We have found the possible region of parameters of the minimal supersymmetric standard model (MSSM) within the bounds from the experimental results of the Higgs mass, the rare decay mode of $b$-quark, the muon $g-2$, the dark matter abundance, and the direct searches for the lighter stop (i.e., one of the supersymmetric partners of top quark) at the LHC. We present numerical results of calculation… ▽ More We have found the possible region of parameters of the minimal supersymmetric standard model (MSSM) within the bounds from the experimental results of the Higgs mass, the rare decay mode of $b$-quark, the muon $g-2$, the dark matter abundance, and the direct searches for the lighter stop (i.e., one of the supersymmetric partners of top quark) at the LHC. We present numerical results of calculations for the one loop effects of supersymmetric particles in the processes of $τたう^+ τたう^-$, $b \overline{b}$, $t \overline{t}$, and $Z h$ production at the ILC by using benchmark points within the possible region of the MSSM parameters. △ Less

Submitted 24 May, 2017; v1 submitted 22 March, 2017; originally announced March 2017.

Comments: 11 pages, 16 figures; Talk presented at the International Workshop on Future Linear Colliders (LCWS2016), Morioka, Japan, 5-9 December 2016. C16-12-05.4

Journal ref: eConf C16-12-05.4, (2017) [arXiv:1703.07671]

Regularization with Numerical Extrapolation for Finite and UV-Divergent Multi-loop Integrals

Authors: E. de Doncker, F. Yuasa, K. Kato, T. Ishikawa, J. Kapenga, O. Olagbemi

Abstract: We give numerical integration results for Feynman loop diagrams such as those covered by Laporta [1] and by Baikov and Chetyrkin [2], and which may give rise to loop integrals with UV singularities. We explore automatic adaptive integration using multivariate techniques from the PARINT package for multivariate integration, as well as iterated in- tegration with programs from the QUADPACK package,… ▽ More We give numerical integration results for Feynman loop diagrams such as those covered by Laporta [1] and by Baikov and Chetyrkin [2], and which may give rise to loop integrals with UV singularities. We explore automatic adaptive integration using multivariate techniques from the PARINT package for multivariate integration, as well as iterated in- tegration with programs from the QUADPACK package, and a trapezoidal method based on a double exponential trans- formation. PARINT is layered over MPI (Message Passing Interface), and incorporates advanced parallel/distributed techniques including load balancing among processes that may be distributed over a cluster or a network/grid of nodes. Results are included for 2-loop vertex and box diagrams and for sets of 2-, 3- and 4-loop self-energy diagrams with or without UV terms. Numerical regularization of integrals with singular terms is achieved by linear and non-linear extrapolation methods. △ Less

Submitted 2 February, 2018; v1 submitted 16 February, 2017; originally announced February 2017.

Comments: replaced by published version

Journal ref: Comput. Phys. Comm. 224 (2018) 164

One loop effects of natural SUSY in third generation fermion production at the ILC

Authors: Yusaku Kouda, Tadashi Kon, Yoshimasa Kurihara, Tadashi Ishikawa, Masato Jimbo, Kiyoshi Kato, Masaaki Kuroda

Abstract: Within the framework of the Minimal Supersymmetric Standard Model, we investigate the 1-loop effects of supersymmetric particles on the third-generation fermion-pair production at the ILC. Three sets of the SUSY parameters are proposed which are consistent with the observed Higgs mass, the muon g-2, the Dark Matter abundance and the decay branching ratio of B meson. We discuss on the possibility o… ▽ More Within the framework of the Minimal Supersymmetric Standard Model, we investigate the 1-loop effects of supersymmetric particles on the third-generation fermion-pair production at the ILC. Three sets of the SUSY parameters are proposed which are consistent with the observed Higgs mass, the muon g-2, the Dark Matter abundance and the decay branching ratio of B meson. We discuss on the possibility of discovering the signals consistent with SUSY as well as of experimentally distinguishing the proposed sets of SUSY parameters. △ Less

Submitted 8 March, 2017; v1 submitted 26 September, 2016; originally announced September 2016.

Journal ref: Prog. Theor. Exp. Phys., 2017-5, 053B02 (2017)

Practical quasi parton distribution functions

Authors: Tomomi Ishikawa, Yan-Qing Ma, Jian-Wei Qiu, Shinsuke Yoshida

Abstract: A completely new strategy to calculate parton distribution functions on the lattice has recently been proposed. In this method, lattice calculable observables, called quasi distributions, are related to normal distributions. The quasi distributions are known to contain power-law UV divergences arise from a Wilson line in the non-local operator, while the normal distributions only have logatithmic… ▽ More A completely new strategy to calculate parton distribution functions on the lattice has recently been proposed. In this method, lattice calculable observables, called quasi distributions, are related to normal distributions. The quasi distributions are known to contain power-law UV divergences arise from a Wilson line in the non-local operator, while the normal distributions only have logatithmic UV divergences. We propose possible method to subtract the power divegence to make the matching of the quasi with the normal distributions well-defined. We also demonstrate the matching of the quasi quark distribution between continuum and lattice implementing the power divergence subtraction. The matching calculations are carried out by one-loop perturbation. △ Less

Submitted 7 September, 2016; originally announced September 2016.

Comments: 28 pages, 7 figures

Report number: RBRC-1199

A development of an accelerator board dedicated for multi-precision arithmetic operations and its application to Feynman loop integrals

Authors: Shinji Motoki, Hiroshi Daisaka, Naohito Nakasato, Tadashi Ishikawa, Fukuko Yuasa, Toshiyuki Fukushige, Atsushi Kawai, Junichiro Makino

Abstract: Higher order corrections in perturbative quantum field theory are required for precise theoretical analysis to investigate new physics beyond the Standard Model. This indicates that we need to evaluate Feynman loop diagram with multi-loop integral which may require multi-precision calculation. We developed a dedicated accelerator system for multi-precision calculation (GRAPE9-MPX). We present perf… ▽ More Higher order corrections in perturbative quantum field theory are required for precise theoretical analysis to investigate new physics beyond the Standard Model. This indicates that we need to evaluate Feynman loop diagram with multi-loop integral which may require multi-precision calculation. We developed a dedicated accelerator system for multi-precision calculation (GRAPE9-MPX). We present performance results of our system for the case of Feynman two-loop box and three-loop selfenergy diagrams with multi-precision. △ Less

Submitted 30 November, 2014; v1 submitted 13 October, 2014; originally announced October 2014.

Comments: 6 pages, 8 figures, 1 table; submitted to the proceedings of the 16th International workshop on Advanced Computing and Analysis Techniques in physics research (ACAT 2014), Prague

Neutral $B$ meson mixings and $B$ meson decay constants with static heavy and domain-wall light quarks

Authors: Yasumichi Aoki, Tomomi Ishikawa, Taku Izubuchi, Christoph Lehner, Amarjit Soni

Abstract: Neutral $B$ meson mixing matrix elements and $B$ meson decay constants are calculated. Static approximation is used for $b$ quark and domain-wall fermion formalism is employed for light quarks. The calculations are carried out on $2+1$ flavor dynamical ensembles generated by RBC/UKQCD Collaborations with lattice spacings $0.086$fm ($a^{-1}\sim 2.3$GeV) and $0.11$fm ($1.7$GeV), and a fixed physical… ▽ More Neutral $B$ meson mixing matrix elements and $B$ meson decay constants are calculated. Static approximation is used for $b$ quark and domain-wall fermion formalism is employed for light quarks. The calculations are carried out on $2+1$ flavor dynamical ensembles generated by RBC/UKQCD Collaborations with lattice spacings $0.086$fm ($a^{-1}\sim 2.3$GeV) and $0.11$fm ($1.7$GeV), and a fixed physical spatial volume of about $(2.7{\rm fm})^3$. In the static quark action, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings, HYP1 and HYP2, and their results are combined in taking the continuum limit. For the matching between the lattice and the continuum theory, one-loop perturbative $O(a)$ improvements are made to reduce discretization errors. As the most important quantity of this work, we obtain SU(3) breaking ratio $ξくしー=1.208(60)$, where the error includes statistical and systematic one. (Uncertainty from infinite $b$ quark mass is not included.) We also find other neutral $B$ meson mixing quantities $f_B\sqrt{\hat{B}_B}=240(22)$MeV, $f_{B_s}\sqrt{\hat{B}_{B_s}}=290(22)$MeV, $\hat{B}_B=1.17(22)$, $\hat{B}_{B_s}=1.22(13)$ and $B_{B_s}/B_B=1.028(74)$, $B$ meson decay constants $f_B=219(17)$MeV, $f_{B_s}=264(19)$MeV and $f_{B_s}/f_B=1.193(41)$, in the static limit of $b$ quark, which do not include infinite $b$ quark mass uncertainty. △ Less

Submitted 29 June, 2015; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: 33 pages, 18 figures, v3: published version

Report number: RBRC-1080

Journal ref: Phys. Rev. D 91, 114505 (2015)

Full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to $e^+e^- \rightarrow e^+e^- γがんま$ at the ILC with GRACE-Loop

Authors: P. H. Khiem, Y. Kurihara, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, N. Nakazawa, Y. Shimizu, T. Ueda, J. A. M. Vermaseren, Y. Yasui

Abstract: By using the GRACE-Loop system, we calculate the full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to the process $e^+e^- \rightarrow e^+e^- γがんま$, which is important for future investigations at the International Linear Collider (ILC). With the GRACE-Loop system, the calculations are checked numerically by three consistency tests: ultraviolet finiteness, infrared finiteness, and gauge-paramete… ▽ More By using the GRACE-Loop system, we calculate the full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to the process $e^+e^- \rightarrow e^+e^- γがんま$, which is important for future investigations at the International Linear Collider (ILC). With the GRACE-Loop system, the calculations are checked numerically by three consistency tests: ultraviolet finiteness, infrared finiteness, and gauge-parameter independence. The results show good numerical stability when quadruple precision is used. In the phenomenological results, we find that the electroweak corrections to the total cross section range from $\sim -4\%$ to $\sim -21\%$ when $\sqrt{s}$ varies from $250$ GeV to $1$ TeV. The corrections also significantly affect the differential cross sections, which are a function of the invariant masses and angles and the final-particle energies. Such corrections will play an important role for the high-precision program at the ILC. △ Less

Submitted 27 August, 2014; v1 submitted 25 March, 2014; originally announced March 2014.

Comments: 15 pages, 7 figures, Talk presented at the International Workshop on Future Linear Colliders (LCWS13), Tokyo, Japan, 11-15 November 2013

Full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to $t \bar{t} γがんま$ and $e^-e^+γがんま$ productions at ILC with GRACE-Loop

Authors: P. H. Khiem, Y. Kurihara, J. Fujimoto, M. Igarashi, T. Ishikawa, T. Kaneko, K. Kato, N. Nakazawa, Y. Shimizu, K. Tobimatsu, T. Ueda, J. A. M. Vermaseren, Y. Yasui

Abstract: The full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to $t \bar{t} γがんま$ and $e^-e^+γがんま$ productions at the International Linear Collider (ILC) are presented in this paper. The computation is performed with the help of GRACE-Loop system. In the physical results, we discuss on the cross section, electroweak corrections, and the top quark forward-backward asymmetry ($A_{FB}$) which are the functio… ▽ More The full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to $t \bar{t} γがんま$ and $e^-e^+γがんま$ productions at the International Linear Collider (ILC) are presented in this paper. The computation is performed with the help of GRACE-Loop system. In the physical results, we discuss on the cross section, electroweak corrections, and the top quark forward-backward asymmetry ($A_{FB}$) which are the function of the center-of-mass energy. △ Less

Submitted 25 March, 2014; originally announced March 2014.

Comments: 4 pages, Proceedings of ICISE2013 - Rencontres du Vietnam, Windows on the Universe. 11-17 Aug 2013. Quy Nhon, Vietnam

Neutral B meson mixing with static heavy and domain-wall light quarks

Authors: Tomomi Ishikawa, Yasumichi Aoki, Taku Izubuchi, Christoph Lehner, Amarjit Soni

Abstract: Neutral B meson mixing matrix elements and B meson decay constants are calculated. Static approximation is used for b quark and domain-wall fermion formalism is employed for light quarks. The calculations are done on 2+1 flavor dynamical ensembles, whose lattice spacings are 0.086 fm and 0.11 fm with a fixed physical spatial volume of about (2.7 fm)^3. In the static quark action, link-smearings ar… ▽ More Neutral B meson mixing matrix elements and B meson decay constants are calculated. Static approximation is used for b quark and domain-wall fermion formalism is employed for light quarks. The calculations are done on 2+1 flavor dynamical ensembles, whose lattice spacings are 0.086 fm and 0.11 fm with a fixed physical spatial volume of about (2.7 fm)^3. In the static quark action, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings and their results are combined in taking a continuum limit. For the matching between the lattice and the continuum theory, one-loop perturbative calculations are used including O(a) improvements to reduce discretization errors. We obtain SU(3) braking ratio ξくしー=1.222(60) in the static limit of b quark. △ Less

Submitted 3 December, 2013; originally announced December 2013.

Comments: 7 pages, 3 figures, talk presented at 31st International Symposium on Lattice Field Theory (LATTICE 2013), July 29 - August 3, 2013, Mainz, Germany

Report number: RBRC-1052

Full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to $e^+e^- \rightarrow t \bar{t} γがんま$ with GRACE-Loop

Authors: P. H. Khiem, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. Shimizu, T. Ueda, J. A. M. Vermaseren, Y. Yasui

Abstract: We present the full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to the process $e^+e^- \rightarrow t \bar{t} γがんま$ at the International Linear Collider (ILC). The computation is performed with the help of the GRACE-Loop system. We present the total cross-section and the top quark forward-backward asymmetry ($A_{FB}$) as a function of the center-of-mass energy and compare them with the process… ▽ More We present the full $\mathcal{O}(αあるふぁ)$ electroweak radiative corrections to the process $e^+e^- \rightarrow t \bar{t} γがんま$ at the International Linear Collider (ILC). The computation is performed with the help of the GRACE-Loop system. We present the total cross-section and the top quark forward-backward asymmetry ($A_{FB}$) as a function of the center-of-mass energy and compare them with the process $e^+e^- \rightarrow t \bar{t}$. We find that the value of $A_{FB}$ in $t \bar{t} γがんま$ production is larger than $A_{FB}$ in $t\bar{t}$ production. It is an important result for the measurement of the top quark forward-backward asymmetry at the ILC. Applying a structure function method, we also subtract the QED correction to gain the genuine weak correction in both the $αあるふぁ$ scheme and the $G_μみゅー$ scheme ($δでるた_{W}^{G_μみゅー}$). We obtain numerical values for $δでるた_{W}^{G_μみゅー}$ which are changing from 2% to -24% when we vary the center-of-mass energy from 360 GeV to 1 TeV. △ Less

Submitted 1 May, 2013; v1 submitted 6 November, 2012; originally announced November 2012.

Comments: 13 pages, 9 figures

Pair-production and Three-body Decay of the Lighter Stop at the ILC in One-loop Order in the MSSM

Authors: M. Jimbo, T. Inoue, T. Jujo, T. Kon, T. Ishikawa, Y. Kurihara, K. Kato, M. Kuroda

Abstract: We have been developing a program package called GRACE/SUSY-loop which automatically calculates the MSSM amplitudes in one-loop order. We present numerical results of calculations for pair-production and three-body decay of the lighter stop ($\widetilde{t}_1$) at the International Linear Collider (ILC) using GRACE/SUSY-loop. Since the distributions of missing transverse energy (MET) depend on mass… ▽ More We have been developing a program package called GRACE/SUSY-loop which automatically calculates the MSSM amplitudes in one-loop order. We present numerical results of calculations for pair-production and three-body decay of the lighter stop ($\widetilde{t}_1$) at the International Linear Collider (ILC) using GRACE/SUSY-loop. Since the distributions of missing transverse energy (MET) depend on mass spectrum of SUSY particles, we consider two scenarios on three-body decay of $\widetilde{t}_1$. In these scenarios, both QCD and EW corrections have positive sign for decay widths and cross sections. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Comments: 6 pages, 5 figures (7 eps files), Prepared for the Proceedings of International Workshop on Future Linear Colliders, LCWS11, 26-30 September 2011, Palacio de Exposiciones y Congresos, Granada, Spain

Report number: CUCP-12-1, KEK-CP-265

Full QED+QCD Low-Energy Constants through Reweighting

Authors: Tomomi Ishikawa, Thomas Blum, Masashi Hayakawa, Taku Izubuchi, Chulwoo Jung, Ran Zhou

Abstract: The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor… ▽ More The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor dynamical domain-wall quark ensembles. △ Less

Submitted 17 August, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: 5 pages, 9 figures, REVTeX 4.1, v2: published version

Report number: RBRC-945

Journal ref: Phys. Rev. Lett. 109, 072002 (2012)

Numerical approach to multi-loop integrals

Authors: K. Kato, E. de Doncker, N. Hamaguchi, T. Ishikawa, T. Koike, Y. Kurihara, Y. Shimizu, F. Yuasa

Abstract: For the calculation of multi-loop Feynman integrals, a novel numerical method, the Direct Computation Method (DCM) is developed. It is a combination of a numerical integration and a series extrapolation. In principle, DCM can handle diagrams of arbitrary internal masses and external momenta, and can calculate integrals with general numerator function. As an example of the performance of DCM, a num… ▽ More For the calculation of multi-loop Feynman integrals, a novel numerical method, the Direct Computation Method (DCM) is developed. It is a combination of a numerical integration and a series extrapolation. In principle, DCM can handle diagrams of arbitrary internal masses and external momenta, and can calculate integrals with general numerator function. As an example of the performance of DCM, a numerical computation of two-loop box diagrams is presented. Further discussion is given on the choice of control parameters in DCM. This method will be an indispensable tool for the higher order radiative correction when it is tested for a wider class of physical parameters and the separation of divergence is done automatically. △ Less

Submitted 30 January, 2012; originally announced January 2012.

Comments: 7 pages; The XXth International Workshop High Energy Physics and Quantum Field Theory, September 24-October 1, 2011, Sochi Russia

Numerical Computation of Two-loop Box Diagrams with Masses

Authors: F. Yuasa, E. de Doncker, N. Hamaguchi, T. Ishikawa, K. Kato, Y. Kurihara, J. Fujimoto, Y. Shimizu

Abstract: A new approach is presented to evaluate multi-loop integrals, which appear in the calculation of cross-sections in high-energy physics. It relies on a fully numerical method and is applicable to a wide class of integrals with various mass configurations. As an example, the computation of two-loop planar and non-planar box diagrams is shown. The results are confirmed by comparisons with other techn… ▽ More A new approach is presented to evaluate multi-loop integrals, which appear in the calculation of cross-sections in high-energy physics. It relies on a fully numerical method and is applicable to a wide class of integrals with various mass configurations. As an example, the computation of two-loop planar and non-planar box diagrams is shown. The results are confirmed by comparisons with other techniques, including the reduction method, and by a consistency check using the dispersion relation. △ Less

Submitted 25 May, 2012; v1 submitted 3 December, 2011; originally announced December 2011.

Comments: 16 pages, 8 figures

Report number: KEK-CP-263

Extrapolation Algorithms for Infrared Divergent Integrals

Authors: Elise de Doncker, Junpei Fujimoto, Nobuyuki Hamaguchi, Tadashi Ishikawa, Yoshimasa Kurihara, Marko Ljucovic, Yoshimitsu Shimizu, Fukuko Yuasa

Abstract: This paper describes applications of extrapolation for the computation of coefficients in an expansion of infrared divergent integrals. An extrapolation procedure is performed with respect to a parameter introduced by dimensional regularization. While this treats typical IR singularities at the boundaries of the integration domain, special care needs to be taken in cases where the integrand is sin… ▽ More This paper describes applications of extrapolation for the computation of coefficients in an expansion of infrared divergent integrals. An extrapolation procedure is performed with respect to a parameter introduced by dimensional regularization. While this treats typical IR singularities at the boundaries of the integration domain, special care needs to be taken in cases where the integrand is singular in the interior of the domain as well as on the boundaries. A double extrapolation is devised for a class of massless vertex integrals. Quadruple precision results are presented, demonstrating high accuracy. The computations are supported by the use of general adaptive integration programs from the QUADPACK package, in iterated integrations with highly singular integrand functions. △ Less

Submitted 17 October, 2011; originally announced October 2011.

Comments: 11 pages, 1 figure, 6 tables, Talk given at 3rd Computational Particle Physics Workshop CPP2010, September 23-25, 2010, KEK Japan

Numerical Approach to Calculation of Feynman Loop Integrals

Authors: F. Yuasa, T. Ishikawa, Y. Kurihara, J. Fujimoto, Y. Shimizu, N. Hamaguchi, E. de Doncker, K. Kato

Abstract: In this paper, we describe a numerical approach to evaluate Feynman loop integrals. In this approach the key technique is a combination of a numerical integration method and a numerical extrapolation method. Since the computation is carried out in a fully numerical way, our approach is applicable to one-, two- and multi-loop diagrams. Without any analytic treatment it can compute diagrams with not… ▽ More In this paper, we describe a numerical approach to evaluate Feynman loop integrals. In this approach the key technique is a combination of a numerical integration method and a numerical extrapolation method. Since the computation is carried out in a fully numerical way, our approach is applicable to one-, two- and multi-loop diagrams. Without any analytic treatment it can compute diagrams with not only real masses but also complex masses for the internal particles. As concrete examples we present numerical results of a scalar one-loop box integral with complex masses and two-loop planar and non-planar box integrals with masses. We discuss the quality of our numerical computation by comparisons with other methods and also propose a self consistency check. △ Less

Submitted 20 September, 2011; originally announced September 2011.

Comments: 8 pages, 5 figures, Talk given at 3rd Computational Particle Physics Workshop CPP2010, September 23-25, 2010, KEK Japan

Slepton NLG (Non-Linear Gauge) in GRACE/SUSY-loop

Authors: Masato Jimbo, Tadashi Ishikawa, Masaaki Kuroda

Abstract: We have been developing a program package called GRACE/SUSY-loop which is for the automatic calculations of the MSSM amplitudes in one-loop order. The non-linear gauge (NLG) fixing conditions play the crucial role in the calculations in one-loop order which contain a large number of Feynman diagrams. We present the recent progress in GRACE/SUSY-loop which is obtained by extending the non-linear ga… ▽ More We have been developing a program package called GRACE/SUSY-loop which is for the automatic calculations of the MSSM amplitudes in one-loop order. The non-linear gauge (NLG) fixing conditions play the crucial role in the calculations in one-loop order which contain a large number of Feynman diagrams. We present the recent progress in GRACE/SUSY-loop which is obtained by extending the non-linear gauge formalism to the slepton sector. △ Less

Submitted 5 December, 2010; v1 submitted 30 November, 2010; originally announced November 2010.

Comments: 8 pages, 3 figures (5 eps files), To appear in Proceedings of CPP2010, 3rd Computational Particle Physics Workshop, 23rd-25th September 2010 at KEK (Tsukuba)

Report number: CUCP-10-3, KEK-CP-248

Present Status of GRACE/SUSY-loop

Authors: M. Jimbo, K. Iizuka, T. Ishikawa, K. Kato, T. Kon, Y. Kurihara, M. Kuroda

Abstract: GRACE/SUSY-loop is a program package for the automatic calculation of the MSSM amplitudes in one-loop order. We present features of GRACE/SUSY-loop, processes calculated using GRACE/SUSY-loop and an extension of the non-linear gauge formalism applied to GRACE/SUSY-loop. GRACE/SUSY-loop is a program package for the automatic calculation of the MSSM amplitudes in one-loop order. We present features of GRACE/SUSY-loop, processes calculated using GRACE/SUSY-loop and an extension of the non-linear gauge formalism applied to GRACE/SUSY-loop. △ Less

Submitted 17 June, 2010; originally announced June 2010.

Comments: 7 pages, 2 figures, Talk presented at LCWS10/ILC10

Report number: CUCP-10-2, KEK-CP-236, KU-PH-006

Systematic study of 1-loop correction on sparticle decay widths using GRACE/SUSY-loop

Authors: K. Iizuka, T. Kon, K. Kato, T. Ishikawa, Y. Kurihara, M. Jimbo, M. Kuroda

Abstract: The 1-loop corrected decay widths of sparticles (charginos, neutralinos, gluino and sfermions) in the framework of the MSSM are calculated systematically using GRACE/SUSY-loop, which is the program package for the automatic calculation of the MSSM amplitudes in the 1-loop order. We present the renormalization scheme used in our system and show some numerical results of decay widths of sfermions… ▽ More The 1-loop corrected decay widths of sparticles (charginos, neutralinos, gluino and sfermions) in the framework of the MSSM are calculated systematically using GRACE/SUSY-loop, which is the program package for the automatic calculation of the MSSM amplitudes in the 1-loop order. We present the renormalization scheme used in our system and show some numerical results of decay widths of sfermions and gluino using the SPS1a' parameter set and other SUSY parameter sets. △ Less

Submitted 25 January, 2010; v1 submitted 15 January, 2010; originally announced January 2010.

Comments: 6 pages, 4 figures, added references, contribution to the proceedings of the 9th International Symposium on Radiative Corrections (RADCOR 2009), October 25-30 2009, Ascona, Switzerland

Report number: KEK-CP-228, CUCP-10-1, KU-PH-004

Journal ref: PoS RADCOR2009:068,2010

Numerical Evaluation of Feynman Integrals by a Direct Computation Method

Authors: F. Yuasa, T. Ishikawa, J. Fujimoto, N. Hamaguchi, E. de Doncker, Y. Shimizu

Abstract: A purely numerical method, Direct ComputationMethod is applied to evaluate Feynman integrals. This method is based on the combination of an efficient numerical integration and an efficient extrapolation. In addition, high-precision arithmetic and parallelization technique can be used in this method if required. We present the recent progress in development of this method and show results such as… ▽ More A purely numerical method, Direct ComputationMethod is applied to evaluate Feynman integrals. This method is based on the combination of an efficient numerical integration and an efficient extrapolation. In addition, high-precision arithmetic and parallelization technique can be used in this method if required. We present the recent progress in development of this method and show results such as one-loop 5-point and two-loop 3-point integrals. △ Less

Submitted 18 April, 2009; originally announced April 2009.

Comments: 5 pages,9 figures, LaTeX, PoS style; ACAT2008, XII International Workshop on Advanced Computing and Analysis Techniques in Physics Research, Erice, Sicily, Italy November 3-7, 2008

Report number: KEK-CP-222

Journal ref: PoS ACAT08:122,2008

Status reports from the GRACE Group

Authors: Y. Yasui, T. Ueda, E. de Doncker, J. Fujimoto, N. Hamaguchi, T. Ishikawa, Y. Shimizu, F. Yuasa

Abstract: We discuss a new approach for the numerical evaluation of loop integrals. The fully numerical calculations of an infrared one-loop vertex and a box diagram are demonstrated. To perform these calculations, we apply an extrapolation method based on the $εいぷしろん$-algorithm. In our approach, the super high precision control in the numerical manipulation is essential to handle the infrared singularity. We… ▽ More We discuss a new approach for the numerical evaluation of loop integrals. The fully numerical calculations of an infrared one-loop vertex and a box diagram are demonstrated. To perform these calculations, we apply an extrapolation method based on the $εいぷしろん$-algorithm. In our approach, the super high precision control in the numerical manipulation is essential to handle the infrared singularity. We adopt a multi-precision library named {\tt HMLib} for the precision control in the calculations. △ Less

Submitted 16 October, 2007; originally announced October 2007.

Comments: 6 pages, 3 figures, Proceedings for the International Workshop LCWS/ILC 2007, DESY in Hamburg, Germany, between May 30, 2007, and June 3, 2007

Journal ref: ECONF C0705302:LOOP04,2007

Precise Numerical Evaluation of the Scalar One-Loop Integrals with the Infrared Divergence

Authors: F. Yuasa, E. de Doncker, J. Fujimoto, N. Hamaguchi, T. Ishikawa, Y. Shimizu

Abstract: We present a new approach for obtaining very precise integration results for infrared vertex and box diagrams, where the integration is carried out directly without performing any analytic integration of Feynman parameters. Using an appropriate numerical integration routine with an extrapolation method, together with a multi-precision library, we have obtained integration results which agree wit… ▽ More We present a new approach for obtaining very precise integration results for infrared vertex and box diagrams, where the integration is carried out directly without performing any analytic integration of Feynman parameters. Using an appropriate numerical integration routine with an extrapolation method, together with a multi-precision library, we have obtained integration results which agree with the analytic results to 10 digits even for such a very small photon mass as $10^{-150}$ GeV in the infrared vertex diagram. △ Less

Submitted 19 September, 2007; v1 submitted 6 September, 2007; originally announced September 2007.

Comments: 10 pages, 2 figures, contribution to the proceedings of the ACAT workshop, Amsterdam, April 23-27, 2007

Report number: KEK-CP-196

Journal ref: PoSACAT2007:087,2007

Light quark masses from unquenched lattice QCD

Authors: CP-PACS/JLQCD Collaborations, :, T. Ishikawa, S. Aoki, M. Fukugita, S. Hashimoto, K-I. Ishikawa, N. Ishizuka, Y. Iwasaki, K. Kanaya, T. Kaneko, Y. Kuramashi, M. Okawa, Y. Taniguchi, N. Tsutsui, A. Ukawa, N. Yamada, T. Yoshié

Abstract: We calculate the light meson spectrum and the light quark masses by lattice QCD simulation, treating all light quarks dynamically and employing the Iwasaki gluon action and the nonperturbatively O(a)-improved Wilson quark action. The calculations are made at the squared lattice spacings at an equal distance a^2~0.005, 0.01 and 0.015 fm^2, and the continuum limit is taken assuming an O(a^2) discr… ▽ More We calculate the light meson spectrum and the light quark masses by lattice QCD simulation, treating all light quarks dynamically and employing the Iwasaki gluon action and the nonperturbatively O(a)-improved Wilson quark action. The calculations are made at the squared lattice spacings at an equal distance a^2~0.005, 0.01 and 0.015 fm^2, and the continuum limit is taken assuming an O(a^2) discretization error. The light meson spectrum is consistent with experiment. The up, down and strange quark masses in the \bar{MS} scheme at 2 GeV are \bar{m}=(m_{u}+m_{d})/2=3.55^{+0.65}_{-0.28} MeV and m_s=90.1^{+17.2}_{-6.1} MeV where the error includes statistical and all systematic errors added in quadrature. These values contain the previous estimates obtained with the dynamical u and d quarks within the error. △ Less

Submitted 21 July, 2008; v1 submitted 16 April, 2007; originally announced April 2007.

Comments: 4 pages, 3 figures, revtex4; v2: contents partly modified, published version

Report number: UTCCS-P-31, UTHEP-542, KEK-CP-193, HUPD-0702, RBRC-666

Journal ref: Phys.Rev.D78:011502,2008