• ### Compression, inversion, and approximate PCA of dense kernel matrices at near-linear computational complexity(1706.02205)

March 22, 2019 math.NA, math.PR, cs.DS, cs.CC
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 homogeneously-distributed 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 fill-in 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.
• ### An efficient open-source implementation to compute the Jacobian matrix for the Newton-Raphson power flow algorithm(1804.06742)

April 18, 2018 stat.CO, cs.CE
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 Newton-Raphson 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 open-source implementations need 3-14x the time to create the Jacobian matrix.
• ### pandapower - an Open Source Python Tool for Convenient Modeling, Analysis and Optimization of Electric Power Systems(1709.06743)

April 18, 2018 cs.CE
pandapower is a Python based, BSD-licensed power system analysis tool aimed at automation of static and quasi-static 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 Newton-Raphson power flow solver formerly based on PYPOWER, which has been accelerated with just-in-time 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 three-winding 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 third-party libraries. pandapower has been successfully applied in several grid studies as well as for educational purposes. A comprehensive, publicly available case-study demonstrates a possible application of pandapower in an automated time series calculation.
• ### Spectroscopic determination of magnetic-field-dependent interactions in an ultracold Yb(3P2)-Li mixture(1708.00215)

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 Mott-insulator 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.
• ### Spin dependent inelastic collisions between metastable state two-electron atoms and ground state alkali-atoms(1705.06841)

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, mJ-independence of collisions involving Li(F=1/2) atoms is found. A systematic mJ-dependence in collisions with Li(F=3/2) atoms, in particular suppressed losses for stretched collisional states, is observed on the other hand. Further, mJ-changing 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.
• ### Image Extrapolation for the Time Discrete Metamorphosis Model - Existence and Applications(1705.04490)

May 12, 2017 math.NA
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.
• ### Collisional stability of localized Yb(${}^3\mathrm{P}_2$) atoms immersed in a Fermi sea of Li(1607.06180)

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 few-body systems separately.
• ### A one-dimensional liquid of fermions with tunable spin(1408.0928)

Aug. 5, 2014 cond-mat.quant-gas
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 one-dimensional (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, strongly-correlated 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 spin-1/2 Luttinger liquid. In the large-N limit, the system exhibits properties of a bosonic spinless liquid. Our results provide a testing ground for many-body theories and may lead to the observation of fundamental 1D effects.
• ### Optimal preparation of quantum states on an atom chip device(1405.6918)

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 Bose-Einstein condensate realized on an atom chip. As a possible application of our scheme, we apply it to improve the sensitivity of an atomic interferometer.
• ### Experimental realization of quantum zeno dynamics(1309.1552)

It is generally impossible to probe a quantum system without disturbing it. However, it is possible to exploit the back-action 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 Bose-Einstein condensate in a five-level 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 two-level 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.
• ### Distribution of Scattering Matrix Elements in Quantum Chaotic Scattering(1304.5284)

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 long-standing problem was the computation of the distribution of the off-diagonal scattering-matrix elements. We report here an exact solution to this problem and present analytical results for systems with preserved and with violated time-reversal 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.