Search

We investigate numerically the joint distribution of magic ($M$) and entanglement ($S$) in $N$-qubit Haar-random quantum states. The distribution $P_N(M,S)$ as well as the marginals become exponentially localized, and centered around the values $\tilde{M_2} \to N-2$ and $\tilde{S} \to N/2$ as $N\to\infty$. Magic and entanglement fluctuations are, however, found to become exponentially uncorrelated. Although exponentially many states with magic $M_2=0$ and entropy $S\approx S_\text{Haar}$ exist, they represent an exponentially small fraction compared to typical quantum states, which are characterized by large magic and entanglement entropy, and uncorrelated magic and entanglement fluctuations., Comment: 6 pages, 7 figures

We study the emergence of complexity in deep random $N$-qubit $T$-gate doped Clifford circuits, as reflected in their spectral properties and in magic generation, characterized by the stabilizer R\'enyi entropy distribution and the non-stabilizing power of the circuit. For pure (undoped) Clifford circuits, a unique periodic orbit structure in the space of Pauli strings implies peculiar spectral correlations and level statistics with large degeneracies. $T$-gate doping induces an exponentially fast transition to chaotic behavior, described by random matrix theory. We compare these complexity indicators with magic generation properties of the Clifford+$T$ ensemble, and determine the distribution of magic, as well as the average non-stabilizing power of the quantum circuit ensemble. In the dilute limit, $N_T \ll N$, magic generation is governed by single-qubit behavior. Magic is generated in approximate quanta, increases approximately linearly with the number of $T$-gates, $N_T$, and displays a discrete distribution for small $N_T$. At $N_T\approx N$, the distribution becomes quasi-continuous, and for $N_T\gg N$ it converges to that of Haar-random unitaries, and averages to a finite magic density, $m_2$, $\lim_{N\to\infty} \langle m_2 \rangle_\text{Haar} = 1$. This is in contrast to the spectral transition, where ${\cal O} (1)$ $T$-gates suffice to remove spectral degeneracies and to induce a transition to chaotic behavior in the thermodynamic limit. Magic is therefore a more sensitive indicator of complexity., Comment: 13 pages, 9 figures

A magnetic impurity in a BCS superconductor induces the formation of a Shiba state and drives a local quantum phase transition. We generalize this concept to a one-dimensional superconductor with fractionalized excitations, where the dominant instability is superconducting. In this framework, conduction electrons fractionalize into gapless charge and gapped spin excitations. We show that magnetic impurity interacts exclusively with the spin degrees of freedom and induces a quantum phase transition. Furthermore, charge excitations influence dynamical observables, giving rise to the phenomenon we term the fractionalized Shiba state. At zero temperature, the tunneling spectrum exhibits universal power-law scaling with an exponent of $-1/2$ at half filling, stemming from the gapless charge modes that form a standard Luttinger liquid. Extending this analysis to finite temperatures reveals that the spectral features retain universal behavior at the critical point., Comment: 6 pages, 4 figures

A magnetic impurity with a larger $S=1$ spin remains partially screened by the Kondo effect when embedded in a metal. However, when placed within an $s$-wave superconductor, the interplay between the superconducting energy gap $\Delta$ and the Kondo temperature $T_K$ induces a quantum phase transition from an underscreened doublet Kondo to an unscreened triplet phase, typically occurring when $\Delta/T_K\approx 1$. We investigate the Kondo compensation of the impurity spin resulting from this partial screening across the quantum phase transition, which together with the spin-spin correlation function serves as a measure of the Kondo cloud's integrity. Deep within the unscreened triplet phase, $\Delta/T_K\gg 1$, the compensation vanishes, signifying complete decoupling of the impurity spin from the environment, while in the partially screened doublet phase, $\Delta/T_K\ll 1$, it asymptotically approaches $1/2$, indicating that half of the spin is screened. Notably, there is a universal jump in the compensation precisely at the phase transition, which we accurately calculate. The spin-spin correlation function exhibits an oscillatory pattern with an envelope function decaying as $\sim 1/x$ at short distances. At larger distances, the superconducting gap induces an exponentially decaying behavior $\sim \exp(-x/\xi_\Delta)$ governed by the superconducting correlation length $\xi_\Delta$, irrespective of the phase, without any distinctive features across the transition. Furthermore, the spectral functions of some relevant operators are evaluated and discussed. In terms of the methods used, a consistent description is provided through the application of multiplicative, numerical and density matrix renormalization group techniques., Comment: 9 pages, 7 figures

Stabilizer entropies and quantum magic have been extensively explored in real-space formulations of quantum systems within the framework of resource theory. However, interesting and transparent physics often emerges in momentum space, such as Cooper pairing. Motivated by this, we investigate the momentum-space structure of Pauli strings and stabilizer entropies in the one-dimensional transverse-field quantum Ising model. By mapping the Ising chain onto momentum-space qubits, where the stabilizer state corresponds to the paramagnetic state, we analyze the evolution of the Pauli string distribution. In the ferromagnetic phase, the distribution is broad, whereas in the paramagnetic phase, it develops a two-peaked structure. We demonstrate that all ferromagnetic states possess the same degree of magic in the thermodynamic limit, while stabilizer entropies are non-analytic at the critical point and vanish with increasing transverse field. The momentum-space approach to quantum magic not only complements its real-space counterpart but also provides advantages in terms of analyzing nonstabilizerness and classical simulability., Comment: 5 pages, 3 figures

We investigate the critical behavior of the Kondo compensation in the presence of a power-law pseudogap in the density of states, $\varrho(\omega)\sim |\omega|^\epsilon$. For $\epsilon<1$, this model exhibits a quantum phase transition from a partially screened doublet ground state to a fully screened many-body singlet ground state with increasing Kondo coupling, $j$. At the critical point, $j_c$, the Kondo compensation is found to scale as $\kappa(j

We propose a numerical method to efficiently compute quantum generating functions (QGF) for a wide class of observables in one-dimensional quantum systems at high temperature. We obtain high-accuracy estimates for the cumulants and reconstruct full counting statistics from the QGF. We demonstrate its potential on spin $S=1/2$ anisotropic Heisenberg chain, where we can reach time scales hitherto inaccessible to state-of-the-art classical and quantum simulations. Our results challenge the conjecture of the Kardar--Parisi--Zhang universality for isotropic integrable quantum spin chains., Comment: 7 pages, 3 figures plus Supporting Information

We investigate a specific limit of the one-dimensional non-Hermitian Hubbard Hamiltonian with complex interactions. In this framework, fermions with different spin quantum numbers are mapped onto two distinct spin species, resulting in two $XY$ spin chains that are coupled through Ising $ZZ$ interaction. The spin ladder model is then examined within the semiclassical limit using a spin-coherent state basis, where the dynamics is governed by a set of coupled Landau-Lifshitz-Gilbert equations. The non-Hermitian interactions in this model generate a spin-transfer torque term. We analyze the system's evolution toward several potential steady states, including a state of decoupled chains that is accessible when each chain has uniform initial conditions. Other possible steady states involve dimerized configurations with decoupled rungs, where the rung spins are either ferromagnetically or antiferromagnetically coupled, depending on the sign of the imaginary interactions. The spin dynamics is then studied in the infinite-temperature limit, which favors dimerized steady states. Despite the decoupled rungs, we observe the formation of ferromagnetic domains along each chain in the steady state. Additionally, we investigate the spin correlation functions and identify signatures of anomalous spin dynamics., Comment: 14 pages, 9 figures

Topological superconductors (SCs) hold great promise for fault-tolerant quantum hardware, however, their experimental realization is very challenging. Recently, superconducting artificial molecules (Andreev molecules) have opened new avenues to engineer topological superconducting materials. In this work, we demonstrate a heteroatomic Andreev molecule, where two normal artificial atoms realized by quantum dots (QDs) are coupled by a superconducting island (SCI). We show that the two normal atoms strongly hybridize and form a 3-electron-based molecular state. Our density matrix renormalization group (DMRG) calculations explain quantitatively the robust binding of electrons. The tunability of the structure allows us to drive a quantum phase transition from an antiferromagnetic Andreev molecular state to a heteroatomic Andreev molecule with ferromagnetically coupled QDs using simple electrical gating.

Information propagation in the one-dimensional infinite temperature Hubbard model with a dissipative particle sink at the end of a semi-infinite chain is studied. In the strongly interacting limit, the two-site mutual information and the operator entanglement entropy exhibit a rich structure with two propagating information fronts and superimposed interference fringes. A classical reversible cellular automaton model quantitatively captures the transport and the slow, classical part of the correlations, but fails to describe the rapidly propagating information jet. The fast quantum jet resembles coherent free particle propagation, with the accompanying long-ranged interference fringes that are exponentially damped by short-ranged spin correlations in the many-body background., Comment: 5 pages, 3 figures, supplemental material included

We focus on the biorthogonal work statistics of the interacting many-body Hatano-Nelson model after switching on the imaginary vector potential. We introduce a generalized Loschmidt echo $G(t)$ utilizing the biorthogonal metric operator. It is well suited for numerical analysis and its Fourier transform yields the probability distribution of work done. The statistics of work displays several universal features, including an exponential decay with the square of both the system size and imaginary vector potential for the probability to stay in the ground state. Additionally, its high energy tail follows a universal power law with exponent $-3$. This originates from the peculiar temporal power law decay of $G(t)$ with a time dependent exponent. The mean and the variance of work scale linearly and logarithmically with system size while all higher cumulants are non-extensive. Our results are relevant for non-unitary field theories as well., Comment: 7 pages, 3 figures

The collective tunneling of a Wigner necklace - a crystal-like state of a small number of strongly interacting electrons confined to a suspended nanotube and subject to a double well potential - is theoretically analyzed and compared with experiments in [Shapir \emph{et al.}, Science {\bf 364}, 870 (2019)]. Density Matrix Renormalization Group computations, exact diagonalization, and instanton theory provide a consistent description of this very strongly interacting system, and show good agreement with experiments. Experimentally extracted and theoretically computed tunneling amplitudes exhibit a scaling collapse. Collective quantum fluctuations renormalize the tunneling, and substantially enhance it as the number of electrons increases., Comment: 10 pages, 9 figures

We explore the Kardar-Parisi-Zhang (KPZ) scaling in the one-dimensional Hubbard model, which exhibits global $SU_c(2)\otimes SU_s(2)$ symmetry at half-filling, for the pseudo-charge and the total spin. We analyze dynamical scaling properties of high temperature charge and spin correlations and transport. At half-filling, we observe a clear KPZ scaling in both charge and spin sectors. Away from half-filling, the $SU_c(2)$ charge symmetry is reduced to $U_c(1)$, while the $SU_s(2)$ symmetry for the total spin is retained. Consequently, transport in the charge sector becomes ballistic, while KPZ scaling is preserved in the spin sector. These findings confirm the link between non-abelian symmetries and KPZ scaling in the presence of integrability. We study two settings of the model: one involving a quench from a bi-partitioned state asymptotically close to the $T=\infty$ equilibrium state of the system, and another where the system is coupled to two markovian reservoirs at the two edges of the chain., Comment: 14 pages, 18 figures

We investigate the decay of spatial correlations of $\mathcal{PT}$-symmetric non-Hermitian one-dimensional models that host higher-order exceptional points. Beyond a certain correlation length, they develop anomalous power-law behavior that indicates strong suppression of correlations in the non-Hermitian setups as compared to the Hermitian ones. The correlation length is also reflected in the entanglement entropy where it marks a change from logarithmic growth at short distance to a constant value at large distance, characteristic of an insulator, despite the spectrum being gapless. Two different families of models are investigated, both having a similar spectrum constrained by particle-hole symmetry. The first model offers an experimentally attractive way to generate arbitrary higher-order exceptional points and represents a non-Hermitian extension of the Dirac Hamiltonian for general spin. At the critical point it displays a decay of the correlations $\sim 1/x^2$ and $1/x^3$ irrespective of the order of the exceptional point. The second model is constructed using unidirectional hopping and displays enhanced suppression of correlations $\sim 1/x^a$, $a\ge 2$ with a power law that depends on the order of the exceptional point., Comment: 19 pages, 11 figures

We investigate the quantum quench dynamics of the interacting Hatano-Nelson model with open boundary conditions using both abelian bosonization and numerical methods. Specifically, we follow the evolution of the particle density and current profile in real space over time by turning the imaginary vector potential on or off in the presence of weak interactions. Our results reveal spatio-temporal Friedel oscillations in the system with light cones propagating ballistically from the open ends, accompanied by local currents of equal magnitude for both switch off and on protocols. Remarkably, the bosonization method accurately accounts for the density and current patterns with a single overall fitting parameter. The continuity equation is satisfied by the long wavelength part of the density and current, despite the non-unitary time evolution when the Hatano-Nelson term is switched on., Comment: 9 pages, 4 figures

We study the dissipative dynamics of a one-dimensional bosonic system described in terms of the bipartite Bose-Hubbard model with alternating gain and loss. This model exhibits the $\mathcal{PT}$ symmetry under some specific conditions and features a $\mathcal{PT}$-symmetry phase transition. It is characterized by an order parameter corresponding to the population imbalance between even and odd sites, similar to the continuous phase transitions in the Hermitian realm. In the noninteracting limit, we solve the problem exactly and compute the parameter dependence of the order parameter. The interacting limit is addressed at the mean-field level, which allows us to construct the phase diagram for the model. We find that both the interaction and dissipation rates induce a $\mathcal{PT}$-symmetry breaking. On the other hand, periodic modulation of the dissipative coupling in time stabilizes the $\mathcal{PT}$-symmetric regime. Our findings are corroborated numerically on a tight-binding chain with gain and loss., Comment: 8 pages, 7 figures

Recent experiments demonstrated that single-particle quantum walks can reveal the topological properties of single-particle states. Here, we generalize this picture to the many-body realm by focusing on multiparticle quantum walks of strongly interacting fermions. After injecting $N$ particles with multiple flavors in the interacting SU$(N)$ Su-Schrieffer-Heeger chain, their multiparticle continuous-time quantum walk is monitored by a variety of methods. We find that the many-body Berry phase in the $N$-body part of the spectrum signals a topological transition upon varying the dimerization, similarly to the single-particle case. This topological transition is captured by the single- and many-body mean chiral displacement during the quantum walk and remains present for strong interaction as well as for moderate disorder. Our predictions are well within experimental reach for cold atomic gases and can be used to detect the topological properties of many-body excitations through dynamical probes., Comment: 7 pages, 3 figures, and Supplemental Material

Particle formation represents a central theme in various branches of physics, often associated to confinement. Here we show that dynamical hadron formation can be spectroscopically detected in an ultracold atomic setting within the most paradigmatic and simplest model of condensed matter physics, the repulsive $\textrm{SU}(N)$ Hubbard model. By starting from an appropriately engineered initial state of the ${{\textrm{SU}(3)}}$ Hubbard model, not only mesons (doublons) but also baryons (trions) are naturally generated during the time evolution. In the strongly interacting limit, baryons become heavy and attract each other strongly, and their residual interaction with mesons generates meson diffusion, as captured by the evolution of the equal time density correlation function. Hadrons remain present in the long time limit, while the system thermalizes to a negative temperature state. Our conclusions extend to a large variety of initial conditions, all spatial dimensions, and for SU($N>2$) Hubbard models., Comment: 7 pages, 5 figures

We perform a comparative study for the magnetization dynamics within linear response theory of one and two dimensional massive Dirac electrons, after switching on either a real (hermitian) or an imaginary (non-hermitian) magnetic field. While hermitian dc magnetic fields polarize the spins in the direction of the external magnetic field, non-hermitian magnetic fields induce only off diagonal response. An imaginary dc magnetic field perpendicular to the mass term induces finite magnetization in the third direction only according to the right hand rule. This can be understood by analyzing the non-hermitian equation of motion of the spin, which becomes analogous to a classical particle in crossed electric and magnetic fields. Therein, the spin expectation value, the mass term and imaginary magnetic field play the role of the classical momentum, magnetic and electric field, respectively. The latter two create a drift velocity perpendicular to them, which gives rise to the off-diagonal component of the dc spin susceptibility, similarly to how the Hall effect develops in the classical description., Comment: 8 pages, 2 figures

We develop a non-Abelian time evolving block decimation (NA-TEBD) approach to study of open systems governed by Lindbladian time evolution, while exploiting an arbitrary number of abelian or non-abelian symmetries. We illustrate this method in a one-dimensional fermionic $SU(2)$ Hubbard model on a semi-infinite lattice with localized particle loss at one end. We observe a ballistic front propagation with strongly renormalized front velocity, and a hydrodynamic current density profile. For large loss rates, a suppression of the particle current is observed, as a result of the quantum Zeno effect. Operator entanglement is found to propagate faster than the depletion profile, preceding the latter., Comment: 11 pages, 8 figures

We consider a PT-symmetric Fermi gas with an exceptional point, representing the critical point between PT-symmetric and symmetry broken phases. The low energy spectrum remains linear in momentum and is identical to that of a hermitian Fermi gas. The fermionic Green's function decays in a power law fashion for large distances, as expected from gapless excitations, albeit the exponent is reduced from $-1$ due to the quantum Zeno effect. In spite of the gapless nature of the excitations, the ground state entanglement entropy saturates to a finite value, independent of the subsystem size due to the non-hermitian correlation length intrinsic to the system. Attractive or repulsive interaction drives the system into the PT-symmetry broken regime or opens up a gap and protects PT-symmetry, respectively. Our results challenge the concept of universality in non-hermitian systems, where quantum criticality can be masked due to non-hermiticity., Comment: 8 pages, 5 figures

The Lindhard function represents the basic building block of many-body physics and accounts for charge response, plasmons, screening, Friedel oscillation, RKKY interaction etc. Here we study its non-Hermitian version in one dimension, where quantum effects are traditionally enhanced due to spatial confinement, and analyze its behavior in various limits of interest. Most importantly, we find that the static limit of the non-Hermitian Lindhard function has no divergence at twice the Fermi wavenumber and vanishes identically for all other wavenumbers at zero temperature. Consequently, no Friedel oscillations are induced by a non-Hermitian, imaginary impurity to lowest order in the impurity potential at zero temperature. Our findings are corroborated numerically on a tight-binding ring by switching on a weak real or imaginary potential. We identify conventional Friedel oscillations or heavily suppressed density response, respectively., Comment: 6 pages, 3 figures

Linear response theory plays a prominent role in various fields of physics and provides us with extensive information about the thermodynamics and dynamics of quantum and classical systems. Here we develop a general theory for the linear response in non-Hermitian systems with non-unitary dynamics and derive a modified Kubo formula for the generalized susceptibility for arbitrary (Hermitian and non-Hermitian) system and perturbation. As an application, we evaluate the dynamical response of a non-Hermitian, one-dimensional Dirac model with imaginary and real masses, perturbed by a time-dependent electric field. The model has a rich phase diagram, and in particular, features a tachyon phase, where excitations travel faster than an effective speed of light. Surprisingly, we find that the dc conductivity of tachyons is finite, and the optical sum rule is exactly satisfied for all masses. Our results highlight the peculiar properties of the Kubo formula for non-Hermitian systems and are applicable for a large variety of settings.

Magnetic impurities embedded in a metal are screened by the Kondo effect, signaled by the formation of an extended correlation cloud, the so-called Kondo or screening cloud. In a superconductor, the Kondo state turns into sub-gap Yu-Shiba-Rusinov (Shiba) states, and a quantum phase transition occurs between screened and unscreened phases once the superconducting energy gap $\Delta$ becomes sufficiently large compared to the Kondo temperature, $T_K$. Here we show that, although the Kondo state does not form in the unscreened phase, the Kondo cloud does exist in both quantum phases. However, while screening is complete in the screened phase, it is only partial in the unscreened phase. Compensation, a quantity introduced to characterize the integrity of the cloud, is universal, and shown to be related to the magnetic impurities' $g$-factor, monitored experimentally by bias spectroscopy., Comment: 5 pages, 4 figures, 3 pages of supplementary material

We study the non-equilibrium dynamics and transport of a PT-symmetric Luttinger liquid (LL) after an interaction quench. The system is prepared in domain wall initial state. After a quantum quench to spatially homogeneous, PT-symmetric LL, the domain wall develops into a flat central region that spreads out ballistically faster than the conventional Lieb-Robinson maximal speed. By evaluating the current inside the regular lightcone, we find a universal conductance $e^2/h$, insensitive to the strength of the PT-symmetric interaction. On the other hand, by repeating the very same time evolution with a hermitian LL Hamiltonian, the conductance is heavily renormalized by the hermitian interaction as $e^2/hK$ with $K$ the LL parameter. Our analytical results are tested numerically, confirming the universality of the conductance in the non-hermitian realm., Comment: 5 pages, 2 figures, 2 pages of supplementary material

We study the Coulomb glass emerging from the interplay of strong interactions and disorder in a model of spinless fermions on the Bethe lattice. In the infinite coordination number limit, strong interactions induce a metallic Coulomb glass phase with a pseudogap structure at the Fermi energy. Quantum and thermal fluctuations both melt this glass and induce a disordered quantum liquid phase. We combine self-consistent diagrammatic perturbation theory with continuous time quantum Monte-Carlo simulations to obtain the complete phase diagram of the electron glass, and to characterize its dynamical properties in the quantum liquid, as well as in the replica symmetry broken glassy phase. Tunneling spectra display an Efros-Shklovskii pseudogap upon decreasing temperatures, but the density of states remains finite at the Fermi energy due to residual quantum fluctuations. Our results bear relevance to the metallic glass phase observed in Si inversion layers., Comment: 16 pages, 13 figures, accepted version

We introduce the notion of non-Abelian tensors, and use them to construct a general non-Abelian time evolving block decimation (NA-TEBD) scheme that uses an arbitrary number of Abelian and non-Abelian symmetries. Our approach increases the speed and memory storage efficiency of matrix product state based computations by several orders of magnitudes, and makes large bond dimensions accessible even on simple desktop architectures. We use it to study post-quench dynamics in the repulsive SU(3) Hubbard model, and to determine the time evolution of various local operators and correlation functions efficiently. Interactions turn algebraic charge relaxation into exponential, and suppress coherent quantum oscillations rapidly., Comment: 10 + 5 pages, 15 figures

We investigate the stability of a Luttinger liquid, upon suddenly coupling it to a dissipative environment. Within the Lindblad equation, the environment couples to local currents and heats the quantum liquid up to infinite temperatures. The single particle density matrix reveals the fractionalization of fermionic excitations in the spatial correlations by retaining the initial non-integer power law exponents, accompanied by an exponential decay in time with interaction dependent rate. The spectrum of the time evolved density matrix is gapped, which collapses gradually as $-\ln(t)$. The von Neumann entropy crosses over from the early time $-t\ln(t)$ behaviour to $\ln(t)$ growth for late times. The early time dynamics is captured numerically by performing simulations on spinless interacting fermions, using several numerically exact methods. Our results could be tested experimentally in bosonic Luttinger liquids., Comment: 8 pages, 3 figures

Semiconductors in the proximity of superconductors have been proposed to support phases hosting Majorana bound states. When the systems undergo a topological phase transition towards the Majorana phase, the spectral gap closes, then reopens, and the quasiparticle band spin polarization is inverted. We focus on two paradigmatic semiconductor-superconductor heterostructures and propose an all-electrical spectroscopic probe sensitive to the spin inversion at the topological transition. Our proposal relies on the indirect coupling of a time-dependent electric field to the electronic spin due to the strong Rashba spin-orbit coupling in the semiconductor. We analyze within linear response theory the dynamical correlation functions and demonstrate that some components of the susceptibility can be used to detect the nontrivial topological phases.

We study a one-dimensional Fermi gas in the presence of dissipative coupling to environment through the Lindblad equation. The dissipation involves energy exchange with the environment and favours the relaxation of electrons to excitations. After switching on the dissipation, the system approaches a steady state, which is described by a generalized Gibbs ensemble. The fermionic single particle density matrix resembles deceivingly to that in a hermitian interaction quench. It decays inversely with the distance for short times due to the fermionic correlations in the initial state, which changes into a non-integer power law decay for late times, representing dissipation induced Luttinger liquid behaviour. However, the crossover between the two regions occurs due to dissipation induced damping, and is unrelated to the propagation of excitations. The velocity of information spreading is set by the dissipative coupling, and differs significantly from the original sound velocity. The thermodynamic entropy grows as $-t\ln t$ initially, and saturates to an extensive value. Our results can be tested experimentally in one-dimensional Dirac systems., Comment: 5 pages, 2 figures + supplementary material

For most chiralities, semiconducting nanotubes display topologically protected end states of multiple degeneracies. We demonstrate using density matrix renormalization group based quantum chemistry tools that the presence of Coulomb interactions induces the formation of robust end spins. These are the close analogues of ferromagnetic edge states emerging in graphene nanoribbons. The interaction between the two ends is sensitive to the length of the nanotube, its dielectric constant, as well as the size of the end spins: for $S=1/2$ end spins their interaction is antiferromagnetic, while for $S>1/2$ it changes from antiferromagnetic to ferromagnetic as the nanotube length increases. The interaction between end spins can be controlled by changing the dielectric constant of the environment, thereby providing a possible platform for two-spin quantum manipulations., Comment: 5 pages, 4 figures + supplementary material

Various promising qubit concepts have been put forward recently based on engineered superconductor (SC) subgap states like Andreev bound states, Majorana zero modes or the Yu-Shiba-Rusinov (Shiba) states. The coupling of these subgap states via a SC strongly depends on their spatial extension and is an essential next step for future quantum technologies. Here we investigate the spatial extension of a Shiba state in a semiconductor quantum dot coupled to a SC for the first time. With detailed transport measurements and numerical renormalization group calculations we find a remarkable more than 50 nm extension of the zero energy Shiba state, much larger than the one observed in very recent scanning tunneling microscopy (STM) measurements. Moreover, we demonstrate that its spatial extension increases substantially in magnetic field., Comment: 11 pages, 7 figures

A Luttinger liquid (LL) describes low energy excitations of many interacting one dimensional systems, and exhibits universal response both in and out of equilibrium. We analyze its behaviour in the non-hermitian realm after quantum quenching to a PT-symmetric LL by focusing on the fermionic single particle density matrix. For short times, we demonstrate the emergence of unique phenomena, characteristic to non-hermitian systems, that correlations propagate faster than the conventional maximal speed, known as the Lieb-Robinson bound. These emergent supersonic modes travel with velocities that are multiples of the conventional light-cone velocity. This behaviour is argued to be generic for correlators in non-hermitian systems. In the long time limit, we find typical LL behaviour, extending the LL universality to the non-equilibrium non-hermitian case. Our analytical results are benchmarked numerically and indicate that the dispersal of quantum information is much faster in non-hermitian systems., Comment: 6+1 pages, 3 figures

While much is known about repulsive quantum impurity models, significantly less attention has been devoted to their attractive counterparts. This motivated us to study the attractive SU(N) Anderson impurity model. While for the repulsive case, the phase diagram features mild N dependence and the ground state is always a Fermi liquid, in the attractive case a Kosterlitz-Thouless charge localization phase transition is revealed for N > 2. Beyond a critical value of attractive interaction an abrupt jump appears in the number of particles at the impurity site, and a singular Fermi liquid state emerges, where the scattering of quasiparticles is found to exhibit power law behavior with fractional power. The capacity diverges exponentially at the quantum critical point, signaling the Kosterlitz-Thouless transition., Comment: 5 pages, 4 figures, includes Supplemental Information

We study quantum quenches in the $S=1$ Heisenberg spin chain and show that the dynamics can be described by the recently developed semi-semiclassical method based on particles propagating along classical trajectories but scattering quantum mechanically. We analyze the non-equilibrium time evolution of the distribution of the total spin in half of the system and compare the predictions of the semi-semiclassical theory with those of a non-Abelian time evolving block decimation (TEBD) algorithm which exploits the SU(2) symmetry. We show that while the standard semiclassical approach using the universal low energy scattering matrix cannot describe the dynamics, the hybrid semiclassical method based on the full scattering matrix gives excellent agreement with the first principles TEBD simulation., Comment: 14 pages, 11 figures

Motivated by recent experiments in ultracold gases, we focus on the properties of the center of mass coordinate of an interacting one dimensional Fermi gas, displaying several distinct phases. While the variance of the center of mass vanishes in insulating phases such as phase separated and charge density wave phases, it remains finite in the metallic phase, which realizes a Luttinger liquid. By combining numerics with bosonization, we demonstrate that the autocorrelation function of the center of mass coordinate is universal throughout the metallic phase. It exhibits persistent oscillations and its short time dynamics reveal important features of the quantum liquid, such as the Luttinger liquid parameter and the renormalized velocity. The full counting statistics of the center of mass follows a normal distribution already for small systems. Our results apply to non-integrable systems as well and are within experimental reach for e.g. carbon nanotubes and cold atomic gases., Comment: 6 pages, 3 figures

The quantum crystal of electrons, predicted more than eighty years ago by Eugene Wigner, is still one of the most elusive states of matter. Here, we present experiments that observe the one-dimensional Wigner crystal directly, by imaging its charge density in real-space. To measure this fragile state without perturbing it, we developed a new scanning probe platform that utilizes a pristine carbon nanotube as a scanning charge perturbation to image, with minimal invasiveness, the many-body electronic density within another nanotube. The obtained images, of few electrons confined in one-dimension, match those of strongly interacting crystals, with electrons ordered like pearls on a necklace. Comparison to theoretical modeling demonstrates the dominance of Coulomb interactions over kinetic energy and the weakness of exchange interactions. Our experiments provide direct evidence for this long-sought electronic state, and open the way for studying other fragile interacting states by imaging their many-body density in real-space.

We propose a scheme to realize the Kondo model with tunable anisotropy using alkaline-earth atoms in an optical lattice. The new feature of our setup is Floquet engineering of interactions using time-dependent Zeeman shifts, that can be realized either using state-dependent optical Stark shifts or magnetic fields. The properties of the resulting Kondo model strongly depend on the anisotropy of the ferromagnetic interactions. In particular, easy-plane couplings give rise to Kondo singlet formation even though microscopic interactions are all ferromagnetic. We discuss both equilibrium and dynamical properties of the system that can be measured with ultracold atoms, including the impurity spin susceptibility, the impurity spin relaxation rate, as well as the equilibrium and dynamical spin correlations between the impurity and the ferromagnetic bath atoms. We analyze the non-equilibrium time evolution of the system using a variational non-Gaussian approach, which allows us to explore coherent dynamics over both short and long timescales, as set by the bandwidth and the Kondo singlet formation, respectively. In the quench-type experiments, when the Kondo interaction is suddenly switched on, we find that real-time dynamics shows crossovers reminiscent of poor man's renormalization group flow used to describe equilibrium systems. For bare easy-plane ferromagnetic couplings, this allows us to follow the formation of the Kondo screening cloud as the dynamics crosses over from ferromagnetic to antiferromagnetic behavior. On the other side of the phase diagram, our scheme makes it possible to measure quantum corrections to the well-known Korringa law describing the temperature dependence of the impurity spin relaxation rate. Theoretical results discussed in our paper can be measured using currently available experimental techniques., Comment: 22 pages, 12 figures

We use a semiclassical approach to study out of equilibrium dynamics and transport in quantum systems with massive quasiparticle excitations having internal quantum numbers. In the universal limit of low energy quasiparticles, the system is described in terms of a classical gas of colored hard-core particles. Starting from an inhomogeneous initial state, in this limit we give analytic expressions for the space and time dependent spin density and spin current profiles. Depending on the initial state, the spin transport is found to be ballistic or diffusive. In the ballistic case we identify a `second front' that moves more slowly than the maximal quasiparticle velocity. Our analytic results also capture the diffusive broadening of this ballistically propagating front. To go beyond the universal limit, we study the effect of non-trivial scattering processes in the $O(3)$ non-linear sigma model by performing Monte Carlo simulations, and observe local equilibration around the second front in terms of the densities of the particle species., Comment: published version

We investigate the spin-resolved transport properties, such as the linear conductance and the tunnel magnetoresistance, of a double quantum dot device attached to ferromagnetic leads and look for signatures of SU(4) symmetry in the Kondo regime. We show that the transport behavior greatly depends on the magnetic configuration of the device, and the spin-SU(2) as well as the orbital and spin-SU(4) Kondo effects become generally suppressed when the magnetic configuration of the leads varies from the antiparallel to the parallel one. Furthermore, a finite spin polarization of the leads lifts the spin degeneracy and drives the system from the SU(4) to an orbital-SU(2) Kondo state. We analyze in detail the crossover and show that the Kondo temperature between the two fixed points has a non-monotonic dependence on the degree of spin polarization of the leads. In terms of methods used, we characterize transport by using a combination of analytical and numerical renormalization group approaches., Comment: 13 pages, 11 figures

The two-channel Kondo impurity model realizes a local non-Fermi liquid state with finite residual entropy. The competition between the two channels drives the system to an impurity quantum critical point. We show that the out-of-time-ordered (OTO) commutator for the impurity spin reveals markedly distinct behaviour depending on the low energy impurity state. For the one channel Kondo model with Fermi liquid ground state, the OTO commutator vanishes for late times, indicating the absence of the butterfly effect. For the two channel case, the impurity OTO commutator is completely temperature independent and saturates quickly to its upper bound 1/4, and the butterfly effect is maximally enhanced. These compare favourably to numerics on spin chain representation of the Kondo model. Our results imply that a large late time value of the OTO commutator does not necessarily diagnose quantum chaos., Comment: 8 pages, 6 figures

We construct a Fermi liquid theory to describe transport in a superconductor-quantum dot- normal metal junction close to the singlet-doublet (parity changing) transition of the dot. Though quasiparticles do not have a definite charge in this chargeless Fermi liquid, in case of particle-hole symmetry, a mapping to the Anderson model unveils a hidden U(1) symmetry and a corresponding pseudo-charge. In contrast to other correlated Fermi-liquids, the back scattering noise reveals an effective charge equal to the charge of Cooper pairs, $e^* = 2e$. In addition,we find a strong suppression of noise when the linear conductance is unitary, even for its non-linear part., Comment: 5 pages, 4 figures and 5 pages of supplementary material with 3 figures

We present a detailed theoretical analysis of the Wigner crystal states in confined semiconducting carbon nanotubes. We show by robust scaling arguments as well as by detailed semi-microscopic calculations that the effective exchange interaction has an SU(4) symmetry, and can reach values even as large as $J\sim 100 {\rm \,K}$ in weakly screened, small diameter nanotubes, close to the Wigner crystal - electron liquid crossover. Modeling the nanotube carefully and analyzing the magnetic structure of the inhomogeneous electron crystal, we recover the experimentally observed 'phase boundaries' of Deshpande and Bockrath [V. V. Deshpande and M. Bockrath, Nature Physics $\mathbf 4$, 314 (2008)]. Spin-orbit coupling only slightly modifies these phase boundaries, but breaks the spin symmetry down to SU(2)$\times$SU(2), and in Wigner molecules it gives rise to interesting excitation spectra, reflecting the underlying SU(4) as well as the residual SU(2)$\times$SU(2) symmetries., Comment: 14 pages, 13 figures

This paper is a corrected version of Phys. Rev. B 95, 165404 (2017), which we have retracted because it contained a trivial but fatal sign error that lead to incorrect conclusions. --- We extend a recently-eveloped Fermi-liquid (FL) theory for the asymmetric single-impurity Anderson model [C. Mora $et al.$, Phys. Rev. B, 92, 075120 (2015)] to the case of an arbitrary local magnetic field. To describe the system's low-lying quasiparticle excitations for arbitrary values of the bare Hamiltonian's model parameters, we construct an effective low-energy FL Hamiltonian whose FL parameters are expressed in terms of the local level's spin-dependent ground-state occupations and their derivatives with respect to level energy and local magnetic field. These quantities are calculable with excellent accuracy from the Bethe Ansatz solution of the Anderson model. Applying this effective model to a quantum dot in a nonequilibrium setting, we obtain exact results for the curvature of the spectral function, $c_A$, describing its leading $\sim\varepsilon^2$ term, and the transport coefficients $c_V$ and $c_T$, describing the leading $\sim V^2$ and $\sim T^2$ terms in the nonlinear differential conductance. A sign change in $c_A$ or $c_V$ is indicative of a change from a local maximum to a local minimum in the spectral function or nonlinear conductance, respectively, as is expected to occur when an increasing magnetic field causes the Kondo resonance to split into two subpeaks. We find that the fields $B_A$, $B_T$ and $B_V$ at which $c_A$, $c_T$ and $c_V$ change sign, respectively, are all of order $T_K$, as expected, with $B_A = B_T = B_V = 0.75073\,T_K$ in the Kondo limit., Comment: 21 pages, 18 figures. Version 1 of this paper, published as Phys. Rev. B 95, 165404 (2017), contained a trivial but fatal sign error that lead to incorrect physical conclusions; it has therefore been retracted. Version 3 is a corrected version with major revisions, additional NRG data and different conclusions; it has been accepted for publication in PRB

We develop a hybrid semiclassical method to study the time evolution of one dimensional quantum systems in and out of equilibrium. Our method handles internal degrees of freedom completely quantum mechanically by a modified time evolving block decimation method, while treating orbital quasiparticle motion classically. We can follow dynamics up to timescales well beyond the reach of standard numerical methods to observe the crossover between pre-equilibrated and locally phase equilibrated states. As an application, we investigate the quench dynamics and phase fluctuations of a pair of tunnel coupled one dimensional Bose condensates. We demonstrate the emergence of soliton-collision induced phase propagation, soliton-entropy production and multistep thermalization. Our method can be applied to a wide range of gapped one-dimensional systems., Comment: published version

When a strongly correlated quantum dot is tunnel-coupled to a superconductor, it leads to the formation of Shiba bound states inside the superconducting gap. They have been measured experimentally in a superconductor-quantum dot-normal lead setup. Side coupling the quantum dot to a topological superconducting wire that supports Majorana bound states at its ends, drastically affects the structure of the Shiba states and induces supplementary in-gap states. The anomalous coupling between the Majorana bound states and the quantum dot gives rise to a characteristic imbalance in the spin resolved spectral functions for the dot operators. These are clear fingerprints for the existence of Majorana fermions and they can be detected experimentally in transport measurements. In terms of methods employed, we have used analytical approaches combined with the numerical renormalization group approach., Comment: 5 pages, 3 figures

We consider a strongly correlated quantum dot, tunnel coupled to two superconducting leads and capacitively coupled to a single mode microwave cavity. When the superconducting gap is the largest energy scale, multiple Shiba states are formed inside the gap. The competition of these states for the ground state signals a quantum phase transition. We demonstrate that photonic measurements can be used to probe such localized Shiba states. Moreover, the quantum phase transition can be pinpointed exactly from the sudden change in the transmission signal. Calculations were performed using the numerical renormalization group approach., Comment: 9 pages, 8 figures

We investigate the effects of spin-orbit interaction (SOI) and valley mixing on the transport and dynamical properties of a carbon nanotube (CNT) quantum dot in the Kondo regime. As these perturbations break the pseudo-spin symmetry in the CNT spectrum but preserve time-reversal symmetry, they induce a finite splitting $\Delta$ between formerly degenerate Kramers pairs. Correspondingly, a crossover from the SU(4) to the SU(2)-Kondo effect occurs as the strength of these symmetry breaking parameters is varied. Clear signatures of the crossover are discussed both at the level of the spectral function as well as of the conductance. In particular, we demonstrate numerically and support with scaling arguments, that the Kondo temperature scales inversely with the splitting $\Delta$ in the crossover regime. In presence of a finite magnetic field, time reversal symmetry is also broken. We investigate the effects of both parallel and perpendicular fields (with respect to the tube's axis), and discuss the conditions under which Kondo revivals may be achieved., Comment: 13 pages, 17 figures

We generalize Nozi\`eres' Fermi-liquid theory for the low-energy behavior of the Kondo model to that of the single-impurity Anderson model. In addition to the electrons' phase shift at the Fermi energy, the low-energy Fermi-liquid theory is characterized by four Fermi-liquid parameters: the two given by Nozi\`eres that enter to first order in the excitation energy, and two additional ones that enter to second order and are needed away from particle-hole symmetry. We express all four parameters in terms of zero-temperature physical observables, namely the local charge and spin susceptibilities and their derivatives with respect to the local level position. We determine these in terms of the bare parameters of the Anderson model using Bethe Ansatz and Numerical Renormalization Group (NRG) calculations. Our low-energy Fermi-liquid theory applies throughout the crossover from the strong-coupling Kondo regime via the mixed-valence regime to the empty-orbital regime. From the Fermi-liquid theory, we determine the conductance through a quantum dot symmetrically coupled to two leads in the regime of small magnetic field, low temperature and small bias voltage, and compute the coefficients of the $\sim B^2$, $\sim T^2$, and $\sim V^2$ terms \textit{exactly} in terms of the Fermi-liquid parameters. The coefficients of $T^2$, $V^2$ and $B^2$ are found to change sign during the Kondo to empty-orbital crossover. The crossover becomes universal in the limit that the local interaction is much larger than the level width. For completeness, we also compute the shot noise and discuss the resulting Fano factor., Comment: 17 pages, 10 figures

We investigate the equilibrium and non-equilibrium transport through a quantum dot in the Kondo regime, embedded between a normal metal and a topological superconductor supporting Majorana bound states at its end points. We find that the Kondo physics is significantly modified by the presence of the Majorana modes. When the Majorana modes are coupled, aside from the Kondo scale $T_K$, a new energy scale $T^*\ll T_K$ emerges, that controls the low energy physics of the system. At low temperatures, the ac-conductance is suppressed for frequencies below $T^*$, while the noise spectrum acquires a $\sim \omega^3$ dependence. At high temperatures, $T \gg T_K$, the regular logarithmic dependence in the differential conductance is also affected. Under non-equilibrium conditions, and in particular in the $\{T, B\}\to 0$ limit, the differential conductance becomes negative. These findings indicate that the changes in transport may serve as clues for detecting the Majorana bound states in such systems. In terms of methods used, we characterize the transport by using a combination of perturbative and renormalization group approaches., Comment: 14 pages, 11 figures