In [MV] the difference qKZ equations were considered modulo a prime number $p$ and a family of polynomial solutions of the qKZ equations modulo $p$ was constructed by an elementary procedure as suitable $p$-approximations of the hypergeometric integrals. In this paper, we study in detail the first family of nontrivial example of the qKZ equations in characteristic $p$. We describe all solutions of these qKZ equations in characteristic $p$ by demonstrating that they all stem from the $p$-hypergeometric solutions. We also prove a Lagrangian property (called the orthogonality property) of the subbundle of the qKZ bundle spanned by the $p$-hypergeometric sections. This paper extends the results of [VV1] on the differential KZ equations to the difference qKZ equations.

This paper introduces a Heisenberg picture approach to Matrix Product States (MPS), offering a rigorous yet intuitive framework to explore their structure and classification. MPS efficiently represent ground states of quantum many-body systems, with infinite MPS (iMPS) capturing long-range correlations and thermodynamic behavior. We classify MPS into projective and non-projective types, distinguishing those with finite correlation structures from those requiring ergodic quantum channels to define a meaningful limit. Using the Markov-Dobrushin inequality, we establish conditions for infinite-volume states and introduce ergodic and mixing MPS. As an application, we analyze the depolarizing MPS, highlighting its lack of finite correlations and the need for an alternative ergodic description. This work deepens the mathematical foundations of MPS and iMPS, providing new insights into entanglement, phase transitions, and quantum dynamics.

We solve the Dry Ten Martini Problem for the unitary almost Mathieu operator with Diophantine frequencies in the non-critical regime.

New non-perturbative results on the eigenvalues of the spheroidal equation are presented. The results, found using an all orders WKB analysis, include a perturbative/non-perturbative (P/NP) relation as well as the first exponential correction to the perturbative series which is valid in certain regions of parameters. The quantum periods are also computed.

This paper solves ``The Dry Ten Martini Problem'' for $C^2$ cosine-type quasiperiodic Schr\"odinger operators with large coupling constants and Diophantine frequencies, a model originally introduced by Sinai in 1987 \cite{sinai}. This shows that the analyticity assumption on the potential is not essential for obtaining a dry Cantor spectrum and can be replaced by a certain geometric condition in the low regularity case. In addition, we prove the homogeneity of the spectrum and the absolute continuity of the integrated density of states (IDS).

This paper contains a review of the theoretical foundations of Clifford algebras, spinors and spinor bundles in the so-called co-frame formalism. A compact index-free notation is introduced, along with a series of identities useful for computations in supergravity theories.

This note examines the BV formulation of $N=1$, $D=4$ supergravity in the first-order Palatini--Cartan framework. Challenges in achieving an off-shell formulation are addressed by introducing corrections to the rank 2 BV action, offering in addition a solid foundation for the study of the theory on manifolds with boundary.

We rigorously construct a variety of orbits for certain delay differential equations, including the electrodynamic equations formulated by Wheeler and Feynman in 1949. These equations involve delays and advances that depend on the trajectory itself, making it unclear how to formulate them as evolution equations in a conventional phase space. Despite their fundamental significance in physics, their mathematical treatment remains limited. Our method applies broadly to various functional differential equations that have appeared in the literature, including advanced/delayed equations, neutral or state-dependent delay equations, and nested delay equations, under appropriate regularity assumptions. Rather than addressing the notoriously difficult problem of proving the existence of solutions for all the initial conditions in a set, we focus on the direct construction of a diverse collection of solutions. This approach is often sufficient to describe physical phenomena. For instance, in certain models, we establish the existence of families of solutions exhibiting symbolic dynamics. Our method is based on the assumption that the system is, in a weak sense, close to an ordinary differential equation (ODE) with "hyperbolic" solutions as defined in dynamical systems. We then derive functional equations to obtain space-time corrections. As a byproduct of the method, we obtain that the solutions constructed depend very smoothly on parameters of the model. Also, we show that many formal approximations currently used in physics are valid with explicit error terms. Several of the relations between different orbits of the ODE persist qualitatively in the full problem.

We study the non-modal stability of black hole spacetimes under linear perturbations. We show that large-amplitude growth can occur at finite time, despite asymptotic decay of linear perturbations. In the example presented, the physical mechanism is a transient form of superradiance, and is qualitatively similar to the transition to turbulence in Navier-Stokes shear flows. As part of the construction we provide a theorem for the positivity of QNM energies, and introduce a truncated-Hamiltonian approach to black hole pseudospectra which does not suffer from convergence issues.

Quantum counterparts of Schrodinger's classical bridge problem have been around for the better part of half a century. During that time, several quantum approaches to this multifaceted classical problem have been introduced. In the present work, we unify, extend, and interpret several such approaches through a classical large deviations perspective. To this end, we consider time-symmetric ensembles that are pre- and post-selected before and after a Markovian experiment is performed. The Schrodinger bridge problem is that of finding the most likely joint distribution of initial and final outcomes that is consistent with obtained endpoint results. The derived distribution provides quantum Markovian dynamics that bridge the observed endpoint states in the form of density matrices. The solution retains its classical structure in that density matrices can be expressed as the product of forward-evolving and backward-evolving matrices. In addition, the quantum Schrodinger bridge allows inference of the most likely distribution of outcomes of an intervening measurement with unknown results. This distribution may be written as a product of forward- and backward-evolving expressions, in close analogy to the classical setting, and in a time-symmetric way. The derived results are illustrated through a two-level amplitude damping example.

This review paper explores the Riccati-type pseudo-potential formulation applied to the quasi-integrable sine-Gordon, KdV, and NLS models. The proposed framework provides a unified methodology for analyzing quasi-integrability properties across various integrable systems, including deformations of the sine-Gordon, Bullough-Dodd, Toda, KdV, pKdV, NLS and SUSY sine-Gordon models. Key findings include the emergence of infinite towers of anomalous conservation laws within the Riccati-type approach and the identification of exact non-local conservation laws in the linear formulations of deformed models. As modified integrable models play a crucial role in diverse fields of non-linear physics-such as Bose-Einstein condensation, superconductivity, gravity models, optics, and soliton turbulence-these results may have far-reaching applications.

A well-developed classification program for 4d $\mathcal{N}=2$ super conformal field theories (SCFTs) leverages Seiberg-Witten geometry on the Coulomb branch of vacua; theories are arranged by increasing $\mathfrak{rank}$, the complex dimension of their Coulomb branch. An alternative organizational scheme focusses on the associated vertex operator algebra (VOA), which is more closely related to the Higgs branch. From the VOA perspective, a natural way to arrange theories is by their ``index of nilpotency'', the smallest integer $\mathfrak{n}$ such that $T^\mathfrak{n} = 0$ in the $C_2$ algebra, where $T$ is the VOA stress tensor. It follows from the Higgs branch reconstruction conjecture that $\mathfrak{n} < \infty$ for any 4d ${\cal N}=2$ SCFT. Extrapolating from several examples, we conjecture that $\mathfrak{n}$ is an RG monotone, $\mathfrak{n}_{\rm IR} \leq \mathfrak{n}_{\rm UV}$. What's more, we find in all cases that $\mathfrak{rank} \leq \mathfrak{n}-1$. Theory ordering by $\mathfrak{n}$ appears thus more refined than ordering by $\mathfrak{rank}$. For example, in the list of $\mathfrak{rank}=1$ theories, the Kodaira SCFTs and $SU(2)$ ${\cal N}=4$ SYM have $\mathfrak{n} =2$, while all others have $\mathfrak{n} >2$.

Hirota's discrete KdV (dKdV) equation is an integrable autonomous partial difference equation on $\mathbb{Z}^2$ that reduces to the Korteweg-de Vries (KdV) equation in a continuum limit. In this paper, we introduce a new non-autonomous version of the dKdV equation. Furthermore, we show that the new equation is integrable and admits discrete Painlev\'e transcendent solutions described by $q$-Painlev\'e equations of $A_J^{(1)}$-surface types ($J=3,4,5,6$).

In strictly axisymmetric configurations of tokamaks, field-line tracing reduces from a three-dimensional ODE system to a two-dimensional one, where Poincar\'e-Bendixson theorem applies and guarantees the nonexistence of chaos. The formulae of functional perturbation theory (FPT) mostly simplify to compact closed-form expressions to allow the computation to finish instantly, which could improve and accelerate the existing plasma control systems by detangling the plasma dynamics from the magnetic topology change. FPT can conveniently calculate how the key geometric objects of magnetic topology: 1. the divertor X-point(s) and the magnetic axis, 2. the last closed flux surface (LCFS) 3. flux surfaces change under perturbation. For example, when the divertor X-point shifts outwards, the LCFS there must expand accordingly, but not necessarily for other places of the LCFS, which could also contract, depending on the perturbation. FPT can not only facilitate adaptive control of plasma, but also enable utilizing as much as possible space in the vacuum vessel by weakening the plasma-wall interaction (PWI) via tuning the eigenvalues of $\mathcal{DP}^m$ of the divertor X-point(s), such that the field line connection lengths in the scrape-off layer (SOL) are long enough to achieve detachment. Increasing flux expansion $f_x$ is another option for detachment and can also be facilitated by FPT. Apart from the edge, FPT can also benefit the understanding of the plasma core. Since the magnetic axis O-point would also shift under perturbation and the shift is known by FPT, the O-point can be controlled without full knowledge of the plasma response, which shall not significantly change the tendency.

The initial-to-final-state inverse problem consists in determining a quantum Hamiltonian assuming the knowledge of the state of the system at some fixed time, for every initial state. This problem was formulated by Caro and Ruiz and motivated by the data-driven prediction problem in quantum mechanics. Caro and Ruiz analysed the question of uniqueness for Hamiltonians of the form $-\Delta + V$ with an electric potential $V = V(\mathrm{t}, \mathrm{x})$ that depends on the time and space variables. In this context, they proved that uniqueness holds in dimension $n \geq 2$ whenever the potentials are bounded and have super-exponential decay at infinity. Although their result does not seem to be optimal, one would expect at least some degree of exponential decay to be necessary for the potentials. However, in this paper, we show that by restricting the analysis to Hamiltonians with time-independent electric potentials, namely $V = V(\mathrm{x})$, uniqueness can be established for bounded integrable potentials exhibiting only super-linear decay at infinity, in any dimension $n \geq 2$. This surprising improvement is possible because, unlike Caro and Ruiz's approach, our argument avoids the use of complex geometrical optics (CGO). Instead, we rely on the construction of stationary states at different energies -- this is possible because the potential does not depend on time. These states will have an explicit leading term, given by a Herglotz wave, plus a correction term that will vanish as the energy grows. Besides the significant relaxation of decay assumptions on the potential, the avoidance of CGO solutions is important in its own right, since such solutions are not readily available in more complicated geometric settings.

In this paper, we address the Wigner distribution and the star exponential function for a time-dependent harmonic oscillator for which the mass and the frequency terms are considered explicitly depending on time. To such an end, we explore the connection between the star exponential, naturally emerging within the context of deformation quantization, and the propagators constructed through the path integral formalism. In particular, the Fourier-Dirichlet expansion of the star exponential implies a distinctive quantization of the Lewis-Riesenfeld invariant. Further, by introducing a judicious time variable, we recovered a time-dependent phase function associated with the Lewis-Riesenfeld construction of the standard Schr\"odinger picture. In particular, we applied our results to the cases of the Caldirola-Kanai and the time-dependent frequency harmonic oscillators, recovering relevant results previously reported in the literature.

We consider the Wigner minor process, i.e. the eigenvalues of an $N\times N$ Wigner matrix $H^{(N)}$ together with the eigenvalues of all its $n\times n$ minors, $H^{(n)}$, $n\le N$. The top eigenvalues of $H^{(N)}$ and those of its immediate minor $H^{(N-1)}$ are very strongly correlated, but this correlation becomes weaker for smaller minors $H^{(N-k)}$ as $k$ increases. For the GUE minor process the critical transition regime around $k\sim N^{2/3}$ was analyzed by Forrester and Nagao (J. Stat. Mech.: Theory and Experiment, 2011) providing an explicit formula for the nontrivial joint correlation function. We prove that this formula is universal, i.e. it holds for the Wigner minor process. Moreover, we give a complete analysis of the sub- and supercritical regimes both for eigenvalues and for the corresponding eigenvector overlaps, thus we prove the decorrelation transition in full generality.

We perform full integration of the stationary axisymmetric Einstein-Maxwell-dilaton-axion (EMDA) theory with and without potential using a recently proposed generalization of Carter's approach to spacetimes beyond type D, allowing the Killing tensor. Crucial to our construction is a new parametrization of the dilaton and axion fields based on the analyticity argument. The general solution in the ungauged case is asymptotically locally flat and contains two more parameters compared to EMDA black holes previously obtained using Harrison transformations. In the gauged case, the general solution is asymptotically AdS and includes flat and hyperbolic topological solutions, as well as generalization of the Kerr-Sen-AdS metric with three additional parameters. Our approach can be applied to more general four-dimensional ungauged and gauged supergravities.

In our recent papers, we studied semiclassical spectral problems for the Bochner-Schr\"odinger operator on a manifold of bounded geometry. We survey some results of these papers in the setting of the magnetic Schr\"odinger operator in the Euclidean space and describe some ideas of the proofs.

In their study of Levin-Wen models [Commun. Math. Phys. 313 (2012) 351-373], Kitaev and Kong proposed a weak Hopf algebra associated with a unitary fusion category $\mathcal{C}$ and a unitary left $\mathcal{C}$-module $\mathcal{M}$, and sketched a proof that its representation category is monoidally equivalent to the unitary $\mathcal{C}$-module functor category $\mathrm{Fun}^{\mathrm{u}}_{\mathcal{C}}(\mathcal{M},\mathcal{M})^\mathrm{rev}$. We give an independent proof of this result without the unitarity conditions. In particular, viewing $\mathcal{C}$ as a left $\mathcal{C} \boxtimes \mathcal{C}^{\mathrm{rev}}$-module, we obtain a quasi-triangular weak Hopf algebra whose representation category is braided equivalent to the Drinfeld center $\mathcal{Z}(\mathcal{C})$. In the appendix, we also compare this quasi-triangular weak Hopf algebra with the tube algebra $\mathrm{Tube}_{\mathcal{C}}$ of $\mathcal{C}$ when $\mathcal{C}$ is pivotal. These two algebras are Morita equivalent by the well-known equivalence $\mathrm{Rep}(\mathrm{Tube}_{\mathcal{C}})\cong\mathcal{Z}(\mathcal{C})$. However, we show that in general there is no weak Hopf algebra structure on $\mathrm{Tube}_{\mathcal{C}}$ such that the above equivalence is monoidal.

This paper discusses how to develop a high-order multiple-relaxation-time lattice Boltzmann (MRT-LB) model for the general d(>=1)-dimensional diagonal-anisotropic diffusion equation. Such an MRT-LB model considers the transformation matrix constructed in a natural way and the DdQ(2d^2+1) lattice structure. A key step in developing the high-order MRT-LB model is to determine the adjustable relaxation parameters and weight coefficients, which are used to eliminate the truncation errors at certain orders of the MRT-LB model, while ensuring the stability of the MRT-LB model. In this work, we first present a unified MRT-LB model for the diagonal-anisotropic diffusion equation. Then, through the direct Taylor expansion, we analyze the macroscopic modified equations of the MRT-LB model up to fourth-order, and further derive the fourth-order consistent conditions of the MRT-LB model. Additionally, we also construct the fourth-order initialization scheme for the present LB method. After that, the condition which guarantees that the MRT-LB model can satisfy the stability structure is explicitly given, and from a numerical perspective, once the stability structure is satisfied, the MRT-LB model must be L^2 stable. In combination with the fourth-order consistent and L^2 stability conditions, the relaxation parameters and weight coefficients of the MRT-LB model can be automatically given by a simple computer code. Finally, we perform numerical simulations of several benchmark problems, and find that the numerical results can achieve a fourth-order convergence rate, which is in agreement with our theoretical analysis. In particular, for the isotropic diffusion equation, we also make a comparison between the fourth-order MRT-LB models with the DdQ(2d^2+1) and DdQ(2d+1) lattice structures, and the numerical results show that the MRT-LB model with the DdQ(2d^2+1) lattice structure is more general.

Several important dynamical systems are in $\mathbb{R}^2$, defined by $(x',y')=(f(x,y),g(x,y))$. A question of fundamental importance is how such systems might behave quantum mechanically. In developing quantum theory, Dirac and others realized that classical Hamiltonian systems can be mapped to their quantum counterparts via canonical quantization. The resulting quantum dynamics is always physical, characterized by completely-positive and trace-preserving evolutions in the Schr\"odinger picture. However, whether non-Hamiltonian systems can be quantized systematically while respecting the same physical requirements has remained a long-standing problem. Here we resolve this question when $f(x,y)$ and $g(x,y)$ are arbitrary polynomials. By leveraging open-systems theory, we prove constructively that every polynomial system admits a physical generator of time evolution in the form of a Lindbladian. We refer to our method as cascade quantization, and demonstrate its power by analyzing several paradigmatic examples of nonlinear dynamics such as bifurcations, noise-activated spiking, and Li\'{e}nard systems. In effect, any classical system whose $f(x,y)$ and $g(x,y)$ are analytic functions can be quantized with arbitrary precision. Crucially, our method is exact. Being free from any approximations, cascade quantization dispenses with simplifying assumptions such as the weakly-nonlinear limit, or semiclassical dynamics in the quantized system -- both of which have been critical in facilitating quantization in the literature. We also highlight the advantages of cascade quantization over the existing proposals, by weighing it against examples from the variational paradigm using Lagrangians, as well as non-variational approaches.

We apply statistical analysis to search for processes responsible for turbulence in physical systems. In our previous studies, we have shown that solar wind turbulence in the inertial range of large magnetohydrodynamic scales exhibits Markov properties. We have recently extended this approach on much smaller kinetic scales. Here we are testing for the Markovian character of stochastic processes in a kinetic regime based on magnetic field and velocity fluctuations in the solar wind, measured onboard the Magnetospheric Multiscale (MMS) mission: behind the bow shock, inside the magnetosheath, and near the magnetopause. We have verified that the Chapman-Kolmogorov necessary conditions for Markov processes is satisfied for local transfer of energy between the magnetic and velocity fields also on kinetic scales. We have confirmed that for magnetic fluctuations, the first Kramers-Moyal coefficient is linear, while the second term is quadratic, corresponding to drift and diffusion processes in the resulting Fokker-Planck equation. It means that magnetic self-similar turbulence is described by generalized Ornstein-Uhlenbeck processes. We show that for the magnetic case, the Fokker-Planck equation leads to the probability density functions of the kappa distributions, which exhibit global universal scale invariance with a linear scaling and lack of intermittency. On the contrary, for velocity fluctuations, higher order Kramers-Moyal coefficients should be taken into account and hence scale invariance is not observed. However, the nonextensity parameter in Tsallis entropy provides a robust measure of the departure of the system from equilibrium. The obtained results are important for a better understanding of the physical mechanism governing turbulent systems in space and laboratory.

We investigate the evolution of anisotropies in Einstein-Gauss-Bonnet theory with a scalar field coupled to the Gauss-Bonnet term. Specifically, we examine the simplest scenario in which the scalar field lacks a kinetic term, and its kinetic contribution arises from an integration by parts of the Gauss-Bonnet scalar. We consider four- and five-dimensional anisotropic spacetimes, focusing on Bianchi I and extended Bianchi I geometries. Our study reveals that the asymptotic solutions correspond to locally symmetric spacetimes where at least two scale factors exhibit analogous behavior or, alternatively, to isotropic configurations where all scale factors evolve identically. Additionally, we discuss the effects of a cosmological constant, finding that the presence of the cosmological constant does not lead to an isotropic universe.

Majority dynamics is a process on a simple, undirected graph $G$ with an initial Red/Blue color for every vertex of $G$. Each day, each vertex updates its color following the majority among its neighbors, using its previous color for tie-breaking. The dynamics achieves \textit{unanimity} if every vertex has the same color after finitely many days, and such color is said to \textit{win}. When $G$ is a $G(n,p)$ random graph, L. Tran and Vu (2019) found a codition in terms of $p$ and the initial difference $2\Delta$ beteween the sizes of the Red and Blue camps, such that unanimity is achieved with probability arbitrarily close to 1. They showed that if $p\Delta^2 \gg1 $, $p\Delta \geq 100$, and $p\geq (1+\varepsilon) n^{-1}\log n$ for a positive constant $\varepsilon$, then unanimity occurs with probability $1 - o(1)$. If $p$ is not extremely small, namely $p > \log^{-1/16} n $, then Sah and Sawhney (2022) showed that the condition $p\Delta^2 \gg 1$ is sufficient. If $n^{-1}\log^2 n \ll p \ll n^{-1/2}\log^{1/4} n$, we show that $p^{3/2}\Delta \gg n^{-1/2}\log n$ is enough. Since this condition holds if $p\Delta \geq 100$ for $p$ in this range, this is an improvement of Tran's and Vu's result. For the closely related problem of finding the optimal condition for $p$ to achieve unanimity when the initial coloring is chosen uniformly at random among all possible Red/Blue assignments, our result implies a new lower bound $p \gg n^{-2/3}\log^{2/3} n$, which improves upon the previous bound of $n^{-3/5}\log n$ by Chakraborti, Kim, Lee and T. Tran (2021).

We study the following Liouville system defined on a flat torus \begin{equation} \left\{ \begin{array}{lr} -\Delta u_i=\sum_{j=1}^n a_{ij}\rho_j\Big(\frac{h_j e^{u_j}}{\int_\Omega h_j e^{u_j}}-1\Big),\nonumber u_j\in H_{per}^1(\Omega)\mbox{ for }i\in I=\{1,\cdots,n\}\nonumber, \end{array} \right. \end{equation} where $h_j\in C^3(\Omega)$, $h_j>0$, $\rho_j>0$ and $u=(u_1,..,u_n)$ is doubly periodic on $\partial\Omega$. The matrix $A=(a_{ij})_{n\times n}$ satisfies certain properties. One central problem about Liouville systems is whether multi-bubble solutions do exist. In this work we present a comprehensive construction of multi-bubble solutions in the most general setting.

It is a folklore belief that metastable wells in low-temperature statistical mechanics models exhibit high-temperature behavior. We prove a rigorous version of this phenomenon in the setting of the exponential random graph model (ERGM) through the lens of concentration of measure. To do this, we first present a new general result deriving concentration inequalities in a metastable well from the metastable mixing of a Markov chain with the appropriate stationary distribution, extending a result of Chatterjee [Cha05] which is suited for more traditional forms of global mixing. We then apply this result to the supercritical (low-temperature) ERGM which was recently proven to exhibit metastable mixing by Bresler, Nagaraj, and Nichani [BNN24], and obtain a novel concentration inequality for Lipschitz observables of the supercritical ERGM conditioned on a large metastable well, answering a question posed by [BNN24]. This extends a result of Ganguly and Nam [GN24] from the subcritical (high-temperature) regime to a metastable well in the supercritical regime, and we are also able to extend the applications of their concentration inequality to these metastable wells. Namely, we obtain an upper bound on the Wasserstein distance between the ERGM conditioned on a metastable well and an appropriate Erd\H{o}s-R\'enyi model, as well as derive a central limit theorem for the count of edges in certain small subcollections of possible edges. Finally, to supplement the mathematical content of the article, we also discuss the results of what appears to be the first simulation study of a metastable well in the supercritical ERGM.

We develop methods for systematic construction of superintegrable polynomials in matrix/eigenvalue models. Our consideration is based on a tight connection of superintegrable property of Gaussian Hermitian model and $W_{1 + \infty}$ algebra in Fock representation. Motivated by this example, we propose a set of assumptions that may allow one to recover superintegrable polynomials. The main two assumptions are box adding/removing rule (Pierri rule) and existence of Hamiltonian for superintegrable polynomials. We detail our method in case of the Gaussian Hermitian model, and then apply it to the cubic Kontsevich model.

The triad refers to embedding of two systems of polynomials, symmetric ones and those of the Baker-Akhiezer type into a power series of the Noumi-Shiraishi type. It provides an alternative definition of Macdonald theory and its extensions. The basic triad is associated with the vector representation of the Ding-Iohara-Miki (DIM) algebra. We discuss lifting this triad to two elliptic generalizations and further to the bi-elliptic triad. At the algebraic level, it corresponds to elliptic and bi-elliptic DIM algebras. This completes the list of polynomials associated with Seiberg-Witten theory with adjoint matter in various dimensions.

We study a random band matrix $H=(H_{xy})_{x,y}$ of dimension $N\times N$ with mean-zero complex Gaussian entries, where $x,y$ belong to the discrete torus $(\mathbb{Z}/\sqrt{N}\mathbb{Z})^{2}$. The variance profile $\mathbb{E}|H_{xy}|^{2}=S_{xy}$ vanishes when the distance between $x,y$ is larger than some band-width parameter $W$ depending on $N$. We show that if the band-width satisfies $W\geq N^{\mathfrak{c}}$ for some $\mathfrak{c}>0$, then in the large-$N$ limit, we have the following results. The first result is a local semicircle law in the bulk down to scales $N^{-1+\varepsilon}$. The second is delocalization of bulk eigenvectors. The third is a quantum unique ergodicity for bulk eigenvectors. The fourth is universality of local bulk eigenvalue statistics. The fifth is a quantum diffusion profile for the associated $T$ matrix. Our method is based on embedding $H$ inside a matrix Brownian motion $H_{t}$ as done in [Dubova-Yang '24] and [Yau-Yin '25] for band matrices on the one-dimensional torus. In this paper, the key additional ingredient in our analysis of $H_{t}$ is a new CLT-type estimate for polynomials in the entries of the resolvent of $H_{t}$.