Source author record

Enrico Rinaldi

Enrico Rinaldi 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

53works
14topics
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

53 published item(s)

preprint2022arXiv

Confinement/deconfinement transition in the D0-brane matrix model -- A signature of M-theory?

We study the confinement/deconfinement transition in the D0-brane matrix model (often called the BFSS matrix model) and its one-parameter deformation (the BMN matrix model) numerically by lattice Monte Carlo simulations. Our results confirm general expectations from the dual string/M-theory picture for strong coupling. In particular, we observe the confined phase in the BFSS matrix model, which is a nontrivial consequence of the M-theory picture. We suggest that these models provide us with an ideal framework to study the Schwarzschild black hole, M-theory, and furthermore, the parameter region of the phase transition between type IIA superstring theory and M-theory. A detailed study of M-theory via lattice Monte Carlo simulations of the D0-brane matrix model might be doable with much smaller computational resources than previously expected.

preprint2022arXiv

Detailed analysis of excited state systematics in a lattice QCD calculation of $g_A$

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations $\gtrsim1$ fm are commonly used to combat these issues. In this work, focusing on $g_A$, we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2--1 fm), combined with a multi-state analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

preprint2022arXiv

On the evolution of the Hubble constant with the SNe Ia Pantheon Sample and Baryon Acoustic Oscillations: a feasibility study for GRB-cosmology in 2030

The difference from 4 to 6 $σ$ in the Hubble constant ($H_0$) between the values observed with the local (Cepheids and Supernovae Ia, SNe Ia) and the high-z probes (CMB obtained by the Planck data) still challenges the astrophysics and cosmology community. Previous analysis has shown that there is an evolution in the Hubble constant that scales as $f(z)=\mathcal{H}_{0}/(1+z)^η$, where $\mathcal{H}_0$ is $H_{0}(z=0)$ and $η$ is the evolutionary parameter. Here, we investigate if this evolution still holds by using the SNe Ia gathered in the Pantheon sample and the BAOs. We assume $H_{0}=70 \textrm{km s}^{-1}\textrm{Mpc}^{-1}$ as the local value and divide the Pantheon into 3 bins ordered in increasing values of redshift. Similar to our previous analysis but varying two cosmological parameters contemporaneously ($H_0$, $Ω_{0m}$ in the $Λ$CDM model and $H_0$, $w_a$ in the $w_{0}w_{a}$CDM model), for each bin we implement a MCMC analysis obtaining the value of $H_0$ $[...]$. Subsequently, the values of $H_0$ are fitted with the model $f(z)$. Our results show that a decreasing trend with $η\sim10^{-2}$ is still visible in this sample. The $η$ coefficient reaches zero in 2.0 $σ$ for the $Λ$CDM model up to 5.8 $σ$ for $w_{0}w_{a}$CDM model. This trend, if not due to statistical fluctuations, could be explained through a hidden astrophysical bias, such as the effect of stretch evolution, or it requires new theoretical models, a possible proposition is the modified gravity theories, $f(R)$ $[...]$. This work is also a preparatory to understand how the combined probes still show an evolution of the $H_0$ by redshift and what is the current status of simulations on GRB cosmology to obtain the uncertainties on the $Ω_{0m}$ comparable with the ones achieved through SNe Ia.

preprint2022arXiv

Optical and X-ray GRB Fundamental Planes as Cosmological Distance Indicators

Gamma-Ray Bursts (GRBs), can be employed as standardized candles, extending the distance ladder beyond Supernovae Type Ia (SNe Ia, $z=2.26$). We standardize GRBs using the 3D fundamental plane relation (the Dainotti relation) among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain $Ω_{\text{M}}= 0.299 \pm 0.009$ assuming a flat $Λ$CDM cosmology with and without correcting GRBs for selection biases and redshift evolution. Using a 3D optical Dainotti correlation, we find this sample is as efficacious in the determination of $Ω_{\text{M}}$ as the X-ray sample. We trimmed our GRB samples to achieve tighter planes to simulate additional GRBs. We determined how many GRBs are needed as standalone probes to achieve a comparable precision on $Ω_{\text{M}}$ to the one obtained by SNe Ia only. We reach the same error measurements derived using SNe Ia in 2011 and 2014 with 142 and 284 simulated optical GRBs, respectively, considering the errorbars on the variables halved. These error limits will be reached in 2038 and in 2047, respectively. Using a doubled sample (obtained by future machine learning approaches allowing a lightcurve reconstruction and the estimates of GRB redhifts when z is unknown) compared to the current sample, with errorbars halved we will reach the same precision as SNe Ia in 2011 and 2014, now and in 2026, respectively. If we consider the current SNe precision, this will be reached with 390 optical GRBs by 2054.

preprint2022arXiv

Running Hubble constant from the SNe Ia Pantheon sample?

The mismatch between different independent measurements of the expansion rate of the Universe is known as the Hubble constant ($H_0$) tension, and it is a serious and pressing problem in cosmology. We investigate this tension considering the dataset from the Pantheon sample, a collection of 1048 Type Ia Supernovae (SNe Ia) with a redshift range $0<z<2.26$. We perform a binned analysis in redshift to study if the $H_0$ tension also occurs in SNe Ia data. Hence, we build equally populated subsamples in three and four bins, and we estimate $H_{0}$ in each bin considering the $Λ$CDM and $w_{0}w_{a}$CDM cosmological models. We perform a statistical analysis via a Markov Chain Monte Carlo (MCMC) method for each bin. We observe that $H_0$ evolves with the redshift, using a fit function $H_{0}(z)=\tilde{H}_{0} (1+z)^{-α}$ with two fitting parameters $α$ and $\tilde{H}_{0}$. Our results show a decreasing behavior of $H_0$ with $α\sim 10^{-2}$ and a consistency with no evolution between 1.2 $σ$ and 2.0 $σ$. Considering the $H_0$ tension, we extrapolate $H_{0}(z)$ until the redshift of the last scattering surface, $z=1100$, obtaining values of $H_0$ consistent in 1 $σ$ with the cosmic microwave background (CMB) measurements by Planck. Finally, we discuss possible $f(R)$ modified gravity models to explain a running Hubble constant with the redshift, and we infer the form of the scalar field potential in the dynamically equivalent Jordan frame.

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

The hyperon spectrum from lattice QCD

Hyperon decays present a promising alternative for extracting $\vert V_{us} \vert$ from lattice QCD combined with experimental measurements. Currently $\vert V_{us} \vert$ is determined from the kaon decay widths and a lattice calculation of the associated form factor. In this proceeding, I will present preliminary work on a lattice determination of the hyperon mass spectrum. I will additionally summarize future goals in which we will calculate the hyperon transition matrix elements, which will provide an alternative means for accessing $\vert V_{us} \vert$. This work is based on a particular formulation of SU(2) chiral perturbation theory for hyperons; determining the extent to which this effective field theory converges is instrumental in understanding the limits of its predictive power, especially since some hyperonic observables are difficult to calculate near the physical pion mass (e.g., hyperon-to-nucleon form factors), and thus the use of heavier than physical pion masses is likely to yield more precise results when combined with extrapolations to the physical point.}

preprint2022arXiv

Using Multivariate Imputation by Chained Equations to Predict Redshifts of Active Galactic Nuclei

Redshift measurement of active galactic nuclei (AGNs) remains a time-consuming and challenging task, as it requires follow up spectroscopic observations and detailed analysis. Hence, there exists an urgent requirement for alternative redshift estimation techniques. The use of machine learning (ML) for this purpose has been growing over the last few years, primarily due to the availability of large-scale galactic surveys. However, due to observational errors, a significant fraction of these data sets often have missing entries, rendering that fraction unusable for ML regression applications. In this study, we demonstrate the performance of an imputation technique called Multivariate Imputation by Chained Equations (MICE), which rectifies the issue of missing data entries by imputing them using the available information in the catalog. We use the Fermi-LAT Fourth Data Release Catalog (4LAC) and impute 24% of the catalog. Subsequently, we follow the methodology described in Dainotti et al. (2021) and create an ML model for estimating the redshift of 4LAC AGNs. We present results which highlight positive impact of MICE imputation technique on the machine learning models performance and obtained redshift estimation accuracy.

preprint2021arXiv

Matrix Model simulations using Quantum Computing, Deep Learning, and Lattice Monte Carlo

Matrix quantum mechanics plays various important roles in theoretical physics, such as a holographic description of quantum black holes. Understanding quantum black holes and the role of entanglement in a holographic setup is of paramount importance for the development of better quantum algorithms (quantum error correction codes) and for the realization of a quantum theory of gravity. Quantum computing and deep learning offer us potentially useful approaches to study the dynamics of matrix quantum mechanics. In this paper we perform a systematic survey for quantum computing and deep learning approaches to matrix quantum mechanics, comparing them to Lattice Monte Carlo simulations. In particular, we test the performance of each method by calculating the low-energy spectrum.

preprint2021arXiv

Near-conformal dynamics in a chirally broken system

Composite Higgs models must exhibit very different dynamics from quantum chromodynamics (QCD) regardless whether they describe the Higgs boson as a dilatonlike state or a pseudo-Nambu-Goldstone boson. Large separation of scales and large anomalous dimensions are frequently desired by phenomenological models. Mass-split systems are well-suited for composite Higgs models because they are governed by a conformal fixed point in the ultraviolet but are chirally broken in the infrared. In this work we use lattice field theory calculations with domain wall fermions to investigate a system with four light and six heavy flavors. We demonstrate how a nearby conformal fixed point affects the properties of the four light flavors that exhibit chiral symmetry breaking in the infrared. Specifically we describe hyperscaling of dimensionful physical quantities and determine the corresponding anomalous mass dimension. We obtain $y_m=1+γ^*= 1.47(5)$ suggesting that $N_f=10$ lies inside the conformal window. Comparing the low energy spectrum to predictions of dilaton chiral perturbation theory, we observe excellent agreement which supports the expectation that the 4+6 mass-split system exhibits near-conformal dynamics with a relatively light $0^{++}$ isosinglet scalar.

preprint2021arXiv

Partial Deconfinement at Strong Coupling on the Lattice

We provide evidence for partial deconfinement -- the deconfinement of a SU($M$) subgroup of the SU($N$) gauge group -- by using lattice Monte Carlo simulations. We take matrix models as concrete examples. By appropriately fixing the gauge, we observe that the $M\times M$ submatrices deconfine. This gives direct evidence for partial deconfinement at strong coupling. We discuss the applications to QCD and holography.

preprint2021arXiv

Two-nucleon S-wave interactions at the $SU(3)$ flavor-symmetric point with $m_{ud}\simeq m_s^{\rm phys}$: a first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multi-particle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a positive definite set of two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance in order to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the $SU(3)$-flavor symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of $a\sim0.086$~fm, lattice volume of $L=48a\simeq4.1$~fm and pion mass $m_π\simeq714$ MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

preprint2020arXiv

$F_K / F_π$ from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of $F_K/F_π$ using Möbius domain-wall fermions computed on gradient-flowed $N_f=2+1+1$ highly-improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from $130 \lesssim m_π\lesssim 400$ MeV, four lattice spacings of $a\sim 0.15, 0.12, 0.09$ and $0.06$ fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the $a\sim0.06$ fm ensemble is helpful, but not necessary to achieve a subpercent determination of $F_K/F_π$. We also include an estimate of the strong isospin breaking corrections and arrive at a final result of $F_{K^\pm}/F_{π^\pm} = 1.1942(45)$ with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of $|V_{us}|/|V_{ud}| = 0.2311(10)$.

preprint2019arXiv

A critical point for bifurcation cascades and featureless turbulence

In this Letter we show that a bifurcation cascade and fully sustained turbulence can share the phase space of a fluid flow system, resulting in the presence of competing stable attractors. We analyse the toroidal pipe flow, which undergoes subcritical transition to turbulence at low pipe curvatures and supercritical transition at high curvatures, as was previously documented. We provide decisive evidence that the nature of the supercritical transition is of Ruelle--Takens type and that, in a narrow range of intermediate curvatures, its dynamics competes with that of sustained turbulence emerging through subcritical transition mechanisms.

preprint2019arXiv

Lattice QCD Determination of $g_A$

The nucleon axial coupling, $g_A$, is a fundamental property of protons and neutrons, dictating the strength with which the weak axial current of the Standard Model couples to nucleons, and hence, the lifetime of a free neutron. The prominence of $g_A$ in nuclear physics has made it a benchmark quantity with which to calibrate lattice QCD calculations of nucleon structure and more complex calculations of electroweak matrix elements in one and few nucleon systems. There were a number of significant challenges in determining $g_A$, notably the notorious exponentially-bad signal-to-noise problem and the requirement for hundreds of thousands of stochastic samples, that rendered this goal more difficult to obtain than originally thought. I will describe the use of an unconventional computation method, coupled with "ludicrously'" fast GPU code, access to publicly available lattice QCD configurations from MILC and access to leadership computing that have allowed these challenges to be overcome resulting in a determination of $g_A$ with 1% precision and all sources of systematic uncertainty controlled. I will discuss the implications of these results for the convergence of $SU(2)$ Chiral Perturbation theory for nucleons, as well as prospects for further improvements to $g_A$ (sub-percent precision, for which we have preliminary results) which is part of a more comprehensive application of lattice QCD to nuclear physics. This is particularly exciting in light of the new CORAL supercomputers coming online, Sierra and Summit, for which our lattice QCD codes achieve a machine-to-machine speed up over Titan of an order of magnitude.

preprint2019arXiv

Thermal phase transition in Yang-Mills matrix model

We study the bosonic matrix model obtained as the high-temperature limit of two-dimensional maximally supersymmetric SU($N$) Yang-Mills theory. So far, no consensus about the order of the deconfinement transition in this theory has been reached and this hinders progress in understanding the nature of the black hole/black string topology change from the gauge/gravity duality perspective. On the one hand, previous works considered the deconfinement transition consistent with two transitions which are of second and third order. On the other hand, evidence for a first order transition was put forward more recently. We perform high-statistics lattice Monte Carlo simulations at large $N$ and small lattice spacing to establish that the transition is really of first order. Our findings flag a warning that the required large-$N$ and continuum limit might not have been reached in earlier publications, and that was the source of the discrepancy. Moreover, our detailed results confirm the existence of a new partially deconfined phase which describes non-uniform black strings via the gauge/gravity duality. This phase exhibits universal features already predicted in quantum field theory.

preprint2018arXiv

Simulating the weak death of the neutron in a femtoscale universe with near-Exascale computing

The fundamental particle theory called Quantum Chromodynamics (QCD) dictates everything about protons and neutrons, from their intrinsic properties to interactions that bind them into atomic nuclei. Quantities that cannot be fully resolved through experiment, such as the neutron lifetime (whose precise value is important for the existence of light-atomic elements that make the sun shine and life possible), may be understood through numerical solutions to QCD. We directly solve QCD using Lattice Gauge Theory and calculate nuclear observables such as neutron lifetime. We have developed an improved algorithm that exponentially decreases the time-to solution and applied it on the new CORAL supercomputers, Sierra and Summit. We use run-time autotuning to distribute GPU resources, achieving 20% performance at low node count. We also developed optimal application mapping through a job manager, which allows CPU and GPU jobs to be interleaved, yielding 15% of peak performance when deployed across large fractions of CORAL.

preprint2016arXiv

Extra-dimensional models on the lattice

In this review we summarize the ongoing effort to study extra-dimensional gauge theories with lattice simulations. In these models the Higgs field is identified with extra-dimensional components of the gauge field. The Higgs potential is generated by quantum corrections and is protected from divergencies by the higher dimensional gauge symmetry. Dimensional reduction to four dimensions can occur through compactification or localization. Gauge-Higgs unification models are often studied using perturbation theory. Numerical lattice simulations are used to go beyond these perturbative expectations and to include non-perturbative effects. We describe the known perturbative predictions and their fate in the strongly-coupled regime for various extra-dimensional models.

preprint2016arXiv

Neutrinoless double beta decay from lattice QCD

While the discovery of non-zero neutrino masses is one of the most important accomplishments by physicists in the past century, it is still unknown how and in what form these masses arise. Lepton number-violating neutrinoless double beta decay is a natural consequence of Majorana neutrinos and many BSM theories, and many experimental efforts are involved in the search for these processes. Understanding how neutrinoless double beta decay would manifest in nuclear environments is key for understanding any observed signals. In these proceedings we present an overview of a set of one- and two-body matrix elements relevant for experimental searches for neutrinoless double beta decay, describe the role of lattice QCD calculations, and present preliminary lattice QCD results.

preprint2016arXiv

Precision lattice test of the gauge/gravity duality at large-$N$

We pioneer a systematic, large-scale lattice simulation of D0-brane quantum mechanics. The large-$N$ and continuum limits of the gauge theory are taken for the first time at various temperatures $0.4 \leq T \leq 1.0$. As a way to directly test the gauge/gravity duality conjecture we compute the internal energy of the black hole directly from the gauge theory and reproduce the coefficient of the supergravity result $E/N^2=7.41T^{14/5}$. This is the first confirmation of the supergravity prediction for the internal energy of a black hole at finite temperature coming directly from the dual gauge theory. We also constrain stringy corrections to the internal energy.

preprint2016arXiv

S-parameter and vector decay constant in QCD with eight fundamental fermions

SU(3) gauge theory with eight massless fundamental fermions seems to be near the conformal boundary, and is a candidate theory of walking technicolor. Along the series of study by LatKMI collaboration using HISQ fermions, S-parameter and vector decay constant, which provide important constraints in the model of electroweak symmetry breaking, are calculated for this theory. Use of various volumes allows a systematic investigation of finite volume effects. A strong sensitivity of the S-parameter to the volume is found.

preprint2016arXiv

Strongly interacting dynamics and the search for new physics at the LHC

We present results for the spectrum of a strongly interacting SU(3) gauge theory with $N_f = 8$ light fermions in the fundamental representation. Carrying out non-perturbative lattice calculations at the lightest masses and largest volumes considered to date, we confirm the existence of a remarkably light singlet scalar particle. We explore the rich resonance spectrum of the 8-flavor theory in the context of the search for new physics beyond the standard model at the Large Hadron Collider (LHC). Connecting our results to models of dynamical electroweak symmetry breaking, we estimate the vector resonance mass to be about 2 TeV with a width of roughly 450 GeV, and predict additional resonances with masses below ~3 TeV.

preprint2016arXiv

Supergravity from D0-brane Quantum Mechanics

The gauge/gravity duality conjecture claims the equivalence between gauge theory and superstring/M-theory. In particular, the one-dimensional gauge theory of D0-branes and type IIA string theory should agree on properties of hot black holes. Type IIA superstring theory predicts the leading $N^2$ behavior of the black hole internal energy to be $E/N^2=a_0T^{14/5}+ a_1T^{23/5}+a_2T^{29/5}+\cdots$ with the supergravity prediction $a_0=7.41$ and unknown coefficients $a_1$, $a_2$, $\ldots$ associated with stringy corrections. In order to test this duality we perform a lattice study of the gauge theory and extract a continuum, large-$N$ value of $a_0=7.4\pm 0.5$---the first direct confirmation of the supergravity prediction at finite temperature---and constrain the stringy corrections ($a_1=-9.7\pm2.2$ and $a_2=5.6\pm1.8$). We also study the sub-leading $1/N^2$ corrections to the internal energy.

preprint2016arXiv

Topological observables in many-flavour QCD

SU(3) gauge theory with eight massless flavours is believed to be walking, while the corresponding twelve- and four-flavour appear IR-conformal and confining respectively. Looking at the simulations performed by the LatKMI collaboration of these theories, we use the topological susceptibility as an additional probe of the IR dynamics. By drawing a comparison with SU(3) pure gauge theory, we see a dynamical quenching effect emerge at larger number of flavours, which is suggestive of emerging near-conformal and conformal behaviour.

preprint2016arXiv

Two-Nucleon Higher Partial-Wave Scattering from Lattice QCD

We present a determination of nucleon-nucleon scattering phase shifts for l >= 0. The S, P, D and F phase shifts for both the spin-triplet and spin-singlet channels are computed with lattice Quantum ChromoDynamics. For l > 0, this is the first lattice QCD calculation using the Luscher finite-volume formalism. This required the design and implementation of novel lattice methods involving displaced sources and momentum-space cubic sinks. To demonstrate the utility of our approach, the calculations were performed in the SU(3)-flavor limit where the light quark masses have been tuned to the physical strange quark mass, corresponding to m_pi = m_K ~ 800 MeV. In this work, we have assumed that only the lowest partial waves contribute to each channel, ignoring the unphysical partial wave mixing that arises within the finite-volume formalism. This assumption is only valid for sufficiently low energies; we present evidence that it holds for our study using two different channels. Two spatial volumes of V ~ (3.5 fm)^3 and V ~ (4.6 fm)^3 were used. The finite-volume spectrum is extracted from the exponential falloff of the correlation functions. Said spectrum is mapped onto the infinite volume phase shifts using the generalization of the Luscher formalism for two-nucleon systems.

preprint2016arXiv

Walking and conformal dynamics in many-flavor QCD

In the search for a realistic walking technicolor model, QCD with many flavors is an attractive candidate. From the series of studies by the LatKMI collaboration, we present updated results of the scaling properties of various hadron spectra, including the (pseudo)scalar, vector, and baryon channels, for $N_f=8$ QCD analyzed with the HISQ action. By comparing these with $N_f=12$ QCD, which has properties consistent with conformality, possible signals of walking dynamics are discussed. We also present a preliminary result of the flavor-singlet pseudoscalar mass in many-flavor QCD.

preprint2015arXiv

Conformality in twelve-flavour QCD

The spectrum of twelve-flavor QCD has been studied in details by the LatKMI collaboration. In this proceeding we present our updated results for the spectrum obtained with the HISQ action at two lattice spacings, several volumes and fermion masses. In particular, we emphasize the existence of a flavor-singlet scalar state parametrically light with respect to the rest of the spectrum, first reported in our paper. This feature is expected to be present for theories in the conformal window, but the lattice calculation of such a state is difficult and requires noise-reduction techniques together with large statistics, in order to evaluate disconnected diagrams. Being able to provide a robust observed connection between a light flavor-singlet scalar and (near-)conformality is an important step towards observing a light composite Higgs boson in walking technicolor theories on the lattice. We also show updated results for the mass anomalous dimension $γ_m$ obtained from various spectral quantities, including the string tension, under the assumption that the theory is inside the conformal window.

preprint2015arXiv

Direct Detection of Stealth Dark Matter through Electromagnetic Polarizability

We calculate the spin-independent scattering cross section for direct detection that results from the electromagnetic polarizability of a composite scalar baryon dark matter candidate -- "Stealth Dark Matter", that is based on a dark SU(4) confining gauge theory. In the nonrelativistic limit, electromagnetic polarizability proceeds through a dimension-7 interaction leading to a very small scattering cross section for dark matter with weak scale masses. This represents a lower bound on the scattering cross section for composite dark matter theories with electromagnetically charged constituents. We carry out lattice calculations of the polarizability for the lightest baryons in SU(3) and SU(4) gauge theories using the background field method on quenched configurations. We find the polarizabilities of SU(3) and SU(4) to be comparable (within about 50%) normalized to the baryon mass, which is suggestive for extensions to larger SU(N) groups. The resulting scattering cross sections with a xenon target are shown to be potentially detectable in the dark matter mass range of about 200-700 GeV, where the lower bound is from the existing LUX constraint while the upper bound is the coherent neutrino background. Significant uncertainties in the cross section remain due to the more complicated interaction of the polarizablity operator with nuclear structure, however the steep dependence on the dark matter mass, $1/m_B^6$, suggests the observable dark matter mass range is not appreciably modified. We briefly highlight collider searches for the mesons in the theory as well as the indirect astrophysical effects that may also provide excellent probes of stealth dark matter.

preprint2015arXiv

Lattice QCD input for axion cosmology

One intriguing BSM particle is the QCD axion, which could simultaneously provide a solution to the Strong CP problem and account for some, if not all, of the dark matter density in the universe. This particle is a pNGB of the conjectured Peccei-Quinn (PQ) symmetry of the Standard Model. Its mass and interactions are suppressed by a heavy symmetry breaking scale, $f_a$, whose value is roughly greater than $10^{9}$ GeV (or, conversely, the axion mass, $m_a$, is roughly less than $10^4\ μ\text{eV}$). The density of axions in the universe, which cannot exceed the relic dark matter density and is a quantity of great interest in axion experiments like ADMX, is a result of the early-universe interplay between cosmological evolution and the axion mass as a function of temperature. The latter quantity is proportional to the second derivative of the QCD free energy with respect to the CP-violating phase, $θ$. However, this quantity is generically non-perturbative and previous calculations have only employed instanton models at the high temperatures of interest (roughly 1 GeV). In this and future works, we aim to calculate the temperature-dependent axion mass at small $θ$ from first-principle lattice calculations, with controlled statistical and systematic errors. Once calculated, this temperature-dependent axion mass is input for the classical evolution equations of the axion density of the universe. Due to a variety of lattice systematic effects at the very high temperatures required, we perform a calculation of the leading small-$θ$ cumulant of the theta vacua on large volume lattices for SU(3) Yang-Mills with high statistics as a first proof of concept, before attempting a full QCD calculation in the future. From these pure glue results, the misalignment mechanism yields the axion mass bound $m_a \geq (14.6\pm0.1) \ μ\text{eV}$ when PQ-breaking occurs after inflation.

preprint2015arXiv

Lattice study for conformal windows of SU(2) and SU(3) gauge theories with fundamental fermions

We present our investigation of SU(2) gauge theory with 8 flavours, and SU(3) gauge theory with 12 flavours. For the SU(2) case, at strong bare coupling, $β\lesssim 1.45$, the distribution of the lowest eigenvalue of the Dirac operator can be described by chiral random matrix theory for the Gaussian symplectic ensemble. Our preliminary result indicates that the chiral phase transition in this theory is of bulk nature. For the SU(3) theory, we use high-precision lattice data to perform the step-scaling study of the coupling, $g_{\rm GF}$, in the Gradient Flow scheme. We carefully examine the reliability of the continuum extrapolation in the analysis, and conclude that the scaling behaviour of this SU(3) theory is not governed by possible infrared conformality at $g_{\rm GF}^{2} \lesssim 6$.

preprint2015arXiv

Lattice study of the scalar and baryon spectra in many-flavor QCD

In the search for a composite Higgs boson in walking technicolor models, many flavor QCD, in particular with $N_f=8$, is an attractive candidate, and has been found to have a composite flavor-singlet scalar as light as the pion. Based on lattice simulations of this theory with the HISQ action, we will present our preliminary results on the scalar decay constant using the fermionic bilinear operator, and on the mass of the lightest baryon state which could be a dark matter candidate. Combining these two results, implications for dark matter direct detection are also discussed.

preprint2015arXiv

Nuclear Parity Violation from Lattice QCD

The electroweak interaction at the level of quarks and gluons are well understood from precision measurements in high energy collider experiments. Relating these fundamental parameters to Hadronic Parity Violation in nuclei however remains an outstanding theoretical challenge. One of the most interesting observables in this respect is the parity violating hadronic neutral current: it is hard to measure in collider experiments and is thus the least constrained observable of the Standard Model. Precision measurements of parity violating transitions in nuclei can help to improve these constraints. In these systems however, the weak interaction is masked by effects of the seven orders of magnitude stronger non-perturbative strong interaction. Therefore, in order to relate experimental measurements of the parity violating pion-nucleon couplings to the fundamental Lagrangian of the SM, these non-perturbative effects have to be well understood. In this paper, we are going to present a Lattice QCD approach for computing the $ΔI{=}2$ parity violating matrix element in proton proton scattering. This process does not involve disconnected diagrams in the isospin symmetric limit and is thus a perfect testbed for studying the feasibility of the more involved calculation of the parity violating pion-nucleon coupling.

preprint2015arXiv

Stealth Dark Matter: Dark scalar baryons through the Higgs portal

We present a new model of "Stealth Dark Matter": a composite baryonic scalar of an $SU(N_D)$ strongly-coupled theory with even $N_D \geq 4$. All mass scales are technically natural, and dark matter stability is automatic without imposing an additional discrete or global symmetry. Constituent fermions transform in vector-like representations of the electroweak group that permit both electroweak-breaking and electroweak-preserving mass terms. This gives a tunable coupling of stealth dark matter to the Higgs boson independent of the dark matter mass itself. We specialize to $SU(4)$, and investigate the constraints on the model from dark meson decay, electroweak precision measurements, basic collider limits, and spin-independent direct detection scattering through Higgs exchange. We exploit our earlier lattice simulations that determined the composite spectrum as well as the effective Higgs coupling of stealth dark matter in order to place bounds from direct detection, excluding constituent fermions with dominantly electroweak-breaking masses. A lower bound on the dark baryon mass $m_B \gtrsim 300$ GeV is obtained from the indirect requirement that the lightest dark meson not be observable at LEP II. We briefly survey some intriguing properties of stealth dark matter that are worthy of future study, including: collider studies of dark meson production and decay; indirect detection signals from annihilation; relic abundance estimates for both symmetric and asymmetric mechanisms; and direct detection through electromagnetic polarizability, a detailed study of which will appear in a companion paper.

preprint2015arXiv

SU(3) gauge theory with four degenerate fundamental fermions on the lattice

As a part of the project studying large $N_f$ QCD, the LatKMI Collaboration has been investigating the SU(3) gauge theory with four fundamental fermions (four-flavor QCD). The main purpose of studying four-flavor QCD is to provide a qualitative comparison to $N_f= 8$, $12$, $16$ QCD; however, a quantitative comparison to real-world QCD is also interesting. To make such comparisons more meaningful, it is desirable to use the same kind of lattice action consistently, so that qualitative difference of different theories are less affected by artifacts of lattice discretization. Here, we adopt the highly-improved staggered quark action with the tree-level Symanzik gauge action (HISQ/tree), which is exactly the same as the setup for our simulations for $SU(3)$ gauge theories with $N_f=8$, $12$ and $16$ fundamental fermions~\cite{Aoki:2013xza, Aoki:2012eq, Aoki:2014oma}. In the next section, we show the fermion mass dependence of $F_π$, $\langle\barψψ\rangle$, $M_π$, $M_ρ$, $M_N$ and their chiral extrapolations. In section 3, preliminary results of the measurement of the mass of the flavor-singlet scalar bound state will be reported.

preprint2015arXiv

Two-nucleon scattering in multiple partial waves

We determine scattering phase shifts for S,P,D, and F partial wave channels in two-nucleon systems using lattice QCD methods. We use a generalization of Luscher's finite volume method to determine infinite volume phase shifts from a set of finite volume ground- and excited-state energy levels on two volumes, V=(3.4 fm)^3 and V=(4.5 fm)^3. The calculations are performed in the SU(3)-flavor limit, corresponding to a pion mass of approximately 800 MeV. From the energy dependence of the phase shifts we are able to extract scattering parameters corresponding to an effective range expansion.

preprint2014arXiv

Composite bosonic baryon dark matter on the lattice: SU(4) baryon spectrum and the effective Higgs interaction

We present the spectrum of baryons in a new SU(4) gauge theory with fundamental fermion constituents. The spectrum of these bosonic baryons is of significant interest for composite dark matter theories. Here, we compare the spectrum and properties of SU(3) and SU(4) baryons, and then compute the dark-matter direct detection cross section via Higgs boson exchange for TeV-scale composite dark matter arising from a confining SU(4) gauge sector. Comparison with the latest LUX results leads to tight bounds on the fraction of the constituent-fermion mass that may arise from electroweak symmetry breaking. Lattice calculations of the dark matter mass spectrum and the Higgs-dark matter coupling are performed on quenched $16^{3} \times 32$, $32^{3} \times 64$, $48^{3} \times 96$, and $64^{3} \times128$ lattices with three different lattice spacings, using Wilson fermions with moderate to heavy pseudoscalar meson masses. Our results lay a foundation for future analytic and numerical study of composite baryonic dark matter.

preprint2014arXiv

Light composite scalar in eight-flavor QCD on the lattice

We present the first observation of a flavor-singlet scalar meson as light as the pion in $N_f=8$ QCD on the lattice, using the Highly Improved Staggered Quark action. Such a light scalar meson can be regarded as a composite Higgs with mass 125 GeV. In accord with our previous lattice results showing that the theory exhibits walking behavior, the light scalar may be a technidilaton, a pseudo Nambu-Goldstone boson of the approximate scale symmetry in walking technicolor.

preprint2014arXiv

Phase Structure Study of SU(2) Lattice Gauge Theory with 8 Flavors

We present the investigation of the strong bare-coupling regime of SU(2) lattice gauge theory with 8 fermion flavors in the fundamental representation. The simulations are performed with unimproved staggered fermions and the plaquette gauge action. One bulk phase transition is observed through the measurement of the plaquette. The results of cold-start and hot-start simulations, as well as the hysteresis study, indicate the order of this transition can be weakly first order. Using the smeared Polyakov loops, and a method inspired by the constraint effective potential, we study the vacuum structure near the confining-deconfining phase transition. The Dirac operator eigenvalue spectrum is investigated, where further analysis is needed to clarify the properties of the chiral phase structure.

preprint2013arXiv

A light composite scalar in eight-flavor QCD on the lattice

In search for a composite Higgs boson (techni-dilaton) in the walking technicolor, we present our preliminary results on the first observation of a light flavor-singlet scalar in a candidate theory for the walking technicolor, the Nf=8 QCD, which was found in our previous paper to have spontaneous chiral symmetry breaking together with remnants of the conformality. Based on simulations with the HISQ-type action on several lattice sizes with various fermion masses, we find evidence of a flavor-singlet scalar meson with mass comparable to that of the Nambu-Goldstone pion in both the small fermion-mass region, where chiral perturbation theory works, and the intermediate fermion-mass region where the hyperscaling relation holds. We further discuss its chiral limit extrapolation in comparison with other states studied in our previous paper: the scalar has a mass much smaller than that of the vector meson, which is compared to the Nambu-Goldstone pion having a vanishing mass in that limit.

preprint2013arXiv

Gluonic observables and the scalar spectrum of twelve-flavour QCD

We measure glueball masses and the string tension in twelve-flavour QCD, aiming at comparing the emerging gluonic spectrum to the mesonic one. When approaching the critical surface at zero quark mass, the hierarchy of masses in the different sectors of the spectrum gives a new handle to determine the existence of an infrared fixed point. We describe the details of our gluonic measurements and the results obtained on a large number of gauge configurations generated with the HISQ action. In particular, we focus on the scalar glueball and its mixing with a flavour-singlet fermionic state, which is lighter than the pseudoscalar (would-be pion) state. The results are interesting in view of a light composite Higgs boson in walking technicolor theories.

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

Light composite scalar in twelve-flavor QCD on the lattice

Based on lattice simulations using highly improved staggered quarks for twelve-flavor QCD with several bare fermion masses, we observe a flavor-singlet scalar state lighter than the pion in the correlators of fermionic interpolating operators. The same state is also investigated using correlators of gluonic interpolating operators. Combined with our previous study, that showed twelve-flavor QCD to be consistent with being in the conformal window, we infer that the lightness of the scalar state is due to infrared conformality. This result shed some light on the possibility of a light composite Higgs boson ("technidilaton") in walking technicolor theories.

preprint2013arXiv

Searching for a continuum 4D field theory arising from a 5D non-abelian gauge theory

The anisotropic 5D SU(2) Yang-Mills model has been widely investigated on the lattice during the last decade. In the case where all dimensions are large in size, it was previously claimed that there is a new phase in the phase diagram, called the Layer phase. In this phase, the gauge fields would be localized on 4D layers. Previous works claim that the phase transition to the Layer phase is of second order, which would allow a continuum limit to be taken. We present the extension of the previous work to large lattices, for which we found a first order phase transition. This leaves the scenario that this 5D theory can be dimensionally reduced to a continuum 4D field theory, doubtful.

preprint2013arXiv

The transition to a layered phase in the anisotropic five-dimensional SU(2) Yang-Mills theory

We extend to large lattices the work of a previous investigation of the phase diagram of the anisotropic five-dimensional SU(2) Yang-Mills model using Monte Carlo simulations in the regime where the lattice spacing in the fifth dimension is larger than in the other four dimensions. We find a first order phase transition between the confining and deconfining phase at the anisotropic parameter point $β_4=2.60$ which was previously claimed to be the critical point at which the order of the transition changes from first to second. We conclude that large lattices are required to establish the first order nature of this line of transitions and consequently that the scenario of dimensional reduction of the five-dimensional theory to a continuum four-dimensional theory via the existence of the so-called "layer phase" is unpromising.

preprint2012arXiv

Light scalar spectrum in extra-dimensional gauge theories

The phase diagram of five-dimensional SU(2) gauge theories with one compactified dimension on anisotropic lattices has a rich structure. In this contribution we show how to control non-perturbatively the scale hierarchy between the cut-off and the compactification scale in the bare parameter space. There exists a set of strong bare couplings where the the five-dimensional lattice theory can be described by an effective four-dimensional theory with a scalar field in the adjoint representation. We present a detailed study of the light scalar spectrum as it arises from the non-perturbative dynamics of the full five-dimensional lattice theory. We also investigate the mixing with scalar glueball states in the attempt to further establish the extra-dimensional nature of light scalar states.

preprint2012arXiv

Light scalars in strongly-coupled extra-dimensional theories

The low-energy dynamics of five-dimensional Yang-Mills theories compactified on S^1 can be described by a four-dimensional gauge theory coupled to a scalar field in the adjoint representation of the gauge group. Perturbative calculations suggest that the mass of this elementary scalar field is protected against power divergences, and is controlled by the size of the extra dimension R. As a first step in the study of this phenomenon beyond perturbation theory, we investigate the phase diagram of a SU(2) Yang-Mills theory in five dimensions regularized on anisotropic lattices and we determine the ratios of the relevant physical scales. The lattice system shows a dimensionally reduced phase where the four-dimensional correlation length is much larger than the size of the extra dimension, but still smaller than the four-dimensional volume. In this region of the bare parameter space, at energies below 1/R, the non-perturbative spectrum contains a \emph{light} scalar state. This state has a mass that is independent of the cut-off, and a small overlap with glueball operators. Our results suggest that light scalar fields can be introduced in a lattice theory using compactified extra dimensions, rather than fine tuning the bare mass parameter.

preprint2012arXiv

Scalar mass corrections from compact extra dimensions on the lattice

We explore the phase diagram of the SU(2) Yang-Mills theory in 5 dimensions by numerical simulations. The lattice system shows a dimensionally-reduced phase where the extra dimension is small compared to the four dimensional correlation length. In the low-energy regime of this phase, the system behaves like a four-dimensional gauge theory coupled to an adjoint scalar field. By tuning the bare parameters of the lattice model, we identify lines of constant physics, and analyse the behaviour of the non-perturbative scalar mass as a function of the compactification and the cut-off scales. The perturbative prediction that the effective theory contains a light particle with a mass that is independent of the cut-off is tested against non-pertubative results.

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.

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.

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.