Search

Valence bond (VB) states as the formation mechanism of Cooper pairs, eventually leading to high-temperature superconductivity, remain a controversial topic. Although various VB-like states find variational relevance in the description of specific spin models and quantum spin liquids, in the realm of many-body fermionic Hamiltonians, the evidence for such states as ground states wave functions remains elusive, challenging the valence-bond pairing mechanism. Here, we present evidence of a VB ground state with pairing tendencies, particularly at finite doping. We achieved this for the generic two-orbital Hubbard model in low dimension, where the VB states can be associated with the presence of the topological order manifested by edge states. Utilizing density-matrix renormalization group calculations, the study reveals key properties akin to those observed in superconductors' phase diagrams, such as pairing restricted to the regime of small but nonzero doping, presence of coherent pairs, and density oscillations in the charge sector.

We study tilted chains of spinless fermions in the presence of the nearest-neighbor density-density interaction for which the noninteracting counterpart displays Stark localization. We demonstrate that the latter two-body interaction can be decomposed into two (orthogonal) parts which, respectively, commute and do not commute with the single-particle Hamiltonian. We derive an explicit form of the noncommuting part that decreases with tilt and describes the nearest-neighbor correlated hopping and the pair-hopping interaction. When the density-density coupling is replaced by the pair-hopping interaction of the same magnitude then the charge dynamics may be faster by a few orders of magnitude than in the original model.

The integration and parallelization of nanopore sensors are essential for improving the throughput of nanopore measurements. Solid-state nanopores traditionally have been used in isolation, which prevents the realization of their full potential in applications. In this study, we present the microfluidic integration of an array of 30 nanopores, which, to our knowledge, is the highest number reported to date. Our microfluidic network was fabricated using high-resolution epoxy photoresists, and the solid-state membranes were bonded through a dry process using complementary surface chemistries. We successfully measured integrated nanopores using external electrodes. This paper discusses the limitations of our methods, particularly concerning microfluidic interfacing and scaling to higher channel counts. Additionally, we present theoretical analysis of current blockades and noise in integrated nanopores, predicting that maintaining low series resistance between the nanopore and electrode is crucial for resolving short events., Comment: 18 pages, 9 figures

The dynamics of disordered two-dimensional systems is much less understood than the dynamics of disordered chains, mainly due to the lack of appropriate numerical methods. We demonstrate that a single-trajectory version of the fermionic truncated Wigner approximation (fTWA) gives unexpectedly accurate results for the dynamics of one-dimensional (1D) systems with moderate or strong disorder. Additionally, the computational complexity of calculations carried out within this approximation is small enough to enable studies of two-dimensional (2D) systems larger than standard fTWA. Using this method, we analyze the dynamics of interacting spinless fermions propagating on disordered 1D and 2D lattices. We find for both spatial dimensions that the imbalance exhibits a universal dependence on the rescaled time $t/\xi_W$, where in 2D the time-scale $\xi_W$ follows a stretched-exponential dependence on disorder strength., Comment: 11 pages, 8 figures

Eigenstate thermalization has been shown to occur for few-body observables in a wide range of nonintegrable interacting models. For intensive observables that are sums of local operators, because of their polynomially vanishing Hilbert-Schmidt norm, weak eigenstate thermalization occurs in quadratic and integrable interacting systems. Here, we unveil a novel weak eigenstate thermalization phenomenon that occurs in quadratic models whose single-particle sector exhibits quantum chaos (quantum-chaotic quadratic models) and in integrable interacting models. In such models, we show that there are few-body observables with a nonvanishing Hilbert-Schmidt norm that are guarrantied to exhibit a polynomially vanishing variance of the diagonal matrix elements, a phenomenon we dub normal weak eigenstate thermalization. For quantum-chaotic quadratic Hamiltonians, we prove that normal weak eigenstate thermalization is a consequence of single-particle eigenstate thermalization, i.e., it can be viewed as a manifestation of quantum chaos at the single-particle level. We report numerical evidence of normal weak eigenstate thermalization for quantum-chaotic quadratic models such as the 3D Anderson model in the delocalized regime and the power-law random banded matrix model, as well as for the integrable interacting spin-$\frac{1}{2}$ XYZ and XXZ models.

We investigate chains of interacting spinless fermions subject to a finite external field $F$ (also called Stark chains) and focus on the regime where the charge thermalization follows the subdiffusive hydrodynamics. First, we study reduced models conserving the dipole moment and derive an explicit Einstein relation which links the subdiffusive transport coefficient with the correlations of the dipolar current. This relation explains why the decay rate, $\Gamma_q$, of the density modulation with wave-vector $q$ shows $q^4$-dependence. In the case of the Stark model, a similar Einstein relation is also derived and tested using various numerical methods. They confirm an exponential reduction of the transport coefficient with increasing $F$. On the other hand, our study of the Stark model indicates that upon increasing $q$ there is a crossover from subdiffusive behavior, $\Gamma_q \propto q^4$, to the normal diffusive relaxation, $\Gamma_q \propto q^2$, at the wave vector $q^*$ which vanishes for $F \to 0$., Comment: 4+3 pages and 4+3 figures

Hilbert space fragmentation is an ergodicity breaking phenomenon, in which Hamiltonian shatters into exponentially many dynamically disconnected sectors. In many fragmented systems, these sectors can be labelled by statistically localized integrals of motion, which are {\em nonlocal} operators. We study the paradigmatic nearest-neighbor pair hopping (PH) model exhibiting the so-called strong fragmentation. We show that this model hosts local integrals of motion (LIOMs), which correspond to frozen density modes with long wavelengths. The latter modes become subdiffusive when longer-range pair hoppings are allowed. Finally, we make a connection with a tilted (Stark) chain. Contrary to the dipole-conserving effective models, the tilted chain is shown to support either Hamiltonian or dipole moment as an LIOM. Numerical results are obtained from a numerical algorithm, in which finding LIOMs is reduced to the data compression problem.

Relaxation rates in nearly integrable systems usually increase quadratically with the strength of the perturbation that breaks integrability. We show that the relaxation rates can be significantly smaller in systems that are integrable along two intersecting lines in the parameter space. In the vicinity of the intersection point, the relaxation rates of certain observables increase with the fourth power of the distance from this point, whereas for other observables one observes standard quadratic dependence on the perturbation. As a result, one obtains exceedingly long-living prethermalization but with a reduced number of the nearly conserved operators. We show also that such a scenario can be realized in spin ladders., Comment: 12 pages, 8 figures (including Supplementary Material)

We study the transport dynamics of an interacting tilted (Stark) chain. We show that the crossover between diffusive and subdiffusive dynamics is governed by $F\sqrt{L}$, where $F$ is the strength of the field, and $L$ is the wave-length of the excitation. While the subdiffusive dynamics persist for large fields, the corresponding transport coefficient is exponentially suppressed with $F$ so that the finite-time dynamics appear almost frozen. We explain the crossover scale between the diffusive and subdiffusive transport by bounding the dynamics of the dipole moment for arbitrary initial state. We also prove its emergent conservation at infinite temperature. Consequently, the studied chain is one of the simplest experimentally realizable models for which numerical data are consistent with the hydrodynamics of fractons., Comment: 8 pages, 4 figures (including Supplementary Material). Changes in the second version: The typo of the first version in the arrow directions of Fig. 1 panel (b) and Fig. S1 (upper panel) is corrected

The random-field spin-1/2 XXZ chains, and the corresponding Anderson insulators of spinless fermions with density-density interaction, have been intensively studied in the context of many-body localization. However, we recently argued [B. Krajewski et al., Phys. Rev. Lett. 129, 260601(2022)] that the two-body density-density interaction in these models is not generic since only a small fraction of this interaction represents a true local perturbation to the Anderson insulator. Here we study ergodicity of strongly disordered Anderson insulator chains choosing other forms of the two-body interaction for which the strength of the true perturbation is of the same order of magnitude as the bare two-body interaction. Focusing on the strong interaction regime, numerical results for the level statistics and the eigenstate thermalization hypothesis are consistent with emergence of ergodicity at arbitrary strong disorder.

We investigate full quantum mechanical evolution of two electrons nonlinearly coupled to quantum phonons and simulate the dynamical response of the system subject to a short spatially uniform optical pulse that couples to dipole-active vibrational modes. Nonlinear electron-phonon coupling can either soften or stiffen the phonon frequency in the presence of electron density. In the former case, an external optical pulse tuned just below the phonon frequency generates attraction between electrons and leads to a long-lived bound state even after the optical pulse is switched off. It originates from a dynamical modification of the self-trapping potential that induces a metastable state. By increasing the pulse frequency, the attractive electron-electron interaction changes to repulsive. Two sequential optical pulses with different frequencies can switch between attractive and repulsive interaction. Finally, we show that the pulse-induced binding of electrons is shown to be efficient also for weakly dispersive optical phonons, in the presence anharmonic phonon spectrum and in two dimensions., Comment: 14 pages

One of the most famous quantum systems with topological properties, the spin $\mathcal{S}=1$ antiferromagnetic Heisenberg chain, is well-known to display exotic $\mathcal{S}=1/2$ edge states. However, this spin model has not been analyzed from the more general perspective of strongly correlated systems varying the electron-electron interaction strength. Here, we report the investigation of the emergence of the Haldane edge in a system of interacting electrons -- the two-orbital Hubbard model -- with increasing repulsion strength $U$ and Hund interaction $J_\mathrm{H}$. We show that interactions not only form the magnetic moments but also form a topologically nontrivial fermionic many-body ground-state with zero-energy edge states. Specifically, upon increasing the strength of the Hubbard repulsion and Hund exchange, we identify a sharp transition point separating topologically trivial and nontrivial ground-states. Surprisingly, such a behaviour appears already at rather small values of the interaction, in a regime where the magnetic moments are barely developed.

We discuss the interplay between many-body localization and spin-symmetry. To this end, we study the time evolution of several observables in the anisotropic t-J model. Like the Hubbard chain, the studied model contains charge and spin degrees of freedom, yet it has smaller Hilbert space and thus allows for numerical studies of larger systems. We compare the field disorder that breaks the Z_2 spin symmetry and a potential disorder that preserves the latter symmetry. In the former case, sufficiently strong disorder leads to localization of all studied observables, at least for the studied system sizes. However, in the case of symmetry-preserving disorder, we observe that odd operators under the Z_2 spin transformation relax towards the equilibrium value at relatively short time scales that grow only polynomially with the disorder strength. On the other hand, the dynamics of even operators and the level statistics within each symmetry sector are consistent with localization. Our results indicate that localization exists within each symmetry sector for symmetry preserving disorder. Odd operators' apparent relaxation is due to their time evolution between distinct symmetry sectors.

In recent years the ergodicity of disordered spin chains has been investigated via extensive numerical studies of the level statistics or the transport properties. However, a clear relationship between these results has yet to be established. We present the relation between the diffusion constant and the energy-level structure, which leads to the Thouless localization criterion. Together with the exponential-like dependence of the diffusion constant on the strength of quasiperiodic or random fields, the Thouless criterion explains the nearly linear drift with the system size of the crossover/transition to the nonergodic regime. Moreover, we show that the Heisenberg spin chain in the presence of the quasiperiodic fields can be well approached via a sequence of simple periodic systems, where diffusion remains finite even at large fields.

We study sample-to-sample fluctuations of the gap ratio in the energy spectra in finite disordered spin chains. The chains are described by the random-field Ising model and the Heisenberg model. We show that away from the ergodic/nonergodic crossover, the fluctuations are correctly captured by the Rosenzweig-Porter (RP) model. However, fluctuations in the microscopic models significantly exceed those in the RP model in the vicinity of the crossover. We show that upon introducing an extension to the RP model, one correctly reproduces the fluctuations in all regimes, i.e., in the ergodic and nonergodic regimes as well as at the crossover between them. Finally, we demonstrate how to reduce the sample-to-sample fluctuations in both studied microscopic models.

The isotropic Heisenberg chain represents a particular case of an integrable many-body system exhibiting superdiffusive spin transport at finite temperatures. Here, we show that this model has distinct properties also at finite magnetization $m\ne0$, even upon introducing the SU(2) invariant perturbations. Specifically, we observe nonmonotonic dependence of the diffusion constant ${\cal D}_0(\Delta)$ on the spin anisotropy $\Delta$, with a pronounced maximum at $\Delta =1$. The latter dependence remains true also in the zero magnetization sector, with superdiffusion at $\Delta=1$ that is remarkably stable against isotropic perturbation (at least in finite-size systems), consistent with recent experiments with cold atoms., Comment: 5+5 pages, 4+5 figures

Thermalization (generalized thermalization) in nonintegrable (integrable) quantum systems requires two ingredients: equilibration and agreement with the predictions of the Gibbs (generalized Gibbs) ensemble. We prove that observables that exhibit eigenstate thermalization in single-particle sector equilibrate in many-body sectors of quantum-chaotic quadratic models. Remarkably, the same observables do not exhibit eigenstate thermalization in many-body sectors (we establish that there are exponentially many outliers). Hence, the generalized Gibbs ensemble is generally needed to describe their expectation values after equilibration, and it is characterized by Lagrange multipliers that are smooth functions of single-particle energies.

Using the truncated Wigner approximation (TWA) we study quench dynamics of two-dimensional lattice systems consisting of interacting spinless fermions with potential disorder. First, we demonstrate that the semiclassical dynamics generally relaxes faster than the full quantum dynamics. We obtain this result by comparing the semiclassical dynamics with exact diagonalization and Lanczos propagation of one-dimensional chains. Next, exploiting the TWA capabilities of simulating large lattices, we investigate how the relaxation rates depend on the dimensionality of the studied system. We show that strongly disordered one-dimensional and two-dimensional systems exhibit a transient, logarithmic-in-time relaxation, which was recently established for one-dimensional chains. Such relaxation corresponds to the infamous $1/f$-noise at strong disorder., Comment: 9 pages, 9 figures

We consider a chain of interacting fermions with random disorder that was intensively studied in the context of many-body localization. We show that only a small fraction of the two-body interaction represents a true local perturbation to the Anderson insulator. While this true perturbation is nonzero at any finite disorder strength W, it decreases with increasing W. This establishes a view that the strongly disordered system should be viewed as a weakly perturbed integrable model, i.e., a weakly perturbed Anderson insulator. As a consequence, the latter can hardly be distinguished from a strictly integrable system in finite-size calculations at large W. We then introduce a rescaled model in which the true perturbation is of the same order of magnitude as the other terms of the Hamiltonian, and show that the system remains ergodic at arbitrary large disorder.

Purely ballistic transport is a rare feature even for integrable models. By numerically studying the Heisenberg chain with the power-law exchange, \mbox{$J\propto1/r^\alpha$}, where $r$ is a distance, we show that for spin anisotropy $\Delta \simeq \exp(-\alpha+2)$ the system exhibits a quasiballistic spin transport and the presence of fermionic excitation which do not decay up to extremely long times $\sim10^3/J$. This conclusion is reached on the base of the dynamics of spin domains, the dynamical spin conductivity, inspecting the matrix elements of the spin-current operator, and by the analysis of most conserved operators. Our results smoothly connects two models where fully ballistic transport is present: free particles with nearest-neighbor hopping and the isotropic Haldane-Shastry model.

The anomalous spin diffusion of the integrable easy-axis Heisenberg chain originates in the ballistic transport of symmetry sectors with nonzero magnetization. Ballistic transport is replaced by normal dissipative transport in all magnetization sectors upon introducing the integrability-breaking perturbations, including external driving. Such behavior implies that the diffusion constant obtained for the integrable model is relevant for the spread of spin excitations but not for the spin conductivity. We present numerical results for closed systems and driven open systems, indicating that the diffusion constant shows a discontinuous variation as the function of perturbation strength.

We numerically study the relaxation of correlation functions in weakly perturbed integrable XXZ chain. The decay of the spin-current and the energy-current correlations at zero magnetization are well described by single, but quite distinct, relaxation rates governed by the square of the perturbation strength $g$. However, at finite magnetization a single correlation function reveals multiple relaxation rates. The result can be understood in terms of multi-scale relaxation scenario, where various relaxation times are linked with various quantities which are conserved in the reference integrable system. On the other hand, the correlations of non-commuting quantities, being conserved at particular anisotropies $\Delta$, decay non-exponentially with characteristic time scale linear in $g$.

We study dynamical correlation functions in the random-field Heisenberg chain, which probes the relaxation times at different length scales. Firstly, we show that the relaxation time associated with the dynamical imbalance (examining the relaxation at the smallest length scale) decreases with disorder much faster than the one determined by the dc conductivity (probing the global response of the system). We argue that the observed dependence of relaxation on the length scale originates from local nonresonant regions. The latter have particularly long relaxation times or remain frozen, allowing for nonzero dc transport via higher-order processes. Based on the numerical evidence, we introduce a toy model that suggests that the nonresonant regions asymptotic dynamics are essential for the proper understanding of the disordered chains with many-body interactions.

We investigate the high-temperature dynamical conductivity $\sigma(\omega)$ in two one-dimensional integrable quantum lattice models: the anisotropic XXZ spin chain and the Hubbard chain. The emphasis is on the metallic regime of both models, where besides the ballistic component, the regular part of conductivity might reveal a diffusive-like transport. To resolve the low-frequency dynamics, we upgrade the microcanonical Lanczos method enabling studies of finite-size systems with up to $L\leq 32$ sites for XXZ spin model with the frequency resolution $\delta \omega \sim 10^{-3} J$. Results for the XXZ chain reveal a fine structure of $\sigma(\omega)$ spectra, which originates from the discontinuous variation of the stiffness, previously found at commensurate values of the anisotropy parameter $\Delta$. Still, we do not find a clear evidence for a diffusive component, at least not for commensurate values of $\Delta$, particularly for $\Delta =0.5$, as well as for $\Delta \to 0$. Similar is the conclusion for the Hubbard model away from half-filling, where the spectra reveal more universal behavior.

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.

We study the ballistic transport in integrable lattice models, i.e., the spin XXZ and Hubbard chains, close to the noninteracting limit. The stiffnesses of spin and charge currents reveal, at high temperatures, a discontinuous reduction (jump) when the interaction is introduced. We show that the jumps are related to the large degeneracy of the parent noninteracting models. These degeneracies are properly captured by the degenerate perturbation calculations which may be performed for large systems. We find that the discontinuities and the quasilocality of the conserved current in this limit can be traced back to the nonlocal character of an effective interaction. From the latter observation we identify a class of observables which show discontinuities in both models. We also argue that the known local conserved quantities are insufficient to explain the stiffnesses in the Hubbard chain in the regime of weak interaction., Comment: updated version of the manuscript

Topological phases of matter are among the most intriguing research directions in Condensed Matter Physics. It is known that superconductivity induced on a topological insulator's surface can lead to exotic Majorana modes, the main ingredient of many proposed quantum computation schemes. In this context, the iron-based high critical temperature superconductors are a promising platform to host such an exotic phenomenon in real condensed-matter compounds. The Coulomb interaction is commonly believed to be vital for the magnetic and superconducting properties of these systems. This work bridges these two perspectives and shows that the Coulomb interaction can also drive a canonical superconductor with orbital degrees of freedom into the topological state. Namely, we show that above a critical value of the Hubbard interaction the system simultaneously develops spiral spin order, a highly unusual triplet amplitude in superconductivity, and, remarkably, Majorana fermions at the edges of the system.

We examine the standard model of many-body localization (MBL), i.e., the disordered chain of interacting spinless fermions, by representing it as the network in the many-body (MB) basis of noninteracting localized Anderson states. By studying eigenstates of the full Hamiltonian, for strong disorders we find that the dynamics is confined up to very long times to disconnected MB clusters in the Fock space. By keeping only resonant contributions and simplifying the quantum problem to rate equations (REs) for MB states, in analogy with percolation problems, the MBL transition is located via the universal cluster distribution and the emergence of the macroscopic cluster. On the ergodic side, our approximate RE approach to the relaxation processes captures well the diffusion transport, as found for the full quantum model. In a broad transient regime, we find an anomalous, i.e., subdiffusivelike, transport, emerging from weak links between MB states., Comment: 7 pages, 9 figures

Low dimensional structures in the non-trivial topological phase can host the in-gap Majorana bound states, identified experimentally as zero-bias peaks in differential conductance. Theoretical methods for studying Majorana modes are mostly based on the bulk-boundary correspondence or exact diagonalization of finite systems via, e.g., Bogoliubov-de Gennes formalism. In this paper, we develop an efficient method for identifying the Majorana in-gap (edge) states via looking for extreme eigenvalues of symmetric matrices. The presented approach is based on the Krylov method and allows for study the spatial profile of the modes as well as the spectrum of the system. The advantage of this method is the calculation cost, which shows linear dependence on the number of lattice sites. The latter problem may be solved for very large clusters of arbitrary shape/geometry. In order to demonstrate the efficiency of our approach, we study two- and three-dimensional clusters described by the Kitaev and Rashba models for which we determine the number of Majorana modes and calculate their spatial structures. Additionally, we discuss the impact of the system size on the physical properties of the topological phase of the magnetic nanoisland deposited on the superconducting surface. In this case, we have shown that the eigenvalues of the in-gap states depend on the length of the system edge., Comment: 15 pages, 15 figures

We study a quantum particle coupled to hard-core bosons and propagating on disordered ladders with $R$ legs. The particle dynamics is studied with the help of rate equations for the boson-assisted transitions between the Anderson states. We demonstrate that for finite $R < \infty$ and sufficiently strong disorder the dynamics is subdiffusive, while the two-dimensional planar systems with $R\to \infty$ appear to be diffusive for arbitrarily strong disorder. The transition from diffusive to subdiffusive regimes may be identified via statistical fluctuations of resistivity. The corresponding distribution function in the diffusive regime has fat tails which decrease with the system size $L$ much slower than $1/\sqrt{L}$. Finally, we present evidence that similar non--Gaussian fluctuations arise also in standard models of many-body localization, i.e., in strongly disordered quantum spin chains.

We study dynamics of a single qubit encoded in two pairs of Majorana modes, whereby each pair is hosted on a trijunction described by the Kitaev model extended by many-body interactions. We demonstrated that the challenging phase-gate may be efficiently implemented via braiding of partially overlapping modes. Although such qubit acquires both geometric and dynamical phases during the braiding protocol, the latter phase may be eliminated if the Majorana modes are hosted by systems with appropriate particle-hole symmetry., Comment: 9 pages, 8 figures

Even though foundations of the eigenstate thermalization hypothesis (ETH) are based on random matrix theory, physical Hamiltonians and observables substantially differ from random operators. One of the major challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we focus on their impact on fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.

A quantum particle propagates subdiffusively on a strongly disordered chain when it is coupled to itinerant hard-core bosons. We establish a generalized Einstein relation (GER) that relates such subdiffusive spread to an unusual time-dependent drift velocity, which appears as a consequence of a constant electric field. We show that GER remains valid much beyond the regime of the linear response. Qualitatively, it holds true up to strongest drivings when the nonlinear field-effects lead to the Stark-like localization. Numerical calculations based on full quantum evolution are substantiated by much simpler rate equations for the boson-assisted transitions between localized Anderson states., Comment: Accepted for publication in Phys. Rev. Lett