All projects

World2vec

World models can achieve strong control performance without faithfully capturing the dynamics of their environment, limiting their potential for mechanistic interpretability and scientific discovery. Here, we introduce world2vec to learn mechanistic world models without reward information directly from observations and actions. World2vec learns latent dynamics that faithfully reproduce the true system over long prediction horizons, with at least one order of magnitude fewer parameters than DreamerV3, TD-MPC2 and LeWM. This enables state of the art performance in controlling continuous environments from pixels and accurate system identification across a benchmark of 71 attractors, demonstrating how world models are amenable to scientific discovery.

The world2vec framework for mechanistic world modeling

A physical state s\boldsymbol s evolves under actions a\boldsymbol a as s˙=f(s,a)\dot{\boldsymbol s}=\mathbf f(\boldsymbol s,\boldsymbol a) and is observed through a nonlinear mixing x=g(s)\boldsymbol x=\mathbf g(\boldsymbol s). World2vec learns an encoder h\mathbf h that maps observations to latent states z=h(x)\boldsymbol z=\mathbf h(\boldsymbol x), and a latent vector field f^\hat{\mathbf f} with z˙=f^(z,a)\dot{\boldsymbol z}=\hat{\mathbf f}(\boldsymbol z,\boldsymbol a). The model is faithful if evolving the encoded state in latent space gives the same trajectory as encoding the evolved physical state:

z^t(T)=z(t)+∫0Tf^(z^t(τ),a(t+τ)) dτ  =!  h(x(t+T))\hat{\boldsymbol z}_t(T)=\boldsymbol z(t)+\int_0^T\hat{\mathbf f}\big(\hat{\boldsymbol z}_t(\tau),\boldsymbol a(t+\tau)\big)\,d\tau\;\overset{!}{=}\;\mathbf h\big(\boldsymbol x(t+T)\big)

The condition must hold for all times tt and horizons T≥0T\geq0. Both sides can be computed from observations and actions alone, without the physical state or the true dynamics. In training, the encoder maps a short history of frames Xn\boldsymbol X_n to the initial latent (a single frame does not show velocities), Euler steps advance it under the recorded actions, and the rollout is matched to the encodings of the observed future:

z^n,0=h(Xn),z^n,j+1=z^n,j+Δt f^(z^n,j,an+j),min⁡h,f^  En[∑j=1Kwj ∥h(Xn+j)−z^n,j∥22]+λ R(h,f^)\begin{gathered} \hat{\boldsymbol z}_{n,0}=\mathbf h(\boldsymbol X_n),\qquad\hat{\boldsymbol z}_{n,j+1}=\hat{\boldsymbol z}_{n,j}+\Delta t\,\hat{\mathbf f}(\hat{\boldsymbol z}_{n,j},\boldsymbol a_{n+j}),\\[4pt] \min_{\mathbf h,\hat{\mathbf f}}\;\mathbb E_n\Big[\sum_{j=1}^{K}w_j\,\big\|\mathbf h(\boldsymbol X_{n+j})-\hat{\boldsymbol z}_{n,j}\big\|_2^2\Big]+\lambda\,\mathcal R(\mathbf h,\hat{\mathbf f}) \end{gathered}

For a single step, this objective matches latent velocities up to a finite difference, as in neural ODEs. Over KK steps, it matches the encoded displacement with the velocities accumulated along the model's own predicted trajectory, so training exposes how small one-step errors compound. We call this principle latent dynamics matching.

A regularizer R\mathcal R prevents collapse, since a constant encoder satisfies the condition trivially. world2vec-dcl normalizes prediction scores against sampled negative embeddings (dynamics contrastive learning); world2vec-sigreg pulls embedding projections towards an isotropic Gaussian (SIGReg from LeJEPA).

Encoder, action-conditioned multistep latent dynamics, dynamics matching with an Euler step, and the SIGReg or DCL regularizers
Fig. 1A · World2vec jointly learns representations and dynamics in latent space. The dynamics model is an action-conditioned neural ODE, which can be parameterized to give an interpretable representation of the dynamics; the regularizer encodes assumptions about the data-generating process.

We measure the quality of the representation with consistency and the quality of the rollouts with LRC and EQA. Consistency is the R² between two independently trained encoders after an optimal affine alignment. LRC (latent rollout consistency) is the R² between latent rollouts and the encodings of the observed future. EQA (equivariance accuracy) is the top-1 retrieval accuracy of the encoded future among sampled negatives; unlike LRC, a collapsed embedding cannot score high on it. Hyperparameters are selected on validation EQA.

Pendulum, Reacher and CartPole with their physical variables
Fig. 2 · Pendulum, Reacher and CartPole from the DeepMind Control Suite, with their physical variables.

Setup: RGB frames from random policies on each environment. A small vision transformer encodes a short history of frames into a low-dimensional latent, and an MLP parameterizes the vector field; the whole model is at least an order of magnitude smaller than the baselines and trains in minutes on a single GPU. Planning uses the cross-entropy method towards an encoded goal, without a reward model or state decoder.

Planning without a reward model or state decoder

Model-predictive control on Pendulum swing-up. At every step, the planner encodes the last frames into a latent, rolls out candidate action sequences with the learned dynamics, refines them with the cross-entropy method, executes the first action of the lowest-cost plan and replans. The cost is the latent distance to the encoded upright frame.

MPC demo: environment frames, imagined future vs. reality, and the planner's candidates in latent space
Loading interactive demo…

One replayed Pendulum episode, planned with world2vec-dcl and with world2vec-sigreg. Averages over seeds are in Table 1 below.

World2vec faithfully models dynamics and enables long-range predictions and control

World2vec predicts faithfully over long horizons and enables model-predictive control from pixels on Pendulum and Reacher, where LeWM, the previous best JEPA-based approach, earns low rewards. Training on multistep rollouts with an Euler-integrated vector field turns one-step baselines that lose track of the system within a second into models that stay consistent with it.

Reward versus EQA after one second for world2vec with DCL and SIGReg, LeWM, TD-MPC2 and DreamerV3
Fig. 1B · Reward versus EQA after a 1 s rollout, averaged across environments.

On Pendulum and Reacher, world2vec-dcl is competitive with DreamerV3 and TD-MPC2, which are trained online with reward, and far ahead of LeWM; world2vec-sigreg solves both tasks as well. CartPole remains an open challenge for reward-free world models, although it can be solved in principle from non-linearly mixed observations (Fig. 6).

DreamerV3 and TD-MPC2 reach high reward but mostly predict poorly in latent space. LeWM and world2vec are trained offline on random-policy data, the one-step baselines are LeWM-like variants within the same framework, and all models are evaluated on the same held-out test episodes. World2vec's representations are also more consistent across training runs than those of the reward-trained models.

Reward
Method# paramsCons.LRCEQAPendulumReacherCartPole
DreamerV317M43.452.3±1.452.7±0.5668±140792±25899±6
TD-MPC25M87.515.5±2.412.1±0.5911±2744±41894±9
LeWM18M99.50.0±0.01.7±0.45±069±4171±5
Baseline (SIGReg ≈ LeWM)300k77.437.4±1.720.3±1.212±1087±6201±4
Baseline (DCL)300k88.351.2±0.811.1±0.23±3104±8154±8
world2vec-sigreg300k96.996.1±0.067.9±1.0434±142570±20221±1
world2vec-dcl300k95.597.7±0.170.7±0.7826±0895±16232±5
Table 1 · Consistency, and LRC and EQA after a 1 s rollout, in % on the test split, averaged over Pendulum, Reacher and CartPole, and reward per environment (dm_control return, maximum 1000; each row's best planner; on Reacher, the average over both inverse-kinematics solutions). Mean ± bootstrap standard error over model seeds, with configurations selected on validation EQA. On Reacher, the baselines and world2vec run at 8.3 Hz, TD-MPC2 and DreamerV3 at 25 Hz, and LeWM, trained on the original LeWM data, at 5 Hz. The baselines predict one step with discrete dynamics; world2vec predicts many steps with Euler integration.
LRC and EQA versus rollout horizon in seconds on Pendulum, Reacher and CartPole for world2vec, LeWM, TD-MPC2 and DreamerV3
Fig. 3 · LRC and EQA versus rollout horizon on Pendulum, Reacher and CartPole; models and selection as in Table 1. Lines show the mean over seeds, bands the bootstrap standard error. The gains are largest on Pendulum and CartPole, which have substantial intrinsic dynamics.
Success
Method25 steps200 stepss / ep.
LeWM81.6±2.452.9±3.88.64±0.12
world2vec-sigreg, 25 Hz85.8±1.072.4±0.22.61±0.06
world2vec-dcl, 25 Hz86.2±1.180.2±1.62.55±0.03
world2vec-sigreg, 8.3 Hz78.9±1.476.0±2.21.08±0.07
world2vec-dcl, 8.3 Hz75.8±2.777.6±2.81.00±0.01
Table 2 · Reacher under the LeWM goal-reaching protocol: success rate (%) for goals 25 or 200 steps ahead (budget kept at 50 steps), and planning time per episode (s, for goals 25 steps ahead). LeWM is trained on its original data. Mean ± bootstrap standard error over training seeds.

Because LeWM earns low rewards in this evaluation, we also test Reacher under LeWM's own goal-reaching protocol. There, world2vec matches LeWM for nearby goals, is clearly better for distant ones and plans several times faster.

A sufficiently large training horizon and latent dynamics matching are essential for latent world modeling

Compared with the one-step baseline, world2vec trains on multistep rollouts and uses a vector field integrated by Euler steps as the predictor. Both components improve the model.

Increasing the training horizon drastically improves both world2vec-dcl and world2vec-sigreg, most of all on CartPole, where one-step training fails entirely. The exact horizon matters little, as long as training goes well beyond single-step prediction.

We test the second component by replacing the Euler step with discrete transitions, z^n,j+1=MLP(z^n,j,an+j)\hat{\boldsymbol z}_{n,j+1}=\mathrm{MLP}(\hat{\boldsymbol z}_{n,j},\boldsymbol a_{n+j}), in which the MLP predicts the next latent state directly. The Euler models are better on consistency, EQA and state decoding for both variants. Discrete dynamics still reach a high LRC when the latent space is as small as the physical state.

The analysis of latent world models found that standard predictors do not capture velocity information. Modeling the latent change rather than the next latent state counteracts this failure.

Embedding geometry qualitatively reflects environment structure

The Reacher arm is described by two angles. During data collection θ₀ performs full rotations, while θ₁ only swings left and right relative to the first arm. The embedding reflects this, with θ₀ encoded with circular boundaries and θ₁ represented linearly with open boundaries.

PCA of Reacher embeddings colored by the two joint angles
Fig. 4 · PCA of Reacher embeddings, colored by θ₀ (left, circular) and θ₁ (right, linear).

Continuous-time ODE parameterization allows flexible sampling

Because world2vec models the rate of change, the integration step can differ between training and planning. A Reacher model trained at one control rate also plans at a finer rate when the Euler step is reduced accordingly, without retraining, and controls the arm at least as well.

Reacher control at the native 10 Hz and retimed to 50 Hz, with planner observations, joint torques and replanning ticks
Fig. 5 · The same world2vec-dcl model, start state and target, planned at the rate it was trained on (top) and retimed to a finer rate by reducing the step size (bottom). Frames are the planner's observations at the time they were taken; below each strip, the two joint torques and a tick at every replanning.

Empirical identifiability of world2vec

No reward-free world model has solved CartPole from raw observations. Following the analysis of latent world models, we replace the pixels with an increasingly non-linear mixing of the state, which makes the representation learning task harder with each layer.

For these objectives, theory predicts representations that are consistent up to an affine transformation and dynamics that are modeled faithfully over time. The state itself is not expected to be recovered linearly, since the noiseless systems lack variability along some dimensions.

As the mixing deepens, LRC and EQA decrease gradually, with world2vec-dcl degrading least, while the one-step baselines collapse as soon as the mixing has more than one layer. Reward declines only slightly, so both world2vec variants still solve CartPole from non-linearly mixed observations, whereas the baselines fail. The gap between EQA here and under pixel observations is modest, which suggests that the result can transfer to pixels with better estimation of the dynamics model.

Consistency, LRC and EQA after one second, and reward, versus mixing depth 1 to 5 on CartPole for world2vec-dcl, world2vec-sigreg and their one-step baselines
Fig. 6 · World2vec across increasingly non-linear observations of the CartPole state. More mixing layers make the representation learning task harder, because the state is transformed more non-linearly. From left: consistency, LRC and EQA after a 1 s rollout, and reward, versus mixing depth; solid lines show world2vec, dashed lines the one-step baselines (purple: DCL, orange: SIGReg).

Setup: the positions and angles of the state, but not its velocities, pass through a randomly initialized, injective MLP into a higher-dimensional observation. An MLP encodes the observation history in place of the vision transformer.

Symbolic regression for understanding latent mechanisms

With a sparse polynomial library as the vector field (SINDy), world2vec-SINDy learns dynamics that can be read as equations. We evaluate it on chaotic attractors from the dysts database, with the states mixed nonlinearly into high-dimensional observations and perturbed by process noise of scale α. Training uses one-step prediction with RK4 integration and no actions.

Consistency, LRC, affine recovery and dynamics recovery on 71 dysts systems, their dependence on process noise, and recovered SprottA equations
Fig. 7 · (a) Consistency, 50-step LRC, affine recovery R² and iDynR²₅₀ for 71 dysts systems at process noise α = 0.1, means over three seeds, ordered by decreasing SINDy LRC; gray circles: identity, blue squares: MLP, red triangles: SINDy dynamics. (b) Medians and interquartile ranges across systems for α ∈ {0, 0.1, 0.5, 1}. (c) SprottA: true states and their optimal linear reconstruction from 50-dimensional observations; learned embeddings (gray) and predicted rollouts (red) in ground-truth coordinates, with the equations of each model below.

World2vec-SINDy improves substantially over an identity predictor and performs comparably to an MLP vector field, and all three give highly consistent ensembles. Process noise adds variation along every state dimension and enables linear state recovery: with more noise, affine recovery and consistency increase for all three models. Noise also caps LRC, so at high noise we judge the dynamics against the true rollouts (iDynR²), where MLP and SINDy stay far above the identity baseline.

Unlike the MLP, SINDy yields equations. For SprottA, a system whose coefficients agree well across seeds, the recovered equations approximate the true system after a least-squares alignment to the true states and removal of small coefficients. Expressing the equations in ground-truth coordinates uses the true states for this alignment; recovering them without the true states, for example with unsupervised symbolic regression, is left for future work.

Read the paper

First page of “world2vec: Mechanistic world models through latent dynamics learning” (opens the PDF)
Preprint, 2026

world2vec: Mechanistic world models through latent dynamics learning

Tobias Schmidt*, Rodrigo González Laiz*, Krisztián Szuppinger, Steffen Schneider