Search

We solve a random matrix model for QCD at finite chemical potential, obtained by generalizing the Stephanov model by modifying the random-matrix integration measure with a one-parameter trace deformation. This allows one to check how important the integration measure is for the qualitative features of random matrix models, as well as to test the robustness and universality of the qualitative picture of the original model. While for a small trace deformation the phase diagram is identical to that of the Stephanov model, for a large deformation an exotic phase with spontaneous charge-conjugation breaking appears., Comment: 18 pages, 7 figures

We consider the problem of making nonparametric inference in multi-dimensional diffusion models from low-frequency data. Statistical analysis in this setting is notoriously challenging due to the intractability of the likelihood and its gradient, and computational methods have thus far largely resorted to expensive simulation-based techniques. In this article, we propose a new computational approach which is motivated by PDE theory and is built around the characterisation of the transition densities as solutions of the associated heat (Fokker-Planck) equation. Employing optimal regularity results from the theory of parabolic PDEs, we prove a novel characterisation for the gradient of the likelihood. Using these developments, for the nonlinear inverse problem of recovering the diffusivity (in divergence form models), we then show that the numerical evaluation of the likelihood and its gradient can be reduced to standard elliptic eigenvalue problems, solvable by powerful finite element methods. This enables the efficient implementation of a large class of statistical algorithms, including (i) preconditioned Crank-Nicolson and Langevin-type methods for posterior sampling, and (ii) gradient-based descent optimisation schemes to compute maximum likelihood and maximum-a-posteriori estimates. We showcase the effectiveness of these methods via extensive simulation studies in a nonparametric Bayesian model with Gaussian process priors. Interestingly, the optimisation schemes provided satisfactory numerical recovery while exhibiting rapid convergence towards stationary points despite the problem nonlinearity; thus our approach may lead to significant computational speed-ups. The reproducible code is available online at https://github.com/MattGiord/LF-Diffusion., Comment: 50 pages, 11 figures, 5 tables

I derive constraints on the Dirac spectrum in the chirally symmetric phase of a gauge theory with two massless fermion flavors. Using only general properties of correlation functions of scalar and pseudoscalar bilinears, I prove that in the chiral limit of vanishing quark mass $m$ all the corresponding susceptibilities must be power series in $m^2$ with finite coefficients, from which spectral constraints follow. I then use them to show that effective breaking of the anomalous $\mathrm{U}(1)_A$ symmetry is allowed in the $\mathrm{SU}(2)_A$ symmetric phase in the chiral limit, and leads to distinctive spectral features: (i) the spectral density must develop a singular $m^4/\lambda$ peak as $m\to 0$, (ii) the two-point eigenvalue correlator of near-zero modes must be singular, and (iii) near-zero modes cannot be localized. Moreover, in the symmetric phase the topological charge distribution must be indistinguishable from that of an ideal gas of instantons and anti-instantons of vanishing density, to leading order in the fermion mass., Comment: 13 pages

We consider the statistical linear inverse problem of making inference on an unknown source function in an elliptic partial differential equation from noisy observations of its solution. We employ nonparametric Bayesian procedures based on Gaussian priors, leading to convenient conjugate formulae for posterior inference. We review recent results providing theoretical guarantees on the quality of the resulting posterior-based estimation and uncertainty quantification, and we discuss the application of the theory to the important classes of Gaussian series priors defined on the Dirichlet-Laplacian eigenbasis and Mat\'ern process priors. We provide an implementation of posterior inference for both classes of priors, and investigate its performance in a numerical simulation study., Comment: 21 Pages, 8 figures, 5 tables. To appear in BAYSM 2023 proceedings

This work studies nonparametric Bayesian estimation of the intensity function of an inhomogeneous Poisson point process in the important case where the intensity depends on covariates, based on the observation of a single realisation of the point pattern over a large area. It is shown how the presence of covariates allows to borrow information from far away locations in the observation window, enabling consistent inference in the growing domain asymptotics. In particular, optimal posterior contraction rates under both global and point-wise loss functions are derived. The rates in global loss are obtained under conditions on the prior distribution resembling those in the well established theory of Bayesian nonparametrics, here combined with concentration inequalities for functionals of stationary processes to control certain random covariate-dependent loss functions appearing in the analysis. The local rates are derived with an ad-hoc study that builds on recent advances in the theory of P\'olya tree priors, extended to the present multivariate setting with a novel construction that makes use of the random geometry induced by the covariates., Comment: 55 pages

We study the effects of taste degeneracy on the continuum scaling of the localization properties of the staggered Dirac operator in high-temperature QCD using numerical simulations on the lattice, focusing in particular on the position of the mobility edge separating localized and delocalized modes at the low end of the spectrum. We find that, if the continuum limit is approached at fixed spatial volume, the restoration of taste symmetry leads to sizeable systematic effects on estimates for the mobility edge obtained from spectral statistics, which become larger and larger as the lattice spacing is decreased. Such systematics, however, are found to decrease if the volume is increased at fixed lattice spacing. We argue that spectral statistics estimate correctly the position of the mobility edge in the thermodynamic limit at fixed spacing, and support this with an independent numerical analysis based directly on the properties of the Dirac eigenvectors, that are unaffected by taste degeneracy. We then provide a theoretical argument justifying the exchange of the thermodynamic and continuum limits when studying localization. This allows us to use spectral statistics to determine the position of the mobility edge, and to obtain a controlled continuum extrapolation of the mobility edge for the first time., Comment: Revised and expanded discussion, includes argument for exchange of thermodynamic and continuum limits for the mobility edge; numerical results unchanged; to appear in PRD; 21 pages, 12 figures

We investigate the connection between localization of low-lying Dirac modes and Polyakov-loop ordering in the lattice $\mathrm{SU}(2)$ Higgs model at finite temperature, probed with the staggered Dirac operator. After mapping out the phase diagram of the model at a fixed temporal extension in lattice units, we study the localization properties of the low-lying modes of the staggered Dirac operator, how these properties change across the various transitions, and how these modes correlate with the gauge and Higgs fields. We find localized low modes in the deconfined and in the Higgs phase, where the Polyakov loop is strongly ordered, but in both cases they disappear as one crosses over to the confined phase. Our findings confirm the general expectations of the "sea/islands" picture, and the more detailed expectations of its refined version concerning the favorable locations of localized modes, also in the presence of dynamical scalar matter., Comment: Revised version, matches published version; 21 pages, 22 figures

Research on the QCD phase diagram with lattice field theory methods is dominated by the use of rooted staggered fermions, as they are the computationally cheapest discretization available. We show that rooted staggered fermions at a nonzero baryochemical potential $\mu_B$ predict a sharp rise in the baryon density at low temperatures and $\mu_B \gtrsim 3 m_\pi/2$, where $m_\pi$ is the Goldstone pion mass. We elucidate the nature of the non-analyticity behind this sharp rise in the density by a comparison of reweighting results with a Taylor expansion of high order. While at first sight this non-analytic behavior becomes apparent at the same position where the pion condensation transition takes place in the phase-quenched theory, the nature of the non-analyticity in the two theories appears to be quite different: While at nonzero isospin density the data are consistent with a genuine thermodynamic (branch-point) singularity, the results at nonzero baryon density point to an essential singularity at $\mu_B=0$. The effect is absent for four flavors of degenerate quarks, where rooting is not used. For the two-flavor case, we show numerical evidence that the magnitude of the effect diminishes on finer lattices. We discuss the implications of this technical complication on future studies of the QCD phase diagram., Comment: 12 pages, 6 figures

I show that spontaneous breaking of vector flavor symmetry on the lattice is impossible in gauge theories with a positive functional-integral measure, for discretized Dirac operators linear in the quark masses, if the corresponding propagator and its commutator with the flavor symmetry generators can be bounded in norm independently of the gauge configuration and uniformly in the volume. Under these assumptions, any order parameter vanishes in the symmetric limit of fermions of equal masses. I show that these assumptions are satisfied by staggered, minimally doubled and Ginsparg-Wilson fermions for positive fermion mass, for any value of the lattice spacing, and so in the continuum limit if this exists. They are instead not satisfied by Wilson fermions, for which spontaneous vector flavor symmetry breaking is known to take place in the Aoki phase. The existence of regularizations unaffected by residual fermion doubling for which the symmetry cannot break spontaneously on the lattice establishes rigorously (at the physicist's level) the impossibility of its spontaneous breaking in the continuum for any number of flavors., Comment: Revised version: title changed; added references on the rooting trick; added discussion on improved and on rooted staggered fermions; 16 pages

We studied integration contour deformations in the chiral random matrix theory of Stephanov with the goal of alleviating the finite-density sign problem. We considered simple ans\"atze for the deformed integration contours, and optimized their parameters. We find that optimization of a single parameter manages to considerably improve on the severity of the sign problem. We show numerical evidence that the improvement achieved is exponential in the degrees of freedom of the system, i.e., the size of the random matrix. We also compare the optimization method with contour deformations coming from the holomorphic flow equations., Comment: 10 pages, 7 figures

Ample numerical evidence from lattice calculations shows a strong connection between the confining properties of gauge theories at finite temperature and the localisation properties of the low-lying Dirac eigenmodes. In this contribution we discuss recent progress on this topic, focussing on results for QCD at imaginary chemical potential ${\mu}_I/T = \pi$ at temperatures above the Roberge-Weiss transition temperature. These confirm the general picture of low modes turning from delocalised to localised at the deconfinement transition, in a previously unexplored setup with a genuine, physical transition in the presence of dynamical fermions. This further supports the use of Dirac eigenmodes as a tool to investigate the mechanisms behind confinement and the deconfinement transition., Comment: 8 pages, 4 figures; contribution to the proceedings of the XVth Quark Confinement and Hadron Spectrum conference (ConfXV), 1-6 August 2022, University of Stavanger, Stavanger (Norway)

We study the localization properties of the eigenmodes of the staggered Dirac operator across the deconfinement transition in finite-temperature $\mathbb{Z}_3$ pure gauge theory on the lattice in 2+1 dimensions. This allows for nontrivial tests of the sea-islands picture of localization, according to which low modes should localize on favorable Polyakov-loop fluctuations in the deconfined phase of a gauge theory. We observe localized low modes in the deconfined phase of the theory, both in the real Polyakov-loop sector, where they are expected, and in the complex Polyakov-loop sectors, where they are not. Our findings expose the limitations of the standard sea-islands picture, and call for its refinement. An improved picture, where spatial hopping terms play a more prominent role, is proposed and found to be in excellent agreement with numerical results., Comment: 22 pages, 10 figures, revised version, to be published in PRD

Besov priors are nonparametric priors that can model spatially inhomogeneous functions. They are routinely used in inverse problems and imaging, where they exhibit attractive sparsity-promoting and edge-preserving features. A recent line of work has initiated the study of their asymptotic frequentist convergence properties. In the present paper, we consider the theoretical recovery performance of the posterior distributions associated to Besov-Laplace priors in the density estimation model, under the assumption that the observations are generated by a possibly spatially inhomogeneous true density belonging to a Besov space. We improve on existing results and show that carefully tuned Besov-Laplace priors attain optimal posterior contraction rates. Furthermore, we show that hierarchical procedures involving a hyper-prior on the regularity parameter lead to adaptation to any smoothness level., Comment: 35 pages, to appear in the Electronic Journal of Statistics

The equation of state of the quark gluon plasma is a key ingredient of heavy ion phenomenology. In addition to the traditional Taylor method, several novel approximation schemes have been proposed with the aim of calculating it at finite baryon density. In order to gain a pragmatic understanding of the limits of these schemes, we compare them to direct results at $\mu_B>0$, using reweighting techniques free from an overlap problem. We use 2stout improved staggered fermions with 8 time-slices and cover the entire RHIC BES range in the baryochemical potential, up to $\mu_B/T=3$., Comment: 7 pages, 3 figures

I discuss the possible effects of a finite density of localised near-zero Dirac modes in the chiral limit of gauge theories with $N_f$ degenerate fermions. I focus in particular on the fate of the massless quasi-particle excitations predicted by the finite-temperature version of Goldstone's theorem, for which I provide an alternative and generalised proof based on a Euclidean ${\rm SU}(N_f)_A$ Ward-Takahashi identity. I show that localised near-zero modes can lead to a divergent pseudoscalar-pseudoscalar correlator that modifies this identity in the chiral limit. As a consequence, massless quasi-particle excitations can disappear from the spectrum of the theory in spite of a non-zero chiral condensate. Three different scenarios are possible, depending on the detailed behaviour in the chiral limit of the ratio of the mobility edge and the fermion mass, which I prove to be a renormalisation-group invariant quantity., Comment: 1+72 pages (40 + appendices and bibliography), 1 figure; to be published in JHEP

We rigorously prove that deep Gaussian process priors can outperform Gaussian process priors if the target function has a compositional structure. To this end, we study information-theoretic lower bounds for posterior contraction rates for Gaussian process regression in a continuous regression model. We show that if the true function is a generalized additive function, then the posterior based on any mean-zero Gaussian process can only recover the truth at a rate that is strictly slower than the minimax rate by a factor that is polynomially suboptimal in the sample size $n$., Comment: 20 pages, to appear in Advances in Neural Information Processing Systems 36 (NeurIPS 2022)

We study the 2+1 dimensional XY model at nonzero chemical potential $\mu$ on deformed integration manifolds, with the aim of alleviating its sign problem. We investigate several proposals for the deformations, and considerably improve on the severity of the sign problem with respect to standard reweighting approaches. We present numerical evidence that the reduction of the sign problem is exponential both in $\mu^2$ and in the spatial volume. We also present a new approach to the optimization procedure based on reweighting, that sensibly reduces its computational cost., Comment: 14 pages, 7 figures

Most lattice studies of hot and dense QCD matter rely on extrapolation from zero or imaginary chemical potentials. The ill-posedness of numerical analytic continuation puts severe limitations on the reliability of such methods. We studied the QCD chiral transition at finite real baryon density with the more direct sign reweighting approach. We simulate up to a baryochemical potential-temperature ratio of $\mu_B/T=2.7$, covering the RHIC Beam Energy Scan range, and penetrating the region where methods based on analytic continuation are unpredictive.This opens up a new window to study QCD matter at finite $\mu_B$ from first principles., Comment: 10 pages, 3 figures; Contribution to the XXXIII International (ONLINE) Workshop on High Energy Physics "Hard Problems of Hadron Physics: Non-Perturbative QCD & Related Quests"; Based on 2108.09213 [hep-lat]. arXiv admin note: substantial text overlap with arXiv:2112.02134

Approaches to finite baryon density lattice QCD usually suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We test a method - sign reweighting - that works directly at finite chemical potential and is yet free from any such uncontrolled systematics: with this approach the only problem is the sign problem itself. In practice the approach involves the generation of configurations with the positive fermionic weights given by the absolute value of the real part of the quark determinant, and a reweighting by a sign. There are only two sectors, +1 and -1 and as long as the average $\left\langle \pm \right\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting has no overlap problem - unlike reweighting from $\mu=0$ - and the results are reliable. We also present results based on this algorithm on the phase diagram of lattice QCD with two different actions: as a first test, we apply the method to calculate the position of the critical endpoint with unimproved staggered fermions at $N_\tau=4$; as a second application, we study the phase diagram with 2stout improved staggered fermions at $N_\tau=6$. This second one is already a reasonably fine lattice - relevant for phenomenology. We demonstrate that the method penetrates the region of the phase diagram where the Taylor and imaginary chemical potential methods lose predictive power., Comment: 9 pages, 4 figures; Contribution to the Proceedings of The 38th International Symposium on Lattice Field Theory, LATTICE2021; Based on 2004.10800 and 2108.09213

The low-lying Dirac modes become localised at the finite-temperature transition in QCD and in other gauge theories, suggesting a general connection between their localisation and deconfinement. The simplest model where this connection can be tested is $\mathbb{Z}_2$ gauge theory in 2+1 dimensions. We show that in this model the low modes in the staggered Dirac spectrum are delocalised in the confined phase and become localised in the deconfined phase. We also show that localised modes correlate with disorder in the Polyakov loop configuration, in agreement with the "sea/island" picture of localisation, and with negative plaquettes. These results further support the conjecture that localisation and deconfinement are closely related., Comment: 8 pages, 5 figures, contribution to the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021

I show that a finite density of near-zero localised Dirac modes can lead to the disappearance of the massless excitations predicted by the finite-temperature version of Goldstone's theorem in the chirally broken phase of a gauge theory., Comment: 8 pages; contribution to the 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We study the localization properties of the low-lying Dirac eigenmodes in QCD at imaginary chemical potential $\hat{\mu}_I=\pi$ at temperatures above the Roberge-Weiss transition temperature $T_{\rm RW}$. We find that modes are localized up to a temperature-dependent "mobility edge" and delocalized above it, and that the mobility edge extrapolates to zero at a temperature compatible with $T_{\rm RW}$. This supports the existence of a strong connection between localization of the low Dirac modes and deconfinement, studied here for the first time in a model with a genuine deconfinement transition in the continuum limit in the presence of dynamical fermions., Comment: 8 pages, 5 figures

State-of-the-art lattice QCD studies of hot and dense strongly interacting matter currently rely on extrapolation from zero or imaginary chemical potentials. The ill-posedness of numerical analytic continuation puts severe limitations on the reliability of such methods. Here we use the more direct sign reweighting method to perform lattice QCD simulation of the QCD chiral transition at finite real baryon density on phenomenologically relevant lattices. This method does not require analytic continuation and avoids the overlap problem associated with generic reweighting schemes, so has only statistical but no uncontrolled systematic uncertainties for a fixed lattice setup. This opens up a new window to study hot and dense strongly interacting matter from first principles. We perform simulations up to a baryochemical potential-temperature ratio of $\mu_B/T=2.5$ covering most of the RHIC Beam Energy Scan range in the chemical potential. We also clarify the connection of the approach to the more traditional phase reweighting method., Comment: 11 pages, 6 figures

It is by now well established that Dirac fermions coupled to non-Abelian gauge theories can undergo an Anderson-type localization transition. This transition affects eigenmodes in the lowest part of the Dirac spectrum, the ones most relevant to the low-energy physics of these models. Here we review several aspects of this phenomenon, mostly using the tools of lattice gauge theory. In particular, we discuss how the transition is related to the finite-temperature transitions leading to the deconfinement of fermions, as well as to the restoration of chiral symmetry that is spontaneously broken at low temperature. Other topics we touch upon are the universality of the transition, and its connection to topological excitations (instantons) of the gauge field and the associated fermionic zero modes. While the main focus is on Quantum Chromodynamics, we also discuss how the localization transition appears in other related models with different fermionic contents (including the quenched approximation), gauge groups, and in different space-time dimensions. Finally, we offer some speculations about the physical relevance of the localization transition in these models., Comment: Invited review for the special issue of Universe "Modern Approaches to Non-Perturbative QCD and other Confining Gauge Theories", ed. D. Antonov. Revised version, 47 pages, 17 figures

We study, by means of numerical lattice simulations, the properties of the reconfinement phase transition taking place in trace deformed $SU(3)$ Yang-Mills theory defined on $\mathbb{R}^3\times S^1$, in which center symmetry is recovered even for small compactification radii. We show, by means of a finite size scaling analysis, that the reconfinement phase transition is first-order, like the usual $SU(3)$ thermal phase transition. We then investigate two different physical phenomena, which are known to characterize the standard confinement/deconfinement phase transition, namely the condensation of thermal magnetic monopoles and the change in the localization properties of the eigenmodes of the Dirac operator. Regarding the latter, we show that the mobility edge signalling the Anderson-like transition in the Dirac spectrum vanishes as one enters the reconfined phase, as it happens in the standard confined phase. Thermal monopoles, instead, show a peculiar behavior: their density decreases going through reconfinement, at odds with the standard thermal theory; nonetheless, they condense at reconfinement, like at the usual confinement transition. The coincidence of monopole condensation and Dirac mode delocalization, even in a framework different from that of the standard confinement transition, suggests the existence of a strict link between them., Comment: 10 pages, 10 figures

We study nonparametric Bayesian models for reversible multi-dimensional diffusions with periodic drift. For continuous observation paths, reversibility is exploited to prove a general posterior contraction rate theorem for the drift gradient vector field under approximation-theoretic conditions on the induced prior for the invariant measure. The general theorem is applied to Gaussian priors and $p$-exponential priors, which are shown to converge to the truth at the minimax optimal rate over Sobolev smoothness classes in any dimension., Comment: 41 pages, 1 figure, to appear in the Annals of Statistics

I show that a finite density of near-zero localised Dirac modes in the chirally broken phase of a gauge theory can lead to the disappearance of the massless excitations predicted by the Goldstone theorem at finite temperature., Comment: Revised version: title changed; error corrected in Eq. 11, conclusions unchanged; discussion improved, comments added about symmetry breaking and relation with Goldstone's theorem; published version, 13 pages

All approaches currently used to study finite baryon density lattice QCD suffer from uncontrolled systematic uncertainties in addition to the well-known sign problem. We formulate and test an algorithm, sign reweighting, that works directly at finite $\mu = \mu_B/3$ and is yet free from any such uncontrolled systematics. With this algorithm the {\em only} problem is the sign problem itself. This approach involves the generation of configurations with the positive fermionic weight $|{\rm Re\; det} D(\mu)|$ where $D(\mu)$ is the Dirac matrix and the signs ${\rm sign} \; ( {\rm Re\; det} D(\mu) ) = \pm 1$ are handled by a discrete reweighting. Hence there are only two sectors, $+1$ and $-1$ and as long as the average $\langle\pm 1\rangle \neq 0$ (with respect to the positive weight) this discrete reweighting by the signs carries no overlap problem and the results are reliable. The approach is tested on $N_t = 4$ lattices with $2+1$ flavors and physical quark masses using the unimproved staggered discretization. By measuring the Fisher (sometimes also called Lee-Yang) zeros in the bare coupling on spatial lattices $L/a = 8, 10, 12$ we conclude that the cross-over present at $\mu = 0$ becomes stronger at $\mu > 0$ and is consistent with a true phase transition at around $\mu_B/T \sim 2.4$., Comment: 13 pages, 17 figures, references added

We perform the first direct determination of the position of the leading singularity of the pressure in the complex chemical potential $\mu_B$ plane in lattice QCD using numerical simulations with 2-stout improved rooted staggered fermions. This provides a direct determination of the radius of convergence of the Taylor expansion of the pressure that does not rely on a finite-order truncation of the expansion. The analyticity issues in the complex $\mu_B$ plane of the grand canonical partition function of QCD with rooted staggered fermions are solved with a careful redefinition of the fermion determinant that makes it a polynomial in the fugacity on any finite lattice, without changing the continuum limit of the observables. By performing a finite volume scaling study at a single coarse lattice spacing, we show that the limiting singularity is not on the real line in the thermodynamic limit, thus showing that the radius of convergence of the Taylor expansion gives a lower bound on the location of a possible phase transition. In the vicinity of the crossover temperature at zero chemical potential, the radius of convergence turns out to be $\mu_B/T \approx 2$ and roughly temperature independent., Comment: 4 pages, 1 figure; Quark Matter 2019 - the XXVIIIth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions

The phase diagram and the location of the critical endpoint (CEP) of lattice QCD with unimproved staggered fermions on a $N_t=4$ lattice was determined fifteen years ago with the multiparameter reweighting method by studying Fisher zeros. We first reproduce the old result with an exact algorithm (not known at the time) and with statistics larger by an order of magnitude. As an extension of the old analysis we introduce stout smearing in the fermion action in order to reduce the finite lattice spacing effects. First we show that increasing the smearing parameter $\rho$ the crossover at $\mu = 0$ gets weaker, i.e., the leading Fisher zero gets farther away from the real axis. Furthermore as the chemical potential is increased the overlap problem gets severe sooner than in the unimproved case, therefore shrinking the range of applicability of the method. Nevertheless certain qualitative features remain, even after introducing the smearing. Namely, at small chemical potentials the Fisher zeros first get farther away from the real axis and later at around $a\mu _q = 0.1 - 0.15$ they start to get closer, i.e., the crossover first gets weaker and later stronger as a function of $\mu$. However, because of the more severe overlap problem the CEP is out of reach with the smeared action., Comment: 8 pages, 5 figures

In typical statistical mechanical systems the grand canonical partition function at finite volume is proportional to a polynomial of the fugacity $e^{\mu/T}$. The zero of this Lee-Yang polynomial closest to the origin determines the radius of convergence of the Taylor expansion of the pressure around $\mu=0$. The computationally cheapest formulation of lattice QCD, rooted staggered fermions, with the usual definition of the rooted determinant, does not admit such a Lee-Yang polynomial. We argue that the radius of convergence is then bounded by the spectral gap of the reduced matrix of the unrooted staggered operator. This is a cutoff effect that potentially affects all estimates of the radius of convergence with the standard staggered rooting. We suggest a new definition of the rooted staggered determinant at finite chemical potential that allows for a definition of a Lee-Yang polynomial, and, therefore of the numerical study of Lee-Yang zeros. We also describe an algorithm to determine the Lee-Yang zeros and apply it to configurations generated with the 2-stout improved staggered action at $N_t = 4$. We perform a finite-volume scaling study of the leading Lee-Yang zeros and estimate the radius of convergence of the Taylor expansion extrapolated to an infinite volume. We show that the limiting singularity is not on the real line, thus giving a lower bound on the location of any possible phase transitions at this lattice spacing. In the vicinity of the crossover temperature at zero chemical potential, the radius of convergence turns out to be $\mu_B/T \approx 2$ and roughly temperature independent. Our simulations are performed at strange quark chemical potential $\mu_s=0$, but the method can be straightforwardly extended to strangeness chemical potential $\mu_S=0$ or strangeness neutrality., Comment: 12 pages, 11 figures; v2: contains corrections to the formulas from the erratum PRD 104, 119901(E) (2021); results unchanged

For $\mathcal{O}$ a bounded domain in $\mathbb{R}^d$ and a given smooth function $g:\mathcal{O}\to\mathbb{R}$, we consider the statistical nonlinear inverse problem of recovering the conductivity $f>0$ in the divergence form equation $$ \nabla\cdot(f\nabla u)=g\ \textrm{on}\ \mathcal{O}, \quad u=0\ \textrm{on}\ \partial\mathcal{O}, $$ from $N$ discrete noisy point evaluations of the solution $u=u_f$ on $\mathcal O$. We study the statistical performance of Bayesian nonparametric procedures based on a flexible class of Gaussian (or hierarchical Gaussian) process priors, whose implementation is feasible by MCMC methods. We show that, as the number $N$ of measurements increases, the resulting posterior distributions concentrate around the true parameter generating the data, and derive a convergence rate $N^{-\lambda}, \lambda>0,$ for the reconstruction error of the associated posterior means, in $L^2(\mathcal{O})$-distance., Comment: 34 pages, to appear in Inverse Problems

We introduce a continuous-time Markov chain describing dynamic allelic partitions which extends the branching process construction of the Pitman sampling formula in Pitman (2006) and the birth-and-death process with immigration studied in Karlin and McGregor (1967), in turn related to the celebrated Ewens sampling formula. A biological basis for the scheme is provided in terms of a population of individuals grouped into families, that evolves according to a sequence of births, deaths and immigrations. We investigate the asymptotic behaviour of the chain and show that, as opposed to the birth-and-death process with immigration, this construction maintains in the temporal limit the mutual dependence among the multiplicities. When the death rate exceeds the birth rate, the system is shown to have reversible distribution identified as a mixture of Pitman sampling formulae, with negative binomial mixing distribution on the population size. The population therefore converges to a stationary random configuration, characterised by a finite number of families and individuals., Comment: 17 pages, to appear in ALEA , Latin American Journal of Probability and Mathematical Statistics

We investigate the QCD phase diagram for nonzero background magnetic fields using first-principles lattice simulations. At the physical point (in terms of quark masses), the thermodynamics of this system is controlled by two opposing effects: magnetic catalysis (enhancement of the quark condensate) at low temperature and inverse magnetic catalysis (reduction of the condensate) in the transition region. While the former is known to be robust and independent of the details of the interactions, inverse catalysis arises as a result of a delicate competition, effective only for light quarks. By performing simulations at different quark masses, we determine the pion mass above which inverse catalysis does not take place in the transition region anymore. Even for pions heavier than this limiting value - where the quark condensate undergoes magnetic catalysis - our results are consistent with the notion that the transition temperature is reduced by the magnetic field. These findings will be useful to guide low-energy models and effective theories of QCD., Comment: Revised version; matches published version; 14 pages, 7 figures

We study different estimators of the radius of convergence of the Taylor series of the pressure in finite density QCD. We adopt the approach in which the radius of convergence is estimated first in a finite volume, and the infinite-volume limit is taken later. This requires an estimator for the radius of convergence that is reliable in a finite volume. Based on general arguments about the analytic structure of the partition function in a finite volume, we demonstrate that the ratio estimator cannot work in this approach, and propose three new estimators, capable of extracting reliably the radius of convergence, which coincides with the distance from the origin of the closest Lee-Yang zero. We also provide an estimator for the phase of the closest Lee-Yang zero, necessary to assess whether the leading singularity is a true critical point. We demonstrate the usage of these estimators on a toy model, namely 4 flavors of unimproved staggered fermions on a small $6^3 \times 4$ lattice, where both the radius of convergence and the Taylor coefficients to any order can be obtained by a direct determination of the Lee-Yang zeros. Interestingly, while the relative statistical error of the Taylor expansion coefficients steadily grows with order, that of our estimators stabilizes, allowing for an accurate estimate of the radius of convergence. In particular, we show that despite the more than 100\% error bars on high-order Taylor coefficients, the given ensemble contains enough information about the radius of convergence., Comment: 14 pages, 11 figures

I study the localisation properties of low Dirac eigenmodes in 2+1 dimensional SU(3) pure gauge theory, both in the low-temperature, confined and chirally-broken phase and in the high-temperature, deconfined and chirally-restored phase, by means of numerical lattice simulations. While these modes are delocalised at low temperature, they become localised at high temperature, up to a critical point in the Dirac spectrum where a BKT-type Anderson transition takes place. All results point to localisation appearing at the deconfinement temperature, and support previous expectations about the close relation between deconfinement, chiral symmetry breaking, and localisation., Comment: 1+33 pages, 16 figures; several minor changes (typos fixed, a few remarks and extra references added); matches published version

I discuss recent results on the relation between the localisation of low-lying Dirac eigenmodes, the restoration of chiral symmetry, and deconfinement in QCD and QCD-like models, providing evidence of a close connection between the three phenomena., Comment: 12 pages, 6 figures. Contribution to the XIIIth Quark Confinement and the Hadron Spectrum conference (Confinement 2018), 31 July - 6 August 2018, Maynooth University (Ireland). v2: improved figures

We consider the statistical inverse problem of recovering an unknown function $f$ from a linear measurement corrupted by additive Gaussian white noise. We employ a nonparametric Bayesian approach with standard Gaussian priors, for which the posterior-based reconstruction of $f$ corresponds to a Tikhonov regulariser $\bar f$ with a reproducing kernel Hilbert space norm penalty. We prove a semiparametric Bernstein-von Mises theorem for a large collection of linear functionals of $f$, implying that semiparametric posterior estimation and uncertainty quantification are valid and optimal from a frequentist point of view. The result is applied to study three concrete examples that cover both the mildly and severely ill-posed cases: specifically, an elliptic inverse problem, an elliptic boundary value problem and the heat equation. For the elliptic boundary value problem, we also obtain a nonparametric version of the theorem that entails the convergence of the posterior distribution to a prior-independent infinite-dimensional Gaussian probability measure with minimal covariance. As a consequence, it follows that the Tikhonov regulariser $\bar f$ is an efficient estimator of $f$, and we derive frequentist guarantees for certain credible balls centred at $\bar{f}$., Comment: 34 pages, to appear in SIAM/ASA Journal on Uncertainty Quantification (JUQ)

We discuss the reliability of available methods to constrain the location of the QCD critical endpoint with lattice simulations. In particular we calculate the baryon fluctuations up to $\chi^B_8$ using simulations at imaginary chemical potentials. We argue that they contain no hint of criticality., Comment: XXVIIth International Conference on Ultrarelativistic Nucleus-Nucleus Collisions (Quark Matter 2018)

We compute the Euclidean correlators of the stress tensor in pure $SU(3)$ Yang-Mills theory at finite temperature at zero and finite spatial momenta with lattice simulations. We perform continuum extrapolations using $N_\tau=10,12,16,20$ lattices with renormalized anisotropy 2. We use these correlators to estimate the shear viscosity of the gluon plasma in the deconfined phase. For $T=1.5T_c$ we obtain $\eta/s=0.17(2)$.

We will discuss the issue of Landau levels of quarks in lattice QCD in an external magnetic field. We will show that in the two-dimensional case the lowest Landau level can be identified unambiguously even if the strong interactions are turned on. Starting from this observation, we will then show how one can define a "lowest Landau level" in the four-dimensional case, and discuss how much of the observed effects of a magnetic field can be explained in terms of it. Our results can be used to test the validity of low-energy models of QCD that make use of the lowest-Landau-level approximation., Comment: 8 pages, 5 figures; to appear in the proceedings of the 35th International Symposium on Lattice Field Theory (Lattice 2017), June 19-24, 2017, Granada, Spain

The Dirac operator in finite temperature QCD is equivalent to the Hamiltonian of an unconventional Anderson model, with on-site noise provided by the fluctuations of the Polyakov lines. The main features of its spectrum and eigenvectors, concerning the density of low modes and their localization properties, are qualitatively reproduced by a toy-model random Hamiltonian, based on an Ising-type spin model mimicking the dynamics of the Polyakov lines. Here we study the low modes of this toy model in the vicinity of the ordering transition of the spin model, and show that at the critical point the spectral density at the origin has a singularity, and the localization properties of the lowest modes change. This provides further evidence of the close relation between deconfinement, chiral transition and localization of the low modes., Comment: 11 pages, 7 figures

The spectrum of the two-dimensional continuum Dirac operator in the presence of a uniform background magnetic field consists of Landau levels, which are degenerate and separated by gaps. On the lattice the Landau levels are spread out by discretization artefacts, but a remnant of their structure is clearly visible (Hofstadter butterfly). If one switches on a non-Abelian interaction, the butterfly structure will be smeared out, but the lowest Landau level (LLL) will still be separated by a gap from the rest of the spectrum. In this talk we discuss how one can define the LLL in QCD and check how well certain physical quantities are approximated by taking into account only the LLL., Comment: 1+6 pages, 2 figures. Contribution to the 34th International Symposium on Lattice Field Theory (Lattice 2016), 24-30 July 2016, University of Southampton, UK

We study the problems of deconfinement, chiral symmetry restoration and localisation of the low Dirac eigenmodes in a toy model of QCD, namely unimproved staggered fermions on lattices of temporal extension $N_T=4$. This model displays a genuine deconfining and chirally-restoring first-order phase transition at some critical value of the gauge coupling. Our results indicate that the onset of localisation of the lowest Dirac eigenmodes takes place at the same critical coupling where the system undergoes the first-order phase transition. This provides further evidence of the close relation between deconfinement, chiral symmetry restoration and localisation of the low modes of the Dirac operator on the lattice., Comment: 1+17 pages, 9 figures; revised introduction, minor changes in figures, new references added; matches published version

We discuss chiral symmetry restoration and eigenmode localisation in finite-temperature QCD by looking at the lattice Dirac operator as a random Hamiltonian. We argue that the features of QCD relevant to both phenomena are the presence of order in the Polyakov line configuration, and the correlations that this induces between spatial links across time slices. This ties the fate of chiral symmetry and of localisation of the lowest Dirac eigenmodes to the confining properties of the theory. We then show numerical results obtained in a QCD-inspired Anderson-like toy model, derived by radically simplifying the QCD dynamics while keeping the important features mentioned above. The toy model reproduces all the important qualitative aspects of chiral symmetry breaking and localisation in QCD, thus supporting the central role played by the confinement/deconfinement transition in triggering both phenomena., Comment: 1+6 pages, 2 figures. Contribution to the 34th International Symposium on Lattice Field Theory (Lattice 2016), 24-30 July 2016, University of Southampton, UK