--- branches/development/src/nonbonded/LJ.cpp	2011/09/14 21:15:17	1629
+++ branches/development/src/nonbonded/LJ.cpp	2012/02/29 20:33:01	1683
@@ -36,7 +36,8 @@
  * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).             
  * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).          
  * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).          
- * [4]  Vardeman & Gezelter, in progress (2009).                        
+ * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
+ * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
  */
 
 #include <stdio.h>
@@ -45,6 +46,7 @@
 #include <cmath>
 #include "nonbonded/LJ.hpp"
 #include "utils/simError.h"
+#include "types/LennardJonesInteractionType.hpp"
 
 namespace OpenMD {
 
@@ -123,51 +125,40 @@ namespace OpenMD {
 
     for (at = atomTypes->beginType(i); at != NULL; 
          at = atomTypes->nextType(i)) {
-      
-      if (at->isLennardJones())
-        addType(at);
+      if (at->isLennardJones()){
+	 addType(at);
+      }
     }
-
     ForceField::NonBondedInteractionTypeContainer* nbiTypes = forceField_->getNonBondedInteractionTypes();
     ForceField::NonBondedInteractionTypeContainer::MapTypeIterator j;
     NonBondedInteractionType* nbt;
+    ForceField::NonBondedInteractionTypeContainer::KeyType keys;
 
+
     for (nbt = nbiTypes->beginType(j); nbt != NULL; 
          nbt = nbiTypes->nextType(j)) {
-      
+
       if (nbt->isLennardJones()) {
-        
-        pair<AtomType*, AtomType*> atypes = nbt->getAtomTypes();
-        
-        GenericData* data = nbt->getPropertyByName("LennardJones");
-        if (data == NULL) {
-          sprintf( painCave.errMsg, "LJ::rebuildMixingMap could not find\n"
-               "\tLennard-Jones parameters for %s - %s interaction.\n", 
-                   atypes.first->getName().c_str(),
-                   atypes.second->getName().c_str());
-          painCave.severity = OPENMD_ERROR;
-          painCave.isFatal = 1;
-          simError(); 
-        }
-    
-        LJParamGenericData* ljData = dynamic_cast<LJParamGenericData*>(data);
-        if (ljData == NULL) {
+        keys = nbiTypes->getKeys(j);
+        AtomType* at1 = forceField_->getAtomType(keys[0]);
+        AtomType* at2 = forceField_->getAtomType(keys[1]);
+
+        LennardJonesInteractionType* ljit = dynamic_cast<LennardJonesInteractionType*>(nbt);
+
+        if (ljit == NULL) {
           sprintf( painCave.errMsg,
-                   "LJ::rebuildMixingMap could not convert GenericData to\n"
-                   "\tLJParam for %s - %s interaction.\n", 
-                   atypes.first->getName().c_str(),
-                   atypes.second->getName().c_str());
+                   "LJ::initialize could not convert NonBondedInteractionType\n"
+                   "\tto LennardJonesInteractionType for %s - %s interaction.\n", 
+                   at1->getName().c_str(),
+                   at2->getName().c_str());
           painCave.severity = OPENMD_ERROR;
           painCave.isFatal = 1;
           simError();          
         }
-        
-        LJParam ljParam = ljData->getData();
 
-        RealType sigma = ljParam.sigma;
-        RealType epsilon = ljParam.epsilon;
-
-        addExplicitInteraction(atypes.first, atypes.second, sigma, epsilon);
+        RealType sigma = ljit->getSigma();
+        RealType epsilon = ljit->getEpsilon();
+        addExplicitInteraction(at1, at2, sigma, epsilon);
       }
     }  
     initialized_ = true;
@@ -219,10 +210,6 @@ namespace OpenMD {
   
   void LJ::addExplicitInteraction(AtomType* atype1, AtomType* atype2, RealType sigma, RealType epsilon){
     
-    // in case these weren't already in the map
-    addType(atype1);
-    addType(atype2);
-
     LJInteractionData mixer;
     mixer.sigma = sigma;
     mixer.epsilon = epsilon;
@@ -240,9 +227,7 @@ namespace OpenMD {
   }
  
   void LJ::calcForce(InteractionData &idat) {
-    
     if (!initialized_) initialize();
-
     map<pair<AtomType*, AtomType*>, LJInteractionData>::iterator it;
     it = MixingMap.find( idat.atypes );