---
title: Non-symmetric Riccati Equation
url: https://www.emergentmind.com/topics/non-symmetric-riccati-equation
type: topic
---

# Non-symmetric Riccati Equation

Non-symmetric Riccati equation denotes a family of quadratic matrix equations in which the unknown matrix enters bilinearly, but the coefficient structure does not enforce the symmetry assumptions of the classical symmetric Riccati theory. In the recent cone-theoretic formulation, the nonsymmetric algebraic Riccati equation is
$$
XBX+DX+XA+C=0,
$$
and a stabilizing solution is one for which \(A_{\rm cl}=A+BX\) and \(D_{\rm cl}=D+XB\) are both stable [2411.14470]. Across transport theory, M-matrix analysis, T-Sylvester generalizations, differential equations, and forward-backward stochastic control, the subject includes minimal nonnegative solutions, cone-preserving solutions, invariant-subspace descriptions, and numerical methods ranging from monotone iteration and Newton–Kleinman procedures to doubling, projection, and low-rank ADI algorithms [1212.6461, 2003.03693, 2604.23208, 2504.14288].

## 1. Principal equation classes

The literature uses several closely related nonsymmetric Riccati models. They differ in sign convention, linear operator, and intended application, but all retain a quadratic term in the unknown matrix and a nonsymmetric coefficient structure.

| Variant | Equation | Setting |
|---|---|---|
| Cone-preserving algebraic Riccati equation | \(XBX+DX+XA+C=0\) | Proper cones and stability [2411.14470] |
| Transport-theory NARE | \(AX+XB-XCX+E=0\) | Large-scale transport theory [1407.2520] |
| M-matrix / linearization form | \(XCX-XD-AX+B=0\) | Transport theory and M-matrix analysis [1109.5006, 1011.1363] |
| Nonsymmetric T-Riccati equation | \(DX+X^TA-X^TBX+C=0\) | T-Sylvester/T-Lyapunov operator [2003.03693] |
| Nonsymmetric differential Riccati equation | \(\dot X=-AX-XD+XSX+Q\) | Large-scale time-dependent problems [1801.01291] |
| Equilibrium Riccati equations | coupled equations for \(P_1(t,s)\) and \(P_2(s)\) | Forward-backward stochastic LQ control [2504.14288] |

A further large-scale generalized NARE is written as
$$
AX\hat E+EX\hat A-EX\hat B C X\hat E+B\hat C=0,
$$
with \(X\in\mathbb R^{n\times\hat n}\), low-rank data, and dimensions \(n,\hat n\gg1\) in the applications emphasized in the low-rank ADI literature [2604.23208].

The T-Riccati equation replaces the standard Sylvester operator by the T-Sylvester operator
$$
\mathcal S_T(X)=DX+X^TA,
$$
and reduces to the T-Lyapunov equation when \(B=0\) [2003.03693]. The differential form introduces a time derivative and is treated by projection onto low-dimensional subspaces, after which a reduced Riccati differential equation is solved [1801.01291].

## 2. Order structures, positivity, and solution concepts

A central distinction in nonsymmetric Riccati theory is the ordering used to define admissible solutions. In the M-matrix literature, the main object is the minimal nonnegative solution \(X_{\min}\), meaning the entrywise smallest nonnegative solution [1109.5006, 1212.6461]. In that framework, the notion of minimality replaces the usual positive-semidefinite ordering of the symmetric theory [1212.6461].

The cone-theoretic formulation replaces entrywise nonnegativity by an arbitrary proper cone \(K\subset\mathbb R^n\). A proper cone is convex, closed, pointed, and has nonempty interior, and it induces the partial order
$$
x\succeq_K y\quad\Longleftrightarrow\quad x-y\in K.
$$
A matrix \(M\) is \(K\)-nonnegative, or cone-preserving, if \(MK\subseteq K\). A solution \(X\) of the Riccati equation is cone-preserving if \(X\succeq_{\pi(K)}0\), where \(\pi(K)\) is the cone of linear operators preserving \(K\) [2411.14470]. When \(K=\mathbb R_+^n\), this reduces to entrywise nonnegativity.

The same framework introduces cross-positivity. If \(K^*\) is the dual cone, then \(M\) is cross-positive on \(K\) if
$$
x\in K,\quad y\in K^*,\quad y^\top x=0
\quad\Longrightarrow\quad
y^\top(Mx)\ge 0.
$$
Equivalently,
$$
M\text{ is cross-positive on }K
\quad\Longleftrightarrow\quad
e^{tM}\succeq_{\pi(K)}0,\ \forall\,t\ge0.
$$
On the nonnegative orthant, cross-positivity means that all off-diagonal entries are nonnegative [2411.14470].

The T-Riccati setting also uses an order structure, but it is expressed directly through sign assumptions: \(B\ge0\), \(C\le0\), and the Kronecker representation of \(\mathcal S_T\) is a nonsingular M-matrix. Under these assumptions, fixed-point and Newton–Kleinman sequences are nondecreasing and converge to the minimal nonnegative solution [2003.03693].

A plausible implication is that entrywise nonnegativity is only one special case of a broader order-theoretic picture. The cone-preserving formulation makes this explicit by extending orthant-based results to arbitrary proper cones [2411.14470].

## 3. Existence theory, stability, and invariant subspaces

One of the sharpest recent existence criteria concerns the block matrix
$$
L=\begin{pmatrix}A&B\\ C&D\end{pmatrix}.
$$
If \(L\) is cross-positive on \(K\times K\) and stable, meaning \(\Re\lambda<0\) for every eigenvalue \(\lambda\) of \(L\), then the Riccati equation
$$
XBX+DX+XA+C=0
$$
admits a unique cone-preserving solution \(X\) such that \(A+BX\) and \(D+XB\) are both stable and cross-positive on \(K\). Conversely, if such a cone-preserving stabilizing solution exists, then \(L\) must be stable [2411.14470]. This result extends the nonnegative orthant case and turns a previously sufficient condition into an equivalence.

The M-matrix literature treats a different but related block structure. With
$$
K=\begin{pmatrix}A&-B\\ -C&D\end{pmatrix},
$$
the algebraic Riccati equation
$$
AX+XD-XBX+C=0
$$
has a minimal nonnegative solution when \(K\) is a nonsingular M-matrix or an irreducible singular M-matrix. A further extension shows that if \(K\) is a reducible singular M-matrix but regular, in the sense that there exists \(v>0\) such that \(Kv>0\), then the minimal nonnegative solution still exists and is unique [1212.6461].

Invariant-subspace formulations are fundamental in both analysis and computation. For the linearizing matrix
$$
\mathcal H=
\begin{bmatrix}
D & -C\\
B & -A
\end{bmatrix},
$$
a matrix \(X\) solves the NARE if and only if
$$
\mathcal H
\begin{bmatrix}
I_n\\ X
\end{bmatrix}
=
\begin{bmatrix}
I_n\\ X
\end{bmatrix}
(D-CX),
$$
so the columns of \(\bigl[I_n;\,X\bigr]\) span an invariant subspace of \(\mathcal H\) [1011.1363]. In the regular M-matrix case, the spectrum of the related linearization splits across the imaginary axis, and the stable eigenvalues coincide with the spectrum of \(D-CX\) [1212.6461].

A recurrent difficulty is the critical or close-to-critical regime. In the critical transport-theory case \((a,c)=(0,1)\), the block matrix \(M\) is singular, factors such as \(I-G_kH_k\) may become singular, and doubling iterations may stall or fail [1109.5006]. The close-to-critical case arises when problematic eigenvalues are small but nonzero, producing ill-conditioning and slow convergence even though the singular critical theory no longer applies directly [1011.1363].

## 4. Iterative methods and structure-preserving algorithms

Monotone fixed-point iteration is the oldest common mechanism across the subject. In the cone-preserving algebraic Riccati problem, one uses
$$
D\,X_{i+1}+X_{i+1}A=-\bigl(X_iBX_i+C\bigr),\qquad X_0=0,
$$
obtaining a monotonically increasing sequence of cone-preserving matrices. A key auxiliary result states that such a sequence converges if it is bounded above in a single vectorial direction \(r\succ_K0\) [2411.14470]. In the T-Riccati setting, the iteration
$$
X_0=0,\qquad
\mathcal S_T(X_{k+1})=X_k^TBX_k-C
$$
is well defined under the M-matrix assumptions and converges to the minimal nonnegative solution when an appropriate nonnegative upper bound exists [2003.03693].

Newton-type methods provide faster local convergence. For the T-Riccati equation, the Fréchet derivative at \(X\) is
$$
\mathcal R_T'(X)[H]=(D-X^TB)H+H^T(A-BX),
$$
and the Newton–Kleinman update solves a T-Sylvester equation whose coefficients are shifted by the current iterate. Under the same sign and M-matrix assumptions, the iterates remain nonnegative, nondecreasing, and bounded above by the minimal nonnegative solution, and they converge monotonically to that solution [2003.03693]. For classical M-matrix algebraic Riccati equations, Newton’s method is locally quadratic when the shifted coefficients remain nonsingular [1212.6461].

The structure-preserving doubling algorithm (SDA) is the main high-order method in transport theory and M-matrix settings. Starting from a Cayley-transformed initialization, the recurrences
$$
\begin{aligned}
E_{k+1}&=E_k(I-G_kH_k)^{-1}E_k,\\
F_{k+1}&=F_k(I-H_kG_k)^{-1}F_k,\\
G_{k+1}&=G_k+E_k(I-G_kH_k)^{-1}G_kF_k,\\
H_{k+1}&=H_k+F_k(I-H_kG_k)^{-1}H_kE_k
\end{aligned}
$$
converge quadratically in the nonsingular M-matrix case [1407.2520, 1109.5006]. In regular singular M-matrix problems, SDA and ADDA remain well defined and quadratically convergent under the simple-zero-eigenvalue assumption on the linearization [1212.6461].

Shift strategies were developed to remove the spectral obstructions behind breakdown and slow convergence. In the critical transport-theory case, one- and two-parameter rank-one shifts move the zero eigenvalues of the block matrix while leaving the minimal nonnegative solution unchanged; the shifted SDA is then guaranteed to converge quadratically with no breakdown [1109.5006]. The subspace-shift technique generalizes this idea to close-to-critical problems by shifting an entire invariant subspace associated with the problematic eigenvalues. If \(V\) and \(U\) span the right and left central invariant subspaces, then
$$
\widehat{\mathcal H}
=
\mathcal H\bigl(I+sV(U^*V)^{-1}U^*\bigr)
$$
scales the selected eigenvalues by \(1+s\), preserves the antistable subspace, and therefore preserves the minimal Riccati solution while improving the convergence rate of doubling iterations [1011.1363].

## 5. Large-scale, low-rank, and projected computation

Large-scale nonsymmetric Riccati problems typically arise when the coefficient matrices are sparse or structured and the solution has low numerical rank. In the transport-theory NARE
$$
AX+XB-XCX+E=0,
$$
the matrices \(B\) and \(C\) are symmetric and low-ranked, while \(A\) and \(E\) are rank-one updates of nonsingular diagonal matrices. This permits \(O(n)\) solves with shifted \(A\) and \(E\) by the Sherman–Morrison–Woodbury formula [1407.2520].

A balancing strategy can simplify the large-scale SDA. With the diagonal scaling \(\Phi=\operatorname{diag}(\sqrt{q_1},\dots,\sqrt{q_n})\), the transformed coefficients satisfy
\[
\widetilde B=\widetilde C^T,\qquad \widetilde A=\widetilde A^T,\qquad \widetilde E=\widetilde E^T,
\]
and carefully chosen low-rank initial factors imply
\[
Q_{1k}=P_{2k},\qquad Q_{2k}=P_{1k},\qquad E_k=E_k^T,\qquad F_k=F_k^T,\qquad H_k=G_k^T.
\]
As a consequence, the modified large-scale SDA reduces the flop operations of the original SDA\(_{ls}\) by about one half while retaining \(O(r^2n)\) complexity and \(O(rn)\) memory per iteration [1407.2520].

Projection methods are prominent for nonsymmetric differential Riccati equations. Extended block Arnoldi bases \(\mathcal V_m\) and \(\mathcal W_m\) are built for \((A,F)\) and \((D,G)\), and the approximation
\[
X_m(t)=\mathcal V_mY_m(t)\mathcal W_m^T
\]
is obtained by solving a reduced Riccati differential equation for \(Y_m(t)\). The reduced problem may be integrated by matrix exponential formulas, BDF schemes, or a two-stage Rosenbrock method [1801.01291]. For transport-theory differential problems, the rank-one perturbation structure of \(A\) and \(D\) again makes the required linear solves cheap via Sherman–Morrison–Woodbury [1801.01291].

The most general low-rank large-scale development in the supplied literature is a low-rank ADI method for the generalized NARE
$$
AX\hat E+EX\hat A-EX\hat B C X\hat E+B\hat C=0.
$$
The approximate solution is stored as
\[
\widetilde X^{(i)}=V^{(i)}\bar X^{(i)}(\widehat W^{(i)})^T,
\]
and each iteration requires two sparse shifted linear solves plus low-rank updates. The residual factors as
\[
R^{(i)}=B_\perp^{(i)}\widehat C_\perp^{(i)},
\]
which supports Petrov–Galerkin residual orthogonality and autonomous shift generation through small projected eigenproblems [2604.23208]. Numerical results include a benchmark example of order \(10^6\), and the method is presented as covering Lyapunov equations, Sylvester equations, and symmetric Riccati equations as special cases [2604.23208].

## 6. Applications, variants, and conceptual boundaries

Transport theory is the most persistent application domain in the supplied literature. Both the classical NARE
\[
XCX-XD-AX+B=0
\]
and the alternative form
\[
AX+XB-XCX+E=0
\]
arise from models with rank-one or low-rank scattering terms and diagonal-plus-rank-one coefficients [1109.5006, 1407.2520]. Large-scale numerical work is organized around this structure, which permits efficient shifts, balancing, and low-rank updates.

The subject also extends beyond transport and queueing models. The nonsymmetric differential Riccati equation is stated to arise in control theory, transport theory, applied probability, and others [1801.01291]. The nonsymmetric T-Riccati equation has applications in macroeconomics and policy dynamics, where the linear part is governed by the T-Sylvester operator rather than the standard Sylvester operator [2003.03693].

A distinct development appears in forward-backward stochastic linear-quadratic control. There, time-inconsistency invalidates the classical dynamic-programming principle, and equilibrium strategies are characterized by a coupled system of matrix-valued, non-local ordinary differential equations with a non-symmetric structure. The unknowns are \(P_1(t,s)\in S^n\) and \(P_2(s)\in\mathbb R^{m\times n}\), and the feedback law depends on the combination
\[
V(s)=P_1(s,s)+P_2(s)^TG_2(s)P_2(s).
\]
Under uniform boundedness, uniform positive-definiteness, and monotonicity assumptions, the equilibrium Riccati system admits a unique global solution and satisfies the a priori bound
\[
\sup_{s\in[0,T]}\|V(s)\|\le C_\delta
\]
[2504.14288]. The same source gives a one-dimensional example in which, without the positivity assumptions, the scalar ERE blows up in finite backward time [2504.14288].

The relation to the symmetric Riccati equation is therefore subtle. In the forward-backward stochastic setting, if \(H=0\), \(M=N=0\), and \(G_2=0\), then \(P_2\equiv0\) and the equilibrium Riccati system reduces to the standard symmetric Riccati ODE [2504.14288]. In the M-matrix setting, the comparison is formulated differently: the nonsymmetric problem has no symmetry such as \(A=D^T\) or \(B=C^T\), and minimality is interpreted componentwise rather than through positive-semidefinite order [1212.6461]. This suggests that the nonsymmetric theory is not a single perturbation of the symmetric case, but a collection of structurally different Riccati problems linked by quadratic matrix nonlinearity, order-preserving mechanisms, and invariant-subspace geometry.

Source: https://www.emergentmind.com/topics/non-symmetric-riccati-equation