From 7d14d9213e5c55786486ce45ec52f5b2a3e46322 Mon Sep 17 00:00:00 2001 From: Yair Fosado <yair.fosado@ed.ac.uk> Date: Sun, 15 Oct 2023 09:31:29 +0100 Subject: [PATCH] Upload New File --- Notebooks/fhLMFT.nb | 1200 +++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 1200 insertions(+) create mode 100644 Notebooks/fhLMFT.nb diff --git a/Notebooks/fhLMFT.nb b/Notebooks/fhLMFT.nb new file mode 100644 index 0000000..d8777bf --- /dev/null +++ b/Notebooks/fhLMFT.nb @@ -0,0 +1,1200 @@ +(* Content-type: application/vnd.wolfram.mathematica *) + +(*** Wolfram Notebook File ***) +(* http://www.wolfram.com/nb *) + +(* CreatedBy='Mathematica 13.0' *) + +(*CacheID: 234*) +(* Internal cache information: +NotebookFileLineBreakTest +NotebookFileLineBreakTest +NotebookDataPosition[ 158, 7] +NotebookDataLength[ 52421, 1192] +NotebookOptionsPosition[ 51895, 1174] +NotebookOutlinePosition[ 52332, 1191] +CellTagsIndexPosition[ 52289, 1188] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ + +Cell[CellGroupData[{ +Cell[BoxData[ + RowBox[{ + RowBox[{"(*", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**"}], "*", "**"}], "**"}], "**"}], "**"}], "**"}], "**"}], + "**"}], "**"}], "**"}], "******)"}], "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "ANALYSIS", " ", "FOR", " ", "TWO", " ", "POLYMER", " ", "SPECIES", " ", + "AND", " ", "NO", " ", "SOLVENT"}], "*)"}], "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**"}], "*", "**"}], "**"}], "**"}], "**"}], "**"}], "**"}], + "**"}], "**"}], "**"}], "******)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"fmix", "[", + RowBox[{"\[CapitalPhi]_", ",", "\[Chi]_", ",", "M_"}], "]"}], ":=", + RowBox[{ + RowBox[{ + FractionBox["\[CapitalPhi]", "M"], "*", + RowBox[{"Log", "[", "\[CapitalPhi]", "]"}]}], "+", + RowBox[{ + FractionBox[ + RowBox[{"(", + RowBox[{"1", "-", "\[CapitalPhi]"}], ")"}], "M"], "*", + RowBox[{"Log", "[", + RowBox[{"1", "-", "\[CapitalPhi]"}], "]"}]}], "+", + RowBox[{"\[Chi]", "*", "\[CapitalPhi]", "*", + RowBox[{"(", + RowBox[{"1", "-", "\[CapitalPhi]"}], ")"}]}]}]}], ";"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "find", " ", "local", " ", "minima", " ", "and", " ", "maximum", " ", + "of", " ", "the", " ", "free", " ", "energy"}], "*)"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", "\[CapitalPhi]"], " ", + RowBox[{"(", + RowBox[{"fmix", "[", + RowBox[{"\[CapitalPhi]", ",", "\[Chi]", ",", "M"}], "]"}], ")"}]}], + ";"}], "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + RowBox[{ + RowBox[{"We", " ", "now", " ", "that", " ", "for", " ", "\[Chi]"}], "<", + SubscriptBox["\[Chi]", "c"]}], ",", " ", + RowBox[{ + RowBox[{ + "there", " ", "is", " ", "always", " ", "a", " ", "minimum", " ", + "located", " ", "at", " ", "\[CapitalPhi]min"}], "=", + RowBox[{"0.5", ".", " ", "However"}]}], ",", " ", + RowBox[{ + RowBox[{"above", " ", + SubscriptBox["\[Chi]", "c"]}], "=", + RowBox[{ + RowBox[{"2", "/", "M"}], " ", "more", " ", "minima", " ", + "develop"}]}]}], "*)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"der", "[", + RowBox[{"\[CapitalPhi]_", ",", "\[Chi]_", ",", "M_"}], "]"}], "=", + RowBox[{ + RowBox[{ + FractionBox["1", "M"], "*", + RowBox[{"Log", "[", + FractionBox["\[CapitalPhi]", + RowBox[{"1", "-", "\[CapitalPhi]"}]], "]"}]}], "+", + RowBox[{"\[Chi]", "*", + RowBox[{"(", + RowBox[{"1", "-", + RowBox[{"2", "*", "\[CapitalPhi]"}]}], ")"}]}]}]}], ";"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + RowBox[{ + RowBox[{ + "Below", " ", "we", " ", "find", " ", "the", " ", "minimum", " ", "of", + " ", "fmix", " ", "for", " ", "a", " ", "given", " ", "value", " ", + "of", " ", "M"}], "=", "M1"}], ",", " ", + RowBox[{"and", " ", "changing", " ", "\[Chi]", " ", + RowBox[{"in", " ", "[", + RowBox[{"0", ",", "3"}], "]"}], " ", "every", " ", "0.1"}]}], "*)"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"M1", "=", "1"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[Chi]lo", "=", "0.0"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[Chi]hi", "=", "5.0"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"delta\[Chi]", "=", "0.05"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Ndata", "=", + RowBox[{ + RowBox[{"Floor", "[", + RowBox[{"(", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]hi", "-", "\[Chi]lo"}], ")"}], "/", "delta\[Chi]"}], + ")"}], "]"}], "+", "1"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + SubscriptBox["\[Chi]", "c"], "=", + RowBox[{"2", "/", "M1"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"Initialize", " ", "lists"}], "*)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[CapitalPhi]min1", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[CapitalPhi]min2", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[CapitalPhi]min3", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"listNsol", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"fullList\[CapitalPhi]min1", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"fullList\[CapitalPhi]min2", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"fullListFmin", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"fullList\[Chi]", "=", + RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", + "\[IndentingNewLine]", + RowBox[{"For", "[", + RowBox[{ + RowBox[{"n", "=", "0"}], ",", + RowBox[{"n", "<", "Ndata"}], ",", + RowBox[{"n", "++"}], ",", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "The", " ", "value", " ", "of", " ", "\[Chi]", " ", "at", " ", "which", + " ", "we", " ", "will", " ", "find", " ", "the", " ", "minima"}], + "*)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"\[Chi]1", "=", + RowBox[{"\[Chi]lo", "+", + RowBox[{"n", "*", "delta\[Chi]"}]}]}], ";", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + RowBox[{ + RowBox[{"We", " ", "will", " ", "solve", " ", "der"}], "==", + RowBox[{ + "0", " ", "in", " ", "the", " ", "interval", " ", "\[CapitalPhi]", + " ", + RowBox[{ + RowBox[{"in", " ", "[", + RowBox[{"0", ",", "1"}], "]"}], ".", " ", "To"}], " ", "do", " ", + "that"}]}], ",", + RowBox[{ + "ask", " ", "Plot", " ", "to", " ", "plot", " ", "the", " ", + "derivative"}], ",", + RowBox[{ + "then", " ", "extract", " ", "the", " ", "list", " ", "of", " ", + "data", " ", + RowBox[{"points", ":"}]}]}], "*)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"(", + RowBox[{"points", "=", + RowBox[{ + RowBox[{"Cases", "[", + RowBox[{ + RowBox[{"Normal", "[", + RowBox[{"Plot", "[", + RowBox[{ + RowBox[{"der", "[", + RowBox[{"\[CapitalPhi]", ",", "\[Chi]1", ",", "M1"}], "]"}], + ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", ",", "0", ",", "1"}], "}"}]}], "]"}], + "]"}], ",", + RowBox[{ + RowBox[{"Line", "[", "l_", "]"}], ":>", "l"}], ",", "Infinity"}], + "]"}], "//", "First"}]}], ")"}], "//", "Shallow"}], ";", + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "Now", " ", "we", " ", "find", " ", "which", " ", "of", " ", "these", + " ", "points", " ", "change", " ", "the", " ", "sign", " ", "of", " ", + "der"}], "*)"}], "\[IndentingNewLine]", + RowBox[{"signChanges", "=", + RowBox[{ + RowBox[{"Position", "[", + RowBox[{ + RowBox[{ + RowBox[{ + RowBox[{"Times", "@@", "##"}], "&"}], "/@", + RowBox[{"Partition", "[", + RowBox[{ + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", "2", ",", "1"}], + "]"}]}], ",", + RowBox[{"x_", "/;", + RowBox[{"x", "<", "0"}]}]}], "]"}], "//", "Flatten"}]}], ";", + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"Now", " ", "use", " ", "find", " ", "root"}], "*)"}], + "\[IndentingNewLine]", + RowBox[{"\[CapitalPhi]min", "=", + RowBox[{ + RowBox[{ + RowBox[{"\[CapitalPhi]", "/.", + RowBox[{"FindRoot", "[", + RowBox[{ + RowBox[{ + RowBox[{"der", "[", + RowBox[{"\[CapitalPhi]", ",", "\[Chi]1", ",", "M1"}], "]"}], "==", + "0"}], ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{"#", ",", "1"}], "]"}], "]"}], ",", + RowBox[{"points", "[", + RowBox[{"[", + RowBox[{ + RowBox[{"#", "+", "1"}], ",", "1"}], "]"}], "]"}]}], "}"}]}], + "]"}]}], "&"}], "/@", "signChanges"}]}], ";", + "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{"\[CapitalPhi]", "->", "0."}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", "->", "1."}], "}"}]}], "}"}], ";", + "\[IndentingNewLine]", "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"fullList\[CapitalPhi]min1", ",", + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", + + RowBox[{"AppendTo", "[", + RowBox[{"fullList\[CapitalPhi]min2", ",", + RowBox[{"1", "-", + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], ";", + "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"fullListFmin", ",", + RowBox[{"fmix", "[", + RowBox[{ + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "1", "]"}], "]"}], ",", "\[Chi]1", ",", "M1"}], + "]"}]}], "]"}], ";", "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"fullList\[Chi]", ",", "\[Chi]1"}], "]"}], ";", + "\[IndentingNewLine]", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"Above", " ", + SubscriptBox["\[Chi]", "c"], " ", "find", " ", "the", " ", "minima", + " ", "and", " ", "maximum"}], "*)"}], "\[IndentingNewLine]", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]1", ">", + RowBox[{"2", "/", "M1"}]}], ",", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"number", " ", "of", " ", "solutions", " ", "found"}], + "*)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Nsol", "=", + RowBox[{"Length", "[", "\[CapitalPhi]min", "]"}]}], ";", + "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"listNsol", ",", "Nsol"}], "]"}], ";", + "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"\[CapitalPhi]min1", ",", + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ";", + "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"\[CapitalPhi]min2", ",", + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "2", "]"}], "]"}]}], "]"}], ";", + "\[IndentingNewLine]", + RowBox[{"AppendTo", "[", + RowBox[{"\[CapitalPhi]min3", ",", + RowBox[{"\[CapitalPhi]min", "[", + RowBox[{"[", "3", "]"}], "]"}]}], "]"}], ";"}], + "\[IndentingNewLine]", ",", "\[IndentingNewLine]", + RowBox[{"Continue", "[", "]"}]}], "]"}]}]}], + RowBox[{"(*", + RowBox[{ + "Close", " ", "if", " ", "statement", " ", "on", " ", "\[Chi]1"}], + "*)"}], "\[IndentingNewLine]", "]"}], + RowBox[{"(*", + RowBox[{"close", " ", "for", " ", "loop"}], "*)"}], "\[IndentingNewLine]", + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**", "**", "**", "**", "**", "**", "**"}], "******)"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "Plot", " ", "the", " ", "values", " ", "of", " ", "\[CapitalPhi]min", + " ", "as", " ", "function", " ", "of", " ", "\[Chi]"}], "*)"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**", "**", "**", "**", "**", "**", "**"}], "******)"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"data1", "=", + RowBox[{"Transpose", "@", + RowBox[{"{", + RowBox[{"fullList\[Chi]", ",", "fullList\[CapitalPhi]min1"}], + "}"}]}]}], ";"}], "\n", + RowBox[{ + RowBox[{"data2", "=", + RowBox[{"Transpose", "@", + RowBox[{"{", + RowBox[{"fullList\[Chi]", ",", "fullList\[CapitalPhi]min2"}], + "}"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**", "**", "**"}], "*****)"}], "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "Plot", " ", "the", " ", "values", " ", "of", " ", "\[Chi]", " ", "as", + " ", "function", " ", "\[CapitalPhi]min"}], "*)"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{ + "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", "**", + "**", "**", "**"}], "*****)"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + RowBox[{"\[Chi]binodal", "[", + RowBox[{"\[CapitalPhi]_", ",", "M_"}], "]"}], ":=", + RowBox[{ + FractionBox["1", + RowBox[{"M", "*", + RowBox[{"(", + RowBox[{ + RowBox[{"2", "*", "\[CapitalPhi]"}], "-", "1"}], ")"}]}]], "*", + RowBox[{"Log", "[", + FractionBox["\[CapitalPhi]", + RowBox[{"1", "-", "\[CapitalPhi]"}]], "]"}]}]}], ";"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Plot", "[", + RowBox[{ + RowBox[{"\[Chi]binodal", "[", + RowBox[{"\[CapitalPhi]", ",", "M1"}], "]"}], ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", ",", "0", ",", "1.0"}], "}"}]}], "]"}], ";"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{ + SubscriptBox["\[PartialD]", "\[CapitalPhi]"], " ", + RowBox[{"(", + RowBox[{"\[Chi]binodal", "[", + RowBox[{"\[CapitalPhi]", ",", "M"}], "]"}], ")"}]}], ";"}], + "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Limit", "[", + RowBox[{ + RowBox[{"\[Chi]binodal", "[", + RowBox[{"\[CapitalPhi]", ",", "M1"}], "]"}], ",", + RowBox[{"x", "->", "0.5"}]}], "]"}], ";"}], "\[IndentingNewLine]", + "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"**", "**", "**", "**", "**", "**"}], "****)"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"Merge", " ", "all", " ", "plots"}], "*)"}], + "\[IndentingNewLine]", + RowBox[{"(*", + RowBox[{"**", "**", "**", "**", "**", "**"}], "****)"}], + "\[IndentingNewLine]", "\[IndentingNewLine]", + RowBox[{"Manipulate", "[", + RowBox[{ + RowBox[{"Row", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"Plot", "[", + RowBox[{ + RowBox[{"fmix", "[", + RowBox[{"\[CapitalPhi]", ",", "\[Chi]", ",", "M1"}], "]"}], ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", ",", "0", ",", "1.0"}], "}"}], ",", + RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", " ", + RowBox[{"AxesLabel", "\[Rule]", + RowBox[{"{", + RowBox[{ + "\"\<\[CapitalPhi]\>\"", ",", + "\"\<\!\(\*SubscriptBox[\(F\), \(mix\)]\)\>\""}], "}"}]}], ",", + RowBox[{"PlotRange", "\[Rule]", + RowBox[{"{", + RowBox[{"Automatic", ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.71"}], ",", "0.57"}], "}"}]}], "}"}]}], ",", + RowBox[{"PlotLegends", "->", + RowBox[{"Placed", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<", + RowBox[{"2", "/", "M1"}]}], ",", " ", + "\"\<\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "==", + RowBox[{"2", "/", "M1"}]}], ",", + "\"\<\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", + ",", " ", + "\"\<\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\""}], + "]"}]}], "]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0.9", ",", "0.9"}], "}"}]}], "]"}]}], ",", + RowBox[{"Epilog", "->", + RowBox[{"{", + RowBox[{ + RowBox[{"PointSize", "[", "Large", "]"}], ",", "Red", ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<=", + RowBox[{"2", "/", "M1"}]}], ",", + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{"0.5", ",", + RowBox[{"fmix", "[", + RowBox[{"0.5", ",", "\[Chi]", ",", "M1"}], "]"}]}], "}"}], + "]"}], ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", ">", + RowBox[{"2", "/", "M1"}]}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fmix", "[", + RowBox[{ + RowBox[{"\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", "\[Chi]", ",", "M1"}], "]"}]}], + "}"}], "]"}], ",", + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"\[CapitalPhi]min3", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fmix", "[", + RowBox[{ + RowBox[{"\[CapitalPhi]min3", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", "\[Chi]", ",", "M1"}], "]"}]}], + "}"}], "]"}]}], "}"}]}], " ", "]"}]}], "]"}], ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<=", + RowBox[{"2", "/", "M1"}]}], ",", + RowBox[{"{", + RowBox[{"Text", "[", + RowBox[{ + RowBox[{"Style", "[", + RowBox[{ + "\"\<\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)\>\"", + ",", "Red"}], "]"}], ",", + RowBox[{"{", + RowBox[{"0.5", ",", + RowBox[{"-", "0.05"}]}], "}"}]}], "]"}], " ", "}"}], ",", + " ", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", ">", + RowBox[{"2", "/", "M1"}]}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"Text", "[", + RowBox[{ + RowBox[{"Style", "[", + RowBox[{ + "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \ +\(1\)], \(*\)]\)\>\"", ",", "Red"}], "]"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"-", "0.05"}]}], "}"}]}], "]"}], ",", + RowBox[{"Text", "[", + RowBox[{ + RowBox[{"Style", "[", + RowBox[{ + "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \ +\(2\)], \(*\)]\)\>\"", ",", "Red"}], "]"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"\[CapitalPhi]min3", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"-", "0.05"}]}], "}"}]}], "]"}]}], "}"}]}], " ", + "]"}]}], "]"}]}], "\[IndentingNewLine]", "}"}]}], ",", + RowBox[{"ImageSize", "\[Rule]", "400"}]}], "]"}], " ", ",", + "\"\< \>\"", ",", + RowBox[{"ListPlot", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"data1", ",", "data2"}], "}"}], ",", + RowBox[{"PlotStyle", "\[Rule]", "Black"}], ",", " ", + RowBox[{"AxesLabel", "\[Rule]", + RowBox[{"{", + RowBox[{ + "\"\<\[Chi]\>\"", ",", + "\"\<\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)\>\""}], + "}"}]}], ",", + RowBox[{"PlotRange", "\[Rule]", + RowBox[{"{", + RowBox[{"Automatic", ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.05"}], ",", "1.05"}], "}"}]}], "}"}]}], ",", + RowBox[{"PlotLegends", "->", + RowBox[{"Placed", "[", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<", + RowBox[{"2", "/", "M1"}]}], ",", " ", + "\"\<\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", + ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "==", + RowBox[{"2", "/", "M1"}]}], ",", + "\"\<\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", + ",", "\"\<\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\ +\""}], "]"}]}], "]"}], ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<=", + RowBox[{"2", "/", "M1"}]}], ",", + "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \ +\(1\)], \(*\)]\)=0.5\\n\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \ +\(2\)], \(*\)]\)=0.5\>\"", ",", + RowBox[{ + "\"\<\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \ +\(1\)], \(*\)]\)=\>\"", "<>", + RowBox[{"ToString", "@", + RowBox[{"SetPrecision", "[", + RowBox[{ + RowBox[{"\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", "2"}], "]"}]}], "<>", + "\"\<\\n\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\),\ + \(2\)], \(*\)]\)=\>\"", "<>", + RowBox[{"ToString", "@", + RowBox[{"SetPrecision", "[", + RowBox[{ + RowBox[{"\[CapitalPhi]min3", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "-", + FractionBox["2", "M1"]}], ")"}], "/", "delta\[Chi]"}], + "]"}], "]"}], "]"}], ",", "4"}], "]"}]}], "<>", + "\"\<\>\""}]}], "]"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"0.15", ",", "0.8"}], "}"}]}], "]"}]}], ",", + RowBox[{"Epilog", "->", + RowBox[{"{", + RowBox[{ + RowBox[{"PointSize", "[", "Large", "]"}], ",", "Red", ",", + RowBox[{"{", + RowBox[{ + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"fullList\[Chi]", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fullList\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}]}], "}"}], "]"}], ",", + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"fullList\[Chi]", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fullList\[CapitalPhi]min2", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}]}], "}"}], "]"}]}], " ", "}"}]}], + "}"}]}], ",", + RowBox[{"ImageSize", "->", "400"}]}], "]"}], ",", "\"\< \>\"", + ",", + RowBox[{"Plot", "[", + RowBox[{ + RowBox[{"\[Chi]binodal", "[", + RowBox[{"\[CapitalPhi]", ",", "M1"}], "]"}], ",", + RowBox[{"{", + RowBox[{"\[CapitalPhi]", ",", "0", ",", "1.0"}], "}"}], ",", + RowBox[{"PlotStyle", "\[Rule]", "Black"}], ",", " ", + RowBox[{"AxesLabel", "\[Rule]", + RowBox[{"{", + RowBox[{ + "\"\<\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)\>\"", ",", + "\"\<\[Chi]\>\""}], "}"}]}], ",", + RowBox[{"PlotRange", "\[Rule]", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.05"}], ",", "1.05"}], "}"}], ",", + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.1"}], ",", + RowBox[{"\[Chi]hi", "+", "0.05"}]}], "}"}]}], "}"}]}], ",", + RowBox[{"PlotLegends", "->", + RowBox[{"Placed", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "<", + RowBox[{"2", "/", "M1"}]}], ",", " ", + "\"\<\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", ",", + RowBox[{"If", "[", + RowBox[{ + RowBox[{"\[Chi]", "==", + RowBox[{"2", "/", "M1"}]}], ",", + "\"\<\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\"", + ",", " ", + "\"\<\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)\>\""}], + "]"}]}], "]"}], "}"}], ",", + RowBox[{"{", + RowBox[{"0.5", ",", "0.9"}], "}"}]}], "]"}]}], ",", + RowBox[{"Epilog", "->", + RowBox[{"{", + RowBox[{ + RowBox[{"PointSize", "[", "Large", "]"}], ",", "Red", ",", + RowBox[{"{", + RowBox[{ + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"fullList\[CapitalPhi]min1", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fullList\[Chi]", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}]}], "}"}], "]"}], ",", + RowBox[{"Point", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"fullList\[CapitalPhi]min2", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}], ",", + RowBox[{"fullList\[Chi]", "[", + RowBox[{"[", + RowBox[{"Round", "[", + RowBox[{ + RowBox[{"(", + RowBox[{"\[Chi]", "/", "delta\[Chi]"}], ")"}], "+", "1"}], + "]"}], "]"}], "]"}]}], "}"}], "]"}]}], " ", "}"}]}], + "}"}]}], ",", + RowBox[{"ImageSize", "->", "400"}]}], "]"}]}], + "\[IndentingNewLine]", "}"}], "]"}], ",", + RowBox[{"{", + RowBox[{ + "\[Chi]", ",", "\[Chi]lo", ",", "\[Chi]hi", ",", "delta\[Chi]", ",", + RowBox[{"Appearance", "->", "\"\<Labeled\>\""}]}], "}"}], ",", + RowBox[{"SaveDefinitions", "\[Rule]", "True"}]}], "]"}], + "\[IndentingNewLine]", "\[IndentingNewLine]"}]}]], "Input", + CellChangeTimes->{{3.905953446829336*^9, 3.905953492125682*^9}, { + 3.9059565843625603`*^9, 3.905956585652594*^9}, {3.905958734391288*^9, + 3.9059587652377853`*^9}, {3.905958815678087*^9, 3.905958897740513*^9}, { + 3.905958929934346*^9, 3.905958942215672*^9}, {3.905958983211403*^9, + 3.905958989100917*^9}, 3.905970817605957*^9, 3.906008275263885*^9, { + 3.906185709523773*^9, 3.9061859379602222`*^9}, {3.906186013052332*^9, + 3.906186014192399*^9}, {3.9061860651259403`*^9, 3.906186314734708*^9}, { + 3.9061864679879093`*^9, 3.906186469545937*^9}, {3.906186520020146*^9, + 3.906186523397768*^9}, {3.906186624892364*^9, 3.906186659071211*^9}, { + 3.906186715929151*^9, 3.906186839377307*^9}, {3.9061869001825333`*^9, + 3.906186913981937*^9}, {3.906187057092443*^9, 3.906187444506194*^9}, { + 3.906187968429184*^9, 3.906187986498597*^9}, {3.906189150233098*^9, + 3.9061891507997913`*^9}, {3.906191972901394*^9, 3.90619217777065*^9}, { + 3.906192234535191*^9, 3.906192245882135*^9}, {3.906192282955419*^9, + 3.9061923907181168`*^9}, {3.9061924247654247`*^9, + 3.9061924301124077`*^9}, {3.90619255628575*^9, 3.906192562752557*^9}, { + 3.9061926976438923`*^9, 3.906192707781095*^9}, {3.90619273841507*^9, + 3.906192793193458*^9}, {3.90619289192037*^9, 3.906193007366138*^9}, { + 3.906193157959475*^9, 3.906193199355279*^9}, {3.906193252022624*^9, + 3.9061933095631447`*^9}, {3.906193345309771*^9, 3.906193386248043*^9}, { + 3.906193422158856*^9, 3.906193447834227*^9}, {3.9061934828694143`*^9, + 3.906193512125957*^9}, {3.906193550870611*^9, 3.906193619320952*^9}, { + 3.906194760665752*^9, 3.9061950056828623`*^9}, {3.906195061852645*^9, + 3.9061952375640373`*^9}, {3.906195271121324*^9, 3.9061952726558104`*^9}, { + 3.906195320935075*^9, 3.906195349846089*^9}, {3.90619538752774*^9, + 3.906195431747758*^9}, {3.9061956221278954`*^9, 3.906195657691861*^9}, { + 3.906195694287621*^9, 3.906195717274535*^9}, {3.906195750594314*^9, + 3.906195910726684*^9}, {3.906195984383193*^9, 3.906196054006936*^9}, { + 3.906196103292314*^9, 3.906196148773584*^9}, {3.906196236312628*^9, + 3.9061962870400667`*^9}, {3.906196344362204*^9, 3.906196371610655*^9}, { + 3.9061967073557367`*^9, 3.906196755785902*^9}, {3.906196810261244*^9, + 3.906196851540406*^9}, {3.906197247135024*^9, 3.906197320278891*^9}, { + 3.906197369961841*^9, 3.906197373199267*^9}, {3.9061974496224747`*^9, + 3.9061974557970533`*^9}, {3.906197496724215*^9, 3.906197567619028*^9}, { + 3.9061976895323267`*^9, 3.906197718988682*^9}, {3.906197759996933*^9, + 3.906197841393283*^9}, {3.9061979263434877`*^9, 3.9061979873695707`*^9}, { + 3.906198020047338*^9, 3.90619802186659*^9}, {3.906198066715631*^9, + 3.906198094341017*^9}, {3.906198141310308*^9, 3.90619817160544*^9}, { + 3.906198244562826*^9, 3.90619825042142*^9}, {3.906198281177109*^9, + 3.906198300289432*^9}, {3.90619837789198*^9, 3.906198422525544*^9}, { + 3.906198497978826*^9, 3.906198498329142*^9}, {3.9061986008503313`*^9, + 3.906198669710833*^9}, {3.906198796518795*^9, 3.906198904104203*^9}, { + 3.906199265626696*^9, 3.906199371046061*^9}, {3.9061994040209503`*^9, + 3.9061994727403812`*^9}, 3.90619957360422*^9, {3.906199628433597*^9, + 3.906199730597781*^9}, {3.9061998670677757`*^9, 3.906199889280673*^9}, { + 3.906199938170035*^9, 3.9061999517432413`*^9}, {3.906201951374099*^9, + 3.906201951714698*^9}, {3.906202326409339*^9, 3.906202414491892*^9}, { + 3.906202453586145*^9, 3.906202609460306*^9}, {3.9062026661628027`*^9, + 3.9062026737128468`*^9}, 3.9062084158027678`*^9, {3.906209293147564*^9, + 3.90620929390555*^9}, 3.90620977065466*^9, {3.906211219133668*^9, + 3.906211229647726*^9}, {3.906211267877417*^9, 3.90621127967336*^9}, { + 3.906213645758916*^9, 3.906213646078842*^9}, {3.906213717222803*^9, + 3.906213719562325*^9}, {3.9062137700123177`*^9, 3.906213770232682*^9}, { + 3.9062138890405493`*^9, 3.906213951744598*^9}, {3.9062139970542927`*^9, + 3.906213997275508*^9}, {3.906214066575585*^9, 3.90621406941665*^9}, { + 3.9062141052473507`*^9, 3.9062141283889437`*^9}, {3.906214159141612*^9, + 3.906214191919896*^9}, {3.90621426039881*^9, 3.906214260542759*^9}, { + 3.906214336934362*^9, 3.9062143589742603`*^9}, {3.906214389827056*^9, + 3.9062143908259277`*^9}}, + CellLabel-> + "In[320]:=",ExpressionUUID->"cef6a996-4bad-4901-8d11-377a6ad0b7cd"], + +Cell[BoxData[ + TagBox[ + StyleBox[ + DynamicModuleBox[{$CellContext`\[Chi]$$ = 0., Typeset`show$$ = True, + Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", + Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = + "\"untitled\"", Typeset`specs$$ = {{ + Hold[$CellContext`\[Chi]$$], 0., 5., 0.05}}, Typeset`size$$ = { + 1270., {130., 135.}}, Typeset`update$$ = 0, Typeset`initDone$$, + Typeset`skipInitDone$$ = False}, + DynamicBox[Manipulate`ManipulateBoxes[ + 1, StandardForm, "Variables" :> {$CellContext`\[Chi]$$ = 0.}, + "ControllerVariables" :> {}, + "OtherVariables" :> { + Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, + Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, + Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, + Typeset`skipInitDone$$}, "Body" :> Row[{ + Plot[ + $CellContext`fmix[$CellContext`\[CapitalPhi], \ +$CellContext`\[Chi]$$, $CellContext`M1], {$CellContext`\[CapitalPhi], 0, 1.}, + PlotStyle -> Blue, + AxesLabel -> { + "\[CapitalPhi]", "\!\(\*SubscriptBox[\(F\), \(mix\)]\)"}, + PlotRange -> {Automatic, {-0.71, 0.57}}, PlotLegends -> Placed[{ + + If[$CellContext`\[Chi]$$ < 2/$CellContext`M1, + "\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + + If[$CellContext`\[Chi]$$ == 2/$CellContext`M1, + "\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + "\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)"]]}, {0.9, + 0.9}], Epilog -> { + PointSize[Large], Red, + If[$CellContext`\[Chi]$$ <= 2/$CellContext`M1, + Point[{0.5, + $CellContext`fmix[ + 0.5, $CellContext`\[Chi]$$, $CellContext`M1]}], + If[$CellContext`\[Chi]$$ > 2/$CellContext`M1, { + Point[{ + Part[$CellContext`\[CapitalPhi]min1, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], + $CellContext`fmix[ + Part[$CellContext`\[CapitalPhi]min1, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], \ +$CellContext`\[Chi]$$, $CellContext`M1]}], + Point[{ + Part[$CellContext`\[CapitalPhi]min3, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], + $CellContext`fmix[ + Part[$CellContext`\[CapitalPhi]min3, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], \ +$CellContext`\[Chi]$$, $CellContext`M1]}]}]], + If[$CellContext`\[Chi]$$ <= 2/$CellContext`M1, { + Text[ + Style[ + "\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)", Red], { + 0.5, -0.05}]}, + If[$CellContext`\[Chi]$$ > 2/$CellContext`M1, { + Text[ + Style[ + "\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(1\)], \ +\(*\)]\)", Red], { + Part[$CellContext`\[CapitalPhi]min1, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], -0.05}], + Text[ + Style[ + "\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(2\)], \ +\(*\)]\)", Red], { + Part[$CellContext`\[CapitalPhi]min3, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], \ +-0.05}]}]]}, ImageSize -> 400], " ", + ListPlot[{$CellContext`data1, $CellContext`data2}, PlotStyle -> + Black, AxesLabel -> { + "\[Chi]", "\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)"}, + PlotRange -> {Automatic, {-0.05, 1.05}}, PlotLegends -> Placed[{ + + If[$CellContext`\[Chi]$$ < 2/$CellContext`M1, + "\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + + If[$CellContext`\[Chi]$$ == 2/$CellContext`M1, + "\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + "\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)"]], + + If[$CellContext`\[Chi]$$ <= 2/$CellContext`M1, + "\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(1\)], \ +\(*\)]\)=0.5\n\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(2\)], \ +\(*\)]\)=0.5", + "\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(1\)], \ +\(*\)]\)=" <> ToString[ + SetPrecision[ + Part[$CellContext`\[CapitalPhi]min1, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], 2]] <> + "\n\!\(\*SuperscriptBox[SubscriptBox[\(\[CapitalPhi]\), \(2\)], \ +\(*\)]\)=" <> ToString[ + SetPrecision[ + Part[$CellContext`\[CapitalPhi]min3, + + Round[($CellContext`\[Chi]$$ - + 2/$CellContext`M1)/$CellContext`delta\[Chi]]], 4]] <> + ""]}, {0.15, 0.8}], Epilog -> { + PointSize[Large], Red, { + Point[{ + Part[$CellContext`fullList\[Chi], + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]], + Part[$CellContext`fullList\[CapitalPhi]min1, + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]]}], + + Point[{ + Part[$CellContext`fullList\[Chi], + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]], + Part[$CellContext`fullList\[CapitalPhi]min2, + + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + + 1]]}]}}, ImageSize -> 400], " ", + Plot[ + $CellContext`\[Chi]binodal[$CellContext`\[CapitalPhi], \ +$CellContext`M1], {$CellContext`\[CapitalPhi], 0, 1.}, PlotStyle -> Black, + AxesLabel -> { + "\!\(\*SuperscriptBox[\(\[CapitalPhi]\), \(*\)]\)", "\[Chi]"}, + PlotRange -> {{-0.05, 1.05}, {-0.1, $CellContext`\[Chi]hi + 0.05}}, + PlotLegends -> Placed[{ + If[$CellContext`\[Chi]$$ < 2/$CellContext`M1, + "\[Chi]<\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + + If[$CellContext`\[Chi]$$ == 2/$CellContext`M1, + "\[Chi]=\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)", + "\[Chi]>\!\(\*SubscriptBox[\(\[Chi]\), \(c\)]\)"]]}, {0.5, + 0.9}], Epilog -> { + PointSize[Large], Red, { + Point[{ + Part[$CellContext`fullList\[CapitalPhi]min1, + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]], + Part[$CellContext`fullList\[Chi], + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]]}], + + Point[{ + Part[$CellContext`fullList\[CapitalPhi]min2, + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + 1]], + Part[$CellContext`fullList\[Chi], + + Round[$CellContext`\[Chi]$$/$CellContext`delta\[Chi] + + 1]]}]}}, ImageSize -> 400]}], + "Specifications" :> {{$CellContext`\[Chi]$$, 0., 5., 0.05, Appearance -> + "Labeled"}}, "Options" :> {}, "DefaultOptions" :> {}], + ImageSizeCache->{1314., {175., 181.}}, + SingleEvaluation->True], + Deinitialization:>None, + DynamicModuleValues:>{}, + Initialization:>({$CellContext`fmix[ + Pattern[$CellContext`\[CapitalPhi], + Blank[]], + Pattern[$CellContext`\[Chi], + Blank[]], + Pattern[$CellContext`M, + Blank[]]] := ($CellContext`\[CapitalPhi]/$CellContext`M) + Log[$CellContext`\[CapitalPhi]] + (( + 1 - $CellContext`\[CapitalPhi])/$CellContext`M) + Log[1 - $CellContext`\[CapitalPhi]] + $CellContext`\[Chi] \ +$CellContext`\[CapitalPhi] ( + 1 - $CellContext`\[CapitalPhi]), $CellContext`\[CapitalPhi][ + Pattern[$CellContext`m, + Blank[]], + Pattern[$CellContext`h, + Blank[]], + Pattern[$CellContext`T, + Blank[]]] := ($CellContext`r[$CellContext`T] $CellContext`m^2)/ + 2 + $CellContext`u $CellContext`m^4 - $CellContext`h $CellContext`m, \ +$CellContext`r[ + Pattern[$CellContext`T, + + Blank[]]] := $CellContext`a ($CellContext`T - $CellContext`Tc), \ +$CellContext`a = 1., $CellContext`Tc = 1., $CellContext`u = + 0.1, $CellContext`M1 = + 1, $CellContext`\[CapitalPhi]min1 = CompressedData[" +1:eJwB7QES/iFib1JlAQAAADwAAAAZve09xG3XP5qERbstI9Q/0jJSVw7G0T8r +ArJV0s/PPxSSeRdRucw/rrQpCiIayj98/c/E+dXHP2/YmqP+2cU/bNOqRtYY +xD8ngM4DnYjCP+YIcry/IcE/i7LtSpO8vz8ctRaI5HK9P0cNb3aiXrs/feCN +BVN5uT9JRmRyh723P1tFuxmlJrY/v6YvLrywtD8GIWY8aFizP9aHe8O3GrI/ +3BKxAxn1sD9MyzRplcqvPxtia2ef0q0/0SAy+8n+qz9qPXzi8kuqP6StkvpN +t6g/tK0F3Fk+pz8ymgVI1t6lP+vZvw+8lqQ/MtT5MDZkoz9DUFfym0WiP6Nk +A9NrOaE/Nse0K0c+oD93KozI3KWeP0NTqU177Jw/5Olh+1NOmz9vvNeXfcmZ +PyPO6yc3XJg/tO01AuQElz/DdplVCMKVP0Vw/RNGkpQ/V8MRNFp0kz9cmP8+ +GmeSP9P7gyByaZE/aUJIMWJ6kD9BdBPp+jGPP8NC1QXQiI0/K84/Oqv3iz/q +uuDEEH2KP0SYQEyfF4k/b1SVwg3Ghz+zAxx7KYeGP6mSLWzUWYU/1FlWmQM9 +hD8fIFKhvS+DPzTxSmwZMYI/EnIl9zxAgT9aVOs4XFyAP2PnEEBwCX8/ekNu +STtxfT9FruHD + "], $CellContext`delta\[Chi] = + 0.05, $CellContext`\[CapitalPhi]min3 = CompressedData[" +1:eJwB7QES/iFib1JlAQAAADwAAAB6IQnhHUnkP8Q9XSJp7uU/meZW1Pgc5z+w +f5NqCwzoP3CbIbqr0eg/SJd1fXd56T+hAMyOgQrqP5pEGVeAieo/U4pUbsr5 +6j/+Xwy/2F3rP+R94xCQt+s/q0milm0I7D8pL/1uo1HsPyQRMrErlOw/8ENO +n9XQ7D83d7MRTwjtP71fyFwrO+0/JQs6euhp7T/fO3P48pTtP8GOkAepvO0/ +pd2J31zh7T9Ls2ypVgPuP95JiQnWIu4/HQxNYBNA7j8nPNjRQFvuP19sViCL +dO4/E7s/YhqM7j9Zpn+bEqLuPwxlBD+Utu4/9mTwnLzJ7j8Ai9pAptvuPz9X +z0Jp7O4/g7NEjRv87j/T2bsZ0QrvPxe5kCWcGO8/sfAkYI0l7z8HN0ETtDHv +P5bAvkYePe8/lFDu39hH7z9JNFO971HvP27gX89tW+8/53FfLl1k7z9aAwgu +x2zvP6/n+260dO8/7L117ix87z9islsUOIPvP2zV57/cie8/wQAXUyGQ7z8J +fOy8C5bvPybxzoKhm+8/q6r1yOeg7z8PkBNa46XvP9lJT66Yqu8/lKaa8Quv +7z+ct3oJQbPvP87TTpo7t+8/BHYjDP+67z/dUhyPjr7vP4zefx/twe8/eT5t +iR3F7z96Dvbh + "], $CellContext`data1 = CompressedData[" +1:eJxd0g1MzHEYB/BzL3W934uh2JSGpSUl79U9pFdl9CLlZpOSiONmLNXktbw0 +XFJX3nrT2yJ1tHSt6Y6U9C7amoUWSVaJ3ovj2W2eZ/vvt+/+z/N5fvvvbxEi +8Q1jMhiMsD+P5vy/ukRpck3li/7P5drssFJTKvK+Vpv/nU2kv1Wbw/dpqp3M +d2jz63pNdRKvi/R/IP4nMt9N9vUQ5wvxeonXR7xvxPtOvAHSP0D8QTI/RPb9 +IN4w8YaJ95N4v4g3QrxR4o0Sb4x448SbIN4k8SaJN0W8aeLNEI8Bw7l5cra2 +nwGmlf1O6kNv0WdA2jnXChvLNtHf8ZWzoMXeN9D4RTPumwXWzJKAlsZG9Jgw +WxY1d3t5PXpMeKCw0l9kVoceE2KmGtVjb2rwPiyQvEvLmu54jh4LIlsL3Nrn +qdFjw9K4BtbthGr02NDDlSqrLJ6hx4bLJf1eScoq9Diw4PGchG5pJXoccDeS +HM/cXYGeDsR2JXL8osrR0wFPt1BpfGUZejqwy7XCNGfJE/R0oeq+7XJlsQI9 +XdCx2Oe4f2cpelxoj49WlZmVoMeFTsEjlunPYvS44P3KISK57iF6ejA/5GD4 +3ZYi9PSg2dx+ona6ED19OOAU83FoSwF6+pBdJBv3KctDTx8URZzOIOdc9AzA +Po3j0f4+Bz0D6OuoMlHezEbPEOzbxlasDs1CzxBctgUOpbtmomcIWaGs1oPr +MtAzgtU1Cmt353voGcGJZVdedubcQc8YXPzyfaK/3ULPGLp790z4bU1Hzxgk +yrcpsbVy9EzAsqFv6RpxKnomoOhfxexmp6DHA9VxeQC3Ohn/Lx64HpkSuslu +oM+DQBXfIfhYEs7zQJw642wWLsN9PGidOGsujbiOHh8iNnnYhZy8hh4fNh4T +9o7bXUWPD6pNbzhNCYno8cG6QbS+cOQyenz4+rRLzYu9hJ4ANqeC9925F9ET +gMQ/udroRTx6AihlLYi2ir+AngDyZTaRbUHn0RNAW9AOOcvpHHpCWGjum1Fp +exa/jxAcBr0iTe3OoCcEgXTxiCOcRk8Iwf59a8XiOPSEsPczD47onUJvNpx0 +Oey54Wis6DeInUEj + "], $CellContext`data2 = CompressedData[" +1:eJxd0nss1XEUAPDret0X7sMWhbp0W6ZpYsuS3UOtdtMIwxZGFo2htZhRovrr +Vh6TEVqi2KiRrscSqbVCU1iYynPe3JdnhNSt79fmnO23785+53zOd7/fEUdc +9otkMhiMyL+P7twZI9LCAl1USHfmr7ZzF2ddvEfv27fz/2cXqv+6nV+K0kUf +6v+2nX/u0MUA8kZQ/Sjyx1D/OJo3iZxp5M0gbw55SuSpkadF9Vrkz6P+BTRv +EXlLyFtC3jLyVpD3E3mryFtF3hryfiFvHXkbyNtA3ibyfiNvC3kMsHhR1bGy +NUzqGZAgZo9ay8aJzwCFe6httHpC+q/dWQ/yJwN7Vq2myDw9SEy7H8Pl0f/K +hGtF4sZn3dPEY0JPTmJqUvIM8ZhQwujIvMmZJffRh3GZZVC6fJZ4+iDxCoj/ +tDZLPANYCeO19IfOEc8APqaO8bMb5ohnANWy0rw4lpJ4hhBivxn/xFdJPEMI +sjlS45BL98gItCe8i3u7lMQzgqNJCoEPS0U8IwjIbwtxcFMRzxgOco+lTEer +iGcMw25XVxdyVcRjgWFWtnFFk4p4LHg6KB8OGaUeC7wUsRWB+mrisWFIJmf3 +2aqJx4bFwbPhIqB7zQH74/3dEKwmHgcO9Abuu5ugJh4HJGXSiN33qMeF0PK0 +ImEp9bhwgWkgza2nHg+aI7QPm9qpx4ONO9+h/Af1eFDf3ukZraSeCWQoTmbs +2aCeCTx/9Nqym6MhnilE1eVIHlpoiGcKNdr94RkSDfFMwdrTQ1TrpCGeGWSV +vD1l464hnhnknFMP95+mHh/OuPg3a3w1ZL/4cMM8vDMumPp8mLgS5hh6kfbz +wV2f5dgaS+fxQTyyHl+bQD0BqJoT1YeuU08Ah18Gm7vepp4AliunWwbk1BNA +K2OXvzibegIwT1U2cfOoJ4TizS+3SoqoJ4SqyuW2qcfUE8J8luj8WBn1hFAn +86l6UEk9IRSUF85zq6kngvSGFLaHgn4fEbj1eBe6NVBPBK4xdrytRuqJoNrP +KivzDfVE4DSUtlf1jnrmkGwXJ7f+oJH+AbtXWW4= + + "], $CellContext`fullList\[Chi] = CompressedData[" +1:eJxF0ElKA1EUBdBK3yfVuAHXEOdVa3AJQnToenQHLsGxI0FBEATBgaKIokjU +qFFjQ/kPeCef/+55n6KW1zZXJ5UoiibRf7a3yuzk4dzNV8Zl9tz386COzI/z +jfUyJ/rT/PCgzBl3bn7BX+qv7F1zN9wtd8fdcw/c1HzKP+qf7D1zM27GvXCv +3Bs35+bcO/fBfXILbsF9cd/cDxcVYR4VwUfFXz2uFGGvUgRW5apctQjv1Lga +V+fqXJ1rcA2uyTW5JtfiWlyba3NtrsN1uC7X5bpcj+txfa7P9bkBN+CG3JAb +ciNuxMX+V8zH+thezCVcwiVcwiVcyqVcyqVcymW+J+MyLuMybqn4BacAMF4= + + "], $CellContext`fullList\[CapitalPhi]min1 = CompressedData[" +1:eJxTTMoPSmVkYGBIZYCBB/ajNHn0l+UrZrIAacm9b22P5F63n9XiultX+Yr9 +JaOgcL5jF+21mTaFXjp/3l5kUqV44M4z9uu2aHIpSZ2yr/l7/sjPq8ft82/M +Wvzv5lH7nMur3K5JHLFXbzjHPLfjkP0zjqI9+xUP2ndveus1ec9+e5mtYh1P +ivbau/Pmly2K221f+6CXNbhyp72nW0pR+97t9tGuuyWXqm2z379MX2/Phi32 +bIppNhkRm+2vtVcf3i61yf6O0EZmya8b7H1Om2ROPbXeXjopO33+pbX2FxWM +fp/8t9o+y7bm0SfvVfZL1k765bt9hf2Wtax3Iu2W2xvNYvW4dm+p/eub+/n3 +TFtib3Tlp4FZymJ754DwT7NdF9kvTmG+nG250N7s+BZtd7sF9uVaPSfuLJ1n +7xy80rf6zRz7Jy8Tfwf7zbbP33N9eu3JmfbK516rm8fMsN/y1pTpCct0+8Nl +M0M5Dk21dy34K+w2aYp9+GFBk6iSyfYxM/7bSaVPsr/8u1mhKHOifaaTh2FS +1QR7xxLhl78M++0PO11lvdDRa699zt5q9fdu+1e7HhwRqO2yd5nh4DNfvNM+ +P2TqId5j7fabmWWqNdvb7FdO0s25EtlqfyUybCazbYu9vELQwr36zfYmH71y +JA2b7IWKVL/bODTaR4W8toiJabBPfi7gUMBZb1/lnOdpXVhrDwDxbBLS + + "], $CellContext`fullList\[CapitalPhi]min2 = CompressedData[" +1:eJztxc9L03EcB+CVRZvzUE1IS61pi8QwjCBpyUstimVoKio0h2tkMHGLyCFa +S6vTqjkx0XKSzXLgipnNH2CadSg1JBvMGJVOUXPN7/fznoVipNipf6Bzz+WR +ai7lFG8QCATFgr+m8P9/O+K5c3R53YcyqWg6WjELV4oqVsvm0Pgt37MSNQ9D +1b0ScZgfV63SvqduPzx1BmN5xXfYBKM1N0IDmFVEFlSbApBl5Onf/wpguShs +0KtawDvjzNba3gV0KFobdEIOhfFr+sfZHApiDnUm1HMIHs9sGf/I4Ui5a1uW +kEde43BhgpzHfvHRSr+Wh09+ZeVHPY/Nltot7f08nkyYfIXTPDJcpe35IQyT +CpPoUyzDz4kzakkqQ/wxrztVybBvPH/PnTIGWRs0O+8yqOxV1u2tDOc3bkJ9 +D8OAJtjcP8Kwevtzqv0LQ8/IWLqWYzC7Tph3rTI8e/gy0h1KuNhdJ2uOIHQG +96rNMkJ0epqkK4lgsb0+GZNCqDvLfN5ThNOHcwcom3A9XD2mUxLmLhclqi4Q +UkKEiUOlBOnUb31XGYEfMLAD1wgHXyjDk28Rlhz+wa8mwpBgR660lhBu5PrF +DYSWtQ83bVaC07E0PP+IsGiRnJtpI3Qrspz3HYQH9qZFcQehurdSlOYiyD2Z +TfJeQnJJXNh6H6EjJ8pS84qQNFm1m39DqIjTmaLfEv4ALbErHQ== + + "], $CellContext`\[Chi]binodal[ + Pattern[$CellContext`\[CapitalPhi], + Blank[]], + Pattern[$CellContext`M, + Blank[]]] := (1/($CellContext`M (2 $CellContext`\[CapitalPhi] - 1))) + Log[$CellContext`\[CapitalPhi]/( + 1 - $CellContext`\[CapitalPhi])], $CellContext`\[Chi]hi = 5.}; + Typeset`initDone$$ = True), + SynchronousInitialization->True, + UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, + UnsavedVariables:>{Typeset`initDone$$}, + UntrackedVariables:>{Typeset`size$$}], "Manipulate", + Deployed->True, + StripOnInput->False], + Manipulate`InterpretManipulate[1]]], "Output", + CellChangeTimes->{ + 3.9062139578721237`*^9, 3.906214003058313*^9, 3.906214071516457*^9, { + 3.906214108229823*^9, 3.9062141307226133`*^9}, {3.9062141620246058`*^9, + 3.906214193906011*^9}, 3.906214262737403*^9, {3.90621434269233*^9, + 3.9062143610871553`*^9}, {3.906214392657584*^9, 3.906214403580596*^9}}, + CellLabel-> + "Out[344]=",ExpressionUUID->"3d3f30c9-a929-4460-9ca9-2ef5435a002f"] +}, Open ]] +}, +WindowSize->{1440, 782}, +WindowMargins->{{0, Automatic}, {Automatic, 0}}, +TaggingRules-><|"TryRealOnly" -> False|>, +FrontEndVersion->"13.0 for Mac OS X x86 (64-bit) (February 4, 2022)", +StyleDefinitions->"Default.nb", +ExpressionUUID->"804de27d-3276-4e9f-9c82-4cfab8c3be4c" +] +(* End of Notebook Content *) + +(* Internal cache information *) +(*CellTagsOutline +CellTagsIndex->{} +*) +(*CellTagsIndex +CellTagsIndex->{} +*) +(*NotebookFileOutline +Notebook[{ +Cell[CellGroupData[{ +Cell[580, 22, 35422, 835, 2216, "Input",ExpressionUUID->"cef6a996-4bad-4901-8d11-377a6ad0b7cd"], +Cell[36005, 859, 15874, 312, 375, "Output",ExpressionUUID->"3d3f30c9-a929-4460-9ca9-2ef5435a002f"] +}, Open ]] +} +] +*) + -- GitLab