Search

The {\it exchange} interaction arising from the particle indistinguishability is of central importance to physics of many-particle quantum systems. Here we study analytically the dynamical generation of quantum entanglement induced by this interaction in an isolated system, namely, an ideal Fermi gas confined in a chaotic cavity, which evolves unitarily from a non-Gaussian pure state. We find that the breakdown of the quantum-classical correspondence of particle motion, via dramatically changing the spatial structure of many-body wavefunction, leads to profound changes of the entanglement structure. Furthermore, for a class of initial states, such change leads to the approach to thermal equilibrium everywhere in the cavity, with the well-known Ehrenfest time in quantum chaos as the thermalization time. Specifically, the quantum expectation values of various correlation functions at different spatial scales are all determined by the Fermi-Dirac distribution. In addition, by using the reduced density matrix (RDM) and the entanglement entropy (EE) as local probes, we find that the gas inside a subsystem is at equilibrium with that outside, and its thermal entropy is the EE, even though the whole system is in a pure state. As a by-product of this work, we provide an analytical solution supporting an important conjecture on thermalization, made and numerically studied by Garrison and Grover in: Phys. Rev. X \textbf{8}, 021026 (2018), and strengthen its statement., Comment: 31 pages, 4 figures

In most quantum sensing schemes, interactions between the constituent particles of the sensor are expected to lead to thermalisation and degraded sensitivity. However, recent theoretical and experimental work has shown that the phenomenon of quantum many-body scarring can slow down, or even prevent thermalisation. We show that scarring can be exploited for quantum sensing that is robust against certain strong interactions. In the ideal case of perfect scars with harmonic energy gaps, the optimal sensing time can diverge despite the strong interactions. We demonstrate the idea with two examples: a spin-1 model with Dzyaloshinskii-Moriya interaction, and a spin-1/2 mixed-field Ising model. We also briefly discuss some non-ideal perturbations, and the addition of periodic controls to suppress their effect on sensing., Comment: 8 + 4 pages

We review recent progress in understanding nearly integrable models within the framework of generalized hydrodynamics (GHD). Integrable systems have infinitely many conserved quantities and stable quasiparticle excitations: when integrability is broken, only a few residual conserved quantities survive, eventually leading to thermalization, chaotic dynamics and conventional hydrodynamics. In this review, we summarize recent efforts to take into account small integrability breaking terms, and describe the transition from GHD to standard hydrodynamics. We discuss the current state of the art, with emphasis on weakly inhomogeneous potentials, generalized Boltzmann equations and collision integrals, as well as bound-state recombination effects. We also identify important open questions for future works., Comment: Review paper, 57 pages, 7 figures. Published in JSTAT special issue 2021 "Emergent Hydrodynamics in Integrable Many-body Systems"

Interaction between quasiparticles of a different nature, such as magnons and phonons in a magnetic medium, leads to the mixing of their properties and the formation of hybrid states in the areas of intersection of individual spectral branches. We recently reported the discovery of a new phenomenon mediated by the magnon-phonon interaction: the spontaneous bottleneck accumulation of magneto-elastic bosons under electromagnetic pumping of pure magnons into a ferrimagnetic yttrium iron garnet film. Here, by studying the transport properties of the accumulated magneto-elastic bosons, we reveal that such accumulation occurs in two frequency-distant groups of quasiparticles: quasi-phonons and quasi-magnons. They propagate with different speeds in different directions relative to the magnetization field. The theoretical model we propose qualitatively describes the double accumulation effect, and the analysis of the two-dimensional spectrum of quasiparticles in the hybridization region allows us to determine the wavevectors and frequencies of each of the groups.

Under the influence of external environments, quantum systems can undergo various different processes, including decoherence and equilibration. We observe that macroscopic objects are both objective and thermal, thus leading to the expectation that both objectivity and thermalisation can peacefully coexist on the quantum regime too. Crucially, however, objectivity relies on distributed classical information that could conflict with thermalisation. Here, we examine the overlap between thermal and objective states. We find that in general, one cannot exist when the other is present. However, there are certain regimes where thermality and objectivity are more likely to coexist: in the high temperature limit, at the non-degenerate low temperature limit, and when the environment is large. This is consistent with our experiences that everyday-sized objects can be both thermal and objective., Comment: Published version: 17 + 3 pages, 4 figures, 1 table

We study an expanding two-fluid model of non-relativistic dark matter and radiation which are allowed to interact during a certain time span and to establish an approximate thermal equilibrium. Such interaction which generates an effective bulk viscous pressure at background level is expected to be relevant for times around the transition from radiation to matter dominance. We quantify the magnitude of this pressure for dark matter particles masses within the range $1 {\rm eV} \lesssim m_{\chi} \lesssim 10 {\rm eV}$ around the matter-radiation equality epoch (i.e., redshift $z_{\rm eq}\sim 3400$) and demonstrate that the existence of a transient bulk viscosity has consequences which may be relevant for addressing current tensions of the standard cosmological model: i) the additional (negative) pressure contribution modifies the expansion rate around $z_{\rm eq}$, yielding a larger $H_0$ value and ii) large-scale structure formation is impacted by suppressing the amplitude of matter overdensity growth via a new viscous friction-term contribution to the M\'esz\'aros effect. As a result, both the $H_0$ and the $S_8$ tensions of the current standard cosmological model are significantly alleviated., Comment: 8 pages/5 figures. V.2: version accepted for publication in PRD

We report the observation of long-lived Floquet prethermal states in a bulk solid composed of dipolar-coupled $^{13}$C nuclei in diamond at room temperature. For precessing nuclear spins prepared in an initial transverse state, we demonstrate pulsed spin-lock Floquet control that prevents their decay over multiple-minute long periods. We observe Floquet prethermal lifetimes $T_2'\approx$90.9s, extended >60,000-fold over the nuclear free induction decay times. The spins themselves are continuously interrogated for $\sim$10min, corresponding to the application of $\approx$5.8M control pulses. The $^{13}$C nuclei are optically hyperpolarized by lattice Nitrogen Vacancy (NV) centers; the combination of hyperpolarization and continuous spin readout yields significant signal-to-noise in the measurements. This allows probing the Floquet thermalization dynamics with unprecedented clarity. We identify four characteristic regimes of the thermalization process, discerning short-time transient processes leading to the prethermal plateau, and long-time system heating towards infinite temperature. This work points to new opportunities possible via Floquet control in networks of dilute, randomly distributed, low-sensitivity nuclei. In particular, the combination of minutes-long prethermal lifetimes and continuous spin interrogation opens avenues for quantum sensors constructed from hyperpolarized Floquet prethermal nuclei., 5 pages, 5 figures. SI: 2 pages, 4 figures

Microcavity exciton-polaritons, known to exhibit non-equilibrium Bose condensation at high critical temperatures, can be also brought in thermal equilibrium with the surrounding medium and form a quantum degenerate Bose-Einstein distribution. It happens when their thermalization time in the regime of positive detunings -- or, alternatively, for high-finesse microcavities -- becomes shorter than their lifetime. Here we present the self-consistent finite-temperature Hartree-Fock-Bogoliubov description for such a system of polaritons, universally addressing the excitation spectrum, momentum-dependent interactions, condensate depletion, and the background population of dark excitons that contribute to the system's chemical potential. Employing the derived expressions, we discuss the implications for the Bogoliubov sound velocity, confirmed by existing experiments, and define the critical temperatures of (quasi-)condensation and the integral particle lifetime dependencies on the detuning. Large positive detunings are shown to provide conditions for the total lifetime reaching nanosecond timescales. This allows realization of thermodynamically-equilibrium polariton systems with Bose-Einstein condensate forming at temperatures as high as tens of Kelvin., Comment: 13 pages, 8 figures

Analyses of thermal diffusivity data on complex insulators and on strongly correlated electron systems hosted in similar complex crystal structures suggest that quantum chaos is a good description for thermalization processes in these systems, particularly in the high-temperature regime where the many phonon bands and their interactions dominate the thermal transport. Here we observe that for these systems diffusive thermal transport is controlled by a universal Planckian timescale τ ∼ ℏ / k B T and a unique velocity v E . Specifically, v E ≈ v p h for complex insulators, and v p h ≲ v E ≪ v F in the presence of strongly correlated itinerant electrons ( v p h and v F are the phonon and electron velocities, respectively). For the complex correlated electron systems we further show that charge diffusivity, while also reaching the Planckian relaxation bound, is largely dominated by the Fermi velocity of the electrons, hence suggesting that it is only the thermal (energy) diffusivity that describes chaos diffusivity.

The generalized amplitude damping (GAD) quantum channel implements the interaction between a qubit and an environment with arbitrary temperature and arbitrary interaction time. Here, we implement a photonic version of the GAD for the case of infinite interaction time (full thermalization). We also show that this quantum channel works as a thermal bath with controlled temperature., Comment: 5 pages, 3 figures

We investigate the relaxation dynamics of open non-integrable quantum many-body systems in the thermodynamic limit by using a tensor-network formalism. We simulate the Lindblad quantum master equation (LQME) of infinite systems by making use of the uniform matrix product operators (MPO) as the ansatz of their density matrices. Furthermore, we establish a method to measure the thermodynamic equivalence between two states described by the uniform MPOs. We numerically show that when an initial state of the LQME is a thermal Gibbs state, a time evolved state is always indistinguishable from a Gibbs state with a time-dependent effective temperature in the weak-dissipation and thermodynamic limit.

We consider weakly ionized plasma where frequent elastic scattering of electrons on neutrals changes the individual acts and the rate of electron–electron collisions significantly. In this case, the kinetics of electron thermalization is very different from that in fully ionized plasma. The colliding electrons do not move freely. They, instead, diffuse because of fast scattering on neutrals. We demonstrate how a proper account of this diffusion enables one to estimate the characteristic time of electron thermalization. We also present a rigorous derivation of the kinetic equation for electrons by using Bogolyubov's method based on Liouville equations for multi-particle distribution functions.

Relaxation to a thermal state is the inevitable fate of non-equilibrium interacting quantum systems without special conservation laws. While thermalization in one-dimensional (1D) systems can often be suppressed by integrability mechanisms, in two spatial dimensions thermalization is expected to be far more effective due to the increased phase space. In this work we propose a general framework for escaping or delaying the emergence of the thermal state in two-dimensional (2D) arrays of Rydberg atoms via the mechanism of quantum scars, i.e. initial states that fail to thermalize. The suppression of thermalization is achieved in two complementary ways: by adding local perturbations or by adjusting the driving Rabi frequency according to the local connectivity of the lattice. We demonstrate that these mechanisms allow to realize robust quantum scars in various two-dimensional lattices, including decorated lattices with non-constant connectivity. In particular, we show that a small decrease of the Rabi frequency at the corners of the lattice is crucial for mitigating the strong boundary effects in two-dimensional systems. Our results identify synchronization as an important tool for future experiments on two-dimensional quantum scars., Comment: 12 pages, 10 figures

Select nuclear spins in a solid isolated from the rest due to local fields reach equilibrium via multielectron spin interactions., Disorder and many body interactions are known to impact transport and thermalization in competing ways, with the dominance of one or the other giving rise to fundamentally different dynamical phases. Here we investigate the spin diffusion dynamics of 13C in diamond, which we dynamically polarize at room temperature via optical spin pumping of engineered color centers. We focus on low-abundance, strongly hyperfine-coupled nuclei, whose role in the polarization transport we expose through the integrated impact of variable radio-frequency excitation on the observable bulk 13C magnetic resonance signal. Unexpectedly, we find good thermal contact throughout the nuclear spin bath, virtually independent of the hyperfine coupling strength, which we attribute to effective carbon-carbon interactions mediated by the electronic spin ensemble. In particular, observations across the full range of hyperfine couplings indicate the nuclear spin diffusion constant takes values up to two orders of magnitude greater than that expected from homo-nuclear spin couplings.

Quantum many-body scars have attracted much interest as a violation of the eigenstate thermalization hypothesis (ETH) due to recent experimental observation in Rydberg atoms and related theoretical studies. In this paper, we construct a model hosting exact Floquet quantum many-body scars, which violate the Floquet version of ETH. We consider two uniformly-driven static Hamiltonians prohibiting neighboring up spins (Rydberg blockade) like the PXP model, and construct a binary drive composed of them. We show that there exists a four-dimensional subspace which completely avoids thermalization to infinite temperature and that any other states, including some special scar states reported in the static PXP model, are vulnerable to heating and relax to infinite temperature. We also construct a more generalized periodic drive composed of time-dependent PXP-type Hamiltonians showing exact Floquet quantum many-body scars and discuss possible experimental realization of the model in Rydberg atoms., Comment: 17 pages, 4 figures

Monolayers of transition metal dichalcogenides (TMDs) have been established in the last years as promising materials for novel optoelectronic devices. However, the performance of such devices is often limited by the dissociation of tightly bound excitons into free electrons and holes. While previous studies have investigated tunneling at large electric fields, we focus in this work on phonon-assisted exciton dissociation that is expected to be the dominant mechanism at small fields. We present a microscopic model based on the density matrix formalism providing access to time- and momentum-resolved exciton dynamics including phonon-assisted dissociation. We track the pathway of excitons from optical excitation via thermalization to dissociation, identifying the main transitions and dissociation channels. Furthermore, we find intrinsic limits for the quantum efficiency and response time of a TMD-based photodetector and investigate their tunability with externally accessible knobs, such as excitation energy, substrate screening, temperature and strain. Our work provides microscopic insights in fundamental mechanisms behind exciton dissociation and can serve as a guide for the optimization of TMD-based optoelectronic devices.

Long-lived dark states, in which an experimentally accessible qubit is not in thermal equilibrium with a surrounding spin bath, are pervasive in solid-state systems. We explain the ubiquity of dark states in a large class of inhomogenous central spin models using the proximity to integrable lines with exact dark eigenstates. At numerically accessible sizes, dark states persist as eigenstates at large deviations from integrability, and the qubit retains memory of its initial polarization at long times. Although the eigenstates of the system are chaotic, exhibiting exponential sensitivity to small perturbations, they do not satisfy the eigenstate thermalization hypothesis. Rather, we predict long relaxation times that increase exponentially with system size. We propose that this intermediate chaotic but non-ergodic regime characterizes mesoscopic quantum dot and diamond defect systems, as we see no numerical tendency towards conventional thermalization with a finite relaxation time., Comment: 14 pages, 15 figures

A promising probe to unmask particle dark matter is to observe its effect on neutron stars, the prospects of which depend critically on whether captured dark matter thermalizes in a timely manner with the stellar core via repeated scattering with the Fermi-degenerate medium. In this work we estimate the timescales for thermalization for multiple scenarios. These include: (a) spin-0 and spin-$\frac{1}{2}$ dark matter, (b) scattering on non-relativistic neutron and relativistic electron targets accounting for the respective kinematics, (c) interactions via a range of Lorentz-invariant structures, (d) mediators both heavy and light in comparison to the typical transfer momenta in the problem. We discuss the analytic behavior of the thermalization time as a function of the dark matter and mediator masses, and the stellar temperature. Finally, we identify parametric ranges where both stellar capture is efficient and thermalization occurs within the age of the universe. For dark matter that can annihilate in the core, these regions indicate parametric ranges that can be probed by upcoming infrared telescopes observing cold neutron stars., Comment: 12 pages revtex4, 7 figures, 2 tables; v2: references added, minor typos fixed

We use fast periodic control to realize finite-time Otto cycles exhibiting quantum advantage. Such periodic modulation of the working medium - bath interaction Hamiltonian during the thermalization strokes can give rise to non-Markovian anti-Zeno dynamics, and corresponding reduction in the thermalization times. Faster thermalization can in turn significantly enhance the power output in engines, or equivalently, the rate of refrigeration in refrigerators. This improvement in performance of dynamically controlled Otto thermal machines arises due to the time-energy uncertainty relation of quantum mechanics., Comment: 12 pages, 13 figures. Close to published version

We introduce a diagnostic for quantum thermalization based on mixed-state entanglement. Specifically, given a pure state on a tripartite system $ABC$, we study the scaling of entanglement negativity between $A$ and $B$. For representative states of self-thermalizing systems, either eigenstates or states obtained by a long-time evolution of product states, negativity shows a sharp transition from an area-law scaling to a volume-law scaling when the subsystem volume fraction is tuned across a finite critical value. In contrast, for a system with quasiparticles, it exhibits a volume-law scaling irrespective of the subsystem fraction. For many-body localized systems, the same quantity shows an area-law scaling for eigenstates, and volume-law scaling for long-time evolved product states, irrespective of the subsystem fraction. We provide a combination of numerical observations and analytical arguments in support of our conjecture. Along the way, we prove and utilize a `continuity bound' for negativity: we bound the difference in negativity for two density matrices in terms of the Hilbert-Schmidt norm of their difference., Comment: 25 pages, 15 figures; typos fixed, references added

We study the effect of a quantum quench between two tunnel coupled Tomonaga-Luttinger liquids (TLLs) with different speed of sound and interaction parameter. The quench dynamics is induced by switching off the tunnelling and letting the two systems evolve independently. We fully diagonalize the problem within a quadratic approximation for the initial tunnelling. Both the case of zero and finite temperature in the initial state are considered. We focus on correlation functions associated with the antisymmetric and symmetric combinations of the two TLLs (relevant for interference measurements), which turn out to be coupled due to the asymmetry in the two systems' Hamiltonians. The presence of different speeds of sound leads to multiple lightcones separating different decaying regimes. In particular, in the large time limit, we are able to identify a prethermal regime where the two-point correlation functions of vertex operators of symmetric and antisymmetric sector can be characterized by two emerging effective temperatures, eventually drifting towards a final stationary regime that we dubbed $quasi-thermal$, well approximated at large scale by a thermal-like state, where these correlators become time independent and are characterized by a unique correlation length. If the initial state is at equilibrium at non-zero temperature $T_0$, all the effective temperatures acquire a linear correction in $T_0$, leading to faster decay of the correlation functions. Such effects can play a crucial role for the correct description of currently running cold atoms experiments., Comment: Accepted version: title and discussion section modified, one appendix removed

The wave turbulence theory predicts that a conservative system of nonlinear waves can exhibit a process of condensation, which originates in the singularity of the Rayleigh-Jeans equilibrium distribution of classical waves. Considering light propagation in a multimode fiber, we show that light condensation is driven by an energy flow toward the higher-order modes, and a bi-directional redistribution of the wave-action (or power) to the fundamental mode and to higher-order modes. The analysis of the near-field intensity distribution provides experimental evidence of this mechanism. The kinetic equation also shows that the wave-action and energy flows can be inverted through a thermalization toward a negative temperature equilibrium state, in which the high-order modes are more populated than low-order modes. In addition, a Bogoliubov stability analysis reveals that the condensate state is stable., Comment: 3 figures

The early Universe may have passed through an extended period of matter-dominated expansion following inflation and prior to the onset of radiation domination. Sub-horizon density perturbations grow gravitationally during such an epoch, collapsing into bound structures if it lasts long enough. The strong analogy between this phase and structure formation in the present-day universe allows the use of N-body simulations and approximate methods for halo formation to model the fragmentation of the inflaton condensate into inflaton halos. For a simple model we find these halos have masses of up to $20\,\mathrm{kg}$ and radii of the order of $10^{-20}\,\mathrm{m}$, roughly $10^{-24}$ seconds after the Big Bang. We find that the N-body halo mass function matches predictions of the mass-Peak Patch method and the Press-Schechter formalism within the expected range of scales. A long matter-dominated phase would imply that reheating and thermalization occurs in a universe with large variations in density, potentially modifying the dynamics of this process. In addition, large overdensities can source gravitational waves and may lead to the formation of primordial black holes., Comment: 12 pages, 8 figures, accepted by PRD

Cosmological evolution and particle creation in $R^2$-modified gravity are considered for the case of the dominant decay of the scalaron into a pair of gauge bosons due to conformal anomaly. It is shown that in the process of thermalization superheavy dark matter with the coupling strength typical for the GUT SUSY can be created. Such dark matter would have the proper cosmological density if the particle mass is close to $10^{12}$ GeV., Comment: 16 pages, 11 figures. Some improvements have been done based on referee's comment

We employ the holographic approach to study the thermalization in the quenched strongly-coupled field theories with very general anisotropic scalings including Lifshitz and hyperscaling violating fixed points. The holographic dual is a Vaidya-like time-dependent geometry where the asymptotic metric has general anisotropic scaling isometries. We find the Ryu-Takanayagi extremal surface and use it to calculate the time-dependent entanglement entropy between a strip region with width $2R$ and its outside region. In the special case with an isotropic metric, we also explore the entanglement entropy for a spherical region of radius $R$. The growth of the entanglement entropy characterizes the thermalization rate after a quench. We study the thermalization process in the early times and late times in both large $R$ and small $R$ limits. The allowed scaling parameter regions are constrained by the null energy conditions as well as the condition for the existence of the Ryu-Takanayagi extremal surfaces. This generalizes the previous works on this subject. All obtained results can be compared with experiments and other methods of probing thermalization., Comment: The published version

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (non-integrable) discrete nonlinear Schr{\"o}dinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS towards AL) is varied, the region in the space of initial energy and norm-densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes., Comment: 11 figures

We theoretically study the out-of-equilibrium dynamics in momentum space of a weakly interacting disordered Bose gas launched with a finite velocity. In the absence of interactions, coherent multiple scattering gives rise to a background of diffusive particles, on top of which a coherent backscattering interference emerges. We revisit this scenario in the presence of interactions, using a diagrammatic quantum transport theory. We find that the dynamics is governed by coupled kinetic equations describing the thermalization of the diffusive and coherent components of the gas. This phenomenon leads to a destruction of coherent backscattering, well described by an exponential relaxation whose rate is controlled by the particle collision time. These predictions are confirmed by numerical simulations., Comment: 16 pages, 12 figures

I review recent progress in thermalization in heavy-ion collisions, with particular emphasis on hydrodynamic attractor, and also report recent progress in hydrodynamic fluctuations., Comment: 8 pages, 3 figures, Proceedings for Quark Matter 2019

In this manuscript, it is presented an overview of the Quark-Gluon Plasma properties measured, so far, using hard probes. We will focus on both quantitative and qualitative properties that have been (or are about to be) measured, making a link between the theoretical description and experimental results. Throughout the manuscript, highlights to some of the conferences' results will be given, but without an extensive overview. A personal opinion of the most important developments and critical problems that need more work in the future is presented at the end., Comment: [v1:] 8 pages, 5 figures, Quark Matter 2019 Proceedings; [v2:] Reference added; Typos corrected

Photon thermalisation and condensation in dye-filled microcavities is a growing area of scientific interest, at the intersection of photonics, quantum optics and statistical physics. We give here a short introduction to the topic, together with an explanation of some of our more important recent results. A key result across several projects is that we have a model based on a detailed physical description which has been used to accurately describe experimental observations. We present a new open-source package in Python called PyPBEC which implements this model. The aim is to enable the reader to readily simulate and explore the physics of photon condensates themselves, so this article also includes a working example code which can be downloaded from the GitHub repository., Comment: This manuscript accompanies the software repository https://github.com/photonbec/PyPBEC and has been submitted to the conference proceedings of Photon 2020

Thermal diffusivity measurements of samples transmitting thermal radiation require adjustments to the data treatment procedures in laser flash analysis. Conventionally, an unconstrained diathermic model is used. Current results show that the alternative coupled radiative-conductive models produce substantially different results -- for instance, at high temperatures in oxide ceramics. However, care must be taken to ensure accurate implementations of each constituent computational technique. The latter are presented in this work., Comment: 39 pages, 13 figures, 6 tables

In this paper, we study the exact dynamics of general open systems interacting with its environment through particle exchanges. The paper includes two main results. First, by taking advantage of the propagating function in the coherent state representation, we solve the exact master equation, whose solution is expressed in terms of the Keldysh nonequilibrium Green functions. Second, in the dynamical perspective, we provide a rigorous thermalization process of open quantum systems., Comment: 12 pages

High temperature virial expansion is a powerful tool in equilibrium statistical mechanics. In this letter we generalize the high temperature virial expansion approach to treat far-from-equilibrium quench dynamics. As an application of our framework, we study the dynamics of a Bose gas quenched from non-interacting to unitarity, and we compare our theoretical results with unexplained experimental results by the Cambridge group [Eigen et al., Nature 563, 221 (2018)]. We show that, during the quench dynamics, the momentum distribution decreases for low-momentum part with $kk^*$, where $k^*$ is a characteristic momentum scale separating the low- and the high-momentum regimes. We determine the universal value of $k^*\lambda$ that agrees perfectly with the experiment, with $\lambda$ being the thermal de Broglie wave length. We also find a jump of the half-way relaxation time across $k^*\lambda$ and the non-monotonic behavior of energy distribution, both of which agree with the experiment. Finally, we address the issue whether the long-time steady state thermalizes or not, and we find that this state does thermalize except for the very high momentum tail with $k\lambda\gg 1$. Our framework can also be applied to quench dynamics in other systems., Comment: 6 pages, 4 figures

Prethermalization, by introducing emergent quasiconserved observables, plays a crucial role in protecting Floquet many-body phases over exponentially long time, while the ultimate fate of such quasiconserved operators can signal thermalization to infinite temperature. To elucidate the properties of prethermal quasiconservation in many-body Floquet systems, here we systematically analyze infinite temperature correlations between observables. We numerically show that the late-time behavior of the autocorrelations unambiguously distinguishes quasiconserved observables from non-conserved ones, allowing to single out a set of linearly-independent quasiconserved observables. By investigating two Floquet spin models, we identify two different mechanism underlying the quasi-conservation law. First, we numerically verify energy quasiconservation when the driving frequency is large, so that the system dynamics is approximately described by a static prethermal Hamiltonian. More interestingly, under moderate driving frequency, another quasiconserved observable can still persist if the Floquet driving contains a large global rotation. We show theoretically how to calculate this conserved observable and provide numerical verification. Having systematically identified all quasiconserved observables, we can finally investigate their behavior in the infinite-time limit and thermodynamic limit, using autocorrelations obtained from both numerical simulation and experiments in solid state nuclear magnetic resonance systems., Comment: 12 pages, 9 figures. arXiv admin note: substantial text overlap with arXiv:1912.05799

We investigate the effects of producing dark matter by Hawking evaporation of primordial black holes (PBHs) in scenarios that may have a second well-motivated dark matter production mechanism, such as freeze-out, freeze-in, or gravitational production. We show that the interplay between PBHs and the alternative sources of dark matter can give rise to model-independent modifications to the required dark matter abundance from each production mechanism, which in turn affect the prospects for dark matter detection. In particular, we demonstrate that for the freeze-out mechanism, accounting for evaporation of PBHs after freeze-out demands a larger annihilation cross section of dark matter particles than its canonical value for a thermal dark matter. For mechanisms lacking thermalization due to a feeble coupling to the thermal bath, we show that the PBH contribution to the dark matter abundance leads to the requirement of an even feebler coupling. Moreover, we show that when a large initial abundance of PBHs causes an early matter-dominated epoch, PBH evaporation alone cannot explain the whole abundance of dark matter today. In this case, an additional production mechanism is required, in contrast to the case when PBHs are formed and evaporate during a radiation-dominated epoch., Comment: 16 pages, 2 figures

A system in thermal equilibrium with a bath will generally be in an athermal state, if the system-bath coupling is strong. In some cases, it will be possible to extract work from that athermal state, after disconnecting the system from the bath. We use this observation to devise a battery charging and storing unit, simply consisting of a system, acting as the battery, and a bath. The charging cycle---connect, let thermalize, disconnect, extract work---requires very little external control and the charged state of the battery, being a part of global thermal equilibrium, can be maintained indefinitely and for free. The efficiency, defined as the ratio of the extractable work stored in the battery and the total work spent on connecting and disconnecting, is always $\leq 1$, which is a manifestation of the second law of thermodynamics. Moreover, coupling, being a resource for the device, is also a source of dissipation: the entropy production per charging cycle is always significant, strongly limiting the efficiency in all coupling strength regimes. We show that our general results also hold for generic microcanonical baths. We illustrate our theory on the Caldeira-Leggett model with a harmonic oscillator (the battery) coupled to a harmonic bath, for which we derive general asymptotic formulas in both weak and ultrastrong coupling regimes, for arbitrary Ohmic spectral densities. We show that the efficiency can be increased by connecting several copies of the battery to the bath. Finally, as a side result, we derive a general formula for Gaussian ergotropy, that is, the maximal work extractable by Gaussian unitary operations from Gaussian states of multipartite continuous-variable systems., Comment: 11+15 pages, 2 figures

If, during the early Universe epoch, the dark matter particle thermalizes in a hidden sector which does not thermalize with the Standard Model thermal bath, its relativistic thermal decoupling can easily lead to the observed relic density, even if the dark matter particle mass is many orders of magnitude heavier than the usual ∼ eV hot relic mass scale. This straightforward scenario simply requires that the temperature of the hidden sector thermal bath is one to five orders of magnitude cooler than the temperature of the Standard Model thermal bath. In this way the resulting relic density turns out to be determined only by the dark matter mass scale and the ratio of the temperatures of both sectors. In a model independent way we determine that this can work for a dark matter mass all the way from ∼1 keV to ∼30 PeV. We also show how this scenario works explicitly in the framework of two illustrative models. One of them can lead to a PeV neutrino flux from dark matter decay of the order of the one needed to account for the high energy neutrinos observed by IceCube., SCOPUS: ar.j, info:eu-repo/semantics/published

Quasi-two-dimensional transition-metal dichalcogenides are a key platform for exploring emergent nanoscale phenomena arising from complex interactions. Access to the underlying degrees-of-freedom on their natural time scales motivates the use of advanced ultrafast probes sensitive to self-organised atomic-scale patterns. Here, we report the first ultrafast investigation of TaTe2, which exhibits unique charge and lattice trimer order characterised by a transition upon cooling from stripe-like chains into a $(3 \times 3)$ superstructure of trimer clusters. Utilising MeV-scale ultrafast electron diffraction, we capture the photo-induced TaTe2 structural dynamics -- exposing a rapid $\approx\!1.4$ ps melting of its low-temperature ordered state followed by recovery via thermalisation into a hot cluster superstructure. Density-functional calculations indicate that the initial quench is triggered by intra-trimer Ta charge transfer which destabilises the clusters, unlike melting of charge density waves in other TaX2 compounds. Our work paves the way for further exploration and ultimately rapid optical and electronic manipulation of trimer superstructures., Comment: Main text and supplementary information: 30 pages, 4 main figures, 24 supplementary figures

Entanglement properties of bipartite unitary operators are studied via their local invariants, namely the entangling power $e_p$ and a complementary quantity, the gate typicality $g_t$. We characterize the boundaries of the set $K_2$ representing all two-qubit gates projected onto the plane $(e_p, g_t)$ showing that the fractional powers of the \textsc{swap} operator form a parabolic boundary of $K_2$, while the other bounds are formed by two straight lines. In this way a family of gates with extreme properties is identified and analyzed. We also show that the parabolic curve representing powers of \textsc{swap} persists in the set $K_N$, for gates of higher dimensions ($N>2$). Furthermore, we study entanglement of bipartite quantum gates applied sequentially $n$ times and analyze the influence of interlacing local unitary operations, which model generic Hamiltonian dynamics. An explicit formula for the entangling power a gate applied $n$ times averaged over random local unitary dynamics is derived for an arbitrary dimension of each subsystem. This quantity shows an exponential saturation to the value predicted by the random matrix theory (RMT), indicating "thermalization" in the entanglement properties of sequentially applied quantum gates that can have arbitrarily small, but nonzero, entanglement to begin with. The thermalization is further characterized by the spectral properties of the reshuffled and partially transposed unitary matrices., 21 pages, 8 figures

The Toda lattice is a nonlinear but integrable system. Here we study the thermalization problem in one-dimensional, perturbed Toda lattices in the thermodynamic limit. We show that the thermalization time, $T_{eq}$, follows a universal law; i.e., $T_{eq}\sim \epsilon^{-2}$, where the perturbation strength, $\epsilon$, characterizes the nonlinear perturbations added to the Toda potential. This universal law applies generally to weak nonlinear lattices due to their equivalence to perturbed Toda systems.

Long range interactions are relevant for a wide range of phenomena in physics where they often present a challenge to theory. In condensed matter, the interplay of Coulomb interaction and disorder remains largely an unsolved problem. In two dimensional films the long-range part of the Coulomb interaction may be screened by a nearby metallic overlay. This technique is employed in this work to present experimental evidence for its effectiveness in limiting the spatial range of the Coulomb interaction. We use this approach to study the effects of the long-range Coulomb interaction on the out-of-equilibrium dynamics of electron-glasses using amorphous indium-oxide films. The results demonstrate that electronic relaxation times, extending over thousands of seconds, do not hinge on the long-range Coulomb interaction nor on the presence of a real gap in the density of states. Rather, they emphasize the dominant role played by disorder in controlling the slow thermalization processes of Anderson insulators taken far from equilibrium., Comment: 9 pages, 11 figures. arXiv admin note: text overlap with arXiv:1805.07052

To explore the structure of the QCD phase diagram in high baryon density domain, several high-energy nuclear collision experiments in a wide range of beam energies are currently performed or planned using many accelerator facilities. In these experiments search for a first-order phase transition and the QCD critical point is one of the most important topics. To find the signature of the phase transition, experimental data should be compared to appropriate dynamical models which quantitatively describe the process of the collisions. In this study we develop a new dynamical model on the basis of the non-equilibrium hadronic transport model JAM and 3+1D hydrodynamics. We show that the new model reproduce well the experimental beam-energy dependence of hadron yields and particle ratio by the partial thermalization of the system in our core-corona approach., Comment: 4 pages, 3 figures; contribution to the proceedings of the 8th International Conference on Quarks and Nuclear Physics (QNP2018), Tsukuba, November 13-17, 2018

In this paper we study the thermal equilibration of small bipartite Bose-Hubbard systems, both quantum mechanically and in mean-field approximation. In particular we consider small systems composed of a single-mode ``thermometer'' coupled to a three-mode ``bath,'' with no additional environment acting on the four-mode system, and test the hypothesis that the thermometer will thermalize if and only if the bath is chaotic. We find that chaos in the bath alone is neither necessary nor sufficient for equilibration in these isolated four-mode systems. The two subsystems can thermalize if the combined system is chaotic even when neither subsystem is chaotic in isolation, and under full quantum dynamics there is a minimum coupling strength between the thermometer and the bath below which the system does not thermalize even if the bath itself is chaotic. We show that the quantum coupling threshold scales like $1/N$ (where $N$ is the total particle number), so that the classical results are obtained in the limit $N\ensuremath{\rightarrow}\ensuremath{\infty}$.

We present new results on the calculation of the dark matter relic abundance within the Higgs induced right-handed (RH) neutrino mixing model, solving the associated density matrix equation. For a benchmark value of the dark matter mass $M_{\rm DM} = 220\,{\rm TeV}$, we show the evolution of the abundance and how this depends on reheat temperature, dark matter lifetime and source RH neutrino mass $M_{\rm S}$, with the assumption $M_{\rm S} < M_{\rm DM}$. We compare the results with those obtained within the Landau-Zener approximation, showing that the latter largely overestimates the final abundance giving some analytical insight. However, we also notice that since in the density matrix formalism the production is non-resonant, this allows source RH neutrino masses below the W boson mass, making dark matter more stable at large mass values. This opens an allowed region for initial vanishing source RH neutrino abundance. For example, for $M_{\rm S} \gtrsim 1\,{\rm GeV}$, we find $M_{\rm DM}\gtrsim 20\,{\rm PeV}$. Otherwise, for $M_{\rm S} > M_W\sim 100\,{\rm GeV}$, one has to assume a thermalisation of the source RH neutrinos prior to the freeze-in of the dark matter abundance. This results into a large allowed range for $M_{\rm DM}$, depending on $M_{\rm S}$. For example, imposing $M_{\rm S} \gtrsim 300\,{\rm GeV}$, allowing also successful leptogenesis, we find $0.5 \lesssim M_{\rm DM}/{\rm TeV} \lesssim 500$. We also discuss in detail leptogenesis with two quasi-degenerate RH neutrinos, showing a case when observed dark matter abundance and matter-antimatter asymmetry are simultaneously reproduced. Finally, we comment on how an initial thermal source RH neutrino abundance can be justified and on how our results suggest that also the interesting case where $M_{\rm DM} < M_{\rm S}$, embeddable in usual high scale two RH neutrino seesaw models, might be viable., 46 pages, 7 figures; v2: appendices B and C added, it matches version to be published in JCAP

Heavy quarks are effective probes of the properties of the QGP created in ultra-relativistic heavy-ion collisions. The ALICE Collaboration measured the non-strange D-meson production in Pb--Pb collisions at $\sqrt{s_{\rm NN\,}}=5.02$ TeV. The in-medium energy loss can be studied via the nuclear modification factor measurement. The measurement of the D-meson elliptic flow, $v_2$, allows us to investigate the participation of the heavy quarks in the collective expansion of the system at low momentum and their possible thermalization in the medium. Furthermore, the Event-Shape Engineering technique is used to measure D-meson elliptic flow in order to study the coupling of the charm quarks to the light quarks of the underlying medium., Comment: 5 pages, 6 figures, The 18th International Conference on Strangeness in Quark Matter (SQM 2019)

We show that a quantum dynamical localization effect can be observed in a generic thermalization process of two weakly-coupled chaotic subsystems. Specifically, our model consists of the minimal experimentally relevant subsystems that exhibit chaos, which are 3-site Bose-Hubbard units. Due to the high dimensionality of the composite 6-site system, the quantum localization effect is weak and can not be resolved merely by the breakdown of quantum-to-classical correspondence. Instead, we adopt an intrinsic definition of localization as the memory of initial conditions, that is not related to the underlying classical dynamics. We discuss the dynamics in the chaotic sea, and in the vicinity of the mobility edge, beyond which ergodization is suppressed., Comment: 15 pages, 10 figures, elaborated version, extra paragraph. This work is dedicated to Shmuel Fishman who passed away on 2-April-2019. Publication date (Phys. Rev. A) is 2-April-2020

We provide evidence that a clean kicked Bose-Hubbard model exhibits a many-body dynamically localized phase. This phase shows ergodicity breaking up to the largest sizes we were able to consider. We argue that this property persists in the limit of large size. The Floquet states violate eigenstate thermalization and then the asymptotic value of local observables depends on the initial state and is not thermal. This implies that the system does not generically heat up to infinite temperature, for almost all the initial states. Differently from many-body localization here the entanglement entropy linearly increases in time. This increase corresponds to space-delocalized Floquet states which are nevertheless localized across specific subsectors of the Hilbert space: In this way the system is prevented from randomly exploring all the Hilbert space and does not thermalize., Comment: 18 pages, 20 figures, revised version published in PRB

As recently manifested , the quench dynamics of isolated quantum systems consisting of a finite number of particles, is characterized by an exponential spreading of wave packets in the many-body Hilbert space. This happens when the inter-particle interaction is strong enough, thus resulting in a chaotic structure of the many-body eigenstates considered in an unperturbed basis. The semi-analytical approach used here, allows one to estimate the rate of the exponential growth as well as the relaxation time, after which the equilibration (thermalization) emerges. The key ingredient parameter in the description of this process is the width $\Gamma$ of the Local Density of States (LDoS) defined by the initially excited state, the number of particles and the interaction strength. In this paper we show that apart from the meaning of $\Gamma$ as the decay rate of survival probability, the width of the LDoS is directly related to the diagonal entropy and the latter can be linked to the thermodynamic entropy of a system equilibrium state emerging after the complete relaxation. The analytical expression relating the two entropies is derived phenomenologically and numerically confirmed in a model of bosons with random two-body interaction, as well as in a deterministic model which becomes completely integrable in the continuous limit., Comment: The original version has been divided in two separate papers, this is the first one

The phenomenon of chiral-symmetry protected thermalization gap in Hermitian photonic systems is counterintuitive as it implies that the photon coherence can be continuously improved by disorders towards an asymptotic limit. We show that the phenomenon disappears in time-independent, non-Hermitian photonic systems even when the chiral symmetry is well preserved. In fact, the degree of thermalization generally increases with the disorder strength, in agreement with intuition. As non-Hermitian characteristics (e.g., weak gain and loss) can be expected in realistic physical situations, the phenomenon of thermalization gap may be observed but only in well controlled experiments with high quality materials., Comment: 6 pages, 3 figures

International audience; Many systems, including biological tissues and foams, are made of highly packed units having high deformability but low compressibility. At two dimensions, these systems offer natural tesselations of plane with fixed density, in which transitions from ordered to disordered patterns are often observed, in both directions. Using a modified Cellular Potts Model algorithm that allows rapid thermalization of extensive systems, we numerically explore the order-disorder transition of monodisperse, two-dimensional cellular systems driven by thermal agitation. We show that the transition follows most of the predictions of Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory developed for melting of 2D solids, extending the validity of this theory to systems with many-body interactions. In particular, we show the existence of an intermediate hexatic phase, which preserves the orientational order of the regular hexagonal tiling, but looses its positional order. In addition to shedding light on the structural changes observed in experimental systems, our study shows that soft cellular systems offer macroscopic systems in which KTHNY melting scenario can be explored, in the continuation of Bragg's experiments on bubble rafts.