-
Experimental control of mode-competition dynamics in a chaotic multimode semiconductor laser for decision making
Authors:
Ryugo Iwami,
Takatomo Mihana,
Kazutaka Kanno,
Makoto Naruse,
Atsushi Uchida
Abstract:
Photonic computing has been widely used to accelerate the computational performance in machine learning. Photonic decision-making is a promising approach that uses photonic computing technologies to solve the multi-armed bandit problem based on reinforcement learning. Photonic decision making using chaotic mode competition dynamics has been proposed. However, the experimental conditions for achiev…
▽ More
Photonic computing has been widely used to accelerate the computational performance in machine learning. Photonic decision-making is a promising approach that uses photonic computing technologies to solve the multi-armed bandit problem based on reinforcement learning. Photonic decision making using chaotic mode competition dynamics has been proposed. However, the experimental conditions for achieving superior decision-making performance have not yet been established. In this study, we experimentally investigate the mode competition dynamics in a chaotic multimode semiconductor laser in the presence of optical feedback and injection. We control chaotic mode-competition dynamics via optical injection, and we found that positive wavelength detuning results in an efficient mode concentration to one of the longitudinal modes with a small optical injection power. We experimentally investigate two-dimensional bifurcation map of the laser dynamics of the total intensity. Complex mixed dynamics are observed in the presence of optical feedback and injection. We experimentally conduct decision making to solve the bandit problem using chaotic mode-competition dynamics. A fast mode concentration property is observed at positive wavelength detuning, which resulted in a fast convergence of the correct decision rate. Our findings could be useful in accelerating the decision-making performance in adaptive optical networks using reinforcement learning.
△ Less
Submitted 27 December, 2023;
originally announced December 2023.
-
Solving multi-armed bandit problems using a chaotic microresonator comb
Authors:
Jonathan Cuevas,
Ryugo Iwami,
Atsushi Uchida,
Kaoru Minoshima,
Naoya Kuse
Abstract:
The Multi-Armed Bandit (MAB) problem, foundational to reinforcement learning-based decision-making, addresses the challenge of maximizing rewards amidst multiple uncertain choices. While algorithmic solutions are effective, their computational efficiency diminishes with increasing problem complexity. Photonic accelerators, leveraging temporal and spatial-temporal chaos, have emerged as promising a…
▽ More
The Multi-Armed Bandit (MAB) problem, foundational to reinforcement learning-based decision-making, addresses the challenge of maximizing rewards amidst multiple uncertain choices. While algorithmic solutions are effective, their computational efficiency diminishes with increasing problem complexity. Photonic accelerators, leveraging temporal and spatial-temporal chaos, have emerged as promising alternatives. However, despite these advancements, current approaches either compromise computation speed or amplify system complexity. In this paper, we introduce a chaotic microresonator frequency comb (chaos comb) to tackle the MAB problem, where each comb mode is assigned to a slot machine. Through a proof-of-concept experiment, we employ 44 comb modes to address an MAB with 44 slot machines, demonstrating performance competitive with both conventional software algorithms and other photonic methods. Further, the scalability of decision making is explored with up to 512 slot machines using experimentally obtained temporal chaos in different time slots. Power-law scalability is achieved with an exponent of 0.96, outperforming conventional software-based algorithms. Moreover, we find that a numerically calculated chaos comb accurately reproduces experimental results, paving the way for discussions on strategies to increase the number of slot machines.
△ Less
Submitted 22 August, 2023; v1 submitted 21 August, 2023;
originally announced August 2023.
-
Chaotic mode-competition dynamics in a multimode semiconductor laser with optical feedback and injection
Authors:
Ryugo Iwami,
Kazutaka Kanno,
Atsushi Uchida
Abstract:
Photonic computing is attracting increasing interest to accelerate information processing in machine learning applications. The mode-competition dynamics of multimode semiconductor lasers is useful for solving the multi-armed bandit problem in reinforcement learning for computing applications. In this study, we numerically evaluate the chaotic mode-competition dynamics in a multimode semiconductor…
▽ More
Photonic computing is attracting increasing interest to accelerate information processing in machine learning applications. The mode-competition dynamics of multimode semiconductor lasers is useful for solving the multi-armed bandit problem in reinforcement learning for computing applications. In this study, we numerically evaluate the chaotic mode-competition dynamics in a multimode semiconductor laser with optical feedback and injection. We observe the chaotic mode-competition dynamics among the longitudinal modes and control them by injecting an external optical signal into one of the longitudinal modes. We define the dominant mode as the mode with the maximum intensity; the dominant-mode ratio for the injected mode increases as the optical injection strength increases. We find that the characteristics of the dominant mode ratio in terms of the optical injection strength are different among the modes owing to the different optical feedback phases. We propose a control technique for the characteristics of the dominant mode ratio by precisely tuning the initial optical frequency detuning between the optical injection signal and injected mode. We also evaluate the relationship between the region for the large dominant mode ratio and injection locking range. The region for the large dominant mode ratio does not correspond to the injection-locking range. This discrepancy results from the complex mode-competition dynamics in multimode semiconductor lasers with both optical feedback and injection. This control technique of chaotic mode-competition dynamics in multimode lasers is promising for applications in reinforcement learning and reservoir computing as photonic artificial intelligence.
△ Less
Submitted 15 November, 2022;
originally announced November 2022.
-
Controlling chaotic itinerancy in laser dynamics for reinforcement learning
Authors:
Ryugo Iwami,
Takatomo Mihana,
Kazutaka Kanno,
Satoshi Sunada,
Makoto Naruse,
Atsushi Uchida
Abstract:
Photonic artificial intelligence has attracted considerable interest in accelerating machine learning; however, the unique optical properties have not been fully utilized for achieving higher-order functionalities. Chaotic itinerancy, with its spontaneous transient dynamics among multiple quasi-attractors, can be employed to realize brain-like functionalities. In this paper, we propose a method fo…
▽ More
Photonic artificial intelligence has attracted considerable interest in accelerating machine learning; however, the unique optical properties have not been fully utilized for achieving higher-order functionalities. Chaotic itinerancy, with its spontaneous transient dynamics among multiple quasi-attractors, can be employed to realize brain-like functionalities. In this paper, we propose a method for controlling the chaotic itinerancy in a multi-mode semiconductor laser to solve a machine learning task, known as the multi-armed bandit problem, which is fundamental to reinforcement learning. The proposed method utilizes ultrafast chaotic itinerant motion in mode competition dynamics controlled via optical injection. We found that the exploration mechanism is completely different from a conventional searching algorithm and is highly scalable, outperforming the conventional approaches for large-scale bandit problems. This study paves the way to utilize chaotic itinerancy for effectively solving complex machine learning tasks as photonic hardware accelerators.
△ Less
Submitted 12 May, 2022;
originally announced May 2022.
-
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
Authors:
Kazuyuki Kanaya,
Mizuki Shirogane,
Shinji Ejiri,
Ryo Iwami,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takashi Umeda
Abstract:
We study the latent heat and the pressure gap between the hot and cold phases at the first-order transition temperature $T=T_c$ of SU(3) Yang-Mills theory, using the small flow-time expansion (SF$t$X) method based on the gradient flow. We first examine alternative procedures in the SFtX method -- the order of the continuum and vanishing flow-time extrapolations. We confirm that the final results a…
▽ More
We study the latent heat and the pressure gap between the hot and cold phases at the first-order transition temperature $T=T_c$ of SU(3) Yang-Mills theory, using the small flow-time expansion (SF$t$X) method based on the gradient flow. We first examine alternative procedures in the SFtX method -- the order of the continuum and vanishing flow-time extrapolations. We confirm that the final results adopting the two orders, as well as other alternatives in which the perturbative order of the matching coefficients and the renormalization scale of the flow scheme are varied, are all consistent with each other. We also confirm $Δp$ is consistent with zero, as expected from the dynamical balance of two phases at $T_c$. For the latent heat in the continuum limit, we find $Δε/T^4 = 1.117(40)$ for the spatial volume $L^3$ corresponding to the aspect ratio $N_s/N_t=T_cL=8$ and $1.349(38)$ for $N_s/N_t=6$. From hysteresis curves, we show that the entropy density in the hot phase is sensitive to the spatial volume, while that in the confined phase is insensitive.
△ Less
Submitted 21 October, 2021;
originally announced October 2021.
-
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
Authors:
Mizuki Shirogane,
Shinji Ejiri,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takashi Umeda
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 January, 2021; v1 submitted 20 November, 2020;
originally announced November 2020.
-
Determination of the endpoint of the first order deconfiniement phase transition in the heavy quark region of QCD
Authors:
Shinji Ejiri,
Shota Itagaki,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Atsushi Kiyohara,
Mizuki Shirogane,
Yusuke Taniguchi,
Takashi Umeda
Abstract:
We study the endpoint of the first order deconfinement phase transition of 2 and 2+1 flavor QCD in the heavy quark region. We perform simulations of quenched QCD and apply the reweighting method to study the heavy quark region. The quark determinant for the reweighting is evaluated by a hopping parameter expansion. To reduce the overlap problem, we introduce an external source term of the Polyakov…
▽ More
We study the endpoint of the first order deconfinement phase transition of 2 and 2+1 flavor QCD in the heavy quark region. We perform simulations of quenched QCD and apply the reweighting method to study the heavy quark region. The quark determinant for the reweighting is evaluated by a hopping parameter expansion. To reduce the overlap problem, we introduce an external source term of the Polyakov loop in the simulation. We study the location of critical point at which the first order phase transition changes to crossover by investigating the histogram of the Polyakov loop and applying the finite-size scaling analysis. We estimate the truncation error of the hopping parameter expansion, and discuss the lattice spacing dependence and the spatial volume dependence in the result of the critical point.
△ Less
Submitted 24 December, 2019;
originally announced December 2019.
-
End point of the first-order phase transition of QCD in the heavy quark region by reweighting from quenched QCD
Authors:
Shinji Ejiri,
Shota Itagaki,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Atsushi Kiyohara,
Mizuki Shirogane,
Takashi Umeda
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 March, 2020; v1 submitted 22 December, 2019;
originally announced December 2019.
-
Equation of state near the first order phase transition point of SU(3) gauge theory using gradient flow
Authors:
Mizuki Shirogane,
Shinji Ejiri,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takashi Umeda
Abstract:
We study energy gap (latent heat) between the hot and cold phases at the first order phase transition point of the SU(3) gauge theory. Performing simulations on lattices with various spatial volumes and lattice spacings, we calculate the energy gap by a method using the Yang-Mills gradient flow and compare it with that by the conventional derivative method.
We study energy gap (latent heat) between the hot and cold phases at the first order phase transition point of the SU(3) gauge theory. Performing simulations on lattices with various spatial volumes and lattice spacings, we calculate the energy gap by a method using the Yang-Mills gradient flow and compare it with that by the conventional derivative method.
△ Less
Submitted 18 December, 2018; v1 submitted 10 November, 2018;
originally announced November 2018.
-
Equation of state in (2+1)-flavor QCD at physical point with improved Wilson fermion action using gradient flow
Authors:
WHOT-QCD Collaboration,
:,
Kazuyuki Kanaya,
Shinji Ejiri,
Ryo Iwami,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takash Umeda
Abstract:
We study the energy-momentum tensor and the equation of state as well as the chiral condensate in (2+1)-flavor QCD at the physical point applying the method of Makino and Suzuki based on the gradient flow. We adopt a nonperturbatively O(a)-improved Wilson quark action and the renormalization group-improved Iwasaki gauge action. At Lattice 2016, we have presented our preliminary results of our stud…
▽ More
We study the energy-momentum tensor and the equation of state as well as the chiral condensate in (2+1)-flavor QCD at the physical point applying the method of Makino and Suzuki based on the gradient flow. We adopt a nonperturbatively O(a)-improved Wilson quark action and the renormalization group-improved Iwasaki gauge action. At Lattice 2016, we have presented our preliminary results of our study in (2+1)-flavor QCD at a heavy u, d quark mass point. We now extend the study to the physical point and perform finite-temperature simulations in the range T \simeq 155--544 MeV (Nt = 4--14 including odd Nt's) at a \simeq 0.09 fm. We show our final results of the heavy QCD study and present some preliminary results obtained at the physical point so far.
△ Less
Submitted 27 October, 2017;
originally announced October 2017.
-
Determination of latent heat at the finite temperature phase transition of SU(3) gauge theory
Authors:
Shinji Ejiri,
Ryo Iwami,
Mizuki Shirogane,
Naoki Wakabayashi,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takashi Umeda
Abstract:
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 $N_t=6,$ 8 and 12 and extrapolate the results to the continuum limit. The energy density and pressure are evaluated by the…
▽ More
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 $N_t=6,$ 8 and 12 and extrapolate the results to the continuum limit. The energy density and pressure are evaluated by the derivative method with nonperturabative anisotropy coefficients. We find that the pressure gap vanishes at all values of $N_t$. 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.$ We also tested a method using the Yang-Mills gradient flow. The preliminary results are consistent with those by the derivative method within the error.
△ Less
Submitted 30 January, 2017;
originally announced January 2017.
-
O(4) scaling analysis in two-flavor QCD at finite temperature and density with improved Wilson quarks
Authors:
T. Umeda,
S. Ejiri,
R. Iwami,
K. Kanaya,
H. Ohno,
A. Uji,
N. Wakabayashi,
S. Yoshida
Abstract:
We study the curvature of the chiral transition/crossover line between the low-temperature hadronic phase and the high-temperature quark-gluon-plasma phase at low densities, performing simulations of two-flavor QCD with improved Wilson quarks. After confirming that the chiral order parameter defined by a Ward-Takahashi identity is consistent with the scaling of the O(4) universality class at zero…
▽ More
We study the curvature of the chiral transition/crossover line between the low-temperature hadronic phase and the high-temperature quark-gluon-plasma phase at low densities, performing simulations of two-flavor QCD with improved Wilson quarks. After confirming that the chiral order parameter defined by a Ward-Takahashi identity is consistent with the scaling of the O(4) universality class at zero chemical potential, we extend the scaling analysis to finite chemical potential to determine the curvature of the chiral transition/crossover line at low densities assuming the O(4) universality. To convert the curvature in lattice units to that of the $T_c(μ_B)$ in physical units, we evaluate the lattice scale applying a gradient flow method. We find $κ=0.0006(1)$ in the chiral limit, which is much smaller than that obtained in (2+1)-flavor QCD with improved staggered quarks.
△ Less
Submitted 8 January, 2017; v1 submitted 30 December, 2016;
originally announced December 2016.
-
Temperature dependence of topological susceptibility using gradient flow
Authors:
Yusuke Taniguchi,
Shinji Ejiri,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Takashi Umeda,
Ryo Iwami,
Naoki Wakabayashi
Abstract:
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 grad…
▽ More
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.
△ Less
Submitted 8 November, 2016;
originally announced November 2016.
-
Equation of state in (2+1)-flavor QCD with gradient flow
Authors:
Kazuyuki Kanaya,
Shinji Ejiri,
Ryo Iwami,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Yusuke Taniguchi,
Takashi Umeda,
Naoki Wakabayashi
Abstract:
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…
▽ More
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.
△ Less
Submitted 29 October, 2016;
originally announced October 2016.
-
Exploring Nf=2+1 QCD thermodynamics from the gradient flow
Authors:
Yusuke Taniguchi,
Shinji Ejiri,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa,
Hiroshi Suzuki,
Takashi Umeda,
Naoki Wakabayashi
Abstract:
The energy-momentum tensor plays an important role in QCD thermodynamics. Its expectation value contains information of the pressure and the energy density as its diagonal part. Further properties like viscosity and specific heat can be extracted from its correlation function. Recently a new method based on the gradient flow was introduced to calculate the energy-momentum tensor on the lattice, an…
▽ More
The energy-momentum tensor plays an important role in QCD thermodynamics. Its expectation value contains information of the pressure and the energy density as its diagonal part. Further properties like viscosity and specific heat can be extracted from its correlation function. Recently a new method based on the gradient flow was introduced to calculate the energy-momentum tensor on the lattice, and has been successfully applied to quenched QCD.
In this paper, we apply the gradient flow method to calculate the energy-momentum tensor in (2+1)-flavor QCD. As the first application of the method with dynamical quarks, we study at a single but fine lattice spacing a=0.07 fm with heavy u and d quarks ($m_π/m_ρ=0.63$) and approximately physical s quark. Performing simulations on lattices with Nt=16 to 4, the temperature range of T=174-697 MeV is covered. We find that the results of the pressure and the energy density by the gradient flow method are consistent with the previous results using the T-integration method at T<280 MeV, while the results show disagreement at T>350 MeV (Nt<8), presumably due to the small-Nt lattice artifact of $O((aT)^2)=O(1/N_t^2)$.
We also apply the gradient flow method to evaluate the chiral condensate taking advantage of the gradient flow method that renormalized quantities can be directly computed avoiding the difficulty of explicit chiral violation with lattice quarks. We compute the renormalized chiral condensate in the MS-bar scheme at renormalization scale $μ=2$ GeV with a high precision to study the temperature dependence of the chiral condensate and its disconnected susceptibility. Even with the Wilson-type quark action, we obtain the chiral condensate and its disconnected susceptibility showing a clear signal of pseudocritical temperature at T~190 MeV related to the chiral restoration crossover.
△ Less
Submitted 20 February, 2019; v1 submitted 6 September, 2016;
originally announced September 2016.
-
Latent heat at the first order phase transition point of SU(3) gauge theory
Authors:
Mizuki Shirogane,
Shinji Ejiri,
Ryo Iwami,
Kazuyuki Kanaya,
Masakiyo Kitazawa
Abstract:
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 e…
▽ More
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.$
△ Less
Submitted 10 May, 2016;
originally announced May 2016.
-
Many flavor approach to study the nature of chiral phase transition of two-flavor QCD
Authors:
Norikazu Yamada,
Shinji Ejiri,
Ryo Iwami
Abstract:
We perform lattice numerical simulations to study the phase transition of QCD at finite temperature to clarify the nature of the transition of massless two flavor QCD. We investigate QCD with two light and Nf heavy quarks instead of two-flavor QCD, and focus on the light quark mass dependence of the critical heavy mass, below which the transition is of first order. The heavy quarks are incorporate…
▽ More
We perform lattice numerical simulations to study the phase transition of QCD at finite temperature to clarify the nature of the transition of massless two flavor QCD. We investigate QCD with two light and Nf heavy quarks instead of two-flavor QCD, and focus on the light quark mass dependence of the critical heavy mass, below which the transition is of first order. The heavy quarks are incorporated into two flavor configurations in the form of the hopping parameter expansion through the reweighting technique. The nature of the transition is identified by the shape of the constraint effective potential at the critical temperature. Our result indicates that the critical heavy mass remains finite in the chiral limit of the two flavors, suggesting the phase transition of massless two-flavor QCD is of second order.
△ Less
Submitted 15 February, 2016;
originally announced February 2016.
-
Critical point search from an extended parameter space of lattice QCD at finite temperature and density
Authors:
Shinji Ejiri,
Ryo Iwami,
Norikazu Yamada
Abstract:
Aiming to understand the phase structure of lattice QCD at nonzero temperature and density, we study the phase transitions of QCD in an extended parameter space, where the number of flavor and quark masses are considered as parameters. Performing simulations of 2 flavor QCD and using the reweighting method, we investigate (2+Nf) flavor QCD at finite density, where two light flavors and Nf massive…
▽ More
Aiming to understand the phase structure of lattice QCD at nonzero temperature and density, we study the phase transitions of QCD in an extended parameter space, where the number of flavor and quark masses are considered as parameters. Performing simulations of 2 flavor QCD and using the reweighting method, we investigate (2+Nf) flavor QCD at finite density, where two light flavors and Nf massive flavors exist. Calculating probability distribution functions, we determine the critical surface terminating first order phase transitions in the parameter space of the light quark mass, the heavy quark mass and the chemical potential. Through the study of the many flavor system, we discuss the phase structure of QCD at finite density.
△ Less
Submitted 18 January, 2016;
originally announced January 2016.
-
Many flavor approach to study the critical point in finite density QCD
Authors:
Ryo Iwami,
Shinji Ejiri,
Norikazu Yamada
Abstract:
We discuss the QCD critical point at finite density through the study of many flavor QCD, in which two light flavors and Nf massive flavors exist. Performing simulations of QCD with two flavors of improved Wilson fermions, we calculate probability distribution functions of many flavor QCD at finite temperature and density. The dynamical effects of massive flavors and the chemical potential are add…
▽ More
We discuss the QCD critical point at finite density through the study of many flavor QCD, in which two light flavors and Nf massive flavors exist. Performing simulations of QCD with two flavors of improved Wilson fermions, we calculate probability distribution functions of many flavor QCD at finite temperature and density. The dynamical effects of massive flavors and the chemical potential are added using the reweighting technique. From the shape of the distribution functions, we determine the critical surface separating the first order transition and crossover regions in the parameter space of the light and massive quark masses and the chemical potentials. It is found that the critical massive quark mass becomes larger as the chemical potential increases in (2+Nf) flavor QCD. The indication to the (2+1) flavor QCD is then discussed.
△ Less
Submitted 28 December, 2015;
originally announced December 2015.
-
Exploring the nature of chiral phase transition in two-flavor QCD using extra heavy quarks
Authors:
Shinji Ejiri,
Ryo Iwami,
Norikazu Yamada
Abstract:
Chiral phase transition of two flavor QCD at finite quark masses is known to be crossover except near the chiral limit, but it can turn to a first order transition when adding many extra flavors. This property is used to explore the nature of the phase transition of massless two flavor QCD using lattice numerical simulations. The extra heavy flavors being incorporated in the form of the hopping pa…
▽ More
Chiral phase transition of two flavor QCD at finite quark masses is known to be crossover except near the chiral limit, but it can turn to a first order transition when adding many extra flavors. This property is used to explore the nature of the phase transition of massless two flavor QCD using lattice numerical simulations. The extra heavy flavors being incorporated in the form of the hopping parameter expansion through the reweighting, the number of the extra flavors and their masses appear only in a single parameter, defined by $h$. We determine the critical value of the parameter, at which the first order and the crossover regions are separated, and examine its dependence on the two flavor mass. The lattice calculations are carried out at $N_t$=4, and show that the critical value does not depend on the two flavor mass in the range we have studied ($0.46 \le m_π/m_ρ\le 0.66$) and appears to remain finite and positive in the chiral limit, suggesting that the phase transition of massless two flavor QCD is of second order.
△ Less
Submitted 22 March, 2016; v1 submitted 19 November, 2015;
originally announced November 2015.
-
Towards the QCD equation of state at the physical point using Wilson fermion
Authors:
T. Umeda,
S. Ejiri,
R. Iwami,
K. Kanaya
Abstract:
We study the (2+1)-flavor QCD at nonzero temperatures using nonperturbatively improved Wilson quarks of the physical masses by the fixed scale approach. We perform physical point simulations at finite temperatures with the coupling parameters which were adopted by the PACS-CS Collaboration in their studies using the reweighting technique. Zero temperature values are obtained on the PACS-CS configu…
▽ More
We study the (2+1)-flavor QCD at nonzero temperatures using nonperturbatively improved Wilson quarks of the physical masses by the fixed scale approach. We perform physical point simulations at finite temperatures with the coupling parameters which were adopted by the PACS-CS Collaboration in their studies using the reweighting technique. Zero temperature values are obtained on the PACS-CS configurations which are open to the public on the ILDG/JLDG. Finite temperature configurations are generated with the RHMC algorithm. The lattice sizes are $32^3 \times N_t$ with $N_t=14$, 13, $\cdots$, 4 which correspond to $T \approx 160$--550 MeV. We present results of some basic observables at these temperatures and the status of our calculation of the equation of state.
△ Less
Submitted 14 November, 2015;
originally announced November 2015.
-
Multipoint reweighting method and its applications to lattice QCD
Authors:
R. Iwami,
S. Ejiri,
K. Kanaya,
Y. Nakagawa,
D. Yamamoto,
T. Umeda
Abstract:
The reweighting method is widely used in numerical studies of QCD, in particular, for the cases in which the conventional Monte-Carlo method cannot be applied directly, e.g., finite density QCD. However, the application range of the reweighing method is restricted due to several problems. One of the most severe problems here is the overlap problem. To solve it, we examine a multipoint reweighting…
▽ More
The reweighting method is widely used in numerical studies of QCD, in particular, for the cases in which the conventional Monte-Carlo method cannot be applied directly, e.g., finite density QCD. However, the application range of the reweighing method is restricted due to several problems. One of the most severe problems here is the overlap problem. To solve it, we examine a multipoint reweighting method in which simulations at several simulation points are combined in the data analyses. We systematically study the applicability and limitation of the multipoint reweighting method in two-flavor QCD at zero density. Measuring histograms of physical quantities at a series of simulation points, we apply the multipoint reweighting method to calculate the meson masses as continuous functions of the gauge coupling $β$ and the hopping parameters $κ$. We then determine lines of constant physics and beta functions, which are needed in a calculation of the equation of state at finite temperature.
△ Less
Submitted 20 November, 2015; v1 submitted 7 August, 2015;
originally announced August 2015.
-
Multipoint reweighting method and beta functions for the calculation of QCD equation of state
Authors:
Ryo Iwami,
S. Ejiri,
K. Kanaya,
Y. Nakagawa,
T. Umeda,
D. Yamamoto
Abstract:
We study a reweighting method aiming at numerical studies of QCD at finite density, in which the conventional Monte-Carlo method cannot be applied directly. One of the most important problems in the reweighting method is the overlap problem. To solve it, we propose to perform simulations at several simulation points and combine their results in the data analyses. In this report, we introduce this…
▽ More
We study a reweighting method aiming at numerical studies of QCD at finite density, in which the conventional Monte-Carlo method cannot be applied directly. One of the most important problems in the reweighting method is the overlap problem. To solve it, we propose to perform simulations at several simulation points and combine their results in the data analyses. In this report, we introduce this multipoint reweighting method and test if the method works well by measuring histograms of physical quantities. Using this method, we calculate the meson masses as continuous functions of the gauge coupling beta and the hopping parameters kappa in QCD at zero density. We then determine lines of constant physics in the (beta, kappa) space and evaluate the derivatives of the lattice spacing with respect to beta and kappa along the lines of constant physics (inverse of the beta functions), which are needed in a calculation of the equation of state.
△ Less
Submitted 26 January, 2015;
originally announced January 2015.
