ML: Diffusion sampling

Machine-learning ideas tagged Diffusion sampling in the ML taxonomy of the Math2NN corpus.

✗ Failed on benchmark 2026

Endpoint-Jacobian diffusion backpropagation

Group W consecutive diffusion or flow-model loss terms and approximate every intermediate parameter Jacobian by a time-weighted interpolation of the Jacobians at the group’s two endpoints. Sum the intermediate upstream signals into two endpoint cotangents, then perform only two full DiT backward passes instead of W. Add a cosine-similarity gate comparing predicted and actual intermediate velocity changes so that groups violating the local-linearity assumption use exact backpropagation.

Useful8/10
Difficulty6/10
Novelty7/10
Paper: JAGG: Jacobian-Aggregated Group Gradient for Efficient GRPO Training of Diffusion Models arXiv:2607.17572
△ Mechanism confirmed, baseline not beaten 2026

Pole-residue latent diffusion decoder

Represent a one-dimensional sharp-feature signal by a small unordered set of complex singularities and residues instead of predicting all grid amplitudes. A transformer diffusion model predicts these tokens, and a differentiable meromorphic decoder evaluates the result directly at arbitrary coordinates, avoiding grid-specific interpolation and preserving discontinuity structure.

Useful8/10
Difficulty6/10
Novelty8/10
Paper: The Singularity Space: A Generative Diffusion Framework for Signal Representation arXiv:2607.10930
✗ Mechanism failed 2026

Recorded-Mesh Neural ODE Backpropagation

Run an adaptive neural ODE solver once to determine accepted step sizes, then train using a regular fixed-length replay of those steps rather than differentiating through adaptive accept/reject logic. The replay can be fused across a batch of trajectories and differentiated with an ordinary reverse sweep, giving the exact discrete gradient of the replayed solver and predictable GPU work.

Useful7/10
Difficulty5/10
Novelty6/10
Paper: GRADSOLVE: fast exact gradients for ODE ensembles on GPUs arXiv:2609.02876
△ Mechanism confirmed, baseline not beaten 2026

Dilation-Matched Metropolized Dynamics

Replace the unstable classical derivative of a discretized rough energy component with a matched dilation quotient derived from its intrinsic scale recursion. Use this field inside kick-drift-kick proposals and apply an exact Metropolis correction, allowing the proposal field to be measurable and nonconservative rather than an exact neural-energy gradient. The experiment should test whether acceptance rates and posterior samples remain stable as the rough-energy resolution increases.

Useful7/10
Difficulty5/10
Novelty8/10
Paper: Posterior Convergence without Force Convergence: Resolution-Stable Sampling for Rough Bayesian Inverse Problems arXiv:2608.18365
△ Mechanism confirmed, baseline not beaten 2026

Tau-leaped parallel discrete Hamiltonian sampler

Approximate the exact event-by-event lifted sampler by drawing independent Poisson jump counts over a short interval and applying compatible discrete moves in parallel. This converts sequential neighbor events into batched GPU-friendly updates while retaining the Hamiltonian rate structure; the step size controls the error-versus-throughput tradeoff.

Useful7/10
Difficulty5/10
Novelty8/10
Paper: Hamiltonian dynamics for sampling on discrete spaces arXiv:2608.17961
✓✓ Beats tuned baseline 2026

Context-free denoiser with analytic quadratic score injection

Train one denoiser only for the nonquadratic residual distribution, then modify the diffusion sampler using an analytically computed quadratic Gaussian context. Changing $K$ at inference changes the target distribution without retraining the denoiser, enabling transfer across temperatures, masses, coupling strengths, and boundary conditions whenever those changes remain quadratic.

Useful7/10
Difficulty6/10
Novelty6/10
Paper: Nuclear Quantum Effects as a Denoising Problem arXiv:2607.19680
Unverified 2026

KL-optimal joint rare-event tilting

Modify a diffusion sampler over a complete scenario trajectory z by exponentially tilting its prior toward a learned failure score s(z), rather than conditioning independently on environment and execution variables. The resulting sampler spends more evaluations in rare failure regions while preserving an explicit importance weight for estimating probabilities under the original distribution.

Useful6/10
Difficulty5/10
Novelty4/10
Paper: Diffusion-Guided Search via Exponential Tilting (DiffTilt): An Application to Falsification of Safety-Critical Systems arXiv:2607.23134
Unverified 2026

Hard Sequential Neural ODE Solver

Partition a long integration interval into M short segments and assign one neural trajectory approximator to each segment. Instead of asking a single network to satisfy the ODE and initial condition over the entire horizon, construct every segment so that its value at the left boundary is exactly the terminal value predicted by the previous segment. This removes interface discontinuities from the optimization problem and should improve long-horizon trajectory accuracy and gradient stability.

Useful6/10
Difficulty5/10
Novelty6/10
Paper: Modeling of an ODE-constrained optimization problem describing tumor dynamics, and numerical approximation via sequential physics-informed neural networks arXiv:2608.18974
Unverified 2026

Two-query SPSA Langevin sampler for black-box neural energies

Replace backpropagated gradients in a Langevin sampler with a simultaneous-perturbation estimate obtained from two noisy evaluations of a neural energy or simulator-defined negative log-density. This enables posterior or latent-space sampling when the energy contains nondifferentiable code, stochastic simulation, discrete operations, or an inaccessible neural-network gradient.

Useful6/10
Difficulty5/10
Novelty8/10
Paper: Zeroth-Order Langevin Monte Carlo via SPSA under Noisy Function Measurements arXiv:2608.07837
Unverified 2026

Hashed Local-Density Particle Layer

Replace an O(N^2) kernel-density interaction in a particle neural SDE or diffusion sampler with a clipped, randomly shifted histogram density estimate. Feed the local estimated density into the particle drift as a multiplicative gain, preserving density-dependent dynamics while evaluating all particles through occupied-cell hashing in expected O(N) time for fixed dimension and number of shifts.

Useful6/10
Difficulty5/10
Novelty8/10
Paper: Density-Dependent McKean--Vlasov Diffusions: Subgaussian Occupancy Bounds and Polynomial Propagation of Chaos arXiv:2607.19583
Unverified 2026

Resolvent response basis for conditioned modules

Compress a module whose output changes with a scalar condition such as diffusion time, temperature, or compute budget by representing its response in a low-rank basis generated by resolvent-like functions. Distinct spectral modes produce rational factors \((1-\tau\lambda_k)^{-1}\), allowing a small number of learned components to approximate a large hypernetwork or condition-dependent parameter table.

Useful5/10
Difficulty6/10
Novelty7/10
Paper: Spectral duality structures and the Fisher--Rao geometry of reset distributions arXiv:2608.15805