
The phenomenon of manybody localised (MBL) systems has attracted significant
interest in recent years, for its intriguing implications from a perspective of
both condensedmatter and statistical physics: they are insulators even at
nonzero temperature and fail to thermalise, violating expectations from
quantum statistical mechanics. What is more, recent seminal experimental
developments with ultracold atoms in optical lattices constituting analog
quantum simulators have pushed manybody localised systems into the realm of
physical systems that can be measured with high accuracy. In this work, we
introduce experimentally accessible witnesses that directly probe distinct
features of MBL, distinguishing it from its Anderson counterpart. We insist on
building our toolbox from techniques available in the laboratory, including
onsite addressing, superlattices, and timeofflight measurements,
identifying witnesses based on fluctuations, densitydensity correlators,
densities, and entanglement. We build upon the theory of out of equilibrium
quantum systems, in conjunction with tensor network and exact simulations,
showing the effectiveness of the tools for realistic models.

We develop a hybrid semiclassical method to study the time evolution of one
dimensional quantum systems in and out of equilibrium. Our method handles
internal degrees of freedom completely quantum mechanically by a modified time
evolving block decimation method, while treating orbital quasiparticle motion
classically. We can follow dynamics up to timescales well beyond the reach of
standard numerical methods to observe the crossover between preequilibrated
and locally phase equilibrated states. As an application, we investigate the
quench dynamics and phase fluctuations of a pair of tunnel coupled one
dimensional Bose condensates. We demonstrate the emergence of solitoncollision
induced phase propagation, solitonentropy production and multistep
thermalization. Our method can be applied to a wide range of gapped
onedimensional systems.

In contrast to lattice systems where powerful numerical techniques such as
matrix product state based methods are available to study the nonequilibrium
dynamics, the nonequilibrium behaviour of continuum systems is much harder to
simulate. We demonstrate here that Hamiltonian truncation methods can be
efficiently applied to this problem, by studying the quantum quench dynamics of
the 1+1 dimensional Ising field theory using a truncated free fermionic space
approach. After benchmarking the method with integrable quenches corresponding
to changing the mass in a free Majorana fermion field theory, we study the
effect of an integrability breaking perturbation by the longitudinal magnetic
field. In both the ferromagnetic and paramagnetic phases of the model we find
persistent oscillations with frequencies set by the lowlying particle
excitations not only for small, but even for moderate size quenches. In the
ferromagnetic phase these particles are the various nonperturbative confined
bound states of the domain wall excitations, while in the paramagnetic phase
the single magnon excitation governs the dynamics, allowing us to capture the
time evolution of the magnetisation using a combination of known results from
perturbation theory and form factor based methods. We point out that the
dominance of low lying excitations allows for the numerical or experimental
determination of the mass spectra through the study of the quench dynamics.

We present a model for an autonomous quantum thermal machine comprised of two
qubits capable of manipulating and even amplifying the local coherence in a
nondegenerate external system. The machine uses only thermal resources,
namely, contact with two heat baths at different temperatures, and the external
system has a nonzero initial amount of coherence. The method we propose allows
for an interconversion between energy, both work and heat, and coherence in an
autonomous configuration working in outofequilibrium conditions. This model
raises interesting questions about the role of fundamental limitations on
transformations involving coherence and opens up new possibilities in the
manipulation of coherence by autonomous thermal machines.

Fixed points for scalar theories in $4\varepsilon$, $6\varepsilon$ and
$3\varepsilon$ dimensions are discussed. It is shown how a large range of
known fixed points for the four dimensional case can be obtained by using a
general framework with two couplings. The original maximal symmetry, $O(N)$, is
broken to various subgroups, both discrete and continuous. A similar discussion
is applied to the six dimensional case. Perturbative applications of the
$a$theorem are used to help classify potential fixed points. At lowest order
in the $\varepsilon$expansion it is shown that at fixed points there is a
lower bound for $a$ which is saturated at bifurcation points.

Rigorous nonequilibrium actions for the manybody problem are usually derived
by means of path integrals combined with a discrete temporal mesh on the
SchwingerKeldysh time contour. The latter suffers from a fundamental
limitation: the initial state on this contour cannot be arbitrary, but
necessarily needs to be described by a noninteracting density matrix, while
interactions are switched on adiabatically. The KostantinovPerel' contour
overcomes these and other limitations, allowing generic initialstate
preparations. In this Article, we apply the technique of the discrete temporal
mesh to rigorously build the nonequilibrium path integral on the
KostantinovPerel' time contour.

Dynamic renormalization group (RG) of fluctuating viscoelastic equations is
investigated to clarify the cause for numerically reported disappearance of
anomalous heat conduction (recovery of Fourier's law) in lowdimensional
momentumconserving systems. RG flow is obtained explicitly for simplified two
model cases: a onedimensional continuous medium under low pressure and
incompressible viscoelastic medium of arbitrary dimensions. Analyses of these
clarify that the inviscid fixed point of contributing the anomalous heat
conduction becomes unstable under the RG flow of nonzero elasticwave speeds.
The dynamic RG analysis further predicts a universal scaling of describing the
crossover between the growth and saturation of observed heat conductivity,
which is confirmed through the numerical experiments of FermiPastaUlam
$\beta$ (FPU$\beta$) lattices.

Fekete's lemma shows the existence of limits in subadditive sequences. This
lemma, and generalisations of it, also have been used to prove the existence of
thermodynamic limits in statistical mechanics. In this paper it is shown that
the two variable supermultiplicative relation $$ p_{n_1}(m_1)\,p_{n_2}(m_2)
\leq p_{n_1+n_2}(m_1 + m_2) $$ together with mild assumptions, imply the
existence of the limit $$ \log \mathcal{P}_\#(\epsilon) = \lim_{n\to\infty}
\frac{1}{n} \log p_n(\lfloor \epsilon n \rfloor).$$ This is a generalisation of
Fekete's lemma. The existence of $\mathcal{P}_\#(\epsilon)$ are proven for
models of adsorbing walks and polygons, and for pulled polygons. In addition,
numerical data are presented estimating the general shape of $\log
\mathcal{P}_\# (\epsilon)$ of models of square lattice selfavoiding walks and
polygons.

The BornMarkov approximation is widely used to study dynamics of open
quantum systems coupled to external baths. Using Keldysh formalism, we show
that the dynamics of a system of bosons (fermions) linearly coupled to
noninteracting bosonic (fermionic) bath falls outside this paradigm if the
bath spectral function has nonanalyticities as a function of frequency. In
this case, we show that the dissipative and noise kernels governing the
dynamics have distinct power law tails. The Green's functions show a short time
"quasi" Markovian exponential decay before crossing over to a power law tail
governed by the nonanalyticity of the spectral function. We study a system of
bosons (fermions) hopping on a one dimensional lattice, where each site is
coupled linearly to an independent bath of noninteracting bosons (fermions).
We obtain exact expressions for the Green's functions of this system which show
power law decay $\sim tt'^{3/2}$. We use these to calculate density and
current profile, as well as unequal time currentcurrent correlators. While the
density and current profiles show interesting quantitative deviations from
Markovian results, the currentcurrent correlators show qualitatively distinct
long time power law tails $tt'^{3}$ characteristic of nonMarkovian
dynamics. We show that the power law decays survive in presence of
interparticle interaction in the system, but the crossover time scale is
shifted to larger values with increasing interaction strength.

We study the longtime behavior of conservative interacting particle systems
in $Z$: the activated random walk model for reactiondiffusion systems and the
stochastic sandpile. We prove that both systems undergo an absorbingstate
phase transition.

Identifying the physical basis of heterosis (or hybrid vigor) has remained
elusive despite over a hundred years of research on the subject. The three main
theories of heterosis are dominance theory, overdominance theory, and epistasis
theory. Kacser and Burns (1981) identified the molecular basis of dominance,
which has greatly enhanced our understanding of its importance to heterosis.
This paper aims to explain how overdominance, and some features of epistasis,
can similarly emerge from the molecular dynamics of proteins. Possessing
multiple alleles at a gene locus results in the synthesis of different
allozymes at reduced concentrations. This in turn reduces the rate at which
each allozyme forms soluble oligomers, which are toxic and must be degraded,
because allozymes coaggregate at low efficiencies. The model developed in this
paper will be used to explain how heterozygosity can impact the metabolic
efficiency of an organism. It can also explain why the viabilities of some
inbred lines seem to decline rapidly at high inbreeding coefficients (F > 0.5),
which may provide a physical basis for truncation selection for heterozygosity.
Finally, the model has implications for the ploidy level of organisms. It can
explain why polyploids are frequently found in environments where severe
physical stresses promote the formation of soluble oligomers. The model can
also explain why complex organisms, which need to synthesize aggregationprone
proteins that contain intrinsically unstructured regions (IURs) and multiple
domains because they facilitate complex protein interaction networks (PINs),
tend to be diploid while haploidy tends to be restricted to relatively simple
organisms.

We investigate the asymptotic state of a periodically driven manybody
quantum system which is weakly coupled to an environment. The combined action
of the modulations and the environment steers the system towards a state being
characterized by a timeperiodic density operator. To resolve this asymptotic
nonequilibrium state at stroboscopic instants of time, we introduce the
dissipative Floquet map, evaluate the stroboscopic density operator as its
eigenelement and elucidate how particle interactions affect properties of the
density operator. We illustrate the idea with a periodically modulated
BoseHubbard dimer and discuss the relations between the interactioninduced
bifurcations in a meanfield dynamics and changes in the characteristics of the
genuine quantum manybody state. We argue that Floquet maps provide insight
into the system relaxation towards its asymptotic state and may help to
understand whether it is possible (or not) to construct a stroboscopic
timeindependent generator mimicking the action of the original timedependent
one.

We consider a system composed of a fixed number of particles with total
energy smaller or equal to some prescribed value. The particles are
noninteracting, indistinguishable and distributed over fixed number of energy
levels. The energy levels are degenerate and degeneracy is a function of the
number of particles. Three cases of the degeneracy function is considered. It
can increase with either the same rate as the number of particles or slower, or
faster. We find useful properties of the entropy of the system and solve the
entropy optimization problem. It turned out, there are several solutions.
Depending on the magnitude of the total energy, the maximum of the entropy can
be in the interior of system's state space or on the boundary. On the boundary
it can have further three cases depending on the degeneracy function. The main
result, Law of Large Numbers yields the most probable state of the system,
which equals to the point of maximum of the entropy. This point can be either
MaxwellBoltzmann statistics or BoseEinstein statistics, or ZipfMandelbort
law. We also find the limiting laws for the fluctuations. These laws are
different for different cases of the entropy's maximum. They can be mixture of
Normal, Exponential and Discrete distributions. Explicit rate of convergence is
provided for all the theorems.

We present a new state function, denotes the energy of thermal motion, it
shows that the energy classification for the internal energy has been
completed. In a new theoretical structure, we present the concept of the heat
partial pressure, rewrite the equations of the first law, by which, the energy
conversion and energy transport can be distinguished by an explicit way, we
redefine the concept of the entropy, and confirm a conversion potential related
to the dissipative compensation, we discuss in detail on the three sources of
irreversibility, and the conversions between the calorimetric entropy and the
configurational entropy, an interesting conclusion shows that the second law
itself has contained the mechanism of evolution.

Internal diffusionlimited aggregation (IDLA) is a stochastic growth model on
a graph $G$ which describes the formation of a random set of vertices growing
from the origin (some fixed vertex) of $G$. Particles start at the origin and
perform simple random walks; each particle moves until it lands on a site which
was not previously visited by other particles. This random set of occupied
sites in $G$ is called the IDLA cluster.
In this paper we consider IDLA on Sierpinski gasket graphs, and show that the
IDLA cluster fills balls (in the graph metric) with probability 1.

For conformal field theories, it is shown how the Ward identity corresponding
to dilatation invariance arises in a Wilsonian setting. In so doing, several
points which are opaque in textbook treatments are clarified. Exploiting the
fact that the Exact Renormalization Group furnishes a representation of the
conformal algebra allows dilatation invariance to be stated directly as a
property of the action, despite the presence of a regulator. This obviates the
need for formal statements that conformal invariance is recovered once the
regulator is removed. Furthermore, the proper subset of conformal primary
fields for which the Ward identity holds is identified for all
dimensionalities.

For nonequilibrium systems in a steady state we present two necessary and
sufficient conditions for the emergence of $q$canonical ensembles, also known
as Tsallis statistics. These conditions are invariance requirements over the
definition of subsystem and environment, and over the joint rescaling of
temperature and energy. Our approach is independent of, but complementary to,
the notions of Tsallis nonextensive statistics and Superstatistics.

A Legendre transform of the recently discovered conformal fixedpoint
equation is constructed, providing an unintegrated equation encoding full
conformal invariance within the framework of the effective average action.

We study the slowing down of particle beams passing through the dusty plasma
with powerlaw kappadistributions. Three plasma components, electrons, ions
and dust particles, can have different kappaparameter. We derive the
deceleration factor (the velocity moment equation) and the slowing down time of
a test particle, and numerically study the slowing down properties of an
electron beam, a proton beam and a dust particle beam, respectively, in the
kappadistributed dusty plasma. We show that the slowing down properties of
particle beams depend strongly on the kappaparameters of the plasma
components, and the dust component plays a dominant role in the slowing down.
And the slowing down also depends on mass and charge of the dust particles in
the dusty plasma. More detailed results are shown in 17 numerical graphs.

We investigate the size scaling of the entanglement entropy (EE) in
nonequilibrium steady states (NESSs) of a onedimensional open quantum system
with a random potential. It models a mesoscopic conductor, composed of a long
quantum wire (QWR) with impurities and two electron reservoirs at zero
temperature. The EE at equilibrium obeys the logarithmic law. However, in NESSs
far from equilibrium the EE grows anomalously fast, obeying the `quasi volume
law,' although the conductor is driven by the zerotemperature reservoirs. This
anomalous behavior arises from both the far from equilibrium condition and
multiple scatterings due to impurities.

The results of the renormalization group are commonly advertised as the
existence of power law singularities near critical points. The classic
predictions are often violated and logarithmic and exponential corrections are
treated on a casebycase basis. We use the mathematics of normal form theory
to systematically group these into universality families of seemingly unrelated
systems united by common scaling variables. We recover and explain the existing
literature and predict the nonlinear generalization for the universal
homogeneous scaling functions. We show that this procedure leads to a better
handling of the singularity even in classic cases and elaborate our framework
using several examples.

The statistics of the smallest eigenvalue of WishartLaguerre ensemble is
important from several perspectives. The smallest eigenvalue density is
typically expressible in terms of determinants or Pfaffians. These results are
of utmost significance in understanding the spectral behavior of
WishartLaguerre ensembles and, among other things, unveil the underlying
universality aspects in the asymptotic limits. However, obtaining exact and
explicit expressions by expanding determinants or Pfaffians becomes impractical
if large dimension matrices are involved. For the real matrices ($\beta=1$)
Edelman has provided an efficient recurrence scheme to work out exact and
explicit results for the smallest eigenvalue density which does not involve
determinants or matrices. Very recently, an analogous recurrence scheme has
been obtained for the complex matrices ($\beta=2$). In the present work we
extend this to $\beta$WishartLaguerre ensembles for the case when exponent
$\alpha$ in the associated Laguerre weight function, $\lambda^\alpha
e^{\beta\lambda/2}$, is a nonnegative integer, while $\beta$ is positive
real. This also gives access to the smallest eigenvalue density of fixed trace
$\beta$WishartLaguerre ensemble, as well as moments for both cases. Moreover,
comparison with earlier results for the smallest eigenvalue density in terms of
certain hypergeometric function of matrix argument results in an effective way
of evaluating these explicitly. Exact evaluations for large values of $n$ (the
matrix dimension) and $\alpha$ also enable us to compare with TracyWidom
density and large deviation results of Katzav and Castillo. We also use our
result to obtain the density of the largest of the proper delay times which are
eigenvalues of the WignerSmith matrix and are relevant to the problem of
quantum chaotic scattering.

Introduced in the late 1960's, the asymmetric exclusion process (ASEP) is an
important model from statistical mechanics which describes a system of
interacting particles hopping left and right on a onedimensional lattice with
open boundaries. It has been known for awhile that there is a tight connection
between the partition function of the ASEP and moments of AskeyWilson
polynomials, a family of orthogonal polynomials which are at the top of the
hierarchy of classical orthogonal polynomials in one variable. On the other
hand, AskeyWilson polynomials can be viewed as a specialization of the
multivariate MacdonaldKoornwinder polynomials (also known as Koornwinder
polynomials), which in turn give rise to the Macdonald polynomials associated
to any classical root system via a limit or specialization. In light of the
fact that Koornwinder polynomials generalize the AskeyWilson polynomials, it
is natural to ask whether one can find a particle model whose partition
function is related to Koornwinder polynomials. In this article we answer this
question affirmatively, by showing that the "homogeneous" Koornwinder moments
at q=t recover the partition function for the twospecies exclusion process. We
also provide a "hook length" formula for Koornwinder moments when q=t=1.

The relationship between anomalous superdiffusive behavior and particle
trapping probability is analyzed on a rocking ratchet potential with spatially
correlated weak disorder. The trapping probability density is shown,
analytically and numerically, to have an exponential form as a function of
space. The trapping processes with a low or no thermal noise are only
transient, but they can last much longer than the characteristic time scale of
the system and therefore might be detected experimentally. Using the result for
the trapping probability we obtain an analytical expression for the number of
wells where a given number of particles are trapped. We have also obtained an
analytical approximation for the secondmoment of the particle distribution
function C2 as a function of time, when trapped particles coexist with constant
velocity untrapped particles. We also use the expression for C2 to characterize
the anomalous superdiffusive motion in the absence of thermal noise for the
transient time.

This work constructs a welldefined and operational form factor expansion in
a model having a massless spectrum of excitations. More precisely, the dynamic
twopoint functions in the massless regime of the XXZ spin1/2 chain are
expressed in terms of properly regularised series of multiple integrals. These
series are obtained by taking, in an appropriate way, the thermodynamic limit
of the finite volume form factor expansions. The series are structured in way
allowing one to identify directly the contributions to the correlator stemming
from the conformaltype excitations on the Fermi surface and those issuing from
the massive excitations (deep holes, particles and bound states). The obtained
form factor series opens up the possibility of a systematic and exact study of
asymptotic regimes of dynamical correlation functions in the massless regime of
the XXZ spin $1/2$ chain. Furthermore, the assumptions on the microscopic
structure of the model's Hilbert space that are necessary so as to write down
the series appear to be compatible with any model  not necessarily integrable
 belonging to the Luttinger liquid universality class. Thus, the present
analysis provides also the phenomenological structure of form factor expansions
in massless models belonging to this universality class.