Source author record

Masudul Haque

Masudul Haque 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

52works
16topics
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

52 published item(s)

preprint2026arXiv

Universality, Robustness, and Limits of the Eigenstate Thermalization Hypothesis in Open Quantum Systems

The eigenstate thermalization hypothesis (ETH) underpins much of our modern understanding of the thermalization of closed quantum many-body systems. Here, we investigate the statistical properties of observables in the eigenbasis of the Lindbladian operator of a Markovian open quantum system. We demonstrate the validity of a Lindbladian ETH ansatz through extensive numerical simulations of several physical models. To highlight the robustness of Lindbladian ETH, we consider what we dub the dilute-click regime of the model, in which one postselects only quantum trajectories with a finite fraction of quantum jumps. The average dynamics are generated by a non-trace-preserving Liouvillian, and we show that the Lindbladian ETH ansatz still holds in this case. On the other hand, the no-click limit is a singular point at which the Lindbladian reduces to a doubled non-Hermitian Hamiltonian and Lindbladian ETH breaks down.

preprint2022arXiv

A bipartite Sachdev-Ye-Kitaev model: Conformal limit and level statistics

We study a bipartite version of the Sachdev-Ye-Kitaev (SYK) model. We show that the model remains solvable in the limit of large-$N$ in the same sense as the original model if the ratio of both flavors is kept finite. The scaling dimensions of the two species can be tuned continuously as a function of the ratio. We also investigate the finite-size spectral properties of the model. We show how the level statistics differs from the original SYK model and infer an additional exchange symmetry in the bipartite model.

preprint2022arXiv

Entanglement of midspectrum eigenstates of chaotic many-body systems: Reasons for deviation from random ensembles

Eigenstates of local many-body interacting systems that are far from spectral edges are thought to be ergodic and close to being random states. This is consistent with the eigenstate thermalization hypothesis and volume-law scaling of entanglement. We point out that systematic departures from complete randomness are generically present in mid-spectrum eigenstates, and focus on the departure of the entanglement entropy from the random-state prediction. We show that the departure is (partly) due to spatial correlations and due to orthogonality to the eigenstates at the spectral edge, which imposes structure on the mid-spectrum eigenstates.

preprint2022arXiv

Quantum Noise as a Symmetry-Breaking Field

We investigate the effect of quantum noise on the measurement-induced quantum phase transition in monitored random quantum circuits. Using the efficient simulability of random Clifford circuits, we find that the transition is broadened into a crossover and that the phase diagram as a function of projective measurements and noise exhibits several distinct regimes. We show that a mapping to a classical statistical mechanics problem accounts for the main features of the random circuit phase diagram. The bulk noise maps to an explicit permutation symmetry breaking coupling; this symmetry is spontaneously broken when the noise is switched off. These results have implications for the realization of entanglement transitions in noisy quantum circuits.

preprint2021arXiv

Diffusive Operator Spreading for Random Unitary Free Fermion Circuits

We study a model of free fermions on a chain with dynamics generated by random unitary gates acting on nearest neighbor bonds and present an exact calculation of time-ordered and out-of-time-ordered correlators. We consider three distinct cases: the random circuit with spatio-temporal disorder (i) with and (ii) without particle number conservation and (iii) the particle non-conserving case with purely temporal disorder. In all three cases, temporal disorder causes diffusive operator spreading and $\sim\sqrt{t}$ entanglement growth. This is in sharp contrast to Anderson localization for the case of static disorder and with the ballistic behavior observed in both the clean case of Hamiltonian evolution and in fully random unitary quantum circuits.

preprint2020arXiv

Few-particle dynamics of fractional quantum Hall lattice models

Considering lattice Hamiltonians designed using conformal field theory to have fractional quantum Hall states as ground states, we study the dynamics of one or two particles on such lattices. Examining the eigenspectrum and dynamics of the single-particle sector, we demonstrate that these Hamiltonians cannot be regarded as describing interacting particles placed on a Chern band, so that the physics is fundamentally different from fractional Chern insulators. The single-particle spectrum is shown to consist of eigenstates localized in shells which have larger radius for larger eigenenergies. This leads to chiral dynamics along reasonably well-defined orbits, both in the single-particle and in the two-particle sectors.

preprint2020arXiv

Gradient descent with momentum --- to accelerate or to super-accelerate?

We consider gradient descent with `momentum', a widely used method for loss function minimization in machine learning. This method is often used with `Nesterov acceleration', meaning that the gradient is evaluated not at the current position in parameter space, but at the estimated position after one step. In this work, we show that the algorithm can be improved by extending this `acceleration' --- by using the gradient at an estimated position several steps ahead rather than just one step ahead. How far one looks ahead in this `super-acceleration' algorithm is determined by a new hyperparameter. Considering a one-parameter quadratic loss function, the optimal value of the super-acceleration can be exactly calculated and analytically estimated. We show explicitly that super-accelerating the momentum algorithm is beneficial, not only for this idealized problem, but also for several synthetic loss landscapes and for the MNIST classification task with neural networks. Super-acceleration is also easy to incorporate into adaptive algorithms like RMSProp or Adam, and is shown to improve these algorithms.

preprint2020arXiv

Non-Hermitian scattering on a tight-binding lattice

We analyze the scattering dynamics and spectrum of a quantum particle on a tight-binding lattice subject to a non-Hermitian (purely imaginary) local potential. The reflection, transmission and absorption coefficients are studied as a function of the strength of this absorbing potential. The system is found to have an exceptional point at a certain strength of the potential. Unusually, all (or nearly all) of the spectrum pairs up into mutually coalescing eigenstate pairs at this exceptional point. At large potential strengths, the absorption coefficient decreases and the effect of the imaginary potential is similar to that of a real potential. We quantify this similarity by utilizing properties of a localized eigenstate.

preprint2020arXiv

Statistics of correlation functions in the random Heisenberg chain

Ergodic quantum many-body systems satisfy the eigenstate thermalization hypothesis (ETH). However, strong disorder can destroy ergodicity through many-body localization (MBL) -- at least in one dimensional systems -- leading to a clear signal of the MBL transition in the probability distributions of energy eigenstate expectation values of local operators. For a paradigmatic model of MBL, namely the random-field Heisenberg spin chain, we consider the full probability distribution of eigenstate correlation functions across the entire phase diagram. We find gaussian distributions at weak disorder, as predicted by pure ETH. At intermediate disorder -- in the thermal phase -- we find further evidence for anomalous thermalization in the form of heavy tails of the distributions. In the MBL phase, we observe peculiar features of the correlator distributions: a strong asymmetry in $S_i^z S_{i+r}^z$ correlators skewed towards negative values; and a multimodal distribution for spin-flip correlators. A quantitative quasi-degenerate perturbation theory calculation of these correlators yields a surprising agreement of the full distribution with the exact results, revealing, in particular, the origin of the multiple peaks in the spin-flip correlator distribution as arising from the resonant and off-resonant admixture of spin configurations. The distribution of the $S_i^zS_{i+r}^z$ correlator exhibits striking differences between the MBL and Anderson insulator cases.

preprint2019arXiv

Many-body quantum dynamics of initially trapped systems due to a Stark potential --- thermalization vs. Bloch oscillations

We analyze the dynamics of an initially trapped cloud of interacting quantum particles on a lattice under a linear (Stark) potential. We reveal a dichotomy: initially trapped interacting systems possess features typical of both many-body-localized and self-thermalizing systems. We consider both fermions ($t$-$V$ model) and bosons (Bose-Hubbard model). For the zero and infinite interaction limits, both systems are integrable: we provide analytic solutions in terms of the moments of the initial cloud shape, and clarify how the recurrent dynamics (many-body Bloch oscillations) depends on the initial state. Away from the integrable points, we identify and explain the time scale at which Bloch oscillations decohere.

preprint2016arXiv

Response of fermions in Chern bands to spatially local quenches

We study the dynamical evolution of Chern-band systems after subjecting them to local quenches. For open-boundary systems, we show for half-filling that the chiral nature of edge states is manifested in the time-dependent chiral response to local density quenches on the edge. In the presence of power-law traps, we show how to mimic the half-filling situation by choosing the appropriate number of fermions depending on the trap size, and explore chiral responses of edges to local quenches in such a configuration. We find that perturbations resulting from the quenches propagate at different group velocities depending on the bulk band gap. Our results provide different routes to check dynamically the non-trivial nature of Chern bands.

preprint2016arXiv

Wavepacket dynamics on Chern band lattices in a trap

The experimental realization of lattices with Chern bands in ultracold-atom and photonic systems has motivated the study of time-dependent phenomena, such as spatial propagation, in lattices with nontrivial topology. We study the dynamics of gaussian wavepackets on the Haldane honeycomb Chern-band lattice model, in the presence of a harmonic trap. We focus on the transverse response to a force, which is due partly to the Berry curvature and partly to the transverse component of the energy band curvature. We evaluate the accuracy of a semiclassical description, which treats the wavepacket as a point particle in both real and momentum space, in reproducing the motion of a realistic wavepacket with finite extent. We find that, in order to accurately capture the wavepacket dynamics, the extent of the wavepacket in momentum space needs to be taken into account. The dynamics is sensitive to the interplay of band dispersion and Berry curvature over the finite region of momentum (reciprocal) space where the wavepacket has support. Moreover, if the wavepacket is prepared with a finite initial momentum, the semiclassical analysis reproduces its motion as long as it has a large overlap with the eigenstates of a single band. The semiclassical description generally improves with increasing real-space size of the wavepacket, as long as the external conditions (e.g., external force) remain uniform throughout the spatial extent of the wavepacket.

preprint2015arXiv

Bose-Hubbard ladder subject to effective magnetic field: quench dynamics in a harmonic trap

Motivated by a recent experiment with optical lattices that has realized a ladder geometry with an effective magnetic field (Atala et al., Nature Physics \textbf{10}, 588 (2014)), we study the dynamics of bosons on a tight-binding two-leg ladder with complex hopping amplitudes. This system displays a quantum phase transition even without interactions. We study the non-equilibrium dynamics without and with interactions, in the presence of a harmonic trapping potential. In particular we consider dynamics induced by quenches of the trapping potential and of the magnitude of the rung hopping. We present a striking "slowing down" effect in the collective mode dynamics near the phase transition. This manifestation of a slowing down phenomenon near a quantum phase transition can be visualized unusually directly: the collective mode dynamics can be followed experimentally in real time and real space by imaging the atomic cloud.

preprint2015arXiv

Dynamics of the impurity screening cloud following quantum quenches of the Resonant Level Model

We present a detailed analysis of the evolution of impurity screening cloud in the Resonant Level Model (RLM) followed by quenches of system parameters. The considered quenches are either local or found to be effectively local. The screening cloud is characterized by impurity-bath correlators and by the entanglement of a block centered around the impurity with the rest of the system. We consider several local quench protocols, such as, quenches from initially decoupled to coupled impurity, quenches between different finite couplings of the impurity, and `detuning quenches' involving changing the onsite impurity potential away from or to the chemical potential of the bath. The relevant correlators and the entanglement display versions of the `light-cone' effect, since the information about the quench travels through the bath at finite speed. At long times (`inside' the light cone), the impurity screening cloud relaxes exponentially to the final equilibrium structure, with the relaxation rate given by the emergent energy scale of impurity screening. Also, snaphots of the time-evolving spatial profile of impurity-bath correlators show exponential dependences on space, with the length scale given by the screening length.

preprint2015arXiv

Inhomogeneous quenches in a fermionic chain: exact results

We consider the non-equilibrium physics induced by joining together two tight binding fermionic chains to form a single chain. Before being joined, each chain is in a many-fermion ground state. The fillings (densities) in the two chains might be the same or different. We present a number of exact results for the correlation functions in the non-interacting case. We present a short-time expansion, which can sometimes be fully resummed, and which reproduces the so-called `light cone' effect or wavefront behavior of the correlators. For large times, we show how all interesting physical regimes may be obtained by stationary phase approximation techniques. In particular, we derive semiclassical formulas in the case when both time and positions are large, and show that these are exact in the thermodynamic limit. We present subleading corrections to the large-time behavior, including the corrections near the edges of the wavefront. We also provide results for the return probability or Loschmidt echo. In the maximally inhomogeneous limit, we prove that it is exactly gaussian at all times. The effects of interactions on the Loschmidt echo are also discussed.

preprint2015arXiv

Nonsmooth and level-resolved dynamics illustrated with the tight binding model

We point out that in the first order time-dependent perturbation theory, the transition probability may behave nonsmoothly in time and have kinks periodically. Moreover, the detailed temporal evolution can be sensitive to the exact locations of the eigenvalues in the continuum spectrum, in contrast to coarse-graining ideas. Underlying this nonsmooth and level-resolved dynamics is a simple equality about the sinc function $\sinc x \equiv \sin x / x$. These physical effects appear in many systems with approximately equally spaced spectra, and is also robust for larger-amplitude coupling beyond the domain of perturbation theory. We use a one-dimensional periodically driven tight-binding model to illustrate these effects, both within and outside the perturbative regime.

preprint2015arXiv

Off-diagonal matrix elements of local operators in many-body quantum systems

In the time evolution of isolated quantum systems out of equilibrium, local observables generally relax to a long-time asymptotic value, governed by the expectation values (diagonal matrix elements) of the corresponding operator in the eigenstates of the system. The temporal fluctuations around this value, response to further perturbations, and the relaxation toward this asymptotic value, are all determined by the off-diagonal matrix elements. Motivated by this non-equilibrium role, we present generic statistical properties of off-diagonal matrix elements of local observables in two families of interacting many-body systems with local interactions. Since integrability (or lack thereof) is an important ingredient in the relaxation process, we analyze models that can be continuously tuned to integrability. We show that, for generic non-integrable systems, the distribution of off-diagonal matrix elements is a gaussian centered at zero. As one approaches integrability, the peak around zero becomes sharper, so that the distribution is approximately a combination of two gaussians. We characterize the proximity to integrability through the deviation of this distribution from a gaussian shape. We also determine the scaling dependence on system size of the average magnitude of off-diagonal matrix elements.

preprint2015arXiv

Overlap distributions for quantum quenches in the anisotropic Heisenberg chain

The dynamics after a quantum quench is determined by the weights of the initial state in the eigenspectrum of the final Hamiltonian, i.e., by the distribution of overlaps in the energy spectrum. We present an analysis of such overlap distributions for quenches of the anisotropy parameter in the one-dimensional anisotropic spin-1/2 Heisenberg model (XXZ chain). We provide an overview of the form of the overlap distribution for quenches from various initial anisotropies to various final ones, using numerical exact diagonalization. We show that if the system is prepared in the antiferromagnetic Néel state (infinite anisotropy) and released into a non-interacting setup (zero anisotropy, XX point) only a small fraction of the final eigenstates gives contributions to the post-quench dynamics, and that these eigenstates have identical overlap magnitudes. We derive expressions for the overlaps, and present the selection rules that determine the final eigenstates having nonzero overlap. We use these results to derive concise expressions for time-dependent quantities (Loschmidt echo, longitudinal and transverse correlators) after the quench. We use perturbative analyses to understand the overlap distribution for quenches from infinite to small nonzero anisotropies, and for quenches from large to zero anisotropy.

preprint2015arXiv

Quantum quench in 2D using the variational Baeriswyl wavefunction

By combining the Baeriswyl wavefunction with equilibrium and time-dependent variational principles, we develop a non-equilibrium formalism to study quantum quenches for two dimensional spinless fermions with nearest-neighbour hopping and repulsion. The variational ground state energy and the short time dynamics agree convincingly with the results of numerically exact simulations. We find that depending on the initial and final interaction strength, the quenched system either exhibits undamped oscillations or relaxes to a time independent steady state. The time averaged expectation value of the CDW order parameter rises sharply when crossing from the steady state regime to the oscillating regime, indicating that the system, being non-integrable, shows signs of thermalization with an effective temperature above or below the equilibrium critical temperature, respectively.

preprint2014arXiv

Entanglement Spectrum and Entanglement Hamiltonian of a Chern insulator with open boundaries

We study the entanglement spectrum of a Chern insulator on a cylinder geometry, with the cut separating the two partitions taken parallel to the cylinder edge, at varying distances from the edge. In contrast to similar studies on a torus, there is only one cut, and hence only one virtual edge mode in the entanglement spectrum. The entanglement spectrum has a gap when the cut is close enough to the physical edge of the cylinder such that the edge mode spatially extends over the cut. This effect is suppressed for parameter choices where the edge mode is sharply localized at the edge. In the extreme case of a perfectly localized edge mode, the entanglement spectrum is gapless even if the smaller partition consists of a single edge row. For the single-row cut, we construct the corresponding entanglement Hamiltonian, which is a one-dimensional tight-binding Hamiltonian with complex long-range hopping and interesting properties. We also study and explain the effect of two different schemes of flux insertion through a ring described by such an entanglement Hamiltonian.

preprint2014arXiv

Finite-size scaling of eigenstate thermalization

According to the eigenstate thermalization hypothesis (ETH), even isolated quantum systems can thermalize because the eigenstate-to-eigenstate fluctuations of typical observables vanish in the limit of large systems. Of course, isolated systems are by nature finite, and the main way of computing such quantities is through numerical evaluation for finite-size systems. Therefore, the finite-size scaling of the fluctuations of eigenstate expectation values is a central aspect of the ETH. In this work, we present numerical evidence that for generic non-integrable systems these fluctuations scale with a universal power law $D^{-1/2}$ with the dimension $D$ of the Hilbert space. We provide heuristic arguments, in the same spirit as the ETH, to explain this universal result. Our results are based on the analysis of three families of models, and several observables for each model. Each family includes integrable members, and we show how the system size where the universal power law becomes visible is affected by the proximity to integrability.

preprint2014arXiv

Global characteristics of all eigenstates of local many-body Hamiltonians: participation ratio and entanglement entropy

In the spectrum of many-body quantum systems, the low-energy eigenstates were the traditional focus of research. The interest in the statistical properties of the full eigenspectrum has grown more recently, in particular in the context of non-equilibrium questions. Wave functions of interacting lattice quantum systems can be characterized either by local observables, or by global properties such as the participation ratio (PR) in a many-body basis or the entanglement between various partitions. We present a study of the PR and of the entanglement entropy (EE) between two roughly equal spatial partitions of the system, in all the eigenfunctions of local Hamiltonians. Motivated by the similarity of the PR and EE - both are generically larger in the bulk and smaller near the edges of the spectrum - we quantitatively analyze the correlation between them. We elucidate the effect of (proximity to) integrability, showing how low-entanglement and low-PR states appear also in the middle of the spectrum as one approaches integrable points. We also determine the precise scaling behavior of the eigenstate-to-eigenstate fluctuations of the PR and EE with respect to system size, and characterize the statistical distribution of these quantities near the middle of the spectrum.

preprint2014arXiv

Modulated trapping of interacting bosons in one dimension

We investigate the response of harmonically confined bosons with contact interactions (trapped Lieb-Liniger gas) to modulations of the trapping strength. We explain the structure of resonances at a series of driving frequencies, where size oscillations and energy grow exponentially. For strong interactions (Tonks-Girardeau gas), we show the effect of resonant driving on the bosonic momentum distribution. The treatment is `exact' for zero and infinite interactions, where the dynamics is captured by a single-variable ordinary differential equation. For finite interactions the system is no longer exactly solvable. For weak interactions, we show how interactions modify the resonant behavior for weak and strong driving, using a variational approximation which adds interactions to the single-variable description in a controlled way.

preprint2014arXiv

Real-space structure of the impurity screening cloud in the Resonant Level Model

We present a detailed investigation of the impurity screening cloud in the Resonant Level Model. The screening is visible in the structure of impurity-bath correlators as a function of distance from the impurity. We characterize the screening cloud through scaling analyses of impurity-bath correlators and entanglement entropies. We devise and study several situations where the screening cloud is destroyed or modified: finite temperatures, energetic detuning of the impurity from the chemical potential, and the situation of an unconventional bath with diverging density of states.

preprint2014arXiv

Repulsive to attractive interaction quenches of 1D Bose gas in a harmonic trap

We consider quantum quenches of harmonically trapped one-dimensional bosons from repulsive to attractive interactions, and the resulting breathing dynamics of the so-called `super-Tonks-Girardeau' (sTG) state. This state is highly excited compared to the ground state of the attractive gas, and is the lowest eigenstate where the particles are not bound or clustered. We analyze the dynamics from a spectral point of view, identifying the relevant eigenstates of the interacting trapped many-body system, and analyzing the nature of these quantum eigenstates. To obtain explicit eigenspectra, we use Hamiltonians with finite-dimensional Hilbert spaces to approximate the Lieb-Liniger system. We employ two very different approximate approaches: an expansion in a truncated single-particle harmonic-trap basis and a lattice (Bose-Hubbard) model. We show how the breathing frequency, identified with the energy difference between the sTG state and another particular eigenstate, varies with interaction.

preprint2013arXiv

Bethe ansatz description of edge-localization in the open-boundary XXZ spin chain

At large values of the anisotropy Δ, the open-boundary Heisenberg spin-1/2 chain has eigenstates displaying localization at the edges. We present a Bethe ansatz description of this `edge-locking' phenomenon in the entire Δ>1 region. We focus on the simplest spin sectors, namely the highly polarized sectors with only one or two overturned spins, i.e., one-particle and two-particle sectors. Edge-locking is associated with pure imaginary solutions of the Bethe equations, which are not commonly encountered in periodic chains. In the one-particle case, at large anisotropies there are two eigenstates with imaginary Bethe momenta, related to localization at the two edges. For any finite chain size, one of the two solutions become real as the anisotropy is lowered below a certain value. For two particles, a richer scenario is observed, with eigenstates having the possibility of both particles locked on the same or different edge, one locked and the other free, and both free either as single magnons or as bound composites corresponding to `string' solutions. For finite chains, some of the edge-locked spins get delocalized at certain values of the anisotropy (`exceptional points'), corresponding to imaginary solutions becoming real. We characterize these phenomena thoroughly by providing analytic expansions of the Bethe momenta for large chains, large anisotropy, and near the exceptional points. In the large-chain limit all the exceptional points coalesce at the isotropic point (Δ=1) and edge-locking becomes stable in the whole Δ>1 region.

preprint2013arXiv

Breathing mode in the Bose-Hubbard chain with a harmonic trapping potential

We investigate the breathing mode of harmonically trapped bosons in an optical lattice at small site occupancies. The Bose-Hubbard model with a trapping potential is used to describe the breathing-mode dynamics initiated through weak quenches of the trap strength. We connect to results for continuum bosons (Lieb-Liniger and Gross-Pitaevskii results) and also present deviations from continuum physics. We take a spectral perspective, identifying the breathing mode frequency with a particular energy gap in the spectrum of the trapped Bose-Hubbard Hamiltonian. We present the low energy eigenspectrum of the trapped many-boson system, and study overlaps of the initial state with eigenstates of the quenched Hamiltonian. There is an intermediate interaction regime, between a "free-boson" limit and a "free-fermion" limit, in which the Bose-Hubbard breathing mode frequency approaches the Gross-Pitaevskii prediction. In addition, we present a striking failure of the time-dependent Gutzwiller approximation for describing breathing modes.

preprint2013arXiv

Entanglement spectrum of the two dimensional Bose-Hubbard model

We study the entanglement spectrum (ES) of the Bose-Hubbard model on the two dimensional square lattice at unit filling, both in the Mott insulating and in the superfluid phase. In the Mott phase, we demonstrate that the ES is dominated by the physics at the boundary between the two subsystems. On top of the boundary-local (perturbative) structure, the ES exhibits substructures arising from one-dimensional dispersions along the boundary. In the superfluid phase, the structure of the ES is qualitatively different, and reflects the spontaneously broken U(1) symmetry of the phase. We attribute the basic low-lying structure to a so-called "tower of states" (TOS) Hamiltonian of the model. We then discuss how these characteristic structures evolve across the superfluid to Mott insulator transition and their influence on the behavior of the entanglement entropies. Finally, we briefly outline the implications of the ES structure on the efficiency of matrix-product-state based algorithms in two dimensions.

preprint2013arXiv

Fine structures in the spectrum of the open-boundary Heisenberg chain at large anisotropies

At large anisotropies, the spectrum of the Heisenberg XXZ spin chain separates into `bands' with energies largely determined by the number of domain walls. The band structure is richer with open boundary conditions: there are more bands and the bands develop intricate fine structures. We characterize and explain these structures and substructures in the open-boundary chain. The fine structures are explained using degenerate perturbation theory. We also present some dynamical consequences of these sub-band structures, through explicit time evolution of the wavefunction from initial states motivated by the fine structure analysis.

preprint2013arXiv

Quantum quenches and work distributions in ultra-low-density systems

We present results on quantum quenches in systems with a fixed number of particles in a large region. We show that the typical differences between local and global quenches present in systems with regular thermodynamic limit are lacking in this low-density limit. In particular, we show that in this limit local quenches may not lead to equilibration to the new ground state, and that global quenches can have power-law work distributions ("edge singularities") typically associated with local quenches for finite-density systems. We also show that this regime allows for large edge singularity exponents beyond that allowed by the constraints of the usual thermodynamic limit. This large-exponent singularity has observable consequences in the time evolution, leading to a distinct intermediate power-law regime in time. We demonstrate these results first using local quantum quenches in a low-density Kondo-like system, and additionally through global and local quenches in Bose-Hubbard, Aubry-Andre, and hard-core boson systems in the low-density regime.

preprint2013arXiv

Single-impurity Kondo physics at extreme particle-hole asymmetry

We study the fate of the Kondo effect with one-dimensional conduction baths at very low densities, such that the system explores the bottom of the conduction band. This can involve either finite low densities, or a small number of fixed conduction electrons in a large system, i.e., the limit of large bath sizes can be taken with either fixed small density or with fixed number. We characterize the Kondo physics for such systems through the energy gain due to Kondo coupling, which is a general analog of the Kondo temperature scale, and through real-space profiles of densities and spin-spin correlation functions.

preprint2013arXiv

Slow interaction ramps in trapped many-particle systems: universal deviations from adiabaticity

For harmonic-trapped atomic systems, we report system-independent non-adiabatic features in the response to interaction ramps. We provide results for several different systems in one, two, and three dimensions: bosonic and fermionic Hubbard models realized through optical lattices, a Bose-Einstein condensate, a fermionic superfluid and a fermi liquid. The deviation from adiabaticity is characterized through the heating or excitation energy produced during the ramp. We find that the dependence of the heat on the ramp time is sensitive to the ramp protocol but has aspects common to all systems considered. We explain these common features in terms of universal dynamics of the system size or cloud radius.

preprint2013arXiv

Viewpoint: Toward Fractional Quantum Hall physics with cold atoms

Viewpoint on Nigel R. Cooper and Jean Dalibard, "Reaching Fractional Quantum Hall States with Optical Flux Lattices", Phys. Rev. Lett. 110, 185301 (2013), and N. Y. Yao, A. V. Gorshkov, C. R. Laumann, A. M. Läuchli, J. Ye, and M. D. Lukin, "Realizing Fractional Chern Insulators in Dipolar Spin Systems", Phys. Rev. Lett. 110, 185302 (2013). Researchers propose new ways to recreate fractional quantum Hall physics using ultracold atoms and molecules.

preprint2012arXiv

Boundary-locality and perturbative structure of entanglement spectra in gapped systems

The entanglement between two parts of a many-body system can be characterized in detail by the entanglement spectrum. Focusing on gapped phases of one-dimensional systems, we show how this spectrum is dominated by contributions from the boundary between the parts. This contradicts the view of an "entanglement Hamiltonian" as a bulk entity. The boundary-local nature of the entanglement spectrum is clarified through its hierarchical level structure, through the combination of two single-boundary spectra to form a two-boundary spectrum, and finally through consideration of dominant eigenfunctions of the entanglement Hamiltonian. We use the boundary-locality to formulate a perturbative scheme for calculating entanglement spectra.

preprint2012arXiv

Entanglement spectrum of the Heisenberg XXZ chain near the ferromagnetic point

We study the entanglement spectrum (ES) of a finite XXZ spin 1/2 chain in the limit Δ-> -1^+ for both open and periodic boundary conditions. At Δ=-1 (ferromagnetic point) the model is equivalent to the Heisenberg ferromagnet and its degenerate ground state manifold is the SU(2) multiplet with maximal total spin. Any state in this so-called "symmetric sector" is an equal weight superposition of all possible spin configurations. In the gapless phase at Δ>-1 this property is progressively lost as one moves away from the Δ=-1 point. We investigate how the ES obtained from the states in this manifold reflects this change, using exact diagonalization and Bethe ansatz calculations. We find that in the limit Δ->-1^+ most of the ES levels show divergent behavior. Moreover, while at Δ=-1 the ES contains no information about the boundaries, for Δ>-1 it depends dramatically on the choice of boundary conditions. For both open and periodic boundary conditions the ES exhibits an elegant multiplicity structure for which we conjecture a combinatorial formula. We also study the entanglement eigenfunctions, i.e. the eigenfunctions of the reduced density matrix. We find that the eigenfunctions corresponding to the non diverging levels mimic the behavior of the state wavefunction, whereas the others show intriguing polynomial structures. Finally we analyze the distribution of the ES levels as the system is detuned away from Δ=-1.

preprint2012arXiv

Itinerant electrons in the Coulomb phase

We study the interplay between magnetic frustration and itinerant electrons. For example, how does the coupling to mobile charges modify the properties of a spin liquid, and does the underlying frustration favor insulating or conducting states? Supported by Monte Carlo simulations, our goal is in particular to provide an analytical picture of the mechanisms involved. The models under considerations exhibit Coulomb phases in two and three dimensions, where the itinerant electrons are coupled to the localized spins via double exchange interactions. Because of the Hund coupling, magnetic loops naturally emerge from the Coulomb phase and serve as conducting channels for the mobile electrons, leading to doping-dependent rearrangements of the loop ensemble in order to minimize the electronic kinetic energy. At low electron density ρ, the double exchange coupling mainly tends to segment the very long loops winding around the system into smaller ones while it gradually lifts the extensive degeneracy of the Coulomb phase with increasing ρ. For higher doping, the results are strongly lattice dependent, displaying loop crystals with a given loop length for some specific values of ρ, which can melt into another loop crystal by varying ρ. Finally, we contrast this to the qualitatively different behavior of analogous models on kagome or triangular lattices.

preprint2012arXiv

Linear quantum quench in the Heisenberg XXZ chain: time dependent Luttinger model description of a lattice system

We study variable-rate linear quenches in the anisotropic Heisenberg (XXZ) chain, starting at the XX point. This is equivalent to switching on a nearest neighbour interaction for hard-core bosons or an interaction quench for free fermions. The physical observables we investigate are: the energy pumped into the system during the quench, the spin-flip correlation function, and the bipartite fluctuations of the z component of the spin in a box. We find excellent agreement between exact numerics (infinite system time-evolving block decimation, iTEBD) and analytical results from bosonization, as a function of the quench time, spatial coordinate and interaction strength. This provides a stringent and much-needed test of Luttinger liquid theory in a non-equilibrium situation.

preprint2012arXiv

Non-equilibrium dynamics in Bose-Hubbard ladders

Motivated by a recent experiment on the non-equilibrium dynamics of interacting bosons in ladder-shaped optical lattices, we report exact calculations on the sweep dynamics of Bose-Hubbard systems in finite two-leg ladders. The sweep changes the energy bias between the legs linearly over a finite time. As in the experiment, we study the cases of [a] the bosons initially all in the lower-energy leg (ground state sweep) and [b] the bosons initially all in the higher-energy leg (inverse sweep). The approach to adiabaticity in the inverse sweep is intricate, as the transfer of bosons is non-monotonic as a function of both sweep time and intra-leg tunnel coupling. Our exact study provides explanations for these non-monotonicities based on features of the full spectrum, without appealing to concepts (e.g., gapless excitation spectrum) that are more appropriate for the thermodynamic limit. We also demonstrate and study Stueckelberg oscillations in the finite-size ladders.

preprint2012arXiv

Spin modulation instabilities and phase separation dynamics in trapped two-component Bose condensates

In the study of trapped two-component Bose gases, a widely used dynamical protocol is to start from the ground state of a one-component condensate and then switch half the atoms into another hyperfine state. The slightly different intra-component and inter-component interactions can then lead to highly nontrivial dynamics. We study and classify the possible subsequent dynamics, over a wide variety of parameters spanned by the trap strength and by the inter- to intra-component interaction ratio. A stability analysis suited to the trapped situation provides us with a framework to explain the various types of dynamics in different regimes.

preprint2012arXiv

Strongly interacting bosons in multi-chromatic potentials supporting mobility edges: localization, quasi-condensation and expansion dynamics

We provide an account of the static and dynamic properties of hard-core bosons in a one-dimensional lattice subject to a multi-chromatic quasiperiodic potential for which the single-particle spectrum has mobility edges. We use the mapping from strongly interacting bosons to weakly interacting fermions, and provide exact numerical results for hard-core bosons in and out of equilibrium. In equilibrium, we find that the system behaves like a quasi-condensate (insulator) depending on whether the Fermi surface of the corresponding fermionic system lies in a spectral region where the single-particle states are delocalized (localized). We also study non-equilibrium expansion dynamics of initially trapped bosons, and demonstrate that the extent of partial localization is determined by the single-particle spectrum.

preprint2012arXiv

Strongly interacting one-dimensional bosons in arbitrary-strength optical lattices: from Bose-Hubbard to sine-Gordon and beyond

We analyze interacting one-dimensional bosons in the continuum, subject to a periodic sinusoidal potential of arbitrary depth. Variation of the lattice depth tunes the system from the Bose-Hubbard limit for deep lattices, through the sine-Gordon regime of weak lattices, to the complete absence of a lattice. Using the Bose-Fermi mapping between strongly interacting bosons and weakly interacting fermions, we derive the phase diagram in the parameter space of lattice depth and chemical potential. This extends previous knowledge from tight-binding (Bose-Hubbard) studies in a new direction which is important because the lattice depth is a readily adjustable experimental parameter. Several other results (equations of state, energy gaps, profiles in harmonic trap) are presented as corollaries to the physics contained in this phase diagram. Generically, both incompressible (gapped) and compressible phases coexist in a trap; this has implications for experimental measurements.

preprint2011arXiv

Crossover from adiabatic to sudden interaction quench in a Luttinger liquid

Motivated by recent experiments on interacting cold atoms, we analyze interaction quenches in Luttinger liquids (LL), where the interaction is ramped from zero to a finite value within a finite time. The fermionic single particle density matrix reveals several regions of spatial and temporal coordinates relative to the quench time, termed as Fermi liquid, sudden quench LL, adiabatic LL regimes, and a LL regime with time dependent exponent. The various regimes can also be observed in the momentum distribution of the fermions, directly accessible through time of flight experiments. Most of our results apply to arbitrary quench protocols.

preprint2010arXiv

Disentangling Entanglement Spectra of Fractional Quantum Hall States on Torus Geometries

We analyze the entanglement spectrum of Laughlin states on the torus and show that it is arranged in towers, each of which is generated by modes of two spatially separated chiral edges. This structure is present for all torus circumferences, which allows for a microscopic identification of the prominent features of the spectrum by perturbing around the thin torus limit.

preprint2010arXiv

Entanglement between particle partitions in itinerant many-particle states

We review `particle partitioning entanglement' for itinerant many-particle systems. This is defined as the entanglement between two subsets of particles making up the system. We identify generic features and mechanisms of particle entanglement that are valid over whole classes of itinerant quantum systems. We formulate the general structure of particle entanglement in many-fermion ground states, analogous to the `area law' for the more usually studied entanglement between spatial regions. Basic properties of particle entanglement are first elucidated by considering relatively simple itinerant models. We then review particle-partitioning entanglement in quantum states with more intricate physics, such as anyonic models and quantum Hall states.

preprint2010arXiv

Entanglement Scaling of Fractional Quantum Hall states through Geometric Deformations

We present a new approach to obtaining the scaling behavior of the entanglement entropy in fractional quantum Hall states from finite-size wavefunctions. By employing the torus geometry and the fact that the torus aspect ratio can be readily varied, we can extract the entanglement entropy of a spatial block as a continuous function of the block boundary. This approach allows us to extract the topological entanglement entropy with an accuracy superior to what is possible on the spherical or disc geometry, where no natural continuously variable parameter is available. Other than the topological information, the study of entanglement scaling is also useful as an indicator of the difficulty posed by fractional quantum Hall states for various numerical techniques.

preprint2010arXiv

Finite-rate quenches of site bias in the Bose-Hubbard dimer

For a Bose-Hubbard dimer, we study quenches of the site energy imbalance, taking a highly asymmetric Hamiltonian to a fully symmetric one. The ramp is carried out over a finite time that interpolates between the instantaneous and adiabatic limits. We provide results for the excess energy of the final state compared to the ground state energy of the final Hamiltonian, as a function of the quench rate. This excess energy serves as the analog of the defect density that is considered in the Kibble-Zurek picture of ramps across phase transitions. We also examine the fate of quantum `self-trapping' when the ramp is not instantaneous.

preprint2010arXiv

Interaction induced fractional Bloch and tunneling oscillations

We study the dynamics of few interacting bosons in a one-dimensional lattice with dc bias. In the absence of interactions the system displays single particle Bloch oscillations. For strong interaction the Bloch oscillation regime reemerges with fractional Bloch periods which are inversely proportional to the number of bosons clustered into a bound state. The interaction strength is affecting the oscillation amplitude. Excellent agreement is found between numerical data and a composite particle dynamics approach. For specific values of the interaction strength a particle will tunnel from the interacting cloud to a well defined distant lattice location.

preprint2009arXiv

Edge-locking and quantum control in highly polarized spin chains

For an open-boundary spin chain with anisotropic Heisenberg (XXZ) interactions, we present states in which a connected block near the edge is polarized oppositely to the rest of the chain. We show that such blocks can be `locked' to the edge of the spin chain, and that there is a hierarchy of edge-locking effects at various orders of the anisotropy. The phenomenon enables dramatic control of quantum state transmission: the locked block can be freed by flipping a single spin or a few spins.

preprint2008arXiv

Survival benefits in mimicry: a quantitative framework

Mimicry is a resemblance between species that benefits at least one of the species. It is a ubiquitous evolutionary phenomenon particularly common among prey species, in which case the advantage involves better protection from predation. We formulate a mathematical description of mimicry among prey species, to investigate benefits and disadvantages of mimicry. The basic setup involves differential equations for quantities representing predator behavior, namely, the probabilities for attacking prey at the next encounter. Using this framework, we present new quantitative results, and also provide a unified description of a significant fraction of the quantitative mimicry literature. The new results include `temporary' mutualism between prey species, and an optimal density at which the survival benefit is greatest for the mimic. The formalism leads naturally to extensions in several directions, such as the evolution of mimicry, the interplay of mimicry with population dynamics, etc. We demonstrate this extensibility by presenting some explorations on spatiotemporal pattern dynamics.

preprint2007arXiv

Trapped fermionic clouds distorted from the trap shape due to many-body effects

We present a general approach for calculating densities and other local quantities of trapped Fermi gases, when the cloud shape is distorted with respect to the trap shape due to global energy considerations. Our approach provides a consistent way to explore physics beyond the local density approximation, if this is necessary due to the distortion. We illustrate this by analyzing in detail experimentally observed distortions in an imbalanced Fermi mixture in an elongated trap. In particular, we demonstrate in that case dramatic deviations from ellipsoidal cloud shapes arising from the competition between surface and bulk energies.

preprint2006arXiv

Pairing of a trapped resonantly-interacting fermion mixture with unequal spin populations

We consider the phase separation of a trapped atomic mixture of fermions with unequal spin populations near a Feshbach resonance. In particular, we determine the density profile of the two spin states and compare with the recent experiments of Partridge et al. (cond-mat/0511752). Overall we find quite good agreement. We identify the remaining discrepancies and pose them as open problems.

preprint2004arXiv

Sweeping a molecular Bose-Einstein condensate across a Feshbach resonance

We consider the dissociation of a molecular Bose-Einstein condensate during a magnetic-field sweep through a Feshbach resonance that starts on the molecular side of the resonance and ends on the atomic side. In particular, we determine the energy distribution of the atoms produced after the sweep. We find that the shape of the energy distribution strongly depends on the rate of the magnetic-field sweep, in a manner that is in good agreement with recent experiments.