Source author record

Biagio Lucini

Biagio Lucini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

62 published item(s)

preprint2023arXiv

Update on $SU(2)$ with one adjoint Dirac flavor

We present an update of our ongoing study of the SU(2) gauge theory with one flavor of Dirac fermion in the adjoint representation. Compared to our previous results we now have data at larger lattice volumes, smaller values of the fermion mass, and also larger values of $β$. We present data for the spectrum of mesons, baryons, glueballs, and the hybrid fermion-glue state, as well as new estimates of the mass anomalous dimension from both finite-size hyperscaling and the Dirac mode number, and discuss the implications of these data for the presence or otherwise of chiral symmetry breaking in this theory.

preprint2022arXiv

Applications of Machine Learning to Lattice Quantum Field Theory

There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment of this approach in the future.

preprint2022arXiv

Composite dynamics in Sp($2N$) gauge theories

Sp($2N$) gauge theories with fermonic matter provide an ideal laboratory to build extensions of the standard model based on novel composite dynamics. Examples include composite Higgs along with top partial compositeness and composite dark matter. Without fermions, their study also complements those based on SU($N_c$) gauge theories with which they share a common sector in the large $N_c=2N$ limit. We report on our recent progress in the numerical studies of Sp($2N$) gauge theories discretised on a four-dimensional Euclidean lattice. In particular, we present preliminary results for the low-lying mass spectra of mesons and chimera baryons in the theories with $N=2$. We also compute the topological susceptibility for various values of $N$, extrapolate the results to the large $N$ limit, and discuss certain universal properties in Yang-Mills theories.

preprint2022arXiv

Inverse Renormalization Group in Quantum Field Theory

We propose inverse renormalization group transformations within the context of quantum field theory that produce the appropriate critical fixed point structure, give rise to inverse flows in parameter space, and evade the critical slowing down effect in calculations pertinent to criticality. Given configurations of the two-dimensional $ϕ^{4}$ scalar field theory on sizes as small as $V=8^{2}$, we apply the inverse transformations to produce rescaled systems of size up to $V'=512^{2}$ which we utilize to extract two critical exponents. We conclude by discussing how the approach is generally applicable to any method that successfully produces configurations from a statistical ensemble and how it can give novel insights into the structure of the renormalization group.

preprint2022arXiv

Investigating the conformal behaviour of SU(2) with one adjoint Dirac flavor

We present a major update on our investigations of SU(2) gauge theory with one Dirac flavor in the adjoint representation on the lattice. In particular we consider larger volumes, as well as four different values of the gauge coupling. We provide results for the spectrum including gluonic, fermionic, and hybrid observables, Polyakov loops, and the anomalous dimension of the fermionic condensate from the Dirac mode number. These data confirm that the theory is close to the lower boundary of the conformal window for adjoint fermions. Our investigations provide important insights regarding the realization of different infrared scenarios that have been conjectured for this theory.

preprint2022arXiv

Lattice studies of the $Sp(4)$ gauge theory with two fundamental and three antisymmetric Dirac fermions

We consider the $Sp(4)$ gauge theory coupled to $N_f=2$ fundamental and $n_f=3$ antisymmetric flavours of Dirac fermions in four dimensions. This theory serves as the microscopic origin for composite Higgs models with $SU(4)/Sp(4)$ coset, supplemented by partial top compositeness. We study numerically its lattice realisation, and couple the fundamental plaquette action to Wilson-Dirac fermions in mixed representations, by adopting a (rational) hybrid Monte Carlo method, to perform non-trivial tests of the properties of the resulting lattice theory. We find evidence of a surface (with boundaries) of first-order bulk phase transitions in the three-dimensional space of bare parameters (one coupling and two masses). Explicit evaluation of the Dirac eigenvalues confirms the expected patterns of global symmetry breaking. After investigating finite volume effects in the weak-coupling phase of the theory, for the largest available lattice we study the mass spectra of the lightest spin-0 and spin-1 flavoured mesons composed of fermions in each representation, and of the lightest half-integer spin composite particle made of fermions in different representations -- the chimera baryon. This work sets the stage for future systematical studies of the non-perturbative dynamics in phenomenologically relevant regions of parameter space.

preprint2022arXiv

Quantitative analysis of phase transitions in two-dimensional XY models using persistent homology

We use persistent homology and persistence images as an observable of three different variants of the two-dimensional XY model in order to identify and study their phase transitions. We examine models with the classical XY action, a topological lattice action, and an action with an additional nematic term. In particular, we introduce a new way of computing the persistent homology of lattice spin model configurations and, by considering the fluctuations in the output of logistic regression and k-nearest neighbours models trained on persistence images, we develop a methodology to extract estimates of the critical temperature and the critical exponent of the correlation length. We put particular emphasis on finite-size scaling behaviour and producing estimates with quantifiable error. For each model we successfully identify its phase transition(s) and are able to get an accurate determination of the critical temperatures and critical exponents of the correlation length.

preprint2022arXiv

Quantum field theories, Markov random fields and machine learning

The transition to Euclidean space and the discretization of quantum field theories on spatial or space-time lattices opens up the opportunity to investigate probabilistic machine learning within quantum field theory. Here, we will discuss how discretized Euclidean field theories, such as the $ϕ^{4}$ lattice field theory on a square lattice, are mathematically equivalent to Markov fields, a notable class of probabilistic graphical models with applications in a variety of research areas, including machine learning. The results are established based on the Hammersley-Clifford theorem. We will then derive neural networks from quantum field theories and discuss applications pertinent to the minimization of the Kullback-Leibler divergence for the probability distribution of the $ϕ^{4}$ machine learning algorithms and other probability distributions.

preprint2022arXiv

The density of states method in Yang-Mills theories and first order phase transitions

Extensions of the standard model that lead to first-order phase transitions in the early universe can produce a stochastic background of gravitational waves, which may be accessible to future detectors. Thermodynamic observables at the transition, such as the latent heat, can be determined by lattice simulations, and then used to predict the expected signatures in a given theory. In lattice calculations, the emergence of metastabilities in proximity of the phase transition may make the precise determination of these observables quite challenging, and may lead to large uncontrolled numerical errors. In this contribution, we discuss as a prototype lattice calculation the first order deconfinement transition that arises in the strong SU(3) Yang-Mills sector. We adopt the novel logarithmic linear relaxation method, which can provide a determination of the density of states of the system with exponential error suppression. Thermodynamic observables can be reconstructed with a controlled error, providing a promising direction for accurate model predictions in the future.

preprint2022arXiv

Towards glueball masses of large-$N$ $\mathrm{SU}(N)$ pure-gauge theories without topological freezing

In commonly used Monte Carlo algorithms for lattice gauge theories the integrated autocorrelation time of the topological charge is known to be exponentially-growing as the continuum limit is approached. This $\mathit{topological}\,\,\textit{freezing}$, whose severity increases with the size of the gauge group, can result in potentially large systematics. To provide a direct quantification of the latter, we focus on $\mathrm{SU}(6)$ Yang--Mills theory at a lattice spacing for which conventional methods associated to the decorrelation of the topological charge have an unbearable computational cost. We adopt the recently proposed $\mathit{parallel}\,\,\mathit{tempering}\,\,\mathit{on}\,\,\mathit{boundary}\,\,\mathit{conditions}$ algorithm, which has been shown to remove systematic effects related to topological freezing, and compute glueball masses with a typical accuracy of $2-5\%$. We observe no sizeable systematic effect in the mass of the first lowest-lying glueball states, with respect to calculations performed at nearly-frozen topological sector.

preprint2021arXiv

Adding machine learning within Hamiltonians: Renormalization group transformations, symmetry breaking and restoration

We present a physical interpretation of machine learning functions, opening up the possibility to control properties of statistical systems via the inclusion of these functions in Hamiltonians. In particular, we include the predictive function of a neural network, designed for phase classification, as a conjugate variable coupled to an external field within the Hamiltonian of a system. Results in the two-dimensional Ising model evidence that the field can induce an order-disorder phase transition by breaking or restoring the symmetry, in contrast with the field of the conventional order parameter which causes explicit symmetry breaking. The critical behavior is then studied by proposing a Hamiltonian-agnostic reweighting approach and forming a renormalization group mapping on quantities derived from the neural network. Accurate estimates of the critical point and of the critical exponents related to the operators that govern the divergence of the correlation length are provided. We conclude by discussing how the method provides an essential step toward bridging machine learning and physics.

preprint2021arXiv

Efficient computations of continuous action densities of states for lattice models

The Logarithmic Linear Relaxation (LLR) algorithm is an efficient method for computing densities of states for systems with a continuous spectrum. A key feature of this method is exponential error reduction, which allows us to evaluate the density of states of a system over hundreds of thousands of orders of magnitude with a fixed level of relative accuracy. As a consequence of exponential error reduction, the LLR method provides a robust alternative to traditional Monte Carlo calculations in cases in which states suppressed by the Boltzmann weight play nevertheless a relevant role, e.g., as transition regions between dominant configuration sets. After reviewing the algorithm, we will show an application in U(1) Lattice Gauge Theory that has enabled us to obtain the most accurate estimate of the critical coupling with modest computational resources, defeating exponential tunneling times between metastable vacua. As a further showcase, we will then present an application of the LLR method to the decorrelation of the topological charge in SU(3) Lattice Gauge Theory near the continuum limit. Finally, we will review in general applications of the LLR algorithm to systems affected by a strong sign problem and discuss the case of the Bose gas at finite chemical potential.

preprint2021arXiv

Machine learning with quantum field theories

The precise equivalence between discretized Euclidean field theories and a certain class of probabilistic graphical models, namely the mathematical framework of Markov random fields, opens up the opportunity to investigate machine learning from the perspective of quantum field theory. In this contribution we will demonstrate, through the Hammersley-Clifford theorem, that the $ϕ^{4}$ scalar field theory on a square lattice satisfies the local Markov property and can therefore be recast as a Markov random field. We will then derive from the $ϕ^{4}$ theory machine learning algorithms and neural networks which can be viewed as generalizations of conventional neural network architectures. Finally, we will conclude by presenting applications based on the minimization of an asymmetric distance between the probability distribution of the $ϕ^{4}$ machine learning algorithms and target probability distributions.

preprint2021arXiv

Sp(4) gauge theories and beyond the standard model physics

We review numerical results for models with gauge group Sp(2N), discussing the glueball spectrum in the large-N limit, the quenched meson spectrum of Sp(4) with Dirac fermions in the fundamental and in the antisymmetric representation and the Sp(4) gauge model with two dynamical Dirac flavours. We also present preliminary results for the meson spectrum in the Sp(4) gauge theory with two fundamental and three antisymmetric Dirac flavours. The main motivation of our programme is to test whether this latter model is viable as a realisation of Higgs compositeness via the pseudo Nambu Goldstone mechanism and at the same time can provide partial top compositeness. In this respect, we report and briefly discuss preliminary results for the mass of the composite baryon made with two fundamental and one antisymmetric fermion (chimera baryon), whose physical properties are highly constrained if partial top compositeness is at work. Our investigation shows that a fully non-perturbative study of Higgs compositeness and partial top compositeness in Sp(4) is within reach with our current lattice methodology.

preprint2020arXiv

Color dependence of tensor and scalar glueball masses in Yang-Mills theories

We report the masses of the lightest spin-0 and spin-2 glueballs obtained in an extensive lattice study of the continuum and infinite volume limits of $Sp(N_c)$ gauge theories for $N_c=2,4,6,8$. We also extrapolate the combined results towards the large-$N_c$ limit. We compute the ratio of scalar and tensor masses, and observe evidence that this ratio is independent of $N_{c}$. Other lattice studies of Yang-Mills theories at the same space-time dimension provide a compatible ratio. We further compare these results to various analytical ones and discuss them in view of symmetry-based arguments related to the breaking of scale invariance in the underlying dynamics, showing that a constant ratio might emerge in a scenario in which the $0^{++}$ glueball is interpreted as a dilaton state.

preprint2020arXiv

Sp(4) gauge theories on the lattice: Quenched fundamental and antisymmetric fermions

We perform lattice studies of meson mass spectra and decay constants of the $Sp(4)$ gauge theory in the quenched approximation. We consider two species of (Dirac) fermions as matter field content, transforming in the 2-index antisymmetric and the fundamental representation of the gauge group, respectively. All matter fields are formulated as Wilson fermions. We extrapolate to the continuum and massless limits, and compare to each other the results obtained for the two species of mesons. In the case of two fundamental and three antisymmetric fermions, the long-distance dynamics is relevant for composite Higgs models. This is the first lattice study of this class of theories. The global $SU(4) \times SU(6)$ symmetry is broken to the $Sp(4) \times SO(6)$ subgroup, and the condensates align with the explicit mass terms present in the lattice formulation of the theory. The main results of our quenched calculations are that, with fermions in the 2-index antisymmetric representation of the group, the masses squared and decay constant squared of all the mesons we considered are larger than the corresponding quantities for the fundamental representation, by factors that vary between $\sim$1.2 and $\sim$2.7. We also present technical results that will be useful for future lattice investigations of dynamical simulations, of composite chimera baryons, and of the approach to large $N$ in the $Sp(2N)$ theories considered. We briefly discuss their high-temperature behaviour, where symmetry restoration and enhancement are expected.

preprint2019arXiv

Free energy of the self-interacting relativistic lattice Bose gas at finite density

The density of state approach has recently been proposed as a potential route to circumvent the sign problem in systems at finite density. In this study, using the Linear Logarithmic Relaxation (LLR) algorithm, we extract the generalised density of states, which is defined in terms of the imaginary part of the action, for the self-interacting relativistic lattice Bose gas at finite density. After discussing the implementation and testing the reliability of our approach, we focus on the determination of the free energy difference between the full system and its phase-quenched counterpart. Using a set of lattices ranging from $4^4$ to $16^4$ , we show that in the low density phase, this overlap free energy can be reliably extrapolated to the thermodynamic limit. The numerical precision we obtain with the LLR method allows us to determine with sufficient accuracy the expectation value of the phase factor, which is used in the calculation of the overlap free energy, down to values of ${\cal O}(10^{-480})$. When phase factor measurements are extended to the dense phase, a change of behaviour of the overlap free energy is clearly visible as the chemical potential crosses a critical value. Using fits inspired by the approximate validity of mean-field theory, which is confirmed by our simulations, we extract the critical chemical potential as the non-analyticity point in the overlap free energy, obtaining a value that is in agreement with other determinations. Implications of our findings and potential improvements of our methodology are also discussed.

preprint2019arXiv

Sp(4) gauge theories on the lattice: $N_f=2$ dynamical fundamental fermions

We perform lattice studies of the gauge theory with Sp(4) gauge group and two flavours of (Dirac) fundamental matter. The global SU(4) symmetry is spontaneously broken by the fermion condensate. The dynamical Wilson fermions in the lattice action introduce a mass that breaks the global symmetry also explicitly. The resulting pseudo-Nambu-Goldstone bosons describe the SU(4)/Sp(4) coset, and are relevant, in the context of physics beyond the Standard Model, for composite Higgs models. We discuss scale setting, continuum extrapolation and finite volume effects in the lattice theory. We study mesonic composite states, which span representations of the unbroken Sp(4) global symmetry, and we measure masses and decay constants of the (flavoured) spin-0 and spin-1 states accessible to the numerical treatment, as a function of the fermion mass. With help from the effective field theory treatment of such mesons, we perform a first extrapolation towards the massless limit. We assess our results by critically comparing to the literature on other models and to the quenched results, and we conclude by outlining future avenues for further exploration. The results of our spectroscopic analysis provide new input data for future phenomenological studies in the contexts of composite Higgs models, and of dark matter models with a strongly coupled dynamical origin.

preprint2018arXiv

Sp(4) gauge theory on the lattice: towards SU(4)/Sp(4) composite Higgs (and beyond)

The Sp(4) gauge theory with two Dirac fundamental flavours provides a candidate for the microscopic origin of composite-Higgs models based on the SU(4)/Sp(4) coset. We employ a combination of two different, complementary strategies for the numerical lattice calculations, based on the Hybrid Monte Carlo and on the Heat Bath algorithms. We perform pure Yang-Mills, quenched computations and exploratory studies with dynamical Wilson fermions. We present the first results in the literature for the spectrum of glueballs of the pure Sp(4) Yang-Mills theory, an EFT framework for the interpretation of the masses and decay constants of the lightest pion, vector and axial-vector mesons, and a preliminary calculation of the latter in the quenched approximation. We show the first numerical evidence of a bulk phase transition in the lattice theory with dynamical Wilson fermions, and perform the technical steps necessary to set up future investigations of the mesonic spectrum of the full theory.

preprint2016arXiv

BSMBench: a flexible and scalable supercomputer benchmark from computational particle physics

Lattice Quantum ChromoDynamics (QCD), and by extension its parent field, Lattice Gauge Theory (LGT), make up a significant fraction of supercomputing cycles worldwide. As such, it would be irresponsible not to evaluate machines' suitability for such applications. To this end, a benchmark has been developed to assess the performance of LGT applications on modern HPC platforms. Distinct from previous QCD-based benchmarks, this allows probing the behaviour of a variety of theories, which allows varying the ratio of demands between on-node computations and inter-node communications. The results of testing this benchmark on various recent HPC platforms are presented, and directions for future development are discussed.

preprint2016arXiv

Form the density-of-states method to finite density quantum field theory

During the last 40 years, Monte Carlo calculations based upon Importance Sampling have matured into the most widely employed method for determinig first principle results in QCD. Nevertheless, Importance Sampling leads to spectacular failures in situations in which certain rare configurations play a non-secondary role as it is the case for Yang-Mills theories near a first order phase transition or quantum field theories at finite matter density when studied with the re-weighting method. The density-of-states method in its LLR formulation has the potential to solve such overlap or sign problems by means of an exponential error suppression. We here introduce the LLR approach and its generalisation to complex action systems. Applications include U(1), SU(2) and SU(3) gauge theories as well as the Z3 spin model at finite densities and heavy-dense QCD.

preprint2016arXiv

Large mass hierarchies from strongly-coupled dynamics

Besides the Higgs particle discovered in 2012, with mass 125 GeV, recent LHC data show tentative signals for new resonances in diboson as well as diphoton searches at high center-of-mass energies (2 TeV and 750 GeV, respectively). If these signals are confirmed (or other new resonances are discovered at the TeV scale), the large hierarchies between masses of new bosons require a dynamical explanation. Motivated by these tentative signals of new physics, we investigate the theoretical possibility that large hierarchies in the masses of glueballs could arise dynamically in new strongly-coupled gauge theories extending the standard model of particle physics. We study lattice data on non-Abelian gauge theories in the (near-)conformal regime as well as a simple toy model in the context of gauge/gravity dualities. We focus our attention on the ratio $R$ between the mass of the lightest spin-2 and spin-0 resonances, that for technical reasons is a particularly convenient and clean observable to study. For models in which (non-perturbative) large anomalous dimensions arise dynamically, we show indications that this mass ratio can be large, with $R>5$. Moreover, our results suggest that $R$ might be related to universal properties of the IR fixed point. Our findings provide an interesting step towards understanding large mass ratios in the non-perturbative regime of quantum field theories with (near) IR conformal behaviour.

preprint2016arXiv

Overcoming strong metastabilities with the LLR method

In previous work, it has been shown that the recently proposed LLR method is very efficient at overcoming strong metastabilities that arise near first-order phase transition points. Here we present a systematic study of the performance of the algorithm near (pseudo-)critical points for $q$-state Potts models with $q$ as large as 20, in two and three dimensions. In particular, we shall focus our study on the ergodicity of the replica exchange step and the underlying physical mechanism. When compared with both analytical and numerical results present in the literature, our determinations of thermodynamic observables (including the order-disorder interface tension at criticality) show an impressive degree of relative accuracy (up to $2.5 \times 10^{-6}$), which confirms the reliability and the efficiency of the proposed approach.

preprint2015arXiv

An efficient algorithm for numerical computations of continuous densities of states

In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in Lattice Gauge Theories and in some Statistical Mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly-constant {\it relative} precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using Compact U(1) Lattice Gauge Theory. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from $8^4$ to $20^4$. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to $18^4$. Robust results for the $20^4$ volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. Other situations where the method is expected to be superior to importance sampling techniques are pointed out.

preprint2015arXiv

Charged hadrons in local finite-volume QED+QCD with C* boundary conditions

In order to calculate QED corrections to hadronic physical quantities by means of lattice simulations, a coherent description of electrically-charged states in finite volume is needed. In the usual periodic setup, Gauss's law and large gauge transformations forbid the propagation of electrically-charged states. A possible solution to this problem, which does not violate the axioms of local quantum field theory, has been proposed by Wiese and Polley, and is based on the use of C* boundary conditions. We present a thorough analysis of the properties and symmetries of QED in isolation and QED coupled to QCD, with C* boundary conditions. In particular we learn that a certain class of electrically-charged states can be constructed in this setup in a fully consistent fashion, without relying on gauge fixing. We argue that this class of states covers most of the interesting phenomenological applications in the framework of numerical simulations. We also calculate finite-volume corrections to the mass of stable charged particles and show that these are much smaller than in non-local formulations of QED.

preprint2015arXiv

Investigating some technical improvements to glueball calculations

We briefly discuss some issues concerning including the pseudoscalar glueball interpolating operators into the variational basis for computing the masses of the eta and etaprime mesons. As a start we present some preliminary results for the correlators from pseudoscalar glueball operators on nf = 2+1+1 twisted mass gauge configurations. Preliminary results for the effect of open boundary conditions on the masses of three glueballs computed in quenched QCD are presented. The statistics of glueball correlators are briefly discussed.

preprint2015arXiv

Numerical results for gauge theories near the conformal window

A novel strong interaction beyond the standard model could provide a dynamical explanation of electroweak symmetry breaking. Experimental results strongly constrain properties of models that realise this mechanism. Whether these constraints are obeyed by any strongly interacting quantum field theory is a non-perturbative problem that needs to be addressed by first-principle calculations. Monte Carlo simulations of lattice regularised gauge theories is a powerful tool that enables us to address this question. Recently various lattice investigations have appeared that have studied candidate models of strongly interacting dynamics beyond the standard model. After a brief review of the main methods and of some recent results, we focus on the analysis of SU(2) gauge theory with one adjoint Dirac fermion flavour, which is shown to have a near-conformal behaviour with an anomalous dimension of order one. The implications of our findings are also discussed.

preprint2015arXiv

Recent results from SU(2) with one adjoint Dirac fermion

We present some results for SU(2) with one adjoint Dirac flavour from lattice studies. Data for the spectroscopy, the static potential, topological charge, and the anomalous dimension of the fermionic condensate are included. Our findings are found to be in- consistent with conventional confining behaviour, instead pointing tentatively towards a theory lying within or very near the onset of the conformal window, with an anomalous dimension of the fermionic condensate of almost 1. Implications of these findings on the building of models of strongly-interacting dynamics beyond the standard model are discussed.

preprint2015arXiv

Topology and glueballs in SU(7) Yang-Mills with open boundary conditions

It is well known that the topology of gauge configurations generated in a Markov Monte-Carlo chain freezes as the continuum limit is approached. The corresponding autocorrelation time increases exponentially with the inverse lattice spacing, affecting the ergodicity of the simulation. In SU(N) gauge theories for large N this problem sets in at much coarser lattice spacings than for N=3. This means that its systematics can be studied on lattices that are smaller in terms of the number of lattice sites. It has been shown that using open boundary conditions in time allows instantons to be created and destroyed, restoring topological mobility and ergodicity. However, with open boundary conditions translational invariance is lost and the influence of spurious states propagating from the boundary into the bulk on physical correlators needs to be carefully evaluated. Moreover, while the total topological charge can be changed, the mobility of instantons across the lattice is still reduced. We consider SU(7) Yang-Mills theory and analyse its topological content in the periodic and open boundary condition cases. We also investigate scalar and pseudo-scalar glueball correlation functions.

preprint2014arXiv

A novel density of state method for complex action systems

Recently, a new and efficient algorithm (the LLR method) has been proposed for computing densities of states in statistical systems and gauge theories. In this talk, we explore whether this novel density of states method can be applied to numerical computations of observables in systems for which the action is complex. To this purpose, we introduce a generalised density of states, in terms of which integrals of oscillating observables can be determined semi-analytically, and we define a strategy to compute it with the LLR method. As a case study, we apply these ideas to the Z(3) spin model at finite density, finding a remarkable agreement of our results for the phase twist with those obtained with the worm algorithm for all explored chemical potentials, including values for which there are cancellations over sixteen orders of magnitude. These findings open new perspectives for dealing with the sign problem on physically more relevant systems.

preprint2014arXiv

Glueballs from the Lattice

Recent numerical calculations of the glueball spectrum in QCD, in SU($N$) Yang-Mills theory in the large-$N$ limit and in candidate theories of strongly interacting dynamics beyond the standard model (in which the lowest-lying scalar plays the role of the Higgs boson) are reviewed and their implications for our theoretical understanding of glueballs in QCD-like theories and in strongly coupled gauge theories with a (near-)conformal dynamics are discussed.

preprint2014arXiv

The density of states approach to dense quantum systems

We develop a first-principle generalised density of state method for studying numerically quantum field theories with a complex action. As a proof of concept, we show that with our approach we can solve numerically the strong sign problem of the $Z_3$ spin model at finite density. Our results are confirmed by standard simulations of the theory dual to the considered model, which is free from a sign problem. Our method opens new perspectives on ab initio simulations of cold dense quantum systems, and in particular of Yang-Mills theories with matter at finite densities, for which Monte Carlo based importance sampling are unable to produce sufficiently accurate results.

preprint2014arXiv

The infrared regime of SU(2) with one adjoint Dirac flavour

SU(2) gauge theory with one Dirac flavour in the adjoint representation is investigated on a lattice. Initial results for the gluonic and mesonic spectrum, static potential from Wilson and Polyakov loops, and the anomalous dimension of the fermionic condensate from the Dirac mode number are presented. The results found are not consistent with conventional confining behaviour, instead tentatively pointing towards a theory lying within or very near the onset of the conformal window, with the anomalous dimension of the fermionic condensate in the range $0.9 \lesssim γ_* \lesssim 0.95$. The implications of our work for building a viable theory of strongly interacting dynamics beyond the standard model are discussed.

preprint2013arXiv

First results for SU(2) Yang-Mills with one adjoint Dirac Fermion

We present a first exploratory study of SU(2) gauge theory with one Dirac flavour in the adjoint representation. We provide initial results for the spectroscopy and the anomalous dimension for the chiral condensate. Our investigation indicates that the theory is conformal or near-conformal, with an anomalous dimension of order one. A discussion of the relevance of these findings in relation to walking technicolor scenarios is also presented.

preprint2013arXiv

Infrared conformality and bulk critical points: SU(2) with heavy adjoint quarks

The lattice phase structure of a gauge theory can be a serious obstruction to Monte Carlo studies of its continuum behaviour. This issue is particularly delicate when numerical studies are performed to determine whether a theory is in a (near-)conformal phase. In this work we investigate the heavy mass limit of the SU(2) gauge theory with Nf=2 adjoint fermions and its lattice phase diagram, showing the presence of a critical point ending a line of first order bulk phase transition. The relevant gauge observables and the low-lying spectrum are monitored in the vicinity of the critical point with very good control over different systematic effects. The scaling properties of masses and susceptibilities open the possibility that the effective theory at criticality is a scalar theory in the universality class of the four-dimensional Gaussian model. This behaviour is clearly different from what is observed for SU(2) gauge theory with two dynamical adjoint fermions, whose (near-)conformal numerical signature is hence free from strong-coupling bulk effects.

preprint2013arXiv

Large-N mesons

We present an update of our project of computing the meson spectrum and decay constants in large-N QCD. The results are obtained in the quenched approximation with the Wilson fermion action for N = 2, 3, 4, 5, 6, 7 and 17 and extrapolated to infinite N. We non-perturbatively determine the renormalization factors for local quark bilinears that are needed to compute the decay constants. We extrapolate our SU(7) results to the continuum limit, employing four different lattice spacings.

preprint2013arXiv

Mesons in large-N QCD

We present the results of a systematic, first-principles study of the spectrum and decay constants of mesons for different numbers of color charges N, via lattice computations. We restrict our attention to states in the non-zero isospin sector, evaluating the masses associated with the ground-state and first excitation in the pseudoscalar, vector, scalar, and axial vector channels. Our results are based on a new set of simulations of four dimensional SU(N) Yang-Mills theories with the number of colors ranging from N=2 to N=17; the spectra and the decay constants are computed in the quenched approximation (which becomes exact in the 't Hooft limit) using Wilson fermions. After discussing the extrapolations to the chiral and large-N limits, we present a comparison of our results to some of the numerical computations and analytical predictions available in the literature - including, in particular, those from holographic computations.

preprint2013arXiv

SU(N) gauge theories at large N

We review the theoretical developments and conceptual advances that stemmed from the generalization of QCD to the limit of a large number of color charges, originally proposed by 't Hooft. Then, after introducing the gauge-invariant non-perturbative formulation of non-Abelian gauge theories on a spacetime lattice, we present a selection of results from recent lattice studies of theories with a different number of colors, and the findings obtained from their extrapolation to the 't Hooft limit. We conclude with a brief discussion and a summary.

preprint2013arXiv

Topology of Minimal Walking Technicolor

We perform a lattice study of the topological susceptibility and instanton size distribution of the $\su{2}$ gauge theory with two adjoint Dirac fermions (also known as Minimal Walking Technicolor), which is known to be in the conformal window. In the theory deformed with a small mass term, by drawing a comparison with the pure gauge theory, we find that topological observables are decoupled from the fermion dynamics. This provides further evidence for the infrared conformality of the theory. A study of the instanton size distribution shows that this quantity can be used to detect the onset of finite size effects.

preprint2013arXiv

Towards a density of states approach for dense matter systems

The density-of-states method (Phys.Rev.Lett. 109 (2012) 111601) features an exponential error suppression and is not restricted to theories with positive probabilistic weight. It is applied to the SU(2) gauge theory at finite densities of heavy quarks. The key ingredient here is the Polyakov line probability distribution, which is obtained of over 80 orders of magnitude. We briefly address whether the exponential error suppression could be sufficient to simulate theories with a strong sign problem.

preprint2012arXiv

Scaling properties of SU(2) gauge theory with mixed fundamental-adjoint action

We study the phase diagram of the SU(2) lattice gauge theory with fundamental-adjoint Wilson plaquette action. We confirm the presence of a first order bulk phase transition and we estimate the location of its end-point in the bare parameter space. If this point is second order, the theory is one of the simplest realizations of a lattice gauge theory admitting a continuum limit at finite bare couplings. All the relevant gauge observables are monitored in the vicinity of the fixed point with very good control over finite-size effects. The scaling properties of the low-lying glueball spectrum are studied while approaching the end-point in a controlled manner.

preprint2012arXiv

SU(N_c) gauge theories at deconfinement

The deconfinement transition in SU($N_c$) Yang--Mills is investigated by Monte Carlo simulations of the gauge theory discretized on a spacetime lattice. We present new results for $ 4 \le N_c \le 8$ (in particular, for $N_c = 5$ and $N_c = 7$), which are analysed together with previously published results. The increased amount of data, the improved statistics and simulations closer to the continuum limit provide us with better control over systematic errors. After performing the thermodynamic limit, numerical results for the ratio of the critical temperature $T_c$ over the square root of the string tension $\sqrtσ$ obtained on lattices with temporal extensions $N_t = 5,6,7,8$ are extrapolated to the continuum limit. The continuum results at fixed $N_c$ are then extrapolated to $N_c = \infty$. We find that our data are accurately described by the formula $T_c/\sqrtσ = 0.5949(17) + 0.458(18)/N_c^2$. Possible systematic errors affecting our calculations are also discussed.

preprint2012arXiv

The density of states in gauge theories

The density of states is calculated for a SU(2) and a compact U(1) lattice gauge theory using a modified version of the Wang-Landau algorithm. We find that the density of states of the SU(2) gauge theory can be reliably calculated over a range of 120,000 orders of magnitude for lattice sizes as big as 20^4. We demonstrate the potential of the algorithm by reproducing the SU(2) average action, its specific heat and the critical couplings of the weak first order transition in U(1).

preprint2011arXiv

Finite volume effects in SU(2) with two adjoint fermions

Many evidences from lattice simulations support the idea that SU(2) with two Dirac flavors in the adjoint representation (also called Minimal Walking Technicolor) is IR conformal. A possible way to see this is through the behavior of the spectrum of the mass-deformed theory. When fermions are massive, a mass-gap is generated and the theory is confined. IR-conformality is recovered in the chiral limit: masses of particles vanish in the chiral limit, while their ratios stay finite. In order to trust this analysis one has to relay on the infinite volume extrapolation. We will discuss the finite volume effects on the mesonic spectrum, investigated by varying the size of the lattice and by changing the boundary conditions for the fields.

preprint2011arXiv

Glueball masses in the large N limit

The lowest-lying glueball masses are computed in SU($N$) gauge theory on a spacetime lattice for constant value of the lattice spacing $a$ and for $N$ ranging from 3 to 8. The lattice spacing is fixed using the deconfinement temperature at temporal extension of the lattice $N_T = 6$. The calculation is conducted employing in each channel a variational ansatz performed on a large basis of operators that includes also torelon and (for the lightest states) scattering trial functions. This basis is constructed using an automatic algorithm that allows us to build operators of any size and shape in any irreducible representation of the cubic group. A good signal is extracted for the ground state and the first excitation in several symmetry channels. It is shown that all the observed states are well described by their large $N$ values, with modest ${\cal O}(1/N^2)$ corrections. In addition spurious states are identified that couple to torelon and scattering operators. As a byproduct of our calculation, the critical couplings for the deconfinement phase transition for N=5 and N=7 and temporal extension of the lattice $N_T=6$ are determined.

preprint2011arXiv

Improved Lattice Spectroscopy of Minimal Walking Technicolor

We present a numerical study of spectroscopic observables in the SU(2) gauge theory with two adjoint fermions using improved source and sink operators. We compare in detail our improved results with previous determinations of masses that used point sources and sinks and we investigate possible systematic effects in both cases. Such comparison enables us to clearly assess the impact of a short temporal extent on the physical picture, and to investigate some effects due to the finite spatial box. While confirming the IR-conformal behaviour of the theory, our investigation shows that in order to make firm quantitative predictions, a better handle on finite size effects is needed.

preprint2011arXiv

k-String tensions and the 1/N expansion

We address the question of whether the large-N expansion in pure SU(N) gauge theories requires that k-string tensions must have a power series expansion in 1/N^2, as in the sine law, or whether 1/N contributions are also allowable, as in Casimir scaling. We find that k-string tensions may, in fact, have 1/N corrections, and consistency with the large-N expansion in the open-string sector depends crucially on an exact cancellation, which we will prove, among terms involving odd powers of 1/N in particular combinations of Wilson loops. It is shown how these cancellations are fulfilled, and consistency with the large-N expansion achieved, in a concrete example, namely, strong-coupling lattice gauge theory with the heat-kernel action. This is a model which has both a 1/N^2 expansion and Casimir scaling of the k-string tensions. Analysis of the closed string channel in this model confirms our conclusions, and provides further insights into the large-N dependence of energy eigenstates and eigenvalues.

preprint2011arXiv

Magnetic monopole plasma phase in (2+1)d compact quantum electrodynamics with fermionic matter

We present the first evidence from lattice simulations that the magnetic monopoles in three dimensional compact quantum electrodynamics (cQED3) with N_f=2 and N_f= 4 four-component fermion flavors are in a plasma phase. The evidence is based mainly on the divergence of the monopole susceptibility (polarizability) with the lattice size at weak gauge couplings. A weak four-Fermi term added to the cQED3 action enabled simulations with massless fermions. The exact chiral symmetry of the interaction terms forbids symmetry breaking lattice discretization counterterms to appear in the theory's effective action. It is also shown that the scenario of a monopole plasma does not depend on the strength of the four-Fermi coupling. Other observables such as the densities of "isolated" dipoles and monopoles and the so-called specific heat show that a crossover from a dense monopole plasma to a dilute monopole gas occurs at strong couplings. The implications of our results on the stability of U(1) spin liquids in two spatial dimensions are also discussed.

preprint2010arXiv

A numerical investigation of orientifold planar equivalence for quenched mesons

We measure on the lattice the quenched pseudoscalar and vector meson masses at a fixed value of the lattice spacing for SU(N) gauge theory with fermions in the adjoint, in the symmetric and in the antisymmetric representations of the gauge group. Simulations are performed for N=3,4,6 in all those representations, with the addition of N=2 for the adjoint representation. We illustrate a strategy for separating the even from the odd-power contributions in 1/N in the masses. Using this technique, we extrapolate the vector mass to the large-N limit in the chiral region and show that at N = infty this mass is the same within errors in all the three representations, as predicted by orientifold planar equivalence. Possible implications of our investigation for studying orientifold planar equivalence in the dynamical case are discussed.

preprint2010arXiv

Confining vs. conformal scenario for SU(2) with adjoint fermions. Gluonic observables

Walking technicolor is a mechanism for electroweak symmetry breaking without Higgs field. The Higgs mechanism is provided by chiral symmetry breaking in the technicolor theory. An essential ingredient is the vicinity to an IR fixed point, which could reconcile technicolor with the electroweak precision tests. SU(2) gauge theory with two Dirac adjoint fermions has been proposed as a candidate for walking technicolor. Understanding whether this theory is confining or IR-conformal is a challenging problem, which can be addressed by means of numerical simulations. We have pointed out that a clean signal for the existence of an IR fixed point in this theory can be obtained by comparing the mesonic and gluonic sectors. We review some technical details of our calculations. Possible systematic errors are discussed.

preprint2010arXiv

Improved Spectroscopy of Minimal Walking Technicolor

We present an improved study of spectroscopic observables in the $SU(2)$ Yang-Mills theory with two adjoint fermions. We make an improvement on the precision of previous results which clarify the scale of finite volume effects present. This analysis adds to the evidence for near-conformal dynamics of this theory, while indicating a preference for a low anomalous mass dimension of the massless theory.

preprint2010arXiv

Mesonic spectroscopy of Minimal Walking Technicolor

We investigate the structure and the novel emerging features of the mesonic non-singlet spectrum of the Minimal Walking Technicolor (MWT) theory. Precision measurements in the nonsinglet pseudoscalar and vector channels are compared to the expectations for an IR-conformal field theory and a QCD-like theory. Our results favor a scenario in which MWT is (almost) conformal in the infrared, while spontaneous chiral symmetry breaking seems less plausible.

preprint2010arXiv

Orientifold Planar Equivalence: The Quenched Meson Spectrum

A numerical study of Orientifold Planar Equivalence is performed in SU(N) Yang-Mills theories for N=2,3,4,6. Quenched meson masses are extracted in the antisymmetric, symmetric and adjoint representations for the pseudoscalar and vector channels. An extrapolation of the vector mass as a function of the pseudoscalar mass to the large-N limit shows that the numerical results agree within errors for the three theories, as predicted by Orientifold Planar Equivalence. As a byproduct of the extrapolation, the size of the corrections up to O(1/N^3) are evaluated. A crucial prerequisite for the extrapolation is the determination of an analytical relationship between the corrections in the symmetric and in the antisymmetric representations, order by order in a 1/N expansion.

preprint2010arXiv

Strongly Interacting Dynamics beyond the Standard Model on a Spacetime Lattice

Strong theoretical arguments suggest that the Higgs sector of the Standard Model of the Electroweak interactions is an effective low-energy theory, with a more fundamental theory that is expected to emerge at an energy scale of the order of the TeV. One possibility is that the more fundamental theory be strongly interacting and the Higgs sector be given by the low-energy dynamics of the underlying theory. We review recent works aimed to determining observable quantities by numerical simulations of strongly interacting theories proposed in the literature for explaining the Electroweak symmetry breaking mechanism. These investigations are based on Monte Carlo simulations of the theory formulated on a spacetime lattice. We focus on the so-called Minimal Walking Technicolour scenario, a SU(2) gauge theory with two flavours of fermions in the adjoint representation. The emerging picture is that this theory has an infrared fixed point that dominates the large distance physics. We shall discuss the first numerical determinations of quantities of phenomenological interest for this theory and analyse future directions of quantitative studies of strongly interacting beyond the Standard Model theories with Lattice techniques. In particular, we report on a finite size scaling determination of the chiral condensate anomalous dimension $γ$, for which we find $0.05 \le γ\le 0.25$.

preprint2010arXiv

The glueball spectrum at large N

The lowest-lying glueball masses are computed in SU($N$) gauge theory on a spacetime lattice for constant value of the lattice spacing $a$ and for $N$ ranging from 3 to 8. The lattice spacing is fixed using the deconfinement temperature at temporal extension of the lattice $N_T = 6$. The calculation is conducted employing in each channel a variational ansatz performed on a large basis of operators that includes also torelon and (for the lightest states) scattering trial functions. This basis is constructed using an automatic algorithm that allows us to build operators of any size and shape in any irreducible representation of the cubic group. A good signal is extracted for the ground state and the first excitation in several symmetry channels. It is shown that all the observed states are well described by their large $N$ values, with modest ${\cal O}(1/N^2)$ corrections. In addition spurious states are identified that couple to torelon and scattering operators.

preprint2010arXiv

The infrared dynamics of Minimal Walking Technicolor

We study the gauge sector of Minimal Walking Technicolor, which is an SU(2) gauge theory with nf=2 flavors of Wilson fermions in the adjoint representation. Numerical simulations are performed on lattices Nt x Ns^3, with Ns ranging from 8 to 16 and Nt=2Ns, at fixed β=2.25, and varying the fermion bare mass m0, so that our numerical results cover the full range of fermion masses from the quenched region to the chiral limit. We present results for the string tension and the glueball spectrum. A comparison of mesonic and gluonic observables leads to the conclusion that the infrared dynamics is given by an SU(2) pure Yang-Mills theory with a typical energy scale for the spectrum sliding to zero with the fermion mass. The typical mesonic mass scale is proportional to, and much larger than this gluonic scale. Our findings are compatible with a scenario in which the massless theory is conformal in the infrared. An analysis of the scaling of the string tension with the fermion mass towards the massless limit allows us to extract the chiral condensate anomalous dimension γ*, which is found to be γ*=0.22+-0.06.

preprint2008arXiv

Orientifold Planar Equivalence: The Chiral Condensate

The recently introduced orientifold planar equivalence is a promising tool for solving non-perturbative problems in QCD. One of the predictions of orientifold planar equivalence is that the chiral condensates of a theory with $N_f$ flavours of Dirac fermions in the symmetric (or antisymmetric) representation and $N_f$ flavours of Majorana fermions in the adjoint representation have the same large $N$ value for any value of the mass of the (degenerate) fermions. Assuming the invariance of the theory under charge conjugation, we prove this statement on the lattice for staggered quenched condensates in SU($N$) Yang-Mills in the large $N$ limit. Then, we compute numerically those quenched condensates for $N$ up to 8. After separating the even from the odd corrections in $1/N$, we are able to show that our data support the equivalence; however, unlike other quenched observables, subleading terms in $1/N$ are needed for describing the data for the symmetric and antisymmetric representation at $N$=3. Possible lessons for the unquenched case are discussed.

preprint2008arXiv

The running of the coupling in SU(N) pure gauge theories

The running of the coupling is studied in SU(4) gauge theory using the Schrödinger functional technique. Up to energies of the order of the square root of the string tension $σ$, the running is found to agree with the two-loop perturbative formula. Relating the perturbative to the non-perturbative regime of the running and converting to the $\msbar$ scheme allows one to extract the ratio $Λ_{\msbar}/\sqrtσ$. The result is then used in combination with similar calculations present in the literature for SU(2) and SU(3) to extract $Λ_{\msbar}/\sqrtσ$ in the large $N$ limit. Our results for $N=3,4$ agree with a recent study of the same quantity performed using the Parisi mean field scheme as an intermediate scheme, while $Λ_{\msbar}/\sqrtσ$ in SU(2) turns out to differ by 2.5%. Possible explanations of this discrepancy are discussed.