
We study the Bell inequality in a holographic model of the casually
disconnected EinsteinPodolskyRosen (EPR) pair. The
ClauserHorneShimonyHolt(CHSH) form of Bell inequality is constructed using
holographic SchwingerKeldysh (SK) correlators. We show that the manifestation
of quantum correlation in Bell inequality can be holographically reproduced
from the classical fluctuations of dual accelerating string in the bulk
gravity. The violation of this holographic Bell inequality supports the
essential quantum property of this holographic model of an EPR pair.

We explore the possibility of detecting entangled photon pairs from cosmic
microwave background or other cosmological sources coming from two patches of
the sky. The measurements use two detectors with different photon polarizer
directions. When two photon sources are separated by a large angle relative to
the earth, such that each detector has only one photon source in its field of
view, a null test of unentangled photons can be performed. The deviation from
this unentangled background is, in principle, the signature of photon
entanglement. To confirm whether the deviation is consistent with entangled
photons, we derive a photon polarization correlation to compare with, similar
to that in a Bell inequality measurement. However, since photon coincidence
measurement cannot be used to discriminate unentangled cosmic photons, it is
unlikely that the correlation expectation value alone can violate Bell
inequality to provide the signature for entanglement.

Using symmetry properties, we determine the mixing pattern of a class of
nonlocal quark bilinear operators containing a straight Wilson line along a
spatial direction. We confirm the previous study that mixing among the lowest
dimensional operators, which have mass dimension equals three, can occur if
chiral symmetry is broken in the lattice action. For higher dimensional
operators, we find that the dimension three operators will always mix with
dimension four operators even if chiral symmetry is preserved. Also, the number
of dimension four operators involved in the mixing is large hence it is
impractical to remove the mixing by the improvement procedure. Our result is
important to determining the Bjorken$x$ dependence parton distribution
functions using the quasidistribution method on a Euclidean lattice. The
requirement of using large hadron momentum in this approach makes the control
of errors from dimension four operators even more important.

We show that the correct way to extract parton distribution functions from
the reduced Ioffetime distribution (RITD), a ratio of the Ioffetime
distribution for a moving hadron and a hadron at rest, is through a
factorization formula. This factorization exists because, at small distances,
forming the ratio does not change the infrared behavior of the numerator, which
is factorizable. We illustrate the effect of such a factorization by applying
it to results in the literature.

We present the first direct latticeQCD calculation of the Bjorken$x$
dependence of the valence quark distribution in the pion. Using largemomentum
effective theory (LaMET), we calculate the boosted pion state with long Wilson
link operators. With lattice data calculated at 310 MeV pion mass, we
illustrate the numerical impact of oneloop matching and meson mass
corrections. Our results are comparable quantitatively with the results
extracted from experimental data as well as from the DysonSchwinger equation.
Future calculations at physical pion mass and larger momentum will be able to
discern discrepancies in various existing analyses.

We present a latticeQCD calculation of the unpolarized isovector parton
distribution function (PDF) using ensembles at the physical pion mass with the
large proton boost momentum $P_z \in \{2.2,3.0\}$ GeV within the framework of
largemomentum effective theory (LaMET). In contrast with our previous
physicalpion PDF result, we increase the statistics significantly, double the
boost momentum, increase the investment in excitedstate contamination
systematics, and switch to $\gamma_t$ operator to avoid mixing with scalar
matrix elements. We used four sourcesink separations in our analysis to
control the systematics associated with excitedstate contamination. The
oneloop LaMET matching corresponding to the new operator is calculated and
applied to our lattice data. We detail the systematics that affect PDF
calculations, providing a guide such that future calculations can further
improve on the precision of lattice PDF. We find our final parton distribution
in reasonable agreement with the PDF provided by the CT14 global analysis.

We present the first latticeQCD calculation of the kaon distribution
amplitude using the largemomentum effective theory (LaMET) approach. The
momentumsmearing technique has been implemented to improve signals at large
meson momenta. We subtract the power divergence due to Wilson line to high
precision using multiple lattice spacings. The kaon structure clearly shows an
asymmetry of the distribution amplitude around $x=1/2$, a clear sign of its
skewness. We also study the leading SU(3) flavor symmetry breaking relations
for the pion, kaon and eta meson distribution amplitudes, and the results are
consistent with the prediction from chiral perturbation theory.

In the framework of quantum chromodynamics (QCD), parton distribution
functions (PDFs) quantify how the momentum and spin of a hadron are divided
among its quark and gluon constituents. Two main approaches exist to determine
PDFs. The first approach, based on QCD factorization theorems, realizes a QCD
analysis of a suitable set of hardscattering measurements, often using a
variety of hadronic observables. The second approach, based on firstprinciple
operator definitions of PDFs, uses lattice QCD to compute directly some
PDFrelated quantities, such as their moments. Motivated by recent progress in
both approaches, in this document we present an overview of latticeQCD and
globalanalysis techniques used to determine unpolarized and polarized proton
PDFs and their moments. We provide benchmark numbers to validate present and
future latticeQCD calculations and we illustrate how they could be used to
reduce the PDF uncertainties in current unpolarized and polarized global
analyses. This document represents a first step towards establishing a common
language between the two communities, to foster dialogue and to further improve
our knowledge of PDFs.

We present lattice results for the isovector unpolarized parton distribution
with nonperturbative RI/MOMscheme renormalization on the lattice. In the
framework of largemomentum effective field theory (LaMET), the full
Bjorken$x$ dependence of a momentumdependent quasidistribution is calculated
on the lattice and matched to the ordinary lightcone parton distribution at
oneloop order, with power corrections included. The important step of RI/MOM
renormalization that connects the lattice and continuum matrix elements is
detailed in this paper. A few consequences of the results are also addressed
here.

We show that the empirical linear relation between the magnitude of the EMC
effect in deep inelastic scattering on nuclei and the short range correlation
scaling factor $a_2$ extracted from highenergy quasielastic scattering at
$x\ge 1$ is a natural consequence of scale separation and derive the
relationship using effective field theory. While the scaling factor $a_2$ is a
ratio of nuclear matrix elements that individually depend on the calculational
scheme, we show that the ratio is independent of this choice. We perform
Green's function Monte Carlo calculations with both chiral and ArgonneUrbana
potentials to verify this and determine the scaling factors for light nuclei.
The resulting values for $^3$He and $^4$He are in good agreement with
experimental values. We also present results for $^9$Be and $^{12}$C extracted
from variational Monte Carlo calculations.

The renormalization of entanglement entropy of quantum field theories is
investigated in the simplest setting with a $\lambda \phi^4$ scalar field
theory. The 3+1 dimensional spacetime is separated into two regions by an
infinitely flat 2dimensional interface. The entanglement entropy of the system
across the interface has an elegant geometrical interpretation using the
replica trick, which requires putting the field theory on a curved spacetime
background. We demonstrate that the theory, and hence the entanglement entropy,
is renormalizable at order $\lambda$ once all the relevant operators up to
dimension4 are included in the action. This exercise has a onetoone
correspondence to entanglement entropy interpretation of the black hole entropy
which suggests that our treatment is sensible. Our study suggests that
entanglement entropy is renormalizable and is a physical quantity.

We present the first lattice results on isovector unpolarized and
longitudinally polarized parton distribution functions (PDFs) at physical pion
mass. The PDFs are obtained using the largemomentum effective field theory
(LaMET) framework where the full Bjorken$x$ dependence of finitemomentum
PDFs, called "quasiPDFs", can be calculated on the lattice. The quasiPDF
nucleon matrix elements are renormalized nonperturbatively in RI/MOMscheme.
However, the recent renormalized quasiPDFs suffer from unphysical oscillations
that alter the shape of the true distribution as a function of Bjorken$x$. In
this paper, we propose two possible solutions to overcome this problem, and
demonstrate the efficacy of the methods on the 2+1+1flavor lattice data at
physical pion mass with lattice spacing 0.09~fm and volume $(5.76\mbox{
fm})^3$.

We show that universality near a critical end point implies a characteristic
relation between third and fourthorder baryon susceptibilities $\chi_3$ and
$\chi_4$, resulting in a bananashaped loop when $\chi_4$ is plotted as a
function of $\chi_3$ along a freezeout line. Including the individual
enhancements of $\chi_3$ and $\chi_4$ near a critical point, these features may
be a consistent set of observations supporting the interpretation of baryon
fluctuations data as arising from criticality.

We present the first latticeQCD calculation of the pion distribution
amplitude using the large momentum effective field theory (LaMET) approach,
which allows us to extract lightcone parton observables from a Euclidean
lattice. The mass corrections needed to extract the pion distribution amplitude
from this approach are calculated to all orders in $m^2_{\pi} /P_z^2$. We also
implement the Wilson line renormalization which is crucial to remove the power
divergences in this approach, and find that it reduces the oscillation at the
end points of the distribution amplitude. Our exploratory result at 310MeV
pion mass favors a singlehump form broader than the asymptotic form of the
pion distribution amplitude.

Lowenergy electronic recoil caused by solar neutrinos in multiton xenon
detectors is an important subject not only because it is a source of the
irreducible background for direct searches of weaklyinteracting massive
particles (WIMPs), but also because it provides a viable way to measure the
solar $pp$ and $^{7}\textrm{Be}$ neutrinos at the precision level of current
standard solar model predictions. In this work we perform $\textit{ab initio}$
manybody calculations for the structure, photoionization, and
neutrinoionization of xenon. It is found that the atomic binding effect yields
a sizable suppression to the neutrinoelectron scattering cross section at low
recoil energies. Compared with the previous calculation based on the free
electron picture, our calculated event rate of electronic recoil in the same
detector configuration is reduced by about $25\%$. We present in this paper the
electronic recoil rate spectrum in the energy window of 100 eV  30 keV with
the standard per ton per year normalization for xenon detectors, and discuss
its implication for low energy solar neutrino detection (as the signal) and
WIMP search (as a source of background).

Recent developments showed that hadron lightcone parton distributions could
be directly extracted from spacelike correlators, known as quasi parton
distributions, in the large hadron momentum limit. Unlike the normal lightcone
parton distribution, a quasi parton distribution contains ultraviolet (UV)
power divergence associated with the Wilson line self energy. We show that to
all orders in the coupling expansion, the power divergence can be removed by a
"mass" counterterm in the auxiliary $z$field formalism, in the same way as the
renormalization of power divergence for an open Wilson line. After adding this
counterterm, the quasi quark distribution is improved such that it contains at
most logarithmic divergences. Based on a simple version of discretized gauge
action, we present the oneloop matching kernel between the improved
nonsinglet quasi quark distribution with a lattice regulator and the
corresponding quark distribution in dimensional regularization.

We study how to engineer holographic models with features of a high
temperature superconductor phase diagram. We introduce a field in the bulk
which provides a tunable "doping" parameter in the boundary theory. By
designing how this field changes the effective masses of other order parameter
fields, desired phase diagrams can be engineered. We give examples of
generating phase diagrams with phase boundaries similar to a superconducting
dome and an antiferromagnetic phase by including two order parameter fields.
We also explore whether the pseudo gap phase can be described without adding
another order parameter field and discuss the potential scaling symmetry
associated with a quantum critical point hidden under the superconducting dome
in this phase diagram.

We present the first latticeQCD calculation of the isovector polarized
parton distribution functions (both helicity and transversity) using the
largemomentum effective field theory (LaMET) approach for direct Bjorken$x$
dependence. We first review the detailed steps of the procedure in the
unpolarized case, then generalize to the helicity and transversity cases. We
also derive a new masscorrection formulation for all three cases. We then
compare the effects of each finitemomentum correction using lattice data
calculated at $M_\pi\approx 310$ MeV. Finally, we discuss the implications of
these results for the poorly known antiquark structure and predict the
seaflavor asymmetry in the transversely polarized nucleon.

We study the nonlinear responses of relativistic chiral matter to the
external fields, such as the electric field ${\bf E}$, gradients of temperature
and chemical potential, ${\bf \nabla} T$ and ${\bf \nabla} \mu$. Using the
kinetic theory with Berry curvature corrections under the relaxation time
approximation, we compute the transport coefficients of possible new electric
currents that are forbidden in usual chirally symmetric matter, but are allowed
in chirally asymmetric matter by parity. In particular, we find a new type of
electric current proportional to ${\bf \nabla} \mu \times {\bf E}$ due to the
interplay between the effects of the Berry curvature and collisions. We also
derive an analogue of the "WiedemannFranz" law specific for anomalous
nonlinear transport in relativistic chiral matter.

The transition magnetic moment of a steriletoactive neutrino conversion
gives rise to not only radiative decay of a sterile neutrino, but also its
nonstandard interaction (NSI) with matter. For sterile neutrinos of keVmass
as dark matter candidates, their decay signals are actively searched for in
cosmic Xray spectra. In this work, we consider the NSI that leads to atomic
ionization, which can be detected by direct dark matter experiments. It is
found that this inelastic scattering process for a nonrelativistic sterile
neutrino has a pronounced enhancement in the differential cross section at
energy transfer about half of its mass, manifesting experimentally as peaks in
the measurable energy spectra. The enhancement effects gradually smear out as
the sterile neutrino becomes relativistic. Using data taken with germanium
detectors that have fine energy resolution in keV and subkeV regimes,
constraints on sterile neutrino mass and its transition magnetic moment are
derived and compared with those from astrophysical observations.

We evaluate the third and fourthorder baryon, charge and strangeness
susceptibilities near a chiral critical point using the NambuJonaLasinio
model. We identify robust qualitative behaviours of the susceptibilities along
hypothetical freezeout lines that agree with previous model studies.
Quantitatively, baryon number fluctuations are the largest in magnitude and
thus offer the strongest signal when freezeout occurs farther away from a
critical point. Charge and strangeness susceptibilities also diverge at a
critical point, but the area where the divergence dominates is smaller, meaning
freezeout must occur closer to a critical point for a signal to be visible in
these observables. In case of strangeness, this is attributable to the
relatively large strange quark mass. Plotting the third and fourthorder
susceptibilities against each other along the freezeout line exhibits clearly
their nonmontonicity and robust features.

Scattering of subGeV dark matter (DM) particles with hydrogen atoms is
studied in this paper. The interactions of DM with electrons and nucleons are
both included and formulated in a general framework based on nonrelativistic
effective field theory. On the assumption of same dark matter coupling
strengths, it is found that DMelectron interactions dominate the inelastic
atomic transitions to discrete excited states and ionization continuum around
the threshold regions, and DMnucleon interactions become more important with
increasing energy and dominate in elastic scattering. The conclusion should
apply, qualitatively, to practical detector species so that electronic and
nuclear contributions in DM scattering processes can be disentangled, while
issues including binding effects and recoil mechanism in manybody systems will
require further detailed calculations.

We calculate the entanglement entropy of the $SU(N)$ YangMills gauge
theories on the lattice under the strong coupling expansion in powers of
$\beta=2N/g^{2}$, where $g$ is the coupling constant. Using the replica method,
our Lagrangian formalism maintains gauge invariance on the lattice. At
$O(\beta^{2})$ and $O(\beta^{3})$, the entanglement entropy is solely
contributed by the central plaquettes enclosing the conical singularity of the
$n$sheeted Riemann surface. The area law emerges naturally to the highest
order $O(\beta^{3})$ of our calculation. The leading $O(\beta)$ term is
negative, which could in principle be canceled by taking into account the
"cosmological constant" living in interface of the two entangled subregions.
This unknown cosmological constant resembles the ambiguity of edge modes in the
Hamiltonian formalism. We further speculate this unknown cosmological constant
can show up in the entanglement entropy of scalar and spinor field theories as
well. Furthermore, it could play the role of a counterterm to absorb the
ultraviolet divergence of entanglement entropy and make entanglement entropy a
finite physical quantity.

We present the first direct lattice calculation of the isovector seaquark
parton distributions using the formalism developed recently by one of the
authors. We use $N_f=2+1+1$ HISQ lattice gauge ensembles (generated by MILC
Collaboration) and clover valence fermions with pion mass 310 MeV. We are able
to obtain the qualitative features of the nucleon sea flavor structure even at
this large pion mass: We observe violation of the Gottfried sum rule,
indicating $\overline{d}(x) > \overline{u}(x)$; the helicity distribution obeys
$\Delta \overline{u}(x) > \Delta \overline{d}(x)$, which is consistent with the
STAR data at large and small leptonic pseudorapidity.

We calculate modelindependently the impact of the critical point on higher
order baryon susceptibilities $\chi_n$, showing how they depend on fluctuations
of the order parameter. Including all tree level diagrams, we find new
contributions to $\chi_4$ equally important to the kurtosis of the order
parameter fluctuations, and we characterize the kurtosis and other nonguassian
moments as functions on the phase diagram. Important features of this analysis
are then confirmed by a GrossNeveu model study with good agreement with other
model studies as well as lattice and experimental data. This suggests the
universality of the characteristic peak in baryon susceptibilities as a signal
of the critical point. We discuss leveraging measurements of different $\chi_n$
to extrapolate the location of the critical point.