
Certain lattice wave systems in translationally invariant settings have one
or more spectral bands that are strictly flat or independent of momentum in the
tight binding approximation, arising from either internal symmetries or
finetuned coupling. These flat bands display remarkable stronglyinteracting
phases of matter. Originally considered as a theoretical convenience useful for
obtaining exact analytical solutions of ferromagnetism, flat bands have now
been observed in a variety of settings, ranging from electronic systems to
ultracold atomic gases and photonic devices. Here we review the design and
implementation of flat bands and chart future directions of this exciting
field.

The microcanonical GrossPitaevskii (aka semiclassical BoseHubbard) lattice
model dynamics is characterized by a pair of energy and norm densities. The
grand canonical Gibbs distribution fails to describe a part of the density
space, due to the boundedness of its kinetic energy spectrum. We define
Poincare equilibrium manifolds and compute the statistics of microcanonical
excursion times off them. The tails of the distribution functions quantify the
proximity of the manybody dynamics to a weaklynonergodic phase, which occurs
when the average excursion time is infinite. We find that a crossover to
weaklynonergodic dynamics takes place inside the nonGibbs phase, being
unnoticed by the largest Lyapunov exponent. In the ergodic part of the
nonGibbs phase, the Gibbs distribution should be replaced by an unknown
modified one. We relate our findings to the corresponding integrable limit,
close to which the actions are interacting through a short range coupling
network.

Linear wave equations on Hamiltonian lattices with translational invariance
are characterized by an eigenvalue band structure in reciprocal space. Flat
band lattices have at least one of the bands completely dispersionless. Such
bands are coined flat bands. Flat bands occur in finetuned networks, and can
be protected by (e.g. chiral) symmetries. Recently a number of such systems
were realized in structured optical systems, excitonpolariton condensates, and
ultracold atomic gases. Flat band networks support compact localized modes.
Local defects couple these compact modes to dispersive states and generate Fano
resonances in the wave propagation. Disorder (i.e. a finite density of defects)
leads to a dense set of Fano defects, and to novel scaling laws in the
localization length of disordered dispersive states. Nonlinearities can
preserve the compactness of flat band modes, along with renormalizing (tuning)
their frequencies. These strictly compact nonlinear excitations induce tunable
Fano resonances in the wave propagation of a nonlinear flat band lattice.

We study nonHermitian photonic lattices that exhibit competition between
conservative and nonHermitian (gain/loss) couplings. A bipartite sublattice
symmetry enforces the existence of nonHermitian flat bands, which are
typically embedded in an auxiliary dispersive band and give rise to
nondiffracting "compact localized states". Band crossings take the form of
nonHermitian degeneracies known as exceptional points. Excitations of the
lattice can produce either diffracting or amplifying behaviors. If the
nonHermitian coupling is finetuned to generate an effective $\pi$ flux, the
lattice spectrum becomes completely flat, a nonHermitian analogue of
AharonovBohm caging in which the magnetic field is replaced by balanced gain
and loss. When the effective flux is zero, the nonHermitian band crossing
points give rise to asymmetric diffraction and anomalous linear amplification.

In the study of subdiffusive wavepacket spreading in disordered KleinGordon
(KG) nonlinear lattices, a central open question is whether the motion
continues to be chaotic despite decreasing densities, or tends to become
quasiperiodic as nonlinear terms become negligible. In a recent study of such
KG particle chains with quartic (4th order) anharmonicity in the onsite
potential, it was shown that $q$Gaussian probability distribution functions of
sums of position observables with $q > 1$ always approach pure Gaussians
($q=1$) in the long time limit and hence the motion of the full system is
ultimately "strongly chaotic". In the present paper, we show that these results
continue to hold even when a sextic (6th order) term is gradually added to the
potential and ultimately prevails over the 4th order anharmonicity, despite
expectations that the dynamics is more "regular", at least in the regime of
small oscillations. Analyzing this system in the subdiffusive energy domain
using $q$statistics, we demonstrate that groups of oscillators centered around
the initially excited one (as well as the full chain) possess strongly chaotic
dynamics and are thus far from any quasiperiodic torus, for times as long as
$t=10^9$.

We analyze the spectrum and eigenstates of a quantum particle in a bipartite
twodimensional tightbinding dice network with short range hopping under the
action of a dc bias. We find that the energy spectrum consists of a periodic
repetition of onedimensional energy band multiplets, with one member in the
multiplet being strictly flat. The corresponding macroscopic degeneracy invokes
eigenstates localized exponentially perpendicular to the dc field direction,
and superexponentially along the dc field direction. We also show that the
band multiplet is characterized by a topological winding number (Zak phase),
which changes abruptly if we vary the dc field strength. These changes are
induced by gap closings between the flat and dispersive bands, and reflect the
number of these closings.

We analyze the stability of a nonHermitian coupler with respect to
modulational inhomogeneous perturbations in the presence of unbalanced gain and
loss. At the paritytime (PT) symmetry point the coupler is unstable. Suitable
symmetry breakings lead to an asymmetric coupler, which hosts nonlinear
supermodes. A subset of those broken symmetry cases finally yields nonlinear
supermodes which are stable against modulational perturbations. The lack of
symmetry requirements will facilitate experimental implementations and relevant
photonics applications.

We study flat bands in bipartite tightbinding networks with discrete
translational invariance. Chiral flat bands with chiral symmetry eigenenergy E
= 0 and host compact localized eigenstates for finite range hopping. For a
bipartite network with a majority sublattice chiral flat bands emerge. We
present a simple generating principle of chiral flat band networks and as a
showcase add to the previously observed cases a number of new potentially
realizable chiral flat bands in various lattice dimensions. Chiral symmetry
respecting network perturbations  including disorder and synthetic magnetic
fields  preserve both the flatband and the modified compact localized states.
Chiral flatbands are spectrally protected by gaps, and pseudogaps in the
presence of disorder due to Griffiths effects.

Certain tight binding lattices host macroscopically degenerate flat spectral
bands. Their origin is rooted in local symmetries of the lattice, with
destructive interference leading to the existence of compact localized
eigenstates. We study the robustness of this localization to disorder in
different classes of flat band lattices in one and two dimensions. Depending on
the flat band class, the flat band states can either be robust, preserving
their strong localization for weak disorder W, or they are destroyed and
acquire large localization lengths $\xi$ that diverge with a variety of
unconventional exponents $\nu$, $\xi \sim 1/W^{\nu}$.

Flat bands (FB) are strictly dispersionless bands in the Bloch spectrum of a
periodic lattice Hamiltonian, recently observed in a variety of photonic and
dissipative condensate networks. FB Hamiltonians are finetuned networks, still
lacking a comprehensive generating principle. We introduce a FB generator based
on local network properties. We classify FB networks through the properties of
compact localized states (CLS) which are exact FB eigenstates and occupy $U$
unit cells. We obtain the complete twoparameter FB family of twoband $d=1$
networks with nearest unit cell interaction and $U=2$. We discover a novel high
symmetry sawtooth chain with identical hoppings in a transverse dc field,
easily accessible in experiments. Our results pave the way towards a complete
description of FBs in networks with more bands and in higher dimensions.

We study the fate of dynamical localization of two quantum kicked rotors with
contact interaction. This interaction mimics experimental realizations with
ultracold atomic gases. Dynamical localization for a single rotor takes place
in momentum space. The contact interaction affects the evolution of the
relative momentum $k$ of a pair of interacting rotors in a nonanalytic way.
Consequently the evolution operator $U$ is exciting large relative momenta with
amplitudes which decay only as a power law $1/k^4$. This is in contrast to the
centerofmass momentum $K$ for which the amplitudes excited by $U$ decay
superexponentially fast. Therefore dynamical localization is preserved for the
centerofmass momentum, but destroyed for the relative momentum for any
nonzero strength of interaction.

We consider the dynamics of noninteracting electrons on a square lattice in
the presence of a magnetic flux {\alpha} and a dc electric field E oriented
along the lattice diagonal. In general, the adiabatic dynamics of an electron
will be characterized by Bloch oscillations in the electrical field direction
and dispersive ballistic transport in the perpendicular direction. For rational
values of {\alpha} and a corresponding discrete set of values of E({\alpha})
vanishing gaps in the spectrum induce a fractionalization of the charge in the
perpendicular direction  while left movers are still performing dispersive
ballistic transport, the complementary fraction of right movers is propagating
in a dispersionless relativistic manner in the opposite direction.
Generalizations to other field directions, lattice symmetries and the possible
probing of the effect with atomic BoseEinstein condensates and photonic
networks are discussed.

We consider the asymmetric active coupler (AAC) consisting of two coupled
dissimilar waveguides with gain and loss. We show that under generic
conditions, not restricted by paritytime symmetry, there exist finitepower,
constantintensity nonlinear supermodes (NS), resulting from the balance
between gain, loss, nonlinearity, coupling and dissimilarity. The system is
shown to possess nonreciprocal dynamics enabling directed power transport and
optical isolation functionality.

Sinusoidal Bloch oscillations appear in band structures exposed to external
fields. LandauZener (LZ) tunneling between different bands is usually a
counteracting effect limiting Bloch oscillations. Here we consider a flat band
network with two dispersive and one flat bands, e.g. for ultracold atoms and
optical waveguide networks. Using external synthetic gauge and gravitational
fields we obtain a perturbed yet gapless band structure with almost flat parts.
The resulting Bloch oscillations consist of two parts  a fast scan through the
nonflat part of the dispersion structure, and an almost complete halt for
substantial time when the atomic wave packet is trapped in the original flat
band part of the unperturbed spectrum, made possible due to LZ tunneling.

Localization of waves by disorder is a fundamental physical problem
encompassing a diverse spectrum of theoretical, experimental and numerical
studies in the context of metalinsulator transitions, the quantum Hall effect,
light propagation in photonic crystals, and dynamics of ultracold atoms in
optical arrays, to name just a few examples. Large intensity light can induce
nonlinear response, ultracold atomic gases can be tuned into an interacting
regime, which leads again to nonlinear wave equations on a mean field level.
The interplay between disorder and nonlinearity, their localizing and
delocalizing effects is currently an intriguing and challenging issue in the
field of lattice waves. In particular it leads to the prediction and
observation of two different regimes of destruction of Anderson localization 
asymptotic weak chaos, and intermediate strong chaos, separated by a crossover
condition on densities. On the other side approximate full quantum interacting
many body treatments were recently used to predict and obtain a novel many body
localization transition, and two distinct phases  a localization phase, and a
delocalization phase, both again separated by some typical density scale. We
will discuss selftrapping, nonergodicity and nonGibbsean phases which are
typical for such discrete models with particle number conservation and their
relation to the above crossover and transition physics. We will also discuss
potential connections to quantum many body theories.

Properly modulated flatband lattices have a divergent density of states at
the flatband energy. Quasiperiodic modulations are known to host a metal
insulator transition already in one space dimension. Their embedding into
flatband geometries consequently allows for a precise engineering and fine
tuning of mobility edges. We obtain analytic expressions for singular mobility
edges for two flatband lattice examples. In particular, we engineer cases with
arbitrarily small energy separations of mobility edge, zeroes, and
divergencies.

Flat band networks are characterized by coexistence of dispersive and flat
bands. Flat bands (FB) are generated by compact localized eigenstates (CLS)
with local network symmetries, based on destructive interference. Correlated
disorder and quasiperiodic potentials hybridize CLS without additional
renormalization, yet with surprising consequencies: (i) states are expelled
from the FB energy $E_{FB}$, (ii) the localization length of eigenstates
vanishes as $\xi \sim 1 / \ln (E E_{FB})$, (iii) the density of states
diverges logarithmically (particlehole symmetry) and algebraically (no
particlehole symmetry), (iv) mobility edge curves show algebraic singularities
at $E_{FB}$. Our analytical results are based on perturbative expansions of the
CLS, and supported by numerical data in one and two lattice dimensions.

We study spreading wave packets in a disordered nonlinear ladder with broken
timereversal symmetry induced by synthetic gauge fields. The model describes
the dynamics of interacting bosons in a disordered and driven optical ladder
within a meanfield approximation. The second moment of the wave packet $m_{2}
= g t^{\alpha}$ grows subdiffusively with the universal exponent $\alpha \simeq
1/3$ similar to the timereversal case. However the prefactor $g$ is strongly
modified by the field strength and shows a nonmonotonic dependence. For a weak
field, the prefactor increases since timereversal enhanced backscattering is
suppressed. For strong fields the spectrum of the linear wave equation reduces
the localization length through the formation of gaps and narrow bands.
Consequently the prefactor for the subdiffusive spreading law is suppressed.

Localization of waves by disorder is a fundamental physical problem
encompassing a diverse spectrum of theoretical, experimental and numerical
studies in the context of metalinsulator transition, quantum Hall effect,
light propagation in photonic crystals, and dynamics of ultracold atoms in
optical arrays. Large intensity light can induce nonlinear response, ultracold
atomic gases can be tuned into an interacting regime, which leads again to
nonlinear wave equations on a mean field level. The interplay between disorder
and nonlinearity, their localizing and delocalizing effects is currently an
intriguing and challenging issue in the field. We will discuss recent advances
in the dynamics of nonlinear lattice waves in random potentials. In the absence
of nonlinear terms in the wave equations, Anderson localization is leading to a
halt of wave packet spreading.
Nonlinearity couples localized eigenstates and, potentially, enables
spreading and destruction of Anderson localization due to nonintegrability,
chaos and decoherence. The spreading process is characterized by universal
subdiffusive laws due to nonlinear diffusion. We review extensive computational
studies for one and twodimensional systems with tunable nonlinearity power.
We also briefly discuss extensions to other cases where the linear wave
equation features localization: AubryAndre localization with quasiperiodic
potentials, WannierStark localization with dc fields, and dynamical
localization in momentum space with kicked rotors.

Macroscopically degenerate flat bands (FB) in periodic lattices host compact
localized states which appear due to destructive interference and local
symmetry. Interference provides a deep connection between the existence of flat
band states (FBS) and the appearance of Fano resonances for wave propagation.
We introduce generic transformations detangling FBS and dispersive states into
lattices of Fano defects. Inverting the transformation, we generate a continuum
of FB models. Our procedure allows us to systematically treat perturbations
such as disorder and explain the emergence of energydependent localization
length scaling in terms of Fano resonances.

Sil, Maiti, and Chakrabarti (SMC) (Phys. Rev. Lett. 101, 076803 (2008),
arXiv:0801.2670) introduce an aperiodic twoleg ladder network composed of
atomic sites with onsite potentials distributed according to a quasiperiodic
AubryAndre potential. SMC claim the existence of multiple mobility edges, i.e.
metalinsulator transitions at multiple values of the Fermi energy. SMC use
numerical calculations of the conductance and den sity of states, and an
analytical result in a limiting case. In the following, we restate in the real
space the change of basis done by SMC. This change of basis reduces their model
to two decoupled chains, each with its own Aubry Andre potential. Each of the
chains has its own critical point, where all states change from metallic to
insulat ing, without any notion of a mobility edge. Since the two chains are
not interacting, their spectra may overlap such that localized states (from one
chain) and delocal ized states (from the other chain) coexist at the same en
ergy. Consequently, the analysis of the density of states and the conductance
alone is misleading, and that is due to the choice of the symmetric
quasiperiodic potential. We briefly outline the possibility to obtain true
mobility edges.

Transport properties of particles and waves in spatially periodic structures
that are driven by external timedependent forces manifestly depend on the
spacetime symmetries of the corresponding equations of motion. A systematic
analysis of these symmetries uncovers the conditions necessary for obtaining
directed transport. In this work we give a unified introduction into the
symmetry analysis and demonstrate its action on the motion in onedimensional
periodic, both in time and space, potentials. We further generalize the
analysis to quasiperiodic drivings, higher space dimensions, and quantum
dynamics. Recent experimental results on the transport of cold and ultracold
atomic ensembles in acdriven optical potentials are reviewed as illustrations
of theoretical considerations.

We study the critical behaviour of Anderson localized modes near intersecting
flat and dispersive bands in the quasionedimensional diamond ladder with weak
diagonal disorder $W$. The localization length $\xi$ of the flat band states
scales with disorder as $\xi \sim W^{\gamma}$, with $\gamma \approx 1.3$, in
contrast to the dispersive bands with $\gamma =2$. A small fraction of
dispersive modes mixed with the flat band states is responsible for the unusual
scaling. Anderson localization is therefore controlled by two different length
scales. Nonlinearity can produce qualitatively different wave spreading
regimes, from enhanced expansion to resonant tunneling and selftrapping.

Do nonlinear waves destroy Anderson localization? Computational and
experimental studies yield subdiffusive nonequilibrium wave packet spreading.
Chaotic dynamics and phase decoherence assumptions are used for explaining the
data. We perform a quantitative analysis of the nonequilibrium chaos
assumption, and compute the time dependence of main chaos indicators  Lyapunov
exponents and deviation vector distributions. We find a slowing down of chaotic
dynamics, which does not cross over into regular dynamics up to the largest
observed time scales, still being fast enough to allow for a thermalization of
the spreading wave packet. Strongly localized chaotic spots meander through the
system as time evolves. Our findings confirm for the first time that
nonequilibrium chaos and phase decoherence persist, fueling the prediction of a
complete delocalization.

Anderson and dynamical localization have been experimentally observed with
ultracold atomic matter. Feshbach resonances are used to efficiently control
the strength of interactions between atoms. This allows to study the
delocalization effect of interactions for localized wave packets. The
delocalization processes are subdiffusive and slow, thereby limiting the
quantitative experimental and numerical analysis. We propose an elegant
solution of the problem by proper ramping the interaction strength in time. We
demonstrate that subdiffusion is speeded up to normal diffusion for interacting
disordered and kicked atomic systems. The door is open to test these
theoretical results experimentally, and to attack similar computational quests
in higher space dimensions