Next Article in Journal
Projected-Reflected Subgradient-Extragradient Method and Its Real-World Applications
Next Article in Special Issue
Quark Cluster Expansion Model for Interpreting Finite-T Lattice QCD Thermodynamics
Previous Article in Journal
Effective Rotor Fault Diagnosis Model Using Multilayer Signal Analysis and Hybrid Genetic Binary Chicken Swarm Optimization
Previous Article in Special Issue
Empirical Consequences of Emergent Mass
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Perturbative RG Analysis of the Condensate Dependence of the Axial Anomaly in the Three-Flavor Linear Sigma Model

Institute of Physics, Eötvös University Budapest, H-1117 Budapest, Hungary
Symmetry 2021, 13(3), 488; https://doi.org/10.3390/sym13030488
Submission received: 20 December 2020 / Revised: 2 March 2021 / Accepted: 13 March 2021 / Published: 16 March 2021
(This article belongs to the Special Issue Chiral Symmetry in Physics)

Abstract

:
Coupling of ‘t Hooft’s determinant term is investigated in the framework of the three-flavor linear sigma model as a function of the chiral condensate. Using perturbation theory around the minimum point of the effective action, we calculate the renormalization group flow of the first field-dependent correction to the coupling of the conventional U A ( 1 ) breaking determinant term. It is found that, at low temperatures, mesonic fluctuations make the anomaly increase when the chiral condensate decreases. As an application, we analyze the effect at the zero temperature nuclear liquid–gas transition.

1. Introduction

The U A ( 1 ) subgroup of approximate U L ( 3 ) × U R ( 3 ) chiral symmetry is anomalously broken in quantum chromodynamics (QCD). QCD is a strongly coupled theory, and as such, most accurate results can be expected to emerge through lattice simulations. These, however, lack the ability to simulate the system at finite density due to the notorius sign problem. To tackle this issue, it is common to build effective models upon chiral symmetry, which are expected to capture essential features of QCD in the low energy regime. Even at zero density, these models are known to provide reasonable results for temperatures below that of the chiral transition [1].
In effective theories, such as the Nambu–Jona–Lasinio or linear sigma models, this is taken into account by ‘t Hooft’s determinant term. Coefficients of operators in the Lagrangian of field theories, including the aforementioned determinant term, are considered to be (coupling) constants, without any field or environment dependence. In the quantum version of the action, however, fluctuations introduce temperature (T), baryochemical potential ( μ B ) and also field dependence as they become coefficient functions. When talking about field dependence of a given coupling, one has in mind the resummation of higher dimensional operators that can reappear when Taylor expanding the coefficient functions in terms of the field variable(s) around a conveniently chosen expansion point.
In QCD, it is well established that the anomalous breaking of U A ( 1 ) symmetry should gradually disappear beyond the critical temperature, as at high T the instanton density causing the anomaly exponentially vanishes [2,3]. At lower temperatures, however, the situation is far from being understood in a satisfactory fashion. One has also great interest in gaining results regarding the anomaly evolution at finite μ B due to the sign problem, as mentioned earlier.
The finite temperature and/or density behavior of the U A ( 1 ) anomaly represents an active direction of research. More conservative results usually argue that the evaporation of the anomaly should follow that of the chiral condensate and thus the U A ( 1 ) symmetry restores around the critical temperature ( T C ) of the chiral transition [4,5,6,7,8]. There are also several arguments and results that indicate that it is visible even beyond T C [9,10,11,12,13]. For example, earlier renormalization group studies indicate that when considering a field-dependent anomaly coefficient, it decreases as a function of the chiral condensate [11,14], and this profile function can also depend explicitly on the tempereture (though the former effect is more dominant). Effective restoration of the anomaly has, e.g., a consequence regarding the order of the chiral transition [15], the axion mass [16], and the fate of η meson, whose mass if substantially drops [17,18,19,20,21,22], in a nuclear medium could lead to the formation of an η -nucleon bound state.
The goal of this paper is to calculate the first field-dependent correction to the coupling of the ‘t Hooft determinant in the effective action perturbatively, i.e., we determine the first Taylor coefficient of the anomaly function. Determining this coefficient (at least qualitatively) allows for obtaining the behavior of the anomaly strength as the function of the chiral condensate. This allows for providing answer for the question whether the aforementioned results on anomaly strengthening [11,14] can be reproduced within a simple perturbative renormalizatrion group setting of the linear sigma model, or determining a full functional dependence of the effective action is necessary. Fluctuations will be included using the functional variant of the renormalization group (FRG) [23], in the so-called local potential approximation (LPA). Even though renormalizable, we think of the model as an effective field theory; therefore, an ultra violet (UV) cutoff is inherently part of the system, which we set to Λ = 1 GeV (we expect the linear sigma model to emerge from QCD around this scale). Our task is to integrate out all fluctuations below Λ .
The paper is organized as follows. In Section 2, we introduce the model and the corresponding method of the FRG. Section 3 is devoted for calculating the effective action and discussing the problem of the expansion point of the Taylor series. After appropriate parametrization of the model, in Section 4, as an application, we show how the anomaly strengthens at the zero temperature nuclear liquid–gas phase transition. Section 5 contains the summary.

2. Model and Method

The model we are working with in this paper is the three-flavor linear sigma model, which is defined via the following Euclidean Lagrangian:
L = Tr ( i M i M ) + m 2 Tr ( M M ) + g 1 Tr ( M M ) 2 + g 2 Tr ( M M Tr ( M M ) ) 2 + a ( det M + det M ) ( h 0 s 0 + h 8 s 8 ) ,
where M contains the meson fields, M = ( s a + i π a ) T a [ T a = λ a / 2 are generators of the U ( 3 ) group with λ a being the Gell-Mann matrices, a = 0 , . . . , 8 ], m 2 is the mass parameter and g 1 , g 2 refer to independent quartic couplings. As discussed in the previous section, the determinant term and the corresponding a parameter is responsible for the U A ( 1 ) anomaly. We also have explicit symmetry breaking terms containing h 0 and h 8 , which represent finite quark masses.
Our main goal is to calculate the quantum effective action, Γ , built upon the theory defined via (1). As announced in the introduction, we think of (1) as an inherently effective model, which is only valid up to the scale Λ = 1 GeV , therefore, one needs to take into account fluctuations with a cutoff Λ . The scale dependent quantum effective action, Γ k , which includes fluctuations with momenta larger than k (i.e., they are integrated out) is defined as
Z k [ J ] = D M D M exp L ( J M + h . c . ) M R k M , Γ k [ M ] = log Z k [ J ] ( J M + h . c . ) M R k M ,
where we omitted matrix indices, J is the conjugate source variable to M, and R k is an appropriately chosen (bilocal) regulator function freezing fluctuations with momenta smaller than k. We note that the integrals can be considered either in direct or Fourier spaces. The Γ k functional, for homogeneous field configurations, obeys the following, so-called flow equation [23]:
k Γ k [ M ] = 1 2 Tr x q ( Γ k ( 2 ) + R k ) 1 ( q ) k R k ( q ) ,
where Γ k ( 2 ) is the second functional derivative of Γ k in a homogeneous background field M, thus the x integral merely gives a spacetime volume factor. We also assumed that the regulator is diagonal in momentum space.
Our aim is to calculate the scale dependent effective action, Γ k , in an approximation that takes into account the evolution of the anomaly at the next-to-leading order, i.e., we wish to determine in Γ k the coefficient of the operator Tr ( M M ) · ( det M + det M ) . Our ansatz for Γ k is as follows:
Γ k = x [ Tr ( i M i M ) + m k 2 Tr ( M M ) + g 1 , k Tr ( M M ) 2 + g 2 , k Tr ( M M Tr ( M M ) ) 2 ( h 0 s 0 + h 8 s 8 ) + a k ( det M + det M ) + a 1 , k Tr ( M M ) · ( det M + det M ) ] .
This is sometimes called the Local Potential Approximation (LPA), where momentum dependence is only introduced into the two point function, via the standard kinetic term in (4). Note that the LPA can be considered as the leading order of the derivative expansion, and there is substantial evidence that these kind of series do converge [23]. As seen in (4), instead of working with a completely general field-dependent potential, we are employing perturbation theory in terms of the small parameter 1 / Λ . That is, by gradually including higher dimensional operators, since their coefficients scale with inverse powers of the scale, the ansatz (4) in the UV can be thought of as a power series in 1 / Λ . We choose (4) to be compatible with (1), but all couplings come with k–dependence. The only exceptions are h 0 and h 8 , as one point couplings do not flow with respect to the scale. Furthermore, notice the new term proportional to a 1 , k , which is key for our purposes to determine the anomaly behavior at k = 0 . First, our task is to calculate Γ k ( 2 ) from (4), then plug it into (3), and identify the individual differential equations for m k 2 , g 1 , k , g 2 , k , a k , and a 1 , k . Finally, these equations need to be integrated from k = 1 GeV to k = 0 to obtain Γ Γ k = 0 . In the ansatz (4), obviously the actual strength of the anomaly is not described by the parameter a, but rather a + a 1 · Tr ( M M ) | min , where we need to evaluate the chiral condensates in the minimum point of the effective action. Therefore, what we are basically after is the relative sign of a 1 to a at k = 0 to decide whether the anomaly strengthens or weakens as the chiral condensate gradually evaporates.
We finally note that there are various choices for the regulator function, R k . In this paper we will stick to R k ( q , p ) R k ( q ) ( 2 π ) 3 δ ( q + p ) = ( k 2 q 2 ) Θ ( k 2 q 2 ) ( 2 π ) 4 δ ( q + p ) , where Θ ( x ) is the step function. This variant has been shown to be the optimal choice for the LPA [24], maximizing the radius of convergence of an amplitude expansion.

3. Calculation of the Effective Action

The first step is to calculate Γ k ( 2 ) . In principle it is a 18 × 18 matrix in the s a π a space, and there is not much hope that one can invert such a complicated expression analytically. Luckily, it is not necessary at all, as in (4) we kept the field dependence up to the order of O ( M 5 ) . By working with a restricted background, Γ k ( 2 ) is easily invertable and by expanding the rhs of (3) in terms of the field variables it still allows for identifying each operator that are being kept in (4).
A convenient choice is to work with M = s 0 T 0 + s 8 T 8 . In such a background, the operators that need to be identified are as follows:
ρ : = Tr ( M M ) = 1 2 ( s 0 2 + s 8 2 ) ,
τ : = Tr ( M M Tr ( M M ) ) 2 = 1 24 s 8 2 ( 8 s 0 2 4 2 s 0 s 8 + s 8 2 ) ,
Δ : = det M + det M = 1 36 ( 2 6 s 0 3 3 6 s 0 s 8 2 2 3 s 8 3 ) .
The Γ k ( 2 ) matrix elements in the scalar sector read
Γ k , s 0 s 0 ( 2 ) = q 2 + m k 2 + g 1 , k ( 3 s 0 2 + s 8 2 ) + 2 3 g 2 , k s 8 2 + 2 3 a k s 0 + a k , 1 5 3 2 3 s 0 3 1 2 6 s 0 s 8 2 1 6 3 s 8 3 ,
Γ k , s 0 s 8 ( 2 ) = 2 g 1 , k s 0 s 8 + g 2 , k 4 3 s 0 s 8 1 2 s 8 2 a k 6 s 8 a 1 , k 1 2 6 s 0 2 s 8 + 1 2 3 s 0 s 8 2 + 1 6 s 8 3 ,
Γ k , s 8 s 8 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + 3 s 8 2 ) + g 2 , k 2 3 s 0 2 2 s 0 s 8 + 1 2 s 8 2 a k 1 6 s 0 + 1 3 s 8 a 1 , k 1 6 6 s 0 3 + 1 2 3 s 0 2 s 8 + 3 2 s 0 s 8 2 + 5 3 3 s 8 3 ,
Γ k , s 1 s 1 ( 2 ) = Γ k , s 2 s 2 ( 2 ) = Γ k , s 3 s 3 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) + g 2 , k 2 3 s 0 2 + 2 s 0 s 8 + 1 6 s 8 2 + a k 1 6 s 0 + 1 3 s 8 + a 1 , k 1 6 6 s 0 3 + 1 2 3 s 0 2 s 8 1 6 s 0 s 8 2 + 1 3 3 s 8 3 ,
Γ k , s 4 s 4 ( 2 ) = Γ k , s 5 s 5 ( 2 ) = Γ k , s 6 s 6 ( 2 ) = Γ k , s 7 s 7 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) + g 2 , k 2 3 s 0 2 1 2 s 0 s 8 + 1 6 s 8 2 a k 1 6 s 0 + 1 2 3 s 8 a 1 , k 1 6 6 s 0 3 + 1 4 3 s 0 2 s 8 + 1 6 s 0 s 8 2 + 5 12 3 s 8 3 ,
while the pseudoscalar components are
Γ k , π 0 π 0 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) 2 3 a k s 0 a 1 , k 1 3 2 3 s 0 3 + 1 2 3 2 s 0 s 8 2 + 1 6 3 s 8 3 ,
Γ k , π 0 π 8 ( 2 ) = g 2 , k 2 3 s 0 s 8 1 3 2 s 8 2 + a k 6 s 8 + a 1 , k 1 2 6 s 0 2 s 8 + 1 2 6 s 8 3 ,
Γ k , π 8 π 8 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) + g 2 , k 2 3 s 0 s 8 + 1 6 s 8 2 + a k 1 6 s 0 + 1 3 s 8 + a 1 , k 5 6 1 6 s 0 3 + 1 2 3 s 0 2 s 8 + 1 3 3 s 8 3 ,
Γ k , π 1 π 1 ( 2 ) = Γ k , π 2 π 2 ( 2 ) = Γ k , π 3 π 3 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) + g 2 , k 2 3 s 0 s 8 1 6 s 8 2 + a k 1 6 s 0 1 3 s 8 + a 1 , k 5 6 1 6 s 0 3 1 2 3 s 0 2 s 8 2 3 3 s 8 3 ,
Γ k , π 4 π 4 ( 2 ) = Γ k , π 5 π 5 ( 2 ) = Γ k , π 6 π 6 ( 2 ) = Γ k , π 7 π 7 ( 2 ) = q 2 + m k 2 + g 1 , k ( s 0 2 + s 8 2 ) + g 2 , k 1 3 2 s 0 s 8 + 5 6 s 8 2 + a k 1 6 s 0 + 1 2 3 s 8 + a 1 , k 5 6 1 6 s 0 3 + 1 4 3 s 0 2 s 8 + 1 12 3 s 8 3 .
Using that k R k ( q ) = 2 k Θ ( k 2 q 2 ) and Γ k ( 2 ) ( q ) + R k ( q ) = Γ k ( 2 ) ( k ) for q < k , from (3) we get
k Γ k = x k 5 32 π 2 Tr Γ k ( 2 ) ( k ) 1 .
Plugging in the matrix elements calculated in the M = s 0 T 0 + s 8 T 8 background, we can expand the rhs of (16) in terms of s 0 and s 8 . After this step, using (5) we identify the ρ , τ and Δ operators as
k Γ k = x k m k 2 · ρ + k g 1 , k · ρ 2 + k g 2 , k · τ + k a k · Δ + k a 1 , k · ρ Δ + . . . ,
where
k m k 2 = k 5 32 π 2 8 ( 15 g 1 , k + 4 g 2 , k ) 3 ( k 2 + m k 2 ) 2 , k g 1 , k = k 5 32 π 2 8 ( 117 g 1 , k 2 + 48 g 1 , k g 2 , k + 16 g 2 , k 2 ) 9 ( k 2 + m k 2 ) 3 , k g 2 , k = k 5 32 π 2 48 g 1 , k g 2 , k + 32 g 2 , k 2 ( k 2 + m k 2 ) 3 , k a k = k 5 32 π 2 8 a k ( 3 g 1 , k 4 g 2 , k ) ( k 2 + m k 2 ) 3 24 a 1 , k ( k 2 + m k 2 ) 2 , k a 1 , k = k 5 32 π 2 32 a k ( 9 g 1 , k 2 + 6 g 1 , k g 2 , k + 2 g 2 , k 2 ) ( k 2 + m k 2 ) 4 + 16 a 1 , k ( 33 g 1 , k 2 g 2 , k ) 3 ( k 2 + m k 2 ) 3 .
Note that we treated the anomaly as perturbation and dropped every term beyond O ( a k , a 1 , k ) . Introducing scale independent variables, from (18) one easily reproduces the well known 1-loop β functions of the couplings in the linear sigma model [25]. Our task now is to solve Equation (18) starting from k = Λ 1 GeV to k = 0 .
Solving (18) all the way down to k = 0 would require m k 2 > 0 throughout the renormalization group flow. Since we wish to obtain phenomenologically reasonable results, the potential has to show spontaneous symmetry breaking. That is to say, when all fluctuations are integrated out, m k = 0 2 has to be negative. However, then there exists a critical scale k crit > 0 , for which all denominators in (18) blow up and the flow equations lose their meaning. The way out is to realize is that one actually has the choice to determine the flow equations in the minimum point of the effective action, Γ k | s 0 , min , s 8 , min , rather than evaluating it in a vanishing background. That is, all renormalization group flows are to be extracted at s 0 , min , s 8 , min . This way one always has a positive definite denominator and the flow equation is valid for any k.
One, therefore, repeats the calculations starting from (16), but this time expands only in terms of s 8 so that the ρ dependence of the parameters can be traced via s 0 . A long but straightforward calculation leads once again to the possibility of identifying the invariants appear in (17), whose coefficients now read as
k m k 2 = k 5 32 π 2 [ 18 g 1 , k ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 2 + 6 g 1 , k ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 2 + 16 ( 3 g 1 , k + 2 g 2 , k ) 3 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 2 + 72 g 1 , k 2 ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 3 + 72 g 1 , k 2 ρ 0 ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 3 + 64 ( 3 g 1 , k + 2 g 2 , k ) 2 ρ 0 9 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 3 ] , k g 1 , k = k 5 32 π 2 [ 36 g 1 , k 2 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 3 + 36 g 1 , k 2 ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 3 + 32 ( 3 g 1 , k + 2 g 2 , k ) 2 9 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 3 ] , k g 2 , k = k 5 32 π 2 [ 6 g 2 , k 2 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 3 9 g 2 , k / 2 ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 2 + 3 g 2 , k ( 6 g 1 , k + g 2 , k ) ρ 0 ( g 2 , k 3 g 1 , k ) ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 2 + 30 g 2 , k 2 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 3 + 3 g 2 , k ( g 2 , k 21 g 1 , k ) / 2 ρ 0 ( g 2 , k 3 g 1 , k ) ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 2 ] , k a k = k 5 32 π 2 [ 36 g 1 , k ( a k + 2 a 1 , k ρ 0 ) ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 3 18 ( a k + a 1 , k ρ 0 ) ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 2 12 g 1 , k ( 3 a k + 10 a 1 , k ρ 0 ) ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 3 6 a k + 10 a 1 , k ρ 0 ρ 0 ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 2 + 16 ( 3 g 1 , k + 2 g 2 , k ) ( 3 a k + a 1 , k ρ 0 ) 3 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 3 + 4 ( 6 a k + a 1 , k ρ 0 ) ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 2 ] ,
k a 1 , k = k 5 32 π 2 [ 36 g 1 , k ( a k + 2 a 1 , k ρ 0 ) ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 3 + 9 a k ρ 0 2 ( k 2 + m k 2 + 2 g 1 , k ρ 0 ) 2 + 12 g 1 , k ( 3 a k + 10 a 1 , k ρ 0 ) ρ 0 ( k 2 + m k 2 + 6 g 1 , k ρ 0 ) 3 16 ( 3 g 1 , k + 2 g 2 , k ) ( 3 a k + a 1 , k ρ 0 ) 3 ρ 0 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 3 12 a k ρ 0 2 ( k 2 + m k 2 + 2 g 1 , k ρ 0 + 4 g 2 , k ρ 0 / 3 ) 2 ] ,
where we have denoted the expansion point by ρ 0 , which is to be set to the value of ρ corresponding to the minimum point of the effective action [note that ρ = Tr ( M M ) / 2 ]. As a side remark, one easily checks that choosing ρ 0 = 0 (19) would lead back to the earlier results, (18). Our task is to integrate the system of equations (19) from k = Λ 1 GeV down to k = 0 with the boundary conditions m Λ 2 = m 2 , g 1 , Λ = g 1 , g 2 , Λ = g 2 , a Λ = a , a 1 , Λ = 0 , where m 2 , g 1 , g 2 , a are such constants that reproduce as accurately as possible the mesonic spectrum in the infrared. Here we used that at the UV scale the coefficient of the operator ρ Δ can be set to zero due to perturbative renormalizability. This might be questionable if the UV scale was not high enough, as being a dimension 5 operator, dimensional analyis suggests that its coefficient, a 1 , is of O ( 1 / Λ ) . Obviously, if the linear sigma model was not an effective theory, and Λ could be sent to infinity, the term in question would not be present. But, in principle the a 1 coupling could be included already in the UV action. Investigation of such a scenario is beyond the scope of this paper.
Before solving the coupled system of equations (19), we need to fix the explicit symmetry breaking terms, i.e., the values for h 0 , h 8 . Instead of h 0 and h 8 , we will work in the nonstrange–strange basis, i.e., h ns = 2 3 h 0 + 1 3 h 8 , h s = 1 3 h 0 2 3 h 8 . The partially conserved axialvector current (PCAC) relations give
m π 2 f π = h ns , m K 2 f K = h ns 2 + h s 2 ,
where m π 2 = δ 2 Γ / δ π i 2 ( q = 0 ) [ i = 1 , 2 , 3 ] and m K 2 = δ 2 Γ / δ π i 2 ( q = 0 ) [ j = 4 , 5 , 6 , 7 ] . Using physical pion and kaon masses, 140 MeV , 494 MeV , respectively, and decay constants, f π = 93 MeV , f K = 113 MeV , one gets
h ns = m π 2 f π ( 122 MeV ) 3 , h s = 1 2 ( 2 m K 2 f K m π 2 f π ) ( 335 MeV ) 3 ,
which is equivalent to
h 0 = 2 3 ( m π 2 f π / 2 + m K 2 f K ) ( 285 MeV ) 3 , h 8 = 2 3 ( m π 2 f π m K 2 f K ) ( 310 MeV ) 3 .
Now we use that Ward identities of chiral symmetry lead to
δ Γ δ s ns ( q = 0 ) = m π 2 s ns h ns , δ Γ δ s s ( q = 0 ) = m K 2 m π 2 2 s ns + m K 2 s s h s .
Combined with (21), this shows that no matter how we choose the remaining parameters m 2 , g 1 , g 2 , a , in the minimum point of the effective action
s ns , min = f π , s s , min = 2 ( f K f π / 2 ) .
That leads to ρ 0 = ( s ns , min 2 + s s , min 2 ) / 2 , and thus we are ready to fix the aforementioned parameters. Solving (19), the values { m 2 , g 1 , g 2 , a } { 0.6835 GeV 2 , 29.7 , 91.5 , 4.4 GeV } lead to the masses of the pion, kaon, η , η as m π 133 MeV , m K 494 MeV , m η 537 MeV , m η 957 MeV , respectively. Note that the strength of the axial anomaly, a, is negative, and it remains so throughout the renormalization group flow. However, after solving (19), one concludes that at k = 0 the coefficient a 1 , k = 0 is positive. That is to say, since the actual strength of the determinant term is A : = a k = 0 + a 1 , k = 0 · ρ 0 , when the chiral condensate evaporates, the absolute value of A becomes larger. This shows that at low temperatures T, where the T dependence of the anomaly parameters is negligible, the anomaly is actually strengthening as chiral symmetry gradually restores. That is one of the main results of the paper.
In what follows, we provide a rough estimate how the anomaly behaves at the zero temperature nuclear liquid–gas transition.

4. Anomaly Strengthening at The Nuclear Liquid–Gas Transition

In this section we apply our results to the zero temperature nuclear liquid–gas transition. We assume that the nucleon field couples to the mesons via Yukawa interaction, L int = g ψ ¯ M 5 ψ , ψ T = ( p , n ) , M 5 = a = ns , 1 , 2 , 3 ( s a + i π a γ 5 ) T a , where the nonstrange generator is T ns = 2 / 3 T 0 + 1 / 3 T 8 , while γ 5 is the fifth Dirac matrix. In principle, one would also need to include the dynamics of an ω vector particle into the system [26,27] that models the repulsive interaction between nucleons, but as we will see in a moment, for our purposes it plays no role.
First, we exploit some of the zero temperature properties of nuclear matter. Note that, in the current model, the nucleon mass entirely originates from the spontaneous breaking of chiral symmetry,
m N ( s ns ) = g Y s ns / 2 ,
and since m N ( f π ) 939 MeV in the vacuum, we arrive at g Y 20.19 . Normal nuclear density, n N 0.17 fm 3 ( 109.131 MeV ) 3 leads to the Fermi momentum, p F , of the nucleons, since at the mean field level, for T = 0 we have
n N = 4 p n F p 2 + m N 2 p F | T = 0 2 3 π 2 p F 3 ,
therefore, p F 267.9 MeV 1.36 fm 1 . This leads to the nonstrange condensate in the liquid phase, s ns , liq , because the Landau mass, which is defined as
M L = p F 2 + m N 2 ( s ns , liq )
is known to be M L 0.8 m N ( f π ) 751.2 MeV , and thus s ns , liq 69.52 MeV [26,27]. This shows that as we increase the chemical potential, the nonstrange chiral condensate, s ns , jumps: f π s ns , liq . This will definitely be accompanied by a jump in the strange condensate, but it has been shown to be significantly smaller [14]. Neglecting the change in s s , the ρ chiral invariant jumps as ( f π 2 + s s , min 2 ) / 2 ( s ns , liq 2 + s s , min 2 ) / 2 . As discussed in the previous section, the anomaly strength is A = a k = 0 + a 1 , k = 0 · ρ , which also jumps accordingly, and the change in A becomes
Δ A = a 1 , k = 0 · Δ ρ ,
where Δ ρ = ( s ns , liq 2 f π 2 ) / 2 . Solving (19) one gets a k = 0 9.05 GeV and a 1 , k = 0 494.5 GeV 1 , therefore, the relative change in the anomaly at the liquid–gas transition is
Δ A A = a 1 , k = 0 · Δ ρ a k = 0 + a 1 , k = 0 · ρ 0 0.2 = 20 % ,
which is in the ballpark of the result of [14]. One can now check how robust this result is with respect to changing the cutoff Λ . A thorough investigation reveals that in a cutoff interval of 0.8 1.5 GeV , while the mass spectrum can be maintained within a few percent error after reparametrization, the Δ A / A ratio is less stable. One finds that the latter is a monotonically decreasing function of the cutoff and varies roughly between 15– 40 % in the above interval. Results show that when going below 1 GeV the cutoff dependence gets stronger, which is understandable, since non-renormalizable operators are absent at the UV scale. That is, if the latter is chosen to be too small, the model cannot provide robust results (more operators would be needed). Going beyond 1.5 2 GeV , in turn, would be physically inappropriate as at those scales quark degrees of freedom would definitely play a crucial role. From these findings it is safe to say is that the relative change of the anomaly strength is of O ( 10 % ) at the transition point.
At this point, we once again wish to emphasize that we have neglected the drop in the strange condensate, and also, the present analysis is based on perturbation theory. In principle higher order operators that break the U A ( 1 ) subgroup should also be resummed, e.g., terms such as Tr ( M M ) n ( det M + det M ) could be of huge importance. The lesson we wish to point out here is that the present, rather simple perturbative calculation can also capture the phenomenon of strengthening anomaly as the chiral condensate (partially) evaporates.

5. Conclusions

In this paper, we investigated how the U A ( 1 ) anomaly behaves as a function of the chiral condensate. We worked with the three-flavor linear sigma model, and calculated the leading correction in a 1 / Λ expansion to the conventional anomaly term caused by quantum fluctuations. We have found that the coefficient of the aforementioned operator, Tr ( M M ) · ( det M + det M ) , causes the actual strength of the anomaly to become larger once the chiral condensate evaporates. For the sake of an example, we demonstrated that at the zero temperature nuclear liquid–gas transition, where (on top of a jump in the nuclear density) the chiral condensate partially restores, the actual strength of the anomaly increases. This could also happen toward the full restoration of chiral symmetry, where quark dynamics also play a significant role. Note that our findings are based solely on calculating mesonic fluctuations, and no instanton effects have been taken into account.
The linear sigma model, being an effective field theory, cannot accommodate instantons as the fundamental model of QCD. Still, there are at least two directions worth exploring in the effective model framework. Recently it has been shown [28] that 3 Q -point interactions are generated by instantons with Q topological charge, which can be embedded into the linear sigma model via [ ( det M ) Q + ( det M ) Q ] operators. Another important direction could be to assign environment dependence even to the bare anomaly coefficient(s) from QCD data and see how these compete against thermal effects caused by mesonic fluctuations.
Finally, we wish to point out that our study calls for an extension via a non-perturbative treatment, where fluctuations are taken into account beyond the O ( a ) order, and the coefficient function of the determinant term is obtained in a functional fashion, rather than at the lowest order of its Taylor series. The aforementioned directions are under progress and will be reported in a separate study.

Funding

This work was supported by the National Research, Development and Innovation Fund under Project No. PD127982, by the János Bolyai Research Scholarship of the Hungarian Academy of Sciences, and by the ÚNKP-20-5 New National Excellence Program of the Ministry for Innovation and Technology from the source of the National Research, Development and Innovation Fund.

Acknowledgments

The author thanks A. Patkos and Zs. Szep for discussions on related topics.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Fukushima, K.; Hatsuda, T. The phase diagram of dense QCD. Rept. Prog. Phys. 2011, 74, 014001. [Google Scholar] [CrossRef] [Green Version]
  2. Hooft, G. Computation of the quantum effects due to a four-dimensional pseudoparticle. Phys. Rev. D 1976, 14, 3432. [Google Scholar]
  3. Gross, D.J.; Pisarski, R.D.; Yaffe, L.G. QCD and instantons at finite temperature. Rev. Mod. Phys. 1981, 53. [Google Scholar] [CrossRef]
  4. Mitter, M.; Schaefer, B.-J. Fluctuations and the axial anomaly with three quark flavors. Phys. Rev. D 2014, 89, 054027. [Google Scholar] [CrossRef] [Green Version]
  5. Cossu, G.; Fukaya, H.; Hashimoto, S.; Noaki, J.; Tomiya, A. On the Axial U(1) Symmetry at Finite Temperature. Available online: http://arxiv.org/abs/1511.05691 (accessed on 12 March 2021).
  6. Ishii, M.; Yonemura, K.; Takahashi, J.; Kouno, J.; Yahiro, M. Determination of U(1)A restoration from pion and a0-meson screening masses: Toward the chiral regime. Phys. Rev. D 2016, 93, 016002. [Google Scholar] [CrossRef] [Green Version]
  7. Nicola, A.G.; de Elvira, J.R. Chiral and U(1)A restoration for the scalar and pseudoscalar meson nonets. Phys. Rev. D 2018, 98, 014020. [Google Scholar] [CrossRef] [Green Version]
  8. Bottaro, S.; Meggiolaro, E. QCD axion and topological susceptibility in chiral effective Lagrangian models at finite temperature. Phys. Rev. D 2020, 102, 014048. [Google Scholar] [CrossRef]
  9. Dick, V.; Karsch, F.; Laermann, E.; Mukherjee, S.; Sharma, S. Microscopic origin of UA(1) symmetry violation in the high temperature phase of QCD. Phys. Rev. D 2015, 91, 094504. [Google Scholar] [CrossRef] [Green Version]
  10. Sharma, S.; Dick, V.; Karsch, F.; Laermann, E.; Mukherjee, S. The topological structures in strongly coupled QGP with chiral fermions on the lattice. Nucl. Phys. A 2016, 956, 793. [Google Scholar] [CrossRef] [Green Version]
  11. Fejos, G.; Hosaka, A. Thermal properties and evolution of the UA(1) factor for 2+1 flavors. Phys. Rev. D 2016, 94, 036005. [Google Scholar] [CrossRef] [Green Version]
  12. Li, X.; Fu, W.; Liu, Y. New insight about the effective restoration of UA(1) symmetry. Phys. Rev. D 2020, 101, 054034. [Google Scholar] [CrossRef] [Green Version]
  13. Braun, J.; Leonhardt, M.; Pawlowski, J.M.; Rosenblüh, D. Chiral and Effective U(1)A Symmetry Restoration in QCD. Available online: http://arxiv.org/abs/2012.06231 (accessed on 12 March 2021).
  14. Fejos, G.; Hosaka, A. Axial anomaly and hadronic properties in a nuclear medium. Phys. Rev. D 2018, 98, 036009. [Google Scholar] [CrossRef] [Green Version]
  15. Pelissetto, A.; Vicari, E. Relevance of the axial anomaly at the finite-temperature chiral transition in QCD. Phys. Rev. D 2013, 88, 105018. [Google Scholar] [CrossRef] [Green Version]
  16. Horvatic, D.; Kekez, D.; Klabucar, D. Temperature dependence of the axion mass in a scenario where the restoration of chiral symmetry drives the restoration of the UA(1) Symmetry. Universe 2019, 10, 208. [Google Scholar] [CrossRef] [Green Version]
  17. Nagahiro, H.; Takizawa, M.; Hirenzaki, S. eta- and eta’-mesic nuclei and UA(1) anomaly at finite density. Phys. Rev. C 2006, 74, 045203. [Google Scholar] [CrossRef] [Green Version]
  18. Sakai, S.; Jido, D. In medium eta’ mass and eta’ N interaction based on chiral effective theory. Phys. Rev. C 2013, 88, 064906. [Google Scholar] [CrossRef] [Green Version]
  19. Sakai, S.; Hosaka, A.; Nagahiro, H. Effect of the final state interaction of eta’N on the eta’ photoproduction off nucleon. Phys. Rev. C 2017, 95, 045206. [Google Scholar] [CrossRef] [Green Version]
  20. Jido, D.; Masutani, H.; Hirenzaki, S. Structure of eta’ mesonic nuclei in a relativistic mean field theory. Prog. Theor. Exp. Phys. 2019, 2019, 053D02. [Google Scholar] [CrossRef]
  21. Horvatic, D.; Kekez, D.; Klabucar, D. Eta’ and eta mesons at high T when the UA(1) and chiral symmetry breaking are tied. Phys. Rev. D 2019, 99, 014007. [Google Scholar] [CrossRef] [Green Version]
  22. Rai, S.K.; Tiwari, V.K. Exploring axial U(1) restoration in a modified 2+1 flavor Polyakov quark meson model. Eur. Phys. J. Plus 2020, 135, 844. [Google Scholar] [CrossRef]
  23. Kopietz, P.; Bartosch, L.; Schütz, F. Introduction to the Functional Renormalization Group; Springer: New York, NY, USA, 2010. [Google Scholar]
  24. Litim, D. Optimisation of the exact renormalisation group. Phys. Lett. B 2000, 486, 92–99. [Google Scholar] [CrossRef] [Green Version]
  25. Pisarski, R.D.; Wilczek, F. Remarks on the chiral phase transition in chromodynamics. Phys. Rev. D 1984, 29, 338. [Google Scholar] [CrossRef]
  26. Floerchinger, S.; Wetterich, C. Chemical freeze-out in heavy ion collisions at large baryon densities. Nucl. Phys. A 2012, 890, 891. [Google Scholar] [CrossRef] [Green Version]
  27. Drews, M.; Weise, W. Functional renormalization group studies of nuclear and neutron matter. Prog. Part. Nucl. Phys. 2017, 93, 69. [Google Scholar] [CrossRef] [Green Version]
  28. Pisarski, R.D.; Rennecke, F. Multi-instanton contributions to anomalous quark interactions. Phys. Rev. D 2020, 101, 114019. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fejős, G. Perturbative RG Analysis of the Condensate Dependence of the Axial Anomaly in the Three-Flavor Linear Sigma Model. Symmetry 2021, 13, 488. https://doi.org/10.3390/sym13030488

AMA Style

Fejős G. Perturbative RG Analysis of the Condensate Dependence of the Axial Anomaly in the Three-Flavor Linear Sigma Model. Symmetry. 2021; 13(3):488. https://doi.org/10.3390/sym13030488

Chicago/Turabian Style

Fejős, Gergely. 2021. "Perturbative RG Analysis of the Condensate Dependence of the Axial Anomaly in the Three-Flavor Linear Sigma Model" Symmetry 13, no. 3: 488. https://doi.org/10.3390/sym13030488

APA Style

Fejős, G. (2021). Perturbative RG Analysis of the Condensate Dependence of the Axial Anomaly in the Three-Flavor Linear Sigma Model. Symmetry, 13(3), 488. https://doi.org/10.3390/sym13030488

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop