Search

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. We propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to non-perturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasi-particle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group (DMRG) simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments., 7 pages, 3 figures plus supplemental materials

We compute ground-state properties of the isotropic, antiferromagnetic Heisenberg model on the sodalite cage geometry. This is a 60-spin spherical molecule with 24 vertex-sharing tetrahedra which can be regarded as a molecular analogue of a capped kagome lattice and which has been synthesized with high-spin rare-earth atoms. Here, we focus on the S=1/2S=1/2 case where quantum effects are strongest. We employ the SU(2)-symmetric density-matrix renormalization group (DMRG). We find a threefold degenerate ground state that breaks the spatial symmetry and that splits up the molecule into three large parts which are almost decoupled from each other. This stands in sharp contrast to the behaviour of most known spherical molecules. On a methodological level, the disconnection leads to ``glassy dynamics’’ within the DMRG that cannot be targeted via standard techniques. In the presence of finite magnetic fields, we find broad magnetization plateaus at 4/5, 3/5, and 1/5 of the saturation, which one can understand in terms of localized magnons, singlets, and doublets which are again nearly decoupled from each other. At the saturation field, the zero-point entropy is S=\ln(182)\approx 5.2S=ln(182)≈5.2 in units of the Boltzmann constant.

The ability to tune the electronic properties of molecular arrays is an important step in the development of molecule-scale electronic devices. However, control over internal device charge distributions by tuning interactions between molecules has proved challenging. Here, we show that gate-tunable charge patterning can occur in one-dimensional molecular arrays on graphene field-effect transistors. One-dimensional molecular arrays are fabricated using an edge-templated self-assembly process that allows organic molecules (F4TCNQ) to be precisely positioned on graphene devices. The charge configurations of the molecular arrays can be reversibly switched between different collective charge states by tuning the graphene Fermi level via a back-gate electrode. Charge pinning at the ends of the molecular arrays allows the charge state of the entire array to be controlled by adding or removing an edge molecule and changing the total number of molecules in an array between odd and even integers. Charge patterns altered in this way propagate down the array in a cascade effect, allowing the array to function as a charge-based molecular shift register. An extended multi-site Anderson impurity model is used to quantitatively explain this behaviour. One-dimensional molecular arrays on graphene field-effect transistors can be reversibly switched between different periodic charge states by tuning the graphene Fermi level via a back-gate electrode and by manipulating individual molecules, allowing them to function as a nanoscale shift register.

This upload contains all code, raw data, processed data and plots used in a publication released under the same name as the name of this upload. The majority of the uploaded data consist of matrix product states (MPS) which have been generated through the DMRG-X algorithm and stored as all further data and plots derive from them. The file system of choice to store the MPSs has been HDF5. By following the steps described in various README files, the reader should be able to reproduce the results shown in the provided plots and also of the linked publication. In addition to that, the reader should also be able to verify whether the processed data has been generated from the raw data, i.e., the MPSs and whether the MPSs have been generated from the provided code. In order to do this, numerical information such as selected parameters or used disorder realizations are provided in the respective HDF5 files, thus allowing perfect reproducibility of all data uploaded., K.D. and C.K. acknowledge support by the Deutsche Forschungsgemeinschaft (DFG) through the Emmy Noether program (Grant No. KA3360/2-1) as well as by "Niedersächsisches Vorab" through the "Quantum- and Nano-Metrology (QUANOMET)" initiative within Project No. P-1. D.M.K. was supported by the DFG via Germany's Excellence Strategy—Cluster of Excellence Matter and Light for Quantum Computing (ML4Q, Project No. EXC 2004/1, Grant No. 390534769). We acknowledge support from the Max Planck–New York City Center for Non-Equilibrium Quantum Phenomena.

The last decade has witnessed an impressive progress in the theoretical understanding of transport properties of clean, one-dimensional quantum lattice systems. Many physically relevant models in one dimension are Bethe-ansatz integrable, including the anisotropic spin-1/2 Heisenberg (also called spin-1/2 XXZ chain) and the Fermi-Hubbard model. Nevertheless, practical computations of, for instance, correlation functions and transport coefficients pose hard problems from both the conceptual and technical point of view. Only due to recent progress in the theory of integrable systems on the one hand and due to the development of numerical methods on the other hand has it become possible to compute their finite temperature and nonequilibrium transport properties quantitatively. Most importantly, due to the discovery of a novel class of quasilocal conserved quantities, there is now a qualitative understanding of the origin of ballistic finite-temperature transport, and even diffusive or super-diffusive subleading corrections, in integrable lattice models. We shall review the current understanding of transport in one-dimensional lattice models, in particular, in the paradigmatic example of the spin-1/2 XXZ and Fermi-Hubbard models, and we elaborate on state-of-the-art theoretical methods, including both analytical and computational approaches. Among other novel techniques, we discuss matrix-product-states based simulation methods, dynamical typicality, and, in particular, generalized hydrodynamics. We will discuss the close and fruitful connection between theoretical models and recent experiments, with examples from both the realm of quantum magnets and ultracold quantum gases in optical lattices., Review article (78 pages, 34 figures). V3: Revised version, additional references

The functional renormalization group (FRG) provides a flexible tool to study correlations in low-dimensional electronic systems. In this paper, we present a novel FRG approach to the steady-state of quantum wires out of thermal equilibrium. Our method is correct up to second order in the two-particle interaction and accounts for inelastic scattering. We combine semi-analytic solutions of the flow equations with MPI parallelization techniques, which allows us to treat systems of up to 60 lattice sites. The equilibrium limit is well-understood and serves as a benchmark. We compute effective distribution functions, the local density of states, and the steady-state current and demonstrate that all of these quantities depend strongly on the choice of the cutoff employed within the FRG. Non-equilibrium is plagued by the lack of physical arguments in favor of a certain cutoff as well as by the appearance of secular higher-order terms which are only partly included in our approach. This demonstrates the inadequacy of a straightforward second-order FRG scheme to study interacting quantum wires out of equilibrium in the absence of a natural cutoff choice.

We employ a novel, unbiased renormalization-group approach to investigate non-equilibrium phase transitions in infinite lattice models. This allows us to address the delicate interplay of fluctuations and ordering tendencies in low dimensions out of equilibrium. We study a prototypical model for the metal to insulator transition of spinless interacting fermions coupled to electronic baths and driven out of equilibrium by a longitudinal static electric field. The closed system features a Berezinskii-Kosterlitz-Thouless transition between a metallic and a charge-ordered phase in the equilibrium limit. We compute the non-equilibrium phase diagram and illustrate a highly non-monotonic dependence of the phase boundary on the strength of the electric field: For small fields, the induced currents destroy the charge order, while at higher electric fields it reemerges due to many-body Wannier-Stark localization physics. Finally, we show that the current in such an interacting non-equilibrium system can counter-intuitively flow opposite to the direction of the electric field. This non-equilibrium steady-state is reminiscent of an equilibrium distribution function with an effective negative temperature.

New journal of physics 22(8), 083039 (2020). doi:10.1088/1367-2630/ab990d, Published by IOP, [London]

This upload contains all code, raw data, processed data and plots used in a publication released under the same name as the name of this upload. By following the steps described in various README files, the reader should be able to reproduce the results shown in the provided plots and also of the linked publication. In addition to that, the reader should also be able to verify whether the processed data has been generated from the raw data and whether the raw data has been generated from the provided code. In order to do this, numerical information such as selected parameters or used disorder realizations are provided in the respective HDF5 files, thus allowing perfect reproducibility of all data uploaded., This work has been supported by the DFG (Grants No. CRC 183, No. FOR 2724, No. EI 519/14-1, and No. EI 519/15-1) and through the Emmy Noether Program (Program No. KA 3360/2-1), the ERC (TAQ), and the Templeton Foundation. This work has also received funding from the European Union's Horizon 2020 Research and Innovation Program under Grant Agreement No. 817482 (PASQuanS).

We report a systematic study of finite-temperature spin transport in quantum and classical one-dimensional magnets with isotropic spin interactions, including both integrable and non-integrable models. Employing a phenomenological framework based on a generalized Burgers' equation in a time-dependent stochastic environment, we identify four different universality classes of spin fluctuations. These comprise, aside from normal spin diffusion, three types of superdiffusive transport: the KPZ universality class and two distinct types of anomalous diffusion with multiplicative logarithmic corrections. Our predictions are supported by extensive numerical simulations on various examples of quantum and classical chains. Contrary to common belief, we demonstrate that even non-integrable spin chains can display a diverging spin diffusion constant at finite temperatures., Comment: Typos corrected, references added and new supplementary material A on the classical simulations

We show how second-order Floquet engineering can be employed to realize systems in which many-body localization coexists with topological properties in a driven system. This allows one to implement and dynamically control a topologically-protected qubit even at high energies. Floquet engineering - the idea that a periodically driven non-equilibrium system can effectively emulate the physics of a different Hamiltonian - is used to simulate an ffective three-body interaction among spins in one dimension, using time-dependent two-body interactions only. In the effective system emulated topology and disorder coexist which provides an intriguing inroad into the interplay of many-body localization, defying our standard understanding of thermodynamics, and topological phases of matter, which are of fundamental and technological importance. We demonstrate explicitly how combining Floquet engineering, topology and many-body localization allows one to harvest the advantages (time-dependent control, topological protection and reduction of heating, respectively) of each of these sub-fields while protecting from their disadvantages (heating, static control parameters and strong disorder)., Comment: 6 pages, 5 figures, cosmetic changes, replaced by published version

Metals in one spatial dimension are described at the lowest energy scales by the Luttinger liquid theory. It is well understood that this free theory, and even interacting integrable models, can support ballistic transport of conserved quantities including energy. In contrast, realistic one-dimensional metals, even without disorder, contain integrability-breaking interactions that are expected to lead to thermalization and conventional diffusive linear response. We argue that the expansion of energy when such a non-integrable Luttinger liquid is locally heated above its ground state shows superdiffusive behavior (i.e., spreading of energy that is intermediate between diffusion and ballistic propagation), by combining an analytical anomalous diffusion model with numerical matrix product state calculations on a specific perturbed spinless fermion chain. Different metals will have different scaling exponents and shapes in their energy spreading, but the superdiffusive behavior is stable and should be visible in time-resolved experiments., v3: 6 pages, 3 figures, text revised and improved

The problem of characterizing low-temperature spin dynamics in antiferromagnetic spin chains has so far remained elusive. We reinvestigate it by focusing on isotropic antiferromagnetic chains whose low-energy effective field theory is governed by the quantum non-linear sigma model. We outline an exact non-perturbative theoretical approach to analyse the low-temperature behaviour in the vicinity of non-magnetized states, and obtain explicit expressions for the spin diffusion constant and the NMR relaxation rate, which we compare with previous theoretical results in the literature. Surprisingly, in SU(2)-invariant spin chains in the vicinity of half-filling we find a crossover from the semi-classical regime to a strongly interacting quantum regime characterized by zero spin Drude weight and diverging spin conductivity, indicating super-diffusive spin dynamics. The dynamical exponent of spin fluctuations is argued to belong to the Kardar-Parisi-Zhang universality class. Furthermore, by employing numerical tDMRG simulations, we find robust evidence that the anomalous spin transport persists also at high temperatures, irrespectively of the spectral gap and integrability of the model., 5 pages + SM, 1 figure in the main text + 3 in the SM. v2: we corrected the typos and updated the figures

We compute the phase diagram of twisted bilayer graphene near the magic angle where the occurrence of flat bands enhances the effects of electron-electron interactions and thus unleashes strongly correlated phenomena. Most importantly, we find a crossover between d+id superconductivity and antiferromagnetic insulating behavior near half filling of the lowest electron band when the temperature is increased. This is consistent with recent experiments. Our results are obtained using unbiased many-body renormalization group techniques combined with a mean-field analysis of the effective couplings. We provide a qualitative understanding by considering the competition between Fermi-surface nesting and van Hove singularities.

We show that any interacting integrable model possesses a class of initial states for which the leading corrections to ballistic transport are subdiffusive rather than diffusive. These initial states are natural to realize experimentally and include the domain wall initial condition that has been the object of much recent scrutiny. Upon performing numerical matrix product state simulations in the spin-$1/2$ XXZ chain, we find that such states can exhibit subdiffusive $t^{1/3}$ scaling of fronts of spin, energy and entanglement entropy across the entire range of anisotropies. This demonstrates that Tracy-Widom scaling is not incompatible with model interactions, as was previously believed., v3: 5+3 pages, 6 figures, text revised

We present a functional renormalization-group approach to interacting topological Green's function invariants with a focus on the nature of transitions. The method is applied to chiral symmetric fermion chains in the Mott limit that can be driven into a Haldane phase. We explicitly show that the transition to this phase is accompanied by a zero of the fermion Green's function. Our results for the phase boundary are quantitatively benchmarked against density matrix renormalization-group data.

The functional renormalization group (FRG) has been used widely to investigate phase diagrams, in particular the one of the two-dimensional Hubbard model. So far, the study of one-dimensional models has not attracted as much attention. We use the FRG to investigate the phases of a one-dimensional spinless tight-binding chain with nearest and next-nearest neighbor interactions at half filling. The phase diagram of this model has already been established with other methods, and phase transitions from a metallic phase to ordered phases take place at intermediate to strong interactions. The model is thus well suited to analyze the potential and the limitations of the FRG in this regime of interactions. We employ flow equations that are exact up to second order in the interaction, which implies that we take into account the frequency dependence of the two-particle vertex as well as the feedback of the dynamic self-energy. For intermediate nearest neighbor interactions, our scheme captures the phase transition from a metallic phase to a charge density wave with alternating occupation. The critical interaction, at which this transition occurs, is underestimated due to our approximations. Similarly, for intermediate next-nearest neighbor interactions, we observe a transition to a charge density wave with occupation pattern ..00110011... We show that taking into account a feedback of the two-particle vertex in the flow equation is essential for the detection of those phases., 15 pages, 13 figures

We study the finite-energy density phase diagram of spinless fermions with attractive interactions in one dimension in the presence of uncorrelated diagonal disorder. Unlike the case of repulsive interactions, a delocalized Luttinger-liquid phase persists at weak disorder in the ground state, which is a well-known result. We revisit the ground-state phase diagram and show that the recently introduced occupation-spectrum discontinuity computed from the eigenspectrum of the one-particle density matrix is noticeably smaller in the Luttinger liquid compared to the localized regions. Moreover, we use the functional renormalization group scheme to study the finite-size dependence of the conductance, which also resolves the existence of the Luttinger liquid and is computationally cheap. Our main results concern the finite-energy density case. Using exact diagonalization and by computing various established measures of the many-body localization-delocalization transition, we argue that the zero-temperature Luttinger liquid smoothly evolves into a finite-energy density ergodic phase without any intermediate phase transition.

Quantum integrable systems, such as the interacting Bose gas in one dimension and the XXZ quantum spin chain, have an extensive number of local conserved quantities that endow them with exotic thermalization and transport properties. We discuss recently introduced hydrodynamic approaches for such integrable systems from the viewpoint of kinetic theory and extend the previous works by proposing a numerical scheme to solve the hydrodynamic equations for finite times and arbitrary locally equilibrated initial conditions. We then discuss how such methods can be applied to describe non-equilibrium steady states involving ballistic heat and spin currents. In particular, we show that the spin Drude weight in the XXZ chain, previously accessible only by rigorous techniques of limited scope or controversial thermodynamic Bethe ansatz arguments, may be evaluated from hydrodynamics in very good agreement with density-matrix renormalization group calculations., Comment: 11+4 pages, published version

We introduce a mechanism for light-induced Floquet engineering of the Fermi surface to dynamically tip the balance between competing instabilities in correlated condensed matter systems in the vicinity of a van-Hove singularity. We first calculate how the Fermi surface is deformed by an off-resonant, high-frequency light field and then determine the impact of this deformation on the ordering tendencies using an unbiased functional renormalization group approach. As a testbed, we investigate Floquet engineering in cuprates driven by light. We find that the $d$-wave superconducting ordering tendency in this system can be strongly enhanced over the Mott insulating one. This gives rise to extended regions of induced $d$-wave superconductivity in the effective phase diagram in the presence of a light field., Comment: 4 pages, 3 figures plus Supplementary Information

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 establish a general connection between ballistic and diffusive transport in systems where the ballistic contribution in canonical ensemble vanishes. A lower bound on the Green-Kubo diffusion constant is derived in terms of the curvature of the ideal transport coefficient, the Drude weight, with respect to the filling parameter. As an application, we explicitly determine the lower bound on the high temperature diffusion constant in the anisotropic spin 1/2 Heisenberg chain for anisotropy parameters $\Delta \geq 1$, thus settling the question whether the transport is sub-diffusive or not. Addi- tionally, the lower bound is shown to saturate the diffusion constant for a certain classical integrable model., Comment: 5+9 pages, 2 figures, v2 features major changes in presentation and additional example

The conventional theory of hydrodynamics describes the evolution in time of chaotic many-particle systems from local to global equilibrium. In a quantum integrable system, local equilibrium is characterized by a local generalized Gibbs ensemble or equivalently a local distribution of pseudo-momenta. We study time evolution from local equilibria in such models by solving a certain kinetic equation, the "Bethe-Boltzmann" equation satisfied by the local pseudo-momentum density. Explicit comparison with density matrix renormalization group time evolution of a thermal expansion in the XXZ model shows that hydrodynamical predictions from smooth initial conditions can be remarkably accurate, even for small system sizes. Solutions are also obtained in the Lieb-Liniger model for free expansion into vacuum and collisions between clouds of particles, which model experiments on ultracold one-dimensional Bose gases., 6+5 pages, published version

Disorder effects are especially pronounced around nodal points in linearly dispersing band structures as present in graphene or Weyl semimetals. Despite the enormous experimental and numerical progress, even a simple quantity like the average density of states cannot be assessed quantitatively by analytical means. We demonstrate how this important problem can be solved employing the functional renormalization group method, and, for the two-dimensional case, we demonstrate excellent agreement with reference data from numerical simulations based on tight-binding models. In three dimensions our analytic results also improve drastically on existing approaches.

We introduce approximate, functional renormalization group based schemes to obtain correlation functions in pure excited eigenstates of large fermionic many-body systems at arbitrary energies. The algorithms are thouroughly benchmarked and their strengths and shortcomings are documented using a one-dimensional interacting tight-binding chain as a prototypical testbed. We study two `toy applications' from the world of Luttinger liquid physics: the survival of power laws in lowly-excited states as well as the spectral function of high-energy `block' excitations which feature several single-particle Fermi edges., Comment: 17 pages, 8 figures; Corrected typos, improving some explanations; results unaffected

We devise a functional renormalization group treatment for a chain of interacting spinless fermions which is correct up to second order in the interaction strength. We treat both inhomogeneous systems in real-space as well as the translational invariant case in a k-space formalism. The strengths and shortcomings of the different schemes as well as technical details of their implementation are discussed. We use the method to study two proof-of-principle problems in the realm of Luttinger liquid physics, namely reflection at interfaces and power laws in the occupation number as a function of crystal momentum., Comment: 9 pages

The ability to understand and control the electronic properties of individual molecules in a device environment is crucial for developing future technologies at the nanometre scale and below. Achieving this, however, requires the creation of three-terminal devices that allow single molecules to be both gated and imaged at the atomic scale. We have accomplished this by integrating a graphene field effect transistor with a scanning tunnelling microscope, thus allowing gate-controlled charging and spectroscopic interrogation of individual tetrafluoro-tetracyanoquinodimethane molecules. We observe a non-rigid shift in the molecule's lowest unoccupied molecular orbital energy (relative to the Dirac point) as a function of gate voltage due to graphene polarization effects. Our results show that electron–electron interactions play an important role in how molecular energy levels align to the graphene Dirac point, and may significantly influence charge transport through individual molecules incorporated in graphene-based nanodevices., The development of single-molecule electronics calls for precise tuning of the electronic properties of individual molecules that go beyond two-terminal control. Here, Wickenburg et al. show gate-tunable switch of charge states of an isolated molecule using a graphene-based field-effect transistor.

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., version as published

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 noninteracting system, where for comparison we compute a quasiexact reference solution via a Monte Carlo sampling of the single-particle levels. Thereafter, we present extensive data obtained for large interacting systems of $L\ensuremath{\sim}100$ sites and large bond dimensions $\ensuremath{\chi}\ensuremath{\sim}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 $S(L)\ensuremath{\sim}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.

One of the most straightforward ways to study thermal properties beyond linear response is to monitor the relaxation of an arbitrarily large left-right temperature gradient ${T}_{L}\ensuremath{-}{T}_{R}$. In one-dimensional systems which support ballistic thermal transport, the local energy currents $\ensuremath{\langle}j(t)\ensuremath{\rangle}$ acquire a nonzero value at long times, and it was recently investigated whether or not this steady state fulfills a simple additive relation $\ensuremath{\langle}j(t\ensuremath{\rightarrow}\ensuremath{\infty})\ensuremath{\rangle}=f({T}_{L})\ensuremath{-}f({T}_{R})$ in integrable models. In this paper, we probe the nonequilibrium dynamics of the Hubbard chain using density matrix renormalization group (DMRG) numerics. We show that the above form provides an effective description of thermal transport in this model; violations are below the finite-time accuracy of the DMRG. As a second setup, we study how an initially equilibrated system radiates into different nonthermal states (such as the vacuum).

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 within any standard matrix product state (MPS) based formulation of the method. Most importantly, we show how to ‘combine’ the Schrödinger and Heisenberg time evolutions of arbitrary pure states |ψ〉|ψ〉 and operators AA in the evaluation of 〈A〉ψ(t)=〈ψ|A(t)|ψ〉〈A〉ψ(t)=〈ψ|A(t)|ψ〉. This includes quantum quenches. The generalization to (non-)thermal mixed state dynamics 〈A〉ρ(t)=Tr[ρA(t)]〈A〉ρ(t)=Tr[ρA(t)] induced by an initial density matrix ρρ is straightforward. In the context of linear response (ground state or finite temperature T>0T>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>0T>0. Finally, we advocate the python programming language as an elegant option for beginners to set up a DMRG code.

Integrable models such as the spin-1/2 Heisenberg chain, the Lieb-Liniger, or the one-dimensional Hubbard model are known to avoid thermalization, which was also demonstrated in several quantum-quench experiments. Another dramatic consequence of integrability is the zero-frequency anomaly in transport coefficients, which results in ballistic finite-temperature transport, despite the presence of strong interactions. While this aspect of nonergodic dynamics has been known for a long time, there has so far not been any unambiguous experimental realization thereof. We make a concrete proposal for the observation of ballistic transport via local quantum-quench experiments in fermionic quantum-gas microscopes. Such an experiment would also unveil the coexistence of ballistic and diffusive transport channels in one and the same system and provide a means of measuring finite-temperature Drude weights. The connection between local quenches and linear-response functions is established via time-dependent Einstein relations.

The rates at which energy and particle densities move to equalize arbitrarily large temperature and chemical potential differences in an isolated quantum system have an emergent thermodynamical description whenever energy or particle current commutes with the Hamiltonian. Concrete examples include the energy current in the 1D spinless fermion model with nearest-neighbor interactions (XXZ spin chain), energy current in Lorentz-invariant theories or particle current in interacting Bose gases in arbitrary dimension. Even far from equilibrium, these rates are controlled by state functions, which we call ``expansion potentials'', expressed as integrals of equilibrium Drude weights. This relation between nonequilibrium quantities and linear response implies non-equilibrium Maxwell relations for the Drude weights. We verify our results via DMRG calculations for the XXZ chain., v2: to appear in PRL

We study the spin and thermal conductivity of spin-$\frac{1}{2}$ ladders and chains at finite temperature, 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-$\frac{1}{2}$ $XXZ$ chains, for which our analysis suggests nonzero 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\ensuremath{\sim}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 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 generalize nonlinear Luttinger liquid theory to describe the dynamics of one-dimensional quantum critical systems at low temperatures. Analyzing density-matrix renormalization group results for the spin autocorrelation function in the XXZ chain we provide, in particular, direct evidence for spin diffusion in sharp contrast to the exponential decay in time predicted by conventional Luttinger liquid theory. Furthermore, we discuss how the frequencies and exponents of the oscillatory contributions from the band edges are renormalized by irrelevant interactions and obtain excellent agreement between our finite temperature nonlinear Luttinger liquid theory and the numerical data.

We study the interplay of interactions and disorder in a one-dimensional fermion lattice coupled adiabatically to infinite reservoirs. We employ both the functional renormalization group (FRG) as well as matrix product state techniques, which serve as an accurate benchmark for small systems. Using the FRG, we compute the length- and temperature-dependence of the conductance averaged over $10^4$ samples for lattices as large as $10^{5}$ sites. We identify regimes in which non-ohmic power law behavior can be observed and demonstrate that the corresponding exponents can be understood by adapting earlier predictions obtained perturbatively for disordered Luttinger liquids. In presence of both disorder and isolated impurities, the conductance has a universal single-parameter scaling form. This lays the groundwork for an application of the functional renormalization group to the realm of many-body localization.

Using a combination of numerically exact and renormalization-group techniques we study the nonequilibrium transport of electrons in an one-dimensional interacting system subject to a quasiperiodic potential. For this purpose we calculate the growth of the mean-square displacement as well as the melting of domain walls. While the system is nonintegrable for all studied parameters, there is no on finite region default of parameters for which we observe diffusive transport. In particular, our model shows a rich dynamical behavior crossing over from superdiffusion to subdiffusion. We discuss the implications of our results for the general problem of many-body localization, with a particular emphasis on the rare region Griffiths picture of subdiffusion., Comment: 6 pages, 5 figures. A more detailed analysis of the dynamical exponents extraction and discussion of the relevant times. Adds a log-derivative for the FRG section

We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is twofold: 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 ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ in the 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 test bed for our numerics, we compute ${\ensuremath{\sigma}}_{r}eg(\ensuremath{\omega})$ for the integrable XXZ spin chain in its gapped phase.

We investigate sudden quenches across the critical point in the transverse field Ising chain with a perturbing non-integrable next-nearest-neighbour interaction. Expressions for the return (Loschmidt) amplitude and associated rate function are derived to linear order in the next-nearest-neighbour coupling. In the thermodynamic limit these quantities exhibit non-analytic behaviour at a set of critical times, a phenomenon referred to as a dynamical quantum phase transition. We quantify the effect of the integrability breaking perturbation on the location and shape of these non-analyticities. Our results agree with those of earlier numerical studies and offer further support for the assertion that the dynamical quantum phase transitions exhibited by this model are a generic feature of its post-quench dynamics and is robust with respect to the inclusion of non-integrable perturbations., 19 pages, 5 figures

We introduce a numerical algorithm to simulate the time evolution of a matrix product state under a long-ranged Hamiltonian. In the effectively one-dimensional representation of a system by matrix product states, long-ranged interactions are necessary to simulate not just many physical interactions but also higher-dimensional problems with short-ranged interactions. Since our method overcomes the restriction to short-ranged Hamiltonians of most existing methods, it proves particularly useful for studying the dynamics of both power-law interacting one-dimensional systems, such as Coulombic and dipolar systems, and quasi two-dimensional systems, such as strips or cylinders. First, we benchmark the method by verifying a long-standing theoretical prediction for the dynamical correlation functions of the Haldane-Shastry model. Second, we simulate the time evolution of an expanding cloud of particles in the two-dimensional Bose-Hubbard model, a subject of several recent experiments., 5 pages + 3 pages appendices, 4 figures