ReadingResearchRadarInvestment frameworkAI Company Primary Sources
Sign in / Sign up
← Back to reading
NewsQuick take2 min readDisconfirmAI Editorial

NVIDIA Says TF32x3 Accelerates Machine Learning Interatomic Potentials While Preserving Energy Conservation, Single-Pass TF32 Fails to Converge

Machine learning interatomic potentials are sensitive to numerical errors in forces, leading most to stick with FP32 inference. NVIDIA reported in a technical blog that BF16x9 and TF32x3—two Tensor Core methods for simulating FP32—can accelerate PET molecular dynamics steps by 1.1–1.7x and approximately 1.5x, respectively, without retraining the model. TF32x3 maintains physical behavior consistent with FP32 in energy conservation and structural relaxation convergence tests, whereas the faster single-pass TF32 fails to meet force convergence standards.

2026-10-29 23:00 GMT+8

Machine learning interatomic potentials (MLIPs) use neural networks to predict atomic energies and forces at a fraction of the cost of quantum mechanical methods like density functional theory, but their precision requirements remain high. Forces are derivatives of energy with respect to atomic positions; in molecular dynamics simulations, forces at every step determine updates to atomic velocities and positions, meaning numerical errors accumulate along the trajectory. Therefore, even though GPUs offer much higher throughput for low-precision matrix multiplication than for FP32, many MLIPs continue to rely on FP32 inference.

In a technical blog post, NVIDIA evaluated two Tensor Core FP32 simulation methods. On Blackwell GPUs, BF16x9 decomposes each FP32 operand into three scaled BF16 components, computes nine pairwise products, and accumulates them in FP32; it can be enabled via the environment variable CUBLAS_EMULATE_SINGLE_PRECISION=1 without modifying model code. On Ampere GPUs, TF32x3 splits each FP32 operand into a leading TF32 component and a residual TF32 component, computing three cross-products (leading × leading, leading × residual, residual × leading) while omitting the smallest residual × residual term. Implemented via CUTLASS's OpMultiplyAddFastF32 operator, inputs and outputs remain in FP32, making the method transparent to the rest of the model.

Sign in / Sign up
ReadingResearchRadarInvestment frameworkAI Company Primary Sources

The real difference lies not in point precision but in physical consistency. The NVIDIA team validated this using three tests: when dragging a CO molecule away from a Pt(111) surface, TF32x3 followed the smooth energy curve of FP32, while single-pass TF32 introduced visible point-to-point fluctuations; in an NVE simulation of a 192-atom water box, TF32x3 showed bounded energy fluctuations comparable to FP32, whereas single-pass TF32 exhibited systematic drift; during structural relaxation of a 539-atom Cantor alloy slab, TF32x3 reached the force convergence standard of 1.0×10⁻³ eV/Å within 218–222 steps (comparable to FP32's 218 steps), while single-pass TF32 failed to converge even after doubling the step limit to 1,000 steps, with forces plateauing around 1.2×10⁻³ eV/Å.

Regarding performance, BF16x9 achieved an end-to-end speedup of 1.7x on a 17,280-atom system (1.1–1.2x on a 192-atom system), with errors falling below the model's own non-determinism levels; TF32x3 accelerated MD steps on A100 GPUs by about 1.5x, with force errors roughly three times the differences between FP32 recompilations. However, the magnitude of acceleration depends on system size and the proportion of total runtime spent on GEMM operations—smaller systems see more limited gains due to higher startup overhead and non-GEMM workloads. TF32x3 performance also relies on kernel tuning specific to matrix shapes and hardware. Neither method requires retraining or changes to the model architecture, allowing them to serve as independent performance layers that can be combined with other software optimizations.

Sources

  • developer.nvidia.comNVIDIA Technical Blog↗
Continue researching →