(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     40672,       1343]
NotebookOptionsPosition[     36503,       1214]
NotebookOutlinePosition[     36862,       1230]
CellTagsIndexPosition[     36819,       1227]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Atomic Clusters", "Title"],

Cell["\<\
APAM 1601
Columbia University\
\>", "Subsubtitle"],

Cell[CellGroupData[{

Cell["Introduction", "Section"],

Cell[TextData[{
 "In this notebook, we continue to explore the use of random walks and the \
Rosenbluth-Teller-Metropolis algorithm to understand complicated physical \
systems in \"thermal equilibrium\".  For this analysis, we will explore the \
three-dimensional structure of clusters of \"atoms\" interacting through the \
Lennard-Jones potential. ",
 ButtonBox["Sir John Edward Lennard-Jones (1894-1954) ",
  BaseStyle->"Hyperlink",
  ButtonData:>{
    URL["http://www.quantum-chemistry-history.com/Le-Jo1Ue.htm"], None}],
 "was a leading theoretical chemist of his time. His simple (yet highly \
nonlinear) model for the potential energy between two atoms or molecules has \
become a much studied system of molecular configuration. The potential energy \
between two atoms is a combination of an ",
 StyleBox["attractive",
  FontSlant->"Italic"],
 " van der Waals force (due to the electric polarization of atoms) and a ",
 StyleBox["repulsive",
  FontSlant->"Italic"],
 " term resulting from the overlap of electron orbits. Today's research \
trying to understand the configurations found in inorganic nanoscience and in \
protein folding continue to make these sorts of models an active area of \
research. "
}], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"lj", "[", "r_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"r", "^", "12"}]}], " ", "-", " ", 
   RowBox[{"2", " ", "/", 
    RowBox[{"r", "^", "6"}]}]}]}]], "Input",
 CellLabel->"In[1]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"lj2", "[", "r2_", "]"}], " ", ":=", " ", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"r2", "^", "6"}]}], " ", "-", " ", 
   RowBox[{"2", " ", "/", 
    RowBox[{"r2", "^", "3", " ", 
     RowBox[{"(*", " ", 
      StyleBox[
       RowBox[{
        RowBox[{"where", " ", "r2"}], " ", "=", " ", 
        RowBox[{"r", " ", ".", " ", "r"}]}],
       FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]}]}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lj", "[", "r", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "0.001`", ",", "2.`"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "1"}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<r\>\"", ",", "\"\<\[Phi](r)\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Lennard-Jones Potential\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[3]:="],

Cell[TextData[{
 "Atoms are attracted to each other, and a minimum energy is found when \
separated at ",
 StyleBox["r",
  FontSlant->"Italic"],
 "~ 1. When the atoms move closer together, a strong potential gradient \
separates them."
}], "Text"],

Cell[CellGroupData[{

Cell["Monte Carlo Simulation Rules", "Subsection"],

Cell[TextData[{
 "The \"rules\" for our numerical simulation are similar to those used for \
many problems in statistical physics (like an ideal ",
 StyleBox["non-interacting",
  FontSlant->"Italic"],
 " atmosphere)\[Dash]with one difference. For particles interacting by way of \
the Lennard-Jones potential, all particles are strongly ",
 StyleBox["interacting with each other",
  FontSlant->"Italic"],
 ". Many large problems in statistical physics assume that particles ",
 StyleBox["only interact with nearest neighbors. ",
  FontSlant->"Italic"],
 "When a particle moves, we need to compute the change in system energy by \
summing over all interactions."
}], "Text"],

Cell[TextData[{
 "We will follow the following set of rules:\nStep 1) Initialize the \
positions of the particles within a sphere of given radius..\nStep 2) Select \
a random particle. Make a potential random step  ",
 StyleBox["in any direction. However, ",
  FontSlant->"Italic"],
 "if the step causes the particle to escapee \nStep 3) Calculate the energy \
",
 StyleBox["change",
  FontSlant->"Italic"],
 " to the entire system if the particle makes the step.\nStep 4a) If the \
energy change is negative, accept the step.\nStep 4b) If the energy change is \
positive, pick a random number between 0 and 1. If the random number is less \
than ",
 StyleBox["Exp[ - (energy change)/kT]",
  FontWeight->"Bold"],
 ", then accept the step.\nStep 4c) Otherwise, reject the step and leave the \
position (and energy) unchanged.\nStep 5) Go to Step 2 and repeat (for a very \
long time....)"
}], "Text"],

Cell[TextData[{
 "These rules comprise the \"Metropolis\" method for computational \
statistical physics. It is problably the most often used computational method \
in physics! See online biography for ",
 ButtonBox["Nick Metropolis",
  BaseStyle->"Hyperlink",
  ButtonData:>{
    URL["http://scienceworld.wolfram.com/biography/Metropolis.html"], None}],
 ". The real insight behind the algorythm was developed by late ",
 ButtonBox["Marshall Rosenbluth",
  BaseStyle->"Hyperlink",
  ButtonData:>{
    URL["http://www.utexas.edu/faculty/council/2004-2005/memorials/rosenbluth/\
rosenbluth.html"], None}],
 ", the father of modern plasma physics. "
}], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"General", "::", "spell1"}], "]"}], ";"}]], "Input",
 CellLabel->"In[4]:="]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["The Potential Energy of a Lennard-Jones Cluster", "Section"],

Cell[TextData[{
 "The total potential energy for ",
 StyleBox["N",
  FontSlant->"Italic"],
 " particles will involve the sum between all pairs of particles:\n\t",
 Cell[BoxData[
  FormBox[
   SubscriptBox["U", "tot"], TraditionalForm]]],
 " = ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    UnderoverscriptBox["\[Sum]", 
     RowBox[{"i", "=", "2"}], "N"], 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"j", "=", "1"}], 
      RowBox[{"i", "-", "1"}]], 
     RowBox[{"\[Phi]", "(", 
      RowBox[{"|", 
       RowBox[{
        SubscriptBox[
         StyleBox["x",
          FontWeight->"Bold",
          FontSlant->"Plain"], "i"], " ", "-", " ", 
        SubscriptBox[
         StyleBox["x",
          FontWeight->"Bold",
          FontSlant->"Plain"], "j"]}], "|"}], ")"}]}]}], TraditionalForm]]],
 "\nwhere \[Phi](",
 Cell[BoxData[
  FormBox["r", TraditionalForm]]],
 ") is the Lennard-Jones potential. For ",
 StyleBox["N",
  FontSlant->"Italic"],
 " particles, this sum has ",
 StyleBox["N",
  FontSlant->"Italic"],
 " (",
 StyleBox["N",
  FontSlant->"Italic"],
 " \[Dash] 1)/2 pairings. (For this reason, we will not be able to simulate \
many particles!) "
}], "Text"],

Cell[TextData[{
 "When only one particle moves, the change in potential energy, \
\[CapitalDelta]",
 StyleBox["U",
  FontSlant->"Italic"],
 ",",
 " requires only those pairings involving the particle that moved. The change \
in energy is\n\t\[CapitalDelta]",
 StyleBox["U",
  FontSlant->"Italic"],
 " = ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    UnderscriptBox["\[Sum]", 
     RowBox[{"i", "\[NotEqual]", "k"}]], 
    RowBox[{"[", 
     RowBox[{
      RowBox[{"\[Phi]", "(", 
       RowBox[{
        SubscriptBox[
         StyleBox["x",
          FontWeight->"Bold",
          FontSlant->"Plain"], "i"], " ", "-", " ", 
        RowBox[{"(", 
         RowBox[{
          SubscriptBox[
           StyleBox["x",
            FontWeight->"Bold",
            FontSlant->"Plain"], "k"], " ", "+", " ", 
          StyleBox[
           RowBox[{"\[CapitalDelta]", 
            StyleBox["x",
             FontWeight->"Bold"]}]]}], ")"}]}], ")"}], " ", "-", " ", 
      RowBox[{"\[Phi]", "(", 
       RowBox[{
        SubscriptBox[
         StyleBox["x",
          FontWeight->"Bold",
          FontSlant->"Plain"], "i"], " ", "-", " ", 
        SubscriptBox[
         StyleBox["x",
          FontWeight->"Bold",
          FontSlant->"Plain"], "k"]}], ")"}]}], "]"}]}], TraditionalForm]]]
}], "Text"],

Cell[TextData[{
 "If \[CapitalDelta]",
 StyleBox["U",
  FontSlant->"Italic"],
 " > 0, the  step is accepted only if a random number between 0 and 1 is less \
than ",
 StyleBox["Exp[-\[CapitalDelta]U / T]",
  FontWeight->"Bold"],
 ". If \[CapitalDelta]",
 StyleBox["U",
  FontSlant->"Italic"],
 " < 0, the step is always accepted."
}], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"Needs", "[", "\"\<BarCharts`\>\"", "]"}], ";", 
    RowBox[{"Needs", "[", "\"\<Histograms`\>\"", "]"}]}], ")"}], 
  ";"}]], "Input",
 CellLabel->"In[5]:="],

Cell[CellGroupData[{

Cell["Initializing the Cluster", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"maxR", " ", "=", " ", "3.0"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
    "all", " ", "particles", " ", "must", " ", "remain", " ", "within", " ", 
     "maxR"}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[6]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"nWalkers", "=", "6"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"walkers", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"RandomReal", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}], ",", "3"}], "]"}], ",", 
     RowBox[{"{", "nWalkers", "}"}]}], "]"}]}], ";"}]}], "Input",
 CellLabel->"In[7]:="],

Cell[BoxData[
 RowBox[{"walkers", " ", "//", " ", "Dimensions"}]], "Input",
 CellLabel->"In[9]:="],

Cell[BoxData[
 RowBox[{"ListPointPlot3D", "[", 
  RowBox[{"walkers", ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"PointSize", "[", "0.02`", "]"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}]}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[10]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Evaluating the total potential", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Phi]lj", "[", "rvec_", "]"}], " ", ":=", " ", 
  RowBox[{"lj2", "[", 
   RowBox[{"rvec", ".", "rvec"}], "]"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{"\[Phi]lj", "[", "r", "]"}], " ", "takes", " ", "a", " ", 
     "vector", " ", "as", " ", "argument"}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[11]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"\[Phi]lj", "[", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "0.5`"}], "}"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "2.`"}], ",", "2.`"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "1"}], "}"}]}]}], "]"}]], "Input",
 CellLabel->"In[12]:="],

Cell[BoxData[
 RowBox[{"\[Phi]lj", "[", 
  RowBox[{
   RowBox[{"walkers", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
   " ", "-", " ", 
   RowBox[{
   "walkers", "\[LeftDoubleBracket]", "3", "\[RightDoubleBracket]"}]}], 
  "]"}]], "Input",
 CellLabel->"In[13]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"potE", "[", "w_List", "]"}], " ", ":=", " ", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"\[Phi]lj", "[", 
       RowBox[{
        RowBox[{"w", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}], 
        " ", "-", " ", 
        RowBox[{"w", "\[LeftDoubleBracket]", "j", "\[RightDoubleBracket]"}]}],
        "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", 
        RowBox[{"i", "-", "1"}]}], "}"}]}], "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "2", ",", " ", "nWalkers"}], "}"}]}], 
   "]"}]}]], "Input",
 CellLabel->"In[14]:="],

Cell[BoxData[
 RowBox[{"potE", "[", "walkers", "]"}]], "Input",
 CellLabel->"In[15]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Evaluating the Cluster Size", "Subsection"],

Cell["\<\
The cluster size can be \"measured\" by the average radius from the \
center-of-mass (COM).\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"avgR", "[", "w_List", "]"}], " ", ":=", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"com", ",", " ", "\[Delta]r"}], "}"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"com", " ", "=", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Plus", " ", "@@", " ", "w"}], ")"}], "/", 
       RowBox[{"Length", "[", "w", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"\[Delta]r", " ", "=", " ", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"#", " ", "-", " ", "com"}], ")"}], "&"}], " ", "/@", " ", 
       "w"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Plus", " ", "@@", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"Sqrt", "[", 
           RowBox[{"#", " ", ".", " ", "#"}], "]"}], "&"}], " ", "/@", " ", 
         "\[Delta]r"}], ")"}]}], "/", 
      RowBox[{"Length", "[", "w", "]"}]}]}]}], "]"}]}]], "Input",
 CellLabel->"In[16]:="],

Cell[BoxData[
 RowBox[{"avgR", "[", "walkers", "]"}]], "Input",
 CellLabel->"In[17]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Evaluating the change in potential", "Subsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.1"}], ";"}]], "Input",
 CellLabel->"In[18]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"s", "[", "]"}], ":=", 
  RowBox[{"RandomReal", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "\[CapitalDelta]x"}], ",", "\[CapitalDelta]x"}], "}"}], 
    ",", "3"}], "]"}]}]], "Input",
 CellLabel->"In[19]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"dpotE", "[", 
   RowBox[{"w_List", ",", "k_Integer", ",", "step_"}], "]"}], " ", ":=", " ", 
  
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"rkOld", ",", "rkNew"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"rkOld", " ", "=", " ", 
      RowBox[{"w", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}]}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"rkNew", " ", "=", " ", 
      RowBox[{"rkOld", " ", "+", " ", "step"}]}], ";", "\[IndentingNewLine]", 
     
     RowBox[{"Sum", "[", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"i", " ", "\[Equal]", " ", "k"}], ",", " ", "0", ",", 
         RowBox[{
          RowBox[{"\[Phi]lj", "[", 
           RowBox[{
            RowBox[{
            "w", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}], " ", 
            "-", " ", "rkNew"}], "]"}], " ", "-", " ", 
          RowBox[{"\[Phi]lj", "[", 
           RowBox[{
            RowBox[{
            "w", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}], " ", 
            "-", " ", "rkOld"}], "]"}]}]}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"i", ",", " ", "1", ",", " ", "nWalkers"}], "}"}]}], 
      "]"}]}]}], "]"}]}]], "Input",
 CellLabel->"In[20]:="],

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"dpotE", "[", 
    RowBox[{"walkers", ",", 
     RowBox[{"RandomInteger", "[", 
      RowBox[{"{", 
       RowBox[{"1", ",", "nWalkers"}], "}"}], "]"}], ",", 
     RowBox[{"s", "[", "]"}]}], "]"}], ",", 
   RowBox[{"{", "5", "}"}]}], "]"}]], "Input",
 CellLabel->"In[21]:="]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["A Random Walker Biased by Lennard-Jones Potential", "Section"],

Cell["\<\
We now have everything needed to define a random walker that can thermally \
agitate particles interacting through the Lennard-Jones potential.\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"step", "[", "w_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"k", ",", "trialStep", ",", "\[CapitalDelta]U"}], "}"}], ",", 
    RowBox[{
     RowBox[{"k", "=", 
      RowBox[{"RandomInteger", "[", 
       RowBox[{"{", 
        RowBox[{"1", ",", "nWalkers"}], "}"}], "]"}]}], ";", 
     RowBox[{"trialStep", "=", 
      RowBox[{"s", "[", "]"}]}], ";", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"#1", ".", "#1"}], "&"}], ")"}], "[", 
         RowBox[{
          RowBox[{"w", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}],
           "+", "trialStep"}], "]"}], ">", 
        SuperscriptBox["maxR", "2"]}], ",", 
       RowBox[{
        RowBox[{"reject", "+=", "1"}], ";", 
        RowBox[{"Return", "[", "w", "]"}]}]}], "]"}], ";", 
     RowBox[{"\[CapitalDelta]U", "=", 
      RowBox[{"dpotE", "[", 
       RowBox[{"w", ",", "k", ",", "trialStep"}], "]"}]}], ";", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"\[CapitalDelta]U", "<", "0"}], ",", 
       RowBox[{"ReplacePart", "[", 
        RowBox[{"w", ",", 
         RowBox[{
          RowBox[{"w", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}],
           "+", "trialStep"}], ",", "k"}], "]"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"RandomReal", "[", "]"}], "<", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"-", 
            FractionBox["\[CapitalDelta]U", "T"]}]]}], ",", 
         RowBox[{"ReplacePart", "[", 
          RowBox[{"w", ",", 
           RowBox[{
            RowBox[{
            "w", "\[LeftDoubleBracket]", "k", "\[RightDoubleBracket]"}], "+", 
            "trialStep"}], ",", "k"}], "]"}], ",", 
         RowBox[{
          RowBox[{"reject", "+=", "1"}], ";", "w"}]}], "]"}]}], "]"}]}]}], 
   "]"}]}]], "Input",
 CellLabel->"In[22]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"T", " ", "=", " ", "0.1"}], ";"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.05"}], ";", " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "adjust", " ", "to", " ", "keep", " ", "rejection", " ", "rate"}], " ", 
     "<", " ", 
     RowBox[{"50", "%"}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"w1", " ", "=", " ", 
   RowBox[{"NestList", "[", 
    RowBox[{"step", ",", " ", "walkers", ",", " ", "10"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"w1", " ", "//", " ", "Dimensions"}]}], "Input",
 CellLabel->"In[23]:="],

Cell[BoxData["reject"], "Input",
 CellLabel->"In[27]:="],

Cell["Run the simulation\[Ellipsis]", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"w1", " ", "=", " ", 
     RowBox[{"NestList", "[", 
      RowBox[{"step", ",", "walkers", ",", "12000"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<Rejection Rate = \>\"", ",", " ", 
      RowBox[{
       RowBox[{"reject", "/", "12000"}], " ", "100.0", " ", "PerCent"}]}], 
     "]"}], ";"}], ")"}], " ", "//", " ", "Timing"}]}], "Input",
 CellLabel->"In[28]:="],

Cell["\<\
Evaluate the total potential energy\[Ellipsis]\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{"potN", " ", "=", " ", 
    RowBox[{"potE", " ", "/@", " ", 
     RowBox[{"w1", "\[LeftDoubleBracket]", 
      RowBox[{"Range", "[", 
       RowBox[{"1", ",", "12000", ",", "nWalkers"}], "]"}], 
      "\[RightDoubleBracket]"}]}]}], ")"}], " ", "//", " ", 
  "Shallow"}]], "Input",
 CellLabel->"In[30]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"potN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Total Potential Energy\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<n/N\>\"", ",", "\"\<\[Sum]\[Phi]\>\""}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[31]:="],

Cell["\<\
Observe how the cluster size changes\[Ellipsis]\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"radN", " ", "=", " ", 
   RowBox[{"avgR", " ", "/@", " ", 
    RowBox[{"w1", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "12000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[32]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"radN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Avg Cluster Radius\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<n/N\>\"", ",", 
      "\"\<\[LeftAngleBracket]\[Delta]r\[RightAngleBracket]\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellLabel->"In[33]:="],

Cell[BoxData[
 RowBox[{"cl1", "=", 
  RowBox[{"ListPointPlot3D", "[", 
   RowBox[{
    RowBox[{"w1", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], ",",
     
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], "}"}]}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}]}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellLabel->"In[34]:="],

Cell[BoxData[
 RowBox[{"cl2", "=", 
  RowBox[{"ListPointPlot3D", "[", 
   RowBox[{
    RowBox[{"w1", "\[LeftDoubleBracket]", "10000", "\[RightDoubleBracket]"}], 
    ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"PointSize", "[", "0.02`", "]"}], "}"}]}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellLabel->"In[35]:="],

Cell["\<\
A annimation of the final state\[Ellipsis]\
\>", "Text"],

Cell[BoxData[
 RowBox[{"ListAnimate", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"ListPointPlot3D", "[", 
     RowBox[{
      RowBox[{"w1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
      ",", 
      RowBox[{"PlotStyle", "\[Rule]", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", 
            FractionBox["maxR", "2"]}], ",", 
           FractionBox["maxR", "2"]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", 
            FractionBox["maxR", "2"]}], ",", 
           FractionBox["maxR", "2"]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", 
            FractionBox["maxR", "2"]}], ",", 
           FractionBox["maxR", "2"]}], "}"}]}], "}"}]}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"12000", "-", 
       RowBox[{"50", " ", "nWalkers"}]}], ",", "12000", ",", "nWalkers"}], 
     "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.4458663789041653`*^9, 3.445866385219451*^9}},
 CellLabel->"In[74]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Ice", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T", " ", "=", " ", "0.05"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.02"}], ";"}]}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "adjust", " ", "to", " ", "keep", " ", "rejection", " ", "rate"}], " ", 
     "<", " ", 
     RowBox[{"50", "%"}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[40]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"w2", " ", "=", " ", 
     RowBox[{"NestList", "[", 
      RowBox[{"step", ",", 
       RowBox[{
       "w1", "\[LeftDoubleBracket]", "12000", "\[RightDoubleBracket]"}], ",", 
       "12000"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<Rejection Rate = \>\"", ",", " ", 
      RowBox[{
       RowBox[{"reject", "/", "12000"}], " ", "100.0", " ", "PerCent"}]}], 
     "]"}], ";"}], ")"}], " ", "//", " ", "Timing"}]}], "Input",
 CellLabel->"In[42]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"potN", " ", "=", " ", 
   RowBox[{"potE", " ", "/@", " ", 
    RowBox[{"w2", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "12000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[44]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"potN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Total Potential Energy\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<n/N\>\"", ",", "\"\<\[Sum]\[Phi]\>\""}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[45]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"radN", " ", "=", " ", 
   RowBox[{"avgR", " ", "/@", " ", 
    RowBox[{"w2", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "12000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[46]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"radN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Avg Cluster Radius\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<n/N\>\"", ",", 
      "\"\<\[LeftAngleBracket]\[Delta]r\[RightAngleBracket]\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellLabel->"In[47]:="],

Cell[BoxData[
 RowBox[{"cl3", "=", 
  RowBox[{"ListPointPlot3D", "[", 
   RowBox[{
    RowBox[{"w2", "\[LeftDoubleBracket]", "12000", "\[RightDoubleBracket]"}], 
    ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"PointSize", "[", "0.02`", "]"}], "}"}]}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", 
          FractionBox["maxR", "2"]}], ",", 
         FractionBox["maxR", "2"]}], "}"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellLabel->"In[48]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Warming", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T", " ", "=", " ", "0.2"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.05"}], ";"}]}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "adjust", " ", "to", " ", "keep", " ", "rejection", " ", "rate"}], " ", 
     "<", " ", 
     RowBox[{"50", "%"}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[49]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"w3", " ", "=", " ", 
     RowBox[{"NestList", "[", 
      RowBox[{"step", ",", 
       RowBox[{
       "w1", "\[LeftDoubleBracket]", "12000", "\[RightDoubleBracket]"}], ",", 
       "6000"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<Rejection Rate = \>\"", ",", " ", 
      RowBox[{
       RowBox[{"reject", "/", "6000"}], " ", "100.0", " ", "PerCent"}]}], 
     "]"}], ";"}], ")"}], " ", "//", " ", "Timing"}]}], "Input",
 CellLabel->"In[51]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"potN", " ", "=", " ", 
   RowBox[{"potE", " ", "/@", " ", 
    RowBox[{"w3", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[53]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"potN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Total Potential Energy\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<n/N\>\"", ",", "\"\<\[Sum]\[Phi]\>\""}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[54]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"radN", " ", "=", " ", 
   RowBox[{"avgR", " ", "/@", " ", 
    RowBox[{"w3", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[55]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"radN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Avg Cluster Radius\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<n/N\>\"", ",", 
      "\"\<\[LeftAngleBracket]\[Delta]r\[RightAngleBracket]\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellLabel->"In[56]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Melting", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T", " ", "=", " ", "0.5"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.1"}], ";"}]}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "adjust", " ", "to", " ", "keep", " ", "rejection", " ", "rate"}], " ", 
     "<", " ", 
     RowBox[{"50", "%"}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[57]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"w4", " ", "=", " ", 
     RowBox[{"NestList", "[", 
      RowBox[{"step", ",", 
       RowBox[{
       "w1", "\[LeftDoubleBracket]", "12000", "\[RightDoubleBracket]"}], ",", 
       "6000"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<Rejection Rate = \>\"", ",", " ", 
      RowBox[{
       RowBox[{"reject", "/", "6000"}], " ", "100.0", " ", "PerCent"}]}], 
     "]"}], ";"}], ")"}], " ", "//", " ", "Timing"}]}], "Input",
 CellLabel->"In[59]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"potN", " ", "=", " ", 
   RowBox[{"potE", " ", "/@", " ", 
    RowBox[{"w4", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[61]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"potN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Total Potential Energy\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<n/N\>\"", ",", "\"\<\[Sum]\[Phi]\>\""}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[62]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"radN", " ", "=", " ", 
   RowBox[{"avgR", " ", "/@", " ", 
    RowBox[{"w4", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[63]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"radN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Avg Cluster Radius\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<n/N\>\"", ",", 
      "\"\<\[LeftAngleBracket]\[Delta]r\[RightAngleBracket]\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellLabel->"In[64]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Liquid", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T", " ", "=", " ", "1.0"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"\[CapitalDelta]x", " ", "=", " ", "0.2"}], ";"}]}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "adjust", " ", "to", " ", "keep", " ", "rejection", " ", "rate"}], " ", 
     "<", " ", 
     RowBox[{"50", "%"}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[65]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"reject", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"w5", " ", "=", " ", 
     RowBox[{"NestList", "[", 
      RowBox[{"step", ",", 
       RowBox[{
       "w1", "\[LeftDoubleBracket]", "12000", "\[RightDoubleBracket]"}], ",", 
       "6000"}], "]"}]}], ";", 
    RowBox[{"Print", "[", 
     RowBox[{"\"\<Rejection Rate = \>\"", ",", " ", 
      RowBox[{
       RowBox[{"reject", "/", "6000"}], " ", "100.0", " ", "PerCent"}]}], 
     "]"}], ";"}], ")"}], " ", "//", " ", "Timing"}]}], "Input",
 CellLabel->"In[67]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"potN", " ", "=", " ", 
   RowBox[{"potE", " ", "/@", " ", 
    RowBox[{"w5", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[69]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"potN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Total Potential Energy\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<n/N\>\"", ",", "\"\<\[Sum]\[Phi]\>\""}], "}"}]}]}], 
  "]"}]], "Input",
 CellLabel->"In[70]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"radN", " ", "=", " ", 
   RowBox[{"avgR", " ", "/@", " ", 
    RowBox[{"w5", "\[LeftDoubleBracket]", 
     RowBox[{"Range", "[", 
      RowBox[{"1", ",", "6000", ",", "nWalkers"}], "]"}], 
     "\[RightDoubleBracket]"}]}]}], ";"}]], "Input",
 CellLabel->"In[71]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"radN", ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Avg Cluster Radius\>\""}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<n/N\>\"", ",", 
      "\"\<\[LeftAngleBracket]\[Delta]r\[RightAngleBracket]\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellLabel->"In[72]:="],

Cell["\<\
A annimation of the final state\[Ellipsis]\
\>", "Text"],

Cell[BoxData[
 RowBox[{"ListAnimate", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"ListPointPlot3D", "[", 
     RowBox[{
      RowBox[{"w5", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
      ",", 
      RowBox[{"PlotStyle", "\[Rule]", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], ",", 
      RowBox[{"PlotRange", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "maxR"}], ",", "maxR"}], "}"}]}], "}"}]}]}], "]"}], 
    ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"6000", "-", 
       RowBox[{"50", " ", "nWalkers"}]}], ",", "6000", ",", "nWalkers"}], 
     "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.4458664336724453`*^9, 3.44586643992553*^9}},
 CellLabel->"In[75]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Summary", "Section"],

Cell[TextData[{
 "We've used the Rosenbluth-Teller-Metropolis algorithm to describe the \
equilibrium potential (or average binding energy) and equilibrium dimensions \
of a small cluster of atoms interacting by way of the Lennard-Jones \
potential. When the temperature of the system is low (",
 StyleBox["T \[LessEqual]",
  FontSlant->"Italic"],
 " 0.2), then the cluster is tightly bound. As the temperature increases, the \
atoms within the cluster separate. Finally, for high temperatures (",
 StyleBox["T ",
  FontSlant->"Italic"],
 "\[GreaterEqual] 1), the atoms are no longer \"clustered\" but appear more \
like un-bound randomly-moving particles or pairs of particles. "
}], "Text"]
}, Closed]]
}, Open  ]]
},
WindowToolbars->{},
WindowSize->{746, 775},
WindowMargins->{{87, Automatic}, {Automatic, 0}},
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (May 21, 2008)",
StyleDefinitions->"TutorialBook.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 32, 0, 51, "Title"],
Cell[625, 25, 60, 3, 41, "Subsubtitle"],
Cell[CellGroupData[{
Cell[710, 32, 31, 0, 86, "Section"],
Cell[744, 34, 1227, 23, 149, "Text"],
Cell[1974, 59, 250, 8, 28, "Input"],
Cell[2227, 69, 466, 14, 28, "Input"],
Cell[2696, 85, 524, 15, 45, "Input"],
Cell[3223, 102, 247, 7, 41, "Text"],
Cell[CellGroupData[{
Cell[3495, 113, 50, 0, 34, "Subsection"],
Cell[3548, 115, 673, 14, 96, "Text"],
Cell[4224, 131, 899, 19, 269, "Text"],
Cell[5126, 152, 659, 15, 59, "Text"],
Cell[5788, 169, 133, 4, 24, "Input"]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell[5970, 179, 66, 0, 54, "Section"],
Cell[6039, 181, 1188, 43, 107, "Text"],
Cell[7230, 226, 1290, 44, 75, "Text"],
Cell[8523, 272, 342, 12, 42, "Text"],
Cell[8868, 286, 220, 7, 26, "Input"],
Cell[CellGroupData[{
Cell[9113, 297, 46, 0, 34, "Subsection"],
Cell[9162, 299, 314, 10, 26, "Input"],
Cell[9479, 311, 459, 15, 64, "Input"],
Cell[9941, 328, 98, 2, 26, "Input"],
Cell[10042, 332, 594, 19, 44, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[10673, 356, 52, 0, 34, "Subsection"],
Cell[10728, 358, 402, 11, 46, "Input"],
Cell[11133, 371, 409, 13, 26, "Input"],
Cell[11545, 386, 277, 8, 26, "Input"],
Cell[11825, 396, 653, 19, 46, "Input"],
Cell[12481, 417, 87, 2, 24, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[12605, 424, 49, 0, 34, "Subsection"],
Cell[12657, 426, 115, 3, 23, "Text"],
Cell[12775, 431, 1006, 29, 81, "Input"],
Cell[13784, 462, 87, 2, 24, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[13908, 469, 56, 0, 34, "Subsection"],
Cell[13967, 471, 118, 3, 24, "Input"],
Cell[14088, 476, 266, 9, 24, "Input"],
Cell[14357, 487, 1302, 35, 98, "Input"],
Cell[15662, 524, 334, 10, 42, "Input"]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell[16045, 540, 68, 0, 54, "Section"],
Cell[16116, 542, 167, 3, 41, "Text"],
Cell[16286, 547, 1970, 55, 207, "Input"],
Cell[18259, 604, 773, 21, 97, "Input"],
Cell[19035, 627, 56, 1, 26, "Input"],
Cell[19094, 630, 45, 0, 23, "Text"],
Cell[19142, 632, 529, 14, 81, "Input"],
Cell[19674, 648, 70, 2, 24, "Text"],
Cell[19747, 652, 355, 10, 26, "Input"],
Cell[20105, 664, 309, 8, 46, "Input"],
Cell[20417, 674, 71, 2, 24, "Text"],
Cell[20491, 678, 301, 8, 26, "Input"],
Cell[20795, 688, 354, 10, 45, "Input"],
Cell[21152, 700, 812, 25, 62, "Input"],
Cell[21967, 727, 876, 27, 82, "Input"],
Cell[22846, 756, 66, 2, 23, "Text"],
Cell[22915, 760, 1185, 34, 151, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[24137, 799, 22, 0, 54, "Section"],
Cell[24162, 801, 459, 15, 43, "Input"],
Cell[24624, 818, 616, 17, 81, "Input"],
Cell[25243, 837, 301, 8, 26, "Input"],
Cell[25547, 847, 309, 8, 46, "Input"],
Cell[25859, 857, 301, 8, 26, "Input"],
Cell[26163, 867, 354, 10, 45, "Input"],
Cell[26520, 879, 876, 27, 82, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[27433, 911, 26, 0, 54, "Section"],
Cell[27462, 913, 458, 15, 43, "Input"],
Cell[27923, 930, 614, 17, 81, "Input"],
Cell[28540, 949, 300, 8, 26, "Input"],
Cell[28843, 959, 309, 8, 46, "Input"],
Cell[29155, 969, 300, 8, 26, "Input"],
Cell[29458, 979, 354, 10, 45, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[29849, 994, 26, 0, 54, "Section"],
Cell[29878, 996, 457, 15, 43, "Input"],
Cell[30338, 1013, 614, 17, 81, "Input"],
Cell[30955, 1032, 300, 8, 26, "Input"],
Cell[31258, 1042, 309, 8, 46, "Input"],
Cell[31570, 1052, 300, 8, 26, "Input"],
Cell[31873, 1062, 354, 10, 45, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[32264, 1077, 25, 0, 54, "Section"],
Cell[32292, 1079, 457, 15, 43, "Input"],
Cell[32752, 1096, 614, 17, 81, "Input"],
Cell[33369, 1115, 300, 8, 26, "Input"],
Cell[33672, 1125, 309, 8, 46, "Input"],
Cell[33984, 1135, 300, 8, 26, "Input"],
Cell[34287, 1145, 354, 10, 45, "Input"],
Cell[34644, 1157, 66, 2, 23, "Text"],
Cell[34713, 1161, 1004, 29, 79, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[35754, 1195, 26, 0, 54, "Section"],
Cell[35783, 1197, 692, 13, 96, "Text"]
}, Closed]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)