pith. sign in

arxiv: 2605.20780 · v1 · pith:IR6ZRZRYnew · submitted 2026-05-20 · 💻 cs.LG · cs.CV

Learning to Think in Physics: Breaking Shortcut Learning in Scientific Diffusion via Representation Alignment

Haozhe Jia , Pengyu Yin , Wenshuo Chen , Shaofeng Liang , Lei Wang , Bowen Tian , Xiucheng Wang , Nanqian Jia
show 1 more author
Yutao Yue
This is my paper

Pith reviewed 2026-05-21 06:20 UTC · model grok-4.3

classification 💻 cs.LG cs.CV
keywords physics-informed diffusion modelsshortcut learningrepresentation alignmentPDE residual supervisionintermediate feature alignmentscientific diffusiondiffusion transformers
0
0 comments X

The pith

REPA-P aligns intermediate features with PDE residuals to break shortcut learning in physics diffusion models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Physics-informed diffusion models typically supervise only their final outputs with physical constraints, allowing intermediate representations to learn shortcuts that break under new boundary conditions. REPA-P addresses this by attaching lightweight 1x1 projection heads to selected intermediate layers, decoding their activations into physical quantities, and applying PDE residual losses to those quantities during training. These heads add no cost at inference since they are discarded afterward. Experiments on four PDE problems show faster convergence, lower residual errors, and stronger robustness to distribution shifts, with benefits seen on both convolutional and transformer diffusion backbones. Supervising just a few layers captures most of the gain and works alongside standard output supervision.

Core claim

By decoding hidden activations from intermediate layers into physical states via simple projection heads and enforcing first-principles PDE residuals on those states, diffusion models avoid learning shortcuts and achieve better adherence to physics across training and out-of-distribution cases.

What carries the argument

Lightweight 1×1 projection heads attached to selected intermediate layers that decode activations into physical quantities for applying PDE residual losses.

If this is right

  • Convergence on PDE tasks accelerates by up to a factor of two.
  • Physics residuals decrease by as much as 66.4 percent.
  • Out-of-distribution robustness improves by up to 49.3 percent.
  • Performance gains appear consistently on U-Net and Diffusion Transformer architectures.
  • Most benefits come from supervising a small number of intermediate layers in addition to the output.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar intermediate supervision could apply to other generative models used for scientific simulation.
  • Embedding physics constraints earlier in the network might reduce reliance on large training datasets.
  • Extending the approach to three-dimensional or time-dependent problems would test its scalability.
  • Combining this with other physics-informed techniques could yield further gains in complex settings.

Load-bearing premise

Lightweight 1x1 projection heads can decode hidden activations into accurate physical quantities for residual computation without introducing major approximation errors.

What would settle it

Running the same experiments but measuring whether the decoded quantities from the projection heads match the true physical fields to high accuracy; if they do not, the residual losses would not enforce the intended constraints.

Figures

Figures reproduced from arXiv: 2605.20780 by Bowen Tian, Haozhe Jia, Lei Wang, Nanqian Jia, Pengyu Yin, Shaofeng Liang, Wenshuo Chen, Xiucheng Wang, Yutao Yue.

Figure 1
Figure 1. Figure 1: Overview of REPA-P. We decode intermediate diffusion features into physical states using lightweight projection heads and enforce PDE and boundary-condition residuals as supervision to align latent representations with valid physics. ing both understanding and generation quality, especially in out-of-distribution scenarios (Bastek et al., 2025). Based on this hypothesis, we propose Physics-Informed Represe… view at source ↗
Figure 2
Figure 2. Figure 2: Darcy flow qualitative comparison (Baseline vs. REPA-P). Top: baseline diffusion; bottom: REPA-P. Each row shows (left→right) the predicted pressure p, the permeability field K, and the PDE residual RMAE(K, p) (log scale). Compared to the baseline, REPA-P produces pressure fields that better respect the structure induced by K and achieves consistently lower residuals, indicating improved satisfaction of th… view at source ↗
Figure 3
Figure 3. Figure 3: Mechanics topology optimization (Baseline vs. REPA-P). Each row shows (left→right) the generated density ρ with CE (%) and mean ρ¯, the equilibrium residual RMAE(ρ, u1, u2) (log scale; lower is better), and the SIMP reference with compliance C and volume limit Vmax (red: displacement BCs; blue: load). REPA-P yields cleaner slender members and lower residuals than the baseline under the same volume constrai… view at source ↗
Figure 4
Figure 4. Figure 4: Training convergence curves on Darcy flow (first 60K of 120K total iterations). Left: Test data loss (log scale). Right: Physics residual error (log scale). REPA-P achieves significantly faster convergence and lower final loss on both metrics compared to the baseline. The shaded regions indicate standard deviation across 3 runs. Best viewed in color. pliance error (CE%), demonstrating that mid-layer align￾… view at source ↗
Figure 5
Figure 5. Figure 5: Physics residual (normalized, log scale) across U-Net layers. Baseline (red) applies physics loss only at output; Ours (blue) applies REPA-P alignment at intermediate layers, achieving 47%-100% reduction [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
read the original abstract

Physics-informed diffusion models typically enforce PDE constraints only on final outputs, leaving intermediate representations unconstrained and prone to shortcut learning under shifted boundary conditions. We introduce **REPA-P**, a teacher-free, architecture-agnostic framework that aligns intermediate features with physical states using first-principles residuals. REPA-P attaches lightweight $1{\times}1$ projection heads to selected layers, decodes hidden activations into physical quantities, and applies PDE residual losses during training. These heads are discarded at inference, introducing **zero overhead**. Across four PDE tasks, including Darcy flow, topology optimization, electrostatic potential, and turbulent channel flow, REPA-P accelerates convergence by up to $2{\times}$, reduces physics residuals by up to $66.4\%$, and improves out-of-distribution robustness by up to $49.3\%$, with consistent gains on both U-Net and Diffusion Transformer backbones. Ablations show that supervising a small set of intermediate layers captures most benefits and complements output-level physics losses. Code is available at [https://github.com/Hxxxz0/REPA-P](https://github.com/Hxxxz0/REPA-P).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes REPA-P, a teacher-free and architecture-agnostic method for physics-informed diffusion models. It attaches lightweight 1×1 projection heads to selected intermediate layers, decodes hidden activations into physical quantities, and applies PDE residual losses to align representations with first-principles physics. This is claimed to mitigate shortcut learning under shifted boundary conditions. Across four PDE tasks (Darcy flow, topology optimization, electrostatic potential, turbulent channel flow), the approach is reported to accelerate convergence by up to 2×, reduce physics residuals by up to 66.4%, and improve out-of-distribution robustness by up to 49.3%, with consistent gains on U-Net and Diffusion Transformer backbones and zero inference overhead after discarding the heads. Ablations indicate that supervising a small set of layers captures most benefits and complements output-level losses.

Significance. If the decoding mechanism proves accurate, REPA-P could meaningfully advance physics-informed generative modeling by constraining internal representations rather than outputs alone, offering a lightweight alternative to teacher-based or output-only supervision. The open code repository and consistent results across architectures and tasks are notable strengths that support reproducibility and potential adoption in scientific computing applications.

major comments (2)
  1. [Method and Experiments] The central mechanism depends on the 1×1 projection heads accurately decoding intermediate activations into physical fields without substantial approximation error or layer-specific bias. In the description of REPA-P and the turbulent channel flow experiments, no direct quantification of decoding fidelity (e.g., reconstruction error against ground-truth fields) is provided; without this, the reported residual reductions and robustness gains may reflect auxiliary supervision rather than genuine physics-aligned representations.
  2. [Ablations] The ablation studies on layer selection demonstrate benefits from intermediate supervision but do not include controls that isolate the contribution of accurate physical decoding versus generic regularization. This leaves open whether the PDE residual losses are enforcing first-principles constraints at the chosen layers or simply adding auxiliary objectives.
minor comments (2)
  1. [Method] The abstract and method would benefit from explicit equations showing how the decoded quantities enter the PDE residual computation, including any assumptions about boundary conditions or discretization.
  2. [Experiments] Consider reporting statistical significance or variance across multiple runs for the quantitative improvements (e.g., convergence speed and residual reductions) to strengthen the empirical claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and describe the revisions we will incorporate to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Method and Experiments] The central mechanism depends on the 1×1 projection heads accurately decoding intermediate activations into physical fields without substantial approximation error or layer-specific bias. In the description of REPA-P and the turbulent channel flow experiments, no direct quantification of decoding fidelity (e.g., reconstruction error against ground-truth fields) is provided; without this, the reported residual reductions and robustness gains may reflect auxiliary supervision rather than genuine physics-aligned representations.

    Authors: We agree that direct quantification of decoding fidelity would provide stronger evidence that the observed gains arise from physics-aligned representations. The current manuscript emphasizes end-task metrics (convergence speed, residual reduction, and OOD robustness) as the primary validation. In the revised version we will add explicit measurements of reconstruction error (e.g., relative L2 or MSE) between the outputs of the 1×1 projection heads and the corresponding ground-truth physical fields for the turbulent channel flow task, and, where computationally feasible, for the other PDE tasks as well. These new results will be presented in a dedicated subsection under Experiments. revision: yes

  2. Referee: [Ablations] The ablation studies on layer selection demonstrate benefits from intermediate supervision but do not include controls that isolate the contribution of accurate physical decoding versus generic regularization. This leaves open whether the PDE residual losses are enforcing first-principles constraints at the chosen layers or simply adding auxiliary objectives.

    Authors: We acknowledge the value of controls that separate physics-specific constraints from generic regularization. Our existing ablations already show that intermediate supervision yields gains beyond output-level physics losses alone. To isolate the role of the first-principles PDE residuals, we will add a new control experiment in the revised manuscript: at the same selected layers we will replace the PDE residual loss with a non-physics auxiliary objective (e.g., an L2 penalty on activations or supervision to random targets) while keeping all other training details identical. The performance difference between this control and the original REPA-P will be reported to clarify whether the first-principles nature of the loss is essential. revision: yes

Circularity Check

0 steps flagged

No circularity: external PDE residuals and empirical baselines keep derivation self-contained

full rationale

The paper introduces REPA-P by attaching 1x1 projection heads to decode intermediate activations and apply first-principles PDE residual losses during training, with heads discarded at inference. These residuals derive from standard physics equations external to the model outputs rather than being fitted or defined in terms of the paper's own predictions. Reported improvements (convergence speed, residual reduction, OOD robustness) are measured against standard output-level baselines on four PDE tasks, not quantities constructed solely from the method's internal choices. No self-citations, uniqueness theorems, or ansatzes from prior author work are invoked as load-bearing steps in the provided text. The central mechanism remains independently verifiable against external physics benchmarks and does not reduce to tautology by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The framework depends on the ability of simple linear projections to map internal features to physically meaningful quantities and on the assumption that supervising a small subset of layers is sufficient to capture most benefits.

free parameters (1)
  • choice of supervised layers
    Selection of which intermediate layers receive projection heads is a design decision that affects performance but is not derived from first principles.
axioms (1)
  • domain assumption Hidden activations in diffusion models can be decoded into physical state variables via 1x1 convolutions with sufficient fidelity for residual computation
    This premise is required for the PDE losses on intermediate features to be meaningful.
invented entities (1)
  • REPA-P projection heads no independent evidence
    purpose: Decode selected layer activations into physical quantities for residual supervision
    New lightweight components introduced by the method and discarded after training.

pith-pipeline@v0.9.0 · 5762 in / 1341 out tokens · 32836 ms · 2026-05-21T06:20:13.309896+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages · 14 internal anchors

  1. [1]

    Chen, W., Jia, H., Lai, S., Wu, K., Xiao, H., Hu, L., and Yue, Y

    URL https:// arxiv.org/abs/2403.14404. Chen, W., Jia, H., Lai, S., Wu, K., Xiao, H., Hu, L., and Yue, Y . Free-T2M: Frequency enhanced text-to-motion diffusion model with consistency loss, 2025a. Chen, W., Li, H., Liang, S., Wang, L., Jia, H., Yuan, K., Wu, J., Tian, B., and Yue, Y . POLARIS: Projection-orthogonal least squares for robust and adaptive inv...

    work page arXiv 2026

  2. [2]

    Diffusion Posterior Sampling for General Noisy Inverse Problems

    URL https://arxiv.org/ abs/2209.14687. Cuomo, S., di Cola, V . S., Giampaolo, F., Rozza, G., Raissi, M., and Piccialli, F. Scientific machine learn- ing through physics-informed neural networks: Where we are and what’s next,

    work page internal anchor Pith review Pith/arXiv arXiv

  3. [3]

    S., Giampaolo, F., Rozza, G., Raissi, M., et al

    URL https://arxiv. org/abs/2201.05624. Esser, P., Kulal, S., Blattmann, A., Entezari, R., M ¨uller, J., Saini, H., Levi, Y ., Lorenz, D., Sauer, A., Boesel, F., Podell, D., Dockhorn, T., English, Z., Lacey, K., Goodwin, A., Marek, Y ., and Rombach, R. Scaling recti- fied flow transformers for high-resolution image synthe- sis,

    work page arXiv

  4. [4]

    URL https: //arxiv.org/abs/2312.06662. Ho, J. and Salimans, T. Classifier-free diffusion guid- ance,

    work page arXiv

  5. [5]

    Classifier-Free Diffusion Guidance

    URL https://arxiv.org/abs/ 2207.12598. Ho, J., Jain, A., and Abbeel, P. Denoising diffusion prob- abilistic models,

    work page internal anchor Pith review Pith/arXiv arXiv

  6. [6]

    Denoising Diffusion Probabilistic Models

    URL https://arxiv.org/ abs/2006.11239. Jia, H., Chen, W., Huang, Z., Wang, L., Xiao, H., Jia, N., Wu, K., Lai, S., Tian, B., and Yue, Y . Physics-informed representation alignment for sparse radio-map reconstruc- tion. InProceedings of the 33rd ACM International Con- ference on Multimedia, MM ’25, pp. 12352–12360. As- sociation for Computing Machinery, 20...

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  7. [7]

    & Song, J

    URL https:// arxiv.org/abs/2201.11793. Li, H., Chen, W., Liang, S., Wang, L., Jia, H., and Yue, Y . Guided path sampling: Steering diffusion models back on track with principled path guidance,

    work page arXiv

  8. [8]

    Fourier Neural Operator for Parametric Partial Differential Equations

    URL https://arxiv.org/abs/ 2010.08895. 10 Learning to Think in Physics Li, Z., Zheng, H., Kovachki, N., Jin, D., Chen, H., Liu, B., Azizzadenesheli, K., and Anandkumar, A. Physics- informed neural operator for learning partial differential equations,

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  9. [9]

    Physics-informed neural operator for learning partial differential equations, 2024

    URL https://arxiv.org/abs/ 2111.03794. Lu, C., Zhou, Y ., Bao, F., Chen, J., Li, C., and Zhu, J. Dpm-solver: A fast ode solver for diffusion probabilistic model sampling in around 10 steps,

    work page arXiv

  10. [10]

    Dpm-solver: A fast ode solver for diffusion probabilistic model sampling in around 10 step s

    URL https: //arxiv.org/abs/2206.00927. Maz´e, F. and Ahmed, F. Diffusion models beat gans on topology optimization,

    work page arXiv

  11. [11]

    org/abs/2208.09591

    URL https://arxiv. org/abs/2208.09591. Nichol, A. and Dhariwal, P. Improved denoising diffusion probabilistic models,

    work page arXiv

  12. [12]

    Improved Denoising Diffusion Probabilistic Models

    URL https://arxiv. org/abs/2102.09672. Ning, M., Li, M., Su, J., Jia, H., Liu, L., Benes, M., Chen, W., Salah, A. A., and Onal Ertugrul, I. DCTdiff: Intrigu- ing properties of image generative modeling in the DCT space. InProceedings of the 42nd International Confer- ence on Machine Learning, volume 267 ofProceedings of Machine Learning Research, pp. 4649...

    work page internal anchor Pith review Pith/arXiv arXiv

  13. [13]

    Scalable Diffusion Models with Transformers

    URL https://arxiv.org/abs/ 2212.09748. Rombach, R., Blattmann, A., Lorenz, D., Esser, P., and Ommer, B. High-resolution image synthesis with latent diffusion models,

    work page internal anchor Pith review Pith/arXiv arXiv

  14. [14]

    High-Resolution Image Synthesis with Latent Diffusion Models

    URL https://arxiv.org/ abs/2112.10752. Ronneberger, O., Fischer, P., and Brox, T. U-net: Convolu- tional networks for biomedical image segmentation,

    work page internal anchor Pith review Pith/arXiv arXiv

  15. [15]

    U-Net: Convolutional Networks for Biomedical Image Segmentation

    URLhttps://arxiv.org/abs/1505.04597. Salimans, T. and Ho, J. Progressive distillation for fast sampling of diffusion models,

    work page internal anchor Pith review Pith/arXiv arXiv

  16. [16]

    Progressive Distillation for Fast Sampling of Diffusion Models

    URL https:// arxiv.org/abs/2202.00512. Sohl-Dickstein, J., Weiss, E. A., Maheswaranathan, N., and Ganguli, S. Deep unsupervised learning using nonequilibrium thermodynamics,

    work page internal anchor Pith review Pith/arXiv arXiv

  17. [17]

    Deep Unsupervised Learning using Nonequilibrium Thermodynamics

    URL https: //arxiv.org/abs/1503.03585. Song, J., Meng, C., and Ermon, S. Denoising diffusion implicit models.CoRR, abs/2010.02502,

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  18. [18]

    Denoising Diffusion Implicit Models

    URL https://arxiv.org/abs/2010.02502. Song, Y ., Sohl-Dickstein, J., Kingma, D. P., Kumar, A., Ermon, S., and Poole, B. Score-based generative model- ing through stochastic differential equations,

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  19. [19]

    Score-Based Generative Modeling through Stochastic Differential Equations

    URL https://arxiv.org/abs/2011.13456. Song, Y ., Dhariwal, P., Chen, M., and Sutskever, I. Con- sistency models,

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  20. [20]

    Consistency Models

    URL https://arxiv.org/ abs/2303.01469. Wang, H., Han, J., Fan, W., Zhang, W., and Liu, H. Phyda: Physics-guided diffusion models for data assimilation in atmospheric systems,

    work page internal anchor Pith review Pith/arXiv arXiv

  21. [21]

    org/abs/2505.12882

    URL https://arxiv. org/abs/2505.12882. Zhang, L., Rao, A., and Agrawala, M. Adding conditional control to text-to-image diffusion models,

    work page arXiv

  22. [22]

    Adding Conditional Control to Text-to-Image Diffusion Models

    URL https://arxiv.org/abs/2302.05543. Zheng, H., Nie, W., Vahdat, A., Azizzadenesheli, K., and Anandkumar, A. Fast sampling of diffusion models via operator learning,

    work page internal anchor Pith review Pith/arXiv arXiv

  23. [23]

    URL https://arxiv.org/ abs/2211.13449. Zhu, Y . and Zabaras, N. Bayesian deep convolutional encoder–decoder networks for surrogate modeling and uncertainty quantification.Journal of Computational Physics, 366:415–447, August

    work page arXiv

  24. [24]

    doi: 10.1016/j.jcp.2018.04.018

    ISSN 0021-9991. doi: 10.1016/j.jcp.2018.04.018. URL http://dx. doi.org/10.1016/j.jcp.2018.04.018. 11 Learning to Think in Physics A. Experimental Details This section provides additional details on the experimental setup, residual computation, and evaluation metrics for each benchmark task. A.1. Darcy Flow Problem Formulation.We study steady two-dimension...

    work page doi:10.1016/j.jcp.2018.04.018 2018

  25. [25]

    Each head consists of a 1×1 convolution mapping from the hidden dimension to 2 output channels (for K and p), followed by bilinear upsampling to the target resolution 64×64

    For REPA-P, we attach 1×1 convolutional projection heads to the bottleneck and selected decoder blocks. Each head consists of a 1×1 convolution mapping from the hidden dimension to 2 output channels (for K and p), followed by bilinear upsampling to the target resolution 64×64 . The mid-layer alignment weight is set to cmid = 0.1 for the main results, whic...

    work page 2025

  26. [26]

    Best results (excluding baseline) inbold.↓: lower is better;↑: higher is better. Darcy Flow Topology Optimization (ID) Charge Turbulence cmid Data↓Phys.↓Phys.↓CE%↓VFE%↓Phys.↓PSNR↑Phys.↓ Baseline 0.0180 0.0260 5.2e-3 9.24 3.38 0.381 37.64 1.91e-3 0.001 0.0156 0.0177 7.8e-3 11.23 4.13 0.245 38.47 1.86e-3 0.005 0.0142 0.01654.5e-3 4.17 3.020.189 37.91 2.04e-...

    work page arXiv