Pre-training neural operators to solve partial differential equations (PDEs) is a frontier research direction: a single pre-trained model can be applied across many different PDE types, avoiding the cost of retraining from scratch for each equation. However, as the variety of PDEs grows, a prominent problem emerges — the solution operators of different PDEs differ dramatically in structure. Diffusion equations produce smooth solutions, advection equations produce sharp ones, and reaction-diffusion equations exhibit complex multi-scale behavior. Existing approaches force ever-larger models to cover this diversity.
In *AOT-POT: Adaptive Operator Transformation for Large-Scale PDE Pre-training*, Lv, Wang, Hao, Wu, Xu, Zhou, Wu, and Zhang propose the opposite idea: rather than enlarging the model, transform the operator. If the solution operator is transformed into a simpler, more uniform form before training, a medium-sized model can cover an entire family of operators.
Method
The key requirement is that the transformation must be adaptive and input-dependent: the same PDE can have differently shaped solution operators under different parameters. The transformation mechanism expands the hidden representation into multiple parallel streams and adaptively aggregates and redistributes the information flow before and after each sublayer. Mixing between streams is achieved through doubly stochastic matrices with Sinkhorn projection.
The authors explicitly connect this to classical numerical analysis: preconditioners have been used in traditional PDE solvers for decades. The contribution is transplanting the same idea into neural networks — replacing "carefully hand-crafted preconditioning for each PDE" with "letting the network learn how to precondition itself."
Results
- Across 12 PDE benchmarks, AOT-POT achieves the best performance with only 3% additional parameters.
- Relative L2 error is reduced by up to 77.6% (average 40.9%).
- After fine-tuning, error drops by 92% on in-domain PDEs and 89% on out-of-domain types unseen during pre-training.
- What are the stability conditions for Sinkhorn-projected doubly stochastic matrices during training — how many iterations are needed to converge?
- How should the number of parallel streams be chosen — does it need to scale with the number of PDE types?
- Does the model still benefit from the transformation on very unusual PDE types, or would more aggressive transformation strategies be required?
Open Questions
References
1. Lv, Q., Wang, H., Hao, Z., et al. (2026). *AOT-POT: Adaptive Operator Transformation for Large-Scale PDE Pre-training*. arXiv:2605.15793 [cs.LG]. 2. Lu, L., et al. (2021). *DeepONet: Learning Nonlinear Operators for Identifying Differential Equations Based on the Universal Approximation Theorem of Operators*. Nature Machine Intelligence. 3. Li, Z., et al. (2021). *Fourier Neural Operator for Parametric Partial Differential Equations*. ICLR.