Source author record

T. Umeda

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

60works
13topics
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

60 published item(s)

preprint2022arXiv

High-power laser experiment forming a supercritical collisionless shock in a magnetized uniform plasma at rest

We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures of plasma density and temperatures, which are compared with one-dimensional particle-in-cell simulations. It is indicated that just after the laser irradiation, the Al plasma is magnetized by a self-generated Biermann battery field, and the plasma slaps the incident N plasma. The compressed external field in the N plasma reflects N ions, leading to counter-streaming magnetized N flows. Namely we identify the edge of the reflected N ions. Such interacting plasmas form a magnetized collisionless shock.

preprint2015arXiv

Low energy spectral and scattering theory for relativistic Schroedinger operators

Spectral and scattering theory at low energy for the relativistic Schroedinger operator are investigated. Some striking properties at thresholds of this operator are exhibited, as for example the absence of 0-energy resonance. Low energy behavior of the wave operators and of the scattering operator are studied, and stationary expressions in terms of generalized eigenfunctions are proved for the former operators. Under slightly stronger conditions on the perturbation the absolute continuity of the spectrum on the positive semi axis is demonstrated. Finally, an explicit formula for the action of the free evolution group is derived. Such a formula, which is well known in the usual Schroedinger case, was apparently not available in the relativistic setting.

preprint2015arXiv

Multipoint reweighting method and beta functions for the calculation of QCD equation of state

We study a reweighting method aiming at numerical studies of QCD at finite density, in which the conventional Monte-Carlo method cannot be applied directly. One of the most important problems in the reweighting method is the overlap problem. To solve it, we propose to perform simulations at several simulation points and combine their results in the data analyses. In this report, we introduce this multipoint reweighting method and test if the method works well by measuring histograms of physical quantities. Using this method, we calculate the meson masses as continuous functions of the gauge coupling beta and the hopping parameters kappa in QCD at zero density. We then determine lines of constant physics in the (beta, kappa) space and evaluate the derivatives of the lattice spacing with respect to beta and kappa along the lines of constant physics (inverse of the beta functions), which are needed in a calculation of the equation of state.

preprint2015arXiv

Multipoint reweighting method and its applications to lattice QCD

The reweighting method is widely used in numerical studies of QCD, in particular, for the cases in which the conventional Monte-Carlo method cannot be applied directly, e.g., finite density QCD. However, the application range of the reweighing method is restricted due to several problems. One of the most severe problems here is the overlap problem. To solve it, we examine a multipoint reweighting method in which simulations at several simulation points are combined in the data analyses. We systematically study the applicability and limitation of the multipoint reweighting method in two-flavor QCD at zero density. Measuring histograms of physical quantities at a series of simulation points, we apply the multipoint reweighting method to calculate the meson masses as continuous functions of the gauge coupling $β$ and the hopping parameters $κ$. We then determine lines of constant physics and beta functions, which are needed in a calculation of the equation of state at finite temperature.

preprint2015arXiv

Towards a spin-ensemble quantum memory for superconducting qubits

This article reviews efforts to build a new type of quantum device, which combines an ensemble of electronic spins with long coherence times, and a small-scale superconducting quantum processor. The goal is to store over long times arbitrary qubit states in orthogonal collective modes of the spin-ensemble, and to retrieve them on-demand. We first present the protocol devised for such a multi-mode quantum memory. We then describe a series of experimental results using NV center spins in diamond, which demonstrate its main building blocks: the transfer of arbitrary quantum states from a qubit into the spin ensemble, and the multi-mode retrieval of classical microwave pulses down to the single-photon level with a Hahn-echo like sequence. A reset of the spin memory is implemented in-between two successive sequences using optical repumping of the spins.

preprint2015arXiv

Towards the QCD equation of state at the physical point using Wilson fermion

We study the (2+1)-flavor QCD at nonzero temperatures using nonperturbatively improved Wilson quarks of the physical masses by the fixed scale approach. We perform physical point simulations at finite temperatures with the coupling parameters which were adopted by the PACS-CS Collaboration in their studies using the reweighting technique. Zero temperature values are obtained on the PACS-CS configurations which are open to the public on the ILDG/JLDG. Finite temperature configurations are generated with the RHMC algorithm. The lattice sizes are $32^3 \times N_t$ with $N_t=14$, 13, $\cdots$, 4 which correspond to $T \approx 160$--550 MeV. We present results of some basic observables at these temperatures and the status of our calculation of the equation of state.

preprint2014arXiv

Histograms in heavy-quark QCD at finite temperature and density

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

preprint2014arXiv

Isotopic identification of engineered nitrogen-vacancy spin qubits in ultrapure diamond

Nitrogen impurities help to stabilize the negatively-charged-state of NV$^-$ in diamond, whereas magnetic fluctuations from nitrogen spins lead to decoherence of NV$^-$ qubits. It is not known what donor concentration optimizes these conflicting requirements. Here we used 10-MeV $^{15}$N$^{3+}$ ion implantation to create NV$^-$ in ultrapure diamond. Optically detected magnetic resonance of single centers revealed a high creation yield of $40\pm3$% from $^{15}$N$^{3+}$ ions and an additional yield of $56\pm3$% from $^{14}$N impurities. High-temperature anneal was used to reduce residual defects, and charge stable NV$^-$, even in a dilute $^{14}$N impurity concentration of 0.06 ppb were created with long coherence times.

preprint2014arXiv

Multi-mode storage and retrieval of microwave fields in a spin ensemble

A quantum memory at microwave frequencies, able to store the state of multiple superconducting qubits for long times, is a key element for quantum information processing. Electronic and nuclear spins are natural candidates for the storage medium as their coherence time can be well above one second. Benefiting from these long coherence times requires to apply the refocusing techniques used in magnetic resonance, a major challenge in the context of hybrid quantum circuits. Here we report the first implementation of such a scheme, using ensembles of nitrogen-vacancy (NV) centres in diamond coupled to a superconducting resonator, in a setup compatible with superconducting qubit technology. We implement the active reset of the NV spins into their ground state by optical pumping and their refocusing by Hahn echo sequences. This enables the storage of multiple microwave pulses at the picoWatt level and their retrieval after up to $35 μ$s, a three orders of magnitude improvement compared to previous experiments.

preprint2013arXiv

Extending spin coherence times of diamond qubits by high temperature annealing

Spins of negatively charged nitrogen-vacancy (NV$^-$) defects in diamond are among the most promising candidates for solid-state qubits. The fabrication of quantum devices containing these spin-carrying defects requires position-controlled introduction of NV$^-$ defects having excellent properties such as spectral stability, long spin coherence time, and stable negative charge state. Nitrogen ion implantation and annealing enable the positioning of NV$^-$ spin qubits with high precision, but to date, the coherence times of qubits produced this way are short, presumably because of the presence of residual radiation damage. In the present work, we demonstrate that a high temperature annealing at 1000$^\circ$C allows 2 millisecond coherence times to be achieved at room temperature. These results were obtained for implantation-produced NV$^-$ defects in a high-purity, 99.99% $^{12}$C enriched single crystal chemical vapor deposited diamond. We discuss these remarkably long coherence times in the context of the thermal behavior of residual defect spins. [Published in Physical Review B {\bf{88}}, 075206 (2013)]

preprint2013arXiv

Scaling properties of the chiral phase transition in the low density region of two-flavor QCD with improved Wilson fermions

We study scaling behavior of a chiral order parameter in the low density region, performing a simulation of two-flavor QCD with improved Wilson quarks. The scaling behavior of the chiral order parameter defined by a Ward-Takahashi identity agrees with the scaling function of the three-dimensional O(4) spin model at zero chemical potential. We extend the scaling study to finite density QCD. Applying the reweighting method and calculating derivatives of the chiral order parameter with respect to the chemical potential, the scaling properties of the chiral phase transition are discussed in the low density region. We moreover calculate the curvature of the phase boundary of the chiral phase transition in the temperature and chemical potential plane assuming the O(4) scaling relation.

preprint2012arXiv

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

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

preprint2012arXiv

Efficiently Engineered Room Temperature Single Photons in Silicon Carbide

We report the first observation of stable single photon sources in silicon carbide (SiC). These sources are extremely bright and operate at room temperature demonstrating that SiC is a viable material in which to realize various quantum information, computation and photonic applications. The maximum single photon count rate detected is 700k counts/s with an inferred quantum efficiency around 70%. The single photon sources are due to intrinsic deep level defects constituted of carbon antisite-vacancy pairs. These are shown to be formed controllably by electron irradiation. The variability of the temporal kinetics of these single defects is investigated in detail.

preprint2012arXiv

Electron spin resonance detected by a superconducting qubit

A new method for detecting the magnetic resonance of electronic spins at low temperature is demonstrated. It consists in measuring the signal emitted by the spins with a superconducting qubit that acts as a single-microwave-photon detector, resulting in an enhanced sensitivity. We implement this new type of electron-spin resonance spectroscopy using a hybrid quantum circuit in which a transmon qubit is coupled to a spin ensemble consisting of NV centers in diamond. With this setup we measure the NV center absorption spectrum at 30mK at an excitation level of \thicksim15\,μ_{B} out of an ensemble of 10^{11} spins.

preprint2012arXiv

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

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

preprint2012arXiv

Finite density QCD phase transition in the heavy quark region

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

preprint2012arXiv

Long coherence time of spin qubits in $^{12}$C enriched polycrystalline CVD diamond

Single defects in diamond and especially negatively charged nitrogen vacancy (NV) centers are very promising quantum systems with wide applications in physics and biology. It was shown that their coherence properties can be strongly improved by growing ultrapure diamond with low concentration of parasitic spins associated with nitrogen electron spins and nuclear spins related to $^{13}$C carbon isotope. Here we report a high quality $^{12}$C-enriched polycrystalline CVD diamond material with properties comparable with single crystals. We find single NVs in the grains of this material, which show extremely long electron spin coherence time $T_2 > 2\,ms$.

preprint2012arXiv

Phase structure of finite temperature QCD in the heavy quark region

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

preprint2012arXiv

Probability distribution functions in the finite density lattice QCD

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

preprint2011arXiv

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

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

preprint2011arXiv

Histogram method in finite density QCD with phase quenched simulations

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

preprint2011arXiv

Hybrid quantum circuit with a superconducting qubit coupled to a spin ensemble

We report the experimental realization of a hybrid quantum circuit combining a superconducting qubit and an ensemble of electronic spins. The qubit, of the transmon type, is coherently coupled to the spin ensemble consisting of nitrogen-vacancy (NV) centers in a diamond crystal via a frequency-tunable superconducting resonator acting as a quantum bus. Using this circuit, we prepare arbitrary superpositions of the qubit states that we store into collective excitations of the spin ensemble and retrieve back later on into the qubit. These results constitute a first proof of concept of spin-ensemble based quantum memory for superconducting qubits.

preprint2011arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

QCD thermodynamics at zero and finite densities with improved Wilson quarks

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

preprint2010arXiv

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

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

preprint2009arXiv

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

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

preprint2009arXiv

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

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

preprint2008arXiv

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

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

preprint2008arXiv

Magnetic and electric screening masses from Polyakov-loop correlations

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

preprint2008arXiv

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

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

preprint2005arXiv

Properties of hadron screening masses at small baryonic density

The properties of hadron screening masses around the deconfinement phase transition at finite baryonic density can be studied by evaluating the Taylor coefficients with respect to the iso-scalar and iso-vector chemical potentials. We simulate 2-flavour lattice QCD with staggered fermions on a 12*12*24*6 lattice with ma = 0.05 and 0.10 and report investigations of nucleon, pseudo-scalar and vector mesons. We present new, strong evidence that in the confining phase, the screening masses at vanishing chemical potential have significant temperature dependence, but the effect of the iso-scalar chemical potential is very small. Above the critical temperature, the second derivative terms of mesons rapidly increase as contrasted to the case of baryon. We also study the responses of the screening masses to an iso-vector chemical potential and discuss some of the issues related to the properties of hadron masses at finite density.

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

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

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

preprint2004arXiv

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

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

preprint2003arXiv

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

Quenched charmonium spectrum

We study charmonium using the standard relativistic formalism in the quenched approximation, on a set of lattices with isotropic lattice spacings ranging from 0.1 to 0.04 fm. We concentrate on the calculation of the hyperfine splitting between eta_c and J/psi, aiming for a controlled continuum extrapolation of this quantity. The splitting extracted from the non-perturbatively improved clover Dirac operator shows very little dependence on the lattice spacing for $a \leq 0.1$ fm. The dependence is much stronger for Wilson and tree-level improved clover operators, but they still yield consistent extrapolations if sufficiently fine lattices, $a \leq 0.07$ fm ($a M(η_c) \leq 1$), are used. Our result for the hyperfine splitting is 77(2)(6) MeV (where Sommer's parameter, r_0, is used to fix the scale). This value remains about 30% below experiment. Dynamical fermions and OZI-forbidden diagrams both contribute to the remainder. Results for the eta_c and J/psi wave functions are also presented.

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

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

Renormalization group flow of SU(3) lattice gauge theory - Numerical studies in a two coupling space

We investigate the renormalization group (RG) flow of SU(3) lattice gauge theory in a two coupling space with couplings $β_{11}$ and $β_{12}$ corresponding to $1\times 1$ and $1\times 2$ loops respectively. Extensive numerical calculations of the RG flow are made in the fourth quadrant of this coupling space, i.e., $β_{11}>0$ and $β_{12}<0$. Swendsen's factor two blocking and the Schwinger-Dyson method are used to find an effective action for the blocked gauge field. The resulting renormalization group flow runs quickly towards an attractive stream which has an approximate line shape. This is numerical evidence of a renormalized trajectory which locates close to the two coupling space. A model flow equation which incorporates a marginal coupling (asymptotic scaling term), an irrelevant coupling and a non-perturbative attraction towards the strong coupling limit reproduces qualitatively the observed features. We further examine the scaling properties of an action which is closer to the attractive stream than the currently used improved actions. It is found that this action shows excellent restoration of rotational symmetry even for coarse lattices with $a \sim 0.3$ fm.

preprint2001arXiv

Screening mass responses to the chemical potential at finite temperature

Responses to chemical potential of the pseudoscalar meson screening mass and the chiral condensate in lattice QCD are investigated. On a $16 \times 8^2 \times \4$ lattice with two flavors of staggered quarks the first and second responses below and above $T_c$ are evaluated. Different behavior in the low and the high temperature phases are observed, which may be explained as a consequence of the chiral symmetry breaking and restoration.