• A black hole described in $SU(N)$ gauge theory consists of $N$ D-branes. By separating one of the D-branes from others and studying the interaction between them, the black hole geometry can be probed. In order to obtain quantitative results, we employ the lattice Monte Carlo simulation. As a proof of the concept, we perform an explicit calculation in the matrix model dual to the black zero-brane in type IIA string theory. We demonstrate this method actually works in the high temperature region, where the stringy correction is large. We argue possible dual gravity interpretations.
  • We study the thermodynamics of the `ungauged' D0-brane matrix model by Monte Carlo simulation. Our results appear to be consistent with the conjecture by Maldacena and Milekhin.
  • Motivated by the relevance of edge state solutions as mediators of transition, we use direct numerical simulations to study the effect of spatially non-uniform viscosity on their energy and stability in minimal channel flows. What we seek is a theoretical support rooted in a fully non-linear framework that explains the modified threshold for transition to turbulence in flows with temperature-dependent viscosity. Consistently over a range of subcritical Reynolds numbers, we find that decreasing viscosity away from the walls weakens the streamwise streaks and the vortical structures responsible for their regeneration. The entire self-sustained cycle of the edge state is maintained on a lower kinetic energy level with a smaller driving force, compared to a flow with constant viscosity. Increasing viscosity away from the walls has the opposite effect. In both cases, the effect is proportional to the strength of the viscosity gradient. The results presented highlight a local shift in the state space of the position of the edge state relative to the laminar attractor with the consequent modulation of its basin of attraction in the proximity of the edge state and of the surrounding manifold. The implication is that the threshold for transition is reduced for perturbations evolving in the neighbourhood of the edge state in case viscosity decreases away from the walls, and vice versa.
  • Studying SU(3) gauge theories with increasing number of light fermions is relevant both for understanding the strong dynamics of QCD and for constructing strongly interacting extensions of the Standard Model (e.g. UV completions of composite Higgs models). In order to contrast these many-flavors strongly interacting theories with QCD, we study the flavor-singlet spectrum as an interesting probe. In fact, some composite Higgs models require the Higgs boson to be the lightest flavor-singlet scalar in the spectrum of a strongly interacting new sector with a well defined hierarchy with the rest of the states. Moreover, introducing many light flavors at fixed number of colors can influence the dynamics of the lightest flavor-singlet pseudoscalar. We present the on-going study of these flavor-singlet channels using multiple interpolating operators on high-statistics ensembles generated by the LatKMI collaboration and we compare results with available data obtained by the Lattice Strong Dynamics collaboration. For the theory with 8 flavors, the two collaborations have generated configurations that complement each others with the aim to tackle the massless limit using the largest possible volumes.
  • We present state-of-the-art results from a lattice QCD calculation of the nucleon axial coupling, $g_A$, using M\"obius Domain-Wall fermions solved on the dynamical $N_f = 2 + 1 + 1$ HISQ ensembles after they are smeared using the gradient-flow algorithm. Relevant three-point correlation functions are calculated using a method inspired by the Feynman-Hellmann theorem, and demonstrate significant improvement in signal for fixed stochastic samples. The calculation is performed at five pion masses of $m_\pi\sim \{400, 350, 310, 220, 130\}$~MeV, three lattice spacings of $a\sim\{0.15, 0.12, 0.09\}$~fm, and we do a dedicated volume study with $m_\pi L\sim\{3.22, 4.29, 5.36\}$. Control over all relevant sources of systematic uncertainty are demonstrated and quantified. We achieve a preliminary value of $g_A = 1.285(17)$, with a relative uncertainty of 1.33\%.
  • Outstanding problems in nuclear physics require input and guidance from lattice QCD calculations of few baryons systems. However, these calculations suffer from an exponentially bad signal-to-noise problem which has prevented a controlled extrapolation to the physical point. The variational method has been applied very successfully to two-meson systems, allowing for the extraction of the two-meson states very early in Euclidean time through the use of improved single hadron operators. The sheer numerical cost of using the same techniques in two-baryon systems has been prohibitive. We present an alternate strategy which offers some of the same advantages as the variational method while being significantly less numerically expensive. We first use the Matrix Prony method to form an optimal linear combination of single baryon interpolating fields generated from the same source and different sink interpolators. Very early in Euclidean time this linear combination is numerically free of excited state contamination, so we coin it a calm baryon. This calm baryon operator is then used in the construction of the two-baryon correlation functions. To test this method, we perform calculations on the WM/JLab iso-clover gauge configurations at the SU(3) flavor symmetric point with m{\pi} $\sim$ 800 MeV --- the same configurations we have previously used for the calculation of two-nucleon correlation functions. We observe the calm baryon removes the excited state contamination from the two-nucleon correlation function to as early a time as the single-nucleon is improved, provided non-local (displaced nucleon) sources are used. For the local two-nucleon correlation function (where both nucleons are created from the same space-time location) there is still improvement, but there is significant excited state contamination in the region the single calm baryon displays no excited state contamination.
  • We report on salient features of a mixed lattice QCD action using valence M\"{o}bius domain-wall fermions solved on the dynamical $N_f=2+1+1$ HISQ ensembles generated by the MILC Collaboration. The approximate chiral symmetry properties of the valence fermions are shown to be significantly improved by utilizing the gradient-flow scheme to first smear the HISQ configurations. The greater numerical cost of the M\"{o}bius domain-wall inversions is mitigated by the highly efficient QUDA library optimized for NVIDIA GPU accelerated compute nodes. We have created an interface to this optimized QUDA solver in Chroma. We provide tuned parameters of the action and performance of QUDA using ensembles with the lattice spacings $a \simeq \{0.15, 0.12, 0.09\}$ fm and pion masses $m_\pi \simeq \{310, 220,130\}$ MeV. We have additionally generated two new ensembles with $a\sim0.12$ fm and $m_\pi\sim\{400, 350\}$ MeV. With a fixed flow-time of $t_{gf}=1$ in lattice units, the residual chiral symmetry breaking of the valence fermions is kept below 10\% of the light quark mass on all ensembles, $m_{res} \lesssim 0.1\times m_l$, with moderate values of the fifth dimension $L_5$ and a domain-wall height $M_5 \leq 1.3$. As a benchmark calculation, we perform a continuum, infinite volume, physical pion and kaon mass extrapolation of $F_{K^\pm}/F_{\pi^\pm}$ and demonstrate our results are independent of flow-time, and consistent with the FLAG determination of this quantity at the level of less than one standard deviation.
  • We entertain the possibility that neutrino masses and dark matter (DM) originate from a common composite dark sector. A minimal effective theory can be constructed based on a dark $SU(3)_D$ interaction with three flavors of massless dark quarks; electroweak symmetry breaking gives masses to the dark quarks. By assigning a $\mathbb Z_2$ charge to one flavor, a stable "dark kaon" can provide a good thermal relic DM candidate. We find that "dark neutrons" may be identified as right handed Dirac neutrinos. Some level of "neutron-anti-neutron" oscillation in the dark sector can then result in non-zero Majorana masses for light Standard Model neutrinos. A simple ultraviolet completion is presented, involving additional heavy $SU(3)_D$-charged particles with electroweak and lepton Yukawa couplings. At our benchmark point, there are "dark pions" that are much lighter than the Higgs and we expect spectacular collider signals arising from the UV framework. This includes the decay of the Higgs boson to $\tau \tau \ell \ell'$, where $\ell$($\ell'$) can be any lepton, with displaced vertices. We discuss the observational signatures of this UV framework in dark matter searches and primordial gravitational wave experiments; the latter signature is potentially correlated with the $H \to \tau \tau \ell \ell'$ decay.
  • Based on the highly improved staggered quark action, we perform lattice simulations of $N_f=8$ QCD and confirm our previous observation of a flavor-singlet scalar meson (denoted as $\sigma$) as light as the pion and various "walking signals" through low-lying spectra, with higher statistics, smaller fermion masses $m_f$, and larger volumes. We measure $M_\pi$, $F_\pi$, $M_\rho$, $M_{a_0}$, $M_{a_1}$, $M_{b_1}$, $M_N$, $M_\sigma$, $F_\sigma$, $\langle \bar{\psi} \psi\rangle$ (both directly and through the GMOR relation), and the string tension. The data are consistent with the spontaneously broken phase of the chiral symmetry, in agreement with the previous results: ratios of the quantities to $M_\pi$ monotonically increase in the smaller $m_f$ region towards the chiral limit similarly to $N_f=4$ QCD, in sharp contrast to $N_f=12$ QCD where the ratios become flattened. The hyperscaling relation holds with roughly a universal value of the anomalous dimension, $\gamma_m \simeq 1$, with a notable exception of $M_\pi$ with $\gamma_m \simeq 0.6$ as in the previous results. This is a salient feature ("walking signal") of $N_f=8$, unlike either $N_f=4$ which has no hyperscaling relation at all, or $N_f=12$ QCD which exhibits universal hyperscaling. We further confirm the previous observation of the light $\sigma$ with mass comparable to the pion in the studied $m_f$ region. In a chiral limit extrapolation of the $\sigma$ mass using the dilaton chiral perturbation theory and also using the simple linear fit, we find the value consistent with the 125 GeV Higgs boson within errors. Our results suggest that the theory could be a good candidate for walking technicolor model, having anomalous dimension $\gamma_m \simeq 1$ and a light flavor-singlet scalar meson as a technidilaton, which can be identified with the 125 GeV composite Higgs in $N_f=8$ one-family model.
  • We report on a lattice QCD calculation of the nucleon axial charge, $g_A$, using M\"{o}bius Domain-Wall fermions solved on the dynamical $N_f=2+1+1$ HISQ ensembles after they are smeared using the gradient-flow algorithm. The calculation is performed with three pion masses, $m_\pi\sim\{310,220,130\}$ MeV. Three lattice spacings ($a\sim\{0.15,0.12,0.09\}$ fm) are used with the heaviest pion mass, while the coarsest two spacings are used on the middle pion mass and only the coarsest spacing is used with the near physical pion mass. On the $m_\pi\sim220$ MeV, $a\sim0.12$ fm point, a dedicated volume study is performed with $m_\pi L \sim \{3.22,4.29,5.36\}$. Using a new strategy motivated by the Feynman-Hellmann Theorem, we achieve a precise determination of $g_A$ with relatively low statistics, and demonstrable control over the excited state, continuum, infinite volume and chiral extrapolation systematic uncertainties, the latter of which remains the dominant uncertainty. Our final determination at 2.6\% total uncertainty is $g_A = 1.278(21)(26)$, with the first uncertainty including statistical and systematic uncertainties from fitting and the second including model selection systematics related to the chiral and continuum extrapolation. The largest reduction of the second uncertainty will come from a greater number of pion mass points as well as more precise lattice QCD results near the physical pion mass.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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_\pi$, $\langle\bar{\psi}\psi\rangle$, $M_\pi$, $M_\rho$, $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.
  • 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, $\beta \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$.
  • 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 $\Delta 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.
  • 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.
  • 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.
  • The Lattice Strong Dynamics collaboration is investigating the properties of a SU(3) gauge theory with $N_f = 8$ light fermions on the lattice. We measure the masses of the lightest pseudoscalar, scalar and vector states using simulations with the nHYP staggered-fermion action on large volumes and at small fermion masses, reaching $M_{\rho}/M_{\pi} \approx 2.2$. The axial-vector meson and the nucleon are also studied for the same range of fermion masses. One of the interesting features of this theory is the dynamical presence of a light flavor-singlet scalar state with $0^{++}$ quantum numbers that is lighter than the vector resonance and has a mass consistent with the one of the pseudoscalar state for the whole fermion mass range explored. We comment on the existence of such state emerging from our lattice simulations and on the challenges of its analysis. Moreover we highlight the difficulties in pursuing simulations in the chiral regime of this theory using large volumes.
  • LatKMI Collaboration discusses the topological insights in many-flavor QCD on the lattice. We explore walking/conformal/confining phase in $N_\mathrm{f}$ = 4, 8 and 12 (in particular $N_\mathrm{f}$ = 8) lattice QCD via the topological charge and susceptibility, eigenvalues and anomalous dimension.