Search

In this short review, we provide an overview of recent progress in deploying advanced characterization techniques to understand the effects of local inhomogeneities in moir\'e heterostructures over multiple length scales. Particular emphasis is placed on correlating the impact of twist angle misalignment, nano-scale disorder, and atomic relaxation on the moir\'e potential and its collective excitations, particularly moir\'e excitons. Finally, we discuss future technological applications leveraging based on moi\'e excitons., Comment: Submitted to Small in response to an invitation for the special issue that celebrates the scientific accomplishments of Prof. P M Ajayan

We recapitulate recent developments of the functional renormalization group (FRG) approach to the steady state of systems out of thermal equilibrium. In particular, we discuss second-order truncation schemes which account for the frequency-dependence of the two particle vertex and which incorporate inelastic processes. Our focus is on two different types of one-dimensional fermion chains: i) infinite, open systems which feature a translation symmetry, and ii) finite systems coupled to left and right reservoirs. In addition to giving a detailed and unified review of the technical derivation of the FRG schemes, we briefly summarize some of the key physical results. In particular, we compute the non-equilibrium phase diagram and analyze the fate of the Berezinskii-Kosterlitz-Thouless transition in the infinite, open system., Comment: Appearing as part of a EPJB collection "Developments in the functional renormalization group approach to correlated electron systems"

Using the time-dependent Ginzburg Landau equations we study vortex motion driven by an applied current in two dimensional superconductors in the presence of a physical boundary. At smaller sourced currents the vortex lattice moves as a whole, with each vortex moving at the same velocity. At larger sourced current, vortex motion is organized into channels, with vortices in channels nearer to the sample edges moving faster than those farther away from sample edges, opposite to the Poiseuille flow of basic hydrodynamics where the velocity is lowest at the boundaries. At intermediate currents, a stick-slip motion of the vortex lattice occurs in which vortices in the channel at the boundary break free from the Abrikosov lattice, accelerate, move past their neighbors and then slow down and reattach to the vortex lattice at which point the stick-slip process starts over. These effects could be observed experimentally, e.g. using fast scanning microscopy techniques.

In superconductors, vortices induced by a magnetic field are nucleated randomly due to some random fluctuations or pinned by impurities or boundaries, impeding the development of vortex based quantum devices. Here, we propose a superconducting structure which allows to nucleate and control vortices on-demand by controlling magnetic fields and currents. Using time-dependent Ginzburg Landau theory, we study a driven vortex motion in two-dimensional Corbino geometries of superconductor-normal metal-superconductor Josephson junctions. We remedy the randomness of nucleation by introducing normal conducting rails to the Corbino disk to guide the nucleation process and motion of vortices towards the junction. We elaborate on the consequences of rail-vortex and vortex-vortex interactions to the quantization of resistance across the junction. Finally, we simulate the nucleations and manipulations of two and four vortices in Corbino networks, and discuss its application to Majorana zero mode braiding operations. Our study provides a potential route towards quantum computation with non-Abelian anyons.

We review recent progress in utilizing ultrafast light-matter interaction to control the macroscopic properties of quantum materials. Particular emphasis is placed on photoinduced phenomena that do not result from ultrafast heating effects but rather emerge from microscopic processes that are inherently nonthermal in nature. Many of these processes can be described as transient modifications to the free-energy landscape resulting from the redistribution of quasiparticle populations, the dynamical modification of coupling strengths and the resonant driving of the crystal lattice. Other pathways result from the coherent dressing of a material's quantum states by the light field. We discuss a selection of recently discovered effects leveraging these mechanisms, as well as the technological advances that led to their discovery. A road map for how the field can harness these nonthermal pathways to create new functionalities is presented., Comment: 36 pages, 12 figures; all authors contributed equally to this work

We study the nonlinear steady-state transport of spinless fermions through a quantum dot with a local two-particle interaction. The dot degree of freedom is in addition coupled to a phonon mode. This setup combines the nonequilibrium physics of the interacting resonant level model and that of the Anderson-Holstein model. The fermion-fermion interaction defies a perturbative treatment. We mainly focus on the antiadiabatic limit, with the phonon frequency being larger than the lead-dot tunneling rate. In this regime also the fermion-boson coupling cannot be treated perturbatively. Our goal is two-fold. We investigate the competing roles of the fermion-fermion and fermion-boson interactions on the emergent low-energy scale $T_{\rm K}$ and show how $T_{\rm K}$ manifests in the transport coefficients as well as the current-voltage characteristics. For small to intermediate interactions, the latter is in addition directly affected by both interactions independently. With increasing fermion-boson interaction the Franck-Condon blockade suppresses the current at small voltages and the emission of phonons leads to shoulders or steps at multiples of the phonon frequency, while the local fermion-fermion interaction implies a negative differential conductance at voltages larger than $T_{\rm K}$. We, in addition, use the model to investigate the limitations of our low-order truncated functional renormalization group approach on the Keldysh contour. In particular, we quantify the role of the broken current conservation., Comment: 13 pages, 9 figures

Quantum geometry has been identified as an important ingredient for the physics of quantum materials and especially of flat-band systems, such as moir\'e materials. On the other hand, the coupling between light and matter is of key importance across disciplines and especially for Floquet and cavity engineering of solids. Here we present fundamental relations between light-matter coupling and quantum geometry of Bloch wave functions, with a particular focus on flat-band and moir\'e materials, in which the quenching of the electronic kinetic energy could allow one to reach the limit of strong light-matter coupling more easily than in highly dispersive systems. We show that, despite the fact that flat bands have vanishing band velocities and curvatures, light couples to them via geometric contributions. Specifically, the intra-band quantum metric allows diamagnetic coupling inside a flat band; the inter-band Berry connection governs dipole matrix elements between flat and dispersive bands. We illustrate these effects in two representative model systems: (i) a sawtooth quantum chain with a single flat band, and (ii) a tight-binding model for twisted bilayer graphene. For (i) we highlight the importance of quantum geometry by demonstrating a nonvanishing diamagnetic light-matter coupling inside the flat band. For (ii) we explore the twist-angle dependence of various light-matter coupling matrix elements. Furthermore, at the magic angle corresponding to almost flat bands, we show a Floquet-topological gap opening under irradiation with circularly polarized light despite the nearly vanishing Fermi velocity. We discuss how these findings provide fundamental design principles and tools for light-matter-coupling-based control of emergent electronic properties in flat-band and moir\'e materials., Comment: 19 pages, 7 figures, including appendix; updated version: typos corrected and references added

The field of quantum simulations in ultra-cold atomic gases has been remarkably successful. In principle it allows for an exact treatment of a variety of highly relevant lattice models and their emergent phases of matter. But so far there is a lack in the theoretical literature concerning the systematic study of the effects of the trap potential as well as the finite size of the systems, as numerical studies of such non periodic, correlated fermionic lattices models are numerically demanding beyond one dimension. We use the recently introduced real-space truncated unity functional renormalization group to study these boundary and trap effects with a focus on their impact on the superconducting phase of the 2D Hubbard model. We find that in the experiments not only lower temperatures need to be reached compared to current capabilities, but also system size and trap potential shape play a crucial role to simulate emergent phases of matter., Comment: 21 pages, 9 Figures

We present a quantitative analysis of two-particle interaction effects in generalized, one-dimensional Aubry-Andr\'e-Harper models with the Fermi energy placed in one of the band gaps. We investigate systems with periodic as well as open boundary conditions; for the latter focusing on the number of edge states and the boundary charge. Both these observables are important for the classification of noninteracting topological systems. In our first class of models the unit cell structure stems from periodically modulated single-particle parameters. In the second it results from the spatial modulation of the two-particle interaction. For both types of models, we find that the single-particle band gaps are renormalized by the interaction in accordance with expectations employing general field theoretical arguments. While interaction induced effective edge states can be found in the local single-particle spectral function close to a boundary, the characteristics of the boundary charge are not modified by the interaction. This indicates that our results for the Rice-Mele and Su-Schriefer-Heeger model [Phys. Rev. B 102, 085122 (2020)] are generic and can be found in lattice models with more complex unit cells as well.

Recent optical conductivity measurements reveal the presence of Hubbard excitons in certain Mott insulators. In light of these results, it is important to revisit the dynamics of these materials to account for excitonic correlations. We investigate time-resolved excitation and relaxation dynamics as a function of temperature in perovskite-type LaVO${}_{3}$ thin films using ultrafast optical pump-probe spectroscopy. LaVO${}_{3}$ undergoes a series of phase transitions at roughly the same critical temperature $T_C\cong 140\ K$, including a second-order magnetic phase transition (PM $\xrightarrow{}$ AFM) and a first-order structural phase transition, accompanied by \textit{C}-type spin order (SO) and \textit{G}-type orbital order (OO). Ultrafast optical pump-probe spectroscopy at 1.6 eV monitors changes in the spectral weight of the Hubbard exciton resonance which serves as a sensitive reporter of spin and orbital fluctuation dynamics. We observe dramatic slowing down of the spin, and orbital dynamics in the vicinity of $T_C\cong 140$ K, reminiscent of a second-order phase transition, despite the (weakly) first-order nature of the transition. We emphasize that since it is spectral weight changes that are probed, the measured dynamics are not reflective of conventional exciton generation and recombination, but are related to the dynamics of Hubbard exciton formation in the presence of a fluctuating many-body environment.

The compelling original idea of a time crystal has referred to a structure that repeats in time as well as in space, an idea that has attracted significant interest recently. While obstructions to realize such structures became apparent early on, focus has shifted to seeing a symmetry breaking in time in periodically driven systems, a property of systems referred to as discrete time crystals. In this work, we introduce Stark time crystals based on a type of localization that is created in the absence of any spatial disorder. We argue that Stark time crystals constitute a phase of matter coming very close to the original idea and exhibit a symmetry breaking in space and time. Complementing a comprehensive discussion of the physics of the problem, we move on to elaborating on possible practical applications and argue that the physical demands of witnessing genuine signatures of many-body localization in large systems may be lessened in such physical systems., Comment: 8 pages, 13 figures, replaced with final version

We present fresh evidence for the presence of discrete quantum time crystals in two spatial dimensions. Discrete time crystals are intricate quantum systems that break discrete time translation symmetry in driven quantum many-body systems undergoing non-equilibrium dynamics. They are stabilized by many-body localization arising from disorder. We directly target the thermodynamic limit using instances of infinite tensor network states and implement disorder in a translationally invariant setting by introducing auxiliary systems at each site. We discuss how such disorder can be realized in programmable quantum simulators: This gives rise to the interesting situation in which a classical tensor network simulation can contribute to devising a blueprint of a quantum simulator featuring pre-thermal time crystalline dynamics, one that will yet ultimately have to be built in order to explore the stability of this phase of matter for long times., Comment: 13 pages, 9 figures, replaced by final version

We investigate the interacting, one-dimensional Rice-Mele model, a prototypical fermionic model of topological properties. To set the stage, we firstly compute the single-particle spectral function, the local density, and the boundary charge in the absence of interactions. The boundary charge is fully determined by bulk properties indicating a bulk-boundary correspondence. In a large parameter regime it agrees with the one obtained from an effective low-energy theory (arXiv:2004.00463). Secondly, we investigate the robustness of our results towards two-particle interactions. To resum the series of leading logarithms for small gaps, which dismantle plain perturbation theory in the interaction, we use an essentially analytical renormalization group approach. It is controlled for small interactions and can directly be applied to the microscopic lattice model. We benchmark the results against numerical density matrix renormalization group data. The main interaction effect in the bulk is a power-law renormalization of the gap with an interaction dependent exponent. The important characteristics of the boundary charge are unaltered and can be understood from the renormalized bulk properties, elevating the bulk-boundary correspondence to the interacting regime. This requires a consistent treatment not only of the low-energy gap renormalization but also of the high-energy band width one. In contrast to low-energy field theories our renormalization group approach also provides the latter. We show that the interaction spoils the relation between the bulk properties and the number of edge states, consistent with the observation that the Rice-Mele model with finite potential modulation does not reveal any zero-energy edge states.

We review the intriguing many-body physics resulting out of the interplay of a single, local impurity and the two-particle interaction in a one-dimensional Fermi system. Even if the underlying homogeneous correlated system is taken to be metallic, this interplay leads to an emergent quantum phase transition between metallic and insulating states. We show that the zero temperature critical point and the universal low-energy physics associated to it, is realized in two different models, the field theoretical local sine-Gordon model and spinless fermions on a lattice with nearest-neighbor hopping and two-particle interaction, as well as in an experimental setup consisting of a highly tunable quantum circuit. Despite the different high-energy physics of the three systems the universal low-energy scaling curves of the conductance as a function of temperature agree up to a very high precision without any free parameter. Overall this provides a convincing example of how emergent universality in complex systems originating from a common underlying quantum critical point establishes a bridge between different fields of physics. In our case between field theory, quantum many-body theory of correlated Fermi systems, and experimental circuit quantum electrodynamics., Comment: version as accepted for publication in EPJST

Experimental advances in the fabrication and characterization of few-layer materials stacked at a relative twist of small angle have recently shown the emergence of flat energy bands. As a consequence electron interactions become relevant, providing inroads into the physics of strongly correlated two-dimensional systems. Here, we demonstrate by combining large scale ab initio simulations with numerically exact strong correlation approaches that an effective one-dimensional system emerges upon stacking two twisted sheets of GeSe, in marked contrast to all Moir\'e systems studied so far. This not only allows to study the necessarily collective nature of excitations in one dimension, but can also serve as a promising platform to scrutinize the crossover from two to one dimension in a controlled setup by varying the twist angle, which provides an intriguing benchmark with respect to theory. We thus establish twisted bilayer GeSe as an intriguing inroad into the strongly correlated physics of low-dimensional systems., Comment: version as published

Many-body localized systems in which interactions and disorder come together defy the expectations of quantum statistical mechanics: In contrast to ergodic systems, they do not thermalize when undergoing nonequilibrium dynamics. What is less clear, however, is how topological features interplay with many-body localized phases as well as the nature of the transition between a topological and a trivial state within the latter. In this work, we numerically address these questions, using a combination of extensive tensor network calculations, specifically DMRG-X, as well as exact diagonalization, leading to a comprehensive characterization of Hamiltonian spectra and eigenstate entanglement properties., Comment: 7 pages, 8 figures, submitted to PRB

Dissipationless charge transport is one of the defining properties of superconductors (SC). The interplay between dimensionality and disorder in determining the onset of dissipation in SCs remains an open theoretical and experimental problem. In this work, we present measurements of the dissipation phase diagrams of SCs in the two dimensional (2D) limit, layer by layer, down to a monolayer in the presence of temperature (T), magnetic field (B), and current (I) in 2H-NbSe2. Our results show that the phase-diagram strongly depends on the SC thickness even in the 2D limit. At four layers we can define a finite region in the I-B phase diagram where dissipationless transport exists at T=0. At even smaller thicknesses, this region shrinks in area. In a monolayer, we find that the region of dissipationless transport shrinks towards a single point, defined by T=B=I=0. In applied field, we show that time-dependent-Ginzburg-Landau (TDGL) simulations that describe dissipation by vortex motion, qualitatively reproduce our experimental I-B phase diagram. Last, we show that by using non-local transport and TDGL calculations that we can engineer charge flow and create phase boundaries between dissipative and dissipationless transport regions in a single sample, demonstrating control over non-equilibrium states of matter., Comment: Manuscript, figures and supplemental information

Using projected entangled-pair states (PEPS) we analyze the localization properties of two-dimensional systems on a square lattice. We compare the dynamics found for three different disorder types: (i) quenched disorder, (ii) sum of two quasi-periodic potentials along both spatial dimensions and (iii) a single quasi-periodic potential rotated with respect to the underlying lattice by a given angle. We establish the rate of loss of information, a quantity measuring the error made while simulating the dynamics, as a good hallmark of localization physics by comparing to entanglement build-up as well as the inverse participation ratio in exactly solvable limits. We find that the disorder strength needed to localize the system increases both with the dimensionality of as well as the interaction strength in the system. The first two cases of potential (i) and (ii) behave similar, while case (iii) requires larger disorder strength to localize., Comment: 4 pages, 4 figures, submitted to PRB(R)

We consider the steady state behavior of observables in the Su-Schrieffer-Heeger model and in the one dimensional transverse field quantum Ising model after a sudden quantum quench of the parameter controlling the gap. In the thermodynamic limit, and for times $t\rightarrow\infty$, we find non-analyticities even in simple local observables as a function of the quench parameter, that is, a non-equilibrium quantum phase transition. We trace the appearance of this non-equilibrium quantum phase transition to an effective metal-insulator transition which occurs on the level of the generalized Gibbs Ensemble (describing the steady-state of the equilibrated system). Studying whether these transitions are robust, we find, in the paradigmatic case of the SSH model, that they persist for both quantum quench protocols of finite duration in time as well as thermal initial states, while they are washed out in the presence of fermion-fermion interactions and for finite system size., Comment: 12 pages, 6 figures

We introduce the Loschmidt amplitude as a powerful tool to perform spectroscopy of generic many-body wave functions and use it to interrogate the wave function obtained after ramping the transverse field quantum Ising model through its quantum critical point. Previous results are confirmed and a more complete understanding of the population of defects and of the effects of magnon-magnon interaction or finite-size corrections is obtained. The influence of quantum coherence is clarified., Comment: 4.5 pages + SM, 5 + 3 figures, accepted version

We study the dynamics arising from a double quantum quench where the parameters of a given Hamiltonian are abruptly changed from being in an equilibrium phase A to a different phase B and back (A$\to$B$\to$A). As prototype models, we consider the (integrable) transverse field Ising as well as the (non-integrable) ANNNI model. The return amplitude features non-analyticities after the first quench through the equilibrium quantum critical point (A$\to$B), which is routinely taken as a signature of passing through a so-called dynamical quantum phase transition. We demonstrate that non-analyticities after the second quench (B$\to$A) can be avoided and reestablished in a recurring manner upon increasing the time $T$ spent in phase B. The system retains an infinite memory of its past state, and one has the intriguing opportunity to control at will whether or not dynamical quantum phase transitions appear after the second quench., Comment: 6 pages, 5 figures

We consider a one-dimensional interacting spinless fermion model, which displays the well-known Luttinger liquid (LL) to charge density wave (CDW) transition as a function of the ratio between the strength of the interaction, $U$, and the hopping, $J$. We subject this system to a spatially uniform drive which is ramped up over a finite time interval and becomes time-periodic in the long time limit. We show that by using a density matrix renormalization group (DMRG) approach formulated for infinite system sizes, we can access the large-time limit even when the drive induces finite heating. When both the initial and long-time states are in the gapless (LL) phase, the final state has power law correlations for all ramp speeds. However, when the initial and final state are gapped (CDW phase), we find a pseudothermal state with an effective temperature that depends on the ramp rate, both for the Magnus regime in which the drive frequency is very large compared to other scales in the system and in the opposite limit where the drive frequency is less than the gap. Remarkably, quantum defects (instantons) appear when the drive tunes the system through the quantum critical point, in a realization of the Kibble-Zurek mechanism., Comment: 4.5 + 5 pages, 5 figures, version as published

We have used a combination of ultrafast coherent phonon spectroscopy, ultrafast thermometry, and time-dependent Landau theory to study the inversion symmetry breaking phase transition at $T_c = 200$ K in the strongly spin-orbit coupled correlated metal Cd$_2$Re$_2$O$_7$. We establish that the structural distortion at $T_c$ is a secondary effect through the absence of any softening of its associated phonon mode, which supports a purely electronically driven mechanism. However, the phonon lifetime exhibits an anomalously strong temperature dependence that decreases linearly to zero near $T_c$. We show that this behavior naturally explains the spurious appearance of phonon softening in previous Raman spectroscopy experiments and should be a prevalent feature of correlated electron systems with linearly coupled order parameters., Comment: 5 pages main text, 5 figures, 7 pages supplementary information

We study correlated quantum wires subject to harmonic modulation of the onsite-potential concentrating on the limit of large times, where the response of the system has synchronized with the drive. We identify the ratio $\Delta\epsilon/\Omega$ of the driving amplitude $\Delta\epsilon$ and the frequency of driving $\Omega$ as the scale determining the crossover from a modified Luttinger liquid picture to a system that behaves effectively like a higher dimensional one. We exemplify this crossover by studying the frequency dependency of the boundary density of state $\rho_{\rm B}(\omega)$ as well as the temperature dependency of the linear conductance $G(T)$ through the wire, if the latter is contacted to leads. Both observables are known to exhibit Luttinger liquid physics without driving given by characteristic power-law suppression as $\omega\to\epsilon_{\rm F}$ (with $\epsilon_{\rm F}$ the Fermi energy) or $T\to 0$, respectively. With driving we find that this suppression is modified from a single power-law to a superposition of an infinite number of power laws. At small $\Delta\epsilon/\Omega\ll 1$ only a few terms of this infinite sum are relevant as the prefactors of higher terms are suppressed exponentially. Thus a picture similar to the equilibrium Luttinger liquid one emerges. Increasing $\Delta\epsilon/\Omega$ an increasing number of power laws contribute to the sum and approaching $\Delta\epsilon/\Omega\gg 1$ the system behaves effectively two dimensional, for which the suppression is wiped out completely., Comment: 13 pages, 8 figures

The properties of prototypical examples of one-dimensional fermionic systems undergoing a sudden quantum quench from a gapless state to a (partially) gapped state are analyzed. By means of a Generalized Gibbs Ensemble analysis or by numerical solutions in the interacting cases, we observe an anomalous, non-monotonic response of steady state correlation functions as a function of the strength of the mechanism opening the gap. In order to interpret this result, we calculate the full dynamical evolution of these correlation functions, which shows a freezing of the propagation of the quench information (light cone) for large quenches. We argue that this freezing is responsible for the non-monotonous behaviour of observables. In continuum non-interacting models, this freezing can be traced back to a Klein-Gordon equation in the presence of a source term. We conclude by arguing in favour of the robustness of the phenomenon in the cases of non-sudden quenches and higher dimensionality., Comment: 5+11 pages, 6 figures

We use a numerical solution of the deterministic TDGL equations to determine the response induced by a probe field in a material quenched into a superconducting state. We characterize differences in response according to whether the probe is applied before, during, or after the phase stiffness has built up to its final steady state value., Comment: 10 pages 12 figures, version as published

We propose a route towards engineering non-thermal states of matter, which show largely unexplored physics. The main idea relies on the adiabatic passage of a thermal ensemble under slow variations of the system Hamiltonian. If the temperature of the initial thermal ensemble is either zero or infinite the ensemble after the passage is a simple thermal one with the same vanishing or infinite temperature. However, for any finite non-zero temperature intriguing non-thermal ensembles can be achieved. We exemplify this in: (a) a single oscillator (b) a dimerized interacting one dimensional chain of spinless fermions, (c) a BCS-type superconductor and (d) the topological Kitaev chain. We solve these models with a combination of methods; either exactly, numerically using the density matrix renormalization group (DMRG) or within an approximate functional renormalization group (FRG) scheme. The designed states show strongly non-thermal behavior in each of the considered models. For example, for the chain of spinless fermions we exemplify how long ranged non-thermal power-law correlations can be stabilized and for the Kitaev chain we elucidate how the non-thermal ensemble can largely alter the transition temperature separating topological and trivial phases., Comment: 11 pages, 9 figures, version as published

We study the expectation values of observables and correlation functions at long times after a global quantum quench. Our focus is on metallic (`gapless') fermionic many-body models and small quenches. The system is prepared in an eigenstate of an initial Hamiltonian, and the time evolution is performed with a final Hamiltonian which differs from the initial one in the value of one global parameter. We first derive general relations between time-averaged expectation values of observables as well as correlation functions and those obtained in an eigenstate of the final Hamiltonian. Our results are valid to linear and quadratic order in the quench parameter g and generalize prior insights in several essential ways. This allows us to develop a phenomenology for the thermalization of local quantities up to a given order in g. Our phenomenology is put to a test in several case studies of one-dimensional models representative of four distinct classes of Hamiltonians: quadratic ones, effectively quadratic ones, those characterized by an extensive set of (quasi-) local integrals of motion, and those for which no such set is known (and believed to be nonexistent). We show that for each of these models, all observables and correlation functions thermalize to linear order in g. The more local a given quantity, the longer the linear behavior prevails when increasing g. Typical local correlation functions and observables for which the term O(g) vanishes thermalize even to order g^2. Our results show that lowest order thermalization of local observables is an ubiquitous phenomenon even in models with extensive sets of integrals of motion., Comment: version as published

We investigate the thermodynamics and transient dynamics of the (unbiased) Ohmic two-state system by exploiting the equivalence of this model to the interacting resonant level model. For the thermodynamics, we show, by using the numerical renormalization group (NRG) method, how the universal specific heat and susceptibility curves evolve with increasing dissipation strength, $\alpha$, from those of an isolated two-level system at vanishingly small dissipation strength, with the characteristic activated-like behavior in this limit, to those of the isotropic Kondo model in the limit $\alpha\to 1^{-}$. For the transient dynamics of the two-level system, $P(t)=\langle \sigma_{z}(t)\rangle$, with initial-state preparation $P(t\leq0)=+1$, we apply the time-dependent extension of the NRG (TDNRG) to the interacting resonant level model, and compare the results obtained with those from the noninteracting-blip approximation (NIBA), the functional renormalization group (FRG), and the time-dependent density matrix renormalization group (TD-DMRG). We demonstrate excellent agreement on short to intermediate time scales between TDNRG and TD-DMRG for $0\lesssim\alpha \lesssim 0.9$ for $P(t)$, and between TDNRG and FRG in the vicinity of $\alpha=1/2$. Furthermore, we quantify the error in the NIBA for a range of $\alpha$, finding significant errors in the latter even for $0.1\leq \alpha\leq 0.4$. We also briefly discuss why the long-time errors in the present formulation of the TDNRG prevents an investigation of the crossover between coherent and incoherent dynamics. Our results for $P(t)$ at short to intermediate times could act as useful benchmarks for the development of new techniques to simulate the transient dynamics of spin-boson problems., Comment: 27 pages, 22 figures; published version

We study the properties of excited states in one-dimensional many-body localized (MBL) systems using a matrix product state algorithm. First, the method is tested for a large disordered non-interacting system, where for comparison we compute a quasi-exact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of L~100 sites and large bond dimensions chi~1700, which allows us to quantitatively analyze the scaling behavior of the entanglement S in the system. The MBL phase is characterized by a logarithmic growth (L)~log(L) over a large scale separating the regimes where volume and area laws hold. We check the validity of the eigenstate thermalization hypothesis. Our results are consistent with the existence of a mobility edge.

We study the zero-frequency current noise of the interacting resonant level model for arbitrary bias voltages using a functional renormalization group approach. For this we extend the existing nonequilibrium scheme by deriving and solving flow equations for the current-vertex functions. On-resonance artificial divergences of the latter found in lowest-order perturbation theory in the two-particle interaction are consistently removed. Away from resonance they are shifted to higher orders. This allows us to gain a comprehensive picture of the current noise in the scaling limit. At high bias voltages, the current noise exhibits a universal power-law decay, whose exponent is to leading order in the interaction identical to that of the current. The effective charge on resonance is analyzed in detail employing properties of the vertex correction. We find that it is only modified to second or higher order in the two-particle interaction., Comment: 15 pages, 16 figures

We report on strong renormalization encountered in periodically driven interacting quantum dots in the non-adiabatic regime. Correlations between lead and dot electrons enhance or suppress the amplitude of driving depending on the sign of the interaction. Employing a newly developed flexible renormalization group based approach for periodic driving to an interacting resonant level we show analytically that the magnitude of this effect follows a power law. Our setup can act as a non-Markovian, single-parameter quantum pump., Comment: 6 pages, 3 figures

We study the thermal conductivity of the one-dimensional Fermi-Hubbard model at finite temperature using a density matrix renormalization group approach. The integrability of this model gives rise to ballistic thermal transport. We calculate the temperature dependence of the thermal Drude weight at half filling for various interactions and moreover, we compute its filling dependence at infinite temperature. The finite-frequency contributions originating from the fact that the energy current is not a conserved quantity are investigated as well. We report evidence that breaking the integrability through a nearest-neighbor interaction leads to vanishing Drude weights and diffusive energy transport. Moreover, we demonstrate that energy spreads ballistically in local quenches with initially inhomogeneous energy density profiles in the integrable case. We discuss the relevance of our results for thermalization in ultra-cold quantum gas experiments and for transport measurements with quasi-one dimensional materials.

We study the spin and thermal conductivity of spin-1/2 ladders at finite temperature. This is relevant for experiments with quantum magnets. Using a state-of-the-art density matrix renormalization group algorithm, we compute the current autocorrelation functions on the real-time axis and then carry out a Fourier integral to extract the frequency dependence of the corresponding conductivities. The finite-time error is analyzed carefully. We first investigate the limiting case of spin-1/2 XXZ chains, for which our analysis suggests non-zero dc-conductivities in all interacting cases irrespective of the presence or absence of spin Drude weights. For ladders, we observe that all models studied are normal conductors with no ballistic contribution. Nonetheless, only the high-temperature spin conductivity of XX ladders has a simple diffusive, Drude-like form, while Heisenberg ladders exhibit a more complicated low-frequency behavior. We compute the dc spin conductivity down to temperatures of the order of T~0.5J, where J is the exchange coupling along the legs of the ladder. We further extract mean-free paths and discuss our results in relation to thermal conductivity measurements on quantum magnets.

We show that to understand the orthogonality catastrophe in the half-filled lattice model of spinless fermions with repulsive nearest neighbor interaction and a local impurity in its Luttinger liquid phase one has to take into account (i) the impurity scaling, (ii) unusual finite size $L$ corrections of the form $\ln(L)/L$, as well as (iii) the renormalization group flow of the umklapp scattering. The latter defines a length scale $L_u$ which becomes exceedingly large the closer the system is to its transition into the charge-density wave phase. Beyond this transition umklapp scattering is relevant in the renormalization group sense. Field theory can only be employed for length scales larger than $L_u$. For small to intermediate two-particle interactions, for which the regime $L > L_u$ can be accessed, and taking into account the finite size corrections resulting from (i) and (ii) we provide strong evidence that the impurity backscattering contribution to the orthogonality exponent is asymptotically given by $1/16$. While further increasing the two-particle interaction leads to a faster renormalization group flow of the impurity towards the cut chain fixed point, the increased bare amplitude of the umklapp scattering renders it virtually impossible to confirm the expected asymptotic value of $1/16$ given the accessible system sizes. We employ the density matrix renormalization group., Comment: 12 pages, 9 figures

We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is two-fold: First, we compute both the charge and the spin current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity sigma(omega) in presence of small integrability-breaking next-nearest neighbor interactions (the extended Hubbard model). The DC conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a testbed for our numerics, we compute sigma(omega) for the integrable XXZ spin chain in its gapped phase.

The equilibrium physics of quantum impurities frequently involves a universal crossover from weak to strong reservoir-impurity coupling, characterized by single-parameter scaling and an energy scale $T_K$ (Kondo temperature) that breaks scale invariance. For the non-interacting resonant level model, the non-equilibrium time evolution of the Loschmidt echo after a local quantum quench was recently computed explicitely [R. Vasseur, K. Trinh, S. Haas, and H. Saleur, Phys. Rev. Lett. 110, 240601 (2013)]. It shows single-parameter scaling with variable $T_K t$. Here, we scrutinize whether similar universal dynamics can be observed in various interacting quantum impurity systems. Using density matrix and functional renormalization group approaches, we analyze the time evolution resulting from abruptly coupling two non-interacting Fermi or interacting Luttinger liquid leads via a quantum dot or a direct link. We also consider the case of a single Luttinger liquid lead suddenly coupled to a quantum dot. We investigate whether the field theory predictions for the universal scaling as well as for the large time behavior successfully describe the time evolution of the Loschmidt echo and the entanglement entropy of microscopic models., Comment: 14 pages, 10 figures

We show that the single-particle spectral properties of gapless one-dimensional Fermi systems in the Luttinger liquid state reached at intermediate times after an abrupt quench of the two-particle interaction are highly indicative of the unusual nonequilibrium nature of this state. The line shapes of the momentum integrated and resolved spectral functions strongly differ from their ground state as well as finite temperature equilibrium counterparts. Using an energy resolution improved version of radio-frequency spectroscopy of quasi one-dimensional cold Fermi gases it should be possible to experimentally identify this nonequilibrium state by its pronounced spectral signatures., Comment: 5 pages, 3 figures

We discuss a few simple modifications to time-dependent density matrix renormalization group (DMRG) algorithms which allow to access larger time scales. We specifically aim at beginners and present practical aspects of how to implement these modifications almost effortlessly within any standard matrix product state (MPS) based formulation of the method. Most importantly, we show how to 'combine' the Schroedinger and Heisenberg time evolutions of arbitrary pure states |psi> and operators A in the evaluation of _psi(t)=. This includes quantum quenches. The generalization (non-)thermal mixed state dynamics _rho(t)=Tr[rho A(t)] induced by an initial density matrix rho is straightforward. In the context of equilibrium (ground state or finite temperature T>0) correlation functions, one can extend the simulation time by a factor of two by 'exploiting time translation invariance', which is efficiently implementable within MPS DMRG. We present a simple analytic argument for why a recently-introduced disentangler succeeds in reducing the effort of time-dependent simulations at T>0. Finally, we advocate the python programming language as an elegant option for beginners to set up a DMRG code.

We reveal several distinct regimes of the relaxation dynamics of a small quantum system coupled to an environment within the plane of the dissipation strength and the reservoir temperature. This is achieved by discriminating between coherent dynamics with damped oscillatory behavior on all time scales, partially coherent behavior being nonmonotonic at intermediate times but monotonic at large ones, and purely monotonic incoherent decay. Surprisingly, elevated temperature can render the system `more coherent' by inducing a transition from the partially coherent to the coherent regime. This provides a refined view on the relaxation dynamics of open quantum systems., Comment: 5 pages, 3 figures

We present a comprehensive study of the spectral and transport properties in the Anderson--Holstein model both in and out of equilibrium using the functional renormalization group (FRG). We show how the previously established machinery of Matsubara and Keldysh FRG can be extended to include the local phonon mode. Based on the analysis of spectral properties in equilibrium we identify different regimes depending on the strength of the electron--phonon interaction and the frequency of the phonon mode. We supplement these considerations with analytical results from the Kondo model. We also calculate the non-linear differential conductance through the Anderson--Holstein quantum dot and find clear signatures of the presence of the phonon mode., Comment: 19 pages, 8 figures

We study the dynamics resulting out of an abrupt change of the two-particle interaction in two models of closed one-dimensional Fermi systems: (a) the field theoretical Tomonaga-Luttinger model and (b) a microscopic lattice model. Using a nonperturbative approach which is controlled for small two-particle interactions we are able to reach large times allowing us to access the properties of the steady state of the lattice model. Comparing those to the exact solution of the full dynamics in the Tomonaga-Luttinger model we provide evidence for universal Luttinger liquid behavior., Comment: 9 pages, 5 figures, 1 table as accepted for publication by PRB

We study linear response and nonequilibrium steady-state thermoelectric transport through a single-level quantum dot tunnel coupled to two reservoirs held at different temperatures as well as chemical potentials. A fermion occupying the dot interacts with those in the reservoirs by a short-ranged two-particle interaction. For parameters for which particles flow against a bias voltage from the hot to the cold reservoir this setup acts as an energy-conversion device with which electrical energy is gained out of waste heat. We investigate how correlations affect its efficiency and output power. In linear response the changes in the thermoelectric properties can be traced back to the interaction induced renormalization of the resonance line shape. In particular, small to intermediate repulsive interactions reduce the maximum efficiency. In nonequilibrium the situation is more complex and we identify a parameter regime in which for a fixed lower bound of the output power the efficiency increases., Comment: 6 pages, 6 figures

The nonequilibrium dynamics of a small quantum system coupled to a dissipative environment is studied. We show that (1) the oscillatory dynamics close to a coherent-to-incoherent transition is surprisingly different from the one of the classical damped harmonic oscillator and that (2) non-markovian memory plays a prominent role in the time evolution after a quantum quench., Comment: 5 pages, 3 figures

We study the steady-state properties as well as the relaxation dynamics of the nonequilibrium interacting resonant level model at finite temperatures. It constitutes the prototype model of a correlated charge fluctuating quantum dot. The two reservoirs are held at different chemical potentials---the difference being the bias voltage---and different temperatures; we discuss the transport through as well as the occupancy of the single level dot. First, we show analytically that in the steady state the reservoir temperatures in competition with the other energy scales act as infrared cutoffs. This is rather intuitive but, depending on the parameter regime under consideration, leads to a surprisingly rich variety of power laws in the current as a function of the temperatures and the bias voltage with different interaction dependent exponents. Next we clarify how finite reservoir temperatures affect the dynamics. They allow to tune the interplay of the two frequencies characterizing the oscillatory part of the time evolution of the model at zero temperature. For the exponentially decaying part we disentangle the contributions of the level-lead hybridization and the temperatures to the decay rates. We identify a coherent-to-incoherent transition in the long time dynamics as the temperature is raised. It occurs at an interaction dependent critical temperature. Finally, taking different temperatures in the reservoirs we discuss the relaxation dynamics of a temperature gradient driven current., Comment: 12 pages, 6 figures, 1 table

We introduce a real time version of the functional renormalization group which allows to study correlation effects on nonequilibrium transport through quantum dots. Our method is equally capable to address (i) the relaxation out of a nonequilibrium initial state into a (potentially) steady state driven by a bias voltage and (ii) the dynamics governed by an explicitly time-dependent Hamiltonian. All time regimes from transient to asymptotic can be tackled; the only approximation is the consistent truncation of the flow equations at a given order. As an application we investigate the relaxation dynamics of the interacting resonant level model which describes a fermionic quantum dot dominated by charge fluctuations. Moreover, we study decoherence and relaxation phenomena within the ohmic spin-boson model by mapping the latter to the interacting resonant level model.

We address the question whether observables of an exactly solvable model of electrons coupled to (optical) phonons relax into large time stationary state values and investigate if the asymptotic expectation values can be computed using a stationary density matrix. Two initial nonequilibrium situations are considered. A sudden quench of the electron-phonon coupling, starting from the noninteracting canonical equilibrium at temperature T in the electron as well as in the phonon subsystems, leads to a rather simple dynamics. A richer time evolution emerges if the initial state is taken as the product of the phonon vacuum and the filled Fermi sea supplemented by a highly excited additional electron. Our model has a natural set of constants of motion, with as many elements as degrees of freedom. In accordance with earlier studies of such type of models we find that expectation values which become stationary can be described by the density matrix of a generalized Gibbs ensemble which differs from that of a canonical ensemble. For the model at hand it appears to be evident that the eigenmode occupancy operators should be used in the construction of the stationary density matrix., Comment: 15 pages, 11 figures, published version