Source author record

S. Aoki

S. Aoki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

hep-lat hep-ex hep-ph physics.ins-det nucl-th hep-th physics.acc-ph

Catalog footprint

What is connected

229works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Measurement of muon neutrino induced charged current interactions without charged pions in the final state using a new T2K off-axis near detector WAGASCI-BabyMIND

We report a flux-integrated cross section measurement of muon neutrino interactions on water and hydrocarbon via charged current reactions without charged pions in the final state with the WAGASCI-BabyMIND detector which was installed in the T2K near detector hall in 2018. The detector is located 1.5$^\circ$ off-axis and is exposed to a more energetic neutrino flux than ND280, another T2K near detector, which is located at a different off-axis position. The total flux-integrated cross section is measured to be $1.26 \pm 0.18\,(stat.+syst.) \times 10^{-39} $ $\mathrm{cm^{2}/nucleon}$ on CH and $1.44 \pm 0.21\,(stat.+syst.) \times 10^{-39} $ $\mathrm{cm^{2}/nucleon}$ on H$_{2}$O. These results are compared to model predictions provided by the NEUT v5.3.2 and GENIE v2.8.0 MC generators and the measurements are compatible with these models. Differential cross sections in muon momentum and cosine of the muon scattering angle are also reported. This is the first such measurement reported with the WAGASCI-BabyMIND detector and utilizes the 2020 and 2021 datasets.

preprint2025arXiv

Results from the T2K experiment on neutrino mixing including a new far detector $μ$-like sample

We have made improved measurements of three-flavor neutrino mixing with 19.7(16.3)$\times 10^{20}$ protons on target in (anti-)neutrino-enhanced beam modes. A new sample of muon-neutrino events with tagged pions has been added at the far detector, as well as new proton and photon-tagged samples at the near detector. Significant improvements have been made to the flux and neutrino interaction modeling. T2K data continue to prefer the normal mass ordering and upper octant of $\sin^2θ_{23}$ with a near-maximal value of the charge-parity violating phase with best-fit values in the normal ordering of $δ_{\scriptscriptstyle\mathrm{CP}}=-2.18\substack{+1.22 \\ -0.47}$, $\sin^2θ_{23}=0.559\substack{+0.018 \\ -0.078}$ and $Δm^2_{32}=(+2.506\substack{+0.039 \\ -0.052})\times 10^{-3}$ eV$^{2}$.

preprint2023arXiv

Updated constraints on sterile neutrino mixing in the OPERA experiment using a new $ν_e$ identification method

This paper describes a new $ν_e$ identification method specifically designed to improve the low-energy ($< 30\,\mathrm{GeV}$) $ν_e$ identification efficiency attained by enlarging the emulsion film scanning volume with the next generation emulsion readout system. A relative increase of 25-70% in the $ν_e$ low-energy region is expected, leading to improvements in the OPERA sensitivity to neutrino oscillations in the framework of the 3 + 1 model. The method is applied to a subset of data where the detection efficiency increase is expected to be more relevant, and one additional $ν_e$ candidate is found. The analysis combined with the $ν_τ$ appearance results improves the upper limit on $\sin^2 2θ_{μe}$ to 0.016 at 90% C.L. in the MiniBooNE allowed region $Δm^2_{41} \sim 0.3\,\mathrm{eV}^2$.

preprint2022arXiv

Measurements of protons and charged pions emitted from $ν_μ$ charged-current interactions on iron at a mean neutrino energy of 1.49$\,$GeV using a nuclear emulsion detector

This study conducted an analysis of muons, protons, and charged pions emitted from $ν_μ$ charged-current interactions on iron using a nuclear emulsion detector. The emulsion detector with a 65$\,$kg iron target was exposed to a neutrino beam corresponding to 4.0$\times$10$^{19}$ protons on target with a mean neutrino energy of 1.49$\,$GeV. The measurements were performed at a momentum threshold of 200 (50)$\,$MeV/$c$ for protons (pions), which are the lowest momentum thresholds attempted up to now. The measured quantities are the multiplicities, emission angles, and momenta of the muons, protons, and charged pions. In addition to these inclusive measurements, exclusive measurements such as the muon-proton emission-angle correlations of specific channels and the opening angle between the protons of CC0$π$2$p$ events were performed. The data were compared to Monte Carlo (MC) predictions and some significant differences were observed. The results of the study demonstrate the capability of detailed measurements of neutrino-nucleus interactions using a nuclear emulsion detector to improve neutrino interaction models.

preprint2022arXiv

Milliarcsecond astrometry for the Galilean moons using stellar occultations

A stellar occultation occurs when a Solar System object passes in front of a star for an observer. This technique allows the determination of sizes and shapes of the occulting body with kilometer precision. Also, this technique constrains the occulting body's positions, albedos, densities, etc. In the context of the Galilean moons, these events can provide their best ground-based astrometry, with uncertainties in the order of 1 mas ($\sim$ 3 km at Jupiter's distance during opposition). We organized campaigns and successfully observed a stellar occultation by Io (JI) in 2021, one by Ganymede (JIII) in 2020, and one by Europa (JII) in 2019, with stations in North and South America. Also, we re-analyzed two previously published events, one by Europa in 2016 and another by Ganymede in 2017. Then, we fit the known 3D shape of the occulting satellite and determine its center of figure. That resulted in astrometric positions with uncertainties in the milliarcsecond level. The positions obtained from these stellar occultations can be used together with dynamical models to ensure highly accurate orbits of the Galilean moons. These orbits can help plan future space probes aiming at the Jovian system, such as JUICE by ESA and Europa Clipper by NASA, and allow more efficient planning of flyby maneuvers.

preprint2022arXiv

Role of axial U(1) anomaly in chiral susceptibility of QCD at high temperature

The chiral susceptibility, or the first derivative of the chiral condensate with respect to the quark mass, is often used as a probe for the QCD phase transition since the chiral condensate is an order parameter of $SU(2)_L \times SU(2)_R$ symmetry breaking. However, the chiral condensate also breaks the axial $U(1)$ symmetry, which is usually not paid attention to as it is already broken by anomaly and apparently gives little impact on the transition. We investigate the susceptibilities in the scalar and pseudoscalar channels in order to quantify how much the axial $U(1)$ breaking contributes to the chiral phase transition. Employing a chirally symmetric lattice Dirac operator, and its eigenmode decomposition, we separate the axial $U(1)$ breaking effects from others. Our result in two-flavor QCD indicates that both of the connected and disconnected chiral susceptibilities are dominated by the axial $U(1)$ breaking at temperatures $T\gtrsim 190$ MeV after the quadratically divergent constant is subtracted.

preprint2022arXiv

Scintillator ageing of the T2K near detectors from 2010 to 2021

The T2K experiment widely uses plastic scintillator as a target for neutrino interactions and an active medium for the measurement of charged particles produced in neutrino interactions at its near detector complex. Over 10 years of operation the measured light yield recorded by the scintillator based subsystems has been observed to degrade by 0.9--2.2\% per year. Extrapolation of the degradation rate through to 2040 indicates the recorded light yield should remain above the lower threshold used by the current reconstruction algorithms for all subsystems. This will allow the near detectors to continue contributing to important physics measurements during the T2K-II and Hyper-Kamiokande eras. Additionally, work to disentangle the degradation of the plastic scintillator and wavelength shifting fibres shows that the reduction in light yield can be attributed to the ageing of the plastic scintillator.

preprint2022arXiv

What is chiral susceptibility probing?

In the early days of QCD, the axial $U(1)$ anomaly was considered as a trigger for the breaking of the $SU(2)_L\times SU(2)_R$ symmetry through topological excitations of gluon fields. However, it has been a challenge for lattice QCD to quantify the effect. In this work, we simulate QCD at high temperatures with chiral fermions. The exact chiral symmetry enables us to separate the contribution from the axial $U(1)$ breaking from others among the susceptibilities in the scalar and pseudoscalar channels. Our result in two-flavor QCD indicates that the chiral susceptibility, which is conventionally used as a probe for $SU(2)_L\times SU(2)_R$ breaking, is actually dominated by the axial $U(1)$ breaking at temperatures $T\ge 165$ MeV.

preprint2021arXiv

Constraint on the Matter-Antimatter Symmetry-Violating Phase in Neutrino Oscillations

The current laws of physics do not explain the observed imbalance of matter and antimatter in the universe. Sakharov proposed that an explanation would require the violation of CP symmetry between matter and antimatter. The only CP violation observed so far is in the weak interactions of quarks, and it is too small to explain the matter-antimatter imbalance of the universe. It has been shown that CP violation in the lepton sector could generate the matter-antimatter disparity through the process called leptogenesis. The quantum mixing of neutrinos, the neutral leptons in the Standard Model, provides a potential source of CP violation through a complex phase dCP, which may have consequences for theoretical models of leptogenesis. This CP violation can be measured in muon neutrino to electron neutrino oscillations and the corresponding antineutrino oscillations, which are experimentally accessible with accelerator-produced beams as established by the T2K experiment. Until now, the value of dCP has not been significantly constrained by neutrino oscillation experiments. Here the T2K collaboration reports a measurement that favors large enhancement of the neutrino oscillation probability, excluding values of dCP which result in a large enhancement of the observed anti-neutrino oscillation probability at three standard deviations (3 sigma). The 3 sigma confidence level interval for dCP, which is cyclic and repeats every 2pi, is [-3.41,-0.03] for the so-called normal mass ordering, and [-2.54,-0.32] for the inverted mass ordering. Our results show an indication of CP violation in the lepton sector. Herein we establish methods for sensitive searches for matter-antimatter asymmetry in neutrino oscillations using accelerator-produced neutrino beams. Future measurements with larger data samples will determine whether the leptonic CP violation is larger than the quark sector CP violation.

preprint2021arXiv

Sensitivity of the SHiP experiment to dark photons decaying to a pair of charged particles

Dark photons are hypothetical massive vector particles that could mix with ordinary photons. The simplest theoretical model is fully characterised by only two parameters: the mass of the dark photon m$_{γ^{\mathrm{D}}}$ and its mixing parameter with the photon, $\varepsilon$. The sensitivity of the SHiP detector is reviewed for dark photons in the mass range between 0.002 and 10 GeV. Different production mechanisms are simulated, with the dark photons decaying to pairs of visible fermions, including both leptons and quarks. Exclusion contours are presented and compared with those of past experiments. The SHiP detector is expected to have a unique sensitivity for m$_{γ^{\mathrm{D}}}$ ranging between 0.8 and 3.3$^{+0.2}_{-0.5}$ GeV, and $\varepsilon^2$ ranging between $10^{-11}$ and $10^{-17}$.

preprint2020arXiv

First combined measurement of the muon neutrino and antineutrino charged-current cross section without pions in the final state at T2K

This paper presents the first combined measurement of the double-differential muon neutrino and antineutrino charged-current cross sections with no pions in the final state on hydrocarbon at the off-axis near detector of the T2K experiment. The data analyzed in this work comprise 5.8$\times$10$^{20}$ and 6.3$\times$10$^{20}$ protons on target in neutrino and antineutrino mode respectively, at a beam energy peak of 0.6 GeV. Using the two measured cross sections, the sum, difference and asymmetry were calculated with the aim of better understanding the nuclear effects involved in such interactions. The extracted measurements have been compared with the prediction from different Monte Carlo generators and theoretical models showing that the difference between the two cross sections have interesting sensitivity to nuclear effects.

preprint2020arXiv

First observation of a tau neutrino charged current interaction with charm production in the Opera experiment

An event topology with two secondary vertices compatible with the decay of short-lived particles was found in the analysis of neutrino interactions in the Opera target. The observed topology is compatible with tau neutrino charged current (CC) interactions with charm production and neutrino neutral current (NC) interactions with $c\overline{c}$ pair production. However, other processes can mimic this topology. A dedicated analysis was implemented to identify the underlying process. A Monte Carlo simulation was developed and complementary procedures were introduced in the kinematic reconstruction. A multivariate analysis technique was used to achieve an optimal separation of signal from background. Most likely, this event is a $ν_τ$ CC interaction with charm production, the tau and charm particle decaying into 1 prong and 2 prongs, respectively. The significance of this observation is evaluated.

preprint2020arXiv

FLAG Review 2019

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to π$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_π$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ as well as those for $D$- and $B$-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $α_s$. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.

preprint2020arXiv

Search for heavy neutrinos with the T2K near detector ND280

This paper reports on the search for heavy neutrinos with masses in the range $140 < M_N < 493$ MeV/c$^2$ using the off-axis near detector ND280 of the T2K experiment. These particles can be produced from kaon decays in the standard neutrino beam and then subsequently decay in ND280. The decay modes under consideration are $N \to \ell^{\pm}_α π^{\mp}$ and $N \to \ell^+_α \ell^-_β ν(\barν)$ ($α,β=e,μ$). A search for such events has been made using the Time Projection Chambers of ND280, where the background has been reduced to less than two events in the current dataset in all channels. No excess has been observed in the signal region. A combined Bayesian statistical approach has been applied to extract upper limits on the mixing elements of heavy neutrinos to electron-, muon- and tau- flavoured currents ($U_e^2$, $U_μ^2$, $U_τ^2$) as a function of the heavy neutrino mass, e.g. $U_e^2 < 10^{-9}$ at $90\%$ C.L. for a mass of $390$ MeV/c$^2$. These constraints are competitive with previous experiments.

preprint2020arXiv

Simultaneous measurement of the muon neutrino charged-current cross section on oxygen and carbon without pions in the final state at T2K

This paper reports the first simultaneous measurement of the double differential muon neutrino charged-current cross section on oxygen and carbon without pions in the final state as a function of the outgoing muon kinematics, made at the ND280 off-axis near detector of the T2K experiment. The ratio of the oxygen and carbon cross sections is also provided to help validate various models' ability to extrapolate between carbon and oxygen nuclear targets, as is required in T2K oscillation analyses. The data are taken using a neutrino beam with an energy spectrum peaked at 0.6 GeV. The extracted measurement is compared with the prediction from different Monte Carlo neutrino-nucleus interaction event generators, showing particular model separation for very forward-going muons. Overall, of the models tested, the result is best described using Local Fermi Gas descriptions of the nuclear ground state with RPA suppression.

preprint2020arXiv

SND@LHC

We propose to build and operate a detector that, for the first time, will measure the process $pp\toνX$ at the LHC and search for feebly interacting particles (FIPs) in an unexplored domain. The TI18 tunnel has been identified as a suitable site to perform these measurements due to very low machine-induced background. The detector will be off-axis with respect to the ATLAS interaction point (IP1) and, given the pseudo-rapidity range accessible, the corresponding neutrinos will mostly come from charm decays: the proposed experiment will thus make the first test of the heavy flavour production in a pseudo-rapidity range that is not accessible by the current LHC detectors. In order to efficiently reconstruct neutrino interactions and identify their flavour, the detector will combine in the target region nuclear emulsion technology with scintillating fibre tracking layers and it will adopt a muon identification system based on scintillating bars that will also play the role of a hadronic calorimeter. The time of flight measurement will be achieved thanks to a dedicated timing detector. The detector will be a small-scale prototype of the scattering and neutrino detector (SND) of the SHiP experiment: the operation of this detector will provide an important test of the neutrino reconstruction in a high occupancy environment.

preprint2020arXiv

T2K measurements of muon neutrino and antineutrino disappearance using $3.13\times 10^{21}$ protons on target

We report measurements by the T2K experiment of the parameters $θ_{23}$ and $Δm^2_{32}$ which govern the disappearance of muon neutrinos and antineutrinos in the three-flavor PMNS neutrino oscillation model at T2K's neutrino energy and propagation distance. Utilizing the ability of the experiment to run with either a mainly neutrino or a mainly antineutrino beam, muon-like events from each beam mode are used to measure these parameters separately for neutrino and antineutrino oscillations. Data taken from $1.49 \times 10^{21}$ protons on target (POT) in neutrino mode and $1.64 \times 10^{21}$ POT in antineutrino mode are used. The best-fit values obtained by T2K were $\sin^2\left(θ_{23}\right)=0.51^{+0.06}_{-0.07} \left(0.43^{+0.21}_{-0.05}\right)$ and $Δm^2_{32}=2.47^{+0.08}_{-0.09} \left(2.50^{+0.18}_{-0.13}\right)$\evmass for neutrinos (antineutrinos). No significant differences between the values of the parameters describing the disappearance of muon neutrinos and antineutrinos were observed. An analysis using an effective two-flavor neutrino oscillation model where the sine of the mixing angle is allowed to take non-physical values larger than 1 is also performed to check the consistency of our data with the three-flavor model. Our data were found to be consistent with a physical value for the mixing angle.

preprint2019arXiv

First Measurement of the Charged Current $\overlineν_μ$ Double Differential Cross Section on a Water Target without Pions in the final state

This paper reports the first differential measurement of the charged-current $\overlineν_μ$ interaction cross section on water with no pions in the final state. The unfolded flux-averaged measurement using the T2K off-axis near detector is given in double differential bins of $μ^+$ momentum and angle. The integrated cross section in a restricted phase space is $σ=\left(1.11\pm0.18\right)\times10^{-38}$ cm$^{2}$ per water molecule. Comparisons with several nuclear models are also presented.

preprint2019arXiv

Measurement of neutrino and antineutrino neutral-current quasielastic-like interactions on oxygen by detecting nuclear de-excitation $γ$-rays

Neutrino- and antineutrino-oxygen neutral-current quasielastic-like interactions are measured at Super-Kamiokande using nuclear de-excitation $γ$-rays to identify signal-like interactions in data from a $14.94 \ (16.35)\times 10^{20}$ protons-on-target exposure of the T2K neutrino (antineutrino) beam. The measured flux-averaged cross sections on oxygen nuclei are $\langle σ_{ν{\rm -NCQE}} \rangle = 1.70 \pm 0.17 ({\rm stat.}) ^{+ {\rm 0.51}}_{- {\rm 0.38}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.82 GeV and $\langle σ_{\barν {\rm -NCQE}} \rangle = 0.98 \pm 0.16 ({\rm stat.}) ^{+ {\rm 0.26}}_{- {\rm 0.19}} ({\rm syst.}) \times 10^{-38} \ {\rm cm^2/oxygen}$ with a flux-averaged energy of 0.68 GeV, for neutrinos and antineutrinos, respectively. These results are the most precise to date, and the antineutrino result is the first cross section measurement of this channel. They are compared with various theoretical predictions. The impact on evaluation of backgrounds to searches for supernova relic neutrinos at present and future water Cherenkov detectors is also discussed.

preprint2019arXiv

Measurement of the muon neutrino charged-current single $π^+$ production on hydrocarbon using the T2K off-axis near detector ND280

We report the measurements of single and double differential cross section of muon neutrino charged-current interactions on carbon with a single positively charged pion in the final state at the T2K off-axis near detector using $5.56\times10^{20}$ protons on target. The analysis uses data control samples for the background subtraction and the cross section signal, defined as a single negatively charged muon and a single positively charged pion exiting from the target nucleus, is extracted using an unfolding method. The model dependent cross section, integrated over the T2K off-axis neutrino beam spectrum peaking at $0.6$~GeV, is measured to be $σ= (11.76 \pm 0.44 \text{(stat)} \pm 2.39 \text{(syst)}) \times 10^{-40} \text{cm}^2$~$\text{nucleon}^{-1}$. Various differential cross sections are measured, including the first measurement of the Adler angles for single charged pion production in neutrino interactions with heavy nuclei target.

preprint2019arXiv

Search for Electron Antineutrino Appearance in a Long-baseline Muon Antineutrino Beam

Electron antineutrino appearance is measured by the T2K experiment in an accelerator-produced antineutrino beam, using additional neutrino beam operation to constrain parameters of the PMNS mixing matrix. T2K observes 15 candidate electron antineutrino events with a background expectation of 9.3 events. Including information from the kinematic distribution of observed events, the hypothesis of no electron antineutrino appearance is disfavored with a significance of 2.40σ and no discrepancy between data and PMNS predictions is found. A complementary analysis that introduces an additional free parameter which allows non-PMNS values of electron neutrino and antineutrino appearance also finds no discrepancy between data and PMNS predictions.

preprint2016arXiv

Chiral behavior of light meson form factors in 2+1 flavor QCD with exact chiral symmetry

We present a study of chiral behavior of light meson form factors in QCD with three flavors of overlap quarks. Gauge ensembles are generated at single lattice spacing 0.12 fm with pion masses down to 300 MeV. The pion and kaon electromagnetic form factors and the kaon semileptonic form factors are precisely calculated using the all-to-all quark propagator. We discuss their chiral behavior using the next-to-next-to-leading order chiral perturbation theory.

preprint2016arXiv

Determination of the muon charge sign with the dipolar spectrometers of the OPERA experiment

The OPERA long-baseline neutrino-oscillation experiment has observed the direct appearance of $ν_τ$ in the CNGS $ν_μ$ beam. Two large muon magnetic spectrometers are used to identify muons produced in the $τ$ leptonic decay and in $ν_μ^{CC}$ interactions by measuring their charge and momentum. Besides the kinematic analysis of the $τ$ decays, background resulting from the decay of charmed particles produced in $ν_μ^{CC}$ interactions is reduced by efficiently identifying the muon track. A new method for the charge sign determination has been applied, via a weighted angular matching of the straight track-segments reconstructed in the different parts of the dipole magnets. Results obtained for Monte Carlo and real data are presented. Comparison with a method where no matching is used shows a significant reduction of up to 40\% of the fraction of wrongly determined charges.

preprint2016arXiv

Measurement of coherent $π^{+}$ production in low energy neutrino-Carbon scattering

We report the first measurement of the flux-averaged cross section for charged current coherent $π^{+}$ production on carbon for neutrino energies less than 1.5 GeV to a restricted final state phase space region in the T2K near detector, ND280. Comparisons are made with predictions from the Rein-Sehgal coherent production model and the model by Alvarez-Ruso {\it et al.}, the latter representing the first implementation of an instance of the new class of microscopic coherent models in a neutrino interaction Monte Carlo event generator. This results contradicts the null results reported by K2K and SciBooNE in a similar neutrino energy region.

preprint2016arXiv

Measurement of double-differential muon neutrino charged-current interactions on C$_8$H$_8$ without pions in the final state using the T2K off-axis beam

We report the measurement of muon neutrino charged-current interactions on carbon without pions in the final state at the T2K beam energy using 5.734$\times10^{20}$ protons on target. For the first time the measurement is reported as a flux-integrated, double-differential cross-section in muon kinematic variables ($\cosθ_μ$, $p_μ$), without correcting for events where a pion is produced and then absorbed by final state interactions. Two analyses are performed with different selections, background evaluations and cross-section extraction methods to demonstrate the robustness of the results against biases due to model-dependent assumptions. The measurements compare favorably with recent models which include nucleon-nucleon correlations but, given the present precision, the measurement does not solve the degeneracy between different models. The data also agree with Monte Carlo simulations which use effective parameters that are tuned to external data to describe the nuclear effects. The total cross-section in the full phase space is $σ= (0.417 \pm 0.047 \text{(syst)} \pm 0.005 \text{(stat)})\times 10^{-38} \text{cm}^2$ $\text{nucleon}^{-1}$ and the cross-section integrated in the region of phase space with largest efficiency and best signal-over-background ratio ($\cosθ_μ>0.6$ and $p_μ> 200$ MeV) is $σ= (0.202 \pm 0.0359 \text{(syst)} \pm 0.0026 \text{(stat)}) \times 10^{-38} \text{cm}^2$ $\text{nucleon}^{-1}$.

preprint2016arXiv

Measurement of Muon Antineutrino Oscillations with an Accelerator-Produced Off-Axis Beam

T2K reports its first measurements of the parameters governing the disappearance of $\barν_μ$ in an off-axis beam due to flavor change induced by neutrino oscillations. The quasimonochromatic $\barν_μ$ beam, produced with a peak energy of 0.6 GeV at J-PARC, is observed at the far detector Super-Kamiokande, 295 km away, where the $\barν_μ$ survival probability is expected to be minimal. Using a dataset corresponding to $4.01 \times 10^{20}$ protons on target, $34$ fully contained $μ$-like events were observed. The best-fit oscillation parameters are $\sin^2 (\barθ_{23}) = 0.45$ and $|Δ\bar{m}^2_{32}| = 2.51 \times 10^{-3}$ eV$^2$ with 68% confidence intervals of 0.38 - 0.64 and 2.26 - 2.80 $\times 10^{-3}$ eV$^2$ respectively. These results are in agreement with existing antineutrino parameter measurements and also with the $ν_μ$ disappearance parameters measured by T2K.

preprint2016arXiv

Proposal for an Extended Run of T2K to $20\times10^{21}$ POT

Recent measurements by the T2K neutrino oscillation experiment indicate that CP violation in neutrino mixing may be observed in the future by long-baseline neutrino oscillation experiments. We propose an extension to the currently approved T2K running from $7.8\times 10^{21}~\mbox{POT}$ to $20\times 10^{21}~\mbox{POT}$, aiming at initial observation of CP violation with 3$\,σ$ or higher significance for the case of maximum CP violation. The program also contains a measurement of mixing parameters, $θ_{23}$ and $Δm^2_{32}$, with a precision of 1.7$^\circ$ or better and 1%, respectively. With accelerator and beamline upgrades, as well as analysis improvements, this program would occur before the next generation of long-baseline neutrino oscillation experiments that are expected to start operation in 2026.

preprint2016arXiv

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

preprint2016arXiv

Sensitivity of the T2K accelerator-based neutrino experiment with an Extended run to $20\times10^{21}$ POT

Recent measurements at the T2K experiment indicate that CP violation in neutrino mixing may be observed in the future by long-baseline neutrino oscillation experiments. We explore the physics program of an extension to the currently approved T2K running of $7.8\times 10^{21}$ protons-on-target to $20\times 10^{21}$ protons-on-target,aiming at initial observation of CP violation with 3$\,σ$ or higher significance for the case of maximum CP violation. With accelerator and beam line upgrades, as well as analysis improvements, this program would occur before the next generation of long-baseline neutrino oscillation experiments that are expected to start operation in 2026.

preprint2015arXiv

$η^\prime$ meson mass from topological charge density correlator in QCD

The flavor-singlet component of the eta prime meson is related to the topological structure of the SU(3) gauge field through the chiral anomaly. We perform a 2+1-flavor lattice QCD calculation and demonstrate that the two-point function of a gluonically defined topological charge density after a short Yang-Mills gradient flow contains the propagation of the eta prime meson, by showing that its mass in the chiral and continuum limit is consistent with the experimental value. The gluonic correlator does not suffer from the contamination of the pion contribution, and the clean signal is obtained at significantly lower numerical cost compared to the conventional method with the quark bilinear operators.

preprint2015arXiv

A facility to Search for Hidden Particles (SHiP) at the CERN SPS

A new general purpose fixed target facility is proposed at the CERN SPS accelerator which is aimed at exploring the domain of hidden particles and make measurements with tau neutrinos. Hidden particles are predicted by a large number of models beyond the Standard Model. The high intensity of the SPS 400~GeV beam allows probing a wide variety of models containing light long-lived exotic particles with masses below ${\cal O}$(10)~GeV/c$^2$, including very weakly interacting low-energy SUSY states. The experimental programme of the proposed facility is capable of being extended in the future, e.g. to include direct searches for Dark Matter and Lepton Flavour Violation.

preprint2015arXiv

Discovery of tau neutrino appearance in the CNGS neutrino beam with the OPERA experiment

The OPERA experiment was designed to search for $ν_μ \rightarrow ν_τ$ oscillations in appearance mode, i.e. by detecting the $τ$-leptons produced in charged current $ν_τ$ interactions. The experiment took data from 2008 to 2012 in the CERN Neutrinos to Gran Sasso beam. The observation of $ν_μ \rightarrow ν_τ$ appearance, achieved with four candidate events in a sub-sample of the data, was previously reported. In this paper, a fifth $ν_τ$ candidate event, found in an enlarged data sample, is described. Together with a further reduction of the expected background, the candidate events detected so far allow assessing the discovery of $ν_μ\rightarrow ν_τ$ oscillations in appearance mode with a significance larger than 5 $σ$.

preprint2015arXiv

Light meson electromagnetic form factors from three-flavor lattice QCD with exact chiral symmetry

We study the chiral behavior of the electromagnetic (EM) form factors of pion and kaon in three-flavor lattice QCD. In order to make a direct comparison of the lattice data with chiral perturbation theory (ChPT), we employ the overlap quark action that has exact chiral symmetry. Gauge ensembles are generated at a lattice spacing of 0.11 fm with four pion masses ranging between M_pi \simeq 290 MeV and 540 MeV and with a strange quark mass m_s close to its physical value. We utilize the all-to-all quark propagator technique to calculate the EM form factors with high precision. Their dependence on m_s and on the momentum transfer is studied by using the reweighting technique and the twisted boundary conditions for the quark fields, respectively. A detailed comparison with SU(2) and SU(3) ChPT reveals that the next-to-next-to-leading order terms in the chiral expansion are important to describe the chiral behavior of the form factors in the pion mass range studied in this work. We estimate the relevant low-energy constants and the charge radii, and find reasonable agreement with phenomenological and experimental results.

preprint2015arXiv

Measurement of the $ν_μ$ CCQE cross section on carbon with the ND280 detector at T2K

The Charged-Current Quasi-Elastic (CCQE) interaction, $ν_{l} + n \rightarrow l^{-} + p$, is the dominant CC process at $E_ν\sim 1$ GeV and contributes to the signal in accelerator-based long-baseline neutrino oscillation experiments operating at intermediate neutrino energies. This paper reports a measurement by the T2K experiment of the $ν_μ$ CCQE cross section on a carbon target with the off-axis detector based on the observed distribution of muon momentum ($p_μ$) and angle with respect to the incident neutrino beam ($θ_μ$). The flux-integrated CCQE cross section was measured to be $(0.83 \pm 0.12) \times 10^{-38}\textrm{ cm}^{2}$ in good agreement with NEUT MC value of ${0.88 \times 10^{-38}} \textrm{ cm}^{2}$. The energy dependence of the CCQE cross section is also reported. The axial mass, $M_A^{QE}$, of the dipole axial form factor was extracted assuming the Smith-Moniz CCQE model with a relativistic Fermi gas nuclear model. Using the absolute (shape-only) $p_μcosθ_μ$ distribution, the effective $M_A^{QE}$ parameter was measured to be ${1.26^{+0.21}_{-0.18} \textrm{ GeV}/c^{2}}$ (${1.43^{+0.28}_{-0.22} \textrm{ GeV}/c^{2}}$).

preprint2015arXiv

Measurement of the $ν_μ$ charged current quasi-elastic cross-section on carbon with the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ charged current quasi-elastic cross-sections on carbon in the T2K on-axis neutrino beam. The measured charged current quasi-elastic cross-sections on carbon at mean neutrino energies of 1.94 GeV and 0.93 GeV are $(11.95\pm 0.19(stat.)_{-1.47}^{+1.82} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$ and $(10.64\pm 0.37(stat.)_{-1.65}^{+2.03} (syst.))\times 10^{-39}\mathrm{cm}^2/\mathrm{neutron}$, respectively. These results agree well with the predictions of neutrino interaction models. In addition, we investigated the effects of the nuclear model and the multi-nucleon interaction.

preprint2015arXiv

Measurement of the Electron Neutrino Charged-current Interaction Rate on Water with the T2K ND280 pi-zero Detector

This paper presents a measurement of the charged current interaction rate of the electron neutrino beam component of the beam above $1.5$~GeV using the large fiducial mass of the T2K $π^0$ detector. The predominant poriton of the $ν_e$ flux ($\sim$85 %) at these energies comes from kaon decays. The measured ratio of the observed beam interaction rate to the predicted rate in the detector with water targets filled is 0.89 $\pm$ 0.08 (stat.) $\pm$ 0.11 (sys.), and with the water targets emptied is 0.90 $\pm$ 0.09 (stat.) $\pm$ 0.13 (sys.). The ratio obtained for the interactions on water only from an event subtraction method is 0.87 $\pm$ 0.33 (stat.) $\pm$ 0.21 (sys.). This is the first measurement of the interaction rate of electron neutrinos on water, which is particularly of interest to experiments with water Cherenkov detectors.

preprint2015arXiv

Measurement of the muon neutrino inclusive charged-current cross section in the energy range of 1-3 GeV with the T2K INGRID detector

We report a measurement of the $ν_μ$-nucleus inclusive charged current cross section (=$σ^{cc}$) on iron using data from exposed to the J-PARC neutrino beam. The detector consists of 14 modules in total, which are spread over a range of off-axis angles from 0$^\circ$ to 1.1$^\circ$. The variation in the neutrino energy spectrum as a function of the off-axis angle, combined with event topology information, is used to calculate this cross section as a function of neutrino energy. The cross section is measured to be $σ^{cc}(1.1\text{ GeV}) = 1.10 \pm 0.15$ $(10^{-38}\text{cm}^2/\text{nucleon})$, $σ^{cc}(2.0\text{ GeV}) = 2.07 \pm 0.27$ $(10^{-38}\text{cm}^2/\text{nucleon})$, and $σ^{cc}(3.3\text{ GeV}) = 2.29 \pm 0.45$ $(10^{-38}\text{cm}^2/\text{nucleon})$, at energies of 1.1, 2.0, and 3.3 GeV, respectively. These results are consistent with the cross section calculated by the neutrino interaction generators currently used by T2K. More importantly, the method described here opens up a new way to determine the energy dependence of neutrino-nucleus cross sections.

preprint2015arXiv

Measurements of neutrino oscillation in appearance and disappearance channels by the T2K experiment with 6.6E20 protons on target

We report on measurements of neutrino oscillation using data from the T2K long-baseline neutrino experiment collected between 2010 and 2013. In an analysis of muon neutrino disappearance alone, we find the following estimates and 68% confidence intervals for the two possible mass hierarchies: Normal Hierarchy: $\sin^2θ_{23}=0.514^{+0.055}_{-0.056}$ and $Δm^2_{32}=(2.51\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ Inverted Hierarchy: $\sin^2θ_{23}=0.511\pm0.055$ and $Δm^2_{13}=(2.48\pm0.10)\times 10^{-3}$ eV$^2$/c$^4$ The analysis accounts for multi-nucleon mechanisms in neutrino interactions which were found to introduce negligible bias. We describe our first analyses that combine measurements of muon neutrino disappearance and electron neutrino appearance to estimate four oscillation parameters and the mass hierarchy. Frequentist and Bayesian intervals are presented for combinations of these parameters, with and without including recent reactor measurements. At 90% confidence level and including reactor measurements, we exclude the region: $δ_{CP}=[0.15,0.83]π$ for normal hierarchy and $δ_{CP}=[-0.08,1.09]π$ for inverted hierarchy. The T2K and reactor data weakly favor the normal hierarchy with a Bayes Factor of 2.2. The most probable values and 68% 1D credible intervals for the other oscillation parameters, when reactor data are included, are: $\sin^2θ_{23}=0.528^{+0.055}_{-0.038}$ and $|Δm^2_{32}|=(2.51\pm0.11)\times 10^{-3}$ eV$^2$/c$^4$.

preprint2015arXiv

Neutrino Oscillation Physics Potential of the T2K Experiment

The observation of the recent electron neutrino appearance in a muon neutrino beam and the high-precision measurement of the mixing angle $θ_{13}$ have led to a re-evaluation of the physics potential of the T2K long-baseline neutrino oscillation experiment. Sensitivities are explored for CP violation in neutrinos, non-maximal $\sin^22θ_{23}$, the octant of $θ_{23}$, and the mass hierarchy, in addition to the measurements of $δ_{CP}$, $\sin^2θ_{23}$, and $Δm^2_{32}$, for various combinations of $ν$-mode and $\barν$-mode data-taking. With an exposure of $7.8\times10^{21}$~protons-on-target, T2K can achieve 1-$σ$ resolution of 0.050(0.054) on $\sin^2θ_{23}$ and $0.040(0.045)\times10^{-3}~\rm{eV}^2$ on $Δm^2_{32}$ for 100\%(50\%) neutrino beam mode running assuming $\sin^2θ_{23}=0.5$ and $Δm^2_{32} = 2.4\times10^{-3}$ eV$^2$. T2K will have sensitivity to the CP-violating phase $δ_{\rm{CP}}$ at 90\% C.L. or better over a significant range. For example, if $\sin^22θ_{23}$ is maximal (i.e $θ_{23}$=$45^\circ$) the range is $-115^\circ<δ_{\rm{CP}}<-60^\circ$ for normal hierarchy and $+50^\circ<δ_{\rm{CP}}<+130^\circ$ for inverted hierarchy. When T2K data is combined with data from the NO$ν$A experiment, the region of oscillation parameter space where there is sensitivity to observe a non-zero $δ_{CP}$ is substantially increased compared to if each experiment is analyzed alone.

preprint2014arXiv

Computation of the electromagnetic pion form factor from lattice QCD in the epsilon regime

We calculate the electromagnetic pion form factor in lattice QCD with 2+1 flavors of the dynamical overlap quarks. Up and down quark masses are set below their physical values so that the system is in the so-called epsilon regime with the small size of our lattice ~ 1.8 fm. The finite volume corrections are generally expected to be ~ 100% in the epsilon regime. We, however, find a way to automatically cancel the dominant part of them. Inserting non-zero momenta and taking appropriate ratios of the two and three point functions, we can eliminate the contribution from the zero-momentum pion mode. Then the remaining finite volume effect is a small perturbation from the non-zero modes. Our lattice data agree with this theoretical prediction and the extracted pion charge radius is consistent with the experiment.

preprint2014arXiv

Evidence for $ν_μ\to ν_τ$ appearance in the CNGS neutrino beam with the OPERA experiment

The OPERA experiment is designed to search for $ν_μ \rightarrow ν_τ$ oscillations in appearance mode i.e. through the direct observation of the $τ$ lepton in $ν_τ$ charged current interactions. The experiment has taken data for five years, since 2008, with the CERN Neutrino to Gran Sasso beam. Previously, two $ν_τ$ candidates with a $τ$ decaying into hadrons were observed in a sub-sample of data of the 2008-2011 runs. Here we report the observation of a third $ν_τ$ candidate in the $τ^-\toμ^-$ decay channel coming from the analysis of a sub-sample of the 2012 run. Taking into account the estimated background, the absence of $ν_μ \rightarrow ν_τ$ oscillations is excluded at the 3.4 $σ$ level.

preprint2014arXiv

Histograms in heavy-quark QCD at finite temperature and density

We study the phase structure of lattice QCD with heavy quarks at finite temperature and density by a histogram method. We determine the location of the critical point at which the first-order deconfining transition in the heavy-quark limit turns into a crossover at intermediate quark masses through a change of the shape of the histogram under variation of coupling parameters. We estimate the effect of the complex phase factor which causes the sign problem at finite density, and show that, in heavy-quark QCD, the effect is small around the critical point. We determine the critical surface in 2+1 flavor QCD in the heavy-quark region at all values of the chemical potential mu including mu=infty.

preprint2014arXiv

Measurement of TeV atmospheric muon charge ratio with the full OPERA data

The OPERA detector, designed to search for $ν_μ \to ν_τ$ oscillations in the CNGS beam, is located in the underground Gran Sasso laboratory, a privileged location to study TeV-scale cosmic rays. For the analysis here presented, the detector was used to measure the atmospheric muon charge ratio in the TeV region. OPERA collected charge-separated cosmic ray data between 2008 and 2012. More than 3 million atmospheric muon events were detected and reconstructed, among which about 110000 multiple muon bundles. The charge ratio $R_μ \equiv N_{μ^+}/N_{μ^-}$ was measured separately for single and for multiple muon events. The analysis exploited the inversion of the magnet polarity which was performed on purpose during the 2012 Run. The combination of the two data sets with opposite magnet polarities allowed minimizing systematic uncertainties and reaching an accurate determination of the muon charge ratio. Data were fitted to obtain relevant parameters on the composition of primary cosmic rays and the associated kaon production in the forward fragmentation region. In the surface energy range 1-20 TeV investigated by OPERA, $R_μ$ is well described by a parametric model including only pion and kaon contributions to the muon flux, showing no significant contribution of the prompt component. The energy independence supports the validity of Feynman scaling in the fragmentation region up to $200$ TeV/nucleon primary energy.

preprint2014arXiv

Measurement of the inclusive $ν_μ$ charged current cross section on iron and hydrocarbon in the T2K on-axis neutrino beam

We report a measurement of the $ν_μ$ inclusive charged current cross sections on iron and hydrocarbon in the T2K on-axis neutrino beam. The measured inclusive charged current cross sections on iron and hydrocarbon averaged over the T2K on-axis flux with a mean neutrino energy of 1.51 GeV are $(1.444\pm0.002(stat.)_{-0.157}^{+0.189}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, and $(1.379\pm0.009(stat.)_{-0.147}^{+0.178}(syst.))\times 10^{-38}\mathrm{cm}^2/\mathrm{nucleon}$, respectively, and their cross section ratio is $1.047\pm0.007(stat.)\pm0.035(syst.)$. These results agree well with the predictions of the neutrino interaction model, and thus we checked the correct treatment of the nuclear effect for iron and hydrocarbon targets in the model within the measurement precisions.

preprint2014arXiv

Measurement of the Inclusive Electron Neutrino Charged Current Cross Section on Carbon with the T2K Near Detector

The T2K off-axis near detector, ND280, is used to make the first differential cross-section measurements of electron neutrino charged current interactions at energies ~1 GeV as a function of electron momentum, electron scattering angle and four-momentum transfer of the interaction. The total flux-averaged $ν_e$ charged current cross-section on carbon is measured to be $1.11\pm0.09~(stat)\pm0.18~(syst)\times10^{-38} cm^2/nucleon$. The differential and total cross-section measurements agree with the predictions of two leading neutrino interaction generators, NEUT and GENIE. The NEUT prediction is $1.23\times10^{-38} cm^2/nucleon$ and the GENIE prediction is $1.08\times10^{-38} cm^2/nucleon$. The total $ν_e$ charged current cross-section result is also in agreement with data from the Gargamelle experiment.

preprint2014arXiv

Measurement of the intrinsic electron neutrino component in the T2K neutrino beam with the ND280 detector

The T2K experiment has reported the first observation of the appearance of electron neutrinos in a muon neutrino beam. The main and irreducible background to the appearance signal comes from the presence in the neutrino beam of a small intrinsic component of electron neutrinos originating from muon and kaon decays. In T2K, this component is expected to represent 1.2% of the total neutrino flux. A measurement of this component using the near detector (ND280), located 280 m from the target, is presented. The charged current interactions of electron neutrinos are selected by combining the particle identification capabilities of both the time projection chambers and electromagnetic calorimeters of ND280. The measured ratio between the observed electron neutrino beam component and the prediction is 1.01+-0.10 providing a direct confirmation of the neutrino fluxes and neutrino cross section modeling used for T2K neutrino oscillation analyses. Electron neutrinos coming from muons and kaons decay are also separately measured, resulting in a ratio with respect to the prediction of 0.68+-0.30 and 1.10+-0.14, respectively.

preprint2014arXiv

Measurement of the muon beam direction and muon flux for the T2K neutrino experiment

The Tokai-to-Kamioka (T2K) neutrino experiment measures neutrino oscillations by using an almost pure muon neutrino beam produced at the J-PARC accelerator facility. The T2K muon monitor was installed to measure the direction and stability of the muon beam which is produced together with the muon neutrino beam. The systematic error in the muon beam direction measurement was estimated, using data and MC simulation, to be 0.28 mrad. During beam operation, the proton beam has been controlled using measurements from the muon monitor and the direction of the neutrino beam has been tuned to within 0.3 mrad with respect to the designed beam-axis. In order to understand the muon beam properties,measurement of the absolute muon yield at the muon monitor was conducted with an emulsion detector. The number of muon tracks was measured to be $(4.06\pm0.05)\times10^4$ cm$^{-2}$ normalized with $4\times10^{11}$ protons on target with 250 kA horn operation. The result is in agreement with the prediction which is corrected based on hadron production data.

preprint2014arXiv

Measurement of the neutrino-oxygen neutral-current interaction cross section by observing nuclear deexcitation $γ$ rays

We report the first measurement of the neutrino-oxygen neutral-current quasielastic (NCQE) cross section. It is obtained by observing nuclear deexcitation $γ$-rays which follow neutrino-oxygen interactions at the Super-Kamiokande water Cherenkov detector. We use T2K data corresponding to $3.01 \times 10^{20}$ protons on target. By selecting only events during the T2K beam window and with well-reconstructed vertices in the fiducial volume, the large background rate from natural radioactivity is dramatically reduced. We observe 43 events in the $4-30$ MeV reconstructed energy window, compared with an expectation of 51.0, which includes an estimated 16.2 background events. The background is primarily nonquasielastic neutral-current interactions and has only 1.2 events from natural radioactivity. The flux-averaged NCQE cross section we measure is $1.55 \times 10^{-38}$ cm$^2$ with a 68\% confidence interval of $(1.22, 2.20) \times 10^{-38}$ cm$^2$ at a median neutrino energy of 630 MeV, compared with the theoretical prediction of $2.01 \times 10^{-38}$ cm$^2$.

preprint2014arXiv

Observation of nu_tau appearance in the CNGS beam with the OPERA experiment

The OPERA experiment is searching for nu_mu -> nu_tau oscillations in appearance mode i.e. via the direct detection of tau leptons in nu_tau charged current interactions. The evidence of nu_mu -> nu_tau appearance has been previously reported with three nu_tau candidate events using a sub-sample of data from the 2008-2012 runs. We report here a fourth nu_tau candidate event, with the tau decaying into a hadron, found after adding the 2012 run events without any muon in the final state to the data sample. Given the number of analysed events and the low background, nu_mu -> nu_tau oscillations are established with a significance of 4.2sigma.

preprint2014arXiv

Precise Measurement of the Neutrino Mixing Parameter θ_{23} from Muon Neutrino Disappearance in an Off-axis Beam

New data from the T2K neutrino oscillation experiment produce the most precise measurement of the neutrino mixing parameter theta_{23}. Using an off-axis neutrino beam with a peak energy of 0.6 GeV and a data set corresponding to 6.57 x 10^{20} protons on target, T2K has fit the energy-dependent nu_mu oscillation probability to determine oscillation parameters. Marginalizing over the values of other oscillation parameters yields sin^2 (theta_{23}) = 0.514 +0.055/-0.056 (0.511 +- 0.055), assuming normal (inverted) mass hierarchy. The best-fit mass-squared splitting for normal hierarchy is Delta m^2_{32} = (2.51 +- 0.10) x 10^{-3} eV^2/c^4 (inverted hierarchy: Delta m^2_{13} = (2.48 +- 0.10) x 10^{-3} eV^2/c^4). Adding a model of multinucleon interactions that affect neutrino energy reconstruction is found to produce only small biases in neutrino oscillation parameter extraction at current levels of statistical uncertainty.

preprint2014arXiv

Procedure for short-lived particle detection in the OPERA experiment and its application to charm decays

The OPERA experiment, designed to perform the first observation of $ν_μ\rightarrow ν_τ$ oscillations in appearance mode through the detection of the $τ$ leptons produced in $ν_τ$ charged current interactions, has collected data from 2008 to 2012. In the present paper, the procedure developed to detect $τ$ particle decays, occurring over distances of the order of 1 mm from the neutrino interaction point, is described in detail. The results of its application to the search for charmed hadrons are then presented as a validation of the methods for $ν_τ$ appearance detection.

preprint2014arXiv

Search for short baseline $ν_e$ disappearance with the T2K near detector

The T2K experiment has performed a search for $ν_e$ disappearance due to sterile neutrinos using $5.9 \times 10^{20}$ protons on target for a baseline of $280 m$ in a neutrino beam peaked at about $500 MeV$. A sample of ν_e CC interactions in the off-axis near detector has been selected with a purity of 63\% and an efficiency of 26\%. The p-value for the null hypothesis is 0.085 and the excluded region at 95\% CL is approximately $sin^2 2 θ_{ee} > 0.3$ for $Δm^2_{eff} > 7 eV^2 / c^4$.

preprint2014arXiv

Spin-Orbit Force from Lattice QCD

We present a first attempt to determine nucleon-nucleon potentials in the parity-odd sector, which appear in 1P1, 3P0, 3P1, 3P2-3F2 channels, in Nf=2 lattice QCD simulations. These potentials are constructed from the Nambu-Bethe-Salpeter wave functions for J^P=0^-, 1^- and 2^-, which correspond to A1^-, T1^- and T2^- + E^- representation of the cubic group, respectively. We have found a large and attractive spin-orbit potential VLS(r) in the isospin-triplet channel, which is qualitatively consistent with the phenomenological determination from the experimental scattering phase shifts. The potentials obtained from lattice QCD are used to calculate the scattering phase shifts in 1P1, 3P0, 3P1 and 3P2-3F2 channels. The strong attractive spin-orbit force and a weak repulsive central force in spin-triplet P-wave channels lead to an attraction in the 3P2 channel, which is related to the P-wave neutron paring in neutron stars.

preprint2014arXiv

Strong coupling constant from vacuum polarization functions in three-flavor lattice QCD with dynamical overlap fermions

We determine the strong coupling constant $α_s$ from a lattice calculation of vacuum polarization functions (VPF) in three-flavor QCD with dynamical overlap fermions. Fitting lattice data of VPF to the continuum perturbative formula including the operator product expansion, we extract the QCD scale parameter $Λ_{\overline{MS}}^{(3)}$. At the $Z$ boson mass scale, we obtain $α_s^{(5)}(M_Z)=0.1181(3)(^{+14}_{-12})$, where the first error is statistical and the second is our estimate of various systematic uncertainties.

preprint2014arXiv

Topology density correlator on dynamical domain-wall ensembles with nearly frozen topological charge

Global topological charge decorrelates very slowly or even freezes in fine lattice simulations. On the other hand, its local fluctuations are expected to survive and lead to the correct physical results as long as the volume is large enough. We investigate this issue on recently generated configurations including dynamical domain-wall fermions at lattice spacings a = 0.08 fm and finer. We utilize the Yang-Mills gradient flow to define the topological charge density operator and calculate its long-distance correlation, through which we propose a new method for extracting the topological susceptibility in a sub-volume. This method takes care of the finite volume correction, which reduces the bias caused by the global topological charge. Our lattice data clearly show a shorter auto-correlation time than that of the naive definition using the whole lattice, and are less sensitive to the global topological history. Numerical results show a clear sea-quark mass dependence, which agrees well with the prediction of chiral perturbation theory.

preprint2013arXiv

Charmed baryons at the physical point in 2+1 flavor lattice QCD

We investigate the charmed baryon mass spectrum using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are tuned to their physical values, reweighted from those employed in the configuration generation. At the physical point, the inverse lattice spacing determined from the $Ω$ baryon mass gives $a^{-1}=2.194(10)$ GeV, and thus the spatial extent becomes $L = 32 a = 2.88(1)$ fm. Our results for the charmed baryon masses are consistent with experimental values, except for the mass of $Ξ_{cc}$, which has been measured by only one experimental group so far and has not been confirmed yet by others. In addition, we report values of other doubly and triply charmed baryon masses, which have never been measured experimentally.

preprint2013arXiv

Large-scale simulations with chiral symmetry

We carry out a comparative study among five-dimensional formulations of chirally symmetric fermions about the algorithmic performance, chiral symmetry violation and topological tunneling to find a computationally inexpensive formulation with good chiral symmetry. With our choice of the lattice action, we have launched large-scale simulations on fine lattices aiming at a precision study of light and heavy quark physics. We report on the comparative study, current status of the large-scale simulations, and preliminary results on the residual quark mass and auto-correlation.

preprint2013arXiv

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-axis Beam

The T2K collaboration reports a precision measurement of muon neutrino disappearance with an off-axis neutrino beam with a peak energy of 0.6 GeV. Near detector measurements are used to constrain the neutrino flux and cross section parameters. The Super-Kamiokande far detector, which is 295 km downstream of the neutrino production target, collected data corresponding to $3.01 \times 10^{20}$ protons on target. In the absence of neutrino oscillations, $205 \pm 17$ (syst.) events are expected to be detected and only 58 muon neutrino event candidates are observed. A fit to the neutrino rate and energy spectrum assuming three neutrino flavors, normal mass hierarchy and $θ_{23}\leq π/4$ yields a best-fit mixing angle $\sin^2(2θ_{23})=1.000$ and mass splitting $|Δm^2_{32}| =2.44 \times 10^{-3}$ eV$^2$/c$^4$. If $θ_{23}\geq π/4$ is assumed, the best-fit mixing angle changes to $\sin^2(2θ_{23})=0.999$ and the mass splitting remains unchanged.

preprint2013arXiv

Measurement of the Inclusive NuMu Charged Current Cross Section on Carbon in the Near Detector of the T2K Experiment

T2K has performed the first measurement of νμ inclusive charged current interactions on carbon at neutrino energies of ~1 GeV where the measurement is reported as a flux-averaged double differential cross section in muon momentum and angle. The flux is predicted by the beam Monte Carlo and external data, including the results from the NA61/SHINE experiment. The data used for this measurement were taken in 2010 and 2011, with a total of 10.8 x 10^{19} protons-on-target. The analysis is performed on 4485 inclusive charged current interaction candidates selected in the most upstream fine-grained scintillator detector of the near detector. The flux-averaged total cross section is <σ_CC>_ϕ=(6.91 +/- 0.13 (stat) +/- 0.84 (syst)) x10^{-39} cm^2/nucleon for a mean neutrino energy of 0.85 GeV.

preprint2013arXiv

New results on $ν_μ\to ν_τ$ appearance with the OPERA experiment in the CNGS beam

The OPERA neutrino experiment is designed to perform the first observation of neutrino oscillations in direct appearance mode in the $ν_μ\to ν_τ$ channel, via the detection of the $τ$-leptons created in charged current $ν_τ$ interactions. The detector, located in the underground Gran Sasso Laboratory, consists of an emulsion/lead target with an average mass of about 1.2 kt, complemented by electronic detectors. It is exposed to the CERN Neutrinos to Gran Sasso beam, with a baseline of 730 km and a mean energy of 17 GeV. The observation of the first $ν_τ$ candidate event and the analysis of the 2008-2009 neutrino sample have been reported in previous publications. This work describes substantial improvements in the analysis and in the evaluation of the detection efficiencies and backgrounds using new simulation tools. The analysis is extended to a sub-sample of 2010 and 2011 data, resulting from an electronic detector-based pre-selection, in which an additional $ν_τ$ candidate has been observed. The significance of the two events in terms of a $ν_μ\to ν_τ$ oscillation signal is of 2.40 $σ$.

preprint2013arXiv

Nucleon strange quark content from N_f=2+1 lattice QCD with exact chiral symmetry

We calculate the strange quark content of the nucleon in 2+1-flavor lattice QCD. Chirally symmetric overlap fermion formulation is used to avoid the contamination from up and down quark contents due to an operator mixing between strange and light scalar operators, \bar{s}s and \bar{u}u+\bar{d}d. At a lattice spacing a=0.112(1) fm, we perform calculations at four values of degenerate up and down quark masses, which cover a range of the pion mass M_pi \simeq 300-540 MeV. We employ two different methods: one is a direct method where we calculate the strange quark content by directly inserting the strange scalar operator. The other is an indirect method where the quark content is extracted from a derivative of the nucleon mass in terms of the strange quark mass. With these two methods we obtain consistent results with each other. Our best estimate f_{T_s}=0.009(15)(16) is in good agreement with our previous studies in two-flavor QCD.

preprint2013arXiv

Overlap/Domain-wall reweighting

We investigate the eigenvalues of nearly chiral lattice Dirac operators constructed with five-dimensional implementations. Allowing small violation of the Ginsparg-Wilson relation, the HMC simulation is made much faster while the eigenvalues are not significantly affected. We discuss the possibility of reweighting the gauge configurations generated with domain-wall fermions to those of exactly chiral lattice fermions.

preprint2013arXiv

Phase shifts in I=2 ππ-scattering from two lattice approaches

We present a lattice QCD study of the phase shift of $I{=}2$ $ππ$ scattering on the basis of two different approaches: the standard finite volume approach by Luscher and the recently introduced HAL QCD potential method. Quenched QCD simulations are performed on lattices with extents $N_s{=}16,24,32,48$ and $N_t{=}128$ as well as lattice spacing $a{\sim}0.115\,\mathrm{fm}$ and a pion mass of $m_π{\sim}940\,\mathrm{MeV}$. The phase shift and the scattering length are calculated in these two methods. In the potential method, the error is dominated by the systematic uncertainty associated with the violation of rotational symmetry due to finite lattice spacing. In Luscher's approach, such systematic uncertainty is difficult to be evaluated and thus is not included in this work. A systematic uncertainty attributed to the quenched approximation, however, is not evaluated in both methods. In case of the potential method, the phase shift can be calculated for arbitrary energies below the inelastic threshold. The energy dependence of the phase shift is also obtained from Luscher's method using different volumes and/or nonrest-frame extension of it. The results are found to agree well with the potential method.

preprint2013arXiv

Search for $ν_μ\rightarrow ν_e$ oscillations with the OPERA experiment in the CNGS beam

A first result of the search for \numu $\rightarrow$ \nue oscillations in the OPERA experiment, located at the Gran Sasso Underground Laboratory, is presented. The experiment looked for the appearance of \nue in the CNGS neutrino beam using the data collected in 2008 and 2009. Data are compatible with the non-oscillation hypothesis in the three-flavour mixing model. A further analysis of the same data constrains the non-standard oscillation parameters $θ_{new}$ and $Δm^2_{new}$ suggested by the LSND and MiniBooNE experiments. For large $Δm^{2}_{new}$ values ($>$0.1 eV$^{2}$), the OPERA 90% C.L. upper limit on sin$^{2}(2θ_{new})$ based on a Bayesian statistical method reaches the value $7.2 \times 10^{-3}$.

preprint2013arXiv

Short Distance Repulsion Among Baryons

We review recent investigations on the short distance behaviors of potentials among baryons, which are formulated through the Nambu-Bethe-Salpeter (NBS) wave function. After explaining the method to define the potentials, we analyze the short distance behavior of the NBS wave functions and the corresponding potentials by combining the operator product expansion and a renormalization group analysis in the perturbation theory of QCD. These analytic results are compared with numerical results obtained in lattice QCD simulations.

preprint2013arXiv

The T2K Neutrino Flux Prediction

The Tokai-to-Kamioka (T2K) experiment studies neutrino oscillations using an off-axis muon neutrino beam with a peak energy of about 0.6 GeV that originates at the J-PARC accelerator facility. Interactions of the neutrinos are observed at near detectors placed at 280 m from the production target and at the far detector -- Super-Kamiokande (SK) -- located 295 km away. The flux prediction is an essential part of the successful prediction of neutrino interaction rates at the T2K detectors and is an important input to T2K neutrino oscillation and cross section measurements. A FLUKA and GEANT3 based simulation models the physical processes involved in the neutrino production, from the interaction of primary beam protons in the T2K target, to the decay of hadrons and muons that produce neutrinos. The simulation uses proton beam monitor measurements as inputs. The modeling of hadronic interactions is re-weighted using thin target hadron production data, including recent charged pion and kaon measurements from the NA61/SHINE experiment. For the first T2K analyses the uncertainties on the flux prediction are evaluated to be below 15% near the flux peak. The uncertainty on the ratio of the flux predictions at the far and near detectors is less than 2% near the flux peak.

preprint2012arXiv

1+1+1 flavor QCD + QED simulation at the physical point

We present the results of 1+1+1 flavor QCD+QED simulation at the physical point, in which the dynamical quark effects in QED and the up-down quark mass difference are incorporated by the reweighting technique. The physical quark masses together with the lattice spacing are determined with $m_{π^+}$, $m_{K^+}$, $m_{K^0}$ and $m_{Ω^-}$ as physical inputs. Calculations are carried out using a set of 2+1 flavor QCD configurations near the physical point generated by the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $β=1.9$ on a $32^3\times 64$ lattice. We evaluate the values of the up, down and strange quark masses individually with non-perturbative QCD renormalization.

preprint2012arXiv

Application of fixed scale approach to static quark free energies in quenched and 2+1 flavor lattice QCD with improved Wilson quark action

Free energies between static quarks and Debye screening masses in the quark-gluon plasma are studied on the basis of Polyakov-line correlations in lattice simulations of 2+1 flavors QCD with the renormalization-group improved gluon action and the $O(a)$-improved Wilson quark action. We perform simulations at $m_{\rm PS}/m_{\rm V} = 0.63$ (0.74) for light (strange) flavors with lattice sizes of $32^3 \times N_t$ with $N_t=4$--12. We adopt the fixed-scale approach, where temperature can be varied without changing the spatial volume and renormalization factor. We find that, at short distance, the free energies of static quarks in color-singlet channel converge to the static-quark potential evaluated from the Wilson-loop at zero-temperature, in accordance with the expected insensitivity of short distance physics to the temperature. At long distance, the free energies of static quarks approach to twice the single-quark free energies, implying that the interaction between static quarks is fully screened. The screening properties can be well described by the screened Coulomb form with appropriate Casimir factor at high temperature. We also discuss a limitation of the fixed-scale approach at high temperature.

preprint2012arXiv

Chiral behavior of kaon semileptonic form factors in lattice QCD with exact chiral symmetry

We calculate the kaon semileptonic form factors in lattice QCD with three flavors of dynamical overlap quarks. Gauge ensembles are generated at pion masses as low as 290 MeV and at a strange quark mass near its physical value. We precisely calculate relevant meson correlators using the all-to-all quark propagator. Twisted boundary conditions and the reweighting technique are employed to vary the momentum transfer and the strange quark mass. We discuss the chiral behavior of the form factors by comparing with chiral perturbation theory and experiments.

preprint2012arXiv

Equation of state in 2+1 flavor QCD with improved Wilson quarks by the fixed scale approach

We study the equation of state in 2+1 flavor QCD with nonperturbatively improved Wilson quarks coupled with the RG-improved Iwasaki glue. We apply the $T$-integration method to nonperturbatively calculate the equation of state by the fixed-scale approach. With the fixed-scale approach, we can purely vary the temperature on a line of constant physics without changing the system size and renormalization constants. Unlike the conventional fixed-$N_t$ approach, it is easy to keep scaling violations small at low temperature in the fixed scale approach. We study 2+1 flavor QCD at light quark mass corresponding to $m_π/m_ρ\simeq 0.63$, while the strange quark mass is chosen around the physical point. Although the light quark masses are heavier than the physical values yet, our equation of state is roughly consistent with recent results with highly improved staggered quarks at large $N_t$.

preprint2012arXiv

Finite density QCD phase transition in the heavy quark region

We extend our previous study of the QCD phase structure in the heavy quark region to non-zero chemical potentials. To identify the critical point where the first order deconfining transition terminates, we study an effective potential defined by the probability distribution function of the plaquette and the Polyakov loop. The reweighting technique is shown to be powerful in evaluating the effective potential in a wide range of the plaquette and Polyakov loop expectation values. We adopt the cumulant expansion to overcome the sign problem in the calculation of complex phase of the quark determinant. We find that the method provides us with an intuitive and powerful way to study the phase structure. We estimate the location of the critical point at finite chemical potential in the heavy quark region.

preprint2012arXiv

First Muon-Neutrino Disappearance Study with an Off-Axis Beam

We report a measurement of muon-neutrino disappearance in the T2K experiment. The 295-km muon-neutrino beam from Tokai to Kamioka is the first implementation of the off-axis technique in a long-baseline neutrino oscillation experiment. With data corresponding to 1.43 10**20 protons on target, we observe 31 fully-contained single muon-like ring events in Super-Kamiokande, compared with an expectation of 104 +- 14 (syst) events without neutrino oscillations. The best-fit point for two-flavor nu_mu -> nu_tau oscillations is sin**2(2 theta_23) = 0.98 and |Δm**2_32| = 2.65 10**-3 eV**2. The boundary of the 90 % confidence region includes the points (sin**2(2 theta_23),|Δm**2_32|) = (1.0, 3.1 10**-3 eV**2), (0.84, 2.65 10**-3 eV**2) and (1.0, 2.2 10**-3 eV**2).

preprint2012arXiv

Measurement of the neutrino velocity with the OPERA detector in the CNGS beam

The OPERA neutrino experiment at the underground Gran Sasso Laboratory has measured the velocity of neutrinos from the CERN CNGS beam over a baseline of about 730 km. The measurement is based on data taken by OPERA in the years 2009, 2010 and 2011. Dedicated upgrades of the CNGS timing system and of the OPERA detector, as well as a high precision geodesy campaign for the measurement of the neutrino baseline, allowed reaching comparable systematic and statistical accuracies. An arrival time of CNGS muon neutrinos with respect to the one computed assuming the speed of light in vacuum of (6.5 +/- 7.4(stat.)((+8.3)(-8.0)sys.))ns was measured corresponding to a relative difference of the muon neutrino velocity with respect to the speed of light (v-c)/c =(2.7 +/-3.1(stat.)((+3.4)(-3.3)(sys.))x10^(-6). The above result, obtained by comparing the time distributions of neutrino interactions and of protons hitting the CNGS target in 10.5 microseconds long extractions, was confirmed by a test performed at the end of 2011 using a short bunch beam allowing to measure the neutrino time of flight at the single interaction level.

preprint2012arXiv

Phase structure of finite temperature QCD in the heavy quark region

We study the quark mass dependence of the finite temperature QCD phase transition in the heavy quark region using an effective potential defined through the probability distribution function of the average plaquette. Performing a simulation of SU(3) pure gauge theory, we first confirm that the distribution function has two peaks indicating that the phase transition is of first order in the heavy quark limit, while the first order transition turns into a crossover as the quark mass decreases from infinity, where the mass dependence of the distribution function is evaluated by the reweighting method combined with the hopping parameter expansion. We determine the endpoint of the first order transition region for N_f=1, 2, 3 and 2+1 cases. The quark mass dependence of the latent heat is also evaluated in the first order transition region.

preprint2012arXiv

Pion form factors in the epsilon regime

We consider how to extract the pion form factors in the epsilon regime. Using the correlators with non-zero momenta and taking appropriate ratios of them, we eliminate the dominant finite volume effect from the zero-momentum pion mode. Our preliminary lattice result for the pion charge radius is consistent with the experiment.

preprint2012arXiv

Probability distribution functions in the finite density lattice QCD

We study the phase structure of QCD at high temperature and density by lattice QCD simulations adopting a histogram method. We try to solve the problems which arise in the numerical study of the finite density QCD, focusing on the probability distribution function (histogram). As a first step, we investigate the quark mass dependence and the chemical potential dependence of the probability distribution function as a function of the Polyakov loop when all quark masses are sufficiently large, and study the properties of the distribution function. The effect from the complex phase of the quark determinant is estimated explicitly. The shape of the distribution function changes with the quark mass and the chemical potential. Through the shape of the distribution, the critical surface which separates the first order transition and crossover regions in the heavy quark region is determined for the 2+1-flavor case.

preprint2012arXiv

The T2K Side Muon Range Detector

The T2K experiment is a long baseline neutrino oscillation experiment aiming to observe the appearance of ν e in a νμ beam. The νμ beam is produced at the Japan Proton Accelerator Research Complex (J-PARC), observed with the 295 km distant Super- Kamiokande Detector and monitored by a suite of near detectors at 280m from the proton target. The near detectors include a magnetized off-axis detector (ND280) which measures the un-oscillated neutrino flux and neutrino cross sections. The present paper describes the outermost component of ND280 which is a side muon range detector (SMRD) composed of scintillation counters with embedded wavelength shifting fibers and Multi-Pixel Photon Counter read-out. The components, performance and response of the SMRD are presented.

preprint2012arXiv

Thermodynamics in 2+1 flavor QCD with improved Wilson quarks by the fixed scale approach

We study thermodynamic properties of 2+1 flavor QCD with improved Wilson quarks coupled with the RG improved Iwasaki glue, using the fixed scale approach. We present the results for the equation of state, renormalized Polyakov loop, and chiral condensate.

preprint2011arXiv

$ρ$ Meson Decay in 2+1 Flavor Lattice QCD

We perform a lattice QCD study of the $ρ$ meson decay from the $N_f=2+1$ full QCD configurations generated with a renormalization group improved gauge action and a non-perturbatively $O(a)$-improved Wilson fermion action. The resonance parameters, the effective $ρ\toππ$ coupling constant and the resonance mass, are estimated from the $P$-wave scattering phase shift for the isospin I=1 two-pion system. The finite size formulas are employed to calculate the phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, $m_π=410\,{\rm MeV}$ ($m_π/m_ρ=0.46$) and $m_π=300\,{\rm MeV}$ ($m_π/m_ρ=0.35$), on a $32^3\times 64$ ($La=2.9\,{\rm fm}$) lattice at the lattice spacing $a=0.091\,{\rm fm}$. We compare our results at these two quark masses with those given in the previous works using $N_f=2$ full QCD configurations and the experiment.

preprint2011arXiv

Charm quark system at the physical point of 2+1 flavor lattice QCD

We investigate the charm quark system using the relativistic heavy quark action on 2+1 flavor PACS-CS configurations previously generated on $32^3 \times 64$ lattice. The dynamical up-down and strange quark masses are set to the physical values by using the technique of reweighting to shift the quark hopping parameters from the values employed in the configuration generation. At the physical point, the lattice spacing equals $a^{-1}=2.194(10)$ GeV and the spatial extent $L=2.88(1)$ fm. The charm quark mass is determined by the spin-averaged mass of the 1S charmonium state, from which we obtain $m_{\rm charm}^{\msbar}(μ= m_{\rm charm}^{\msbar}) = 1.260(1)(6)(35)$ GeV, where the errors are due to our statistics, scale determination and renormalization factor. An additional systematic error from the heavy quark is of order $α_s^2 f(m_Q a)(a Λ_{QCD})$, which is estimated to be a percent level if the factor $f(m_Q a)$ analytic in $m_Q a$ is of order unity. Our results for the charmed and charmed-strange meson decay constants are $f_D=226(6)(1)(5)$ MeV, $f_{D_s}=257(2)(1)(5)$ MeV, again up to the heavy quark errors of order $α_s^2 f(m_Q a)(a Λ_{QCD})$. Combined with the CLEO values for the leptonic decay widths, these values yield $|V_{cd}| = 0.205(6)(1)(5)(9)$, $|V_{cs}| = 1.00(1)(1)(3)(3)$, where the last error is on account of the experimental uncertainty of the decay widths.

preprint2011arXiv

Charmonium spectral functions with the variational method in zero and finite temperature lattice QCD

We propose a method to evaluate spectral functions on the lattice based on a variational method. On a lattice with a finite spatial extent, spectral functions consist of discrete spectra only. Adopting a variational method, we calculate the locations and the heights of spectral functions at low-lying discrete spectra. We first test the method in the case of analytically solvable free Wilson quarks at zero and finite temperatures and confirm that the method well reproduces the analytic results for low-lying spectra. We find that we can systematically improve the results by increasing the number of trial states. We then apply the method to calculate the charmonium spectral functions for S and P-wave states at zero-temperature in quenched QCD and compare the results with those obtained using the conventional maximum entropy method (MEM). The results for the ground state are consistent with the location and the area of the first peak in spectral functions from the MEM, while the variational method leads to a mass which is closer to the experimental value for the first excited state. We also investigate the temperature dependence of the spectral functions for S-wave states below and above $T_c$. We obtain no clear evidences for dissociation of $J/ψ$ and $η_c$ up to 1.4$T_c$.

preprint2011arXiv

Chiral interpolation in a finite volume

A simulation of lattice QCD at (or even below) the physical pion mass is feasible on a small lattice size of \sim 2 fm. The results are, however, subject to large finite volume effects. In order to precisely understand the chiral behavior in a finite volume, we develop a new computational scheme to interpolate the conventional epsilon and p regimes within chiral perturbation theory. In this new scheme, we calculate the two-point function in the pseudoscalar channel, which is described by a set of Bessel functions in an infra-red finite way as in the epsilon regime, while chiral logarithmic effects are kept manifest as in the p regime. The new ChPT formula is compared to our 2+1- flavor lattice QCD data near the physical up and down quark mass, mud \sim 3 MeV on an L \sim 1.8 fm lattice. We extract the pion mass = 99(4) MeV, from which we attempt a chiral "interpolation" of the observables to the physical point.

preprint2011arXiv

Determination of the chiral condensate from QCD Dirac spectrum on the lattice

We calculate the chiral condensate of QCD with 2, 2+1 and 3 flavors of sea quarks. Lattice QCD simulations are performed employing dynamical overlap fermions with up and down quark masses covering a range between 3 and 100 MeV. On L ~ 1.8-1.9 fm lattices at a lattice spacing ~ 0.11 fm, we calculate the eigenvalue spectrum of the overlap-Dirac operator. By matching the lattice data with the analytical prediction from chiral perturbation theory at the next-to-leading order, the chiral condensate in the massless limit of up and down quarks is determined.

preprint2011arXiv

Histogram method in finite density QCD with phase quenched simulations

We propose a new approach to finite density QCD based on a histogram method with phase quenched simulations at finite chemical potential. Integrating numerically the derivatives of the logarithm of the quark determinant with respect to the chemical potential, we calculate the reweighting factor and the complex phase of the quark determinant. The complex phase is handled with a cumulant expansion to avoid the sign problem. We examine the applicability of this method.

preprint2011arXiv

Indication of Electron Neutrino Appearance from an Accelerator-produced Off-axis Muon Neutrino Beam

The T2K experiment observes indications of $ν_μ\rightarrow ν_e$ appearance in data accumulated with $1.43\times10^{20}$ protons on target. Six events pass all selection criteria at the far detector. In a three-flavor neutrino oscillation scenario with $|Δm_{23}^2|=2.4\times10^{-3}$ eV$^2$, $\sin^2 2θ_{23}=1$ and $\sin^2 2θ_{13}=0$, the expected number of such events is 1.5$\pm$0.3(syst.). Under this hypothesis, the probability to observe six or more candidate events is 7$\times10^{-3}$, equivalent to 2.5$σ$ significance. At 90% C.L., the data are consistent with 0.03(0.04)$<\sin^2 2θ_{13}<$ 0.28(0.34) for $δ_{\rm CP}=0$ and a normal (inverted) hierarchy.

preprint2011arXiv

Kaon semileptonic form factors in QCD with exact chiral symmetry

We report on our calculation of the kaon semileptonic form factors in Nf=2+1 lattice QCD. Chiral symmetry is exactly preserved by using the overlap quark action for a straightforward comparison with chiral perturbation theory (ChPT). We simulate three pion masses down to 290 MeV at a single lattice spacing of 0.11 fm and at a strange quark mass very close to its physical value. The form factors near zero momentum transfer are precisely calculated by using the all-to-all propagator and twisted boundary conditions. We compare the normalizations and slopes of the form factors with ChPT and experiments.

preprint2011arXiv

Measurement of charm production in neutrino charged-current interactions

The nuclear emulsion target of the CHORUS detector was exposed to the wide-band neutrino beam of the CERN SPS of 27 GeV average neutrino energy from 1994 to 1997. In total about 100000 charged-current neutrino interactions with at least one identified muon were located in the emulsion target and fully reconstructed, using newly developed automated scanning systems. Charmed particles were searched for by a program recognizing particle decays. The observation of the decay in nuclear emulsion makes it possible to select a sample with very low background and minimal kinematical bias. 2013 charged-current interactions with a charmed hadron candidate in the final state were selected and confirmed through visual inspection. The charm production rate induced by neutrinos relative to the charged-current cross-section is measured to be sigma(nu_mu N -> mu- C X)/sigma(CC) = (5.75 +-0.32 stat +-0.30 syst)%. The charm production cross-section as a function of the neutrino energy is also obtained. The results are in good agreement with previous measurements. The charm-quark hadronization produces the following charmed hadrons with relative fractions (in %): f_Dzero = 43.7+-4.5, f_Lambda_c^plus = 19.2+-4.2, f_Dplus = 25.3+-4.2, and f_D_splus = 11.8+-4.7.

preprint2011arXiv

Measurements of the T2K neutrino beam properties using the INGRID on-axis near detector

Precise measurement of neutrino beam direction and intensity was achieved based on a new concept with modularized neutrino detectors. INGRID (Interactive Neutrino GRID) is an on-axis near detector for the T2K long baseline neutrino oscillation experiment. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The module has a sandwich structure of iron target plates and scintillator trackers. INGRID directly monitors the muon neutrino beam profile center and intensity using the number of observed neutrino events in each module. The neutrino beam direction is measured with accuracy better than 0.4 mrad from the measured profile center. The normalized event rate is measured with 4% precision.

preprint2011arXiv

Nucleon strange quark content from two-flavor lattice QCD with exact chiral symmetry

Strange quark content of the nucleon is calculated in dynamical lattice QCD employing the overlap fermion formulation. For this quantity, exact chiral symmetry guaranteed by the Ginsparg-Wilson relation is crucial to avoid large contamination due to a possible operator mixing with $\bar{u}u+\bar{d}d$. Gauge configurations are generated with two dynamical flavors on a 16^3 x 32 lattice at a lattice spacing a \simeq 0.12fm. We directly calculate the relevant three-point function on the lattice including a disconnected strange quark loop utilizing the techniques of all-to-all quark propagator and low-mode averaging. Our result f_{T_s} = 0.032(8)(22), is in good agreement with our previous indirect estimate using the Feynman-Hellmann theorem.

preprint2011arXiv

Scaling behavior of chiral phase transition in two-flavor QCD with improved Wilson quarks at finite density

We study scaling behavior of a chiral order parameter performing a simulation of two-flavor QCD with improved Wilson quarks. It has been shown that the scaling behavior of the chiral order parameter defined by a Ward-Takahashi identity agrees with the scaling function of the three-dimensional O(4) spin model at zero chemical potential. We extend the scaling study to finite density QCD. Calculating derivatives of the chiral order parameter with respect to the chemical potential in two-flavor QCD, the scaling property of chiral phase transition is discussed in the low density region. We moreover calculate the curvature of the phase boundary of the chirl phase transition in the temperature and chemical potential plane assuming the O(4) scaling relation.

preprint2011arXiv

Two-photon decay of $π^0$ from two-flavor lattice QCD

We study the correction to the radiative $π^0$ decay width due to finite light quark mass. Using lattice QCD with the overlap fermion formulation, we calculate the three-point function of the form $<PV_μV_ν>$ in the (Euclidean) momentum space, which corresponds to the $π^0\rightarrow γ^*γ^*$ amplitude. To fit the lattice data, we use two different modifications of vector meson dominance (VMD) ansatz. One is a combined form of VMD with the next-to-leading order (NLO) chiral perturbation theory (ChPT), and the other is a resummed form of pion-loop diagrams. We extract one of the low energy constants in NLO ChPT, and estimate $π^0\rightarrow γγ$ decay width including finite mass correction.

preprint2011arXiv

Workshop on Precision Measurements of alphas

These are the proceedings of the "Workshop on Precision Measurements of alphas" held at the Max-Planck-Institute for Physics, Munich, February 9-11, 2011. The workshop explored in depth the determination of alphas(mZ) in the MS-bar scheme from the key categories where high precision measurements are currently being made, including DIS and global PDF fits, tau-decays, electroweak precision observables and Z-decays, event-shapes, and lattice QCD. These proceedings contain a short summary contribution from the speakers, as well as the lists of authors, conveners, participants, and talks.

preprint2010arXiv

An application of the variational analysis to calculate the meson spectral functions

We present a new method to calculate meson spectral functions (SPFs) on the lattice based on a variational method. Because, on a finite volume lattice, the meson SPFs have discrete spectra only, a suitable way to extract such discrete signals is needed. Using a variational method, we can calculate several discrete quantities such as the position and the area of spectral peaks for low-lying states. Moreover data accuracy can be improved by increasing the number of basis functions. In this report, we first confirm our method in the free quark case and show that our method works well. Then, we apply the method to a quenched lattice QCD simulation and calculate the charmonium SPFs for S and P-waves at zero temperature. Our results for the ground state are well consistent with the position and the area of the lowest peaks of charmonium SPFs calculated by the conventional maximum entropy method. For first excited states, the signals may be reliablly extracted with our method because the charmonium mass converges to a value close to the experimental one when the number of basis functions is increased. We also investigate the SPFs for S-wave charmonia at below and above $T_c$. Our results suggest that $J/ψ$ and $η_c$ may survive up to 1.4$T_c$.

preprint2010arXiv

Calculation of $ρ$ meson decay width from the PACS-CS configurations

We present preliminary results on the $ρ$ meson decay width from $N_f=2+1$ full QCD configurations generated by PACS-CS Collaboration. The decay width is estimated from the $P$-wave scattering phase shift for the isospin $I=1$ two-pion system. The finite size formula presented by Lüscher in the center of mass frame and its extension to non-zero total momentum frame by Rummukainen and Gottlieb are employed for the calculations of the phase shift. Our calculations are carried out at $m_π=410\ {\rm MeV}$ ($m_π/m_ρ=0.46$) and $a=0.091\ {\rm fm}$ on a $32^3\times 64$ ($La=2.9 {\rm fm}$) lattice.

preprint2010arXiv

Determination of the chiral condensate from 2+1-flavor lattice QCD

We perform a precise calculation of the chiral condensate in QCD using lattice QCD with 2+1 flavors of dynamical overlap quarks. Up and down quark masses cover a range between 3 and 100 MeV on a 16^3x48 lattice at a lattice spacing around 0.11 fm. At the lightest sea quark mass, the finite volume system on the lattice is in the epsilon-regime. By matching the low-lying eigenvalue spectrum of the Dirac operator with the prediction of chiral perturbation theory at the next-to-leading order, we determine the chiral condensate in 2+1-flavor QCD with strange quark mass fixed at its physical value as Sigma (MS-bar at 2 GeV) = [242(04)(^+19_-18}) MeV}]^3, where the errors are statistical and systematic, respectively.

preprint2010arXiv

Electric and Magnetic Screening Masses at Finite Temperature from Generalized Polyakov-Line Correlations in Two-flavor Lattice QCD

Screenings of the quark-gluon plasma in electric and magnetic sectors are studied on the basis of generalized Polyakov-line correlation functions in lattice QCD simulations with two flavors of improved Wilson quarks. Using the Euclidean-time reflection ($\R$) and the charge conjugation ($\Ca$), electric and magnetic screening masses are extracted in a gauge invariant manner. Long distance behavior of the standard Polyakov-line correlation in the quark-gluon plasma is found to be dictated by the magnetic screening. Also, ratio of the two screening masses agrees with that obtained from the dimensionally-reduced effective field theory and the ${\cal N}=4$ supersymmetric Yang-Mills theory.

preprint2010arXiv

Energy dependence of nucleon-nucleon potentials in lattice QCD

Recently a new approach to calculate the nuclear potential from lattice QCD has been proposed. In the approach the nuclear potential is constructed from Bethe-Salpeter (BS) wave functons through the Schroedinger equation. The procedure leads to non-local but energy independent potential, which can be expanded in terms of local functions. In several recent applications of this method, local potentials, which correspond to the leading order (LO) terms of the expansion, are calculated from the BS wave function at E~0 MeV, where E is the center of mass energy. It is therefore important to check the validity of the LO approximation obtained at E~0. In this report, in order to check how well the LO approximation for the NN potentials works, we compare the LO potentials determined from the BS wave function at E~45 MeV with those at E~0 MeV in quenched QCD. We find that the difference of the LO potentials between two energies are not found wihin the statistical errors. This shows that the LO approximation for the potential is valid at low energies to describe the NN interactions.

preprint2010arXiv

EOS in 2+1 flavor QCD with improved Wilson quarks by the fixed-scale approach

We present the status of our study on the equation of state in 2+1 flavor QCD with non-perturbatively improved Wilson quarks coupled with the RG improved glue. We apply the T-integration method to non-perturbatively calculate the equation of state by the fixed-scale approach.

preprint2010arXiv

Equation of State and Heavy-Quark Free Energy at Finite Temperature and Density in Two Flavor Lattice QCD with Wilson Quark Action

We study the equation of state at finite temperature and density in two-flavor QCD with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 \times 4$ lattice. Along the lines of constant physics at $m_{\rm PS}/m_{\rm V} = 0.65$ and 0.80, we compute the second and forth derivatives of the grand canonical partition function with respect to the quark chemical potential $μ_q = (μ_u+μ_d)/2$ and the isospin chemical potential $μ_I = (μ_u-μ_d)/2$ at vanishing chemical potentials, and study the behaviors of thermodynamic quantities at finite $μ_q$ using these derivatives for the case $μ_I=0$. In particular, we study density fluctuations at none-zero temperature and density by calculating the quark number and isospin susceptibilities and their derivatives with respect to $μ_q$. To suppress statistical fluctuations, we also examine new techniques applicable at low densities. We find a large enhancement in the fluctuation of quark number when the density increased near the pseudo-critical temperature, suggesting a critical point at finite $μ_q$ terminating the first order transition line between hadronic and quark gluon plasma phases. This result agrees with the previous results using staggered-type quark actions qualitatively. Furthermore, we study heavy-quark free energies and Debye screening masses at finite density by measuring the first and second derivatives of these quantities for various color channels of heavy quark-quark and quark-anti-quark pairs. The results suggest that, to the leading order of $μ_q$, the interaction between two quarks becomes stronger at finite densities, while that between quark and anti-quark becomes weaker.

preprint2010arXiv

Kaon-Nucleon potential from lattice QCD

We study the KN interactions in the I(J^π)=0(1/2^-) and 1(1/2^-) channels and associated exotic state Θ^+ from 2+1 flavor full lattice QCD simulation for relatively heavy quark mass corresponding to m_π=871 MeV. The s-wave KN potentials are obtained from the Bethe-Salpeter wave function by using the method recently developed by HAL QCD (Hadrons to Atomic nuclei from Lattice QCD) Collaboration. Potentials in both channels reveal short range repulsions: Strength of the repulsion is stronger in the I=1 potential, which is consistent with the prediction of the Tomozawa-Weinberg term. The I=0 potential is found to have attractive well at mid range. From these potentials, the $KN$ scattering phase shifts are calculated and compared with the experimental data.

preprint2010arXiv

Light meson form factors in N_f=2+1 QCD with dynamical overlap quarks

We report on our calculation of pion and kaon form factors in three-flavor QCD using the overlap quark action. Gauge ensembles are generated on a 16^3 \times 48 lattice at a lattice spacing of 0.11 fm with pion masses down to 310 MeV. Connected and disconnected meson correaltors are calculated using the all-to-all quark propagator. We present our preliminary analysis on the chiral behavior of the electromagnetic and scalar form factors as well as a comparison of the shape of the K to pi form factors with experiment.

preprint2010arXiv

Measurement of the atmospheric muon charge ratio with the OPERA detector

The OPERA detector at the Gran Sasso underground laboratory (LNGS) was used to measure the atmospheric muon charge ratio in the TeV energy region. We analyzed 403069 atmospheric muons corresponding to 113.4 days of livetime during the 2008 CNGS run. We computed separately the muon charge ratio for single and for multiple muon events in order to select different energy regions of the primary cosmic ray spectrum and to test the charge ratio dependence on the primary composition. The measured charge ratio values were corrected taking into account the charge-misidentification errors. Data have also been grouped in five bins of the "vertical surface energy". A fit to a simplified model of muon production in the atmosphere allowed the determination of the pion and kaon charge ratios weighted by the cosmic ray energy spectrum.

preprint2010arXiv

Non-perturbative renormalization of quark mass in Nf=2+1 QCD with the Schroedinger functional scheme

We present an evaluation of the quark mass renormalization factor for Nf=2+1 QCD. The Schroedinger functional scheme is employed as the intermediate scheme to carry out non-perturbative running from the low energy region, where renormalization of bare mass is performed on the lattice, to deep in the high energy perturbative region, where the conversion to the renormalization group invariant mass or the MS-bar scheme is safely carried out. For numerical simulations we adopted the Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy regions and three lattice spacings to take the continuum limit at each scale. The regularization independent step scaling function of the quark mass for the Nf=2+1 QCD is obtained in the continuum limit. Renormalization factors for the pseudo scalar density and the axial vector current are also evaluated for the same action and the bare couplings as two recent large scale Nf=2+1 simulations; previous work of the CP-PACS/JLQCD collaboration, which covered the up-down quark mass range heavier than $m_π\sim 500$ MeV and that of PACS-CS collaboration for much lighter quark masses down to $m_π=155$ MeV. The quark mass renormalization factor is used to renormalize bare PCAC masses in these simulations.

preprint2010arXiv

Nucleon strange quark content in 2+1-flavor QCD

We calculate the strange quark content of the nucleon directly from its disconnected three-point function in N_f=2+1 QCD. Chiral symmetry is crucial to avoid a possibly large contamination due to operator mixing, and is exactly preserved by employing the overlap quark action. We also use the all-to-all quark propagator and the low-mode averaging technique in order to accurately calculate the relevant nucleon correlator. Our preliminary result extrapolated to the physical point is f_{T_s} = 0.013(12)(16). This is in good agreement with our previous estimate in N_f=2 QCD as well as those from our indirect calculations using the Feynman-Hellmann theorem.

preprint2010arXiv

Physical Point Simulation in 2+1 Flavor Lattice QCD

We present the results of the physical point simulation in 2+1 flavor lattice QCD with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action at $β=1.9$ on a $32^3 \times 64$ lattice. The physical quark masses together with the lattice spacing is determined with $m_π$, $m_K$ and $m_Ω$ as physical inputs. There are two key algorithmic ingredients to make possible the direct simulation at the physical point: One is the mass-preconditioned domain-decomposed HMC algorithm to reduce the computational cost. The other is the reweighting technique to adjust the hopping parameters exactly to the physical point. The physics results include the hadron spectrum, the quark masses and the pseudoscalar meson decay constants. The renormalization factors are nonperturbatively evaluated with the Schr{ö}dinger functional method. The results are compared with the previous ones obtained by the chiral extrapolation method.

preprint2010arXiv

QCD thermodynamics at zero and finite densities with improved Wilson quarks

The WHOT-QCD Collaboration is pushing forward lattice studies of QCD at finite temperatures and densities using improved Wilson quarks. We first present results on QCD at zero and finite densities with two flavors of degenerate quarks (N_F=2 QCD) adopting the conventional fixed-Nt approach. We then report on the status of a study of N_F=2+1 QCD adopting a fixed-scale approach armed with the T-integration method which we have developed.

preprint2010arXiv

The order of the deconfinement phase transition in a heavy quark mass region

We study the quark mass dependence of the QCD phase transition by an effective potential defined through the distribution function of observables. As a test of the method, we study the first order deconfinement phase transition in the heavy quark mass limit and its fate at lighter quark masses. We confirm that the distribution function for the plaquette has two peaks indicating that the phase transition is of first order in the heavy quark limit. We then study the quark mass dependence of the distribution function by a reweighting method combined with the hopping parameter expansion. We find that the first order transition turns into a crossover as the quark mass decreases. We determine the critical point for the cases of $N_f$=1, 2, 3 and 2+1. We find that the probability distribution function provides us with a powerful tool to study the order of transitions.

preprint2009arXiv

$π^0$ to two-photon decay in lattice QCD

We calculate the neutral pion ($π^0$) to off-shell two photon ($γ^*γ^*$) transition form factor in lattice QCD. The transition form factor can be extracted from the three-point function of the form (axial-vector)-(vector)-(vector) as a function of off-shell two-photon momentum. Since the axial-anomaly plays an important role in the $π^0\toγγ$ decay process, we employ the overlap fermion, which preserves the exact chiral symmetery on the lattice. After extrapolating to the chiral and the vanishing photon momentum limit with a fit function based on vector meson dominance (VMD) model, we find that the Adler-Bell-Jackiw anomaly is correctly reproduced.

preprint2009arXiv

Calculation of nucleon strange quark content with dynamical overlap quarks

We calculate the nucleon strange quark content directly from disconnected three-point functions. Numerical simulations are carried out in two-flavor QCD using the overlap quark action with up and down quark masses down to a fifth of the physical strange quark mass. To improve the statistical accuracy, we calculate the nucleon two-point functions with the low-mode averaging technique, whereas the all-to-all quark propagator is used for the disconnected quark loop. We obtain the y parameter, which is the ratio of the strange and light quark contents, y = 0.024(45) at the physical point. This is in a good agreement with our earlier calculation from the nucleon spectrum through the Feynman-Hellmann theorem.

preprint2009arXiv

Chiral properties of light mesons with $N_f=2+1$ overlap fermions

We present an update of the light meson spectrum with $N_f$=2+1 overlap fermions on a $16^3\times 48$ lattice at five different up and down quark masses and two strange quark masses. Based on our experience with the previous simulation with $N_f=2$, we carry out the chiral extrapolation with the prediction of the chiral perturbation theory at the next-to-next-to leading order. We also check the consistency of our analysis by using alternative chiral extrapolation with a reduced theory in which the strange quark mass is integrated out.

preprint2009arXiv

Flavor-singlet mesons in N_f=2+1 QCD with dynamical overlap quarks

We report on our study of flavor-singlet mesons in three-flavor QCD with dynamical overlap quarks. Gauge ensembles are generated on a 16^3 x 48 lattice at a lattice spacing of 0.10 fm with the strange quark masses around its physical value m_s,phys and up and down quark masses down to m_s,phys/5. Connected and disconnected meson correlators are calculated using the all-to-all quark propagator. We present our preliminary results on the spectrum of flavor-singlet pseudoscalar and vector mesons.

preprint2009arXiv

Heavy-quark free energy at finite temperature with 2+1 flavors of improved Wilson quarks in fixed scale approach

The free energy between a static quark and an antiquark is studied by using the color-singlet Polyakov-line correlation at finite temperature. We perform simulations on $32^3 \times 12$, 10, 8, 6, 4 lattices in the high temperature phase with the RG-improved gluon action and 2+1 flavors of the clover-improved Wilson quark action. Since the simulations are based on the fixed scale approach that the temperature can be varied without changing the spatial volume and renormalization factor, it is possible to investigate temperature dependence of the heavy-quark free energy without any adjustment of the overall constant. We find that, the heavy-quark free energies at short distance converge to the heavy-quark potential evaluated from the Wilson-loop operator at zero temperature, in accordance with the expected insensitivity of short distance physics to the temperature. At long distance, the heavy-quark free energies approach to twice the single-quark free energies, implying that the interaction between heavy quarks is screened. The Debye screening mass obtained from the long range behavior of the heavy-quark free energy is compared with results of the thermal perturbation theory and those of $N_f=2$ and $N_f=0$ lattice simulations.

preprint2009arXiv

Lambda-nucleon force from lattice QCD

We study the $Λ$-nucleon ($ΛN$) force by using lattice QCD. The Bethe-Salpeter amplitude is calculated for the lowest scattering state of the $ΛN$ so as to obtain the $ΛN$ potential. The numerical calculation is twofold: (i) Full lattice QCD by using 2+1 flavor PACS-CS gauge configurations with, $β=1.9$, corresponding to the lattice spacing of $a=0.0907(13)$ fm, on a $32^3\times 64$ lattice. A set of parameter $(κ_{ud},κ_s)=(0.13770,0.13640)$ is used, which corresponds to $m_π\approx 300$ MeV and $m_K\approx 594$ MeV. The spatial lattice volume corresponds to (2.86 fm)$^3$. (ii) Quenched lattice QCD with $β=5.7$, the lattice spacing of $a=0.1416(9)$ fm, on the $32^3\times48$ lattice. Two sets of hopping parameters $(κ_{ud},κ_s)=(0.1665,0.1643),(0.1670,0.1643)$ are used. The spatial lattice volume is (4.5 fm)$^3$. For the full QCD, we find that the $Λp$ has a relatively strong (weak) repulsive core in the $^1S_0$ ($^3S_1$) channel at short distance, while the potential has slight attractive region at medium distance. The lowest scattering energy in the finite lattice volume is calculated; Slightly negative values obtained in both spin channels. For the quenched QCD, we find that the results are qualitatively in agreement with those in the full QCD calculation.

preprint2009arXiv

Nuclear forces from quenched and 2+1 flavor lattice QCD using the PACS-CS gauge configurations

Two of recent progress in lattice QCD approach to nuclear force are reported. (i) Tensor force from quenched lattice QCD: By truncating the derivative expansion of inter-nucleon potential to the strictly local terms, we obtain central force V_C(r) and tensor force V_T(r) separately from s-wave and d-wave components of Bethe-Salpeter wave function for two nucleon state with J^P=1^+. Numerical calculation is performed with quenched QCD on 32^4 lattice using the standard plaquette action at beta=5.7 with the standard Wilson quark action with kappa=0.1640, 0.1665, 0.1678. Preliminary results show that the depths of the resulting tensor force amount to 20 to 40 MeV, which is enhanced in the light quark mass region. (ii) Nuclear force from 2+1 flavor QCD with PACS-CS gauge configuration: Preliminary full QCD results are obtained by using 2+1 flavor gauge configurations generated by PACS-CS collaboration. The resulting potential has the midium range attraction of about 30 MeV similar to the preceding quenched calculations. However, the repulsive core at short distance is significantly stronger than the corresponding quenched QCD result.

preprint2009arXiv

Nucleon sigma term and strange quark content in 2+1-flavor QCD with dynamical overlap fermions

We study the sigma term and the strange quark content of nucleon in 2+1-flavor QCD with dynamical overlap fermions. We analyze the lattice data of nucleon mass taken at two different strange quark masses with five values of up and down quark masses each. Using the reweighting technique, we study the strange quark mass dependence of the nucleon and extract the strange quark content.

preprint2009arXiv

Pion form factors from two-flavor lattice QCD with exact chiral symmetry

We calculate pion vector and scalar form factors in two-flavor lattice QCD and study the chiral behavior of the vector and scalar radii <r^2>_{V,S}. Numerical simulations are carried out on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with quark masses down to \sim m_s/6, where m_s is the physical strange quark mass. Chiral symmetry, which is essential for a direct comparison with chiral perturbation theory (ChPT), is exactly preserved in our calculation at finite lattice spacing by employing the overlap quark action. We utilize the so-called all-to-all quark propagator in order to calculate the scalar form factor including the contributions of disconnected diagrams and to improve statistical accuracy of the form factors. A detailed comparison with ChPT reveals that the next-to-next-to-leading-order contributions to the radii are essential to describe their chiral behavior in the region of quark mass from m_s/6 to m_s/2. Chiral extrapolation based on two-loop ChPT yields <r^2>_V=0.409(23)(37)fm and <r^2>_S=0.617(79)(66)fm, which are consistent with phenomenological analysis. We also present our estimates of relevant low-energy constants.

preprint2009arXiv

Precise determination of the strong coupling constant in Nf=2+1 lattice QCD with the Schrödinger functional scheme

We present an evaluation of the running coupling constant for Nf=2+1 QCD. The Schroedinger functional scheme is used as the intermediate scheme to carry out non-perturbative running from the low energy region, where physical scale is introduced, to deep in the high energy perturbative region, where conversion to the MS-bar scheme is safely performed. Possible systematic errors due to the use of perturbation theory occur only in the conversion from three-flavor to four-flavor running coupling constant near the charm mass threshold, where higher order terms beyond 5th order in the $β$ function may not be negligible. For numerical simulations we adopted Iwasaki gauge action and non-perturbatively improved Wilson fermion action with the clover term. Seven renormalization scales are used to cover from low to high energy region and three lattice spacings to take the continuum limit at each scale. A physical scale is introduced from the previous Nf=2+1 simulation of the CP-PACS/JL-QCD collaboration, which covered the up-down quark mass range heavier than $m_π\sim 500$ MeV.

preprint2009arXiv

Scintillator counters with WLS fiber/MPPC readout for the side muon range detector (SMRD)of the T2K experiment

The T2K neutrino experiment at J-PARC uses a set of near detectors to measure the properties of an unoscillated neutrino beam and neutrino interaction cross-sections. One of the sub-detectors of the near-detector complex, the side muon range detector (SMRD), is described in the paper. The detector is designed to help measure the neutrino energy spectrum, to identify background and to calibrate the other detectors. The active elements of the SMRD consist of 0.7 cm thick extruded scintillator slabs inserted into air gaps of the UA1 magnet yokes. The readout of each scintillator slab is provided through a single WLS fiber embedded into a serpentine shaped groove. Two Hamamatsu multi-pixel avalanche photodiodes (MPPC's) are coupled to both ends of the WLS fiber. This design allows us to achieve a high MIP detection efficiency of greater than 99%. A light yield of 25-50 p.e./MIP, a time resolution of about 1 ns and a spatial resolution along the slab better than 10 cm were obtained for the SMRD counters.

preprint2009arXiv

Towards the equation of state in 2+1 flavor QCD with improved Wilson quarks in the fixed scale approach

We report on the status of our study towards the equation of state in 2+1 flavor QCD with improved Wilson quarks. To reduce the computational cost which is quite demanding for Wilson-type quarks, we adopt the fixed scale approach, i.e. the temperature T is varied by N_t at fixed lattice spacing. Since the conventional integral method to obtain the pressure is inapplicable at a fixed scale, we adopt the "T-integral method", to calculate the pressure non-perturbatively. Reduction of the computational cost of T=0 simulations thus achieved is indispensable to study EOS in QCD with dynamical quarks.

preprint2009arXiv

Universality of the Nf=2 Running Coupling

We investigate universality of the Nf=2 running coupling in the Schödinger functional scheme, by calculating the step scaling function in lattice QCD with the renorm alization group (RG) improved gauge action at both weak(u=0.9796) and strong(u=3.3340) couplings, where u=\bar{g}^2_SF with \bar{g}_SF being the running coupling in this scheme. In our main calculations, we use the tree-level value for O(a) improvement coefficients of boundary gauge fields. In addition we employ the 1-loop value for them in order to see how scaling behaviours are affected by them. In the continuum limit, the step scaling function obtained from the RG improved gauge actions agrees with the previous result obtained from the plaquette action within errors at both couplings, though errors of our result are larger. Combined fits using all data with the RG improved action as well as the plaquette action reduce errors in the continuum limit by 2% at the weak coupling and 22% at the strong coupling.

preprint2008arXiv

2+1 flavor lattice QCD simulation with $O(a)$-improved Wilson quarks

We present simulation details and results for the light hadron spectrum in N f = 2 + 1 lattice QCD with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing of 0.09 fm on a (2.9fm)^3 box using the PACS-CS computer. We employ the Luscher's domain-decomposed HMC algorithm with several improvements to reduce the degenerate up-down quark mass toward the physical value. So far the resulting pseudoscalar meson mass is ranging from 702MeV down to 156MeV. We discuss on the stability and the efficiency of the algorithm. The light harden spectrum extrapolated at the physical point is compared with the experimental values. We also present the values of the quark masses and the pseudoscalar meson decay constants.

preprint2008arXiv

2+1 Flavor Lattice QCD toward the Physical Point

We present the first results of the PACS-CS project which aims to simulate 2+1 flavor lattice QCD on the physical point with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Numerical simulations are carried out at the lattice spacing of $a=0.0907(13)$fm on a $32^3\times 64$ lattice with the use of the DDHMC algorithm to reduce the up-down quark mass. Further algorithmic improvements make possible the simulation whose ud quark mass is as light as the physical value. The resulting PS meson masses range from 702MeV down to 156MeV, which clearly exhibit the presence of chiral logarithms. An analysis of the PS meson sector with SU(3) ChPT reveals that the NLO corrections are large at the physical strange quark mass. In order to estimate the physical ud quark mass, we employ the SU(2) chiral analysis expanding the strange quark contributions analytically around the physical strange quark mass. The SU(2) LECs ${\bar l}_3$ and ${\bar l}_4$ are comparable with the recent estimates by other lattice QCD calculations. We determine the physical point together with the lattice spacing employing $m_π$, $m_K$ and $m_Ω$ as input. The hadron spectrum extrapolated to the physical point shows an agreement with the experimental values at a few % level of statistical errors, albeit there remain possible cutoff effects. We also find that our results of $f_π=134.0(4.2)$MeV, $f_K=159.4(3.1)$MeV and $f_K/f_π=1.189(20)$ with the perturbative renormalization factors are compatible with the experimental values. For the physical quark masses we obtain $m_{\rm ud}^\msbar=2.527(47)$MeV and $m_{\rm s}^\msbar=72.72(78)$MeV extracted from the axial-vector Ward-Takahashi identity with the perturbative renormalization factors.

preprint2008arXiv

Charm quark system in 2+1 flavor lattice QCD using the PACS-CS configurations

We study heavy-heavy and heavy-light quark systems for charm with a relativistic heavy quark action in 2+1 flavor lattice QCD. Configurations are generated by the PACS-CS Collaboration at the lattice spacing is $a=0.09$ fm with the lattice size of $32^3\times 64$ employing the $O(a)$-improved Wilson quark action and the Iwasaki gauge action. We present preliminary results for the charmonium spectrum and the $D$ and $D_s$ meson decay constants evaluated at 3.5 MeV$< m_{\rm ud}<$ 12 MeV with $m_{\rm s}$ around the physical value. We investigate the dynamical quark mass dependences of the hyperfine and the orbital splittings. The decay constants are compared with the recent experimental values.

preprint2008arXiv

Equation of state at finite density in two-flavor QCD with improved Wilson quarks

We study the equation of state in two-flavor QCD at finite temperature and density. Simulations are made with the RG-improved gluon action and the clover-improved Wilson quark action. Along the lines of constant physics for $m_{\rm PS}/m_{\rm V} = 0.65$ and 0.80, we compute the derivatives of the quark determinant with respect to the quark chemical potential $μ_q$ up to the fourth order at $μ_q=0$. We adopt several improvement techniques in the evaluation. We study thermodynamic quantities and quark number susceptibilities at finite $μ_q$ using these derivatives. We find enhancement of the quark number susceptibility at finite $μ_q$, in accordance with previous observations using staggered-type quarks. This suggests the existence of a nearby critical point.

preprint2008arXiv

Light meson spectrum with $N_f=2+1$ dynamical overlap fermions

We report on a numerical simulation with 2+1 dynamical flavors of overlap fermions. We calculate pseudo-scalar masses and decay constants on a $16^3\times 48 \times (0.11 {\rm fm})^4$ lattice at five different up and down quark masses and two strange quark masses. The lightest pion mass corresponds to $\approx 310$ MeV. We also study the validity of the chiral perturbation theory using the results of the numerical simulation with two dynamical flavors and conclude that the one-loop formulae cannot be directly applied in the strange quark mass region. We therefore extrapolate our 2+1-flavor results to the chiral limit by fitting the data to the two-loop formulae of the chiral perturbation theory.

preprint2008arXiv

Light quark masses from unquenched lattice QCD

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.

preprint2008arXiv

Magnetic and electric screening masses from Polyakov-loop correlations

Screening properties of the quark gluon plasma are studied from Polyakov-loop correlation in lattice QCD simulations with two flavors of improved Wilson quarks at temperatures $T/\Tpc \simeq 1$--4 where $\Tpc$ is the pseudocritical temperature. Using the Euclidean-time reflection symmetry and the charge conjugation symmetry, we introduce various types of Polyakov-loop correlation functions and extract screening masses in magnetic and electric sectors. We find that the temperature dependence of the screening masses are well described by the weak coupling expansion. We also find that a ratio of the screening masses in the electric sector to the magnetic sector shows qualitative agreement with a prediction from the dimensionally-reduced effective field theory and the N=4 supersymmetric Yang-Mills theory at $1.3 < T/\Tpc < 3$.

preprint2008arXiv

Pion vector and scalar form factors with dynamical overlap quarks

We calculate the pion vector and scalar form factors in two-flavor QCD. Gauge configurations are generated with dynamical overlap quarks on a 16^3 x 32 lattice at a lattice spacing of 0.12 fm with sea quark masses down to a sixth of the physical strange quark mass. Contributions of disconnected diagrams to the scalar form factor is calculated employing the all-to-all quark propagators. We present a detailed comparison of the vector and scalar radii with chiral perturbation theory to two loops.

preprint2008arXiv

SU(2) and SU(3) chiral perturbation theory analyses on meson and baryon masses in 2+1 flavor lattice QCD

We investigate the quark mass dependence of meson and baryon masses obtained from 2+1 flavor dynamical quark simulations performed by the PACS-CS Collaboration. With the use of SU(2) and SU(3) chiral perturbation theories up to NLO, we examine the chiral behavior of the pseudoscalar meson masses and the decay constants in terms of the degenerate up-down quark mass ranging form 3 MeV to 24 MeV and two choices of the strange quark mass around the physical value. We discuss the convergence of the SU(2) and SU(3) chiral expansions and present the results for the low energy constants. We find that the SU(3) expansion is not convergent at NLO for the physical strange quark mass. The chiral behavior of the nucleon mass is also discussed based on the SU(2) heavy baryon chiral perturbation theory up to NNLO. Our results show that the expansion is well behaved only up to m_pi^2 ~ 0.2 GeV^2.

preprint2008arXiv

Thermodynamics of SU(3) gauge theory at fixed lattice spacing

We study thermodynamics of SU(3) gauge theory at fixed scales on the lattice, where we vary temperature by changing the temporal lattice size N_t=(Ta_t)^{-1}. In the fixed scale approach, finite temperature simulations are performed on common lattice spacings and spatial volumes. Consequently, we can isolate thermal effects in observables from other uncertainties, such as lattice artifact, renormalization factor, and spatial volume effect. Furthermore, in the EOS calculations, the fixed scale approach is able to reduce computational costs for zero temperature subtraction and parameter search to find lines of constant physics, which are demanding in full QCD simulations. As a test of the approach, we study the thermodynamics of the SU(3) gauge theory on isotropic and anisotropic lattices. In addition to the equation of state, we calculate the critical temperature and the static quark free energy at a fixed scale.

preprint2008arXiv

Topological susceptibility in (2+1)-flavor lattice QCD with overlap fermion

We determine the topological susceptibility χ_t in the topologically-trivial sector generated by lattice simulations of N_f = 2+1 QCD with overlap Dirac fermion, on a 16^3 x 48 lattice with lattice spacing ~ 0.11 fm, for five sea quark masses m_q ranging from m_s/6 to m_s (where m_s is the physical strange quark mass). The χ_t is extracted from the plateau (at large time separation) of the 2-point and 4-point time-correlation functions of the flavor-singlet pseudoscalar meson η', which arises from the finite size effect due to fixed topology. In the small m_q regime, our result of χ_t agrees with the chiral effective theory. Using the formula χ_t = Σ(m_u^{-1} + m_d^{-1} + m_s^{-1})^{-1} by Leutwyler-Smilga, we obtain the chiral condensate Σ^{MSbar}(2 GeV) = [249(4)(2) MeV]^3.

preprint2007arXiv

2+1 flavor light hadron spectrum and quark masses with the O(a) improved Wilson-clover quark formalism

We present a summary of results of the joint CP-PACS and JLQCD project toward a 2+1 flavor full QCD simulation with the O(a)-improved Wilson quark formalism and the Iwasaki gauge action. Configurations were generated during 2002-2005 at three lattice spacings, a~0.076, 0.100 and 0.122 fm, keeping the physical volume constant at (2.0fm)^3. Up and down quark masses are taken in the range m_{PS}/m_V~0.6-0.78. We have completed the analysis for the light meson spectrum and quark masses in the continuum limit using the full configuration set. The predicted meson masses reproduce experimental values in the continuum limit at a 1% level. The average up and down, and strange quark masses turn out to be m_{ud}^{\bar{MS}}(μ=2 GeV)=3.50(14)({}^{+26}_{-15}) MeV and m_s^{\bar{MS}}(μ=2 GeV)=91.8(3.9)({}^{+6.8}_{-4.1}) MeV. We discuss our future strategy toward definitive results on hadron spectroscopy with the Wilson-clover formalism.

preprint2007arXiv

Application of chiral perturbation theory to 2+1 flavor lattice QCD with O(a)-improved Wilson quarks

We apply chiral perturbation theory to the pseudoscalar meson mass and decay constant data obtained in the PACS-CS Project toward 2+1 flavor lattice QCD simulations with the O(a)-improved Wilson quarks. We examine the existence of chiral logarithms in the quark mass range from m_{ud}=47 MeV down to 6 MeV on a (2.8 fm)^3 box with the lattice spacing a=0.09 fm. Several low energy constants are determined. We also discuss the magnitude of finite size effects based on chiral perturbation theory.

preprint2007arXiv

First Nonperturbative Test of a Relativistic Heavy Quark Action in Quenched Lattice QCD

We perform a numerical test of a relativistic heavy quark(RHQ) action, recently proposed by Tsukuba group, in quenched lattice QCD at $a\simeq 0.1$ fm. With the use of the improvement parameters previously determined at one-loop level for the RHQ action, we investigate a restoration of rotational symmetry for heavy-heavy and heavy-light meson systems around the charm quark mass. We focused on two quantities, the meson dispersion relation and the pseudo-scalar meson decay constants. It is shown that the RHQ action significantly reduces the discretization errors due to the charm quark mass. We also calculate the S-state hyperfine splittings for the charmonium and charmed-strange mesons and the $D_s$ meson decay constant. The remaining discretization errors in the physical quantities are discussed.

preprint2007arXiv

Heavy-Quark Free Energy, Debye Mass, and Spatial String Tension at Finite Temperature in Two Flavor Lattice QCD with Wilson Quark Action

We study Polyakov loop correlations and spatial Wilson loop at finite Temperature in two-flavor QCD simulations with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 \times 4$ lattice. From the line of constant physics at $m_{\rm PS}/m_{\rm V}=0.65$ and 0.80, we extract the heavy-quark free energies, the effective running coupling $g_{\rm eff}(T)$ and the Debye screening mass $m_D(T)$ for various color channels of heavy quark--quark and quark--anti-quark pairs above the critical temperature. The free energies are well approximated by the screened Coulomb form with the appropriate Casimir factors at high temperature. The magnitude and the temperature dependence of the Debye mass are compared to those of the next-to-leading order thermal perturbation theory and to a phenomenological formula in terms of $g_{\rm eff}(T)$. We make a comparison between our results with the Wilson quark action and the previous results with the staggered quark action. The spatial string tension is also studied in the high temperature phase and is compared to the next-to-next-leading order prediction in an effective theory with dimensional reduction.

preprint2007arXiv

Lattice QCD Calculation of the $ρ$ Meson Decay Width

We present a lattice QCD calculation of the $ρ$ meson decay width via the $P$-wave scattering phase shift for the I=1 two-pion system. Our calculation uses full QCD gauge configurations for $N_f=2$ flavors generated using a renormalization group improved gauge action and an improved Wilson fermion action on a $12^3\times24$ lattice at $m_π/m_ρ=0.41$ and the lattice spacing $1/a=0.92 {\rm GeV}$. The phase shift calculated with the use of the finite size formula for the two-pion system in the moving frame shows a behavior consistent with the existence of a resonance at a mass close to the vector meson mass obtained in spectroscopy. The decay width estimated from the phase shift is consistent with the experiment, when the quark mass is scaled to the realistic value.

preprint2007arXiv

Light hadron spectrum with 2+1 flavor dynamical $O(a)$-improved Wilson quarks

We present preliminary results for the light harden spectrum in $N_f=2+1$ lattice QCD obtained with the nonperturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at $β=1.90$ on a $32^3 \times 64$ lattice using the PACS-CS computer. We employ Lüscher's domain-decomposed HMC algorithm to reduce the up-down quark masses toward the physical value. The pseudoscalar meson masses range from 730 MeV down to 210 MeV. We compare the light harden spectrum extrapolated to the physical point with the experimental values.

preprint2007arXiv

Non-perturbative improvement of the axial current with three dynamical flavors and the Iwasaki gauge action

We perform a non-perturbative determination of the improvement coefficient c_A to remove O(a) discretization errors in the axial vector current in three-flavor lattice QCD with the Iwasaki gauge action and the standard O$(a)$-improved Wilson quark action. An improvement condition with a good sensitivity to c_A is imposed at constant physics. Combining our results with the perturbative expansion, c_A is now known rather precisely for 1/a \gtrsim 1.6 GeV.

preprint2007arXiv

Two-flavor lattice QCD in the epsilon-regime and chiral Random Matrix Theory

The low-lying eigenvalue spectrum of the QCD Dirac operator in the epsilon-regime is expected to match with that of chiral Random Matrix Theory (ChRMT). We study this correspondence for the case including sea quarks by performing two-flavor QCD simulations on the lattice. Using the overlap fermion formulation, which preserves exact chiral symmetry at finite lattice spacings, we push the sea quark mass down to \sim 3 MeV on a 16^3\times 32 lattice at a lattice spacing a \simeq 0.11 fm. We compare the low-lying eigenvalue distributions and find a good agreement with the analytical predictions of ChRMT. By matching the lowest-lying eigenvalue we extract the chiral condensate, Σ(2 GeV)[MSbar] = [251(7)(11) MeV]^3, where errors represent statistical and higher order effects in the epsilon expansion. We also calculate the eigenvalue distributions on the lattices with heavier sea quarks at two lattice spacings. Although the epsilon expansion is not applied for those sea quarks, we find a reasonable agreement of the Dirac operator spectrum with ChRMT. The value of Sigma, after extrapolating to the chiral limit, is consistent with the estimate in the epsilon-regime.

preprint2006arXiv

$ρ$ meson decay from the lattice

We present preliminary results on the $ρ$ meson decay width estimated from the scattering phase shift of the I=1 two-pion system. The phase shift is calculated by the finite size formula for non-zero total momentum frame (the moving frame) derived by Rummukainen and Gottlieb, using the $N_f=2$ improved Wilson fermion action at $m_π/m_ρ=0.41$ and $L=2.53 {\rm fm}$.

preprint2006arXiv

2+1 Flavor Lattice QCD with Luescher's Domain-Decomposed HMC Algorithm

We report on a study of 2+1 flavor lattice QCD with the $O(a)$-improved Wilson quarks on a $16^3\times 32$ lattice at the lattice spacing $1/a\approx 2$GeV employing Luescher's domain-decomposed HMC(LDDHMC) algorithm. This is dedicated to a preliminary study for the PACS-CS project which plans to complete the Wilson-clover $N_f=2+1$ program lowering the up-down quark masses close to the physical values as much as possible. We focus on three issues: (i) how light quark masses we can reach with LDDHMC, (ii) efficiency of the algorithm compared with the conventional HMC, (iii) parameter choice for the production runs on PACS-CS.

preprint2006arXiv

An application of the UV-filtering preconditioner to the Polynomial Hybrid Monte Carlo algorithm

We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the $O(a)$-improved Wilson action at $β=5.2, c_{\mathrm{sw}}=2.02, M_π/M_ρ\sim 0.8$ and 0.7 on a $16^3\times 48$ lattice. We find that the UV-filtering reduces the magnitude of the molecular dynamics force from the pseudo fermion by a factor 3 by tuning the UV-filter parameter. Combining with the multi-time scale molecular dynamics integrator we achieve a factor 2 improvement.

preprint2006arXiv

An estimate of the eta and eta-prime meson masses in Nf=2+1 lattice QCD

Masses of the eta and eta-prime mesons are estimated in Nf=2+1 lattice QCD with the non-perturbatively O(a) improved Wilson quark action and the Iwasaki RG-improved gluon action, using CP-PACS/JLQCD configurations on a 16^3 x 32 lattice at beta=1.83 (lattice spacing is 0.122 fm). We apply a stochastic noise estimator technique combined with smearing method to evaluate correlators among flavor SU(2) singlet pseudoscalar operators and strange pseudoscalar operators for 10 combinations of up/down and strange quark masses. The correlator matrix is then diagonalized to identify signals for mass eigenstates. Masses of the ground state and the first excited state extrapolated to the physical point are m_eta= 0.545(16) GeV and m_eta-prime= 0.871(46) GeV, being close to the experimental values of the eta and eta-prime masses.

preprint2006arXiv

Kaon $B$-parameters for Generic $ΔS=2$ Four-Quark Operators in Quenched Domain Wall QCD

We present a study of $B$-parameters for generic $ΔS=2$ four-quark operators in domain wall QCD. Our calculation covers all the $B$-parameters required to study the neutral kaon mixing in the standard model (SM) and beyond it. We evaluate one-loop renormalization factors of the operators employing the plaquette and Iwasaki gauge actions. Numerical simulations are carried out in quenched QCD with both gauge actions on $16^3\times 32\times 16$ and $24^3\times 32\times 16$ at the lattice spacing $1/a\approx 2$GeV. We investigate the relative magnitudes of the non-SM $B$-parameters to the SM one, which are compared with the previous results obtained with the overlap and the clover quark actions.

preprint2006arXiv

Nonperturbative $O(a)$ improvement of the Wilson quark action with the RG-improved gauge action using the Schrödinger functional method

We perform a nonperturbative determination of the $O(a)$-improvement coefficient $c_{\rm SW}$ and the critical hopping parameter $κ_c$ for $N_f$=3, 2, 0 flavor QCD with the RG-improved gauge action using the Schrödinger functional method. In order to interpolate $c_{\rm SW}$ and $κ_c$ as a function of the bare coupling, a wide range of $β$ from the weak coupling region to the moderately strong coupling points used in large-scale simulations is studied. Corrections at finite lattice size of $O(a/L)$ turned out to be large for the RG-improved gauge action, and hence we make the determination at a size fixed in physical units using a modified improvement condition. This enables us to avoid $O(a)$ scaling violations which would remain in physical observables if $c_{\rm SW}$ determined for a fixed lattice size $L/a$ is used in numerical simulations.

preprint2005arXiv

I=2 Pion Scattering Length from Two-Pion Wave Functions

We calculate the two-pion wave function in the ground state of the I=2 $S$-wave system and find the interaction range between two pions, which allows us to examine the validity of the necessary condition for the finite-volume method for the scattering length proposed by Lüscher. We work in the quenched approximation employing a renormalization group improved gauge action for gluons and an improved Wilson action for quarks at $1/a=1.207(12) {\rm GeV}$ on $16^3 \times 80$, $20^3 \times 80$ and $24^3 \times 80$ lattices. We conclude that the necessary condition is satisfied within the statistical errors for the lattice sizes $L\ge 24$ ($3.92 {\rm fm}$) when the quark mass is in the range that corresponds to $m_π^2 = 0.273-0.736 {\rm GeV}^2$. We obtain the scattering length with a smaller statistical error from the wave function than from the two-pion time correlator.

preprint2005arXiv

Light hadron spectrum and quark masses in 2+1 flavor QCD

CP-PACS and JLQCD collaborations are carrying out a joint project of the 2+1 flavor full QCD simulation. Gauge configurations are generated for the non-perturbatively $O(a)$-improved Wilson quark action and the Iwasaki gauge action using PHMC algorithm at three lattice spacings, $a\sim 0.076$, 0.010 and 0.122 fm, with a fixed physical volume $(2.0 fm)^3$. We present analysis for the light meson spectrum and quark masses in the continuum limit, which are determined using data obtained from the simulations at the two coarser lattices. Our simulations reproduce experimental values of meson masses. The ud and strange quark masses turn out to be $m_{ud}^{\bar{MS}}(μ=2 GeV)=3.34(23) MeV$ and $m_s^{\bar{MS}}(μ=2 GeV)=86.7(5.9) MeV$. We also show preliminary results at our finest lattice spacing for which simulations are still being continued.

preprint2005arXiv

Non-perturbative O(a)-improvement of Wilson quark action in three-flavor QCD with plaquette gauge action

We perform a non-perturbative determination of the O(a)-improvement coefficient c_SW for the Wilson quark action in three-flavor QCD with the plaquette gauge action. Numerical simulations are carried out in a range of β=12.0-5.2 on a single lattice size of 8^3x16 employing the Schrödinger functional setup of lattice QCD. As our main result, we obtain an interpolation formula for c_SW and the critical hopping parameter K_c as a function of the bare coupling. This enables us to remove O(a) scaling violation from physical observables in future numerical simulation in the wide range of β. Our analysis with a perturbatively modified improvement condition for c_SW suggests that finite volume effects in c_SW are not large on the 8^3x16 lattice. We investigate N_f dependence of c_SW by additional simulations for N_f=4, 2 and 0 at β=9.6. As a preparatory step for this study, we also determine c_SW in two-flavor QCD at β=5.2. At this β, several groups carried out large-scale calculations of the hadron spectrum, while no systematic determination of c_SW has been performed.

preprint2005arXiv

The PACS-CS Project

We describe our plan to develop a large-scale cluster system with a peak speed of 14.3Tflops for lattice QCD at the Center for Computational Sciences, University of Tsukuba, as a successor to the current 0.6Tflops CP-PACS computer. The system consist of 2560 nodes connected by a 16x16x10 three-dimensional hyper crossbar network. Each node has a single low-voltage 2.8GHz Xeon processor and 2GBytes of memory with 6.4GBytes/sec bandwidth, and 160 GBytes of disk in RAID1 mode. The network link in each of the three directions is made of dual Gigabit Ethernet with the peak throughput of 250MByte/sec. Hence each node has an aggregate network bandwidth of 750MByte/sec. The system will run under Linux and SCore, and an extension of the PM driver is developed for the network. The system will be developed jointly with Hitachi Limited. The installation is scheduled in the first quarter of Japanese Fiscal 2006 (April-June 2006) and the start of operation is expected in July 2006.

preprint2004arXiv

A scaling study of the step scaling function in SU(3) gauge theory with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the renormalization-group improved Iwasaki gauge action and the perturbatively improved Lüscher-Weisz gauge action. We confirm that the step scaling functions from the improved gauge actions agree with that previously obtained from the plaquette action within errors in the continuum limit at both weak and strong coupling regions. We also investigate how different choices of boundary counter terms for the improved gauge actions affect the scaling behavior. In the extrapolation to the continuum limit, we observe that the cut off dependence becomes moderate for the Iwasaki action, if a perturbative reduction of scaling violations is applied to the simulation results. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

A scaling study of the step scaling function of quenched QCD with improved gauge actions

We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the Iwasaki gauge action and the Luescher-Weisz gauge action. In particular, we test the choice of boundary counter terms and apply a perturbative procedure for removal of lattice artifacts for the simulation results in the extrapolation procedure. We confirm the universality of the step scaling functions at both weak and strong coupling regions. We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality.

preprint2004arXiv

Bulk first-order phase transition in three-flavor lattice QCD with $O(a)$-improved Wilson fermion action at zero temperature

Three-flavor QCD simulation with the $O(a)$-improved Wilson fermion action is made employing an exact fermion algorithm developed for odd number of quark flavors. For the plaquette gauge action, an unexpected first-order phase transition is found in the strong coupling regime ($β\lesssim$ 5.0) at relatively heavy quark masses ($m_{\mathrm{PS}}/m_{\mathrm{V}}\sim$ 0.74--0.87). Strong metastability persists on a large lattice of size $12^3\times 32$, which indicates that the transition has a bulk nature. The phase gap becomes smaller toward weaker couplings and vanishes at $β\simeq 5.0$, which corresponds to a lattice spacing $a\simeq$ 0.1 fm. The phase transition is not found if the improved gauge actions are employed. Our results imply that realistic simulations of QCD with three flavors of dynamical Wilson-type fermions at lattice spacings in the range $a=$ 0.1--0.2 fm require use of improved gauge actions. Possible origins of the phase transition is discussed.

preprint2004arXiv

Chiral extrapolations with small sea quark mass data in two-flavor lattice QCD

We present results on the light hadron spectrum and quark mass in two-flavor QCD calculated with small sea quark masses down to $m_{PS}/m_{V}=0.35$. The configurations are generated using the RG improved gauge and tadpole-improved clover quark action at $β=1.8$, where $a^{-1} \simeq 1$ GeV. Chiral extrapolations are made using not only polynomials and ChPT in the continuum but also formulae of Wilson chiral perturbation theory (WChPT) including $O(a^2)$ chiral breaking terms. We examine the viability of WChPT and its influence on quark masses.

preprint2004arXiv

I=2 $ππ$ Scattering Phase Shift with two Flavors of $O(a)$ Improved Dynamical Quarks

We present a lattice QCD calculation of phase shift including the chiral and continuum extrapolations in two-flavor QCD. The calculation is carried out for I=2 S-wave $ππ$ scattering. The phase shift is evaluated for two momentum systems, the center of mass and laboratory systems, by using the finite volume method proposed by Lüscher in the center of mass system and its extension to general systems by Rummukainen and Gottlieb. The measurements are made at three different bare couplings $β= 1.80$, 1.95 and 2.10 using a renormalization group improved gauge and a tadpole improved clover fermion action, and employing a set of configurations generated for hadron spectroscopy in our previous work. The illustrative values we obtain for the phase shift in the continuum limit are $δ$(deg.) $= - 3.50(64)$, $ - 9.5(30)$ and $ - 16.9(64)$ for $\sqrt{s}({\rm GeV})$ $=0.4$, $ 0.6$ and $ 0.8$, which are consistent with experiment.

preprint2004arXiv

I=2 Pion Scattering Length from Two-Pion Wave Function

We present a report on a calculation of scattering length for I=2 $S$-wave two-pion system from two-pion wave function. Calculations are made with an RG-improved action for gluons and improved Wilson action for quarks at $a^{-1}=1.207(12) {\rm GeV}$ on $16^3 \times 80$, $20^3 \times 80$ and $24^3 \times 80$ lattices. We investigate the validity of necessary condition for application of Lüscher's formula through the wave function. We find that the condition is satisfied for lattice volumes $L\ge 3.92 {\rm fm}$ for the quark mass range $m_π^2 = 0.273-0.736 {\rm GeV}^2$. We also find that the scattering length can be extracted with a smaller statistical error from the wave function than with a time correlation function used in previous studies.

preprint2004arXiv

Lattice QCD calculation of the proton decay matrix element in the continuum limit

We present a quenched lattice QCD calculation of the αand βparameters of the proton decay matrix element. The simulation is carried out using the Wilson quark action at three values of the lattice spacing in the range a\approx 0.1-0.064 fm to study the scaling violation effect. We find only mild scaling violation when the lattice scale is determined by the nucleon mass. We obtain in the continuum limit, |α(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |β(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with αand βin a relatively opposite sign, where the first error is statistical and the second is due to the uncertainty in the determination of the physical scale.

preprint2004arXiv

Light hadron spectroscopy in two-flavor QCD with small sea quark masses

We extend the study of the light hadron spectrum and the quark mass in two-flavor QCD to smaller sea quark mass, corresponding to $m_{PS}/m_{V}=0.60$--0.35. Numerical simulations are carried out using the RG-improved gauge action and the meanfield-improved clover quark action at $β=1.8$ ($a = 0.2$ fm from $ρ$ meson mass). We observe that the light hadron spectrum for small sea quark mass does not follow the expectation from chiral extrapolations with quadratic functions made from the region of $m_{PS}/m_{V}=0.80$--0.55. Whereas fits with either polynomial or continuum chiral perturbation theory (ChPT) fails, the Wilson ChPT (WChPT) that includes $a^2$ effects associated with explicit chiral symmetry breaking successfully fits the whole data: In particular, WChPT correctly predicts the light quark mass spectrum from simulations for medium heavy quark mass, such as $m_{PS}/m_V \simgt 0.5$. Reanalyzing the previous data %at $m_{PS}/m_{V}=0.80$--0.55 with the use of WChPT, we find the mean up and down quark mass being smaller than the previous result from quadratic chiral extrapolation by approximately 10%, $m_{ud}^{\bar{\rm MS}}(μ=2 {GeV}) = 3.11(17)$ [MeV] in the continuum limit.

preprint2004arXiv

Non-perturbative renormalization of meson decay constants in quenched QCD for a renormalization group improved gauge action

Renormalization constants ($Z$-factors) of vector and axial-vector currents are determined non-perturbatively in quenched QCD for a renormalization group improved gauge action and a tadpole improved clover quark action using the Schrödinger functional method. Non-perturbative values of $Z$-factors turn out to be smaller than one-loop perturbative values by $O(15%)$ at lattice spacing of $a^{-1}\approx$ 1 GeV. The pseudoscalar and vector meson decay constants calculated with the non-perturbative $Z$-factors show a much better scaling behavior compared to previous results obtained with tadpole improved one-loop $Z$-factors. In particular, the non-perturbative $Z$-factors normalized at infinite physical volume show that scaling violation of the decay constants are within about 10% up to the lattice spacing $a^{-1}\sim 1$ GeV. The continuum estimates obtained from data in the range $a^{-1}=$ 1 -- 2 GeV agree with those determined from finer lattices ($a^{-1}\sim 2-4$ GeV) with the standard action.

preprint2003arXiv

An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion

We propose an exact simulation algorithm for lattice QCD with dynamical Kogut-Susskind fermion in which the N_f-flavor fermion operator is defined as the N_f/4-th root of the Kogut-Susskind (KS) fermion operator. The algorithm is an extension of the Polynomial Hybrid Monte Carlo (PHMC) algorithm to KS fermions. The fractional power of the KS fermion operator is approximated with a Hermitian Chebyshev polynomial, with which we can construct an algorithm for any number of flavors. The error which arises from the approximation is corrected by the Kennedy-Kuti noisy Metropolis test. Numerical simulations are performed for the two-flavor case for several lattice parameters in order to confirm the validity and the practical feasibility of the algorithm. In particular tests on a 16^4 lattice with a quark mass corresponding to m_{PS}/m_V ~ 0.68 are successfully accomplished. We conclude that our algorithm provides an attractive exact method for dynamical QCD simulations with KS fermions.

preprint2003arXiv

Calculation of Non-Leptonic Kaon Decay Amplitudes from $K\toπ$ Matrix Elements in Quenched Domain-Wall QCD

We explore application of the domain wall fermion formalism of lattice QCD to calculate the $K\toππ$ decay amplitudes in terms of the $K\toπ$ and $K\to 0$ hadronic matrix elements through relations derived in chiral perturbation theory. Numerical simulations are carried out in quenched QCD using domain-wall fermion action for quarks and an RG-improved gauge action for gluons on a $16^3\times 32\times 16$ and $24^3\times 32\times 16$ lattice at $β=2.6$ corresponding to the lattice spacing $1/a\approx 2$GeV. Quark loop contractions which appear in Penguin diagrams are calculated by the random noise method, and the $ΔI=1/2$ matrix elements which require subtractions with the quark loop contractions are obtained with a statistical accuracy of about 10%. We confirm the chiral properties required of the $K\toπ$ matrix elements. Matching the lattice matrix elements to those in the continuum at $μ=1/a$ using the perturbative renormalization factor to one loop order, and running to the scale $μ=m_c=1.3$ GeV with the renormalization group for $N_f=3$ flavors, we calculate all the matrix elements needed for the decay amplitudes. With these matrix elements, the $ΔI=3/2$ decay amplitude shows a good agreement with experiment in the chiral limit. The $ΔI=1/2$ amplitude, on the other hand, is about 50--60% of the experimental one even after chiral extrapolation. In view ofthe insufficient enhancement of the $ΔI=1/2$ contribution, we employ the experimental values for the real parts of the decay amplitudes in our calculation of $ε'/ε$. We find that the $ΔI=3/2$ contribution is larger than the $ΔI=1/2$ contribution so that $ε'/ε$ is negative and has a magnitude of order $10^{-4}$. Possible reasons for these unsatisfactory results are discussed.

preprint2003arXiv

Continuum limit of proton decay matrix elements in quenched lattice QCD

We present a lattice QCD calculation of the parameters αand βwhich are necessary in the theoretical estimation of the proton lifetime in grand unified theories (GUTs) using chiral lagrangian approach. The simulation is carried out using the Wilson quark action at three gauge coupling constants in the quenched approximation. We obtain |α(2GeV)|=0.0091(08)(^{+10}_{-19})GeV^3 and |β(2GeV)|=0.0098(08)(^{+10}_{-20})GeV^3 in the continuum limit where the first error is statistical and the second one is due to scale setting.

preprint2003arXiv

Flavor Singlet Meson Mass in the Continuum Limit in Two-Flavor Lattice QCD

We present results for the mass of the eta-prime meson in the continuum limit for two-flavor lattice QCD, calculated on the CP-PACS computer, using a renormalization-group improved gauge action, and Sheikholeslami and Wohlert's fermion action with tadpole-improved csw. Correlation functions are measured at three values of the coupling constant beta corresponding to the lattice spacing a approx. 0.22, 0.16, 0.11 fm and for four values of the quark mass parameter kappa corresponding to mpi over mrho approx. 0.8, 0.75, 0.7 and 0.6. For each beta, kappa pair, 400-800 gauge configurations are used. The two-loop diagrams are evaluated using a noisy source method. We calculate eta-prime propagators using local sources, and find that excited state contributions are much reduced by smearing. A full analysis for the smeared propagators gives metaprime=0.960(87)+0.036-0.248 GeV, in the continuum limit, where the second error represents the systematic uncertainty coming from varying the functional form for chiral and continuum extrapolations.

preprint2003arXiv

Heavy quark expansion parameters from lattice NRQCD

We present a lattice QCD calculation of the heavy quark expansion parameters $μ_π^2$ and $μ_G^2$ for heavy-light mesons and heavy-light-light baryons. The calculation is carried out on a 20$^3\times$48 lattice at $β$ = 6.0 in the quenched approximation, using the lattice NRQCD action for heavy quarks. We obtain the parameters $μ_π^2$ and $μ_G^2$ in two different methods: a direct calculation of the matrix elements and an indirect calculation through the mass spectrum, and confirm that the both methods give consistent results. We also discuss an application to the lifetime ratios.

preprint2003arXiv

I=2 pion-pion scattering phase shift in the continuum limit calculated with two-flavor full QCD

We present a calculation of the scattering phase shift for the I=2 S-wave pion-pion system in the continuum limit with two-flavor full QCD. Calculations are made at three lattice spacings, using the finite volume method of Lüscher in the center of mass frame, and its extension to the laboratory frame.

preprint2003arXiv

Light hadron spectrum in three-flavor QCD with O(a)-improved Wilson quark action

We report on a calculation of the light hadron spectrum and quark masses in three-flavor dynamical QCD using the non-perturbatively O(a)-improved Wilson quark action and a renormalization-group improved gauge action. Simulations are carried out on a 16^3 \times 32 lattice at β=1.9, where a^{-1} \simeq 2GeV, with 6 ud quark masses corresponding to m_{pi}/m_{rho} \simeq 0.64-0.77 and 2 s quark masses close to the physical value. We observe that the inclusion of dynamical strange quark brings the lattice QCD meson spectrum to good agreement with experiment. Dynamical strange quarks also lead to a reduction of the uds quark masses by about 15%.

preprint2003arXiv

Non-perturbative calculation of $Z_V$ and $Z_A$ in domain-wall QCD on a finite box

We report on a non-perturbative evaluation of the renormalization factors for the vector and axial-vector currents, $Z_V$ and $Z_A$, in the quenched domain-wall QCD (DWQCD) with plaquette and renormalization group improved gauge actions. We take the Dirichlet boundary condition for both gauge and domain-wall fermion fields on the finite box, and introduce the flavor-chiral Ward-Takahashi identities to calculate the renormalization factors. As a test of the method, we numerically confirm the expected relation that $Z_V \simeq Z_A$ in DWQCD. Employing two different box sizes for the numerical simulations at several values of the gauge coupling constant $g^2$ and the domain-wall height $M$, we extrapolate $Z_V$ to the infinite volume to remove $a/L$ errors. We finally give the interpolation formula of $Z_V$ in the infinite volume as a function of $g^2$ and $M$.

preprint2003arXiv

Non-Perturbative Determination of $c_{\rm SW}$ in Three-flavor Dynamical QCD

We present a fully non-perturbative determination of the $O(a)$ improvement coefficient $c_{\rm SW}$ in three-flavor dynamical QCD for the RG improved as well as the plaquette gauge actions, using the Schrödinger functional scheme. Results are compared with one-loop estimates at weak gauge coupling.

preprint2003arXiv

Non-perturbative renormalization of vector and axial vector currents in quenched QCD for a renormalization group improved gauge action

Renormalization constants of vector ($Z_V$) and axial-vector ($Z_A$) currents are determined non-perturbatively in quenched QCD for an RG-improved gauge action and a tadpole-improved clover quark action using the Schrödinger functional method. Meson decay constants $f_ρ$ and $f_π$ show much better scaling when $Z_V$ and $Z_A$ estimated for infinite physical volume are used instead of $Z$-factors from tadpole-improved one-loop perturbation theory.

preprint2003arXiv

Study of finite volume effects in the non-perturbative determination of $\csw$ with the SF method in full three-flavor lattice QCD

The non-perturbative $\csw$ determined by the Schrödinger functional (SF) method with the RG-improved gauge action in dynamical $N_f=3$ QCD shows a finite volume effect when the numerical simulations are carried out at a constant lattice size $L/a$. We remove the unwanted finite volume effect by keeping physical lattice extent $L$ at a constant. The details of the method and the result obtained for non-perturbative $\csw$ with a constant $L$ are reported.

preprint2003arXiv

Two flavors of dynamical quarks on anisotropic lattices

We report on our study of two-flavor full QCD on anisotropic lattices using $O(a)$-improved Wilson quarks coupled with an RG-improved glue. The bare gauge and quark anisotropies corresponding to the renormalized anisotropy $ξ=a_s/a_t = 2$ are determined as functions of $β$ and $κ$, which covers the region of spatial lattice spacings $a_s\approx 0.28$--0.16 fm and $m_{PS}/m_V\approx 0.6$--0.9. The calibrations of the bare anisotropies are performed with the Wilson loop and the meson dispersion relation at 4 lattice cutoffs and 5--6 quark masses. Using the calibration results we calculate the meson mass spectrum and the Sommer scale $r_0$. We confirm that the values of $r_0$ calculated for the calibration using pseudo scalar and vector meson energy momentum dispersion relation coincide in the continuum limit within errors. This work serves to lay ground toward studies of heavy quark systems and thermodynamics of QCD including the extraction of the equation of state in the continuum limit using Wilson-type quark actions.

preprint2002arXiv

Charmonium Spectrum from Quenched Anisotropic Lattice QCD

We present a detailed study of the charmonium spectrum using anisotropic lattice QCD. We first derive a tree-level improved clover quark action on the anisotropic lattice for arbitrary quark mass. The heavy quark mass dependences of the improvement coefficients, i.e. the ratio of the hopping parameters $ζ=K_t/K_s$ and the clover coefficients $c_{s,t}$, are examined at the tree level. We then compute the charmonium spectrum in the quenched approximation employing $ξ= a_s/a_t = 3$ anisotropic lattices. Simulations are made with the standard anisotropic gauge action and the anisotropic clover quark action at four lattice spacings in the range $a_s$=0.07-0.2 fm. The clover coefficients $c_{s,t}$ are estimated from tree-level tadpole improvement. On the other hand, for the ratio of the hopping parameters $ζ$, we adopt both the tree-level tadpole-improved value and a non-perturbative one. We calculate the spectrum of S- and P-states and their excitations. The results largely depend on the scale input even in the continuum limit, showing a quenching effect. When the lattice spacing is determined from the $1P-1S$ splitting, the deviation from the experimental value is estimated to be $\sim$30% for the S-state hyperfine splitting and $\sim$20% for the P-state fine structure. Our results are consistent with previous results at $ξ= 2$ obtained by Chen when the lattice spacing is determined from the Sommer scale $r_0$. We also address the problem with the hyperfine splitting that different choices of the clover coefficients lead to disagreeing results in the continuum limit.

preprint2002arXiv

Dynamical fermions on anisotropic lattices

preprint2002arXiv

Exploring QCD at small sea quark masses with improved Wilson-type quarks

We explore the region of small sea quark masses below $m_{PS}/m_V=0.5$ in two-flavor QCD using a mean-field improved clover quark action and an RG-improved gauge action at $a \simeq 0.2$ fm on $12^3 \times 24$ and $16^3 \times 24$ lattices. We find that instability of the standard BiCGStab algorithm at small quark masses can be mostly removed by the BiCGStab(DS-$L$) algorithm, which employs $L$-th minimal residual polynomials with a dynamical selection of $L$. We also find singular spikes of $ΔH$ in the HMC algorithm at moderate values of $Δτ$. Nature of the spike is studied. We also study finite-size effects and chiral properties of meson masses.

preprint2002arXiv

I=2 Pion Scattering Length with the Wilson Fermion

The calculation of the I=2 pion scattering length in quenched lattice QCD is revisited. The calculation is carried out with the Wilson fermion action employing Lüscher's finite size scaling method at $β=5.9$, 6.1, and 6.3 corresponding to the range of lattice spacing $a\simeq 0.12 - 0.07$ fm. We obtain in the continuum limit $a_0m_π= -0.0410(69)$, which is consistent with the prediction of chiral perturbation theory $a_0m_π=-0.0444(10)$.

preprint2002arXiv

I=2 Pion Scattering Phase Shift with Wilson Fermions

We present results of phase shift for I=2 $S$-wave $ππ$ system with the Wilson fermions in the quenched approximation. The finite size method proposed by Lüscher is employed, and calculations are carried out at $β=5.9$ ($a^{-1}=1.934(16)$ GeV from $m_ρ$) on $24^3 \times 60$, $32^3 \times 60$, and $48^3 \times 60$ lattices.

preprint2002arXiv

I=2 Pion Scattering Phase Shift with Wilson Fermions

We present a lattice QCD calculation of the scattering phase shift for the I=2 $S$-wave two-pion system using the finite size method proposed by Lüscher. We work in the quenched approximation employing the standard plaquette action at $β=5.9$ for gluons and the Wilson fermion action for quarks. The phase shift is extracted from the energy eigenvalues of the two-pion system, which are obtained by a diagonalization of the pion 4-point function evaluated for a set of relative spatial momenta. In order to change momentum of the two-pion system, calculations are carried out on $24^3\times 60$, $32^3\times 60$, and $48^3\times 60$ lattices. The phase shift is successfully calculated over the momentum range $0 < p^2 < 0.3 {\rm GeV}^2$.

preprint2002arXiv

Light Hadron Spectrum and Quark Masses from Quenched Lattice QCD

We present details of simulations for the light hadron spectrum in quenched QCD carried out on the CP-PACS parallel computer. Simulations are made with the Wilson quark action and the plaquette gauge action on 32^3x56 - 64^3x112 lattices at four lattice spacings (a \approx 0.1-0.05 fm) and the spatial extent of 3 fm. Hadronic observables are calculated at five quark masses (m_{PS}/m_V \approx 0.75 - 0.4), assuming the u and d quarks being degenerate but treating the s quark separately. We find that the presence of quenched chiral singularities is supported from an analysis of the pseudoscalar meson data. We take m_π, m_ρand m_K (or m_ϕ) as input. After chiral and continuum extrapolations, the agreement of the calculated mass spectrum with experiment is at a 10% level. In comparison with the statistical accuracy of 1-3% and systematic errors of at most 1.7% we have achieved, this demonstrates a failure of the quenched approximation for the hadron spectrum: the meson hyperfine splitting is too small, and the octet masses and the decuplet mass splittings are both smaller than experiment. Light quark masses are calculated using two definitions: the conventional one and the one based on the axial-vector Ward identity. The two results converge toward the continuum limit, yielding m_{ud}=4.29(14)^{+0.51}_{-0.79} MeV. The s quark mass depends on the strange hadron mass chosen for input: m_s = 113.8(2.3)^{+5.8}_{-2.9} MeV from m_K and m_s = 142.3(5.8)^{+22.0}_{-0} MeV from m_ϕ, indicating again a failure of the quenched approximation. We obtain Λ_{\bar{MS}}^{(0)}= 219.5(5.4) MeV. An O(10%) deviation from experiment is observed in the pseudoscalar meson decay constants.

preprint2002arXiv

Polynomial Hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors

We present a polynomial hybrid Monte Carlo (PHMC) algorithm for lattice QCD with odd numbers of flavors of O(a)-improved Wilson quark action. The algorithm makes use of the non-Hermitian Chebyshev polynomial to approximate the inverse square root of the fermion matrix required for an odd number of flavors. The systematic error from the polynomial approximation is removed by a noisy Metropolis test for which a new method is developed. Investigating the property of our PHMC algorithm in the N_f=2 QCD case, we find that it is as efficient as the conventional HMC algorithm for a moderately large lattice size (16^3 times 48) with intermediate quark masses (m_{PS}/m_V ~ 0.7-0.8). We test our odd-flavor algorithm through extensive simulations of two-flavor QCD treated as an N_f = 1+1 system, and comparing the results with those of the established algorithms for N_f=2 QCD. These tests establish that our PHMC algorithm works on a moderately large lattice size with intermediate quark masses (16^3 times 48, m_{PS}/m_V ~ 0.7-0.8). Finally we experiment with the (2+1)-flavor QCD simulation on small lattices (4^3 times 8 and 8^3 times 16), and confirm the agreement of our results with those obtained with the R algorithm and extrapolated to a zero molecular dynamics step size.

preprint2001arXiv

Charmonium spectrum from quenched QCD on anisotropic lattices

We present our final results of the charmonium spectrum in quenched QCD on anisotropic lattices. Simulations are made with the plaquette gauge action and a tadpole improved clover quark action employing $ξ= a_s/a_t = 3$. We calculate the spectrum of S- and P-states and their excitation, and study the scaling behavior of mass splittings. Comparison is made with the experiment and previous lattice results. The issue of hyperfine splitting for different choices of the clover coefficients obtained by Klassen is discussed.

preprint2001arXiv

Chiral property of domain-wall fermion from eigenvalues of 4D Wilson-Dirac Operator

We investigate a chiral property of the domain-wall fermion (DWF) system using the four-dimensional hermitian Wilson-Dirac operator $H_W$. A formula expressing the Ward-Takahashi identity quark mass $m_{5q}$ with eigenvalues of this operator is derived, which well explains the $N_5$ dependence of $m_{5q}$ observed in previous numerical simulations. We further discuss the chiral property of DWF in the large volume in terms of the spectra of $H_W$.

preprint2001arXiv

Equation of state for pure SU(3) gauge theory on anisotropic lattices

We present results for the equation of state for pure SU(3) gauge theory obtained on anisotropic lattices with the anisotropy $ξ\equiv a_s/a_t = 2$. The pressure and energy density are calculated on $N_t / ξ= 4, 5$ and 6 lattices with the integral method. They are found to satisfy the leading $1/N_t^2$ scaling from our coarsest lattice $N_t/ξ=4$. This enables us to carry out well controlled continuum extrapolations. We find that the pressure and energy density agree with those obtained using the isotropic plaquette action, but have smaller and more reliable errors.

preprint2001arXiv

Equation of state in finite-temperature QCD with two flavors of improved Wilson quarks

We present results of a first study of equation of state in finite-temperature QCD with two flavors of Wilson-type quarks. Simulations are made on lattices with temporal size $N_t=4$ and 6, using an RG-improved action for the gluon sector and a meanfield-improved clover action for the quark sector. The lines of constant physics corresponding to fixed values of the ratio $m_{\rm PS}/m_{\rm V}$ of the pseudo-scalar to vector meson masses at zero temperature are determined, and the beta functions which describe the renormalization-group flow along these lines are calculated. Using these results, the energy density and the pressure are calculated as functions of temperature along the lines of constant physics in the range $m_{\rm PS}/m_{\rm V} = 0.65$--0.95. The quark mass dependence in the equation of state is found to be small for $m_{\rm PS}/m_{\rm V} \simlt 0.8$. Comparison of results for $N_t=4$ and $N_t=6$ lattices show significant scaling violation present in the $N_t=4$ results. At high temperatures the results for $N_t=6$ are quite close to the continuum Stefan-Boltzmann limit, suggesting the possibility of a precise continuum extrapolation of thermodynamic quantities from simulations at $N_t\simgt 6$.

preprint2001arXiv

Hybrid Quarkonia with Dynamical Sea Quarks

We present a dynamical lattice calculation with 2 flavours for bottomonium states with an additional gluonic excitation. Using improved actions for the quarks and gauge fields at a lattice spacing of $a \approx 0.1$ fm, we find 10.977(61)(62) GeV for the energy of the lowest lying $b\bar bg$-hybrid, where the first error is statistical and the second denotes the systematic uncertainty due to the determination of scale. In a parallel quenched simulation we demonstrate explicitly that vacuum polarisation effects are less than 10% of the splitting with the ground state.

preprint2001arXiv

I=2 Pion Scattering Length and Phase Shift with Wilson Fermions

We present preliminary results of scattering length and phase shift for I=2 S-wave $ππ$ system with the Wilson fermions in the quenched approximation. The finite size method presented by Lüscher is employed, and calculations are carried out at $β=5.9$ on a $24^3\times 60$ and $32^3\times 60$ lattice.

preprint2001arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 GeV)= 0.5746(61)(191)$, where the first error is statistical and the second error represents an estimate of scaling violation and ${\cal O}(α^2)$ errors in the renormalization factor added in quadrature, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is consistent, albeit somewhat small, with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for light quark masses are also reported.

preprint2001arXiv

Light Hadron Spectroscopy with Two Flavors of Dynamical Quarks on the Lattice

We present results of a numerical calculation of lattice QCD with two degenerate flavors of dynamical quarks, identified with up and down quarks, and with a strange quark treated in the quenched approximation. The lattice action and simulation parameters are chosen with a view to carrying out an extrapolation to the continuum limit as well as chiral extrapolations. Gauge configurations are generated with a renormalization-group improved gauge action and a mean field improved clover quark action at three values of $β$ and four sea quark masses. The sizes of lattice are chosen so that the physical spatial size is kept constant. Hadron masses, light quark masses and meson decay constants are measured at five valence quark masses. We also carry out complementary quenched simulations with the same improved actions. The quenched spectrum from this analysis agrees well in the continuum limit with the one of our earlier work using the standard action. We find the two-flavor full QCD meson masses in the continuum limit to be much closer to experimental meson masses than those from quenched QCD. We take these results as manifestations of sea quark effects in two-flavor full QCD. For baryon masses full QCD values for strange baryons are in agreement with experiment, while they differ increasingly with decreasing strange quark content, resulting in a nucleon mass higher than experiment. The pattern suggests finite size effects as a possible origin for this deviation. For light quark masses in the continuum limit we obtain values which are reduced by about 25% compared to the values in quenched QCD. We also present results for decay constants where large scaling violations obstruct a continuum extrapolation. Need for a non-perturbative estimate of renormalization factors is discussed.

preprint2001arXiv

Maximum entropy analysis of hadron spectral functions and excited states in quenched lattice QCD

Employing the maximum entropy method we extract the spectral functions from meson correlators at four lattice spacings in quenched QCD with the Wilson quark action. We confirm that the masses and decay constants, obtained from the position and the area of peaks, agree well with the results from the conventional exponential fit. For the first excited state, we obtain $m_{π_1} = 660(590)$ MeV, $m_{ρ_1} = 1540(570)$ MeV, and $f_{ρ_1} = 0.085(36)$ in the continuum limit.

preprint2001arXiv

Non-perturbative renormalization in domain-wall QCD by Schrödinger functional scheme

We formulate and numerically test the Schrödinger functional scheme for domain-wall QCD. We then apply it to a non-perturbative calculation of the renormalization factors for vector and axial-vector currents in quenched domain-wall QCD with plaquette and renormalization group improved gauge actions at $a^{-1}\simeq 2$ GeV.

preprint2001arXiv

Spectral function and excited states in lattice QCD with maximum entropy method

We apply the maximum entropy method to extract the spectral functions for pseudoscalar and vector mesons from hadron correlators previously calculated at four different lattice spacings in quenched QCD with the Wilson quark action. We determine masses and decay constants for the ground and excited states of the pseudoscalar and vector channels from position and area of peaks in the spectral functions. We obtain the results, $m_{π_1} = 660(590)$ MeV and $m_{ρ_1} = 1540(570)$ MeV for the masses of the first excited state masses, in the continuum limit of quenched QCD. We also find unphysical states which have infinite mass in the continuum limit, and argue that they are bound states of two doublers of the Wilson quark action. If the interpretation is correct, this is the first time that the state of doublers is identified in lattice QCD numerical simulations.

preprint2001arXiv

Topological Susceptibility in Lattice QCD with Two Flavors of Dynamical Quarks

We present a study of the topological susceptibility in lattice QCD with two degenerate flavors of dynamical quarks. The topological charge is measured on gauge configurations generated with a renormalization group improved gauge action and a mean field improved clover quark action at three values of $β=6/g^2$, corresponding to lattice spacings of $a \approx 0.22$, 0.16 and 0.11 fm, with four sea quark masses at each $β$. The study is supplemented by simulations of pure SU(3) gauge theory with the same gauge action at 5 values of $β$ with lattice spacings 0.09 fm$\simlt a \simlt$0.27 fm. We employ a field theoretic definition of the topological charge together with cooling. For the topological susceptibility in the continuum limit of pure SU(3) gauge theory we obtain $χ_t^{1/4} = 197^{+13}_{-16}$ MeV where the error shows statistical and systematic ones added in quadrature. In full QCD $χ_t$ at heavy sea quark masses is consistent with that of pure SU(3) gauge theory. A decrease of $χ_t$ toward light quark masses, as predicted by the anomalous Ward-Takahashi identity for U(1) chiral symmetry, becomes clearer for smaller lattice spacings. The cross-over in the behavior of $χ_t$ from heavy to light sea quark masses is discussed.

preprint2000arXiv

Calculation of $K\to π$ matrix elements in quenched domain-wall QCD

We explore the possibility for an evaluation of non-leptonic $ΔS=1$ $K$ decay amplitudes through the calculation of $K\toπ$ matrix elements using domain-wall QCD. The relation between the physical $K\toππ$ matrix elements and $K\toπ$ matrix elements deduced from chiral perturbation theory is recapitulated. Quenched numerical simulations are performed on an $16^3\times 32\times 16$ lattice at a lattice spacing $a^{-1}\approx 2$GeV for both the standard plaquette gauge action and a renormalization-group improved gauge action, and reasonable signals for $K\toπ$ matrix elements are obtained. Preliminary results are reported on the $K\toππ$ matrix elements, and results from two actions are compared.

preprint2000arXiv

Decay constants of B and D mesons from improved relativistic lattice QCD with two flavours of sea quarks

We present a calculation of the B and D meson decay constants in lattice QCD with two (Nf=2) flavours of light dynamical quarks, using an O(a)-improved Wilson action for both light and heavy quarks and a renormalization-group improved gauge action. Simulations are made at three values of lattice spacing a=0.22, 0.16, 0.11 fm and four values of sea quark mass in the range m_PS/m_V \~= 0.8-0.6. Our estimate for the continuum values of the decay constants are fBd = 208(10)(11) MeV, fBs = 250(10)(13)(^{+8}_{-0}) MeV, fDd = 225(14)(14) MeV, fDs = 267(13)(17)(^{+10}_{-0}) MeV for Nf=2 where the statistical and systematic errors are separately listed, and the third error for fBs and fDs show uncertainty of determination of strange quark mass. We also carry out a set of quenched simulations using the same action to make a direct examination of sea quark effects. Taking the ratio of results for Nf=2 and Nf=0, we obtain fb^{Nf=2}/fb^{Nf=0} = 1.11(6), fbs^{Nf=2}/fbs^{Nf=0} = 1.14(5), fd^{Nf=2}/\fd^{Nf=0} = 1.03(6), fds^{Nf=2}/\fds^{Nf=0} = 1.07(5). They show a 10-15% increase in the Nf=2 results over those of Nf=0 for the B meson decay constants, while evidence for such a trend is statistically less clear for the D meson decay constants.

preprint2000arXiv

Dynamical Quark Effects on Light Quark Masses

We present results for light quark masses from a systematic lattice QCD study with two degenerate flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 and the lattice spacing $a=0.22$--0.11 fm. In the continuum limit we find $m_{ud}^{\bar{MS}}(2 {\rm GeV})=3.44^{+0.14}_{-0.22}$ MeV using the $π$ and $ρ$ meson masses as physical input, and $m_s^{\bar{MS}}(2 {\rm GeV})=88^{+4}_{-6}$ MeV or $90^{+5}_{-11}$ MeV with the $K$ or $ϕ$ meson mass as additional input. The quoted errors represent statistical and systematic combined, the latter including those from continuum and chiral extrapolations, and from renormalization factors. Compared to quenched results, two flavors of dynamical quarks reduce quark masses by about 25%.

preprint2000arXiv

Eigenvalues of the hermitian Wilson-Dirac operator and chiral properties of the domain-wall fermion

Chiral properties of QCD formulated with the domain-wall fermion (DWQCD) are studied using the anomalous quark mass m_{5q} and the spectrum of the 4-dimensional Wilson-Dirac operator. Numerical simulations are made with the standard plaquette gauge action and a renormalization-group improved gauge action. Results are reported on the density of zero eigenvalue obtained with the accumulation method, and a comparison is made with the results for m_{5q}.

preprint2000arXiv

Full QCD Light Hadron Spectrum and Quark Masses: Final Results from CP-PACS

We present the final results of the CP-PACS calculation of the light hadron spectrum and quark masses with two flavors of dynamical quarks. Simulations are made with a renormalization-group improved gauge action and a mean-field improved clover quark action for sea quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 and the lattice spacing $a=0.22$--0.11 fm. For the meson spectrum in the continuum limit a clearly improved agreement with experiment is observed compared to the quenched case, demonstrating the importance of sea quark effects. For light quark masses we obtain $m_{ud}^{\bar{MS}}(2GeV)=3.44^{+0.14}_{-0.22}$ MeV and $m_s^{\bar{MS}}(2GeV)=88^{+4}_{-6}$ MeV ($K$-input) and $m_s^{\bar{MS}}(2GeV)=90^{+5}_{-11}$ MeV ($ϕ$-input), which are reduced by about 25% compared to the values in quenched QCD.

preprint2000arXiv

Kaon B parameter from quenched domain-wall QCD

We report on a calculation of $B_K$ with domain wall fermion action in quenched QCD. Simulations are made with a renormalization group improved gauge action at $β=2.6$ and 2.9 corresponding to $a^{-1}\approx 2$GeV and 3GeV. Effects due to finite fifth dimensional size $N_5$ and finite spatial size $N_σ$ are examined in detail. Matching to the continuum operator is made perturbatively at one loop order. We obtain $B_K(μ= 2 {GeV})= 0.5746(61)$, where the error is statistical only, as an estimate of the continuum value in the $\msbar$ scheme with naive dimensional regularization. This value is smaller but consistent with $B_K(μ= 2 {GeV})= 0.628(42)$ obtained by the JLQCD Collaboration using the Kogut-Susskind quark action. Results for strange quark mass are also reported.

preprint2000arXiv

Perturbative Renormalization Factors of Baryon Number Violating Operators for Improved Quark and Gauge Actions in Lattice QCD

We calculate one-loop renormalization factors of three-quark operators, which appear in the low energy effective Lagrangian of the nucleon decay, for $O(a)$-improved quark action and gauge action including six-link loops. This calculation is required to predict the hadronic nucleon decay matrix elements in the continuum regularization scheme using lattice QCD. We present detailed numerical results of the one-loop coefficients for general values of the clover coefficients employing the several improved gauge actions in the Symanzik approach and in the Wilson's renormalization group approach. The magnitudes of the one-loop coefficients for the improved gauge actions show sizable reduction compared to those for the plaquette action.

preprint2000arXiv

Phase structure and critical temperature of two-flavor QCD with a renormalization group improved gauge action and clover improved Wilson quark action

We study the finite-temperature phase structure and the transition temperature of QCD with two flavors of dynamical quarks on a lattice with the temporal size $N_t=4$, using a renormalization group improved gauge action and the Wilson quark action improved by the clover term. The region of a parity-broken phase is identified, and the finite-temperature transition line is located on a two-dimensional parameter space of the coupling ($β=6/g^2$) and hopping parameter $K$. Near the chiral transition point, defined as the crossing point of the critical line of the vanishing pion mass and the line of finite-temperature transition, the system exhibits behavior well described by the scaling exponents of the three-dimensional O(4) spin model. This indicates a second-order chiral transition in the continuum limit. The transition temperature in the chiral limit is estimated to be $T_c = 171(4)$ MeV.

preprint2000arXiv

Quenched charmonium spectrum on anisotropic lattices

We present the results of quenched charmonium spectrum for S- and P-states, obtained by a relativistic heavy quark method on anisotropic lattices. Simulations are carried out using the standard plaquette gauge action and a meanfield-improved clover quark action at $a_t^{-1} = 3$--6 GeV with the renormalized anisotropy fixed to $ξ\equiv a_s/a_t =3$. We study the scaling of our fine and hyperfine mass splittings, and compare with previous results.

preprint2000arXiv

Sea Quark Effects on Quarkonia

We study the effects of two dynamical sea quarks on the spectrum of heavy quarkonia. Within the non-relativistic approach to Lattice QCD we found sizeable changes to the hyperfine splitting, but we could not observe any changes for the fine structure. We also investigated the scaling behaviour of our results for several different lattice spacings.

preprint1999arXiv

$B \to πl\barν$ Form Factors with NRQCD Heavy Quark and Clover Light Quark Actions

We report results on semileptonic $B\toπl\barν$ decay form factors near $q^2_{\rm max}$ using NRQCD heavy quark and clover light quark actions and currents improved through $O(αa)$. An inconsistency with the soft pion relation $f^0(q^2_{\rm max})= f_B/f_π$ found in a previous work is confirmed, and a possible solution with nonperturbative renormalization is discussed. We find that $f^+(q^2)$ is well described by the $B^*$ pole near $q^2_{\rm max}$, and its $1/M_B$ scaling is also consistent with the prediction of the pole dominance model.

preprint1999arXiv

Comparative Study of full QCD Hadron Spectrum and Static Quark Potential with Improved Actions

We investigate effects of action improvement on the light hadron spectrum and the static quark potential in two-flavor QCD for $a^{-1} \approx 1$ GeV and $m_{PS}/m_V = 0.7-0.9$. We compare a renormalization group improved action with the plaquette action for gluons, and the SW-clover action with the Wilson action for quarks. We find a significant improvement in the hadron spectrum by improving the quark action, while the gluon improvement is crucial for a rotationally invariant static potential. We also explore the region of light quark masses corresponding to $m_{PS}/m_V \geq 0.4$ on a 2.7 fm lattice using the improved gauge and quark action. A flattening of the potential is not observed up to 2 fm.

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with renormalization group improved action

A lattice study of the equation of state for pure SU(3) gauge theory using a renormalization-group (RG) improved action is presented. The energy density and pressure are calculated on a $16^3\times 4$ and a $32^3\times 8$ lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.

preprint1999arXiv

Equation of state for pure SU(3) gauge theory with RG improved action

We present results for the equation of state for pure SU(3) gauge theory obtained with a renormalization-group (RG) improved action. The energy density and pressure are calculated on a $16^3\times 4$ and a $32^3\times 8$ lattice employing the integral method. Extrapolating the results to the continuum limit, we find the energy density and pressure to be in good agreement with those obtained with the standard plaquette action within the error of 3-4%.

preprint1999arXiv

Equation of state in finite-temperature QCD with improved Wilson quarks

We study finite-temperature phase transition and equation of state for two-flavor QCD at $N_t=4$ using an RG-improved gauge action and a meanfield-improved clover quark action. The pressure is computed using the integral method. The O(4) scaling of chiral order parameter is also examined.

preprint1999arXiv

Eta meson mass and topology in QCD with two light flavors

We present results for the mass of the flavor singlet meson calculated on two-flavor full QCD configurations generated by the CP-PACS full QCD project. We also investigate topological charge fluctuations and their dependence on the sea quark mass.

preprint1999arXiv

Heavy Quarkonia from Anisotropic and Isotropic Lattices

We report on recent results for the spectrum of heavy quarkonia. Using coarse and anisotropic lattices we achieved an unprecedented control over statistical and systematic errors for higher excited states such as exotic hybrid states. In a parallel study on isotropic lattices we also investigate the effect of two dynamical flavours on the spin structure of charmonium and bottomonium for several symmetric lattices.

preprint1999arXiv

Heavy-light decay constants from clover heavy quark action in QCD with two flavors of dynamical quarks

We present results on an analysis of the decay constants f_B and f_Bs with two flavours of sea quark. The calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge coupling, with m_pi/m_rho ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with heavy quarks whose mass is in the range of the b-quark. A detailed comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.

preprint1999arXiv

Heavy-light decay constants from relativistic Nf=2,0 QCD

We present results on an analysis of the decay constants fB and fBs with two flavours of sea quark. The calculation has been carried out on 3 different bare gauge couplings and 4 sea quark masses at each gauge coupling, with m_pi/m_rho ranging from 0.8 to 0.6. We employ the Fermilab formalism to perform calculations with heavy quarks whose mass is in the range of the b-quark. A comparison with a quenched calculation using the same action is made to elucidate the effects due to the sea quarks.

preprint1999arXiv

Heavy-light spectrum and decay constant from NRQCD with two flavors of dynamical quarks

We report on a study of B mesons on N_f = 2 full QCD configurations using an RG-improved gauge action, NRQCD heavy quark action and tadpole-improved clover light quark action. Results on the heavy-light spectrum and the decay constants from 16^3x32 lattices at a^{-1} ~ 1.5 GeV are presented, and compared with quenched results obtained with the same action combination at matching lattice spacings.

preprint1999arXiv

Hybrid Quarkonia on Asymmetric Lattices

We report on a study of heavy hybrid states using the NRQCD approach on coarse and asymmetric lattices, where we discard vacuum polarisation effects and neglect all spin-correction terms. We find a clear hybrid signal on all our lattices ($a_s= 0.15 ... 0.47$ fm). We have studied in detail the lattice spacing artefacts, finite volume effects and mass dependence. Within the above approximations we predict the hybrid excitation in Charmonium to be 1.323(13) GeV above its ground state. The bottomonium hybrid was found to be 1.542(8) GeV above its ground state.

preprint1999arXiv

I=2 Pion Scattering Length with Wilson Fermions

We present results for I=2 pion scattering length with the Wilson fermions in the quenched approximation. The finite size method presented by Lüscher is employed, and calculations are carried out at $β=5.9$, 6.1, and 6.3. In the continuum limit, we obtain a result in reasonable agreement with the experimental value.

preprint1999arXiv

Lattice QCD with Domain-Wall Fermions

We study the quenched lattice QCD using domain-wall fermions at $β=6.0$. Behaviors of both pion mass and the explicit breaking term in the axial Ward-Takahashi identity support the existence of the chiral zero modes. We observe a good agreement between the pion decay constants $f_π$ from both the conserved axial current and the local current perturbatively renormalized at 1-loop. Finally the possible existence of the parity broken phase is also examined in this model.

preprint1999arXiv

Light hadron spectrum and quark masses in QCD with two flavors of dynamical quarks

We present updated results of the CP-PACS calculation of the light hadron spectrum in $N_{\rm f}=2$ full QCD. Simulations are made with an RG-improved gauge action and a tadpole-improved clover quark action for sea quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 and the lattice spacing $a=0.22$--0.09 fm. A comparison of the full QCD spectrum with new quenched results, obtained with the same improved action, shows clearly the existence of sea quark effects in vector meson masses. Results for light quark masses in $N_{\rm f}=2$ QCD are also presented.

preprint1999arXiv

Non-perturbative determination of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action

We report results of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action, employing the Reguralization Independent scheme (RI) of Martinelli et al. to non-perturbatively evaluate the renormalization factor relating the bare quark mass on the lattice to that in the continuum. Calculations are carried out at β=6.0, 6.2, and 6.4, from which we find $m^{\bar{MS}}_{ud} (2 GeV)= 4.23(29) MeV$ for the average up and down quark mass and, with the $ϕ$ meson mass as input, $m^{\bar{MS}}_{s} (2 GeV)= 129(12) MeV$ for the strange mass in the continuum limit. These values are about 20% larger than those obtained with the one-loop perturbative renormalization factor.

preprint1999arXiv

Nucleon decay matrix elements with the Wilson quark action: an update

We present preliminary results of a new lattice computation of hadronic matrix elements of baryon number violating operators which appear in the low-energy effective Lagrangian of (SUSY-)Grand Unified Theories. The contribution of irrelevant form factor which has caused an underestimate of the matrix elements in previous studies is subtracted in this calculation. Our results are 2$\sim$4 times larger than the most conservative values often employed in phenomenological analyses of nucleon decay with specific GUT models.

preprint1999arXiv

Pion decay constant in quenched QCD with Kogut-Susskind quarks

We present a non-perturbative calculation for the pion decay constant with quenched Kogut-Susskind quarks. Numerical simulations are carried out at $β= 6.0$ and 6.2 with various operators extending over all flavors. The renormalization correction is applied for each flavor by computing non-perturbative renormalization constants, and it is compared with a perturbative calculation. We also study the behavior of $f_π$ in the continuum limits for both non-perturbative and perturbative calculations. The results in the continuum limit is also discussed.

preprint1999arXiv

Quenched Light Hadron Spectrum

We present results of a large-scale simulation for the flavor non-singlet light hadron spectrum in quenched lattice QCD with the Wilson quark action. Hadron masses are calculated at four values of lattice spacing in the range $a \approx$ 0.1 - 0.05 fm on lattices with a physical extent of 3 fm at five quark masses corresponding to $m_π/m_ρ\approx 0.75$ - 0.4. The calculated spectrum in the continuum limit shows a systematic deviation from experiment, though the magnitude of deviation is contained within 11%. Results for decay constants and light quark masses are also reported.

preprint1999arXiv

Quenched QCD with domain-wall fermions on coarse lattices

We investigate the existence of chiral zero modes at a^{-1} \simeq 1 GeV in quenched domain-wall QCD. Simulations are carried out for the plaquette and an RG-improved gauge actions on a 12^3x24xN_s lattice with N_s=10-50. We find that the pion mass in the chiral limit remains non-vanishing as N_s\to\infty for both gauge actions. Possible origins of this non-vanishing pion mass are discussed.

preprint1998arXiv

Full QCD light hadron spectrum from the CP-PACS

We report on an on-going two-flavor full QCD study on CP-PACS using an RG-improved gauge action and a tadpole-improved SW quark action. Runs are made for three lattice spacings $a^{-1}\approx 0.9$, 1.3, and 2.5 GeV on $12^3\times24$, $16^3\times32$, and $24^3\times48$ lattices. Four sea quark masses having $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.6 are simulated, for each of which hadron masses are evaluated for valence quark masses corresponding to $m_{\rm PS}/m_{\rm V} \approx 0.8$--0.5. Results for hadron and light quark masses are presented and compared with those obtained in quenched QCD.

preprint1998arXiv

Non-perturbative renormalization factors of bilinear quark operators for Kogut-Susskind fermions and light quark masses in quenched QCD

Light quark masses are computed for Kogut-Susskind fermions by evaluating non-perturbatively the renormalization factor for bilinear quark operators. Calculations are carried out in the quenched approximation at β=6.0, 6.2, and 6.4. For the average up and down quark mass we find $m_{\bar MS}(2 GeV)= 4.15(27) MeV$ in the continuum limit, which is significantly larger than $3.51(20) MeV$ ($q^*=1/a$) or $3.40(21) MeV$ ($q^*=π/a$) obtained with the one-loop perturbative renormalization factor.

preprint1998arXiv

Quenched Light Hadron Spectrum with the Wilson Quark Action: Final Results from CP-PACS

We report the final results of the CP-PACS calculation for the quenched light hadron spectrum with the Wilson quark action. Our data support the presence of quenched chiral singularities, and this motivates us to use mass formulae based on quenched chiral perturbation theory in order to extrapolate hadron masses to the physical point. Hadron masses and decay constants in the continuum limit show unambiguous systematic deviations from experiment. We also report the results for light quark masses.

preprint1998arXiv

The static quark potential in full QCD

We report results on the static quark potential in two-flavor full QCD. The calculation is performed for three values of lattice spacing $a^{-1}\approx 0.9, 1.3$ and 2.5 GeV on $12^3{\times}24, 16^3{\times}32$ and $24^3{\times}48$ lattices respectively, at sea quark masses corresponding to $m_π/m_ρ\approx 0.8-0.6$. An RG-improved gauge action and a tadpole-improved SW clover quark action are employed. We discuss scaling of $m_ρ/\sqrtσ$ and effects of dynamical quarks on the potential.

preprint1997arXiv

B meson decay constant with the Wilson and Clover heavy quark actions

We present results of our quenched study of the B meson decay constant obtained with a parallel set of simulations with the Wilson and Clover actions at β=5.9, 6.1 and 6.3. Systematic errors associated with the large b-quark mass are analyzed within the Fermilab non-relativistic formalism. As our best estimate in the continuum limit we obtain $f_B=163 \pm 16 MeV$ and $f_{B_s}=175 \pm 18 MeV$ with the Clover action.

preprint1997arXiv

CP-PACS Result for the Quenched Light Hadron Spectrum

The quenched hadron spectrum in the continuum obtained with the Wilson quark action in recent simulations on the CP-PACS is presented. Results for the light quark masses and the QCD scale parameter are reported.

preprint1997arXiv

Hadron spectroscopy and static quark potential in full QCD: A comparison of improved actions on the CP-PACS

We present first results from a full QCD calculation on the CP-PACS, comparing various actions at $a^{-1} \sim 1 GeV$ and $m_π/m_ρ\approx 0.7$--0.9. We use the plaquette and a renormalization group improved action for the gluons, and the Wilson and the SW-Clover action for quarks. We find that significant improvements in the hadron spectrum results from improving the quarks, while the gluon improvement is required for a rotationally invariant static potential. An ongoing effort towards exploring the chiral limit in full QCD is described.

preprint1997arXiv

Kaon B parameter from quenched Lattice QCD

We present results of a large-scale simulation for the Kaon B parameter $B_K$ in quenched lattice QCD with the Kogut-Susskind quark action. Calculating $B_K$ at 1% statistical accuracy for seven values of lattice spacing in the range $a\approx 0.24-0.04$ fm on lattices up to $56^3\times 96$, we verify a quadratic $a$ dependence of $B_K$ theoretically predicted. Strong indications are found that, with our level of accuracy, $α_{\bar{MS}}(1/a)^2$ terms arising from our one-loop matching procedure have to be included in the continuum extrapolation. We present $B_K$(NDR, 2 GeV)=0.628(42) as our final value, as obtained by a fit including the $α_{\bar{MS}}(1/a)^2$ term.

preprint1997arXiv

Scaling Study of the Two-Flavor Chiral Phase Transition with the Kogut-Susskind Quark Action in Lattice QCD

We report on a study of two-flavor finite-temperature chiral phase transition employing the Kogut-Susskind quark action and the plaquette gluon action in lattice QCD for a lattice with $N_t=4$ temporal size. Hybrid R simulations of $10^4$ trajectories are made at quark masses of $m_q=0.075, 0.0375, 0.02, 0.01$ in lattice units for the spatial sizes $8^3, 12^3$ and $16^3$. The spatial size dependence of various susceptibilities confirm the previous conclusion of the absence of a phase transition down to $m_q=0.02$. At $m_q=0.01$ an increase of susceptibilities is observed up to the largest volume $16^3$ explored in the present work. We argue, however, that this increase is likely to be due to an artifact of too small a lattice size and it cannot be taken to be the evidence for a first-order transition. Analysis of critical exponents estimated from the quark mass dependence of susceptibilities shows that they satisfy hyperscaling consistent with a second-order transition located at $m_q=0$. The exponents obtained from larger lattice, however, deviate significantly from both those of O(2), which is the exact symmetry group of the Kogut-Susskind action at finite lattice spacing, and those of O(4) expected from an effective sigma model analysis in the continuum limit.

preprint1996arXiv

Tensor Charge of the Nucleon on the Lattice

Tensor charge of the nucleon, which will be measured in Drell-Yan processes in polarized proton-proton collisions at RHIC, are studied in quenched lattice QCD simulation. On the 16$^3\times$ 20 lattice with $β=5.7$, connected parts of the tensor charge are determined with small statistical error, while the disconnected parts are found to be small with relatively large error bars. Flavor-singlet tensor charge ($δΣ= δu + δd + δs$) is not suppressed as opposed to the flavor singlet axial charge ($ΔΣ= Δu + Δd + Δs$).

S. Aoki

What is connected

Connect this record

See the researcher in context

Building this map preview

229 published item(s)

Measurement of muon neutrino induced charged current interactions without charged pions in the final state using a new T2K off-axis near detector WAGASCI-BabyMIND

Results from the T2K experiment on neutrino mixing including a new far detector $μ$-like sample

Updated constraints on sterile neutrino mixing in the OPERA experiment using a new $ν_e$ identification method

Measurements of protons and charged pions emitted from $ν_μ$ charged-current interactions on iron at a mean neutrino energy of 1.49$\,$GeV using a nuclear emulsion detector

Milliarcsecond astrometry for the Galilean moons using stellar occultations

Role of axial U(1) anomaly in chiral susceptibility of QCD at high temperature

Scintillator ageing of the T2K near detectors from 2010 to 2021

What is chiral susceptibility probing?

Constraint on the Matter-Antimatter Symmetry-Violating Phase in Neutrino Oscillations

Sensitivity of the SHiP experiment to dark photons decaying to a pair of charged particles

First combined measurement of the muon neutrino and antineutrino charged-current cross section without pions in the final state at T2K

First observation of a tau neutrino charged current interaction with charm production in the Opera experiment

FLAG Review 2019

Search for heavy neutrinos with the T2K near detector ND280

Simultaneous measurement of the muon neutrino charged-current cross section on oxygen and carbon without pions in the final state at T2K

SND@LHC

T2K measurements of muon neutrino and antineutrino disappearance using $3.13\times 10^{21}$ protons on target

First Measurement of the Charged Current $\overlineν_μ$ Double Differential Cross Section on a Water Target without Pions in the final state

Measurement of neutrino and antineutrino neutral-current quasielastic-like interactions on oxygen by detecting nuclear de-excitation $γ$-rays

Measurement of the muon neutrino charged-current single $π^+$ production on hydrocarbon using the T2K off-axis near detector ND280

Search for Electron Antineutrino Appearance in a Long-baseline Muon Antineutrino Beam

Chiral behavior of light meson form factors in 2+1 flavor QCD with exact chiral symmetry

Determination of the muon charge sign with the dipolar spectrometers of the OPERA experiment

Measurement of coherent $π^{+}$ production in low energy neutrino-Carbon scattering

Measurement of double-differential muon neutrino charged-current interactions on C$_8$H$_8$ without pions in the final state using the T2K off-axis beam

Measurement of Muon Antineutrino Oscillations with an Accelerator-Produced Off-Axis Beam

Proposal for an Extended Run of T2K to $20\times10^{21}$ POT

Review of lattice results concerning low-energy particle physics

Sensitivity of the T2K accelerator-based neutrino experiment with an Extended run to $20\times10^{21}$ POT

$η^\prime$ meson mass from topological charge density correlator in QCD

A facility to Search for Hidden Particles (SHiP) at the CERN SPS

Discovery of tau neutrino appearance in the CNGS neutrino beam with the OPERA experiment

Light meson electromagnetic form factors from three-flavor lattice QCD with exact chiral symmetry

Measurement of the $ν_μ$ CCQE cross section on carbon with the ND280 detector at T2K

Measurement of the $ν_μ$ charged current quasi-elastic cross-section on carbon with the T2K on-axis neutrino beam

Measurement of the Electron Neutrino Charged-current Interaction Rate on Water with the T2K ND280 pi-zero Detector

Measurement of the muon neutrino inclusive charged-current cross section in the energy range of 1-3 GeV with the T2K INGRID detector

Measurements of neutrino oscillation in appearance and disappearance channels by the T2K experiment with 6.6E20 protons on target

Neutrino Oscillation Physics Potential of the T2K Experiment

Computation of the electromagnetic pion form factor from lattice QCD in the epsilon regime

Evidence for $ν_μ\to ν_τ$ appearance in the CNGS neutrino beam with the OPERA experiment

Histograms in heavy-quark QCD at finite temperature and density

Measurement of TeV atmospheric muon charge ratio with the full OPERA data

Measurement of the inclusive $ν_μ$ charged current cross section on iron and hydrocarbon in the T2K on-axis neutrino beam

Measurement of the Inclusive Electron Neutrino Charged Current Cross Section on Carbon with the T2K Near Detector

Measurement of the intrinsic electron neutrino component in the T2K neutrino beam with the ND280 detector

Measurement of the muon beam direction and muon flux for the T2K neutrino experiment

Measurement of the neutrino-oxygen neutral-current interaction cross section by observing nuclear deexcitation $γ$ rays

Observation of nu_tau appearance in the CNGS beam with the OPERA experiment

Precise Measurement of the Neutrino Mixing Parameter θ_{23} from Muon Neutrino Disappearance in an Off-axis Beam

Procedure for short-lived particle detection in the OPERA experiment and its application to charm decays

Search for short baseline $ν_e$ disappearance with the T2K near detector

Spin-Orbit Force from Lattice QCD

Strong coupling constant from vacuum polarization functions in three-flavor lattice QCD with dynamical overlap fermions

Topology density correlator on dynamical domain-wall ensembles with nearly frozen topological charge

Charmed baryons at the physical point in 2+1 flavor lattice QCD

Large-scale simulations with chiral symmetry

Measurement of Neutrino Oscillation Parameters from Muon Neutrino Disappearance with an Off-axis Beam

Measurement of the Inclusive NuMu Charged Current Cross Section on Carbon in the Near Detector of the T2K Experiment

New results on $ν_μ\to ν_τ$ appearance with the OPERA experiment in the CNGS beam

Nucleon strange quark content from N_f=2+1 lattice QCD with exact chiral symmetry

Overlap/Domain-wall reweighting

Phase shifts in I=2 ππ-scattering from two lattice approaches

Search for $ν_μ\rightarrow ν_e$ oscillations with the OPERA experiment in the CNGS beam

Short Distance Repulsion Among Baryons

The T2K Neutrino Flux Prediction

1+1+1 flavor QCD + QED simulation at the physical point

Application of fixed scale approach to static quark free energies in quenched and 2+1 flavor lattice QCD with improved Wilson quark action

Chiral behavior of kaon semileptonic form factors in lattice QCD with exact chiral symmetry

Equation of state in 2+1 flavor QCD with improved Wilson quarks by the fixed scale approach

Finite density QCD phase transition in the heavy quark region

First Muon-Neutrino Disappearance Study with an Off-Axis Beam

Measurement of the neutrino velocity with the OPERA detector in the CNGS beam

Phase structure of finite temperature QCD in the heavy quark region