
As manybody Floquet theory becomes more popular, it is important to find
ways to connect theory with experiment. Theoretical calculations can have a
periodic driving field that is always on, but experiment cannot. Hence, we need
to know how long a driving field is needed before the system starts to look
like the periodically driven Floquet system. We answer this question here for
noninteracting band electrons in the infinitedimensional limit by studying the
properties of the system under pulsed driving fields and illustrating how they
approach the Floquet limit. Our focus is on determining the minimal pulse
lengths needed to recover the qualitative and semiquantitative Floquet theory
results.

In a large transverse field, there is an energy cost associated with flipping
spins along the axis of the field. This penalty can be employed to relate the
transversefield Ising model in a large field to the XY model in no field (when
measurements are performed at the right stroboscopic times). We describe the
details for how this relationship works and, in particular, we also show under
what circumstances it fails. We examine wavefunction overlap between the two
models and observables, such as spinspin Green's functions. In general, the
mapping is quite robust at short times, but will ultimately fail if the run
time becomes too long. There is also a tradeoff between the length of time one
can run a simulation out to and the time jitter of the stroboscopic
measurements that must be balanced when planning to employ this mapping.

We present a theory for charge and heat transport parallel to the interfaces
of a multilayer (ML) in which the interfacing gives rise the redistribution of
the electronic charges. The ensuing electrical field couples selfconsistently
to the itinerant electrons, so that the properties of the ML crucially depend
on an interplay between the onsite Coulomb forces and the long range
electrostatic forces. The ML is described by the FalicovKimball model and the
selfconsistent solution is obtained by iterating simultaneously the DMFT and
the Poisson equations. This yields the reconstructed charge profile, the
electrical potential, the planar density of states, the transport function, and
the transport coefficients of the device.
We find that a heterostructure built of two MottHubbard insulators exhibits,
in a large temperature interval, a linear conductivity and a large
temperatureindependent thermopower. The charge and energy currents are
confined to the central part of the ML. Our results indicate that correlated
multilayers have the potential for applications; by tuning the band shift and
the Coulomb correlation on the central planes, we can bring the chemical
potential in the immediate proximity of the MottHubbard gap edge and optimize
the transport properties of the device. In such a heterostructure, a small gate
voltage can easily induce a MI transition. This switching does not involve the
diffusion of electrons over macroscopic distances and it is much faster than in
ordinary semiconductors. Furthermore, the right combination of strongly
correlated materials with small ZT can produce, theoretically at least, a
heterostructure with a large ZT.

One of the challenges for fermionic cold atom experiments in optical lattices
is to cool the systems to low enough temperature that they can form quantum
degenerate ordered phases. In particular, there has been significant work in
trying to find the antiferromagnetic phase transition of the Hubbard model in
three dimensions, without success. Here, we attack this problem from a
different angle by enhancing the ordering temperature via an increase in the
degeneracy of the atomic species trapped in the optical lattice. In addition to
developing the general theory, we also discuss some potential systems where one
might be able to achieve these results experimentally.

We examine the exact equation of motion for the relaxation of populations of
strongly correlated electrons after a nonequilibrium excitation by a pulsed
field, and prove that the populations do not change when the Green's functions
have no average time dependence. We show how the average time dependence enters
into the equation of motion to lowest order and describe what governs the
relaxation process of the electron populations in the longtime limit. While
this result may appear, on the surface, to be required by any steadystate
solution, the proof is nontrivial, and provides new critical insight into how
nonequilibrium populations relax, which goes beyond the assumption that they
thermalize via a simple relaxation rate determined by the imaginary part of the
selfenergy, or that they can be described by a quasiequilibrium condition
with a FermiDirac distribution and a timedependent temperature. We also
discuss the implications of this result to approximate theories, which may not
satisfy the exact relation in the equation of motion.

We examine the time evolution of an asymmetric Hubbard dimer, which has a
different onsite interaction on the two sites. The Hamiltonian has a
timedependent hopping term, which can be employed to describe an electric
field (which creates a Hamiltonian with complex matrix elements), or it can
describe a modulation of the lattice (which has real matrix elements). By
examining the symmetries under spin and pseudospin, we show that the former
case involves at most a 3 x 3 blockit can be mapped onto the time evolution
of a timeindependent Hamiltonian, so the dynamics can be evaluated
analytically and exactly (by solving a nontrivial cubic equation). We also show
that the latter case reduces to at most 2 x 2 blocks, and hence the time
evolution for a single Trotter step can be determined exactly, but the time
evolution generically requires a Trotter product.

In this review, we develop the formalism employed to describe
chargedensitywave insulators in pump/probe experiments using ultra short
driving pulses of light. The theory emphasizes exact results in the simplest
model for a chargedensity wave insulator (given by a noninteracting systems
with two bands and a gap) and by employing nonequilibrium dynamical meanfield
theory to solve the FalicovKimball model in its ordered phase. We show both
how to develop the formalism and how the solutions behave. Care is taken to
describe the details behind these calculations and to show how to verify their
accuracy via sumrule constraints.

We describe different approximations associated with employing a constant
matrix element for the coupling of light to multiband electrons in the context
of timeresolved angleresolved photoemission spectroscopy (TRARPES). In
particular, we demonstrate that the constant matrix approximation only holds
for one of the coordinate representations, and changing to other bases,
requires including nonconstant corrections to the matrix element. We also
discuss some simplifying approximations, where a constant matrix element is
employed in multiple bases, and the consequences of this further approximation
(especially with respect to the TRARPES signal no longer being nonnegative).
We also discuss issues related to gauge invariance of the final spectra.

Recent developments in the techniques of ultrafast pumpprobe photoemission
have made possible the search for collective modes in strongly correlated
systems out of equilibrium. Including inelastic scattering processes and a
retarded interaction, we simulate time and angle resolved photoemission
spectroscopy (trARPES) to study the amplitude mode of a dwave superconductor,
a collective mode excited through the nonlinear lightmatter coupling to the
pump pulse. We find that the amplitude mode oscillations of the dwave order
parameter occur in phase at a single frequency that is twice the
quasisteadystate maximum gap size after pumping. We comment on the necessary
conditions for detecting the amplitude mode in trARPES experiments.

We review recent work on the theory for pump/probe photoemission spectroscopy
of electronphonon mediated superconductors in both the normal and the
superconducting states. We describe the formal developments that allow one to
solve the MigdalEliashberg theory in nonequilibrium for an ultrashort laser
pumping field, and explore the solutions which illustrate the relaxation as
energy is transferred from electrons to phonons. We focus on exact results
emanating from sum rules and approximate numerical results which describe rules
of thumb for relaxation processes. In addition, in the superconducting state,
we describe how Higg's oscillations can be excited due to the nonlinear
coupling with the electric field and how pumping the system can enhance
superconductivity.

We determine the exact timeresolved photoemission spectroscopy for a nesting
driven chargedensitywave (described by the spinless FalicovKimball model
within dynamical meanfield theory). The pumpprobe experiment involves two
light pulses: the first is an ultrashort intense pump pulse that excites the
system into nonequilibrium, and the second is a lower amplitude higher
frequency probe pulse that photoexcites electrons. We examine three different
cases: the strongly correlated metal, the quantumcritical charge density wave
and the critical Mott insulator. Our results show that the quantum critical
charge density wave has an ultra efficient relaxation channel that allows
electrons to be deexcited during the pump pulse, resulting in little net
excitation. In contrast, the metal and the Mott insulator show excitations that
are closer to what one expects from these systems. In addition, the pump field
produces spectral band narrowing, peak sharpening, and a spectral gap
reduction, all of which rapidly return to their field free values after the
pump is over.

We develop a generalized gradient expansion of the inhomogeneous dynamical
meanfield theory method for determining properties of ultracold atoms in a
trap. This approach goes beyond the wellknown local density approximation and
at higher temperatures, in the normal phase, it shows why the local density
approximation works so well, since the local density and generalized gradient
approximations are essentially indistinguishable from each other (and from the
exact solution within full inhomogeneous dynamical meanfield theory). But
because the generalized gradient expansion only involves nearestneighbor
corrections, it does not work as well at low temperatures, when the systems
enter into ordered phases. This is primarily due to the problem that ordered
phases often satisfy some global constraints which determine the spatial
ordering pattern, and the local density and generalized gradient approximations
are not able to impose those kinds of constraints; they also overestimate the
tendency to order. The theory is applied to phase separation of different mass
fermionic mixtures represented by the FalicovKimball model and to determining
the entropy per particle of a fermionic system represented by the Hubbard
model. The generalized gradient approximation is a useful diagnostic for the
accuracy of the local density approximationwhen both methods agree, they are
likely accurate, when they disagree, neither is likely to be correct.

In complex materials various interactions play important roles in determining
the material properties. Angle Resolved Photoelectron Spectroscopy (ARPES) has
been used to study these processes by resolving the complex single particle
self energy $\Sigma(E)$ and quantifying how quantum interactions modify bare
electronic states. However, ambiguities in the measurement of the real part of
the self energy and an intrinsic inability to disentangle various contributions
to the imaginary part of the self energy often leave the implications of such
measurements open to debate. Here we employ a combined theoretical and
experimental treatment of femtosecond timeresolved ARPES (trARPES) and show
how measuring the population dynamics using trARPES can be used to separate
electronboson interactions from electronelectron interactions. We demonstrate
the analysis of a welldefined electronboson interaction in the unoccupied
spectrum of the cuprate Bi$_{2}$Sr$_{2}$CaCu$_{2}$O$_{8+x}$ characterized by an
excited population decay time constant $\tau_{QP}$ that maps directly to a
discrete component of the equilibrium self energy not readily isolated by
static ARPES experiments.

Recent work proposed a variant on Ramsey interferometry for coupled
spin$1/2$ systems that directly measures the retarded spinspin Green's
function. We expand on that work by investigating nonequilibrium retarded
spinspin Green's functions within the transversefield Ising model. We derive
the lowest four spectral moments to understand the shorttime behavior and we
employ a Lehmannlike representation to determine the spectral behavior. We
simulate a Ramsey protocol for a nonequilibrium quantum spin system that
consists of a coherent superposition of the ground state and diabatically
excited higherenergy states via a temporally ramped transverse magnetic field.
We then apply the Ramsey spectroscopy protocol to the final Hamiltonian, which
has a constant transverse field. The shorttime behavior directly relates to
LiebRobinson bounds for the transport of manybody correlations, while the
longtime behavior relates to the excitation spectra of the Hamiltonian.
Compressive sensing is employed in the data analysis to efficiently extract
that spectra.

We model the bangbang optimization protocol as a shortcut to adiabaticity in
the groundstate preparation of an iontrapbased quantum simulator. Compared
to a locally adiabatic evolution, the bangbang protocol produces a somewhat
lower groundstate probability, but its implementation is so much simpler than
the locally adiabatic approach, that it remains an excellent choice to use for
maximizing groundstate preparation in systems that cannot be solved with
conventional computers. We describe how one can optimize the shortcut and
provide specific details for how it can be implemented with current
iontrapbased quantum simulators.

Using the KadanoffBaymKeldysh formalism, we employ nonequilibrium dynamical
meanfield theory to exactly solve for the nonlinear response of an
electronmediated chargedensitywaveordered material. We examine both the dc
current and the order parameter of the conduction electrons as the ordered
system is driven by the electric field. Although the formalism we develop
applies to all models, for concreteness, we examine the chargedensitywave
phase of the FalicovKimball model, which displays a number of anomalous
behaviors including the appearance of subgap density of states as the
temperature increases. These subgap states should have a significant impact on
transport properties, particularly the nonlinear response of the system to a
large dc electric field.

The goal of adiabatic groundstate preparation is to start a simple quantum
system in its ground state and adiabatically evolve the Hamiltonian to a
complex one, maintaining the ground state throughout the evolution. In
iontrapbased quantum simulations, coherence times are too short to allow for
adiabatic evolution for large chains, so the system evolves diabatically,
creating excitations to higher energy states. Because the probability for
diabatic excitation depends exponentially on the excitation energy and because
the thermal distribution also depends exponentially on the excitation energy,
we investigate whether the diabatic excitation can create a thermal
distribution; as this could serve as an alternative for creating thermal states
of complex quantum systems without requiring contact with a heat bath. In this
work, we explore this relationship and determine situations where diabatic
excitation can approximately create such a thermal state.

Theoretical models of spins coupled to bosons provide a simple setting for
studying a broad range of important phenomena in manybody physics, from
virtually mediated interactions to decoherence and thermalization. In many
atomic, molecular, and optical systems, such models also underlie the most
successful attempts to engineer strong, longranged interactions for the
purpose of entanglement generation. Especially when the coupling between the
spins and bosons is strongsuch that it cannot be treated
perturbativelythe properties of such models are extremely challenging to
calculate theoretically. Here, exact analytical expressions for nonequilibrium
spinspin correlation functions are derived for a specific model of spins
coupled to bosons. The spatial structure of the coupling between spins and
bosons is completely arbitrary, and thus the solution can be applied to systems
in any number of dimensions. The explicit and nonperturbative inclusion of the
bosons enables the study of entanglement generation (in the form of spin
squeezing) even when the bosons are driven strongly and nearresonantly, and
thus provides a quantitative view of the breakdown of adiabatic elimination
that inevitably occurs as one pushes towards the fastest entanglement
generation possible. The solution also helps elucidate the effect of finite
temperature on spin squeezing. The model considered is relevant to a variety of
atomic, molecular, and optical systems, such as atoms in cavities or trapped
ions. As an explicit example, the results are used to quantify phonon effects
in trapped ion quantum simulators, which are expected to become increasingly
important as these experiments push towards larger numbers of ions.

One of the goals in quantum simulation is to adiabatically generate the
ground state of a complicated Hamiltonian by starting with the ground state of
a simple Hamiltonian and slowly evolving the system to the complicated one. If
the evolution is adiabatic and the initial and final ground states are
connected due to having the same symmetry, then the simulation will be
successful. But in most experiments, adiabatic simulation is not possible
because it would take too long, and the system has some level of diabatic
excitation. In this work, we quantify the extent of the diabatic excitation
even if we do not know {\it a priori} what the complicated ground state is.
Since many quantum simulator platforms, like trapped ions, can measure the
probabilities to be in a product state, we describe techniques that can employ
these measurements to estimate the probability of being in the ground state of
the system after the diabatic evolution. These techniques do not require one to
know any properties about the Hamiltonian itself, nor to calculate its
eigenstate properties. All the information is derived by analyzing the
productstate measurements as functions of time.

One of the challenges with quantum simulation in ion traps is that the
effective spinspin exchange couplings are not uniform across the lattice. This
can be particularly important in Penning trap realizations where the presence
of an ellipsoidal boundary at the edge of the trap leads to dislocations in the
crystal. By adding an additional anharmonic potential to better control
interion spacing, and a triangular shaped rotating wall potential to reduce the
appearance of dislocations, one can achieve better uniformity of the ionic
positions. In this work, we calculate the axial phonon frequencies and the
spinspin interactions driven by a spindependent optical dipole force, and
discuss what effects the more uniform ion spacing has on the spin simulation
properties of Penning trap quantum simulators. Indeed, we find the spinspin
interactions behave more like a power law for a wide range of parameters.

Generating big data pervades much of physics. But some problems, which we
call extreme data problems, are too large to be treated within big data
science. The nonequilibrium quantum manybody problem on a lattice is just such
a problem, where the Hilbert space grows exponentially with system size and
rapidly becomes too large to fit on any computer (and can be effectively
thought of as an infinitesized data set). Nevertheless, much progress has been
made with computational methods on this problem, which serve as a paradigm for
how one can approach and attack extreme data problems. In addition, viewing
these physics problems from a computerscience perspective leads to new
approaches that can be tried to solve them more accurately and for longer
times. We review a number of these different ideas here.

The control of physical properties of solids with short laser pulses is an
intriguing prospect of ultrafast materials science. Continuouswave
highfrequency laser driving with circular polarization was predicted to induce
a lightmatter coupled new state possessing a quasistatic band structure with
an energy gap and a quantum Hall effect, coined "Floquet topological
insulator". Whereas the envisioned Floquet topological insulator requires well
separated Floquet bands and therefore highfrequency pumping, a natural
followup question regards the creation of Floquetlike states in graphene with
realistic pump laser pulses. Here we predict that with short lowfrequency
laser pulses attainable in pumpprobe experiments, states with local spectral
gaps at the Dirac points and novel pseudospin textures can be achieved in
graphene using circular light polarization. We demonstrate that time and
angleresolved photoemission spectroscopy can track these states by measuring
sizeable energy gaps and quasiFloquet energy bands that form on femtosecond
time scales. By analyzing Floquet energy level crossings and snapshots of
pseudospin textures near the Dirac points, we identify transitions to new
states with optically induced nontrivial changes of sublattice mixing that can
lead to Berry curvature corrections of electrical transport and magnetization.

Ultracold mixtures of different atomic species have great promise for
realizing novel manybody phenomena. In a binary mixture of femions with a
large mass difference and repulsive interspecies interactions, a disordered
Mott insulator phase can occur. This phase displays an incompressible total
density, although the relative density remain compressible. We use
strongcoupling and Monte Carlo calculations to show that this phase exists for
a broad parameter region for ultracold gases confined in a harmonic trap on a
threedimensional optical lattice, for experimentally accessible values of the
trap parameters.

The Coulomb repulsion between ions in a linear Paul trap give rise to
anharmonic terms in the potential energy when expanded about the equilibrium
positions. We examine the effect of these anharmonic terms on the accuracy of a
quantum simulator made from trapped ions. To be concrete, we consider a linear
chain of $\text{Yb}^{171+}$ ions stabilized close to the zigzag transition. We
find that for typical experimental temperatures, frequencies change by no more
than a factor of $0.01\%$ due to the anharmonic couplings. Furthermore, shifts
in the effective spinspin interactions (driven by a spindependent optical
dipole force) also tend to be small for detunings to the blue of the transverse
centerofmass frequency. However, detuning the spin interactions near other
frequencies can lead to nonnegligible anharmonic contributions to the effective
spinspin interactions. We also examine an odd behavior exhibited by the
harmonic spinspin interactions for a range of intermediate detunings, where
nearest neighbor spins with a larger spatial separation on the ion chain
interact more strongly than nearest neighbors with a smaller spatial
separation.

We develop the theory to describe the equilibrium ion positions and phonon
modes for a trapped ion quantum simulator in an oblate Paul trap that creates
twodimensional Coulomb crystals in a triangular lattice. By coupling the
internal states of the ions to laser beams propagating along the symmetry axis,
we study the effective Ising spinspin interactions that are mediated via the
axial phonons and are less sensitive to ion micromotion. We find that the axial
mode frequencies permit the programming of Ising interactions with inverse
power law spinspin couplings that can be tuned from uniform to $r^{3}$ with
DC voltages. Such a trap could allow for interesting new geometrical
configurations for quantum simulations on moderately sized systems including
frustrated magnetism on triangular lattices or AharonovBohm effects on ion
tunneling. The trap also incorporates periodic boundary conditions around loops
which could be employed to examine time crystals.