-
Impact of local CP-odd domain in hot QCD on axionic domain-wall interpretation for NANOGrav 15-year Data
Authors:
Linlin Huang,
Yuanyuan Wang,
He-Xu Zhang,
Shinya Matsuzaki,
Hiroyuki Ishida,
Mamiya Kawaguchi,
Akio Tomiya
Abstract:
We argue that the axionic domain-wall with a QCD bias may be incompatible with the NANOGrav 15-year data on a stochastic gravitational wave (GW) background, when the domain wall network collapses in the hot-QCD induced local CP-odd domain. This is due to the drastic suppression of the QCD bias set by the QCD topological susceptibility in the presence of the CP-odd domain with nonzero $θ$ parameter…
▽ More
We argue that the axionic domain-wall with a QCD bias may be incompatible with the NANOGrav 15-year data on a stochastic gravitational wave (GW) background, when the domain wall network collapses in the hot-QCD induced local CP-odd domain. This is due to the drastic suppression of the QCD bias set by the QCD topological susceptibility in the presence of the CP-odd domain with nonzero $θ$ parameter of order one which the QCD sphaleron could generate. We quantify the effect on the GW signals by working on a low-energy effective model of Nambu-Jona-Lasinio type in the mean field approximation. We find that only at $θ=π$, the QCD bias tends to get significantly large enough due to the criticality of the thermal CP restoration, which would, however, give too big signal strengths to be consistent with the NANOGrav 15-year data and would also be subject to the strength of the phase transition at the criticality.
△ Less
Submitted 23 May, 2024; v1 submitted 17 March, 2024;
originally announced March 2024.
-
Sparse modeling approach to extract spectral functions with covariance of Euclidean-time correlators of lattice QCD
Authors:
Junichi Takahashi,
Hiroshi Ohno,
Akio Tomiya
Abstract:
We present our sparse modeling study to extract spectral functions from Euclidean-time correlation functions. In this study covariance between different Euclidean times of the correlation function is taken into account, which was not done in previous studies. In order to check applicability of the method, we firstly test it with mock data which imitate possible charmonium spectral functions. Then,…
▽ More
We present our sparse modeling study to extract spectral functions from Euclidean-time correlation functions. In this study covariance between different Euclidean times of the correlation function is taken into account, which was not done in previous studies. In order to check applicability of the method, we firstly test it with mock data which imitate possible charmonium spectral functions. Then, we extract spectral functions from correlation functions obtained from lattice QCD at finite temperature.
△ Less
Submitted 26 November, 2023;
originally announced November 2023.
-
Equivariant Transformer is all you need
Authors:
Akio Tomiya,
Yuki Nagai
Abstract:
Machine learning, deep learning, has been accelerating computational physics, which has been used to simulate systems on a lattice. Equivariance is essential to simulate a physical system because it imposes a strong induction bias for the probability distribution described by a machine learning model. This reduces the risk of erroneous extrapolation that deviates from data symmetries and physical…
▽ More
Machine learning, deep learning, has been accelerating computational physics, which has been used to simulate systems on a lattice. Equivariance is essential to simulate a physical system because it imposes a strong induction bias for the probability distribution described by a machine learning model. This reduces the risk of erroneous extrapolation that deviates from data symmetries and physical laws. However, imposing symmetry on the model sometimes occur a poor acceptance rate in self-learning Monte-Carlo (SLMC). On the other hand, Attention used in Transformers like GPT realizes a large model capacity. We introduce symmetry equivariant attention to SLMC. To evaluate our architecture, we apply it to our proposed new architecture on a spin-fermion model on a two-dimensional lattice. We find that it overcomes poor acceptance rates for linear models and observe the scaling law of the acceptance rate as in the large language models with Transformers.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
Self-learning Monte Carlo with equivariant Transformer
Authors:
Yuki Nagai,
Akio Tomiya
Abstract:
Machine learning and deep learning have revolutionized computational physics, particularly the simulation of complex systems. Equivariance is essential for simulating physical systems because it imposes a strong inductive bias on the probability distribution described by a machine learning model. However, imposing symmetry on the model can sometimes lead to poor acceptance rates in self-learning M…
▽ More
Machine learning and deep learning have revolutionized computational physics, particularly the simulation of complex systems. Equivariance is essential for simulating physical systems because it imposes a strong inductive bias on the probability distribution described by a machine learning model. However, imposing symmetry on the model can sometimes lead to poor acceptance rates in self-learning Monte Carlo (SLMC). Here, we introduce a symmetry equivariant attention mechanism for SLMC, which can be systematically improved. We evaluate our architecture on a spin-fermion model (\textit{i.e.}, double exchange model) on a two-dimensional lattice. Our results show that the proposed method overcomes the poor acceptance rates of linear models and exhibits a similar scaling law to large language models, with model quality monotonically increasing with the number of layers. Our work paves the way for the development of more accurate and efficient Monte Carlo algorithms with machine learning for simulating complex physical systems.
△ Less
Submitted 2 July, 2024; v1 submitted 20 June, 2023;
originally announced June 2023.
-
Use of Schwinger-Dyson equation in constructing an approximate trivializing map
Authors:
Peter Boyle,
Taku Izubuchi,
Luchang Jin,
Chulwoo Jung,
Christoph Lehner,
Nobuyuki Matsumoto,
Akio Tomiya
Abstract:
We construct an approximate trivializing map by using a Schwinger-Dyson equation. The advantage of this method is that: (1) The basis for the flow kernel can be chosen arbitrarily by hand. (2) It can be applied to the general action of interest. (3) The coefficients in the kernel are determined by lattice estimates of the observables, which does not require analytic calculations beforehand. We per…
▽ More
We construct an approximate trivializing map by using a Schwinger-Dyson equation. The advantage of this method is that: (1) The basis for the flow kernel can be chosen arbitrarily by hand. (2) It can be applied to the general action of interest. (3) The coefficients in the kernel are determined by lattice estimates of the observables, which does not require analytic calculations beforehand. We perform the HMC with the effective action obtained by the Schwinger-Dyson method, and show that we can have better control of the effective action than the known $t$-expansion construction. However, the algorithmic overhead is still large and overwhelming the gain though faster decorrelation is observed for long-range observables in some cases. This contribution reports the preliminary results of this attempt.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
GomalizingFlow.jl: A Julia package for Flow-based sampling algorithm for lattice field theory
Authors:
Akio Tomiya,
Satoshi Terasaki
Abstract:
GomalizingFlow.jl: is a package to generate configurations for quantum field theory on the lattice using the flow based sampling algorithm in Julia programming language. This software serves two main purposes: to accelerate research of lattice QCD with machine learning with easy prototyping, and to provide an independent implementation to an existing public Jupyter notebook in Python/PyTorch. Goma…
▽ More
GomalizingFlow.jl: is a package to generate configurations for quantum field theory on the lattice using the flow based sampling algorithm in Julia programming language. This software serves two main purposes: to accelerate research of lattice QCD with machine learning with easy prototyping, and to provide an independent implementation to an existing public Jupyter notebook in Python/PyTorch. GomalizingFlow.jl implements, the flow based sampling algorithm, namely, RealNVP and Metropolis-Hastings test for two dimension and three dimensional scalar field, which can be switched by a parameter file. HMC for that theory also implemented for comparison. This package has Docker image, which reduces effort for environment construction. This code works both on CPU and NVIDIA GPU.
△ Less
Submitted 4 October, 2023; v1 submitted 18 August, 2022;
originally announced August 2022.
-
Fate of the first-order chiral phase transition in QCD: Implications for dark QCD studied via a Nambu-Jona-Lasinio model
Authors:
Yuanyuan Wang,
Mamiya Kawaguchi,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
The first-order nature of the chiral phase transition in QCD-like theories can play crucial roles to address a dark side of the Universe, where the created out-of equilibrium is essential to serve as cosmological and astrophysical probes such as gravitational wave productions, which have extensively been explored. This interdisciplinary physics is built based on a widely-accepted conjecture that t…
▽ More
The first-order nature of the chiral phase transition in QCD-like theories can play crucial roles to address a dark side of the Universe, where the created out-of equilibrium is essential to serve as cosmological and astrophysical probes such as gravitational wave productions, which have extensively been explored. This interdisciplinary physics is built based on a widely-accepted conjecture that the thermal chiral phase transition in QCD-like theories with massless (light) three flavors is of first order. We find that such a first order feature may not hold, when ordinary or dark quarks are externally coupled to a weak enough background field of photon or dark photon (which we collectively call a ``magnetic" field). We assume that a weak ``magnetic" background field could be originated from some ``magnetogenesis" in the early Universe. We work on a Nambu-Jona-Lasinio model which can describe the chiral phase transition in a wide class of QCD-like theories. We show that in the case with massless (light) three flavors, the first-order feature goes away when $2 f_π^2 \lesssim eB ( \ll (4 πf_π)^2)$, where $eB$ is the ``magnetic" field strength and $f_π$ the pion decay constant at the vacuum. This disappearance is the generic consequence of the presence of the ``magnetically" induced scale anomaly and the ``magnetic" catalysis for the chiral symmetry breaking, and would impact or constrain modeling dark QCD coupled to an external ``magnetic" field.
△ Less
Submitted 9 November, 2022; v1 submitted 8 August, 2022;
originally announced August 2022.
-
New aspect of chiral and axial breaking in QCD
Authors:
Chuan-Xin Cui,
Mamiya Kawaguchi,
Jin-Yang Li,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
Violation of the $U(1)$ axial symmetry in QCD is stricter than the chiral $SU(2)$ breaking, simply because of the presence of the quantum axial anomaly. If the QCD gauge coupling is sent to zero, the strength of the $U(1)$ axial breaking coincides with that of the chiral $SU(2)$ breaking, which we shall in short call an axial-chiral coincidence. This coincidence is trivial since QCD then becomes a…
▽ More
Violation of the $U(1)$ axial symmetry in QCD is stricter than the chiral $SU(2)$ breaking, simply because of the presence of the quantum axial anomaly. If the QCD gauge coupling is sent to zero, the strength of the $U(1)$ axial breaking coincides with that of the chiral $SU(2)$ breaking, which we shall in short call an axial-chiral coincidence. This coincidence is trivial since QCD then becomes a non-interacting theory. Actually, there exists another limit in the QCD parameter space, where an axial-chiral coincidence occurs even with nonzero QCD gauge coupling, that can be dubbed a nontrivial coincidence: it is the case with the massive light quarks $(m_l\neq 0)$ and the massless strange quark ($m_s=0$), due to the flavor-singlet nature of the topological susceptibility. This coincidence is robust and tied to the anomalous chiral Ward-Takahashi identity, which is operative even at hot QCD. This implies that the chiral $SU(2)$ symmetry is restored simultaneously with the $U(1)$ axial symmetry at high temperatures. This simultaneous restoration is independent of $m_l (\neq 0)$, hence is irrespective to the order of the chiral phase transition. In this paper, we discuss how the real-life QCD can be evolved from the nontrivial chiral-axial coincidence limit, by working on a Nambu-Jona-Lasinio model with the $U(1)$ axial anomaly contribution properly incorporated. It is shown that at high temperatures the large differences between the restorations of the chiral $SU(2)$ symmetry and the $U(1)$ axial symmetry for two light quarks and a sufficiently large current mass for the strange quark is induced by a significant interference of the topological susceptibility. Thus the deviation from the nontrivial coincidence, which is monitored by the strange quark mass controlling the topological susceptibility, provides a new way of understanding the chiral $SU(2)$ and $U(1)$ axial breaking in QCD.
△ Less
Submitted 10 March, 2024; v1 submitted 25 May, 2022;
originally announced May 2022.
-
Schwinger model at finite temperature and density with beta VQE
Authors:
Akio Tomiya
Abstract:
We investigate a quantum gauge theory at finite temperature and density using a variational algorithm for near-term quantum devices. We adapt $β$-VQE to evaluate thermal and quantum expectation values and study the phase diagram for massless Schwinger model along with the temperature and density. By compering the exact variational free energy, we find the variational algorithm work for $T>0$ and…
▽ More
We investigate a quantum gauge theory at finite temperature and density using a variational algorithm for near-term quantum devices. We adapt $β$-VQE to evaluate thermal and quantum expectation values and study the phase diagram for massless Schwinger model along with the temperature and density. By compering the exact variational free energy, we find the variational algorithm work for $T>0$ and $μ>0$ for the Schwinger model. No significant volume dependence of the variational free energy is observed in $μ/g \in[0, 1.4]$. We calculate the chiral condensate and take the continuum extrapolation. As a result, we obtain qualitative picture of the phase diagram for massless Schwinger model.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
Extensively parallelizable chiral fermion
Authors:
Yuki Nagai,
Akio Tomiya
Abstract:
Chiral symmetry is a key to investigating quantum physics, from condensed matter to particle physics. We propose a novel way of realizing a chiral fermion, known as the overlap-Dirac operator, without explicitly calculating the low modes of the Wilson-Dirac operator. We introduce a projection operator inspired by the Sakurai-Sugiura method and formulate the exact sign function and overlap-Dirac op…
▽ More
Chiral symmetry is a key to investigating quantum physics, from condensed matter to particle physics. We propose a novel way of realizing a chiral fermion, known as the overlap-Dirac operator, without explicitly calculating the low modes of the Wilson-Dirac operator. We introduce a projection operator inspired by the Sakurai-Sugiura method and formulate the exact sign function and overlap-Dirac operator with a contour-integral form. Like the Sakurai-Sugiura method, the proposing method is multi-scale parallelizable, which fits the multi-core/multi-GPGPU paradigm. We confirm that the quality of chiral symmetry realized with the proposed method is sufficient for double precision. We evaluate the strong scaling of the proposing method.
△ Less
Submitted 4 April, 2022;
originally announced April 2022.
-
QCD knows new quarks
Authors:
Chuan-Xin Cui,
Hiroyuki Ishida,
Mamiya Kawaguchi,
Jin-Yang Li,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
We find that a big gap between indicators for the breaking strengths of the global chiral SU(2) and U(1) axial symmetries in QCD of the Standard Model (SM) can be interpreted as a new fine-tuning problem. This may thus imply calling for a class of Beyond the SM, which turns out to favor having a new chiral symmetry and the associated massless new quark insensitive to the chiral SU(2) symmetry for…
▽ More
We find that a big gap between indicators for the breaking strengths of the global chiral SU(2) and U(1) axial symmetries in QCD of the Standard Model (SM) can be interpreted as a new fine-tuning problem. This may thus imply calling for a class of Beyond the SM, which turns out to favor having a new chiral symmetry and the associated massless new quark insensitive to the chiral SU(2) symmetry for the lightest up and down quarks, so that the fine-tuning is relaxed. Our statistical estimate shows that QCD of the SM is by more than 300 standard deviations off the desired parameter space, which is free from the fine-tuning, and the significance will be greater as the lattice measurements on the QCD hadron observables get more accurate. As one viable candidate, we introduce a dark QCD model with massless new quarks, which can survive current experimental, cosmological, and astrophysical limits, and also leave various phenomenological and cosmological consequences, to be probed in the future. This is a new indication from QCD, which gives a new avenue to deeper understand QCD, and provides a new guideline to consider going beyond the SM.
△ Less
Submitted 25 October, 2023; v1 submitted 27 December, 2021;
originally announced December 2021.
-
HMC with Normalizing Flows
Authors:
Sam Foreman,
Taku Izubuchi,
Luchang Jin,
Xiao-Yong Jin,
James C. Osborn,
Akio Tomiya
Abstract:
We propose using Normalizing Flows as a trainable kernel within the molecular dynamics update of Hamiltonian Monte Carlo (HMC). By learning (invertible) transformations that simplify our dynamics, we can outperform traditional methods at generating independent configurations. We show that, using a carefully constructed network architecture, our approach can be easily scaled to large lattice volume…
▽ More
We propose using Normalizing Flows as a trainable kernel within the molecular dynamics update of Hamiltonian Monte Carlo (HMC). By learning (invertible) transformations that simplify our dynamics, we can outperform traditional methods at generating independent configurations. We show that, using a carefully constructed network architecture, our approach can be easily scaled to large lattice volumes with minimal retraining effort. The source code for our implementation is publicly available online at https://github.com/nftqcd/fthmc.
△ Less
Submitted 2 December, 2021;
originally announced December 2021.
-
Correlated Dirac Eigenvalues and Axial Anomaly in Chiral Symmetric QCD
Authors:
Heng-Tong Ding,
Sheng-Tai Li,
Swagato Mukherjee,
Akio Tomiya,
Xiao-Dan Wang,
Yu Zhang
Abstract:
We investigate the Dirac eigenvalue spectrum ($ρ(λ,m_l)$) to study the microscopic origin of axial anomaly in high temperature phase of QCD. We propose novel relations between the derivatives ($\partial^n ρ(λ,m_l)/\partial m_l^n$) of the Dirac eigenvalue spectrum with respect to the quark mass ($m_l$) and the $(n+1)$-point correlations among the eigenvalues ($λ$) of the massless Dirac operator. Ba…
▽ More
We investigate the Dirac eigenvalue spectrum ($ρ(λ,m_l)$) to study the microscopic origin of axial anomaly in high temperature phase of QCD. We propose novel relations between the derivatives ($\partial^n ρ(λ,m_l)/\partial m_l^n$) of the Dirac eigenvalue spectrum with respect to the quark mass ($m_l$) and the $(n+1)$-point correlations among the eigenvalues ($λ$) of the massless Dirac operator. Based on these relations, we present lattice QCD results for $\partial^n ρ(λ,m_l)/\partial m_l^n$ ($n=1, 2, 3$) with $m_l$ corresponding to pion masses $m_π=160-55$ MeV, and at a temperature of about 1.6 times the chiral phase transition temperature. Calculations were carried out using (2+1)-flavors of highly improved staggered quarks and the tree-level Symanzik gauge action with the physical strange quark mass, three lattice spacings $a=0.12, 0.08, 0.06$ fm, and lattices having aspect ratios $4-9$. We find that $ρ(λ\to0,m_l)$ develops a peaked structure. This peaked structure, which arises due to non-Poisson correlations within the infrared part of the Dirac eigenvalue spectrum, becomes sharper as $a\to0$, and its amplitude is proportional to $m_l^2$. After continuum and chiral extrapolations, we find that the axial anomaly remains manifested in two-point correlation functions of scalar and pseudo-scalar mesons in the chiral limit. We demonstrate that the behavior of $ρ(λ\to0,m_l)$ is responsible for it.
△ Less
Submitted 7 December, 2021; v1 submitted 1 December, 2021;
originally announced December 2021.
-
New interpretation of chiral phase transition: Violation of trilemma in QCD
Authors:
Chuan-Xin Cui,
Jin-Yang Li,
Shinya Matsuzaki,
Mamiya Kawaguchi,
Akio Tomiya
Abstract:
We find that the chiral phase transition (chiral crossover) in QCD at physical point is triggered by big imbalance among three fundamental quantities essential for the QCD vacuum structure: susceptibility functions for the chiral symmetry, axial symmetry, and the topological charge. The balance, dobbed the QCD trilemma, is unavoidably violated when one of the magnitudes among them is highly domina…
▽ More
We find that the chiral phase transition (chiral crossover) in QCD at physical point is triggered by big imbalance among three fundamental quantities essential for the QCD vacuum structure: susceptibility functions for the chiral symmetry, axial symmetry, and the topological charge. The balance, dobbed the QCD trilemma, is unavoidably violated when one of the magnitudes among them is highly dominated, or suppressed. Based on a three-flavor Nambu-Jona-Lasinio model, we explicitly evaluate the amount of violation of the QCD trilemma at physical point, and show that the violation takes place not only at vacuum, but even in a whole temperature regime including the chiral crossover epoch. This work confirms and extends the suggestion recently reported from lattice QCD with 2 flavors on dominance of the axial and topological susceptibilities left in the chiral susceptibility at high temperatures. It turns out that the imbalance is essentially due to the flavor symmetry violation of the lightest three flavors, and the flavor breaking specifically brings enhancement of the axial anomaly contribution in the chiral order parameter, while the the strength of the axial breaking and the transition rate of the topological charge are fairly insensitive to the flavor symmetry. The violation of QCD trilemma and its flavor dependence can be tested by lattice simulations with 2 + 1 flavors in the future, and would also give a new guiding principle to explore the flavor dependence of the chiral phase transition, such as the Columbia plot, including possible extension with external fields.
△ Less
Submitted 17 June, 2022; v1 submitted 10 June, 2021;
originally announced June 2021.
-
Gauge covariant neural network for 4 dimensional non-abelian gauge theory
Authors:
Yuki Nagai,
Akio Tomiya
Abstract:
Quantum-Chromo dynamics (QCD) is a fundamental theory for quarks and gluons, which describes both the sub-atomic world and the history of our universe. The simulation for QCD on a lattice (lattice QCD) is one of the most challenging computational tasks. Recently, machine-learning techniques have been applied to solve various problems in lattice QCD. We propose gauge covariant neural networks and t…
▽ More
Quantum-Chromo dynamics (QCD) is a fundamental theory for quarks and gluons, which describes both the sub-atomic world and the history of our universe. The simulation for QCD on a lattice (lattice QCD) is one of the most challenging computational tasks. Recently, machine-learning techniques have been applied to solve various problems in lattice QCD. We propose gauge covariant neural networks and their training rule for lattice QCD, which can treat realistic quarks and gluons in four dimensions. We find that the smearing procedure can be regarded as extended versions of residual neural networks with fixed parameters. To show the applicability of our neural networks, we develop the self-learning hybrid Monte-Carlo for two-color QCD, where results are consistent with the results of the Hybrid Monte Carlo.
△ Less
Submitted 20 March, 2023; v1 submitted 22 March, 2021;
originally announced March 2021.
-
Detecting scale anomaly in chiral phase transition of QCD: new critical endpoint pinned down
Authors:
Mamiya Kawaguchi,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
Violation of scale symmetry, scale anomaly, being a radical concept in quantum field theory, is of importance to comprehend the vacuum structure of QCD, and should potentially contribute to the chiral phase transition in thermal QCD, as well as the chiral and U(1) axial symmetry. Though it should be essential, direct evidence of scale anomalies has never been observed in the chiral phase transitio…
▽ More
Violation of scale symmetry, scale anomaly, being a radical concept in quantum field theory, is of importance to comprehend the vacuum structure of QCD, and should potentially contribute to the chiral phase transition in thermal QCD, as well as the chiral and U(1) axial symmetry. Though it should be essential, direct evidence of scale anomalies has never been observed in the chiral phase transition. We propose a methodology to detect a scale anomaly in the chiral phase transition, which is an electromagnetically induced scale anomaly: apply a weak magnetic field background onto two-flavor massless QCD with an extremely heavy strange quark, first observe the chiral crossover; second, adjusting the strange quark mass to be smaller and smaller, observe the second-order chiral phase transition, and then the first-order one in the massless-three flavor limit. Thus, the second-order chiral phase transition, observed as the evidence of the quantum scale anomaly, is a new critical endpoint. It turns out that this electromagnetic scale anomaly gets most operative in the weak magnetic field regime, rather than a strong field region. We also briefly address accessibility of lattice QCD, a prospected application to dense matter system, and implications to astrophysical observations, such as gravitational wave productions provided from thermomagnetic QCD-like theories.
△ Less
Submitted 8 December, 2021; v1 submitted 10 February, 2021;
originally announced February 2021.
-
QCD phase structure in strong magnetic fields
Authors:
H. -T. Ding,
S. -T. Li,
Q. Shi,
A. Tomiya,
X. -D. Wang,
Y. Zhang
Abstract:
In this proceedings we discuss the natural connection between the reduction of neutral pion mass in the vacuum, and the magnetic catalysis as well as the reduction of transition temperature in the external magnetic field. We also present the first results on fluctuations of and correlations among conserved charges in strong magnetic fields from lattice QCD computations.
In this proceedings we discuss the natural connection between the reduction of neutral pion mass in the vacuum, and the magnetic catalysis as well as the reduction of transition temperature in the external magnetic field. We also present the first results on fluctuations of and correlations among conserved charges in strong magnetic fields from lattice QCD computations.
△ Less
Submitted 9 November, 2020;
originally announced November 2020.
-
Correlated Dirac Eigenvalues and Axial Anomaly in Chiral Symmetric QCD
Authors:
H. -T. Ding,
S. -T. Li,
Swagato Mukherjee,
A. Tomiya,
X. -D. Wang,
Y. Zhang
Abstract:
We introduce novel relations between the derivatives [$\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$] of the Dirac eigenvalue spectrum [$ρ(λ, m_{l})$] with respect to the light sea quark mass ($m_{l}$) and the $(n+1)$-point correlations among the eigenvalues ($λ$) of the massless Dirac operator. Using these relations we present lattice QCD results for $\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$ (…
▽ More
We introduce novel relations between the derivatives [$\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$] of the Dirac eigenvalue spectrum [$ρ(λ, m_{l})$] with respect to the light sea quark mass ($m_{l}$) and the $(n+1)$-point correlations among the eigenvalues ($λ$) of the massless Dirac operator. Using these relations we present lattice QCD results for $\partial^{n}ρ(λ, m_{l})/\partial m_{l}^n$ ($n=1, 2, 3$) for $m_{l}$ corresponding to pion masses $m_π=160-55$ MeV, and at a temperature of about 1.6 times the chiral phase transition temperature. Calculations were carried out using (2+1) flavors of highly improved staggered quarks with the physical value of strange quark mass, three lattice spacings $a=0.12, 0.08, 0.06$ fm, and lattices having aspect ratios $4-9$. We find that $ρ(λ\to0, m_{l})$ develops a peaked structure. This peaked structure arises due to non-Poisson correlations within the infrared part of the Dirac eigenvalue spectrum, becomes sharper as $a\to0$, and its amplitude is proportional to $m_{l}^2$. We demonstrate that this $ρ(λ\to0,m_l)$ is responsible for the manifestations of axial anomaly in two-point correlation functions of light scalar and pseudoscalar mesons. After continuum and chiral extrapolations we find that axial anomaly remains manifested in two-point correlation functions of scalar and pseudoscalar mesons in the chiral limit.
△ Less
Submitted 24 February, 2021; v1 submitted 28 October, 2020;
originally announced October 2020.
-
Self-learning Monte-Carlo for non-abelian gauge theory with dynamical fermions
Authors:
Yuki Nagai,
Akinori Tanaka,
Akio Tomiya
Abstract:
In this paper, we develop the self-learning Monte-Carlo (SLMC) algorithm for non-abelian gauge theory with dynamical fermions in four dimensions to resolve the autocorrelation problem in lattice QCD. We perform simulations with the dynamical staggered fermions and plaquette gauge action by both in HMC and SLMC for zero and finite temperature to examine the validity of SLMC. We confirm that SLMC ca…
▽ More
In this paper, we develop the self-learning Monte-Carlo (SLMC) algorithm for non-abelian gauge theory with dynamical fermions in four dimensions to resolve the autocorrelation problem in lattice QCD. We perform simulations with the dynamical staggered fermions and plaquette gauge action by both in HMC and SLMC for zero and finite temperature to examine the validity of SLMC. We confirm that SLMC can reduce autocorrelation time in non-abelian gauge theory and reproduces results from HMC. For finite temperature runs, we confirm that SLMC reproduces correct results with HMC, including higher-order moments of the Polyakov loop and the chiral condensate. Besides, our finite temperature calculations indicate that four flavor QC${}_2$D with $\hat{m} = 0.5$ is likely in the crossover regime in the Colombia plot.
△ Less
Submitted 22 October, 2020;
originally announced October 2020.
-
Chiral properties of (2+1)-flavor QCD in strong magnetic fields at zero temperature
Authors:
H. -T. Ding,
S. -T. Li,
A. Tomiya,
X. -D. Wang,
Y. Zhang
Abstract:
We present lattice QCD results for masses and magnetic polarizabilities of light and strange pseudoscalar mesons, chiral condensates, decay constants of neutral pion, and neutral kaon in the presence of background magnetic fields with $eB$ ranging up to around 3.35 GeV$^2$ ($\sim70~M_π^2$) in the vacuum. The computations were carried out in (2+1)-flavor QCD mostly on $32^3 \times 96$ lattices usin…
▽ More
We present lattice QCD results for masses and magnetic polarizabilities of light and strange pseudoscalar mesons, chiral condensates, decay constants of neutral pion, and neutral kaon in the presence of background magnetic fields with $eB$ ranging up to around 3.35 GeV$^2$ ($\sim70~M_π^2$) in the vacuum. The computations were carried out in (2+1)-flavor QCD mostly on $32^3 \times 96$ lattices using the highly improved staggered quark action with $M_π \approx $ 220 MeV at zero temperature. We find that the masses of neutral pseudoscalar mesons monotonously decrease as the magnetic field strength grows and then saturate at a nonzero value, while there exists a nonmonotonous behavior of charged pion and kaon masses in the magnetic field. We observe a $qB$ scaling of the up and down quark flavor components of neutral pion mass, neutral pion decay constant as well as the quark chiral condensates at 0.05 $\lesssim eB\lesssim$ 3.35 GeV$^2$. We show that the correction to the Gell-Mann-Oakes-Renner relation involving the neutral pion is less than 6% and the correction for the relation involving neutral kaon is less than 30% at $eB\lesssim$ 3.35 GeV$^2$. We also derive the Ward-Takahashi identities for QCD in the magnetic field in the continuum formulation including the relation between integrated neutral pseudoscalar meson correlators and chiral condensates.
△ Less
Submitted 6 July, 2021; v1 submitted 2 August, 2020;
originally announced August 2020.
-
Chiral phase structure of three flavor QCD in a background magnetic field
Authors:
Heng-Tong Ding,
Christian Schmidt,
Akio Tomiya,
Xiao-Dan Wang
Abstract:
We investigate the chiral phase structure of three flavor QCD in a background $U(1)$ magnetic field using the standard staggered action and the Wilson plaquette gauge action. We perform simulations on lattices with a temporal extent of $N_τ=4$ and four spatial extents of $N_σ= 8,16, 20$ and 24. We choose a smaller-than-physical quark mass in lattice spacing as $am = 0.030$ such that there exists a…
▽ More
We investigate the chiral phase structure of three flavor QCD in a background $U(1)$ magnetic field using the standard staggered action and the Wilson plaquette gauge action. We perform simulations on lattices with a temporal extent of $N_τ=4$ and four spatial extents of $N_σ= 8,16, 20$ and 24. We choose a smaller-than-physical quark mass in lattice spacing as $am = 0.030$ such that there exists a crossover transition at vanishing magnetic fields, and adopt two values of magnetic field strength in lattice spacing $a \sqrt{ e{B}}\simeq 1.5$ and 2. We find that the transition becomes stronger in the presence of a background magnetic field, and turns into a first order as seen from the volume scaling of the order parameter susceptibility as well as the metastable states in the time history of the chiral condensate. On the other hand, the chiral condensate and transition temperature always increase with $B$ even within the regime of a first order phase transition. This suggests that the discrepancy in the behavior of chiral condensates and transition temperature as a function of $B$ between earlier lattice studies using larger-than-physical pion masses with standard staggered fermions and those using physical pions with improved staggered fermions is mainly due to lattice cutoff effects.
△ Less
Submitted 18 August, 2020; v1 submitted 23 June, 2020;
originally announced June 2020.
-
Analysis on nonperturbative flavor violation at chiral crossover criticality in QCD
Authors:
Mamiya Kawaguchi,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
We discuss the violation of quark-flavor symmetry at high temperatures, induced from nonperturbative thermal loop corrections and axial anomaly, based on a three-flavor linear-sigma model including an axial-anomaly induced-flavor breaking term. We employ a nonperturbative analysis following the Cornwall-Jackiw-Tomboulis formalism, and show that the model undergoes a chiral crossover with a pseudo-…
▽ More
We discuss the violation of quark-flavor symmetry at high temperatures, induced from nonperturbative thermal loop corrections and axial anomaly, based on a three-flavor linear-sigma model including an axial-anomaly induced-flavor breaking term. We employ a nonperturbative analysis following the Cornwall-Jackiw-Tomboulis formalism, and show that the model undergoes a chiral crossover with a pseudo-critical temperature, consistently with lattice observations. We find following features regarding the flavor breaking eminent around and above the pseudo-critical temperature: i) up-and down-quark condensates drop faster than the strange quark's toward the criticality, but still keep nonzero value even going far above the critical temperature; ii) the introduced anomaly-related flavor-breaking effect acts as a catalyzer toward the chiral restoration, and reduces the amount of flavor breaking in the up, down and strange quark condensates; iii) a dramatic deformation for the meson flavor mixing structure is observed, in which the anomaly-induced favor breaking is found to be almost irrelevant; iv) the meson spectroscopy gets corrected by the net nonperturbative flavor breaking effects, where the scalar meson mass hierarchy (inverse mass hierarchy) is significantly altered by the presence of the anomaly-related flavor breaking; v) the topological susceptibility significantly gets the contribution from the surviving strange quark condensate, which cannot be dictated by the chiral perturbation theory, and deviates from the dilute instanton gas prediction. There the anomaly-induced flavor breaking plays a role of the destructive interference for the net flavor violation; vi) the U(1)_A breaking is enhanced by the strange quark condensate, which may account for the tension in the effective U(1)_A restoration observed on lattices with two flavors and 2+1 flavors near the chiral limit.
△ Less
Submitted 4 January, 2021; v1 submitted 14 May, 2020;
originally announced May 2020.
-
Nonperturbative Flavor Breaking in Topological Susceptibility at Chiral Crossover
Authors:
Mamiya Kawaguchi,
Shinya Matsuzaki,
Akio Tomiya
Abstract:
We demonstrate that the QCD topological susceptibility nonperturbatively gets a significant contribution signaled by flavor-nonuniversal quark condensates at around the pseudo-critical temperature of the chiral crossover. It implies a remarkable flavor breaking in the axial anomaly as well as the QCD theta vacuum in high temperature QCD, which are almost flavor universal in the vacuum. A nontrivia…
▽ More
We demonstrate that the QCD topological susceptibility nonperturbatively gets a significant contribution signaled by flavor-nonuniversal quark condensates at around the pseudo-critical temperature of the chiral crossover. It implies a remarkable flavor breaking in the axial anomaly as well as the QCD theta vacuum in high temperature QCD, which are almost flavor universal in the vacuum. A nontrivial flavor breaking is triggered by nonperturbative thermal loop corrections at around the chiral crossover, which is different from the trivial flavor violation just scaled by the quark mass ratio, observed at asymptotically high temperatures. This critical flavor violation cannot be dictated by the chiral perturbation theory with that lattice QCD usually compares, or the dilute instanton gas approximation based on that its astrophysical implications have conventionally been made. This would give an impact on the thermal history and the cosmological evolution of QCD axion including the estimate of the relic abundance as a dark matter candidate.
△ Less
Submitted 4 January, 2021; v1 submitted 25 March, 2020;
originally announced March 2020.
-
Meson masses in external magnetic fields with HISQ fermions
Authors:
Heng-Tong Ding,
Sheng-Tai Li,
Swagato Mukherjee,
Akio Tomiya,
Xiao-Dan Wang
Abstract:
We studied the temporal correlation function of mesons in the pseudo-scalar channel in (2+1)-flavor QCD in the presence of external magnetic fields at zero temperature. The simulations were performed on $32^3 \times 96$ lattices using the Highly Improved Staggered Quarks (HISQ) action with $m_π \approx $ 230 MeV. The strength of magnetic fields $|eB|$ ranges from 0 to around 3.3 GeV$^2$ (…
▽ More
We studied the temporal correlation function of mesons in the pseudo-scalar channel in (2+1)-flavor QCD in the presence of external magnetic fields at zero temperature. The simulations were performed on $32^3 \times 96$ lattices using the Highly Improved Staggered Quarks (HISQ) action with $m_π \approx $ 230 MeV. The strength of magnetic fields $|eB|$ ranges from 0 to around 3.3 GeV$^2$ ($\sim 60 m_π^2$). We found that the masses of neutral pseudo-scalar particles, e.g. neutral pion and kaon, monotonouslly decrease as the magnetic field grows and then saturate at a nonzero value. It is observed that heavier neutral pseudo-scalars are less affected by magnetic fields. Moreover, we found a non-monotonous behavior of charged pion and kaon mass in magnetic field for the first time. In the case of small magnetic field (0 $\leq~|eB| \lesssim$ 0.3 GeV$^2~\sim 6m_π^2$ ) the mass of charged pseudo-scalar grows with magnetic field and can be well described by the Lowest Landau Level approximation, while for $|eB|$ larger than 0.3 GeV$^2$ the mass starts to decrease. The possible connection between $|eB|$ dependences of neutral pion mass and the decreasing behavior of pseudo-critical temperature in magnetic field is discussed. Due to the nonzero value of neutral pion mass our simulation indicates that the superconducting phase of QCD does not exist in the current window of magnetic field.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
Dirac Eigenvalue spectrum of $N_f$=2+1 QCD towards the chiral limit using HISQ fermions
Authors:
Heng-Tong Ding,
Olaf Kaczmarek,
Frithjof Karsch,
Sheng-Tai Li,
Swagato Mukherjee,
Akio Tomiya,
Yu Zhang
Abstract:
We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number to determine the eigenvalue spectrum. Simulations of (2 + 1)-flavor QCD are performed using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ$ = 8 lattices with aspect ratios $N_σ/N_τ$ ranging from 5 to 7. The strange quark mass is fixed to its physical value $m_{s}^{\rm phy}$, and…
▽ More
We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number to determine the eigenvalue spectrum. Simulations of (2 + 1)-flavor QCD are performed using the Highly Improved Staggered Quarks (HISQ/tree) action on $N_τ$ = 8 lattices with aspect ratios $N_σ/N_τ$ ranging from 5 to 7. The strange quark mass is fixed to its physical value $m_{s}^{\rm phy}$, and the light quark masses $m_{l}$ are varied from $m_{s}^{\rm phy}/40$ to $m_{s}^{\rm phy}/160$ which correspond to pion mass $m_π$ ranging from 110 MeV to 55 MeV in the continuum limit. We compute the chiral condensate and $χ_π - χ_δ$ through the eigenvalue spectrum obtained from the the eigenvalue filtering method. We compare these results with those obtained from a direct calculation of the observables which involves inversions of the fermion matrix using the stochastic "noise vector" method. We find that these approaches yield consistent results. Furthermore, we also investigate the quark mass and temperature dependences of the Dirac eigenvalue density at zero eigenvalues to gain more insights about the $U_A(1)$ symmetry breaking in QCD.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
Classically Emulated Digital Quantum Simulation of the Schwinger Model with Topological Term via Adiabatic State Preparation
Authors:
Bipasha Chakraborty,
Masazumi Honda,
Taku Izubuchi,
Yuta Kikuchi,
Akio Tomiya
Abstract:
We perform a digital quantum simulation of a gauge theory with a topological term in Minkowski spacetime, which is practically inaccessible by standard lattice Monte Carlo simulations. We focus on $1+1$ dimensional quantum electrodynamics with the $θ$-term known as the Schwinger model. We construct the true vacuum state of a lattice Schwinger model using adiabatic state preparation which, in turn,…
▽ More
We perform a digital quantum simulation of a gauge theory with a topological term in Minkowski spacetime, which is practically inaccessible by standard lattice Monte Carlo simulations. We focus on $1+1$ dimensional quantum electrodynamics with the $θ$-term known as the Schwinger model. We construct the true vacuum state of a lattice Schwinger model using adiabatic state preparation which, in turn, allows us to compute an expectation value of the fermion mass operator with respect to the vacuum. Upon taking a continuum limit we find that our result in massless case agrees with the known exact result. In massive case, we find an agreement with mass perturbation theory in small mass regime and deviations in large mass regime. We estimate computational costs required to take a reasonable continuum limit. Our results imply that digital quantum simulation is already useful tool to explore non-perturbative aspects of gauge theories with real time and topological terms.
△ Less
Submitted 22 April, 2022; v1 submitted 2 January, 2020;
originally announced January 2020.
-
Entanglement after Quantum Quenches in Lifshitz Scalar Theories
Authors:
Keun-Young Kim,
Mitsuhiro Nishida,
Masahiro Nozaki,
Minsik Seo,
Yuji Sugimoto,
Akio Tomiya
Abstract:
We study the time evolution of the entanglement entropy after quantum quenches in Lifshitz free scalar theories, with the dynamical exponent $z>1$, by using the correlator method. For quantum quenches we consider two types of time-dependent mass functions: end-critical-protocol (ECP) and cis-critical-protocol (CCP). In both cases, at early times the entanglement entropy is independent of the subsy…
▽ More
We study the time evolution of the entanglement entropy after quantum quenches in Lifshitz free scalar theories, with the dynamical exponent $z>1$, by using the correlator method. For quantum quenches we consider two types of time-dependent mass functions: end-critical-protocol (ECP) and cis-critical-protocol (CCP). In both cases, at early times the entanglement entropy is independent of the subsystem size. After a critical time ($t_c$), the entanglement entropy starts depending on the subsystem size significantly. This critical time $t_c$ for $z = 1$ in the fast ECP and CCP has been explained well by the fast quasi-particle of the quasi-particle picture. However, we find that for $z > 1$ this explanation does not work and $t_c$ is delayed. We explain why $t_c$ is delayed for $z>1$ based on the quasiparticle picture: in essence, it is due to the competition between the fast and slow quasiparticles. At late times, in the ECP, the entanglement entropy slowly increases while, in the CCP, it is oscillating with a well defined period by the final mass scale, independently of $z$. We give an interpretation of this phenomena by the correlator method. As $z$ increases, the entanglement entropy increases, which can be understood by long-range interactions due to $z$.
△ Less
Submitted 30 January, 2020; v1 submitted 13 June, 2019;
originally announced June 2019.
-
Phase structure of three flavor QCD in external magnetic fields using HISQ fermions
Authors:
Akio Tomiya,
Heng-Tong Ding,
Xiao-Dan Wang,
Yu Zhang,
Swagato Mukherjee,
Christian Schmidt
Abstract:
We study the phase structure of QCD with three degenerate flavors in the external magnetic fields using highly improved staggered quarks (HISQ). The simulations are performed on $16^3\times 6$ lattice. In order to investigate the quark mass dependence of the chiral transition we choose the values of the bare quark masses 0.015 and 0.0009375 in the lattice unit, corresponding to $m_π=320$ MeV and 8…
▽ More
We study the phase structure of QCD with three degenerate flavors in the external magnetic fields using highly improved staggered quarks (HISQ). The simulations are performed on $16^3\times 6$ lattice. In order to investigate the quark mass dependence of the chiral transition we choose the values of the bare quark masses 0.015 and 0.0009375 in the lattice unit, corresponding to $m_π=320$ MeV and 80 MeV in the continuum limit. We found no indication of a first order phase transition in the current window of quark masses and external magnetic fields. Unlike to the case with standard staggered fermions inverse magnetic catalysis is always observed at about the critical temperature. The microscopic origin of this phenomena are further discussed by looking into the Dirac eigenvalue spectrum.
△ Less
Submitted 2 April, 2019;
originally announced April 2019.
-
Phase transition encoded in neural network
Authors:
Kouji Kashiwa,
Yuta Kikuchi,
Akio Tomiya
Abstract:
We discuss an aspect of neural networks for the purpose of phase transition detection. To this end, we first train the neural network by feeding Ising/Potts configurations with labels of temperature so that it can predict the temperature of input. We do not explicitly supervise whether the configurations are in ordered/disordered phase. Nevertheless, we can identify the critical temperature from t…
▽ More
We discuss an aspect of neural networks for the purpose of phase transition detection. To this end, we first train the neural network by feeding Ising/Potts configurations with labels of temperature so that it can predict the temperature of input. We do not explicitly supervise whether the configurations are in ordered/disordered phase. Nevertheless, we can identify the critical temperature from the parameters (weights and biases) of trained neural network. We attempt to understand how temperature-supervised neural networks capture the information of phase transition by paying attention to what quantities they learn. Our detailed analyses reveal that they learn different physical quantities depending on how well they are trained. Main observation in this study is how the weights in the trained neural-network can have information of the phase transition in addition to temperature.
△ Less
Submitted 26 June, 2019; v1 submitted 4 December, 2018;
originally announced December 2018.
-
Deep Learning and Holographic QCD
Authors:
Koji Hashimoto,
Sotaro Sugishita,
Akinori Tanaka,
Akio Tomiya
Abstract:
We apply the relation between deep learning (DL) and the AdS/CFT correspondence to a holographic model of QCD. Using a lattice QCD data of the chiral condensate at a finite temperature as our training data, the deep learning procedure holographically determines an emergent bulk metric as neural network weights. The emergent bulk metric is found to have both a black hole horizon and a finite-height…
▽ More
We apply the relation between deep learning (DL) and the AdS/CFT correspondence to a holographic model of QCD. Using a lattice QCD data of the chiral condensate at a finite temperature as our training data, the deep learning procedure holographically determines an emergent bulk metric as neural network weights. The emergent bulk metric is found to have both a black hole horizon and a finite-height IR wall, so shares both the confining and deconfining phases, signaling the cross-over thermal phase transition of QCD. In fact, a quark antiquark potential holographically calculated by the emergent bulk metric turns out to possess both the linear confining part and the Debye screening part, as is often observed in lattice QCD. From this we argue the discrepancy between the chiral symmetry breaking and the quark confinement in the holographic QCD. The DL method is shown to provide a novel data-driven holographic modeling of QCD, and sheds light on the mechanism of emergence of the bulk geometries in the AdS/CFT correspondence.
△ Less
Submitted 27 September, 2018;
originally announced September 2018.
-
Deep Learning and AdS/CFT
Authors:
Koji Hashimoto,
Sotaro Sugishita,
Akinori Tanaka,
Akio Tomiya
Abstract:
We present a deep neural network representation of the AdS/CFT correspondence, and demonstrate the emergence of the bulk metric function via the learning process for given data sets of response in boundary quantum field theories. The emergent radial direction of the bulk is identified with the depth of the layers, and the network itself is interpreted as a bulk geometry. Our network provides a dat…
▽ More
We present a deep neural network representation of the AdS/CFT correspondence, and demonstrate the emergence of the bulk metric function via the learning process for given data sets of response in boundary quantum field theories. The emergent radial direction of the bulk is identified with the depth of the layers, and the network itself is interpreted as a bulk geometry. Our network provides a data-driven holographic modeling of strongly coupled systems. With a scalar $φ^4$ theory with unknown mass and coupling, in unknown curved spacetime with a black hole horizon, we demonstrate our deep learning (DL) framework can determine them which fit given response data. First, we show that, from boundary data generated by the AdS Schwarzschild spacetime, our network can reproduce the metric. Second, we demonstrate that our network with experimental data as an input can determine the bulk metric, the mass and the quadratic coupling of the holographic model. As an example we use the experimental data of magnetic response of a strongly correlated material Sm$_{0.6}$Sr$_{0.4}$MnO$_3$. This AdS/DL correspondence not only enables gravity modeling of strongly correlated systems, but also sheds light on a hidden mechanism of the emerging space in both AdS and DL.
△ Less
Submitted 22 February, 2018;
originally announced February 2018.
-
Entanglement Spreading and Oscillation
Authors:
Mitsuhiro Nishida,
Masahiro Nozaki,
Yuji Sugimoto,
Akio Tomiya
Abstract:
We study dynamics of quantum entanglement in smooth global quenches with a finite rate, by computing the time evolution of entanglement entropy in 1 + 1 dimensional free scalar theory with time-dependent masses which start from a nonzero value at early time and either crosses or approaches zero. The time-dependence is chosen so that the quantum dynamics is exactly solvable. If the quenches asympto…
▽ More
We study dynamics of quantum entanglement in smooth global quenches with a finite rate, by computing the time evolution of entanglement entropy in 1 + 1 dimensional free scalar theory with time-dependent masses which start from a nonzero value at early time and either crosses or approaches zero. The time-dependence is chosen so that the quantum dynamics is exactly solvable. If the quenches asymptotically approach a critical point at late time, the early-time and late-time entropies are proportional to the time and subsystem size respectively. Their proportionality coefficients are determined by scales: in a fast limit, an initial correlation length; in a slow limit, an effective scale defined when adiabaticity breaks down. If the quenches cross a critical point, the time evolution of entropy is characterized by the scales: the initial correlation length in the fast limit and the effective correlation length in the slow limit. The entropy oscillates, and the entanglement oscillation comes from a coherence between right-moving and left-moving waves if we measure the entropy after time characterized by the quench rate. The periodicity of the late-time oscillation is consistent with the periodicity of the oscillation of zero modes which are zero-momentum spectra of two point functions of a fundamental field and its conjugate momentum.
△ Less
Submitted 12 January, 2018; v1 submitted 28 December, 2017;
originally announced December 2017.
-
Towards reduction of autocorrelation in HMC by machine learning
Authors:
Akinori Tanaka,
Akio Tomiya
Abstract:
In this paper we propose new algorithm to reduce autocorrelation in Markov chain Monte-Carlo algorithms for euclidean field theories on the lattice. Our proposing algorithm is the Hybrid Monte-Carlo algorithm (HMC) with restricted Boltzmann machine. We examine the validity of the algorithm by employing the phi-fourth theory in three dimension. We observe reduction of the autocorrelation both in sy…
▽ More
In this paper we propose new algorithm to reduce autocorrelation in Markov chain Monte-Carlo algorithms for euclidean field theories on the lattice. Our proposing algorithm is the Hybrid Monte-Carlo algorithm (HMC) with restricted Boltzmann machine. We examine the validity of the algorithm by employing the phi-fourth theory in three dimension. We observe reduction of the autocorrelation both in symmetric and broken phase as well. Our proposing algorithm provides consistent central values of expectation values of the action density and one-point Green's function with ones from the original HMC in both the symmetric phase and broken phase within the statistical error. On the other hand, two-point Green's functions have slight difference between one calculated by the HMC and one by our proposing algorithm in the symmetric phase. Furthermore, near the criticality, the distribution of the one-point Green's function differs from the one from HMC. We discuss the origin of discrepancies and its improvement.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Chiral phase transition of three flavor QCD with nonzero magnetic field using standard staggered fermions
Authors:
Akio Tomiya,
Heng-Tong Ding,
Swagato Mukherjee,
Christian Schmidt,
Xiao-Dan Wang
Abstract:
Lattice simulations for (2+1)-flavor QCD with external magnetic field demonstrated that the quark mass is one of the important parameters responsible for the (inverse) magnetic catalysis. We discuss the dependences of chiral condensates and susceptibilities, the Polyakov loop on the magnetic field and quark mass in three degenerate flavor QCD. The lattice simulations are performed using standard s…
▽ More
Lattice simulations for (2+1)-flavor QCD with external magnetic field demonstrated that the quark mass is one of the important parameters responsible for the (inverse) magnetic catalysis. We discuss the dependences of chiral condensates and susceptibilities, the Polyakov loop on the magnetic field and quark mass in three degenerate flavor QCD. The lattice simulations are performed using standard staggered fermions and the plaquette action with spatial sizes Ns = 16 and 24 and a fixed temporal size Nt = 4. The value of the quark masses are chosen such that the system undergoes a first order chiral phase transition and crossover with zero magnetic field. We find that in light mass regime, the quark chiral condensate undergoes magnetic catalysis in the whole temperature region and the phase transition tend to become stronger as the magnetic field increases. In crossover regime, deconfinement transition temperature is shifted by the magnetic field when quark mass ma is less than 0.4. The lattice cutoff effects are also discussed.
△ Less
Submitted 8 November, 2017;
originally announced November 2017.
-
(De)Constructing magnetized dimensions
Authors:
Yoshiyuki Tatsuta,
Akio Tomiya
Abstract:
We provide an origin of family replications in the standard model of particle physics by constructing renormalizable, asymptotically free, four dimensional local gauge theories that dynamically generate the fifth and sixth dimensions with magnetic fluxes.
We provide an origin of family replications in the standard model of particle physics by constructing renormalizable, asymptotically free, four dimensional local gauge theories that dynamically generate the fifth and sixth dimensions with magnetic fluxes.
△ Less
Submitted 15 March, 2017;
originally announced March 2017.
-
Quantum Quench and Scaling of Entanglement Entropy
Authors:
Pawel Caputa,
Sumit R. Das,
Masahiro Nozaki,
Akio Tomiya
Abstract:
Global quantum quench with a finite quench rate which crosses critical points is known to lead to universal scaling of correlation functions as functions of the quench rate. In this work, we explore scaling properties of the entanglement entropy of a subsystem in a harmonic chain during a mass quench which asymptotes to finite constant values at early and late times and for which the dynamics is e…
▽ More
Global quantum quench with a finite quench rate which crosses critical points is known to lead to universal scaling of correlation functions as functions of the quench rate. In this work, we explore scaling properties of the entanglement entropy of a subsystem in a harmonic chain during a mass quench which asymptotes to finite constant values at early and late times and for which the dynamics is exactly solvable. When the initial state is the ground state, we find that for large enough subsystem sizes the entanglement entropy becomes independent of size. This is consistent with Kibble-Zurek scaling for slow quenches, and with recently discussed "fast quench scaling" for quenches fast compared to physical scales, but slow compared to UV cutoff scales.
△ Less
Submitted 28 May, 2017; v1 submitted 14 February, 2017;
originally announced February 2017.
-
Evidence of effective axial U(1) symmetry restoration at high temperature QCD
Authors:
A. Tomiya,
G. Cossu,
S. Aoki,
H. Fukaya,
S. Hashimoto,
T. Kaneko,
J. Noaki
Abstract:
We study the axial U(1) symmetry at finite temperature in two-flavor lattice QCD. Employing the Mobius domain-wall fermions, we generate gauge configurations slightly above the critical temperature Tc with different lattice sizes L = 2-4 fm. Our action allows frequent topology tunneling while keeping good chiral symmetry close enough to that of overlap fermions. This allows us to recover full chir…
▽ More
We study the axial U(1) symmetry at finite temperature in two-flavor lattice QCD. Employing the Mobius domain-wall fermions, we generate gauge configurations slightly above the critical temperature Tc with different lattice sizes L = 2-4 fm. Our action allows frequent topology tunneling while keeping good chiral symmetry close enough to that of overlap fermions. This allows us to recover full chiral symmetry by an overlap/domain-wall reweighting. Above the phase transition, a strong suppression of the low-lying modes is observed in both of overlap and domain-wall Dirac spectra. We, however, find a sizable violation of the Ginsparg-Wilson relation in the Mobius domain-wall Dirac eigenmodes, which dominates the signals of the axial U(1) symmetry breaking near the chiral limit. We also find that the use of overlap fermion only in the valence sector is dangerous since it suffers from the artifacts due to partial quenching. Reweighting the Mobius domain-wall fermion determinant to that of the overlap fermion, we observe the axial U(1) breaking to vanish in the chiral limit, which is stable against the changes of the lattice volume and lattice spacing.
△ Less
Submitted 26 October, 2017; v1 submitted 6 December, 2016;
originally announced December 2016.
-
Detection of phase transition via convolutional neural network
Authors:
Akinori Tanaka,
Akio Tomiya
Abstract:
We design a Convolutional Neural Network (CNN) which studies correlation between discretized inverse temperature and spin configuration of 2D Ising model and show that it can find a feature of the phase transition without teaching any a priori information for it. We also define a new order parameter via the CNN and show that it provides well approximated critical inverse temperature. In addition,…
▽ More
We design a Convolutional Neural Network (CNN) which studies correlation between discretized inverse temperature and spin configuration of 2D Ising model and show that it can find a feature of the phase transition without teaching any a priori information for it. We also define a new order parameter via the CNN and show that it provides well approximated critical inverse temperature. In addition, we compare the activation functions for convolution layer and find that the Rectified Linear Unit (ReLU) is important to detect the phase transition of 2D Ising model.
△ Less
Submitted 24 June, 2017; v1 submitted 28 September, 2016;
originally announced September 2016.
-
On the axial $U(1)$ symmetry at finite temperature
Authors:
JLQCD Collaboration,
Guido Cossu,
Hidenori Fukaya,
Shoji Hashimoto,
Jun-ichi Noaki,
Akio Tomiya
Abstract:
We study the $U(1)_A$ anomaly in two-flavor lattice QCD at finite temperature with the Möbius domain-wall Dirac operator. We generate gauge configurations in the temperature range $(0.9, 1.2) T_c$ on different physical volumes, $L=$ 2--4 fm, and lattice spacings. We measure the difference of the susceptibilities of the flavor non-singlet scalar ($χ_δ$) and pseudoscalar ($χ_π$) mesons. They are rel…
▽ More
We study the $U(1)_A$ anomaly in two-flavor lattice QCD at finite temperature with the Möbius domain-wall Dirac operator. We generate gauge configurations in the temperature range $(0.9, 1.2) T_c$ on different physical volumes, $L=$ 2--4 fm, and lattice spacings. We measure the difference of the susceptibilities of the flavor non-singlet scalar ($χ_δ$) and pseudoscalar ($χ_π$) mesons. They are related by an axial $U(1)$ transformation and the difference vanishes if the axial $U(1)$ symmetry is respected. We identify the source of axial $U(1)$ symmetry breaking at finite temperature in the lowest eigenmodes, for the observable $χ_π- χ_δ$. We then reweight the Möbius domain-wall fermion partition function to that of the overlap-Dirac operator to fully recover chiral symmetry. Our data show a significant discrepancy in the results coming from the Möbius domain-wall valence quarks, the overlap valence quarks on our DWF configurations and the reweighted ones that have full chiral symmetry. After recovering full chiral symmetry we conclude that the difference $χ_π- χ_δ$ shows a suppression in the chiral limit that is compatible with an effective restoration of $U(1)_A$ at $T \gtrsim T_c$ in the scalar meson channels.
△ Less
Submitted 18 November, 2015;
originally announced November 2015.
-
Violation of chirality of the Möbius domain-wall Dirac operator from the eigenmodes
Authors:
JLQCD collaboration,
Guido Cossu,
Hidenori Fukaya,
Shoji Hashimoto,
Akio Tomiya
Abstract:
We investigate the effects of the violation of the Ginsparg-Wilson (GW) relation in the Möbius domain-wall fermion formulation on the lattice with finite fifth dimension. Using a decomposion in terms of the eigenmodes of its four-dimensional effective Dirac operator, we isolate the GW-violating terms for various physical quantities including the residual mass and the meson susceptibilities relevan…
▽ More
We investigate the effects of the violation of the Ginsparg-Wilson (GW) relation in the Möbius domain-wall fermion formulation on the lattice with finite fifth dimension. Using a decomposion in terms of the eigenmodes of its four-dimensional effective Dirac operator, we isolate the GW-violating terms for various physical quantities including the residual mass and the meson susceptibilities relevant for the effective restoration of the axial U(1) symmetry at finite temperature. Numerical result shows that the GW-violating effect is more significant, or even overwhelming, for the quantities that are dominated by the low-lying eigenmodes.
△ Less
Submitted 26 October, 2015;
originally announced October 2015.
-
Effective lattice action for the configurations smeared by the Wilson flow
Authors:
Aya Kagimura,
Akio Tomiya,
Ryo Yamamura
Abstract:
We investigate a trajectory for the Wilson flow in the theory space. For this purpose, we determine the coefficient of the plaquette and rectangular terms in the action for the configurations defined by the solution of the Wilson flow. The demon method regarded as one of the inverse Monte Carlo methods is used for the determination of them. Starting from the conventional Wilson plaquette action of…
▽ More
We investigate a trajectory for the Wilson flow in the theory space. For this purpose, we determine the coefficient of the plaquette and rectangular terms in the action for the configurations defined by the solution of the Wilson flow. The demon method regarded as one of the inverse Monte Carlo methods is used for the determination of them. Starting from the conventional Wilson plaquette action of quenched QCD, we find that the coefficient of the plaquette grows while that of the rectangular tends to negative with the development of the flow as the known improved actions. We also find that the trajectory forms a straight line in the two-coupling theory space.
△ Less
Submitted 27 August, 2015; v1 submitted 20 August, 2015;
originally announced August 2015.
-
Effects of near-zero Dirac eigenmodes on axial U(1) symmetry at finite temperature
Authors:
Akio Tomiya,
Guido Cossu,
Hidenori Fukaya,
Shoji Hashimoto,
Junichi Noaki
Abstract:
We study the axial U(1)A symmetry of Nf = 2 QCD at finite temperature using the Dirac eigenvalue spectrum. The gauge configurations are generated employing the Mobius domain-wall fermion action on 16^3x8 and 32^3x8 lattices. The physical spatial size of these lattices is around 2 fm and 4 fm, respectively, and the simulated temperature is around 200 MeV, which is slightly above the critical temper…
▽ More
We study the axial U(1)A symmetry of Nf = 2 QCD at finite temperature using the Dirac eigenvalue spectrum. The gauge configurations are generated employing the Mobius domain-wall fermion action on 16^3x8 and 32^3x8 lattices. The physical spatial size of these lattices is around 2 fm and 4 fm, respectively, and the simulated temperature is around 200 MeV, which is slightly above the critical temperature of the chiral phase transition. Although the Mobius domain-wall Dirac operator is expected to have a good chiral symmetry and our data actually show small values of the residual mass, we observe significant violation of the Ginsparg-Wilson relation for the low- lying eigenmodes of the Mobius domain-wall Dirac operator. Using the reweighting technique, we compute the overlap-Dirac operator spectrum on the same set of configurations and find a significant difference of the spectrum between the two Dirac operators for the low-lying eigenvalues. The overlap-Dirac spectrum shows a gap from zero, which is insensitive to the spacial volume.
△ Less
Submitted 23 December, 2014;
originally announced December 2014.
-
Axial U(1) symmetry at finite temperature with Möbius domain-wall fermions
Authors:
JLQCD Collaboration,
Guido Cossu,
Hidenori Fukaya,
Shoji Hashimoto,
Takashi Kaneko,
Jun-ichi Noaki,
Akio Tomiya
Abstract:
We investigate the axial $U(1)$ symmetry restoration at finite temperature in two flavor QCD. We employ the Möbius domain-wall formalism that is designed to achieve good chiral symmetry. We show the measurements of a difference of meson susceptibilities, sensitive to the $U(1)_A$ symmetry breaking. The signal is dominated by zero and near-zero modes. By reweighting the measure to that of overlap f…
▽ More
We investigate the axial $U(1)$ symmetry restoration at finite temperature in two flavor QCD. We employ the Möbius domain-wall formalism that is designed to achieve good chiral symmetry. We show the measurements of a difference of meson susceptibilities, sensitive to the $U(1)_A$ symmetry breaking. The signal is dominated by zero and near-zero modes. By reweighting the measure to that of overlap fermions we find a suppression of the $U(1)_A$ breaking effects above the chiral transition temperature.
△ Less
Submitted 17 December, 2014;
originally announced December 2014.
-
Determination of the mass anomalous dimension for $N_f=12$ and $N_f=9$ SU($3$) gauge theories
Authors:
Etsuko Itou,
Akio Tomiya
Abstract:
We show the numerical simulation result for the mass anomalous dimension of the SU($3$) gauge theory coupled to $N_f = 12$ fundamental fermions. We use two independent methods, namely the step scaling method and the hyperscaling method of the Dirac mode number, to determine the anomalous dimension in the vicinity of the infrared fixed point of the theory. We show the continuum extrapolations keepi…
▽ More
We show the numerical simulation result for the mass anomalous dimension of the SU($3$) gauge theory coupled to $N_f = 12$ fundamental fermions. We use two independent methods, namely the step scaling method and the hyperscaling method of the Dirac mode number, to determine the anomalous dimension in the vicinity of the infrared fixed point of the theory. We show the continuum extrapolations keeping the renormalized coupling constant as a reference in both analyses.
Furthermore, some recent works seems to suggest the lower boundary of the conformal window of the SU($3$) gauge theory exists between $N_f=8$ and $10$. We also briefly report our new project, in which the numerical simulation of the SU($3$) gauge theory coupled to $N_f=9$ fundamental fermions has been performed.
△ Less
Submitted 5 November, 2014;
originally announced November 2014.
-
Symmetry breaking caused by large R-charge
Authors:
Akinori Tanaka,
Akio Tomiya,
Takuya Shimotani
Abstract:
We discuss the gauge symmetry breaking via the Hosotani mechanism by using exact results on supersymmetric gauge theories based on the localization method. We use the theories on S^2 x S^1 Euclidean space, and study how the effective potential for the Wilson line phase varies by running an imaginary chemical potential. In order to break the symmetry, we find that large R-charge is necessary. With…
▽ More
We discuss the gauge symmetry breaking via the Hosotani mechanism by using exact results on supersymmetric gauge theories based on the localization method. We use the theories on S^2 x S^1 Euclidean space, and study how the effective potential for the Wilson line phase varies by running an imaginary chemical potential. In order to break the symmetry, we find that large R-charge is necessary. With such large R-charge, we study the phase structure of the theory. In addition, we observed that a finite size effect on our curved space when we take R-charge is not so large.
△ Less
Submitted 23 December, 2014; v1 submitted 30 April, 2014;
originally announced April 2014.
