Source author record

Teiji Kunihiro

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

54works
16topics
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

54 published item(s)

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

Entropy production in longitudinally expanding Yang-Mills field with use of Husimi function$-$semiclassical approximation

We investigate the possible thermalization process of the highly occupied and weakly coupled Yang-Mills fields expanding along the beam axis through an evaluation of the entropy, particle number, and pressure anisotropy. The time evolution of the system is calculated by solving the equation of motion for the Wigner function in the semiclassical approximation with initial conditions mimicking the glasma. For the evaluation of the entropy, we adopt the Husimi-Wehrl (HW) entropy, which is obtained by using the Husimi function, a positive semidefinite quantum distribution function given by smearing the Wigner function. By numerical calculations at $g=0.1$ and $0.2$, the entropy production is found to occur together with the particle creation in two distinct stages: In the first stage, the particle number and the entropy at low longitudinal momenta grow rapidly. In the second stage, the particle number and the entropy of higher longitudinal momentum modes show slower increase. The pressure anisotropy remains in our simulation and implies that the system is still out-of-equilibrium.

preprint2020arXiv

Shear viscosity of classical fields in scalar theory

We investigate the shear viscosity of massless classical scalar fields in the $ϕ^4$ theory on a lattice by using the Green-Kubo formula. Based on the scaling property of the classical field, the shear viscosity is represented using a scaling function. Equilibrium expectation value of the time-correlation function of the energy-momentum tensor is evaluated as the ensemble average of the classical field configurations, whose time evolution is obtained by solving the classical equation of motion starting from the initial condition in thermal equilibrium. It is found that there are two distinct damping time scales in the time-correlation function, which is found to show damped oscillation behavior in the early stage around a slow monotonous decay with an exponential form, and the slow decay part is found to dominate the shear viscosity in the massless classical field theory. This kind of slow decay is also known to exist in the molecular dynamics simulation, then it may be a generic feature of dense matter.

preprint2020arXiv

Shear viscosity of classical Yang-Mills field

We investigate the shear viscosity $η$ of the classical Yang-Mills (CYM) field on a lattice by using the Green-Kubo formula, where the shear viscosity is calculated from the time-correlation function of the energy-momentum tensor in equilibrium. Dependence of the shear viscosity $η(g,T)$ on the coupling $g$ and temperature $T$ is represented by a scaling function $f_η(g^2T)$ as $η(g,T)=Tf_η(g^2T)$ due to the scaling-invariant property of the CYM. The explicit functional form of $f_η(g^2T)$ is successfully determined from the calculated shear viscosity: It turns out that $η(g,T)$ of the CYM field is proportional to $1/g^{1.10-1.88}$ at weak coupling, which is a weaker dependence on $g$ than that in the leading-order perturbation theory but consistent with that of the "anomalous viscosity" $η\propto 1/g^{1.5}$ under the strong disordered field. The obtained shear viscosity is also found to be roughly consistent with that estimated through the analysis of the anisotropy of the pressure of the CYM dynamics in the expanding geometry with recourse to a hydrodynamic equation.

preprint2016arXiv

Effective field theory and the scattering process for magnons in the ferromagnet, antiferromagnet, and ferrimagnet

We discuss that a low-energy effective Lagrangian relying on SO(3) $\rightarrow$ SO(2) is applicable for a ferrimagnet as well as a ferromagnet and an antiferromagnet. The analysis of the particle states shows that there exist not only massless modes with the dispersion relations $ω\propto |\bm{k}|,\, |\bm{k}|^2$, i.e., the so-called type-I and type-II Nambu-Goldstone modes, respectively, but also gapped modes with $ω\propto m^2+|\bm{k}|^2$. We clarify how the coefficients of the terms with one time derivative and those with two time derivatives in the effective Lagrangian determine the order parameters specifying whether the system is in a ferromagnetic, antiferromagnetic or ferrimagnetic state; we stress that the gapped mode related to the spontaneous symmetry breaking appears only in the ferrimagnetic system and not in the ferromagnetic and antiferromagnetic systems. We also establish the power counting scheme and calculate the scattering amplitudes and thereby the scattering lengths between the two Nambu-Goldstone bosons. We show that the scattering length of the gapped mode is finite and proportional to the gap. This characteristic property of the gapped NG mode can be used to discriminate it from gapped excitations which originate in other mechanisms. Finally, we study the effects of the explicit symmetry breaking that are given by an external magnetic field and a single-ion anisotropy, and show that the external magnetic fields do not have any effects on the scattering amplitudes in all the spin systems as was known for the ferromagnet system. In contrast, the anisotropy does affect the scattering amplitudes, the phase shift, and the scattering length except for spin 1/2 systems. This result supports the possibility of the Efimov effect in spin systems discussed in previous studies.

preprint2016arXiv

Entropy production from chaoticity in Yang-Mills field theory with use of the Husimi function

We investigate possible entropy production in Yang-Mills (YM) field theory by using a quantum distribution function called Husimi function $f_{\rm H}(A, E, t)$ for YM field, which is given by a coarse graining of Wigner function and non-negative. We calculate the Husimi-Wehrl (HW) entropy $S_{\rm HW}(t)=-{\rm Tr}f_H \log f_H$ defined as an integral over the phase-space, for which two adaptations of the test-particle method are used combined with Monte-Carlo method. We utilize the semiclassical approximation to obtain the time evolution of the distribution functions of the YM field, which is known to show a chaotic behavior in the classical limit. We also make a simplification of the multi-dimensional phase-space integrals by making a product ansatz for the Husimi function, which is found to give a 10-20 per cent over estimate of the HW entropy for a quantum system with a few degrees of freedom. We show that the quantum YM theory does exhibit the entropy production, and that the entropy production rate agrees with the sum of positive Lyapunov exponents or the Kolmogorov-Sinai entropy, suggesting that the chaoticity of the classical YM field causes the entropy production in the quantum YM theory.

preprint2016arXiv

Functional renormalization group analysis of the soft mode at the QCD critical point

We make an intensive investigation of the soft mode at the quantum chromodynamics (QCD) critical point on the basis of the functional renormalization group (FRG) method in the local potential approximation. We calculate the spectral functions $ρ_{σ, π}(ω,\, p)$ in the scalar ($σ$) and pseudoscalar ($π$) channels beyond the random phase approximation in the quark--meson model. At finite baryon chemical potential $μ$ with a finite quark mass, the baryon-number fluctuation is coupled to the scalar channel and the spectral function in the $σ$ channel has a support not only in the time-like ($ω\,>\,p$) but also in the space-like ($ω\,<\, p$) regions, which correspond to the mesonic and the particle--hole phonon excitations, respectively. We find that the energy of the peak position of the latter becomes vanishingly small with the height being enhanced as the system approaches the QCD critical point, which is a manifestation of the fact that the phonon mode is the {\em soft mode} associated with the second-order transition at the QCD critical point, as has been suggested by some authors. Moreover, our extensive calculation of the spectral function in the $(ω, p)$ plane enables us to see that the mesonic and phonon modes have the respective definite dispersion relations $ω_{σ.{\rm ph}}(p)$, and it turns out that $ω_σ(p)$ crosses the light-cone line into the space-like region, and then eventually merges into the phonon mode as the system approaches the critical point more closely. This implies that the sigma-mesonic mode also becomes soft at the critical point. We also provide numerical stability conditions that are necessary for obtaining the accurate effective potential from the flow equation.

preprint2016arXiv

Mesoscopic dynamics of fermionic cold atoms ---Quantitative analysis of transport coefficients and relaxation times---

We give a quantitative analysis of the dynamical properties of fermionic cold atomic gases in normal phase, such as the shear viscosity, heat conductivity, and viscous relaxation times, using the novel microscopic expressions derived by the renormalization group (RG) method, where the Boltzmann equation is faithfully solved to extract the hydrodynamics without recourse to any ansatz. In particular, we examine the quantum statistical effects, temperature dependence, and scattering-length dependence of the transport coefficients and the viscous relaxation times. The numerical calculation shows that the relation $τ_π=η/P$, which is derived in the relaxation-time approximation (RTA) and is used in most of the literature, turns out to be satisfied quite well, while the similar relation for the viscous relaxation time $τ_J$ of the heat conductivity is satisfied only approximately with a considerable error.

preprint2016arXiv

Novel picture of the soft modes at the QCD critical point based on the FRG method

We investigate the soft mode at the QCD critical point (CP) on the basis of the functional renormalization group. We calculate the spectral functions in the meson channels in the two-flavor quark--meson model. Our result shows that the energy of the peak position of the particle--hole mode in the sigma channel becomes vanishingly small as the system approaches the QCD CP, which is a manifestation of the softening of the phonon mode. We also extract the dispersion curves of the mesonic and the phonon mode, a hydrodynamic mode which leads to a finding that the dispersion curve of the sigma-mesonic mode crosses the light-cone into the space-like momentum region, and then eventually merges into the phonon mode as the system approaches further close to the CP. This may suggest that the sigma-mesonic mode forms the soft mode together with the hydrodynamic mode at the CP.

preprint2016arXiv

Possible relevance of the softening of the sigma meson to $η$ decay into 3$π$ in the nuclear medium

We investigate the role of the softening of the scalar-isoscalar (sigma) meson in the $η\rightarrowπ^+π^-π^0$ and $3π^0$ decay widths in the symmetric nuclear medium using a linear sigma model. Our calculation shows that these decay widths in the nuclear medium increase by up to a factor of four to ten compared with those in the free space mainly depending on the mass of the sigma meson in the free space, which is an input parameter of the model. The enhancements are considerable even at a half of the normal nuclear density. Thus, the $η$ decay into $3π$ could be a possible new probe for the chiral restoration in the nuclear medium. We find that the density dependence of the $η\rightarrow3π^0$ decay is moderate in comparison with that of $η\rightarrowπ^+π^-π^0$, although the former width is greater than the latter one at a given density. This is because the softening of the sigma meson causes cancellation of the terms appearing from the Bose symmetry in the $η\rightarrow3π^0$ decay. The difference between the density dependences should be helpful for experimental confirmation of the findings of the present study.

preprint2016arXiv

Second-order hydrodynamics for fermionic cold atoms: Detailed analysis of transport coefficients and relaxation times

We give a detailed derivation of the second-order (local) hydrodynamics for Boltzmann equation with an external force by using the renormalization group method. In this method, we solve the Boltzmann equation faithfully to extract the hydrodynamics without recourse to any ansatz. Our method leads to microscopic expressions of not only all the transport coefficients that are of the same form as those in Chapman-Enskog method but also those of the viscous relaxation times $τ_i$ that admit physically natural interpretations. As an example, we apply our microscopic expressions to calculate the transport coefficients and the relaxation times of the cold fermionic atoms in a quantitative way, where the transition probability in the collision term is given explicitly in terms of the $s$-wave scattering length $a_s$. We thereby discuss the quantum statistical effects, temperature dependence, and scattering-length dependence of the first-order transport coefficients and the viscous relaxation times: It is shown that as the temperature is lowered, the transport coefficients and the relaxation times increase rapidly because Pauli principle acts effectively. On the other hand, as $a_s$ is increased, these quantities decrease and become vanishingly small at unitarity because of the strong coupling. The numerical calculation shows that the relation $τ_π=η/P$, which is derived in the relaxation-time approximation and used in most of literature without almost any foundation, turns out to be satisfied quite well, while the similar relation for the relaxation time $τ_J$ of the heat conductivity is satisfied only approximately with a considerable error.

preprint2016arXiv

Spectral functions in functional renormalization group approach -- analysis of the collective soft modes at the QCD critical point --

We first review the method to calculate the spectral functions in the functional renormalization group (FRG) approach, which has been recently developed. We also provide the numerical stability conditions given by the present authors for a generic nonlinear evolution equation that are necessary for obtaining the accurate effective potential from the flow equation in the FRG. As an interesting example, we report the recent calculation of the spectral functions of the mesonic and particle-hole excitations using a chiral effective model of Quantum Chromodynamics (QCD); we extract the dispersion relations from them and try to reveal the nature of the soft modes at the QCD critical point (CP) where the phase transition is second order. Our result shows that a clear development and the softening of the phonon mode in the space-like region as the system approaches the CP; furthermore it turns out that the sigma mesonic mode once in the time-like region gets to merge with the phonon mode in the close vicinity of the CP, implying a novel possibility about the nature of the soft mode of the QCD CP.

preprint2015arXiv

Causal hydrodynamics from kinetic theory by doublet scheme in renormalization-group method

We develop a general framework in the renormalization-group (RG) method for extracting a mesoscopic dynamics from an evolution equation by incorporating some excited (fast) modes as additional components to the invariant manifold spanned by zero modes. We call this framework the doublet scheme. The validity of the doublet scheme is first tested and demonstrated by taking the Lorenz model as a simple three-dimensional dynamical system; it is shown that the two-dimensional reduced dynamics on the attractive manifold composed of the would-be zero and a fast modes are successfully obtained in a natural way. We then apply the doublet scheme to construct causal hydrodynamics as a mesoscopic dynamics of kinetic theory, i.e., the Boltzmann equation, in a systematic manner with no ad-hoc assumption. It is found that our equation has the same form as thirteen-moment causal hydrodynamic equation, but the microscopic formulae of the transport coefficients and relaxation times are different. In fact, in contrast to the Grad equation, our equation leads to the same expressions for the transport coefficients as given by the Chapman-Enskog expansion method and suggests novel formulae of the relaxation times expressed in terms of relaxation functions which allow a natural physical interpretation of the relaxation times. Furthermore, our theory nicely gives the explicit forms of the distribution function and the thirteen hydrodynamic variables in terms of the linearized collision operator, which in turn clearly suggest the proper ansatz forms of them to be adopted in the method of moments.

preprint2015arXiv

Derivation of second-order relativistic hydrodynamics for reactive multi-component systems

We derive the second-order hydrodynamic equation for reactive multi-component systems from the relativistic Boltzmann equation. In the reactive system, particles can change their species under the restriction of the imposed conservation laws during the collision process. Our derivation is based on the renormalization group (RG) method, in which the Boltzmann equation is solved in an organized perturbation method as faithfully as possible and possible secular terms are resummed away by a suitable setting of the initial value of the distribution function. The microscopic formulae of the relaxation times and the lengths are explicitly given as well as those of the transport coefficients for the reactive multi-component system. The resultant hydrodynamic equation with these formulae has nice properties that it satisfies the positivity of the entropy production rate and the Onsager's reciprocal theorem, which ensure the validity of our derivation.

preprint2015arXiv

Entropy production in quantum Yang-Mills mechanics in semi-classical approximation

We discuss thermalization of isolated quantum systems by using the Husimi-Wehrl entropy evaluated in the semiclassical treatment. The Husimi-Wehrl entropy is the Wehrl entropy obtained by using the Husimi function for the phase space distribution. The time evolution of the Husimi function is given by smearing the Wigner function, whose time evolution is obtained in the semiclassical approximation. We show the efficiency and usefullness of this semiclassical treatment in describing entropy production of a couple of quantum mechanical systems, whose classical counter systems are known to be chaotic. We propose two methods to evaluate the time evolution of the Husimi-Wehrl entropy, the test-particle method and the two-step Monte-Carlo method. We demonstrate the characteristics of the two methods by numerical calculations, and show that the simultaneous application of the two methods ensures the reliability of the results of the Husimi-Wehrl entropy at a given time.

preprint2015arXiv

Lattice QCD study of four-quark components of the isosinglet scalar mesons: Significance of disconnected diagrams

We study the possible significance of four-quark states in the iso-singlet scalar mesons ($J^{PC}=0^{++}$, $I=0$) by performing two-flavor full lattice QCD simulations on an $8^3 \times 16$ lattice using the improved gauge action and the clover-improved Wilson quark action. In particular, we evaluate the propagators of molecular and tetra-quark operators together with singly disconnected diagrams. In the computation of the singly disconnected diagrams we employ the $Z_2$-noise method with the truncated eigenmode approach. We show that the quark loops given by the disconnected diagrams play an essential role in propagators of tetraquark and molecular operators.

preprint2015arXiv

Parametric instability of classical Yang-Mills fields in a color magnetic background

We investigate instabilities of classical Yang-Mills fields in a time-dependent spatially homogeneous color magnetic background field in a non-expanding geometry for elucidating the earliest stage dynamics of ultra-relativistic heavy-ion collisions. The background gauge field configuration considered in this article is spatially homogeneous and temporally periodic, and is alluded by Berges-Scheffler-Schlichting-Sexty (BSSS). We discuss the whole structure of instability bands of fluctuations around the BSSS background gauge field on the basis of Floquet theory, which enables us to discuss the stability in a systematic way. We find various instability bands on the $(p_z, p_T)$-plane. These instability bands are caused by parametric resonance despite the fact that the momentum dependence of the growth rate for $|\mathbf{p}| \leq \sqrt{B}$ is similar to a Nielsen-Olesen instability. Moreover, some of instability bands are found to emerge not only in the low momentum but also in the high momentum region; typically of the order of the saturation momentum as $|\mathbf{p}| \sim \sqrt{B} \sim Q_{\rm s}$.

preprint2015arXiv

Parametric Instability of Classical Yang-Mills Fields in an Expanding Geometry

We investigate the instability of classical Yang-Mills field in an expanding geometry under a color magnetic background field within the linear regime. We consider homogeneous, boost-invariant and time-dependent color magnetic fields simulating the glasma configuration. We introduce the conformal coordinates which enable us to map an expanding problem approximately into a nonexpanding problem. We find that the fluctuations with finite longitudinal momenta can grow exponentially due to parametric instability. Fluctuations with finite transverse momenta can also show parametric instability, but their momenta are restricted to be small. The most unstable modes start to grow exponentially in the early stage of the dynamics and they may affect the thermalization in heavy-ion collisions.

preprint2015arXiv

Possible medium effects on η-π^0 mixing and η\rightarrow π^0 π^+ π^- decay in the asymmetric nuclear matter

It is known that the decay process η\rightarrowπ^0π^+π^- in free space is possible due to the isospin-symmetry breaking in Quantum Chromodynamics (QCD), i.e., the small mass difference between u and d quarks: The small width is intimately related with the mixing property between ηand π^0 mesons for which the chiral anomaly play a role. In asymmetric nuclear matter such as heavy nuclei, isospin-symmetry breaking is large and it is expected that the mixing property of the mesons changes significantly and the above-mentioned decay width of the ηcreated in such a medium may be enhanced. This is an intriguing possibility to revealing a medium effect on a hadron in nuclei. We apply the in-medium chiral perturbation theory to estimate the modification of the η-π^0 mixing angle and the partial width of ηπ^0π^+π^- in asymmetric nuclear matter as a function of the isospin asymmetry. We find that these quantities can be greatly enhanced in neutron-rich matter, which should be deteactable in experiment.

preprint2015arXiv

Relativistic Causal Hydrodynamics Derived from Boltzmann Equation: a novel reduction theoretical approach

We derive the second-order hydrodynamic equation and the microscopic formulae of the relaxation times as well as the transport coefficients systematically from the relativistic Boltzmann equation. Our derivation is based on a novel development of the renormalization-group method, a powerful reduction theory of dynamical systems, which has been applied successfully to derive the non-relativistic second-order hydrodynamic equation Our theory nicely gives a compact expression of the deviation of the distribution function in terms of the linearized collision operator, which is different from those used as an ansatz in the conventional fourteen-moment method. It is confirmed that the resultant microscopic expressions of the transport coefficients coincide with those derived in the Chapman-Enskog expansion method. Furthermore, we show that the microscopic expressions of the relaxation times have natural and physically plausible forms. We prove that the propagating velocities of the fluctuations of the hydrodynamical variables do not exceed the light velocity, and hence our seconder-order equation ensures the desired causality. It is also confirmed that the equilibrium state is stable for any perturbation described by our equation.

preprint2015arXiv

Study of entropy production in Yang-Mills theory with use of Husimi function

Understanding the thermalization process in a pure quantum system is a challenge in theoretical physics. In this work, we explore possible thermalization mechanism in Yang-Mills(YM) theory by using a positive semi-definite quantum distribution function called a Husimi function which is given by a coarse graining of the Wigner function within the minimum uncertainty. Then entropy is defined in terms of the Husimi function, which is called the Husimi-Wehrl(HW) entropy. We propose two numerical methods to calculate the HW entropy. We find that it is feasible to apply the semi-classical approximation with the use of classical YM equation. It should be noted that the semi-classical approximation is valid in the systems of physical interest including the early stage of heavy-ion collisions. Using a product ansatz for the Husimi function, which is confirmed to reproduce the HW entropy within 20% error (overestimate) for a few-body quantum system, we succeed in a numerical evaluation of HW entropy of YM fields and show that it surely has a finite value and increases in time.

preprint2015arXiv

The $η$ decay into 3$π$ in asymmetric nuclear medium

We explore how the $η$-$π^0$ mixing angle and the $η$ meson decay into $π^{+}π^{-}π^0$ and 3$π^{0}$ are modified in the nuclear medium on the basis of the in-medium chiral effective theory with the isospin asymmetry $α$ varied, where $α\equiv δρ/ρ$ with $δρ=ρ_n-ρ_p$ and $ρ=ρ_n+ρ_p$. We find that the larger the isospin asymmetry $δρ$ and smaller the total density $ρ$, the more enhanced the mixing angle. We show that the decay width in the nuclear medium has an additional density dependence which cannot be renormalized into that of the mixing angle: The additional term originates from the vertex proportional to a low energy constant $c_1$, which only comes into play in the nuclear medium but not in the free space. It turns out that the resultant density effect on the decay widths overwhelms that coming from the isospin asymmetry, and the higher the $ρ$, the more enhanced the decay widths; the width for the $π^{+}π^{-}π^0$ decay is enhanced with a factor two to three at the normal density $ρ_0$ with a minor increase due to $δρ$, while that for the 3$π^0$ decay shows only a small increase of around 10 percent even at $ρ_0$. We mention the possible relevance of the partial restoration of chiral symmetry to the unexpected density effect on the decay widths in the nuclear 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

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.

preprint2014arXiv

Time evolution of gluon coherent state and its von Neumann entropy in heavy-ion collisions

We propose a new prescription for evaluating a von Neumann entropy in the initial stage of high-energy heavy-ion collisions utilizing the time evolution of classical Yang-Mills (CYM) field: The von Neumann entropy is computed for the quantum coherent states constructed so as to give the classical gluon fields as the expectation values. The entropy is to be liberated when the complete decoherence is achieved. As a demonstration, the time evolution of the CYM dynamics is solved with an initial condition which mimics the Glasma state, though in a non-expanding geometry; the Glasma state is characterized by the longitudinal color-electric and -magnetic fields with gluon fields' fluctuations around it. We find that the initial longitudinal fluctuations of the fields play essential roles for the entropy production in two ways: First, the field fluctuations at $t=0$ themselves act as a source of the von Neumann entropy prepared before the time evolution. Second, the initial fluctuations triggers field instabilities, and hence the larger the strength of them, the more the entropy production at later time.

preprint2013arXiv

Entropy production in classical Yang-Mills theory from Glasma initial conditions

We study the thermalization process in classical Yang-Mills (CYM) field theory starting from noisy glasma-like initial conditions by investigating the initial-value sensitivity of trajectories. Kunihiro et al. linked entropy generation to the Kolmogorov-Sinai entropy, which gives the entropy production rate in classical chaotic systems, calculated numerically for CYM fields starting from purely random initial field configurations. In contrast, we here study glasma-like initial conditions. For small random fluctuations we obtain qualitatively similar results while no entropy increase is observed when such fluctuations are absent. We analyze the intermediate time Lyapunov spectrum for several time windows and calculate the Kolmogorov-Sinai entropy. We find a large number of positive Lyapunov exponents at the early stages of time evolution. Also for later times their number is a sizeable fraction of the total number of degrees of freedom. The spectrum of positive Lyapunov exponents at first changes rapidly, but then stabilizes, indicating that the dynamics of the gauge fields approaches a steady state. Thus we conclude that also for glasma-like initial conditions a significant amount of entropy is produced by classical gluon field dynamics.

preprint2013arXiv

Neutrino spectral density at electroweak scale temperature

Motivated by the scenario of resonant leptogenesis in which lepton number creation in the electroweak-scale is relevant, we investigate the spectral properties and possible collective nature of the standard model neutrinos at electroweak scale temperature (T). We adopt the R_xi gauge fixing, which includes the unitary gauge as a limiting case, and allows us to study the broken as well as the restored phases of the gauge symmetry in a unified way. We show that the spectral density of the neutrino has a three-peak structure in the low-momentum region due to the scattering with the thermally excited particles (i.e., Landau damping) when T becomes comparable to the weak-boson masses in the plasma. The three peaks are identified with a novel ultrasoft mode, the usual quasiparticle, and antiplasmino modes. Varying the gauge-fixing parameter, we show that the three-peak structure appears independently of the gauge fixing and thus has a physical significance. We discuss possible implications of the neutrino spectral density obtained in the present work on particle cosmology, in particular in the context of resonant leptogenesis.

preprint2012arXiv

Derivation of relativistic hydrodynamic equations consistent with relativistic Boltzmann equation by renormalization-group method

We review our work on the application of the renormalization-group method to obtain first- and second-order relativistic hydrodynamics of the relativistic Boltzmann equation (RBE) as a dynamical system, with some corrections and new unpublished results. For the first-order equation, we explicitly obtain the distribution function in the asymptotic regime as the invariant manifold of the dynamical system, which turns out to be nothing but the matching condition defining the energy frame, i.e., the Landau-Lifshitz one. It is argued that the frame on which the flow of the relativistic hydrodynamic equation is defined must be the energy frame, if the dynamics should be consistent with the underlying RBE. A sketch is also given for derivation of the second-order hydrodynamic equation, i.e., extended thermodynamics, which is accomplished by extending the invariant manifold so that it is spanned by excited modes as well as the zero modes (hydrodynamic modes) of the linearized collision operator. On the basis of thus constructed resummed distribution function, we propose a novel ansatz for the functional form to be used in Grad moment method; it is shown that our theory gives the same expressions for the transport coefficients as those given in the Chapman-Enskog theory as well as the novel expressions for the relaxation times and lengths allowing natural interpretation.

preprint2012arXiv

New forms of non-relativistic and relativistic hydrodynamic equations as derived by the renormalization-group method - possible functional ansatz in the moment method consistent with Chapman-Enskog theory -

After a brief account of the derivation of the first-order relativistic hydrodynamic equation as a construction of the invariant manifold of relativistic Boltzmann equation, we give a sketch of derivation of the second-order hydrodynamic equation (extended thermodynamics) both in the nonrelativistic and relativistic cases. We show that the resultant equation suggests a novel ansatz for the functional form to be used in Grad moment method, which turns out to give the same expressions for the transport coefficients as those given in the Chapman-Enskog theory as well as the novel expressions for the relaxation times and lenghts allowing natural physical interpretaion.

preprint2012arXiv

Uniqueness of Landau-Lifshitz Energy Frame in Relativistic Dissipative Hydrodynamics

We show that the relativistic dissipative hydrodynamic equation derived from the relativistic Boltzmann equation by the renormalization-group method uniquely leads to the one in the energy frame proposed by Landau and Lifshitz, provided that the macroscopic-frame vector, which defines the local rest frame of the fluid velocity, is independent of the momenta of constituent particles, as it should. We argue that the relativistic hydrodynamic equations for viscous fluids must be defined on the energy frame if it is consistent with the underlying relativistic kinetic equation.

preprint2011arXiv

Axial Anomaly, Mismatched Fermi Surfaces and Vector Interaction in Dense Neutral Quark Matter

We report on effects of the q-$\bar{\rm q}$ vector interaction and/or the U(1)$_A$-anomaly-induced chiral-diquark coupling on the charge-neutral quark matter in $β$-equilibrium. We show that when the vector coupling is absent, there can appear a cross-over region sandwiched by two critical points in the intermediate temperature ($T$) region, while the phase transition in the low-$T$ region including zero temperature keeps being first order until the strength of the anomaly term is increased to have a critical value. On the other hand, when the vector coupling is also present, there appears a crossover region in the low-$T$ area including zero temperature with a new critical point, as was first demonstrated by Kitazawa et al and the present authors without and with the charge-neutral condition, respectively. We remark that the possible chromomagnetic instability is suppressed and can be even completely absent owing to the enhanced diquark coupling due to the anomaly term and/or by the vector interaction.

preprint2011arXiv

First-principle Derivation of Stable First-order Generic-frame Relativistic Dissipative Hydrodynamic Equations from Kinetic Theory by Renormalization-group Method

We derive first-order relativistic dissipative hydrodynamic equations (RDHEs) from relativistic Boltzmann equation (RBE) on the basis of the renormalization-group (RG) method. We introduce a macroscopic-frame vector (MFV) to specify the local rest frame (LRF) on which the macroscopic dynamics is described. The five hydrodynamic modes are identified with the same number of the zero modes of the linearized collision operator, i.e., the collision invariants. After defining the inner product in the function space spanned by the distribution function, the higher-order terms, which give rise to the dissipative effects, are constructed so that they are orthogonal to the zero modes in terms of the inner product: Here, any ansatz's, such as the so-called conditions of fit used in the standard methods in an ad-hoc way, are not necessary. We elucidate that the Burnett term dose not affect the RDHEs owing to the very nature of the hydrodynamic modes as the zero modes. Applying the RG equation, we obtain the RDHE in a generic LRF specified by the MFV, as the coarse-grained and covariant equation. Our generic RDHE reduces to RDHEs in various LRFs, including the energy and particle LRFs with a choice of the MFV. We find that our RDHE in the energy LRF coincides with that of Landau and Lifshitz, while the derived RDHE in the particle LRF is slightly different from that of Eckart, owing to the presence of the dissipative internal energy. We prove that the Eckart equation can not be compatible with the underlying RBE. The proof is made on the basis of the observation that the orthogonality condition to the zero modes coincides with the ansatz's posed on the dissipative parts of the energy-momentum tensor and the particle current in the phenomenological RDHEs. We also present an analytic proof that all of our RDHEs have a stable equilibrium state owing to the positive definiteness of the inner product.

preprint2011arXiv

Multiple critical point structure for chiral phase transition induced by charge neutrality and vector interaction

The combined effect of the repulsive vector interaction and the positive electric chemical potential on the chiral phase transition is investigated by considering neutral color superconductivity. Under the charge-neutrality constraint, the chiral condensate, diquark condensate and quark number densities are obtained in two-plus-one-flavor Nambu-Jona-Lasinio model with the so called Kobayashi-Maskawa-'t Hooft term. We demonstrate that multiple chiral critical-point structures always exist in the Nambu-Jona-Lasinio model within the self-consistent mean-field approximation, and that the number of chiral critical points can vary from zero to four, which is dependent on the magnitudes of vector interaction and the diquark coupling.

preprint2011arXiv

Renormalized Linear Kinetic Theory as Derived from Quantum Field Theory - a novel diagrammatic method for computing transport coefficients -

We propose a novel diagrammatic method for computing transport coefficients in relativistic quantum field theory. The self-consistent equation for summing the diagrams with pinch singularities has a form of a linearized kinetic equation as usual, but our formalism enables us to incorporate higher-order corrections of the coupling systematically for the first time. Furthermore, it is clarified that the higher-order corrections are nicely summarized into that of the vertex function, spectral function, and collision term. We identify the diagrams up to the next-to-next-leading order corrections in the weak coupling expansion of phi^4 theory, which is a difficult task in kinetic approaches and other diagrammatic methods.

preprint2011arXiv

Roles of axial anomaly on neutral quark matter with color superconducting phase

We investigate effects of the axial anomaly term with a chiral-diquark coupling on the phase diagram within a two-plus-one-flavor Nambu-Jona-Lasinio (NJL) model under the charge-neutrality and $β$-equilibrium constraints. We find that when such constraints are imposed, the new anomaly term plays a quite similar role as the vector interaction does on the phase diagram, which the present authors clarified in a previous work. Thus, there appear several types of phase structures with multiple critical points at low temperature $T$, although the phase diagrams with intermediate-$T$ critical point(s) are never realized without these constraints even within the same model Lagrangian. This drastic change is attributed to an enhanced interplay between the chiral and diquark condensates due to the anomaly term at finite temperature; the u-d diquark coupling is strengthened by the relatively large chiral condensate of the strange quark through the anomaly term, which in turn definitely leads to the abnormal behavior of the diquark condensate at finite $T$, inherent to the asymmetric quark matter. We note that the critical point from which the crossover region extends to zero temperature appears only when the strength of the vector interaction is larger than a critical value. We also show that the chromomagnetic instability of the neutral asymmetric homogenous two-flavor color superconducting(2CSC) phase is suppressed and can be even completely cured by the enhanced diquark coupling due to the anomaly term and/or by the vector interaction.

preprint2011arXiv

Ultrasoft Fermionic Mode in Yukawa Theory at High Temperature

We explore whether an ultrasoft fermionic mode exists at extremely high temperature in Yukawa theory with massless fermion (coupling constant is g). We find that the fermion propagator has a pole at ω= \pm p/3-iζ, for ultrasoft momentum p, where ζ\sim g^4T ln 1/g, and the residue is Z \sim g^2. It is shown that one needs to take into account the asymptotic masses and the damping rate of hard particles to get a sensible result for such an ultrasoft fermionic mode; possible vertex correction turns out unnecessary for the scalar coupling in contrast to the gauge coupling.

preprint2011arXiv

Ultrasoft Fermionic Modes at High Temperature

A possible collective fermionic excitation in the ultrasoft energy-momentum region is examined in Yukawa model with scalar coupling and quantum electrodynamics (QED) with g being coupling constant at extremely high temperature T where the fermion mass is negligible. We analytically sum up the ladder diagrams for the vertex correction in the leading order in QED, which is not necessary in the Yukawa model, and find that the fermion pole exists at ω= \pm p/3-iζwith ultrasoft momentum p both for the Yukawa model and QED; ζis the sum of the damping rates of fermion and boson with hard momenta. We also obtain the expression of the residue of the pole, which is as small as of order g^2. We show that the fermion propagator and the vertex function satisfy the Ward-Takahashi identity in QED. Thus we establish the existence of an ultrasoft fermionic mode at extremely high temperature, which was originally called phonino and was suggested in the context of supersymmetry and its breaking at finite T. We discuss the possible origin of such an ultrasoft fermionic mode without recourse to supersymmetry. The case of QCD is briefly mentioned.

preprint2010arXiv

Nature of the sigma meson as revealed by its softening process

The pi pi scattering is studied in two-flavor chiral models with a finite pion mass to investigate the nature of the sigma meson which is observed as the lowest scalar-isoscalar resonance. We compare several models with different origins of the sigma meson, such as the chiral partner of the pion and the dynamically generated pi pi molecule. We find that the dynamically generated sigma meson exhibits a novel pattern of the threshold enhancement reflecting the s-wave nature of the resonance, which is qualitatively different from the softening of the chiral partner introduced as a bare field. This behavior around the threshold energy region is universal as far as the bare sigma pole stays away from the threshold throughout the symmetry restoration process. On the other hand, for m_pi=0, the dynamically generated sigma behaves similarly to the chiral partner in the symmetry restoration limit, implying the possibility of the dynamically generated chiral partner.

preprint2010arXiv

Second-order Relativistic Hydrodynamic Equations for Viscous Systems; how does the dissipation affect the internal energy?

We derive the second-order dissipative relativistic hydrodynamic equations in a generic frame with a continuous parameter from the relativistic Boltzmann equation. We present explicitly the relaxation terms in the energy and particle frames. Our results show that the viscosities are frame-independent but the relaxation times are generically frame-dependent. We confirm that the dissipative part of the energy-momentum tensor in the particle frame satisfies $δT^μ_μ= 0$ obtained for the first-order equation before, in contrast to the Eckart choice $u_μδT^{μν} u_ν= 0$ adopted as a matching condition in the literature. We emphasize that the new constraint $δT^μ_μ= 0$ can be compatible with the phenomenological derivation of hydrodynamics based on the second law of thermodynamics.

preprint2010arXiv

Spectral Function of Fermion Coupled with Massive Vector Boson at Finite Temperature in Gauge Invariant Formalism

We investigate spectral properties of a fermion coupled with a massive gauge boson with a mass m at finite temperature (T) in the perturbation theory. The massive gauge boson is introduced as a U(1) gauge boson in the Stueckelberg formalism with a gauge parameter α. We find that the fermion spectral function has a three-peak structure for T \sim m irrespective of the choice of the gauge parameter, while it tends to have one faint peak at the origin and two peaks corresponding to the normal fermion and anti-plasmino excitations familiar in QED in the hard thermal loop approximation for T \gg m. We show that our formalism successfully describe the fermion spectral function in the whole T region with the correct high-T limit except for the faint peak at the origin, although some care is needed for choice of the gauge parameter for T \gg m. We clarify that for T \sim m, the fermion pole is almost independent of the gauge parameter in the one-loop order, while for T \gg m, the one-loop analysis is valid only for α\ll 1/g where g is the fermion-boson coupling constant, implying that the one-loop analysis can not be valid for large gauge parameters as in the unitary gauge.

preprint2009arXiv

Toward Identifying the QCD Critical Point: attenuation of the sound mode around the critical point

Motivated by the fact that QCD critical point (CP) belongs to the same universality class as the liquid-gas transition, the dynamical density fluctuations around the CP is analyzed using relativistic fluid dynamics for a viscous fluid. It is shown that relativistic effects on the dynamical structure factor of the density fluctuation appears only in the width of the sound and thermal modes through a modification of the transport/kinetic coefficients. It is found that the mechanical density fluctuation which is enhanced by the relativistic effects are attenuated, whereas the entropy fluctuation in turn becomes more prominent around the QCD CP. This attenuation of the density fluctuation around the QCD CP may imply that if the suppression or even total disappearance of Mach cone is observed as the incident energy of the relativistic heavy-ion collisionsis lowered, it can be a signal of the existence of the QCD critical point.

preprint2009arXiv

Towards a Theory of Entropy Production in the Little and Big Bang

We propose a broadly applicable formalism for the description of coarse grained entropy production in quantum mechanical processes. Our formalism is based on the Husimi transform of the quantum state, which encodes the notion that information about any quantum state is limited by the experimental resolution. We show in two analytically tractable cases (the decay of an unstable vacuum state and reheating after cosmic inflation) that the growth rate of the Wehrl entropy associated with the Husimi function approaches the classical Kolmogorov-Sinai entropy. We also discuss various possible applications of our formalism, including the production of entropy in the early stages of a relativistic heavy ion collision.

preprint2009arXiv

Two-flavor lattice QCD study of the axial charges of N(1535) and N(1650)

We show the first lattice QCD results on the axial charges $g_A^{N^*N^*}$ of $N^*(1535)$ and $N^*(1650)$. The measurements are performed with two flavors of dynamical quarks employing the renormalization-group improved gauge action at $β$=1.95 and the mean-field improved clover quark action with the hopping parameters, $κ$=0.1375, 0.1390 and 0.1400. We construct 2$\times$2 correlation matrices and diagonalize them so that the signals of $N^*(1535)$ and $N^*(1650)$ are properly separated. Wraparound contributions in the correlator, which can be another source of signal contaminations, are eliminated by imposing the Dirichlet boundary condition in the temporal direction. The axial charge of $N^*(1535)$ is found to take small values as $g_A^{N^*N^*}\sim {\mathcal O}(0.1)$, whereas that of $N^*(1650)$ is approximately 0.5, which is almost independent of quark masses and consistent with the predictions by the naive nonrelativistic quark model.

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.

preprint2005arXiv

Application of the Renormalization-group Method to the Reduction of Transport Equations

We first give a comprehensive review of the renormalization group method for global and asymptotic analysis, putting an emphasis on the relevance to the classical theory of envelopes and on the importance of the existence of invariant manifolds of the dynamics under consideration. We clarify that an essential point of the method is to convert the problem from solving differential equations to obtaining suitable initial (or boundary) conditions:The RG equation determines the slow motion of the would-be integral constants in the unperturbative solution on the invariant manifold. The RG method is applied to derive the Navier-Stokes equation from the Boltzmann equation, as an example of the reduction of dynamics. We work out to obtain the transport coefficients in terms of the one-body distribution function.

preprint2003arXiv

Chiral Restoration and the Scalar and Vector Correlations in Hot and Dense Matter

First, it is pointed out that hadron/nuclear physics based on QCD should be regarded as ``condensed matter physics'' of the QCD vacuum. We indicate that phase shift analyses which respect chiral symmetry (ChS), analyticity and crossing symmetry of the scattering amplitude show the $σ$ meson pole in the s-channel in the low mass region as well as the $ρ$ meson pole in the $t$-channel in the \pipi scattering in the scalar channel. We review recent developments in exploring possible precursory phenomena of partial restoration of \chis in nuclear medium by examining the spectral function in the scalar and the vector channels. We emphasize that the wave function renormalization of the pion in the medium plays an essential role to induce the decrease of the pion decay constant as the order parameter of chiral transition. An emphasis is also put on the importance to examine the scalar and vector channels simultaneously for exploring the possible restoration of chiral symmetry.

preprint2002arXiv

Lattice study of "f$_{0}$(600) or $σ$"

We investigate the propagator of "f$_{0}$(600) or the $σ$" by the full-QCD simulation with Wilson fermions. We calculate the mesonic correlator in the I=0, $J^P=0^{+}$ channel on the $8^{3} \times 16$ lattice. Plaquet action and Wilson fermion action are adopted. A coupling constant $β$ is set to 4.8 and three kinds of hopping parameter, $κ$=0.1846, 0.1874 and 0.1891 are assayed. The disconnected diagram in the propagator is evaluated through taking average over 500 or 1000 Z2 noise. Simulations with the larger hopping parameter provide us with less noisy results. Though the statistics is not yet enough, our results indicate the existence of a pole with a mass in almost the same order as that of the $ρ$.

preprint1997arXiv

Renormalization-group Resummation of Divergent Series of the Perturbative Wave Function of Quantum Anharmonic Oscillator

The renormalization group method is applied for obtaining the asymptotic form of the wave function of the quantum anharmonic oscillator by resumming the perturbation series. It is shown that the resumed series is the cumulant of the naive perturbation series. Working out up to the sixth order and performing a further resummation proposed by Bender and Bettencourt, we find that the agreement with the WKB result becomes worse in the higher orders than the fourth at which the agreement is the best.

preprint1997arXiv

Some Phenomenological Properties of the Chiral Transition in QCD

On the basis of the NJL model as an effective theory of QCD and analogies with condensed matter physics, we extract simple physical pictures of the chiral phase transition at finite temperature $T$ and/or chemical potential $μ$ and predict interesting phenomena associated with the phase transition which may be seen in relativistic heavy-ion collisions to be performed in RHIC or LHC. We indicate that the macroscopic equation describing non-equilibrium phenomena of the chiral transition like DCC might be given by a diffusion-like equation rather than a wave equation especially when the baryon density is present.

preprint1995arXiv

CHIRAL SYMMETRY AND $U_A(1)$ ANOMALY IN AN EFFECTIVE THEORY OF QCD

We show on the basis of an effective theory of QCD that a wide variety of observables in the hadron world is governed by the chiral symmetry together with an interplay between the axial anomaly and the explicit symmetry breaking due to the current quark mass. We also discuss the nature of the chiral transition at finite temperature and related dynamical phenomena using the effective Lagrangian. Some phenomenological implications of the small vector coupling (``vector limit'') at high temperatures are suggested.

preprint1995arXiv

Chiral Symmetry and Scalar Meson in Hadron and Nuclear Physics

After giving a short introduction to the Nambu-Jona-Lasinio model with an anomaly term, we show the importance of the scalar-scalar correlation in the low-energy hadron dynamics, which correlation may be summarized by a scalar-isoscalar meson, the sigma meson. The discussion is based on the chiral quark model with the sigma-meson degrees of freedom. Possible experiments are proposed to produce the elusive meson in a nucleus and detect it. In relation to a precursory soft mode for the chiral transition, the reason is clarified why the dynamic properties of the superconductor may be described by the diffusive time-dependent Ginzburg-Landau (TDGL) equation. We indicate the chiral symmetry plays a significant role also in nuclei; one may say that the stability of nuclei is due to the chiral symmetry of QCD.