Stochastic Optimal Control provides a unified mathematical framework for solving complex decision-making problems, encompassing paradigms such as maximum entropy reinforcement learning(RL) and imitation learning(IL). However, conventional parametric policies often struggle to represent the multi-modality of the solutions. Though diffusion-based policies are aimed at recovering the multi-modality, they lack an explicit probability density, which complicates policy-gradient optimization. To bridge this gap, we propose MePoly, a novel policy parameterization based on polynomial energy-based models. MePoly provides an explicit, tractable probability density, enabling exact entropy maximization. Theoretically, we ground our method in the classical moment problem, leveraging the universal approximation capabilities for arbitrary distributions. Empirically, we demonstrate that MePoly effectively captures complex non-convex manifolds and outperforms baselines in performance across diverse benchmarks.
Figures & tables
Figure 1: Conceptual comparison of policy parameterizations on a non-convex action manifold. (A) A unimodal Gaussian concentrates mass around a single mean, leading to limited expressivity and mode collapse. (B) Diffusion can represent complex supports via iterative sampling but lacks a tractable likelihood and requires multi-step generation. (C) MePoly yields an expressive, explicit density that can conform to the manifold while retaining tractable log-probabilities and entropies for learning.
Figure 2: Trajectory samples on multi-modal navigation tasks: Each column is a different environment. Each row corresponds to various methods: MePoly (ours), PPO (Gaussian), PPO (Gaussian Mixture Model), and FPO (Flow-Matching). In every panel, we visualize multiple rollouts from the same start state (white circle) under identical environment layouts: walls/obstacles are shown in gray, goal regions in green, and death/unsafe regions in red (see legend). Trajectories are color-coded by time (purple → yellow), and black markers indicate terminal positions. MePoly consistently produces diverse, distinct feasible routes that cover multiple classes (e.g., different goals and passages through slits/around obstacles), whereas baselines often collapse to a single mode or fail to represent alternative valid solutions in highly non-convex environments.
Figure 3: Myth of Mepoly: The LHS shows that Mepoly asymptotically approaches the sample distribution as the order increases, which validates Theorem 3.4 and Corollary 3.5 . The RHS replaces Legendre Polynomials with standard Polynomials, samples become noticeably blurrier and less faithful, indicating that orthogonal bases are critical for stable and accurate manifold approximation.
Figure 4: Sample quality comparison: MePoly (ours) generates samples that closely follow the underlying curved support and preserve the global topology. In contrast, the baselines fail to capture the non-convex and multi-modal structures, resulting in significant geometric mismatch and a collapse to unimodal or fragmented approximations.
Method
PushCube
PickCube
StackCube
PushT
MePoly
0.901/0.018
0.945/0.025
0.446/0.049
0.279/0.020
DP
0.856/0.015
0.992/0.002
0.410/0.001
0.298/0.013
Table 1: Comparative analysis of MePoly performance. Top : ManiSkill manipulation tasks demonstration. Bottom : Quantitative success rate comparison on ManiSkill benchmarks against the Diffusion Policy (DP). Results are aggregated over three random seeds; for each seed, we report the peak success rate achieved during training, and summarize performance as mean/std across seeds.
Efficient exploration is a central problem in reinforcement learning and is often formalized as maximizing the entropy of the state-action occupancy measure. While unconstrained maximum-entropy exploration is relatively well understood, real-world exploration is often constrained by safety, resource, or imitation requirements. This constrained setting is particularly challenging because entropy maximization lacks additive structure, rendering Bellman-equation-based methods inapplicable. Moreover, scalable approaches require policy parameterization, inducing non-convexity in both the objective and the constraints. To our knowledge, the only prior model-free policy-gradient approach for this setting under general policy parameterization is due to Ying et al. (2025). Unfortunately, their guarantees are limited to weak regret and ergodic averages, which do not imply that the final output is a single deployable policy that is near-optimal and nearly feasible. In this work we take a different approach to this problem, and propose Policy Gradient Penalty (PGP) method, a single-loop policy-space method that enforces general convex occupancy-measure constraints via quadratic-penalty regularization. PGP constructs pseudo-rewards that yield gradient estimates of the penalized objective, subsequently exploiting the classical Policy Gradient Theorem. We further establish the regularity of the penalized objective, providing the smoothness properties needed to justify the convergence of PGP. Leveraging hidden convexity and strong duality, we then establish global last-iterate convergence guarantees, attaining an ε-optimal constrained entropy value with ε bounded constraint violation despite policy-induced non-convexity. We validate PGP through ablations on a grid-world benchmark and further demonstrate scalability on two challenging continuous-control tasks.
Florian Wolf, Ilyas Fatkhullin, Niao He
The Computing & Mathematical Sciences Department, California Institute of Technology, Pasadena, CA. · Department of Computer Science, ETH Zurich, Switzerland. · 3ETH AI Center, ETH Zurich, Switzerland.
Balancing policy expressiveness with the exploration-exploitation trade-off is a core challenge in online Reinforcement Learning (RL). While Stochastic Differential Equation (SDE)-based diffusion policies can represent complex, multimodal action distributions, they suffer from two critical limitations: their stochastic reverse processes render entropy intractable (necessitating heuristic exploration), and computing policy gradients through long denoising chains is expensive and unstable. In this work, we show that ODE-based flow matching inherently resolves these issues by enabling both simulation-free policy optimization and tractable entropy computation. Building on this, we introduce Flow Matching Policy Optimization with Mirror Descent and Entropy Constraints (FMER). Our framework exploits this insight in three ways. First, we derive an analytic entropy objective that corrects for the density distortion caused by the tanh transformation (mapping an unbounded latent space to bounded actions), thereby facilitating entropy-constrained policy optimization. Second, we theoretically establish that minimizing an advantage-weighted conditional flow matching loss acts as a simulation-free surrogate for policy mirror descent. This steers the velocity field toward high-value regions while entirely avoiding backpropagation through the ODE solver. Finally, we dynamically tune the mirror descent temperature based on the effective sample size to enforce a robust trust region during training. Empirical evaluations demonstrate that FMER achieves superior performance on the challenging sparse-reward FrankaKitchen environment, while maintaining competitive results across standard dense-reward MuJoCo benchmarks. Code is available at https://github.com/gtanaelle/FMER.git.
Maximum entropy reinforcement learning (MaxEnt RL) has become a standard framework for sequential decision making, yet its standard Gaussian policy parameterization is inherently unimodal, limiting its ability to model complex multimodal action distributions. This limitation has motivated increasing interest in generative policies based on diffusion and flow matching as more expressive alternatives. However, incorporating such policies into MaxEnt RL is challenging for two main reasons: the likelihood and entropy of continuous-time generative policies are generally intractable, and multi-step sampling introduces both long-horizon backpropagation instability and substantial inference latency. To address these challenges, we propose Truncated Rectified Flow Policy (TRFP), a framework built on a hybrid deterministic-stochastic architecture. This design makes entropy-regularized optimization tractable while supporting stable training and effective one-step sampling through gradient truncation and flow straightening. Empirical results on a toy multigoal environment and 10 MuJoCo benchmarks show that TRFP captures multimodal behavior effectively, outperforms strong baselines on most benchmarks under standard sampling, and remains highly competitive under one-step sampling.