Source author record

S. R. Clark

S. R. Clark 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

34works
11topics
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

34 published item(s)

preprint2016arXiv

Beyond mean-field bistability in driven-dissipative lattices: bunching-antibunching transition and quantum simulation

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.

preprint2016arXiv

Non-linear quantum-classical scheme to simulate non-equilibrium strongly correlated fermionic many-body dynamics

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.

preprint2016arXiv

Thermometry of ultracold atoms via non-equilibrium work distributions

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.

preprint2016arXiv

Transport enhancement from incoherent coupling between one-dimensional quantum conductors

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.

preprint2015arXiv

An optically stimulated superconducting-like phase in K3C60 far above equilibrium Tc

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.

preprint2015arXiv

Capturing exponential variance using polynomial resources: applying tensor networks to non-equilibrium stochastic processes

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}^{-βW}$, motivated by Jarzynski's equality, with $W$ the work done quenching from equilibrium at inverse temperature $β$, is exactly and efficiently captured by tensor networks.

preprint2015arXiv

Coexistence of energy diffusion and local thermalization in nonequilibrium XXZ spin chains with integrability breaking

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.

preprint2015arXiv

Dephasing enhanced transport in nonequilibrium strongly correlated quantum systems

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.

preprint2015arXiv

Hubbard model for atomic impurities bound by the vortex lattice of a rotating BEC

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

preprint2015arXiv

Pressure dependent relaxation in the photo-excited Mott insulator ETF2TCNQ: Influence of hopping and correlations on quasiparticle recombination rates

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.

preprint2015arXiv

Proposed parametric cooling of bilayer cuprate superconductors by terahertz excitation

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.

preprint2015arXiv

Quantum interference between charge excitation paths in a solid state Mott insulator

The competition between electron localization and de-localization in Mott insulators underpins the physics of strongly-correlated electron systems. Photo-excitation, which re-distributes charge between sites, can control this many-body process on the ultrafast timescale. To date, time-resolved studies have been performed in solids in which other degrees of freedom, such as lattice, spin, or orbital excitations come into play. However, the underlying quantum dynamics of bare electronic excitations has remained out of reach. Quantum many-body dynamics have only been detected in the controlled environment of optical lattices where the dynamics are slower and lattice excitations are absent. By using nearly-single-cycle near-IR pulses, we have measured coherent electronic excitations in the organic salt ET-F2TCNQ, a prototypical one-dimensional Mott Insulator. After photo-excitation, a new resonance appears on the low-energy side of the Mott gap, which oscillates at 25 THz. Time-dependent simulations of the Mott-Hubbard Hamiltonian reproduce the oscillations, showing that electronic delocalization occurs through quantum interference between bound and ionized holon-doublon pairs.

preprint2015arXiv

Spatially resolved ultrafast magnetic dynamics launched at a complex-oxide hetero-interface

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.

preprint2015arXiv

THz-Frequency Modulation of the Hubbard U in an Organic Mott Insulator

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.

preprint2015arXiv

Total correlations of the diagonal ensemble herald the many-body localization transition

The intriguing phenomenon of many-body localization (MBL) has attracted significant interest recently, but a complete characterization is still lacking. In this work, we introduce the total correlations, a concept from quantum information theory capturing multi-partite correlations, to the study of this phenomenon. We demonstrate that the total correlations of the diagonal ensemble provides a meaningful diagnostic tool to pin-down, probe, and better understand the MBL transition and ergodicity breaking in quantum systems. In particular, we show that the total correlations has sub-linear dependence on the system size in delocalized, ergodic phases, whereas we find that it scales extensively in the localized phase developing a pronounced peak at the transition. We exemplify the power of our approach by means of an exact diagonalization study of a Heisenberg spin chain in a disordered field.

preprint2014arXiv

What is a quantum simulator?

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.

preprint2013arXiv

Ab initio derivation of Hubbard models for cold atoms in optical lattices

We derive ab initio local Hubbard models for several optical lattice potentials of current interest, including the honeycomb and Kagomé 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/.

preprint2013arXiv

Extracting quantum work statistics and fluctuation theorems by single qubit interferometry

We propose an experimental scheme to verify the quantum non-equilibrium fluctuation relations using current technology. Specifically, we show that the characteristic function of the work distribution for a non-equilibrium quench of a general quantum system can be extracted from Ramsey interferometry of a single probe qubit. Our scheme paves the way for the full characterisation of non-equilibrium processes in a variety of complex quantum systems ranging from single particles to many-body atomic systems and spin chains. We demonstrate our idea using a time-dependent quench of the motional state of a trapped ion, where the internal pseudo-spin provides a convenient probe qubit.

preprint2013arXiv

Frozen photons in Jaynes Cummings arrays

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.

preprint2013arXiv

Heat transport in the $XXZ$ spin chain: from ballistic to diffusive regimes and dephasing enhancement

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.

preprint2012arXiv

Capturing the re-entrant behaviour of one-dimensional Bose-Hubbard model

The Bose Hubbard model (BHM) is an archetypal quantum lattice system exhibiting a quantum phase transition between its superfluid (SF) and Mott-insulator (MI) phase. Unlike in higher dimensions the phase diagram of the BHM in one dimension possesses regions in which increasing the hopping amplitude can result in a transition from MI to SF and then back to a MI. This type of re-entrance is well known in classical systems like liquid crystals yet its origin in quantum systems is still not well understood. Moreover, this unusual re-entrant character of the BHM is not easily captured in approximate analytical or numerical calculations. Here we study in detail the predictions of three different and widely used approximations; a multi-site mean-field decoupling, a finite-sized cluster calculation, and a real-space renormalization group (RG) approach. It is found that mean-field calculations do not reproduce re-entrance while finite-sized clusters display a precursor to re-entrance. Here we show for the first time that RG does capture the re-entrant feature and constitutes one of the simplest approximation able to do so. The differing abilities of these approaches reveals the importance of describing short-ranged correlations relevant to the kinetic energy of a MI in a particle-number symmetric way.

preprint2012arXiv

Deconstructing the Hubbard Hamiltonian by Ultrafast Quantum Modulation Spectroscopy in Solid-state Mott Insulators

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.

preprint2012arXiv

Non-equilibrium many-body effects in driven nonlinear resonator arrays

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/γ\sim 20$ which makes our results implementable with current designs in Circuit QED and with near future photonic crystal set ups.

preprint2012arXiv

Re-entrance and entanglement in the one-dimensional Bose-Hubbard model

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.

preprint2012arXiv

Solving search problems by strongly simulating quantum circuits

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}$.

preprint2011arXiv

Algebraically contractible topological tensor network states

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.

preprint2011arXiv

Breathing oscillations of a trapped impurity in a Bose gas

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.

preprint2011arXiv

Capturing long range correlations in two-dimensional quantum lattice systems using correlator product states

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.

preprint2011arXiv

Impurity transport through a strongly interacting bosonic quantum gas

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.

preprint2010arXiv

Dynamical simulations of classical stochastic systems using matrix product states

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.

preprint2010arXiv

Entanglement consumption of instantaneous nonlocal quantum measurements

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.

preprint2010arXiv

Phonon resonances in atomic currents through Bose-Fermi mixtures in optical lattices

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.

preprint2009arXiv

Exact matrix product solutions in the Heisenberg picture of an open quantum spin chain

In recent work Hartmann et al [Phys. Rev. Lett. 102, 057202 (2009)] demonstrated that the classical simulation of the dynamics of open 1D quantum systems with matrix product algorithms can often be dramatically improved by performing time evolution in the Heisenberg picture. For a closed system this was exemplified by an exact matrix product operator solution of the time-evolved creation operator of a quadratic fermi chain with a matrix dimension of just two. In this work we show that this exact solution can be significantly generalized to include the case of an open quadratic fermi chain subjected to master equation evolution with Lindblad operators that are linear in the fermionic operators. Remarkably even in this open system the time-evolution of operators continues to be described by matrix product operators with the same fixed dimension as that required by the solution of a coherent quadratic fermi chain for all times. Through the use of matrix product algorithms the dynamical behaviour of operators in this non-equilibrium open quantum system can be computed with a cost that is linear in the system size. We present some simple numerical examples which highlight how useful this might be for the more detailed study of open system dynamics. Given that Heisenberg picture simulations have been demonstrated to offer significant accuracy improvements for other open systems that are not exactly solvable our work also provides further insight into how and why this advantage arises.

preprint2007arXiv

Generation of twin Fock states via transition from a two-component Mott insulator to a superfluid

We propose the dynamical creation of twin Fock states, which exhibit Heisenberg limited interferometric phase sensitivities, in an optical lattice. In our scheme a two-component Mott insulator with two bosonic atoms per lattice site is melted into a superfluid. This process transforms local correlations between hyperfine states of atom pairs into multi-particle correlations extending over the whole system. The melting time does not scale with the system size which makes our scheme experimentally feasible.