energiesdomains.nb

(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 10.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     35175,        820]
NotebookOptionsPosition[     34372,        789]
NotebookOutlinePosition[     34708,        804]
CellTagsIndexPosition[     34665,        801]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["Need to worry about \[OpenCurlyDoubleQuote]metastable\
\[CloseCurlyDoubleQuote] states...", "Text",
 CellChangeTimes->{{3.659711137467695*^9, 3.6597111451376343`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"diameter", "=", 
   RowBox[{"{", 
    RowBox[{
    "5.0", ",", "5.5", ",", "6.0", ",", "6.5", ",", "7.0", ",", "7.5", ",", 
     "8.0", ",", "8.5", ",", "9.0", ",", "9.5", ",", "10.0", ",", "10.5", ",",
      "11.0", ",", "11.5", ",", "12.0", ",", "12.5", ",", "13.0", ",", "13.5",
      ",", "14.0", ",", "14.5", ",", "15.0", ",", "15.5", ",", "16.0", ",", 
     "16.5", ",", "17.0", ",", "17.5", ",", "18.0", ",", "18.5", ",", "19.0", 
     ",", "19.5", ",", "20", ",", "21.0", ",", "22.0", ",", "23.0", ",", 
     "24.0", ",", "30"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.659817771736498*^9, 3.6598178169345503`*^9}, {
  3.659819009281932*^9, 3.659819010350019*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"volume", " ", "=", " ", 
   RowBox[{
    FractionBox[
     RowBox[{"4", " ", "\[Pi]", " "}], "3"], 
    SuperscriptBox[
     RowBox[{"(", 
      FractionBox["diameter", "2"], ")"}], "3"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.6598179127013187`*^9, 3.6598179259152308`*^9}, {
  3.659818299633399*^9, 3.6598183019554367`*^9}, {3.659819001725017*^9, 
  3.659819015016045*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"domainEnergy", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
    "0.13331", ",", "0.47125", ",", "0.278", ",", "0.56", ",", "0.8878", ",", 
     "0.67669", ",", "1.847", ",", "1.8687", ",", "1.775", ",", "1.51", ",", 
     "1.99", ",", "1.93", ",", "4.1", ",", "4.71", ",", "4.33", ",", "6.82", 
     ",", "7.75", ",", "6.59", ",", "9.25", ",", "11.12", ",", "10.94", ",", 
     "10.1", ",", "15.74", ",", "17.7", ",", "10.32"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGCQAWIQrafzdG6oy2tH478xC0B0au+/syD6mX/jORD96MyK
OyD6pUbIQxBdEJ0yPQxI91gfBdN7WAXngeh/uU1gmsVly2IQ7eMWtgxEhxzv
WQWiRfZMANMNn3ZvBdFvnP5vA9HWIc3HQLSbjPRxEJ0V33YJRPPPErkMohWk
tt4B6/8adBdML2t/CqJbY2+D6fMZj/XDgbTy6R9gWudtvE8EkD4X3AOmn/Nq
hoPoGbVuYDpibdMEED3x2A4wfWPyrhkg+snTK2B6km+NSiSQbsmcDqb3h2Rp
gei1vrPAtLR2vBmIvmjYCKZPaRTagOjX22eC6c1cr31BdOwhXj8QbaJcI9YD
Ui+8Ckw/4zimAaKlfr4A0zdOZBuDaEnLFjB9JKnBAUSbbZ4PppkM9vmB6IQi
5QAQvdlfIhJEPy2wBtMABM3VRA==
  "]],

Cell[BoxData[
 RowBox[{
  RowBox[{"electrostaticEnergy", "=", " ", 
   RowBox[{"{", 
    RowBox[{
    "0.115", ",", "0.1839", ",", "0.2216", ",", "0.3455", ",", "0.33525", ",",
      "0.69943", ",", "0.7445", ",", "0.7264", ",", "0.835", ",", "1.19", ",",
      "1.248", ",", "1.06", ",", "2.66", ",", "2.43", ",", "2.89", ",", 
     "2.51", ",", "3.35", ",", "3.09", ",", "3.44", ",", "4.02", ",", "5.71", 
     ",", "5.08", ",", "6.06", ",", "7.52", ",", "5.39"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGCQAmIQvV3i1+1Ql9eO034euQuitxRdfwiijQq8H4Fogylu
M8KANNOuTDBdcfHdPBC9ZmHsfBB9bcLsJSD6gNCZpSBaubx7NYjuWt0Lppt/
rjsOor9etzoBov/ITLoMoqX2918B0f7lV++C6IOybvdAdHzX22cgmlV19XMQ
vSNivUE4kFbIPQGmqz7o+UaA9F9m8wfRSWeeRoDokx6MkSDagaluIojmy10B
pkP9mmeCaP+9s8G0zyVl1Uggfem6E5i20jPXBtHTjILB9LV7XuYg2uZbLpje
oZxhC6K97h0G02v2ePmB6C8rOPxB9BKxORI9QPrZ7z1gWsnssyaI/uUorQ2i
F5zYZQKiowr/gulP1fcdQDQHS6ITiJaKSQgA0Xfcm8F0+v/ySBB9QbYmCkQD
AFrzyJw=
  "]],

Cell[BoxData[
 RowBox[{
  RowBox[{"totalEnergy", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "5.585"}], ",", 
     RowBox[{"-", "6.175"}], ",", 
     RowBox[{"-", "7.4179"}], ",", 
     RowBox[{"-", "9.41"}], ",", 
     RowBox[{"-", "11.46"}], ",", 
     RowBox[{"-", "13.5"}], ",", 
     RowBox[{"-", "22.29"}], ",", 
     RowBox[{"-", "22.56"}], ",", 
     RowBox[{"-", "22.202"}], ",", 
     RowBox[{"-", "28.87"}], ",", 
     RowBox[{"-", "36.36"}], ",", 
     RowBox[{"-", "37.25"}], ",", 
     RowBox[{"-", "37.38"}], ",", 
     RowBox[{"-", "46.83"}], ",", 
     RowBox[{"-", "51.3"}], ",", 
     RowBox[{"-", "57.18"}], ",", 
     RowBox[{"-", "71.29"}], ",", 
     RowBox[{"-", "83.18"}], ",", 
     RowBox[{"-", "88.223"}], ",", 
     RowBox[{"-", "95.4"}], ",", 
     RowBox[{"-", "103.56"}], ",", 
     RowBox[{"-", "116.24"}], ",", 
     RowBox[{"-", "126.49"}], ",", 
     RowBox[{"-", "142.6"}], ",", 
     RowBox[{"-", "159.95"}]}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.659705974597969*^9, 3.659705975024694*^9}, {
   3.659706020905076*^9, 3.659706028129653*^9}, {3.659706137552537*^9, 
   3.659706140960836*^9}, {3.659707610952242*^9, 3.6597076188176117`*^9}, {
   3.6597076505127907`*^9, 3.659707655136298*^9}, {3.659707709942719*^9, 
   3.659707744839323*^9}, {3.6597079719903584`*^9, 3.659707982439522*^9}, {
   3.6597080820007963`*^9, 3.6597080866579647`*^9}, {3.659708156379036*^9, 
   3.659708160199201*^9}, {3.659708225909191*^9, 3.659708236147286*^9}, {
   3.659708816603302*^9, 3.6597088229446898`*^9}, {3.659711105414178*^9, 
   3.659711108832261*^9}, {3.659711185011084*^9, 3.659711187183123*^9}, {
   3.659711639634636*^9, 3.659711647149457*^9}, {3.65971170705018*^9, 
   3.659711714623616*^9}, {3.659712818263404*^9, 3.659712823334811*^9}, {
   3.659712872420162*^9, 3.6597128775549583`*^9}, {3.659712966793507*^9, 
   3.659712968447443*^9}, {3.65971301864258*^9, 3.6597130225645323`*^9}, {
   3.65971315391183*^9, 3.659713160792383*^9}, 3.659817159895822*^9, {
   3.6598172245718737`*^9, 3.659817246809388*^9}, {3.6598173064414682`*^9, 
   3.6598173319134693`*^9}, {3.6598173858786707`*^9, 
   3.6598174159558773`*^9}, {3.659817493091663*^9, 3.659817508852892*^9}, {
   3.659817616919261*^9, 3.659817618569886*^9}, {3.659817686966467*^9, 
   3.659817710822912*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Length", "[", "domainEnergy", "]"}]], "Input",
 CellChangeTimes->{{3.6597078141377573`*^9, 3.6597078184295263`*^9}}],

Cell[BoxData["25"], "Output",
 CellChangeTimes->{3.659707818707921*^9, 3.659708103043898*^9, 
  3.659708784379023*^9, 3.6597111114059467`*^9, 3.65971119002752*^9, 
  3.6597116503997*^9, 3.659711717895895*^9, 3.659712790530643*^9, 
  3.659712827173398*^9, 3.659712971482122*^9, 3.659713024750319*^9, 
  3.6597131641183357`*^9, 3.659816939309144*^9, 3.659817249631159*^9, 
  3.659817622771188*^9, 3.6598177139538527`*^9, 3.6598178842800407`*^9, 
  3.6598190184051933`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"s", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"volume", "[", 
           RowBox[{"[", "k", "]"}], "]"}], "*", 
          SuperscriptBox["10", 
           RowBox[{"-", "3"}]]}], ",", 
         RowBox[{"domainEnergy", "[", 
          RowBox[{"[", "k", "]"}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", "25"}], "}"}]}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Thick", ",", "Blue"}], "}"}]}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "700"}], ",", 
     RowBox[{"Frame", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], 
     ",", "\[IndentingNewLine]", 
     RowBox[{"FrameTicksStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{"Black", ",", "18"}], "]"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{"{", 
        RowBox[{"Black", ",", 
         RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
         RowBox[{"FontSize", "\[Rule]", "14"}], ",", 
         RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], ",", 
         RowBox[{"FontColor", "\[Rule]", "Black"}]}], "}"}], "]"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Style", "[", 
         RowBox[{
         "\"\<Volume [\!\(\*SuperscriptBox[\(10\), \
\(3\)]\)\!\(\*SuperscriptBox[\(nm\), \(3\)]\)]\>\"", ",", 
          RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], 
          ",", 
          RowBox[{"FontColor", "\[Rule]", "Black"}], ",", 
          RowBox[{"FontSize", "\[Rule]", "22"}]}], "]"}], ",", 
        RowBox[{"Style", "[", 
         RowBox[{"\"\<E [arb units]\>\"", ",", 
          RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], 
          ",", 
          RowBox[{"FontColor", "\[Rule]", "Black"}], ",", 
          RowBox[{"FontSize", "\[Rule]", "22"}]}], "]"}]}], "}"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{"LineLegend", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Directive", "[", 
             RowBox[{"Blue", ",", "DotDashed", ",", 
              RowBox[{"Thickness", "[", "0.0025", "]"}]}], "]"}], ",", 
            RowBox[{"Directive", "[", 
             RowBox[{"Red", ",", 
              RowBox[{"Thickness", "[", "0.0025", "]"}]}], "]"}]}], "}"}], 
          ",", 
          RowBox[{"(*", 
           RowBox[{
            RowBox[{
            "can", " ", "use", " ", "linelegend", " ", "with", " ", "a", " ", 
             "frame"}], ",", " ", 
            RowBox[{"use", " ", "Directive", " ", "on", " ", "the", " ", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"Directive", "[", 
                RowBox[{"Black", ",", 
                 RowBox[{"Dashing", ".."}]}], "]"}], ",", "Black", ",", 
               "\"\<0\>\"", ",", "\"\<10\>\""}], "}"}], " ", "to", " ", "get",
              " ", "hte", " ", "proper", " ", "dashing"}], ",", " ", 
            RowBox[{
            "then", " ", "use", " ", "text", " ", "annotations", " ", "with", 
             " ", "the", " ", "color", " ", "for", " ", "strains"}]}], "*)"}],
           "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Style", "[", 
             RowBox[{"\"\<\!\(\*SubscriptBox[\(E\), \(wall\)]\)\>\"", ",", 
              RowBox[{"FontSize", "\[Rule]", "15"}]}], "]"}], ",", 
            "\[IndentingNewLine]", 
            RowBox[{"Style", "[", 
             RowBox[{
             "\"\<\!\(\*SubscriptBox[\(E\), \(electrostatic\)]\) ; \
\[CapitalDelta]E \[EqualTilde]  7 * \!\(\*SuperscriptBox[\(10\), \
\(-5\)]\)\>\"", ",", 
              RowBox[{"FontSize", "\[Rule]", "15"}]}], "]"}]}], 
           "\[IndentingNewLine]", "}"}]}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0.418", ",", "0.82"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0.8", ",", "0.3"}], "}"}]}], "}"}]}], "]"}]}]}], "]"}]}], 
  ";", 
  RowBox[{"t", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"volume", "[", 
           RowBox[{"[", "k", "]"}], "]"}], "*", 
          SuperscriptBox["10", 
           RowBox[{"-", "3"}]]}], ",", 
         RowBox[{"electrostaticEnergy", "[", 
          RowBox[{"[", "k", "]"}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", "25"}], "}"}]}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Thick", ",", "Red"}], "}"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"u", "=", " ", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"volume", "[", 
           RowBox[{"[", "k", "]"}], "]"}], "*", 
          SuperscriptBox["10", 
           RowBox[{"-", "3"}]]}], ",", 
         RowBox[{"totalEnergy", "[", 
          RowBox[{"[", "k", "]"}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"k", ",", "1", ",", "25"}], "}"}]}], "]"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "700"}], ",", 
     RowBox[{"Frame", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], 
     ",", "\[IndentingNewLine]", 
     RowBox[{"FrameTicksStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{"Black", ",", "18"}], "]"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",", "\[IndentingNewLine]", 
     RowBox[{"BaseStyle", "\[Rule]", 
      RowBox[{"Directive", "[", 
       RowBox[{"{", 
        RowBox[{"Black", ",", 
         RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}], ",", 
         RowBox[{"FontSize", "\[Rule]", "14"}], ",", 
         RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], ",", 
         RowBox[{"FontColor", "\[Rule]", "Black"}]}], "}"}], "]"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Style", "[", 
         RowBox[{
         "\"\<Volume [\!\(\*SuperscriptBox[\(10\), \
\(3\)]\)\!\(\*SuperscriptBox[\(nm\), \(3\)]\)]\>\"", ",", 
          RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], 
          ",", 
          RowBox[{"FontColor", "\[Rule]", "Black"}], ",", 
          RowBox[{"FontSize", "\[Rule]", "22"}]}], "]"}], ",", 
        RowBox[{"Style", "[", 
         RowBox[{"\"\<Total E [arb units]\>\"", ",", 
          RowBox[{"FontFamily", "\[Rule]", "\"\<Latin Modern Roman\>\""}], 
          ",", 
          RowBox[{"FontColor", "\[Rule]", "Black"}], ",", 
          RowBox[{"FontSize", "\[Rule]", "22"}]}], "]"}]}], "}"}]}]}], 
    "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.659707759814535*^9, 3.6597077916092653`*^9}, {
  3.659707919885138*^9, 3.65970792512683*^9}, {3.659708015616788*^9, 
  3.6597080207577333`*^9}, {3.659816943640567*^9, 3.659817059674006*^9}, {
  3.6598174342012873`*^9, 3.659817434770774*^9}, {3.659817733769067*^9, 
  3.659817736717772*^9}, {3.659817930327495*^9, 3.6598180856394053`*^9}, {
  3.65981812116486*^9, 3.659818163168996*^9}, {3.659818195192202*^9, 
  3.659818235213449*^9}, {3.659819501273551*^9, 3.659819561415345*^9}, {
  3.6598207147008257`*^9, 3.659820754385738*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"{", 
   RowBox[{"s", ",", "t"}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.659707793235854*^9, 3.659707802344597*^9}, {
  3.659707933545024*^9, 3.659707944004035*^9}, {3.6597079929963083`*^9, 
  3.65970801454191*^9}, {3.65981946036095*^9, 3.6598194968986588`*^9}}],

Cell[BoxData[
 TagBox[
  GraphicsBox[{{{{}, {{}, 
       {RGBColor[0, 0, 1], PointSize[0.012833333333333334`], Thickness[Large],
         PointBox[CompressedData["
1:eJxTTMoPSmViYGCQBGIQnVqsaxF0cIP98QzZUF+Bg/b7tCKVxXy32WvG9B/6
qnHP/vf9mNxPn/fYz8gTaj5w6qK9W//3vRfTD9m/Cdwh1/r6of1ni+Rbn/4c
s+ctXNN9O+ON/Z5NZ/ZZuZ2x3/W94GTh8qf2OgvZzztrXLSPSrG+79/71/6Q
WPPV6ClX7PUmLPhh+Oyv/Z8Yt54ruTfs08Dgj/0Wsf08srvvQO3/YQ9x3wP7
64sLbLmu/7fn2eH8Iif1sf3DKpF17g//Qd37DKpfwOHlHJWy1uqXUPVCDhD3
v7EPAjk3UNBhY9Idr1MxH+w9wAZIO0D888meAQzkHUrjosNWz/liD7Y+RsoB
4r9vEPkGJYfy9iwBl66f9kDDgTaoOUD8+wfqHlWHYuWnq4uM/tsbg4GKA8T/
jA7VIOkqfQcBm8truS4yOYClNxs6QMKDxWEJyLmLVRwAmde7eg==
         
         "]]}, {}}, {}}, {{}, {{}, 
       {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], Thickness[Large],
         PointBox[CompressedData["
1:eJxTTMoPSmViYGCQBGIQnVqsaxF0cIN9oS3X9cUFe+33aUUqi/lus7/EcI+J
s+u4/e/7MbmfPu+xXz8ltSM59oy9W//3vRfTD9m/VDPkWCNzzf6zRfKtT3+O
2d/+WZe1p+Sq/Z5NZ/ZZuZ2xP/L6F+vuxGf2OgvZzztrXLTPff575cdLz+0P
iTVfjZ5yxX4if5XZarvn9n9i3Hqu5N6w3yHX+jpwxyv7LWL7eWR337EHucaW
67M9xH0P7LXbxW6e+/7ZnmeH84uc1Mf23zRi+g99/QB17zN7j4dVIuvcWR1e
zlEpa61+CfUPswPE/W/s5UHGy7E7bEy643Uq5oM9UDFQB4sDxD+f7M+eAQEu
h9K46LDVc75A3cPhAPHfN3uw9lZuh/L2LAGXrp9Q/QIOEP/+sQc797qYQ7Hy
09VFRv/tg8AGiDhA/M/oAHaOrYSDgM3ltVwXmRwg+uUcIOHB4gD0DNBHog4A
6SDERA==
         "]]}, {}}, {}}}, InsetBox[
     TemplateBox[{
      StyleBox["\"\\!\\(\\*SubscriptBox[\\(E\\), \\(wall\\)]\\)\"", FontSize -> 
        15, StripOnInput -> False],StyleBox[
       "\"\\!\\(\\*SubscriptBox[\\(E\\), \\(electrostatic\\)]\\) ; \
\[CapitalDelta]E \[EqualTilde]  7 * \\!\\(\\*SuperscriptBox[\\(10\\), \
\\(-5\\)]\\)\"", FontSize -> 15, StripOnInput -> False]},
      "LineLegend",
      DisplayFunction->(FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{0, Small, Small, Small}], 
                    Thickness[0.045]], {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{0, Small, Small, Small}], 
                    Thickness[0.045]], 
                    PointBox[
                    NCache[{
                    Scaled[{
                    Rational[1, 2], 
                    Rational[1, 2]}]}, {
                    Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0], 
                    Thickness[0.045]], {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0], 
                    Thickness[0.045]], 
                    PointBox[
                    NCache[{
                    Scaled[{
                    Rational[1, 2], 
                    Rational[1, 2]}]}, {
                    Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ),
      Editable->True,
      InterpretationFunction:>(RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0, 0, 1], 
                    RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> 
                    True, FrameStyle -> RGBColor[0., 0., 0.6666666666666666], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], "RGBColor[0, 0, 1]"], Appearance -> 
                    None, BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0, 0, 1]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0, 0, 1], Editable -> False, Selectable -> 
                    False], ",", 
                   RowBox[{"Dashing", "[", 
                    RowBox[{"{", 
                    RowBox[{"0", ",", "Small", ",", "Small", ",", "Small"}], 
                    "}"}], "]"}], ",", 
                   RowBox[{"Thickness", "[", "0.0025`", "]"}]}], "]"}], ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[1, 0, 0], 
                    RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> 
                    True, FrameStyle -> RGBColor[0.6666666666666666, 0., 0.], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], "RGBColor[1, 0, 0]"], Appearance -> 
                    None, BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[1, 0, 0]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[1, 0, 0], Editable -> False, Selectable -> 
                    False], ",", 
                   RowBox[{"Thickness", "[", "0.0025`", "]"}]}], "]"}]}], 
             "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2}], "}"}], ",", 
           RowBox[{"Joined", "\[Rule]", 
             RowBox[{"{", "False", "}"}]}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", 
             RowBox[{"{", 
               RowBox[{"{", 
                 RowBox[{"False", ",", "Automatic"}], "}"}], "}"}]}]}], 
         "]"}]& )], Scaled[{0.418, 0.82}], ImageScaled[{0.8, 0.3}],
     BaseStyle->{FontSize -> Larger},
     FormatType->StandardForm]},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{0, 0},
   BaseStyle->Directive[{
      GrayLevel[0], FontWeight -> "Bold", FontSize -> 14, FontFamily -> 
      "Latin Modern Roman", FontColor -> GrayLevel[0]}],
   DisplayFunction->Identity,
   Frame->{{True, False}, {True, False}},
   FrameLabel->{{
      FormBox[
       StyleBox[
       "\"E [arb units]\"", FontFamily -> "Latin Modern Roman", FontColor -> 
        GrayLevel[0], FontSize -> 22, StripOnInput -> False], 
       TraditionalForm], None}, {
      FormBox[
       StyleBox[
       "\"Volume [\\!\\(\\*SuperscriptBox[\\(10\\), \
\\(3\\)]\\)\\!\\(\\*SuperscriptBox[\\(nm\\), \\(3\\)]\\)]\"", FontFamily -> 
        "Latin Modern Roman", FontColor -> GrayLevel[0], FontSize -> 22, 
        StripOnInput -> False], TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   FrameTicksStyle->Directive[
     GrayLevel[0], 18],
   GridLines->{None, None},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImageSize->700,
   Method->{},
   PlotRange->{{0, 2.572440784514442}, {0, 17.7}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.02], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}],
  InterpretTemplate[Legended[
    Graphics[{{{}, {{{}, {
          Hue[0.67, 0.6, 0.6], 
          Directive[
           PointSize[0.012833333333333334`], 
           AbsoluteThickness[1.6], FontColor -> GrayLevel[0], 
           Thickness[Large], 
           RGBColor[0, 0, 1]], 
          Point[CompressedData["
1:eJxTTMoPSmViYGCQBGIQnVqsaxF0cIP98QzZUF+Bg/b7tCKVxXy32WvG9B/6
qnHP/vf9mNxPn/fYz8gTaj5w6qK9W//3vRfTD9m/Cdwh1/r6of1ni+Rbn/4c
s+ctXNN9O+ON/Z5NZ/ZZuZ2x3/W94GTh8qf2OgvZzztrXLSPSrG+79/71/6Q
WPPV6ClX7PUmLPhh+Oyv/Z8Yt54ruTfs08Dgj/0Wsf08srvvQO3/YQ9x3wP7
64sLbLmu/7fn2eH8Iif1sf3DKpF17g//Qd37DKpfwOHlHJWy1uqXUPVCDhD3
v7EPAjk3UNBhY9Idr1MxH+w9wAZIO0D888meAQzkHUrjosNWz/liD7Y+RsoB
4r9vEPkGJYfy9iwBl66f9kDDgTaoOUD8+wfqHlWHYuWnq4uM/tsbg4GKA8T/
jA7VIOkqfQcBm8truS4yOYClNxs6QMKDxWEJyLmLVRwAmde7eg==
           
           "]]}, {}}}, {}}, {{}, {{{}, {
          Hue[0.67, 0.6, 0.6], 
          Directive[
           PointSize[0.012833333333333334`], 
           AbsoluteThickness[1.6], 
           Thickness[Large], 
           RGBColor[1, 0, 0]], 
          Point[CompressedData["
1:eJxTTMoPSmViYGCQBGIQnVqsaxF0cIN9oS3X9cUFe+33aUUqi/lus7/EcI+J
s+u4/e/7MbmfPu+xXz8ltSM59oy9W//3vRfTD9m/VDPkWCNzzf6zRfKtT3+O
2d/+WZe1p+Sq/Z5NZ/ZZuZ2xP/L6F+vuxGf2OgvZzztrXLTPff575cdLz+0P
iTVfjZ5yxX4if5XZarvn9n9i3Hqu5N6w3yHX+jpwxyv7LWL7eWR337EHucaW
67M9xH0P7LXbxW6e+/7ZnmeH84uc1Mf23zRi+g99/QB17zN7j4dVIuvcWR1e
zlEpa61+CfUPswPE/W/s5UHGy7E7bEy643Uq5oM9UDFQB4sDxD+f7M+eAQEu
h9K46LDVc75A3cPhAPHfN3uw9lZuh/L2LAGXrp9Q/QIOEP/+sQc797qYQ7Hy
09VFRv/tg8AGiDhA/M/oAHaOrYSDgM3ltVwXmRwg+uUcIOHB4gD0DNBHog4A
6SDERA==
           "]]}, {}}}, {}}}, {
     DisplayFunction -> Identity, PlotRangePadding -> {{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.02], 
         Scaled[0.05]}}, AxesOrigin -> {0, 0}, 
      PlotRange -> {{0, 2.572440784514442}, {0, 17.7}}, DisplayFunction -> 
      Identity, AspectRatio -> GoldenRatio^(-1), Axes -> {True, True}, 
      AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, BaseStyle -> Directive[{
         GrayLevel[0], FontWeight -> "Bold", FontSize -> 14, FontFamily -> 
         "Latin Modern Roman", FontColor -> GrayLevel[0]}], DisplayFunction :> 
      Identity, Frame -> {{True, False}, {True, False}}, FrameLabel -> {{
         Style[
         "E [arb units]", FontFamily -> "Latin Modern Roman", FontColor -> 
          GrayLevel[0], FontSize -> 22], None}, {
         Style[
         "Volume [\!\(\*SuperscriptBox[\(10\), \
\(3\)]\)\!\(\*SuperscriptBox[\(nm\), \(3\)]\)]", FontFamily -> 
          "Latin Modern Roman", FontColor -> GrayLevel[0], FontSize -> 22], 
         None}}, FrameTicks -> {{Automatic, Automatic}, {
        Automatic, Automatic}}, FrameTicksStyle -> Directive[
        GrayLevel[0], 18], GridLines -> {None, None}, GridLinesStyle -> 
      Directive[
        GrayLevel[0.5, 0.4]], ImageSize -> 700, Method -> {}, 
      PlotRange -> {{0, 2.572440784514442}, {0, 17.7}}, PlotRangeClipping -> 
      True, PlotRangePadding -> {{
         Scaled[0.02], 
         Scaled[0.02]}, {
         Scaled[0.02], 
         Scaled[0.05]}}, Ticks -> {Automatic, Automatic}}], 
    Placed[
     Unevaluated[
      LineLegend[{
        Directive[
         RGBColor[0, 0, 1], 
         Dashing[{0, Small, Small, Small}], 
         Thickness[0.0025]], 
        Directive[
         RGBColor[1, 0, 0], 
         Thickness[0.0025]]}, {
        Style["\!\(\*SubscriptBox[\(E\), \(wall\)]\)", FontSize -> 15], 
        Style[
        "\!\(\*SubscriptBox[\(E\), \(electrostatic\)]\) ; \[CapitalDelta]E \
\[EqualTilde]  7 * \!\(\*SuperscriptBox[\(10\), \(-5\)]\)", FontSize -> 15]}, 
       Joined -> {False}, LabelStyle -> {}, LegendLayout -> "Column", 
       LegendMarkers -> {{False, Automatic}}]], {{0.418, 0.82}, {0.8, 0.3}}, 
     Identity]]& ],
  AutoDelete->True,
  Editable->True,
  SelectWithContents->False,
  Selectable->True]], "Output",
 CellChangeTimes->{{3.659707797730761*^9, 3.659707802726017*^9}, {
   3.659707926933926*^9, 3.6597079442995777`*^9}, {3.659707989734167*^9, 
   3.659708022948165*^9}, 3.659708104936406*^9, 3.6597081652405043`*^9, 
   3.6597082474613037`*^9, 3.659711114315867*^9, 3.659711191155529*^9, 
   3.659711652209559*^9, 3.65971171953549*^9, 3.659712829071719*^9, 
   3.659712973070711*^9, 3.659713026557864*^9, 3.659713165514721*^9, {
   3.659816973320565*^9, 3.659817024924039*^9}, 3.659817061929729*^9, 
   3.659817251410819*^9, 3.659817336965542*^9, 3.6598174196342497`*^9, 
   3.659817519907035*^9, 3.65981762483326*^9, {3.6598177156418247`*^9, 
   3.659817739069887*^9}, {3.659817974048779*^9, 3.659817993588952*^9}, 
   3.659818087521462*^9, 3.659818132837584*^9, 3.6598181652407017`*^9, 
   3.659818205420801*^9, 3.6598182376717978`*^9, 3.65981902053049*^9, {
   3.659819484484672*^9, 3.659819546295197*^9}, {3.659820744801161*^9, 
   3.6598207555333357`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["u"], "Input",
 CellChangeTimes->{3.6597081119657097`*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0, 1, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
     1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGCQBGIQnVqsaxF0cIP99cUFtlzhYgf2aUUqi/luszcGgc0S
B37fj8n99HmPfew/519vV8secOv/vvdi+iH7oB1yra8vKh34bJF869OfY/Zv
AkECagf2bDqzz8rtjD0DGGgf0FnIft5Z46I9F8gCL7MDh8Sar0ZPuWLff+ir
Rky/2YE/MW49V3Jv2B88tdB1m7HZgS1i+3lkd9+xl299Hbjjns0BiPse2K9z
f1gloud0gGeH84uc1McQ8xc4Qd37zL7QFmjDBqcDL+eolLVWv4TYl+UOdf8b
+zQQWOZ5YGPSHa9TMR8g/p3uA/XPJ/vDIOdcCjxQGhcdtnrOF4h/TodA/ffN
vvkA0IF8YQfK27MEXLp+2s+aCQQ3w6H+/WO/BGjc9ceRB4qVn64uMvoP8R9/
LNT/jA5g/vz4AwI2l9dyXWRyAIfv5URoeLA4gN33L/kAADw6xo4=
      "]]}, {}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  BaseStyle->Directive[{
     GrayLevel[0], FontWeight -> "Bold", FontSize -> 14, FontFamily -> 
     "Latin Modern Roman", FontColor -> GrayLevel[0]}],
  DisplayFunction->Identity,
  Frame->{{True, False}, {True, False}},
  FrameLabel->{{
     FormBox[
      StyleBox[
      "\"Total E [arb units]\"", FontFamily -> "Latin Modern Roman", 
       FontColor -> GrayLevel[0], FontSize -> 22, StripOnInput -> False], 
      TraditionalForm], None}, {
     FormBox[
      StyleBox[
      "\"Volume [\\!\\(\\*SuperscriptBox[\\(10\\), \
\\(3\\)]\\)\\!\\(\\*SuperscriptBox[\\(nm\\), \\(3\\)]\\)]\"", FontFamily -> 
       "Latin Modern Roman", FontColor -> GrayLevel[0], FontSize -> 22, 
       StripOnInput -> False], TraditionalForm], None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  FrameTicksStyle->Directive[
    GrayLevel[0], 18],
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImageSize->700,
  Method->{},
  PlotRange->{{0, 2.572440784514442}, {-159.95, 0}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.659708112269355*^9, 3.659708167166286*^9, 3.659711115941337*^9, 
   3.659711192824354*^9, 3.659711654382949*^9, 3.6597117211783733`*^9, 
   3.659712830390296*^9, 3.659712974424768*^9, 3.6597130284844103`*^9, 
   3.6597131669543867`*^9, 3.6598169757951803`*^9, 3.659817063581402*^9, 
   3.659817253676845*^9, 3.6598173392135973`*^9, {3.659817422203857*^9, 
   3.659817436725481*^9}, 3.65981752229186*^9, 3.659817626684482*^9, 
   3.659817718156196*^9, 3.659818022532641*^9, {3.659818134510556*^9, 
   3.659818150115535*^9}, 3.6598182069965677`*^9, 3.659818239825873*^9, 
   3.659819022912571*^9, 3.6598195633588743`*^9}]
}, Open  ]]
},
WindowSize->{1855, 1056},
WindowMargins->{{-7, Automatic}, {Automatic, 0}},
FrontEndVersion->"10.1 for Linux x86 (64-bit) (March 23, 2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 173, 2, 31, "Text"],
Cell[734, 24, 715, 13, 55, "Input"],
Cell[1452, 39, 413, 11, 57, "Input"],
Cell[1868, 52, 992, 20, 32, "Input"],
Cell[2863, 74, 974, 20, 32, "Input"],
Cell[3840, 96, 2322, 48, 55, "Input"],
Cell[CellGroupData[{
Cell[6187, 148, 140, 2, 32, "Input"],
Cell[6330, 152, 470, 7, 32, "Output"]
}, Open  ]],
Cell[6815, 162, 7743, 187, 368, "Input"],
Cell[CellGroupData[{
Cell[14583, 353, 316, 6, 32, "Input"],
Cell[14902, 361, 16712, 358, 471, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31651, 724, 71, 1, 32, "Input"],
Cell[31725, 727, 2631, 59, 457, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)