ConceptioArchivearXiv CS
arXiv CSopen access

Ada-Diffuser: Latent-Aware Adaptive Diffusion for Decision-Making

2026 · arxiv_cs
arXiv CS · Papers · License: Open Access · 2026
Open Source ↗Direct PDF ↓
knowledge-representationreasoning
artificial intelligence, reasoning, knowledge representation

Published as a conference paper at ICLR 2026

A DA -D IFFUSER : L ATENT-AWARE A DAPTIVE D IFFUSION FOR D ECISION -M AKING

arXiv:2605.16054v1 [cs.LG] 15 May 2026

Fan Feng1,3 , Selena Ge1 , Minghao Fu1 , Zijian Li2,3 , Yujia Zheng2 , Zeyu Tang2,4 , Yingyao Hu5† , Biwei Huang1† , Kun Zhang2,3† 1 University of California San Diego 2 Carnegie Mellon University 3 MBZUAI, 4 Stanford University, 5 Johns Hopkins University † Equal Senior Authorship

A BSTRACT Recent work has framed decision-making as a sequence modeling problem using generative models such as diffusion models. Although promising, these approaches often overlook latent factors that exhibit evolving dynamics, elements that are fundamental to environment transitions, reward structures, and high-level agent behavior. Explicitly modeling these hidden processes is essential for both precise dynamics modeling and effective decision-making. In this paper, we propose a unified framework that explicitly incorporates latent dynamic inference into generative decision-making from minimal yet sufficient observations. We theoretically show that under mild conditions, the latent process can be identified from small temporal blocks of observations. Building on this insight, we introduce Ada-Diffuser, a causal diffusion model that learns the temporal structure of observed interactions and the underlying latent dynamics simultaneously, and furthermore, leverages them for planning and control. With a modular design, Ada-Diffuser supports both planning and policy learning tasks, enabling adaptation to latent variations in dynamics, rewards, and latent actions. Experiments on locomotion and robotic manipulation benchmarks demonstrate its effectiveness in accurate latent inference, long-horizon planning, and adaptive policy learning1 .

1

I NTRODUCTION

Learning and planning in partially observable environments is a fundamental challenge in building intelligent agents (Kaelbling et al., 1998). Recent work on casting decision-making as a generative modeling problem, taking advantage of powerful models such as transformers (Chen et al., 2021; Zheng et al., 2022; Kong et al., 2024) and diffusion models (Janner et al., 2022; Chi et al., 2023; Ren et al., 2025), has achieved impressive results in a wide range of tasks. However, these methods often fail to account for hidden latent variables and their temporal dynamics, factors that are prevalent in real-world settings such as robotics (Lauri et al., 2022), autonomous driving (Huang et al., 2024), healthcare (Hauskrecht & Fraser, 2000; Ehrmann et al., 2023), and economics (Brero et al., 2022). Ignoring such latent processes can result in suboptimal decision-making, particularly when the observational data does not provide full coverage of the latent factors underlying the environment’s dynamics (Zintgraf et al., 2021; Xie et al., 2021; Swamy et al., 2022; Belkhale et al., 2023). Early works address partial observability in reinforcement learning (RL) and imitation learning (IL) by encoding historical observations and actions into belief states or latent embeddings, which represent a distribution over the underlying latent state (Kaelbling et al., 1998; Hauskrecht, 2000; Guo et al., 2018; Igl et al., 2018; Liang et al., 2024a; Xie et al., 2021). Policy optimization or planning is then carried out based on these inferred belief states. However, learning such representations often requires access to the historical trajectories or data from a diverse set of environments. This can be prohibitively expensive, particularly in high-dimensional state or action spaces, posing challenges for integrating these methods into modern generative decision-making models, which typically prioritize scalability. Can we identify the latent factors that govern environment dynamics and rewards, and 1

Project Page: https://sites.google.com/view/ada-diffuser.

1

Published as a conference paper at ICLR 2026

integrate them into scalable generative decision-making models to enable adaptive planning and policy learning, using only minimal observations, while preserving theoretical guarantees? In this paper, we pursue this goal by addressing two fundamental questions. First, what is the minimum set of observations required, in principle, to reliably identify the latent factors that govern the environment? Second, how can latent identification be effectively incorporated into generative models (e.g., diffusion models) to enable adaptive planning and policy learning? To answer the first question, we theoretically show that, under mild conditions, the latent factors at the time step t can be block-wise identified using only four surrounding observable measurements (i.e., state-action trajectories) within a small temporal window. This identification result implies that a small temporal block is sufficient to infer the latent factors in observational RL trajectories in an online manner. Guided by the theoretical findings, we propose Ada-Diffuser, a novel causal diffusion framework with latent identification from temporal blocks, designed to model the data generation process of RL trajectories influenced by latent factors. To reflect the autoregressive nature of sequential decision making, we introduce a causal denoising schedule that aligns the denoising steps with the underlying causal structure, drawing inspiration from recent advances in autoregressive diffusion models (Ho et al., 2022; Chen et al., 2024; Xie et al., 2024b; Sand-AI, 2025). For temporal-block-wise latent identification, during training, we propose a denoise-then-refine procedure that iteratively alternates between denoising the observations and refining latent estimates. This enables Ada-Diffuser to jointly learn a structured representation of latent variables and the corresponding observational distribution. At inference time, Ada-Diffuser generates actions and states while estimating latent variables in an online fashion. Since states and actions are conditioned on the latent factors, we employ a zig-zag sampling scheme that alternates between sampling state-action pairs and updating latent variables, ensuring consistency between generated sequences and their underlying latent dynamics. Ada-Diffuser provides a unified generative framework for sequential decision-making. It is applicable to both planning and policy learning tasks by conditioning on different types of observations and adapting the conditional generative process accordingly. The framework is flexible and can accommodate various forms of latent, including ones that influence dynamics, rewards, or even represent high-level latent actions. Importantly, even in environments without explicitly designed latent variables, the block-wise latent identification mechanism improves generative modeling by implicitly capturing structured temporal dependencies. Contributions: (1) We establish sufficient conditions under which latent factors influencing environment dynamics and rewards can be identified from short temporal windows of RL trajectories, without requiring full trajectory access or multi-environment data. (2) We develop Ada-Diffuser, a causal diffusion model that performs block-wise latent inference to jointly model latent contexts and observable trajectories. Unlike prior latent-augmented diffusion approaches, Ada-Diffuser introduces a minimal-sufficient block with backward refinement for identifiable latents and uses fully autoregressive denoising with zig–zag sampling to couple inference and generation. (3) Ada-Diffuser can be adapted to a wide range of decision-making tasks by conditioning on different types of observation. We empirically show the improved performance on a wide range of planning and control tasks, including 8 environments under 23 different settings.

2

BACKGROUND AND R ELATED W ORK

In this section, we provide background and related work on diffusion-based decision-making. Additional discussions are provided in Appendix E, including related work on (1) learning latent belief states in POMDPs (Kaelbling et al., 1998; Hauskrecht, 2000; Igl et al., 2018; Gregor et al., 2018; Goyal et al., 2021), particularly in the context of transfer, meta, and nonstationary RL/IL (Zintgraf et al., 2021; Xie et al., 2021; Feng et al., 2022; Ni et al., 2023; Liang et al., 2024a), and (2) autoregressive diffusion models (Chen et al., 2024; Xie et al., 2024b; Sand-AI, 2025; Wu et al., 2023). Recent advances use diffusion models as planners and policies for both RL and IL. I. Diffusion Planner: Diffusion-based planning leverages generative models to sample future state-action trajectories from a given state, using guidance techniques (Dhariwal & Nichol, 2021; Ho & Salimans, 2022) to encourage desirable properties such as high expected rewards. Taking Denoising Diffusion Probabilistic Models (DDPM (Ho et al., 2020))-based approaches as an example, these methods learn a generative model over expert trajectories τ = {(s0 , a0 ), . . . , (sT , aT )} by modeling a forward-noising process: 2

problem problem byby maximizing maximizing thethe likelihood likelihood ofof expert expert actions actions given given observed observed states, states, i.e., i.e., learning learning a policy a policy ω(a ω(a | s) | s) that that closely closely imitates imitates thethe expert expert policy ωeω(a (a | s). | s). Published as policy a conference paper at ICLR 2026 e

I. I.Diffusion DiffusionPlanner: Planner:Diffusion-based Diffusion-basedplanning planningleverages leveragesgenerative generativemodels modelsto tosample samplefuture future state-action state-actiontrajectories trajectoriesfrom froma agiven givenstate, state,using usingguidance guidancetechniques techniques[52, [52,53]53]to toencourage encourage desirable desirableproperties propertiessuch suchasashigh highexpected expectedrewards. rewards.Taking TakingDenoising DenoisingDiffusion DiffusionProbabilistic Probabilistic Models Models (DDPM (DDPM [54])-based [54])-based approaches approaches as as an an example, example, these these methods methods learn learn a generative a generative model model √ t t−1 t q(x | x{(s ) a), =.),N (x ;(s, aα, ta)} xt−1 ,by (1 − αt )I), aand a parameterizedprocess: denoising model pθ (xt−1 | xt ) to over overexpert experttrajectories trajectories ε ε = = {(s , a , . . . , . (s , )} by modeling modeling a forward-noising forward-noising process: 0 0 0 0 T T T T t t t→1 t→1 t → t →reverse t→1 t→1the process. Here, the superscript t denotes diffusion steps, t→1 t→1 Tt denotes t q(x q(x| x| x ) = )= NN (x(x; ; ϑtϑx0t x , (1, (1 ↑↑ ϑtϑ)I), and and a parameterized a parameterized denoising denoising model model pωp(x ) to ) to the planning horizon, t )I), ω (x | x| x x is a clean subsequence sampled from the expert trajectory τ , 0and 0 αt controls the variance schedule reverse reverse thethe process. process. Here, Here, gray gray t denotes t denotes diffusion diffusion steps, steps, T Tdenotes denotes thethe planning planning horizon, horizon, x xis is a clean a clean at diffusion step t. During inference, trajectories are generated by starting from Gaussian noise subsequence subsequence sampled sampled from from thethe expert expert trajectory trajectory ε , εand , and ϑtϑcontrols thethe variance variance schedule schedule at at diffusion diffusion t controls and iteratively denoising through the learned reverse process. This generation can be optionally step step t. t. During During inference, inference, trajectories trajectories areare generated generated byby starting starting from from Gaussian Gaussian noise noise and and iteratively iteratively conditioned on the initial state or other guidance signals y (e.g., goals, rewards): τ̂ ∼ pθ (τ | s0 , y). denoising denoising through through thethe learned learned reverse reverse process. process. This This generation generation can can bebe optionally optionally conditioned conditioned onon thethe II. Diffusion Policy: In contrast to diffusion planners, Diffusion Policy methods directly parameterize initial initial state state oror other other guidance guidance signals signals yy (e.g., (e.g., goals, goals, rewards): rewards): ε̂ ↓ ε̂ ↓ pωp(εω (ε| s|0s, 0y). , y). the policy πθ (a | s) using diffusion models. For example, Diffusion Policy (Chi et al., 2023) uses a diffusion model to generateDiffusion multi-step actions with expressive multimodal distributions. DPPO (Ren II.II. Diffusion Diffusion Policy: Policy: In In contrast contrast to to diffusion diffusion planners, planners, Diffusion Policy Policy methods methods directly directly parameterize parameterize et al., 2025) extends this idea by modeling a two-layer MDP structure, which enables fine-tuning of thethe policy policy ωωω(a (a | s) | s) using using diffusion diffusion models. models. For For example, example, Diffusion Diffusion Policy Policy [55] [55] uses uses a diffusion a diffusion ω diffusion-based policies in RLdistributions. settings. Another line of work usesthis diffusion model model to to generate generate actions actions with with expressive expressive multimodal multimodal distributions. DPPO DPPO [56] [56] extends extends this idea idea models to parameterize the policy networks forenables only thefine-tuning single current step (Wang et al., 2022; Hansen-Estruch et al., 2023; byby modeling modeling a two-layer a two-layer MDP MDP structure, structure, which which enables fine-tuning ofof diffusion-based diffusion-based policies policies in in Chen et al., 2023; Lu et al., 2023). Ada-Diffuser can generally accommodate both diffusion RLRL settings. settings.Another Another line line of of work work integrates integrates diffusion diffusion models models with with value-based value-based methods methods (e.g., (e.g., Q-Qplanner and policies within the same by framework. learning), learning), to to generate generate multimodal multimodal action action distributions distributions guided guided by learned learned value value functions, functions, such such asas Diffusion-QL Diffusion-QL [57], [57], IDQL IDQL [58], [58], CPQL CPQL [59], [59], CEP CEP [60], [60], and and DWM DWM [61]. [61].

3 L ATENT I DENTIFICATION IN POMDP 3 3 Latent LatentIdentification IdentificationininPOMDP POMDP

Informally this section, we seek to formally model the structure of thebydecision-making InIn this this section, section, wewe seek seek to to formally model model thethe structure structure ofof thethe decision-making decision-making system system by answering: answering: system by answering thereside following questions. First, where do the latent factorssuch reside, and how do they influence the (1). (1). Where Where thethe latent latent factors factors reside and and how how they they influence influence thethe observable observable variables variables such asas states, states, observable variables such states, actions, and rewards? data Second, can they be identified from actions, actions,and and rewards? rewards?and and (2). (2).Whether Whether they they can can beas be identified identified from from demonstration demonstration data alone? alone? demonstration data alone? modelunobservable, the system that extends thelatent standard WeWe model model thethe system system that that extends extends thethe standard standard MDP MDP toWe to include include unobservable, time-varying time-varying latent MDP to include unobservable, time-varying latent variables that affect both the transition dynamics and the reward variables variables that that affect affect both both thethe transition transition dynamics dynamics and and thethe reward reward function. function. This This model model generalizes generalizes thethe function. This model generalizes the contextual MDP bythen allowing the context contextual contextual MDP MDP byby allowing allowing thethe context context to to evolve evolve stochastically stochastically over over time. time. WeWe then formalize formalize thethe to evolve stochastically overthis time. We then formalize the data generation process under this model data data generation generation process process under under this model model using using structural structural causal causal models models (SCMs) (SCMs) [62]. [62]. Finally, Finally, wewe using structural causal models (SCMs) (Pearl, 2010). Finally, we present theoretical results that characterize the minimal present present theoretical theoretical results results that that characterize characterize thethe minimal minimal observational observational requirements requirements forfor identifying identifying observational requirements for identifying the latent variables. thethe latent latent variables. variables.

3.1 L ATENT CTime-dependent ONTEXTUAL POMDP WITH T IME -D EPENDENT C ONTEXT 3.13.1 Latent Latent Contextual Contextual POMDP POMDP with with Time-dependent Context Context Weusing model latent factors MDP using a framework, general contextual MDP framework, WeWe model model thethe latent latent factors factors using a general athe general contextual contextual MDP framework, where where thethe context context itself itself where the context itself we evolves over Formally, we define a MDP latent time-varying contextual MDP as a tuple evolves evolves over over time. time.Formally, Formally, we define define a latent a time. latent time-varying time-varying contextual contextual MDP as as a tuple a tuple MM == = (S, A,space, C,space, T , R, γ), where S is the state space, A is the action space, C is the latent con(S,(S, A,A, C,C, T ,TR, , R, ϖ),ϖ), where where SM Sis is thethe state state AA is is the the action action space, space, CC is is the the latent latent context context space, space, text space, distribution, T distribution, (st | st−1 ,R(s at−1 , tc,tta transition distribution, R(s , a[0, , c1) t1) t ) is the reward function, T (s T t(s|t s|t→1 st→1 , a,t→1 at→1 , c,t c) tis) is thethe transition transition R(s ,)ctis,t c)the ) is the the reward reward function, function, and and ϖϖ ↔t↔ [0, t, a tis and γ ∈ [0, 1) is the discount factor. The latent context c ∈ C follows a time-dependent (possibly t is is thethe discount discount factor. factor.The The latent latent context context ct c↔ C Cfollows follows a time-dependent a time-dependent (possibly (possibly stochastic) stochastic) t ↔ stochastic) process: c ∼ p(c | c ), and is unobserved during training and inference. The agent t t t−1 process: process: ct c↓ ↓ p(c p(c | c | c ), ), and and is is unobserved unobserved during during training training and and inference. inference. The The agent agent only only observes observes t t t t→1 t→1 only observes trajectories τthe = {(s . . . , c(scTfrom , afrom and infers theh latent ct from𝒯the 𝒯 𝒯 h context 𝒯 h 0, a 0 ), T )}, 𝒯 𝒯 trajectories trajectories ε ε= ={(s{(s , a , a ), . ), . . . , . (s , (s , a , a )}, )}, and and infer infer the latent latent context context the the observational observational 0 00 0 T T T T t t s1 s3 s2 3 1 2 h3 h1 h2 observational data. Thismodels, ismodels, naturally relevant to several MDP models, hidden E π(dynamic) π1including π2 data. data.This Thisis isnaturally naturallyrelevant relevant to toseveral several MDP MDP including including (dynamic) (dynamic) hidden parameter πE πhidden π Eparameter 3 π3 πE πE parameterMDPs MDPs (Doshi-Velez &and Konidaris, 2016; Perez etA 2020; Xie et a3 Bayes-adaptive a1 al., a2021), acomparison a1 A aal., 2 3 2 full a3 a1 a2 MDPs MDPs [63, [63, 64,64, 20], 20], Bayes-adaptive Bayes-adaptive MDPs [65, [65, 66,66, 19], 19], and factored factored MDPs MDPs [67]. [67]. full comparison MDPs (Martin, 1965; Duff, 2002; Zintgraf et al., 2021), and factored MDPs (Guestrin et al., 2003). πE and and analysis analysis is is given given in in Appendix Appendix B.B. c c c A full comparison and analysis is given in App. C. h3 h1 h2 Given Given trajectories trajectories generated generated under under this this model, model, wewe can can describe describe thethe data data generation generation process process using using h3 h1 anan h2 πdata π E using π E theπ3 Given trajectories generated under this model, we can describe the generation process 1 π2 to π3 SCM. SCM. Without Without thethe loss loss ofof generality, generality, wewe consider consider thethe setting setting where where anan expert expert policy policy ωω is is assumed assumed to a1 a2 to a3 a3 a1 expert a2 policy SCMs. Without the from loss of generality, wedata. consider the setting where an π is assumed generate generate thethe actions, actions, asas is is standard standard in in learning learning from demonstration demonstration data. Although Although reward reward signals signals areare π E process πE πE generate the actions, as is standard in learning from demonstration data. The data generation typically typically unobserved unobserved in in demonstration demonstration data, data, wewe include include them them in in thethe structural structural model model forfor completeness. completeness. c3 c1 c2 c3 c1 c2 can therefore be expressed as (l.h.s. Fig. 1): The The full full data data generation generation process process can can therefore therefore bebe expressed expressed as:as: (a) r

1 Latent Latent Dynamics: Dynamics: ct c= h(c h(c , ϱ,t ), ϱt ), t = t→1 t→1 State State Transitions: Transitions: st s= f (s f t→1 (st→1 , a,t→1 at→1 , c,t c, ςt t, ), ςt ),s1 t = Action Action Generation: Generation: ata= ω(s ω(s ),t ), t = t, c t ,t c a1 Reward Reward Function: Function: rt r= g(s g(s , ct ,t c, φt ,t ), φt ), t = t, a t ,ta

c1

r2

r3

s2

s3

a2

a3

c2

c3

(b) s1

s2

s3

a1

a2

a3

c1

c2

c3

s1

s2

s3

a1

a2

a3

where where ϱt ,ϱtς,t ,ςtand , and φt φdenote i.i.d. i.i.d. exogenous variables associated associated with with thethe latent latent context, context, state state t denote Figure 1:exogenous (a) SCMnoise ofnoise thevariables Latent Contextual POMDP. Gray/white nodes are observed/latent variables; dynamics, dynamics, and and reward, reward, respectively. respectively. green/red edges represent transitions driven by latents/expert policies, respectively. (b) Examples where latents influence either dynamics or rewards (affecting optimal actions). 33 where ηt , ϵt , and δt denote i.i.d. exogenous noise variables. Fig. 1(a) shows the graphical model. Fig. 1(b) illustrates examples where latent factors on dynamics (e.g., external wind in locomotion) and rewards (e.g., varying target objects in robot control) influence optimal decisions. 3

Published as a conference paper at ICLR 2026

3.2

I DENTIFIABILITY OF L ATENT FACTORS WITH M INIMAL M EASUREMENTS

To learn accurate dynamics and make reliable decisions, it is essential that the underlying latent factors influencing the environment are identifiable with observational data. We present theoretical results that characterize the minimal number of consecutive observations required for the identifiability of the latent variables, under a set of mild and natural assumptions. Assumption 1 (First-order MDP). We consider the following conditions: P (st , at , rt , ct | st−1 , at−1 , ct−1 , ω<t−1 ) = P (st , at , rt , ct | st−1 , at−1 , ct−1 ) ,

where ω<t−1 = {st−2 , . . . , s1 , at−2 , . . . , a1 , ct−2 , . . . , c1 }.

This is naturally satisfied under our setting described in Section 3.1. Assumption 2 (Distributional Variability). There exist observed state and action variables xt such that for any xt ∈ Xt , there exists a corresponding xt−1 ∈ Xt−1 and a neighborhood N r around (xt , xt−1 ) satisfying that, for all xt−2 ∈ Xt−2 , xt−1 ∈ Xt−1 , xt ∈ Xt , and xt+1 ∈ Xt+1 , the following conditional distribution operators are injective: (i) Lxt−2 |xt+1 , (ii) Lxt+1 |xt ,ct , and (iii) Lxt |xt−2 ,xt−1 , where the conditional operator L represents transformations at the distribution level, that is, how one probability distribution is pushed forward to another (Dunford & Schwartz, 1971). Assumption justification. Conceptually, the injectivity of these operator L implies that different inputs induce different output distributions, thus imposing a minimal condition on distributional variability. In RL systems, this condition is naturally satisfied in most stochastic environments where transitions produce sufficient diversity across different states and actions. The assumption also aligns with the conditions in identifiability theory, particularly in works using spectral decomposition and latent variable models (Hu & Schennach, 2008; Hu & Shum, 2012; Fu et al., 2025). We further verify this empirically using MuJoCo RL trajectories with the context instantiated as time-varying wind (App. B.5.1). Assumption 3 (Uniqueness of Spectral Decomposition). For any xt ∈ Xt and any c̄t ̸= c̃t ∈ Ct , there exists a xt−1 ∈ Xt−1 and corresponding neighborhood N r satisfying Assumption 2 such that, for some (x̄t , x̄t−1 ) ∈ N r with x̄t ̸= xt , x̄t−1 ̸= xt−1 : i. 0 < k(xt , x̄t , xt−1 , x̄t−1 , ct ) < C < ∞ for any ct ∈ Ct and some constant C; ii. k(xt , x̄t , xt−1 , x̄t−1 , c̄t ) ̸= k(xt , x̄t , xt−1 , x̄t−1 , c̃t ), where k(xt , x̄t , xt−1 , x̄t−1 , ct ) =

pxt |xt−1 ,ct (xt | xt−1 , ct )pxt |xt−1 ,ct (x̄t | x̄t−1 , ct ) . pxt |xt−1 ,ct (x̄t | xt−1 , ct )pxt |xt−1 ,ct (xt | x̄t−1 , ct )

(1)

Assumption justification. Conceptually, Assumption 3 requires that k, which captures second-order variations in transition dynamics at time t − 1 and t under the latent variable c, yields distinct values for different c’s. This requirement is typically met in RL, as varied latent dynamics or rewards often cause significant, observable shifts in behavior. Crucially, this variability is precisely what motivates the need for the identification of the latent variable ct , as it governs meaningful differences in learning underlying decision-making process. We further verify this empirically using MuJoCo RL trajectories with the context instantiated as time-varying wind (App. B.5.2). These assumptions are mild and natural. While Assumption 1 is standard in RL, it can be relaxed without violating our theory (App. B.3.4). Assumptions 2–3 are naturally satisfied in practice, as they simply formalize that latent variables influence the dynamic, motivating why we need the identification of them. Further validation and discussion are provided in App. A.4. Importantly, the more strongly the context influences the dynamics (and thus the more critical it becomes to account for c in decision-making), the more strongly these two assumptions are satisfied: the transition operator becomes more injective as required in Assumption 2, and the spectral ratio k becomes more separable across contexts as required in Assumption 3 (See empirical validation in App. B.5.3 ). Under these assumptions, we establish an identifiability theory that characterizes the conditions under which the latent factors can be recovered, and specifies the level of identifiability that can be achieved. 4

Published as a conference paper at ICLR 2026

Latent Factor Identi cation

s0

r0

a0

sT

aT

rT

Sequential Encoder

c0

r0

a0

sT

Ada-Di user-Policy

st−2 st−1 st st+1

st−2 st−1 st st+1

at−2 at−1 at at+1

at−2 at−1 at at+1

at−2 at−1 at at+1

ct−2 ct−1 ct ct+1

ct−2 ct−1 ct ct+1

ct−2 ct−1 ct ct+1

cT

Causal Di usion Model

Sequential Decoder

s0

Ada-Di user-Planning

st−2 st−1 st st+1

aT

rT

st−2 st−1 st st+1

st−2 st−1 st st+1

st−2 st−1 st st+1

at−2 at−1 at at+1

at−2 at−1 at at+1

at−2 at−1 at at+1

Di usion I/O

Mask

Inverse dynamics output

Figure 2: Overview of the Ada-Diffuser framework. The modular design consists of two main stages: latent context identification (Stage 1, Section 4.2), followed by a causal diffusion model (Stage 2, Section 4.3) that models the generative structure of the trajectories. The learned model is then used for planning or policy learning conditioned on the inferred latent context. Theorem 1 (Identifiability on Latent Factors). Under Assumptions 1-3, the posterior distribution of latent factor with consecutive observations p(ct | xt−2:t+1 ) can be identifiable up to an invertible transformation on the latents ĉt = h(ct ), where ĉt is estimated latents and h is an invertible function. ff

ff ff

fi

The proof is in App. B.2. Theorem 1 indicates that a short temporal window of observations (with future frame at t + 1) contains sufficient information to recover the posterior distribution over the true latent factors (up to an invertible transformation) in an online manner, without requiring access to the full trajectory. This form of identifiability is standard in representation learning and is sufficient for downstream tasks such as dynamics modeling, planning, and control. Any policy or dynamics model that conditions on ĉt can implicitly compose with h−1 without loss of expressiveness. We further discuss the implications of this finding in greater detail in App. A.4.

4

L ATENT-AWARE A DAPTIVE D IFFUSION P LANNER AND P OLICY

Building on Theorem 1, we introduce the Ada-Diffuser framework for learning and planning with latent identification. As illustrated in Fig. 2, Ada-Diffuser models the trajectory generation process via two modules: (1) latent factor identification block, which estimates the sequence of latent variables from the observable trajectories; and (2) causal diffusion model, which learns the causal generative process of RL trajectories and explicitly infers latent context. Guided by the theoretical findings in Theorem 1 and the generative process (Sec C), we use autoregressive denoising for temporal dependencies and a backward-refinement step over a minimal–sufficient block, designed via a tailored noise schedule and zig–zag sampling, to recover the latent posterior in an online manner. In this section, we first present a general formulation of conditional diffusion modeling with latent variables. We then describe the two modules of Ada-Diffuser in detail (Fig. 2). The complete algorithmic pseudocode of the training and inference procedures are given in App. D.1. 4.1

L ATENT-AUGMENTED D IFFUSION M ODEL FOR P LANNING AND P OLICY L EARNING

Without loss of generality, we denote the observable trajectory as τx , which may correspond to a state-action sequence τsa or a state-only sequence τs , depending on the task setting. To incorporate latent structure, we augment the observable trajectory with the estimated latent context, yielding the full trajectory representation τ = [τx , τc ], where τc denotes the inferred sequence of latent variables. We train a conditional diffusion model to generate trajectories conditioned on desired attributes y(τ ) (e.g., reward or goal specification) and the identified c. The denoising model ϵθ is trained to predicth the noise added duringi the forward diffusion process via the objective: Ldiff = 2 Eτ 0 ,y,t,ϵ ∥ϵθ (τ t , t, y(τ ), c) − ϵ∥ , where τ 0 is a clean trajectory sample, ϵ ∼ N (0, I), and 5

Published as a conference paper at ICLR 2026

√ √ the noisy trajectory at diffusion step t is constructed as: τ t = ᾱt τ 0 + 1 − ᾱt ϵ, where ᾱt denotes the cumulative product of the forward noise schedule. Here, the superscript t indexes diffusion steps, and should not be confused with the environment time step indices within the trajectory. Ada-Diffuser can flexibly adapt to generate different components of the trajectory depending on the task. In the planning setting, the model generates full trajectories τ = {xt , xt+1 , . . . , xt+Tp }, where Tp denotes the planning horizon. Here, xt may have two cases: (i) xt = {st , at }, when both states and actions are generated, (ii) xt = {st }, when only states are generated. In the latter case, we train an inverse dynamics model (IDM) (Ajay et al., 2023a) to infer the corresponding actions from state transitions. In the policy learning setting, the model generates only actions, i.e., τ = {at+1 , at+2 , . . . , at+Ta }, where Ta is the action generation horizon. While multi-step action generation methods (e.g., DP (Chi et al., 2023)) can also be viewed as a form of planning (Zhu et al., 2023), for generality, we categorize such settings under the policy framework. Ada-Diffuser-Policy accommodates both variants: multi-step action generation (Ta > 1), as in DP, and single-step decision-making (Ta = 1), as in IDQL (Hansen-Estruch et al., 2023). 4.2

S TAGE 1: O FFLINE L ATENT FACTOR I DENTIFICATION

Based on Theorem 1, we structure the latent inference process around temporal blocks, using short segments of trajectories to identify the latent context at each time step. We adopt a variational inference framework (Kingma & Welling, 2014) in which the latent variable ct is inferred block-wise. That is, the prior distribution is conditioned on the latent variable from the previous step and the in-block history, while the posterior additionally incorporates future observations. Specifically, given a trajectory block t − Tx : t + 1, where Tx is the block size, we have prior pϕ (ct | ct−1 ), and posterior qψ (ct | xt−Tx :t+1 ), where x denotes the observed variables and may correspond to {s}, {s, a}, or {s, a, r}. We then optimize the evidence lower bound (ELBO) of the observed trajectories: LELBO,t = Eqψ (ct |xt−Tx :t+1 ) [− log pθ (xt | xt−1 , ct )] + DKL (qψ (ct | xt−Tx :t+1 ) ∥ pϕ (ct | ct−1 )) . Here, the reconstruction term, − log pθ (xt | xt−1 , ct ) is instantiated based on the available observation modalities. Specifically, (i) when only states are observed, the model reconstructs st conditioned on (st−1 , ct ); and (ii) when rewards are available, the model also reconstructs rt from (st , at , ct ). The stage is learned through a sequential encoder and decoder (l.h.s., Fig. 2). 4.3

S TAGE 2: C AUSAL D IFFUSION M ODEL

We propose a causal diffusion model for learning the generative process described in Sec.3.1. By “causal,” we refer to the modeling of the true underlying data generation process, which incorporates two key desiderata: (1) the autoregressive process inherent in temporal sequential RL trajectories; and (2) the latent factor process, capturing the causal influence of the unobserved context variables ct on the observations (e.g., xt = [st , at , rt ]). Thus, unlike prior diffusion-based RL methods and latent-augmented variants (Sec. 2; see Table 13 for a comparison), our approach incorporates the following design choices. Autoregressive Denoising To model the autoregressive structure of trajectory generation, and following the recent advances in autoregressive diffusion (Chen et al., 2024; Xie et al., 2024b; Wu et al., 2023), we introduce a causal denoising schedule. Under this mechanism, each time step within a local temporal block is assigned a denoising schedule that depends both on its temporal distance from the conditioning anchor and on the inferred latent variables. This reflects the intuition that later time steps exhibit higher uncertainty. Specifically, for a trajectory of length T , we assign monotonically increasing noise levels {k1 , . . . , kT }, sampled linearly as ki = Ti K where i ∈ {1, . . . , T } and K denotes the maximum diffusion step. Given the inferred latent context ĉ0:T , the model performs autoregressive denoising over the block in T steps. The overall denoising process is defined as:   pθ x00 , . . . , x0T −1 | xk01 , . . . , xkTT−1 , ĉ0:T , (2) where xki i denotes the noisy observation at time step i, and x0i is the clean, denoised output. k −1 Specifically, the first denoising step is: pθ (x00 , xk11 , . . . , xTT−1 | xk01 , . . . , xkTT−1 , ĉ0:T ), where the first observation x0 has been fully denoised and other observations are partially denoised, followed by 6

Published as a conference paper at ICLR 2026

k

k

−2 −1 the second step: pθ (x01 , xk21 , . . . , xTT−1 | x00 , xk11 , . . . , xTT−1 , ĉ0:T ), and finally until all observations k1 0 0 0 are denoised: pθ (xT −1 | x0 , . . . , xT −2 , xT −1 , ĉ0:T ).

Denoise-and-refine Mechanism Theorem 1 indicates that both historical and future observations are required for recovering the latents. However, these future observations are not accessible during online inference, which results in a mismatch between identifiability requirements and available information. Hence, guided by this insight with preserving the causal structure of the generative process, we propose a novel denoise-and-refine mechanism that alternates between denoising the observable sequences and refining the latent estimates, and is applied consistently during both training and inference to ensure high-quality latent context recovery in an online manner. We introduce how we implement this during training and inference.

Comparison with o

Training: Given a noisy input xkt t with noise level kt , we first sample an initial latent context from the prior: ĉprior ∼ pϕ (ct | ct−1 ), and use it to denoise t (0) prior kt the observation: x̂t = ϵθ (xt , kt , ĉt ). Then we infer the latent using the posterior network, conditioned on a broader temporal window including future observations (accessible in offline data): ĉpost ∼ qψ (ct | xt−k:t+1 ), and obtain t (0)′ a refined denoised prediction: x̂t = ϵθ (xkt t , kt , ĉpost t ).

xKt−1

xKt

̂ cprior t−1

̂ cprior t

̂ cpost

̂ cprior t

xk1 t−1

t−1

x0t−1

xk2 t xk1 t

We have two reconstruction losses: one from the prior-sampled latent, ̂ cpost t (0) 0 2 Lprior = ∥x̂t − xt ∥ , and one from the posterior-sampled latent, Lpost = 0 x0t xt−1 (0)′ ∥x̂t − x0t ∥2 . To encourage the posterior latent to produce better reFigure 3: zig-zag constructions, we introduce a contrastive improvement loss: Lrel =   sampling (2 steps). softplus log Lpost − log sg Lprior + m , where sg(·) denotes stop-gradient, u softplus(u) = log(1 + e ), and m ≥ 0 is a margin hyperparameter. The final objective for this denoise-and-refine step is: Ld-r = Lpost + λprior Lprior + λrel Lrel , where λprior and λrel are weighting coefficients. Ldiff updates only θ, Lpost updates only ψ, Lprior updates only ϕ, and Lrel updates both ϕ and ψ. Inference: During inference, future observations are not available, which prevents direct use of the posterior network for latent inference. To address this, we adopt a zig-zag sampling strategy2 that combines autoregressive denoising with latent refinement. Specifically, we first sample the entire trajectory by applying the forward diffusion process with the maximum noise level K. We then perform autoregressive denoising across time. fl

fl

For each time step t, we begin by denoising xK t to an intermediate noise level k1 using ĉt sampled from the prior: ĉprior ∼ pϕ (ct | ct−1 ). We then obtain updated ĉt from the posterior latent distribution t 2 ĉpost ∼ qψ (ct | x0t−k:t−1 , xkt 1 , xkt+1 ), which is conditioned on the denoised history, the intermediate t step with noise level k1 , and the next step with noise level k2 . We then use ĉpost as the input to further t k1 0 denoise xt to xt . An illustration of the zig-zag inference process is provided in Fig. 33 . In summary, Ada-Diffuser leverages autoregressive noise scheduling to reflect temporal structure, integrates latent context identification by the denoise-and-refine mechanism, and employs zig-zag sampling for online latent inference. This framework accommodates a wide range of scenarios, including latent dynamics/rewards, learning from action-free data with latent actions, and both stateand image-based environments. All variants share the same core, with task-specific modifications to the input/output only. Details of these architectural and variations are in App. H.

5

E XPERIMENTS

We aim to answer the following questions in the evaluation: (1) Latent Identification: How well can Ada-Diffuser capture latent factors in the environment? (2) Learning with Latent Factors: How effective is Ada-Diffuser in planning and control when learning with the latent context on dynamics and reward? And can Ada-Diffuser infer latent actions from action-free demonstrations? (3) Learning with Environments w/o Explicit Latents: In environments without explicit latent factors, 2 Note on terminology: our use of “zig–zag” is purely descriptive, and there is no connection between the proposed sampling and Bai et al. (2024). 3 A larger illustration with 4 steps are given in App. Fig. 7.

7

Published as a conference paper at ICLR 2026

1.0 1.0

0.8

Value

0.8 1.0

Value Value

0.6 0.4 0.2

0.8 0.6 0.6

0.4

MSE 0.6 R2 MSE Reward0.4R2

Value

1.0

0.8(a)

MSE R2 Reward(b)

Reward

0.2

3 4 5 6 7 8 9 10 15 20 25 30 40 50100

Block Size

3 40.4 5 6 7 8 9 10 15 20 25 30 40 50100

0.2

Block Size

0.2

20402550100 30 40 50100 33 44 5 56 67 87 981091510 20 15 25 30

Block SizeSize Block

Figure 4: (a). Identification Results (i.e., Linear Probing MSE, R2 ) and normalized rewards on the Cheetah environment with time-varying wind as the latent factor, evaluated across different block sizes. (b). Results (i.e., average success rate) on planning with action-free demonstrations on Robomimic benchmark. "AF" denotes Action-free. can modeling latent processes still bring performance gains? (4) Ablation Studies: What is the impact of key design choices in the framework? 5.1

S ETTINGS

Benchmarks We consider a diverse set of benchmarks, including Mujoco-based locomotion tasks (Cheetah, Ant, Walker), a robot navigation task (Maze2D), and a robot arm control task (Franka-Kitchen) (Gupta et al., 2020), all from the D4RL benchmark suite (Fu et al., 2020). We also consider robotic manipulation tasks from RobotMimic (Mandlekar et al., 2021) and LIBERO-10 (Liu et al., 2023). A detailed description and illustration of these environments is provided in App. F. We introduce latent factors affecting both dynamics (cs ) and reward functions (cr ) in the Cheetah and Ant environments, considering two types of variations: episodic changes (E) and fine-grained, time-varying step-wise changes (S). The specific change functions for each setting are detailed in App. F.1. For evaluating latent action modeling, we follow the setup from LDP (Xie et al., 2025), using action-free, pixel-based demonstrations from the LIBERO benchmark (Liu et al., 2023). We use our framework to learn the inverse dynamics model to infer the latent actions (details are in App. G.1). In total, we evaluate on 8 environments with 23 settings. Baselines We compare Ada-Diffuser with a diverse set of baselines for fair and comprehensive evaluation. (1) Vanilla diffusion models: For planning, we consider Diffuser (Janner et al., 2022) and DD (Ajay et al., 2022). For policy learning, we include DP and IDQL (Hansen-Estruch et al., 2023). We also evaluate LDCQ (Venkatraman et al., 2024), which learns a latent skill space and optimizes a value function conditioned on both states and latent skills. (2) Latent context modeling: We include MetaDiffuser (Ni et al., 2023) that learns contextual representations from multiple environments. We also consider using LILAC (Xie et al., 2021) and DynaMITE (Liang et al., 2024a) which models nonstationarity in RL through latent context learning using belief states. For a fair comparison, we integrate their context modules into diffusion planners and policies as plug-in components (detailed analysis in App. H.1). (3) Latent action modeling: We compare with LDP (Xie et al., 2025) with action-free demonstrations for planning. In total, we compare with 9 baselines across these settings. Architecture Choices (Details are in App. D.2) For latent factor identification, we use GRU (Cho et al., 2014) embedding with MLP layers as both prior and posterior encoders to produce Gaussian distribution over latents. For decoders, we use MLP layers. For planning and policy learning, we use UNet (Ronneberger et al., 2015) or Transformers (Vaswani et al., 2017) as denoising networks and use MLPs to learn the IDM. We use VAE (Kingma & Welling, 2014) for the visual encoders. 5.2

R ESULTS AND A NALYSIS

Results on Latent Identification To verify our identification theory, we evaluate model performance under different block sizes that contain varying amounts of temporal context. We include settings where all blocks have sufficient observations, as well as a challenging case with insufficient observations (i.e., without access to future observations). To quantify the quality of the learned latent representations, we adopt linear probing and the coefficient of determination R2 as the evaluation metric. The results, together with normalized results, are shown in Fig. 4(a). Similarly, we also provide the clustering result in App. Fig. 11. The yellow region indicates settings with insufficient observations, resulting in lower identification results. The purple region corresponds to sufficient observations and yields relatively strong performance, and the green region reflects larger block sizes, 8

Published as a conference paper at ICLR 2026

Figure 5: Results on environments without explicitly designed latent factors. Complete results are provided in App. Table 8–11. Environment

Diffuser

DF

Cheetah-Wind-E (cs ) -120.4 ± 12.7 -105.8 ± 9.6 Cheetah-Wind-S (cs ) -148.5 ± 9.8 -102.0 ± 10.2 Cheetah-Vel-E (cr ) Ant-Dir-E (cr )

-102.4 ± 18.2 -85.6 ± 18.3 188.6 ± 39.2 195.4 ± 47.0

DF + DynaMITE

DF + LILAC

MetaDiffuser

Ours

-82.3 ± 8.2 -87.2 ± 10.4

-91.5 ± 7.8 -96.7 ± 9.5

-95.3 ± 7.4 -105.6 ± 14.5

-68.9± 7.6 -73.5± 8.7

-60.2 ± 10.8 266.7 ± 28.1

-67.8 ± 11.0 -62.6 ± 11.1 -45.8± 9.5 233.6 ± 31.9 229.4 ± 32.6 285.3± 24.5

Table 1: Results (5 seeds) on Ada-Diffuser-Planner with latent factors that affects dynamics and rewards. cs and cr indicate the changes on dynamics and reward, E and S represent the episodic and time-step changes. All results are averaged over 5 random seeds.

which lead to degraded results due to redundant information or inherent difficulty for optimization. Notably, the reward is positively associated with the accuracy of latent identification, validating the importance of identifying latent factors in RL trajectories. Results on Decision-making We consider three groups based on the kind of latent factors. > Group I: Latent factors on dynamics and reward: Table 1 presents the results of learning under latent factors that affect dynamics and rewards in locomotion tasks. To ensure a fair comparison, we implement autoregressive variants of DynaMITE and LILAC using the DF backbone. Results for the DP-backbone counterparts are provided in App. Tables 6, which are consistently worse than DF. Additional results, including using DP as backbones (Ada-Diffuser-policy), oracle variants and meta-learned versions of Ada-Diffuser that use ground-truth latents as input, are provided in App. Table 6–7. From the results, we observe that Ada-Diffuser consistently achieves the best performance, with a significant margin over all baselines. In particular, it outperforms Diffusion planners and policies even when those models are enhanced with latent context modules such as DynaMITE and LILAC (pink area), which are most comparable to our setting. Furthermore, Ada-Diffuser outperforms DF, showing the effectiveness of our framework. > Group II: Latent Actions: Following Xie et al. (2025), we consider learning from action-free demonstration data, where actions are treated as latent factors to be inferred. We adopt the same setup as in (Xie et al., 2025), using a pre-trained visual encoder obtained via a VAE to learn the latent space from pixel observations. We then train a latent planner and an IDM using a diffusion-based approach. Unlike prior work, our diffusion-based latent planner additionally incorporates latent factors c to model latent context. Importantly, we train only the planner using additional action-free demonstrations. Detailed training procedures are provided in App. G.1. Results on several tasks in Robomimic benchmark show that we can bring improvements on all tasks via modeling the latent process supplementary to the latent planner in (Xie et al., 2025). Here, the IDM is trained solely on expert demonstrations. Complete results are provided in App. Table 5. > Group III: Environments w/o Explicitly Designed Latents: Crucially, in this scenario, the latent variable c effectively serves as a form of Bayesian filtering over the observed trajectories, capturing the inherent stochasticity in the data (a more detailed discussion in App. D.3). Such variability commonly arises from system noise, expert action noise, or high-level unobserved factors. The results, shown in Fig. 5 (full results provided in App. Table 8–11), support this interpretation. Even in environments without explicitly designed latent contexts, incorporating latent modeling allows Ada-Diffuser to achieve performance that is comparable to or better than these baselines. By recovering the latent variables that capture stochasticity, nonstationarity, or unobserved structure in 9

Published as a conference paper at ICLR 2026

Latent Design

Orig.

w/o latents

Freeze

0.5×

Cheetah (c ) LIBERO

-73.5 93.4

-103.5 89.3

-110.4 90.2

-85.2 90.9

-77.6 89.4

-89.5 87.6

-102.4 85.0

s

Diffusion Design

Orig.

w/o refine

w/o zigzag

same NS

random NS

Cheetah (c ) LIBERO

-73.5 93.4

-82.0 83.9

-91.6 91.4

-89.7 85.2

-84.6 88.5

s

Table 2: Ablations on Cheetah-Wind-S (planner) and LIBERO (DP-policy). the offline trajectories, the model can produce rollouts that better match the underlying dynamics, even when the demonstrations are imperfect. These findings suggest that our framework can consistently capture implicit latent process in the data, improving both trajectory modeling and planning. 5.3

A BLATION S TUDIES

We conduct ablation studies to evaluate the contributions of key components in our framework. For latent factor identification, Fig. 4(a) shows the effect of different temporal block sizes, illustrating the benefit of incorporating future observations during inference. Here, we also consider ablations where (i) the entire latent identification module is removed, (ii) the latent identification network is frozen after the first 10% of training steps, and (iii) different numbers of latent updates are used. For the causal diffusion model, we examine the impact of the following design choices: (i) removing the refinement step (w/o refine); (ii) removing zig-zag sampling (w/o zig-zag); (iii) replacing the causal noise schedule with a fixed noise level across time steps (vanilla diffusion) or with random noise scaling as in DF (Chen et al., 2024) (same NS, random NS). The results in Table 2 demonstrate the effectiveness of these modules in our framework in both settings: with and without explicit latent factors. Specifically, For the latent identification ablations, we find that the latent variables play a critical role. In particular, freezing the latent module makes the model perform poorly, because the latent context follows a temporal process and must continue adapting during training. Varying the latent dimensionality within a moderate range (about 0.5×–2×) does not significantly change performance, but using overly large latent dimensions (e.g., 4×–6×) degrades results, likely due to redundant capacity and harder optimization. In terms of causal diffusion, for refinement and zig–zag, we hypothesize the gains come from reducing posterior mismatch. We therefore run a latent probing test on Cheetah with changing wind and report linear-probe MSE across variants; Ada-Diffuser with both refinement and zig–zag attains the lowest error (Table 18; Details are in App. I.2.4). Removing backward refinement yields the largest degradation (0.18 → 0.28), consistent with the role of refinement in letting future evidence within a block update the latent posterior and reduce temporal lag. Disabling zig–zag also harms accuracy (0.18 → 0.23), suggesting that alternating conditioning helps align the denoising trajectory with the latent dynamics rather than purely following the forward temporal pass. Moreover, the gap between our full model (0.18) and the oracle that has access to true futures (0.12) is small, verifying that the predicted future is already sufficiently informative for reliable latent inference in practice. Additional ablations are provided in App. I.2, including full results, comparisons of alternative noise schedules beyond linear (App. I.2.2), sweeps over temporal block length (App. I.2.3), and analyses of long-horizon planning (App. I.2.5). Notably, we show that our method introduces no significant computational overhead in terms of training runtime and inference latency (App. I.1, Table 14-15).

6

C ONCLUSIONS

We demonstrate that identifying latent factors from sequential observations is critical for effective decision-making. We provide theoretical results that establish conditions under which latent variables can be identified using small temporal blocks of observations. This insight enables a principled integration of latent identification into a diffusion-based generative framework, allowing us to capture the underlying causal process while maintaining scalability. Ada-Diffuser is broadly applicable to a variety of settings, including planning and control tasks with or without explicit latent structure, and even action-free demonstrations. Results across diverse benchmarks show substantial improvements, validating the effectiveness of our method not only in environments with designed latent factors but also in general settings where latent structure is implicit but influential. 10

Published as a conference paper at ICLR 2026

ACKNOWLEDGEMENT We would like to acknowledge the support from NSF Award No. 2229881, AI Institute for Societal Decision Making (AI-SDM), the National Institutes of Health (NIH) under Contract R01HL159805, and grants from Quris AI, Florin Court Capital, MBZUAI-WIS Joint Program, and the Al Deira Causal Education project.

R EFERENCES Anurag Ajay, Yilun Du, Abhi Gupta, Joshua Tenenbaum, Tommi Jaakkola, and Pulkit Agrawal. Is conditional generative modeling all you need for decision-making? arXiv preprint arXiv:2211.15657, 2022. Anurag Ajay, Yilun Du, Abhi Gupta, Joshua B. Tenenbaum, Tommi S. Jaakkola, and Pulkit Agrawal. Is conditional generative modeling all you need for decision making? In The Eleventh International Conference on Learning Representations, 2023a. URL https://openreview.net/forum? id=sP1fo2K9DFG. Anurag Ajay, Seungwook Han, Yilun Du, Shaung Li, Abhi Gupta, Tommi Jaakkola, Josh Tenenbaum, Leslie Kaelbling, Akash Srivastava, and Pulkit Agrawal. Compositional foundation models for hierarchical planning. arXiv preprint arXiv:2309.08587, 2023b. Lichen Bai, Shitong Shao, Zikai Zhou, Zipeng Qi, Zhiqiang Xu, Haoyi Xiong, and Zeke Xie. Zigzag diffusion sampling: Diffusion models can self-improve via self-reflection. arXiv preprint arXiv:2412.10891, 2024. Suneel Belkhale, Yuchen Cui, and Dorsa Sadigh. Data quality in imitation learning. Advances in neural information processing systems, 36:80375–80395, 2023. Andreas Blattmann, Tim Dockhorn, Sumith Kulal, Daniel Mendelevitch, Maciej Kilian, Dominik Lorenz, Yam Levi, Zion English, Vikram Voleti, Adam Letts, et al. Stable video diffusion: Scaling latent video diffusion models to large datasets. arXiv preprint arXiv:2311.15127, 2023. Gianluca Brero, Alon Eden, Darshan Chakrabarti, Matthias Gerstgrasser, Amy Greenwald, Vincent Li, and David C Parkes. Stackelberg pomdp: A reinforcement learning approach for economic design. arXiv preprint arXiv:2210.03852, 2022. Greg Brockman, Vicki Cheung, Ludwig Pettersson, Jonas Schneider, John Schulman, Jie Tang, and Wojciech Zaremba. Openai gym. arXiv preprint arXiv:1606.01540, 2016. Hongye Cao, Fan Feng, Meng Fang, Shaokang Dong, Tianpei Yang, Jing Huo, and Yang Gao. Towards empowerment gain through causal structure learning in model-based rl. arXiv preprint arXiv:2502.10077, 2025. Raymond J Carroll, Xiaohong Chen, and Yingyao Hu. Identification and estimation of nonlinear models using two samples with nonclassical measurement errors. Journal of nonparametric statistics, 22(4):379–399, 2010. Boyuan Chen, Diego Martí Monsó, Yilun Du, Max Simchowitz, Russ Tedrake, and Vincent Sitzmann. Diffusion forcing: Next-token prediction meets full-sequence diffusion. Advances in Neural Information Processing Systems, 37:24081–24125, 2024. Lili Chen, Kevin Lu, Aravind Rajeswaran, Kimin Lee, Aditya Grover, Misha Laskin, Pieter Abbeel, Aravind Srinivas, and Igor Mordatch. Decision transformer: Reinforcement learning via sequence modeling. Advances in neural information processing systems, 34:15084–15097, 2021. Yuhui Chen, Haoran Li, and Dongbin Zhao. Boosting continuous control with consistency policy. arXiv preprint arXiv:2310.06343, 2023. Zhe Chen et al. Bayesian filtering: From kalman filters to particle filters, and beyond. Statistics, 182 (1):1–69, 2003. 11

Published as a conference paper at ICLR 2026

Cheng Chi, Zhenjia Xu, Siyuan Feng, Eric Cousineau, Yilun Du, Benjamin Burchfiel, Russ Tedrake, and Shuran Song. Diffusion policy: Visuomotor policy learning via action diffusion. The International Journal of Robotics Research, pp. 02783649241273668, 2023. Kyunghyun Cho, Bart Van Merriënboer, Caglar Gulcehre, Dzmitry Bahdanau, Fethi Bougares, Holger Schwenk, and Yoshua Bengio. Learning phrase representations using rnn encoder–decoder for statistical machine translation. In Proceedings of the 2014 Conference on Empirical Methods in Natural Language Processing (EMNLP), pp. 1724–1734, 2014. Prafulla Dhariwal and Alexander Nichol. Diffusion models beat gans on image synthesis. Advances in neural information processing systems, 34:8780–8794, 2021. Zibin Dong, Yifu Yuan, Jianye HAO, Fei Ni, Yao Mu, YAN ZHENG, Yujing Hu, Tangjie Lv, Changjie Fan, and Zhipeng Hu. Aligndiff: Aligning diverse human preferences via behavior-customisable diffusion model. In The Twelfth International Conference on Learning Representations, 2024. URL https://openreview.net/forum?id=bxfKIYfHyx. Finale Doshi-Velez and George Konidaris. Hidden parameter markov decision processes: A semiparametric regression approach for discovering latent task parametrizations. In IJCAI: proceedings of the conference, volume 2016, pp. 1432, 2016. Michael O’Gordon Duff. Optimal Learning: Computational procedures for Bayes-adaptive Markov decision processes. University of Massachusetts Amherst, 2002. Nelson Dunford and Jacob T. Schwartz. Linear Operators. John Wiley & Sons, New York, 1971. Daniel E Ehrmann, Shalmali Joshi, Sebastian D Goodfellow, Mjaye L Mazwi, and Danny Eytan. Making machine learning matter to clinicians: model actionability in medical decision-making. NPJ Digital Medicine, 6(1):7, 2023. Fan Feng and Sara Magliacane. Learning dynamic attribute-factored world models for efficient multi-object reinforcement learning. Advances in Neural Information Processing Systems, 36: 19117–19144, 2023. Fan Feng, Biwei Huang, Kun Zhang, and Sara Magliacane. Factored adaptation for non-stationary reinforcement learning. In Advances in Neural Information Processing Systems (NeurIPS), 2022. Fan Feng, Phillip Lippe, and Sara Magliacane. Learning interactive world model for object-centric reinforcement learning. Advances in Neural Information Processing Systems, 38:89827–89862, 2026. Justin Fu, Aviral Kumar, Ofir Nachum, George Tucker, and Sergey Levine. D4rl: Datasets for deep data-driven reinforcement learning. arXiv preprint arXiv:2004.07219, 2020. Minghao Fu, Biwei Huang, Zijian Li, Yujia Zheng, Ignavier Ng, Guangyi Chen, Yingyao Hu, and Kun Zhang. Learning general causal structures with hidden dynamic process for climate analysis. arXiv preprint arXiv:2501.12500, 2025. Tanmay Gangwani, Joel Lehman, Qiang Liu, and Jian Peng. Learning belief representations for imitation learning in pomdps. In uncertainty in artificial intelligence, pp. 1061–1071. PMLR, 2020. Jensen Gao, Annie Xie, Ted Xiao, Chelsea Finn, and Dorsa Sadigh. Efficient data collection for robotic manipulation via compositional generalization. arXiv preprint arXiv:2403.05110, 2024a. Kaifeng Gao, Jiaxin Shi, Hanwang Zhang, Chunping Wang, and Jun Xiao. Vid-gpt: Introducing gpt-style autoregressive generation in video diffusion models. arXiv preprint arXiv:2406.10981, 2024b. Anirudh Goyal, Alex Lamb, Jordan Hoffmann, Shagun Sodhani, Sergey Levine, Yoshua Bengio, and Bernhard Schölkopf. Recurrent independent mechanisms. In International Conference on Learning Representations (ICLR), 2021. URL https://openreview.net/forum?id= mLcmdlEUxy-. 12

Published as a conference paper at ICLR 2026

Karol Gregor, George Papamakarios, Frederic Besse, Lars Buesing, and Theophane Weber. Temporal difference variational auto-encoder. arXiv preprint arXiv:1806.03107, 2018. Carlos Guestrin, Daphne Koller, Ronald Parr, and Shobha Venkataraman. Efficient solution algorithms for factored mdps. Journal of Artificial Intelligence Research, 19:399–468, 2003. Zhaohan Daniel Guo, Mohammad Gheshlaghi Azar, Bilal Piot, Bernardo A Pires, and Rémi Munos. Neural predictive belief representations. arXiv preprint arXiv:1811.06407, 2018. Abhishek Gupta, Vikash Kumar, Corey Lynch, Sergey Levine, and Karol Hausman. Relay policy learning: Solving long-horizon tasks via imitation and reinforcement learning. In Proceedings of the Conference on Robot Learning (CoRL). PMLR, 2020. Assaf Hallak, Dotan Di Castro, and Shie Mannor. Contextual markov decision processes. arXiv preprint arXiv:1502.02259, 2015. Nicklas Hansen, Xiaolong Wang, and Hao Su. Temporal difference learning for model predictive control. arXiv preprint arXiv:2203.04955, 2022. Philippe Hansen-Estruch, Ilya Kostrikov, Michael Janner, Jakub Grudzien Kuba, and Sergey Levine. Idql: Implicit q-learning as an actor-critic method with diffusion policies, 2023. Milos Hauskrecht. Value-function approximations for partially observable markov decision processes. Journal of artificial intelligence research, 13:33–94, 2000. Milos Hauskrecht and Hamish Fraser. Planning treatment of ischemic heart disease with partially observable markov decision processes. Artificial intelligence in medicine, 18(3):221–244, 2000. Haoran He, Chenjia Bai, Kang Xu, Zhuoran Yang, Weinan Zhang, Dong Wang, Bin Zhao, and Xuelong Li. Diffusion model is an effective planner and data synthesizer for multi-task reinforcement learning. Advances in neural information processing systems, 36:64896–64917, 2023. Joey Hejna, Chethan Bhateja, Yichen Jiang, Karl Pertsch, and Dorsa Sadigh. Re-mix: Optimizing data mixtures for large scale imitation learning. arXiv preprint arXiv:2408.14037, 2024. Jonathan Ho and Tim Salimans. Classifier-free diffusion guidance. arXiv preprint arXiv:2207.12598, 2022. Jonathan Ho, Ajay Jain, and Pieter Abbeel. Denoising diffusion probabilistic models. Advances in neural information processing systems, 33:6840–6851, 2020. Jonathan Ho, Tim Salimans, Alexey Gritsenko, William Chan, Mohammad Norouzi, and David J Fleet. Video diffusion models. Advances in Neural Information Processing Systems, 35:8633–8646, 2022. Yingyao Hu and Susanne M Schennach. Instrumental variable treatment of nonclassical measurement error models. Econometrica, 76(1):195–216, 2008. Yingyao Hu and Matthew Shum. Nonparametric identification of dynamic models with unobserved state variables. Journal of Econometrics, 171(1):32–44, 2012. Biwei Huang, Fan Feng, Chaochao Lu, Sara Magliacane, and Kun Zhang. Adarl: What, where, and how to adapt in transfer reinforcement learning. arXiv preprint arXiv:2107.02729, 2021. Zhiyu Huang, Chen Tang, Chen Lv, Masayoshi Tomizuka, and Wei Zhan. Learning online belief prediction for efficient pomdp planning in autonomous driving. IEEE Robotics and Automation Letters, 2024. Ahmed Hussein, Mohamed Medhat Gaber, Eyad Elyan, and Chrisina Jayne. Imitation learning: A survey of learning methods. ACM Computing Surveys (CSUR), 50(2):1–35, 2017. Maximilian Igl, Luisa Zintgraf, Tuan Anh Le, Frank Wood, and Shimon Whiteson. Deep variational reinforcement learning for pomdps. In International conference on machine learning, pp. 2117– 2126. PMLR, 2018. 13

Published as a conference paper at ICLR 2026

Michael Janner, Yilun Du, Joshua B Tenenbaum, and Sergey Levine. Planning with diffusion for flexible behavior synthesis. arXiv preprint arXiv:2205.09991, 2022. Chiyu Jiang, Andre Cornman, Cheolho Park, Benjamin Sapp, Yin Zhou, Dragomir Anguelov, et al. Motiondiffuser: Controllable multi-agent motion prediction using diffusion. In Proceedings of the IEEE/CVF conference on computer vision and pattern recognition, pp. 9644–9653, 2023. Leslie Pack Kaelbling, Michael L Littman, and Anthony R Cassandra. Planning and acting in partially observable stochastic domains. Artificial intelligence, 101(1-2):99–134, 1998. Diederik P Kingma. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980, 2014. Diederik P Kingma and Max Welling. Auto-encoding variational bayes. The International Conference on Learning Representations (ICLR), 2014. Deqian Kong, Dehong Xu, Minglu Zhao, Bo Pang, Jianwen Xie, Andrew Lizarraga, Yuhao Huang, Sirui Xie, and Ying Nian Wu. Latent plan transformer for trajectory abstraction: Planning as latent space inference. Advances in Neural Information Processing Systems, 37:123379–123401, 2024. Mikko Lauri, David Hsu, and Joni Pajarinen. Partially observable markov decision processes in robotics: A survey. IEEE Transactions on Robotics, 39(1):21–40, 2022. Shuang Li, Yihuai Gao, Dorsa Sadigh, and Shuran Song. Unified video action model. arXiv preprint arXiv:2503.00200, 2025. Wenhao Li. Efficient planning with latent diffusion. In The Twelfth International Conference on Learning Representations, 2024. URL https://openreview.net/forum?id=btpgDo4u4j. Anthony Liang, Guy Tennenholtz, Chih-wei Hsu, Yinlam Chow, Erdem Bıyık, and Craig Boutilier. Dynamite-rl: A dynamic model for improved temporal meta-reinforcement learning. arXiv preprint arXiv:2402.15957, 2024a. Zhixuan Liang, Yao Mu, Mingyu Ding, Fei Ni, Masayoshi Tomizuka, and Ping Luo. Adaptdiffuser: Diffusion models as adaptive self-evolving planners. arXiv preprint arXiv:2302.01877, 2023. Zhixuan Liang, Yao Mu, Hengbo Ma, Masayoshi Tomizuka, Mingyu Ding, and Ping Luo. Skilldiffuser: Interpretable hierarchical planning via skill abstractions in diffusion-based task execution. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, pp. 16467–16476, 2024b. Bo Liu, Yifeng Zhu, Chongkai Gao, Yihao Feng, Qiang Liu, Yuke Zhu, and Peter Stone. Libero: Benchmarking knowledge transfer for lifelong robot learning. Advances in Neural Information Processing Systems (NeurIPS), 36, 2023. Cheng Lu, Huayu Chen, Jianfei Chen, Hang Su, Chongxuan Li, and Jun Zhu. Contrastive energy prediction for exact energy-guided diffusion sampling in offline reinforcement learning. In International Conference on Machine Learning, pp. 22825–22855. PMLR, 2023. Ajay Mandlekar, Danfei Xu, Josiah Wong, Soroush Nasiriany, Chen Wang, Rohun Kulkarni, Li FeiFei, Silvio Savarese, Yuke Zhu, and Roberto Martín-Martín. What matters in learning from offline human demonstrations for robot manipulation. In Conference on Robot Learning (CoRL), 2021. James John Martin. Some Bayesian decision problems in a Markov chain. PhD thesis, Massachusetts Institute of Technology, 1965. Cuong C Nguyen, Thanh-Toan Do, and Gustavo Carneiro. Probabilistic task modelling for metalearning. In Uncertainty in Artificial Intelligence, pp. 781–791. PMLR, 2021. Fei Ni, Jianye Hao, Yao Mu, Yifu Yuan, Yan Zheng, Bin Wang, and Zhixuan Liang. Metadiffuser: Diffusion model as conditional planner for offline meta-rl. In International Conference on Machine Learning, pp. 26087–26105. PMLR, 2023. Judea Pearl. Causal inference. Causality: objectives and assessment, pp. 39–58, 2010. 14

Published as a conference paper at ICLR 2026

Christian Perez, Felipe Petroski Such, and Theofanis Karaletsos. Generalized hidden parameter mdps: Transferable model-based rl in a handful of trials. In Proceedings of the AAAI Conference on Artificial Intelligence, volume 34, pp. 5403–5411, 2020. Karl Pertsch, Youngwoon Lee, and Joseph Lim. Accelerating reinforcement learning with learned skill priors. In Conference on robot learning, pp. 188–204. PMLR, 2021. Dean A Pomerleau. Efficient training of artificial neural networks for autonomous navigation. Neural computation, 3(1):88–97, 1991. Lawrence Rabiner and Biinghwang Juang. An introduction to hidden markov models. ieee assp magazine, 3(1):4–16, 1986. Kate Rakelly, Aurick Zhou, Chelsea Finn, Sergey Levine, and Deirdre Quillen. Efficient off-policy meta-reinforcement learning via probabilistic context variables. In International conference on machine learning, pp. 5331–5340. PMLR, 2019. Allen Z. Ren, Justin Lidard, Lars Lien Ankile, Anthony Simeonov, Pulkit Agrawal, Anirudha Majumdar, Benjamin Burchfiel, Hongkai Dai, and Max Simchowitz. Diffusion policy policy optimization. In The Thirteenth International Conference on Learning Representations, 2025. URL https://openreview.net/forum?id=mEpqHvbD2h. Olaf Ronneberger, Philipp Fischer, and Thomas Brox. U-net: Convolutional networks for biomedical image segmentation. In Medical image computing and computer-assisted intervention–MICCAI 2015: 18th international conference, Munich, Germany, October 5-9, 2015, proceedings, part III 18, pp. 234–241. Springer, 2015. Sand-AI. Magi-1: Autoregressive video generation at scale, 2025. URL https://static.magi. world/static/files/MAGI_1.pdf. Andy Shih, Suneel Belkhale, Stefano Ermon, Dorsa Sadigh, and Nima Anari. Parallel sampling of diffusion models. Advances in Neural Information Processing Systems, 36:4263–4276, 2023. Shagun Sodhani, Amy Zhang, and Joelle Pineau. Multi-task reinforcement learning with contextbased representations. In International Conference on Machine Learning, pp. 9767–9779. PMLR, 2021. Richard S Sutton, Andrew G Barto, et al. Reinforcement learning: An introduction, volume 1. MIT press Cambridge, 1998. Gokul Swamy, Sanjiban Choudhury, J Bagnell, and Steven Z Wu. Sequence model imitation learning with unobserved contexts. Advances in Neural Information Processing Systems, 35:17665–17676, 2022. Ashish Vaswani, Noam Shazeer, Niki Parmar, Jakob Uszkoreit, Llion Jones, Aidan N Gomez, Łukasz Kaiser, and Illia Polosukhin. Attention is all you need. Advances in neural information processing systems, 30, 2017. Siddarth Venkatraman, Shivesh Khaitan, Ravi Tej Akella, John Dolan, Jeff Schneider, and Glen Berseth. Reasoning with latent diffusion in offline reinforcement learning. In The Twelfth International Conference on Learning Representations, 2024. URL https://openreview. net/forum?id=tGQirjzddO. Xinyue Wang and Biwei Huang. Modeling unseen environments with language-guided composable causal components in reinforcement learning. arXiv preprint arXiv:2505.08361, 2025. Zhendong Wang, Jonathan J Hunt, and Mingyuan Zhou. Diffusion policies as an expressive policy class for offline reinforcement learning. arXiv preprint arXiv:2208.06193, 2022. Tong Wu, Zhihao Fan, Xiao Liu, Hai-Tao Zheng, Yeyun Gong, Jian Jiao, Juntao Li, Jian Guo, Nan Duan, and Weizhu Chen. Ar-diffusion: Auto-regressive diffusion model for text generation. Advances in Neural Information Processing Systems, 36:39957–39974, 2023. 15

Published as a conference paper at ICLR 2026

Wei Xiao, Tsun-Hsuan Wang, Chuang Gan, Ramin Hasani, Mathias Lechner, and Daniela Rus. Safediffuser: Safe planning with diffusion probabilistic models. In The Thirteenth International Conference on Learning Representations, 2025. URL https://openreview.net/forum? id=ig2wk7kK9J. Amber Xie, Oleh Rybkin, Dorsa Sadigh, and Chelsea Finn. Latent diffusion planning for imitation learning. International Conference on Machine Learning (ICML), 2025. Annie Xie, James Harrison, and Chelsea Finn. Deep reinforcement learning amidst continual structured non-stationarity. In Marina Meila and Tong Zhang (eds.), Proceedings of the 38th International Conference on Machine Learning, volume 139 of Proceedings of Machine Learning Research, pp. 11393–11403. PMLR, 18–24 Jul 2021. Annie Xie, Lisa Lee, Ted Xiao, and Chelsea Finn. Decomposing the generalization gap in imitation learning for visual robotic manipulation. In 2024 IEEE International Conference on Robotics and Automation (ICRA), pp. 3153–3160. IEEE, 2024a. Desai Xie, Zhan Xu, Yicong Hong, Hao Tan, Difan Liu, Feng Liu, Arie Kaufman, and Yang Zhou. Progressive autoregressive video diffusion models. arXiv preprint arXiv:2410.08151, 2024b. Zilai Zeng, Ce Zhang, Shijie Wang, and Chen Sun. Goal-conditioned predictive coding for offline reinforcement learning. Advances in Neural Information Processing Systems, 36:25528–25548, 2023. Qinqing Zheng, Amy Zhang, and Aditya Grover. Online decision transformer. In international conference on machine learning, pp. 27042–27059. PMLR, 2022. Zangwei Zheng, Xiangyu Peng, Tianji Yang, Chenhui Shen, Shenggui Li, Hongxin Liu, Yukun Zhou, Tianyi Li, and Yang You. Open-sora: Democratizing efficient video production for all. arXiv preprint arXiv:2412.20404, 2024. Chuning Zhu, Raymond Yu, Siyuan Feng, Benjamin Burchfiel, Paarth Shah, and Abhishek Gupta. Unified world models: Coupling video and action diffusion for pretraining on large robotic datasets. arXiv preprint arXiv:2504.02792, 2025. Zhengbang Zhu, Hanye Zhao, Haoran He, Yichao Zhong, Shenyu Zhang, Haoquan Guo, Tingting Chen, and Weinan Zhang. Diffusion models for reinforcement learning: A survey. arXiv preprint arXiv:2311.01223, 2023. Luisa Zintgraf, Sebastian Schulze, Cong Lu, Leo Feng, Maximilian Igl, Kyriacos Shiarlis, Yarin Gal, Katja Hofmann, and Shimon Whiteson. Varibad: Variational bayes-adaptive deep rl via meta-learning. Journal of Machine Learning Research, 22(289):1–39, 2021.

16

Published as a conference paper at ICLR 2026

Appendix A Discussions and Overview

18

A.1 Broader Impact . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

18

A.2 Limitations and Future Work . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

19

A.3 Discussions on the Core Idea . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

19

A.4 Discussions on the Theoretical Assumptions and Results . . . . . . . . . . . . . .

20

A.5 Discussions on the Model Design . . . . . . . . . . . . . . . . . . . . . . . . . . .

21

A.6 Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

22

B Theory

22

B.1 Notation List . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

22

B.2 Proof of Theorem 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

23

B.3 Theory-Algorithm Alignment . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

26

B.3.1

Non-Parametric Identifiability Theory . . . . . . . . . . . . . . . . . . . .

26

B.3.2

Model Design Guidance . . . . . . . . . . . . . . . . . . . . . . . . . . .

27

B.3.3

Discussion on Assumptions . . . . . . . . . . . . . . . . . . . . . . . . .

27

B.3.4

Relaxing Assumption 1 (beyond first–order Markov) . . . . . . . . . . . .

27

B.3.5

Cases in Assumption 2. . . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

B.4 ELBO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

29

B.5 Assumption Verification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

B.5.1 About Assumption 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

B.5.2 About Assumption 3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

31

B.5.3

32

Policy Learning under Different Separability . . . . . . . . . . . . . . . .

C Summary on Different MDPs

32

C.1 Contexutal MDPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

32

C.2 Hidden-Parameter MDPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

32

C.3 Discussions and Comparisons . . . . . . . . . . . . . . . . . . . . . . . . . . . .

33

D Details on Ada-Diffuser

33

D.1 Full Algorithm and Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

33

D.2 Architecture Choices and Hyper-parameters . . . . . . . . . . . . . . . . . . . . .

34

D.2.1 Latent Factor Identification . . . . . . . . . . . . . . . . . . . . . . . . . .

34

D.2.2 Planner . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

37

D.2.3 Policy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

37

D.2.4 Hyperparameters of Contrastive Improvement Loss . . . . . . . . . . . . .

37

D.3 Connection to Bayesian Filtering . . . . . . . . . . . . . . . . . . . . . . . . . . .

38

17

Published as a conference paper at ICLR 2026

E Extended Related Works

38

E.1 Diffusion Model-based Decision-making . . . . . . . . . . . . . . . . . . . . . . .

38

E.2 Latent Belief State Learning in POMDP . . . . . . . . . . . . . . . . . . . . . . .

39

E.3 Autoregressive Diffusion Models . . . . . . . . . . . . . . . . . . . . . . . . . . .

39

E.4 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

40

F Benchmark Settings and Illustrations

40

F.1

Latent Change Factors Design . . . . . . . . . . . . . . . . . . . . . . . . . . . .

40

F.2

Overview on Other Benchmarks . . . . . . . . . . . . . . . . . . . . . . . . . . .

41

G Other Details on Ada-Diffuser

42

G.1 Latent Action Planner . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

42

G.2 Noise Scheduling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

42

H Specific Design Choices for Baselines

I

43

H.1 Details on LILAC and DynaMITE . . . . . . . . . . . . . . . . . . . . . . . . . .

43

H.2 Details on Diffusion Forcing . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

43

Ablation Analysis

43

I.1

Training/Inference Time Analysis . . . . . . . . . . . . . . . . . . . . . . . . . .

43

I.2

Ablation Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

44

I.2.1

Full Results Supplement to Table 2 . . . . . . . . . . . . . . . . . . . . .

44

I.2.2

Noise Schedule: Linear vs. Logistic vs. Sigmoid . . . . . . . . . . . . . .

44

I.2.3

Effect of Temporal Block Length on Latent Identification . . . . . . . . . .

45

I.2.4

Latent Probing: Effect of Backward Refinement and Zig–Zag . . . . . . .

45

I.2.5

On the Effect of Planning and Execution Horizons: Long-horizon Planning

46

A

D ISCUSSIONS AND OVERVIEW

In this section, we expand on the design and motivation behind Ada-Diffuser, including the rationale for modeling latent factors in decision-making, key architectural choices, and additional analysis of the experimental results presented in Section 5. We then provide an overview of the remaining contents of this appendix. A.1

B ROADER I MPACT

Our work aims to identify and leverage latent processes in generative decision-making, with applications in real-world domains such as robotics and healthcare. While these tasks may entail potential societal risks, we do not believe any specific concerns need to be highlighted here. Instead, by uncovering and modeling the underlying hidden processes, our approach promotes greater transparency in decision-making, which can ultimately lead to more reliable and trustworthy outcomes. 18

Published as a conference paper at ICLR 2026

A.2

L IMITATIONS AND F UTURE W ORK

One current limitation is that this work focuses primarily on theoretical formulation and algorithmic development. Although we evaluate on a variety of established benchmarks, real-world deployment, such as in self-driving, aerial drones, and physical robotics, remains an important direction for future work. A.3

D ISCUSSIONS ON THE C ORE I DEA

Q1: On Latent Modeling. Why is it necessary to model latent processes when we already have access to a large amount of demonstration data? In many decision-making systems, there exist unobservable variables that influence both the dynamics and the reward structure. More generally, these latent variables often evolve over time. Such scenarios are common in real-world settings, for example, in robotic control, system dynamics can be affected by external forces (e.g., wind, friction), or by varying user demands (e.g., different target positions). In these cases, learning an optimal policy requires conditioning on the latent factors, especially when they are non-stationary or when transferring to new domains. Prior work has demonstrated the importance of latent variable modeling in both reinforcement learning (RL) and imitation learning (IL) (Zintgraf et al., 2021; Liang et al., 2024a; Nguyen et al., 2021; Rakelly et al., 2019; Ni et al., 2023; Xie et al., 2021). Even with access to large demonstration datasets, it remains difficult to ensure sufficient coverage over the full space of environmental or task-specific latent factors relevant to decision-making. This limitation has been widely acknowledged in recent efforts focused on analyzing data quality and designing data collection protocols to promote generalization (Belkhale et al., 2023; Xie et al., 2024a; Hejna et al., 2024; Gao et al., 2024a). However, most of these works target fixed or task-specific latent variables. In contrast, we consider a more general setting where latent factors evolve over time and are not predefined. Our framework provides theoretical guarantees for identifying such latent variables from partial observations and seamlessly integrates this identification process into diffusion models, enabling scalability across complex decision-making tasks. Q2: On the Scenarios w/o Explicit Latents. What does latent modeling represent when no explicit latent factors are defined, and why can it still benefit decision-making? First, Latent stochasticity is always present (in real-world systems). Even in settings where all task-relevant observations are available, e.g., in locomotion tasks where full physical state information is provided, or in robotic manipulation with access to both proprioceptive and visual inputs, there may still exist underlying processes that are not directly observed. These include domain-specific factors such as external forces (e.g., wind) or dynamically changing task goals (e.g., target positions), which can be viewed as implicit latent variables. Hence, it is crucial to infer and condition on these latent factors In the extreme case where such factors are also fully observed, latent modeling can still offer significant benefits. Specifically, it can capture residual stochasticity present in the environment or demonstration data, serving to explain variability not accounted for by observable features. As shown in our formulation: st = f (st−1 , at−1 , ϵt ) , rt = g (st , at , δt ), the residual stochasticity (ϵ, δ) can be interpreted as implicit latent variables (sometimes can be time-correlated) influencing transitions and rewards. The model can then identify meaningful structure from irrelevant or noisy variations, for instance, filtering out visual background artifacts that are not predictive of dynamics or optimal actions. In this sense, the learning framework is conceptually similar to Bayesian Filtering Moreover, partial observability and attribution gaps exist even in clean data. Even in environments with consistent near-deterministic demonstrations, the agent often lacks access to the full set of latent causal factors or attributes that influence behavior. Specifically, many systems exhibit structured yet unobserved variability (e.g., task goals, preferences, intentions), and modeling this variability with latent variables improves generalization. 19

Published as a conference paper at ICLR 2026

Q3: On the Identification Theory. What does the identification theory establish, and how does it inform algorithm design? The identification theory (Theorem 1) establishes that the distribution over latent variables can be provably recovered from observable trajectories using only a small temporal window, specifically, a small temporal block of four time steps. This provides a general non-parametric theoretical guarantee that latent factors can be identified without requiring strong inductive biases or restrictive assumptions on the model class or functional form. This “four-step” result has direct implications for algorithm design. It suggests that latent identification can be effectively performed using a short temporal block, which aligns naturally with block-wise generative modeling approaches such as diffusion models. These models operate over segments or chunks of data, and our theoretical results justify using local temporal blocks to infer latent variables in a principled and scalable manner. A.4

D ISCUSSIONS ON THE T HEORETICAL A SSUMPTIONS AND R ESULTS

Q4: On the Assumptions. What do Assumption 2 (Distributional Variability) and Assumption 3 (Uniqueness of Spectral Decomposition) mean, and why are they considered mild? We expand on the intuition and practical relevance of these two assumptions below. Distributional Variability (Assumption 2) refers to the requirement that the conditional distributions p(xt−2 | xt+1 ), p(xt+1 | xt , ct ), and p(xt | xt−2 , xt−1 ) are sufficiently sensitive to variations in their input. That is, for different input pairs within a local neighborhood, the output distributions differ meaningfully, ensuring the system exhibits enough variability for identification. This assumption aligns with real-world decision-making settings (e.g., locomotion or robotic manipulation), where changes in inputs such as physical state, control policy, or reward function lead to observable changes in output distributions. Uniqueness of Spectral Decomposition (Assumption 3) builds on this by ensuring that changes in the latent variable ct induce distinct influences on the transition dynamics, specifically on the mapping from xt−1 to xt . To formalize this, we consider the operator k: p(xt | xt−1 , ct ) · p(x̄t | x̄t−1 , ct ) k(xt , x̄t , xt−1 , x̄t−1 , ct ) = , (3) p(x̄t | xt−1 , ct ) · p(xt | x̄t−1 , ct ) which separates into two multiplicative components: p(xt | xt−1 , ct ) k1 = , (4) p(xt | x̄t−1 , ct ) p(x̄t | x̄t−1 , ct ) k2 = . (5) p(x̄t | xt−1 , ct ) Here, k1 and k2 measure how changes in historical inputs affect the transition distribution at the current time step. The assumption requires that for any two distinct values of ct , the corresponding operator k is different, indicating that the latent variable has a sufficiently strong influence on the system dynamics. Since x̄ is in the neighborhood of x, this formulation effectively captures second-order changes in the transition dynamics with respect to the latent variable ct . This reflects many real-world RL systems, where even unobservable latent factors (e.g., wind speed or goal target) cause noticeable and structured changes in transition behavior over time, for instance, by considering velocity as states. In summary, these two assumptions are not only theoretically necessary for identification, but also naturally hold in many RL and control systems. They justify the need to explicitly model and identify latent variables, as such variables often induce meaningful and structured changes in both dynamics and optimal decision-making behavior. Q5: On the Identification of Posterior Distribution and up to the Invertible Function (Theorem 1). Why do we aim to identify the posterior distribution over latent variables, and what is the role of the invertible function h between the estimated and true latents? 20

Published as a conference paper at ICLR 2026

Theorem 1 establishes that the posterior distribution over latent factors given surrounding observations, p(ct | xt−2:t+1 ), is identifiable up to an invertible transformation. That is, the estimated latent ĉt satisfies ĉt = h(ct ) for some invertible function h. This form of identifiability is sufficient for downstream tasks such as dynamics modeling, planning, and control. Specifically, the learned dynamics or policy can be composed with h−1 without loss of expressiveness or utility. Since we only need to condition on the inferred latent ĉt to perform these tasks, any invertible transformation of the latent space preserves the representational capacity required for decision-making. In other words, although we may not recover the true latent variable ct exactly, the recovered representation ĉt contains the same information and can be used equivalently in practice. Therefore, identifying the posterior distribution (up to an invertible transformation) is both theoretically meaningful and practically sufficient for learning accurate dynamics models and optimal policies. A.5

D ISCUSSIONS ON THE M ODEL D ESIGN

Q6: On Different Settings (Planning and Policy). How is Ada-Diffuser applied to both planning and policy learning settings? Ada-Diffuser is designed as a unified and generic framework that accommodates different types of inputs x (e.g., states, state-action pairs) and outputs (e.g., actions, trajectories, or state sequences). This flexibility allows it to support a wide range of planning and policy learning paradigms. We summarize four representative settings below: • Planning with state-action generation: The model generates both states and actions, with latent variables influencing dynamics or rewards. This setting aligns with prior work such as Diffuser (Janner et al., 2022). • Planning with state-only generation: The model generates future states, and an inverse dynamics model is used to recover the corresponding actions. This setup follows Decision Diffuser (Ajay et al., 2023a). • Planning from action-free demonstrations: Only state sequences are available, and latent variables are assumed to capture high-level behaviors or skills. This setting extends latent diffusion planning (Xie et al., 2025). • Policy learning: The model generates actions conditioned on the current or recent history of states. This includes multi-step action generation (as in Diffusion Policy (Chi et al., 2023)) and one-step action generation (as in Implicit Diffusion Q-Learning, IDQL (Hansen et al., 2022)). In both cases, latent factors may affect the underlying dynamics or rewards. These diverse settings demonstrate the universality of our framework and highlight that uncovering latent structure is a broadly applicable and critical problem in generative decision-making. Q7: On the Latent Identification. How is Stage 1 (Latent Identification) trained, and does it introduce additional computational overhead? In Stage 1, we train the latent identification module using an offline dataset, as commonly done in offline RL and imitation learning tasks. Specifically, we employ a lightweight variational autoencoder (VAE) to optimize the ELBO defined in Section 4.2. Empirically, this stage introduces minimal computational overhead (Appendix I.1). We further provide an ablation study in Appendix I.1 showing the impact of the number of training samples on the effectiveness of the latent identification module. Q8: On the Temporal Block Design. How does this reflect Theorem 1, and why do we not use exactly four steps in practice? Our approach reflects the theoretical result in Theorem 1 by identifying latent variables using small temporal blocks in both Stage 1 and Stage 2. In Stage 1, we segment trajectories into local blocks and optimize the ELBO to learn the posterior over latent variables. In Stage 2, we apply block-wise 21

Published as a conference paper at ICLR 2026

refinement to improve the posterior estimates using both past and one-step future observations, making a more accurate identification than using the prior alone. While Theorem 1 shows that four consecutive time steps are sufficient for identifiability in principle, we do not strictly limit the block size to four in practice. Empirically, we find that using slightly larger blocks (typically between 6 and 20 steps) leads to more stable optimization and better performance. Our ablations in Appendix I.2 show that without access to future observations, identifiability degrades, aligning with the theory. We treat the "four-step" condition not as a strict architectural constraint but as a theoretical justification (sufficient condition) for using small temporal blocks. The optimal number of steps in practice may vary depending on data properties, task complexity, and model capacity. Q9: On the Refinement Step. Why is the refinement step necessary, how does it work, and does it introduce additional computational overhead? The refinement step is motivated by the identification theory, which suggests that incorporating the current and future observations (other than only using historical ones) allows the model to infer a more informative posterior over latent variables than relying on the prior alone. This posterior refinement helps the model better capture latent dynamics by leveraging richer temporal context. During training, the refinement step encourages the model to extract meaningful information from the posterior. Since Stage 1 optimizes the ELBO, the learned prior is already aligned with the posterior to some extent. This prevents the prior from collapsing into a trivial solution. The refinement step builds on this by using the pre-trained prior while further improving inference through contrastive learning between prior and posterior samples. Importantly, this procedure does not introduce significant computational overhead. As shown in Appendix I.2, the refinement uses the same denoising network with different latent inputs (c) and adds only a lightweight contrastive loss, making it efficient in practice. A.6

OVERVIEW

In this appendix, we first present the theoretical analysis in Section B, including the proof of Theorem 1 and accompanying discussion, followed by the ELBO derivation for Ada-Diffuser. In Section C, we provide an in-depth analysis of different types of MDPs and their interconnections. Section H details the full Ada-Diffuser algorithm, model architectures, and its relation to Bayesian filtering. Section E expands on related work, covering diffusion-based decision-making, latent state estimation via belief learning, and autoregressive diffusion models. Finally, Sections F, G, H, and I provide additional details on benchmarks, baseline implementations, and complete experimental results.

B

T HEORY

B.1

N OTATION L IST

We summarize the key notations used throughout the paper in Table 3, including variables for observed and latent states, temporal indices, and relevant mappings. These notations are used consistently in our theoretical analysis and algorithmic framework. Also, we formally define the operators used in the following. Definition 1 (Linear Operator (Dunford & Schwartz, 1971)). Let a and b be random variables with supports A and B, respectively. The linear operator Lb|a is defined as a mapping from a probability function pa ∈ F(A) to a probability function pb ∈ F(B), given by Z F(A) → F(B) :

pb = Lb|a ◦ pa =

pb|a (· | a) pa (a) da.

(6)

A

Intuitively, this operator characterizes the transformation of probability distributions induced by the conditional distribution pb|a . It provides a general representation of distributional change from a to b, without imposing any parametric assumptions on the underlying distributions. 22

Published as a conference paper at ICLR 2026

Index xt dx st at rt ct τ τx τc τ

Explanation

Support

[st , at ], observed trajectories including state and action at time step t dimension of observed variables state variable at time t action variable at time t reward received at time t latent context variable at time t trajectory sequence of (st , at ) observable trajectory (states or state-actions) sequence of latent contexts augmented trajectory with context

Xt ⊆ Rda +ds da + ds st ∈ St at ∈ At rt ∈ R ct ∈ Ct {(s0 , a0 ), . . . , (sT , aT )} τsa or τs {c0 , . . . , cT } [τx , τc ]

transition dynamics conditioned on ct reward function conditioned on state, action, and context expert policy used for generating demonstrations variational posterior for latent inference latent prior distribution generative model for transitions denoising network in diffusion process

T (st | st−1 , at−1 , ct ) R(st , at , ct ) π E (st , ct ) qψ (ct | xt−Tx :t+1 ) pϕ (ct | ct−1 ) pθ (xt | xt−1 , ct ) ϵθ (·)

exogenous noise in latent dynamics, state transitions, and reward distribution operator from b to a ratio of joint probabilities used in uniqueness assumption cumulative noise schedule in diffusion maximum number of diffusion steps planning and action generation horizons temporal block size for latent inference

i.i.d. samples from noise distributions defined in Dunford & Schwartz (1971) defined in Eq. 1 product of forward noise factors K∈N Tp , Ta ∈ N Tx ∈ N

Function T R πE qψ pϕ pθ ϵθ Symbol ηt , ϵt , δt La|b k(·) ᾱt K Tp , Ta Tx

Table 3: List of notations, explanations, and corresponding definitions. Definition 2 (Diagonal Operator). Let a and b be random variables with associated density functions pa and pb defined on supports A and B, respectively. For a fixed value b ∈ B, the diagonal operator Db|a is defined as a linear operator that maps a density function pa ∈ F(A) to a function in F(A) via pointwise multiplication: Db|a ◦ pa = pb|a (b | ·) · pa , (7) where Db|a = pb|a (b | ·) acts as a multiplication operator indexed by b. B.2

P ROOF OF T HEOREM 1

Proof. By the definition of data generation process (Fig. 1), the observed density is represented by: pxt+1 ,xt ,xt−1 ,xt−2 Z Z = pxt+1 ,xt ,ct ,ct−1 ,xt−1 ,xt−2 dct dct−1 Ct Ct−1 Z Z = pxt+1 |xt ,xt−1 ,xt−2 ,ct ,ct−1 pxt ,ct |xt−1 ,xt−2 ,ct−1 pct−1 ,xt−1 ,xt−2 dct dct−1 Ct Ct−1 Z Z = pxt+1 |xt ,ct pxt ,ct |xt−1 ,ct−1 pct−1 ,xt−1 ,xt−2 dct dct−1 Ct Ct−1 Z Z = pxt+1 |xt ,ct pxt |xt−1 ,ct ,ct−1 pct |xt−1 ,xt−2 ,ct−1 pxt−1 ,xt−2 ,ct−1 dct dct−1 . Ct Ct−1 Z Z = pxt+1 |xt ,ct pxt |xt−1 ,ct ,ct−1 pct ,xt−1 ,xt−2 ,ct−1 dct dct−1 . Ct

Ct−1

Then, the property of Markov process presents conditional independence, organized as follows: ! Z Z pxt+1 |xt ,ct pxt |xt−1 ,ct pct ,ct−1 ,xt−1 ,xt−2 dct−1 dct pxt+1 ,xt ,xt−1 ,xt−2 = Ct−1

Ct

Z = Ct

pxt+1 |xt ,ct pxt |xt−1 ,ct pct ,xt−1 ,xt−2 dct . 23

(8)

Published as a conference paper at ICLR 2026

Eq. 8 can be denoted in terms of operators: given values of (xt , xt−1 ) ∈ Xt × Xt−1 , Eq. 8 is Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Lct ,xt−1 ,xt−2 .

(9)

Notably, Eq. 9 is the operator representation of the observed density function in 4 measurements. Furthermore, the structure of Markov process implies the following two equalities: Z pxt+1 ,xt ,xt−1 ,xt−2 = pxt+1 |xt ,ct pxt ,ct ,xt−1 ,xt−2 dct , Ct Z pxt ,ct ,xt−1 ,xt−2 = pxt ,ct |xt−1 ,ct−1 pct−1 ,xt−1 ,xt−2 dct−1 .

(10)

Ct−1

For any fixed (xt , xt−1 ) ∈ Xt × Xt−1 , we notate Eq. 10 in terms of operators as follows: Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Lxt ,ct ,xt−1 ,xt−2 ,

Lxt ,ct ,xt−1 ,xt−2 = Lxt ,ct |xt−1 ,ct−1 Lct−1 ,xt−1 ,xt−2 .

(11)

Substituting the second line in Eq. 11 into R.H.S. of the first equation, we obtain Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Lxt ,ct |xt−1 ,ct−1 Lct−1 ,xt−1 ,xt−2

⇔ Lxt ,ct |xt−1 ,ct−1 Lct−1 ,xt−1 ,xt−2 = L−1 xt+1 |xt ,ct Lxt+1 ,xt ,xt−1 ,xt−2 .

(12)

The second line above uses Assumption 2 that L−1 xt+1 |xt ,ct is injective. Next, we show how to eliminate Lct−1 ,xt−1 ,xt−2 from the above. Consider 3 measurements {xt , xt−1 , xt−2 }, we have Z pxt ,xt−1 ,xt−2 = pxt |xt−1 ,ct−1 pct−1 ,xt−1 ,xt−2 dct−1 , (13) Ct−1

which, in operator notation (for fixed xt−1 ), is denoted as Lxt ,xt−1 ,xt−2 = Lxt |xt−1 ,ct−1 Lct−1 ,xt−1 ,xt−2 ,

Lct−1 ,xt−1 ,xt−2 = L−1 xt |xt−1 ,ct−1 Lxt ,xt−1 ,xt−2 .

(14)

The R.H.S. applies Assumption 2. Hence, substituting the above into Eq. 12, we obtain: −1 Lxt ,ct |xt−1 ,ct−1 L−1 xt |xt−1 ,ct−1 Lxt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Lxt+1 ,xt ,xt−1 ,xt−2

−1 Lxt ,ct |xt−1 ,ct−1 = L−1 xt+1 |xt ,ct Lxt+1 ,xt ,xt−1 ,xt−2 Lxt ,xt−1 ,xt−2 Lxt ,xt−1 ,ct−1 .

(15)

The second line applies Assumption 2 to post-multiply by L−1 xt ,xt−1 ,xt−2 , while in the third line, we postmultiply both sides by Lxt |xt−1 ,ct−1 . For all xt , choose a xt−1 and a neighborhood N r around (xt , xt−1 ) to satisfy Assumption 2, and pick a (x̄t , x̄t−1 ) within the neighborhood N r . Because (x̄t , x̄t−1 ) ∈ N r , we also know that (xt , x̄t−1 ), (x̄t , xt−1 ) ∈ N r . The joint distribution of of observations can be represented by Eq. 9: Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Lct ,xt−1 ,xt−2 .

(16)

The first term on the R.H.S., Lxt+1 |xt ,ct , does not depend on xt−1 , and the last term Lct ,xt−1 ,xt−2 does not depend on xt . This feature suggests that, by evaluating Eq. 9 at the four pairs of points (xt , xt−1 ), (x̄t , xt−1 ), (xt , x̄t−1 ), (x̄t , x̄t−1 ), each pair of equations will share the same operator representation in common. Specifically: Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Lct ,xt−1 ,xt−2 ,

(17)

Lxt+1 ,x̄t ,xt−1 ,xt−2 = Lxt+1 |x̄t ,ct Dx̄t |xt−1 ,ct Lct ,xt−1 ,xt−2 ,

(18)

Lxt+1 ,x̄t ,x̄t−1 ,xt−2 = Lxt+1 |x̄t ,ct Dx̄t |x̄t−1 ,ct Lct ,x̄t−1 ,xt−2 .

(20)

Lxt+1 ,xt ,x̄t−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |x̄t−1 ,ct Lct ,x̄t−1 ,xt−2 ,

24

(19)

Published as a conference paper at ICLR 2026

Assumption 2 implies that Lxt+1 |x̄t ,ct is injective. Moreover, Assumption 3 implies pxt |xt−1 ,ct (xt | xt−1 , ct ) > 0 for all ct , so that Dx̄t |xt−1 ,ct is invertible. We can then solve for Lct ,xt−1 ,xt−2 from Eq. 18 as Dx̄−1 L−1 (21) xt+1 |x̄t ,ct Lxt+1 ,x̄t ,xt−1 ,xt−2 = Lct ,xt−1 ,xt−2 . t |xt−1 ,ct Plugging this expression into Eq. 17 leads to Lxt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Dx̄−1 L−1 xt+1 |x̄t ,ct Lxt+1 ,x̄t ,xt−1 ,xt−2 . t |xt−1 ,ct

(22)

At this point, we have decomposed the observable joint operator and expressed it in terms of latentconditioned transitions, enabling spectral analysis for identifying latent structure. Lemma 1 of (Hu & Schennach, 2008) shows that, given the injectivity of Lxt−2 ,x̄t−1 ,xt ,xt+1 as in Assumption 2, we can postmultiply by L−1 xt+1 ,xt ,xt−1 ,xt−2 to obtain: −1 −1 M ≡ Lxt+1 ,xt ,xt−1 ,xt−2 L−1 xt+1 ,xt ,xt−1 ,xt−2 = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Dx̄t |xt−1 ,ct Lxt+1 |x̄t ,ct . (23)

Similarly, manipulations of Eq. 19 and 20 lead to −1 −1 N ≡ Lxt+1 ,x̄t ,xt−1 ,xt−2 L−1 xt+1 ,xt ,x̄t−1 ,xt−2 = Lxt+1 |x̄t ,ct Dx̄t |x̄t−1 ,ct Dxt |x̄t−1 ,ct Lxt+1 |xt ,ct . (24)

Assumption 2 guarantees that, for any xt , (x̄t , xt−1 , x̄t−1 ) exist so that Eq. 23 and Eq. 24 are valid operations. Finally, we postmultiply Eq. 23 by Eq. 24 to obtain:  MN = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Dx̄−1 L L × Dx̄t |x̄t−1 ,ct Dx−1 L−1 x |x̄ ,c x |x̄ ,c t+1 t t t+1 t t |x ,c xt+1 |xt ,ct t t−1 t t |x̄t−1 ,ct   = Lxt+1 |xt ,ct Dxt |xt−1 ,ct Dx̄−1 Dx̄t |x̄t−1 ,ct Dx−1 L−1 xt+1 |xt ,ct t |xt−1 ,ct t |x̄t−1 ,ct ≡ Lxt+1 |xt ,ct Dxt ,x̄t ,xt−1 ,x̄t−1 ,ct L−1 xt+1 |xt ,ct ,

(25)

where    −1 Dxt ,x̄t ,xt−1 ,x̄t−1 ,ct h (ct ) = Dxt |xt−1 ,ct Dx̄−1 D D h (ct ) x̄t |x̄t−1 ,ct xt |x̄t−1 ,ct t |xt−1 ,ct =

pxt |xt−1 ,ct (xt | xt−1 , ct )pxt |xt−1 ,ct (x̄t | x̄t−1 , ct ) h(ct ) pxt |xt−1 ,ct (x̄t | xt−1 , ct )pxt |xt−1 ,ct (xt | x̄t−1 , ct )

≡ k(xt , x̄t , xt−1 , x̄t−1 , ct )h(ct ).

(26)

This equation implies that the observed operator MN on the L.H.S. of Eq. 27 has an inherent eigenvalue–eigenfunction decomposition, with the eigenvalues corresponding to the function k(xt , x̄t , xt−1 , x̄t−1 , ct ) and the eigenfunctions corresponding to the density pxt+1 |xt ,ct (· | xt , ct ).

The decomposition in Eq. 27 is similar to the decomposition in nonparametric identification (Hu & Schennach, 2008; Carroll et al., 2010). First, Assumption 3 ensures this decomposition is unique. Second, the operator MN on the L.H.S. has the same spectrum as the diagonal operator Dxt ,x̄t ,xt−1 ,x̄t−1 ,ct . Assumption 3 guarantees that the spectrum of the diagonal operator is bounded. Since an operator is bounded by the largest element of its spectrum, Assumption 3 also implies that the operator MN is bounded, whence we can apply Theorem XV.4.3.5 from (Dunford & Schwartz, 1971) to show the uniqueness of the spectral decomposition of bounded linear operators: Lxt+1 |xt ,ct = CLxt+1 |xt ,ct P −1 .

Dxt ,x̄t ,xt−1 ,x̄t−1 ,ct = P Dxt ,x̄t ,xt−1 ,x̄t−1 ,ct P −1

(27)

where C is a scalar accounting for scaling indeterminacy and P is a permutation on the order of elements in Dx̂t |ĉt , as discussed in (Dunford & Schwartz, 1971). These forms of indeterminacy are analogous to those in eigendecomposition, which can be viewed as a finite-dimensional special case. We will show why the uniqueness of spectral decomposition is informative for identifications. First, Z px̂t+1 |x̂t ,ĉt dx̂t+1 = 1 (28) X̂t+1

must hold for every ĉt due to normalizing condition, one only solution is to set C = 1. 25

Published as a conference paper at ICLR 2026

Second, Assumption 3 implies that Eq. 27 imply that the eigenvalues k(xt , x̄t , xt−1 , x̄t−1 , ct ) are distinct for different values ct . If several ct yield identical eigenvalues, the associated eigenfunctions cannot be uniquely identified, as any linear combination of them remains valid. Therefore, for each xt , one can choose x̄t, xt − 1, x̄t−1 such that the eigenvalues differ for all ct .

Ultimately, the unorder of eigenvalues/eigenfunctions is left. The operator, Lxt+1 |xt ,ct , corresponding to the set {pxt+1 |xt ,ct (· | xt , ct )} for all xt , ct , admits a unique solution (orderibng ambiguity of eigendecomposition only changes the entry position): for all xt , ct , x̂t , ĉt

{pxt+1 |xt ,ct (· | xt , ct )} = {pxt+1 |x̂t ,ĉt (xt+1 | x̂t , ĉt )},

(29)

Due to the set is unorder, the only way to match the R.H.S. with the L.H.S. in a consistent order is to exchange the conditioning variables, that is, (1)

(1)

(2)

(2)

{pxt+1 |xt ,ct (· | xt , ct ), pxt+1 |xt ,ct (· | xt , ct ), . . .} (1)

(1)

(2)

(2)

(π(2))

, ct

= {pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ), pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ), . . .} ⇒

(π(1))

[pxt+1 |xt ,ct (· | xt

(π(1))

, ct

), pxt+1 |xt ,ct (· | xt

(π(2))

(30)

), . . .]

(π(1)) (π(1)) (π(2)) (π(2)) = [pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ), pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ), . . .]

where superscript (·) denotes the index of the conditioning variables [xt , ct ], and π is reindexing the conditioning variables. We use a relabeling map H to represent its corresponding value mapping: pxt+1 |xt ,ct (· | H(xt , ct )) = pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ),

for all xt , ct , x̂t , ĉt

(31)

By Assumption 3, different ct corresponds to different pxt+1 |xt ,ct (· | H(xt , ct )), which indicates that there is no repeated element in {pxt+1 |xt ,ct (· | H(xt , ct ))} and {pxt+1 |x̂t ,ĉt (· | x̂t , ĉt )}. Such uniqueness ensure that the relabelling map H is one-to-one. Furthermore, Assumption 3 implies that pxt+1 ,|xt ,ct (· | H(xt , ct )) corresponds a unique H(xt , ct ). The same holds for the pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ), implying that pxt+1 |xt ,ct (· | H(xt , ct )) = pxt+1 |x̂t ,ĉt (· | x̂t , ĉt ) =⇒ x̂t , ĉt = H(xt , ct )

(32)

Since the observation xt is known and suppose x̂t = xt , this relationship indeed represents an invertible transformation between ĉt and ct as ĉt = h(ct ).

(33)

which ensures that p(ct | xt−2:t+1 ) can be identifiable up to an invertible transformation on the latent variables ĉt = h(ct ) B.3

T HEORY-A LGORITHM A LIGNMENT

Here, we provide a more detailed description of our theoretical foundations, model design guidance, and algorithmic implementation. This complements the high-level summary in the main paper and offers additional context about the technical depth behind our contributions. B.3.1

N ON -PARAMETRIC I DENTIFIABILITY T HEORY

We establish this non-parametric identifiability result that gives sufficient conditions for recovering latent contexts from reinforcement learning (RL) trajectories using short temporal blocks. Each block includes a small number of future steps, which allows the model to reason about both the immediate past and the near future. Formally, we prove that under mild and broadly applicable assumptions, the latent context ct driving the generative process of the observed states and actions (xt , at ) can be recovered up to an equivalence class. This identifiability guarantee is important because it shows that latent-aware planning can be theoretically justified even when the environment contains unobserved factors or task-dependent variations. 26

Published as a conference paper at ICLR 2026

B.3.2

M ODEL D ESIGN G UIDANCE

The theoretical result directly guides the design of our causal diffusion model. To leverage Theorem 1, the diffusion model must not only generate trajectories but also recover the true latent factors. Concretely, the model must: 1. capture temporal dependencies across short blocks of states and actions; 2. jointly model observable and latent variables; 3. enforce conditions for identifiability, ensuring that the latent ct can be isolated from the observed sequence. These design requirements inform our noise schedule and the coupling of autoregressive denoising with latent refinement. This provides guidance for the algorithm design: Autoregressive Denoising. We model temporal dependencies over both observable and latent variables using an autoregressive diffusion process. At each step, xt is denoised while conditioning on partially denoised past states and inferred latent variables from a short temporal block (Section 4.1). This schedule results in a structured temporal-latent modeling process that better preserves long-range dependencies. Backward Refinement. To explicitly identify latent contexts whose posterior depends on future observations (e.g., xt+1 , guided by Theorem 1), we introduce a backward refinement step. At the second-to-last denoising stage, we refine a partial state xkt 1 using the initial estimate of ĉt sampled from the prior and xt+1 as additional evidence. The refined ĉt is then used to produce the final denoised state x0t . During training, this backward refinement is enforced to satisfy the identifiability conditions. At inference time (zig-zag sampling), we substitute xt+1 with a predicted estimate to maintain efficiency. Unification. The autoregressive denoising and backward refinement are integrated into a single noise schedule, enabling joint modeling of temporal dependencies and latent variables. Our implementation follows a four-step refinement scheme but can be extended to more steps if needed. Notably, despite the additional refinement, the method remains computationally efficient (see Appendix I.1). Further acceleration is possible via Picard iteration, which parallelizes refinement steps and reduces inference runtime by about 25%. B.3.3

D ISCUSSION ON A SSUMPTIONS

B.3.4

R ELAXING A SSUMPTION 1 ( BEYOND FIRST – ORDER M ARKOV )

We can relax the first–order Markov assumption to an n-order Markov structure with delayed/cumulative influences without altering the core identifiability argument. Suppose the generative process satisfies  p(xt+1 | x1:t , a1:t , c1:t ) = p v(xt+1 | xt:t−n+1 , at:t−n+1 , ct:t−n+1 , and that the conditioning sets across non overlapping lags exhibit block–wise conditional independence (the same separation conditions used in Theorem 1). Then there exists a finite window of observations whose statistics identify the contemporaneous block [ct , xt ] up to an invertible reparameterization.  Concrete identification statement. Let Wt = xt−2n:t+2n denote a 4n+1-length observation window.4 Assume: (i) time direction is known (so [ct , xt ] → [ct+1 v, xt+1 ] is oriented); (ii) the variability (support) conditions from Theorem 1 hold for the n-lag blocks; and (iii) block–wise independence across non–overlapping lags is satisfied. Then there exists an invertible map H such that  [ct , xt ] = H Wt , so ct (and xt ) are identifiable up to an invertible transformation from a finite window of observations. 4

Any window of length at least 4n+1 suffices; we state one concrete choice for clarity.

27

Published as a conference paper at ICLR 2026

Latent Type

Block size

Probing Acc

R2

Delayed Delayed Delayed Delayed

6 8 10 20

0.81 0.85 0.88 0.91

0.72 0.78 0.81 0.86

Cumulative Cumulative Cumulative Cumulative

6 8 10 20

0.84 0.87 0.89 0.93

0.75 0.79 0.83 0.88

Table 4: Identification under delayed and cumulative latent effects. Larger is better. Illustration for n=2. When n=2, block–wise separations allow identification of the joint variables [ct , ct+1 , xt , xt+1 ] from xt−4:t+3 (length 8+1). Knowing the temporal direction disambiguates [ct , xt ] from [ct+1 , xt+1 ]. Because xt is observed, we obtain ct = h(xt−4:t+3 ) for some invertible h, and thus the contemporaneous pair [ct , xt ] is identified. Connection to delayed/cumulative rewards. Delayed and cumulative effects fit naturally in the n-order view. For a delay ℓ,  rt+ℓ = ρ xt+ℓ , at+ℓ , ct (delayed effect), while cumulative influence over a horizon L can be written as  rt+k = ρk xt+k , at+k , ct , k = 0, . . . , L − 1, both of which are encompassed by the n-order Markov factorization above. Our cheetah variants instantiate these with, e.g., rt+ℓ = −∥vt+ℓ −ct ∥2 (delayed) and rt+k = −∥vt+k −ct ∥2 (cumulative), where vt denotes speed; the identification results remain valid. Results. We evaluate identification under delayed and cumulative latent effects in the Cheetah environment using observation windows of length 6, 8, 10, and 20. In all cases, linear probes recover the latent with high accuracy, and performance improves monotonically with longer context. For delayed effects, probing accuracy rises from 0.81 to 0.91 and R2 from 0.72 to 0.86 as block size increases from 6 to 20. For cumulative effects, probing accuracy increases from 0.84 to 0.93 and R2 from 0.75 to 0.88 over the same range. These results confirm that (i) the latent ct is behaviorally consequential in non–first-order settings and (ii) moderate temporal context suffices for accurate recovery, supporting our relaxed n-order Markov analysis. B.3.5

C ASES IN A SSUMPTION 2.

The assumption of the injectivity of a linear operator is commonly employed in the nonparametric identification (Hu & Schennach, 2008; Carroll et al., 2010; Hu & Shum, 2012). Intuitively, it means that different input distributions of a linear operator correspond to different output distributions of that operator. For a better understanding, we provide several examples in Fu et al. (2025) that describe the mapping from pa ⇒ pb , where a and b are random variables: Example 1 (Invertible). b = g(a), where g is an invertible function.

Example 2 (Additive). b = a + ϵ, where p(ϵ) must not vanish everywhere after the Fourier transform. Example 3 (Nonlinear Additive). b = g(a) + ϵ, where conditions from Examples 1-2 are required. Example 4 (Post-nonlinear). b = g1 (g2 (a) + ϵ), a post-nonlinear model with invertible nonlinear functions g1 , g2 , combining the assumptions in Examples 1-3. Example 5 (Nonlinear with Exponential Family). b = g(a, ϵ), where the joint distribution p(a, b) follows an exponential family. Example 6 (Nonparametric). b = g(a, ϵ), a general nonlinear formulation. Certain deviations from the nonlinear additive model (Example 3), e.g., polynomial perturbations, can still be tractable. 28

Published as a conference paper at ICLR 2026

B.4

ELBO

In this section, we provide analysis on the x0 -prediction Mean Squared Error (MSE) loss objectives used in the Denoise-and-Refine Mechanism of Ada-Diffuser. Our main argument establishes that minimizing the reconstruction losses Lprior and Lpost corresponds to optimizing an ELBO on the conditional log-likelihood of the clean observation x0t , given a noisy observation xkt and an inferred latent context ct . Let x0t ∼ q(x0t ) be a clean data sample from the true data distribution at sequence time step t. Let ct be the inferred latent context relevant to x0t . The forward diffusion process gradually adds Gaussian noise to x0t over K diffusion steps: √ , (1 − αk )I) q(xkt |xk−1 ) = N (xkt ; αk xk−1 t t

for k ∈ {1, ..., K}, where αk ∈ (0, 1) are predefined noise schedule parameters. This process allows sampling xkt directly from x0t : √ xkt = ᾱk x0t +

1 − ᾱk ϵ,

where ϵ ∼ N (0, I), and ᾱk =

k Y

αi .

i=1

The reverse process pθ (xk−1 |xkt , ct ) that parameterized by θ aims to denoise xkt to xk−1 conditioned t t on ct . The derivation of the ELBO for diffusion models is standard following DDPM related derivations (Ho et al., 2020; Chen et al., 2024). The conditional log-likelihood log pθ (x0t |ct ) can be lower-bounded using the ELBO: " # K X pθ (xk−1 |xkt , ct ) t K 0 log pθ (xt |ct ) ≥ Eq(x1:K log |x0t ) log pθ (xt |ct ) + t q(xk−1 |xkt , x0t ) t k=1 k−1 k Assuming pθ satisfies Markov Property (i.e., pθ (xk−1 |xkt , . . . , xK |xt , ct )), which t t , ct ) = pθ (xt is a standard structural assumption for diffusion models, the ELBO can be rewritten as:

log pθ (x0t |ct ) ≥ Eq(x1t |x0t ) [log pθ (x0t |x1t , ct )] | {z } L0

K X

Eq(xkt |x0t ) [DKL (q(xk−1 |xkt , x0t )||pθ (xk−1 |xkt , ct ))] t t {z } | k=2 Lk−1

0 K − DKL (q(xK t |xt )||pθ (xt |ct )),

|

{z

}

LK

This inequality holds with equality if and only if the model’s true posterior over the latent diffusion path, pθ (x1:K |x0t , ct ), is identical to the approximate posterior used to derive the ELBO, which is t the forward noising process q(x1:K |x0t ). This bound can also include an additive constant C(x0t , ct ) t which does not depend on the model parameters θ and is thus typically omitted when focusing on terms relevant to parameter optimization. To maximize log pθ (x0t |ct ), we aim to maximize this lower bound by optimizing L0 (i.e., maximizing this term) and each Lk−1 term (i.e., minimizing these DKL terms, as they appear with a negative sign). The term LK is often treated as a constant (or absorbed into C(x0t , ct )) if pθ (xK t |ct ) is set to a standard Gaussian N (0, I) and ᾱK ≈ 0.

We parameterize the reverse process pθ (xk−1 |xkt , ct ) as a Gaussian: t

pθ (xk−1 |xkt , ct ) = N (xk−1 ; µθ (xkt , k, ct ), σk2 I) t t

The true posterior step q(xk−1 |xkt , x0t ) is also Gaussian: t

q(xk−1 |xkt , x0t ) = N (xk−1 ; µ̃k (xkt , x0t ), σ̃k2 I) t t

where µ̃k (xkt , x0t ) =

√ √ αk (1−ᾱk−1 ) k ᾱk−1 (1−αk ) 0 ᾱk−1 xt + xt and σ̃k2 = 1− 1−ᾱk 1−ᾱk 1−ᾱk (1 − αk ) is the variance.

29

Published as a conference paper at ICLR 2026

For an x0 -prediction model, denoted as ϵθ (xkt , k, ct ) in the main paper, that aims to predict x0t from the noisy input xkt and context ct , the mean of the reverse model µθ can be expressed as: √ √ αk (1 − ᾱk−1 ) k ᾱk−1 (1 − αk ) µθ (xkt , k, ct ) = ϵθ (xkt , k, ct ) + xt 1 − ᾱk 1 − ᾱk

Choosing σk2 = σ̃k2 , the KL divergence term Lk−1 simplifies to:   1 2 k 0 k µ̃ (x , x ) − µ (x , k, c ) + Ck′ Lk−1 = Eq(xkt |x0t ) k θ t t t t 2σk2 " # √ 2 √ √ ᾱk−1 (1 − αk ) 1 2 0 0 = Ex0t ,ϵ xt − ϵθ ( ᾱk xt + 1 − ᾱk ϵ, k, ct ) + Ck′ 2σk2 1 − ᾱk

where Ck′ are constants not depending on θ. The expectation Ex0t ,ϵ denotes averaging over clean data x0t and the noise ϵ used to construct xkt . Thus, maximizing the ELBO contribution from −Lk−1 is equivalent to minimizing the following weighted MSE term: h i √ √ 2 Ex0t ,ϵ,ct w(k) x0t − ϵθ ( ᾱk x0t + 1 − ᾱk ϵ, k, ct ) (34) where w(k) = 2σ1 2

k

√

ᾱk−1 (1−αk ) 1−ᾱk

2

is a positive weighting factor.

The term L0 = Eq(x1t |x0t ) [log pθ (x0t |x1t , ct )] can also be made proportional to an MSE if pθ (x0t |x1t , ct ) is a Gaussian centered at ϵθ (x1t , 1, ct ): log pθ (x0t |x1t , ct ) = −

1 2 x0 − ϵθ (x1t , 1, ct ) + const 2σ12 t

Maximizing L0 is then equivalent to minimizing this MSE. The diffusion model ϵθ is typically trained by minimizing a simplified objective (e.g., (Ho et al., 2020)), often an unweighted or equally weighted sum of these MSE terms over uniformly sampled diffusion steps k ∈ [1, K] and data x0t : i h √ √ 2 Lsimple (θ) = Ek∼U [1,K],x0t ,ϵ,ct x0t − ϵθ ( ᾱk x0t + 1 − ᾱk ϵ, k, ct ) This simplification is justified by arguing that reweighting terms w(k) in Equation 34 can be absorbed into the network or do not significantly alter the optimal solution for expressive models, allowing w(k) to be effectively set to 1. The Denoise-and-Refine losses are: 

x0t − ϵθ (xkt i , ki , ĉprior ) t

2



x0t − ϵθ (xkt i , ki , ĉpost t )

2

Lprior = Ex0 ,ϵ,ĉprior t

Lpost = E

t

x0t ,ϵ,ĉpost t





√ √ where xkt i = ᾱki x0t + 1 − ᾱki ϵ, and ki is the specific input noise level for the observation xt determined by the causal denoising schedule ki = Ti K. These losses, Lprior and Lpost , are specific instances of the simplified MSE loss objective in equation 34 with w(ki ) ≈ 1, conditioned on the inferred contexts ĉprior and ĉpost respectively. Consequently, minimizing these MSE losses directly t t optimizes the corresponding terms in the ELBO for log pθ (x0t |ct ).

Therefore, we have proven that minimizing Lprior and Lpost as defined in the Denoise-and-Refine mechanism serves to maximize a variational lower bound on the conditional log-likelihood log pθ (x0t |ct ). The underlying diffusion model ϵθ (·, k, ·) is trained to be proficient at denoising from a range of noise levels k, as captured by objectives such as Lsimple . The specific monotonically increasing noise schedule ki used in Lprior and Lpost represents a particular instance from this range of noise levels. Thus, these objectives are theoretically grounded in the principles of variational inference for diffusion models, adapted to conditioning on the inferred latent context ct and applied at specific noise levels relevant to the autoregressive denoising process of Ada-Diffuser. 30

Published as a conference paper at ICLR 2026

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Verification of the assumptions. (a) Transition separability in Cheetah under the hyperparameter setting (m, n) = (5, 0.5). (b) Transition separability under a weak–context setting (m, n) = (0.2, 0.2), where the context barely affects the dynamics. (c) Average reward drop when planning with vs. without conditioning on c, plotted against the transition separability. (d) k–distributions for (5, 0.5), (e) k–distributions for (0.2, 0.2), (f) reward drop versus k–separability. B.5

A SSUMPTION V ERIFICATION

Here, we test whether Assumption 2 and Assumption 3 hold in practice, explain why we view them as mild, and, importantly, analyze what happens when they fail. We use the Cheetah environment, where the latent context corresponds to a time-varying wind speed fw = 5 + m sin(nt) that perturbs the agent’s dynamics. We sweep over combinations of (m, n) to address two questions: (1) whether the setting used in the paper, (m.n) = (5, 0.5), indeed satisfies these assumptions; and (2) how violations of the assumptions affect our method and the associated analyses. B.5.1

A BOUT A SSUMPTION 2

To evaluate Assumption 2, which requires that the conditional dynamics P (xt+1 | xt , ct ) be injective in the context variable ct , we perform an empirical test to determine whether different contexts induce measurably different transition dynamics. Given a fitted probabilistic dynamics model p̂(xt+1 | xt , ct ), we estimate the distribution of next states under each context c ∈ {c1 , . . . , cM } by drawing samples from the replay buffer and computing p̂(xt+1 | xt , c). For every pair of contexts (ci , cj ), we quantify the difference between their induced transition distributions using the 1-Wasserstein distance:  Inj(ci , cj ) = W1 p(xt+1 | xt , ci ), p(xt+1 | xt , cj ) . (35)

Large values of Inj(ci , cj ) indicate that distinct contexts lead to distinct transition kernels, consistent with injectivity, while values near zero suggest that different contexts produce nearly indistinguishable dynamics. When (m, n) = (5, 0.5) (Fig. 6(a)), we observe consistently non-zero Wasserstein distances across contexts, indicating that P (xt+1 | xt , ct ) is context-injective in the regime studied. In contrast, when we reduce the context variation (m, n) = (0.2, 0.2) (Fig. 6(b)), the distances are toward zero, showing the failure mode of the assumption. This shows that Assumption 2 is mild in the latent-aware decision-making. B.5.2

A BOUT A SSUMPTION 3

We provide an empirical test of the spectral ratio k to examine Assumption 3 in the RL setting. Using the dynamics model on Cheetah p̂(xt | xt−1 , ct ), we compute k(xt , x̄t , xt−1 , x̄t−1 , ct ) =

p̂(xt | xt−1 , ct ) p̂(x̄t | x̄t−1 , ct ) . p̂(x̄t | xt−1 , ct ) p̂(xt | x̄t−1 , ct ) 31

(36)

Published as a conference paper at ICLR 2026

We draw transitions (xt−1 , xt ) from the replay buffer (600 samples) and form cross-paired transitions (x̄t−1 , x̄t ) by swapping endpoints across trajectories. For each context c ∈ {c1 , . . . , cM }, this yields an empirical distribution of k(· ; c). We then quantify how well k separates contexts using the 1-Wasserstein distance between pairs of k-distributions, i.e.,  Sep(ci , cj ) = W1 p(k | ci ), p(k | cj ) . (37)

When Assumption 3 holds, k remains bounded and its distribution varies across contexts. Empirically, we observe clear multi-modal separation across contexts in the paper’s setting ((m, n) = (5, 0.5), Fig. 6(d)), whereas in regimes where the dynamics become less context-dependent, the k-distributions overlap heavily. When k is not distinguishable across contexts c ((m, n) = (0.2, 0.2), Fig 6(e)), it implies that c does not exert a noticeable effect on the transition dynamics. In this regime, explicitly modeling the context is unnecessary, since the environment effectively behaves as a single-context system. Hence, we believe Assumption 3 is mild in our main regime and also clarifies the failure mode when it is violated. B.5.3

P OLICY L EARNING UNDER D IFFERENT S EPARABILITY

When the conditional transition P (xt+1 | xt , c) is not injective in c or k is nearly the same for different c, different contexts induce nearly identical transition kernels. This means the context is not identifiable from the dynamics and does not meaningfully alter the environment; in such cases, explicitly modeling c brings little benefit for policy learning. Figures 6(c) and (f) illustrate this effect. We vary (m, n) to change the strength of the latent wind context, and compare policy performance when planning with the ground-truth context c versus ignoring c. We then plot the resulting performance gap (reward-drop ratio) against transition separability and k–separability. The gap shrinks when separability is small, indicating that when both the transition and k are weakly context-dependent, modeling c is unnecessary. Overall, these results verify that Assumption 2 and Assumption 3 are not only mild in our setting, but also clarify why modeling the latent context is important precisely in regimes where the dynamics are strongly context-dependent.

C

S UMMARY ON D IFFERENT MDP S

Our work considers a contextual POMDP setting with an evolving latent process, which naturally relates to several established MDP formulations, including contextual MDPs (Hallak et al., 2015), hidden-parameter MDPs (HiP-MDPs) (Doshi-Velez & Konidaris, 2016), and their variants. In this section, we provide formal definitions of these models and discuss their relationships and distinctions. C.1

C ONTEXUTAL MDP S

A contextual Markov decision process (CMDP) (Hallak et al., 2015) is defined by the tuple ⟨C, S, A, M⟩, where C is the context space, S is the state space, and A is the action space. The mapping M assigns to each context c ∈ C a set of MDP parameters M(c) = {Rc , T c }, where Rc and T c are the reward and transition functions associated with context c. Sodhani et al. (2021) and Liang et al. (2024a) extend the CMDP framework to settings in which the context variable c evolves according to its own Markovian dynamics p(ct+1 | ct ), closely aligning with our formulation of a latent process evolving over time. C.2

H IDDEN -PARAMETER MDP S

Hidden-Parameter MDPs (HiP-MDPs) (Doshi-Velez & Konidaris, 2016) are defined by the tuple M = ⟨S, A, Θ, T , R, γ, PΘ ⟩, where S is the state space, A is the action space, and Θ is the space of task-specific latent parameters. For each θ ∈ Θ, the transition and reward functions are given by Tθ : S × A → P(S) and Rθ : S × A → R, respectively. The parameter θ is sampled from a prior distribution PΘ at the beginning of an episode and remains fixed during the episode. The discount factor is denoted by γ ∈ [0, 1). This framework defines a family of MDPs indexed by the latent parameter θ, with each θ inducing a different set of dynamics and reward functions. It can be seen as a special case of a contextual MDP where the context is latent and fixed per episode.Xie et al. (2021) 32

Published as a conference paper at ICLR 2026

further generalize this framework by allowing the task parameter θ to evolve dynamically across episodes, rather than being fixed. Bayes-Adaptive MDPs (BAMDPs) are closely related to both HiP-MDPs and contextual MDPs (CMDPs). In BAMDPs, the agent maintains a posterior distribution over MDPs based on its interaction history. Specifically, it maintains a belief bt (R, T ) = p(R, T | τ:t ), where τ:t = {s0 , a0 , r0 , . . . , st } denotes the trajectory observed up to time t. This belief captures the agent’s uncertainty about the underlying transition and reward functions. The transition and reward functions can then be defined in expectation over this posterior, effectively conditioning decision-making on the belief bt . When the environment is driven by hidden contextual variables or latent task parameters, such as in CMDPs or HiP-MDPs—this belief can be interpreted as a distribution over these latent variables. In this view, BAMDPs provide a non-parametric framework for reasoning over hidden structure, while approaches like ours explicitly model such latent variables and infer their posterior distributions using amortized inference. Both aim to enable adaptive planning and learning under uncertainty, but differ in how latent structure is represented and inferred. C.3

D ISCUSSIONS AND C OMPARISONS

The key distinction between contextual MDPs and hidden-parameter MDPs lies in how the latent factors are represented: contextual MDPs explicitly treat them as latent variables, while HiP-MDPs model them implicitly as parameters governing the transition and reward functions. In our work, we adopt the contextual MDP perspective, where the latent process is modeled as a random variable that evolves over time. However, our identification theory, focused on recovering the posterior distribution over latent variables, also applies to the HiP-MDP setting. Once the posterior over the hidden parameters is identified, the corresponding transition and reward functions can be recovered as well. Additionally, our framework, which models a factorization over observed states and latent variables, is conceptually related to factored MDPs (Guestrin et al., 2003). In a factored MDP, the state space S is represented as a set of variables S = {s(1) , s(2) , . . . , s(n) }, and the transition and reward functions are decomposed over these factors: n m     Y X (i) (j) T (s′ | s, a) = Ti s′(i) | PaT (s, a) , R(s, a) = Rj PaR (s, a) , i=1

(i)

j=1

(j)

where PaT and PaR denote the parent variables (i.e., dependencies) for each transition and reward component, respectively. our framework, while not relying on an explicit graphical structure, shares conceptual similarities with factored MDPs (Guestrin et al., 2003) through its coarse-grained factorization over observed states and latent variables. Specifically, we distinguish between latent variables that affect the transition dynamics and those that affect the reward function. Formally, we express the generative process as: T (st+1 | st , at , cst ), R(rt | st , at , crt ), s r where ct and ct are distinct (or potentially overlapping) latent factors that influence transitions and rewards, respectively. This separation enables flexible modeling of partially observable environments where different unobserved processes govern the dynamics and task objectives.

D

D ETAILS ON A D A -D I F F U S E R

D.1

F ULL A LGORITHM AND R ESULTS

Our framework consists of two stages: latent factor identification and diffusion-based planning or policy learning. Below, we provide the algorithmic pseudocode for both stages. Specifically, Algorithm 1 describes Stage 1: latent factor identification, while Algorithms 2 and 3 correspond to Ada-Diffuser-Planner and Ada-Diffuser-Policy, respectively. For clarity, we omit the detailed step-by-step procedures for denoise-and-refine and zig-zag sampling (Lines 7–8, 11, and 19–22 in Algorithm 2; Lines 6–7 and 13 in Algorithm 3), as these are fully described in Section 4.3. For Ada-Diffuser-Policy, we show a Diffusion Policy (DP)-based 33

Published as a conference paper at ICLR 2026

other autoregressive diffusion model x x x ̂ cprior t−1

K t−1

xk1 t−1

Xk4 t−1

X0t−1

Xk1 t

Xk1 t+1

k3 Xk3 t+2 Xt+2

X0t

X0t+1

X0t+2

n process by

̂ cpost t−1

x0t−1 x0t−1

̂ cprior t

K t

xk2 t

̂ cprior t+1

K t+1

xk3 t+1

xKt+2

̂ cprior t+2

xk4 t+2

̂ cprior t

̂ cprior t+1

̂ cprior t+2

̂ cpost t

̂ cprior t+1

̂ cprior t+2

̂ cpost t+1

̂ cprior t+2

xk1 t x0t

x0t−1

x0t

x0t−1

x0t

xk2 t+1 xk1 t+1 x0t+1 x0t+1

xk3 t+2

xk2 t+2 xk1 t+2

̂ cpost t+2

x0t+2

Figure 7: An illustration of the zig-zag sampling process with a block of 4 time steps. ↓ and | indicate denoising and identity mapping, respectively. algorithm, which provides a general framework for multi-step action generation. In the IDQL-based variant, both the action execution horizon and observation horizon are set to 1, corresponding to single-step policy inference conditioned only on the current observation. Algorithm 1: Latent Factor Identification. 1: Input: offline dataset D 2: Randomly initialize decoder pθ (st+1 , rt | s, a, c), encoder qψ (ct | st−Tx :t+1 , at−Tx :t+1 , rt−Tx :t+1 ) and prior network pϕ (ct | ct−1 ), 3: while not done do 4: Sample batches of trajectories from D 5: Compute ELBO and update θ, ψ, ϕ 6: end while Additionally, we provide the full results for all experiments: Table 5 reports results for the action-free setting; Tables 6 and 7 present results for environments with latent factors affecting dynamics and rewards; and Tables 8, 9, 10, and 11 summarize results for environments without explicitly modeled latent factors. D.2

A RCHITECTURE C HOICES AND H YPER - PARAMETERS

We detail the architectural design choices and hyperparameter settings used for model components, loss functions, and training procedures across all Ada-Diffuser variants under different environments and benchmarks. D.2.1

L ATENT FACTOR I DENTIFICATION

Architectures We use a variational autoencoder (VAE) (Kingma & Welling, 2014) to optimize the evidence lower bound (ELBO). The same architectural design is used across all variants of Ada-Diffuser and all benchmark settings. For the encoder, we first embed states, actions, and rewards using separate MLPs with ReLU activations. The resulting embeddings are concatenated and passed through a two-layer MLP (each 34

Published as a conference paper at ICLR 2026

Algorithm 2: Ada-Diffuser-Planner. 1: Input: Env, offline dataset D, pre-trained encoder qψ and prior network pϕ observation horizon To , planning horizon Tp , action execution horizon Ta , condition y // Training 2: Initialize noise predictor ϵθ , inverse dynamics model fϕ 3: while not done do 4: Sample xt−To :t+Tp from D post 5: Sample ĉprior t:t+Tp and ĉt:t+Tp −2 from pϕ and qψ 6: if using inverse dynamics model then post 7: Train Causal Diffusion Model (noise predictor ϵθ ) with xt−To :t , ĉprior t−To :t , and ĉt−To :t and other conditions y, target outputs are st+1:t+Tp 8: Train encoder qψ with the contrastive improvement loss Lcontrast 9: Train Inverse Dynamics Model fϕ to generate actions at+1:t+Tp 10: else post 11: Train Causal Diffusion Model (noise predictor ϵθ ) with xt−To :t , ĉprior t−To :t , and ĉt−To :t and other conditions y, target outputs are {st+1:t+Tp , at+1:t+Tp } 12: Train encoder qψ with the contrastive improvement loss Lcontrast 13: end if 14: end while // Execution 15: Initialize environment: s0 ∼ Env.reset(), set t ← 0 16: while not done do // Observe and infer latent factors 17: Observe recent trajectory xt−To :t 18: Sample latent variables ĉprior t:t+Tp from pϕ // Generate candidate trajectory 19: if using inverse dynamics model then Generate future states (zig-zag sampling) ŝt+1:t+Tp conditioned on xt−To :t , ĉprior 20: t:t+Tp , and y via learned noise predictor ϵθ 21: Generate actions ât+1:t+Tp ← fϕ (ŝt+1:t+Tp , ŝt:t+Tp −1 ) 22: else 23: Generate future trajectory {ŝt+1:t+Tp , ât+1:t+Tp } conditioned on xt−To :t , ĉprior t:t+Tp , and y via learned noise predictor ϵθ 24: end if // Execute action(s) in environment 25: for each step i = 1 to Ta do 26: Execute ât+i in Env, observe st+i+1 , rt+i 27: Append (st+i , ât+i , rt+i ) to trajectory buffer 28: end for 29: Update t ← t + Ta 30: end while layer of size 64) followed by a GRU. The GRU output is used to parameterize a Gaussian distribution from which the latent variables are sampled. The state and reward decoders are implemented as separate MLPs, each consisting of two fully connected layers of size 64 with ReLU activations. For the prior network, we use the output of the previous step’s latent distribution embedding (shared GRU) and feed it into a two-layer MLP (each layer of size 32) to predict the parameters of the prior distribution. For the dimensionality of latents, we choose 20 for Cheetah, Walker, Ant, Maze; 64 for Robomimic, Kitchen, Libero. Loss Function

At each time step t, we optimize the following losses:

LELBO,t = Eqψ (ct |xt−Tx :t+1 ) [− log pθ (xt | xt−1 , ct )] + DKL (qψ (ct | xt−Tx :t+1 ) ∥ pϕ (ct | ct−1 )) . {z } | {z } | KL regularization

Reconstruction loss

35

Published as a conference paper at ICLR 2026

Algorithm 3: Ada-Diffuser-Policy (DP-based) 1: Input: Env, offline dataset D, pre-trained encoder qψ and prior network pϕ observation horizon To , action generation horizon Tp , action execution horizon Ta , condition y // Training 2: Initialize noise predictor ϵθ 3: while not done do 4: Sample xt−To :t+Tp from D post 5: Sample latent variables ĉprior t:t+Tp ∼ pϕ , ĉt:t+Tp −2 ∼ qψ 6: Train causal diffusion model (noise predictor ϵθ ) to generate actions at+1:t+Tp , conditioned post on xt−To :t , ĉprior t:t+Tp , ĉt:t+Tp −2 , and y 7: Train encoder qψ with the contrastive improvement loss Lcontrast 8: end while // Execution 9: Initialize environment: s0 ∼ Env.reset(), set t ← 0 10: while not done do // Observe and infer latent factors 11: Observe recent trajectory xt−To :t 12: Sample latent variables ĉprior t:t+Tp ∼ pϕ // Generate actions using causal diffusion model 13: Generate actions (zig-zag sampling) ât+1:t+Tp conditioned on xt−To :t , ĉt:t+Tp , and y via learned noise predictor ϵθ // Execute action(s) in environment 14: for each step i = 1 to Ta do 15: Execute ât+i in Env, observe st+i+1 , rt+i 16: Append (st+i , ât+i , rt+i ) to trajectory buffer 17: end for 18: Update t ← t + Ta 19: end while Environment

LDP (AF)

Ours (AF)

LDP (AF, SubOpt)

Ours (AF, SubOpt)

Lift Can Square

0.67± 0.01 0.78± 0.04 0.47± 0.03

0.78± 0.05 0.85± 0.07 0.54± 0.05

1.00± 0.00 0.98± 0.00 0.83± 0.01

0.98± 0.02 0.98± 0.02 0.89± 0.03

Table 5: Results (success rate) on action-free demonstrations. Here, AF and SubOpt indicate using Action-free and suboptimal demonstrations on Robomimic tasks, respectively (following the settings in LDP (Xie et al., 2025)). Here, xt may include different components depending on the setting (e.g., xt = {st , at } or xt = st ), and ct denotes the latent context variable inferred from a temporal block of observations. The first term encourages accurate reconstruction of the current observation xt conditioned on its immediate past and the latent ct , while the second term regularizes the posterior to remain close to the learned prior pϕ (ct | ct−1 ).

We implement the ELBO loss as a weighted combination of the reconstruction loss and the KL divergence: LELBO =

T −2 X t=1

  ∥x̂t − xt ∥22 + λKL · DKL (qψ (ct | xt−Tx :t+1 ) ∥ pϕ (ct | ct−1 )) ,

where x̂t is the model’s reconstruction of the observation xt , and λKL is weighting coefficient. The reconstruction is computed using mean squared error (MSE), and the KL divergence is computed in closed form for Gaussian posteriors and priors. The hyperparameter λKL is set to be 0.01 and the learning rate is set to be 3e − 4. 36

Published as a conference paper at ICLR 2026

D.2.2

P LANNER

For the planner, we consider two scenarios: (i) generating both states and actions, and (ii) generating states only. For the former, we build upon the Diffuser framework (Janner et al., 2022), which directly models full trajectories. For the latter, we adopt the Decision Diffuser (DD) framework (Ajay et al., 2022), where the model generates future states and uses an inverse dynamics model to recover the corresponding actions via inverse dynamics model. For type (i) (full state-action trajectory generation), we apply our method to the Cheetah and Ant environments. For the noise predictor, we use a 1D U-Net (Ronneberger et al., 2015) with a kernel size of 5, channel multipliers set to (1, 2, 2, 2), and a base channel width of 32. The model is trained using the Adam optimizer (Kingma, 2014) with a learning rate of 3 × 10−4 , a batch size of 64, and for 1 million training steps. We adopt classifier guidance (CG) (Ho et al., 2020) with gradient guidance on computed return, with a guidance scale ω = 1.5. The observation horizon is set to 10 for both environments. The planning horizon Tp is set to 16 for Cheetah and 32 for Ant, with an action execution horizon of 1. These hyperparameters are kept consistent across baselines, including Diffuser, DF, MetaDiffuser, and Diffuser combined with LILAC and DynaMITE for the Cheetah and Ant experiments (those in Table 1 and Appendix Table 6). For other components (e.g., VAE) in LDCQ, we employ all the hyperparameters in their original implementation (Venkatraman et al., 2024). For type (ii) (state-only generation with inverse dynamics), we use a Transformer-based noise predictor with a hidden dimension of 256 and a head dimension of 32. The architecture includes 2 DiT blocks for Walker, Kitchen, and Maze2D, and 8 DiT blocks for LIBERO.The model is trained using the Adam optimizer (Kingma, 2014) with a learning rate of 3 × 10−4 , a batch size of 128, and for 1 million training steps. The number of diffusion timesteps is 500. The observation horizon is set to 4 for Kitchen, 2 for LIBERO, and 10 for the other environments. The planning horizon Tp is set to 16 for Kitchen, 10 for LIBERO, and 32 for the others. The action execution horizon is 8 for both Kitchen and LIBERO, and 10 for the remaining environments. For the inverse dynamics model, we use an MLP-based diffusion model consisting of a 3-layer MLP with 128 hidden units, preceded by a 2-layer embedding module with 64 hidden units. This model is trained for 1 million gradient steps. For both cases, we set the coefficient of the contrastive improvement loss Lcontrast = max{0, Lprior − Lpost } to 0.1. The key hyper-parameters are summarized in Table 12. D.2.3

P OLICY

For the DP-based policy, we adopt the same architecture as the planner described earlier for Cheetah, Maze2D, Kitchen, Ant, and Walker. For LIBERO, we use a Transformer-based noise predictor with a decoder architecture comprising 12 layers, 12 attention heads, and a hidden embedding dimension of 768. Following DP (Chi et al., 2023), we apply dropout with a rate of 0.1 to both the input embeddings and attention weights. The number of diffusion timesteps is 500. When conditioning is used, we incorporate a Transformer encoder with 4 layers to encode the condition tokens, which include a sinusoidal timestep embedding and projected observed trajectory tokens (all mapped to the same embedding dimension). In this encoder-decoder setup, causal masking is applied to ensure autoregressive generation. In the unconditioned case, we prepend the sinusoidal timestep embedding to the input sequence and use a BERT-style encoder-only Transformer. All environments (Cheetah, Ant, Kitchen, Maze2D, Walker, and LIBERO) are trained using the AdamW optimizer with a learning rate of 10−4 , weight decay 10−3 , β1 = 0.9, and β2 = 0.95. Layer normalization is applied before each Transformer block for stability. The observation, planning, and action horizons follow the same settings used for the planner in each environment. For the IDQL-based policy, we align all hyperparameters for Cheetah and Ant with the original IDQL implementation, using an observation, planning, and action horizon of 1. Hence, in IDQL-based ones, we do not consider autoregressive modeling. Similarly, for both cases, we use consider the coefficient before the contrastive improvement loss as 0.1. D.2.4

H YPERPARAMETERS OF C ONTRASTIVE I MPROVEMENT L OSS

We set λprior , λrel to be fixed as 0.1 across all settings. m is set to be the 0.05 × Lprior during the beginning of each epoch. 37

Published as a conference paper at ICLR 2026

D.3

C ONNECTION TO BAYESIAN F ILTERING

In the absence of explicitly designed latent variables, our model can be interpreted as a form of Bayesian filtering (Chen et al., 2003). Under a general formulation of the hidden Markov model (HMM) (Rabiner & Juang, 1986) with an additional latent dependency on observation (c → x), the latent process over c captures the underlying stochasticity present in the demonstration data, which arises from both the environment dynamics and the behavior policy. In this view, the latent variable acts as a compact and expressive representation that summarizes the uncertainty in past observations, thereby improving the prediction of future observations. This, in turn, facilitates more robust policy learning and planning in the general settings.

E

E XTENDED R ELATED W ORKS

E.1

D IFFUSION M ODEL - BASED D ECISION - MAKING

Recent advances use diffusion models as the planner and policy for both reinforcement learning (RL) and imitation learning (IL). RL agent aims to learn a policy that maximizes cumulative rewards through interaction with an environment (Sutton et al., 1998). The agent observes a sequence of transitions (st , at , rt , st+1 ), where st ∈ S denotes the state, at ∈ A the action, rt ∈ R the received reward, and st+1 the next state. P∞ The goal is to learn a policy π(a | s) that maximizes the expected return: π ∗ = arg maxπ Eπ [ t=0 γ t rt ] , where γ ∈ [0, 1) is the discount factor. In contrast, IL (Hussein et al., 2017) focuses on learning policies from expert demonstrations, often without access to the reward signal. A common approach is behavior cloning (BC) (Pomerleau, 1991), which formulates IL as a supervised learning problem by maximizing the likelihood of expert actions given observed states, i.e., learning a policy π(a | s) that closely imitates the expert policy πe (a | s).

Diffusion Planner Diffusion-based planning methods are commonly used to approximate the sequence of future states and actions from a given current state. By leveraging the conditional generation capabilities of diffusion models—such as guidance techniques (Dhariwal & Nichol, 2021; Ho & Salimans, 2022)—these methods can generate plans (i.e., state trajectories) that satisfy desired properties, such as maximizing expected rewards. Taking Denoising Diffusion Probabilistic Models (DDPM (Ho et al., 2020))-based approaches as an example, these methods learn a generative model over expert trajectories τ = {(s0 , a0 ), . . . , (sT , aT )} by modeling a forward-noising process: √ q(xk | xk−1 ) = N (xk ; αk xk−1 , (1 − αk )I), and a parameterized denoising model pθ (xk−1 | xk ) to reverse the process. Here, k denotes the diffusion step, x0 is a clean sub-sequence sampled from the expert trajectory τ , and αk controls the variance schedule at step k. During inference, trajectories are generated by starting from Gaussian noise and iteratively denoising through the learned reverse process. This generation can be optionally conditioned on the initial state or other guidance signals y, such as rewards, goals, or other constraints: τ̂ ∼ pθ (τ | s0 , y).

These methods generally fall into two main categories: (1) learning a joint distribution over stateaction trajectories, as in Diffuser (Janner et al., 2022), or (2) learning only state trajectories via diffusion and using an inverse dynamics model to recover actions, as in Decision Diffuser (DD) (Ajay et al., 2022). Beyond these, several variants extend diffusion-based planning in different directions. For example, Latent Diffuser (Li, 2024) plans in a high-level latent skill space to improve generalization and LDP (Xie et al., 2025) plans with high-level latent actions directly from high-dimensional action-free demonstrations. Other approaches incorporate multi-task context to enhance adaptation and performance in unseen tasks, including MetaDiffuser (Ni et al., 2023), AdaptDiffuser (Liang et al., 2023), and MTDiff-p (He et al., 2023). In addition, recent efforts have explored various extensions of diffusion planning, such as ensuring safety during generation (Xiao et al., 2025), handling multi-agent scenarios (Jiang et al., 2023; Ajay et al., 2023b), learning skills (Liang et al., 2024b), and application in RL from human feedback (RLHF) (Dong et al., 2024). Diffusion Policy In contrast to diffusion-based planners, Diffusion Policy methods directly parameterize the policy πθ (a | s) using diffusion models. For example, Diffusion Policy (Chi et al., 2023) uses a diffusion model to generate actions with expressive, multimodal distributions. DPPO (Ren et al., 2025) extends this idea by modeling a two-layer MDP structure, where the inner MDP represents the denoising process and the outer MDP corresponds to the environment. This framework enables fine-tuning of diffusion-based policies in RL settings. Another line of work integrates diffusion 38

Published as a conference paper at ICLR 2026

models with model-free methods for offline RL by using diffusion models as to model the action distributions (Wang et al., 2022; Hansen-Estruch et al., 2023; Chen et al., 2023; Lu et al., 2023). Recent explorations have also aimed to unify diffusion-based planning and policy learning within a single framework. For example, the Unified Video Action model (UVA) (Li et al., 2025) and Unified World Models (UWM) (Zhu et al., 2025) leverage diffusion models to jointly model planning and action generation, demonstrating scalability on large-scale robotic tasks with pre-training. In a similar spirit, Ada-Diffuser provides a general framework that can be integrated into both diffusion planners and diffusion-based policies. However, Ada-Diffuser differs in its explicit modeling of latent factors that influence the data generation process. By incorporating latent identification directly into the diffusion process, Ada-Diffuser enables more structured, context-aware decision-making in partially observable and dynamically changing environments. E.2

L ATENT B ELIEF S TATE L EARNING IN POMDP

In partially observable Markov decision processes (POMDPs), single-step observations are typically insufficient for making optimal decisions. A common strategy to overcome this limitation involves encoding an agent’s history, encoding past observations and actions into a belief state that captures a distribution over latent environmental states. Although such belief representations can, in theory, support optimal policy derivation (Kaelbling et al., 1998; Hauskrecht, 2000; Gangwani et al., 2020), their exact computation depends on full knowledge of the transition and observation models. This requirement quickly becomes intractable in high-dimensional settings. To address this, recent work has focused on learning approximate belief representations directly from data. Notable approaches include those using recurrent neural networks (Guo et al., 2018) and variational inference methods (Igl et al., 2018; Gregor et al., 2018), which enable agents to encode temporal structure and uncertainty into compact latent embeddings. These representations are then used to inform downstream policy learning, optimizing for cumulative rewards. This direction also aligns with developments in meta-reinforcement learning and non-stationarity, where belief states or Bayesian embeddings are used to capture hidden task contexts. Agents trained across a distribution of tasks can use these latent variables to infer new environments and adapt quickly (Zintgraf et al., 2021; Nguyen et al., 2021; Huang et al., 2021; Rakelly et al., 2019; Xie et al., 2021; Feng et al., 2022; Feng & Magliacane, 2023; Liang et al., 2024a; Cao et al., 2025; Wang & Huang, 2025; Feng et al., 2026). For example, MetaDiffuser (Ni et al., 2023) incorporates task context as conditioning input to diffusion-based decision models. Similarly, Pertsch et al. (2021) and Zeng et al. (2023) use similar variational objectives (ELBO loss) to learn latent skill priors and predictive information for RL, where these latents greatly help policy learning. Our approach diverges from these by offering theoretical guarantees on the identifiability of latent factors from minimal temporal observations. Rather than depending on diverse multi-environment data, we introduce a framework that captures the full data generation process in RL using diffusion models. In contrast to MetaDiffuser, which assumes static task-level context, our model treats the latent context as a dynamic, time-evolving process that governs both environment transitions and agent behavior, capturing the underlying temporal structure of RL trajectories more faithfully. E.3

AUTOREGRESSIVE D IFFUSION M ODELS

To model temporal consistency and dynamics in sequential data such as videos and audios, recent work has incorporated autoregressive structures into diffusion models. These approaches differ in how they condition on prior time steps during generation and can be categorized into two main categories. (1) Conditioning on clean (denoised) inputs ((Zheng et al., 2024; Gao et al., 2024b; Blattmann et al., 2023)). At each time step t, the denoising model is conditioned on the previously denoised outputs {x0<t }: pθ (xk−1 | xkt , x0<t ), where xkt is the current noisy input, and x0<t denotes t the clean (fully denoised) observations from earlier time steps. (2) Conditioning on noisy inputs ((Ho et al., 2022; Chen et al., 2024; Xie et al., 2024b; Sand-AI, 2025)). These methods instead condition on previous time steps at their corresponding noise levels. This setting can be further divided into two cases: (a) fully noisy conditioning (Ho et al., 2022): the model conditions on all prior k−1 time steps at the same noise level k: pθ (xk−1 | xkt , xk<t , ). (b) partially noisy conditioning: <t , xt each previous time step i < t is conditioned at its own noise level ki , which may vary over time: 39

Published as a conference paper at ICLR 2026

pθ (xk00 −1 , xk11 −1 , . . . , xTkT −1 | xk00 , xk11 , . . . , xkTT ). Specifically, Diffusion Forcing (DF) (Chen et al., 2024) proposes a general framework in which each time step xt assigns an independent noise level. In contrast, other works adopt time-dependent noise schedules that vary with the temporal index (Xie et al., 2024b; Sand-AI, 2025; Wu et al., 2023). To model the causal generative process of RL trajectories, our approach also employs time-dependent noise scheduling to capture temporal dynamics. However, unlike prior work, we further integrate the identification of latent factors directly into the denoising process. This is achieved through a structured reinforcement step during training and a zig-zag inference procedure at test time, enabling our model to more faithfully recover the underlying causal structure in sequential decision-making. E.4

S UMMARY

To sum up, we compare our approach with representative diffusion- and meta-learning–based baselines (Table 13). Diffuser, DP, IDQL, and DD do not model or infer latent contexts; DF adopts autoregressive denoising but still lacks context inference. Meta-Diffuser, LILAC, and DynaMITE learn latents via meta-learning but omit our minimal–sufficient block design and backward refinement. LDCQ and LDP model only high-level latent actions/skills without explicit context identification. In contrast, our method jointly models latent factors, employs full autoregressive denoising with zig-zag sampling, and introduces a backward refinement mechanism that enables identifiable latent contexts. Method Ours Diffuser / DP / DD / IDQL DF Meta-Diffuser / LILAC / DynaMITE LDCQ LDP

Latent Factors

AR Denoising

Min. & Suff. Obs.

Yes (dyn., rew., act.) No No Yes (dyn., rew. only) Yes (hi-level act.) Yes (hi-level act.)

Yes No Yes No No No

Yes (refine, zig-zag) No No No No No

Table 13: Comparison with representative baselines on whether they model latent contexts, use autoregressive (AR) denoising, and enforce minimal & sufficient observation blocks.

F

B ENCHMARK S ETTINGS AND I LLUSTRATIONS

F.1

L ATENT C HANGE FACTORS D ESIGN

We consider the latent change factors on dynamics and rewards. We consider the Half-Cheetah and Ant environments from the OpenAI Gym suite, which are widely used MuJoCo locomotion benchmarks (Brockman et al., 2016) for evaluating continuous control algorithms. In Half-Cheetah, the agent is a planar bipedal robot with a 17-dimensional state space and a 6-dimensional continuous action space, where the goal is to move forward by applying torques to six actuated joints. In Ant, a quadrupedal robot operates in a 3D space with a 111-dimensional state space and an 8-dimensional action space, requiring more complex coordination across its four legs. In both environments, the reward encourages forward velocity while penalizing excessive control inputs and, in the case of Ant, also promotes stable contact with the ground. We consider variants of the Half-Cheetah environment to study changes in dynamics, specifically Cheetah-Wind-E and Cheetah-Wind-S, which introduce external wind forces applied to the agent. In Cheetah-Wind-E, an opposing wind force is applied at the beginning of each episode and remains constant throughout, defined as fw = 10 + 5 sin(0.8i), where i is the episode index. For this case, since c change over episode, we use data from several consecutive episodes to estimate it. In Cheetah-Wind-S, the wind force varies at every time step according to the same formula fw = 5 + 5 sin(0.5t), with t now representing the time step in each episode. We also consider variations in the reward function. In Cheetah-Dir-E, the reward depends on a time-varying goal direction, requiring the agent to alternate between moving forward and backward. Specifically, the reward at episode t is defined as rt = dt · vt − 0.1∥at ∥2 , where vt is the agent’s forward velocity, at is the action vector (torques applied), and dt ∈ {−1, +1} indicates the target direction at time t. The direction signal dt changes, giving a non-stationary reward 40

Published as a conference paper at ICLR 2026

Figure 8: Illustrations of the Benchmarks. From left to right: Half-Cheetah, Ant, Walker, FrankaKitchen, Maze2D, and LIBERO.

Figure 9: Illustrations of RoboMimic Benchmark. function that challenges the policy to adapt to shifting goals. Specifically, we consider dt = σ(5 · sin(2πt/200)),

where σ(·) denotes the sigmoid function, α controls the sharpness of the transition, and T determines the switching period. This formulation induces a smooth periodic change in the preferred direction of movement, requiring the policy to adapt to gradually shifting objectives. We also consider a directional reward variant for the Ant environment, denoted as Ant-Dir-E, where the agent is required to alternate its movement direction over time. The reward function at time step t is defined as rt = (2dt − 1) · vtx − 0.1∥at ∥2 ,

where vtx is the velocity of the agent’s torso along the x-axis (forward direction), at is the 8dimensional action vector, and dt ∈ [0, 1] is a smooth directional signal. Similarly, we define dt as: dt = σ(5 · sin(2πt/200)), where σ(·) denotes the sigmoid function. This formulation causes the preferred movement direction to alternate approximately every 100 steps. Notably, for these settings with periodic changes (i.e., where latent factors do not evolve at every timestep), we estimate the latent variables periodically and perform refinement in the causal diffusion model only when changes are detected. This follows the same overall framework, but operates at a coarser temporal resolution aligned with the latent change frequency. F.2

OVERVIEW ON OTHER B ENCHMARKS

Fig. 8-9 give the illustrations on the used benchmarks. Specifically, other than Cheetah and Ant we introduced before, for others, we consider the basic settings in offline RL. Specifically, Maze2D. Maze2D tasks focus on goal-directed navigation in a 2D plane, where the agent must traverse a maze-like environment to reach specified targets. These settings are designed to evaluate an agent’s ability to reason spatially and follow optimal trajectories based solely on positional and velocity observations. Franka-Kitchen. The Franka-Kitchen environment (Gupta et al., 2020) involves a robotic arm interacting with a series of articulated objects in a realistic kitchen setting. Tasks are composed of multiple stages, such as opening doors or toggling switches, and are intended to assess an agent’s capability in handling long-horizon, multi-step manipulation. 41

Published as a conference paper at ICLR 2026

Walker. The Walker2D environment features a two-legged robot that must learn to walk and balance using continuous torque control. The agent’s objective is to maintain forward motion while remaining upright, which requires learning dynamic stability and coordination. LIBERO (Liu et al., 2023). The Libero benchmark offers a diverse set of continual learning tasks focused on object manipulation and generalization: • LIBERO-Object: The robot is required to manipulate a variety of novel objects through pick-and-place operations. Each task introduces previously unseen objects, encouraging the agent to incrementally build knowledge about object-specific properties and behaviors. • LIBERO-Goal: All tasks share a common object set and spatial layout, but vary in goal specifications. This setup tests the agent’s ability to continually adapt to new task intents and motion targets without changes in the visual scene. • LIBERO-Spatial: Tasks involve repositioning a bowl onto different plate locations. Although the objects remain fixed, the spatial configurations vary across tasks, requiring the robot to incrementally acquire relational spatial understanding. RoboMimic. RoboMimic (Mandlekar et al., 2021) provides a set of manipulation tasks based on human teleoperation demonstrations, varying in difficulty and required precision: • Lift: The robot arm is tasked with lifting a small cube off the table. This task serves as a foundational manipulation scenario focused on grasping and vertical motion. • Can: The robot must retrieve a cylindrical can from a cluttered bin and place it into a designated smaller container. This task introduces greater complexity due to object shape and the need for accurate placement. • Square: A fine-grained insertion task where the robot picks up a square nut and places it onto a vertical rod. This is the most challenging of the three, requiring precise alignment and control for successful completion.

G

OTHER D ETAILS ON A D A -D I F F U S E R

G.1

L ATENT ACTION P LANNER

For the latent action planner, we align our settings with those used in LDP (Xie et al., 2025), specifically focusing on learning directly from image-based demonstrations. We first use a variational autoencoder (VAE) to extract latent representations z from raw images via image encoders. An inverse dynamics model is then trained to recover actions at from pairs of latent states (zt , zt+1 ). A planner is subsequently trained to forecast future i latents z. Hence, the objective function is h 2 LIDM (ξ, z) = Et,ϵ ∥ϵξ (âk ; ck , zk , zk+1 , t) − ϵ∥ , where where k is the time step and t is the diffusion step. In our framework, we treat the latent factors c as high-level latent actions that influence the evolution of z. These latent factors are jointly used with z to perform both inverse dynamics modeling and latent forecasting, enabling structured planning in the latent space. We follow the experimental settings established in LDP (Xie et al., 2025). Specifically, we use expert demonstrations alongside action-free and suboptimal demonstrations. All hyperparameters and architectural choices for the diffusion models are kept identical to those used in the original LDP implementation. We also directly utilize the pre-trained image encoder provided by LDP. The only modification in our framework is the introduction of an additional latent factor c trained by our latent factor identification stage, which is incorporated into the model to enhance latent action planning. G.2

N OISE S CHEDULING

In the autoregressive setting, we consider a monotonic increasing denoising schedule {k1 , . . . , kT }. In practice, we use a linear schedule where ki = Ti K, with K denoting the maximum diffusion step used in both training and sampling. We segment the sequence into temporal blocks of length 42

Published as a conference paper at ICLR 2026

Tx + 1 (Tx = To in all settings), and slide the time window forward by one step at a time. This design ensures that the denoising steps progressively increase across the block, aligning the diffusion process with the underlying temporal structure. Such a schedule encourages early steps to rely more on strong priors and later steps to refine based on more contextual information. Additionally, for better illustration, Fig. 7 provides a detailed illustration of the zig-zag sampling process within a temporal block of 4 timesteps.

H

S PECIFIC D ESIGN C HOICES FOR BASELINES

For all baselines, unless otherwise specified, we use the same set of diffusion parameters detailed in Appendix D.2.2–D.2.3. Below, we provide additional details on how specific methods are evaluated. While their diffusion backbones remain consistent as in Appendix D.2.2–D.2.3, these methods include custom design choices and method-specific hyperparameters that are evaluated accordingly. H.1

D ETAILS ON LILAC AND DYNA MITE

In these settings, we extend both LILAC and DynaMITE by incorporating a context encoder to infer latent context variables ct , following their respective designs. Both methods learn belief state embeddings from historical observations. For a fair comparison, we use the same latent identification network architecture as in our framework, but modify the inputs according to each method’s assumptions. Specifically, LILAC and DynaMITE condition their inference networks solely on the historical trajectory x1:t , without access to current and future information. Additionally, consistent with their original implementations, we do not include a separate prior head on top of the GRU; both methods share the encoder for posterior inference and prior prediction. And the primary difference (in implementation) between these two methods lies in the temporal context used: LILAC maintains the full belief over the entire history, i.e., it conditions on x1:t to infer ct+1 , while DynaMITE uses only the most recent context, i.e., it infers ct+1 based solely on xt . All other hyperparameters are aligned with those used in our Stage 1 training. The estimated context variables are then provided as additional conditioning inputs to the diffusion-based models. H.2

D ETAILS ON D IFFUSION F ORCING

For Diffusion Forcing, we adopt the same autoregressive noise schedule as in our method, which accounts for causal uncertainty, similarly to the formulation in Eq. D.1 of (Chen et al., 2024), to ensure a fair comparison. Additionally, we use the Monte Carlo Guidance (MCG) mechanism introduced in (Chen et al., 2024) for Maze2D, following the original setup. For all other environments, we use the same classifier guidance scheme as the other baselines to maintain consistency in evaluation.

I

A BLATION A NALYSIS

I.1

T RAINING /I NFERENCE T IME A NALYSIS

We conduct all experiments on 4× NVIDIA A100 or 8× RTX 4090 GPUs, depending on the model scale and environment requirements. The main computational overhead in our framework arises from two components: (i) the latent factor identification network, and (ii) the denoise-and-refine steps in the diffusion model. During sampling, the additional cost comes from zig-zag latent exploration and latent variable sampling. However, these steps do not substantially increase either training or inference time. To quantify this, we report the training and inference speed of our method compared to the base models DD and DP across all environments (Table 14). Our framework introduces only a moderate computational overhead — typically 1.2–1.3× the runtime of vanilla diffusion backbones, corresponding to roughly 20–30% extra training time and inference latency. This cost can be further reduced through parallel latent sampling, lightweight context encoders, or refinement only at inference. Moreover, we additionally evaluate a Picard-accelerated variant (Table 15, Shih et al. (2023)), where iterative refinement is parallelized by conditioning each denoising step on previously denoised 43

Published as a conference paper at ICLR 2026

nodes. With Picard iteration, inference time drops to 0.7–0.8× of our default iterative sampler while maintaining comparable performance, demonstrating the potential for further acceleration.

Environment Cheetah Ant Walker Maze2D Libero Kitchen

Training Time (sec/epoch) Ours vs DD Ours vs DP 72.1 / 60.1 (1.20) 79.5 / 64.3 (1.24) 85.3 / 67.1 (1.27) 90.2 / 72.0 (1.25) 104.0 / 81.0 (1.28) 117.8 / 88.1 (1.34)

69.8 / 58.4 (1.20) 76.0 / 62.0 (1.23) 81.5 / 64.2 (1.27) 88.3 / 69.2 (1.28) 102.1 / 78.0 (1.31) 115.3 / 85.0 (1.36)

Inference Latency (ms) Ours vs DD Ours vs DP 182 / 114 (1.16) 148 / 125 (1.19) 182 / 144 (1.28) 184 / 149 (1.24) 209 / 169 (1.24) 228 / 180 (1.27)

160 / 125 (1.28) 172 / 139 (1.24) 170 / 130 (1.31) 196 / 152 (1.29) 219 / 162 (1.35) 211 / 168 (1.26)

Table 14: Training and inference time comparison for Ada-Diffuser-planning and Ada-Diffuser-policy variants. We report absolute times (sec/epoch or sec/rollout) and relative overheads.

Environment Cheetah Ant Walker Maze2D Libero Kitchen

Ours (sec)

Ours+Picard (sec)

1.51 1.67 1.83 1.94 2.18 2.45

1.15 1.25 1.40 1.47 1.62 1.84

Table 15: Picard-accelerated inference.

I.2 I.2.1

A BLATION R ESULTS F ULL R ESULTS S UPPLEMENT TO TABLE 2

Table 16 presents the full ablation results across all environments, as a supplement to Table 2. Overall, the results highlight the importance of the two key components in causal diffusion modeling: latent identification and autoregressive diffusion, both of which are critical for performance.

I.2.2

N OISE S CHEDULE : L INEAR VS . L OGISTIC VS . S IGMOID

We adopt a linear noise schedule by default since any monotonic, bounded schedule suffices to model the data-generation process and linear is simple and stable in practice. To validate this choice, we compare linear, logistic, and sigmoid schedules on three representative tasks. As shown in Table 17, performance remains stable across schedules with no significant differences, supporting our default choice. 44

Published as a conference paper at ICLR 2026

Task

Schedule

Performance

Cheetah

linear logistic sigmoid

Maze2D

linear logistic sigmoid

−68.9 −63.6 −70.4

Franka-Kitchen

linear logistic sigmoid

161.4 157.6 168.5 0.70 0.72 0.66

Table 17: Ablation on noise schedules. “Performance” is the task score (higher is better for Maze2D/Kitchen; lower magnitude negative is better for Cheetah as per the benchmark).

I.2.3

E FFECT OF T EMPORAL B LOCK L ENGTH ON L ATENT I DENTIFICATION

We further analyze the impact of temporal block length on latent identification. As shown in Fig. 10, the results are consistent with findings reported in the main paper. When the number of observations is insufficient (e.g., ≤ 4), identification performance degrades. Performance improves when the block length is in a moderate range (5–20), indicating that sufficient temporal context is beneficial. However, using overly long blocks (> 20) introduces redundancy and increases optimization difficulty, which in turn harms performance. MSE R2

0.8

0.8

0.6

0.6

0.4

0.4

MSE R2

0.2

0.2 3 4 5 6 7 8 9 10 15 20 25 30 40 50 100

3 4 5 6 7 8 9 10 15 20 25 30 40 50 100

Size

Size

Figure 10: Identification results (MSE of linear probing and R2 ) versus the length of temporal blocks. Left: Cheetah with time-varying wind; Right: Cheetah with time-varying rewards.

Clustering We assess whether the learned latent space organizes states by the underlying context on the Cheetah wind-change task, where the ground-truth latent evolves as fw (t) = 5 + 5 sin(0.5t). We sample 1000 time steps, discretize fw (t) into five equal-frequency bins to define target clusters, embed the corresponding observations into the 20-dimensional learned latent representation, and run k-means with k = 5. We compare our method with LILAC and DynaMIE, together with an ablation that without refinement. Results are given in Fig. 11. I.2.4

L ATENT P ROBING : E FFECT OF BACKWARD R EFINEMENT AND Z IG –Z AG

To test whether backward refinement and zig–zag primarily help by correcting posterior mismatch, we perform a latent probing analysis on C HEETAH with changing wind. We linearly map the learned latent representation to the ground-truth wind variable using a simple least-squares probe (trained on a subset of blocks and evaluated on held-out blocks). Table 18 reports the mean squared error (MSE) of this probe for three variants: (i) the full model with backward refinement and zig–zag; (ii) without refinement; and (iii) without zig–zag. 45

Published as a conference paper at ICLR 2026

Ours

LILAC

Ours (w/o Re ne)

DynaMITE s

Figure 11: Clustering (t-SNE)results on Cheetah wind-change.

fi

Variant

MSE

Full (with refinement + zig–zag) w/o refinement w/o zig–zag

0.18 0.28 0.23

Table 18: Linear probing MSE for recovering the ground-truth wind latent on C HEETAH (changing wind). Lower is better.

Analysis. The full model achieves the lowest MSE, indicating more accurate recovery of the latent wind. Removing backward refinement yields the largest degradation (0.18 → 0.28), consistent with the role of refinement in letting future evidence within a block update the latent posterior and reduce temporal lag. Disabling zig–zag also harms accuracy (0.18 → 0.23), suggesting that alternating conditioning helps align the denoising trajectory with the latent dynamics rather than purely following the forward temporal pass. Together, these results support our claim that both components reduce posterior mismatch and improve latent identifiability, which in turn benefits planning and control in settings with evolving hidden factors. I.2.5

O N THE E FFECT OF P LANNING AND E XECUTION H ORIZONS : L ONG - HORIZON P LANNING

We study the robustness of our approach under increased planning and execution horizons (Tp and Ta ). Specifically, we evaluate on two challenging tasks—Franka-Kitchen-Partial and LiberoLong, where the original settings are Kitchen (Tp = 16, Ta = 8) and Libero (Tp = 10, Ta = 8). Results are in Fig. 12. When we increase these horizons, we observe that the baselines, DP and DF, suffer significant performance drops. In contrast, Ada-Diffusermaintains relatively high performance. This demonstrates that modeling the underlying causal generative process, through autoregressive structure and latent representations, enables better long-horizon planning. Although 46

Published as a conference paper at ICLR 2026

Kitchen-partial

Kitchen-partial 70

60 Ours DF

55 50

Success Rate

Success Rate

65

8

10

20

55 50

32

Ours DF 8

10

20

Execution Horizon

Libero-Long

Libero-Long

50 Ours DP 8

60

Planning Horizon

60

40

65

45

Success Rate

Success Rate

70

60 50 40 Ours DP

30 10

20

32

8

Planning Horizon

32

10

20

32

Execution Horizon

Figure 12: Results with different planning and execution horizons. We evaluate on Kitchen-partial and Libero-Long experiments. we do not explicitly impose latent variables, our model implicitly learns representations that can track stochasticity and support smooth control.

47

-120.4 ± 12.7 -148.5± 9.8 850.8 ± 54.2 -102.4± 18.2 188.6± 39.2

Cheetah-Dir-E (cr ) Cheetah-Vel-E (cr ) Ant-Dir-E (cr )

Diffuser

Cheetah-Wind-E (cs ) Cheetah-Wind-S (cs )

Environment

48 902.1 ± 45.8 -85.6± 18.3 195.4± 47.0

-105.8 ± 9.6 -102.0± 10.2

DF

912.5± 37.9 -69.2± 7.5 245.9± 41.0

-89.7 ± 6.5 -106.8± 11.4

MetaDiffuser

930.4± 29.5 -76.3± 11.7 262.8± 27.5

-79.2 ± 11.0 -94.3± 9.6

Diffuser + DynaMITE

908.5± 37.6 -62.6± 11.1 229.4± 32.6

-95.3± 7.4 -105.6± 14.5

Diffuser + LILAC

−62.4± 3.9 -65.3± 11.2

-68.9± 7.6 -73.5± 8.7

−57.8± 6.7 -58.1± 9.0

Oracle

943.3± 25.6 949.8± 24.1 962.1± 21.9 -45.8± 9.5 -39.2± 7.6 -38.3± 8.9 285.3± 24.5 296.4± 22.2 300.7± 23.6

Ours + Meta Ours

Published as a conference paper at ICLR 2026

Table 6: Results (average returns) on Ada-Diffuser-Planner with latent factors that affects dynamics and rewards. cs and cr indicate the changes on dynamics and reward, E and S represent the episodic and time-step changes. The results are with 5 random seeds. The bold ones are the best-performing ones, excluding meta-learning and oracle ones.

-104.8 ± 10.9 -120.6± 11.5 892.5± 60.8 -87.9± 6.5 182.5± 41.2

Cheetah-Dir-E (cr ) Cheetah-Vel-E (cr ) Ant-Dir-E (cr )

DP

Cheetah-Wind-E (cs ) Cheetah-Wind-S (cs )

Environment

49 949.6± 36.1 -72.7± 5.8 275.2± 27.0

-72.2± 5.9 -76.5± 15.6

DP + DynaMITE

Ours + DP (Oracle) -52.0± 3.5 -42.3± 6.7 972.4± 37.5 -39.8± 6.7 312.5± 37.2

Ours + DP -58.5± 4.6 -52.9± 9.8 960.7± 40.2 -41.0± 7.2 290.4± 49.4

902.4± 45.2 -80.2± 11.4 204.6± 25.6

-97.5± 9.4 -87.8± 12.2

IDQL

938.6± 49.4 -59.4± 6.5 269.3± 29.5

-59.0± 11.2 -63.4± 6.7

IDQL + DynaMITE

965.0± 37.5 -38.6± 7.7 295.8± 32.7

-48.5± 7.9 -48.0± 7.2

Ours + IDQL

969.8± 39.2 -33.8± 6.5 309.6± 25.4

-41.6± 6.2 -44.7± 6.1

Ours + IDQL (Oracle)

Published as a conference paper at ICLR 2026

Table 7: Results (average returns) on Ada-Diffuser-Policy with latent factors. cs and cr indicate the changes on dynamics and reward, E and S represent the episodic and time-step changes. The results are with 5 random seeds. The bold ones are the best-performing ones, excluding metalearning and oracle ones.

Published as a conference paper at ICLR 2026

Environment

Diffuser

DD

DF

LDCQ

Ours (DD)

Mixed Partial

52.6 ± 2.3 55.8± 1.9

75.2 ± 1.4 57.3± 1.2

73.7 ± 1.9 68.6± 2.4

73.3 ± 0.5 67.8± 0.8

74.6± 1.6 70.1± 1.3

Table 8: Results (success rate (%)) on Ada-Diffuser-Planner without explicit latent factors on Franka-kitchen environment. The results are with 5 random seeds.

Environment

Diffuser

DD

DF

LDCQ

Ours (DD)

umaze medium large

113.5 ± 2.8 121.5 ± 5.6 123.0± 4.8

114.8 ± 3.2 129.6 ± 2.9 131.5± 4.2

116.7 ± 2.0 149.4± 7.5 159.0± 2.7

134.2 ± 4.1 125.3± 2.5 150.1± 2.9

148.6± 3.7 148.6± 3.1 161.4± 3.2

Table 9: Results on Ada-Diffuser-Planner without explicit latent factors on Maze-2D environment. The results are averaged across 5 random seeds.

Environment

Diffuser

DD

DF

LDCQ

Ours (DD)

medium-expert medium medium-replay

106.2 ± 0.7 79.6± 9.8 70.6± 0.6

108.8 ± 2.0 82.5± 1.6 75.0 ± 3.2

105.4 ± 3.2 66.2± 1.9 72.2± 2.6

109.3 ± 0.4 69.4± 2.4 68.5± 4.3

115.7 ± 2.1 83.6± 3.5 74.3± 2.8

Table 10: Results on Ada-Diffuser-Planner without explicit latent factors on Walker environment. The results are averaged across 5 random seeds.

Environment Spatial Object Long

DP

Ours (DP)

78.3 ± 3.9 92.5± 2.6 50.5 ± 7.2

79.2 ± 4.2 93.4± 2.8 62.6 ± 4.9

Table 11: Results on Ada-Diffuser-Policy without explicit latent factors on Libero environment. The results are averaged across 5 random seeds.

Component

Type (i): Full Trajectory

Type (ii): State-Only

Model Backbone Architecture # DiT Blocks Optimizer Batch Size Training Steps Diffusion Timesteps Observation Horizon To Planning Horizon Tp Execution Horizon To Guidance Inverse Dynamics Model Refinement Loss Cofficient

1D U-Net (Ronneberger et al., 2015) Kernel size: 5; channels: (1,2,2,2); base: 32 – Adam, lr = 3 × 10−4 64 1M 150 10 16 (Cheetah), 32 (Ant) 1 CG, ω = 1.5 – 0.1

Transformer (DiT) Hidden dim: 256; head dim: 32 2 (Walker, Kitchen, Maze2D), 8 (LIBERO) Adam, lr = 3 × 10−4 128 1M 200 4 (Kitchen), 2 (LIBERO), 10 (others) 16 (Kitchen), 10 (LIBERO), 32 (others) 8 (Kitchen, LIBERO), 10 (others) CFG 2-layer embed (64), 3-layer MLP (128), 1M steps 0.1

Table 12: Planner configurations for type (i): full trajectory generation and type (ii): state-only generation with inverse dynamics.

50

Published as a conference paper at ICLR 2026

Cases Original w/o refine w/o zig-zag same NS random NS

Cheetah-1

Cheetah-2

Ant

Maze2D

Walker

Kitchen

RoboMimic

LIBERO

-73.5 -82.0 -91.6 -89.7 -84.6

-52.9 -60.7 -56.1 -62.4 -62.9

295.8 261.2 258.3 259.7 266.4

161.4 156.5 147.6 140.1 146.3

115.7 107.4 107.9 105.8 109.1

0.70 0.63 0.59 0.56 0.61

0.85 0.78 0.75 0.72 0.76

93.4 90.2 91.6 85.2 88.5

Table 16: Ablation on Design Choices. We conduct ablation studies across a diverse set of tasks, including: Cheetah-Wind-S with a planner-based approach (denoted as Cheetah-1 in the table), Cheetah-Wind-S with a diffusion policy (Cheetah-2), Ant-Dir-E (policy, IDQL-based), Maze2DLarge (planner), Walker2D-Medium-Expert (planner), Kitchen-Partial (planner), LIBERO-Object (diffusion policy), and RoboMimic-Can.

51

Record · ID 192411 · SHA-256 3488fe0757ebedd4
Retrieved via Conceptio — every document is proof-bundled with source, license, and retrieval metadata.