Simulating charging characteristics of lithium iron phosphate by electro-ionic optimization on a quantum annealer

The rapid evolution of quantum computing hardware opens up new avenues in the simulation of energy materials. Today's quantum annealers are able to tackle complex combinatorial optimization problems. A formidable challenge of this type is posed by materials with site-occupational disorder for which atomic arrangements with a low, or lowest, energy must be found. In this article, a method is presented for the identification of the correlated ground-state distribution of both lithium ions and redox electrons in lithium iron phosphate (LFP), a widely employed cathode material in lithium-ion batteries. The point-charge Coulomb energy model employed correctly reproduces the LFP charging characteristics. As is shown, grand-canonical transformation of the energy cost function makes the combinatorial distribution problem solvable on quantum annealing (QA) hardware. The QA output statistics follow a pseudothermal behavior characterized by a problem-dependent effective sampling temperature, which has bearings on the estimated scaling of the QA performance with system size. This work demonstrates the potential of quantum computation for the joint simulation of the electronic and ionic structure in energy materials.