
Dense kernel matrices $\Theta \in \mathbb{R}^{N \times N}$ obtained from
point evaluations of a covariance function $G$ at locations $\{ x_{i} \}_{1
\leq i \leq N}$ arise in statistics, machine learning, and numerical analysis.
For covariance functions that are Green's functions of elliptic boundary value
problems and homogeneouslydistributed sampling points, we show how to identify
a subset $S \subset \{ 1 , \dots , N \}^2$, with $\# S = O ( N \log (N)
\log^{d} ( N /\epsilon ) )$, such that the zero fillin incomplete Cholesky
factorisation of the sparse matrix $\Theta_{ij} 1_{( i, j ) \in S}$ is an
$\epsilon$approximation of $\Theta$. This factorisation can provably be
obtained in complexity $O ( N \log( N ) \log^{d}( N /\epsilon) )$ in space and
$O ( N \log^{2}( N ) \log^{2d}( N /\epsilon) )$ in time; we further present
numerical evidence that $d$ can be taken to be the intrinsic dimension of the
data set rather than that of the ambient space. The algorithm only needs to
know the spatial configuration of the $x_{i}$ and does not require an analytic
representation of $G$. Furthermore, this factorization straightforwardly
provides an approximate sparse PCA with optimal rate of convergence in the
operator norm. Hence, by using only subsampling and the incomplete Cholesky
factorization, we obtain, at nearly linear complexity, the compression,
inversion and approximate PCA of a large class of covariance matrices. By
inverting the order of the Cholesky factorization we also obtain a solver for
elliptic PDE with complexity $O ( N \log^{d}( N /\epsilon) )$ in space and $O (
N \log^{2d}( N /\epsilon) )$ in time.

Power flow calculations for systems with a large number of buses, e.g. grids
with multiple voltage levels, or time series based calculations result in a
high computational effort. A common power flow solver for the efficient
analysis of power systems is the NewtonRaphson algorithm. The main
computational effort of this method results from the linearization of the
nonlinear power flow problem and solving the resulting linear equation. This
paper presents an algorithm for the fast linearization of the power flow
problem by creating the Jacobian matrix directly in CRS format. The increase in
speed is achieved by reducing the number of iterations over the nonzero
elements of the sparse Jacobian matrix. This allows to efficiently create the
Jacobian matrix without having to approximate the problem. A comparison of the
calculation time of three power grids shows that comparable opensource
implementations need 314x the time to create the Jacobian matrix.

pandapower is a Python based, BSDlicensed power system analysis tool aimed
at automation of static and quasistatic analysis and optimization of balanced
power systems. It provides power flow, optimal power flow, state estimation,
topological graph searches and short circuit calculations according to IEC
60909. pandapower includes a NewtonRaphson power flow solver formerly based on
PYPOWER, which has been accelerated with justintime compilation. Additional
enhancements to the solver include the capability to model constant current
loads, grids with multiple reference nodes and a connectivity check. The
pandapower network model is based on electric elements, such as lines, two and
threewinding transformers or ideal switches. All elements can be defined with
nameplate parameters and are internally processed with equivalent circuit
models, which have been validated against industry standard software tools. The
tabular data structure used to define networks is based on the Python library
pandas, which allows comfortable handling of input and output parameters. The
implementation in Python makes pandapower easy to use and allows comfortable
extension with thirdparty libraries. pandapower has been successfully applied
in several grid studies as well as for educational purposes. A comprehensive,
publicly available casestudy demonstrates a possible application of pandapower
in an automated time series calculation.

We present experimental results on the inelastic and elastic interspecies
interactions between ytterbium (Yb) in the metastable ${}^3\mathrm{P}_2$ state
loaded into a deep optical lattice and spin polarized lithium (Li) in its
ground state. Focusing on the $m_J = 0$ magnetic sublevel of
Yb(${}^3\mathrm{P}_2$), bias magnetic fields between 20 G and 800 G are
investigated and significantly enhanced inelastic collision rates with high
magnetic fields are found. In addition, by direct spectroscopy of the Yb
Mottinsulator immersed in the Li Fermi gas an upper boundary of the background
scattering length of the Yb(${}^3\mathrm{P}_2,
m_J=0$)Li(${}^2\mathrm{S}_{1/2}, F=1/2, m_F=+1/2$) system is estimated,
revealing the absence of useful Feshbach resonances. These observations are
qualitatively consistent with the theoretical calculations.

Experimentally the spin dependence of inelastic collisions between ytterbium
(Yb) in the metastable 3P0 state and lithium (Li) in the Li ground state
manifold is investigated at low magnetic fields. Using selective excitation all
magnetic sublevels mJ of 174Yb(3P0) are accessed and four of the six lowest
lying magnetic sublevels of 6Li are prepared by optical pumping. On the one
hand, mJindependence of collisions involving Li(F=1/2) atoms is found. A
systematic mJdependence in collisions with Li(F=3/2) atoms, in particular
suppressed losses for stretched collisional states, is observed on the other
hand. Further, mJchanging processes are found to be of minor relevance. The
span of observed inelastic collision rates is between 1*10^{11} cm^3/s and
40*10^{11} cm^3/s, and a possible origin of the observed behavior is
discussed.

The space of images can be equipped with a Riemannian metric measuring both
the cost of transport of image intensities and the variation of image
intensities along motion lines. The resulting metamorphosis model was
introduced and analyzed by Trouv\'e and Younes, and a variational time
discretization for the geodesic interpolation was proposed by Berkels et al. In
this paper, this time discrete model is expanded and an image extrapolation via
a discrete exponential map is consistently derived for the variational time
discretization. For a given weakly differentiable initial image and an initial
image variation, the exponential map allows to compute a discrete geodesic
extrapolation path in the space of images. It is shown that a time step of this
shooting method can be formulated in the associated deformations only. For
sufficiently small time steps local existence and uniqueness are proved using a
suitable fixed point formulation and the implicit function theorem. A spatial
Galerkin discretization with cubic splines on coarse meshes for the deformation
and piecewise bilinear finite elements on fine meshes for the image intensities
are used to derive a fully practical algorithm. Different applications
underline the efficiency and stability of the proposed approach.

We establish an experimental method for a detailed investigation of inelastic
collisional properties between ytterbium (Yb) in the metastable
${}^3\mathrm{P}_2$ state and ground state lithium (Li). By combining an optical
lattice and a direct excitation to the ${}^3\mathrm{P}_2$ state we achieve high
selectivity on the collisional partners. Using this method we determine
inelastic loss coefficients in collisions between
$^{174}$Yb(${}^3\mathrm{P}_2$) with magnetic sublevels of $m_J=0$ and $2$ and
ground state $^6$Li to be $(4.4\pm0.3)\times10^{11}~\mathrm{cm}^3/\mathrm{s}$
and $(4.7\pm0.8)\times10^{11}~\mathrm{cm}^3/\mathrm{s}$, respectively. Absence
of spin changing processes in Yb(${}^3\mathrm{P}_2$)Li inelastic collisions at
low magnetic fields is confirmed by inelastic loss measurements on the $m_J=0$
state. We also demonstrate that our method allows us to look into loss
processes in fewbody systems separately.

Correlations in systems with spin degree of freedom are at the heart of
fundamental phenomena, ranging from magnetism to superconductivity. The effects
of correlations depend strongly on dimensionality, a striking example being
onedimensional (1D) electronic systems, extensively studied theoretically over
the past fifty years. However, the experimental investigation of the role of
spin multiplicity in 1D fermions  and especially for more than two spin
components  is still lacking. Here we report on the realization of 1D,
stronglycorrelated liquids of ultracold fermions interacting repulsively
within SU(N) symmetry, with a tunable number N of spin components. We observe
that static and dynamic properties of the system deviate from those of ideal
fermions and, for N>2, from those of a spin1/2 Luttinger liquid. In the
largeN limit, the system exhibits properties of a bosonic spinless liquid. Our
results provide a testing ground for manybody theories and may lead to the
observation of fundamental 1D effects.

Atom chips provide compact and robust platforms towards practical quantum
technologies. A quick and faithful preparation of arbitrary input states for
these systems is crucial but represents a very challenging experimental task.
This is especially difficult when the dynamical evolution is noisy and
unavoidable setup imperfections have to be considered. Here, we experimentally
prepare with very high small errors different internal states of a Rubidium
BoseEinstein condensate realized on an atom chip. As a possible application of
our scheme, we apply it to improve the sensitivity of an atomic interferometer.

It is generally impossible to probe a quantum system without disturbing it.
However, it is possible to exploit the backaction of quantum measurements and
strong couplings to tailor and protect the coherent evolution of a quantum
system. This is a profound and counterintuitive phenomenon known as quantum
Zeno dynamics (QZD). Here we demonstrate QZD with a rubidium BoseEinstein
condensate in a fivelevel Hilbert space. We harness measurements and strong
couplings to dynamically disconnect different groups of quantum states and
constrain the atoms to coherently evolve inside a twolevel subregion. In
parallel to the foundational importance due to the realization of a dynamical
superselection rule and the theory of quantum measurements, this is an
important step forward in protecting and controlling quantum dynamics and,
broadly speaking, quantum information processing.

Scattering is an important phenomenon which is observed in systems ranging
from the micro to macroscale. In the context of nuclear reaction theory the
Heidelberg approach was proposed and later demonstrated to be applicable to
many chaotic scattering systems. To model the universal properties,
stochasticity is introduced to the scattering matrix on the level of the
Hamiltonian by using random matrices. A longstanding problem was the
computation of the distribution of the offdiagonal scatteringmatrix elements.
We report here an exact solution to this problem and present analytical results
for systems with preserved and with violated timereversal invariance. Our
derivation is based on a new variant of the supersymmetry method. We also
validate our results with scattering data obtained from experiments with
microwave billiards.