BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ken-Ichi Ishikawa

Ken-Ichi Ishikawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phhep-thnucl-exnucl-thhep-exmath.NAMathematical SoftwareNumerical Analysisphysics.comp-ph

Catalog footprint

What is connected

22works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Nucleon isovector tensor charge from lattice QCD with physical light quarks

We present preliminary results for the axial, scalar and tensor charges of the nucleon measured in 2+1 flavor QCD with the physical light quarks ($m_π=135$ MeV). Our simulations are carried out with gauge configurations generated by the PACS Collaboration with the stout-smeared $O(a)$ improved Wilson fermions and Iwasaki gauge action at a single lattice spacing of $0.085\ (\mathrm{fm})$. There are two lattice ensembles of the PACS gauge configurations, which have physical lattice sizes over $(10\ \mathrm{fm})^4$ and $(5\ \mathrm{fm})^4$, respectively. We compute the nucleon three-point correlation functions in the axial, scalar, and tensor channels. For the renormalization, we use the Rome-Southampton method as the intermediate scheme in order to evaluate the renormalization constants for the scalar and tensor currents in fully nonperturbative manner. We then evaluate the renormalized values of the scalar and tensor charges ($g_S$ and $g_T$) in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV with a help of the continuum perturbation theory for the matching between two schemes. We compare our preliminary results of $g_S$ and $g_T$ with those of other collaboration results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Scale setting for large-$N$ SUSY Yang-Mills on the lattice

In this paper we study the large $N$ limit of four-dimensional Supersymmetric Yang-Mills on the lattice using twisted reduced models. We have generated configurations with dynamical massive gluinos and various lattice 't Hooft couplings, and verified that the Pfaffian remains positive. We have determined the lattice spacing in terms of various observables obtaining compatible results. Extrapolating to the massless gluino limit we obtain the lattice spacing dependence on the bare couplings for the supersymmetric theory. The observed dependence goes along the expected behaviour predicted by perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The lower moments of nucleon structure functions in lattice QCD with physical quark masses

We present results for the nucleon structure functions and form factors obtained from 2+1 flavor lattice QCD with physical light quark masses ($m_π=135$ MeV) in a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared ${\mathscr{O}}(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β=1.82$ and $2.00$ corresponding to lattice spacings of $0.085$ fm and $0.064$ fm respectively. The lower moments of structure functions, $\langle x \rangle_{u-d}$ and $\langle x \rangle_{Δu - Δd}$ given by the twist-2 operators being properly renormalized, are evaluated in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV only at $β=1.82$, since the renormalization factors at $β=2.00$ have not yet determined nonperturbatively in the RI/MOM scheme. Instead, at two lattice spacings, we evaluate appropriate ratios of $g_{A}/g_{V}$ and $\langle x \rangle_{u-d}/\langle x \rangle_{Δu -Δd}$, which are not renormalized in the continuum limit. These quantities thus can be directly compared with the experimental data without the renormalization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$K_{l3}$ form factors at the physical point on (10.9 fm)$^3$ volume

We present the calculation of the $K_{l3}$ form factors with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on a large volume of (10.9 fm)$^3$ at one lattice spacing of $a = 0.085$ fm. We extract the form factors from 3-point functions with three different time separations between the source and sink operators to confirm suppression of excited state contributions. The form factors are calculated in very close to the zero momentum transfer, $q^2 = 0$, thanks to the large volume, so that stable interpolations to $q^2 = 0$ are carried out. Using our form factors, we obtain the form factor at $q^2 = 0$, $f_+(0) = 0.9603(16)(^{+14}_{\ -4})(44)(19)(1)$, where the first, second, and fifth errors are statistical, systematic errors from fit functions and the isospin breaking effect, respectively. The third and fourth errors denote the finite lattice spacing effects estimated from the renormalization factor and contribution beyond the leading order SU(3) chiral perturbation theory (ChPT). The result of $f_+(0)$ yields the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, $|V_{us}| = 0.2255(13)(4)$, where the first error comes from our calculation and the second from the experiment. This value is consistent with the ones determined from the unitarity of the CKM matrix and the $K_{l2}$ decay within one standard deviation, while it is slightly larger than recent lattice calculations by at most 1.5 $σ$. Furthermore, we evaluate the shape of the form factors and the phase space integral from our results. We confirm that those results are consistent with the experiment, and also $|V_{us}|$ determined with our phase space integral agrees with the one in the above.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$K_{l3}$ form factors in $N_f = 2+1$ QCD at physical point on large volume

We present our results of the $K_{l3}$ form factors on the volume whose spatial extent is more than $L=$10 fm, with the physical pion and kaon masses using the stout-smearing clover $N_f = 2+1$ quark action and Iwasaki gauge action at $a^{-1}\approx2.3$ GeV. The $K_{l3}$ form factor at zero momentum transfer is obtained from fit based on the next-to-leading (NLO) formula in SU(3) chiral perturbation theory. We estimate systematic errors of the form factor, mainly coming from the finite lattice spacing effect. We also determine the value of $|V_{us}|$ by combining our result with the experiment and check the consistency with the standard model prediction. The result is compared with the previous lattice calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A thick-restart Lanczos type method for Hermitian $J$-symmetric eigenvalue problems

A thick-restart Lanczos type algorithm is proposed for Hermitian $J$-symmetric matrices. Since Hermitian $J$-symmetric matrices possess doubly degenerate spectra or doubly multiple eigenvalues with a simple relation between the degenerate eigenvectors, we can improve the convergence of the Lanczos algorithm by restricting the search space of the Krylov subspace to that spanned by one of each pair of the degenerate eigenvector pairs. We show that the Lanczos iteration is compatible with the $J$-symmetry, so that the subspace can be split into two subspaces that are orthogonal to each other. The proposed algorithm searches for eigenvectors in one of the two subspaces without the multiplicity. The other eigenvectors paired to them can be easily reconstructed with the simple relation from the $J$-symmetry. We test our algorithm on randomly generated small dense matrices and a sparse large matrix originating from a quantum field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon form factors and root-mean-square radii on a (10.8 fm$)^4$ lattice at the physical point

We present the nucleon form factors and root-mean-square (RMS) radii measured on a (10.8 fm$)^4$ lattice at the physical point. We compute the form factors at small momentum transfer region in $q^2\le 0.102$ GeV$^2$ with the standard plateau method choosing four source-sink separation times $t_{\rm sep}$ from 0.84 to 1.35 fm to examine the possible excited state contamination. We obtain the electric and magnetic form factors and their RMS radii for not only the isovector channel but also the proton and neutron ones without the disconnected diagram. We also obtain the axial-vector coupling and the axial radius from the axial-vector form factor. We find that these three form factors do not show large $t_{\rm sep}$ dependence in our lattice setup. On the other hand, the induced pseudoscalar and pseudoscalar form factors show the clear effects of the excited state contamination, which affect the generalized Goldberger-Treiman relation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Towards higher order numerical stochastic perturbation computation applied to the twisted Eguchi-Kawai model

We have evaluated perturbation coefficients of Wilson loops up to $O(g^8)$ for the four-dimensional twisted Eguchi-Kawai model using the numerical stochastic perturbation theory (NSPT) in arXiv:1902.09847. In this talk we present a progress report on the higher order calculation up to $O(g^{63})$, for which we apply a fast Fourier transformation (FFT) based convolution algorithm to the multiplication of polynomial matrices in the NSPT aiming for higher order calculation. We compare two implementations with the CPU-only version and the GPU version of the FFT based convolution algorithm, and find a factor 9 improvement on the computational speed of the NSPT algorithm with SU($N=225$) at $O(g^{31})$. The perturbation order dependence of the computational time, we investigate it up to $O(g^{63})$, shows a mild scaling behavior on the truncation order.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A performance evaluation of CCS QCD Benchmark on the COMA (Intel(R) Xeon Phi$^{TM}$, KNC) system

The most computationally demanding part of Lattice QCD simulations is solving quark propagators. Quark propagators are typically obtained with a linear equation solver utilizing HPC machines. The CCS QCD Benchmark is a benchmark program solving the Wilson-Clover quark propagator, and is developed at the Center for Computational Sciences (CCS), University of Tsukuba. We optimized the benchmark program for a \Intel \XeonPhi (Knights Corner, KNC) system named "COMA (PACS-IX)" at CCS Tsukuba under the Intel Parallel Computing Center program. A single precision BiCGStab solver with the overlapped Restricted Additive Schwarz (RAS) preconditioner was implemented using SIMD intrinsics, OpenMP and MPI in the offload mode. With the reverse-offloading technique, we could reduce the communication and offloading overheads. We observed a performance of $\sim 200$ GFlops sustained for the Wilson-Clover hopping matrix multiplication on the lattice sizes larger than $24^3\times 32$ on a sinlge card of the COMA system. A good weak scaling perofmace was observed on the local lattice sizes larger than $24^3\times 32$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical determination of the $Λ$-parameter in SU(3) gauge theory from the twisted gradient flow coupling

We estimate the $Λ$-parameter in the $\overline{\mathrm{MS}}$ scheme for the SU(3) pure gauge theory with the twisted gradient flow method non-perturbatively. We obtain $Λ_{\overline{\mathrm{MS}}}/\sqrtσ=0.527(13)(10)$ and $r_{0}Λ_{\overline{\mathrm{MS}}}=0.605(15)(5)$ which are consistent with the known values.This demonstrates the validity of the present method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Study of quark mass dependence of binding energy for light nuclei in 2+1 flavor lattice QCD

We investigate the formation of light nuclei with the nuclear mass number less than or equal to four in 2+1 flavor QCD using a non-perturbative improved Wilson quark and Iwasaki gauge actions. The quark mass is decreased from our previous work to the one corresponding to the pion mass of 0.30 GeV. In each multi-nucleon channel, the energy shift of the ground state relative to the assembly of free nucleons is calculated on two volumes, whose spatial extents are 4.3 fm and 5.8 fm. From the volume dependence of the energy shift, we distinguish a bound state of multi nucleons from an attractive scattering state. We find that all the ground states measured in this calculation are bound states. As in the previous studies at larger $m_π$, our result indicates that at $m_π= 0.30$ GeV the effective interaction between nucleons in the light nuclei is relatively stronger than the one in nature, since the results for the binding energies are larger than the experimental values and a bound state appears in the dineutron channel, which is not observed in experiment. Possible sources of systematic error in this calculation are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A construction of the Schrödinger Functional for Möbius Domain Wall Fermions

We construct the Schrödinger Functional (SF) setup for the Möbius domain wall fermions (MDWF). The method is an extension of the method proposed by Takeda for the standard domain wall fermion. In order to fulfill the requirement that the lattice Dirac operator with the SF boundary obeys the Lüscher's universality argument: the lattice chiral fermion with the SF boundary condition breaks the chiral symmetry at the temporal boundary, we impose the parity symmetry with respect to the fifth-direction on the MDWF operator. This additional symmetry restricts the choice of the parameter of the MDWF so that the optimal parameter from the Zolotarev optimal approximation cannot be applied. We introduce a modified parameter set having the fifth-dimensional parity symmetry. We investigate the MDWF with the SF boundary by observing eigenvalues of the Hermitian operator and the Ginsparg-Wilson relation violation at the tree-level. We compare the computational cost with that of the standard DWF with the SF scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical evaluation of vacuum polarization tensor in constant external magnetic fields

Hattori-Itakura have recently derived the full Landau-level summation form for the photon vacuum polarization tensor in constant external magnetic fields at the one-loop level. The Landau-level summation form is essential when the photon momentum exceeds the threshold of the pair creation of charged particles in a magnetic field stronger than the squared mass of the charged particle. The tensor has three different form factors depending on the tensor direction with respect to the external magnetic field. The renormalization is nontrivial because these form factors are expressed in terms of double or triple summation forms. We give a numerical UV subtraction method which can be applied to numerically evaluate the form factors in constant external magnetic fields. We numerically investigate the photon vacuum polarization tensor in the form of the Landau-level summation and estimate the systematic errors coming from truncation of the Landau-level summation in a parameter region realized in heavy ion collision experiments. We find that the error is practically controllable at an $O(10^{-2})$ level for electrons and muons in strong magnetic fields expected in heavy ion collisions in the experimentally feasible kinematic parameter regions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice study on chiral dynamics of two-color six-flavors QCD

The electroweak symmetry breaking and origin of masses may be attributed to the breakdown of chiral symmetry due to a strong gauge dynamics. Among several candidates of such gauge systems, we focus on two-color QCD with $N_f=6$ massless Dirac fermions in the fundamental representation, and study on whether the dynamics of this gauge system trigger chiral symmetry breaking or not by simulating with Wilson fermions on lattices up to $L/a=32$. We show the result for the quark mass dependence and the volume dependence of some quantities such as the mass of the lightest pseudoscalar meson, decay constant, and give the three evidences supporting the absence of chiral symmetry breaking in the six-flavor theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multi-nucleon bound states in $N_f=2+1$ lattice QCD

We report on our on-going effort to calculate the properties of light nuclei directly from quarks and gluons based on lattice QCD. After briefly introducing our motivations and aims, we describe our strategy of fixing the strange quark mass at its physical value and approaching the physical point for the up and down quark masses step by step from the region of heavy quark masses. A successful calculation for the pion mass of 0.51GeV is reviewed, and the status for a lighter pion mass of 0.30GeV is reported.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing reweighting method for truncated Overlap fermions

The lattice QCD simulation with the lattice chiral symmetry is very attractive, however, it is difficult to maintain the symmetry at a modest numerical computation cost. A candidate to reduce the computational cost during the configuration generation with the HMC algorithm is to relax the requirement of the chiral symmetry and to use the reweighing method recovering the symmetry at the measurement phase. In this talk, we presented the reweighing method to restore the chiral symmetry of the truncated overlap fermion operator. In order to avoid the large discrepancy between the truncated overlap operator and the exact overlap operator, we split the reweighting factor into several steps gradually increasing the order of truncation. We investigated the truncation dependence of the reweighting factor on a set of quenched $8^3\times 32$ lattice configurations generated with the DBW2 gauge action. We found that a large fluctuation on the reweighting factor between a high-order truncated overlap operator and the exact overlap operator on a couple of configurations. The origin of the large fluctuation seems to be due to a small eigenvalue of the overlap kernel on these configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bound states of multi-nucleon channels in N_f=2+1 lattice QCD

We calculate the energies for multi-nucleon ground states with the nuclear mass number less than or equal to 4 in 2+1 flavor QCD at the lattice spacing of a = 0.09 fm employing a relatively heavy quark mass corresponding to m_pi = 0.51 GeV. We investigate the volume dependence of the energy shift of the ground state and the state of free nucleons to distinguish a bound state from attractive scattering states. From the investigation we conclude that ^4He, ^3He, deuteron and dineutron are bound at m_pi = 0.51 GeV. We compare their binding energies with those in our quenched studies and also with some recent investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Helium nuclei, deuteron and dineutron in 2+1 flavor lattice QCD

We calculate the binding energies for multi-nucleon bound states with the nuclear mass number less than or equal to 4 in 2+1 flavor QCD at the lattice spacing of a = 0.09 fm employing a relatively heavy quark mass corresponding to m_pi = 0.51 GeV. To distinguish a bound state from attractive scattering states, we investigate the volume dependence of the energy shift between the ground state and the state of free nucleons by changing the spatial extent of the lattice from 2.9 fm to 5.8 fm. We conclude that ^4He, ^3He, deuteron and dineutron are bound at m_pi = 0.51 GeV. We compare their binding energies with those in our quenched studies and also with several previous investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Electromagnetic form factor of pion from N_f=2+1 dynamical flavor QCD

We present a calculation of the electromagnetic form factor of the pion in $N_f=2+1$ flavor lattice QCD. Calculations are made on the PACS-CS gauge field configurations generated using Iwasaki gauge action and Wilson-clover quark action on a $32^3\times64$ lattice volume with the lattice spacing estimated as $a=0.0907(13)$ fm at the physical point. Measurements of the form factor are made using the technique of partially twisted boundary condition to reach small momentum transfer as well as periodic boundary condition with integer momenta. Additional improvements including random wall source techniques and a judicious choice of momenta carried by the incoming and outgoing quarks are employed for error reduction. Analyzing the form factor data for the pion mass at $M_π\approx 411$ MeV and 296 MeV, we find that the NNLO SU(2) chiral perturbation theory fit yields $<r^2>=0.441 \pm 0.046 {\rm fm}^2$ for the pion charge radius at the physical pion mass. Albeit the error is quite large, this is consistent with the experimental value of $0.452\pm 0.011 {\rm fm}^2$. Below $M_π\approx 300$ MeV, we find that statistical fluctuations in the pion two- and three-point functions become too large to extract statistically meaningful averages on a $32^3$ spatial volume. We carry out a sample calculation on a $64^4$ lattice with the quark masses close to the physical point, which suggests that form factor calculations at the physical point become feasible by enlarging lattice sizes to $M_πL\approx 4$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Domain Decomposition method on GPU cluster

Pallalel GPGPU computing for lattice QCD simulations has a bottleneck on the GPU to GPU data communication due to the lack of the direct data exchanging facility. In this work we investigate the performance of quark solver using the restricted additive Schwarz (RAS) preconditioner on a low cost GPU cluster. We expect that the RAS preconditioner with appropriate domaindecomposition and task distribution reduces the communication bottleneck. The GPU cluster we constructed is composed of four PC boxes, two GPU cards are attached to each box, and we have eight GPU cards in total. The compute nodes are connected with rather slow but low cost Gigabit-Ethernet. We include the RAS preconditioner in the single-precision part of the mixedprecision nested-BiCGStab algorithm and the single-precision task is distributed to the multiple GPUs. The benchmarking is done with the O(a)-improved Wilson quark on a randomly generated gauge configuration with the size of $32^4$. We observe a factor two improvment on the solver performance with the RAS precoditioner compared to that without the preconditioner and find that the improvment mainly comes from the reduction of the communication bottleneck as we expected.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Recent algorithm and machine developments for lattice QCD

I review recent machine trends and algorithmic developments for dynamical lattice QCD simulations with the HMC algorithm for Wilson-type fermions. The topics include the trend toward multi-core processors and general purpose GPU (GPGPU) computing, and improvements on the quark determinant preconditioning, molecular dynamics integrator, and quark solvers. I also discuss the prospect on the use of these techniques on the forthcoming petaflops machines.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Application of heavy-quark effective theory to lattice QCD: II. Radiative corrections to heavy-light currents

We apply heavy-quark effective theory to separate long- and short-distance effects of heavy quarks in lattice gauge theory. In this approach, the inverse heavy-quark mass and the lattice spacing are treated as short distances, and their effects are lumped into short-distance coefficients. We show how to use this formalism to match lattice gauge theory to continuum QCD, order by order in the heavy-quark expansion. In this paper, we focus on heavy-light currents. In particular, we obtain one-loop results for the matching factors of lattice currents, needed for heavy-quark phenomenology, such as the calculation of heavy-light decay constants, and heavy-to-light transition form factors. Results for the Brodsky-Lepage-Mackenzie scale $q^*$ are also given.

0 citations0 reviewsNo code linkedOpen access