Double-scaled bosonic and fermionic embedded ensembles, complex SYK, and the dual Hilbert space

Abstract: We derive the density of states along with the $2$- and $4$-point functions of embedded ensembles for both fermions and bosons in the double-scaled limit. It is shown the models are equivalent to the double-scaled Sachdev-Ye-Kitaev model, expanding the double-scaled universality class to include both fermionic and bosonic systems. The models can be solved by introducing the Wick product of non-commuting Gaussian random variables. We show that deriving the Wick product is sufficient for computing the density of states, and properties of the Wick product can be used to compute $n$-point functions directly in the energy basis. In this context, the Wick product is equivalent to normal ordering of $q$-oscillators, which leads to the duality between moments of double-scaled models and expectation values in the chord Hilbert space. By considering operator probes as a second set of oscillators, we extend the duality to compute $n$-point functions. Embedded ensembles are equivalent to complex SYK at fixed charge, and we show working directly with embedded ensembles streamlines the derivations.

• The paper establishes an analytic framework linking double-scaled bosonic and fermionic ensembles with DS-SYK models by deriving key spectral and correlation functions.
• It employs Wick product formalism and q-Hermite polynomials to obtain the density of states and explicit 2- and 4-point functions in the energy basis.
• The study demonstrates a deep duality between EGUE ensembles and q-oscillator Hilbert spaces, offering new insights into holographic duality in many-body chaos.

Double-Scaled Embedded Ensembles, Complex SYK, and Dual Hilbert Space: An Essay

Introduction and Context

The paper “Double-scaled bosonic and fermionic embedded ensembles, complex SYK, and the dual Hilbert space” (2604.14522) presents an analytic treatment of embedded random matrix ensembles (EGUE) for fermionic and bosonic systems in the double-scaled limit. The focus is on the derivation of spectral and correlation functions, establishing equivalence with double-scaled Sachdev-Ye-Kitaev (SYK) models, and elucidating dualities with Hilbert spaces of $q$-oscillators. This research systematically generalizes the universality class associated with SYK, previously restricted to fermionic systems, to encompass bosonic models—addressing long-standing questions about the holographic properties of bosonic many-body chaos.

Embedded Ensembles: Double-Scaled Limit and Chord Diagram Universality

The classical random matrix ensembles (GOE, GUE, GSE) effectively capture quantum chaotic behavior but suffer from pathological features in the context of many-body systems, specifically unphysical full-body interactions and the Wigner semicircle law for density of states. Embedded ensembles refine the description by introducing finite-body interactions, thereby ensuring sparsity and realistic physics. For a $p$-body EGUE, the Hilbert space dimension and statistical structure differ for fermions and bosons; the double-scaled limit (large $N$ and $m$, with fixed $m/N$ and $p^2/m$) is crucial for analytic tractability.

The authors rigorously demonstrate that moments of EGUE Hamiltonians, for both particle statistics, are represented as sums over chord diagrams with intersection weights determined by the parameters

$$q_{\pm} = \exp\left\{-\frac{p^2}{m}\frac{1}{1\pm\frac{m}{N}}\right\}$$

($-$ for fermions, $+$ for bosons). This is technically achieved through moment methods and unitary group decompositions, in contrast to previous probabilistic arguments used for SYK and spin glass models. The result confirms exact equivalence with DS-SYK models: the chord diagram combinatorics, and thus spectral and correlation structures, are universal within the double-scaled regime.

Wick Product Formalism and the Density of States

A central innovation of the paper is the introduction of a Wick product for non-commuting Gaussian random variables, shown to coincide with $q$-Hermite polynomials. The authors derive recursive relations for Wick polynomials, exploiting their orthogonality properties to obtain the density of states as a $p$0-normal distribution:

$p$1

where $p$2 and $p$3 encodes the intersection parameter. Limiting cases reproduce physically relevant distributions: for $p$4 (few-body/dilute limit), a Gaussian emerges; for $p$5 (many-body/full interaction), the semicircle law is recovered.

Crucially, the Wick product reformulation subsumes combinatorial aspects previously tied to chord diagrams, aligning analytic and diagrammatic approaches and providing streamlined methods for correlation function derivation.

Exact $p$6- and $p$7-Point Functions in Energy Basis

The paper advances techniques to compute $p$8-point functions directly in the energy basis, rooted in Wick product properties. The $p$9-point function (strength density), essential for transition probabilities and ETH analysis, is derived as a bivariate $N$0-normal distribution:

$N$1

with closed form for the underlying sum in terms of $N$2-Hermite polynomials and Pochhammer symbols.

The $N$3-point function, relevant for out-of-time-order correlators and quantum chaos diagnostics, is resolved via systematic contractions of Wick products, yielding expressions consistent with those obtained through chord diagram enumeration and transfer matrix methods. The analytic apparatus enables explicit evaluation and limiting behaviors, unifying previous diagrammatic and moment-based computations.

Dual Hilbert Space: $N$4-Oscillator Normal Ordering and Holographic Interpretation

By mapping moments and correlation functions to ground-state expectation values of transfer matrices in spaces of $N$5-oscillators, the paper elucidates deep duality relations. Normal ordering of transfer matrix powers reproduces $N$6-Hermite polynomials, establishing:

$N$7

This duality, central to holographic interpretations, extends to $N$8-point functions via a second set of oscillators for operator probes, matching chord Hilbert space approaches from prior SYK literature. The formalism underscores the bulk-boundary correspondence in AdS$N$9 settings, with transfer matrices representing gravitational Hamiltonians in the bulk.

Comparison to Complex SYK at Fixed Charge

The EGUE is explicitly shown to be equivalent to complex SYK (cSYK) in a fixed charge sector. The authors clarify differences in normalization and derivation between their moment methods and previous chord diagram analyses for cSYK. Working directly within fixed particle spaces avoids nontrivial contour integral projections and additional scaling factors, demonstrating analytic efficiency and conceptual clarity. This has practical implications for computations involving bosonic systems, where prior chord diagram analyses lack generalizability.

Implications and Future Directions

The formal equivalence of embedded ensembles (with both statistics) to DS-SYK models in the double-scaled limit enriches the universality class associated with quantum chaos and holographic duality. The paper’s analytic framework—combining Wick products, $m$0-Hermite polynomials, and dual Hilbert space mapping—provides scalable methods for spectral and dynamical function computation. The results directly impact the analysis of ETH, many-body chaos, and AdS$m$1 gravity.

Group-theoretical extensions (e.g., to $m$2 and $m$3 Racah-Wigner structures), perturbative genus expansions, and exploration of strength densities beyond the double-scaled regime represent compelling directions for further research. The techniques developed here are poised to facilitate broader investigations into the microscopic origins of quantum chaotic behavior, the spectrum of operator dynamics in many-body systems, and the robustness of holographic principles across particle statistics.

Conclusion

The study resolves longstanding technical questions regarding the universality and analytic solvability of embedded ensembles in the double-scaled limit, encompassing both fermionic and bosonic cases. The equivalence with DS-SYK and the rigorous dual Hilbert space mapping underpin advances in the analytic characterization of quantum chaotic spectra and correlation functions. The methodological innovations and clarifications provided establish new standards for computational tractability and theoretical generality in the study of many-body quantum chaos and holography.