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

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

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     31764,        907]
NotebookOptionsPosition[     29030,        860]
NotebookOutlinePosition[     29367,        875]
CellTagsIndexPosition[     29324,        872]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["ECMD sr-LDA functional [Paziani et al. PRB 73 (2006) 155111]", "Title",
 CellChangeTimes->{{3.7717615442945957`*^9, 3.7717615745751762`*^9}, {
  3.7717616555034122`*^9, 3.771761658578969*^9}, {3.771761889119907*^9, 
  3.7717618899605417`*^9}},ExpressionUUID->"cc8f10f6-dee0-4be6-9b95-\
a2b645508bc1"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubsuperscriptBox["\[Epsilon]", "c", "LR"], "[", 
   RowBox[{
    SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], "]"}], "\[Equal]", 
  RowBox[{
   RowBox[{
    SuperscriptBox["\[Epsilon]", "LR"], "[", 
    RowBox[{
     SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], "]"}], "-", 
   RowBox[{
    SubscriptBox["t", "s"], "[", 
    RowBox[{
     SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], "]"}], "-", 
   RowBox[{
    SubsuperscriptBox["\[Epsilon]", "x", "LR"], "[", 
    RowBox[{
     SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7717614891109467`*^9, 
  3.771761526234874*^9}},ExpressionUUID->"f09b5b1d-fcfa-4277-b055-\
8d013d150e51"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubscriptBox["t", "s"], "[", 
   RowBox[{
    SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], "]"}], "=", 
  RowBox[{
   FractionBox[
    RowBox[{"3", 
     SubsuperscriptBox["k", "F", "2"]}], "10"], 
   RowBox[{
    SubscriptBox["\[Phi]", "5"], "[", "\[Zeta]", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.771761595434334*^9, 
  3.77176162310323*^9}},ExpressionUUID->"2adb86e5-1613-4560-afb2-\
c874c93162cc"],

Cell[BoxData[
 RowBox[{
  SubscriptBox["k", "F"], "\[Equal]", 
  RowBox[{
   FractionBox["1", 
    RowBox[{"\[Alpha]", " ", 
     SubscriptBox["r", "s"]}]], "\t", "\[Alpha]"}], "\[Equal]", 
  SuperscriptBox[
   RowBox[{"(", 
    FractionBox["4", 
     RowBox[{"9", "\[Pi]"}]], ")"}], 
   RowBox[{"1", "/", "3"}]]}]], "Input",
 CellChangeTimes->{{3.771761629813285*^9, 3.771761651073*^9}, 
   3.771761892459922*^9},ExpressionUUID->"e471ce7c-fc48-482e-8bea-\
3e9e0179671c"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubscriptBox["\[Epsilon]", "x"], "[", 
   RowBox[{
    SubscriptBox["r", "s"], ",", "\[Zeta]", ",", "\[Mu]"}], "]"}], "\[Equal]", 
  FractionBox[
   RowBox[{
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "+", "\[Zeta]"}], ")"}], 
      RowBox[{"4", "/", "3"}]], 
     RowBox[{
      SubscriptBox["f", "x"], "[", 
      RowBox[{
       SubscriptBox["r", "s"], ",", 
       RowBox[{"\[Mu]", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "+", "\[Zeta]"}], ")"}], 
         RowBox[{
          RowBox[{"-", "1"}], "/", "3"}]]}]}], "]"}]}], "+", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "-", "\[Zeta]"}], ")"}], 
      RowBox[{"4", "/", "3"}]], 
     RowBox[{
      SubscriptBox["f", "x"], "[", 
      RowBox[{
       SubscriptBox["r", "s"], ",", 
       RowBox[{"\[Mu]", 
        SuperscriptBox[
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Zeta]"}], ")"}], 
         RowBox[{
          RowBox[{"-", "1"}], "/", "3"}]]}]}], "]"}]}]}], "2"]}]], "Input",
 CellChangeTimes->{{3.771761705400154*^9, 
  3.771761802293672*^9}},ExpressionUUID->"d1a8b6ea-55ae-4a5e-8673-\
b9734def02ca"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubscriptBox["f", "x"], "[", 
   RowBox[{
    SubscriptBox["r", "s"], ",", "\[Mu]"}], "]"}], "\[Equal]", 
  RowBox[{
   RowBox[{"-", 
    FractionBox["\[Mu]", "\[Pi]"]}], 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"2", "y"}], "-", 
        RowBox[{"4", 
         SuperscriptBox["y", "3"]}]}], ")"}], 
      RowBox[{"Exp", "[", 
       RowBox[{"-", 
        FractionBox["1", 
         RowBox[{"4", 
          SuperscriptBox["y", "2"]}]]}], "]"}]}], "-", 
     RowBox[{"3", "y"}], "+", 
     RowBox[{"4", 
      SuperscriptBox["y", "3"]}], "+", 
     RowBox[{
      SqrtBox["\[Pi]"], 
      RowBox[{"Erf", "[", 
       FractionBox["1", 
        RowBox[{"2", "y"}]], "]"}]}]}], ")"}], "\t", "y"}], "\[Equal]", 
  FractionBox[
   RowBox[{"\[Mu]", " ", "\[Alpha]", " ", 
    SubscriptBox["r", "s"]}], "2"]}]], "Input",
 CellChangeTimes->{{3.771761806462194*^9, 
  3.7717619043772717`*^9}},ExpressionUUID->"8257931c-d5a2-4fa5-bfcf-\
242ed662fad5"],

Cell[CellGroupData[{

Cell["Analytic representation of the correlation energy", "Section",
 CellChangeTimes->{{3.7717619451626053`*^9, 3.7717619604379997`*^9}, {
  3.7717628547850637`*^9, 
  3.771762856014763*^9}},ExpressionUUID->"da7db347-895a-4445-adf0-\
185e3e7ae4bc"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Epsilon]cLR", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_", ",", "\[Mu]_"}], "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{
      RowBox[{
       SuperscriptBox[
        RowBox[{
         SubscriptBox["\[Phi]", "2"], "[", "\[Zeta]", "]"}], "3"], 
       RowBox[{"Q", "[", 
        FractionBox[
         RowBox[{"\[Mu]", " ", 
          SqrtBox["rs"]}], 
         RowBox[{
          SubscriptBox["\[Phi]", "2"], "[", "\[Zeta]", "]"}]], "]"}]}], "+", 
      RowBox[{
       RowBox[{"a1", "[", 
        RowBox[{"rs", ",", "\[Zeta]"}], "]"}], 
       SuperscriptBox["\[Mu]", "3"]}], "+", 
      RowBox[{
       RowBox[{"a2", "[", 
        RowBox[{"rs", ",", "\[Zeta]"}], "]"}], 
       SuperscriptBox["\[Mu]", "4"]}], "+", 
      RowBox[{
       RowBox[{"a3", "[", 
        RowBox[{"rs", ",", "\[Zeta]"}], "]"}], 
       SuperscriptBox["\[Mu]", "5"]}], "+", 
      RowBox[{
       RowBox[{"a4", "[", 
        RowBox[{"rs", ",", "\[Zeta]"}], "]"}], 
       SuperscriptBox["\[Mu]", "6"]}], "+", 
      RowBox[{
       RowBox[{"a5", "[", 
        RowBox[{"rs", ",", "\[Zeta]"}], "]"}], 
       SuperscriptBox["\[Mu]", "8"]}]}], 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{
         SuperscriptBox[
          RowBox[{"b0", "[", "rs", "]"}], "2"], 
         SuperscriptBox["\[Mu]", "2"]}]}], ")"}], "4"]], "/.", 
    RowBox[{"\[Alpha]", "\[Rule]", 
     SuperscriptBox[
      RowBox[{"(", 
       FractionBox["4", 
        RowBox[{"9", "\[Pi]"}]], ")"}], 
      RowBox[{"1", "/", "3"}]]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.771761963382757*^9, 3.77176203397927*^9}, {
   3.771762098820107*^9, 3.771762142266923*^9}, {3.771762856723877*^9, 
   3.7717628631754007`*^9}, {3.771778456441491*^9, 3.7717784599917917`*^9}, {
   3.77177856335856*^9, 3.7717785984922667`*^9}, {3.771778747941729*^9, 
   3.771778770864193*^9}, {3.771779376690545*^9, 3.771779376831264*^9}, {
   3.771779500268244*^9, 3.771779535341432*^9}, 3.771781649390931*^9, {
   3.771782124540922*^9, 3.771782126643386*^9}, 3.771836550551532*^9},
 CellLabel->"In[1]:=",ExpressionUUID->"2424c22e-2944-42fb-a926-6931d6b950db"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    SubscriptBox["\[Phi]", "n_"], "[", "\[Zeta]_", "]"}], ":=", 
   FractionBox[
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "+", "\[Zeta]"}], ")"}], 
      RowBox[{"n", "/", "3"}]], "+", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"1", "-", "\[Zeta]"}], ")"}], 
      RowBox[{"n", "/", "3"}]]}], "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.771761663679653*^9, 3.771761699635745*^9}, {
   3.771778463326877*^9, 3.7717784697876587`*^9}, 3.7718365533847017`*^9},
 CellLabel->"In[2]:=",ExpressionUUID->"5e7071ef-4b10-4195-9ab5-72086ae9360d"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Q", "[", "x_", "]"}], ":=", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"2", 
         RowBox[{"Log", "[", "2", "]"}]}], "-", "2"}], 
       SuperscriptBox["\[Pi]", "2"]], 
      RowBox[{"Log", "[", 
       FractionBox[
        RowBox[{"1", "+", 
         RowBox[{"a", " ", "x"}], "+", 
         RowBox[{"b", " ", 
          SuperscriptBox["x", "2"]}], "+", 
         RowBox[{"c", " ", 
          SuperscriptBox["x", "3"]}]}], 
        RowBox[{"1", "+", 
         RowBox[{"a", " ", "x"}], "+", 
         RowBox[{"d", " ", 
          SuperscriptBox["x", "2"]}]}]], "]"}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"a", "\[Rule]", "5.84605"}], ",", 
       RowBox[{"c", "\[Rule]", "3.91744"}], ",", 
       RowBox[{"d", "\[Rule]", "3.44851"}], ",", 
       RowBox[{"b", "->", 
        RowBox[{"3.44851", "-", 
         FractionBox[
          RowBox[{"3", "\[Pi]", " ", "\[Alpha]"}], 
          RowBox[{
           RowBox[{"4", 
            RowBox[{"Log", "[", "2", "]"}]}], "-", "4"}]]}]}]}], "}"}]}]}], 
   ";"}], "\t"}]], "Input",
 CellChangeTimes->{{3.771762166139291*^9, 3.7717621744945383`*^9}, {
   3.7717622073456*^9, 3.771762234218986*^9}, {3.771762346733646*^9, 
   3.771762414996784*^9}, {3.771762878996976*^9, 3.771762881748227*^9}, {
   3.7717784794941*^9, 3.771778523118412*^9}, 3.771836554184534*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"431fbccd-665c-479c-825d-74e37905a953"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"a1", "[", 
     RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"4", 
      SuperscriptBox[
       RowBox[{"b0", "[", "rs", "]"}], "6"], 
      RowBox[{"C3", "[", 
       RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], "+", 
     RowBox[{
      SuperscriptBox[
       RowBox[{"b0", "[", "rs", "]"}], "8"], 
      RowBox[{"C5", "[", 
       RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}]}]}], ";"}], 
  "\t"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"a2", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"4", 
     SuperscriptBox[
      RowBox[{"b0", "[", "rs", "]"}], "6"], 
     RowBox[{"C2", "[", 
      RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], "+", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"b0", "[", "rs", "]"}], "8"], 
     RowBox[{"C4", "[", 
      RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], "+", 
    RowBox[{"6", 
     SuperscriptBox[
      RowBox[{"b0", "[", "rs", "]"}], "4"], "ec"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"a3", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"b0", "[", "rs", "]"}], "8"], 
    RowBox[{"C3", "[", 
     RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"a4", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     SuperscriptBox[
      RowBox[{"b0", "[", "rs", "]"}], "8"], 
     RowBox[{"C2", "[", 
      RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], "+", 
    RowBox[{"4", 
     SuperscriptBox[
      RowBox[{"b0", "[", "rs", "]"}], "6"], "ec"}]}]}], ";", "\t", 
  RowBox[{
   RowBox[{"a5", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"b0", "[", "rs", "]"}], "8"], "ec"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.771762952203158*^9, 3.771762957939657*^9}, {
  3.771762998580441*^9, 3.771763110694063*^9}, {3.7717785319110537`*^9, 
  3.771778539064945*^9}, {3.771778710580977*^9, 3.7717787447614594`*^9}, {
  3.771778782889411*^9, 3.771778885987258*^9}, {3.7717793542622547`*^9, 
  3.771779361369142*^9}, {3.771780821220508*^9, 3.771780824966153*^9}, {
  3.771781783385194*^9, 3.771781793877317*^9}, {3.771836556016198*^9, 
  3.7718365584308977`*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"68450462-8882-467d-85de-668cca5b67df"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"C2", "[", 
     RowBox[{"\[Zeta]_", ",", "rs_"}], "]"}], ":=", 
    RowBox[{"-", 
     FractionBox[
      RowBox[{"3", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         SuperscriptBox["\[Zeta]", "2"]}], ")"}], 
       RowBox[{"gc", "[", "rs", "]"}]}], 
      RowBox[{"8", " ", 
       SuperscriptBox["rs", "3"]}]]}]}], ";"}], "\t"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"C3", "[", 
     RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"1", "-", 
        SuperscriptBox["\[Zeta]", "2"]}], ")"}]}], 
     FractionBox[
      RowBox[{"g", "[", "rs", "]"}], 
      RowBox[{
       SqrtBox[
        RowBox[{"2", "\[Pi]"}]], " ", 
       SuperscriptBox["rs", "3"]}]]}]}], ";"}], "\t"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"C4", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{"9", 
      RowBox[{"c4", "[", 
       RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], 
     RowBox[{"64", 
      SuperscriptBox["rs", "3"]}]]}]}], ";", "\t", 
  RowBox[{
   RowBox[{"C5", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{"-", 
    FractionBox[
     RowBox[{"9", 
      RowBox[{"c5", "[", 
       RowBox[{"rs", ",", "\[Zeta]"}], "]"}]}], 
     RowBox[{"40", 
      SqrtBox[
       RowBox[{"2", "\[Pi]"}]], 
      SuperscriptBox["rs", "3"]}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"gc", "[", "rs_", "]"}], ":=", 
   RowBox[{
    RowBox[{"g", "[", "rs", "]"}], "-", 
    FractionBox["1", "2"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g", "[", "rs_", "]"}], ":=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.7317", "-", "0.752411"}], ")"}], "rs"}], "+", " ", 
      RowBox[{"0.0819306", 
       SuperscriptBox["rs", "2"]}], "-", 
      RowBox[{"0.0127713", 
       SuperscriptBox["rs", "3"]}], "+", 
      RowBox[{"0.00185898", 
       SuperscriptBox["rs", "4"]}]}], ")"}], 
    RowBox[{
     RowBox[{"Exp", "[", 
      RowBox[{
       RowBox[{"-", "0.752411"}], "rs"}], "]"}], "/", "2"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.771763124087843*^9, 3.771763239782792*^9}, {
  3.771778546103352*^9, 3.7717785482662897`*^9}, {3.771778611303747*^9, 
  3.771778616852309*^9}, {3.7717788897975807`*^9, 3.7717789557867327`*^9}, {
  3.7717790425591583`*^9, 3.771779055944982*^9}, {3.771779133030603*^9, 
  3.771779133243347*^9}, {3.771779255825417*^9, 3.77177928375259*^9}, {
  3.771779409854127*^9, 3.7717794246671057`*^9}, {3.771779529748591*^9, 
  3.7717795315356283`*^9}, {3.771779735811041*^9, 3.771779874568338*^9}, {
  3.7718365594479647`*^9, 3.771836570743126*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"76293f49-0f9c-448d-9fa4-7a7821619364"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c4", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      FractionBox[
       RowBox[{"1", "+", 
        SuperscriptBox["\[Zeta]", "2"]}], "2"], ")"}], 
     RowBox[{"gpp", "[", 
      RowBox[{"rs", 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox["2", 
          RowBox[{"1", "+", "\[Zeta]"}]], ")"}], 
        RowBox[{"1", "/", "3"}]]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      FractionBox[
       RowBox[{"1", "-", 
        SuperscriptBox["\[Zeta]", "2"]}], "2"], ")"}], 
     RowBox[{"gpp", "[", 
      RowBox[{"rs", 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox["2", 
          RowBox[{"1", "-", "\[Zeta]"}]], ")"}], 
        RowBox[{"1", "/", "3"}]]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox["\[Zeta]", "2"]}], ")"}], 
     RowBox[{"D2", "[", "rs", "]"}]}], "-", 
    FractionBox[
     RowBox[{
      SubscriptBox["\[Phi]", "8"], "[", "\[Zeta]", "]"}], 
     RowBox[{"5", 
      SuperscriptBox["\[Alpha]", "2"], 
      SuperscriptBox["rs", "2"]}]]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.77176327770026*^9, 3.771763352798826*^9}, {
   3.771763509128861*^9, 3.7717635428199377`*^9}, {3.7717786252373247`*^9, 
   3.7717786421091213`*^9}, {3.77177867451056*^9, 3.771778705676589*^9}, {
   3.771778959529389*^9, 3.7717789703788347`*^9}, 3.771779326066607*^9, {
   3.7717799304393873`*^9, 3.771779934301042*^9}, 3.771836571928794*^9},
 CellLabel->"In[13]:=",ExpressionUUID->"7b438198-7fff-4469-8902-d63008ab1edd"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c5", "[", 
    RowBox[{"rs_", ",", "\[Zeta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      FractionBox[
       RowBox[{"1", "+", 
        SuperscriptBox["\[Zeta]", "2"]}], "2"], ")"}], 
     RowBox[{"gpp", "[", 
      RowBox[{"rs", 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox["2", 
          RowBox[{"1", "+", "\[Zeta]"}]], ")"}], 
        RowBox[{"1", "/", "3"}]]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      FractionBox[
       RowBox[{"1", "-", 
        SuperscriptBox["\[Zeta]", "2"]}], "2"], ")"}], 
     RowBox[{"gpp", "[", 
      RowBox[{"rs", 
       SuperscriptBox[
        RowBox[{"(", 
         FractionBox["2", 
          RowBox[{"1", "-", "\[Zeta]"}]], ")"}], 
        RowBox[{"1", "/", "3"}]]}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox["\[Zeta]", "2"]}], ")"}], 
     RowBox[{"D3", "[", "rs", "]"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{
  3.7717634987152567`*^9, {3.771763552292667*^9, 3.771763554424954*^9}, {
   3.7717786290210857`*^9, 3.771778702617523*^9}, {3.771778974783494*^9, 
   3.771778978908309*^9}, {3.771779127030455*^9, 3.771779128290526*^9}, {
   3.7717793291115313`*^9, 3.771779329972632*^9}, {3.771779932467745*^9, 
   3.7717799359878817`*^9}, 3.771836573567986*^9},
 CellLabel->"In[14]:=",ExpressionUUID->"bb62c6dd-19f1-49b0-93ed-7a72089471a5"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"b0", "[", "rs_", "]"}], ":=", 
    RowBox[{"0.78949", "rs"}]}], ";"}], "\t"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"D2", "[", "rs_", "]"}], ":=", 
    RowBox[{
     FractionBox[
      RowBox[{"Exp", "[", 
       RowBox[{
        RowBox[{"-", "0.547"}], "rs"}], "]"}], 
      SuperscriptBox["rs", "2"]], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "0.388"}], "rs"}], "+", 
       RowBox[{"0.676", 
        SuperscriptBox["rs", "2"]}]}], ")"}]}]}], ";"}], 
  "\t"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"D3", "[", "rs_", "]"}], ":=", 
   RowBox[{
    FractionBox[
     RowBox[{"Exp", "[", 
      RowBox[{
       RowBox[{"-", "0.31"}], "rs"}], "]"}], 
     SuperscriptBox["rs", "3"]], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "4.95"}], "rs"}], "+", 
      SuperscriptBox["rs", "2"]}], ")"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.771763567753668*^9, 3.7717635711197557`*^9}, {
  3.7717636288569927`*^9, 3.771763659653902*^9}, {3.771764117143701*^9, 
  3.771764146373032*^9}, {3.7717642369975147`*^9, 3.771764249484836*^9}, {
  3.771778793333074*^9, 3.771778793734756*^9}, {3.771778908427208*^9, 
  3.771778922603982*^9}, {3.771778999237069*^9, 3.771779031517296*^9}, {
  3.771779119216881*^9, 3.771779125350945*^9}, {3.771836574383712*^9, 
  3.7718365786632423`*^9}},
 CellLabel->"In[15]:=",ExpressionUUID->"1b3202c0-1ab4-429c-8dba-a60935c07597"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"gpp", "[", "rs_", "]"}], ":=", 
   RowBox[{
    FractionBox[
     SuperscriptBox["2", 
      RowBox[{"5", "/", "3"}]], 
     RowBox[{"5", 
      SuperscriptBox["\[Alpha]", "2"], 
      SuperscriptBox["rs", "2"]}]], 
    FractionBox[
     RowBox[{"1", "-", 
      RowBox[{"0.02267", "rs"}]}], 
     RowBox[{"1", "+", 
      RowBox[{"0.4319", " ", "rs"}], "+", 
      RowBox[{"0.04", 
       SuperscriptBox["rs", "2"]}]}]]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.771764158041356*^9, 3.771764207830522*^9}, {
   3.771778645829051*^9, 3.7717786706143827`*^9}, {3.7717789848615627`*^9, 
   3.771778994476886*^9}, 3.771779335374116*^9, {3.771779511382004*^9, 
   3.77177951806966*^9}, {3.771779938347245*^9, 3.771779938762389*^9}, 
   3.7718365796717587`*^9},
 CellLabel->"In[18]:=",ExpressionUUID->"30aeeec3-f7a5-4138-bd8c-548696d2f752"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Epsilon]cLR", "[", 
   RowBox[{"1.0", ",", "0.2", ",", "0.5"}], "]"}], "/.", 
  RowBox[{"ec", "\[Rule]", 
   RowBox[{
    RowBox[{"-", "5.8963985436182129"}], " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "2"}]]}]}]}]], "Input",
 CellChangeTimes->{{3.771780839941061*^9, 3.77178084339225*^9}, {
  3.7717818058410683`*^9, 3.771781807442444*^9}, {3.7718366739261637`*^9, 
  3.7718367162615347`*^9}, {3.771836930343273*^9, 3.771836948258225*^9}},
 CellLabel->"In[30]:=",ExpressionUUID->"408e19ad-23db-4b1a-8d84-8acca9835fae"],

Cell[BoxData[
 RowBox[{"-", "0.02003816512886262`"}]], "Output",
 CellChangeTimes->{
  3.7717821311862507`*^9, {3.7718365235145597`*^9, 3.77183653965664*^9}, {
   3.771836584520483*^9, 3.77183659732627*^9}, {3.771836677187216*^9, 
   3.7718367167452602`*^9}, {3.771836933204549*^9, 3.7718369486178427`*^9}},
 CellLabel->"Out[30]=",ExpressionUUID->"0b40f4ad-6233-4b3d-8a80-07aca06e731a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Epsilon]cLR", "[", 
   RowBox[{"rs", ",", "z", ",", "mu"}], "]"}], "//", "FortranForm"}]], "Input",\

 CellChangeTimes->{{3.77178165542163*^9, 3.771781675040906*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"7d031781-f5af-40e6-b079-1bc0fdb02f2b"],

Cell["\<\
        (0.15092912957740764*ec*mu**8*rs**8 + \
(0.030106005566307166*mu**5*rs**5*
     -       (1 + 0.020711000000000035*rs + 0.0819306*rs**2 - 0.0127713*rs**3 \
+ 0.00185898*rs**4)*(-1 + z**2))/E**(0.752411*rs) + 
     -    mu**6*(0.9685895783716283*ec*rs**6 - (0.05659842359152786*rs**8*(1 \
- rs**2)*
     -          (-0.5 + (1 + 0.020711000000000035*z + 0.0819306*z**2 - \
0.0127713*z**3 + 0.00185898*z**4)/(2.*E**(0.752411*z))))/z**3) + 
     -    mu**3*((0.1932056675843198*rs**3*(1 + 0.020711000000000035*rs + \
0.0819306*rs**2 - 0.0127713*rs**3 + 0.00185898*rs**4)*(-1 + z**2))/
     -        E**(0.752411*rs) - 0.013547702504838225*rs**5*(((-4.95*rs + \
rs**2)*(1 - z**2))/(E**(0.31*rs)*rs**3) + 
     -          (3*1.5**0.3333333333333333*Pi**0.6666666666666666*(1 - \
0.028562410201116776*rs*(1/(1 - z))**0.3333333333333333)*(1 - z**2))/
     -           (10.*rs**2*(1 + 0.5441599014495957*rs*(1/(1 - \
z))**0.3333333333333333 + 0.06349604207872797*rs**2*(1/(1 - \
z))**0.6666666666666666)*
     -             (1/(1 - z))**0.6666666666666666) + \
(3*1.5**0.3333333333333333*Pi**0.6666666666666666*(1 + z**2)*
     -             (1 - 0.028562410201116776*rs*(1/(1 + \
z))**0.3333333333333333))/
     -           (10.*rs**2*(1/(1 + z))**0.6666666666666666*(1 + \
0.5441599014495957*rs*(1/(1 + z))**0.3333333333333333 + 
     -               0.06349604207872797*rs**2*(1/(1 + \
z))**0.6666666666666666)))) + 
     -    mu**4*(2.330975903948103*ec*rs**4 - (0.3632210918893606*rs**6*(1 - \
rs**2)*
     -          (-0.5 + (1 + 0.020711000000000035*z + 0.0819306*z**2 - \
0.0127713*z**3 + 0.00185898*z**4)/(2.*E**(0.752411*z))))/z**3 - 
     -       0.02122440884682295*rs**5*(((-0.388*rs + 0.676*rs**2)*(1 - \
z**2))/(E**(0.547*rs)*rs**2) + 
     -          (3*1.5**0.3333333333333333*Pi**0.6666666666666666*(1 - \
0.028562410201116776*rs*(1/(1 - z))**0.3333333333333333)*(1 - z**2))/
     -           (10.*rs**2*(1 + 0.5441599014495957*rs*(1/(1 - \
z))**0.3333333333333333 + 0.06349604207872797*rs**2*(1/(1 - \
z))**0.6666666666666666)*
     -             (1/(1 - z))**0.6666666666666666) + \
(3*1.5**0.3333333333333333*Pi**0.6666666666666666*(1 + z**2)*
     -             (1 - 0.028562410201116776*rs*(1/(1 + \
z))**0.3333333333333333))/
     -           (10.*rs**2*(1/(1 + z))**0.6666666666666666*(1 + \
0.5441599014495957*rs*(1/(1 + z))**0.3333333333333333 + 
     -               0.06349604207872797*rs**2*(1/(1 + \
z))**0.6666666666666666)) - 
     -          (3*1.5**0.3333333333333333*Pi**0.6666666666666666*((1 - \
z)**2.6666666666666665 + (1 + z)**2.6666666666666665))/(20.*rs**2))) + 
     -    (((1 - z)**0.6666666666666666 + (1 + z)**0.6666666666666666)**3*(-2 \
+ 2*Log(2))*
     -       Log((1 + (31.33952*mu**3*rs**1.5)/((1 - z)**0.6666666666666666 + \
(1 + z)**0.6666666666666666)**3 + 
     -           (29.798101530536215*mu**2*rs)/((1 - z)**0.6666666666666666 + \
(1 + z)**0.6666666666666666)**2 + 
     -           (11.6921*mu*Sqrt(rs))/((1 - z)**0.6666666666666666 + (1 + \
z)**0.6666666666666666))/
     -         (1 + (13.79404*mu**2*rs)/((1 - z)**0.6666666666666666 + (1 + \
z)**0.6666666666666666)**2 + 
     -           (11.6921*mu*Sqrt(rs))/((1 - z)**0.6666666666666666 + (1 + \
z)**0.6666666666666666))))/(8.*Pi**2))/(1 + 0.6232944601*mu**2*rs**2)**4\
\>", "Output",
 CellChangeTimes->{{3.771781666009734*^9, 3.7717816755992613`*^9}, 
   3.771781811362009*^9, 3.7717821333483133`*^9, 3.771836600647401*^9},
 CellLabel->
  "Out[20]//FortranForm=",ExpressionUUID->"e8c89351-e3be-455d-b993-\
40a6da00af8b"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["ECMD part", "Section",
 CellChangeTimes->{{3.771764286765162*^9, 
  3.7717642894676533`*^9}},ExpressionUUID->"1774903d-cd4a-427c-913b-\
6d03354a574f"],

Cell[BoxData[
 RowBox[{
  SubscriptBox["\[CapitalDelta]", 
   RowBox[{"LR", "-", "SR"}]], "=", 
  FractionBox[
   RowBox[{
    RowBox[{
     SubscriptBox["\[Delta]", "2"], " ", 
     SuperscriptBox["\[Mu]", "2"]}], "+", 
    RowBox[{
     SubscriptBox["\[Delta]", "3"], " ", 
     SuperscriptBox["\[Mu]", "3"]}], "+", 
    RowBox[{
     SubscriptBox["\[Delta]", "4"], " ", 
     SuperscriptBox["\[Mu]", "4"]}], "+", 
    RowBox[{
     SubscriptBox["\[Delta]", "5"], " ", 
     SuperscriptBox["\[Mu]", "5"]}], "+", 
    RowBox[{
     SubscriptBox["\[Delta]", "6"], " ", 
     SuperscriptBox["\[Mu]", "6"]}]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{
       SubsuperscriptBox["d", "0", "2"], 
       SuperscriptBox["\[Mu]", "2"]}]}], ")"}], "4"]]}]], "Input",
 CellChangeTimes->{{3.771764296622439*^9, 3.77176430804951*^9}, {
  3.771764345532308*^9, 
  3.771764387349852*^9}},ExpressionUUID->"a08e663f-6de1-46b7-a344-\
1b47fcd2e96e"],

Cell[BoxData[
 RowBox[{
  SubscriptBox["\[Delta]", "3"], "\[Equal]", 
  RowBox[{
   RowBox[{"4", 
    SubsuperscriptBox["d", "0", "6"], 
    SubscriptBox[
     OverscriptBox["C", "~"], "3"]}], "+", 
   RowBox[{
    SubsuperscriptBox["d", "0", "8"], 
    SubscriptBox[
     OverscriptBox["C", "~"], "5"], "\t", 
    SubscriptBox["\[Delta]", "4"]}]}], "\[Equal]", 
  RowBox[{
   RowBox[{"4", 
    SubsuperscriptBox["d", "0", "6"], 
    SubscriptBox["C", "2"]}], "+", 
   RowBox[{
    SubsuperscriptBox["d", "0", "8"], 
    SubscriptBox["C", "4"], "\t", 
    SubscriptBox["\[Delta]", "5"]}]}], "\[Equal]", 
  RowBox[{
   SubsuperscriptBox["d", "0", "8"], 
   SubscriptBox[
    OverscriptBox["C", "~"], "3"], "\t", 
   SubscriptBox["\[Delta]", "6"]}], "\[Equal]", 
  RowBox[{
   SubsuperscriptBox["d", "0", "8"], 
   SubscriptBox["C", "2"]}]}]], "Input",
 CellChangeTimes->{{3.771764399150547*^9, 
  3.77176451268692*^9}},ExpressionUUID->"760031e2-bb91-4ce7-8b6e-\
8566604d6f6c"],

Cell[BoxData[
 RowBox[{
  SubscriptBox[
   OverscriptBox["C", "~"], "3"], "\[Equal]", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{"1", "-", 
      SuperscriptBox["\[Zeta]", "2"]}], ")"}]}], 
   FractionBox[
    RowBox[{
     RowBox[{"g", "[", 
      RowBox[{"0", ",", 
       SubscriptBox["r", "s"], ",", "0"}], "]"}], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", 
        SqrtBox["2"]}], "-", "1"}], ")"}]}], 
    RowBox[{"2", 
     SqrtBox["\[Pi]"], 
     SubsuperscriptBox["r", "s", "3"]}]], "\t", 
   SubscriptBox[
    OverscriptBox["C", "~"], "5"]}], "\[Equal]", 
  RowBox[{"-", 
   FractionBox[
    RowBox[{"3", 
     RowBox[{
      SubscriptBox["c", "5"], "[", 
      RowBox[{
       SubscriptBox["r", "s"], ",", "\[Zeta]"}], "]"}], 
     RowBox[{"(", 
      RowBox[{"3", "-", 
       SqrtBox["2"]}], ")"}]}], 
    RowBox[{"20", 
     SqrtBox[
      RowBox[{"2", "\[Pi]"}]], 
     SubsuperscriptBox["r", "s", "3"]}]]}]}]], "Input",
 CellChangeTimes->{{3.771764516235433*^9, 
  3.771764586923218*^9}},ExpressionUUID->"be645cf4-80f8-46a4-a854-\
7e7f3cf3adaa"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   SubscriptBox["\[Delta]", "2"], "[", 
   SubscriptBox["r", "s"], "]"}], "\[Equal]", 
  RowBox[{"0.073867", 
   SubsuperscriptBox["r", "s", 
    RowBox[{"3", "/", "2"}]], "\t", 
   RowBox[{
    SubscriptBox["d", "0"], "[", 
    RowBox[{
     SubscriptBox["r", "s"], ",", "\[Zeta]"}], "]"}]}], "\[Equal]", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"0.70605", "+", 
     RowBox[{"0.12927", 
      SuperscriptBox["\[Zeta]", "2"]}]}], ")"}], 
   SubscriptBox["r", "s"]}]}]], "Input",
 CellChangeTimes->{{3.771764590484206*^9, 
  3.771764650992613*^9}},ExpressionUUID->"2bcddbd2-7cce-438d-a1b4-\
20e94d4eb82b"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1300, 1391},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 306, 4, 98, "Title",ExpressionUUID->"cc8f10f6-dee0-4be6-9b95-a2b645508bc1"],
Cell[889, 28, 762, 22, 32, "Input",ExpressionUUID->"f09b5b1d-fcfa-4277-b055-8d013d150e51"],
Cell[1654, 52, 453, 14, 51, "Input",ExpressionUUID->"2adb86e5-1613-4560-afb2-c874c93162cc"],
Cell[2110, 68, 471, 14, 50, "Input",ExpressionUUID->"e471ce7c-fc48-482e-8bea-3e9e0179671c"],
Cell[2584, 84, 1207, 40, 52, "Input",ExpressionUUID->"d1a8b6ea-55ae-4a5e-8673-b9734def02ca"],
Cell[3794, 126, 1033, 35, 52, "Input",ExpressionUUID->"8257931c-d5a2-4fa5-bfcf-242ed662fad5"],
Cell[CellGroupData[{
Cell[4852, 165, 249, 4, 67, "Section",ExpressionUUID->"da7db347-895a-4445-adf0-185e3e7ae4bc"],
Cell[5104, 171, 2180, 58, 69, "Input",ExpressionUUID->"2424c22e-2944-42fb-a926-6931d6b950db"],
Cell[7287, 231, 631, 17, 50, "Input",ExpressionUUID->"5e7071ef-4b10-4195-9ab5-72086ae9360d"],
Cell[7921, 250, 1501, 41, 53, "Input",ExpressionUUID->"431fbccd-665c-479c-825d-74e37905a953"],
Cell[9425, 293, 2439, 73, 115, "Input",ExpressionUUID->"68450462-8882-467d-85de-668cca5b67df"],
Cell[11867, 368, 2917, 88, 210, "Input",ExpressionUUID->"76293f49-0f9c-448d-9fa4-7a7821619364"],
Cell[14787, 458, 1627, 46, 53, "Input",ExpressionUUID->"7b438198-7fff-4469-8902-d63008ab1edd"],
Cell[16417, 506, 1439, 41, 53, "Input",ExpressionUUID->"bb62c6dd-19f1-49b0-93ed-7a72089471a5"],
Cell[17859, 549, 1518, 44, 112, "Input",ExpressionUUID->"1b3202c0-1ab4-429c-8dba-a60935c07597"],
Cell[19380, 595, 883, 23, 52, "Input",ExpressionUUID->"30aeeec3-f7a5-4138-bd8c-548696d2f752"],
Cell[CellGroupData[{
Cell[20288, 622, 563, 12, 30, "Input",ExpressionUUID->"408e19ad-23db-4b1a-8d84-8acca9835fae"],
Cell[20854, 636, 386, 6, 34, "Output",ExpressionUUID->"0b40f4ad-6233-4b3d-8a80-07aca06e731a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21277, 647, 282, 6, 30, "Input",ExpressionUUID->"7d031781-f5af-40e6-b079-1bc0fdb02f2b"],
Cell[21562, 655, 3550, 64, 531, "Output",ExpressionUUID->"e8c89351-e3be-455d-b993-40a6da00af8b"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[25161, 725, 156, 3, 67, "Section",ExpressionUUID->"1774903d-cd4a-427c-913b-6d03354a574f"],
Cell[25320, 730, 968, 30, 57, "Input",ExpressionUUID->"a08e663f-6de1-46b7-a344-1b47fcd2e96e"],
Cell[26291, 762, 975, 31, 33, "Input",ExpressionUUID->"760031e2-bb91-4ce7-8b6e-8566604d6f6c"],
Cell[27269, 795, 1090, 39, 65, "Input",ExpressionUUID->"be645cf4-80f8-46a4-a854-7e7f3cf3adaa"],
Cell[28362, 836, 640, 20, 33, "Input",ExpressionUUID->"2bcddbd2-7cce-438d-a1b4-20e94d4eb82b"]
}, Open  ]]
}, Open  ]]
}
]
*)