---
title: LQR-Schrödinger Bridge
url: https://www.emergentmind.com/topics/lqr-schrodinger-bridge
type: topic
---

# LQR-Schrödinger Bridge

The LQR-Schrödinger Bridge problem synthesizes the classical Schrödinger bridge—where the task is to stochastically steer probability flows between prescribed marginals—with the full structure of the linear quadratic regulator (LQR), where the cost functional is quadratic in both control and state. This unification yields a versatile paradigm for exact distributional steering in continuous and discrete time, accommodating linear system dynamics, general quadratic costs, and non-Gaussian initial and terminal distributions. It establishes deep connections between stochastic control, optimal transport, reaction-diffusion partial differential equations (PDEs), and integrable Markov kernels, with rigorous analytic and algorithmic tractability.

## 1. Formulation and Problem Statement

Consider a controlled linear diffusion in $\mathbb{R}^n$ over a finite time horizon $[0,T]$. The state evolves according to
\[
dx_t = A(t)x_t\,dt + B(t)u_t\,dt + \sqrt{\epsilon}\,B(t)\,dw_t,
\]
where $A(t)$ is an $n\times n$ matrix, $B(t)$ is $n\times m$ and $(A,B)$ is uniformly controllable, and $w_t$ is an $m$-dimensional standard Wiener process. For $\epsilon>0$, the dynamics include noise. The initial and final state distributions are prescribed, typically as zero-mean Gaussians $x_0\sim\mathcal{N}(0,\Sigma_0)$ and $x_T\sim\mathcal{N}(0,\Sigma_1)$ for suitable covariance matrices.

The cost to be minimized is
\[
J(u) = \mathbb{E}\left[\int_0^T (u_t^\top R(t)u_t + x_t^\top Q(t)x_t)\,dt\right],
\]
where $R(t)$ is positive definite and $Q(t)$ symmetric (not necessarily positive definite). The control $u_t$ belongs to the set of adapted processes of finite energy.

The task: Find $u_t$ steering the distribution of $x_t$ from $\Sigma_0$ to $\Sigma_1$ (or, in general, between $\mu_0$ and $\mu_1$ with finite second moments), minimizing $J(u)$ subject to the dynamics above [1608.03622][2504.15753].

## 2. Coupled Riccati Equations and LQR Structure

Under Gaussian endpoint distributions, optimal controls are linear state feedbacks:
\[
u_t^* = -R(t)^{-1}B(t)^\top\Pi(t)x_t,
\]
where $\Pi(t)$ solves the matrix Riccati ODE
\[
-\dot{\Pi}(t) = A(t)^\top\Pi(t) + \Pi(t)A(t) - \Pi(t)B(t)R(t)^{-1}B(t)^\top\Pi(t) + Q(t),
\]
subject to two-point boundary conditions involving auxiliary matrix trajectories. The dual equation, in terms of the optimal state covariance $M(t)$, gives a coupled forward Riccati:
\[
\dot{M}(t) = A(t)M(t) + M(t)A(t)^\top + \epsilon B(t)B(t)^\top - M(t)B(t)R(t)^{-1}B(t)^\top M(t),
\]
with $M(0)$ and $M(T)$ involving $\Sigma_0$ and $\Sigma_1$. The unique solution is determined by projection onto a $2n$-dimensional Hamiltonian flow and selecting branches to keep Gramian matrices nonsingular [1608.03622].

For general, non-Gaussian endpoint laws, the structure holds with dynamics and Sinkhorn recursions performed in functional space [2504.15753].

## 3. Markov Kernels, Reaction–Diffusion, and Exact Solvability

The core technical advancement is the construction of exact Markov kernels (Green’s functions) for the underlying reaction-diffusion PDEs encoding both advection/diffusion and quadratic killing (state cost). For time-invariant $(A,B,Q)$, the kernel takes the form
\[
\kappa(t_0, x;\, t, y) = \frac{1}{(2\pi)^{n/2}\sqrt{\det M(t_0, t)}}\;
\exp\left(-\frac{1}{2}(y - \Phi(t, t_0)x)^\top M(t_0, t)^{-1}(y - \Phi(t, t_0)x) - \frac{1}{2}x^\top P_{t_0}x + \frac{1}{2}y^\top P_t y\right),
\]
with $P_t$ the solution of the backward Riccati ODE, $\Phi$ the state transition of the closed-loop system, and $M$ the associated controllability Gramian under optimal feedback [2504.15753]. For the scalar case with $A=0$, $B=I$, $Q\succ0$, the kernel reduces to a non-radial Gaussian, revealing the connection to the Mehler kernel and quantum harmonic oscillator propagators [2406.00503][2407.15245].

Exactly solvable cases enable contractive Sinkhorn-type iterations for coupling arbitrary distributions at endpoints, with the kernel parameterizing the effective cost-to-go/distance functional arising from deterministic optimal control.

## 4. Schrödinger System, Sinkhorn Recursion, and Duality

The time-marginal dynamics and scaling functions underlying the optimal process satisfy a pair of coupled reaction–diffusion PDEs,
\[
\begin{aligned}
\partial_t \widehat\varphi &= (L - q)\,\widehat\varphi, \\
\partial_t \varphi &= (L^* + q)\,\varphi,
\end{aligned}
\]
where $q(x,t) = \frac{1}{2}x^\top Q_t x$ and $L$ is the generator of the uncontrolled process. The product $\rho(x, t) = \widehat\varphi(x, t)\varphi(x, t)$ gives the optimally steered density [2504.15753][2406.00503].

Iterating the backward and forward projections with product normalization realizes a continuous-time Sinkhorn recursion in function space. For Gaussian or Gaussian-mixture approximations, all intermediate function evaluations, including the projections, may be computed in closed form. These recursions converge exponentially fast in Hilbert’s projective metric [2406.00503][2506.17273].

## 5. Discrete-Time Formulation and Generalizations

The problem admits an exact discrete-time analogue. Given trajectories $\{x_0, x_1, ..., x_N\}$, with stepwise cost
\[
J(x_{0:N}) = \sum_{t=0}^{N-1}\left[x_t^T Q_t x_t + (x_{t+1} - x_t)^T R_t (x_{t+1}-x_t)\right] + x_N^T Q_N x_N,
\]
the optimal path probability factorizes via (dynamic programming) Gibbs–Kantorovich potentials $\phi_0$ and $\phi_N$ as
\[
p^*(x_{0:N}) \propto \phi_0(x_0)\left[\prod_{t=0}^{N-1} r_t(x_t, x_{t+1})\right]\phi_N(x_N),
\]
with $r_t(x_t, x_{t+1})$ the (non-normalized) reference measure [2506.17273]. Dual Riccati recursions for the cost-to-go potentials enable efficient forward and backward passes, resulting in Markovian transition kernels and explicit Gaussian (or generalized Gaussian) marginals.

The discrete-time LQR–Schrödinger bridge is extensible to non-Gaussian endpoint marginals via iterative proportional fitting, and generalizes Bures–Wasserstein transport by introducing time-indexed quadratic potentials that effect “curvature” in the path geometry [2506.17273].

## 6. Connections: Quantum Mechanics, Optimal Transport, and Beyond

The forward–backward PDE system underlying the LQR-Schrödinger bridge, under Wick rotation, is equivalent to the Schrödinger equation for the quantum harmonic oscillator with quadratic potential. The fundamental solution is the Mehler kernel, matching the propagator for the quantum oscillator. The closed-form transition kernel structure provides a unified framework that encompasses classical heat kernel diffusion, stochastic control, quantum mechanics, and entropy-regularized optimal transport [2406.00503][2407.15245].

In the zero-noise limit ($\epsilon \to 0$), the problem converges to a deterministic boundary-value optimal transport problem with general quadratic cost:
\[
\inf_{u, x} \int_0^T \left[u(t)^\top R(t) u(t) + x(t)^\top Q(t) x(t)\right] dt,
\]
recovering the quadratic-cost mass transport with specified marginals. The solution path is then entirely determined by the boundary conditions and the LQR structure [1608.03622].

## 7. Numerical Methods and Implementation Remarks

The solution pipeline for the LQR–Schrödinger bridge in both continuous and discrete time is characterized by:
- Backward integration of the matrix Riccati ODE (or recursion) for the value/cost-to-go matrices.
- Forward computation of the closed-loop state transition and Gramian.
- Construction of the Markov kernel as a parameterized non-homogeneous Gaussian (or quasi-Gaussian) in $(x,y)$.
- Iterative Sinkhorn-type proportional fitting in function space to enforce general (including non-Gaussian) endpoint distributions.

For practical purposes, each recursion (per time-step) is $O(d^3)$, and convergence is empirically rapid; typically, a small number of iterations is sufficient. For mixture-of-Gaussians approximations, all required integrals decompose into explicit formulas. In non-quadratic or non-Gaussian contexts, numerical quadrature or Monte Carlo methods can be used to approximate projections, leveraging the explicit form of the kernels [2504.15753][2406.00503][2506.17273].

---

**Key References:**
- Chen–Georgiou–Pavon, "Optimal steering of a linear stochastic system to a final probability distribution, Part III" [1608.03622]
- Blom, Chen, Georgiou, Pavon, "Markov Kernels, Distances and Optimal Control" [2504.15753]
- Chen, Georgiou, Pavon, "Schrödinger Bridge with Quadratic State Cost is Exactly Solvable" [2406.00503]
- Chen, Pavon, "Weyl Calculus and Exactly Solvable Schrödinger Bridges with Quadratic State Cost" [2407.15245]
- Pavon et al., "The LQR-Schrödinger Bridge" [2506.17273]

Source: https://www.emergentmind.com/topics/lqr-schrodinger-bridge