Source author record

Masakiyo Kitazawa

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

46works
8topics
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

46 published item(s)

preprint2023arXiv

Critical point in heavy-quark region of QCD on fine lattices

We perform a finite-size scaling analysis of the critical point in the heavy-quark region of QCD at nonzero temperature. Our previous analysis on the Binder cumulant at $N_t=4$ is extended to finer lattices with $N_t=6$ and $8$. The aspect ratio is also extended up to $15$ to suppress the non-singular contribution. High-precision analysis of the Binder cumulant is realized by an efficient Monte-Carlo simulation with the hopping-parameter expansion (HPE). Effects of higher-order terms in the HPE are incorporated by the reweighting method.

preprint2022arXiv

Anomalous enhancement of dilepton production as a precursor of color superconductivity

We compute the modification of the photon self-energy due to dynamical diquark fluctuations developed near the critical temperature of the color superconductivity through the Aslamasov-Larkin, Maki-Thompson and density of states terms, which are responsible for the paraconductivity in metals at vanishing energy and momentum. It is shown that the rate has a significant enhancement at low invariant-mass region over a rather wide range of temperature in the normal phase. This enhancement is worth exploration in the relativistic heavy-ion collisions, which may thereby reveal the significance of the diquark fluctuations in dense quark matter.

preprint2022arXiv

Scope and convergence of the hopping parameter expansion in finite temperature QCD with heavy quarks around the critical point

Hopping parameter expansion is a useful tool to investigate heavy dynamical quarks in lattice QCD, while the range of its applicability has been sometimes questioned. We study the convergence and the valid range of the hopping parameter expansion in the determination of the critical point (critical quark mass) of QCD with heavy quarks at finite temperature and density. On lattices with sufficiently large spatial extent, the terms in the hopping parameter expansion are classified into Wilson loop terms and Polyakov-type loop terms. We first study the case of the worst convergence in which all the gauge link variables are unit matrices and thus the Wilson loops and the Polyakov-type loops get their maximum values. We perform explicit calculation up to more than 100th order of the hopping parameter expansion. We show that the hopping parameter expansion is convergent up to the chiral limit of free Wilson quarks. We then perform a Monte-Carlo simulation to measure correlation among Polyakov-type loop terms up to the 20th order of the hopping parameter expansion. In previous studies, strong correlation between the leading order Polyakov loop term and the next-to-leading order bent Polyakov loop terms was reported and used to construct an effective theory to incorporate the next-to-leading order effect by a shift of the leading order coupling parameter. We establish that the strong correlation among Polyakov-type loop terms holds also at higher orders of the hopping parameter expansion, and extend the effective theory to incorporate higher-order effects up to high orders. Using the effective theory, we study the truncation error of the hopping parameter expansion. We find that the previous next-to-leading order result of the critical point for $N_t=4$ are well reliable. For $N_t \ge 6$, we need to incorporate higher-order effects in the effective theory.

preprint2021arXiv

Classifying Topological Charge in SU(3) Yang-Mills Theory with Machine Learning

We apply a machine learning technique for identifying the topological charge of quantum gauge configurations in four-dimensional SU(3) Yang-Mills theory. The topological charge density measured on the original and smoothed gauge configurations with and without dimensional reduction is used as inputs for the neural networks (NN) with and without convolutional layers. The gradient flow is used for the smoothing of the gauge field. We find that the topological charge determined at a large flow time can be predicted with high accuracy from the data at small flow times by the trained NN; for example, the accuracy exceeds $99\%$ with the data at $t/a^2\le0.3$. High robustness against the change of simulation parameters is also confirmed with a fixed physical volume. We find that the best performance is obtained when the spatial coordinates of the topological charge density are fully integrated out in preprocessing, which implies that our convolutional NN does not find characteristic structures in multi-dimensional space relevant for the determination of the topological charge.

preprint2021arXiv

Finite-size scaling around the critical point in the heavy quark region of QCD

Finite-size scaling is investigated in detail around the critical point in the heavy-quark region of nonzero temperature QCD. Numerical simulations are performed with large spatial volumes up to the aspect ratio $N_s/N_t=12$ at a fixed lattice spacing with $N_t=4$. We show that the Binder cumulant and the distribution function of the Polyakov loop follow the finite-size scaling in the $Z(2)$ universality class for large spatial volumes with $N_s/N_t \ge 9$, while, for $N_s/N_t \le 8$, the Binder cumulant becomes inconsistent with the $Z(2)$ scaling. To realize the large-volume simulations in the heavy-quark region, we adopt the hopping parameter expansion for the quark determinant: We generate gauge configurations using the leading order action including the Polyakov loop term for $N_t=4$, and incorporate the next-to-leading order effects in the measurements by the multipoint reweighting method. We find that the use of the leading-order configurations is crucially effective in suppressing the overlapping problem in the reweighting and thus reducing the statistical errors.

preprint2021arXiv

Latent heat and pressure gap at the first-order deconfining phase transition of SU(3) Yang-Mills theory using the small flow-time expansion method

We study latent heat and the pressure gap between the hot and cold phases at the first-order deconfining phase transition temperature of the SU(3) Yang-Mills theory. Performing simulations on lattices with various spatial volumes and lattice spacings, we calculate the gaps of the energy density and pressure using the small flow-time expansion (SFtX) method. We find that the latent heat $Δε$ in the continuum limit is $Δε/T^4 = 1.117 \pm 0.040$ for the aspect ratio $N_s/N_t=8$ and $1.349 \pm 0.038$ for $N_s/N_t=6$ at the transition temperature $T=T_c$. We also confirm that the pressure gap is consistent with zero, as expected from the dynamical balance of two phases at $T_c$. From hysteresis curves of the energy density near $T_c$, we show that the energy density in the (metastable) deconfined phase is sensitive to the spatial volume, while that in the confined phase is insensitive. Furthermore, we examine the effect of alternative procedures in the SFtX method - the order of the continuum and the vanishing flow-time extrapolations, and also the renormalization scale and higher-order corrections in the matching coefficients. We confirm that the final results are all very consistent with each other for these alternatives.

preprint2020arXiv

A Study of Stress-Tensor Distribution around Flux Tube in Abelian-Higgs Model

We study the stress-tensor distribution around the flux tube in static quark and anti-quark systems based on the momentum conservation and the Abelian-Higgs (AH) model. We first investigate constraints on the stress-tensor distribution from the momentum conservation and show that the effect of boundaries plays a crucial role to describe the structure of the flux tube in SU(3) Yang-Mills theory which has measured recently on the lattice. We then study the distributions of the stress tensor and energy density around the magnetic vortex with and without boundaries in the AH model, and compare them with the distributions in SU(3) Yang-Mills theory based on the dual superconductor picture. It is shown that a wide parameter range of the AH model is excluded by a comparison with the lattice results in terms of the stress tensor.

preprint2020arXiv

Calculation of PCAC mass with Wilson fermion using gradient flow

We calculate the PCAC mass for $(2+1)$ flavor full QCD with Wilson-type quarks. We adopt the Small Flow-time eXpansion (SFtX) method based on the gradient flow which provides us a general way to compute correctly renormalized observables even if the relevant symmetries for the observable are broken explicitly due to the lattice regularization, such as the Poincáre and chiral symmetries. Our calculation is performed on heavy $u, d$ quarks mass ($m_π/m_ρ\simeq0.63$) and approximately physical $s$ quark mass with fine lattice $a \simeq 0.07$~fm. The results are compared with those computed with the Schrödinger functional method.

preprint2020arXiv

Critical fluctuations in a dynamically expanding heavy-ion collision

For the discovery of the QCD critical point it is crucial to develop dynamical models of the fluctuations of the net-baryon number that can be embedded in simulations of heavy-ion collisions. In this proceeding, we study the dynamical formation of the critical fluctuations of the net-baryon number near the QCD critical point and their survival in the late stages in an expanding system. The stochastic diffusion equation with a non-linear free energy functional is employed for describing the evolution of conserved-charge fluctuations along trajectories in the crossover and first-order transition regions near the QCD critical point.

preprint2020arXiv

Distribution of Energy-Momentum Tensor around a Static Quark in the Deconfined Phase of SU(3) Yang-Mills Theory

Energy momentum tensor (EMT) characterizes the response of the vacuum as well as the thermal medium under the color electromagnetic fields. We define the EMT by means of the gradient flow formalism and study its spatial distribution around a static quark in the deconfined phase of SU(3) Yang-Mills theory on the lattice. Although no significant difference can be seen between the EMT distributions in the radial and transverse directions except for the sign, the temporal component is substantially different from the spatial ones near the critical temperature $T_c$. This is in contrast to the prediction of the leading-order thermal perturbation theory. The lattice data of the EMT distribution also indicate the thermal screening at long distance and the perturbative behavior at short distance.

preprint2020arXiv

End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD

We study the end point of the first-order deconfinement phase transition in two and 2+1 flavor QCD in the heavy quark region of the quark mass parameter space. We determine the location of critical point at which the first-order deconfinement phase transition changes to crossover, and calculate the pseudo-scalar meson mass at the critical point. Performing quenched QCD simulations on lattices with the temporal extents Nt=6 and 8, the effects of heavy quarks are determined using the reweighting method. We adopt the hopping parameter expansion to evaluate the quark determinants in the reweighting factor. We estimate the truncation error of the hopping parameter expansion by comparing the results of leading and next-to-leading order calculations, and study the lattice spacing dependence as well as the spatial volume dependence of the result for the critical point. The overlap problem of the reweighting method is also examined. Our results for Nt=4 and 6 suggest that the critical quark mass decreases as the lattice spacing decreases and increases as the spatial volume increases.

preprint2020arXiv

Issues with Search for Critical Point in QCD with Relativistic Heavy Ion Collisions

A systematic search for a critical point in the phase diagram of QCD matter is underway at the Relativistic Heavy Ion Collider (RHIC) and is planned at several future facilities. Its existence, if confirmed, and its location will greatly enhance our understanding of QCD. In this note we emphasize several important issues that are often not fully recognized in theoretical interpretations of experimental results relevant to the critical point search. We discuss ways in which our understanding on these issues can be improved.

preprint2019arXiv

Classifying topological sector via machine learning

We employ a machine learning technique for an estimate of the topological charge $Q$ of gauge configurations in SU(3) Yang-Mills theory in vacuum. As a first trial, we feed the four-dimensional topological charge density with and without smoothing into the convolutional neural network and train it to estimate the value of $Q$. We find that the trained neural network can estimate the value of $Q$ from the topological charge density at small flow time with high accuracy. Next, we perform the dimensional reduction of the input data as a preprocessing and analyze lower dimensional data by the neural network. We find that the accuracy of the neural network does not have statistically-significant dependence on the dimension of the input data. From this result we argue that the neural network does not find characteristic features responsible for the determination of $Q$ in the higher dimensional space.

preprint2019arXiv

Dynamically Integrated Transport Approach for High-Energy Nuclear Collisions at High Baryon Density

To explore the structure of the QCD phase diagram in high baryon density domain, several high-energy nuclear collision experiments in a wide range of beam energies are currently performed or planned using many accelerator facilities. In these experiments search for a first-order phase transition and the QCD critical point is one of the most important topics. To find the signature of the phase transition, experimental data should be compared to appropriate dynamical models which quantitatively describe the process of the collisions. In this study we develop a new dynamical model on the basis of the non-equilibrium hadronic transport model JAM and 3+1D hydrodynamics. We show that the new model reproduce well the experimental beam-energy dependence of hadron yields and particle ratio by the partial thermalization of the system in our core-corona approach.

preprint2016arXiv

Efficient formulas for efficiency correction of cumulants

We derive formulas which connect cumulants of particle numbers observed with efficiency losses with the original ones based on the binomial model. These formulas can describe the case with multiple efficiencies in a compact form. Compared with the presently suggested ones based on factorial moments, these formulas would drastically reduce the numerical cost for efficiency corrections when the order of the cumulant and the number of different efficiencies are large. The efficiency correction with realistic $p_T$-dependent efficiency would be carried out with the aid of these formulas.

preprint2016arXiv

Equation of state in (2+1)-flavor QCD with gradient flow

The energy-momentum tensor and equation of state are studied in finite-temperature (2+1)-flavor QCD with improved Wilson quarks using the method proposed by Makino and Suzuki based on the gradient flow. We find that the results of the gradient flow are consistent with the previous results using the $T$-integration method at $T \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 280$ MeV ($N_t\mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}}10$), while a disagreement is found at $T \mathrel{\rlap{\raise 0.511ex \hbox{$>$}}{\lower 0.511ex \hbox{$\sim$}}} 350$ MeV ($N_t \mathrel{\rlap{\raise 0.511ex \hbox{$<$}}{\lower 0.511ex \hbox{$\sim$}}} 8$) presumably due to the small-$N_t$ lattice artifact. We also report on the results on the renormalized chiral condensate and its disconnected susceptibility using the method of Hieda and Suzuki. The results show a clear signal of the expected chiral restoration crossover even with Wilson-type quarks which violate the chiral symmetry explicitly.

preprint2016arXiv

Fluctuations of conserved charges in relativistic heavy ion collisions: An introduction

Bulk fluctuations of conserved charges measured by event-by-event analysis in relativistic heavy ion collisions are observables which are believed to carry significant amount of information on the hot medium created by the collisions. Active studies have been done recently experimentally, theoretically, and on the lattice. In particular, non-Gaussianity of the fluctuations has acquired much attention recently. In this review, we give a pedagogical introduction to these issues, and survey recent developments in this field of research. Starting from the definition of cumulants, basic concepts in fluctuation physics, such as thermal fluctuations in statistical mechanics and time evolution of fluctuations in diffusive systems, are described. Phenomena which are expected to occur in finite temperature and/or density QCD matter and their measurement by event-by-event analyses are also elucidated.

preprint2016arXiv

Latent heat at the first order phase transition point of SU(3) gauge theory

We calculate the energy gap (latent heat) and pressure gap between the hot and cold phases of the SU(3) gauge theory at the first order deconfining phase transition point. We perform simulations around the phase transition point with the lattice size in the temporal direction Nt=6, 8 and 12 and extrapolate the results to the continuum limit. We also investigate the spatial volume dependence. The energy density and pressure are evaluated by the derivative method with non-perturabative anisotropy coefficients. We adopt a multi-point reweighting method to determine the anisotropy coefficients. We confirm that the anisotropy coefficients approach the perturbative values as Nt increases. We find that the pressure gap vanishes at all values of Nt when the non-perturbative anisotropy coefficients are used. The spatial volume dependence in the latent heat is found to be small on large lattices. Performing extrapolation to the continuum limit, we obtain $ Δε/T^4 = 0.75 \pm 0.17 $ and $ Δ(ε-3 p)/T^4 = 0.623 \pm 0.056.$

preprint2016arXiv

Non-perturbative production rate of photons with a lattice quark propagator: effect of vertex correction

We analyze the production rate of photons from the thermal medium above the deconfinement temperature with a quark propagator obtained from a lattice QCD numerical simulation. The photon-quark vertex is determined gauge-invariantly, so as to satisfy the Ward-Takahashi identity. The obtained photon production rate shows a suppression compared to perturbative results.

preprint2016arXiv

Temperature dependence of topological susceptibility using gradient flow

We study temperature dependence of the topological susceptibility with the $N_{f}=2+1$ flavors Wilson fermion. We have two major interests in this paper. One is a comparison of gluonic and fermionic definitions of the topological susceptibility. Two definitions are related by the chiral Ward-Takahashi identity but their coincidence is highly non-trivial for the Wilson fermion. By applying the gradient flow both for the gauge and quark fields we find a good agreement of these two measurements. The other is a verification of a prediction of the dilute instanton gas approximation at low temperature region $T_{pc}< T<1.5T_{pc}$, for which we confirm the prediction that the topological susceptibility decays with power $χ_{t}\propto(T/T_{pc})^{-8}$ for three flavors QCD.

preprint2016arXiv

Thermal blurring of event-by-event fluctuations provoked by rapidity conversion

We study the effect of thermal blurring caused by the use of (momentum-space) rapidity as a proxy of coordinate-space rapidity in experimental measurements of conserved charge fluctuations in relativistic heavy ion collisions. In theoretical studies assuming statistical mechanics, calculated fluctuations are those in a spatial volume. Experiments, on the other hand, can measure fluctuations only in a momentum-space in the final state. In a standard argument to compare experimental results for a momentum space with theoretical studies for a coordinate space, rapidities of particles are implicitly regarded as equivalent to their coordinate-space rapidity. We show that the relation of two fluctuations is significantly altered by the existence of the thermal motion, i.e. thermal blurring. We discuss that the thermal blurring can be regarded as a part of the diffusion process, and the effect can be understood by studying the rapidity window dependences of fluctuations. Centrality dependence of the thermal blurring effect is also discussed.

preprint2016arXiv

Understanding experimentally-observed fluctuations

We discuss two topics on the experimental measurements of fluctuation observables in relativistic heavy-ion collisions. First, we discuss the effects of the thermal blurring, i.e. the blurring effect arising from the experimental measurement of fluctuations in momentum space as a proxy of the thermal fluctuations defined in coordinate space, on higher order cumulants. Second, we discuss the effect of imperfect efficiency of detectors on the measurement of higher order cumulants. We derive effective formulas which can carry out the correction of this effect for higher order cumulants based on the binomial model.

preprint2015arXiv

Determination of Reference Scales for Wilson Gauge Action from Yang--Mills Gradient Flow

A parametrization of the lattice spacing ($a$) in terms of the bare coupling ($β$) for the SU(3) Yang--Mills theory with the Wilson gauge action is given in a wide range of~$β$. The Yang--Mills gradient flow with respect to the flow time~$t$ for the dimensionless observable, $t\frac{d}{dt}t^2\langle E(t)\rangle$, is utilized to determine the parametrization. With fine lattice spacings ($6.3\leβ\le7.5$) and large lattice volumes ($N_{\rm s}=64$--$128$), the discretization and finite-volume errors are significantly reduced to the same level as the statistical error.

preprint2015arXiv

Dilepton production spectrum above Tc with a lattice quark propagator

The dilepton production rate from the deconfined medium is analyzed with the photon self-energies constructed from quark propagators obtained by lattice numerical simulation for two values of temperature $T=1.5T_{\rm c}$ and $3T_{\rm c}$ above the critical temperature $T_{\rm c}$. The photon self-energy is calculated by the Schwinger-Dyson equation with the lattice quark propagtor and a vertex function determined so as to satisfy the Ward-Takahashi identity. The obtained dilepton production rate at zero momentum exhibits divergences reflecting van Hove singularity, and is significantly enhanced around $ω\simeq T$ compared with the rate obtained by the perturbative analysis.

preprint2015arXiv

Rapidity window dependences of higher order cumulants and diffusion master equation

We study the rapidity window dependences of higher order cumulants of conserved charges observed in relativistic heavy ion collisions. The time evolution and the rapidity window dependence of the non-Gaussian fluctuations are described by the diffusion master equation. Analytic formulas for the time evolution of cumulants in a rapidity window are obtained for arbitrary initial conditions. We discuss that the rapidity window dependences of the non-Gaussian cumulants have characteristic structures reflecting the non-equilibrium property of fluctuations, which can be observed in relativistic heavy ion collisions with the present detectors. It is argued that various information on the thermal and transport properties of the hot medium can be revealed experimentally by the study of the rapidity window dependences, especially by the combined use, of the higher order cumulants. Formulas of higher order cumulants for a probability distribution composed of sub-probabilities, which are useful for various studies of non-Gaussian cumulants, are also presented.

preprint2015arXiv

Thermodynamics and reference scale of SU(3) gauge theory from gradient flow on fine lattices

We study the parametrization of lattice spacing and thermodynamics of SU(3) gauge theory on the basis of the Yang-Mills gradient flow on fine lattices. The lattice spacing of the Wilson gauge action is determined over a wide range $6.3\leβ\le7.5$ with high accuracy. The measurements of the flow time and lattice spacing dependences of the expectation values of the energy-momentum tensor are performed on fine lattices.

preprint2015arXiv

Thermodynamics of SU(3) Gauge Theory from Gradient Flow

A novel method to study the bulk thermodynamics in lattice gauge theory is proposed on the basis of the Yang-Mills gradient flow with a fictitious time t. The energy density (epsilon) and the pressure (P) of SU(3) gauge theory at fixed temperature are calculated directly on 32^3 x (6,8,10) lattices from the thermal average of the well-defined energy-momentum tensor (T_{mu nu}^R(x)) obtained by the gradient flow. It is demonstrated that the continuum limit can be taken in a controlled manner from the t-dependence of the flowed data.

preprint2014arXiv

Effects of global charge conservation on time evolution of cumulants of conserved charges in relativistic heavy ion collisions

We investigate the effect of the global charge conservation on the cumulants of conserved charges observed in relativistic heavy ion collisions in a finite rapidity window, $Δη$, with a special emphasis on the time evolution of fluctuations in the hadronic medium. It is argued that the experimental result of the net-electric charge fluctuation observed by ALICE does not receive effects from the global charge conservation, because of the finite diffusion distance of charged particles in the hadronic stage. We emphasize that the magnitude of the effect of the global charge conservation can be estimated experimentally by combining the information on the $Δη$ dependences of various cumulants of conserved charges, similarly to other dynamical properties of the hot medium.

preprint2014arXiv

Emergence of soft quark excitations by the coupling with a soft mode of the QCD critical point

We study the quark spectrum at nonzero temperature and density near the critical point (CP) of the chiral phase transition incorporating effects of the scalar- and pseudoscalar-density fluctuations in a chiral effective model with a nonzero current quark mass. It is known that the soft mode associated with the second-order transition at the CP lies in the spacelike region of the scalar-density fluctuation. We find that the soft mode influences the quark spectrum significantly near the CP, resulting in the shift of the quasiquark peak. Effects of the composite stable pions on the quark spectrum near the CP are also discussed.

preprint2014arXiv

Higher order fluctuations of conserved charges in heavy ion collisions

We investigate the effect of the global charge conservation on the cumulants of conserved charges in relativistic heavy ion collisions in a finite rapidity window by studying the time evolution of cumulants in the hadronic medium. It is argued that the global charge conservation does not affect the experimental result of the net-electric charge fluctuation observed by ALICE, because of the finite diffusion distance of charged particles in the hadronic stage.

preprint2014arXiv

Measurement of thermodynamics using gradient flow

We analyze bulk thermodynamics and correlation functions of the energy-momentum tensor in pure Yang-Mills gauge theory using the energy-momentum tensor defined by the gradient flow and small flow time expansion. Our results on thermodynamic observables are consistent with those obtained by the conventional integral method. The analysis of the correlation function of total energy supports the energy conservation. It is also addressed that these analyses with gradient flow require less statistics compared with the previous methods. All these results suggest that the energy-momentum tensor can be successfully defined and observed on the lattice with moderate numerical costs with the gradient flow.

preprint2014arXiv

Possible generation of anomalously soft quark excitations at nonzero temperature: Nonhyperbolic dispersion of parapion and van Hove singularity

We study the quark spectrum at finite temperature near and above the pseudocritical temperature of the chiral phase transition incorporating the effects of the collective modes with the quantum number of the sigma (parasigma) and pion (parapion) in a chiral effective model with a nonzero current quark mass. Below the pion zero-binding temperature where the pionic modes are bound, the quark self-energy has van Hove singularity induced by the scattering of quarks with the composite bound pions with a nonhyperbolic dispersion curve. This singularity is found to cause a drastic change in the quark spectrum from that in the mean field picture near the pseudocritical temperature: The quark spectrum has an unexpected sharp peak at an energy considerably lower than the constituent quark mass, while the spectrum approaches the mean field one at high temperatures. We clarify that the emergence of this anomalous structure of the quark spectral function originates from the composite nature of the pionic modes with a non-Lorentz invariant dispersion relation in the medium at finite temperature.

preprint2013arXiv

Nonequilibrium time evolution of higher order cumulants of conserved charges and event-by-event analysis

We investigate the time evolution of higher order cumulants of conserved charges in a volume with the diffusion master equation. Applying the result to the diffusion of non-Gaussian fluctuations in the hadronic stage of relativistic heavy ion collisions, we show that the fourth-order cumulant of net-electric charge at LHC energy is suppressed compared with the recently observed second-order cumulant at ALICE, if the higher order cumulants at hadronization are suppressed compared with their values in the hadron phase in equilibrium. The significance of the experimental information on the rapidity window dependence of various cumulants in investigating the history of the dynamical evolution of the hot medium created in relativistic heavy ion collisions is emphasized.

preprint2013arXiv

On the effect of secondary protons on baryon and proton number cumulants in event-by-event analysis

We investigate the effects of secondary (knockout) protons, which constitute about 20% of the observed protons at STAR, on the higher order cumulants of proton and baryon numbers measured by event-by-event analyses in relativistic heavy ion collisions. We argue that the contribution of this background effect on the cumulants is expressed by a simple formula, and that hence their effects can be removed in the experimental analysis. It is discussed that this background effect has non-negligible contribution to recently observed proton number cumulants at STAR, especially the third-order one, and that the removal of this effect is crucial to investigate the thermodynamical properties of the primordial hot medium appropriately.

preprint2013arXiv

Shear viscosity to relaxation time ratio in SU(3) lattice gauge theory

We evaluate the ratio of the shear viscosity to the relaxation time of the shear flux above but near the critical temperature $T_c$ in SU(3) gauge theory on the lattice. The ratio is related to Kubo's canonical correlation of the energy-momentum tensor in Euclidean space with the relaxation time approximation and an appropriate regularization. Using this relation, the ratio is evaluated by direct measurements of the Euclidean observables on the lattice. We obtained the ratio with reasonable statistics for the range of temperature $1.3T_c \lesssim T \lesssim 4T_c$. We also found that the characteristic speed of the transverse plane wave in gluon media is almost constant, $v \simeq 0.5$, for $T \gtrsim 1.5T_c$, which is compatible with the causality in the second order dissipative hydrodynamics.

preprint2012arXiv

Relation between baryon number fluctuations and experimentally observed proton number fluctuations in relativistic heavy ion collisions

We explore the relation between proton and nucleon number fluctuations in the final state in relativistic heavy ion collisions. It is shown that the correlations between the isospins of nucleons in the final state are almost negligible over a wide range of collision energy. This leads to a factorization of the distribution function of the proton, neutron, and their antiparticles in the final state with binomial distribution functions. Using the factorization, we derive formulas to determine nucleon number cumulants, which are not direct experimental observables, from proton number fluctuations, which are experimentally observable in event-by-event analyses. With a simple treatment for strange baryons, the nucleon number cumulants are further promoted to the baryon number ones. Experimental determination of the baryon number cumulants makes it possible to compare various theoretical studies on them directly with experiments. Effects of nonzero isospin density on this formula are addressed quantitatively. It is shown that the effects are well suppressed over a wide energy range.

preprint2012arXiv

Revealing baryon number fluctuations from proton number fluctuations in relativistic heavy ion collisions

Baryon number cumulants are invaluable tools to diagnose the primordial stage of heavy ion collisions if they can be measured. In experiments, however, proton number cumulants have been measured as substitutes. In fact, proton number fluctuations are further modified in the hadron phase and different from those of baryon number. We show that the isospin distribution of nucleons at kinetic freezeout is binomial and factorized. This leads to formulas that express the baryon number cumulants sorely in terms of proron number fluctuations, which are experimentally observable.

preprint2012arXiv

Thermal mass and dispersion relations of quarks in the deconfined phase of quenched QCD

Temporal quark correlation functions are analyzed in quenched lattice QCD for two values of temperature above the critical temperature (Tc) for deconfinement, T=1.5Tc and 3Tc. A two-pole ansatz for the quark spectral function is used to determine the bare quark mass and the momentum dependence of excitation spectra on large lattices of size up to 128^3x16. The dependence of the quark correlator on these parameters as well as the finite volume dependence of the excitation energies are analyzed in detail in order to examine the reliability of our analysis. Our results suggest the existence of quasi-particle peaks in the quark spectrum. We furthermore find evidence that the dispersion relation of the plasmino mode has a minimum at non-zero momentum even in the non-perturbative region near Tc. We also elaborate on the enhancement of the quark correlator near the chiral limit which is observed at T=1.5Tc on about half of the gauge configurations. We attribute this to the presence of near zero-modes of the fermion matrix that are associated with non-trivial topology of the gauge configurations.

preprint2010arXiv

Spectral properties of quarks above Tc -- thermal mass, dispersion relation, and self-energy --

Spectral properties of quarks above the critical temperature for deconfinement are analyzed in quenched lattice QCD on lattices of size 128^3x16. We study quark spectral function in energy and momentum space, focusing on the values of the thermal mass and the dispersion relations of normal and plasmino modes at nonzero momentum, as well as their spatial volume dependence. Our numerical result suggests that the dispersion relation of the plasmino mode has a minimum at nonzero momentum even near the critical temperature. The quark self-energy is also analyzed by using the analyticy of the inverse propagator, which is found to be consistent with the spectral function estimated by the two-pole ansatz.

preprint2009arXiv

Evaluation of Israel-Stewart parameters in lattice gauge theory

Navier-Stokes equations are known as hydrodynamic equations which take account of effects of dissipations. There are, however, problems in the relativistic Navier-Stokes equations, i.e. the equations violate causality. Israel-Stewart equations, which evade the problems of Navier-Stokes equations by introducing new parameters, such as the relaxation times, have recently been used in describing the space-time evolution of the quark-gluon plasma produced in high energy heavy ion collisions. The viscosities and the relaxation times are related to each other by imposing entropy constraints on the system. According to Boltzmann-Einstein principle, the probability distribution of the fluctuation in the energy-momentum tensor is related to the entropy of the system. Applying this principle to the entropy in Israel-Stewart theory, one can obtain the ratios of the viscosities to the relaxation times. We evaluate the ratios of the viscosities to the relaxation times in SU(3) lattice gauge theory.

preprint2009arXiv

Quark propagator at finite temperature and finite momentum in quenched lattice QCD

We present an analysis of the quark spectral function above and below the critical temperature for deconfinement performed at zero and non-zero momentum in quenched lattice QCD using clover improved Wilson fermions in Landau gauge. It is found that the temporal quark correlation function in the deconfined phase near the critical temperature is well reproduced by a two-pole ansatz for the spectral function. This indicates that excitation modes of the quark field have small decay rates. The bare quark mass and momentum dependence of the spectral function is analyzed with this ansatz. In the chiral limit we find that the quark spectral function has two collective modes corresponding to the normal and plasmino excitations in the high temperature limit. Over a rather wide temperature range in the deconfined phase the pole mass of these modes at zero momentum, which corresponds to the thermal mass of the quark, is approximately proportional to temperature. With increasing bare quark masses the plasmino mode gradually disappears and the spectral function is dominated by a single pole. We also discuss quasi-particle properties of heavy quarks in the deconfined phase. In the confined phase, it is found that the pole ansatz for the spectral function fails completely.

preprint2009arXiv

Third moments of conserved charges as probes of QCD phase structure

The third moments of conserved charges, the baryon and electric charge numbers, and energy, as well as their mixed moments, carry more information on the state around the QCD phase boundary than previously proposed fluctuation observables and higher order moments. In particular, their signs give plenty of information on the location of the state created in relativistic heavy ion collisions in the temperature and baryon chemical potential plane. We demonstrate this with an effective model.

preprint2009arXiv

Third moments of conserved charges in QCD phase diagram

We point out that the third moments of conserved charges, the baryon and electric charge numbers, and energy, as well as their mixed moments, change their signs around the QCD phase boundary in the temperature and baryon chemical potential plane. These signs can be measured in relativistic heavy ion collisions, and will give clear information on the phase structure of QCD and the state of the system in the early stage of relativistic heavy ion collisions. The behaviors of these moments on the temperature axis and at small quark chemical potential can be analyzed in lattice QCD simulations. We emphasize that the third moments obtained on the lattice, together with the experimental results, will provide a deep understanding about the QCD phase diagram and the location of the state created in heavy ion collisions.

preprint2007arXiv

Chiral transition and mesonic excitations for quarks with thermal masses

We study the effect of a thermal quark mass, m_T, on the chiral phase transition and mesonic excitations in the light quark sector at finite temperature in a simple chirally-symmetric model. We show that while nonzero m_T lowers the chiral condensate, the chiral transition remains of second order. It is argued that the mesonic excitations have large decay rate at energies below 2m_T, owing to the Landau damping of the quarks and the van Hove singularities of the collective modes.

preprint2007arXiv

Further Analysis of Excitations of Quarks at Finite Temperature -- Mass Effect and Pole Structure --

We calculate the spectral function of the massive quark at finite temperature (T) using a Yukawa model and show that the peak in the negative energy region among the three-peaks found in a previous work for the massless quark is largely suppressed. To explore the underlying mechanism of this behavior, we also investigate the pole structure of the retarded Green function of the quark. We will show the result only for the massless quark. We find the residues of the poles corresponding the three-peaks are all comparable at $T\sim m_b$. We also show that the multi-peak structure of the quark spectra is well described in the pole approximation which indicates that the quasi-paricle picture is valid in this $T$ region.

preprint2007arXiv

Spectrum of soft mode with thermal mass of quarks above critical temperature

We study effects of a thermal-quark mass and Landau damping of quarks on the chiral phase transition and its soft modes at finite temperature. For this purpose we employ a simple model with the quark propagator obtained in the hard thermal loop approximation. We show that the chiral phase transition is second order even if quark have a finite thermal mass, and the thermal mass suppresses the chiral condensate. We argue that mesonic spectra have a large width due to scattering between gluons and quarks, and the van Hove singularity at threshold.