Time-dependent currents of 1D bosons in an optical lattice

We analyse the time-dependence of currents in a 1D Bose gas in an optical lattice. For a 1D system, the stability of currents induced by accelerating the lattice exhibits a broad crossover as a function of the magnitude of the acceleration, and the strength of the inter-particle interactions. This differs markedly from mean-field results in higher dimensions. Using the infinite Time Evolving Block Decimation algorithm, we characterise this crossover by making quantitative predictions for the time-dependent behaviour of the currents and their decay rate. We also compute the time-dependence of quasi-condensate fractions which can be measured directly in experiments. We compare our results to calculations based on phase-slip methods, finding agreement with the scaling as the particle density increases, but with significant deviations near unit filling.Comment: 19 pages, 10 figure

Fault-Tolerant Dissipative Preparation of Atomic Quantum Registers with Fermions

We propose a fault tolerant loading scheme to produce an array of fermions in an optical lattice of the high fidelity required for applications in quantum information processing and the modelling of strongly correlated systems. A cold reservoir of Fermions plays a dual role as a source of atoms to be loaded into the lattice via a Raman process and as a heat bath for sympathetic cooling of lattice atoms. Atoms are initially transferred into an excited motional state in each lattice site, and then decay to the motional ground state, creating particle-hole pairs in the reservoir. Atoms transferred into the ground motional level are no longer coupled back to the reservoir, and doubly occupied sites in the motional ground state are prevented by Pauli blocking. This scheme has strong conceptual connections with optical pumping, and can be extended to load high-fidelity patterns of atoms.Comment: 12 pages, 7 figures, RevTex

Entanglement growth in quench dynamics with variable range interactions

Studying entanglement growth in quantum dynamics provides both insight into the underlying microscopic processes and information about the complexity of the quantum states, which is related to the efficiency of simulations on classical computers. Recently, experiments with trapped ions, polar molecules, and Rydberg excitations have provided new opportunities to observe dynamics with long-range interactions. We explore nonequilibrium coherent dynamics after a quantum quench in such systems, identifying qualitatively different behavior as the exponent of algebraically decaying spin-spin interactions in a transverse Ising chain is varied. Computing the build-up of bipartite entanglement as well as mutual information between distant spins, we identify linear growth of entanglement entropy corresponding to propagation of quasiparticles for shorter range interactions, with the maximum rate of growth occurring when the Hamiltonian parameters match those for the quantum phase transition. Counter-intuitively, the growth of bipartite entanglement for long-range interactions is only logarithmic for most regimes, i.e., substantially slower than for shorter range interactions. Experiments with trapped ions allow for the realization of this system with a tunable interaction range, and we show that the different phenomena are robust for finite system sizes and in the presence of noise. These results can act as a direct guide for the generation of large-scale entanglement in such experiments, towards a regime where the entanglement growth can render existing classical simulations inefficient.Comment: 17 pages, 7 figure

Dissipation-induced d-Wave Pairing of Fermionic Atoms in an Optical Lattice

We show how dissipative dynamics can give rise to pairing for two-component fermions on a lattice. In particular, we construct a "parent" Liouvillian operator so that a BCS-type state of a given symmetry, e.g. a d-wave state, is reached for arbitrary initial states in the absence of conservative forces. The system-bath couplings describe single-particle, number conserving and quasi-local processes. The pairing mechanism crucially relies on Fermi statistics. We show how such Liouvillians can be realized via reservoir engineering with cold atoms representing a driven dissipative dynamics.Comment: 5 pages, 3 figures. Replaced with the published versio

A Novel Approach to Study Highly Correlated Nanostructures: The Logarithmic Discretization Embedded Cluster Approximation

This work proposes a new approach to study transport properties of highly correlated local structures. The method, dubbed the Logarithmic Discretization Embedded Cluster Approximation (LDECA), consists of diagonalizing a finite cluster containing the many-body terms of the Hamiltonian and embedding it into the rest of the system, combined with Wilson's idea of a logarithmic discretization of the representation of the Hamiltonian. The physics associated with both one embedded dot and a double-dot side-coupled to leads is discussed in detail. In the former case, the results perfectly agree with Bethe ansatz data, while in the latter, the physics obtained is framed in the conceptual background of a two-stage Kondo problem. A many-body formalism provides a solid theoretical foundation to the method. We argue that LDECA is well suited to study complicated problems such as transport through molecules or quantum dot structures with complex ground states.Comment: 17 pages, 13 figure

Spin-charge separation in two-component Bose-gases

We show that one of the key characteristics of interacting one-dimensional electronic quantum systems, the separation of spin and charge, can be observed in a two-component system of bosonic ultracold atoms even close to a competing phase separation regime. To this purpose we determine the real-time evolution of a single particle excitation and the single-particle spectral function using density-matrix renormalization group techniques. Due to efficient bosonic cooling and good tunability this setup exhibits very good conditions for observing this strong correlation effect. In anticipation of experimental realizations we calculate the velocities for spin and charge perturbations for a wide range of parameters

Modulation spectroscopy with ultracold fermions in an optical lattice

We propose an experimental setup of ultracold fermions in an optical lattice to determine the pairing gap in a superfluid state and the spin ordering in a Mott-insulating state. The idea is to apply a periodic modulation of the lattice potential and to use the thereby induced double occupancy to probe the system. We show by full time-dependent calculation using the adaptive time dependent density-matrix renormalization group method that the position of the peak in the spectrum of the induced double occupancy gives the pairing energy in a superfluid and the interaction energy in a Mott-insulator, respectively. In the Mott-insulator we relate the spectral weight of the peak to the spin ordering at finite temperature using perturbative calculations

High order non-unitary split-step decomposition of unitary operators

We propose a high order numerical decomposition of exponentials of hermitean operators in terms of a product of exponentials of simple terms, following an idea which has been pioneered by M. Suzuki, however implementing it for complex coefficients. We outline a convenient fourth order formula which can be written compactly for arbitrary number of noncommuting terms in the Hamiltonian and which is superiour to the optimal formula with real coefficients, both in complexity and accuracy. We show asymptotic stability of our method for sufficiently small time step and demonstrate its efficiency and accuracy in different numerical models.Comment: 10 pages, 4 figures (5 eps files) Submitted to J. of Phys. A: Math. Ge

Don't break a leg: Running birds from quail to ostrich prioritise leg safety and economy in uneven terrain

Cursorial ground birds are paragons of bipedal running that span a 500-fold mass range from quail to ostrich. Here we investigate the task-level control priorities of cursorial birds by analysing how they negotiate single-step obstacles that create a conflict between body stability (attenuating deviations in body motion) and consistent leg force–length dynamics (for economy and leg safety). We also test the hypothesis that control priorities shift between body stability and leg safety with increasing body size, reflecting use of active control to overcome size-related challenges. Weight-support demands lead to a shift towards straighter legs and stiffer steady gait with increasing body size, but it remains unknown whether non-steady locomotor priorities diverge with size. We found that all measured species used a consistent obstacle negotiation strategy, involving unsteady body dynamics to minimise fluctuations in leg posture and loading across multiple steps, not directly prioritising body stability. Peak leg forces remained remarkably consistent across obstacle terrain, within 0.35 body weights of level running for obstacle heights from 0.1 to 0.5 times leg length. All species used similar stance leg actuation patterns, involving asymmetric force–length trajectories and posture-dependent actuation to add or remove energy depending on landing conditions. We present a simple stance leg model that explains key features of avian bipedal locomotion, and suggests economy as a key priority on both level and uneven terrain. We suggest that running ground birds target the closely coupled priorities of economy and leg safety as the direct imperatives of control, with adequate stability achieved through appropriately tuned intrinsic dynamics

Mean-Field Interacting Boson Random Point Fields in Weak Harmonic Traps

A model of the mean-field interacting boson gas trapped by a weak harmonic potential is considered by the \textit{boson random point fields} methods. We prove that in the Weak Harmonic Trap (WHT) limit there are two phases distinguished by the boson condensation and by a different behaviour of the local particle density. For chemical potentials less than a certain critical value, the resulting Random Point Field (RPF) coincides with the usual boson RPF, which corresponds to a non-interacting (ideal) boson gas. For the chemical potentials greater than the critical value, the boson RPF describes a divergent (local) density, which is due to \textit{localization} of the macroscopic number of condensed particles. Notice that it is this kind of transition that observed in experiments producing the Bose-Einstein Condensation in traps