diff --git a/include/cantera/thermo/SurfPhase.h b/include/cantera/thermo/SurfPhase.h
index 22b43a5f7..96e6b176b 100644
--- a/include/cantera/thermo/SurfPhase.h
+++ b/include/cantera/thermo/SurfPhase.h
@@ -493,6 +493,8 @@ public:
         m_press = p;
     }
 
+    virtual void getPureGibbs(doublereal* g) const;
+
     //!  Returns the vector of nondimensional
     //!  Gibbs Free Energies of the reference state at the current temperature
     //!  of the solution and the reference pressure for the species.
diff --git a/src/thermo/SurfPhase.cpp b/src/thermo/SurfPhase.cpp
index 4f30ba6bb..67f7bceba 100644
--- a/src/thermo/SurfPhase.cpp
+++ b/src/thermo/SurfPhase.cpp
@@ -186,6 +186,12 @@ void SurfPhase::setParameters(int n, doublereal* const c)
     setSiteDensity(c[0]);
 }
 
+void SurfPhase::getPureGibbs(doublereal* g) const
+{
+    _updateThermo();
+    copy(m_mu0.begin(), m_mu0.end(), g);
+}
+
 void SurfPhase::getGibbs_RT(doublereal* grt) const
 {
     _updateThermo();