Chiral Limit Research Articles

Holst-MacDowell-Mansouri action for (extended) supergravity with boundaries and super Chern-Simons theory

In this article, the Cartan geometric approach toward (extended) supergravity in the presence of boundaries will be discussed. In particular, based on new developments in this field, we will derive the Holst variant of the MacDowell-Mansouri action for mathcal{N} = 1 and mathcal{N} = 2 pure AdS supergravity in D = 4 for arbitrary Barbero-Immirzi parameters. This action turns out to play a crucial role in context of boundaries in the framework of supergravity if one imposes supersymmetry invariance at the boundary. For the mathcal{N} = 2 case, it follows that this amounts to the introduction of a θ-topological term to the Yang-Mills sector which explicitly depends on the Barbero-Immirzi parameter. This shows the close connection between this parameter and the θ-ambiguity of gauge theory.We will also discuss the chiral limit of the theory, which turns out to possess some very special properties such as the manifest invariance of the resulting action under an enlarged gauge symmetry. Moreover, we will show that demanding supersymmetry invariance at the boundary yields a unique boundary term corresponding to a super Chern-Simons theory with OSp( mathcal{N} |2) gauge group. In this context, we will also derive boundary conditions that couple boundary and bulk degrees of freedom and show equivalence to the results found in the D’Auria-Fré approach in context of the non-chiral theory. These results provide a step towards of quantum description of supersymmetric black holes in the framework of loop quantum gravity.

The role of the UA(1) breaking term in dynamical chiral symmetry breaking of chiral effective theories

Abstract The spontaneous breaking of chiral symmetry is examined by chiral effective theories, such as the linear $\sigma$ model and the Nambu–Jona-Lasinio (NJL) model. We discuss the properties of the sigma meson regarded as the quantum fluctuation of the chiral condensate when the chiral symmetry is spontaneously broken, mainly by the U$_{A}$(1) anomaly. We derive a mass relation among the SU(3) flavor singlet mesons, $\eta_{0}$ and $\sigma_{0}$, in the linear $\sigma$ model to be satisfied for the anomaly-driven symmetry breaking in the chiral limit, and find that it is also supported in the NJL model. With the explicit breaking of chiral symmetry, we show that both of the chiral effective models reproducing the observed physical quantities suggest that the $\sigma_{0}$ meson should have a mass smaller than $\sim$800 MeV when anomaly-driven symmetry breaking takes place.

Contribution of the QCD Θ -term to the nucleon electric dipole moment

We present a calculation of the contribution of the $\Theta$-term to the neutron and proton electric dipole moments using seven 2+1+1-flavor HISQ ensembles. We also estimate the topological susceptibility for the 2+1+1 theory to be $\chi_Q = (66(9)(4) \rm MeV)^4$ in the continuum limit at $M_\pi = 135$ MeV. The calculation of the nucleon three-point function is done using Wilson-clover valence quarks. The CP-violating form factor $F_3$ is calculated by expanding in small $\Theta$. We show that lattice artifacts introduce a term proportional to $a$ that does not vanish in the chiral limit, and we include this in our chiral-continuum fits. A chiral perturbation theory analysis shows that the $N(0) \pi(0)$ state should provide the leading excited state contribution, and we study the effect of such a state. Detailed analysis of the contributions to the neutron and proton electric dipole moment using two strategies for removing excited state contamination are presented. Using the excited state spectrum from fits to the two-point function, we find $d_n^\Theta$ is small, $|d_n^\Theta| \lesssim 0.01 \overline \Theta e$ fm, whereas for the proton we get $|d_p^\Theta| \sim 0.02 \overline \Theta e$ fm. On the other hand, if the dominant excited-state contribution is from the $N \pi$ state, then $|d_n^\Theta|$ could be as large as $0.05 \overline \Theta e$ fm and $|d_p^\Theta| \sim 0.07 \overline \Theta e$ fm. Our overall conclusion is that present lattice QCD calculations do not provide a reliable estimate of the contribution of the $\Theta$-term to the nucleon electric dipole moments, and a factor of ten higher statistics data are needed to get better control over the systematics and possibly a $3\sigma$ result.

Momentum Space Quantum Monte Carlo on Twisted Bilayer GrapheneSupported by the RGC of Hong Kong SAR of China (Grant Nos. 17303019, 17301420, and AoE/P-701/20), the National Key Research and Development Program of China (Grant No. 2016YFA0300502), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB33000000 and XDB28000000), the National Natural Science Foundation of China (Grant

We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG). The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation when either the two valleys or the two spin degrees of freedom are considered. By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point (CNP). Then, with the exact Green’s functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings.

Fully coupled functional equations for the quark sector of QCD

We present a comprehensive study of the quark sector of $2+1$ flavour QCD, based on a self-consistent treatment of the coupled system of Schwinger-Dyson equations for the quark propagator and the full quark-gluon vertex. The individual form factors of the quark-gluon vertex are expressed in a special tensor basis obtained from a set of gauge-invariant operators. The sole external ingredient used as input to our equations is the Landau gauge gluon propagator with $2+1$ dynamical quark flavours, obtained from studies with Schwinger-Dyson equations, the functional renormalisation group approach, and large volume lattice simulations. The appropriate renormalisation procedure required in order to self-consistently accommodate external inputs stemming from other functional approaches or the lattice is discussed in detail, and the value of the gauge coupling is accurately determined at two vastly separated renormalisation group scales. Our analysis establishes a clear hierarchy among the vertex form factors. We identify only three dominant ones, in agreement with previous results. The components of the quark propagator obtained from our approach are in excellent agreement with the results from Schwinger-Dyson equations, the functional renormalisation group, and lattice QCD simulation, a simple benchmark observable being the chiral condensate in the chiral limit, which is computed as $(245\,\textrm{MeV})^3$. The present approach has a wide range of applications, including the self-consistent computation of bound-state properties and finite temperature and density physics, which are briefly discussed.

Twisted bilayer graphene. III. Interacting Hamiltonian and exact symmetries

We derive the explicit Hamiltonian of twisted bilayer graphene (TBG) with Coulomb interaction projected into the flat bands, and study the symmetries of the Hamiltonian. First, we show that all projected TBG Hamiltonians can be written as Positive Semidefinite Hamiltonian, the first example of which was found in [PRL 122, 246401]. We then prove that the interacting TBG Hamiltonian exhibits an exact U(4) symmetry in the exactly flat band (nonchiral-flat) limit. We further define, besides a first chiral limit where the AA stacking hopping is zero, a new second chiral limit where the AB/BA stacking hopping is zero. In the first chiral-flat limit (or second chiral-flat limit) with exactly flat bands, the TBG is enhanced to have an exact U(4)$\times$U(4) symmetry, whose generators are different between the two chiral limits. While in the first chiral limit and in the non-chiral case these symmetries have been found in [PRX 10, 031034] for the $8$ lowest bands, we here prove that they are valid for projection into any $8 n_\text{max}$ particle-hole symmetric TBG bands, with $n_\text{max}>1$ being the practical case for small twist angles $<1^\circ$. Furthermore, in the first or second chiral-nonflat limit without flat bands, an exact U(4) symmetry still remains. We also elucidate the link between the U(4) symmetry presented here and the similar but different U(4) of [PRL 122, 246401]. Furthermore, we show that our projected Hamiltonian can be viewed as the normal-ordered Coulomb interaction plus a Hartree-Fock term from passive bands, and exhibits a many-body particle-hole symmetry which renders the physics symmetric around charge neutrality. We also provide an efficient parameterization of the interacting Hamiltonian. The existence of two chiral limits, with an enlarged symmetry suggests a possible duality of the model yet undiscovered.

Twisted bilayer graphene. VI. An exact diagonalization study at nonzero integer filling

Using exact diagonalization, we study the projected Hamiltonian with Coulomb interaction in the 8 flat bands of first magic angle twisted bilayer graphene. Employing the U(4) (U(4)$\times$U(4)) symmetries in the nonchiral (chiral) flat band limit, we reduced the Hilbert space to an extent which allows for study around $\nu=\pm 3,\pm2,\pm1$ fillings. In the first chiral limit $w_0/w_1=0$ where $w_0$ ($w_1$) is the $AA$ ($AB$) stacking hopping, we find that the ground-states at these fillings are extremely well-described by Slater determinants in a so-called Chern basis, and the exactly solvable charge $\pm1$ excitations found in [arXiv:2009.14200] are the lowest charge excitations up to system sizes $8\times8$ (for restricted Hilbert space) in the chiral-flat limit. We also find that the Flat Metric Condition (FMC) used in [arXiv:2009.11301,2009.11872,2009.12376,2009.13530,2009.14200] for obtaining a series of exact ground-states and excitations holds in a large parameter space. For $\nu=-3$, the ground state is the spin and valley polarized Chern insulator with $\nu_C=\pm1$ at $w_0/w_1\lesssim0.9$ (0.3) with (without) FMC. At $\nu=-2$, we can only numerically access the valley polarized sector, and we find a spin ferromagnetic phase when $w_0/w_1\gtrsim0.5t$ where $t\in[0,1]$ is the factor of rescaling of the actual TBG bandwidth, and a spin singlet phase otherwise, confirming the perturbative calculation [arXiv:2009.13530]. The analytic FMC ground state is, however, predicted in the intervalley coherent sector which we cannot access [arXiv:2009.13530]. For $\nu=-3$ with/without FMC, when $w_0/w_1$ is large, the finite-size gap $\Delta$ to the neutral excitations vanishes, leading to phase transitions. Further analysis of the ground state momentum sectors at $\nu=-3$ suggests a competition among (nematic) metal, momentum $M_M$ ($\pi$) stripe and $K_M$-CDW orders at large $w_0/w_1$.

Twisted bilayer graphene. II. Stable symmetry anomaly

We show that the entire continuous model of twisted bilayer graphene (TBG) (and not just the two active bands) with particle-hole symmetry is anomalous and hence incompatible with a lattice model. Previous works, e.g., [Phys. Rev. Lett. 123, 036401], [Phys. Rev. X 9, 021013], [Phys. Rev. B 99, 195455], and others [1-4] found that the two flat bands in TBG possess a fragile topology protected by the $C_{2z}T$ symmetry. [Phys. Rev. Lett. 123, 036401] also pointed out an approximate particle-hole symmetry ($\mathcal{P}$) in the continuous model of TBG. In this work, we numerically confirm that $\mathcal{P}$ is indeed a good approximation for TBG and show that the fragile topology of the two flat bands is enhanced to a $\mathcal{P}$-protected stable topology. This stable topology implies $4l+2$ ($l\in\mathbb{N}$) Dirac points between the middle two bands. The $\mathcal{P}$-protected stable topology is robust against arbitrary gap closings between the middle two bands the other bands. We further show that, remarkably, this $\mathcal{P}$-protected stable topology, as well as the corresponding $4l + 2$ Dirac points, cannot be realized in lattice models that preserve both $C_{2z}T$ and $\mathcal{P}$ symmetries. In other words, the continuous model of TBG is anomalous and cannot be realized on lattices. Two other topology related topics, with consequences for the interacting TBG problem, i.e., the choice of Chern band basis in the two flat bands and the perfect metal phase of TBG in the so-called second chiral limit, are also discussed.

Twisted bilayer graphene. IV. Exact insulator ground states and phase diagram

We derive the exact insulator ground states of the projected Hamiltonian of magic-angle twisted bilayer graphene (TBG) flat bands with Coulomb interactions in various limits, and study the perturbations away from these limits. We define the (first) chiral limit where the AA stacking hopping is zero, and a flat limit with exactly flat bands. In the chiral-flat limit, the TBG Hamiltonian has a U(4)$\times$U(4) symmetry, and we find that the exact ground states at integer filling $-4\le \nu\le 4$ relative to charge neutrality are Chern insulators of Chern numbers $\nu_C=4-|\nu|,2-|\nu|,\cdots,|\nu|-4$, all of which are degenerate. This confirms recent experiments where Chern insulators are found to be competitive low-energy states of TBG. When the chiral-flat limit is reduced to the nonchiral-flat limit which has a U(4) symmetry, we find $\nu=0,\pm2$ has exact ground states of Chern number $0$, while $\nu=\pm1,\pm3$ has perturbative ground states of Chern number $\nu_C=\pm1$, which are U(4) ferromagnetic. In the chiral-nonflat limit with a different U(4) symmetry, different Chern number states are degenerate up to second order perturbations. In the realistic nonchiral-nonflat case, we find that the perturbative insulator states with Chern number $\nu_C=0$ ($0<|\nu_C|<4-|\nu|$) at integer fillings $\nu$ are fully (partially) intervalley coherent, while the insulator states with Chern number $|\nu_C|=4-|\nu|$ are valley polarized. However, for $0<|\nu_C|\le4-|\nu|$, the fully intervalley coherent states are highly competitive (0.005meV/electron higher). At nonzero magnetic field $|B|>0$, a first-order phase transition for $\nu=\pm1,\pm2$ from Chern number $\nu_C=\text{sgn}(\nu B)(2-|\nu|)$ to $\nu_C=\text{sgn}(\nu B)(4-|\nu|)$ is expected, which agrees with recent experimental observations. Lastly, the TBG Hamiltonian reduces into an extended Hubbard model in the stabilizer code limit.

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k·p two-band model

We investigate the Twisted Bilayer Graphene (TBG) model to obtain an analytic understanding of its energetics and wavefunctions needed for many-body calculations. We provide an approximation scheme which first elucidates why the BM $K_M$-point centered calculation containing only $4$ plane-waves provides a good analytical value for the first magic angle. The approximation scheme also elucidates why most many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a $\Gamma_M$-point centered model of 6 plane-waves, we analytically understand the small $\Gamma_M$-point gap between the active and passive bands in the isotropic limit $w_0=w_1$. Furthermore, we analytically calculate the group velocities of passive bands in the isotropic limit, and show that they are \emph{almost} doubly degenerate, while no symmetry forces them to be. Furthermore, away from $\Gamma_M$ and $K_M$ points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite it is defined only by vanishing $K_M$-point Dirac velocity. We derive analytically a connected "magic manifold" $w_1=2\sqrt{1+w_0^2}-\sqrt{2+3w_0^2}$, on which the bands remain extremely flat as $w_0$ is tuned between the isotropic ($w_0=w_1$) and chiral ($w_0=0$) limits. We analytically show why going away from the isotropic limit by making $w_0$ less (but not larger) than $w_1$ increases the $\Gamma_M$- point gap between active and passive bands. Finally, perturbatively, we provide an analytic $\Gamma_M$ point $k\cdot p$ $2$-band model that reproduces the TBG band structure and eigenstates in a certain $w_0,w_1$ parameter range. Further refinement of this model suggests a possible faithful $2$-band $\Gamma_M$ point $k\cdot p$ model in the full $w_0, w_1$ parameter range.

Twisted symmetric trilayer graphene: Single-particle and many-body Hamiltonians and hidden nonlocal symmetries of trilayer moiré systems with and without displacement field

We derive the Hamiltonian for trilayer moir\'e systems with the Coulomb interaction projected onto the bands near the charge neutrality point. Motivated by the latest experimental results, we focus on the twisted symmetric trilayer graphene (TSTG) with a mirror-symmetry with respect to the middle layer. We provide a full symmetry analysis of the non-interacting Hamiltonian with a perpendicular displacement field coupling the band structure made otherwise of the twisted bilayer graphene (TBG) and the high velocity Dirac fermions, and we identify a hidden non-local symmetry of the problem. In the presence of this displacement field, we construct an approximate single-particle model, akin to the tripod model for TBG, capturing the essence of non-interacting TSTG. We also derive more quantitative perturbation schemes for the low-energy physics of TSTG with displacement field, obtaining the corresponding eigenstates. This allows us to obtain the Coulomb interaction Hamiltonian projected in the active band TSTG wavefunctions and derive the full many-body Hamiltonian of the system. We also provide an efficient parameterization of the interacting Hamiltonian. Finally, we show that the discrete symmetries at the single-particle level promote the $\mathrm{U} \left( 2 \right) \times \mathrm{U} \left( 2 \right)$ spin-valley symmetry to enlarged symmetry groups of the interacting problem under different limits. The interacting part of the Hamiltonian exhibits a large $\mathrm{U} \left( 4 \right) \times \mathrm{U} \left( 4 \right) \times \mathrm{U} \left( 4 \right) \times \mathrm{U} \left( 4 \right)$ symmetry in the chiral limit. Moreover, by identifying a new symmetry which we dub spatial many-body charge conjugation, we show that the physics of TSTG is symmetric around charge neutrality.

Study of the axial U(1) anomaly at high temperature with lattice chiral fermions

We investigate the axial $U(1)$ anomaly of two-flavor QCD at temperatures 190--330 MeV. In order to preserve precise chiral symmetry on the lattice, we employ the M\"obius domain-wall fermion action as well as overlap fermion action implemented with a stochastic reweighting technique. Compared to our previous studies, we reduce the lattice spacing to 0.07 fm, simulate larger multiple volumes to estimate finite size effect, and take more than four quark mass points, including one below physical point to investigate the chiral limit. We measure the topological susceptibility, axial $U(1)$ susceptibility, and examine the degeneracy of $U(1)$ partners in meson/baryon correlators. All the data above the critical temperature indicate that the axial $U(1)$ violation is consistent with zero within statistical errors. The quark mass dependence suggests disappearance of the $U(1)$ anomaly at a rate comparable to that of the $SU(2{)}_{L}\ifmmode\times\else\texttimes\fi{}SU(2{)}_{R}$ symmetry breaking.

Analysis of nonperturbative flavor violation at chiral crossover criticality in QCD

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.

Clear correlation between monopoles and the chiral condensate in SU(3) QCD

We study spontaneous chiral-symmetry breaking in SU(3) QCD in terms of the dual superconductor picture for quark confinement in the maximally Abelian (MA) gauge, using lattice QCD Monte Carlo simulations with four different lattices of $16^4$, $24^4$, $24^3\times 6$ at $\beta=6.0$ (i.e., the spacing $a \simeq$ 0.1 fm), and $32^4$ at $\beta=6.2$ (i.e., $a \simeq$ 0.075 fm), at the quenched level. First, in the confinement phase, we find Abelian dominance and monopole dominance in the MA gauge for the chiral condensate in the chiral limit,using the two different methods of i) the Banks-Casher relation with the Dirac eigenvalue density and ii) finite quark-mass calculations with the quark propagator and its chiral extrapolation. In the high-temperature deconfined phase, the chiral restoration is observed also for the Abelian and the monopole sectors. Second, we investigate local correlation between the chiral condensate and monopoles, which topologically appear in the MA gauge. We find that the chiral condensate locally takes a quite large value near monopoles. As an interesting possibility, the strong magnetic field around monopoles is responsible to chiral symmetry breaking in QCD, similarly to the magnetic catalysis.

Competition of inhomogeneous chiral phases and two-flavor color superconductivity in the NJL model

We study the phase structure of the two-flavor Nambu--Jona-Lasinio (NJL) model in the chiral limit, extending a previous study of the competition of an inhomogeneous chiral phase and a two-flavor color-superconducting (2SC) phase [1, 2]. There, an analytic expression for the dispersion relations for quasiparticle excitations in the presence of both a particular inhomogeneous chiral condensate, the so-called chiral density wave (CDW), and a homogeneous 2SC condensate was found. In this work we show how to determine the dispersion relations for arbitrary modulations of the chiral condensate in the presence of a homogeneous 2SC condensate, if the dispersion relations in the absence of color superconductivity are known. In our calculations, we employ two different Ans\"atze for the inhomogeneous chiral condensate, the CDW as well as the real-kink crystal (RKC). Depending on the value of the diquark coupling we find a region of the phase diagram where the inhomogeneous chiral and the 2SC condensates coexist, confirming results of Refs. [1, 2]. Decreasing the diquark coupling favors the inhomogeneous phase over the coexistence phase. On the other hand, increasing the diquark coupling leads to a larger 2SC phase, while the inhomogeneous chiral and the coexistence phases become smaller. In agreement with previous studies the RKC Ansatz is energetically preferred over the CDW Ansatz. Both Ans\"atze lead to a qualitatively similar phase diagram, however the coexistence phase is smaller for the RKC Ansatz.

Correlated Dirac Eigenvalues and Axial Anomaly in Chiral Symmetric QCD.

We introduce novel relations between the derivatives [∂^{n}ρ(λ,m_{l})/∂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 ∂^{n}ρ(λ,m_{l})/∂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.06fm, and lattices having aspect ratios 4-9. We find that ρ(λ→0,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→0, and its amplitude is proportional to m_{l}^{2}. We demonstrate that this ρ(λ→0,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.

Predictions for the beauty meson spectrum

We predict the spectrum of the four $1S$ and eight $1P$ nonstrange and strange states in the beauty meson family in the context of effective field theory. By the union of heavy quark effective theory and chiral perturbation theory, the mass formalisms for the heavy-light mesons are defined. Our analysis uses mass expressions involving, for the first time, the full leading self-energy corrections and leading power corrections to the heavy quark and chiral limits. The counterterms present in these expressions are fitted using available experimental and theoretical information on the charmed meson masses and couplings. The results from charm sector are used to make theoretical expectations on the analog beauty meson spectrum. The observed spectrum of the ground state, $B^{(*)}_{(s)}$, and excited, $B_{(s)1}$ and $B^*_{(s)2}$, beauty mesons are well reproduced in our theoretical calculations. The excited scalar, $B^*_{(s)0}$, and axial-vector, $B^\prime_{(s)1}$, beauty mesons have not yet been discovered. Hopefully our predictions may provide valuable clues to further experimental exploration of these missing resonances.

Dislocations under gradient flow and their effect on the renormalized coupling

Non-zero topological charge is prohibited in the chiral limit of gauge-fermion systems because any instanton would create a zero mode of the Dirac operator. On the lattice, however, the geometric $Q_\text{geom}=\langle F{\tilde F}\rangle /32\pi^2$ definition of the topological charge does not necessarily vanish even when the gauge fields are smoothed for example with gradient flow. Small vacuum fluctuations (dislocations) not seen by the fermions may be promoted to instanton-like objects by the gradient flow. We demonstrate that these artifacts of the flow cause the gradient flow renormalized gauge coupling to increase and run faster. In step-scaling studies such artifacts contribute a term which increases with volume. The usual $a/L\to 0$ continuum limit extrapolations can hence lead to incorrect results. In this paper we investigate these topological lattice artifacts in the SU(3) 10-flavor system with domain wall fermions and the 8-flavor system with staggered fermions. Both systems exhibit nonzero topological charge at the strong coupling, especially when using Symanzik gradient flow. We demonstrate how this artifact impacts the determination of the renormalized gauge coupling and the step-scaling $\beta$ function.

Sensitivity of the Polyakov loop and related observables to chiral symmetry restoration

While the Polyakov loop is an order parameter of the deconfinement transition in the heavy quark mass regime of QCD, its sensitivity to the deconfinement of light, dynamical quarks in QCD is not apparent. On the other hand, the quark mass dependence of the Polyakov loop is sensitive to the appearance of a chiral phase transition. Using lattice QCD calculations in the staggered fermion discretization scheme at finite values of the lattice spacing, $aT = 1/8$, we show here, for the first time, that the Polyakov loop expectation value, and the heavy quark free energy extracted from it, behave like energy-like observables in the vicinity of the chiral phase transition temperature $T_c$. Consistent with scaling behavior of energy-like observables in the 3-$d$, O(2) universality class, the quark mass derivatives diverge in the chiral limit at $T_c$ while the temperature derivatives stay finite. The latter will develop a characteristic spike at $T_c$. This, however, may be resolved only in calculations with quark masses being two orders of magnitude smaller than those currently accessible in lattice QCD calculations.

Axial UA(1) Anomaly: A New Mechanism to Generate Massless Bosons

Prior to the establishment of QCD as the correct theory describing hadronic physics, it was realized that the essential ingredients of the hadronic world at low energies are chiral symmetry and its spontaneous breaking. Spontaneous symmetry breaking is a non-perturbative phenomenon, and, thanks to massive QCD simulations on the lattice, we have at present a good understanding of the vacuum realization of the non-abelian chiral symmetry as a function of the physical temperature. As far as the UA(1) anomaly is concerned, and especially in the high temperature phase, the current situation is however far from satisfactory. The first part of this article is devoted to reviewing the present status of lattice calculations, in the high temperature phase of QCD, of quantities directly related to the UA(1) axial anomaly. In the second part, some recently suggested interesting physical implications of the UA(1) anomaly in systems where the non-abelian axial symmetry is fulfilled in the vacuum are analyzed. More precisely it is argued that, if the UA(1) symmetry remains effectively broken, the topological properties of the theory can be the basis of a mechanism, other than Goldstone’s theorem, to generate a rich spectrum of massless bosons at the chiral limit.