Noninteracting Limit Research Articles

Localization, multifractality, and many-body localization in periodically kicked quasiperiodic lattices

We study the combined effect of quasiperiodic disorder, driven and interaction in the periodically kicked Aubry-Andr\'{e} model. In the non-interacting limit, by analyzing the quasienergy spectrum statistics, we verify the existence of a dynamical localization transition in the high-frequency region, whereas the spectrum statistics becomes intricate in the low-frequency region due to the emergence of the extended/localized-to-multifractal edges in the quasienergy spectrum, which separate the multifractal states from the extended (localized) states. When the interaction is introduced, we find the periodically kicked incommensurate potential can lead to a transition from ergodic to many-body-localization phase in the high-frequency region. However, the many-body localization phase vanishes in the low-frequency region even for strong quasiperiodic disorder. Our studies demonstrate that the periodically kicked Aubry-Andr\'{e} model displays rich dynamical phenomena and the driving frequency plays an important role in the formation of many-body localization in addition to the disorder strength.

Magnetic excitations, nonclassicality, and quantum wake spin dynamics in the Hubbard chain

Recent work has demonstrated that quantum Fisher information (QFI), a witness of multipartite entanglement, and magnetic Van Hove correlations $G(r,t)$, a probe of local real-space real-time spin dynamics, can be successfully extracted from inelastic neutron scattering on spin systems through accurate measurements of the dynamical spin structure factor $S(k,\ensuremath{\omega})$. Here we apply theoretically these ideas to the half-filled Hubbard chain with nearest-neighbor hopping, away from the strong-coupling limit. This model has nontrivial redistribution of spectral weight in $S(k,\ensuremath{\omega})$ going from the noninteracting limit ($U=0$) to strong coupling ($U\ensuremath{\rightarrow}\ensuremath{\infty}$), where it reduces to the Heisenberg quantum spin chain. We use the density matrix renormalization group to find $S(k,\ensuremath{\omega})$, from which QFI is then calculated. We find that QFI grows with $U$. With realistic energy resolution it becomes capable of witnessing bipartite entanglement above $U=2.5$ (in units of the hopping), where it also changes slope. This point is also proximate to slope changes of the bandwidth $W(U)$ and the half-chain von Neumann entanglement entropy. We compute $G(r,t)$ by Fourier transforming $S(k,\ensuremath{\omega})$. The results indicate a crossover in the short-time short-distance dynamics at low $U$ characterized by ferromagnetic light-cone wavefronts, to a Heisenberg-type behavior at large $U$ featuring antiferromagnetic light cones and spatially period-doubled antiferromagnetism. We find this crossover has largely been completed by $U=3$. Our results thus provide evidence that, in several aspects, the strong-coupling limit of the Hubbard chain is reached qualitatively already at a relatively modest interaction strength. We discuss experimental candidates for observing the $G(r,t)$ dynamics found at low $U$.

Temporal Entanglement, Quasiparticles, and the Role of Interactions

In quantum many-body dynamics admitting a description in terms of noninteracting quasiparticles, the Feynman-Vernon influence matrix (IM), encoding the effect of the system on the evolution of its local subsystems, can be analyzed exactly. For discrete dynamics, the temporal entanglement (TE) of the corresponding IM satisfies an area law, suggesting the possibility of an efficient representation of the IM in terms of matrix-product states. A natural question is whether integrable interactions, preserving stable quasiparticles, affect the behavior of the TE. While a simple semiclassical picture suggests a sublinear growth in time, one can wonder whether interactions may lead to violations of the area law. We address this problem by analyzing quantum quenches in a family of discrete integrable dynamics corresponding to the real-time Trotterization of the interacting XXZ Heisenberg model. By means of an analytical solution at the dual-unitary point and numerical calculations for generic values of the system parameters, we provide evidence that, away from the noninteracting limit, the TE displays a logarithmic growth in time, thus violating the area law. Our findings highlight the nontrivial role of interactions, and raise interesting questions on the possibility to efficiently simulate the local dynamics of interacting integrable systems.

Competition of Quantum Anomalous Hall States and Charge Density Wave in a Correlated Topological Model

We investigate the topological phase transition driven by non-local electronic correlations in a realistic quantum anomalous Hall model consisting of dxy –d x 2 – y 2 orbitals. Three topologically distinct phases defined in the non-interacting limit evolve to different charge density wave phases under correlations. Two conspicuous conclusions were obtained: The topological phase transition does not involve gap-closing and the dynamical fluctuations significantly suppress the charge order favored by the next nearest neighbor interaction. Our study sheds light on the stability of topological phase under electronic correlations, and we demonstrate a positive role played by dynamical fluctuations that is distinct to all previous studies on correlated topological states.

Driven-Dissipative Criticality within the Discrete Truncated Wigner Approximation.

We present an approach to the numerical simulation of open quantum many-body systems based on the semiclassical framework of the discrete truncated Wigner approximation. We establish a quantum jump formalism to integrate the quantum master equation describing the dynamics of the system, which we find to be exact in both the noninteracting limit and the limit where the system is described by classical rate equations. We apply our method to simulation of the paradigmatic dissipative Ising model, where we are able to capture the critical fluctuations of the system beyond the level of mean-field theory.

Many body density of states in the edge of the spectrum: non-interacting limit

In noninteracting limit, the density of states (dos) of a many body system can be expressed as a convolution of the single body dos of its subunits. We use the formulation to derive, in the edge of the spectrum, a differential equation for the ensemble averaged many body dos that is relatively easier to solve. Our analysis, based on the systems in which the subunits can be modelled by a Gaussian or Wishart random matrix ensemble, indicates that a rescaling of energy by the number of subunits leaves the many body dos in a mathematically invariant form.

Deformation of localized states and state transitions in systems of randomly hopping interacting fermions

We numerically study the random-hopping fermions (the Cruetz ladder) with repulsion and investigate how the interactions deform localized eigenstates by means of the one particle-density matrix (OPDM). The ground state exhibits resurgence of localization from the compact localized state to strong-repulsion-induced localization. On the other hand, excited states in the middle of the spectrum tend to extend by the repulsion. The transition property obtained by numerical calculations of the OPDM is deeply understood by studying a solvable model in which local integrals of motion (LIOMs) are obtained explicitly. The present work clarifies the utility of the OPDM and also how compact-support LIOMs in non-interacting limit are deformed by the repulsion.

Hubbard model for spin-1 Haldane chains

The Haldane phase for antiferromagnetic spin-1 chains is a celebrated topological state of matter, featuring gapped excitations and fractional spin-1/2 edge states. Here, we provide numerical evidence that this phase can be realized with a Hubbard model at half filling, where each $s=1$ spin is stored in a four-site fermionic structure. We find that the noninteracting limit of our proposed model describes a one-dimensional topological insulator, and we conjecture it to be adiabatically connected to the Haldane phase. Our work suggests a route to build spin-1 networks using Hubbard model quantum simulators.

Fermion dispersion renormalization in a two-dimensional semi-Dirac semimetal

We present a non-perturbative study of the quantum many-body effects caused by the long-range Coulomb interaction in a two-dimensional semi-Dirac semimetal. This kind of semimetal may be realized in deformed graphene and a class of other realistic materials. In the non-interacting limit, the dispersion of semi-Dirac fermion is linear in one direction and quadratic in the other direction. When the impact of Coulomb interaction is taken into account, such a dispersion can be significantly modified. To reveal the correlation effects, we first obtain the exact self-consistent Dyson-Schwinger equation of the full fermion propagator and then extract the momentum dependence of the renormalized fermion dispersion from the numerical solutions. Our results show that the fermion dispersion becomes linear in two directions. These results are compared to previous theoretical works on semi-Dirac semimetals.

Dynamics of Strongly Interacting Fermi Gases with Time-Dependent Interactions: Consequence of Conformal Symmetry.

In this Letter, we investigate the effects of a time-dependent, short-ranged interaction on the long-time expansion dynamics of Fermi gases. We show that the effects of the interaction on the dynamics is dictated by how it changes under a conformal transformation, and derive an explicit criterion for the relevancy of time-dependent interactions near both the strongly and noninteracting scale invariant limits. In addition, we show that it is possible to engineer interactions that give rise to nonexponential thermalization dynamics in trapped Fermi gases. To supplement the symmetry analysis, we perform hydrodynamic simulations to show that the moment of inertia of the trapped gas indeed follows a universal time dependence that is determined jointly by the conformal symmetry and time-dependent scattering length a(t). Our results should also be relevant to the dynamics of other systems that are nearly scale invariant and that are governed by a nonrelativistic conformal symmetry.

Localization crossover and subdiffusive transport in a classical facilitated network model of a disordered interacting quantum spin chain

We consider the random-field Heisenberg model, a paradigmatic model for many-body localization (MBL), and add a Markovian dephasing bath coupled to the Anderson orbitals of the model's non-interacting limit. We map this system to a classical facilitated hopping model that is computationally tractable for large system sizes, and investigate its dynamics. The classical model exhibits a robust crossover between an ergodic (thermal) phase and a frozen (localized) phase. The frozen phase is destabilized by thermal subregions (bubbles), which thermalize surrounding sites by providing a fluctuating interaction energy and so enable off-resonance particle transport. Investigating steady state transport, we observe that the interplay between thermal and frozen bubbles leads to a clear transition between diffusive and subdiffusive regimes. This phenomenology both describes the MBL system coupled to a bath, and provides a classical analogue for the many-body localization transition in the corresponding quantum model, in that the classical model displays long local memory times. It also highlights the importance of the details of the bath coupling in studies of MBL systems coupled to thermal environments.

Phenomenology of Spectral Functions in Disordered Spin Chains at Infinite Temperature.

Studies of disordered spin chains have recently experienced a renewed interest, inspired by the question to which extent the exact numerical calculations comply with the existence of a many-body localization phase transition. For the paradigmatic random field Heisenberg spin chains, many intriguing features were observed when the disorder is considerable compared to the spin interaction strength. Here, we introduce a phenomenological theory that may explain some of those features. The theory is based on the proximity to the noninteracting limit, in which the system is an Anderson insulator. Taking the spin imbalance as an exemplary observable, we demonstrate that the proximity to the local integrals of motion of the Anderson insulator determines the dynamics of the observable at infinite temperature. In finite interacting systems our theory quantitatively describes its integrated spectral function for a wide range of disorders.

Floquet engineering and nonequilibrium topological maps in twisted trilayer graphene

Motivated by the recent experimental realization of twisted trilayer graphene and the observed superconductivity that is associated with its flat bands at specific angles, we study trilayer graphene under the influence of different forms of light in the noninteracting limit. Specifically, we study four different types of stacking configurations with a single twisted layer. In all four cases, we study the impact of circularly polarized light and longitudinal light coming from a waveguide. We derive effective time-independent Floquet Hamiltonians and review light-induced changes to the band structure. For circularly polarized light, we find band flattening effects as well as band gap openings. We emphasize that there is a rich band topology, which we summarize in Chern number maps that are different for all four studied lattice configurations. The case of a so-called ABC stacking with top layer twist is especially rich and shows a different phase diagram depending on the handedness of the circularly polarized light. Consequently, we propose an experiment where this difference in typologies could be captured via optical conductivity measurements. In contrast for the case of longitudinal light that is coming from a waveguide, we find that the band structure is very closely related to the equilibrium one but the magic angles can be tuned in situ by varying the intensity of the incident beam of light.

Connecting topological Anderson and Mott insulators in disordered interacting fermionic systems

The topological Anderson and Mott insulators are two phases that have so far been separately and widely explored beyond topological band insulators. Here we combine the two seemingly different topological phases into a system of spin-1/2 interacting fermionic atoms in a disordered optical lattice. We find that the topological Anderson and Mott insulators in the noninteracting and clean limits can be adiabatically connected without gap closing in the phase diagram of our model. Lying between the two phases, we uncover a disordered correlated topological insulator, which is induced from a trivial band insulator by the combination of disorder and interaction, as the generalization of topological Anderson insulators to the many-body interacting regime. The phase diagram is determined by computing various topological properties and confirmed by unsupervised and automated machine learning. We develop an approach to provide a unified and clear description of topological phase transitions driven by interaction and disorder. The topological phases can be detected from disorder/interaction induced edge excitations and charge pumping in optical lattices.

Flat-band ferromagnetism and spin waves in the Haldane-Hubbard model

We study the flat-band ferromagnetic phase of the Haldane-Hubbard model on a honeycomb lattice within a bosonization scheme for flat-band Chern insulators, focusing on the calculation of the spin-wave excitation spectrum. We consider the Haldane-Hubbard model with the noninteracting lower bands in a nearly-flat band limit, previously determined for the spinless model, and at 1/4-filling of its corresponding noninteracting limit. Within the bosonization scheme, the Haldane-Hubbard model is mapped into an effective interacting boson model, whose quadratic term allows us to determine the spin-wave spectrum at the harmonic approximation. We show that the excitation spectrum has two branches with a Goldstone mode and Dirac points at center and at the K and K' points of the first Brillouin zone, respectively. We also consider the effects on the spin-wave spectrum due to an energy offset in the on-site Hubbard repulsion energies and due to the presence of an staggered on-site energy term, both quantities associated with the two triangular sublattices. In both cases, we find that an energy gap opens at the K and K' points. Moreover, we also find some evidences for an instability of the flat-band ferromagnetic phase in the presence of the staggered on-site energy term. We provide some additional results for the square lattice topological Hubbard model previous studied within the bosonization formalism and comment on the differences between the bosonization scheme implementation for the correlated Chern insulators on both square and honeycomb lattices.

Higher-order topological Mott insulator on the pyrochlore lattice

We provide the first unbiased evidence for a higher-order topological Mott insulator in three dimensions by numerically exact quantum Monte Carlo simulations. This insulating phase is adiabatically connected to a third-order topological insulator in the noninteracting limit, which features gapless modes around the corners of the pyrochlore lattice and is characterized by a {mathbb {Z}}_{4} spin-Berry phase. The difference between the correlated and non-correlated topological phases is that in the former phase the gapless corner modes emerge only in spin excitations being Mott-like. We also show that the topological phase transition from the third-order topological Mott insulator to the usual Mott insulator occurs when the bulk spin gap solely closes.

General connection between time-local and time-nonlocal perturbation expansions

There exist two canonical approaches to describe open quantum systems by a time-evolution equation: the Nakajima-Zwanzig quantum master equation, featuring a time-nonlocal memory kernel $\mathcal{K}$, and the time-convolutionless equation with a time-local generator $\mathcal{G}$. These key quantities have recently been shown to be connected by an exact fixed-point relation [Phys. Rev. X 11, 021041 (2021)]. Here we show that this implies a recursive relation between their perturbative expansions, allowing a series for the kernel $\mathcal{K}$ to be translated directly into a corresponding series for the more complicated generator $\mathcal{G}$. This leads to an elegant way of computing the generator using well-developed, standard memory-kernel techniques for strongly interacting open systems. Moreover, it allows for an unbiased comparison of time-local and time-nonlocal approaches independent of the particular technique chosen to calculate expansions of $\mathcal{K}$ and $\mathcal{G}$ (Nakajima-Zwanzig projections, real-time diagrams, etc.). We illustrate this for leading and next-to-leading order calculations of $\mathcal{K}$ and $\mathcal{G}$ for the single impurity Anderson model using both the bare expansion in the system-environment coupling and a more advanced renormalized series. We compare the different expansions obtained, quantify the legitimacy of the generated dynamics (complete positivity) and benchmark with the exact result in the non-interacting limit.

Performance of reservoir discretizations in quantum transport simulations.

Quantum transport simulations often use explicit, yet finite, electronic reservoirs. These should converge to the correct continuum limit, albeit with a trade-off between discretization and computational cost. Here, we study this interplay for extended reservoir simulations, where relaxation maintains a bias or temperature drop across the system. Our analysis begins in the non-interacting limit, where we parameterize different discretizations to compare them on an even footing. For many-body systems, we develop a method to estimate the relaxation that best approximates the continuum by controlling virtual transitions in Kramers turnover for the current. While some discretizations are more efficient for calculating currents, there is little benefit with regard to the overall state of the system. Any gains become marginal for many-body, tensor network simulations, where the relative performance of discretizations varies when sweeping other numerical controls. These results indicate that typical reservoir discretizations have little impact on numerical costs for certain computational tools. The choice of a relaxation parameter is nonetheless crucial, and the method we develop provides a reliable estimate of the optimal relaxation for finite reservoirs.

Gross-Neveu-Heisenberg criticality from competing nematic and antiferromagnetic orders in bilayer graphene

We study the phase diagram of an effective model of competing nematic and antiferromagnetic orders of interacting electrons on the Bernal-stacked honeycomb bilayer, as relevant for bilayer graphene. In the noninteracting limit, the model features a semimetallic ground state with quadratic band touching points at the Fermi level. Taking the effects of short-range interactions into account, we demonstrate the existence of an extended region in the mean-field phase diagram characterized by coexisting nematic and antiferromagnetic orders. By means of a renormalization group approach, we reveal that the quantum phase transition from nematic to coexistent nematic-antiferromagnetic orders is continuous and characterized by emergent Lorentz symmetry. It falls into the $(2+1)$-dimensional relativistic Gross-Neveu-Heisenberg quantum universality class, which has recently been much investigated in the context of interacting Dirac systems in two spatial dimensions. The coexistence-to-antiferromagnetic transition, by contrast, turns out to be weakly first order as a consequence of the absence of continuous spatial rotational symmetry on the honeycomb bilayer. Implications for experiments in bilayer graphene are discussed.

Pair-correlation ansatz for the ground state of interacting bosons in an arbitrary one-dimensional potential

We derive and describe a very accurate variational scheme for the ground state of the system of a few ultra-cold bosons confined in one-dimensional traps of arbitrary shapes. It is based on assumption that all inter-particle correlations have two-body nature. By construction, the proposed ansatz is exact in the noninteracting limit, exactly encodes boundary conditions forced by contact interactions, and gives full control on accuracy in the limit of infinite repulsions. We show its efficiency in a whole range of intermediate interactions for different external potentials. Our results manifest that for generic non-parabolic potentials mutual correlations forced by interactions cannot be captured by distance-dependent functions.