diff --git a/MaterialLib/Adsorption/Adsorption.cpp b/MaterialLib/Adsorption/Adsorption.cpp
index 7a2f2624d1e8a17999b0d5634e89bdfb4ed3af91..a5f257acd0c548445d71bc2569dbf866d94e6434 100644
--- a/MaterialLib/Adsorption/Adsorption.cpp
+++ b/MaterialLib/Adsorption/Adsorption.cpp
@@ -118,29 +118,6 @@ double AdsorptionReaction::getReactionRate(const double p_Ads, const double T_Ad
                                       // this the rate in terms of loading!
 }
 
-void AdsorptionReaction::getDReactionRate(const double p_Ads, const double T_Ads,
-                                     const double M_Ads, const double /*loading*/,
-                                     std::array<double, 3> &dqdr) const
-{
-    const double A = getPotential(p_Ads, T_Ads, M_Ads);
-    const double p_S = getEquilibriumVapourPressure(T_Ads);
-    const double dAdT = MaterialLib::PhysicalConstant::IdealGasConstant * log(p_S/p_Ads) / (M_Ads*1.e3);
-    const double dAdp = - MaterialLib::PhysicalConstant::IdealGasConstant * T_Ads / M_Ads / p_Ads;
-
-    const double W = characteristicCurve(A);
-    const double dWdA = dCharacteristicCurve(A);
-
-    const double rho_Ads = getAdsorbateDensity(T_Ads);
-    const double drhodT = - rho_Ads * getAlphaT(T_Ads);
-
-    dqdr = std::array<double, 3>{{
-        rho_Ads*dWdA*dAdp,
-        drhodT*W + rho_Ads*dWdA*dAdT,
-        -k_rate
-    }};
-}
-
-
 // Evaluate adsorbtion potential A
 double AdsorptionReaction::getPotential(const double p_Ads, double T_Ads, const double M_Ads) const
 {
diff --git a/MaterialLib/Adsorption/Adsorption.h b/MaterialLib/Adsorption/Adsorption.h
index 7f2c51739f1ce9c8168a2bd856acf0d92a1ab5c9..8f139642f51d3b878f335dac3035af84d053f2ae 100644
--- a/MaterialLib/Adsorption/Adsorption.h
+++ b/MaterialLib/Adsorption/Adsorption.h
@@ -39,17 +39,6 @@ public:
     double getReactionRate(const double p_Ads, const double T_Ads,
                            const double M_Ads,
                            const double loading) const override;
-    /**
-     * @brief get_d_reaction_rate
-     * @param p_Ads
-     * @param T_Ads
-     * @param M_Ads
-     * @param loading
-     * @param dqdr array containing the differentials wrt: p, T, C
-     */
-    virtual void getDReactionRate(const double p_Ads, const double T_Ads,
-                                     const double M_Ads, const double loading,
-                                     std::array<double, 3>& dqdr) const;
 
 protected:
     virtual double getAdsorbateDensity(const double T_Ads) const = 0;