Quantum Computations on Fusion Blanket Molten Salts

Abstract

Molten salts such as FLiBe (2LiF--BeF$_2$) are leading blanket materials for breeding and recovering tritium in fusion reactors. Predicting tritium speciation requires accurate electronic ground-state energies for representative molten-salt clusters, a demanding task for correlated electronic-structure methods. Here we report the first application of heterogeneous quantum--classical computing to tritium binding in FLiBe. Clusters drawn from ab initio molecular dynamics are partitioned by an embedded-wavefunction (EWF) method into atom-centered fragments, and the largest fragments are solved on IBM quantum hardware using extended sample-based quantum diagonalization (ext-SQD). Across nine clusters, the heterogeneous quantum--classical workflow reproduces fragment ground-state energies with agreement to full configuration interaction within 0.7~kcal/mol and a mean absolute deviation of 0.3~kcal/mol. In contrast, fragmented and unfragmented conformational energy differences and tritium binding energies differ by 12~kcal/mol and 110~kcal/mol on average, respectively, identifying fragment construction rather than fragment solution as the dominant source of algorithmic bias. To the best of our knowledge, this is the first such demonstration for a charged ionic system and in particular an inorganic molten salt, where electrostatic and polarization effects make the accurate treatment of electronic correlation particularly challenging. These results also identify areas of future research towards an accurate and scalable quantum--classical workflow to compute free-energy estimates of tritium speciation in fusion blankets.