---
title: Joint Splitting Probabilities
url: https://www.emergentmind.com/topics/joint-splitting-probabilities
type: topic
---

# Joint Splitting Probabilities

Joint splitting probabilities denote probability assignments that resolve an observed outcome into contributions indexed by additional labels that are not part of the final marginal event itself. In the cited literature, the phrase appears in several distinct but structurally related settings: complex weak-measurement decompositions for non-commuting observables, joint probability densities for how projection weight is distributed across channels or subspaces, transportation policies that split fixed marginals across pairs of indices, and first-passage objects that retain both exit identity and auxiliary hidden-state information [1303.0078][2109.01166][2508.07386].

## 1. Domain-specific meanings and common structure

Across the surveyed papers, the term is attached to several mathematically different objects rather than to a single standardized definition. What recurs is a decomposition of a marginal probability, density, or transport mass into contributions labeled by auxiliary variables, channels, observables, or internal states.

| Setting | Joint object | Characteristic relation |
|---|---|---|
| Weak-measurement quantum mechanics | \(p(a,b)\), \(p(m|a,b)\) | \(p(m)=\sum_{a,b}p(m|a,b)\,p(a,b)\) |
| Random orthonormal states | \(P_{NK}(t_1,\dots,t_R)\) | \(t_\xi=\sum_{j=1}^K |w_{j\xi}|^2\) |
| Transportation / joint tables | \(T_{ij}\) | \(\sum_j T_{ij}=p_i,\ \sum_i T_{ij}=q_j\) |
| Hidden-state first passage | \(\Pi(x_{\rm exit},y_{\rm exit}\mid x_0,y_0)\) | backward equation on \((x_0,y_0)\) |
| Active matter in confinement | \(\pi_\lambda(z,w)\) | \(p_v(w)\partial_z\pi_\lambda(z,w)=\rho(z,-w)\) |
| Rare-event control | \(q_B,\bar q_B\) | \(\lambda_i^*=2k_B T\,\partial_{r_i}\ln q_B\) |

Observable probabilities are then recovered by summing, integrating, or marginalizing over the auxiliary labels. In quantum weak-measurement formalisms those labels are incompatible observables; in random-matrix theory they are subspace projection weights; in first-passage theory they can be boundary identity, exit state, or number of jumps; in transportation problems they are pairwise allocations of prescribed marginals [1303.0078][1510.05333][2109.01166][2508.07386].

## 2. Complex joint splitting in quantum weak measurements

In the weak-measurement formulation of Hofmann, the starting point is a complex joint probability for two generally non-commuting observables with eigenstates \(\{|a)\}\) and \(\{|b)\}\). For a state described by density operator \(\hat p\), weak measurement of \(|a)(a|\) followed by a precise measurement in the \(|b)\)-basis defines
\[
p(a,b)=(b|a)(a|\hat p|b).
\]
This quantity is generally complex, is operationally reconstructed from weak values conditioned on post-selection of \(|b)\), and completely characterizes the quantum state through
\[
\hat p=\sum_{a,b} p(a,b)\,\frac{|a)(b|}{(b|a)}.
\]
The probability of a different measurement outcome \(|m)\) is then written in the splitting form
\[
p(m)\equiv (m|\hat p|m)=\sum_{a,b} p(m|a,b)\,p(a,b),
\qquad
p(m|a,b)=\frac{(b|m)(m|a)}{(b|a)}.
\]
Here \(p(a,b)\) carries the state dependence, while \(p(m|a,b)\) is state-independent and determined solely by overlaps between eigenstates [1303.0078].

This construction preserves the algebraic structure of classical probability splitting while abandoning positivity and, in general, reality. The objects \(p(a,b)\) and \(p(m|a,b)\) can be complex, and there is no classical sample space of simultaneous sharp values for \(a\) and \(b\). Nonetheless, the marginals reproduce standard Born probabilities, and the full state is recoverable from the complex joint distribution. In this sense, joint splitting probabilities become quasi-probabilistic weights attached to incompatible measurement labels rather than frequencies of jointly realized events [1303.0078].

The same paper emphasizes that the phase of the complex conditional probability encodes an action-like quantity,
\[
S(a,b,m)=\hbar\,\operatorname{Arg}\bigl(p(m|a,b)\bigr),
\]
and identifies this phase as the source of “logical tension” among incompatible properties. A probability is therefore “split” into complex contributions whose coherent recombination yields the real positive marginal \(p(m)\) [1303.0078].

## 3. Deterministic basis transformations and the exclusion of precise non-contextual quantum joints

A related formulation represents the density operator as a complex joint probability
\[
\rho(a,b)=\langle b|a\rangle\,\langle a|\hat\rho|b\rangle
\]
for any pair of observables with non-zero mutual overlap. Transformations to a different joint representation are expressed by complex conditional probabilities,
\[
\rho(c,b)=\sum_a p(c|a,b)\,\rho(a,b),
\qquad
p(c|a,b)=\frac{\langle b|c\rangle\langle c|a\rangle}{\langle b|a\rangle},
\]
and reversibility is encoded by
\[
\sum_c p(a'|c,b)\,p(c|a,b)=\delta_{a,a'}.
\]
In this framework, determinism is not a hidden-variable assignment of simultaneous values, but the state-independent reversibility of transformations between complete joint representations. Coarse-graining drives these complex conditional probabilities toward sharply peaked classical relations, so classical reality emerges as an approximation to the transformation laws of quantum determinism [1107.0536].

This interpretation has an immediate conceptual limit. “Excluding joint probabilities from quantum theory” studies two non-commuting observables in a two-dimensional Hilbert space and proves that there is no precise joint probability that applies for any quantum state and is consistent with quantum imprecise probabilities. The argument uses lower and upper probability operators \(\underline{\omega}(P,Q)\) and \(\overline{\omega}(P,Q)\), with a precise joint probability required to lie between the corresponding lower and upper probabilities for all states. In that sense, a universal, non-contextual, everywhere positive precise joint probability is excluded even for two non-commuting observables in dimension two. Context-dependent joint probabilities remain admissible, for example those associated with sequential measurement schemes, but they are still constrained by the imprecise bounds [1803.06722].

Taken together, these papers fix a persistent ambiguity. Complex joint splitting rules can be exact and state-complete as representations, but they are not classical joint probabilities. Conversely, demanding a precise positive non-contextual joint probability is too strong: the imprecise-probability analysis shows that such an object is not available in general [1107.0536][1803.06722].

## 4. Random orthonormal states, projection weights, and channel splitting

In random-matrix theory, joint splitting probabilities arise as joint probability densities for how projection weight is distributed across a fixed subspace by one or several random orthonormal states. Let \(w\in O(N)\) or \(U(N)\) be Haar distributed, and let
\[
t_\xi=\sum_{j=1}^K |w_{j\xi}|^2,\qquad \xi=1,\dots,R,
\]
be the projection probabilities of \(R\) orthonormal columns onto a \(K\)-dimensional subspace. The central object is the joint density
\[
P_{NK}(t_1,\dots,t_R).
\]
For a single vector, the distribution is Beta:
\[
t\sim \mathrm{Beta}\!\left(\frac K2,\frac{N-K}{2}\right)
\]
in the orthogonal case and
\[
t\sim \mathrm{Beta}(K,N-K)
\]
in the unitary case. The mean is \(K/N\) in both symmetry classes, while the variances differ by the usual orthogonal–unitary prefactor shift [1510.05333].

For several orthonormal states, the \(t_\xi\) are correlated because the underlying columns are orthogonal. The paper derives determinant-based integral representations for general \(R\), and explicit two-point formulas for specific \((N,K)\). In the unitary case with \(N=4\) and \(K=2\),
\[
P_{42}(t_1,t_2)=12
\begin{cases}
t_1 t_2, & t_1+t_2<1,\\
(1-t_1)(1-t_2), & t_1+t_2>1,
\end{cases}
\]
whereas the orthogonal case with the same \((N,K)\) involves a complete elliptic integral and develops a square-root singularity as \(t_1+t_2\to 1\) [1510.05333].

The same mathematics describes partial conductances in a two-terminal mesoscopic scattering problem. If \(t_\xi\) denotes the partial conductance associated with an incoming mode, then \(g=\sum_\xi t_\xi\) is the Landauer–Büttiker conductance, and the joint density \(P_{NK}(t_1,\dots,t_{N-K})\) becomes a joint splitting law for conductance across channels. In this usage, “joint splitting” refers to the statistical allocation of total norm or conductance among correlated channels rather than to first-passage competition or quasi-probabilistic observables [1510.05333].

## 5. Joint tables, transportation policies, and entropic geometry

A distinct classical use of the phrase appears in the study of joint probabilities with fixed marginals. Here the unknown is a joint matrix \(T=(T_{ij})\) satisfying
\[
\sum_j T_{ij}=p_i,\qquad \sum_i T_{ij}=q_j,\qquad T_{ij}\ge 0,
\]
and the entries are interpreted either as joint probabilities, transportation policies, or pixel intensities. In the transportation interpretation, \(T_{ij}\) is the fraction of supply at site \(i\) shipped to demand site \(j\); “joint splitting probabilities” are therefore feasible ways of splitting prescribed marginal masses across pairs \((i,j)\). Additional size constraints \(x_n\in[a_n,b_n]\) and cost constraints \(\sum_n W_n x_n=c\) are incorporated by a linear system \(Ax=y(c)\) [2109.01166].

The paper solves this constrained splitting problem by maximum entropy in the mean. The unknown \(x\) is represented as the expectation of a random vector \(\xi\) supported on the box \(\prod_n[a_n,b_n]\), and the relative entropy with respect to a product reference measure on the box endpoints is maximized under the mean constraint \(A\,E_P[\xi]=y(c)\). The resulting MEM solution is
\[
x_n^*=
\frac{a_n e^{-a_n(A^t\lambda^*)_n}+b_n e^{-b_n(A^t\lambda^*)_n}}
{e^{-a_n(A^t\lambda^*)_n}+e^{-b_n(A^t\lambda^*)_n}},
\]
and in the special case \(a_n=0\), \(b_n=1\),
\[
x_n^*=\frac{1}{1+e^{(A^t\lambda^*)_n}}.
\]
This gives an entropic family of admissible joint splittings compatible with marginals and linear expectations [2109.01166].

The same work endows the space of such splittings with a Hessian Riemannian geometry. In the “pixel space” of probabilities, the induced diagonal metric is
\[
G_n(x_n)=\frac{1}{(x_n-a_n)(b_n-x_n)},
\]
which diverges at the boundary of the feasible box. Geodesics and distances defined by this metric provide a geometry of alternative joint splittings, while a stepwise reduction of the cost constraint yields an interior-point-like approach to the minimum-cost transportation problem [2109.01166].

## 6. Hidden states, first passage, and active-matter confinement

For two-dimensional Markov processes \(\{X(t),Y(t)\}\), a direct joint analogue of splitting probabilities records both the exit boundary of \(X\) and the value of the hidden variable \(Y\) at exit. With \(X\) confined to an interval \(\mathbb{L}\) and \(Y\) taking values in \(\mathbb{Y}\), the joint splitting probability is
\[
\Pi(x_{\rm exit},y_{\rm exit}\mid x_0,y_0),
\]
the probability that \(X\), initialized at \(x_0\), first exits \(\mathbb{L}\) through \(x_{\rm exit}\) and that \(Y\), initialized at \(y_0\), equals \(y_{\rm exit}\) at exit time. It satisfies the stationary backward equation
\[
0=\mathcal{L}_B^{(X)}\Pi+\mathcal{L}_B^{(Y)}\Pi
\]
with absorbing or Robin boundary conditions, and conditional splitting probabilities are then defined by Bayes’ theorem,
\[
\Pi(y_{\rm exit}\mid x_{\rm exit};x_0,y_0)
=
\frac{\Pi(x_{\rm exit},y_{\rm exit}\mid x_0,y_0)}
{\sum_{y'}\Pi(x_{\rm exit},y'\mid x_0,y_0)}.
\]
The paper develops generic spectral formulas for decoupled Brownian \(X\) with an autonomous hidden state \(Y\), and explicit coupled calculations for run-and-tumble motion, diffusion in an intermittent piecewise-linear potential, and diffusion with stochastic resetting [2508.07386].

A complementary first-passage reformulation appears for run-and-tumble particles confined between two walls. There, the process is recast as a jump process with jump kernel \(G(\Delta x)\), and the splitting probability
\[
\pi(x)=\Pr\{\text{hit the right wall at }x=L\text{ before the left wall at }x=0\}
\]
satisfies the integral equation
\[
\pi(x)=\int_0^L dx'\,\pi(x')\,G(x'-x)+\int_L^\infty dx'\,G(x'-x).
\]
Its derivative gives the stationary bulk density,
\[
\rho(x)=\frac{d\pi}{dx},
\]
and the boundary values encode adsorption,
\[
\pi(0)=\frac{f_A}{2},\qquad \pi(L)=1-\frac{f_A}{2}.
\]
Iterating the integral equation yields a decomposition
\[
\pi(0)=\sum_{n=1}^\infty p_n,
\]
where \(p_n\) is the probability of hitting the opposite wall for the first time in exactly \(n\) jumps without intermediate return. In one dimension,
\[
f_A=\frac{2}{2+\lambda},
\qquad
\pi(x)=\frac12+(1-f_A)\left(\frac{x}{L}-\frac12\right),
\]
with \(\lambda=L/(\tau v_0)\) [2410.06393].

The phase-space formulation of confined active matter sharpens this relation. For RTP, ABP, and AOUP dynamics in dimensionless variables \((z,w)\), the velocity-resolved splitting probability \(\pi_\lambda(z,w)\) and the stationary phase-space density \(\rho(z,w)\) obey the exact correspondence
\[
p_v(w)\,\partial_z\pi_\lambda(z,w)=\rho(z,-w),
\qquad
p_v(w)\,\pi_\lambda(0,w)=f_0(-w)=f_\lambda(w).
\]
After integration over velocities,
\[
\rho(z)=\partial_z\pi_\lambda(z),
\qquad
\pi_\lambda(0)=f_w.
\]
Thus the stationary spatial density is literally the derivative of the splitting probability, while the wall-adsorbed fraction equals the wall splitting probability. The correspondence holds in arbitrary spatial dimension for Markovian active dynamics in phase space [2606.28709].

## 7. Memory effects, optimal control, and monitored quantum competition

For non-Markovian Gaussian first-passage problems, splitting probabilities depend on post-first-passage path statistics rather than only on geometry. In the two-target one-dimensional setting with absorbing points at \(0\) and \(L\), the splitting probabilities \(\pi_1(x_0)\) and \(\pi_2(x_0)\) satisfy
\[
\pi_1(x_0)+\pi_2(x_0)=1,
\]
but the key identity is
\[
x_0=\lim_{t\to\infty}\big[\pi_1\,\mu_1(t)+\pi_2\,\mu_2(t)\big],
\]
where \(\mu_i(t)\) are conditional mean trajectories after first passage to target \(i\) if motion is continued. For fractional Brownian motion with \(\psi(t)=\kappa t^{2H}\), the small-\(x_0\) behavior is
\[
\pi_2(x_0)\simeq A_H\left(\frac{x_0}{L}\right)^{1/H-1}.
\]
The paper derives self-consistent equations for \(\mu_1\) and \(\mu_2\) and validates the resulting splitting probabilities in viscoelastic-fluid experiments, where marked deviations from the Brownian law \(\pi_2=x_0/L\) are directly traced to out-of-equilibrium post-first-passage trajectories [2402.05005].

In stochastic control of rare reactive events, the splitting probability reappears as the committor \(\bar q_B(\mathbf r)\), the probability to hit metastable set \(B\) before \(A\). A time-dependent analogue \(q_B(\mathbf r,\tau)\) solves a generalized bridge problem, and in the bistable spectral approximation
\[
q_B(\mathbf r,\tau)=\bar q_B(\mathbf r)e^{-\mu_2\tau}+\bar p_B\big(1-e^{-\mu_2\tau}\big).
\]
The optimal controller that generates a bridge process statistically identical to the original reactive trajectory ensemble is
\[
\lambda_i^*(\mathbf r,\tau)=2k_B T\,\partial_{r_i}\ln q_B(\mathbf r,\tau).
\]
This makes the committor gradient not only a reaction-coordinate object but also a force field that guarantees reactivity while preserving the statistics of the transition path ensemble [2402.05414].

A further quantum extension appears in monitored continuous-time quantum walks with two targets. There the splitting probabilities
\[
P_\alpha(\psi_0)=\sum_{n=1}^\infty |\phi_n^{(\alpha)}|^2,\qquad \alpha=L,R,
\]
describe which target is detected first under repeated monitoring. By introducing symmetric and antisymmetric target superpositions \(\ket{d_\pm}=(\ket{x_L}\pm\ket{x_R})/\sqrt{2}\), the two-target problem maps onto a pair of single-target detection problems. For large systems, sampling times below
\[
\tau_c=\frac{2\pi}{\Delta E}
\]
produce a universal regime
\[
P_L(x_0)\approx P_R(x_0)\approx \frac12,
\]
independent of bulk initial condition and sampling time, whereas \(\tau>\tau_c\) yields a nonuniversal regime with pronounced peaks and dips as functions of \(\tau\) and \(x_0\). At resonant sampling times, dark states can make \(P_L+P_R<1\) [2601.16111].

A neighboring asymptotic theory for symmetric continuous jump processes shows that even the simplest boundary transmission problem retains microscopic information: the transmission probability from the boundary scales as
\[
\pi_{0,\underline{x}}(0)\sim
\frac{2^{\mu-1}}{\sqrt{\pi}}\,
\Gamma\!\left(\frac{1+\mu}{2}\right)
\left(\frac{a_\mu}{x}\right)^{\mu/2},
\]
which explains why the naive continuous-limit prediction \(\pi_{0,\underline{x}}(0)=0\) fails for genuine jump dynamics [2201.13179]. This suggests a broader pattern: whenever splitting observables are refined by hidden labels, incompatible observables, or microscopic jump structure, the resulting “joint” object carries information absent from the marginal splitting probability alone.

Source: https://www.emergentmind.com/topics/joint-splitting-probabilities