We show how one can deterministically generate photonic matrix product states with high bond and physical dimensions with an atomic array if one has access to a Rydberg-blockade mechanism. We develop both a quantum gate and an optimal control approach to universally control the system and analyze the photon retrieval efficiency of atomic arrays. Comprehensive modeling of the system shows that our scheme is capable of generating a large number of entangled photons. We further develop a multi-port photon emission approach that can efficiently distribute entangled photons into free space in several directions, which can become a useful tool in future quantum networks.
Quantum entanglement involving coherent superpositions of macroscopically distinct states is among the most striking features of quantum theory, but its realization is challenging, since such states are extremely fragile. Using a programmable quantum simulator based on neutral atom arrays with interactions mediated by Rydberg states, we demonstrate the deterministic generation of Schrodinger cat states of the Greenberger-Horne-Zeilinger (GHZ) type with up to 20 qubits. Our approach is based on engineering the energy spectrum and using optimal control of the many-body system. We further demonstrate entanglement manipulation by using GHZ states to distribute entanglement to distant sites in the array, establishing important ingredients for quantum information processing and quantum metrology.
A powerful method to interface quantum light with matter is to propagate the light through an ensemble of atoms. Recently, a number of such interfaces have emerged, most prominently Rydberg ensembles, that enable strong nonlinear interactions between propagating photons. A largely open problem is whether these systems produce exotic many-body states of light and developing new tools to study propagation in the large photon number limit is highly desirable. Here, we provide a method based on a spin model that maps quasi one-dimensional (1D) light propagation to the dynamics of an open 1D interacting spin system, where all photon correlations are obtained from those of the spins. The spin dynamics in turn are numerically solved using the toolbox of matrix product states. We apply this formalism to investigate vacuum induced transparency, wherein the different photon number components of a pulse propagate with number-dependent group velocity and separate at output.
Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials, and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.
We study thermal states of strongly interacting quantum spin chains and prove that those can be represented in terms of convex combinations of matrix product states. Apart from revealing new features of the entanglement structure of Gibbs states our results provide a theoretical justification for the use of Whites algorithm of minimally entangled typical thermal states. Furthermore, we shed new light on time dependent matrix product state algorithms which yield hydrodynamical descriptions of the underlying dynamics.
We demonstrate that the optimal states in lossy quantum interferometry may be efficiently simulated using low rank matrix product states. We argue that this should be expected in all realistic quantum metrological protocols with uncorrelated noise and is related to the elusive nature of the Heisenberg precision scaling in presence of decoherence.