Search

Coherent driving has established itself as a powerful tool for guiding a many-body quantum system into a desirable, coherent non-equilibrium state. A thermodynamically large system will, however, almost always saturate to a featureless infinite temperature state under continuous driving and so the optical manipulation of many-body systems is considered feasible only if a transient, prethermal regime exists, where heating is suppressed. Here we show that, counterintuitively, in a broad class of lattices Floquet heating can actually be an advantageous effect. Specifically, we prove that the maximum entropy steady states which form upon driving the ground state of the Hubbard model on unbalanced bi-partite lattices possess uniform off-diagonal long-range order which remains finite even in the thermodynamic limit. This creation of a `hot' condensate can occur on \textit{any} driven unbalanced lattice and provides an understanding of how heating can, at the macroscopic level, expose and alter the order in a quantum system. We discuss implications for recent experiments observing emergent superconductivity in photoexcited materials., Comment: 10 pages, 3 figures

In triangular lattice structures, spatial anisotropy and frustration can lead to rich equilibrium phase diagrams with regions containing complex, highly entangled states of matter. In this work we study the driven two-rung triangular Hubbard model and evolve these states out of equilibrium, observing how the interplay between the driving and the initial state unexpectedly shuts down the particle-hole excitation pathway. This restriction, which symmetry arguments fail to predict, dictates the transient dynamics of the system, causing the available particle-hole degrees of freedom to manifest uniform long-range order. We discuss implications of our results for a recent experiment on photo-induced superconductivity in ${\rm \kappa - (BEDT-TTF)_{2}Cu[N(CN)_{2}]Br}$ molecules., Comment: Main Text: 7 Pages, 4 Figures, Supplementary: 4 Pages, 3 Figures

Far and mid infrared optical pulses have been shown to induce non-equilibrium unconventional orders in complex materials, including photo-induced ferroelectricity in quantum paraelectrics, magnetic polarization in antiferromagnets and transient superconducting correlations in the normal state of cuprates and organic conductors. In the case of non-equilibrium superconductivity, femtosecond drives have generally resulted in electronic properties that disappear immediately after excitation, evidencing a state that lacks intrinsic rigidity. Here, we make use of a new optical device to drive metallic K$_3$C$_{60}$ with mid-infrared pulses of tunable duration, ranging between one picosecond and one nanosecond. The same superconducting-like optical properties observed over short time windows for femtosecond excitation are shown here to become metastable under sustained optical driving, with lifetimes in excess of ten nanoseconds. Direct electrical probing becomes possible at these timescales, yielding a vanishingly small resistance. Such a colossal positive photo-conductivity is highly unusual for a metal and, when taken together with the transient optical conductivities, it is rather suggestive of metastable light-induced superconductivity., Comment: 14 pages, 5 figures, supplementary materials

Superconductivity in organic conductors is often tuned by the application of chemical or external pressure. With this type of tuning, orbital overlaps and electronic bandwidths are manipulated, whilst the properties of the molecular building blocks remain virtually unperturbed.Here, we show that the excitation of local molecular vibrations in the charge-transfer salt $\kappa-(BEDT-TTF)_2Cu[N(CN)_2]Br$ induces a colossal increase in carrier mobility and the opening of a superconducting-like optical gap. Both features track the density of quasi-particles of the equilibrium metal, and can be achieved up to a characteristic coherence temperature $T^* \approxeq 50 K$, far higher than the equilibrium transition temperature $T_C = 12.5 K$. Notably, the large optical gap achieved by photo-excitation is not observed in the equilibrium superconductor, pointing to a light induced state that is different from that obtained by cooling. First-principle calculations and model Hamiltonian dynamics predict a transient state with long-range pairing correlations, providing a possible physical scenario for photo-molecular superconductivity., Comment: 11 pages, 7 figures, supplementary material

The out-of-equilibrium dynamics of quantum systems is one of the most fascinating problems in physics, with outstanding open questions on issues such as relaxation to equilibrium. An area of particular interest concerns few-body systems, where quantum and thermal fluctuations are expected to be especially relevant. In this contribution, we present numerical results demonstrating the impact of conserved quantities (or 'charges') in the outcomes of out-of-equilibrium measurements starting from realistic equilibrium states on a few-body system implementing the Dicke model., Comment: 12 pages, 1 fig. Contribution to Proceedings of the 24th European Conference on Few-Body Problems in Physics (EFB24). Matches journal version published under CC BY 4.0

The $t$-$J$ model is a standard model of strongly correlated electrons, often studied in the context of high-$T_c$ superconductivity. However, most studies of this model neglect three-site terms, which appear at the same order as the superexchange $J$. As these terms correspond to pair-hopping, they are expected to play an important role in the physics of superconductivity when doped sufficiently far from half-filling. In this paper we present a density matrix renormalisation group study of the one-dimensional $t$-$J$ model with the pair hopping terms included. We demonstrate that that these additional terms radically change the one-dimensional ground state phase diagram, extending the superconducting region at low fillings, while at larger fillings, superconductivity is completely suppressed. We explain this effect by introducing a simplified effective model of repulsive hardcore bosons., Comment: 10 Pages, 7 figures

The non-equilibrium dynamics of quantum many-body systems is one of the most fascinating problems in physics. Open questions range from how they relax to equilibrium to how to extract useful work from them. A critical point lies in assessing whether a system has conserved quantities (or 'charges'), as these can drastically influence its dynamics. Here we propose a general protocol to reveal the existence of charges based on a set of exact relations between out-of-equilibrium fluctuations and equilibrium properties of a quantum system. We apply these generalised quantum fluctuation relations to a driven quantum simulator, demonstrating their relevance to obtain unbiased temperature estimates from non-equilibrium measurements. Our findings will help guide research on the interplay of quantum and thermal fluctuations in quantum simulation, in studying the transition from integrability to chaos and in the design of new quantum devices., Comment: v3: substantially improved presentation and discussion, matches published version; accompanied by numerical data in figures

Recent experiments performed on cuprates and alkali-doped fullerides have demonstated that key signatures of superconductivity can be induced above the equilibrium critical temperature by optical modulation. These observations in disparate physical systems may indicate a general underlying mechanism. Multiple theories have been proposed, but these either consider specific features, such as competing instabilities, or focus on conventional BCS-type superconductivity. Here we show that periodic driving can enhance electron pairing in strongly-correlated systems. Focusing on the strongly-repulsive limit of the doped Hubbard model, we investigate in-gap, spatially inhomogeneous, on-site modulations. We demonstrate that such modulations substantially reduce electronic hopping, while simultaneously sustaining super-exchange interactions and pair hopping via driving-induced virtual charge excitations. We calculate real-time dynamics for the one-dimensional case, starting from zero and finite temperature initial states, and show that enhanced singlet--pair correlations emerge quickly and robustly in the out-of-equilibrium many-body state. Our results reveal a fundamental pairing mechanism that might underpin optically induced superconductivity in some strongly correlated quantum materials., Comment: 14 pages, 11 figures

We investigate cold bosonic impurity atoms trapped in a vortex lattice formed by condensed bosons of another species. We describe the dynamics of the impurities by a bosonic Hubbard model containing occupation-dependent parameters to capture the effects of strong impurity-impurity interactions. These include both a repulsive direct interaction and an attractive effective interaction mediated by the BEC. The occupation dependence of these two competing interactions drastically affects the Hubbard model phase diagram, including causing the disappearance of some Mott lobes, Comment: 5 pages main text, 18 pages supplemental material

In the present work we investigate the existence of multiple nonequilibrium steady states in a coherently driven XY lattice of dissipative two-level systems. A commonly used mean-field ansatz, in which spatial correlations are neglected, predicts a bistable behavior with a sharp shift between low- and high-density states. In contrast one-dimensional matrix product methods reveal these effects to be artifacts of the mean-field approach, with both disappearing once correlations are taken fully into account. Instead, a bunching-antibunching transition emerges. This indicates that alternative approaches should be considered for higher spatial dimensions, where classical simulations are currently infeasible. Thus we propose a circuit QED quantum simulator implementable with current technology to enable an experimental investigation of the model considered.

We propose a non-linear, hybrid quantum-classical scheme for simulating non-equilibrium dynamics of strongly correlated fermions described by the Hubbard model in a Bethe lattice in the thermodynamic limit. Our scheme implements non-equilibrium dynamical mean field theory (DMFT) and uses a digital quantum simulator to solve a quantum impurity problem whose parameters are iterated to self-consistency via a classically computed feedback loop where quantum gate errors can be partly accounted for. We analyse the performance of the scheme in an example case., Comment: 15 pages, 5 figures

Estimating the temperature of a cold quantum system is difficult. Usually, one measures a well-understood thermal state and uses that prior knowledge to infer its temperature. In contrast, we introduce a method of thermometry that assumes minimal knowledge of the state of a system and is potentially non-destructive. Our method uses a universal temperature-dependence of the quench dynamics of an initially thermal system coupled to a qubit probe that follows from the Tasaki-Crooks theorem for non-equilibrium work distributions. We provide examples for a cold-atom system, in which our thermometry protocol may retain accuracy and precision at subnanokelvin temperatures., Comment: Updated to published version. 6 pages plus 11 pages of supplemental material, and some numerical data

The control of non-equilibrium phenomena in complex solids is an important research frontier, encompassing new effects like light induced superconductivity. Here, we show that coherent optical excitation of molecular vibrations in the organic conductor K3C60 can induce a non-equilibrium state with the optical properties of a superconductor. A transient gap in the real part of the optical conductivity and a low-frequency divergence of the imaginary part are measured for base temperatures far above equilibrium Tc=20 K. These findings underscore the role of coherent light fields in inducing emergent order., Comment: 40 pages, 23 figures

Static strain in complex oxide heterostructures has been extensively used to engineer electronic and magnetic properties at equilibrium. In the same spirit, deformations of the crystal lattice with light may be used to achieve functional control across hetero-interfaces dynamically. Here, by exciting large amplitude infrared-active vibrations in a LaAlO3 substrate we induce magnetic order melting in a NdNiO3 film across a hetero-interface. Femtosecond Resonant Soft X-ray Diffraction is used to determine the spatial and temporal evolution of the magnetic disordering. We observe a magnetic melt front that grows from the substrate interface into the film, at a speed that suggests electronically driven propagation. Light control and ultrafast phase front propagation at hetero-interfaces may lead to new opportunities in optomagnetism, for example by driving domain wall motion to transport information across suitably designed devices., Comment: 35 pages, 8 Figures (both incl. Supplement)

We propose and analyze a scheme for parametrically cooling bilayer cuprates based on the selective driving of a $c$-axis vibrational mode. The scheme exploits the vibration as a transducer making the Josephson plasma frequencies time-dependent. We show how modulation at the difference frequency between the intra- and interbilayer plasmon substantially suppresses interbilayer phase fluctuations, responsible for switching $c$-axis transport from a superconducting to resistive state. Our calculations indicate that this may provide a viable mechanism for stabilizing non-equilibrium superconductivity even above $T_c$, provided a finite pair density survives between the bilayers out of equilibrium., Comment: 4 pages + 7 page supplemental

We introduce a detector that selectively probes the phononic excitations of a cold Bose gas. The detector is composed of a single impurity atom confined by a double-well potential, where the two lowest eigenstates of the impurity form an effective probe qubit that is coupled to the phonons via density-density interactions with the bosons. The system is analogous to a two-level atom coupled to photons of the radiation field. We demonstrate that tracking the evolution of the qubit populations allows probing both thermal and coherent excitations in targeted phonon modes. The targeted modes are selected in both energy and momentum by adjusting the impurity's potential. We show how to use the detector to observe coherent density waves and to measure temperatures of the Bose gas down to the nano-Kelvin regime. We analyze how our scheme could be realized experimentally, including the possibility of using an array of multiple impurities to achieve greater precision from a single experimental run., Comment: 11+4 pages, 7 figures

We present the path-sum formulation for exact statistical inference of marginals on Gaussian graphical models of arbitrary topology. The path-sum formulation gives the covariance between each pair of variables as a branched continued fraction of finite depth and breadth. Our method originates from the closed-form resummation of infinite families of terms of the walk-sum representation of the covariance matrix. We prove that the path-sum formulation always exists for models whose covariance matrix is positive definite: i.e.~it is valid for both walk-summable and non-walk-summable graphical models of arbitrary topology. We show that for graphical models on trees the path-sum formulation is equivalent to Gaussian belief propagation. We also recover, as a corollary, an existing result that uses determinants to calculate the covariance matrix. We show that the path-sum formulation formulation is valid for arbitrary partitions of the inverse covariance matrix. We give detailed examples demonstrating our results.

We present the path-sum formulation for $\mathsf{OE}[\mathsf{H}](t',t)=\mathcal{T}\,\text{exp}\big(\int_{t}^{t'}\!\mathsf{H}(\tau)\,d\tau\big)$, the time-ordered exponential of a time-dependent matrix $\mathsf{H}(t)$. The path-sum formulation gives $\mathsf{OE}[\mathsf{H}]$ as a branched continued fraction of finite depth and breadth. The terms of the path-sum have an elementary interpretation as self-avoiding walks and self-avoiding polygons on a graph. Our result is based on a representation of the time-ordered exponential as the inverse of an operator, the mapping of this inverse to sums of walks on graphs and the algebraic structure of sets of walks. We give examples demonstrating our approach. We establish a super-exponential decay bound for the magnitude of the entries of the time-ordered exponential of sparse matrices. We give explicit results for matrices with commonly encountered sparse structures.

In this work we analyze the simultaneous emergence of diffusive energy transport and local thermalization in a nonequilibrium one-dimensional quantum system, as a result of integrability breaking. Specifically, we discuss the local properties of the steady state induced by thermal boundary driving in a XXZ spin chain with staggered magnetic field. By means of efficient large-scale matrix product simulations of the equation of motion of the system, we calculate its steady state in the long-time limit.We start by discussing the energy transport supported by the system, finding it to be ballistic in the integrable limit and diffusive when the staggered field is finite. Subsequently, we examine the reduced density operators of neighboring sites and find that for large systems they are well approximated by local thermal states of the underlying Hamiltonian in the nonintegrable regime, even for weak staggered fields. In the integrable limit, on the other hand, this behavior is lost, and the identification of local temperatures is no longer possible. Our results agree with the intuitive connection between energy diffusion and thermalization., Comment: 13 pages, 13 figures

Estimating the expected value of an observable appearing in a non-equilibrium stochastic process usually involves sampling. If the observable's variance is high, many samples are required. In contrast, we show that performing the same task without sampling, using tensor network compression, efficiently captures high variances in systems of various geometries and dimensions. We provide examples for which matching the accuracy of our efficient method would require a sample size scaling exponentially with system size. In particular, the high variance observable $\mathrm{e}^{-\beta W}$, motivated by Jarzynski's equality, with $W$ the work done quenching from equilibrium at inverse temperature $\beta$, is exactly and efficiently captured by tensor networks., Comment: 7 pages, 3 figures, including supplemental material

We use midinfrared pulses with stable carrier-envelope phase offset to drive molecular vibrations in the charge transfer salt ET-F2TCNQ, a prototypical one-dimensional Mott insulator. We find that the Mott gap, which is probed resonantly with 10 fs laser pulses, oscillates with the pump field. This observation reveals that molecular excitations can coherently perturb the electronic on-site interactions (Hubbard U) by changing the local orbital wave function. The gap oscillates at twice the frequency of the vibrational mode, indicating that the molecular distortions couple quadratically to the local charge density.

Quantum simulators are devices that actively use quantum effects to answer questions about model systems and, through them, real systems. Here we expand on this definition by answering several fundamental questions about the nature and use of quantum simulators. Our answers address two important areas. First, the difference between an operation termed simulation and another termed computation. This distinction is related to the purpose of an operation, as well as our confidence in and expectation of its accuracy. Second, the threshold between quantum and classical simulations. Throughout, we provide a perspective on the achievements and directions of the field of quantum simulation., Comment: 13 pages, 2 figures

We characterise and study dynamical localisation of a finite interacting quantum many-body system. We present explicit bounds on the disorder strength required for the onset of localisation of the dynamics of arbitrary ensemble of sites of the XYZ spin-1/2 model. We obtain these results using a novel form of the fractional moment criterion, which we establish, together with a generalisation of the self-avoiding walk representation of the system Green's functions, called path-sums. These techniques are not specific to the XYZ model and hold in a much more general setting. We further present bounds for two observable quantities in the localised regime: the magnetisation of any sublattice of the system as well as the linear magnetic response function of the system. We confirm our results through numerical simulations., Comment: 35 pages; 5 figures

We study the non-equilibrium transport properties of a highly anisotropic two-dimensional lattice of spin-1/2 particles governed by a Heisenberg XXZ Hamiltonian. The anisotropy of the lattice allows us to approximate the system at finite temperature as an array of incoherently coupled one-dimensional chains. We show that in the regime of strong intrachain interactions, the weak interchain coupling considerably boosts spin transport in the driven system. Interestingly, we show that this enhancement increases with the length of the chains, which is related to superdiffusive spin transport. We describe the mechanism behind this effect, compare it to a similar phenomenon in single chains induced by dephasing, and explain why the former is much stronger.

Femtosecond relaxation of photo-excited quasiparticles in the one dimensional Mott insulator ET-F2TCNQ are measured as a function of external pressure, which is used to tune the electronic structure. By fitting the static optical properties and measuring femtosecond decay times at each pressure value, we correlate the relaxation rates with the electronic bandwidth t and on the intersite correlation energy V. The scaling of relaxation times with microscopic parameters is different than for metals and semiconductors. The competition between localization and delocalization of the Mott-Hubbard exciton dictates the efficiency of the decay, as exposed by a fit based on the solution of the time-dependent extended Hubbard Hamiltonian., Comment: 24 pages, 7 figures, final version including supplementary materials

We study the origin of "frozen" states in coupled Jaynes-Cummings-Hubbard arrays in the presence of losses. For the case of half the array initially populated with photons while the other half is left empty we show the emergence of self-localized photon or "frozen" states for specific values of the local atom-photon coupling. We analyze the dynamics in the quantum regime and discover important additional features appear not captured by a semiclassical treatment, which we analyze for different array sizes and filling fractions. We trace the origin of this interaction-induced photon "freezing" to the suppression of excitation of propagating modes in the system at large interaction strengths. We discuss in detail the possibility to experimentally probe the relevant transition by analyzing the emitted photon correlations. We find a strong signature of the effect in the emitted photons statistics., Comment: 5 pages, 6 figures. Contribution to the JOSA B special issue celebrating Jaynes-Cummings physics

In this work we study the heat transport in an XXZ spin-1/2 Heisenberg chain with homogeneous magnetic field, incoherently driven out of equilibrium by reservoirs at the boundaries. We focus on the effect of bulk dephasing (energy-dissipative) processes in different parameter regimes of the system. The non-equilibrium steady state of the chain is obtained by simulating its evolution under the corresponding Lindblad master equation, using the time evolving block decimation method. In the absence of dephasing, the heat transport is ballistic for weak interactions, while being diffusive in the strongly-interacting regime, as evidenced by the heat-current scaling with the system size. When bulk dephasing takes place in the system, diffusive transport is induced in the weakly-interacting regime, with the heat current monotonically decreasing with the dephasing rate. In contrast, in the strongly-interacting regime, the heat current can be significantly enhanced by dephasing for systems of small size.

We derive ab initio local Hubbard models for several optical lattice potentials of current interest, including the honeycomb and Kagom\'{e} lattices, verifying their accuracy on each occasion by comparing the interpolated band structures against the originals. To achieve this, we calculate the maximally-localized generalized Wannier basis by implementing the steepest-descent algorithm of Marzari and Vanderbilt [N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 (1997)] directly in one and two dimensions. To avoid local minima we develop an initialization procedure that is both robust and requires no prior knowledge of the optimal Wannier basis. The MATLAB code that implements our full procedure is freely available online at http://ccpforge.cse.rl.ac.uk/gf/project/mlgws/., Comment: 13 pages, 14 figures

A key insight from recent studies is that noise, such as dephasing, can improve the efficiency of quantum transport by suppressing coherent single-particle interference effects. However, it is not yet clear whether dephasing can enhance transport in an interacting many-body system. Here, we address this question by analyzing the transport properties of a boundary driven spinless fermion chain with nearest-neighbor interactions subject to bulk dephasing. The many-body nonequilibrium stationary state is determined using large-scale matrix product simulations of the corresponding quantum master equation. We find dephasing enhanced transport only in the strongly interacting regime, where it is shown to induce incoherent transitions bridging the gap between bound dark states and bands of mobile eigenstates. The generic nature of the transport enhancement is illustrated by a simple toy model, which contains the basic elements required for its emergence. Surprisingly, the effect is significant even in the linear response regime of the full system, and it is predicted to exist for any large and finite chain. The response of the system to dephasing also establishes a signature of an underlying nonequilibrium phase transition between regimes of transport degradation and enhancement. The existence of this transition is shown not to depend on the integrability of the model considered. As a result, dephasing enhanced transport is expected to persist in more realistic nonequilibrium strongly correlated systems., Comment: 15 pages, 17 figures

Most available theories for correlated electron transport are based on the Hubbard Hamiltonian. In this effective theory, renormalized hopping and interaction parameters only implicitly incorporate the coupling of correlated charge carriers to microscopic degrees of freedom. Unfortunately, no spectroscopy can individually probe such renormalizations, limiting the applicability of Hubbard models. We show here that the role of each individual degree of freedom can be made explicit by using a new experimental technique, which we term 'quantum modulation spectroscopy' and we demonstrate here in the one-dimensional Mott insulator ET-F2TCNQ. We explore the role on the charge hopping of two localized molecular modes, which we drive with a mid infrared optical pulse. Sidebands appear in the modulated optical spectrum, and their linshape is fitted with a model based on the dynamic Hubbard Hamiltonian. A striking asymmetry between the renormalization of doublons and holons is revealed. The concept of quantum modulation spectroscopy can be used to systematically deconstruct Hubbard Hamiltonians in many materials, exposing the role of any mode, electronic or magnetic, that can be driven to large amplitude with a light field., Comment: 20 pages, 5 figures

Simulating quantum circuits using classical computers lets us analyse the inner workings of quantum algorithms. The most complete type of simulation, strong simulation, is believed to be generally inefficient. Nevertheless, several efficient strong simulation techniques are known for restricted families of quantum circuits and we develop an additional technique in this article. Further, we show that strong simulation algorithms perform another fundamental task: solving search problems. Efficient strong simulation techniques allow solutions to a class of search problems to be counted and found efficiently. This enhances the utility of strong simulation methods, known or yet to be discovered, and extends the class of search problems known to be efficiently simulable. Relating strong simulation to search problems also bounds the computational power of efficiently strongly simulable circuits; if they could solve all problems in $\mathrm{P}$ this would imply the collapse of the complexity hierarchy $\mathrm{P} \subseteq \mathrm{NP} \subseteq # \mathrm{P}$., Comment: 9 pages, 3 figures

Re-entrance is a novel feature where the phase boundaries of a system exhibit a succession of transitions between two phases A and B, like A-B-A-B, when just one parameter is varied monotonically. This type of re-entrance is displayed by the 1D Bose Hubbard model between its Mott insulator (MI) and superfluid phase as the hopping amplitude is increased from zero. Here we analyse this counter-intuitive phenomenon directly in the thermodynamic limit by utilizing the infinite time-evolving block decimation algorithm to variationally minimize an infinite matrix product state (MPS) parameterized by a matrix size chi. Exploiting the direct restriction on the half-chain entanglement imposed by fixing chi, we determined that re-entrance in the MI lobes only emerges in this approximate when chi >= 8. This entanglement threshold is found to be coincident with the ability an infinite MPS to be simultaneously particle-number symmetric and capture the kinetic energy carried by particle-hole excitations above the MI. Focussing on the tip of the MI lobe we then applied, for the first time, a general finite-entanglement scaling analysis of the infinite order Kosterlitz-Thouless critical point located there. By analysing chi's up to a very moderate chi = 70 we obtained an estimate of the KT transition as t_KT = 0.30 +/- 0.01, demonstrating the how a finite-entanglement approach can provide not only qualitative insight but also quantitatively accurate predictions., Comment: 12 pages, 8 figures

We study the non-equilibrium behavior of optically driven dissipative coupled resonator arrays. Assuming each resonator is coupled with a two-level system via a Jaynes-Cummings interaction, we calculate the many-body steady state behavior of the system under coherent pumping and dissipation. We propose and analyze the many-body phases using experimentally accessible quantities such as the total excitation number, the emitted photon spectra and photon coherence functions for different parameter regimes. In parallel, we also compare and contrast the expected behavior of this system assuming the local nonlinearity in the cavities is generated by a generic Kerr effect rather than a Jaynes-Cummings interaction. We find that the behavior of the experimentally accessible observables produced by the two models differs for realistic regimes of interactions even when the corresponding nonlinearities are of similar strength. We analyze in detail the extra features available in the Jaynes-Cummings-Hubbard (JCH) model originating from the mixed nature of the excitations and investigate the regimes where the Kerr approximation would faithfully match the JCH physics. We find that the latter is true for values of the light-matter coupling and losses beyond the reach of current technology. Throughout the study we operate in the weak pumping, fully quantum mechanical regime where approaches such as mean field theory fail, and instead use a combination of quantum trajectories and the time evolving block decimation algorithm to compute the relevant steady state observables. In our study we have assumed small to medium size arrays (from 3 up to 16 sites) and values of the ratio of coupling to dissipation rate $g/\gamma \sim 20$ which makes our results implementable with current designs in Circuit QED and with near future photonic crystal set ups., Comment: 22 pages, 6 figures

We show that the series of all walks between any two vertices of any (possibly weighted) directed graph $\mathcal{G}$ is given by a universal continued fraction of finite depth and breadth involving the simple paths and simple cycles of $\mathcal{G}$. A simple path is a walk forbidden to visit any vertex more than once. We obtain an explicit formula giving this continued fraction. Our results are based on an equivalent to the fundamental theorem of arithmetic: we demonstrate that arbitrary walks on $\mathcal{G}$ factorize uniquely into nesting products of simple paths and simple cycles, where nesting is a product operation between walks that we define. We show that the simple paths and simple cycles are the prime elements of the set of all walks on $\mathcal{G}$ equipped with the nesting product. We give an algorithm producing the prime factorization of individual walks, and obtain a recursive formula producing the prime factorization of sets of walks. Our results have already found applications in machine learning, matrix computations and quantum mechanics., Comment: Updated with links between nesting and loop-erasing. Still under review (!)

We introduce the method of path-sums which is a tool for exactly evaluating a function of a discrete matrix with possibly non-commuting entries, based on the closed-form resummation of infinite families of terms in the corresponding Taylor series. If the matrix is finite, our approach yields the exact result in a finite number of steps. We achieve this by combining a mapping between matrix powers and walks on a weighted directed graph with a universal graph-theoretic result on the structure of such walks. We present path-sum expressions for a matrix raised to a complex power, the matrix exponential, matrix inverse, and matrix logarithm. We show that the quasideterminants of a matrix can be naturally formulated in terms of a path-sum, and present examples of the application of the path-sum method. We show that obtaining the inversion height of a matrix inverse and of quasideterminants is an NP-complete problem., Comment: 23 pages, light version submitted to SIAM Journal on Matrix Analysis and Applications (SIMAX). A separate paper with the graph theoretic results is available at: arXiv:1202.5523v1. Results for matrices over division rings will be published separately as well

Motivated by a recent experiment [J. Catani et al., arXiv:1106.0828v1 preprint, 2011], we study breathing oscillations in the width of a harmonically trapped impurity interacting with a separately trapped Bose gas. We provide an intuitive physical picture of such dynamics at zero temperature, using a time-dependent variational approach. In the Gross-Pitaevskii regime we obtain breathing oscillations whose amplitudes are suppressed by self trapping, due to interactions with the Bose gas. Introducing phonons in the Bose gas leads to the damping of breathing oscillations and non-Markovian dynamics of the width of the impurity, the degree of which can be engineered through controllable parameters. Our results reproduce the main features of the impurity dynamics observed by Catani et al. despite experimental thermal effects, and are supported by simulations of the system in the Gross-Pitaevskii regime. Moreover, we predict novel effects at lower temperatures due to self-trapping and the inhomogeneity of the trapped Bose gas., Comment: 7 pages, 3 figures

We consider the excitation of electrons in semiconducting carbon nanotubes by photons from the evanescent field created by a subwavelength-diameter optical fiber. The strongly changing evanescent field of such nanofibers requires dropping the dipole approximation. We show that this leads to novel effects, especially a high dependence of the photon absorption on the relative orientation and geometry of the nanotube-nanofiber setup in the optical and near infrared domain. In particular, we calculate photon absorption probabilities for a straight nanotube and nanofiber depending on their relative angle. Nanotubes orthogonal to the fiber are found to perform much better than parallel nanotubes when they are short. As the nanotube gets longer the absorption of parallel nanotubes is found to exceed the orthogonal nanotubes and approach 100% for extremely long nanotubes. In addition, we show that if the nanotube is wrapped around the fiber in an appropriate way the absorption is enhanced. We find that optical and near infrared photons could be converted to excitations with efficiencies that may exceed 90%. This may provide opportunities for future photodetectors and we discuss possible setups., Comment: 14 pages, 14 figures

We show that transient supersolid quantum states of Rydberg-excitations can be created dynamically from a Mott insulator of ground state atoms in a 2D optical-lattices by irradiating it with short laser pulses. The structure of these supersolids is tunable via the choice of laser parameters. We calculate first, second and fourth order correlation functions as well as the pressure to characterize the supersolid states. Our study is based on the development of a general theoretical tool for obtaining the dynamics of strongly interacting quantum systems whose initial state is accurately known. We show that this method allows to accurately approximate the evolution of quantum systems analytically with a number of operations growing polynomially., Comment: 2 figures

We adapt the bialgebra and Hopf relations to expose internal structure in the ground state of a Hamiltonian with $Z_2$ topological order. Its tensor network description allows for exact contraction through simple diagrammatic rewrite rules. The contraction property does not depend on specifics such as geometry, but rather originates from the non-trivial algebraic properties of the constituent tensors. We then generalise the resulting tensor network from a spin-1/2 lattice to a class of exactly contractible states on spin-S degrees of freedom, yielding the most efficient tensor network description of finite Abelian lattice gauge theories. We gain a new perspective on these states as examples of two-dimensional quantum states with algebraically contractible tensor network representations. The introduction of local perturbations to the network is shown to reduce the von Neumann entropy of string-like regions, creating an unentangled sub-system within the bulk in a certain limit. We also show how perturbations induce finite-range correlations in this system. This class of tensor networks is readily translated onto any lattice, and we differentiate between the physical consequences of bipartite and non-bipartite lattices on the properties of the corresponding quantum states. We explicitly show this on the hexagonal, square, kagome and triangular lattices., Comment: 23 pages, 6 figures

We study the suitability of correlator product states for describing ground-state properties of two-dimensional spin models. Our ansatz for the many-body wave function takes the form of either plaquette or bond correlator product states and the energy is optimized by varying the correlators using Monte Carlo minimization. For the Ising model we find that plaquette correlators are best for estimating the energy while bond correlators capture the expected long-range correlations and critical behavior of the system more faithfully. For the antiferromagnetic Heisenberg model, however, plaquettes outperform bond correlators at describing both local and long-range correlations because of the substantially larger number of local parameters they contain. These observations have quantitative implications for the application of correlator product states to other more complex systems, and give important heuristic insights: in particular the necessity of carefully tailoring the choice of correlators to the system considered, its interactions and symmetries., Comment: 14 pages, 13 figures, clarifying amendments made, typos corrected, legend of Figure 8 corrected, references added

Using near-exact numerical simulations we study the propagation of an impurity through a one-dimensional Bose lattice gas for varying bosonic interaction strengths and filling factors at zero temperature. The impurity is coupled to the Bose gas and confined to a separate tilted lattice. The precise nature of the transport of the impurity is specific to the excitation spectrum of the Bose gas which allows one to measure properties of the Bose gas non-destructively, in principle, by observing the impurity; here we focus on the spatial and momentum distributions of the impurity as well as its reduced density matrix. For instance we show it is possible to determine whether the Bose gas is commensurately filled as well as the bandwidth and gap in its excitation spectrum. Moreover, we show that the impurity acts as a witness to the cross-over of its environment from the weakly to the strongly interacting regime, i.e., from a superfluid to a Mott insulator or Tonks-Girardeau lattice gas and the effects on the impurity in both of these strongly-interacting regimes are clearly distinguishable. Finally, we find that the spatial coherence of the impurity is related to its propagation through the Bose gas, giving an experimentally controllable example of noise-enhanced quantum transport., Comment: 11 pages, 7 figures

We consider 2D networks composed of nodes initially linked by two-qubit mixed states. In these networks we develop a global error correction scheme that can generate distance-independent entanglement from arbitrary network geometries using rank two states. By using this method and combining it with the concept of percolation we also show that the generation of long distance entanglement is possible with rank three states. Entanglement percolation and global error correction have different advantages depending on the given situation. To reveal the trade-off between them we consider their application on networks containing pure states. In doing so we find a range of pure-state schemes, each of which has applications in particular circumstances: For instance, we can identify a protocol for creating perfect entanglement between two distant nodes. However, this protocol can not generate a singlet between any two nodes. On the other hand, we can also construct schemes for creating entanglement between any nodes, but the corresponding entanglement fidelity is lower., Comment: 10 pages, 9 figures, 1 table

Arrays of atoms trapped in optical lattices are appealing as storage media for photons, since motional dephasing of the atoms is eliminated. The regular lattice is also associated with band structure in the dispersion experienced by incident photons. Here we study the influence of this band structure on the efficiency of quantum memories based on electromagnetically induced transparency (EIT) and on Raman absorption. We observe a number of interesting effects, such as both reduced and superluminal group velocities, enhanced atom-photon coupling and anomalous transmission. These effects are ultimately deleterious to the memory efficiency, but they are easily avoided by tuning the optical fields away from the band edges., Comment: 6 pages, 5 figures, accepted for publication in PRA

We present an analysis of Bose-Fermi mixtures in optical lattices for the case where the lattice potential of the fermions is tilted and the bosons (in the superfluid phase) are described by Bogoliubov phonons. It is shown that the Bogoliubov phonons enable hopping transitions between fermionic Wannier-Stark states; these transitions are accompanied by energy dissipation into the superfluid and result in a net atomic current along the lattice. We derive a general expression for the drift velocity of the fermions and find that the dependence of the atomic current on the lattice tilt exhibits negative differential conductance and phonon resonances. Numerical simulations of the full dynamics of the system based on the time-evolving block decimation algorithm reveal that the phonon resonances should be observable under the conditions of a realistic measuring procedure., Comment: 8 pages, 5 figures

We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of 1D quantum systems, to simulate the time-evolution of non-equilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time-evolution is efficiently simulated by repeatedly updating and approximately re-factorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible, even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size., Comment: 16 pages, 10 figures

Relativistic causality has dramatic consequences on the measurability of nonlocal variables and poses the fundamental question of whether it is physically meaningful to speak about the value of nonlocal variables at a particular time. Recent work has shown that by weakening the role of the measurement in preparing eigenstates of the variable it is in fact possible to measure all nonlocal observables instantaneously by exploiting entanglement. However, for these measurement schemes to succeed with certainty an infinite amount of entanglement must be distributed initially and all this entanglement is necessarily consumed. In this work we sharpen the characterisation of instantaneous nonlocal measurements by explicitly devising schemes in which only a finite amount of the initially distributed entanglement is ever utilised. This enables us to determine an upper bound to the average consumption for the most general cases of nonlocal measurements. This includes the tasks of state verification, where the measurement verifies if the system is in a given state, and verification measurements of a general set of eigenstates of an observable. Despite its finiteness the growth of entanglement consumption is found to display an extremely unfavourable exponential of an exponential scaling with either the number of qubits needed to contain the Schmidt rank of the target state or total number of qubits in the system for an operator measurement. This scaling is seen to be a consequence of the combination of the generic exponential scaling of unitary decompositions combined with the highly recursive structure of our scheme required to overcome the no-signalling constraint of relativistic causality., Comment: 32 pages and 14 figures. Updated to published version

We study the generation of singlets in quantum networks with nodes initially sharing a finite number of partially entangled bipartite mixed states. We prove that singlets between arbitrary nodes in such networks can be created if and only if the initial states connecting the nodes have a particular form. We then generalize the method of entanglement percolation, previously developed for pure states, to mixed states of this form. As part of this, we find and compare different distillation protocols necessary to convert groups of mixed states shared between neighboring nodes of the network into singlets. In addition, we discuss protocols that only rely on local rules for the efficient connection of two remote nodes in the network via entanglement swapping. Further improvements of the success probability of singlet generation are developed by using particular forms of `quantum preprocessing' on the network. This includes generalized forms of entanglement swapping and we show how such strategies can be embedded in regular and hierarchical quantum networks., Comment: 17 pages, 21 figures

Quantum memories, capable of controllably storing and releasing a photon, are a crucial component for quantum computers and quantum communications. So far, quantum memories have operated with bandwidths that limit data rates to MHz. Here we report the coherent storage and retrieval of sub-nanosecond low intensity light pulses with spectral bandwidths exceeding 1 GHz in cesium vapor. The novel memory interaction takes place via a far off-resonant two-photon transition in which the memory bandwidth is dynamically generated by a strong control field. This allows for an increase in data rates by a factor of almost 1000 compared to existing quantum memories. The memory works with a total efficiency of 15% and its coherence is demonstrated by directly interfering the stored and retrieved pulses. Coherence times in hot atomic vapors are on the order of microsecond - the expected storage time limit for this memory., Comment: 13 pages, 5 figures