Researcher profile

Matteo Giordano

Matteo Giordano contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2026arXiv

A nonparametric Bayesian analysis of independent and identically distributed observations of covariate-driven Poisson processes

An important task in the statistical analysis of inhomogeneous point processes is to investigate the influence of a set of covariates on the point-generating mechanism. In this article, we consider the nonparametric Bayesian approach to this problem, assuming that $n$ independent and identically distributed realizations of the point pattern and the covariate random field are available. In many applications, different covariates are often vastly diverse in physical nature, resulting in anisotropic intensity functions whose variations along distinct directions occur at different smoothness levels. To model this scenario, we employ hierarchical prior distributions based on multi-bandwidth Gaussian processes. We prove that the resulting posterior distributions concentrate around the ground truth at optimal rate as $n\to\infty$, and achieve automatic adaptation to the anisotropic smoothness. Posterior inference is concretely implemented via a Metropolis-within-Gibbs Markov chain Monte Carlo algorithm that incorporates a dimension-robust sampling scheme to handle the functional component of the proposed nonparametric Bayesian model. Our theoretical results are supported by extensive numerical simulation studies. Further, we present an application to the analysis of a Canadian wildfire dataset.

preprint2022arXiv

Lattice simulations of the QCD chiral transition at real $μ_B$

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 $μ_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 $μ_B$ from first principles.

preprint2022arXiv

Localisation of Dirac eigenmodes and confinement in gauge theories: the Roberge-Weiss transition

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 $μ_I/T = π$ 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.

preprint2022arXiv

Localisation of Dirac modes in gauge theories and Goldstone's theorem at finite temperature

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.

preprint2022arXiv

Radius of convergence in lattice QCD at finite $μ_B$ with rooted staggered fermions

In typical statistical mechanical systems the grand canonical partition function at finite volume is proportional to a polynomial of the fugacity $e^{μ/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 $μ=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 $μ_B/T \approx 2$ and roughly temperature independent. Our simulations are performed at strange quark chemical potential $μ_s=0$, but the method can be straightforwardly extended to strangeness chemical potential $μ_S=0$ or strangeness neutrality.

preprint2021arXiv

Lattice simulations of the QCD chiral transition at real baryon density

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 $μ_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.

preprint2021arXiv

Localization properties of Dirac modes at the Roberge-Weiss phase transition

We study the localization properties of the low-lying Dirac eigenmodes in QCD at imaginary chemical potential $\hatμ_I=π$ 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.

preprint2020arXiv

A reversible allelic partition process and Pitman sampling formula

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.

preprint2020arXiv

Bernstein-von Mises theorems and uncertainty quantification for linear inverse problems

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}$.

preprint2020arXiv

Consistency of Bayesian inference with Gaussian process priors in an elliptic inverse problem

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^{-λ}, λ>0,$ for the reconstruction error of the associated posterior means, in $L^2(\mathcal{O})$-distance.

preprint2020arXiv

New approach to lattice QCD at finite density; results for the critical end point on coarse lattices

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 $μ= μ_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(μ)|$ where $D(μ)$ is the Dirac matrix and the signs ${\rm sign} \; ( {\rm Re\; det} D(μ) ) = \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 $μ= 0$ becomes stronger at $μ> 0$ and is consistent with a true phase transition at around $μ_B/T \sim 2.4$.

preprint2020arXiv

Reconfinement, localization and thermal monopoles in $SU(3)$ trace-deformed Yang-Mills theory

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.

preprint2020arXiv

The effect of stout smearing on the phase diagram from multiparameter reweigthing in lattice QCD

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 $ρ$ the crossover at $μ= 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μ_q = 0.1 - 0.15$ they start to get closer, i.e., the crossover first gets weaker and later stronger as a function of $μ$. However, because of the more severe overlap problem the CEP is out of reach with the smeared action.

preprint2020arXiv

Towards a reliable lower bound on the location of the critical endpoint

We perform the first direct determination of the position of the leading singularity of the pressure in the complex chemical potential $μ_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 $μ_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 $μ_B/T \approx 2$ and roughly temperature independent.

preprint2019arXiv

Magnetic catalysis and inverse catalysis for heavy pions

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.