Towards the simulation of transition-metal oxides of the cathode battery materials using VQE methods

Variational quantum eigensolver (VQE) is a hybrid quantum-classical technique that leverages noisy intermediate scale quantum (NISQ) hardware to obtain the minimum eigenvalue of a model Hamiltonian. VQE has so far been used to simulate condensed matter systems as well as quantum chemistry of small molecules. In this work, we employ VQE methods to obtain the ground-state energy of LiCoO$_2$, a candidate transition metal oxide used for battery cathodes. We simulate Li$_2$Co$_2$O$_4$ and Co$_2$O$_4$ gas-phase models, which represent the lithiated and delithiated states during the discharge and the charge of the Li-ion battery, respectively. Computations are performed using a statevector simulator with a single reference state for three different trial wavefunctions: unitary coupled-cluster singles and doubles (UCCSD), unitary coupled-cluster generalized singles and doubles (UCCGSD) and k-unitary pair coupled-cluster generalized singles and doubles (k-UpCCGSD). The resources in terms of circuit depth, two-qubit entangling gates and wavefunction parameters are analyzed. We find that the k-UpCCGSD with k=5 produces results similar to UCCSD but at a lower cost. Finally, the performance of VQE methods is benchmarked against the classical wavefunction-based methods, such as coupled-cluster singles and doubles (CCSD) and complete active space configuration interaction (CASCI). Our results show that VQE methods quantitatively agree with the results obtained from CCSD. However, the comparison against the CASCI results clearly suggests that advanced trial wavefunctions are likely necessary to capture the multi-reference characteristics as well as the correlations emerging from high-level electronic excitations.