
A black hole described in $SU(N)$ gauge theory consists of $N$ Dbranes. By
separating one of the Dbranes 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 zerobrane 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' D0brane 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 nonuniform viscosity on their energy and stability in minimal
channel flows. What we seek is a theoretical support rooted in a fully
nonlinear framework that explains the modified threshold for transition to
turbulence in flows with temperaturedependent 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 selfsustained 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 manyflavors strongly
interacting theories with QCD, we study the flavorsinglet spectrum as an
interesting probe. In fact, some composite Higgs models require the Higgs boson
to be the lightest flavorsinglet 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 flavorsinglet pseudoscalar. We present
the ongoing study of these flavorsinglet channels using multiple
interpolating operators on highstatistics 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 stateoftheart results from a lattice QCD calculation of the
nucleon axial coupling, $g_A$, using M\"obius DomainWall fermions solved on
the dynamical $N_f = 2 + 1 + 1$ HISQ ensembles after they are smeared using the
gradientflow algorithm. Relevant threepoint correlation functions are
calculated using a method inspired by the FeynmanHellmann 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 signaltonoise problem which has prevented a
controlled extrapolation to the physical point. The variational method has been
applied very successfully to twomeson systems, allowing for the extraction of
the twomeson states very early in Euclidean time through the use of improved
single hadron operators. The sheer numerical cost of using the same techniques
in twobaryon 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 twobaryon correlation functions.
To test this method, we perform calculations on the WM/JLab isoclover 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
twonucleon correlation functions. We observe the calm baryon removes the
excited state contamination from the twonucleon correlation function to as
early a time as the singlenucleon is improved, provided nonlocal (displaced
nucleon) sources are used. For the local twonucleon correlation function
(where both nucleons are created from the same spacetime 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 domainwall 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 gradientflow scheme to first smear the HISQ configurations. The
greater numerical cost of the M\"{o}bius domainwall 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
flowtime 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 domainwall 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 flowtime, 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
"neutronantineutron" oscillation in the dark sector can then result in
nonzero 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
flavorsinglet scalar meson (denoted as $\sigma$) as light as the pion and
various "walking signals" through lowlying 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
flavorsinglet scalar meson as a technidilaton, which can be identified with
the 125 GeV composite Higgs in $N_f=8$ onefamily model.

We report on a lattice QCD calculation of the nucleon axial charge, $g_A$,
using M\"{o}bius DomainWall fermions solved on the dynamical $N_f=2+1+1$ HISQ
ensembles after they are smeared using the gradientflow 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 FeynmanHellmann 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 nucleonnucleon scattering phase shifts for l
>= 0. The S, P, D and F phase shifts for both the spintriplet and spinsinglet
channels are computed with lattice Quantum ChromoDynamics. For l > 0, this is
the first lattice QCD calculation using the Luscher finitevolume formalism.
This required the design and implementation of novel lattice methods involving
displaced sources and momentumspace 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 finitevolume 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 finitevolume 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 twonucleon systems.

While the discovery of nonzero 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 numberviolating 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 twobody 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 extradimensional
gauge theories with lattice simulations. In these models the Higgs field is
identified with extradimensional 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. GaugeHiggs
unification models are often studied using perturbation theory. Numerical
lattice simulations are used to go beyond these perturbative expectations and
to include nonperturbative effects. We describe the known perturbative
predictions and their fate in the stronglycoupled regime for various
extradimensional models.

The gauge/gravity duality conjecture claims the equivalence between gauge
theory and superstring/Mtheory. In particular, the onedimensional gauge
theory of D0branes 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
temperatureand constrain the stringy corrections ($a_1=9.7\pm2.2$ and
$a_2=5.6\pm1.8$). We also study the subleading $1/N^2$ corrections to the
internal energy.

We pioneer a systematic, largescale lattice simulation of D0brane 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 nonperturbative 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
8flavor 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 fourflavour appear IRconformal 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 nearconformal 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
flavorsinglet pseudoscalar mass in manyflavor 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
(fourflavor QCD). The main purpose of studying fourflavor QCD is to provide a
qualitative comparison to $N_f= 8$, $12$, $16$ QCD; however, a quantitative
comparison to realworld 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
highlyimproved staggered quark action with the treelevel 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 flavorsinglet 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
highprecision lattice data to perform the stepscaling 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 nonperturbative strong interaction.
Therefore, in order to relate experimental measurements of the parity violating
pionnucleon couplings to the fundamental Lagrangian of the SM, these
nonperturbative 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 pionnucleon coupling.

We determine scattering phase shifts for S,P,D, and F partial wave channels
in twonucleon 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 excitedstate 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 flavorsinglet 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 staggeredfermion action on large volumes and at
small fermion masses, reaching $M_{\rho}/M_{\pi} \approx 2.2$. The axialvector
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 flavorsinglet 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 manyflavor 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.