diff --git a/src/FMT.cpp b/src/FMT.cpp
index b6484e8..7a773ec 100644
--- a/src/FMT.cpp
+++ b/src/FMT.cpp
@@ -545,17 +545,6 @@ void FMT::add_second_derivative(const vector<DFT_FFT> &v, vector<DFT_Vec> &d2F,
 	} 
            
     }
-  // Now get EOS contributions
-  for(int s = 0;s<Nspecies;s++)
-    {
-      FMT_Species_EOS *species = dynamic_cast<FMT_Species_EOS*>(allSpecies[s]);
-      if(species)
-	species->add_second_derivative(v[s], d2F[s], this);
-    }
-
-
-
-  
 }
 
 // Adds contribution to F_{I,I+J} = sum_{a,b} sum_K dV Phi_{ab}(K) w_a(K-I) w_b(K-I-J)
diff --git a/src/FMT_Species_EOS.cpp b/src/FMT_Species_EOS.cpp
index f4f8ab8..76f5b9d 100644
--- a/src/FMT_Species_EOS.cpp
+++ b/src/FMT_Species_EOS.cpp
@@ -187,5 +187,8 @@ void FMT_Species_EOS::add_second_derivative(const DFT_FFT &v, DFT_Vec &d2F, cons
   //the fmt weights have all necessary normalization factors built in, including dV
   bConjugate = true;
   convolute_eos_eta_weight_with(Lambda, result, bConjugate);
+  cout << "d2F.get(12) = " << d2F.get(12) << endl;
+  cout << "result.Real().get(12) = " << result.Real().get(12) << endl;
   d2F.IncrementBy_Scaled_Vector(result.Real(),dV);
+  cout << "d2F.get(12) = " << d2F.get(12) << endl;  
 }