diff --git a/include/cantera/thermo/LatticeSolidPhase.h b/include/cantera/thermo/LatticeSolidPhase.h
index 754046bc3..d8e4bc637 100644
--- a/include/cantera/thermo/LatticeSolidPhase.h
+++ b/include/cantera/thermo/LatticeSolidPhase.h
@@ -643,17 +643,6 @@ public:
      */
     void setLatticeMoleFractionsByName(int n, std::string x);
 
-    //! Return a changeable reference to the calculation manager
-    //! for species reference-state thermodynamic properties
-    /*!
-     *  This routine returns the calculation manager for the sublattice
-     *
-     * @param k   Species id. The default is -1, meaning return the default
-     *
-     * @internal
-     */
-    virtual SpeciesThermo& speciesThermo(int k = -1);
-
 #ifdef H298MODIFY_CAPABILITY
 
     //! Modify the value of the 298 K Heat of Formation of one species in the phase (J kmol-1)
diff --git a/src/thermo/LatticeSolidPhase.cpp b/src/thermo/LatticeSolidPhase.cpp
index 6473b767c..30557d036 100644
--- a/src/thermo/LatticeSolidPhase.cpp
+++ b/src/thermo/LatticeSolidPhase.cpp
@@ -668,31 +668,6 @@ void LatticeSolidPhase::setParametersFromXML(const XML_Node& eosdata)
 
 }
 //====================================================================================================================
-// Return a changeable reference to the calculation manager
-// for species reference-state thermodynamic properties
-/*
- *
- * @param k   Species id. The default is -1, meaning return the default
- *
- * @internal
- */
-SpeciesThermo& LatticeSolidPhase::speciesThermo(int k)
-{
-    return *m_spthermo;
-    /*
-    int kk;
-    if (k >= 0) {
-      for (int n = 0; n < m_nlattice; n++) {
-    if (lkstart_[n+1] < k) {
-      kk = k - lkstart_[n];
-      return m_lattice[n]->speciesThermo(kk);
-    }
-      }
-    }
-    return m_lattice[0]->speciesThermo(-1);
-    */
-}
-//====================================================================================================================
 
 #ifdef H298MODIFY_CAPABILITY