Source author record

M. Fukugita

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

97works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

97 published item(s)

preprint2009arXiv

Luminosity Functions of Type Ia Supernovae and their Host Galaxies from the Sloan Digital Sky Survey

The sample of 137 low-redshift type Ia supernovae with 0.05 < z < 0.3 obtained from the SDSS-II Supernova Survey for the southern equatorial stripe of 300 square degrees is used to derive the luminosity functions of type Ia supernovae and of their host galaxies in the gri passbands. We show that the luminosity function of type Ia supernova host galaxies matches well with that of galaxies in the general field, suggesting that the occurrence of type Ia supernovae does not favour a particular type of galaxies but is predominantly proportional to the luminosity of galaxies. The only evidence that points to possible correlation between the supernova rate and star formation activity is that the supernova rate in late-type galaxies is higher than that in early-type galaxies by 31+/-35%. The sample contains 8 type Ia supernovae whose host galaxies were not identified, but it is shown that their occurrence is consistent with them occurred in low luminous galaxies beyond the survey. The luminosity function of type Ia supernovae is approximately Gaussian with the full-width half maximum being a factor of 1.4 in luminosity. The Gaussian distribution becomes tighter if the ratio of extinction to reddening, R_V, is lower than the characteristic value for the Milky Way and if luminosity is corrected for the light curve shape. The colour excess is ~0.07 mag which is significantly smaller than reddening expected for field galaxies. This colour excess does not vary with the distance of the supernovae from the centre of the host galaxy to 15 kpc. This suggests that the major part of the colour excess appears to be either intrinsic or reddening that arises in the immediate environment of supernova, rather than interstellar reddening in host galaxies.

preprint2009arXiv

Properties of Disks and Bulges of Spiral and Lenticular Galaxies in the Sloan Digital Sky Survey

A bulge-disk decomposition is made for 737 spiral and lenticular galaxies drawn from a SDSS galaxy sample for which morphological types are estimated. We carry out the bulge-disk decomposition using the growth curve fitting method. It is found that bulge properties, effective radius, effective surface brightness, and also absolute magnitude, change systematically with the morphological sequence; from early to late types, the size becomes somewhat larger, and surface brightness and luminosity fainter. In contrast disks are nearly universal, their properties remaining similar among disk galaxies irrespective of detailed morphologies from S0 to Sc. While these tendencies were often discussed in previous studies, the present study confirms them based on a large homogeneous magnitude-limited field galaxy sample with morphological types estimated. The systematic change of bulge-to-total luminosity ratio, $B/T$, along the morphological sequence is therefore not caused by disks but mostly by bulges. It is also shown that elliptical galaxies and bulges of spiral galaxies are unlikely to be in a single sequence. We infer the stellar mass density (in units of the critical mass density) to be $Ω=$0.0021 for spheroids, i.e., elliptical galaxies plus bulges of spiral galaxies, and $Ω=$0.00081 for disks.

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.

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

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

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.

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.

preprint2006arXiv

The Sloan Digital Sky Survey Monitor Telescope Pipeline

The photometric calibration of the Sloan Digital Sky Survey (SDSS) is a multi-step process which involves data from three different telescopes: the 1.0-m telescope at the US Naval Observatory (USNO), Flagstaff Station, Arizona (which was used to establish the SDSS standard star network); the SDSS 0.5-m Photometric Telescope (PT) at the Apache Point Observatory (APO), New Mexico (which calculates nightly extinctions and calibrates secondary patch transfer fields); and the SDSS 2.5-m telescope at APO (which obtains the imaging data for the SDSS proper). In this paper, we describe the Monitor Telescope Pipeline, MTPIPE, the software pipeline used in processing the data from the single-CCD telescopes used in the photometric calibration of the SDSS (i.e., the USNO 1.0-m and the PT). We also describe transformation equations that convert photometry on the USNO-1.0m u'g'r'i'z' system to photometry the SDSS 2.5m ugriz system and the results of various validation tests of the MTPIPE software. Further, we discuss the semi-automated PT factory, which runs MTPIPE in the day-to-day standard SDSS operations at Fermilab. Finally, we discuss the use of MTPIPE in current SDSS-related projects, including the Southern u'g'r'i'z' Standard Star project, the u'g'r'i'z' Open Star Clusters project, and the SDSS extension (SDSS-II).

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.

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.

preprint2004arXiv

Radiative Corrections to Neutrino Reactions off Proton and Deuteron

Radiative corrections are calculated for antineutrino proton quasielastic scattering, neutrino deuteron scattering, and the asymmetry of polarised neutron beta decay from which $G_{A}/G_{V}$ is determined. A particular emphasis is given to the constant parts that are usually absorbed into the coupling constants, and thereby those that appear in the processes that concern us are unambiguously tied among each other.

preprint2004arXiv

The Galaxy-mass Correlation Function Measured from Weak Lensing in the SDSS

We present galaxy-galaxy lensing measurements over scales 0.025 to 10 Mpc/h in the Sloan Digital Sky Survey. Using a flux-limited sample of 127,001 lens galaxies with spectroscopic redshifts and mean luminosity <L> = L_* and 9,020,388 source galaxies with photometric redshifts, we invert the lensing signal to obtain the galaxy-mass correlation function xi_{gm}. We find xi_{gm} is consistent with a power-law, xi_{gm} = (r/r_0)^{-gamma}, with best-fit parameters gamma = 1.79 +/- 0.06 and r_0 = (5.4+/-0.7)(0.27/Omega_m)^{1/gamma} Mpc/h. At fixed separation, the ratio xi_{gg}/xi_{gm} = b/r where b is the bias and r is the correlation coefficient. Comparing to the galaxy auto-correlation function for a similarly selected sample of SDSS galaxies, we find that b/r is approximately scale independent over scales 0.2-6.7 Mpc/h, with mean <b/r> = (1.3+/-0.2)(Omega_m/0.27). We also find no scale dependence in b/r for a volume limited sample of luminous galaxies (-23.0 < M_r < -21.5). The mean b/r for this sample is <b/r>_{Vlim} = (2.0+/-0.7)(Omega_m/0.27). We split the lens galaxy sample into subsets based on luminosity, color, spectral type, and velocity dispersion, and see clear trends of the lensing signal with each of these parameters. The amplitude and logarithmic slope of xi_{gm} increases with galaxy luminosity. For high luminosities (L ~5 L_*), xi_{gm} deviates significantly from a power law. These trends with luminosity also appear in the subsample of red galaxies, which are more strongly clustered than blue galaxies.

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

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

An exact Polynomial Hybrid Monte Carlo algorithm for dynamical Kogut-Susskind fermions

We present a polynomial Hybrid Monte Carlo (PHMC) algorithm as an exact simulation algorithm with dynamical Kogut-Susskind fermions. The algorithm uses a Hermitian polynomial approximation for the fractional power of the KS fermion matrix. The systematic error from the polynomial approximation is removed by the Kennedy-Kuti noisy Metropolis test so that the algorithm becomes exact at a finite molecular dynamics step size. We performed numerical tests with $N_f$$=$2 case on several lattice sizes. We found that the PHMC algorithm works on a moderately large lattice of $16^4$ at $β$$=$5.7, $m$$=$0.02 ($m_{\mathrm{PS}}/m_{\mathrm{V}}$$\sim$0.69) with a reasonable computational time.

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

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 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.

preprint2002arXiv

The Broad-band Optical Properties of Galaxies with Redshifts 0.0 < z < 0.2

Using photometry and spectroscopy of 144,609 galaxies from the Sloan Digital Sky Survey, we present bivariate distributions of pairs of seven galaxy properties: four optical colors, surface brightness, radial profile shape as measured by the Sersic index, and absolute magnitude. In addition, we present the dependence of local galaxy density (smoothed on 8 h^{-1} Mpc scales) on all of these properties. Several classic, well-known relations among galaxy properties are evident at extremely high signal-to-noise ratio: the color-color relations of galaxies, the color-magnitude relations, the magnitude-surface brightness relation, and the dependence of density on color and absolute magnitude. We show that most of the i-band luminosity density in the universe is in the absolute magnitude and surface brightness ranges used. Some of the relationships between parameters, in particular the color--magnitude relations, show stronger correlations for exponential galaxies and concentrated galaxies taken separately than for all galaxies taken together. We provide a simple set of fits of the dependence of galaxy properties on luminosity for these two sets of galaxies.

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

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

Lepton Flavour Mixing Matrix and CP Violation from Neutrino Oscillation Experiments

The measurement of the charged-current 8B solar neutrino reaction on deuterium at the Sudbury Neutrino Observatory confirms the neutrino oscillation hypothesis for the solar neutrino problem, and the result favours the solution with large neutrino mixing angles. We demonstrate that the current neutrino oscillation data (including atmospheric and reactor neutrinos) are sufficient to construct the lepton flavour mixing matrix with a reasonable accuracy. We also infer the maximum size of CP violation effects consistent with the current neutrino oscillation experiments.

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

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

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.

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

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

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

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

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

Metal Abundance in the Solar Interior

It is shown that the metal abundance in the solar interior is constrained from the current solar neutrino experiment, even if one allows neutrino conversion in the sun due to neutrino oscillation. The result shows that the metal abundance in the interior should be within the range 0.4 - 1.4 times that in the surface, supporting the idea that the sun formed by a contraction of a gas cloud with an almost homogeneous composition.

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.

preprint1995arXiv

Hadron Scattering Lengths in Lattice QCD

Lattice QCD calculation of s-wave hadron scattering lengths in the channels π-π, π-N, K-N, {\bar K}-N and N-N is carried out in the quenched QCD at $β=6/g^2=5.7$. A variant of the method of wall sourceis developed for this purpose, which reduces the computer time by a factor L^3 on an L^3xT lattice compared to the conventional point source method and avoids the Fierz mixing problem. A version of the method in which gauge configurations are not fixed to any gauge can be extended to calculate disconnected quark loop con- tributions in hadron two- and three-point functions. An analytical estimate of statistical errors for this method is worked out, and the magnitude of errors without and with gauge fixing is compared for the case of π-πfour-point functions calculated with the KS quark action. For π-πscattering both I=0 and 2 scattering lengths are evaluated using the KS and Wilson quark actions on a 12^3x20 lattice. For the same size of lattice, π-N, K-N and {\bar K}-N scattering lenghts are calculated with the Wilson quark action. For the π-πand π-N cases simulation results are consistent with the predictions of current algebra and PCAC within one to two standard deviations up to quite heavy quark masses corresponding to $m_π/m_ρ\approx 0.74$, while for the K-N and {\bar K}-N cases the agreement is within a factor of two. For N-N scat- tering simulations with the Wilson action on a 20^4 lattice with heavy quarks with $m_π/m_ρ\approx 0.74-0.95$, where the deuteron is expected to become unbound from a phenomenological study with one-boson exchange potentials, show that the nucleon-nucleon force is attractive for both spin triplet and singlet channels, and that the scattering lengths are substantially larger compared to those for the π-πand π-N cases even for such heavy quarks.

preprint1994arXiv

An Exploratory Study of Nucleon-Nucleon Scattering Lengths in Lattice QCD

An exploratory study is made of the nucleon-nucleon $s$-wave scattering lengths in quenched lattice QCD with the Wilson quark action. The $π$-$N$ and $π$-$π$ scattering lengths are also calculated for comparison. The calculations are made with heavy quarks corresponding to $m_π/m_ρ\approx 0.73-0.95$. The results show that the $N$-$N$ system has an attractive force in both spin-singlet and triplet channels, with their scattering lengths significantly larger than those for the $π$-$N$ and $π$-$π$ cases, a trend which is qualitatively consistent with the experiment. Problems toward a more realistic calculation for light quarks are discussed.

preprint1993arXiv

Reheating during Hierarchical Clustering in the Universe Dominated by the Cold Dark Matter

We investigate reheating of the universe by early formation of stars and quasars in the hierarchical clustering scheme of cold dark matter scenario, with perturbation fluctuations normalized by the COBE data. It is found that ionizing uv flux from OB stars with the abundance given by the standard initial mass function is strong enough to ionize the universe from $z\approx 30$ to the present epoch, if 1--2\% of the collapsed baryons go into stars. This lessens significantly the CMB anisotropies at a small angular scale. Reionization also increases the Jeans mass to $M_{\rm luminous}\approx 10^9M_\odot$ for $z\lsim 10$, which leads to a cut-off of the luminosity function of normal galaxies on a faint side. A strong uv flux is expected at $z\approx 2-5$, and the null result of the Gunn-Peterson test is naturally explained. Early star formation also results in production of heavy elements, and the observational metal abundance sets a strong constraint on the photon energy injection into the intergalactic space.