Single-particle Eigenstates Research Articles

Single-quasiparticle eigenstate thermalization.

Quadratic Hamiltonians that exhibit single-particle quantum chaos are called quantum-chaotic quadratic Hamiltonians. One of their hallmarks is single-particle eigenstate thermalization introduced in Łydżba etal. [Phys. Rev. B 104, 214203 (2021)2469-995010.1103/PhysRevB.104.214203], which describes statistical properties of matrix elements of observables in single-particle eigenstates. However, the latter has been studied only in quantum-chaotic quadratic Hamiltonians that obey the U(1) symmetry. Here, we focus on quantum-chaotic quadratic Hamiltonians that break the U(1) symmetry and, hence, their "single-particle" eigenstates are actually single-quasiparticle excitations introduced on the top of a many-body state. We study their wave functions and matrix elements of one-body observables, for which we introduce the notion of single-quasiparticle eigenstate thermalization. Focusing on spinless fermion Hamiltonians in three dimensions with local hopping, pairing, and on-site disorder, we also study the properties of disorder-induced near zero modes, which give rise to a sharp peak in the density of states at zero energy. Finally, we numerically show equilibration of observables in many-body eigenstates after a quantum quench. We argue that the latter is a consequence of single-quasiparticle eigenstate thermalization, in analogy to the U(1) symmetric case from Łydżba etal. [Phys. Rev. Lett. 131, 060401 (2023)0031-900710.1103/PhysRevLett.131.060401].

Designer quantum states on a fractal substrate: Compact localization, flat bands and the edge modes

Compact localized single particle eigenstates on a deterministic fractal substrate, modelled by a triangular Sierpinski gasket of arbitrarily large size, are unravelled and examined analytically. We prescribe an exact real space renormalization group (RSRG) decimation scheme within a tight binding formalism to discern these states, and argue that the number of such states can be infinite if the fractal substrate is enlarged to its thermodynamic limit. Interestingly, these localized states turn out to populate the non-dispersive, flat bands in a periodic array of Sierpinski gasket motifs, however large they may be. Our results match and corroborate the recently observed compact localized, flat band states engineered on two dimensional photonic waveguide networks with a fractal geometry, and provide a whole subset of them, which, in principle, should be observable in fractal photonic lattice experiments.

Spectral and dynamical characters of 1D incommensurate optical lattices with [formula omitted]-symmetry

In this paper, we investigate the spectral and dynamical properties of a 1D incommensurate optical lattice, which displays a Parity-time (PT) symmetry described by the Non-Hermitian Aubry–André potential. We show through the spectral analysis of the eigenvalues and the associated eigenfunctions that, taking into account a non-Hermitian variant of the site energy modulated by a cosine form, leads the system to undergo a transition from the unbroken-PT phase to the broken-PT phase, which corresponds here to the delocalized-to-localized phase transition. Our findings also indicate that the critical point of transition can be selectively adjusted and consequently, the single-particle eigenstates taken in the metallic regime of the original Aubry–André model thus move from a conducting state to an insulating one upon an increase of the complex phase parameter. Furthermore, we also investigate the dynamical features of the system’s wave function according to different parameters.

Induced position-dependent spin–orbit coupling in atomic Bose–Einstein condensate

In this paper, we study-induced spin–orbit coupling (SOC) in one-dimensional ultracold quantum gases composed of atoms of different species. One species is subjected to an equal combination of Rashba and Dresselhaus SOC generated by Raman transition. The two species interact with each other through spin-independent and spin-exchange contact interactions. The spin-exchange interaction introduces a locking pattern of the momentum and spin degrees of freedom for the species without direct SOC. For small spin-independent interactions, the two species overlap and the induced SOC is effective. For large spin-independent interactions, however, the two species keep separate from each other and the induced SOC is negligible. We propose to generate inhomogeneous SOC by exploiting density engineering technique applied to the directly spin–orbit coupled Bose–Einstein condensate. When the effective potential is of Gaussian type, the single-particle eigenenergies and eigenstates are calculated by exact diagonalization method. For general cases, mean-field ground states are obtained by numerically searching for the minimum of energy functional. With the density engineering technique, it is possible to produce hybrid structures of plane wave, stripe and/or zero-momentum phases.

Anomalous Josephson effect and rectification in junctions between Floquet topological superconductors

Periodically driven Kitaev chains are known to exhibit novel Floquet Majorana fermions and anomalous Floquet end modes. The fact that `quasienergy in Floquet systems is periodic' poses a difficulty in defining the ground state in periodically driven systems. To overcome this problem, we start with the ground state of the undriven Kitaev chain and gradually switch on the periodic driving in chemical potential over a timescale $\tau$. The extent to which the single particle eigenstates of the undriven system get distributed among the Floquet states upon driving can be characterized by inverse participation ratio. In a Josephson junction between two periodically driven Kitaev chains not differing in phase of the pair potentials but differing in phases of the driving potentials, a net average current flows from one superconductor to the other. We term such a current anomalous current. Further, we study current phase relation in junctions between two periodically driven superconductors and find that the system exhibits nonequilibrium Josephson diode effect. The Floquet Majorana end modes and anomalous Floquet end modes whenever present contribute significantly to the anomalous current and the diode effect. Further, the anomalous current and the nonequilibrium Josephson diode effect survive when the periodic driving is adiabatically switched on.

Effect of spin–orbit coupling in one-dimensional quasicrystals with power-law hopping

In the one-dimensional quasiperiodic Aubry-Andr\'{e}-Harper Hamiltonian with nearest-neighbor hopping, all single-particle eigenstates undergo a phase transition from ergodic to localized states at a critical disorder strength $W_c/t = 2.0$. There is no mobility edge in this system. However, in the presence of power-law hopping having the form $1/r^a$, beyond a critical disorder strength mobility edge appears for $a > 1$, while, for $0< a\leq 1$, a multifractal edge separates the extended and the multifractal states. In both these limits, depending on the strength of the disorder, lowest $\beta^s L$ states are delocalized. We have found that, in the presence of the spin-orbit coupling, the critical disorder strength is always larger irrespective of the value of the parameter $a$. Furthermore, we demonstrate that for $a\leq 1$, in the presence of spin-orbit coupling, there exists multiple multifractal edges, and the energy spectrum splits up into alternative bands of delocalized and multifractal states. Moreover, the location of the multifractal edges are generally given by the fraction $(\beta^s \pm \beta^m)$. The qualitative behavior of the energy spectrum remains unaffected for $a > 1$. However, in contrast to the previously reported results, we find that in this limit, similar to the other case, multiple mobility edges can exist with or without the spin-orbit coupling.

Exceptional spectrum and dynamic magnetization

A macroscopic effect can be induced by a local non-Hermitian term in a many-body system, when it manifests simultaneously level coalescence of a full real degeneracy spectrum, leading to exceptional spectrum. In this paper, we propose a family of systems that support such an intriguing property. It is generally consisted of two arbitrary identical Hermitian sub-lattices in association with unidirectional couplings between them. We show exactly that all single-particle eigenstates coalesce in pairs even only single unidirectional coupling appears. It means that all possible initial states obey the exceptional dynamics, resulting in some macroscopic phenomena, which never appears in a Hermitian system. As an application, we study the dynamic magnetization induced by complex fields in an itinerant electron system. It shows that an initial saturated ferromagnetic state at half-filling can be driven into its opposite state according to the dynamics of high-order exceptional point. Any Hermitian quench term cannot realize a steady opposite saturated ferromagnetic state. Numerical simulations for the dynamical processes of magnetization are performed for several representative situations, including lattice dimensions, global random and local impurity distributions. It shows that the dynamic magnetization processes exhibit universal behavior.

Tight-binding billiards.

Recent works have established universal entanglement properties and demonstrated validity of single-particle eigenstate thermalization in quantum-chaotic quadratic Hamiltonians. However, a common property of all quantum-chaotic quadratic Hamiltonians studied in this context so far is the presence of random terms that act as a source of disorder. Here we introduce tight-binding billiards in two dimensions, which are described by noninteracting spinless fermions on a disorder-free square lattice subject to curved open (hard-wall) boundaries. We show that many properties of tight-binding billiards match those of quantum-chaotic quadratic Hamiltonians: The average entanglement entropy of many-body eigenstates approaches the random matrix theory predictions and one-body observables in single-particle eigenstates obey the single-particle eigenstate thermalization hypothesis. On the other hand, a degenerate subset of single-particle eigenstates at zero energy (i.e., the zero modes) can be described as chiral particles whose wave functions are confined to one of the sublattices.

Symmetry-Protected Transport through a Lattice with a Local Particle Loss.

We study particle transport through a chain of coupled sites connected to free-fermion reservoirs at both ends, subjected to a local particle loss. The transport is characterized by calculating the conductance and particle density in the steady state using the Keldysh formalism for open quantum systems. In addition to a reduction of conductance, we find that transport can remain (almost) unaffected by the loss for certain values of the chemical potential in the lattice. We show that this "protected" transport results from the spatial symmetry of single-particle eigenstates. At a finite voltage, the density profile develops a drop at the lossy site, connected to the onset of nonballistic transport.

Reduced Scaling of Optimal Regional Orbital Localization via Sequential Exhaustion of the Single-Particle Space.

Wannier functions have become a powerful tool in the electronic structure calculations of extended systems. The generalized Pipek-Mezey Wannier functions exhibit appealing characteristics (e.g., reaching an optimal localization and the separation of the σ–π orbitals) compared with other schemes. However, when applied to giant nanoscale systems, the orbital localization suffers from a large computational cost overhead if one is interested in localized states in a small fragment of the system. Herein, we present a swift, efficient, and robust approach for obtaining regionally localized orbitals of a subsystem within the generalized Pipek-Mezey scheme. The proposed algorithm introduces a reduced work space and sequentially exhausts the entire orbital space until the convergence of the localization functional. It tackles systems with ∼10000 electrons within 0.5 h with no loss in localization quality compared to the traditional approach. Regionally localized orbitals with a higher extent of localization are obtained via judiciously extending the subsystem’s size. Exemplifying on large bulk and a 4 nm wide slab of diamond with an NV– center, we demonstrate the methodology and discuss how the choice of the localization region affects the excitation energy of the defect. Furthermore, we show how the sequential algorithm is easily extended to stochastic methodologies that do not provide individual single-particle eigenstates. It is thus a promising tool to obtain regionally localized states for solving the electronic structure problems of a subsystem embedded in giant condensed systems.

Automorphic Bloch theorems for hyperbolic lattices

Hyperbolic lattices are a new form of synthetic quantum matter in which particles effectively hop on a discrete tessellation of two-dimensional (2D) hyperbolic space, a non-Euclidean space of uniform negative curvature. To describe the single-particle eigenstates and eigenenergies for hopping on such a lattice, a hyperbolic generalization of band theory was previously constructed, based on ideas from algebraic geometry. In this hyperbolic band theory, eigenstates are automorphic functions, and the Brillouin zone is a higher-dimensional torus, the Jacobian of the compactified unit cell understood as a higher-genus Riemann surface. Three important questions were left unanswered: whether a band theory can be expected to hold for a non-Euclidean lattice, where translations do not generally commute; whether a formal Bloch theorem can be rigorously established; and whether hyperbolic band theory can describe finite lattices realized in an experiment. In the present work, we address all three questions simultaneously. By formulating periodic boundary conditions for finite but arbitrarily large lattices, we show that a generalized Bloch theorem can be rigorously proved but may or may not involve higher-dimensional irreducible representations (irreps) of the nonabelian translation group, depending on the lattice geometry. Higher-dimensional irreps correspond to points in a moduli space of higher-rank stable holomorphic vector bundles, which further generalizes the notion of Brillouin zone beyond the Jacobian. For a large class of finite lattices, only 1D irreps appear, and the hyperbolic band theory previously developed becomes exact.

Single-particle eigenstate thermalization in quantum-chaotic quadratic Hamiltonians

We study the matrix elements of local and nonlocal operators in the single-particle eigenstates of two paradigmatic quantum-chaotic quadratic Hamiltonians; the quadratic Sachdev-Ye-Kitaev (SYK2) model and the three-dimensional Anderson model below the localization transition. We show that they display eigenstate thermalization for normalized observables. Specifically, we show that the diagonal matrix elements exhibit vanishing eigenstate-to-eigenstate fluctuations, and a variance proportional to the inverse Hilbert space dimension. We also demonstrate that the ratio between the variance of the diagonal and the off-diagonal matrix elements is $2$, as predicted by the random matrix theory. We study distributions of matrix elements of observables and establish that they need not be Gaussian. We identify the class of observables for which the distributions are Gaussian.

Entanglement in many-body eigenstates of quantum-chaotic quadratic Hamiltonians

In a recent Letter [Phys. Rev. Lett. 125, 180604 (2020)], we introduced a closed-form analytic expression for the average bipartite von Neumann entanglement entropy of many-body eigenstates of random quadratic Hamiltonians. Namely, of Hamiltonians whose single-particle eigenstates have random coefficients in the position basis. A paradigmatic Hamiltonian for which the expression is valid is the quadratic Sachdev-Ye-Kitaev (SYK2) model in its Dirac fermion formulation. Here we show that the applicability of our result is much broader. Most prominently, it is also relevant for local Hamiltonians such as the three-dimensional (3D) Anderson model at weak disorder. Moreover, it describes the average entanglement entropy in Hamiltonians without particle-number conservation, such as the SYK2 model in the Majorana fermion formulation and the 3D Anderson model with additional terms that break particle-number conservation. We extend our analysis to the average bipartite second R\'enyi entanglement entropy of eigenstates of the same quadratic Hamiltonians, which is derived analytically and tested numerically. We conjecture that our results for the entanglement entropies of many-body eigenstates apply to quadratic Hamiltonians whose single-particle eigenstates exhibit quantum chaos, to which we refer as quantum-chaotic quadratic Hamiltonians.

Readout of quasiperiodic systems using qubits

We develop a theoretical scheme to perform a readout of the properties of a quasiperiodic system by coupling it to one or two qubits. We show that the decoherence dynamics of a single qubit coupled via a pure dephasing type term to a one-dimensional quasiperiodic system with a potential given by the Andr\'e-Aubry-Harper (AAH) model and its generalized versions (GAAH model) is sensitive to the nature of the single-particle eigenstates (SPEs). More specifically, we can use the non-Markovianity of the qubit dynamics as quantified by the backflow of information to clearly distinguish the localized, delocalized, and mixed regimes with a mobility edge of the AAH and GAAH models and evidence the transition between them. By attaching two qubits at distinct sites of the system, we demonstrate that the transport property of the quasiperiodic system is encoded in the scaling of the threshold time to develop correlations between the qubits with the distance between the qubits. This scaling can also be used to distinguish and infer different regimes of transport such as ballistic, diffusive, and no transport engendered by SPEs that are delocalized, critical, and localized respectively. In addition, the localization length of the SPEs can also be gleaned from the exponential decay of correlations at long times as a function of distance between qubits. When there is a mobility edge allowing the coexistence of different kinds of SPEs in the spectrum, such as the coexistence of localized and delocalized states in the GAAH models, we find that the transport behavior and the scaling of the threshold time with qubit separation are governed by the fastest spreading states.

Interplay of disorder and interactions in a flat-band supporting diamond chain

We systematically study the effect of disorder and interactions on a quasi-one dimensional diamond chain possessing flat bands. Disorder localizes all the single particle eigenstates, while at low disorder strengths we obtain weak flat-band based localization (FBL), at high disorder strengths, we see conventional Anderson localization (AL). The compactly localized (CL) eigenstates of flat bands show a persisting oscillatory recurrence in the study of single-particle wavepacket dynamics. For low disorder a damped oscillatory recurrence behavior is observed which is absent for high disorder. Non-interacting many particle fermion states also follow the same trend except showing a delocalizing tendency at intermediate disorder due to the fermionic statistics in the system. As interactions are switched on, for the finite-sizes that we are able to study, a non-ergodic `mixed phase' is observed at low disorder which is separated from the MBL phase at high disorder by a thermal phase at intermediate-disorder. A study of many-body nonequilibrium dynamics reinforces these findings.

Holographic unitary renormalization group for correlated electrons - I: A tensor network approach

We present a unified framework for the renormalisation of the Hamiltonian and eigenbasis of a system of correlated electrons, unveiling thereby the interplay between electronic correlations and many-particle entanglement. For this, we extend substantially the unitary renormalization group (URG) scheme introduced in Refs. [1–3]. We recast the RG as a discrete flow of the Hamiltonian tensor network, i.e., the collection of various 2n-point scattering vertex tensors comprising the Hamiltonian. The renormalisation progresses via unitary transformations that block diagonalizes the Hamiltonian iteratively via the disentanglement of single-particle eigenstates. This procedure incorporates naturally the role of quantum fluctuations. The RG flow equations possess a non-trivial structure, displaying a feedback mechanism through frequency-dependent dynamical self-energies and correlation energies. The interplay between various UV energy scales enables the coupled RG equations to flow towards a stable fixed point in the IR. The effective Hamiltonian at the IR fixed point generically has a reduced parameter space, as well as number of degrees of freedom, compared to the microscopic Hamiltonian. Importantly, the vertex RG flows are observed to govern the RG flow of the tensor network that denotes the coefficients of the many-particle eigenstates. The RG evolution of various many-particle entanglement features of the eigenbasis are, in turn, quantified through the coefficient tensor network. In this way, we show that the URG framework provides a microscopic understanding of holographic renormalisation: the RG flow of the vertex tensor network generates a eigenstate coefficient tensor network possessing a many-particle entanglement metric. We find that the eigenstate tensor network accommodates sign factors arising from fermion exchanges, and that the IR fixed point reached generically involves a trivialisation of the fermion sign factor. Several results are presented for the emergence of composite excitations in the neighbourhood of a gapless Fermi surface, as well as for the condensation phenomenon involving the gapping of the Fermi surface.

Interaction-induced topological properties of two bosons in flat-band systems

In flat-band systems, destructive interference leads to the localization of non-interacting particles and forbids their motion through the lattice. However, in the presence of interactions the overlap between neighbouring single-particle localized eigenstates may enable the propagation of bound pairs of particles. In this work, we show how these interaction-induced hoppings can be tuned to obtain a variety of two-body topological states. In particular, we consider two interacting bosons loaded into the orbital angular momentum $l=1$ states of a diamond-chain lattice, wherein an effective $\pi$ flux may yield a completely flat single-particle energy landscape. In the weakly-interacting limit, we derive effective single-particle models for the two-boson quasiparticles which provide an intuitive picture of how the topological states arise. By means of exact diagonalization calculations, we benchmark these states and we show that they are also present for strong interactions and away from the strict flat-band limit. Furthermore, we identify a set of doubly localized two-boson flat-band states that give rise to a special instance of Aharonov-Bohm cages for arbitrary interactions.

Study of counterintuitive transport properties in the Aubry-André-Harper model via entanglement entropy and persistent current

The single particle eigenstates of the Aubry-Andr\'e-Harper model are known to show a delocalization-localization transition at a finite strength of the quasi-periodic disorder. In this work, we point out that an intimate relationship exists between the sub-band structure of the spectrum and transport properties of the model. To capture the transport properties we have not only used a variety of single-particle measures like inverse participation ratio, and von Neumann entropy, but also many-particle measures such as persistent current and its variance, and many body entanglement entropy. The many-particle measures are very sensitive to the sub-band structure of the spectrum. Even in the delocalized phase, surprisingly the entanglement entropy is substantially suppressed when the Fermi level is in the band gaps whereas the persistent current is vanishingly small for the same locations of the Fermi level. The entanglement entropy seems to follow area-law exclusively for these special locations of Fermi level or filling fractions of free fermions. A study of the standard deviation of persistent current offers further distinguishing features for the special fillings. In the delocalized phase, the standard deviation vs. mean persistent current curves are discontinuous for the non-special values of filling fractions and continuous (closed) for the special values of filling fractions whereas in the localized phase, these curves become straight lines for both types of filling fractions. We have also discussed how the results depend on the system size. Our results, specially on the persistent current, can potentially be tested experimentally using the present day set-ups based on ultra-cold atoms.

Localization and interactions in topological and nontopological bands in two dimensions

A two-dimensional electron gas in a high magnetic field displays macroscopically degenerate Landau levels, which can be split into Hofstadter subbands by means of a weak periodic potential. By carefully engineering such a potential, one can precisely tune the number, bandwidths, bandgaps and Chern character of these subbands. This allows a detailed study of the interplay of disorder, interaction and topology in two dimensional systems. We first explore the physics of disorder and single-particle localization in subbands derived from the lowest Landau level, that nevertheless may have a topological nature different from that of the entire lowest Landau level. By projecting the Hamiltonian onto subbands of interest, we systematically explore the localization properties of single-particle eigenstates in the presence of quenched disorder. We then introduce electron-electron interactions and investigate the fate of many-body localization in subbands of varying topological character.

Effective two-mode description of a few ultra-cold bosons in a double-well potential

We present a construction of an improved two-mode model for modeling the dynamics of interacting ultra-cold bosons confined in a one-dimensional double well trap. Unlike in the typically used two-mode model based on the lowest single-particle eigenstates of the external potential, the improved model uses a basis of properly chosen effective wave functions originating in the many-body model. Accuracy of the improved model is examined and it is shown that within a certain limit of inter-particle interaction strength, the model recovers an exact evolution of the wells' populations much more closely than the traditional two-mode model.