(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 60082, 2041]*) (*NotebookOutlinePosition[ 60891, 2071]*) (* CellTagsIndexPosition[ 60847, 2067]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["The Potts model on the square lattice", "Title", TextAlignment->Center, TextJustification->0], Cell["\<\ V. Adamchik S. Finch R. Ziff\ \>", "Names", TextAlignment->Center, TextJustification->0, FontFamily->"Helvetica", FontSize->14], Cell[CellGroupData[{ Cell["Introduction", "Section"], Cell["\<\ The Potts model encompasses a number of problems in statistical \ physics and lattice theory. It generalizes the Ising model so that each spin \ can have more than two values. It includes the ice-vertex and bond \ percolation models as special cases. It also is related to graph-coloring \ problems: if we color each vertex in a lattice with one of q colors, in how \ many ways do we have exactly n pairs of adjacent vertices colored alike? We \ won't attempt to survey or even summarize the field, instead referring \ interested readers to [1].\ \>", "Text"], Cell["\<\ Here is one special case of the Potts model. We will define a \ certain polynomial (also called the Whitney-Tutte polynomial; see [2,3,4]) \ over the finite planar square lattice. The summation is taken over all \ subgraphs G of the lattice, the subgraphs being formed by removing bonds from \ the lattice but leaving all the sites. Note that many of these subgraphs \ will be disconnected, i.e., have multiple components. Consider the polynomial\ \ \>", "Text"], Cell[BoxData[ \(TraditionalForm\`W(z)\ = \[Sum]\+G z\^\(\ c(G) + n(G)\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`c(G)\)]], " is the number of connected clusters of bonds plus the number of isolated \ vertices, and ", Cell[BoxData[ \(TraditionalForm\`n(G)\)]], " is the nullity or cyclomatic number of ", Cell[BoxData[ \(TraditionalForm\`G\)]], ", i.e., the number of independent cycles in ", Cell[BoxData[ \(TraditionalForm\`G\)]], ". Euler's relation gives a simple way to compute ", Cell[BoxData[ \(TraditionalForm\`n(G)\)]], ":" }], "Text"], Cell[BoxData[ \(TraditionalForm\`n(G) = \(-v\) + c(G) + e(G)\)], "Text", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`v\)]], " is the total number of lattice vertices and ", Cell[BoxData[ \(TraditionalForm\`e(G)\)]], " is the number of bonds in ", Cell[BoxData[ \(TraditionalForm\`G\)]], ". For the infinite square lattice ", Cell[BoxData[ \(TraditionalForm\`\((v\ \[Rule] \ \[Infinity])\)\)]], ", it can be proved (see [2]) that" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(lim\+\(v \[Rule] \[Infinity]\)\[ThinSpace]\(log\ \(W(z)\)\)\/v\), "=", RowBox[{\(1\/\(4\ y\)\), " ", RowBox[{ SubsuperscriptBox["\[Integral]", AdjustmentBox[\(-\[Infinity]\), BoxMargins->{{0, 0}, {-0.395062, 0.395062}}, BoxBaselineShift->0.395062], "\[Infinity]"], \(\(\(sech(\(\[Pi]\ x\)\/\(2\ y\))\)\ \(log(\(cosh(x) - cos(2\ y)\)\/\(cosh(x) - 1\))\)\) \[DifferentialD]x\)}]}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`cos(y) = z\/2\)]], ". This remarkable integral formula is one of the few exact formulas \ associated with the Potts model, and it will be the focus of this paper. \ Define, for our purposes," }], "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(y) = \(1\/\(2\ y\)\) \(\[Integral]\_0 \%\[Infinity]\(\( sech(\(\[Pi]\ x\)\/\(2\ y\))\)\ \(log(\(cosh(x) - cos(2\ y)\)\/\(cosh(x) - 1\))\)\) \[DifferentialD]x\)\)], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["\<\ Observe that the integral can be rewritten in the following \ alternative form (after performing integration by parts):\ \>", "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(y) = \(\(4\ \(\(sin\^2\)(y)\)\)\/\[Pi]\) \(\[Integral]\_\(\ 0 \)\%\[Infinity]\(\(\(\( tan\^\(-1\)\)( tanh(\(\[Pi]\ x\)\/\(4\ y\)))\)\ \(coth(x\/2)\)\ \)\/\(cosh(x) - cos(2\ y)\)\) \[DifferentialD]x\)\)], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell[TextData[{ "Here is a graph of ", Cell[BoxData[ \(TraditionalForm\`\(P\_2\)(y)\)]], " on the interval ", Cell[BoxData[ FormBox[ RowBox[{"y", "\[Element]", RowBox[{"(", RowBox[{"0", ",", " ", FormBox[\(4\ \[Pi]\), "TraditionalForm"]}], ")"}]}], TraditionalForm]]] }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[ \(4\ Sin[y]\^2\)\/\[Pi]\ NIntegrate[ \(ArcTan[Tanh[\(\[Pi]\ x\)\/\(4\ y\)]]\ Coth[x\/2]\)\/\(Cosh[x] - Cos[2\ y]\), {x, 0, \[Infinity]}], {y, 0.0001, 4\ \[Pi]}]\)], "Input", CellLabel->"In[61]:="], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.0757881 0.0147151 0.375774 [ [.17539 .00222 -3 -9 ] [.17539 .00222 3 0 ] [.32696 .00222 -3 -9 ] [.32696 .00222 3 0 ] [.47854 .00222 -3 -9 ] [.47854 .00222 3 0 ] [.63011 .00222 -3 -9 ] [.63011 .00222 3 0 ] [.78169 .00222 -6 -9 ] [.78169 .00222 6 0 ] [.93327 .00222 -6 -9 ] [.93327 .00222 6 0 ] [.01131 .10866 -24 -4.5 ] [.01131 .10866 0 4.5 ] [.01131 .2026 -18 -4.5 ] [.01131 .2026 0 4.5 ] [.01131 .29655 -24 -4.5 ] [.01131 .29655 0 4.5 ] [.01131 .39049 -6 -4.5 ] [.01131 .39049 0 4.5 ] [.01131 .48443 -24 -4.5 ] [.01131 .48443 0 4.5 ] [.01131 .57838 -18 -4.5 ] [.01131 .57838 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .17539 .01472 m .17539 .02097 L s [(2)] .17539 .00222 0 1 Mshowa .32696 .01472 m .32696 .02097 L s [(4)] .32696 .00222 0 1 Mshowa .47854 .01472 m .47854 .02097 L s [(6)] .47854 .00222 0 1 Mshowa .63011 .01472 m .63011 .02097 L s [(8)] .63011 .00222 0 1 Mshowa .78169 .01472 m .78169 .02097 L s [(10)] .78169 .00222 0 1 Mshowa .93327 .01472 m .93327 .02097 L s [(12)] .93327 .00222 0 1 Mshowa .125 Mabswid .0617 .01472 m .0617 .01847 L s .0996 .01472 m .0996 .01847 L s .13749 .01472 m .13749 .01847 L s .21328 .01472 m .21328 .01847 L s .25117 .01472 m .25117 .01847 L s .28907 .01472 m .28907 .01847 L s .36486 .01472 m .36486 .01847 L s .40275 .01472 m .40275 .01847 L s .44064 .01472 m .44064 .01847 L s .51643 .01472 m .51643 .01847 L s .55433 .01472 m .55433 .01847 L s .59222 .01472 m .59222 .01847 L s .66801 .01472 m .66801 .01847 L s .7059 .01472 m .7059 .01847 L s .7438 .01472 m .7438 .01847 L s .81958 .01472 m .81958 .01847 L s .85748 .01472 m .85748 .01847 L s .89537 .01472 m .89537 .01847 L s .97116 .01472 m .97116 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .10866 m .03006 .10866 L s [(0.25)] .01131 .10866 1 0 Mshowa .02381 .2026 m .03006 .2026 L s [(0.5)] .01131 .2026 1 0 Mshowa .02381 .29655 m .03006 .29655 L s [(0.75)] .01131 .29655 1 0 Mshowa .02381 .39049 m .03006 .39049 L s [(1)] .01131 .39049 1 0 Mshowa .02381 .48443 m .03006 .48443 L s [(1.25)] .01131 .48443 1 0 Mshowa .02381 .57838 m .03006 .57838 L s [(1.5)] .01131 .57838 1 0 Mshowa .125 Mabswid .02381 .0335 m .02756 .0335 L s .02381 .05229 m .02756 .05229 L s .02381 .07108 m .02756 .07108 L s .02381 .08987 m .02756 .08987 L s .02381 .12745 m .02756 .12745 L s .02381 .14624 m .02756 .14624 L s .02381 .16502 m .02756 .16502 L s .02381 .18381 m .02756 .18381 L s .02381 .22139 m .02756 .22139 L s .02381 .24018 m .02756 .24018 L s .02381 .25897 m .02756 .25897 L s .02381 .27776 m .02756 .27776 L s .02381 .31533 m .02756 .31533 L s .02381 .33412 m .02756 .33412 L s .02381 .35291 m .02756 .35291 L s .02381 .3717 m .02756 .3717 L s .02381 .40928 m .02756 .40928 L s .02381 .42807 m .02756 .42807 L s .02381 .44686 m .02756 .44686 L s .02381 .46564 m .02756 .46564 L s .02381 .50322 m .02756 .50322 L s .02381 .52201 m .02756 .52201 L s .02381 .5408 m .02756 .5408 L s .02381 .55959 m .02756 .55959 L s .02381 .59716 m .02756 .59716 L s .02381 .61595 m .02756 .61595 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02382 .60332 m .02499 .6033 L .02606 .60326 L .0273 .60319 L .02846 .60308 L .03054 .60283 L .03279 .60244 L .03527 .60189 L .03791 .60115 L .04262 .59946 L .0475 .59721 L .05206 .59464 L .06245 .5871 L .07305 .57705 L .08275 .56577 L .10459 .53325 L .1443 .44937 L .18249 .3392 L .22313 .18986 L .24193 .10913 L .25182 .06354 L .25722 .03769 L .25986 .02482 L .26101 .01915 L .26226 .01648 L .26358 .02287 L .26483 .02876 L .26758 .0412 L .27341 .0653 L .28367 .10108 L .29326 .1281 L .30384 .15191 L .31404 .16976 L .32362 .18257 L .33373 .19242 L .33941 .19646 L .34463 .19929 L .34937 .20115 L .35202 .20192 L .35443 .20245 L .35552 .20264 L .35668 .2028 L .35795 .20294 L .35912 .20303 L .36021 .20309 L .36123 .20311 L .36235 .2031 L .36355 .20306 L .36422 .20302 L .36493 .20297 L Mistroke .3662 .20284 L .36772 .20263 L .3691 .2024 L .37181 .20181 L .37437 .20109 L .37912 .19936 L .38425 .19692 L .39378 .1909 L .40423 .18216 L .4253 .15826 L .44485 .12914 L .46578 .0911 L .48563 .04887 L .49107 .03629 L .49405 .02922 L .49684 .02248 L .49799 .01966 L .49921 .01667 L .49985 .01509 L .50054 .01606 L .50176 .01903 L .50719 .03167 L .52811 .07181 L .54709 .09799 L .5572 .10847 L .56674 .11633 L .57681 .12265 L .58248 .12536 L .5877 .12733 L .59239 .12869 L .59498 .12927 L .59742 .12972 L .5997 .13005 L .6009 .13018 L .60217 .1303 L .60324 .13038 L .60442 .13045 L .6055 .13049 L .6065 .13051 L .60771 .13052 L .60903 .13049 L .60969 .13047 L .6104 .13044 L .61169 .13036 L .61297 .13026 L .61418 .13014 L .61644 .12986 L .61911 .12942 L .62158 .12892 L .62716 .12746 L Mistroke .63182 .12588 L .63684 .12384 L .6459 .11926 L .66653 .10473 L .68768 .08427 L .70777 .05993 L .7181 .04565 L .72782 .03114 L .73207 .02449 L .73438 .02078 L .73655 .01726 L .73768 .0154 L .73889 .01603 L .73958 .01715 L .74024 .01821 L .74146 .02017 L .74599 .02715 L .76567 .05321 L .77585 .06408 L .78665 .07382 L .79612 .08091 L .80658 .08724 L .81608 .09169 L .825 .09478 L .8297 .09599 L .83484 .09701 L .83744 .0974 L .84021 .09772 L .84155 .09785 L .84282 .09795 L .84395 .09802 L .8452 .09808 L .84636 .09812 L .84744 .09814 L .84859 .09815 L .84982 .09814 L .85096 .09812 L .85205 .09809 L .85302 .09804 L .85408 .09798 L .8565 .0978 L .85908 .09753 L .86371 .09686 L .86909 .09579 L .87403 .09452 L .88517 .09072 L .895 .08629 L .90536 .08057 L .9434 .05078 L .97619 .01472 L Mistroke Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", CellLabel->"From In[61]:=", ImageSize->{288, 177.9375}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgo`03?00bo`03?00_o`04?00_o`04?00_o`Lo00ol00c`04?l0 01co00<<04_o00_l00c`04Ol001co00_l00c`05Ol00c`0_l00c`06_l00c`09ol00c`0=Ol902oo0005l`03?006o`03?006o`03?006ldKo00o`03?00do`03 ?00Vo`03?00Fo`03?01fo`007?l00c`0A_l00c`03ol00c`0<_l00c`0:?l00``04ol00``0N?l001co 00?l000[o00"], ImageRangeCache->{{{0, 359}, {221.375, 0}} -> {-1.09682, -0.111942, 0.048533, 0.009788}}], Cell[BoxData[ FormBox[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False], TraditionalForm]], "Output", CellLabel->"Out[61]="] }, Open ]], Cell[TextData[{ "Although it isn't known whether the function ", Cell[BoxData[ \(TraditionalForm\`\(P\_2\)(y)\)]], " has a closed-form expression for all values of ", Cell[BoxData[ \(TraditionalForm\`y\)]], ", it can be evaluated explicitly for any ", Cell[BoxData[ \(TraditionalForm\`y\)]], " which is a rational multiple of \[Pi]. We consider several particular \ cases of ", Cell[BoxData[ \(TraditionalForm\`\(P\_2\)(y)\)]], " in the next section. In addition to the integrals (1) and (2), we'll need \ another representation for ", Cell[BoxData[ FormBox[ FormBox[\(\(P\_2\)(y)\), "TraditionalForm"], TraditionalForm]]], " which can be obtained by performing the logarithmic substitution in (1):" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(y)\), "TraditionalForm"], "=", \(\(2 q\)\/\(\[Pi]\ \)\ \(\[Integral]\_0\%1 x\^\(q - 1\)\/\(x\^\(2\ q\) + 1\)\ \(log(\(x\^\(4 p\) - 2\ x\^\(2 p\)\ \(cos(\(2\ p\ \[Pi]\)\/q)\) + 1 \)\/\((1 - x\^\(2 p\))\)\^2)\) \[DifferentialD]x\)\)}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell[TextData[{ "where ", Cell[BoxData[ \(TraditionalForm\`y = \(\[Pi]\ p\)\/q\)]], ", and p and ", "q", " are positive integers." }], "Text"] }, Closed]], Cell[CellGroupData[{ Cell["Particular Cases", "Section"], Cell["\<\ The following cases had been evaluated previously (see [2]), though \ no derivation of them had been given.\ \>", "Text"], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \ 0\)]]], "Subsubsection"], Cell[TextData[{ "We shall start with the representation (1). Changing the variable of \ integration ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] \(2\ y\ z\)\/\[Pi]\)]], ", we obtain:" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(y)\), "TraditionalForm"], "=", \(\(1\/\[Pi]\) \(\[Integral]\_0\%\[Infinity]\( sech(z)\)\ \(log(\(cosh(\(2\ y\ z\)\/\[Pi]) - cos(2\ y)\)\/\(cosh(\(2\ y\ z\)\/\[Pi]) - 1\))\) \[DifferentialD]z\)\)}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell["Expanding the integrand at y = 0, we get:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(0)\), "TraditionalForm"], "=", \(\(1\/\[Pi]\) \(\[Integral]\_0 \%\[Infinity]\(\( sech(z)\)\ \(log(1 + \[Pi]\^2\/z\^2)\)\) \[DifferentialD]z\)\)}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell["or", "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(0)\), "TraditionalForm"], "=", \(\(\(2\ \)\/\[Pi]\) \(\[Integral]\_0\%1 \(\( log(\(log\^2\)(x) + \[Pi]\^2)\)\/\(x\^2 + 1\)\) \[DifferentialD]x\) - \n\t\t \(\(4\ \)\/\[Pi]\) \(\[Integral]\_0\%1\(\( log(log(1\/x))\)\/\(x\^2 + 1\)\) \[DifferentialD]x\)\)}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "Both integrals can be found in ", "Gradshteyn and Ryzhik, ", StyleBox["Table of Integrals, Series and Products", FontSlant->"Italic"], ":" }], "Text"], Cell[BoxData[ \(TraditionalForm \`\[Integral]\_0\%1\(\( log(log(1\/x))\)\/\(x\^2 + 1\)\) \[DifferentialD]x = \(\ \[Pi]\)\/2\ \(log(\(\(\[CapitalGamma](3\/4)\ \)\/\(\[CapitalGamma](1\/4)\)\) \@\(2\ \[Pi]\))\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell[BoxData[ \(TraditionalForm \`\[Integral]\_0\%1\(\( log(\(log\^2\)(x) + \[Pi]\^2)\)\/\(x\^2 + 1\)\) \[DifferentialD]x = \[Pi]\ \(log( \(\@\(\[Pi]\/2\)\ \(\[CapitalGamma](1\/4)\)\)\/\(2\ \(\[CapitalGamma](3\/4)\)\))\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell["Therefore,", "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(0)\), "TraditionalForm"], "=", \(4\ \(log( \(\[CapitalGamma](1\/4)\)\/\(2 \( \[CapitalGamma](3\/4)\)\)\ ) \)\)}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \[Pi]\/4\)]]], "Subsubsection"], Cell["From representation (3) with p = 1, q = 4 we have:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ FormBox[\(\(P\_2\)(\[Pi]\/4)\), "TraditionalForm"], "=", \(8\/\(\[Pi]\ \)\ \(\[Integral]\_0\%1 x\^3\/\(1 + x\^8\)\ \(log(\(1 + x\^4\)\/\((1 - x\^2)\)\^2)\) \[DifferentialD]x\) = \(4\/\(\[Pi]\ \)\ \(\[Integral]\_0\%1 x\/\(1 + x\^4\)\ \(log(\(1 - x\^4\)\/\(\((1 - x\^2)\) \((1 - x)\)\^2\))\) \[DifferentialD]x\) = \n\t\t 2\/\(\[Pi]\ \)\ \(\[Integral]\_0\%1\( log(1 + x)\)\/\(1 + x\^2\)\ \[DifferentialD]x\)\t - 8\/\(\[Pi]\ \)\ \(\[Integral]\_0\%1 x\/\(1 + x\^4\)\ \(log(1 - x)\) \[DifferentialD]x\)\)\)}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["\<\ We evaluate each of the integrals separately. The first integral \ is:\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(2\/\[Pi]\ \(\[Integral]\_0\%1\( Log[1 + x]\/\(x\^2 + 1\)\) \[DifferentialD]x\)\)], "Input", CellLabel->"In[5]:="], Cell[BoxData[ \(TraditionalForm \`\(\(-\[ImaginaryI]\)\ \(\(log\^2\)(\(-1\) - \[ImaginaryI])\) + \[ImaginaryI]\ \(\(log\^2\)(\(-1\) + \[ImaginaryI])\) + 2\ \[Pi]\ \(log(2)\)\)\/\(2\ \[Pi]\)\)], "Output", CellLabel->"Out[5]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(TraditionalForm\`FullSimplify[%]\)], "Input", CellLabel->"In[6]:=", FormatType->StandardForm], Cell[BoxData[ \(TraditionalForm\`\(log(2)\)\/4\)], "Output", CellLabel->"Out[6]="] }, Open ]], Cell[TextData[{ "Changing the variable of integration ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] 1\/z\)]], " in the second integral, we observe that:" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\[Integral]\_0\%1 x\/\(1 + x\^4\)\ \(log(1 - x)\) \[DifferentialD]x\), "=", RowBox[{ RowBox[{"-", FractionBox[ TagBox["C", Catalan], "4"]}], "+", \(\[Integral]\_1\%\[Infinity] z\/\(1 + z\^4\)\ \(log(z - 1)\) \[DifferentialD]z\)}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["\<\ where C is Catalan's constant. Then taking into account that:\ \>", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\[Integral]\_0\%\[Infinity] x\/\(1 + x\^4\)\ \(log(1 - x)\) \[DifferentialD]x\), "=", RowBox[{\(-\(1\/16\)\), " ", RowBox[{"(", RowBox[{ RowBox[{"4", " ", TagBox["C", Catalan]}], "-", \(2\ \[ImaginaryI]\ \[Pi]\^2\), "-", \(\[Pi]\ \(log(2)\)\), "+", \(4\ \[Pi]\ \(log(1 + \@2)\)\)}], ")"}]}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["we immediately arrive at:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\[Integral]\_0\%1 x\/\(1 + x\^4\)\ \(log(1 - x)\) \[DifferentialD]x\), "=", RowBox[{\(-\(1\/32\)\), " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", TagBox["C", Catalan]}], "-", \(\[Pi]\ \(log(2)\)\), "+", \(4\ \[Pi]\ \(log(1 + \@2)\)\)}], ")"}]}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0], Cell["Finally,", "Text"], Cell[BoxData[ FormBox[ RowBox[{\(\(P\_2\)(\[Pi]\/4)\), "=", RowBox[{ FractionBox[ RowBox[{"2", " ", TagBox["C", Catalan]}], "\[Pi]"], "+", \(log(1 + \@2)\)}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \[Pi]\/3\)]]], "Subsubsection"], Cell["p = 1, q = 3. From (2), we have:", "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(\[Pi]\/3) = \(\(\ 6\)\/\[Pi]\) \(\[Integral]\_0\%1\(\( x\^2\ \)\/\(x\^6 + 1\)\) \(log(\(x\^4 + x\^2 + 1\)\/\((1 - x\^2)\)\^2)\) \[DifferentialD]x\)\)], "Equation", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "Since ", Cell[BoxData[ \(TraditionalForm\`x\^4 + x\^2 + 1 = \(1 - x\^6\)\/\(1 - x\^2\)\)]], ", we have after further simplifications:" }], "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(\[Pi]\/3) = \(\(2\ \)\/\[Pi]\) \(\[Integral]\_0\%1\(\( log(1 - x\^2)\)\/\(x\^2 + 1\)\) \[DifferentialD]x\) - 18\/\[Pi]\ \(\[Integral]\_0\%1\(\( x\^2\ \(log(1 - x\^2)\)\)\/\(x\^6 + 1\)\) \[DifferentialD]x\)\)], "Equation", TextAlignment->Center, TextJustification->0], Cell["The first integral is:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[ 2\/\[Pi]\ \(\[Integral]\_0\%1\( Log[1 - x\^2]\/\(1 + x\^2\)\) \[DifferentialD]x\)]\)], "Input", CellLabel->"In[9]:="], Cell[BoxData[ \(TraditionalForm \`\(\[Pi]\ \(log(2)\) - 4\ \[ImaginaryI]\ \(\(Li\_2\)(1\/2 - \[ImaginaryI]\/2)\) + 4\ \[ImaginaryI]\ \(\(Li\_2\)(1\/2 + \[ImaginaryI]\/2)\)\)\/\(4\ \[Pi]\)\)], "Output", CellLabel->"Out[9]="] }, Open ]], Cell["Taking into account that:", "Text"], Cell[BoxData[{ FormBox[ RowBox[{\(\(Li\_2\) \((1\/2 - \[ImaginaryI]\/2)\)\), "=", RowBox[{\(\[Pi]\^2\/48\), "-", RowBox[{"\[ImaginaryI]", " ", TagBox["C", Catalan]}], "-", \(1\/2\ \(\(log\^2\)(1\/2 + \[ImaginaryI]\/2)\)\)}]}], TraditionalForm], FormBox[ RowBox[{\(\(Li\_2\) \((1\/2 + \[ImaginaryI]\/2)\)\), "=", RowBox[{\(\[Pi]\^2\/48\), "+", RowBox[{"\[ImaginaryI]", TagBox["C", Catalan]}], " ", "-", \(1\/2\ \(\(log\^2\)(1\/2 - \[ImaginaryI]\/2)\)\)}]}], TraditionalForm]}], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["we arrive at:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\(\(2\ \)\/\[Pi]\) \(\[Integral]\_0\%1\(\( log(1 - x\^2)\)\/\(1 + x\^2\)\) \[DifferentialD]x\)\), "=", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"2", " ", TagBox["C", Catalan]}], "\[Pi]"]}], "+", \(\(log(2)\)\/2\)}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["In the second integral", "Text"], Cell[BoxData[ \(TraditionalForm \`18\/\[Pi]\ \(\[Integral]\_0\%1\(\( x\^2\ \(log(1 - x\^2)\)\)\/\(1 + x\^6\)\) \[DifferentialD]x\)\)], "Equation", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "we change the variable ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] 1\/z\)]], ":" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\[Integral]\_0\%1\(\( x\^2\ \(log(1 - x\^2)\)\)\/\(1 + x\^6\)\) \[DifferentialD]x\), " ", "=", " ", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"2", " ", TagBox["C", Catalan]}], "9"]}], "+", \(\[Integral]\_1 \%\[Infinity]\(\( z\^2\ \(log(z\^2 - 1)\)\)\/\(1 + z\^6\)\) \[DifferentialD]z\)}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0], Cell["Taking into account that:", "Text"], Cell[BoxData[ \(TraditionalForm \`\[Integral]\_0 \%\[Infinity]\(\( x\^2\ \(log(1 - x\^2)\)\)\/\(x\^6 + 1\)\) \[DifferentialD]x\ = \(\[ImaginaryI]\ \[Pi]\^2\)\/12 - \[Pi]\/6\ \ \(log(2)\)\)], "Equation",\ TextAlignment->Center, TextJustification->0], Cell[" it then follows that:", "Text"], Cell[BoxData[ FormBox[ RowBox[{" ", RowBox[{ \(\[Integral]\_0\%1\(\( x\^2\ \(log(1 - x\^2)\)\)\/\(x\^6 + 1\)\) \[DifferentialD]x\), "=", RowBox[{ RowBox[{"-", FractionBox[ TagBox["C", Catalan], "9"]}], "-", \(\[Pi]\/12\ \(log(2)\)\)}]}]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0], Cell["Finally,", "Text"], Cell[BoxData[ \(TraditionalForm\`\(P\_2\)(\[Pi]\/3) = \ 2\ \(log(2)\)\)], "Equation", TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \[Pi]\/2\)]]], "Subsubsection"], Cell[TextData[{ "If ", Cell[BoxData[ \(TraditionalForm\`y = \[Pi]\/2\)]], ", then ", Cell[BoxData[ \(TraditionalForm\`p = 1\)]], " and ", Cell[BoxData[ \(TraditionalForm\`q\ = \ 2\)]], ". From (2) we have:" }], "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(\[Pi]\/2) = \(\(\ 4\)\/\[Pi]\) \(\[Integral]\_0\%\[Infinity]\(\( tan\^\(-1\)\)(tanh(x\/2))\)\ \(csch(x)\) \[DifferentialD]x\)\)], "Equation", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "Changing the variable ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] 2\ \(\(tanh\^\(-1\)\)(z)\)\)]], ", the integral is simplified to:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(4\/\[Pi]\ \(\[Integral]\_0\%1\( ArcTan[z]\/z\) \[DifferentialD]z\)\)], "Input", CellLabel->"In[10]:="], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], "\[Pi]"], TraditionalForm]], "Output", CellLabel->"Out[10]="] }, Open ]], Cell["Thus,", "Text"], Cell[BoxData[ FormBox[ RowBox[{\(\(P\_2\)(\[Pi]\/2)\), "=", FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], "\[Pi]"]}], TraditionalForm]], "Equation", TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \(2\ \[Pi]\)\/3\)]]], "Subsubsection"], Cell["p = 2, q = 3. We shall start with representation (1):", "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(\(2 \[Pi]\)\/3) = \(\(6\/\[Pi]\) \(\[Integral]\_0\%1 x\^2\/\(x\^6 + 1\)\ \(log(\(x\^8 + x\^4 + 1\)\/\((1 - x\^4)\)\^2)\) \[DifferentialD]x\) = \n\n\t\t\t \(2\/\[Pi]\) \(\[Integral]\_0\%1\( log(1 - x\^4)\)\/\(x\^2 + 1\)\ \[DifferentialD]x\) - \(18\/\[Pi]\) \(\[Integral]\_0\%1 x\^2\/\(x\^6 + 1\)\ \(log(1 - x\^4)\) \[DifferentialD]x\)\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[ 2\/\[Pi]\ \(\[Integral]\_0\%1\( Log[1 - x\^4]\/\(x\^2 + 1\)\) \[DifferentialD]x\)]\)], "Input", CellLabel->"In[23]:="], Cell[BoxData[ \(TraditionalForm \`\(\[Pi]\ \(log(2)\) - 2\ \[ImaginaryI]\ \(\(Li\_2\)(1\/2 - \[ImaginaryI]\/2)\) + 2\ \[ImaginaryI]\ \(\(Li\_2\)(1\/2 + \[ImaginaryI]\/2)\)\)\/\[Pi]\)], "Output", CellLabel->"Out[23]="] }, Open ]], Cell["Taking into account both formulas in (4), we obtain:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\(2\/\[Pi]\) \(\[Integral]\_0\%1\( log(1 - x\^4)\)\/\(x\^2 + 1\)\ \[DifferentialD]x\)\), "=", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], "\[Pi]"]}], "+", \(\(3\ \(log(2)\)\)\/2\)}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["The second integral is:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\(18\/\[Pi]\) \(\[Integral]\_0\%1 x\^2\/\(x\^6 + 1\)\ \(log(1 - x\^4)\) \[DifferentialD]x\)\), "=", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], "\[Pi]"]}], " ", "+", \(3\/4\ \(log(81\/64)\)\)}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell["since", "Text"], Cell[BoxData[ \(TraditionalForm\`\(\ \[Integral]\_0\%\[Infinity]\(\( x\^2\ \(log(1 - x\^4)\)\)\/\(x\^6 + 1\)\) \[DifferentialD]x = \(\[ImaginaryI]\ \[Pi]\^2\)\/12 + \[Pi]\/12\ \ \(log(81\/64)\)\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell["and", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\[Integral]\_0\%1 x\^2\/\(x\^6 + 1\)\ \(log(1 - x\^4)\) \[DifferentialD]x\), "=", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], "9"]}], "+", \(\[Integral]\_1\%\[Infinity] x\^2\/\(x\^6 + 1\)\ \(log(x\^4 - 1)\) \[DifferentialD]x\)}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell["Therefore,", "Text"], Cell[BoxData[ \(TraditionalForm\`\(P\_2\)(\(2 \[Pi]\)\/3) = 3 \( log(4\/3)\)\)], "Text",\ TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \[Pi]\)]]], "Subsubsection"], Cell["It immediately follows from (2) that", "Text"], Cell[BoxData[ FormBox[ RowBox[{\(\(P\_2\)(k\ \[Pi]) = 0\), ",", " ", FormBox[\(k = 1\), "TraditionalForm"], ",", FormBox["2", "TraditionalForm"], ",", "..."}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0] }, Closed]], Cell[CellGroupData[{ Cell[TextData[Cell[BoxData[ \(TraditionalForm\`y = \(3 \[Pi]\)\/2\)]]], "Subsubsection"], Cell["From representation (2) we have", "Text"], Cell[BoxData[ \(TraditionalForm \`\(P\_2\)(\(3 \[Pi]\)\/2) = \(\(\ 4\)\/\[Pi]\) \((\[Integral]\_0\%1\(\(\( tan\^\(-1\)\)(z)\)\/z\) \[DifferentialD]z + 2 \(\[Integral]\_0\%1 \(\( z\ \(\(tan\^\(-1\)\)(z)\)\)\/\(z\^2 + 3\)\) \[DifferentialD]z\) - 6 \(\[Integral]\_0\%1 \(\( z\ \(\(tan\^\(-1\)\)(z)\)\)\/\(3\ z\^2 + 1\)\) \[DifferentialD]z\))\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell["All these integrals can be easily evaluated:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[ 4\/\[Pi]\ \(( \[Integral]\_0\%1\( ArcTan[z]\/z\) \[DifferentialD]z + 2\ \(\[Integral]\_0\%1\(\( z\ ArcTan[z]\)\/\(z\^2 + 3\)\) \[DifferentialD]z\) - 6\ \(\[Integral]\_0\%1\(\( z\ ArcTan[z]\)\/\(3\ z\^2 + 1\)\) \[DifferentialD]z\))\) // TrigToExp]\)], "Input", CellLabel->"In[24]:="], Cell[BoxData[ FormBox[ RowBox[{\(1\/\[Pi]\), RowBox[{"(", RowBox[{ RowBox[{"4", " ", TagBox["C", Catalan]}], "-", \(4\ \[ImaginaryI]\ \[Pi]\^2\), "-", \(\[Pi]\ \(log(7 - 4\ \@3)\)\), "+", \(2\ \[Pi]\ \(log(2 + \@3)\)\), "-", \(4\ \[ImaginaryI]\ \(\(Li\_2\)(\(-\[ImaginaryI]\)\ \@\(7 - 4\ \@3\))\)\), "+", \(4\ \[ImaginaryI]\ \(\(Li\_2\)(\[ImaginaryI]\ \@\(7 - 4\ \@3\))\)\)}], ")"}]}], TraditionalForm]], "Output", CellLabel->"Out[24]="] }, Open ]], Cell["Taking into account that:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\(Li\_2\)(\[ImaginaryI]\ \@\(7 - 4\ \@3\)) - \(Li\_2\)(\(-\[ImaginaryI]\)\ \@\(7 - 4\ \@3\))\), "=", RowBox[{\(1\/6\), " ", "\[ImaginaryI]", " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", TagBox["C", Catalan]}], "-", \(\[Pi]\ \(log(2 + \@3)\)\)}], ")"}]}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0], Cell["we obtain:", "Text"], Cell[BoxData[ FormBox[ RowBox[{\(\(P\_2\)(\(3 \[Pi]\)\/2)\), "=", RowBox[{\(2\/3\ \(log(2 + \@3)\)\), "-", FractionBox[ RowBox[{"4", " ", TagBox["C", Catalan]}], \(3\ \[Pi]\)]}]}], TraditionalForm]], "Text", TextAlignment->Center, TextJustification->0] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["First Derivative", "Section"], Cell[TextData[{ "The Potts model includes the bond percolation model ([5]) as a special \ case. In [2], Temperley and Lieb derived a closed-form expression for the \ limiting mean cluster density at the critical probability ", Cell[BoxData[ \(TraditionalForm\`1\/2\)]], ": " }], "Text"], Cell[BoxData[ FormBox[ RowBox[{\(\(K\_b\)(1\/2)\), " ", "=", RowBox[{ RowBox[{\(-\(1\/16\)\), "-", RowBox[{\(\(csc(y)\)\/4\), " ", FractionBox[\(\[PartialD]\(W(y)\)\), \(\[PartialD]y\), MultilineFunction->None]}]}], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}]}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["\<\ We will show here that formula (5) can be further simplified \ to\ \>", "Text"], Cell[BoxData[ \(TraditionalForm\`\(K\_b\)(1\/2)\ = \(3\ \@3\)\/2 - 41\/16\)], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["As far as we know, the exact expression in (6) is new.", "Text"], Cell[TextData[{ StyleBox["Proof.", FontWeight->"Bold", FontSlant->"Italic"], " We shall proceed with representation (2). Differentiating, we obtain" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\(\(P\_2\)(y)\)\), \(\[PartialD]y\), MultilineFunction->None], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(\(\(8\ \@3\ \)\/\(\[Pi]\ \)\) \(\[Integral]\_0 \%\[Infinity]\(\(\(\( tan\^\(-1\)\)(tanh(\(3\ x\)\/4))\)\ \(sinh(x)\)\)\/\((2\ \(cosh(x)\) + 1)\)\^2\) \[DifferentialD]x\) - \(27\/\[Pi]\^2\) \(\[Integral]\_0\%\[Infinity] x\ \(cosh(x\/2)\)\ \(csch(3\ x)\) \[DifferentialD]x\)\)}], TraditionalForm]], "NumberedEquation",\ FormatType->StandardForm], Cell[TextData[{ "The second integral in (7) is almost trivial and can be easily evaluated \ by ", StyleBox["Mathematica", FontSlant->"Italic"], " " }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(27\/\[Pi]\^2\ \(\[Integral]\_0\%\[Infinity]\( x\ Cosh[x\/2]\ Csch[3\ x]\) \[DifferentialD]x\)\)], "Input", CellLabel->"In[8]:="], Cell[BoxData[ FormBox[ FractionBox[ RowBox[{"3", " ", RowBox[{"(", RowBox[{ TagBox[ RowBox[{"\[Zeta]", "(", RowBox[{ TagBox["2", (Editable -> True)], ",", TagBox[\(5\/12\), (Editable -> True)]}], ")"}], InterpretTemplate[ Zeta[ #, #2]&]], "+", TagBox[ RowBox[{"\[Zeta]", "(", RowBox[{ TagBox["2", (Editable -> True)], ",", TagBox[\(7\/12\), (Editable -> True)]}], ")"}], InterpretTemplate[ Zeta[ #, #2]&]]}], ")"}]}], \(4\ \[Pi]\^2\)], TraditionalForm]], "Output", CellLabel->"Out[8]="] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(FullSimplify[%]\)], "Input", CellLabel->"In[9]:="], Cell[BoxData[ \(TraditionalForm\`6 - 3\ \@3\)], "Output", CellLabel->"Out[9]="] }, Open ]], Cell[TextData[{ "In the first integral in (7), we change the variable of integration ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] \(-2\)\ \(log(z)\)\)]], ". We obtain:" }], "Text"], Cell[BoxData[ \(TraditionalForm \`\[Integral]\_0 \%\[Infinity]\(\(\(\( tan\^\(-1\)\)(tanh(\(3\ x\)\/4))\)\ \(sinh(x)\)\)\/\((2\ \(cosh(x)\) + 1)\)\^2\) \[DifferentialD]x = \[Integral]\_0\%1 \(\( z\ \((z\^4 - 1)\)\ \(\(tan\^\(-1\)\)(\(z\^3 - 1\)\/\(z\^3 + 1\))\)\)\/\(( z\^4 + z\^2 + 1)\)\^2\) \[DifferentialD]\)], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell[TextData[{ "Then using ", StyleBox["Mathematica", FontSlant->"Italic"], ", we get:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(8\ \@3\)\/\[Pi]\ \(\[Integral]\_0\%1 \(\( z\ \((z\^4 - 1)\)\ ArcTan[\(z\^3 - 1\)\/\(z\^3 + 1\)]\)\/\((z\^4 + z\^2 + 1)\)\^2 \) \[DifferentialD]z\)\)], "Input", CellLabel->"In[10]:="], Cell[BoxData[ \(TraditionalForm\`\(-\@3\)\ \((\(-2\) + \@3)\)\)], "Output", CellLabel->"Out[10]="] }, Open ]], Cell["Thus we have found that:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\(\(P\_2\)(y)\)\), \(\[PartialD]y\), MultilineFunction->None], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(5\ \@3 - 9\)}], " "}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["which completes the proof of identity (5).", "Text"] }, Closed]], Cell[CellGroupData[{ Cell["Second Derivative", "Section"], Cell["\<\ Just as the first derivative has an interpretation as a mean value, \ the second derivative can be regarded as a measure of \"fluctuations\" in \ the number of clusters corresponding to the bond percolation model.\ \>", "Text"], Cell["We will prove here that ", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(\((\(cot(y)\)\ d\/dy)\)\^2\ \(\(P\_2\)(y)\)\), SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(\(-\(1\/2\)\)\ \((25 - 8\ \@3)\) + 18\/\[Pi]\)}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell["which, again, we believe is a new result.", "Text"], Cell[TextData[{ StyleBox["Proof", FontWeight->"Bold", FontSlant->"Italic"], ". First of all, we observe that" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{\(\((\(cot(y)\)\ d\/dy)\)\^2\ \(\(P\_2\)(y)\)\), SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", RowBox[{ RowBox[{\(-\(20\/3\)\), "+", \(4\ \@3\), "+", " ", RowBox[{\(1\/3\), " ", FractionBox[\(\[PartialD]\^2\(\( P\_2\)(y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None]}]}], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}]}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "taking into account identity (5). Differentiating the integral on the \ right side of formula (2) twice with respect to ", Cell[BoxData[ \(TraditionalForm\`y\)]], ", we obtain:" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\^2\(\( P\_2\)(y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(\(-\(16\/\[Pi]\)\) \(\[Integral]\_0 \%\[Infinity]\(\(\ \(\(tan\^\(-1\)\)(tanh(\(3\ x\)\/4))\)\ \((7\ \(sinh(x)\) + sinh(2\ x))\)\)\/\(\ \((2\ \(cosh(x)\) + 1)\)\^3\)\) \[DifferentialD]x\) - \(\(72\ \@3\)\/\(\[Pi]\^2\ \)\) \(\[Integral]\_0 \%\[Infinity]\(\(\ x\ \(sinh(x\/2)\)\)\/\(\((2\ \(cosh(x)\) - 1)\)\ \((2\ \(cosh(x)\) + 1)\)\^2\)\) \[DifferentialD]x\) + \(81\/\[Pi]\^3\) \(\[Integral]\_0 \%\[Infinity]\(\(\ x\ \(csch(x\/2)\)\)\/\(\ 2\ \(cosh(2\ x)\) + 1\)\) \[DifferentialD]x\) - \(243\/\(4\ \[Pi]\^3\)\) \(\[Integral]\_0 \%\[Infinity]\(\(\ x\^2\ \(sech(x\/2)\)\)\/\(\ \((1 - 2\ \(cosh(x)\))\)\^2\)\) \[DifferentialD]x\)\)}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0], Cell[TextData[{ "Changing the variable of integration ", Cell[BoxData[ \(TraditionalForm\`x \[Rule] \(-2\)\ \(log(z)\)\)]], ", we arrive at:" }], "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\^2\(\( P\_2\)(y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(\(-\(16\/\[Pi]\)\) \(\[Integral]\_0\%1 \(\( z\ \((z\^8 + 7\ z\^6 - 7\ z\^2 - 1)\)\ \(\(tan\^\(-1\)\)(\(z\^3 - 1\)\/\(z\^3 + 1\))\)\)\/\(( z\^4 + z\^2 + 1)\)\^3\) \[DifferentialD]z\) + \(\(144\ \@3\)\/\(\[Pi]\^2\ \)\) \(\[Integral]\_0\%1 \(\( z\^4\ \((1 - z\^2)\)\ \(log(z)\)\)\/\(\((z\^4 - z\^2 + 1)\)\ \((z\^4 + z\^2 + 1)\)\^2\)\) \[DifferentialD]z\) + \(648\/\[Pi]\^3\) \(\[Integral]\_0\%1 \(\(\ z\^4\ \(log(z)\)\)\/\(\((z\^2 - 1)\)\ \((z\^8 + z\^4 + 1)\)\)\) \[DifferentialD]z\) - \(972\/\[Pi]\^3\) \(\[Integral]\_0\%1 \(\(\ z\^4\ \(\(log\^2\)(z)\)\)\/\(\((z\^2 + 1)\)\ \((z\^4 - z\^2 + 1)\)\^2\)\) \[DifferentialD]z\)\)}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm], Cell[TextData[{ "All these integrals can be evaluated straightforwardly by ", StyleBox["Mathematica", FontSlant->"Italic"], ". \nRemarkably, it turns out the first integral (which looks more arduous \ than the others) is easily done:" }], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\[Integral]\_0\%1 \(\( z\ \((z\^8 + 7\ z\^6 - 7\ z\^2 - 1)\)\ ArcTan[\(z\^3 - 1\)\/\(z\^3 + 1\)]\)\/\((z\^4 + z\^2 + 1)\)\^3 \) \[DifferentialD]z\)], "Input", CellLabel->"In[4]:="], Cell[BoxData[ \(TraditionalForm\`\[Pi]\/16\)], "Output", CellLabel->"Out[4]="] }, Open ]], Cell[TextData[{ "Running the other three integrals straightforwardly in ", StyleBox["Mathematica ", FontSlant->"Italic"], "yields large expressions containing dilogarithm functions. It requires \ some effort to simplify these expressions. However, we can streamline the \ integration procedure if we observe that the denominators of the integrands \ are products of cyclotomic polynomials. It is well known that, if ", Cell[BoxData[ \(TraditionalForm\`f(z)\)]], " is cyclotomic, then by definition it divides some ", Cell[BoxData[ \(TraditionalForm\`1 - z\^n\)]], ". Hence the following integrand, for example, can be rewritten as" }], "Text"], Cell[BoxData[ \(TraditionalForm \`\(z\^4\ \(log(z)\)\)\/\(\((z\^2 - 1)\)\ \((z\^8 + z\^4 + 1)\)\) = \(\((1 - z)\)\^2\ \((z + 1)\)\^2\ \((z\^2 + 1)\)\^3\ \(log(z)\)\)\/\(3\ \((1 - z\^12)\)\) - \(log(z)\)\/\(3\ \((1 - z\^2)\)\)\)], "Text", TextAlignment->Center, TextJustification->0], Cell["Then, taking into account that:", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\[Integral]\_0\%1\(\( z\^k\ Log[z]\)\/\(1 - z\^p\)\) \[DifferentialD]z\)], "Input"], Cell[BoxData[ FormBox[ RowBox[{"If", "[", RowBox[{\(Re(k) > \(-1\) \[And] Re(p) > 0\), ",", RowBox[{"-", FractionBox[ RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(\(k + 1\)\/p\), ")"}], \(p\^2\)]}], ",", \(\[Integral]\_0\%1\(\( z\^k\ \(log(z)\)\)\/\(1 - z\^p\)\) \[DifferentialD]z\)}], "]"}], TraditionalForm]], "Output", CellLabel->"Out[1]="] }, Open ]], Cell["we have:", "Text"], Cell[BoxData[ FormBox[ RowBox[{ \(\(648\/\[Pi]\^3\) \(\[Integral]\_0\%1 \(\(\ z\^4\ \(log(z)\)\)\/\(\((z\^2 - 1)\)\ \((z\^8 + z\^4 + 1)\)\)\) \[DifferentialD]z\)\), "=", RowBox[{\(-\(3\/\[Pi]\)\), "-", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(1\/12\), ")"}]}], \(4\ \[Pi]\^3\)], "-", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(11\/12\), ")"}]}], \(4\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(1\/3\), ")"}]}], \(2\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(2\/3\), ")"}]}], \(2\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"3", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(1\/6\), ")"}]}], \(2\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"3", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(5\/6\), ")"}]}], \(2\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(5\/12\), ")"}]}], \(4\ \[Pi]\^3\)], "+", FractionBox[ RowBox[{"9", " ", RowBox[{ SuperscriptBox[ TagBox["\[Psi]", PolyGamma], \((1)\)], "(", \(7\/12\), ")"}]}], \(4\ \[Pi]\^3\)]}]}], TraditionalForm]], "Text"], Cell["\<\ Performing further simplifications of polygamma functions, we \ obtain:\ \>", "Text"], Cell[BoxData[ \(TraditionalForm \`\(648\/\[Pi]\^3\) \(\[Integral]\_0\%1 \(\(\ z\^4\ \(log(z)\)\)\/\(\((z\^2 - 1)\)\ \((z\^8 + z\^4 + 1)\)\)\) \[DifferentialD]z\) = \(9\ \((1 - 2\ \@3)\)\)\/\[Pi]\)], "Text", TextAlignment->Center, TextJustification->0], Cell["\<\ Finally, repeating the same procedure for other integrals in the \ formula (13), we arrive at:\ \>", "Text"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FractionBox[\(\[PartialD]\^2\(\( P\_2\)(y)\)\), \(\[PartialD]y\^2\), MultilineFunction->None], SubscriptBox[ StyleBox["\[VerticalSeparator]", SpanMinSize->3], \(y = \[Pi]\/3\)]}], "=", \(\(-\(\(35\ \)\/2\)\) + 54\/\[Pi]\)}], TraditionalForm]], "NumberedEquation", TextAlignment->Center, TextJustification->0, FormatType->StandardForm] }, Closed]], Cell[CellGroupData[{ Cell["References", "Section"], Cell[TextData[{ "1. F. Y. Wu, ", StyleBox["The Potts model", FontSlant->"Italic"], ", Reviews of Modern Physics, 54 (1982) 235-268.." }], "Reference"], Cell[TextData[{ "2. H. N. V. Temperley and E. H. Lieb, ", StyleBox[ "Relations between the 'percolation' and 'colouring' problem and other \ graph-theoretical problems associated with regular planar lattices; some \ exact results for the 'percolation' problem", FontSlant->"Italic"], ", Proc. Royal Soc. London A 322 (1971) 251-280." }], "Reference"], Cell[TextData[{ "3. H. N. V. Temperley, ", StyleBox["Graph Theory and Applications", FontSlant->"Italic"], ", Ellis Horwood Ltd. 1981" }], "Reference"], Cell[TextData[{ "4. N. Biggs, ", StyleBox["Algebraic Graph Theory", FontSlant->"Italic"], ", Cambridge Univ. Press, 2nd ed. 1993." }], "Reference"], Cell[TextData[{ "5. S. Finch, ", StyleBox["Random percolation constants", FontSlant->"Italic"], ", HTML essay at World Wide Web URL \ http://www.mathsoft.com/asolve/constant/rndprc/rndprc.html, \npart of \ Favorite Mathematical Constants collection (1996)." }], "Reference"] }, Closed]] }, Open ]] }, FrontEndVersion->"NeXT 3.0", ScreenRectangle->{{0, 1053}, {0, 832}}, WindowToolbars->{}, Evaluator->"Local", WindowSize->{662, 675}, WindowMargins->{{Automatic, 161}, {Automatic, 5}}, PrintingCopies->1, PrintingPageRange->{1, 10}, ShowCellLabel->False, Magnification->1.25, StyleDefinitions -> "ArticleClassic.nb" ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1731, 51, 103, 2, 123, "Title"], Cell[1837, 55, 145, 8, 67, "Names"], Cell[CellGroupData[{ Cell[2007, 67, 31, 0, 66, "Section"], Cell[2041, 69, 571, 9, 70, "Text"], Cell[2615, 80, 474, 8, 70, "Text"], Cell[3092, 90, 137, 3, 70, "Text"], Cell[3232, 95, 543, 18, 70, "Text"], Cell[3778, 115, 125, 3, 70, "Text"], Cell[3906, 120, 433, 14, 70, "Text"], Cell[4342, 136, 640, 15, 70, "Text"], Cell[4985, 153, 275, 7, 70, "Text"], Cell[5263, 162, 367, 10, 70, "NumberedEquation"], Cell[5633, 174, 143, 3, 70, "Text"], Cell[5779, 179, 419, 11, 70, "NumberedEquation"], Cell[6201, 192, 362, 12, 70, "Text"], Cell[CellGroupData[{ Cell[6588, 208, 267, 7, 70, "Input"], Cell[6858, 217, 14415, 520, 70, 5950, 410, "GraphicsData", "PostScript", "Graphics"], Cell[21276, 739, 192, 5, 70, "Output"] }, Open ]], Cell[21483, 747, 781, 21, 70, "Text"], Cell[22267, 770, 515, 13, 70, "NumberedEquation"], Cell[22785, 785, 158, 7, 70, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[22980, 797, 35, 0, 33, "Section"], Cell[23018, 799, 131, 3, 70, "Text"], Cell[CellGroupData[{ Cell[23174, 806, 81, 1, 70, "Subsubsection"], Cell[23258, 809, 208, 6, 70, "Text"], Cell[23469, 817, 423, 11, 70, "Text"], Cell[23895, 830, 57, 0, 70, "Text"], Cell[23955, 832, 350, 10, 70, "Text"], Cell[24308, 844, 18, 0, 70, "Text"], Cell[24329, 846, 504, 13, 70, "Text"], Cell[24836, 861, 179, 6, 70, "Text"], Cell[25018, 869, 316, 8, 70, "Text"], Cell[25337, 879, 332, 8, 70, "Text"], Cell[25672, 889, 26, 0, 70, "Text"], Cell[25701, 891, 301, 9, 70, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[26039, 905, 86, 1, 70, "Subsubsection"], Cell[26128, 908, 66, 0, 70, "Text"], Cell[26197, 910, 890, 21, 70, "Text"], Cell[27090, 933, 94, 3, 70, "Text"], Cell[CellGroupData[{ Cell[27209, 940, 153, 4, 70, "Input"], Cell[27365, 946, 253, 5, 70, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27655, 956, 117, 3, 70, "Input"], Cell[27775, 961, 88, 2, 70, "Output"] }, Open ]], Cell[27878, 966, 177, 5, 70, "Text"], Cell[28058, 973, 487, 14, 70, "Text"], Cell[28548, 989, 88, 3, 70, "Text"], Cell[28639, 994, 574, 15, 70, "Equation"], Cell[29216, 1011, 41, 0, 70, "Text"], Cell[29260, 1013, 489, 14, 70, "Equation"], Cell[29752, 1029, 24, 0, 70, "Text"], Cell[29779, 1031, 319, 10, 70, "Equation"] }, Closed]], Cell[CellGroupData[{ Cell[30135, 1046, 86, 1, 70, "Subsubsection"], Cell[30224, 1049, 48, 0, 70, "Text"], Cell[30275, 1051, 305, 8, 70, "Equation"], Cell[30583, 1061, 175, 5, 70, "Text"], Cell[30761, 1068, 389, 10, 70, "Equation"], Cell[31153, 1080, 38, 0, 70, "Text"], Cell[CellGroupData[{ Cell[31216, 1084, 172, 4, 70, "Input"], Cell[31391, 1090, 260, 6, 70, "Output"] }, Open ]], Cell[31666, 1099, 41, 0, 70, "Text"], Cell[31710, 1101, 724, 19, 70, "NumberedEquation"], Cell[32437, 1122, 29, 0, 70, "Text"], Cell[32469, 1124, 491, 15, 70, "Equation"], Cell[32963, 1141, 38, 0, 70, "Text"], Cell[33004, 1143, 221, 6, 70, "Equation"], Cell[33228, 1151, 123, 5, 70, "Text"], Cell[33354, 1158, 555, 15, 70, "Equation"], Cell[33912, 1175, 41, 0, 70, "Text"], Cell[33956, 1177, 290, 8, 70, "Equation"], Cell[34249, 1187, 38, 0, 70, "Text"], Cell[34290, 1189, 440, 13, 70, "Equation"], Cell[34733, 1204, 24, 0, 70, "Text"], Cell[34760, 1206, 138, 3, 70, "Equation"] }, Closed]], Cell[CellGroupData[{ Cell[34935, 1214, 86, 1, 70, "Subsubsection"], Cell[35024, 1217, 253, 11, 70, "Text"], Cell[35280, 1230, 270, 7, 70, "Equation"], Cell[35553, 1239, 175, 5, 70, "Text"], Cell[CellGroupData[{ Cell[35753, 1248, 127, 3, 70, "Input"], Cell[35883, 1253, 182, 6, 70, "Output"] }, Open ]], Cell[36080, 1262, 21, 0, 70, "Text"], Cell[36104, 1264, 262, 8, 70, "Equation"] }, Closed]], Cell[CellGroupData[{ Cell[36403, 1277, 93, 1, 70, "Subsubsection"], Cell[36499, 1280, 69, 0, 70, "Text"], Cell[36571, 1282, 575, 14, 70, "Text"], Cell[CellGroupData[{ Cell[37171, 1300, 173, 4, 70, "Input"], Cell[37347, 1306, 248, 6, 70, "Output"] }, Open ]], Cell[37610, 1315, 68, 0, 70, "Text"], Cell[37681, 1317, 485, 15, 70, "Text"], Cell[38169, 1334, 39, 0, 70, "Text"], Cell[38211, 1336, 473, 14, 70, "Text"], Cell[38687, 1352, 21, 0, 70, "Text"], Cell[38711, 1354, 286, 7, 70, "Text"], Cell[39000, 1363, 19, 0, 70, "Text"], Cell[39022, 1365, 510, 14, 70, "Text"], Cell[39535, 1381, 26, 0, 70, "Text"], Cell[39564, 1383, 144, 4, 70, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[39745, 1392, 83, 1, 70, "Subsubsection"], Cell[39831, 1395, 52, 0, 70, "Text"], Cell[39886, 1397, 283, 8, 70, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[40206, 1410, 93, 1, 70, "Subsubsection"], Cell[40302, 1413, 47, 0, 70, "Text"], Cell[40352, 1415, 543, 13, 70, "Text"], Cell[40898, 1430, 60, 0, 70, "Text"], Cell[CellGroupData[{ Cell[40983, 1434, 408, 8, 70, "Input"], Cell[41394, 1444, 601, 15, 70, "Output"] }, Open ]], Cell[42010, 1462, 41, 0, 70, "Text"], Cell[42054, 1464, 482, 13, 70, "Text"], Cell[42539, 1479, 26, 0, 70, "Text"], Cell[42568, 1481, 330, 9, 70, "Text"] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[42947, 1496, 35, 0, 33, "Section"], Cell[42985, 1498, 298, 7, 70, "Text"], Cell[43286, 1507, 532, 14, 70, "NumberedEquation"], Cell[43821, 1523, 89, 3, 70, "Text"], Cell[43913, 1528, 182, 5, 70, "NumberedEquation"], Cell[44098, 1535, 71, 0, 70, "Text"], Cell[44172, 1537, 172, 5, 70, "Text"], Cell[44347, 1544, 763, 18, 70, "NumberedEquation"], Cell[45113, 1564, 169, 6, 70, "Text"], Cell[CellGroupData[{ Cell[45307, 1574, 168, 4, 70, "Input"], Cell[45478, 1580, 840, 23, 70, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[46355, 1608, 72, 2, 70, "Input"], Cell[46430, 1612, 85, 2, 70, "Output"] }, Open ]], Cell[46530, 1617, 193, 5, 70, "Text"], Cell[46726, 1624, 511, 13, 70, "NumberedEquation"], Cell[47240, 1639, 110, 5, 70, "Text"], Cell[CellGroupData[{ Cell[47375, 1648, 256, 6, 70, "Input"], Cell[47634, 1656, 104, 2, 70, "Output"] }, Open ]], Cell[47753, 1661, 40, 0, 70, "Text"], Cell[47796, 1663, 474, 13, 70, "NumberedEquation"], Cell[48273, 1678, 58, 0, 70, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[48368, 1683, 36, 0, 33, "Section"], Cell[48407, 1685, 242, 5, 70, "Text"], Cell[48652, 1692, 40, 0, 70, "Text"], Cell[48695, 1694, 417, 11, 70, "NumberedEquation"], Cell[49115, 1707, 57, 0, 70, "Text"], Cell[49175, 1709, 133, 5, 70, "Text"], Cell[49311, 1716, 725, 19, 70, "NumberedEquation"], Cell[50039, 1737, 218, 6, 70, "Text"], Cell[50260, 1745, 1365, 32, 70, "NumberedEquation"], Cell[51628, 1779, 165, 5, 70, "Text"], Cell[51796, 1786, 1375, 32, 70, "NumberedEquation"], Cell[53174, 1820, 256, 6, 70, "Text"], Cell[CellGroupData[{ Cell[53455, 1830, 239, 5, 70, "Input"], Cell[53697, 1837, 84, 2, 70, "Output"] }, Open ]], Cell[53796, 1842, 673, 14, 70, "Text"], Cell[54472, 1858, 314, 6, 70, "Text"], Cell[54789, 1866, 47, 0, 70, "Text"], Cell[CellGroupData[{ Cell[54861, 1870, 112, 2, 70, "Input"], Cell[54976, 1874, 523, 13, 70, "Output"] }, Open ]], Cell[55514, 1890, 24, 0, 70, "Text"], Cell[55541, 1892, 2279, 63, 70, "Text"], Cell[57823, 1957, 95, 3, 70, "Text"], Cell[57921, 1962, 310, 8, 70, "Text"], Cell[58234, 1972, 118, 3, 70, "Text"], Cell[58355, 1977, 494, 15, 70, "NumberedEquation"] }, Closed]], Cell[CellGroupData[{ Cell[58886, 1997, 29, 0, 33, "Section"], Cell[58918, 1999, 160, 5, 70, "Reference"], Cell[59081, 2006, 362, 8, 70, "Reference"], Cell[59446, 2016, 161, 5, 70, "Reference"], Cell[59610, 2023, 157, 5, 70, "Reference"], Cell[59770, 2030, 284, 7, 70, "Reference"] }, Closed]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)