Eigenstate Thermalization, Random Matrix Theory, and Behemoths.

A nonlocal observation operator has been developed to assimilate GPS radio occultation (RO) refractivity with WRF 3DVAR. For simplicity, in the past GPS RO refractivity was often assimilated using a local observation operator with the assumption that the GPS RO observation was representative of a model local point. Such an operator did not take into account the effects of horizontal inhomogeneity on the derived GPS RO refractivity. In order to more accurately model the observables, Sokolovskiy et al. (2005a) developed a nonlocal observation operator; which would take into account the effects of horizontal inhomogeneity on GPS RO measurements. This nonlocal observation operator calculates the integrated amount of the model refractivity along the ray paths centered at the perigee points. For comparative purposes, the nonlocal observation operator can be simplified by limiting the length of integration near the RO point. This is called the “local operator variant”, which is equivalent to the original local operator except that the original one is performed with fixed tangent points at observation levels. For computational efficiency, assimilation using both the nonlocal operator and local operator variant now is performed with smear tangent points at the mean height of each model vertical level. In this study, the statistics of observation errors using both local and nonlocal operators were estimated based on WRF simulations. The observation errors produced by the nonlocal operator are about two times smaller than those generated by the local operator and in agreement with Sokolovskiy et al. (2005b). Each of the three operators is used to assimilate GPS RO refractivity soundings from the FORMOSAT-3/COSMIC mission using the WRF 3DVAR system. The WRF model then is applied to simulate Typhoons Kaemi (July 2006) which struck Taiwan with significant torrential rainfall. The analysis increments produced by the nonlocal operator and local operator variant are quite similar in horizontal and vertical distributions; whereas, the former is slightly stretched along the ray's direction, as a result of the longer integration length. The simulated typhoon tracks prior to landfall are quite similar for the three operators. Both the nonlocal operator and local operator variant improve the detoured track after landfall as predicted by the local operator. The nonlocal operator outperforms the two local operators in rainfall prediction at later times. The performances of the nonlocal operator in general are promising and can replace the local operator at no marked cost of computational efficiency.

We test the eigenstate thermalization hypothesis (ETH) in 1 + 1 -dimensional SU(2) lattice gauge theory (LGT) with one flavor of dynamical fermions. Using the loop-string-hadron framework of the LGT with a bosonic cutoff, we exactly diagonalize the Hamiltonian for finite size systems and calculate matrix elements (MEs) in the eigenbasis for both local and nonlocal operators. We analyze different indicators to identify the parameter space for quantum chaos at finite lattice sizes and investigate how the ETH behavior emerges in both the diagonal and off-diagonal MEs. Our investigations allow us to study various timescales of thermalization and the emergence of random matrix behavior, and highlight the interplays of the several diagnostics with each other. Furthermore, from the off-diagonal MEs, we extract a smooth function that is closely related to the spectral function for both local and nonlocal operators. We find numerical evidence of the spectral gap and the memory peak in the nonlocal operator case. Finally, we investigate aspects of subsystem ETH in the lattice gauge theory and identify certain features in the subsystem reduced density matrix that are unique to gauge theories.

Complexity of dynamics is at the core of quantum many-body chaos and exhibits a hierarchical feature: higher-order complexity implies more chaotic dynamics. Conventional ergodicity in thermalization processes is a manifestation of the lowest order complexity, which is represented by the eigenstate thermalization hypothesis (ETH) stating that individual energy eigenstates are thermal. Here, we propose a higher-order generalization of the ETH, named the $ k $-ETH ($ k=1,2,\dots $), to quantify higher-order complexity of quantum many-body dynamics at the level of individual energy eigenstates, where the lowest order ETH (1-ETH) is the conventional ETH. As a non-trivial contribution of the higher-order ETH, we show that the $ k $-ETH with $ k\geq 2 $ implies a universal behavior of the $ k $th Renyi entanglement entropy of individual energy eigenstates. In particular, the Page correction of the entanglement entropy originates from the higher-order ETH, while as is well known, the volume law can be accounted for by the 1-ETH. We numerically verify that the 2-ETH approximately holds for a nonintegrable system, but does not hold in the integrable case. To further investigate the information-theoretic feature behind the $ k $-ETH, we introduce a concept named a partial unitary $ k $-design (PU $ k $-design), which is an approximation of the Haar random unitary up to the $ k $th moment, where partial means that only a limited number of observables are accessible. The $ k $-ETH is a special case of a PU $ k $-design for the ensemble of Hamiltonian dynamics with random-time sampling. In addition, we discuss the relationship between the higher-order ETH and information scrambling quantified by out-of-time-ordered correlators. Our framework provides a unified view on thermalization, entanglement entropy, and unitary $ k $-designs, leading to deeper characterization of higher-order quantum complexity.

We study the transport through a resonant level coupled to two leads with the latter being described by Wigner’s random matrices. By taking appropriate thermodynamic limit before taking the long time limit, we obtain the stationary current as a function of voltage bias. The I–V curve is similar to that of single impurity Anderson model. On the other hand, the current matrix and initial density matrix in our model look like random matrices in the eigenbasis of Hamiltonian. They satisfy the description of eigenstate thermalization hypothesis (ETH) and nonequilibrium steady state hypothesis (NESSH), respectively. A statistical formula of current has been derived based on ETH and NESSH (Wang 2017 J. Stat. Mech. 093105). We check this formula in our model and find it to predict the stationary current to a high precision. The shape of I–V curve is explained by the peak structure in the characteristic function of NESSH, which is reminiscent of the transmission coefficient.

We derive the eigenstate thermalization hypothesis (ETH) from a random matrix Hamiltonian by extending the model introduced by Deutsch (1991 Phys. Rev. A 43 2046). We approximate the coupling between a subsystem and a many-body environment by means of a random Gaussian matrix. We show that a common assumption in the analysis of quantum chaotic systems, namely the treatment of eigenstates as independent random vectors, leads to inconsistent results. However, a consistent approach to the ETH can be developed by introducing an interaction between random wave-functions that arises as a result of the orthonormality condition. This approach leads to a consistent form for off-diagonal matrix elements of observables. From there we obtain the scaling of time-averaged fluctuations of generic observables with system size for which we calculate an analytic form in terms of the inverse participation ratio. The analytic results are compared to exact diagonalizations of a quantum spin chain for different physical observables in multiple parameter regimes.

There is a dichotomy in the nonequilibrium dynamics of quantum many body systems. In the presence of integrability, expectation values of local operators equilibrate to values described by a generalized Gibbs ensemble, which retains extensive memory about the initial state of the system. On the other hand, in generic systems such expectation values relax to stationary values described by the thermal ensemble, fixed solely by the energy of the state. At the heart of understanding this dichotomy is the eigenstate thermalization hypothesis (ETH): individual eigenstates in nonintegrable systems are thermal, in the sense that expectation values agree with the thermal prediction at a temperature set by the energy of the eigenstate. In systems where ETH is violated, thermalization can be avoided. Thus establishing the range of validity of ETH is crucial in understanding whether a given quantum system thermalizes. Here we study a simple model with confinement, the quantum Ising chain with a longitudinal field, in which ETH is violated. Despite an absence of integrability, there exist rare (nonthermal) states that persist far into the spectrum. These arise as a direct consequence of confinement: pairs of particles are confined, forming new `meson' excitations whose energy can be extensive in the system size. We show that such states are nonthermal in both the continuum and in the low-energy spectrum of the corresponding lattice model. We highlight that the presence of such states within the spectrum has important consequences, with certain quenches leading to an absence of thermalization and local observables evolving anomalously.

In this thesis, we study systems of nonlinear partial differential equations that can be used to describe phase separation. This process can be modeled from both a macroscopic and microscopic point of view. Typically, the macroscopic model contains a local differential operator, which accounts for short-range interactions. In the microscopic model, however, this operator is then replaced by a nonlocal operator, which describes long-range interactions between the particles. This is done by convolution integrals weighted by suitable interaction kernels. In the first part, we analyze the asymptotic behavior of the nonlocal operator in the case, when the corresponding interaction kernel concentrates around the origin suitably. We prove that the nonlocal operator converges to a local differential operator as the parameter in the interaction kernel is sent to zero. More precisely, we even provide concrete rates of convergence depending on this parameter. The analysis is done in the case of sufficiently smooth bounded domains and a large class of interaction kernels. Indeed, the kernel can be either anisotropic or isotropic as well as of $W^{1,1}$-regularity or even more singular. The convergence is shown with respect to the $L^p$-norm, where $p\in[1,\infty)$. The proof is based on localization and perturbation arguments. In the case of regular kernels, we also prove convergence on the torus with respect to the $L^2$-norm. This proof is based on methods from Fourier analysis. In the second part, we focus on models, where two different phases are separated by a diffuse interface. We intend to apply the results from the first part to show that solutions to the nonlocal model converge to a solution to the local model as the parameter in the nonlocal operator is sent to zero. More precisely, we prove strong convergence for the sequence of solutions along with certain rates of convergence. This is done using an energy method. We consider the following models: - Cahn-Hilliard equation - Allen-Cahn equation - Cahn-Hilliard model for tumor growth - Cahn-Hilliard/Navier-Stokes system Finally, we also investigate the case, where both the parameter in the nonlocal operator and the thickness of the diffuse interface are sent to zero suitably. The analysis is done for the nonlocal Allen--Cahn equation on the torus. We first prove a well-posedness result and suitable error estimates for the solutions. Then, we combine the convergence result from the first part together with existing results for the sharp interface limit of the local Allen--Cahn equation.

Studies of Eigenstate Thermalization Hypothesis (ETH) in two-dimensional CFTs call for calculation of the expectation values of local operators in highly excited energy eigenstates. This can be done efficiently by representing zero modes of these operators in terms of the Virasoro algebra generators. In this paper we present a pedagogical introduction explaining how this calculation can be performed analytically or using computer algebra. We illustrate the computation of zero modes by a number of examples and list explicit expressions for all local operators from the vacuum family with the dimension of less or equal than eight. Finally, we derive an explicit expression for the quantum KdV generator Q7 in terms of the Virasoro algebra generators. The obtained results can be used for quantitative studies of ETH at finite value of central charge.

We say of an isolated macroscopic quantum system in a pure state psi that it is in macroscopic thermal equilibrium (MATE) if psi lies in or close to a suitable subspace mathcal {H}_textrm{eq} of Hilbert space. It is known that every initial state psi _0 will eventually reach and stay there most of the time (“thermalize”) if the Hamiltonian is non-degenerate and satisfies the appropriate version of the eigenstate thermalization hypothesis (ETH), i.e., that every eigenvector is in MATE. Tasaki recently proved the ETH for a certain perturbation H_theta ^textrm{fF} of the Hamiltonian H_0^textrm{fF} of Ngg 1 free fermions on a one-dimensional lattice. The perturbation is needed to remove the high degeneracies of H_0^textrm{fF}. Here, we first point out that also for degenerate Hamiltonians all psi _0 thermalize if the ETH holds, i.e., if every eigenbasis lies in MATE, and we prove that this is the case for H_0^textrm{fF}. Inspired by the fact that there is one eigenbasis of H_0^textrm{fF} for which MATE can be proved more easily than for the others, with smaller error bounds, and also in higher spatial dimensions, we show for any given H_0 that the existence of one eigenbasis in MATE implies quite generally that most eigenbases of H_0 lie in MATE. We also show that, as a consequence, after adding a small generic perturbation, H=H_0+lambda V with lambda ll 1, for most perturbations V the perturbed Hamiltonian H satisfies ETH and all states thermalize.

We investigate the eigenstate thermalization hypothesis (ETH) of highly excited descendant states in two-dimensional large central charge c conformal field theory. We use operator product expansion of twist operators to calculate the short interval expansions of entanglement entropy and relative entropy for an interval of length ℓ up to order ℓ12. Using these results to ensure ETH of a heavy state when compared with the canonical ensemble state up to various orders of c, we get the constraints on the expectation values of the first few quasiprimary operators in the vacuum conformal family at the corresponding order of c. Similarly, we also obtain the constraints from the expectation values of the first few Korteweg-de Vries charges. We check these constraints for some types of special descendant excited states. Among the descendant states we consider, we find that at most only the leading order ones of the ETH constraints can be satisfied for the descendant states that are slightly excited on top of a heavy primary state. Otherwise, the ETH constraints are violated for the descendant states that are heavily excited on top of a primary state.

We show that there are well separated families of quantum expanders with asymptotically the maximal cardinality allowed by a known upper bound. This has applications to the “growth” of certain operator spaces: It implies asymptotically sharp estimates for the growth of the multiplicity of M_N -spaces needed to represent (up to a constant C>1 ) the M_N -version of the n -dimensional operator Hilbert space OH_n as a direct sum of copies of M_N . We show that, when C is close to 1, this multiplicity grows as \exp{\beta n N^2} for some constant \beta>0 . The main idea is to relate quantum expanders with “smooth” points on the matricial analogue of the Euclidean unit sphere. This generalizes to operator spaces a classical geometric result on n -dimensional Hilbert space (corresponding to N=1). In an appendix, we give a quick proof of an inequality (related to Hastings's previous work) on random unitary matrices that is crucial for this paper.

The structure of this survey is the following: at first we repeat in Part I the first 20 years old proof of the strong Elliptic law for random matrices Then in Part II we prove the strong Elliptic law for random matrices Ξ n of the general form, i.e. when their diagonal entries have nonzero expectations, and when we do not require the existence of the probability densities of the entries of random matrices, but instead of Lyapunov condition we require the boundedness of the absolute moments of the entries of the order 4 + δ , δ &gt; 0. In this case the Elliptic law means that the support of the accompanying spectral density of eigenvalues looks like the picture of several galaxies made by telescope. If the distances between the centers of these galaxies are large enough we have several almost elliptical galaxies. We also present the sand clock law density for the sum of diagonal complex matrix and the symmetric random matrix belonging the domain of attraction of the Semicircular law. After this we prove in the Part III the weak Elliptic law, when all entries can have nonzero expectations. At the end of survey we give the proof of the Elliptic law for the Unitary random matrices from the so called Class N 12 of random Unitary matrices and give the over turned stool law density for the sum of diagonal complex matrix A n and where Ξ n belongs to the domain of attraction of the Circular law. We also establish in Part III for Elliptic Law at the first time the convergency rate and the Central Limit Theorem. These statements are based on the VICTORIA -transform of random matrix which is the abbreviation of the following words: Very Important Computational Transformation Of Random Independent Arrays. We follow the main strategy of the theory of limit theorems of the probability theory, i.e. we try to solve the problem of description of all limits of normalized spectral functions where λ k ( A n Ξ n B n + C n ) are eigenvalues of the matrix A n Ξ n B n + C n , A n , B n , and C n are nonrandom matrices, under general (as only possible) conditions on the entries of random matrices Ξ n , χ is the indicator function. We emphasize that the spectral theory of Hermitian random matrices is rather profound. For example, in 1975 in [3] V. Girko proved the general stochastic canonical equation for ACE (Asymptotically Constant Entries)-symmetric matrices: Assume that for any n , the random entries of a symmetric matrix are independent and they are asymptotically constant entries (ACE), i.e., for any ε &gt; 0, and τ &gt; 0 is an arbitrary constant, and that, for every 0 ≤ u ≤ 1 and 0 ≤ v ≤ 1, where the symbol denotes the weak convergence of distribution when n → ∞, is a nondecreasing function with bounded variation in z and continuous in u and v in the domain 0 ≤ u , v ≤ 1. Then, with probability one, for almost all x , where λ k (Ξ n × n ) are eigenvalues, F ( x ) is a distribution function whose Stieltjes transform satisfies the relation G α ( x , y , t ), as a function of x , is a distribution function satisfying the regularized stochastic canonical equation K 3 [3, 23] at the points x of continuity, is a random real functional whose Laplace transform of one-dimensional distribution is equal to The integrand <jats:inli

We consider a realistic nonequilibrium protocol, where a quantum system in thermal equilibrium is suddenly subjected to an external force. Due to this force, the system is driven out of equilibrium and the expectation values of certain observables acquire a dependence on time. Eventually, upon switching off the external force, the system unitarily evolves under its own Hamiltonian and, as a consequence, the expectation values of observables equilibrate towards specific constant long-time values. Summarizing our main results, we show that, in systems which violate the eigenstate thermalization hypothesis (ETH), this long-time value exhibits an intriguing dependence on the strength of the external force. Specifically, for weak external forces, i.e., within the linear response regime, we show that expectation values thermalize to their original equilibrium values, despite the ETH being violated. In contrast, for stronger perturbations beyond linear response, the quantum system relaxes to some nonthermal value which depends on the previous nonequilibrium protocol. While we present theoretical arguments which underpin these results, we also numerically demonstrate our findings by studying the real-time dynamics of two low-dimensional quantum spin models.

Integrable systems do not obey the strong eigenstate thermalization hypothesis (ETH), which has been proposed as a mechanism of thermalization in isolated quantum systems. It has been suggested that an integrable system reaches a steady state described by a generalized Gibbs ensemble (GGE) instead of thermal equilibrium. We prove that a generalized version of the strong ETH holds for noninteracting integrable systems with translation invariance. Our generalized ETH states that any pair of energy eigenstates with similar values of local conserved quantities looks similar with respect to local observables, such as local correlations. This result tells us that an integrable system relaxes to a GGE for any initial state that has subextensive fluctuations of macroscopic local conserved quantities. Contrary to the previous derivations of the GGE, it is not necessary to assume the cluster decomposition property for an initial state.

We study the nonequiliribium dynamics of atom-molecule Bose gases in a double-well potential. In this system, the internal atom-molecule tunneling has significant influence on the dynamics. We investigate the regularity of dynamics by studying the level statistics of the quantum system. We find that chaotic energy eigenstates arise from the competition between the interwell and the atom-molecule internal tunnelings. Furthermore, we show that the physical quantities relax to the microcanonical averages in the full-quantum dynamics when the system is chaotic. This thermalization is caused by the verification of the eigenstate thermalization hypothesis (ETH). We show numerically that the onset of ETH occurs simultaneously with that of chaos. In addition, we show that the energy eigenstates become exponentially localized states simultaneously with the onset of chaos.