Search

We investigate the nature of the chiral phase transition in the massless two-flavor QCD using the renormalization group improved gauge action and the Wilson quark action on $32^3\times 16$, $24^3\times 12$, and $16^3\times 8$ lattices. We calculate the spacial and temporal propagators of the iso-triplet mesons in the pseudo-scalar ($PS$), scalar ($S$), vector ($V$) and axial-vector ($AV$) channels on the lattice of three sizes. We first verify that the RG scaling is excellently satisfied for all cases. This is consistent with the claim that the chiral phase transition is second order. Then we compare the spacial and temporal effective masses between the axial partners, i.e. $PS$ vs $S$ and $V$ vs $AV$, on each of the three size lattices. We find the effective masses of all of six cases for the axial partners agree remarkably. This is consistent with the claim that at least $Z_4$ subgroup of the $U_A(1)$ symmetry in addition to the $SU_A(2)$ symmetry is recovered at the chiral phase transition point., Comment: 4 figures, 1 table. arXiv admin note: text overlap with arXiv:1704.03134

We investigate the nature of the chiral phase transition in the massless two-flavor QCD using the renormalization group improved gauge action and the Wilson quark action on $32^3\times 16$, $24^3\times 12$, and $16^3\times 8$ lattices. Based on the renormalization group equation, we derive the scaling relation for the effective masses of mesons at the chiral phase transition point. If the chiral phase transition is second order, the effective masses as a function of the rescaled time/space do not depend on the lattice size and show the universal behavior. We find that our numerical simulations on the three sizes of lattices are excellently on the scaling curves, which is consistent with the second order phase transition., Comment: 9 pages, 4 figures, 2 tables

We investigate the properties of quarks and gluons above the chiral phase transition temperature $T_c,$ using the RG improved gauge action and the Wilson quark action with two degenerate quarks mainly on a $32^3\times 16$ lattice. In the one-loop perturbation theory, the thermal ensemble is dominated by the gauge configurations with effectively $Z(3)$ center twisted boundary conditions, making the thermal expectation value of the spatial Polyakov loop take a non-trivial $Z(3)$ center. This is in agreement with our lattice simulation of high temperature QCD. We further observe that the temporal propagator of massless quarks at extremely high temperature $\beta=100.0 \, (T \simeq10^{58} T_c)$ remarkably agrees with the temporal propagator of free quarks with the $Z(3)$ twisted boundary condition for $t/L_t \geq 0.2$, but differs from that with the $Z(3)$ trivial boundary condition. As we increase the mass of quarks $m_q$, we find that the thermal ensemble continues to be dominated by the $Z(3)$ twisted gauge field configurations as long as $m_q \le 3.0 \, T$ and above that the $Z(3)$ trivial configurations come in. The transition is essentially identical to what we found in the departure from the conformal region in the zero-temperature many-flavor conformal QCD on a finite lattice by increasing the mass of quarks. We argue that the behavior is consistent with the renormalization group analysis at finite temperature., Comment: 16 pages, 9 figures; 4 tables, an appendix added

We investigate SU(3) gauge theories in four dimensions with Nf fundamental fermions, on a lattice using the Wilson fermion. Clarifying the vacuum structure in terms of Polyakov loops in spatial directions and properties of temporal propagators using a new method "local analysis", we conjecture that the "conformal region" exists together with the confining region and the deconfining region in the phase structure parametrized by beta and K, both in the cases of the large Nf QCD within the conformal window (referred as Conformal QCD) with an IR cutoff and small Nf QCD at T/Tc>1 with Tc being the chiral transition temperature (referred as High Temperature QCD). Our numerical simulation on a lattice of the size 16^3 x 64 shows the following evidence of the conjecture. In the conformal region we find the vacuum is the nontrivial Z(3) twisted vacuum modified by non-perturbative effects and temporal propagators of meson behave at large t as a power-law corrected Yukawa-type decaying form. The transition from the conformal region to the deconfining region or the confining region is a sharp transition between different vacua and therefore it suggests a first order transition both in Conformal QCD and in High Temperature QCD. Within our fixed lattice simulation, we find that there is a precise correspondence between Conformal QCD and High Temperature QCD in the temporal propagators under the change of the parameters Nf and T/Tc respectively. In particular, we find the correspondence between Conformal QCD with Nf = 7 and High Temperature QCD with Nf=2 at T ~ 2 Tc being in close relation to a meson unparticle model. From this we estimate the anomalous mass dimension gamma* = 1.2 (1) for Nf=7. We also show that the asymptotic state in the limit T/Tc --> infty is a free quark state in the Z(3) twisted vacuum., Comment: 48 pages, 47 figures and 7tables; two subsections added on "in relation to mass deformed CFT" and "finite size scaling"; a version to be published in PRD

We give a new perspective on the properties of quarks and gluons at finite temperature T in N_f = 2 ~ 6 QCD. We point out the existence of an IR fixed point for the gauge coupling constant at T>T_c (T_c is the chiral phase transition temperature). Based on this observation we predict theoretically and verify numerically that the correlation functions of a meson G(t) at T/T_c > 1 decay with a power-law corrected Yukawa-type decaying form, G(t)=c exp(-m t)/t^alpha in the "conformal region" defined by m < c Lambda_IR, where Lambda_IR is the IR cutoff, m is the characteristic scale of the spectrum in the meson cannel and c is a constant of order 1. The decaying form is the characteristics of conformal theories with an IR cutoff. We discuss in detail how the resulting hyper scaling relation of physical observables may modify the existing argument about the order of the chiral phase transition in the N_f=2 case., Comment: 7 pages, 5figures, 1 table; revised to make the argument more clear with new figures and a table

We give a new perspective on the dynamics of conformal theories realized in the SU(N) gauge theory, when the number of flavors N_f is within the conformal window. Motivated by the RG argument on conformal theories with a finite IR cutoff \Lambda_{IR}, we conjecture that the propagator of a meson G_H(t) on a lattice behaves at large t as a power-law corrected Yukawa-type decaying form G_H(t) = c_H \exp{(-m_H t)}/t^{\alpha_H} instead of the exponentially decaying form c_H\exp{(-m_H t)}, in the small quark mass region where m_H \le c \Lambda_{IR}: m_H is the mass of the ground state hadron in the channel H and c is a constant of order 1. The transition between the "conformal region" and the "confining region" is a first order transition. Our numerical results verify the predictions for the N_f=7 case and the N_f=16 case in the SU(3) gauge theory with the fundamental representation., Comment: 4 pages, 4 figures; minor modifications as published

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., Comment: 4 pages, 3 figures, revtex4; v2: contents partly modified, published version

We discuss a possibility that the Neutron Electric Dipole Moment (NEDM) can be calculated in lattice QCD simulations in the presence of the CP violating $\theta$ term. In this paper we measure the energy difference between spin-up and spin-down states of the neutron in the presence of an uniform and static external electric field. We first test this method in quenched QCD with the RG improved gauge action on a $16^3\times 32$ lattice at $a^{-1}\simeq$ 2 GeV, employing two different lattice fermion formulations, the domain-wall fermion and the clover fermion for quarks, at relatively heavy quark mass $(m_{PS}/m_V \simeq 0.85)$. We obtain non-zero values of NEDM from calculations with both fermion formulations. We next consider some systematic uncertainties of our method for NEDM, using $24^3\times 32$ lattice at the same lattice spacing only with the clover fermion. We finally investigate the quark mass dependence of NEDM and observe a non-vanishing behavior of NEDM toward the chiral limit. We interpret this behavior as a manifestation of the pathology in the quenched approximation., Comment: LaTeX2e, 51 pages, 43 figures, uses revtex4 and graphicx, References and comments added, typos corrected, accepted by PRD

We report on a direct lattice calculation of the neutron EDM(NEDM) using the external electric field method in both quenched and full QCD. In quenched QCD, we use a $24^3\times 32$ lattice at $\beta=2.6$ with the Iwasaki gauge action and the clover fermion action to examine the viability of this method. In particular we investigate possible effects of violation of the periodic boundary condition of the external electric field on the NEDM signal. We also study the quark mass dependence of NEDM in quenched QCD, and observe that NEDM seems to remain non-zero toward the chiral limit because of the quenched artifact. In 2-flavor full QCD we employ configurations generated by the CP-PACS collaboration on a $24^3\times 48$ lattice at $\beta=2.1$ with the same gluon and quark actions as in the quenched case. Since the number of configurations is limited, we employ 8 different source points per one configuration and take an average over them. Our preliminary result at three quark masses ($m_{PS}/m_V\simeq 0.81, 0.76, 0.69$) indicates that non-zero value for NEDM can be obtained in full QCD. Statistical errors, however, are still too large to show the theoretically expected behaviour for NEDM in full QCD that it vanishes in the chiral limit., Comment: 7 pages, 15 figures, a talk given at the XXIV International Symposium on Lattice Field Theory (Lattice 2006, Hadronic Interactions and Structure), Tucson, Arizona, USA, July 23-28, 2006

We carry out a feasibility study toward a lattice QCD calculation of the neutron electric dipole moment (NEDM) in the presence of the $\theta$ term using two different approaches. In the first method, we calculate the CP-odd electromagnetic form factor $F_3$, which becomes the NEDM in the zero momentum transfer limit. At the first order in $\theta$, we derive a formula connecting the lattice three-point function to the CP-odd electromagnetic form factor. In the second method we directly extract the NEDM from the energy difference between spin-up and spin-down neutron states in the presence of a constant electric field, without expanding a small but non-zero $\theta$. We test both approaches numerically, employing the domain-wall quark action with the RG improved gauge action in quenched QCD at $a^{-1}\simeq 2$ GeV on a $16^3\times 32\times 16$ lattice, and further applying the second method to the clover quark action at a similar lattice spacing and nucleon mass. We obtain good signals from both approaches. In particular the second method works well with both fermion formulations., Comment: 6 pages, 7 eps-figures, talk presented at Lattice 2005 (Chiral fermions), Trinity College, Dublin, Ireland, 25-30 Jul 2005

We carry out a feasibility study for the lattice QCD calculation of the neutron electric dipole moment (NEDM) in the presence of the $\theta$ term. We develop the strategy to obtain the nucleon EDM from the CP-odd electromagnetic form factor $F_3$ at small $\theta$, in which NEDM is given by $\lim_{q^2\to 0}\theta F_3(q^2)/(2m_N)$ where $q$ is the momentum transfer and $m_N$ is the nucleon mass. We first derive a formula which relates $F_3$, a matrix element of the electromagnetic current between nucleon states, with vacuum expectation values of nucleons and/or the current. In the expansion of $\theta$, the parity-odd part of the nucleon-current-nucleon three-point function contains contributions not only from the parity-odd form factors but also from the parity-even form factors multiplied by the parity-odd part of the nucleon two-point function, and therefore the latter contribution must be subtracted to extract $F_3$. We then perform an explicit lattice calculation employing the domain-wall quark action with the RG improved gauge action in quenched QCD at $a^{-1}\simeq 2$ GeV on a $16^3\times 32\times 16$ lattice. At the quark mass $m_f a =0.03$, corresponding to $m_\pi/m_\rho \simeq 0.63$, we accumulate 730 configurations, which allow us to extract the parity-odd part in both two- and three-point functions. Employing two different Dirac $\gamma$ matrix projections, we show that a consistent value for $F_3$ cannot be obtained without the subtraction described above. We obtain $F_3(q^2\simeq 0.58 \textrm{GeV}^2)/(2m_N) =$ $-$0.024(5) $e\cdot$fm for the neutron and $F_3(q^2\simeq 0.58 \textrm{GeV}^2)/(2m_N) =$ 0.021(6) $e\cdot$fm for the proton., Comment: LaTeX2e, 43 pages, 42 eps figures, uses revtex4 and graphicx, comments added and typos corrected, final version to appear in Phys. Rev. D

We present a quenched lattice QCD calculation of the \alpha and \beta 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, |\alpha(NDR,2GeV)|=0.0090(09)(^{+5}_{-19})GeV^3 and |\beta(NDR,2GeV)|=0.0096(09)(^{+6}_{-20})GeV^3 with \alpha and \beta 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., Comment: 4 pages, 3 figures

We investigate the phase structure of lattice QCD for the general number of flavors in the parameter space of gauge coupling constant and quark mass, employing the one-plaquette gauge action and the standard Wilson quark action. Performing a series of simulations for the number of flavors $N_F=6$--360 with degenerate-mass quarks, we find that when $N_F \ge 7$ there is a line of a bulk first order phase transition between the confined phase and a deconfined phase at a finite current quark mass in the strong coupling region and the intermediate coupling region. The massless quark line exists only in the deconfined phase. Based on these numerical results in the strong coupling limit and in the intermediate coupling region, we propose the following phase structure, depending on the number of flavors whose masses are less than $\Lambda_d$ which is the physical scale characterizing the phase transition in the weak coupling region: When $N_F \ge 17$, there is only a trivial IR fixed point and therefore the theory in the continuum limit is free. On the other hand, when $16 \ge N_F \ge 7$, there is a non-trivial IR fixed point and therefore the theory is non-trivial with anomalous dimensions, however, without quark confinement. Theories which satisfy both quark confinement and spontaneous chiral symmetry breaking in the continuum limit exist only for $N_F \le 6$., Comment: RevTeX, 20 pages, 43 PS figures

We present an unquenched lattice calculation for the $B^0-\bar{B}^0$ transition amplitude. The calculation, carried out at an inverse lattice spacing $1/a$ = 2.22(4) GeV, incorporates two flavors of dynamical quarks described by the fully $O(a)$-improved Wilson fermion action and heavy quarks described by NRQCD. A particular attention is paid to the uncertainty that arises from the chiral extrapolation, especially the effect of pion loops, for light quarks, which we find could be sizable for the leptonic decay of the $B_d$ meson, whereas it is small for the $B_s$ meson and the $B$ parameters. We obtain $f_{B_d}=191(10)(^{+12}_{-22})$ MeV, $f_{B_s}/f_{B_d}=1.13(3)(^{+13}_{-\ 2})$, $B_{B_d}(m_b)=0.836(27)(^{+56}_{-62})$, $B_{B_s}/B_{B_d}=1.017(16)(^{+56}_{-17})$, $\xi=1.14(3)(^{+13}_{- 2})$, where the first error is statistical, and the second is systematic, including uncertainties due to chiral extrapolation, finite lattice spacing, heavy quark expansion and perturbative operator matching., Comment: 4 pages, 3 figures; version to appear in PRL

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 $\zeta=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 $\xi = 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 $\zeta$, 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 $\xi = 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., Comment: 43 pages, 49 eps figures, revtex; minor changes, version to appear in Physical Review D

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., Comment: 4 pages, 2 figures, minor corrections and added references

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., Comment: Latex 2.09 with PRD style files and epsf, 40 pages, 22 figures

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 $\pi$ and $\rho$ 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 $\phi$ 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%., Comment: 4 pages, revtex, 3 eps figures. Double statistics at beta=2.1, error estimates elaborated. Accepted for publication in Phys. Rev. Lett

We present a model-independent calculation of hadron matrix elements for all dimension-six operators associated with baryon number violating processes using lattice QCD. The calculation is performed with the Wilson quark action in the quenched approximation at $\beta=6/g^2=6.0$ on a $28^2\times 48\times 80$ lattice. Our results cover all the matrix elements required to estimate the partial lifetimes of (proton,neutron)$\to$($\pi,K,\eta$) +(${\bar \nu},e^+,\mu^+$) decay modes. We point out the necessity of disentangling two form factors that contribute to the matrix element; previous calculations did not make the separation, which led to an underestimate of the physical matrix elements. With a correct separation, we find that the matrix elements have values 3-5 times larger than the smallest estimates employed in phenomenological analyses of the nucleon decays, which could give strong constraints on several GUT models. We also find that the values of the matrix elements are comparable with the tree-level predictions of chiral lagrangian., Comment: 53 pages, 18 eps figures

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., Comment: Contribution to EPS 99, Tampere, Finland. LateX2e file, 2 pages long. Requires hep99 style files (included) and epsf. 3 Encapsulated Postscript figures. (Changed Author list.)

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%., Comment: 17 pages, 14 figures, revtex

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., Comment: 7 pages, Latex, 3 figures

The finite temperature transition of QCD with two degenerate light quarks is studied on the lattice with a renormalization group improved gauge action and the Wilson quark action. We have made simulations on an $8^3\times 4$ lattice near the chiral transition point. It is shown that the chiral condensate which is the order parameter of chiral symmetry satisfies remarkably a scaling relation with the exponents of the three dimensional O(4) Heisenberg model. This indicates that the chiral transition in two-flavor QCD is of second order in the continuum limit., Comment: 13 pages, 6 figures

We investigate and elucidate the phase structure of QCD for general number of flavors $N_F$ with Wilson quarks, varying $N_F$ from 2 up to 300. Based on numerical results combined with the result of the perturbation theory we propose the following picture: When $N_F \ge 17$, there is only a trivial fixed point and therefore the theory in the continuum limit is trivial. On the other hand, when $16 \ge N_F \ge 7$, there is a non-trivial fixed point and therefore the theory is non-trivial with anomalous dimensions, however, without quark confinement. Theories which satisfy both quark confinement and spontaneous chiral symmetry breaking in the continuum limit exist only for $N_F \le 6$. We also discuss the structure of the deconfining phase at finite temperatures for the small number of flavors such as $N_F=2$ and 3, through a systematic study of it for general number of flavors., Comment: LaTeX, 7 pages, 10 PS figures, Talk presented at LATTICE96(poster) and LATTICE96(finite temperature)