Coupled Differential Riccati Equations

• Coupled Differential Riccati Equations (CDREs) are systems of Riccati differential equations interconnected through modes, players, or mean-field interactions.
• They are applied in linear-quadratic differential games, regime-switching models, and PDE control to analyze stability, optimality, and feedback gains.
• Methodologies include invariant subspace and geometric analyses, operator techniques, and residue calculus to reduce and solve the coupled dynamics.

Searching arXiv for recent and foundational papers on coupled differential Riccati equations and closely related Riccati frameworks. Coupled Differential Riccati Equations (CDREs) denote a family of Riccati-type differential systems in which the Riccati unknown is coupled through modes, players, mean fields, or subsystem interactions rather than evolving as an isolated equation. In the narrow control-theoretic sense, CDREs are multiple matrix Riccati differential equations coupled across indices, as in Markov regime-switching stochastic differential games; in broader usage, the literature also includes globally coupled scalar complex Riccati equations, single operator-valued differential Riccati equations attached to coupled PDE systems on product spaces, and state-dependent generalized differential Riccati equations derived from invariant Lagrangian subbundles (Li et al., 11 Sep 2025, Pazó et al., 5 Mar 2025, Acquistapace et al., 2020, Schaft, 2015).

1. Terminological scope and structural variants

The literature does not use the expression “CDRE” in a single uniform sense. In standard control usage, it often refers to multiple matrix Riccati equations $P_i(t)$ coupled across indices $i$, for example in Markov jump systems, zero-sum games, or multi-player LQ problems. Other works use closely related language for coupled scalar Riccati dynamics, operator Riccati equations for coupled infinite-dimensional plants, or generalized differential Riccati equations with side conditions (Schaft, 2015).

Setting	Riccati unknown	Source of coupling
Regime-switching SLQ games	$P_T(t,i)$, $\Sigma(t,i)$	Markov-generator terms $\sum_j \pi_{ij}P(t,j)$ or $\sum_k \lambda_{ik}\Sigma(s,k)$
Globally coupled complex Riccati ensembles	$z_j$, reduced $Z,Q,\xi$	Coefficients depend self-consistently on the mean field $Z$
Coupled PDE control systems	$P(t)\in\mathcal L(Y)$	Coupling lies in the underlying product-space dynamics
Geometric differential Riccati framework	$i$0	Tangent/cotangent coupling and invariance with respect to $i$1
Harmonic-coupled Riccati equations	$i$2	Discrete-time harmonic means across a network

A useful distinction is between coupling of the Riccati unknowns and coupling of the controlled plant. The operator-theoretic PDE paper studies a single DRE on a product Hilbert space, so the underlying hyperbolic/parabolic dynamics are coupled even though the Riccati side is not presented as several matrix equations (Acquistapace et al., 2020). By contrast, the regime-switching game papers formulate one Riccati matrix per mode and couple those matrices explicitly through the generator of the Markov chain (Li et al., 11 Sep 2025, Wu et al., 2024).

This suggests that “CDRE” is best treated as a family-resemblance term. The common feature is not a unique canonical formula, but the presence of Riccati differential structure together with a nontrivial coupling mechanism.

2. Matrix CDREs in differential games

A standard finite-dimensional interpretation appears in linear-quadratic differential games. For a two-player closed-loop Nash game over a finite horizon, the time-dependent value matrices $i$3 and feedback gains $i$4 satisfy coupled differential Riccati equations of the form

$i$5

together with the time-varying stationarity relation obtained from the algebraic block condition for $i$6. The paper "Numerical Method for a Class of Algebraic Riccati Equations" studies only the infinite-horizon stationary counterpart, but it explicitly interprets its coupled algebraic Riccati equations as the time-invariant steady state of such a finite-horizon coupled Riccati flow (Li et al., 2020).

The most explicit stochastic game formulation in the supplied literature is the regime-switching zero-sum SLQ setting. For each mode $i$7, the finite-horizon Riccati matrix $i$8 satisfies

$i$9

where the Markov jump contribution $P_T(t,i)$0 is part of $P_T(t,i)$1. This term is the source of coupling across modes: when Itô’s formula is applied to $P_T(t,i)$2, the jumps of $P_T(t,i)$3 force the equation in mode $P_T(t,i)$4 to depend on all $P_T(t,i)$5 (Li et al., 11 Sep 2025).

Under the paper’s assumptions (A1) uniform convexity-concavity and (A2) $P_T(t,i)$6-stability of $P_T(t,i)$7, the finite-horizon CDREs admit a unique strongly regular solution, with

$P_T(t,i)$8

The saddle feedback is then

$P_T(t,i)$9

the finite-horizon value function is quadratic,

$\Sigma(t,i)$0

and the infinite-horizon limit is the coupled algebraic Riccati system

$\Sigma(t,i)$1

The convergence of $\Sigma(t,i)$2 to $\Sigma(t,i)$3 is exponential: $\Sigma(t,i)$4 and the same estimate holds for the gains $\Sigma(t,i)$5. These estimates are the analytic input for the turnpike theorem in that paper (Li et al., 11 Sep 2025).

3. Stackelberg and hierarchical CDRE systems

A distinct but related hierarchy arises in the zero-sum stochastic linear-quadratic Stackelberg differential game with Markovian regime switching. Here the follower first solves a regime-coupled Riccati equation

$\Sigma(t,i)$6

which yields the rational reaction

$\Sigma(t,i)$7

The leader’s problem is then reduced to a backward SLQ problem whose coefficients already depend on the follower Riccati solution $\Sigma(t,i)$8 (Wu et al., 2024).

After the transformation

$\Sigma(t,i)$9

and the Riccati ansatz

$\sum_j \pi_{ij}P(t,j)$0

the leader’s coupled differential Riccati equation becomes

$\sum_j \pi_{ij}P(t,j)$1

with

$\sum_j \pi_{ij}P(t,j)$2

The coupling is therefore layered: across regimes via $\sum_j \pi_{ij}P(t,j)$3, across hierarchy through the dependence on the follower solution $\sum_j \pi_{ij}P(t,j)$4, and across forward-backward variables through the decoupling relation $\sum_j \pi_{ij}P(t,j)$5 (Wu et al., 2024).

The solvability theory is correspondingly partial. Under (H1)–(H4), the leader’s reduced BSLQ problem has a unique optimal control and the matrix $\sum_j \pi_{ij}P(t,j)$6 is shown to be positive definite. For the leader CDRE itself, the paper proves a conditional inverse-limit construction through a family $\sum_j \pi_{ij}P(t,j)$7 and obtains a direct solvability theorem only under the one-dimensional assumption (H5), where the CDRE can be rewritten as the Riccati equation of a standard regime-switching forward SLQ problem (Wu et al., 2024).

The examples illustrate both explicit and implicit solvability. In a scalar two-regime example,

$\sum_j \pi_{ij}P(t,j)$8

while a second scalar example is rewritten into the one-dimensional FSLQ-type Riccati form and solved numerically (Wu et al., 2024).

4. Mean-field coupled scalar complex Riccati systems

A different use of coupled Riccati dynamics appears in the paper "Low Dimensional Dynamics of Globally Coupled Complex Riccati Equations: Exact Firing-rate Equations for Spiking Neurons with Clustered Substructure" (Pazó et al., 5 Mar 2025). The basic microscopic system is an ensemble of $\sum_j \pi_{ij}P(t,j)$9 heterogeneous complex ODEs

$\sum_k \lambda_{ik}\Sigma(s,k)$0

with coefficients depending self-consistently on the global mean field

$\sum_k \lambda_{ik}\Sigma(s,k)$1

The paper is explicit that this is a globally coupled system of scalar complex Riccati equations, not a matrix Riccati equation (Pazó et al., 5 Mar 2025).

In the thermodynamic limit $\sum_k \lambda_{ik}\Sigma(s,k)$2, the state is described by a conditional density $\sum_k \lambda_{ik}\Sigma(s,k)$3 satisfying a continuity equation. The reduction proceeds in two stages. First, the exact invariant density ansatz

$\sum_k \lambda_{ik}\Sigma(s,k)$4

reduces the PDE to evolution equations

$\sum_k \lambda_{ik}\Sigma(s,k)$5

Second, under the restriction

$\sum_k \lambda_{ik}\Sigma(s,k)$6

together with Lorentzian heterogeneity

$\sum_k \lambda_{ik}\Sigma(s,k)$7

the mean-field integral closes by residue calculus: $\sum_k \lambda_{ik}\Sigma(s,k)$8 The resulting finite-dimensional system is

$\sum_k \lambda_{ik}\Sigma(s,k)$9

The paper describes this as an exact invariant-manifold reduction plus residue-calculus closure, closely related in philosophy to Ott–Antonsen/Lorentzian-ansatz reductions but specialized to the Möbius/Riccati structure of complex Riccati dynamics. It also records important limitations: heterogeneity is placed only in the real part of $z_j$0, analyticity assumptions are mathematically delicate because $z_j$1 enters through $z_j$2, and if the pole-selection sign changes during evolution, “our approach breaks down” (Pazó et al., 5 Mar 2025).

The neuroscience application is clustered QIF dynamics. After a first reduction within each subpopulation and the special internal coupling $z_j$3 with $z_j$4, the cluster variable

$z_j$5

satisfies

$z_j$6

so each cluster becomes one Riccati unit in a higher-level globally coupled ensemble. In the asymptotic regime, $z_j$7, and the dynamics collapse to two-dimensional firing-rate equations for global observables (Pazó et al., 5 Mar 2025).

5. Operator-valued and geometric generalizations

In infinite-dimensional PDE control, the relevant Riccati equation may be a single operator DRE associated with a coupled plant rather than a finite family of coupled matrix equations. The paper "Uniqueness for Riccati equations with unbounded operator coefficients" studies the control system

$z_j$8

with $z_j$9 generating a $Z,Q,\xi$0-semigroup on a separable Hilbert space $Z,Q,\xi$1 and $Z,Q,\xi$2. The finite-horizon Riccati equation is formulated weakly as

$Z,Q,\xi$3

and the infinite-horizon ARE is

$Z,Q,\xi$4

The central result is uniqueness in carefully defined classes $Z,Q,\xi$5 and $Z,Q,\xi$6, not unrestricted weak uniqueness. The paper is explicit that it concerns a single operator-valued DRE/ARE attached to a coupled PDE system, not a classical finite-dimensional system of several coupled Riccati ODEs (Acquistapace et al., 2020).

A geometric generalization is developed in "A geometric approach to differential Hamiltonian systems and differential Riccati equations". For the nonlinear control system

$Z,Q,\xi$7

the Riccati object is a state-dependent symmetric matrix field $Z,Q,\xi$8 defining a Lagrangian graph subbundle

$Z,Q,\xi$9

Invariance of that subbundle under the differential Hamiltonian system yields the generalized differential Riccati equation

$Z$0

together with the additional invariance equations

$Z$1

The paper states explicitly that this is not standard coupled differential Riccati equations in the conventional multi-equation sense; rather, it is a geometric generalized differential Riccati framework closely related to control contraction metrics and incremental stability/stabilizability (Schaft, 2015).

6. Solvability themes, stationary limits, and common distinctions

Several recurrent themes organize the theory. One is the passage from finite-horizon differential systems to stationary algebraic limits. In the two-player Nash setting, the 2020 CARE paper treats the stationary nonlinear algebraic fixed point that one would expect as the limit of time-dependent coupled Riccati dynamics, while explicitly not analyzing backward integration, time discretization, or stability of the CDRE flow itself (Li et al., 2020). In the zero-sum regime-switching game paper, this limiting passage is fully developed: $Z$2 converges exponentially to the stabilizing CARE solution $Z$3, and the gain convergence $Z$4 drives the turnpike estimate for the optimal state and controls (Li et al., 11 Sep 2025).

A second theme is that exact reduction is often tied to an invariant structure. For globally coupled complex Riccati equations, the key object is an invariant density manifold and a residue-calculus closure, with exact asymptotic dynamics when $Z$5 (Pazó et al., 5 Mar 2025). For geometric differential Riccati equations, the decisive structure is an invariant Lagrangian subbundle of $Z$6 (Schaft, 2015). For operator Riccati equations with unbounded coefficients, the decisive structure is a uniqueness class reflecting gain regularity on a fractional domain $Z$7 rather than an unrestricted strong operator identity (Acquistapace et al., 2020).

A third theme is that solvability can be partial and representation-dependent. The leader CDRE in the Stackelberg problem is solved generally only through a conditional inverse-limit construction, whereas a direct theorem is obtained under the one-dimensional hypothesis (H5) (Wu et al., 2024). This contrasts with the regime-switching zero-sum game, where strong regularity and stabilizing solvability are part of the main finite- and infinite-horizon theory (Li et al., 11 Sep 2025).

A common misconception is to treat all coupled Riccati equations as continuous-time CDREs. The paper "Harmonic-Copuled Riccati Equations and its Applications in Distributed Filtering" is explicit that its HCRE system is discrete-time and algebraic: the recursion

$Z$8

converges to a unique fixed point, but there are no actual differential Riccati equations in that paper (Qian et al., 2022). Similarly, the PDE uniqueness paper concerns operator DREs for coupled systems rather than a finite-dimensional family of mutually coupled Riccati ODEs (Acquistapace et al., 2020).

Taken together, these works show that CDREs are best understood through their coupling mechanism. In stochastic games, the coupling is typically through modes, players, and stationarity operators. In mean-field complex Riccati systems, it is through self-consistent coefficients driven by the population mean. In PDE control, it is the coupled plant that lifts a single Riccati equation to a product-space problem. In geometric formulations, coupling occurs between state, tangent, and cotangent dynamics. This suggests that the unifying content of CDRE theory is not a single canonical equation, but a class of Riccati differential structures whose coefficients, domains, or unknowns are constrained by a coupled dynamical architecture.