Source author record

J. Ignacio Cirac

J. Ignacio Cirac appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

130works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

130 published item(s)

preprint2026arXiv

Efficient Characterization of Coherent and Correlated Low-Degree Noise in Layers of Gates

We present a quantum process-tomography protocol based on a low-degree ansatz for the quantum channel, i.e. when it can be expressed as a fixed-degree polynomial in terms of Pauli operators. We demonstrate how to perform tomography of such channels with a logarithmic amount of effort relative to the size of the system, by employing random state preparation and measurements in the Pauli basis. We extend the applicability of the protocol to channels consisting of a layer of quantum gates with a polylogarithmic number of non-Clifford gates, followed by a low-degree noise channel. Rather than inverting the layer of quantum gates on the hardware-which would introduce additional errors-we instead carry out the inversion in classical postprocessing, while adding to the sample complexity a factor at most polynomial in system size. Numerical simulations support our theoretical findings and demonstrate the feasibility of our method.

preprint2026arXiv

Quantum Circuits for Matrix-Product Unitaries

Matrix-product unitaries (MPUs) are many-body unitary operators that, as a consequence of their tensor-network structure, preserve the entanglement area law in 1D systems. However, it is unknown how to implement an MPU as a quantum circuit since the individual tensors describing the MPU are not unitary. In this Letter, we show that a large class of MPUs can be implemented with a polynomial-depth quantum circuit. For an $N$-site MPU built from a repeated bulk tensor with open boundary, we explicitly construct a quantum circuit of polynomial depth $T = O(N^α)$ realizing the MPU, where the constant $α$ depends only on the bulk and boundary tensor and not the system size $N$. We show that this class includes nontrivial unitaries that generate long-range entanglement and, in particular, contains a large class of unitaries constructed from representations of $C^*$-weak Hopf algebras. Furthermore, we also adapt our construction to nonuniform translationally-varying MPUs and show that they can be implemented by a circuit of depth $O(N^β \, \mathrm{poly}\, D)$ where $β\le 1 + \log_2 \sqrt{D}/ s_{\min}$, with $D$ being the bond dimension and $s_{\min}$ the smallest nonzero Schmidt value of the normalized Choi state corresponding to the MPU.

preprint2022arXiv

Generation of photonic tensor network states with Circuit QED

We propose a circuit QED platform and protocol to generate microwave photonic tensor network states deterministically. We first show that using a microwave cavity as ancilla and a transmon qubit as emitter is a good platform to produce photonic matrix product states. The ancilla cavity combines a large controllable Hilbert space with a long coherence time, which we predict translates into a high number of entangled photons and states with a high bond dimension. Going beyond this paradigm, we then consider a natural generalization of this platform, in which several cavity-qubit pairs are coupled to form a chain. The photonic states thus produced feature a two-dimensional entanglement structure and can be interpreted as $\textit{radial plaquette}$ projected entangled pair states [Wei, Malz, and Cirac, Phys. Rev. Lett. 128, 010607 (2022)], which include many paradigmatic states, such as the broad class of isometric tensor network states, graph states, and string-net states.

preprint2022arXiv

Locality optimization for parent Hamiltonians of Tensor Networks

Tensor Network states form a powerful framework for both the analytical and numerical study of strongly correlated phases. Vital to their analytical utility is that they appear as the exact ground states of associated parent Hamiltonians, where canonical proof techniques guarantee a controlled ground space structure. Yet, while those Hamiltonians are local by construction, the known techniques often yield complex Hamiltonians which act on a rather large number of spins. In this paper, we present an algorithm to systematically simplify parent Hamiltonians, breaking them down into any given basis of elementary interaction terms. The underlying optimization problem is a semidefinite program, and thus the optimal solution can be found efficiently. Our method exploits a degree of freedom in the construction of parent Hamiltonians -- the excitation spectrum of the local terms -- over which it optimizes such as to obtain the best possible approximation. We benchmark our method on the AKLT model and the Toric Code model, where we show that the canonical parent Hamiltonians (acting on 3 or 4 and 12 sites, respectively) can be broken down to the known optimal 2-body and 4-body terms. We then apply our method to the paradigmatic Resonating Valence Bond (RVB) model on the kagome lattice. Here, the simplest previously known parent Hamiltonian acts on all the 12 spins on one kagome star. With our optimization algorithm, we obtain a vastly simpler Hamiltonian: We find that the RVB model is the exact ground state of a parent Hamiltonian whose terms are all products of at most four Heisenberg interactions, and whose range can be further constrained, providing a major improvement over the previously known 12-body Hamiltonian.

preprint2022arXiv

Long-range electron-electron interactions in quantum dot systems and applications in quantum chemistry

Long-range interactions play a key role in several phenomena of quantum physics and chemistry. To study these phenomena, analog quantum simulators provide an appealing alternative to classical numerical methods. Gate-defined quantum dots have been established as a platform for quantum simulation, but for those experiments the effect of long-range interactions between the electrons did not play a crucial role. Here we present the first detailed experimental characterization of long-range electron-electron interactions in an array of gate-defined semiconductor quantum dots. We demonstrate significant interaction strength among electrons that are separated by up to four sites, and show that our theoretical prediction of the screening effects matches well the experimental results. Based on these findings, we investigate how long-range interactions in quantum-dot arrays may be utilized for analog simulations of artificial quantum matter. We numerically show that about ten quantum dots are sufficient to observe binding for a one-dimensional $H_2$-like molecule. These combined experimental and theoretical results pave the way for future quantum simulations with quantum dot arrays and benchmarks of numerical methods in quantum chemistry.

preprint2022arXiv

Matrix product operator algebras I: representations of weak Hopf algebras and projected entangled pair states

Matrix Product Operators (MPOs) are tensor networks representing operators acting on 1D systems. They model a wide variety of situations, including communication channels with memory effects, quantum cellular automata, mixed states in 1D quantum systems, or holographic boundary models associated to 2D quantum systems. A scenario where MPOs have proven particularly useful is to represent algebras of non-trivial symmetries. Concretely, the boundary of both symmetry protected and topologically ordered phases in 2D quantum systems exhibit symmetries in the form of MPOs. In this paper, we develop a theory of MPOs as representations of algebraic structures. We establish a dictionary between algebra and MPO properties which allows to transfer results between both setups, covering the cases of pre-bialgebras, weak bialgebras, and weak Hopf algebras. We define the notion of pulling-through algebras, which abstracts the minimal requirements needed to define topologically ordered 2D tensor networks from MPO algebras. We show, as one of our main results, that any semisimple pivotal weak Hopf algebra is a pulling-trough algebra. We demonstrate the power of this framework by showing that they can be used to construct Kitaev's quantum double models for Hopf algebras solely from an MPO representation of the Hopf algebra, in the exact same way as MPO symmetries obtained from fusion categories can be used to construct Levin-Wen string-net models, and to explain all their topological features; it thus allows to describe both Kitaev and string-net models on the same formal footing.

preprint2022arXiv

Preparation and verification of tensor network states

We consider a family of tensor network states defined on regular lattices that come with a natural definition of an adiabatic path to prepare them. This family comprises relevant classes of states, such as injective Matrix Product and Projected Entangled-Pair States, and some corresponding to classical spin models. We show how uniform lower bounds to the gap of the parent Hamiltonian along the adiabatic trajectory can be efficiently computed using semi-definite programming. This allows one to check whether the adiabatic preparation can be performed efficiently with a scalable effort. We also derive a set of observables whose expectation values can be easily determined and that form a complete set, in the sense that they uniquely characterize the state. We identify a subset of those observables which can be efficiently computed if one has access to the quantum state and local measurements, and analyze how they can be used in verification procedures.

preprint2022arXiv

Sequential generation of projected entangled-pair states

We introduce plaquette projected entangled-pair states, a class of states in a lattice that can be generated by applying sequential unitaries acting on plaquettes of overlapping regions. They satisfy area-law entanglement, possess long-range correlations, and naturally generalize other relevant classes of tensor network states. We identify a subclass that can be more efficiently prepared in a radial fashion and that contains the family of isometric tensor network states [M. P. Zaletel and F. Pollmann, Phys. Rev. Lett. 124, 037201 (2020)]. We also show how this subclass can be efficiently prepared using an array of photon sources.

preprint2022arXiv

Spin-Holstein models in trapped-ion systems

In this work, we highlight how trapped-ion quantum systems can be used to study generalized Holstein models, and benchmark expensive numerical calculations. We study a particular spin-Holstein model that can be implemented with arrays of ions confined by individual microtraps, and that is closely related to the Holstein model of condensed matter physics, used to describe electron-phonon interactions. In contrast to earlier proposals, we focus on simulating many-electron systems and inspect the competition between charge-density wave order, fermion pairing and phase separation. In our numerical study, we employ a combination of complementary approaches, based on non-Gaussian variational ansatz states and matrix product states, respectively. We demonstrate that this hybrid approach outperforms standard density-matrix renormalization group calculations.

preprint2021arXiv

Chemistry of a light impurity in a Bose-Einstein condensate

In ultracold atomic gases, a unique interplay arises between phenomena known from condensed matter physics, few-body physics and chemistry. Similar to an electron in a solid, an impurity in an ultracold gas can get dressed by excitations from the medium, forming a quasiparticle called the polaron. We study how dressing of an impurity leads to a modification of its chemical reactivity. Using a Gaussian state variational method in the frame of the impurity, we demonstrate that three-body correlations lead to an instability of the polaron. This instability is connected to an Efimov resonance, but shifted to smaller interactions by many-body effects, showing that polaron formation stimulates Efimov physics and the associated chemistry.

preprint2021arXiv

Convergence guarantees for discrete mode approximations to non-Markovian quantum baths

Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. The non-Markovian environment is often approximated by discrete bosonic modes, thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been previously proposed, the resulting approximation lacks rigorous and general convergence guarantees. In this letter, we show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that this approximation error typically falls off polynomially with the number of modes. Our results lend rigor to classical and quantum algorithms for approximating non-Markovian dynamics.

preprint2021arXiv

From Efimov Physics to the Bose Polaron using Gaussian States

Since the Efimov effect was introduced in 1970, a detailed theoretical understanding of Efimov physics has been developed in the few-body context. However, it has proven to be challenging to describe the role Efimov-type correlations play in many-body systems such as quenched or collapsing Bose-Einstein condensates (BECs). To study the impact the Efimov effect can have in such scenarios, we consider a light impurity immersed in a weakly interacting BEC, forming a Bose polaron. In this case, the higher-order correlations are localized around the impurity, making it more feasible to develop a theoretical description. Specifically, we employ a Gaussian state variational Ansatz in the reference frame of the impurity, capable of both capturing the Efimov effect and the formation of the polaron cloud. We find that the Efimov effect leads to a cooperative binding of bosons to the impurity and the formation of a many-body bound state. As a result, the polaron is not the ground state, but rendered a metastable excited state which can decay into these Efimov clusters. While this decay is slow for small interaction strengths, it becomes more prominent as the attractive scattering length increases, up to the point where the polaron becomes completely unstable. This critical scattering length can be interpreted as a many-body shifted Efimov resonance, where the scattering of two excitations of the bath with the polaron can lead to bound state formation. Compared to the few-body case, the resonance is shifted to smaller attractive scattering lengths due to the participation of the polaron cloud in the cooperative binding process. This corresponds to an intriguing scenario of polaron-assisted chemistry, where many-body effects lead to enhanced signal of the chemical recombination process, which can be directly probed in state-of-the-art experiments.

preprint2021arXiv

Symmetries and local transformations of translationally invariant Matrix Product States

We determine the local symmetries and local transformation properties of translationally invariant matrix product states (MPS). We focus on physical dimension $d=2$ and bond dimension $D=3$ and use the procedure introduced in D. Sauerwein et al., Phys. Rev. Lett. 123, 170504 (2019) to determine all (including non--global) symmetries of those states. We identify and classify the stochastic local transformations (SLOCC) that are allowed among MPS. We scrutinize two very distinct sets of MPS and show the big diversity (also compared to the case $D=2$) occurring in both, their symmetries and the possible SLOCC transformations. These results reflect the variety of local properties of MPS, even if restricted to translationally invariant states with low bond dimension. Finally, we show that states with non-trivial local symmetries are of measure zero for $d = 2$ and $D > 3$.

preprint2020arXiv

Efficient description of many-body systems with Matrix Product Density Operators

Matrix Product States form the basis of powerful simulation methods for ground state problems in one dimension. Their power stems from the fact that they faithfully approximate states with a low amount of entanglement, the "area law". In this work, we establish the mixed state analogue of this result: We show that one-dimensional mixed states with a low amount of entanglement, quantified by the entanglement of purification, can be efficiently approximated by Matrix Product Density Operators (MPDOs). In combination with results establishing area laws for thermal states, this helps to put the use of MPDOs in the simulation of thermal states on a formal footing.

preprint2020arXiv

Exact dynamics in dual-unitary quantum circuits

We consider the class of dual-unitary quantum circuits in $1+1$ dimensions and introduce a notion of ``solvable'' matrix product states (MPSs), defined by a specific condition which allows us to tackle their time evolution analytically. We provide a classification of the latter, showing that they include certain MPSs of arbitrary bond dimension, and study analytically different aspects of their dynamics. For these initial states, we show that while any subsystem of size $\ell$ reaches infinite temperature after a time $t\propto \ell$, irrespective of the presence of conserved quantities, the light-cone of two-point correlation functions displays qualitatively different features depending on the ergodicity of the quantum circuit, defined by the behavior of infinite-temperature dynamical correlation functions. Furthermore, we study the entanglement spreading from such solvable initial states, providing a closed formula for the time evolution of the entanglement entropy of a connected block. This generalizes recent results obtained in the context of the self-dual kicked Ising model. By comparison, we also consider a family of non-solvable initial mixed states depending on one real parameter $β$, which, as $β$ is varied from zero to infinity, interpolate between the infinite temperature density matrix and arbitrary initial pure product states. We study analytically their dynamics for small values of $β$, and highlight the differences from the case of solvable MPSs.

preprint2020arXiv

Fermionic quantum cellular automata and generalized matrix product unitaries

We study matrix product unitary operators (MPUs) for fermionic one-dimensional (1D) chains. In stark contrast with the case of 1D qudit systems, we show that (i) fermionic MPUs do not necessarily feature a strict causal cone and (ii) not all fermionic Quantum Cellular Automata (QCA) can be represented as fermionic MPUs. We then introduce a natural generalization of the latter, obtained by allowing for an additional operator acting on their auxiliary space. We characterize a family of such generalized MPUs that are locality-preserving, and show that, up to appending inert ancillary fermionic degrees of freedom, any representative of this family is a fermionic QCA and viceversa. Finally, we prove an index theorem for generalized MPUs, recovering the recently derived classification of fermionic QCA in one dimension. As a technical tool for our analysis, we also introduce a graded canonical form for fermionic matrix product states, proving its uniqueness up to similarity transformations.

preprint2020arXiv

How much entanglement is needed to reduce the energy variance?

We explore the relation between the entanglement of a pure state and its energy variance for a local one dimensional Hamiltonian, as the system size increases. In particular, we introduce a construction which creates a matrix product state of arbitrarily small energy variance $δ^2$ for $N$ spins, with bond dimension scaling as $\sqrt{N} D_0^{1/δ}$, where $D_0>1$ is a constant. This implies that a polynomially increasing bond dimension is enough to construct states with energy variance that vanishes with the inverse of the logarithm of the system size. We run numerical simulations to probe the construction on two different models, and compare the local reduced density matrices of the resulting states to the corresponding thermal equilibrium. Our results suggest that the spatially homogeneous states with logarithmically decreasing variance, which can be constructed efficiently, do converge to the thermal equilibrium in the thermodynamic limit, while the same is not true if the variance remains constant.

preprint2020arXiv

Nondestructive photon counting in waveguide QED

Number-resolving single-photon detectors represent a key technology for a host of quantum optics protocols, but despite significant efforts, state-of-the-art devices are limited to few photons. In contrast, state-dependent atom counting in arrays can be done with extremely high fidelity up to hundreds of atoms. We show that in waveguide QED, the problem of photon counting can be reduced to atom counting, by entangling the photonic state with an atomic array in the collective number basis. This is possible as the incoming photons couple to collective atomic states and can be achieved by engineering a second decay channel of an excited atom to a metastable state. Our scheme is robust to disorder and finite Purcell factors, and its fidelity increases with atom number. Analyzing the state of the re-emitted photons, we further show that if the initial atomic state is a symmetric Dicke state, dissipation engineering can be used to implement a nondestructive photon-number measurement, in which the incident state is scattered into the waveguide unchanged. Our results generalize to related platforms, including superconducting qubits.

preprint2020arXiv

Quantum East model: localization, non-thermal eigenstates and slow dynamics

We study in detail the properties of the quantum East model, an interacting quantum spin chain inspired by simple kinetically-constrained models of classical glasses. Through a combination of analytics, exact diagonalization and tensor-network methods we show the existence of a transition, from a fast to a slow thermalization regime, which manifests itself throughout the spectrum. On the slow side, by exploiting the localization of the ground state and the form of the Hamiltonian, we explicitly construct a large (exponential in size) number of non-thermal states which become exact finite-energy-density eigenstates in the large-size limit, as expected for a true phase transition. A ``super-spin'' generalization allows us to find a further large class of area-law states proved to display very slow relaxation. These states retain memory of their initial conditions for extremely long times. Our numerical analysis reveals that the localization properties are not limited to the ground state and that} many eigenstates have large overlap with product states and can be approximated well by matrix product states at arbitrary energy densities. The mechanism that induces localization to the ground state, and hence the non-thermal behavior of the system, can be extended to a wide range of models including a number of simple spin chains. We discuss implications of our results for slow thermalization and non-ergodicity more generally in disorder-free systems with constraints and we give numerical evidence that these results may be extended to two dimensional systems.

preprint2020arXiv

Simulating 2+1d $\mathbb{Z}_3$ lattice gauge theory with iPEPS

We simulate a zero-temperature pure $\mathbb{Z}_3$ Lattice Gauge Theory in 2+1 dimensions by using an iPEPS (Infinite Projected Entangled-Pair State) ansatz for the ground state. Our results are therefore directly valid in the thermodynamic limit. They clearly show two distinct phases separated by a phase transition. We introduce an update strategy that enables plaquette terms and Gauss-law constraints to be applied as sequences of two-body operators. This allows the use of the most up-to-date iPEPS algorithms. From the calculation of spatial Wilson loops we are able to prove the existence of a confined phase. We show that with relatively low computational cost it is possible to reproduce crucial features of gauge theories. We expect that the strategy allows the extension of iPEPS studies to more general LGTs.

preprint2020arXiv

Wigner Crystals in Two-Dimensional Transition-Metal Dichalcogenides: Spin Physics and Readout

Wigner crystals are prime candidates for the realization of regular electron lattices under minimal requirements on external control and electronics. However, several technical challenges have prevented their detailed experimental investigation and applications to date. We propose an implementation of two-dimensional electron lattices for quantum simulation of Ising spin systems based on self-assembled Wigner crystals in transition-metal dichalcogenides. We show that these semiconductors allow for minimally invasive all-optical detection schemes of charge ordering and total spin. For incident light with optimally chosen beam parameters and polarization, we predict a strong dependence of the transmitted and reflected signals on the underlying lattice periodicity, thus revealing the charge order inherent in Wigner crystals. At the same time, the selection rules in transition-metal dichalcogenides provide direct access to the spin degree of freedom via Faraday rotation measurements.

preprint2019arXiv

Classification of Matrix-Product Unitaries with Symmetries

We prove that matrix-product unitaries (MPUs) with on-site unitary symmetries are completely classified by the (chiral) index and the cohomology class of the symmetry group $G$, provided that we can add trivial and symmetric ancillas with arbitrary on-site representations of $G$. If the representations in both system and ancillas are fixed to be the same, we can define symmetry-protected indices (SPIs) which quantify the imbalance in the transport associated to each group element and greatly refines the classification. These SPIs are stable against disorder and measurable in interferometric experiments. Our results lead to a systematic construction of two-dimensional Floquet symmetry-protected topological (SPT) phases beyond the standard classification, and thus shed new light on understanding nonequilibrium phases of quantum matter.

preprint2019arXiv

Probing thermalization through spectral analysis with matrix product operators

We combine matrix product operator techniques with Chebyshev polynomial expansions and present a method that is able to explore spectral properties of quantum many-body Hamiltonians. In particular, we show how this method can be used to probe thermalization of large spin chains without explicitly simulating their time evolution, as well as to compute full and local densities of states. The performance is illustrated with the examples of the Ising and PXP spin chains. For the non-integrable Ising chain, our findings corroborate the presence of thermalization for several initial states, well beyond what direct time-dependent simulations have been able to achieve so far.

preprint2019arXiv

Ultrafast molecular dynamics in terahertz-STM experiments: Theoretical analysis using Anderson-Holstein model

We analyze ultrafast tunneling experiments in which electron transport through a localized orbital is induced by a single cycle THz pulse. We include both electron-electron and electron-phonon interactions on the localized orbital using the Anderson-Holstein model and consider two possible filling factors, the singly occupied Kondo regime and the doubly occupied regime relevant to recent experiments with a pentacene molecule. Our analysis is based on variational non-Gaussian states and provides the accurate description of the degrees of freedom at very different energies, from the high microscopic energy scales to the Kondo temperature $T_K$. To establish the validity of the new method we apply this formalism to study the Anderson model in the Kondo regime in the absence of coupling to phonons. We demonstrate that it correctly reproduces key properties of the model, including the screening of the impurity spin, formation of the resonance at the Fermi energy, and a linear conductance of $2e^2/h$. We discuss the suppression of the Kondo resonance by the electron-phonon interaction on the impurity site. When analyzing THz STM experiments we compute the time dependence of the key physical quantities, including current, the number of electrons on the localized orbital, and the number of excited phonons. We find long-lived oscillations of the phonon that persist long after the end of the pulse. We compare the results for the interacting system to the non-interacting resonant level model.

preprint2016arXiv

Dissipative Long-Range Entanglement Generation between Electronic Spins

We propose a scheme for deterministic generation and long-term stabilization of entanglement between two electronic spin qubits confined in spatially separated quantum dots. Our approach relies on an electronic quantum bus, consisting either of quantum Hall edge channels or surface acoustic waves, that can mediate long-range coupling between localized spins over distances of tens of micrometers. Since the entanglement is actively stabilized by dissipative dynamics, our scheme is inherently robust against noise and imperfections.

preprint2016arXiv

Edge states for the Kalmeyer-Laughlin wave function

We study lattice wave functions obtained from the SU(2)$_1$ Wess-Zumino-Witten conformal field theory. Following Moore and Read's construction, the Kalmeyer-Laughlin fractional quantum Hall state is defined as a correlation function of primary fields. By an additional insertion of Kac-Moody currents, we associate a wave function to each state of the conformal field theory. These wave functions span the complete Hilbert space of the lattice system. On the cylinder, we study global properties of the lattice states analytically and correlation functions numerically using a Metropolis Monte Carlo method. By comparing short-range bulk correlations, numerical evidence is provided that the states with one current operator represent edge states in the thermodynamic limit. We show that the edge states with one Kac-Moody current of lowest order have a good overlap with low-energy excited states of a local Hamiltonian, for which the Kalmeyer-Laughlin state approximates the ground state. For some states, exact parent Hamiltonians are derived on the cylinder. These Hamiltonians are SU(2) invariant and nonlocal with up to four-body interactions.

preprint2016arXiv

Lattice effects on Laughlin wave functions and parent Hamiltonians

We investigate lattice effects on wave functions that are lattice analogues of bosonic and fermionic Laughlin wave functions with number of particles per flux $ν=1/q$ in the Landau levels. These wave functions are defined analytically on lattices with $μ$ particles per lattice site, where $μ$ may be different than $ν$. We give numerical evidence that these states have the same topological properties as the corresponding continuum Laughlin states for different values of $q$ and for different fillings $μ$. These states define, in particular, particle-hole symmetric lattice Fractional Quantum Hall states when the lattice is half-filled. On the square lattice it is observed that for $q\leq 4$ this particle-hole symmetric state displays the topological properties of the continuum Laughlin state at filling fraction $ν=1/q$, while for larger $q$ there is a transition towards long-range ordered anti-ferromagnets. This effect does not persist if the lattice is deformed from a square to a triangular lattice, or on the Kagome lattice, in which case the topological properties of the state are recovered. We then show that changing the number of particles while keeping the expression of these wave functions identical gives rise to edge states that have the same correlations in the bulk as the reference lattice Laughlin states but a different density at the edge. We derive an exact parent Hamiltonian for which all these edge states are ground states with different number of particles. In addition this Hamiltonian admits the reference lattice Laughlin state as its unique ground state of filling factor $1/q$. Parent Hamiltonians are also derived for the lattice Laughlin states at other fillings of the lattice, when $μ\leq 1/q$ or $μ\geq 1-1/q$ and when $q=4$ also at half-filling.

preprint2016arXiv

Projected Entangled Pair States with non-Abelian gauge symmetries: an SU(2) study

Over the last years, Projected Entangled Pair States have demonstrated great power for the study of many body systems, as they naturally describe ground states of gapped many body Hamiltonians, and suggest a constructive way to encode and classify their symmetries. The PEPS study is not only limited to global symmetries, but has also been extended and applied for local symmetries, allowing to use them for the description of states in lattice gauge theories. In this paper we discuss PEPS with a local, SU(2) gauge symmetry, and demonstrate the use of PEPS features and techniques for the study of a simple family of many body states with a non-Abelian gauge symmetry. We present, in particular, the construction of fermionic PEPS able to describe both two-color fermionic matter and the degrees of freedom of an SU(2) gauge field with a suitable truncation.

preprint2016arXiv

Rapid adiabatic preparation of injective PEPS and Gibbs states

We propose a quantum algorithm for many-body state preparation. It is especially suited for injective PEPS and thermal states of local commuting Hamiltonians on a lattice. We show that for a uniform gap and sufficiently smooth paths, an adiabatic runtime and circuit depth of $O(\operatorname{polylog}N)$ can be achieved for $O(N)$ spins. This is an almost exponential improvement over previous bounds. The total number of elementary gates scales as $O(N\operatorname{polylog}N)$. This is also faster than the best known upper bound of $O(N^2)$ on the mixing times of Monte Carlo Markov chain algorithms for sampling classical systems in thermal equilibrium.

preprint2016arXiv

Systematic construction of density functionals based on matrix product state computations

We propose a systematic procedure for the approximation of density functionals in density functional theory that consists of two parts. First, for the efficient approximation of a general density functional, we introduce an efficient ansatz whose non-locality can be increased systematically. Second, we present a fitting strategy that is based on systematically increasing a reasonably chosen set of training densities. We investigate our procedure in the context of strongly correlated fermions on a one-dimensional lattice in which we compute accurate training densities with the help of matrix product states. Focusing on the exchange-correlation energy, we demonstrate how an efficient approximation can be found that includes and systematically improves beyond the local density approximation. Importantly, this systematic improvement is shown for target densities that are quite different from the training densities.

preprint2016arXiv

The multi-flavor Schwinger model with chemical potential - Overcoming the sign problem with Matrix Product States

During recent years there has been an increasing interest in the application of matrix product states, and more generally tensor networks, to lattice gauge theories. This non-perturbative method is sign problem free and has already been successfully used to compute mass spectra, thermal states and phase diagrams, as well as real-time dynamics for Abelian and non-Abelian gauge models. In previous work we showed the suitability of the method to explore the zero-temperature phase structure of the multi-flavor Schwinger model at non-zero chemical potential, a regime where the conventional Monte Carlo approach suffers from the sign problem. Here we extend our numerical study by looking at the spatially resolved chiral condensate in the massless case. We recover spatial oscillations, similar to the theoretical predictions for the single-flavor case, with a chemical potential dependent frequency and an amplitude approximately given by the homogeneous zero density condensate value.

preprint2016arXiv

Ultrashort Pulses for Far-Field Nanoscopy

We show that ultrashort pulses can be focused, in a particular instant, to a spot size given by the wavelength associated with its spectral width. For attosecond pulses this spot size is within the nanometer scale. Then we show that a two-level system can be left excited after interacting with an ultrashort pulse whose spectral width is larger than the transition frequency, and that the excitation probability depends not on the field amplitude but on the field intensity. The latter makes the excitation profile have the same spot size as the ultrashort pulse. This unusual phenomenon is caused by quantum electrodynamics in the ultrafast light-matter interaction regime since the usually neglected counterrotating terms describing the interaction with the free electromagnetic modes are crucial for making the excitation probability nonzero and depend on the field intensity. These results suggest that a train of coherent attosecond pulses could be used to excite fluorescent markers with nanoscale resolution. The detection of the light emitted after fluorescence -or any other method used to detect the excitation- could then lead to a new scheme for far-field light nanoscopy.

preprint2015arXiv

Chiral topological spin liquids with projected entangled pair states

Topological chiral phases are ubiquitous in the physics of the Fractional Quantum Hall Effect. Non-chiral topological spin liquids are also well known. Here, using the framework of projected entangled pair states (PEPS), we construct a family of chiral spin liquids on the square lattice which are generalized spin-1/2 Resonating Valence Bond (RVB) states obtained from deformed local tensors with $d+i\, d$ symmetry. On a cylinder, we construct four topological sectors with even or odd number of spinons on the boundary and even or odd number of ($\mathbb{Z}_2$) fluxes penetrating the cylinder which, we argue, remain orthogonal in the limit of infinite perimeter. The analysis of the transfer matrix provides evidence of short-range (long-range) triplet (singlet) correlations as for the critical (non-chiral) RVB state. The Entanglement Spectrum exhibits chiral edge modes, which we confront to predictions of Conformal Field Theory, and the corresponding Entanglement Hamiltonian is shown to be long ranged.

preprint2015arXiv

Efficient variational diagonalization of fully many-body localized Hamiltonians

We introduce a unitary matrix-product operator (UMPO) based variational method that approximately finds all the eigenstates of fully many-body localized (fMBL) one-dimensional Hamiltonians. The computational cost of the variational optimization scales linearly with system size for a fixed bond dimension of the UMPO ansatz. We demonstrate the usefulness of our approach by considering the Heisenberg chain in a strongly disordered magnetic field for which we compare the approximation to exact diagonalization results.

preprint2015arXiv

Exact parent Hamiltonians of bosonic and fermionic Moore-Read states on lattices and local models

We introduce a family of strongly-correlated spin wave functions on arbitrary spin-1/2 and spin-1 lattices in one and two dimensions. These states are lattice analogues of Moore-Read states of particles at filling fraction 1/q, which are non-Abelian Fractional Quantum Hall states in 2D. One parameter enables us to perform an interpolation between the continuum limit, where the states become continuum Moore-Read states of bosons (odd q) and fermions (even q), and the lattice limit. We show numerical evidence that the topological entanglement entropy stays the same along the interpolation for some of the states we introduce in 2D, which suggests that the topological properties of the lattice states are the same as in the continuum, while the 1D states are critical states. We then derive exact parent Hamiltonians for these states on lattices of arbitrary size. By deforming these parent Hamiltonians, we construct local Hamiltonians that stabilize some of the states we introduce in 1D and in 2D.

preprint2015arXiv

Fermionic Projected Entangled Pair States and Local U(1) Gauge Theories

Tensor networks, and in particular Projected Entangled Pair States (PEPS), are a powerful tool for the study of quantum many body physics, thanks to both their built-in ability of classifying and studying symmetries, and the efficient numerical calculations they allow. In this work, we introduce a way to extend the set of symmetric PEPS in order to include local gauge invariance and investigate lattice gauge theories with fermionic matter. To this purpose, we provide as a case study and first example, the construction of a fermionic PEPS, based on Gaussian schemes, invariant under both global and local U(1) gauge transformations. The obtained states correspond to a truncated U(1) lattice gauge theory in 2 + 1 dimensions, involving both the gauge field and fermionic matter. For the global symmetry (pure fermionic) case, these PEPS can be studied in terms of spinless fermions subject to a p-wave superconducting pairing. For the local symmetry (fermions and gauge fields) case, we find confined and deconfined phases in the pure gauge limit, and we discuss the screening properties of the phases arising in the presence of dynamical matter.

preprint2015arXiv

Gauging quantum states: from global to local symmetries in many-body systems

We present an operational procedure to transform global symmetries into local symmetries at the level of individual quantum states, as opposed to typical gauging prescriptions for Hamiltonians or Lagrangians. We then construct a compatible gauging map for operators, which preserves locality and reproduces the minimal coupling scheme for simple operators. By combining this construction with the formalism of projected entangled-pair states (PEPS), we can show that an injective PEPS for the matter fields is gauged into a G-injective PEPS for the combined gauge-matter system, which potentially has topological order. We derive the corresponding parent Hamiltonian, which is a frustration free gauge theory Hamiltonian closely related to the Kogut-Susskind Hamiltonian at zero coupling constant. We can then introduce gauge dynamics at finite values of the coupling constant by applying a local filtering operation. This scheme results in a low-parameter family of gauge invariant states of which we can accurately probe the phase diagram, as we illustrate by studying a Z2 gauge theory with Higgs matter.

preprint2015arXiv

Multi-photon Scattering Theory and Generalized Master Equations

We develop a scattering theory to investigate the multi-photon transmission in a one-dimensional waveguide in the presence of quantum emitters. It is based on a path integral formalism, uses displacement transformations, and does not require the Markov approximation. We obtain the full time-evolution of the global system, including the emitters and the photonic field. Our theory allows us to compute the transition amplitude between arbitrary initial and final states, as well as the S-matrix of the asymptotic in- and out- states. For the case of few incident photons in the waveguide, we also re-derive a generalized master equation in the Markov limit. We compare the predictions of the developed scattering theory and that with the Markov approximation. We illustrate our methods with five examples of few-photon scattering: (i) by a two-level emitter, (ii) in the Jaynes-Cummings model; (iii) by an array of two-level emitters; (iv) by a two-level emitter in the half-end waveguide; (v) by an array of atoms coupled to Rydberg levels. In the first two, we show the application of the scattering theory in the photon scattering by a single emitter, and examine the correctness of our theory with the well-known results. In the third example, we analyze the condition of the Markov approximation for the photon scattering in the array of emitters. In the forth one, we show how a quantum emitter can generate entanglement of out-going photons. Finally, we highlight the interplay between the phenomenon of electromagnetic-induced transparency and the Rydberg interaction, and show how this results in a rich variety of possibilities in the quantum statistics of the scattering photons.

preprint2015arXiv

Non-Abelian string breaking phenomena with Matrix Product States

Using matrix product states, we explore numerically the phenomenology of string breaking in a non-Abelian lattice gauge theory, namely 1+1 dimensional SU(2). The technique allows us to study the static potential between external heavy charges, as traditionally explored by Monte Carlo simulations, but also to simulate the real-time dynamics of both static and dynamical fermions, as the latter are fully included in the formalism. We propose a number of observables that are sensitive to the presence or breaking of the flux string, and use them to detect and characterize the phenomenon in each of these setups.

preprint2015arXiv

Quantum dynamics of propagating photons with strong interactions: a generalized input-output formalism

There has been rapid development of systems that yield strong interactions between freely propagating photons in one dimension via controlled coupling to quantum emitters. This raises interesting possibilities such as quantum information processing with photons or quantum many-body states of light, but treating such systems generally remains a difficult task theoretically. Here, we describe a novel technique in which the dynamics and correlations of a few photons can be exactly calculated, based upon knowledge of the initial photonic state and the solution of the reduced effective dynamics of the quantum emitters alone. We show that this generalized "input-output" formalism allows for a straightforward numerical implementation regardless of system details, such as emitter positions, external driving, and level structure. As a specific example, we apply our technique to show how atomic systems with infinite-range interactions and under conditions of electromagnetically induced transparency enable the selective transmission of correlated multi-photon states.

preprint2015arXiv

Quantum Simulations of Lattice Gauge Theories using Ultracold Atoms in Optical Lattices

Can high energy physics be simulated by low-energy, non-relativistic, many-body systems, such as ultracold atoms? Such ultracold atomic systems lack the type of symmetries and dynamical properties of high energy physics models: in particular, they manifest neither local gauge invariance nor Lorentz invariance, which are crucial properties of the quantum field theories which are the building blocks of the standard model of elementary particles. However, it turns out, surprisingly, that there are ways to configure atomic system to manifest both local gauge invariance and Lorentz invariance. In particular, local gauge invariance can arise either as an effective, low energy, symmetry, or as an "exact" symmetry, following from the conservation laws in atomic interactions. Hence, one could hope that such quantum simulators may lead to new type of (table-top) experiments, that shall be used to study various QCD phenomena, as the confinement of dynamical quarks, phase transitions, and other effects, which are inaccessible using the currently known computational methods. In this report, we review the Hamiltonian formulation of lattice gauge theories, and then describe our recent progress in constructing quantum simulation of Abelian and non-Abelian lattice gauge theories in 1+1 and 2+1 dimensions using ultracold atoms in optical lattices.

preprint2015arXiv

Slowest local operators in quantum spin chains

We numerically construct slowly relaxing local operators in a nonintegrable spin-1/2 chain. Restricting the support of the operator to $M$ consecutive spins along the chain, we exhaustively search for the operator that minimizes the Frobenius norm of the commutator with the Hamiltonian. We first show that the Frobenius norm bounds the time scale of relaxation of the operator at high temperatures. We find operators with significantly slower relaxation than the slowest simple "hydrodynamic" mode due to energy diffusion. Then, we examine some properties of the nontrivial slow operators. Using both exhaustive search and tensor network techniques, we find similar slowly relaxing operators for a Floquet spin chain; this system is hydrodynamically "trivial", with no conservation laws restricting their dynamics. We argue that such slow relaxation may be a generic feature following from locality and unitarity.

preprint2015arXiv

Universal Quantum Transducers based on Surface Acoustic Waves

We propose a universal, on-chip quantum transducer based on surface acoustic waves in piezo-active materials. Because of the intrinsic piezoelectric (and/or magnetostrictive) properties of the material, our approach provides a universal platform capable of coherently linking a broad array of qubits, including quantum dots, trapped ions, nitrogen-vacancy centers or superconducting qubits. The quantized modes of surface acoustic waves lie in the gigahertz range, can be strongly confined close to the surface in phononic cavities and guided in acoustic waveguides. We show that this type of surface acoustic excitations can be utilized efficiently as a quantum bus, serving as an on-chip, mechanical cavity-QED equivalent of microwave photons and enabling long-range coupling of a wide range of qubits.

preprint2015arXiv

Variational Matrix Product Operators for the Steady State of Dissipative Quantum Systems

We present a new variational method, based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations for several dissipative spin models over a wide range of parameters illustrate the performance of the method and show that indeed the stationary state is often well described by a MPO of very moderate dimensions.

preprint2014arXiv

Algorithms for finite Projected Entangled Pair States

Projected Entangled Pair States (PEPS) are a promising ansatz for the study of strongly correlated quantum many-body systems in two dimensions. But due to their high computational cost, developing and improving PEPS algorithms is necessary to make the ansatz widely usable in practice. Here we analyze several algorithmic aspects of the method. On the one hand, we quantify the connection between the correlation length of the PEPS and the accuracy of its approximate contraction, and discuss how purifications can be used in the latter. On the other, we present algorithmic improvements for the update of the tensor that introduce drastic gains in the numerical conditioning and the efficiency of the algorithms. Finally, the state-of-the-art general PEPS code is benchmarked with the Heisenberg and quantum Ising models on lattices of up to $21 \times 21$ sites.

preprint2014arXiv

Approximating Gibbs states of local Hamiltonians efficiently with PEPS

We analyze the error of approximating Gibbs states of local quantum spin Hamiltonians on lattices with Projected Entangled Pair States (PEPS) as a function of the bond dimension ($D$), temperature ($β^{-1}$), and system size ($N$). First, we introduce a compression method in which the bond dimension scales as $D=e^{O(\log^2(N/ε))}$ if $β<O(\log (N))$. Second, building on the work of Hastings [Phys. Rev. B 73, 085115 (2006)], we derive a polynomial scaling relation, $D=\left(N/ε\right)^{O(β)}$. This implies that the manifold of PEPS forms an efficient representation of Gibbs states of local quantum Hamiltonians. From those bounds it also follows that ground states can be approximated with $D=N^{O(\log(N))}$ whenever the density of states only grows polynomially in the system size. All results hold for any spatial dimension of the lattice.

preprint2014arXiv

Chiral projected entangled-pair state with topological order

We show that projected entangled-pair states (PEPS) can describe chiral topologically ordered phases. For that, we construct a simple PEPS for spin-1/2 particles in a two-dimensional lattice. We reveal a symmetry in the local projector of the PEPS that gives rise to the global topological character. We also extract characteristic quantities of the edge conformal field theory using the bulk-boundary correspondence.

preprint2014arXiv

Construction of spin models displaying quantum criticality from quantum field theory

We provide a method for constructing finite temperature states of one-dimensional spin chains displaying quantum criticality. These models are constructed using correlators of products of quantum fields and have an analytical purification. Their properties can be investigated by Monte-Carlo simulations, which enable us to study the low-temperature phase diagram and to show that it displays a region of quantum criticality. The mixed states obtained are shown to be close to the thermal state of a simple nearest neighbour Hamiltonian.

preprint2014arXiv

Frustration free gapless Hamiltonians for Matrix Product States

For every Matrix Product State (MPS) one can always construct a so-called parent Hamiltonian. This is a local, frustration free, Hamiltonian which has the MPS as ground state and is gapped. Whenever that parent Hamiltonian has a degenerate ground state (the so-called non-injective case), we construct another 'uncle' Hamiltonian which is local and frustration free but gapless, and its spectrum is $\R^+$. The construction is obtained by linearly perturbing the matrices building up the state in a random direction, and then taking the limit where the perturbation goes to zero. For MPS where the parent Hamiltonian has a unique ground state (the so-called injective case) we also build such uncle Hamiltonian with the same properties in the thermodynamic limit.

preprint2014arXiv

Long-Distance Transfer and Routing of Static Magnetic Fields

We show how the static magnetic field of a finite source can be transferred and routed to arbitrary long distances. This is achieved by using transformation optics, which results in a device made of a material with a highly anisotropic magnetic permeability. We show that a simplified version of the device, made by a superconducting-ferromagnet hybrid, also leads to an excellent transfer of the magnetic field. The latter is demonstrated with a proof-of-principle experiment where a ferromagnet tube coated with a superconductor improves the transfer of static magnetic fields with respect to conventional methods by a 400\% factor over distances of 14cm.

preprint2014arXiv

Nuclear Spin Dynamics in Double Quantum Dots: Multi-Stability, Dynamical Polarization, Criticality and Entanglement

We theoretically study the nuclear spin dynamics driven by electron transport and hyperfine interaction in an electrically-defined double quantum dot (DQD) in the Pauli-blockade regime. We derive a master-equation-based framework and show that the coupled electron-nuclear system displays an instability towards the buildup of large nuclear spin polarization gradients in the two quantum dots. In the presence of such inhomogeneous magnetic fields, a quantum interference effect in the collective hyperfine coupling results in sizable nuclear spin entanglement between the two quantum dots in the steady state of the evolution. We investigate this effect using analytical and numerical techniques, and demonstrate its robustness under various types of imperfections.

preprint2014arXiv

Optical lattice implementation scheme of a bosonic topological model with fermionic atoms

We present a scheme to implement a Fermi-Hubbard-like model in ultracold atoms in optical lattices and analyze the topological features of its ground state. In particular, we show that the ground state for appropriate parameters has a large overlap with a lattice version of the bosonic Laughlin state at filling factor one half. The scheme utilizes laser assisted and normal tunnelling in a checkerboard optical lattice. The requirements on temperature, interactions, and hopping strengths are similar to those needed to observe the Néel antiferromagnetic ordering in the standard Fermi-Hubbard model in the Mott insulating regime.

preprint2014arXiv

Quantum simulation of the Schwinger model: A study of feasibility

We analyze some crucial questions regarding the practical feasibility of quantum simulation for lattice gauge models. Our analysis focuses on two models suitable for the quantum simulation of the Schwinger Hamiltonian, or QED in 1+1 dimensions, which we investigate numerically using tensor networks. In particular, we explore the effect of representing the gauge degrees of freedom with finite-dimensional systems and show that the results converge rapidly; thus even with small dimensions it is possible to obtain a reasonable accuracy. We also discuss the time scales required for the adiabatic preparation of the interacting vacuum state and observe that for a suitable ramping of the interaction the required time is almost insensitive to the system size and the dimension of the physical systems. Finally, we address the possible presence of noninvariant terms in the Hamiltonian that is realized in the experiment and show that for low levels of noise it is still possible to achieve a good precision for some ground-state observables, even if the gauge symmetry is not exact in the implemented model.

preprint2014arXiv

Robustness of quantum memories based on Majorana zero modes

We analyze the rate at which quantum information encoded in zero-energy Majorana modes is lost in the presence of perturbations. We show that information can survive for times that scale exponentially with the size of the chain both in the presence of quenching and time-dependent quadratic dephasing perturbations, even when the latter have spectral components above the system's energy gap. The origin of the robust storage, namely the fact that a sudden quench affects in the same way both parity sectors of the original spectrum, is discussed, together with the memory performance at finite temperatures and in the presence of particle exchange with a bath.

preprint2014arXiv

RVB superconductors with fermionic projected entangled pair states

We construct a family of simple fermionic projected entangled pair states (fPEPS) on the square lattice with bond dimension $D=3$ which are exactly hole-doped resonating valence bond (RVB) wavefunctions with short-range singlet bonds. Under doping the insulating RVB spin liquid evolves immediately into a superconductor with mixed $d+is$ pairing symmetry whose pair amplitude grows as the square-root of the doping. The relative weight between $s$-wave and $d$-wave components can be controlled by a single variational parameter $c$. We optimize our ansatz w.r.t. $c$ for the frustrated $t-J_1-J_2$ model (including both nearest and next-nearest neighbor antiferromagnetic interactions $J_1$ and $J_2$, respectively) for $J_2\simeq J_1/2$ and obtain an energy very close to the infinite-PEPS state (using full update optimization and same bond dimension). The orbital symmetry of the optimized RVB superconductor has predominant d-wave character, although we argue a residual (complex s-wave) time reversal symmetry breaking component should always be present. Connections of the results to the physics of superconducting cuprates and pnictides are outlined.

preprint2014arXiv

Symmetries and boundary theories for chiral Projected Entangled Pair States

We investigate the topological character of lattice chiral Gaussian fermionic states in two dimensions possessing the simplest descriptions in terms of projected entangled-pair states (PEPS). They are ground states of two different kinds of Hamiltonians. The first one, $\mathcal H_\mathrm{ff}$, is local, frustration-free, and gapless. It can be interpreted as describing a quantum phase transition between different topological phases. The second one, $\mathcal H_\mathrm{fb}$ is gapped, and has hopping terms scaling as $1/r^3$ with the distance $r$. The gap is robust against local perturbations, which allows us to define a Chern number for the PEPS. As for (non-chiral) topological PEPS, the non-trivial topological properties can be traced down to the existence of a symmetry in the virtual modes that are used to build the state. Based on that symmetry, we construct string-like operators acting on the virtual modes that can be continuously deformed without changing the state. On the torus, the symmetry implies that the ground state space of the local parent Hamiltonian is two-fold degenerate. By adding a string wrapping around the torus one can change one of the ground states into the other. We use the special properties of PEPS to build the boundary theory and show how the symmetry results in the appearance of chiral modes, and a universal correction to the area law for the zero Rényi entropy.

preprint2014arXiv

The temperature dependence of the chiral condensate in the Schwinger model with Matrix Product States

We present our recent results for the tensor network (TN) approach to lattice gauge theories. TN methods provide an efficient approximation for quantum many-body states. We employ TN for one dimensional systems, Matrix Product States, to investigate the 1-flavour Schwinger model. In this study, we compute the chiral condensate at finite temperature. From the continuum extrapolation, we obtain the chiral condensate in the high temperature region consistent with the analytical calculation by Sachs and Wipf.

preprint2014arXiv

Unifying Projected Entangled Pair States contractions

The approximate contraction of a Projected Entangled Pair States (PEPS) tensor network is a fundamental ingredient of any PEPS algorithm, required for the optimization of the tensors in ground state search or time evolution, as well as for the evaluation of expectation values. An exact contraction is in general impossible, and the choice of the approximating procedure determines the efficiency and accuracy of the algorithm. We analyze different previous proposals for this approximation, and show that they can be understood via the form of their environment, i.e. the operator that results from contracting part of the network. This provides physical insight into the limitation of various approaches, and allows us to introduce a new strategy, based on the idea of clusters, that unifies previous methods. The resulting contraction algorithm interpolates naturally between the cheapest and most imprecise and the most costly and most precise method. We benchmark the different algorithms with finite PEPS, and show how the cluster strategy can be used for both the tensor optimization and the calculation of expectation values. Additionally, we discuss its applicability to the parallelization of PEPS and to infinite systems (iPEPS).

preprint2013arXiv

Dissipative Dynamics and Phase Transitions in Fermionic Systems

We study abrupt changes in the dynamics and/or steady state of fermionic dissipative systems produced by small changes of the system parameters. Specifically, we consider open fermionic systems whose dynamics is described by master equations that are quadratic (and, under certain conditions, quartic) in creation and annihilation operators. We analyze both phase transitions in steady state, as well as "dynamical transitions". The latter are characterized by abrupt changes in the rate at which the system asymptotically approaches the steady state. We illustrate our general findings with relevant examples of fermionic (and, equivalently, spin) systems, and show that they can be realized in ion chains.

preprint2013arXiv

Dissipative spin chains: Implementation with cold atoms and steady-state properties

We propose a quantum optical implementation of a class of dissipative spin systems, including the XXZ and Ising model, with ultra-cold atoms in optical lattices. Employing the motional degree of freedom of the atoms and detuned Raman transitions we show how to obtain engineerable dissipation and a tunable transversal magnetic field, enabling the study of the dynamics and steady-states of dissipative spin models. As an example of effects made accessible this way, we consider small spin chains and weak dissipation and show by numerical simulation that steady-state expectation values display pronounced peaks at certain critical system parameters. We show that this effect is related to degeneracies in the Hamiltonian and derive a sufficient condition for its occurrence.

preprint2013arXiv

Field-induced superfluids and Bose liquids in Projected Entangled Pair States

In two-dimensional incompressible quantum spin liquids, a large enough magnetic field generically induces "doping" of polarized S=1 triplons or S=1/2 spinons. We review a number of cases such as spin-3/2 AKLT or spin-1/2 Resonating Valence Bond (RVB) liquids where the Projected Entangled Pair States (PEPS) framework provides very simple and comprehensive pictures. On the bipartite honeycomb lattice, simple PEPS can describe Bose condensed triplons (AKLT) or spinons (RVB) superfluids with transverse staggered (Néel) magnetic order. On the Kagome lattice, doping the RVB state with deconfined spinons or triplons (i.e. spinon bound pairs) yields uncondensed Bose liquids preserving U(1) spin-rotation symmetry. We find that spinon (triplon) doping destroys (preserves) the topological Z_2 symmetry of the underlying RVB state. We also find that spinon doping induces longer range interactions in the entanglement Hamiltonian, suggesting the emergence of (additive) log-corrections to the entanglement entropy.

preprint2013arXiv

Fractional quantum Hall states in lattices: Local models and physical implementation

The fractional quantum Hall (FQH) effect is one of the most striking phenomena in condensed matter physics. It is described by a simple Laughlin wavefunction and has been thoroughly studied both theoretically and experimentally. In lattice systems, however, such an effect has not been observed, there are few simple models displaying it, and only few mechanisms leading to it are known. Here we propose a new way of constructing lattice Hamiltonians with local interactions and FQH like ground states. In particular, we obtain a spin 1/2 model with a bosonic Laughlin like ground state, displaying a variety of topological features. We also demonstrate how such a model naturally emerges out of a Fermi-Hubbard like model at half filling, in which the kinetic energy part possesses bands with nonzero Chern number, and we show how this model can be implemented in an optical lattice setup with present technology.

preprint2013arXiv

Ground States of Fermionic lattice Hamiltonians with Permutation Symmetry

We study the ground states of lattice Hamiltonians that are invariant under permutations, in the limit where the number of lattice sites, N -> \infty. For spin systems, these are product states, a fact that follows directly from the quantum de Finetti theorem. For fermionic systems, however, the problem is very different, since mode operators acting on different sites do not commute, but anti-commute. We construct a family of fermionic states, \cal{F}, from which such ground states can be easily computed. They are characterized by few parameters whose number only depends on M, the number of modes per lattice site. We also give an explicit construction for M=1,2. In the first case, \cal{F} is contained in the set of Gaussian states, whereas in the second it is not. Inspired by that constructions, we build a set of fermionic variational wave functions, and apply it to the Fermi-Hubbard model in two spatial dimensions, obtaining results that go beyond the generalized Hartree-Fock theory.

preprint2013arXiv

Lattice Laughlin States of Bosons and Fermions at Filling Fractions $1/q$

We introduce a two-parameter family of strongly-correlated wave functions for bosons and fermions in lattices. One parameter, $q$, is connected to the filling fraction. The other one, $η$, allows us to interpolate between the lattice limit ($η=1$) and the continuum limit ($η\to 0^+$) of families of states appearing in the context of the fractional quantum Hall effect or the Calogero-Sutherland model. We give evidence that the main physical properties along the interpolation remain the same. Finally, in the lattice limit, we derive parent Hamiltonians for those wave functions and in 1D, we determine part of the low energy spectrum.

preprint2013arXiv

Matrix Product States for Lattice Field Theories

The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system. Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method. Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques. While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle. Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics. We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques. Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner and to accurately extract the mass spectrum. In this work we extend the analysis to the determination of the chiral condensate, both for massless and massive fermions. The method allows for accurate finite size and continuum limit extrapolations and produces remarkably precise results, thus showing the feasibility of these techniques for gauge theory problems.

preprint2013arXiv

Optomechanics assisted with a qubit: From dissipative state preparation to many-body physics

We propose and analyze nonlinear optomechanical protocols that can be implemented by adding a single atom to an optomechanical cavity. In particular, we show how to engineer the environment in order to dissipatively prepare the mechanical oscillator in a superposition of Fock states with fidelity close to one. Furthermore, we discuss how a single atom in a cavity with several mechanical oscillators can be exploited to realize nonlinear many-body physics by stroboscopically driving the mechanical oscillators. We show how to prepare non-classical many-body states by either applying coherent protocols or engineering dissipation. The analysis of the protocols is carried out using a perturbation theory for degenerate Liouvillians and numerical tools. Our results apply to other systems where a qubit is coupled to a mechanical oscillator via a bosonic mode, e.g., in cavity quantum electromechanics.

preprint2013arXiv

Purifications of multipartite states: limitations and constructive methods

We analyze the description of quantum many-body mixed states using matrix product states and operators. We consider two such descriptions: (i) as a matrix product density operator of bond dimension D, and (ii) as a purification that is written as a matrix product state of bond dimension D'. We show that these descriptions are inequivalent in the sense that D' cannot be upper bounded by D only. Then we provide two constructive methods to obtain (ii) out of (i). The sum of squares (sos) polynomial method scales exponentially in the number of different eigenvalues, and its approximate version is formulated as a Semidefinite Program, which gives efficient approximate purifications whose D' only depends on D. The eigenbasis method scales quadratically in the number of eigenvalues, and its approximate version is very efficient for rapidly decaying distributions of eigenvalues. Our results imply that a description of mixed states which is both efficient and locally positive semidefinite does not exist, but that good approximations do.

preprint2013arXiv

Quantum simulations of gauge theories with ultracold atoms: local gauge invariance from angular momentum conservation

Quantum simulations of High Energy Physics, and especially of gauge theories, is an emerging and exciting direction in quantum simulations. However, simulations of such theories, compared to simulations of condensed matter physics, must satisfy extra restrictions, such as local gauge and Lorentz invariance. In this paper we discuss these special requirements, and present a new method for quantum simulation of lattice gauge theories using ultracold atoms. This method allows to include local gauge invariance as a fundamental symmetry of the atomic Hamiltonian, arising from natural atomic interactions and conservation laws (and not as a property of a low energy sector). This allows us to implement elementary gauge invariant interactions for three lattice gauge theories: compact QED (U(1)), SU(N) and Z_N, which can be used to build quantum simulators in 1+1 dimensions. We also present a new loop method, which uses the elementary interactions as building blocks in the effective construction of quantum simulations for d+1 dimensional lattice gauge theories (d>1), without having to use Gauss's law as a constraint, as in previous proposals. We discuss in detail the quantum simulation of 2+1 dimensional compact QED and provide a numerical proof of principle. The simplicity of the already gauge invariant elementary interactions of this model suggests it may be useful for future experimental realizations.

preprint2013arXiv

Steady-State Entanglement in the Nuclear Spin Dynamics of a Double Quantum Dot

We propose a scheme for the deterministic generation of steady-state entanglement between the two nuclear spin ensembles in an electrically defined double quantum dot. Due to quantum interference in the collective coupling to the electronic degrees of freedom, the nuclear system is actively driven into a two-mode squeezed-like target state. The entanglement build-up is accompanied by a self-polarization of the nuclear spins towards large Overhauser field gradients. Moreover, the feedback between the electronic and nuclear dynamics leads to multi-stability and criticality in the steady-state solutions.

preprint2013arXiv

Topological order in PEPS: Transfer operator and boundary Hamiltonians

We study the structure of topological phases and their boundaries in the Projected Entangled Pair States (PEPS) formalism. We show how topological order in a system can be identified from the structure of the PEPS transfer operator, and subsequently use these findings to analyze the structure of the boundary Hamiltonian, acting on the bond variables, which reflects the entanglement properties of the system. We find that in a topological phase, the boundary Hamiltonian consists of two parts: A universal non-local part which encodes the nature of the topological phase, and a non-universal part which is local and inherits the symmetries of the topological model, which helps to infer the structure of the boundary Hamiltonian and thus possibly of the physical edge modes.

preprint2012arXiv

A cold-atom quantum simulator for SU(2) Yang-Mills lattice gauge theory

Non-abelian gauge theories play an important role in the standard model of particle physics, and unfold a partially unexplored world of exciting physical phenomena. In this letter, we suggest a realization of a non-abelian lattice gauge theory - SU(2) Yang-Mills in 1+1 dimensions, using ultracold atoms. Remarkably, and in contrast to previous proposals, in our model gauge invariance is a direct consequence of angular momentum conservation and thus is fundamental and robust. Our proposal may serve as well as a starting point for higher dimensional realizations.

preprint2012arXiv

Calculus of continuous matrix product states

We discuss various properties of the variational class of continuous matrix product states, a class of ansatz states for one-dimensional quantum fields that was recently introduced as the direct continuum limit of the highly successful class of matrix product states. We discuss both attributes of the physical states, e.g. by showing in detail how to compute expectation values, as well as properties intrinsic to the representation itself, such as the gauge freedom. We consider general translation non-invariant systems made of several particle species and derive certain regularity properties that need to be satisfied by the variational parameters. We also devote a section to the translation invariant setting in the thermodynamic limit and show how continuous matrix product states possess an intrinsic ultraviolet cutoff. Finally, we introduce a new set of states which are tangent to the original set of continuous matrix product states. For the case of matrix product states, this construction has recently proven relevant in the development of new algorithms for studying time evolution and elementary excitations of quantum spin chains. We thus lay the foundation for similar developments for one-dimensional quantum fields.

preprint2012arXiv

Dissipative Phase Transition in Central Spin Systems

We investigate dissipative phase transitions in an open central spin system. In our model the central spin interacts coherently with the surrounding many-particle spin environment and is subject to coherent driving and dissipation. We develop analytical tools based on a self-consistent Holstein-Primakoff approximation that enable us to determine the complete phase diagram associated with the steady states of this system. It includes first and second-order phase transitions, as well as regions of bistability, spin squeezing and altered spin pumping dynamics. Prospects of observing these phenomena in systems such as electron spins in quantum dots or NV centers coupled to lattice nuclear spins are briefly discussed.

preprint2012arXiv

Gaussian Matrix Product States

We introduce Gaussian Matrix Product States (GMPS), a generalization of Matrix Product States (MPS) to lattices of harmonic oscillators. Our definition resembles the interpretation of MPS in terms of projected maximally entangled pairs, starting from which we derive several properties of GMPS, often in close analogy to the finite dimensional case: We show how to approximate arbitrary Gaussian states by MPS, we discuss how the entanglement in the bonds can be bounded, we demonstrate how the correlation functions can be computed from the GMPS representation, and that they decay exponentially in one dimension, and finally relate GMPS and ground states of local Hamiltonians.

preprint2012arXiv

Laughlin spin liquid states on lattices obtained from conformal field theory

We propose a set of spin system wavefunctions that are very similar to lattice versions of the Laughlin states. The wavefunction are conformal blocks of conformal field theories, and for filling factor ν=1/2 we provide a parent Hamiltonian, which is valid for any even number of spins and is at the same time a 2D generalization of the Haldane-Shastry model. We also demonstrate that the Kalmeyer-Laughlin state is reproduced as a particular case of this model. Finally, we discuss various properties of the spin states and point out several analogies to known results for the Laughlin states.

preprint2012arXiv

Master equation approach to optomechanics with arbitrary dielectrics

We present a master equation describing the interaction of light with dielectric objects of arbitrary sizes and shapes. The quantum motion of the object, the quantum nature of light, as well as scattering processes to all orders in perturbation theory are taken into account. This formalism extends the standard master equation approach to the case where interactions among different modes of the environment are considered. It yields a genuine quantum description, including a renormalization of the couplings and decoherence terms. We apply this approach to analyze cavity cooling of the center-of-mass mode of large spheres. Furthermore, we derive an expression for the steady-state phonon numbers without relying on resolved-sideband or bad-cavity approximations.

preprint2012arXiv

Quantum simulation of small-polaron formation with trapped ions

We propose a quantum simulation of small-polaron physics using a one-dimensional system of trapped ions acted upon by off-resonant standing waves. This system, envisioned as an array of microtraps, in the single-excitation case allows the realization of the anti-adiabatic regime of the Holstein model. We show that the strong excitation-phonon coupling regime, characterized by the formation of small polarons, can be reached using realistic values of the relevant system parameters. Finally, we propose measurements of the quasiparticle residue and the average number of phonons in the ground state, experimental probes validating the polaronic character of the phonon-dressed excitation.

preprint2012arXiv

Quantum states on Harmonic lattices

We investigate bosonic Gaussian quantum states on an infinite cubic lattice in arbitrary spatial dimensions. We derive general properties of such states as ground states of quadratic Hamiltonians for both critical and non-critical cases. Tight analytic relations between the decay of the interaction and the correlation functions are proven and the dependence of the correlation length on band gap and effective mass is derived. We show that properties of critical ground states depend on the gap of the point-symmetrized rather than on that of the original Hamiltonian. For critical systems with polynomially decaying interactions logarithmic deviations from polynomially decaying correlation functions are found. Moreover, we provide a generalization of the matrix product state representation for Gaussian states and show that properties hold analogously to the case of finite dimensional spin systems.

preprint2012arXiv

Resonating valence bond states in the PEPS formalism

We study resonating valence bond (RVB) states in the Projected Entangled Pair States (PEPS) formalism. Based on symmetries in the PEPS description, we establish relations between the toric code state, the orthogonal dimer state, and the SU(2) singlet RVB state on the kagome lattice: We prove the equivalence of toric code and dimer state, and devise an interpolation between the dimer state and the RVB state. This interpolation corresponds to a continuous path in Hamiltonian space, proving that the RVB state is the four-fold degenerate ground state of a local Hamiltonian on the (finite) kagome lattice. We investigate this interpolation using numerical PEPS methods, studying the decay of correlation functions, the change of overlap, and the entanglement spectrum, none of which exhibits signs of a phase transition.

preprint2012arXiv

Robust entanglement generation by reservoir engineering

Following a recent proposal [C. Muschik et. al., Phys. Rev. A 83, 052312 (2011)], engineered dissipative processes have been used for the generation of stable entanglement between two macroscopic atomic ensembles at room temperature [H. Krauter et. al., Phys. Rev. Lett. 107, 080503 (2011)]. This experiment included the preparation of entangled states which are continuously available during a time interval of one hour. Here, we present additional material, further-reaching data and an extension of the theory developed in [C. Muschik et. al., Phys. Rev. A 83, 052312 (2011)]. In particular, we show how the combination of the entangling dissipative mechanism with measurements can give rise to a substantial improvement of the generated entanglement in the presence of noise.

preprint2012arXiv

Simulating Compact Quantum Electrodynamics with ultracold atoms: Probing confinement and nonperturbative effects

Recently, there has been much interest in simulating quantum field theory effects of matter and gauge fields. In a recent work [Phys. Rev. Lett. 107, 275301 (2011)] a method for simulating compact Quantum Electrodynamics (cQED) using Bose-Einstein condensates has been suggested. We suggest an alternative approach, which relies on single atoms in an optical lattice, carrying 2l+1 internal levels, which converges rapidly to cQED as l increases. That enables the simulation of cQED in 2+1 dimensions in both the weak and the strong coupling regimes, hence allowing to probe confinement as well as other nonperturbative effects of the theory. We provide an explicit construction for the case l=1 which is sufficient for simulating the effect of confinement between two external static charges.

preprint2012arXiv

Superradiance-like Electron Transport through a Quantum Dot

We theoretically show that intriguing features of coherent many-body physics can be observed in electron transport through a quantum dot (QD). We first derive a master equation based framework for electron transport in the Coulomb-blockade regime which includes hyperfine (HF) interaction with the nuclear spin ensemble in the QD. This general tool is then used to study the leakage current through a single QD in a transport setting. We find that, for an initially polarized nuclear system, the proposed setup leads to a strong current peak, in close analogy with superradiant emission of photons from atomic ensembles. This effect could be observed with realistic experimental parameters and would provide clear evidence of coherent HF dynamics of nuclear spin ensembles in QDs.

preprint2012arXiv

Tensor network techniques for the computation of dynamical observables in 1D quantum spin systems

We analyze the recently developed folding algorithm [Phys. Rev. Lett. 102, 240603 (2009)] to simulate the dynamics of infinite quantum spin chains, and relate its performance to the kind of entanglement produced under the evolution of product states. We benchmark the accomplishments of this technique with respect to alternative strategies using Ising Hamiltonians with transverse and parallel fields, as well as XY models. Additionally, we evaluate its ability to find ground and thermal equilibrium states.

preprint2012arXiv

Topological and Entanglement Properties of Resonating Valence Bond wavefunctions

We examine in details the connections between topological and entanglement properties of short-range resonating valence bond (RVB) wave functions using Projected Entangled Pair States (PEPS) on kagome and square lattices on (quasi-)infinite cylinders with generalized boundary conditions (and perimeters with up to 20 lattice spacings). Making use of disconnected topological sectors in the space of dimer lattice coverings, we explicitly derive (orthogonal) "minimally entangled" PEPS RVB states. For the kagome lattice, we obtain, using the quantum Heisenberg antiferromagnet as a reference model, the finite size scaling of the energy separations between these states. In particular, we extract two separate (vanishing) energy scales corresponding (i) to insert a vison line between the two ends of the cylinder and (ii) to pull out and freeze a spin at either end. We also investigate the relations between bulk and boundary properties and show that, for a bipartition of the cylinder, the boundary Hamiltonian defined on the edge can be written as a product of a highly non-local projector with an emergent (local) su(2)-invariant one-dimensional (superfluid) t--J Hamiltonian, which arises due to the symmetry properties of the auxiliary spins at the edge. This multiplicative structure, a consequence of the disconnected topological sectors in the space of dimer lattice coverings, is characteristic of the topological nature of the states. For minimally entangled RVB states, it is shown that the entanglement spectrum, which reflects the properties of the edge modes, is a subset (half for kagome RVB) of the spectrum of the local Hamiltonian, providing e.g. a simple argument on the origin of the topological entanglement entropy S0=-ln 2 of Z2 spin liquids. We propose to use these features to probe topological phases in microscopic Hamiltonians and some results are compared to existing DMRG data.

preprint2012arXiv

Variational matrix product ansatz for dispersion relations

A variational ansatz for momentum eigenstates of translation invariant quantum spin chains is formulated. The matrix product state ansatz works directly in the thermodynamic limit and allows for an efficient implementation (cubic scaling in the bond dimension) of the variational principle. Unlike previous approaches, the ansatz includes topologically non-trivial states (kinks, domain walls) for systems with symmetry breaking. The method is benchmarked using the spin-1/2 XXZ antiferromagnet and the spin-1 Heisenberg antiferromagnet and we obtain surprisingly accurate results.

preprint2011arXiv

Adiabatic Preparation of a Heisenberg Antiferromagnet Using an Optical Superlattice

We analyze the possibility to prepare a Heisenberg antiferromagnet with cold fermions in optical lattices, starting from a band insulator and adiabatically changing the lattice potential. The numerical simulation of the dynamics in 1D allows us to identify the conditions for success, and to study the influence that the presence of holes in the initial state may have on the protocol. We also extend our results to two-dimensional systems.

preprint2011arXiv

Applying the variational principle to (1+1)-dimensional quantum field theories

We extend the recently introduced continuous matrix product state (cMPS) variational class to the setting of (1+1)-dimensional relativistic quantum field theories. This allows one to overcome the difficulties highlighted by Feynman concerning the variational procedure applied to relativistic theories, and provides a new way to regularize quantum field theories. A fermionic version of the continuous matrix product state is introduced which is manifestly free from fermion doubling and sign problems. We illustrate the power of the formalism with the simulation of free massive Dirac fermions, the Gross-Neveu model, and the Casimir effect. We find that cMPS can capture chiral symmetry breaking with absolute scaling of the chiral parameter, and that boundary effects can be accommodated with modest computational effort.

preprint2011arXiv

Entanglement distillation by dissipation and continuous quantum repeaters

Even though entanglement is very vulnerable to interactions with the environment, it can be created by purely dissipative processes. Yet, the attainable degree of entanglement is profoundly limited in the presence of noise sources. We show that distillation can also be realized dissipatively, such that a highly entanglement steady state is obtained. The schemes put forward here display counterintuitive phenomena, such as improved performance if noise is added to the system. We also show how dissipative distillation can be employed in a continuous quantum repeater architecture, in which the resources scale polynomially with the distance.

preprint2011arXiv

Entanglement generated by dissipation and steady state entanglement of two macroscopic objects

Entanglement is a striking feature of quantum mechanics and an essential ingredient in most applications in quantum information. Typically, coupling of a system to an environment inhibits entanglement, particularly in macroscopic systems. Here we report on an experiment, where dissipation continuously generates entanglement between two macroscopic objects. This is achieved by engineering the dissipation using laser- and magnetic fields, and leads to robust event-ready entanglement maintained for 0.04s at room temperature. Our system consists of two ensembles containing about 10^{12} atoms and separated by 0.5m coupled to the environment composed of the vacuum modes of the electromagnetic field. By combining the dissipative mechanism with a continuous measurement, steady state entanglement is continuously generated and observed for up to an hour.

preprint2011arXiv

Entanglement spectrum and boundary theories with projected entangled-pair states

In many physical scenarios, close relations between the bulk properties of quantum systems and theories associated to their boundaries have been observed. In this work, we provide an exact duality mapping between the bulk of a quantum spin system and its boundary using Projected Entangled Pair States (PEPS). This duality associates to every region a Hamiltonian on its boundary, in such a way that the entanglement spectrum of the bulk corresponds to the excitation spectrum of the boundary Hamiltonian. We study various specific models, like a deformed AKLT [1], an Ising-type [2], and Kitaev's toric code [3], both in finite ladders and infinite square lattices. In the latter case, some of those models display quantum phase transitions. We find that a gapped bulk phase with local order corresponds to a boundary Hamiltonian with local interactions, whereas critical behavior in the bulk is reflected on a diverging interaction length of the boundary Hamiltonian. Furthermore, topologically ordered states yield non-local Hamiltonians. As our duality also associates a boundary operator to any operator in the bulk, it in fact provides a full holographic framework for the study of quantum many-body systems via their boundary.

preprint2011arXiv

Gapless Hamiltonians for the toric code using the PEPS formalism

We study Hamiltonians which have Kitaev's toric code as a ground state, and show how to construct a Hamiltonian which shares the ground space of the toric code, but which has gapless excitations with a continuous spectrum in the thermodynamic limit. Our construction is based on the framework of Projected Entangled Pair States (PEPS), and can be applied to a large class of two-dimensional systems to obtain gapless "uncle Hamiltonians".

preprint2011arXiv

Hawking Radiation on an Ion Ring in the Quantum Regime

This paper discusses a recent proposal for the simulation of acoustic black holes with ions. The ions are rotating on a ring with an inhomogeneous, but stationary velocity profile. Phonons cannot leave a region, in which the ion velocity exceeds the group velocity of the phonons, as light cannot escape from a black hole. The system is described by a discrete field theory with a nonlinear dispersion relation. Hawking radiation is emitted by this acoustic black hole, generating entanglement between the inside and the outside of the black hole. We study schemes to detect the Hawking effect in this setup.

preprint2011arXiv

Majorana Fermions in Equilibrium and Driven Cold Atom Quantum Wires

We introduce a new approach to create and detect Majorana fermions using optically trapped 1D fermionic atoms. In our proposed setup, two internal states of the atoms couple via an optical Raman transition---simultaneously inducing an effective spin-orbit interaction and magnetic field---while a background molecular BEC cloud generates s-wave pairing for the atoms. The resulting cold atom quantum wire supports Majorana fermions at phase boundaries between topologically trivial and nontrivial regions, as well as `Floquet Majorana fermions' when the system is periodically driven. We analyze experimental parameters, detection schemes, and various imperfections.

preprint2011arXiv

Quantum spin Hamiltonians for the SU(2)_k WZW model

We propose to use null vectors in conformal field theories to derive model Hamiltonians of quantum spin chains and corresponding ground state wave function(s). The approach is quite general, and we illustrate it by constructing a family of Hamiltonians whose ground states are the chiral correlators of the SU(2)_k WZW model for integer values of the level k. The simplest example corresponds to k=1 and is essentially a nonuniform generalization of the Haldane-Shastry model with long-range exchange couplings. At level k=2, we analyze the model for N spin 1 fields. We find that the Renyi entropy and the two-point spin correlator show, respectively, logarithmic growth and algebraic decay. Furthermore, we use the null vectors to derive a set of algebraic, linear equations relating spin correlators within each model. At level k=1, these equations allow us to compute the two-point spin correlators analytically for the finite chain uniform Haldane-Shastry model and to obtain numerical results for the nonuniform case and for higher-point spin correlators in a very simple way and without resorting to Monte Carlo techniques.

preprint2011arXiv

Time-dependent variational principle for quantum lattices

We develop a new algorithm based on the time-dependent variational principle applied to matrix product states to efficiently simulate the real- and imaginary time dynamics for infinite one-dimensional quantum lattice systems. This procedure: (1) is argued to be optimal; (2) does not rely on the Trotter decomposition and thus has no Trotter error; (3) explicitly preserves all symmetries and conservation laws; and (4) has low computational complexity. The algorithm is illustrated using both imaginary time and real-time examples.

preprint2011arXiv

Topologically Protected Quantum State Transfer in a Chiral Spin Liquid

Topology plays a central role in ensuring the robustness of a wide variety of physical phenomena. Notable examples range from the robust current carrying edge states associated with the quantum Hall and the quantum spin Hall effects to proposals involving topologically protected quantum memory and quantum logic operations. Here, we propose and analyze a topologically protected channel for the transfer of quantum states between remote quantum nodes. In our approach, state transfer is mediated by the edge mode of a chiral spin liquid. We demonstrate that the proposed method is intrinsically robust to realistic imperfections associated with disorder and decoherence. Possible experimental implementations and applications to the detection and characterization of spin liquid phases are discussed.

preprint2011arXiv

Unforgeable Noise-Tolerant Quantum Tokens

The realization of devices which harness the laws of quantum mechanics represents an exciting challenge at the interface of modern technology and fundamental science. An exemplary paragon of the power of such quantum primitives is the concept of "quantum money". A dishonest holder of a quantum bank-note will invariably fail in any forging attempts; indeed, under assumptions of ideal measurements and decoherence-free memories such security is guaranteed by the no-cloning theorem. In any practical situation, however, noise, decoherence and operational imperfections abound. Thus, the development of secure "quantum money"-type primitives capable of tolerating realistic infidelities is of both practical and fundamental importance. Here, we propose a novel class of such protocols and demonstrate their tolerance to noise; moreover, we prove their rigorous security by determining tight fidelity thresholds. Our proposed protocols require only the ability to prepare, store and measure single qubit quantum memories, making their experimental realization accessible with current technologies.

preprint2011arXiv

Violation of the area law and long range correlations in infinite matrix product states

We propose to construct a family of states of one-dimensional spin chains by replacing the finite dimensional matrices in matrix product states by infinite dimensional operators constructed from bosonic annihilation and creation operators. The resulting states are demonstrated to violate the area law and to exhibit long range correlations. In addition, we propose an efficient way to prepare the states experimentally, in which the spins interact sequentially with an ancilla system.

preprint2010arXiv

A quantum interface between light and nuclear spins in quantum dots

The coherent coupling of flying photonic qubits to stationary matter-based qubits is an essential building block for quantum communication networks. We show how such a quantum interface can be realized between a traveling-wave optical field and the polarized nuclear spins in a singly charged quantum dot strongly coupled to a high-finesse optical cavity. By adiabatically eliminating the electron a direct effective coupling is achieved. Depending on the laser field applied, interactions that enable either write-in or read-out are obtained.

preprint2010arXiv

Cold atom simulation of interacting relativistic quantum field theories

We demonstrate that Dirac fermions self-interacting or coupled to dynamic scalar fields can emerge in the low energy sector of designed bosonic and fermionic cold atom systems. We illustrate this with two examples defined in two spacetime dimensions. The first one is the self-interacting Thirring model. The second one is a model of Dirac fermions coupled to a dynamic scalar field that gives rise to the Gross-Neveu model. The proposed cold atom experiments can be used to probe spectral or correlation properties of interacting quantum field theories thereby presenting an alternative to lattice gauge theory simulations.

preprint2010arXiv

Complete devil's staircase and crystal--superfluid transitions in a dipolar XXZ spin chain: A trapped ion quantum simulation

Systems with long-range interactions show a variety of intriguing properties: they typically accommodate many meta-stable states, they can give rise to spontaneous formation of supersolids, and they can lead to counterintuitive thermodynamic behavior. However, the increased complexity that comes with long-range interactions strongly hinders theoretical studies. This makes a quantum simulator for long-range models highly desirable. Here, we show that a chain of trapped ions can be used to quantum simulate a one-dimensional model of hard-core bosons with dipolar off-site interaction and tunneling, equivalent to a dipolar XXZ spin-1/2 chain. We explore the rich phase diagram of this model in detail, employing perturbative mean-field theory, exact diagonalization, and quasiexact numerical techniques (density-matrix renormalization group and infinite time evolving block decimation). We find that the complete devil's staircase -- an infinite sequence of crystal states existing at vanishing tunneling -- spreads to a succession of lobes similar to the Mott-lobes found in Bose--Hubbard models. Investigating the melting of these crystal states at increased tunneling, we do not find (contrary to similar two-dimensional models) clear indications of supersolid behavior in the region around the melting transition. However, we find that inside the insulating lobes there are quasi-long range (algebraic) correlations, opposed to models with nearest-neighbor tunneling which show exponential decay of correlations.

preprint2010arXiv

Dissipatively driven entanglement of two macroscopic atomic ensembles

Up to date, the life time of experimentally demonstrated entangled states has been limited, due to their fragility under decoherence and dissipation. Therefore, they are created under strict isolation conditions. In contrast, new approaches harness the coupling of the system to the environment, which drives the system into the desired state. Following these ideas, we present a robust method for generating steady state entanglement between two distant atomic ensembles. The proposed scheme relies on the interaction of the two atomic systems with the common vacuum modes of the electromagnetic field which act as an engineered environment. We develop the theoretical framework for two level systems including dipole-dipole interactions and complement these results by considering the implementation in multi-level ground states.

preprint2010arXiv

Emerging Bosons with Three-Body Interactions from Spin-1 Atoms in Optical Lattices

We study two many-body systems of bosons interacting via an infinite three-body contact repulsion in a lattice: a pairs quasi-condensate induced by correlated hopping and the discrete version of the Pfaffian wavefunction. We propose to experimentally realise systems characterized by such interaction by means of a proper spin-1 lattice Hamiltonian: spin degrees of freedom are locally mapped into occupation numbers of emerging bosons, in a fashion similar to spin-1/2 and hardcore bosons. Such a system can be realized with ultracold spin-1 atoms in a Mott Insulator with filling-factor one. The high versatility of these setups allows us to engineer spin-hopping operators breaking the SU(2) symmetry, as needed to approximate interesting bosonic Hamiltonians with three-body hardcore constraint. For this purpose we combine bichromatic spin-independent superlattices and Raman transitions to induce a different hopping rate for each spin orientation. Finally, we illustrate how our setup could be used to experimentally realize the first setup, i.e. the transition to a pairs quasi-condensed phase of the emerging bosons. We also report on a route towards the realization of a discrete bosonic Pfaffian wavefunction and list some open problems to reach this goal.

preprint2010arXiv

Fermionic Projected Entangled Pair States

We introduce a family of states, the fPEPS, which describes fermionic systems on lattices in arbitrary spatial dimensions. It constitutes the natural extension of another family of states, the PEPS, which efficiently approximate ground and thermal states of spin systems with short-range interactions. We give an explicit mapping between those families, which allows us to extend previous simulation methods to fermionic systems. We also show that fPEPS naturally arise as exact ground states of certain fermionic Hamiltonians. We give an example of such a Hamiltonian, exhibiting criticality while obeying an area law.

preprint2010arXiv

Generalized Hartree-Fock Theory for Interacting Fermions in Lattices: Numerical Methods

We present numerical methods to solve the Generalized Hartree-Fock theory for fermionic systems in lattices, both in thermal equilibrium and out of equilibrium. Specifically, we show how to determine the covariance matrix corresponding to the Fermionic Gaussian state that optimally approximates the quantum state of the fermions. The methods apply to relatively large systems, since their complexity only scales quadratically with the number of lattice sites. Moreover, they are specially suited to describe inhomogenous systems, as those typically found in recent experiments with atoms in optical lattices, at least in the weak interaction regime. As a benchmark, we have applied them to the two-dimensional Hubbard model on a 10x10 lattice with and without an external confinement.

preprint2010arXiv

Ground-state properties of the spin-1/2 antiferromagnetic Heisenberg model on the triangular lattice: A variational study based on entangled-plaquette states

We study, on the basis of the general entangled-plaquette variational ansatz, the ground-state properties of the spin-1/2 antiferromagnetic Heisenberg model on the triangular lattice. Our numerical estimates are in good agreement with available exact results and comparable, for large system sizes, to those computed via the best alternative numerical approaches, or by means of variational schemes based on specific (i.e., incorporating problem dependent terms) trial wave functions. The extrapolation to the thermodynamic limit of our results for lattices comprising up to N=324 spins yields an upper bound of the ground-state energy per site (in units of the exchange coupling) of $-0.5458(2)$ [$-0.4074(1)$ for the XX model], while the estimated infinite-lattice order parameter is $0.3178(5)$ (i.e., approximately 64% of the classical value).

preprint2010arXiv

Hawking Radiation from an Acoustic Black Hole on an Ion Ring

In this article we propose to simulate acoustic black holes with ions in rings. If the ions are rotating with a stationary and inhomogeneous velocity profile, regions can appear where the ion velocity exceeds the group velocity of the phonons. In these regions phonons are trapped like light in black holes, even though we have a discrete field theory and a nonlinear dispersion relation. We study the appearance of Hawking radiation in this setup and propose a scheme to detect it.

preprint2010arXiv

Matrix Product State and mean field solutions for one-dimensional systems can be found efficiently

We consider the problem of approximating ground states of one-dimensional quantum systems within the two most common variational ansatzes, namely the mean field ansatz and Matrix Product States. We show that both for mean field and for Matrix Product States of fixed bond dimension, the optimal solutions can be found in a way which is provably efficient (i.e., scales polynomially). This implies that the corresponding variational methods can be in principle recast in a way which scales provably polynomially. Moreover, our findings imply that ground states of one-dimensional commuting Hamiltonians can be found efficiently.

preprint2010arXiv

Modified spin-wave theory with ordering vector optimization II: Spatially anisotropic triangular lattice and $J_1J_2J_3$ model with Heisenberg interactions

We study the ground state phases of the $S=1/2$ Heisenberg quantum antiferromagnet on the spatially anisotropic triangular lattice and on the square lattice with up to next-next-nearest neighbor coupling (the $J_1J_2J_3$ model), making use of Takahashi's modified spin-wave (MSW) theory supplemented by ordering vector optimization. We compare the MSW results with exact diagonalization and projected-entangled-pair-states calculations, demonstrating their qualitative and quantitative reliability. We find that MSW theory correctly accounts for strong quantum effects on the ordering vector of the magnetic phases of the models under investigation: in particular collinear magnetic order is promoted at the expenses of non-collinear (spiral) order, and several spiral states which are stable at the classical level, disappear from the quantum phase diagram. Moreover, collinear states and non-collinear ones are never connected continuously, but they are separated by parameter regions in which MSW breaks down, signaling the possible appearance of a non-magnetic ground state. In the case of the spatially anisotropic triangular lattice, a large breakdown region appears also for weak couplings between the chains composing the lattice, suggesting the possible occurrence of a large non-magnetic region continuously connected with the spin-liquid state of the uncoupled chains.

preprint2010arXiv

Optically Levitating Dielectrics in the Quantum Regime: Theory and Protocols

We provide a general quantum theory to describe the coupling of light with the motion of a dielectric object inside a high finesse optical cavity. In particular, we derive the total Hamiltonian of the system as well as a master equation describing the state of the center of mass mode of the dielectric and the cavity field mode. In addition, a quantum theory of elasticity is used in order to study the coupling of the center of mass motion with internal vibrational excitations of the dielectric. This general theory is applied to the recent proposal of using an optically levitating nanodielectric as a cavity optomechanical system [Romero-Isart et al. NJP 12, 033015 (2010), Chang et al. PNAS 107, 1005 (2010)]. On this basis, we also design a light-mechanics interface to prepare non-Gaussian states of the mechanical motion, such as quantum superpositions of Fock states. Finally, we introduce a direct mechanical tomography scheme to probe these genuine quantum states by time of flight experiments.

preprint2010arXiv

Scalable Architecture for a Room Temperature Solid-State Quantum Information Processor

The realization of a scalable quantum information processor has emerged over the past decade as one of the central challenges at the interface of fundamental science and engineering. Much progress has been made towards this goal. Indeed, quantum operations have been demonstrated on several trapped ion qubits, and other solid-state systems are approaching similar levels of control. Extending these techniques to achieve fault-tolerant operations in larger systems with more qubits remains an extremely challenging goal, in part, due to the substantial technical complexity of current implementations. Here, we propose and analyze an architecture for a scalable, solid-state quantum information processor capable of operating at or near room temperature. The architecture is applicable to realistic conditions, which include disorder and relevant decoherence mechanisms, and includes a hierarchy of control at successive length scales. Our approach is based upon recent experimental advances involving Nitrogen-Vacancy color centers in diamond and will provide fundamental insights into the physics of non-equilibrium many-body quantum systems. Additionally, the proposed architecture may greatly alleviate the stringent constraints, currently limiting the realization of scalable quantum processors.

preprint2010arXiv

Simulating quantum-optical phenomena with cold atoms in optical lattices

We propose a scheme involving cold atoms trapped in optical lattices to observe different phenomena traditionally linked to quantum-optical systems. The basic idea consists of connecting the trapped atomic state to a non-trapped state through a Raman scheme. The coupling between these two types of atoms (trapped and free) turns out to be similar to that describing light-matter interaction within the rotating-wave approximation, the role of matter and photons being played by the trapped and free atoms, respectively. We explain in particular how to observe phenomena arising from the collective spontaneous emission of atomic and harmonic oscillator samples such as superradiance and directional emission. We also show how the same setup can simulate Bose-Hubbard Hamiltonians with extended hopping as well as Ising models with long-range interactions. We believe that this system can be realized with state of the art technology.

preprint2010arXiv

Simulating two- and three-dimensional frustrated quantum systems with string-bond states

Simulating frustrated quantum magnets is among the most challenging tasks in computational physics. We apply String-Bond States, a recently introduced ansatz which combines Tensor Networks with Monte Carlo based methods, to the simulation of frustrated quantum systems in both two and three dimensions. We compare our results with existing results for unfrustrated and two-dimensional systems with open boundary conditions, and demonstrate that the method applies equally well to the simulation of frustrated systems with periodic boundaries in both two and three dimensions.

preprint2010arXiv

Strong and weak thermalization of infinite non-integrable quantum systems

When a non-integrable system evolves out of equilibrium for a long time, local observables are expected to attain stationary expectation values, independent of the details of the initial state. However, intriguing experimental results with ultracold gases have shown no thermalization in non-integrable settings, triggering an intense theoretical effort to decide the question. Here we show that the phenomenology of thermalization in a quantum system is much richer than its classical counterpart. Using a new numerical technique, we identify two distinct thermalization regimes, strong and weak, occurring for different initial states. Strong thermalization, intrinsically quantum, happens when instantaneous local expectation values converge to the thermal ones. Weak thermalization, well-known in classical systems, happens when local expectation values converge to the thermal ones only after time averaging. Remarkably, we find a third group of states showing no thermalization, neither strong nor weak, to the time scales one can reliably simulate.

preprint2010arXiv

Toward Quantum Superposition of Living Organisms

The most striking feature of quantum mechanics is the existence of superposition states, where an object appears to be in different situations at the same time. The existence of such states has been tested with small objects, like atoms, ions, electrons and photons, and even with molecules. More recently, it has been possible to create superpositions of collections of photons, atoms, or Cooper pairs. Current progress in optomechanical systems may soon allow us to create superpositions of even larger objects, like micro-sized mirrors or cantilevers, and thus to test quantum mechanical phenomena at larger scales. Here we propose a method to cool down and create quantum superpositions of the motion of sub-wavelength, arbitrarily shaped dielectric objects trapped inside a high--finesse cavity at a very low pressure. Our method is ideally suited for the smallest living organisms, such as viruses, which survive under low vacuum pressures, and optically behave as dielectric objects. This opens up the possibility of testing the quantum nature of living organisms by creating quantum superposition states in very much the same spirit as the original Schrödinger's cat "gedanken" paradigm. We anticipate our essay to be a starting point to experimentally address fundamental questions, such as the role of life and consciousness in quantum mechanics.

preprint2009arXiv

Ground-State Properties of Quantum Many-Body Systems: Entangled-Plaquette States and Variational Monte Carlo

We propose a new ansatz for the ground-state wave function of quantum many-body systems on a lattice. The key idea is to cover the lattice with plaquettes and obtain a state whose configurational weights can be optimized by means of a Variational Monte Carlo algorithm. Such a scheme applies to any dimension, without any "sign" instability. We show results for various two dimensional spin models (including frustrated ones). A detailed comparison with available exact results, as well as with variational methods based on different ansatzs is offered. In particular, our numerical estimates are in quite good agreement with exact ones for unfrustrated systems, and compare favorably to other methods for frustrated ones.

preprint2009arXiv

Infinite matrix product states, Conformal Field Theory and the Haldane-Shastry model

We generalize the Matrix Product States method using the chiral vertex operators of Conformal Field Theory and apply it to study the ground states of the XXZ spin chain, the J1-J2 model and random Heisenberg models. We compute the overlap with the exact wave functions, spin-spin correlators and the Renyi entropy, showing that critical systems can be described by this method. For rotational invariant ansatzs we construct an inhomogenous extension of the Haldane-Shastry model with long range exchange interactions.

preprint2009arXiv

Interfacing nuclear spins in quantum dots to cavity or traveling-wave fields

We show how to realize a quantum interface between optical fields and the polarized nuclear spins in a singly charged quantum dot, which is strongly coupled to a high-finesse optical cavity. An effective direct coupling between cavity and nuclear spins is obtained by adiabatically eliminating the (far detuned) excitonic and electronic states. The requirements needed to map qubit and continuous variable states of cavity or traveling-wave fields to the collective nuclear spin are investigated: For cavity fields, we consider adiabatic passage processes to transfer the states. It is seen that a significant improvement in cavity lifetimes beyond present-day technology would be required for a quantum interface. We then turn to a scheme which couples the nuclei to the output field of the cavity and can tolerate significantly shorter cavity lifetimes. We show that the lifetimes reported in the literature and the recently achieved nuclear polarization of ~90% allow both high-fidelity read-out and write-in of quantum information between the nuclear spins and the output field. We discuss the performance of the scheme and provide a convenient description of the dipolar dynamics of the nuclei for highly polarized spins, demonstrating that this process does not affect the performance of our protocol.

preprint2009arXiv

Modified spin-wave theory with ordering vector optimization I: frustrated bosons on the spatially anisotropic triangular lattice

We investigate a system of frustrated hardcore bosons, modeled by an XY antiferromagnet on the spatially anisotropic triangular lattice, using Takahashi's modified spin-wave (MSW) theory. In particular we implement ordering vector optimization on the ordered reference state of MSW theory, which leads to significant improvement of the theory and accounts for quantum corrections to the classically ordered state. The MSW results at zero temperature compare favorably to exact diagonalization (ED) and projected entangled-pair state (PEPS) calculations. The resulting zero-temperature phase diagram includes a 1D quasi-ordered phase, a 2D Neel ordered phase, and a 2D spiraling ordered phase. We have strong indications that the various ordered or quasi-ordered phases are separated by spin-liquid phases with short-range correlations, in analogy to what has been predicted for the Heisenberg model on the same lattice. Within MSW theory we also explore the finite-temperature phase diagram. We find that the zero-temperature long-range-ordered phases turn into quasi-ordered phases (up to a Berezinskii-Kosterlitz-Thouless temperature), while zero-temperature quasi-ordered phases become short-range correlated at finite temperature. These results show that modified spin-wave theory is very well suited for describing ordered and quasi-ordered phases of frustrated XY spins (or, equivalently, of frustrated lattice bosons) both at zero and finite temperatures. While MSW theory, just as other theoretical methods, cannot describe spin-liquid phases, its breakdown provides a fast method for singling out Hamiltonians which may feature these intriguing quantum phases. We thus suggest a tool for guiding our search for interesting systems whose properties are necessarily studied with a physical quantum simulator.

preprint2008arXiv

Dividing Quantum Channels

We investigate the possibility of dividing quantum channels into concatenations of other channels, thereby studying the semigroup structure of the set of completely-positive trace-preserving maps. We show the existence of 'indivisible' channels which can not be written as non-trivial products of other channels and study the set of 'infinitesimal divisible' channels which are elements of continuous completely positive evolutions. For qubit channels we obtain a complete characterization of the sets of indivisible and infinitesimal divisible channels. Moreover, we identify those channels which are solutions of time-dependent master equations for both positive and completely positive evolutions. For arbitrary finite dimension we prove a representation theorem for elements of continuous completely positive evolutions based on new results on determinants of quantum channels and Markovian approximations.

preprint2007arXiv

Quantum processing photonic states in optical lattices

The mapping of photonic states to collective excitations of atomic ensembles is a powerful tool which finds a useful application in the realization of quantum memories and quantum repeaters. In this work we show that cold atoms in optical lattices can be used to perform an entangling unitary operation on the transferred atomic excitations. After the release of the quantum atomic state, our protocol results in a deterministic two qubit gate for photons. The proposed scheme is feasible with current experimental techniques and robust against the dominant sources of noise.

preprint2007arXiv

The computational complexity of PEPS

We determine the computational power of preparing Projected Entangled Pair States (PEPS), as well as the complexity of classically simulating them, and generally the complexity of contracting tensor networks. While creating PEPS allows to solve PP problems, the latter two tasks are both proven to be #P-complete. We further show how PEPS can be used to approximate ground states of gapped Hamiltonians, and that creating them is easier than creating arbitrary PEPS. The main tool for our proofs is a duality between PEPS and postselection which allows to use existing results from quantum compexity.

preprint2006arXiv

Entanglement Percolation in Quantum Networks

Quantum networks are composed of nodes which can send and receive quantum states by exchanging photons. Their goal is to facilitate quantum communication between any nodes, something which can be used to send secret messages in a secure way, and to communicate more efficiently than in classical networks. These goals can be achieved, for instance, via teleportation. Here we show that the design of efficient quantum communication protocols in quantum networks involves intriguing quantum phenomena, depending both on the way the nodes are displayed, and the entanglement between them. These phenomena can be employed to design protocols which overcome the exponential decrease of signals with the number of nodes. We relate the problem of establishing maximally entangled states between nodes to classical percolation in statistical mechanics, and demonstrate that quantum phase transitions can be used to optimize the operation of quantum networks.

preprint2006arXiv

Quantum emulsion: a glassy phase of bosonic mixtures in optical lattices

We numerically investigate mixtures of two interacting bosonic species with unequal parameters in one-dimensional optical lattices. In large parameter regions full phase segregation is seen to minimize the energy of the system, but the true ground state is masked by an exponentially large number of metastable states characterized by microscopic phase separation. The ensemble of these quantum emulsion states, reminiscent of emulsions of immiscible fluids, has macroscopic properties analogous to those of a Bose glass, namely a finite compressibility in absence of superfluidity. Their metastability is probed by extensive quantum Monte Carlo simulations generating a rich correlated stochastic dynamics. The tuning of the repulsion of one of the two species via a Feshbach resonance drives the system through a quantum phase transition to the superfluid state.

preprint2002arXiv

On the effective size of certain "Schrödinger cat'' like states

Several experiments and experimental proposals for the production of macroscopic superpositions naturally lead to states of the general form $|ϕ_1>^{\otimes N}+|ϕ_2>^{\otimes N}$, where the number of subsystems $N$ is very large, but the states of the individual subsystems have large overlap, $|łϕ_1|ϕ_2 \r|^2=1-ε^2$. We propose two different methods for assigning an effective particle number to such states, using ideal Greenberger--Horne--Zeilinger (GHZ)-- states of the form $|0\r^{\otimes n}+|1\r^{\otimes n}$ as a standard of comparison. The two methods are based on decoherence and on a distillation protocol respectively. Both lead to an effective size $n$ of the order of $N ε^2$.