---
title: 'CNEMF-MDPs: Conditional Non-Exchangeable MDPs'
url: https://www.emergentmind.com/topics/conditional-non-exchangeable-mean-field-mdps-cnemf-mdps
type: topic
---

# CNEMF-MDPs: Conditional Non-Exchangeable MDPs

Conditional Non-Exchangeable Mean Field MDPs (CNEMF-MDPs) are infinite-horizon Markov Decision Processes with a continuum of heterogeneous agents interacting through a common noise, without assuming exchangeability. The framework is introduced in both a strong formulation and a label-state formulation, and the control problem is lifted to a standard MDP defined on the Wasserstein space of probability measures over the product of the label and state spaces. In this formulation, the label space represents agent heterogeneity, the state space is the individual state space, and a fixed distribution specifies the population of agent labels. Within this lifted state space, the value function is characterized as the unique fixed point of an appropriate Bellman operator. Subsequent work studies operators on constrained measure spaces of labeled conditional distributions and uses DeepONet-type branch-trunk architectures to approximate the resulting infinite-dimensional Bellman operators, thereby connecting CNEMF-MDPs to neural operator learning on $\mathcal M_\lambda$ [2603.00009] [2603.21683].

## 1. Strong and weak formulations

The canonical CNEMF-MDP setup begins with a label space $I=[0,1]$ endowed with a fixed (Lebesgue) law $\lambda$, an individual state space $X$ that is a compact Polish space with metric $d$, and an action space $A$ that is compact Polish with metric $d_A$. Heterogeneity is encoded by the label $u\in I$, and the population is not assumed to be exchangeable because labels enter the dynamics and reward explicitly. The stochastic environment contains idiosyncratic noises $\varepsilon_t^u$, $u\in I$, i.i.d. with law $\lambda_\varepsilon$, and common noise $\varepsilon_t^0$, i.i.d. with law $\lambda_{\varepsilon 0}$, independent of the idiosyncratic noises. Initial information is represented by $\Gamma^u$, independent across $u$ and from the noises, with $\xi^u=\xi_0(u,\Gamma^u)$ [2603.00009].

In the strong formulation, admissible controls are randomized open-loop families $\alpha=(\alpha_t^u)$ indexed by $u\in I$ and $t\in\mathbb N$,
$$
\alpha_t^u
=
a_t\bigl(u,\Gamma^u,(\varepsilon_s^u)_{s\le t},(\varepsilon_s^0)_{s\le t},Z_t^u\bigr),
$$
where $Z_t^u\sim \mathrm{Unif}[0,1]$ are independent randomizers. The state dynamics are
$$
X_{t+1}^u
=
F\bigl(u,X_t^u,\alpha_t^u,\mu_t^0,\varepsilon_{t+1}^u,\varepsilon_{t+1}^0\bigr),
$$
with conditional joint law
$$
\mu_t^0(dv,dx,da)=P^0_{(X_t^v,\alpha_t^v)}(dx,da)\,dv.
$$
The corresponding discounted reward is
$$
J_S(\xi^\cdot,\alpha)
=
\int_I
E\Bigl[\sum_{t\ge 0}\beta^t f(u,X_t^u,\alpha_t^u,\mu_t^0)\Bigr]\,du,
\qquad \beta\in(0,1).
$$

The weak, or label-state, formulation introduces a representative label $U\sim\lambda$ and a single process $(X_t,\alpha_t)$. One writes
$$
X_{t+1}
=
F\bigl(U,X_t,\alpha_t,\mu_t^0,\varepsilon_{t+1},\varepsilon_{t+1}^0\bigr),
\qquad
\mu_t^0(du,dx,da)=P^0_{(U,X_t,\alpha_t)}(du,dx,da),
$$
with admissible controls
$$
\alpha_t
=
a_t\bigl(U,\Gamma,(\varepsilon_s)_{s\le t},(\varepsilon_s^0)_{s\le t},Z_t\bigr),
$$
and payoff
$$
J_W(\xi,\alpha)=E\Bigl[\sum \beta^t f(U,X_t,\alpha_t,\mu_t^0)\Bigr].
$$
A key theorem shows that $V_S(\xi^\cdot)=V_W(\xi)$ whenever the laws of $\xi^u$ and $\xi|U=u$ coincide. This establishes equivalence between the strong and weak formulations and identifies the control problem with a law-invariant object indexed by the initial label-state distribution [2603.00009].

## 2. Lifted Wasserstein state space and Bellman equation

The defining structural move in CNEMF-MDPs is the lift from agent-level dynamics to a measure-valued MDP. In the general formulation, the lifted state is
$$
\mu_t=\mathrm{Law}^0(U,X_t)\in P_\lambda(I\times X),
$$
and the lifted control is
$$
\nu_t=\mathrm{Law}^0(U,X_t,\alpha_t)\in \tilde A:=P(I\times X\times A),
\qquad
\mathrm{pr}_{1,2}\nu_t=\mu_t.
$$
The transition kernel is
$$
\widehat F(\mu,\nu,e^0)
=
(u,x)\mapsto (u,F(u,x,a,\nu,\varepsilon,e^0))_\#(\nu\otimes\lambda_\varepsilon),
$$
so that $\mu_{t+1}=\widehat F(\mu_t,\nu_t,\varepsilon_{t+1}^0)$. The stage cost is
$$
\hat f(\mu,\nu)
=
\int_{I\times X\times A} f(u,x,a,\nu)\,\nu(du,dx,da),
$$
and the infinite-horizon value satisfies
$$
\widetilde V(\mu)
=
\sup_{(\nu_t)}
E\Bigl[\sum \beta^t \hat f(\mu_t,\nu_t)\Bigr].
$$
The Bellman operator on bounded measurable $V:P_\lambda(I\times X)\to\mathbb R$ is
$$
(TV)(\mu)
=
\sup_{\nu:\mathrm{pr}_{1,2}\nu=\mu}
\Bigl\{
\hat f(\mu,\nu)+\beta E\bigl[V(\widehat F(\mu,\nu,\varepsilon^0))\bigr]
\Bigr\}.
$$
Under Lipschitz-type assumptions on $F$ and $f$ in $(u,x,\mu)$, $T$ is a $\beta$-contraction in sup-norm, so the value function is the unique fixed point $V^*=TV^*$ [2603.00009].

A more specialized representation appears when the individual state space is $\mathbb R^d$ and the first marginal is prescribed. Let $I$ be a label space endowed with a reference measure $\lambda$, and define
$$
\mathcal M_\lambda
:=
\{\mu\in\mathcal P_2(I\times\mathbb R^d): \mathrm{pr}_1{}_\#\mu=\lambda\}.
$$
By disintegration,
$$
\mu_t(du,dx)=\lambda(du)\,\mu_t^u(dx),
$$
so $\mu_t^u\in\mathcal P_2(\mathbb R^d)$ is the conditional law of the $u$-agent. In this form, a stationary Markov feedback is a measurable map $a:I\times\mathbb R^d\to A$, and the transition operator $\mathcal T:\mathcal M_\lambda\times(I\times\mathbb R^d\to A)\to\mathcal M_\lambda$ is given by
$$
\mu_{t+1}(du,dy)
=
\lambda(du)\int_{x\in\mathbb R^d} K(u,x,a(u,x),dy)\,\mu_t^u(dx).
$$
Equivalently,
$$
\mu_{t+1}(du,dy)
=
\int_{I\times\mathbb R^d\times A}\delta_{(u,y)}\,K(u,x,a(u,x),dy)\,\mu_t(du,dx).
$$
The reward functional is
$$
R(\mu,a)
=
\int_{I\times\mathbb R^d} r(u,x,\mu,a(u,x))\,\mu(du,dx),
$$
and the Bellman equation takes the form
$$
V(\mu)
=
\sup_{a:I\times\mathbb R^d\to A}
\{R(\mu,a)+\gamma V(\mathcal T(\mu,a))\},
$$
or, in finite horizon,
$$
V_t(\mu)
=
\sup_a\{R(\mu,a)+V_{t+1}(\mathcal T(\mu,a))\}.
$$
This measure-valued Bellman equation makes explicit that the state variable is the full labeled conditional distribution, not a symmetric empirical law [2603.21683].

## 3. Disintegration, randomization, and measurable selection

The non-exchangeable character of CNEMF-MDPs is encoded by the label marginal constraint and by disintegration. Elements of $\mathcal M_\lambda$ correspond to families of labeled conditional distributions $\{\mu^u\}_{u\in I}$ through
$$
\mu(du,dx)=\lambda(du)\,\mu^u(dx).
$$
This is more restrictive than an arbitrary probability measure on $I\times\mathbb R^d$ and more informative than an exchangeable law on states alone, because the prescribed first marginal keeps the population of labels fixed while allowing label-dependent state laws and controls. A plausible implication is that the mean field in a CNEMF-MDP is best viewed as a conditional family indexed by labels rather than as a single unlabeled population distribution [2603.21683].

Randomization is structurally important. In the precursor CMKV-MDP analysis with common noise and open-loop controls, the lifted MDP on the space of probability measures requires relaxed controls of the form
$$
\nu_t(dx,da)=\mu_t(dx)\,\hat\alpha_t(x,da),
$$
or, equivalently, actions in $\mathcal P(\mathcal X\times A)$. Because the planner must “randomize” over the population to match an arbitrary $\hat a(x)\in\mathcal P(A)$ when the $\sigma$-field on the population is not already rich, one must allow randomized feedback policies
$$
\mathrm{mfa}:\mathcal P(\mathcal X)\times \mathcal X\times[0,1]\to A.
$$
That analysis also proves a measurable optimal coupling lemma: for compact Polish $(\mathcal X,d)$ there exists
$$
\zeta:\mathcal P(\mathcal X)\times\mathcal P(\mathcal X)\times\mathcal X\times[0,1]\to\mathcal X
$$
measurable such that if $X\sim\mu$ and $U\sim\mathrm{Unif}[0,1]$ independent of $X$, then $X'=\zeta(\mu,\mu',X,U)$ satisfies $X'\sim\mu'$ and
$$
E[d(X,X')]=W_1(\mu,\mu').
$$
This measurable-coupling construction is the technical device used to obtain jointly measurable selections and $\varepsilon$-optimal randomized feedback controls [1912.07883].

The CNEMF-MDP framework inherits this emphasis on randomization. In the 2026 formulation, admissible controls are randomized open-loop, and from an $\varepsilon$-optimal lifted randomized feedback $\nu^\varepsilon(\mu)$ one constructs strong and weak randomized feedback controls achieving the same performance. In the related non-exchangeable mean-field control problem with controlled interactions, relaxed admissible laws are the closure, in the Wasserstein topology, of strong feedback-law-valued pairs, and one has
$$
V_{\rm MFC}(\nu)=V_{\rm MFC}^R(\nu),
$$
together with existence of an optimal relaxed control under mild continuity and Lipschitz assumptions on $(b,L,g)$ [2603.00009] [2511.00288].

## 4. Quantitative propagation of chaos and finite-population control

A major contribution of the CNEMF-MDP theory is a quantitative comparison between the lifted infinite-agent problem and a finite $N$-agent MDP with common noise. In the finite system, the state is $x\in X^N$, the control is $a\in A^N$, and transitions are
$$
x'^i
=
F_N\bigl(i/N,x^i,a^i,\mu_{xa}^N,\varepsilon^{i/N},\varepsilon^0\bigr),
$$
with empirical law
$$
\mu_{xa}^N=\frac1N\sum \delta_{(i/N,x^i,a^i)}.
$$
The reward is
$$
R_N(x,a)=\frac1N\sum f_N(i/N,x^i,a^i,\mu_{xa}^N),
$$
and the Bellman operator on $W:X^N\to\mathbb R$ is
$$
(T_NW)(x)=\sup_a\{R_N(x,a)+\beta E[W(x')]\}.
$$
Under parallel Lipschitz assumptions, $T_N$ is a $\beta$-contraction with fixed point $V_N^*$ [2603.00009].

To compare finite and infinite systems, for any lifting $V:P_\lambda(I\times X)\to\mathbb R$ one defines
$$
W(x)=V(\mu_N^\lambda[x]),
$$
where $\mu_N^\lambda[x]$ is the labeled empirical measure sending $i/N$ to $x^i$. The finite-infinite comparison yields, for matching feedbacks,
$$
\bigl|\widetilde T^aV(x)-T_N\widetilde V(x)\bigr|
\le
C\Bigl(M_N^\gamma+\varepsilon_N^f+(\varepsilon_N^F)^\gamma+E\bigl[W_1(\mathrm{Law}^0(\cdots),\mu_N^\lambda[x])\bigr]^\gamma\Bigr),
$$
where
$$
M_N=E[W_1(\text{empirical},\mu)]\lesssim N^{-1/2}.
$$
By choosing a minimizing coupling, one deduces the explicit quantitative bound
$$
|V_N(x)-V^*(\mu_N^\lambda[x])|
\le
C\bigl(M_N^\gamma+\varepsilon_N^f+(\varepsilon_N^F)^\gamma\bigr).
$$
The same analysis yields a near-optimal policy construction: from an $\varepsilon$-optimal lifted feedback $\nu^\varepsilon(\mu)$ one builds an $N$-agent feedback
$$
\pi^N(x)\approx \arg\min W_1\bigl(\nu^\varepsilon(\mu_N^\lambda[x]),\mu_N^\lambda[x,\pi]\bigr),
$$
producing an
$$
O\bigl(\varepsilon+M_N^\gamma+\varepsilon_N^f+(\varepsilon_N^F)^\gamma\bigr)
$$
-optimal $N$-agent policy [2603.00009].

These estimates formalize the propagation of chaos in a non-exchangeable setting with common noise. They also justify the mean-field approximation operationally: $\varepsilon$-optimal policies for the limit CNEMF-MDP can be converted into near-optimal policies for finite heterogeneous systems. The paper identifies large-scale networked systems with heterogeneity, including power grids, financial networks, and graphon games, as application domains for this finite-to-infinite control transfer [2603.00009].

## 5. Neural operator approximation on $\mathcal M_\lambda$

The operator-learning perspective treats the Bellman map and related value operators as functions on the constrained measure space $\mathcal M_\lambda$. The central approximation theorem states, in paraphrased form, that if
$$
V:I\times\mathbb R^d\times\mathcal M_\lambda\to\mathbb R
$$
is continuous with mild growth and $\rho$ is any probability law over $\mathcal M_\lambda$, then for every $\varepsilon>0$ there exist $J\in\mathbb N$, cylindrical features $\phi_1,\dots,\phi_J\in C(I\times\mathbb R^d)$, a branch network $B:\mathbb R^J\to\mathbb R^r$, and trunk networks $T_1,\dots,T_r:I\times\mathbb R^d\to\mathbb R$ such that, with
$$
\Phi_J(\mu)=\bigl(\langle \phi_1,\mu\rangle,\dots,\langle \phi_J,\mu\rangle\bigr),
$$
one has
$$
\int_{\mathcal M_\lambda}
E_{(u,x)\sim\mu}
\Bigl[
\Bigl|
V(u,x,\mu)-\sum_{k=1}^r T_k(u,x)B_k(\Phi_J(\mu))
\Bigr|^2
\Bigr]
\,\rho(d\mu)
\le \varepsilon.
$$
The proof combines cylindrical approximations of probability measures with a DeepONet-type branch-trunk neural architecture. Its two stated ingredients are: first, cylindrical approximation on compact subsets of $\mathcal M_\lambda$ by functions of the form
$$
(u,x,\mu)\mapsto \sum_{k=1}^r f_k(\Phi_J(\mu))\,g_k(u,x),
$$
and second, finite-dimensional universal approximation of each $f_k$ and $g_k$ by small feed-forward nets [2603.21683].

In practical form, the learned operator is written
$$
\widehat V(\mu)(u,x)
=
\sum_{k=1}^r T_k(u,x;\theta_T)\,B_k(\Phi_J(\mu);\theta_B),
$$
or equivalently
$$
\widehat V(\mu)(u,x)=\sum_{k=1}^r \alpha_k(\mu)\,\phi_k(u,x),
$$
with $\alpha=B(\Phi_J(\mu))$ and $\phi=(T_1,\dots,T_r)$. The trunk net takes $(u,x)$ as input and outputs a scalar basis function, while the branch net takes the $J$-vector of moments $\Phi_J(\mu)$ and outputs the coefficients. This is the specific DeepONet branch-trunk architecture proposed for conditional mean-field operators on labeled distributions [2603.21683].

The sampling strategy for training is explicit. One fixes a non-atomic reference $\nu\in\mathcal P_2(\mathbb R^d)$, samples a random transport map $T(u,\cdot):\mathbb R^d\to\mathbb R^d$ in some parametrized family so that for each $u$, $T(u,\cdot)_\#\nu$ is the desired $\mu^u$, draws i.i.d. samples $(u_i,y_i)\sim\lambda\otimes\nu$, and sets $x_i=T(u_i,y_i)$. The empirical law of $(u_i,x_i)$ then lies in $\mathcal M_\lambda$. Repeating this procedure yields training measures $\{\mu^{(m)}\}$, and targets can be computed either by solving a reference discretization of the Bellman equation, or by Monte-Carlo simulation followed by regression. As a proof of concept, the paper applies this machinery to a CNEMF-MDP with controlled SDE / Euler-scheme dynamics and approximates the Bellman operator
$$
\mathcal B[V](\mu)(u,x)
:=
\max_{a\in A}
\Bigl\{
r(u,x,\mu,a)+\gamma E_Z[V(\mathcal T(\mu,a),u,x')]
\Bigr\}.
$$
Training uses $M\approx 10^4$ sampled laws $\mu^{(m)}$, and for each $\mu^{(m)}$, $N\approx 5\times 10^4$ particles to estimate temporal or dynamic-programming targets. The reported diagnostics are mean-squared rollout error versus the number $r$ of trunk-branch sensors, error as a function of the number $J$ of moments, and computational times on a single GPU. In all cases the experiments observe rapid decay of approximation error as $r,J$ increase, confirming the universal-approximation theory in practice [2603.21683].

## 6. Controlled interactions, related models, and interpretive issues

A related non-exchangeable mean-field control problem extends the CNEMF viewpoint by making the interaction structure itself controllable. In that formulation, each agent carries a private state $x\in\mathbb R^d$ and a label $u\in[0,1]$, and there are two kinds of controls: a local action $a\in A_{\rm reg}$ and a pairwise interaction action $r\in A_{\rm int}$. Heterogeneity is encoded by a structural kernel
$$
W:[0,1]^2\to[0,1],
$$
typically of graphon type. In the finite-$n$ system, this kernel is approximated by a step-kernel
$$
W^n(u,v)
=
\sum_{i,j=1}^n
\xi_{ij}^n\,
1_{((i-1)/n,i/n]}(u)\,
1_{((j-1)/n,j/n]}(v),
$$
which converges in cut-norm to $W$. The finite-agent dynamics are McKean-Vlasov-type SDEs involving empirical measures $M_t^{1,i,n}$ and $M_t^{2,i,n}$ that record outgoing and incoming controlled interactions, while the mean-field limit replaces sums by integrals and empirical laws by a flow $(\mu_t,\Lambda_t)$ [2511.00288].

This controlled-interaction model proves several properties that align with the broader CNEMF program. Under mild continuity and Lipschitz assumptions on $(b,L,g)$, the set of relaxed admissible laws is compact and convex in law, relaxed and strong formulations coincide in value, and an optimal relaxed control exists. If $\|f\circ W^n-f\circ W\|_{\Box}\to 0$ for all bounded Lipschitz $f$, then asymptotically optimal finite-$n$ controls admit subsequences whose empirical laws converge weakly to relaxed mean-field controls, and conversely any Lipschitz mean-field control can be lifted to finite systems with convergent empirical laws and payoffs. Together these imply
$$
\lim_{n\to\infty}V_n(\nu^n)=V_{\rm MFC}(\nu),
$$
as well as convergence of $\varepsilon_n$-optimal $n$-player controls to mean-field optimal controls. The corresponding Bellman-type equation is written on the space of measures and kernels:
$$
V(\mu,W)
=
\sup_{\Lambda\in\mathcal A(\mu)}
\Bigl\{\langle L(t,\cdot),\Lambda\rangle+V(\mu',W')\Bigr\}.
$$
Under standard Lipschitz and boundedness assumptions on $(b,L,g)$ and compactness of $A_{\rm reg}\times A_{\rm int}$, $V(\mu,W)$ is continuous in $\mu$ in Wasserstein distance and in $W$ in cut-norm [2511.00288].

Several interpretive points follow directly from these results. A common misconception is that heterogeneity can be absorbed into the state while preserving the standard exchangeable mean-field structure. The CNEMF-MDP framework is introduced precisely for systems “without assuming exchangeability,” and the fixed label distribution is part of the state constraint in $P_\lambda(I\times X)$ or $\mathcal M_\lambda$ [2603.00009]. Another misconception is that pure deterministic feedback always suffices. The common-noise mean-field MDP literature emphasizes the crucial role of relaxed controls and randomization hypotheses, and the non-exchangeable extensions retain that feature through lifted randomized feedbacks and relaxed law-valued controls [1912.07883]. Finally, the cited work identifies continuous-time extensions, learning in CNEMF-MDPs, and relaxed controls under weaker regularity as future directions, while the operator-learning results suggest a concrete computational route for high-dimensional Bellman problems on constrained Wasserstein spaces [2603.00009] [2603.21683].

Source: https://www.emergentmind.com/topics/conditional-non-exchangeable-mean-field-mdps-cnemf-mdps