Bug fixes from:

http://www.cfd-online.com/Forums/openfoam-bugs/64379-bugs-moleculedynamics-
model.html

http://www.cfd-online.com/Forums/openfoam-bugs/64383-fix-taking-reference-
temporary-object-tetherpotentiallist-c.html

src/lagrangian/molecularDynamics/molecule/interactionLists/directInteractionList/directInteractionList.C

@@ -119,8 +119,6 @@ void Foam::directInteractionList::buildDirectInteractionList
 
                     const label cellO(mesh.faceOwner()[f]);
 
-                    const label cellN(mesh.faceNeighbour()[f]);
-
                     forAll(pCells, pC)
                     {
                         const label cellI(pCells[pC]);
@@ -156,6 +154,8 @@ void Foam::directInteractionList::buildDirectInteractionList
                             // boundary faces will not have neighbour
                             // information
 
+                            const label cellN(mesh.faceNeighbour()[f]);
+
                             if (cellN > cellI)
                             {
                                 if

src/lagrangian/molecularDynamics/potential/pairPotential/pairPotentialList/pairPotentialList.C

@@ -225,7 +225,7 @@ const Foam::pairPotential& Foam::pairPotentialList::pairPotentialFunction
 
 bool Foam::pairPotentialList::rCutMaxSqr(const scalar rIJMagSqr) const
 {
-    if (rIJMagSqr <= rCutMaxSqr_)
+    if (rIJMagSqr < rCutMaxSqr_)
     {
         return true;
     }
@@ -243,7 +243,7 @@ bool Foam::pairPotentialList::rCutSqr
     const scalar rIJMagSqr
 ) const
 {
-    if (rIJMagSqr <= rCutSqr (a, b))
+    if (rIJMagSqr < rCutSqr (a, b))
     {
         return true;
     }