---
title: Constrained Dyson Brownian Motion
url: https://www.emergentmind.com/topics/constrained-dyson-brownian-motion
type: topic
---

# Constrained Dyson Brownian Motion

Searching arXiv for the main paper and closely related work on constrained Dyson Brownian motion, interlacing minors, and geometric/ensemble-constrained variants.
Constrained Dyson Brownian motion denotes variants of Dyson’s eigenvalue diffusion in which the particle system is restricted by additional geometric, algebraic, or state-space constraints beyond the standard non-collision condition. A canonical and concrete instance is the joint evolution of the spectra of two consecutive principal minors of a matrix undergoing Dyson’s Brownian motion: the two spectra must satisfy an interlacing relation, and their coupled dynamics form a Markov diffusion on the corresponding Gelfand–Tsetlin cone [1007.0220]. In this setting, the constraint is neither an external hard wall nor an imposed reflection rule; it is inherited from matrix minor structure and is expressed simultaneously through the allowed state space, singular cross-covariances, and the stationary density [1007.0220].

## 1. Dyson Brownian motion and the meaning of constraint

Dyson’s original construction starts from matrix-valued Ornstein–Uhlenbeck dynamics on \(H_n^{(\beta)}\), the space of real symmetric, complex Hermitian, or quaternionic self-dual Hermitian matrices for \(\beta=1,2,4\). The independent real parameters of the matrix evolve as Ornstein–Uhlenbeck processes, and the invariant measure is the Gaussian ensemble \(e^{-\frac{\beta}{2}\operatorname{Tr} B^2}\,dB\) [1007.0220]. The ordered eigenvalues \(\lambda_1(t),\dots,\lambda_n(t)\) then satisfy the Dyson SDE
\[
d\lambda_i(t) = \left( -\lambda_i(t) + \sum_{j\neq i} \frac{1}{\lambda_i(t)-\lambda_j(t)} \right)dt + \sqrt{\frac{2}{\beta}\,db_{ii}(t), \qquad i=1,\dots,n,
\]
whose drift combines Ornstein–Uhlenbeck confinement with logarithmic repulsion [1007.0220].

In this basic form, Dyson Brownian motion already contains a constraint: eigenvalues are ordered and almost surely non-colliding. The more specific notion of constrained Dyson Brownian motion arises when further restrictions are imposed on the admissible configurations or on the geometry of the process. The most explicit example in the literature represented here is the interlacing constraint between the spectra of consecutive principal minors [1007.0220]. Other works treat different forms of constraint, including confinement by an external potential \(V\) in \(\beta\)-Dyson Brownian motion [1707.02700], killing on a subdomain in the colliding regime [2504.08467], confinement to the unit circle in unitary Dyson Brownian motion [2504.16551], and confinement to a Jordan curve in the plane [2603.05316]. These are related but distinct mechanisms.

For the single-level Dyson process, the generator on eigenvalues is
\[
\mathcal{A}_\lambda = \sum_{i=1}^n\left( \frac1\beta \frac{\partial^2}{\partial\lambda_i^2} + \left(-\lambda_i+\sum_{j\neq i}\frac{1}{\lambda_i-\lambda_j}\right) \frac{\partial}{\partial\lambda_i} \right),
\]
and the invariant density is the Gaussian \(\beta\)-ensemble
\[
\pi_\lambda(\lambda)\,d\lambda = C_{n,\beta}^{-1}e^{-\frac12\sum_i\lambda_i^2} \prod_{i<j}|\lambda_j-\lambda_i|^\beta \,d\lambda
\]
[1007.0220]. Constrained variants modify this baseline by changing the support, the coupling structure, or both.

## 2. Interlacing constraints from consecutive principal minors

Let \(B^{(N)}(t)\) denote the top-left \(N\times N\) principal minor of a matrix \(B(t)\) undergoing Dyson Ornstein–Uhlenbeck dynamics, and let \(B^{(N-1)}(t)\) be the consecutive smaller principal minor. If
\[
\lambda^{(N)}(t)=\lambda_1^{(N)}(t)<\cdots<\lambda_N^{(N)}(t), \qquad
\mu^{(N-1)}(t)=\mu_1^{(N-1)}(t)<\cdots<\mu_{N-1}^{(N-1)}(t),
\]
then classical interlacing of a Hermitian matrix with its principal minor gives
\[
\lambda_1 \le \mu_1 \le \lambda_2 \le \mu_2 \le \cdots \le \lambda_{n-1} \le \mu_{n-1} \le \lambda_n
\]
[1007.0220].

This defines the joint state space
\[
\{(\lambda,\mu)\in\mathbb{R}^n\times \mathbb{R}^{n-1} : \lambda_1<\cdots<\lambda_n,\ \mu_1<\cdots<\mu_{n-1},\ \lambda_i\le\mu_i\le\lambda_{i+1}\},
\]
which is a two-row Gelfand–Tsetlin cone [1007.0220]. The constraint is intrinsic: it is not added to the diffusion after the fact, but inherited from the algebraic relation between consecutive minors.

The central theorem of Adler, Nordenstam, and van Moerbeke is that the joint process
\[
(\lambda(t), \mu(t)) = \left(\spec B^{(n)}(t), \ \spec B^{(n-1)}(t)\right)
\]
is itself a Markov diffusion for \(\beta=1,2,4\), and the formulas extend to arbitrary \(\beta>0\) at the level of SDE and generator [1007.0220]. This gives a fully explicit realization of Dyson Brownian motion under interlacing constraints.

The mixed Vandermonde structure is central. Writing
\[
P_n(x)=\prod_{\alpha=1}^n(x-\lambda_\alpha),\qquad
P_{n-1}(x)=\prod_{\beta=1}^{n-1}(x-\mu_\beta),
\]
and
\[
\Delta_n(\lambda,\mu) = \prod_{i=1}^n\prod_{j=1}^{n-1}(\lambda_i-\mu_j)
= \prod_i P_{n-1}(\lambda_i)
= \prod_j P_n(\mu_j),
\]
the interlacing condition controls the sign and singular structure of these factors [1007.0220]. A plausible implication is that the geometry of the constrained state space is best understood through these mixed Vandermonde products rather than through a simple reflecting-boundary picture.

## 3. Two-level Markov diffusion: SDEs, generator, and boundary behavior

The lower level \(\mu\) evolves as an ordinary Dyson Ornstein–Uhlenbeck diffusion of size \(n-1\):
\[
d\mu_\beta(t) = \left( -\mu_\beta(t) + \sum_{\gamma\neq \beta}\frac{1}{\mu_\beta(t)-\mu_\gamma(t)} \right)dt + \sqrt{\frac{2}{\beta}\,db_\beta(t),\quad 1\le\beta\le n-1,
\]
whereas the upper level \(\lambda\) has the same deterministic Dyson drift at size \(n\), but its noise is coupled to \(\mu\) through the full bordered-matrix structure [1007.0220]. The coupling enters through covariances, not through an extra drift term:
\[
\frac{d\lambda_\alpha d\mu_\beta}{dt}
= -\frac{2}{\beta}\frac{1}{(\lambda_\alpha-\mu_\beta)^2}
\frac{P_{n-1}(\lambda_\alpha)P_n(\mu_\beta)}{P_n'(\lambda_\alpha)P_{n-1}'(\mu_\beta)}.
\]
Thus the deterministic drift in \(d\lambda_\alpha\) is pure Dyson at level \(n\), independent of \(\mu\), while the interlacing constraint is enforced stochastically through the quadratic variation structure [1007.0220].

At generator level, the coupled process decomposes as
\[
\mathcal{A}^\top = \mathcal{A}_\lambda^\top + \mathcal{A}_\mu^\top + \mathcal{A}_{\lambda\mu}^\top,
\]
with cross-term
\[
\mathcal{A}_{\lambda\mu}^\top := -\frac{2}{\beta} \sum_{i = 1}^n\sum_{j=1}^{n-1} \frac{\partial}{\partial\lambda_i} \frac{\partial}{\partial\mu_j} \left( \frac{1}{(\lambda_i-\mu_j)^2} \frac{P_{n-1}(\lambda_i)}{P_n'(\lambda_i)} \frac{P_n(\mu_j)}{P_{n-1}'(\mu_j)} \right)
\]
[1007.0220]. This singular second-order operator is the analytic expression of constrained Dyson motion for two consecutive minors.

The process does not hit the forbidden non-interlacing region. Near a boundary face, the gaps
\[
g_i^{(L)}(t)=\mu_i(t)-\lambda_i(t)\ge0, \qquad
g_i^{(R)}(t)=\lambda_{i+1}(t)-\mu_i(t)\ge0
\]
behave like square-root diffusions with positive drift away from zero:
\[
d(\mu_i-\lambda_i) =F_i(\lambda,\mu)\,dt +\sqrt{\mu_i-\lambda_i} \sum_{k\le\ell}\alpha_{k\ell}\,d\tilde b_{k\ell},
\]
\[
d(\lambda_{i+1}-\mu_i) =\hat F_i(\lambda,\mu)\,dt +\sqrt{\lambda_{i+1}-\mu_i} \sum_{k\le\ell}\hat\alpha_{k\ell}\,d\tilde b_{k\ell},
\]
with \(F_i,\hat F_i>0\) at the boundary [1007.0220]. This excludes a simple reflection interpretation. The paper explicitly notes that this two-level diffusion is not Warren’s interlaced Dyson process, whose boundary mechanism is different [1007.0220].

A common misconception is that interlacing-constrained Dyson systems must be implemented by reflecting Brownian motions at the interlacing faces. The two-level minor process shows otherwise: the constraint may instead be encoded by a singular covariance structure and Bessel-type repulsion from the boundary [1007.0220].

## 4. Transition law, invariant measure, and the role of \(\beta\)

For one level, the transition density involves the Harish-Chandra/Itzykson–Zuber integral \(F_n^{(\beta)}\) and the usual Vandermonde factor \(|\Delta_n(\lambda)|^\beta\) [1007.0220]. For two consecutive minors, the transition density is more elaborate and contains three distinct Vandermonde-type contributions:
\[
|\Delta_n(\lambda)\,\Delta_{n-1}(\mu)|\,|\Delta_n(\lambda,\mu)|^{\frac{\beta}{2}-1},
\]
together with a generalized HCIZ-type integral \(G_{n-1}^{(\beta)}\) and integrals over sphere variables that arise from angular degrees of freedom of the border vectors [1007.0220].

The invariant density is
\[
I(\lambda,\mu)\,d\lambda\,d\mu
= \hat Z_{n,\beta}^{-1}\,\bigl(\operatorname{vol}(S^{\beta-1})\bigr)^{2(n-1)}
\exp\!\left(-\frac{\beta}{2}\sum_{i=1}^n\lambda_i^2\right)
|\Delta_n(\lambda)\,\Delta_{n-1}(\mu)|\,|\Delta_n(\lambda,\mu)|^{\frac{\beta}{2}-1}\, d\lambda\,d\mu
\]
[1007.0220]. The support is still the interlacing domain, so the constraint survives even when the mixed factor has exponent zero.

The case \(\beta=2\) is special because
\[
|\Delta_n(\lambda,\mu)|^{\frac{\beta}{2}-1}=|\Delta_n(\lambda,\mu)|^0=1,
\]
so the mixed interaction disappears from the stationary density, although not from the support or from the dynamics [1007.0220]. This often invites overinterpretation. The disappearance of the cross Vandermonde factor at equilibrium does not mean the two levels decouple; the coupling remains present in the admissible state space and in the stochastic cross-variation.

The construction is explicit for \(\beta=1,2,4\), where a matrix model exists, and Corollary 3 extends the SDE and generator to arbitrary \(\beta>0\) by taking the same covariance structure and generator as definition [1007.0220]. In that sense, constrained two-level Dyson diffusion has a genuine \(\beta\)-version beyond the classical matrix cases.

A related but different notion of constrained Dyson motion appears in \(\beta\)-Dyson Brownian motion with confining potential \(V\):
\[
d\lambda_{i,t}
= \sqrt{\frac{2}{\beta}\,dZ_{i,t}
-\biggl(
\frac{V'(\lambda_{i,t})}{2}
- \sum_{j\neq i} \frac{1}{\lambda_{i,t} - \lambda_{j,t}}
\biggr)\,dt,
\]
whose invariant law is the general \(\beta\)-ensemble
\[
\rho^{\beta}_N(\lambda)
= \frac{1}{Z^{\beta}_N}
\exp\!\Biggl\{ -\frac{\beta}{2} \biggl(
-\sum_{1 \leq i \neq j \leq N} \log|\lambda_i - \lambda_j|
+ \sum_{i=1}^N V(\lambda_i)
\biggr) \Biggr\}
\]
[1707.02700]. Here the “constraint” is confinement by \(V\), rather than interlacing between levels.

## 5. Markovianity breaks at three consecutive minors

A sharp structural feature of the consecutive-minor construction is that two levels are Markovian but three consecutive levels are not, at least for \(\beta=2,4\) [1007.0220]. If
\[
(\lambda,\mu,\nu) = (\spec B^{(n)}, \spec B^{(n-1)}, \spec B^{(n-2)}),
\]
then for generic initial conditions the triple process is not Markovian [1007.0220].

The mechanism is generator-theoretic. If a projection \(Y_t=\varphi(X_t)\) of a diffusion \(X_t\) is Markovian, then the full generator must map functions of \(Y\) into functions of \(Y\). For three consecutive minors, Adler, Nordenstam, and van Moerbeke exhibit matrix observables whose evolution under the full Dyson generator depends on information not recoverable from \((\lambda,\mu,\nu)\) alone [1007.0220]. In the \(n=3,\beta=2\) case, the missing information can be seen explicitly through phases of off-diagonal entries that affect determinants but are not encoded in the three spectra [1007.0220].

This establishes a precise boundary for interlacing-constrained Markovian Dyson dynamics in the principal-minor setting: two consecutive rows of the Gelfand–Tsetlin pattern suffice, three do not [1007.0220]. A plausible implication is that the angular variables eliminated successfully in the two-level reduction cease to be eliminable at the three-level stage.

This non-Markovianity is important conceptually because it prevents an overly naive extrapolation from two-level formulas to full interlacing arrays. The two-level process is explicit and Markovian; higher-level constrained dynamics require different constructions.

## 6. Related constrained variants and broader geometric formulations

The principal-minor model is one important realization of constrained Dyson Brownian motion, but the broader literature represented here shows that “constraint” appears in several mathematically distinct forms.

In tridiagonal models for \(\beta\)-Dyson Brownian motion with confining potential \(V\), the eigenvalues follow the constrained log-gas SDE while the matrix is represented by a Jacobi diffusion whose entries evolve according to explicit SDEs derived through Lanczos tridiagonalization [1707.02700]. For the quadratic potential \(V(x)=2x^2\), the eigenvalue distribution is “constrained to lie roughly in \((-\sqrt{n},\sqrt{n})\)” and the large-\(n\) analysis describes the dynamics of bounded principal submatrices after centering by \(D\sqrt n\) [1707.02700]. This is a matrix-level realization of confinement-constrained Dyson flow rather than an interlacing constraint.

A different constrained regime is Dyson Brownian motion with collisions and killing on a domain \(\mathscr U\subset\overline{\mathscr O}\), where
\[
\mathscr O := \{x\in\mathbb{R}^N : x^1 < x^2 < \dots < x^N\},
\]
and the process is defined via a multivalued SDE with interaction potential
\[
V_I(x)
=
-\gamma \sum_{1\le i<j\le N}\ln(x^j-x^i)
\]
on \(\mathscr O\) [2504.08467]. For \(0<\gamma<1/2\), collisions occur almost surely in finite time, yet the singular drift remains integrable through collisions and the killed process in \(\mathscr U\) admits a unique quasi-stationary distribution under suitable assumptions [2504.08467]. This is a constraint by killing rather than by interlacing.

Geometric constraints also arise on compact or curved configuration spaces. For unitary Dyson Brownian motion, the particles are eigenangles on the circle and satisfy
\[
d\lambda_t^i = \frac{1}{N}\sum_{j\ne i}\cot\left(\frac{\lambda_t^i-\lambda_t^j}{2}\right)\,dt + \frac{2}{\sqrt N}\,dB_t^i
\]
in the unitary scaling considered by Bertucci and Pesce [2504.16551]. Here the constraint is that eigenvalues live on \(\mathbb T\), and the macroscopic limit is studied through a viscosity equation for the primitive of the empirical measure [2504.16551]. Likewise, Dyson Brownian motion has been constructed on a rectifiable Jordan curve \(\Gamma\), where the Gibbs density is
\[
\varrho_{\beta,N}(\boldsymbol{z}) = \frac{1}{Z_{\beta,N}(\gamma)} \prod_{1\le i \ne j \le N} |z_i - z_j|^{\beta/2}
\]
and the arc-length parametrization induces an SDE with drift
\[
b_i(x) = \frac{\beta}{2}\sum_{j\neq i} \operatorname{Re} \left\{
\frac{\gamma'(x_i)}{\gamma(x_i) - \gamma(x_j)} \right\}
\]
[2603.05316]. This is a geometric confinement constraint on a nontrivial one-dimensional manifold.

A further geometric perspective appears in the extrinsic construction of Dyson Brownian motion via projection of Brownian motion onto tangent and normal directions of isospectral group orbits in \(H(n)\) [2210.11347]. In that formulation, \(\beta\) controls anisotropy between tangential and normal noise, and in the limit \(\beta=\infty\) the eigenvalues evolve deterministically by Coulomb repulsion while the isospectral orbits move by minus one half times mean curvature [2210.11347]. This suggests that constrained Dyson motion can often be understood as stochastic dynamics near a foliation by group orbits.

## 7. Significance and conceptual interpretation

The two-consecutive-minor diffusion of Adler, Nordenstam, and van Moerbeke provides a particularly sharp and explicit example of constrained Dyson Brownian motion because the entire constrained structure is available in closed form: SDE, generator, transition density, and invariant measure [1007.0220]. The process may be viewed as two Dyson diffusions whose admissible configurations lie in a two-row Gelfand–Tsetlin cone and whose fluctuations are coupled by a singular cross-operator [1007.0220].

Three complementary descriptions of the constraint are equivalent in this setting. First, there is state-space restriction: only interlacing configurations are allowed. Second, there is stochastic enforcement: cross-covariances become singular as interlacing faces are approached and induce Bessel-type repulsion from the boundary. Third, there is stationary encoding: the invariant measure carries levelwise Vandermonde factors together with a mixed power of \(|\Delta_n(\lambda,\mu)|\) [1007.0220].

This suggests a general editorial shorthand, *Editor’s term*: “constraint by inherited geometry,” for those Dyson systems where the restriction is induced by an underlying matrix, group-orbit, or manifold structure rather than by externally imposed boundary conditions. In the principal-minor model, the inherited geometry is the interlacing of spectra of consecutive minors [1007.0220]. In Jordan-curve and unitary models, it is the ambient curve or circle [2603.05316, 2504.16551]. In tridiagonal and \(V\)-constrained models, it is the confining potential and the associated \(\beta\)-ensemble [1707.02700].

The sharp failure of Markovianity at three consecutive minors is equally significant. It shows that not every natural constrained spectral projection of matrix Dyson dynamics closes as a diffusion [1007.0220]. This delineates the scope of explicit constrained Dyson theories and clarifies why two-level interlacing diffusions occupy a special place: they are rich enough to display nontrivial constrained geometry, yet simple enough to remain Markovian.

In summary, constrained Dyson Brownian motion is not a single model but a family of Dyson-type diffusions modified by structural restrictions. Among them, the spectra of two consecutive principal minors under Dyson Ornstein–Uhlenbeck dynamics furnish a paradigmatic finite-\(n\) example: a Markov diffusion on an interlacing cone, with pure Dyson drift on each level, singular noise coupling between levels, explicit transition laws, and a stationary density reflecting both Vandermonde repulsion and interlacing geometry [1007.0220].

Source: https://www.emergentmind.com/topics/constrained-dyson-brownian-motion