where µ ref O is the reference oxygen chemical potential, such that the first and last values in the sequence correspond to relaxed vacancy formation energies. Therefore, ∆H p0 = ∆H s and ∆H pn+1 = ∆H e and are path independent. The forward and reverse migration activation energies can then be computed,