Source author record

D. M. Kennes

D. M. Kennes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2021arXiv

Anti-Poiseuille Flow: Increased Vortex Velocity at Superconductor Edges

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.

preprint2021arXiv

Vortex control in superconducting Corbino geometry networks

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.

preprint2020arXiv

Influence of spin and orbital fluctuations on Mott-Hubbard exciton dynamics in LaVO${}_{3}$ Thin Films

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.

preprint2020arXiv

Loschmidt-amplitude wave function spectroscopy and the physics of dynamically driven phase transitions

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.

preprint2020arXiv

One-dimensional flat bands in twisted bilayer germanium selenide

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é 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.

preprint2020arXiv

The interacting Rice-Mele model: bulk and boundaries

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.

preprint2019arXiv

Absence of dissipationless transport in clean 2D superconductors

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.

preprint2019arXiv

Entanglement and spectra in topological many-body localized phases

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.

preprint2019arXiv

Universality at work -- the local sine-Gordon model, lattice fermions, and quantum circuits

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.

preprint2016arXiv

Current noise of the interacting resonant level model

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.

preprint2016arXiv

Entanglement scaling of excited states in large one-dimensional many-body localized systems

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.

preprint2016arXiv

Renormalization in periodically driven quantum dots

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.

preprint2016arXiv

The Ohmic two-state system from the perspective of the interacting resonant level model: Thermodynamics and transient dynamics

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, $α$, 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 $α\to 1^{-}$. For the transient dynamics of the two-level system, $P(t)=\langle σ_{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α\lesssim 0.9$ for $P(t)$, and between TDNRG and FRG in the vicinity of $α=1/2$. Furthermore, we quantify the error in the NIBA for a range of $α$, finding significant errors in the latter even for $0.1\leq α\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.

preprint2016arXiv

Thermal conductivity of the one-dimensional Fermi-Hubbard model

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.

preprint2015arXiv

Spin and thermal conductivity of quantum spin ladders

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.

preprint2014arXiv

Dynamical regimes of dissipative quantum systems

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.

preprint2014arXiv

Extending the range of real time density matrix renormalization group simulations

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 <A>_psi(t)=<psi|A(t)|psi>. This includes quantum quenches. The generalization (non-)thermal mixed state dynamics <A>_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.

preprint2014arXiv

Renormalization group flows in one-dimensional lattice models: impurity scaling, umklapp scattering and the orthogonality catastrophe

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.

preprint2014arXiv

Spectral properties of one-dimensional Fermi systems after an interaction quench

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.

preprint2014arXiv

Transport properties of the one-dimensional Hubbard model at finite temperature

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.

preprint2014arXiv

Universal quench dynamics of interacting quantum impurity systems

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.

preprint2013arXiv

Efficiency and power of a thermoelectric quantum dot device

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.

preprint2013arXiv

Functional renormalization group study of the Anderson--Holstein model

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.

preprint2013arXiv

Luttinger liquid properties of the steady state after a quantum quench

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.

preprint2013arXiv

Oscillatory dynamics and non-markovian memory in dissipative quantum systems

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.

preprint2013arXiv

The interacting resonant level model in nonequilibrium: finite temperature effects

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.

preprint2011arXiv

A renormalization group approach to time dependent transport through correlated quantum dots

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.

preprint2010arXiv

Relaxation dynamics of an exactly solvable electron-phonon 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.