(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 311850, 7553] NotebookOptionsPosition[ 309832, 7479] NotebookOutlinePosition[ 310224, 7496] CellTagsIndexPosition[ 310181, 7493] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Define spaces for nuclear indices, define indices, and define operators \ themselves\ \>", "Section"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"nocc", " ", "=", " ", RowBox[{"Append", "[", RowBox[{"occ", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"nvirt", "=", RowBox[{"Append", "[", RowBox[{"virt", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"nothervirt", "=", RowBox[{"Append", "[", RowBox[{"othervirt", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"nallvirt", "=", RowBox[{"Append", "[", RowBox[{"allvirt", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"nany", " ", "=", " ", RowBox[{"Append", "[", RowBox[{"any", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"nallany", " ", "=", " ", RowBox[{"Append", "[", RowBox[{"allany", ",", RowBox[{"particleType", "[", "nucleus", "]"}]}], "]"}]}], ";"}], "\n"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"i", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"j", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"k", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"l", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"m", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"n", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"o", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"w", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "occ"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"II", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"JJ", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"KK", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"LL", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"MM", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"NN", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"OO", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"WW", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nocc"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"a", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"b", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"c", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"d", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"e", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"f", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "virt"}], "]"}]}], ";"}], "\n"}], "\n", RowBox[{ RowBox[{"AA", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"BB", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"CC", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"DD", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"EE", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"FF", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nvirt"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"p", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "any"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"q", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "any"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"r", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "any"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"s", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "any"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[Kappa]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[Kappa]\>\"", ",", "allany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[Lambda]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[Lambda]\>\"", ",", "allany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[Mu]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[Mu]\>\"", ",", "allany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[Nu]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[Nu]\>\"", ",", "allany"}], "]"}]}], ";"}], "\n"}], "\n", RowBox[{ RowBox[{"PP", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"QQ", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"RR", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"SS", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\\"", ",", "nany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[CapitalKappa]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[CapitalKappa]\>\"", ",", "nallany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[CapitalLambda]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[CapitalLambda]\>\"", ",", "nallany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"\[CapitalMu]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[CapitalMu]\>\"", ",", "nallany"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"\[CapitalNu]", "=", RowBox[{"createParticleIndex", "[", RowBox[{"\"\<\[CapitalNu]\>\"", ",", "nallany"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"Aia", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "a", "}"}], ",", RowBox[{"{", "i", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Abj", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "j", "}"}], ",", RowBox[{"{", "b", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Aijab", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Acdmn", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"m", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"c", ",", "d"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Aijkabc", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "c"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "j", ",", "k"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"Adefmno", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"m", ",", "n", ",", "o"}], "}"}], ",", RowBox[{"{", RowBox[{"d", ",", "e", ",", "f"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"AIA", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "AA", "}"}], ",", RowBox[{"{", "II", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"ABJ", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "JJ", "}"}], ",", RowBox[{"{", "BB", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"AIJAB", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"AA", ",", "BB"}], "}"}], ",", RowBox[{"{", RowBox[{"II", ",", "JJ"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"ACDMN", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"MM", ",", "NN"}], "}"}], ",", RowBox[{"{", RowBox[{"CC", ",", "DD"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"AIJKABC", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"AA", ",", "BB", ",", "CC"}], "}"}], ",", RowBox[{"{", RowBox[{"II", ",", "JJ", ",", "KK"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"ADEFMNO", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"MM", ",", "NN", ",", "OO"}], "}"}], ",", RowBox[{"{", RowBox[{"DD", ",", "EE", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"AiJaB", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "BB"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "JJ"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"AcDmN", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"m", ",", "NN"}], "}"}], ",", RowBox[{"{", RowBox[{"c", ",", "DD"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"AijKabC", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", "CC"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "j", ",", "KK"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"AdeFmnO", "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"m", ",", "n", ",", "OO"}], "}"}], ",", RowBox[{"{", RowBox[{"d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ SubscriptBox["G", "\[Mu]\[CapitalNu]\[Kappa]\[CapitalLambda]"], "=", RowBox[{"createSQM", "[", RowBox[{"\"\\"", ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "\[CapitalNu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Kappa]", ",", "\[CapitalLambda]"}], "}"}], ",", "nonsymm"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ SubscriptBox["A", "\[Mu]\[CapitalNu]\[Kappa]\[CapitalLambda]"], "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"\[Mu]", ",", "\[CapitalNu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Kappa]", ",", "\[CapitalLambda]"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"H2en", "=", RowBox[{ SubscriptBox["G", "\[Mu]\[CapitalNu]\[Kappa]\[CapitalLambda]"], SubscriptBox["A", "\[Mu]\[CapitalNu]\[Kappa]\[CapitalLambda]"]}]}], ";"}], "\n"}], "\n", RowBox[{ RowBox[{ SubscriptBox["G", "\[Mu]\[Nu]\[Kappa]\[Lambda]"], "=", RowBox[{"createSQM", "[", RowBox[{"\"\<\!\(\*OverscriptBox[\(g\), \(_\)]\)\>\"", ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "\[Nu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Kappa]", ",", "\[Lambda]"}], "}"}], ",", "antisymm"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ SubscriptBox["A", "\[Mu]\[Nu]\[Kappa]\[Lambda]"], "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"\[Mu]", ",", "\[Nu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Kappa]", ",", "\[Lambda]"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"H2ee", "=", RowBox[{ RowBox[{"1", "/", "4"}], " ", SubscriptBox["G", "\[Mu]\[Nu]\[Kappa]\[Lambda]"], SubscriptBox["A", "\[Mu]\[Nu]\[Kappa]\[Lambda]"]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ SubscriptBox["G", "\[CapitalMu]\[CapitalNu]\[CapitalKappa]\[CapitalLambda]"], "=", RowBox[{"createSQM", "[", RowBox[{"\"\<\!\(\*OverscriptBox[\(g\), \(_\)]\)\>\"", ",", RowBox[{"{", RowBox[{"\[CapitalMu]", ",", "\[CapitalNu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[CapitalKappa]", ",", "\[CapitalLambda]"}], "}"}], ",", "antisymm"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ SubscriptBox["A", "\[CapitalMu]\[CapitalNu]\[CapitalKappa]\[CapitalLambda]"], "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", RowBox[{"\[CapitalMu]", ",", "\[CapitalNu]"}], "}"}], ",", RowBox[{"{", RowBox[{"\[CapitalKappa]", ",", "\[CapitalLambda]"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"H2nn", "=", RowBox[{ RowBox[{"1", "/", "4"}], " ", SubscriptBox["G", "\[CapitalMu]\[CapitalNu]\[CapitalKappa]\[CapitalLambda]"], SubscriptBox["A", "\[CapitalMu]\[CapitalNu]\[CapitalKappa]\[CapitalLambda]"]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ SubscriptBox["T", "\[Mu]\[Nu]"], "=", RowBox[{"createSQM", "[", RowBox[{"\"\\"", ",", RowBox[{"{", "\[Mu]", "}"}], ",", RowBox[{"{", "\[Nu]", "}"}], ",", "antisymm"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ SubscriptBox["A", "\[Mu]\[Nu]"], "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "\[Mu]", "}"}], ",", RowBox[{"{", "\[Nu]", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"H1e", "=", " ", RowBox[{ SubscriptBox["T", "\[Mu]\[Nu]"], "*", SubscriptBox["A", "\[Mu]\[Nu]"]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{ SubscriptBox["T", "\[CapitalMu]\[CapitalNu]"], "=", RowBox[{"createSQM", "[", RowBox[{"\"\\"", ",", RowBox[{"{", "\[CapitalMu]", "}"}], ",", RowBox[{"{", "\[CapitalNu]", "}"}], ",", "antisymm"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ SubscriptBox["A", "\[CapitalMu]\[CapitalNu]"], "=", RowBox[{"createSQS", "[", RowBox[{ RowBox[{"{", "\[CapitalMu]", "}"}], ",", RowBox[{"{", "\[CapitalNu]", "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"H1n", "=", " ", RowBox[{ SubscriptBox["T", "\[CapitalMu]\[CapitalNu]"], "*", SubscriptBox["A", "\[CapitalMu]\[CapitalNu]"]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"H1", "=", RowBox[{"H1e", "+", "H1n"}]}], ";"}], "\n", RowBox[{ RowBox[{"H2", " ", "=", " ", RowBox[{"H2ee", "+", "H2nn", "+", "H2en"}]}], ";"}], "\n", RowBox[{ RowBox[{"H1", "+", "H2"}], "//", "TraditionalForm"}]}], "Input"], Cell[BoxData[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"\[Mu]\"\>", "\<\"\[Nu]\"\>"], Subsuperscript["T", "\[Mu]", "\[Nu]"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[ATilde]\"\>", "\<\"\[Nu]\"\>", \ "\<\"\[Mu]\"\>"], Subsuperscript["\[ATilde]", "\[Nu]", "\[Mu]"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"\[CapitalMu]\"\>", \ "\<\"\[CapitalNu]\"\>"], Subsuperscript["T", "\[CapitalMu]", "\[CapitalNu]"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[ATilde]\"\>", "\<\"\[CapitalNu]\"\>", "\<\"\ \[CapitalMu]\"\>"], Subsuperscript["\[ATilde]", "\[CapitalNu]", "\[CapitalMu]"], Editable->False], TraditionalForm]}], "+", RowBox[{ FractionBox["1", "4"], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"\[Mu]\[Nu]\"\>", "\<\"\[Kappa]\[Lambda]\"\>"], Subsuperscript[ "\!\(\*OverscriptBox[\(g\), \(_\)]\)", "\[Mu]\[Nu]", "\[Kappa]\[Lambda]"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[ATilde]\"\>", "\<\"\[Kappa]\[Lambda]\"\>", "\<\ \"\[Mu]\[Nu]\"\>"], Subsuperscript["\[ATilde]", "\[Kappa]\[Lambda]", "\[Mu]\[Nu]"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"\[Mu]\[CapitalNu]\"\>", "\<\"\ \[Kappa]\[CapitalLambda]\"\>"], Subsuperscript["g", "\[Mu]\[CapitalNu]", "\[Kappa]\[CapitalLambda]"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[ATilde]\"\>", \ "\<\"\[Kappa]\[CapitalLambda]\"\>", "\<\"\[Mu]\[CapitalNu]\"\>"], Subsuperscript[ "\[ATilde]", "\[Kappa]\[CapitalLambda]", "\[Mu]\[CapitalNu]"], Editable->False], TraditionalForm]}], "+", RowBox[{ FractionBox["1", "4"], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"\[CapitalMu]\[CapitalNu]\"\>", "\<\"\[CapitalKappa]\ \[CapitalLambda]\"\>"], Subsuperscript[ "\!\(\*OverscriptBox[\(g\), \(_\)]\)", "\[CapitalMu]\[CapitalNu]", "\[CapitalKappa]\[CapitalLambda]"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[ATilde]\"\>", "\<\"\[CapitalKappa]\ \[CapitalLambda]\"\>", "\<\"\[CapitalMu]\[CapitalNu]\"\>"], Subsuperscript[ "\[ATilde]", "\[CapitalKappa]\[CapitalLambda]", "\[CapitalMu]\[CapitalNu]"], Editable->False], TraditionalForm]}]}], TraditionalForm]], "Output", CellChangeTimes->{3.405108659123886*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Compute", " ", SubscriptBox[ RowBox[{ SubscriptBox["CISD", "e"], "[", "T", "]"}], "n"], " ", "matrix", " ", "elements"}]], "Section"], Cell["\<\ Comments: 1) only matrix elements unique to this method (including 2+1 mixed \ triples,i.e. double electronic and single nuclear excitations) are shown 2) all matrix elements between the reference and the 2+1 mixed triples are \ zero\ \>", "Text"], Cell[BoxData[{ RowBox[{ RowBox[{"dim", "=", RowBox[{"deltaIndex", "[", RowBox[{"i", ",", "m"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"djn", "=", RowBox[{"deltaIndex", "[", RowBox[{"j", ",", "n"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"din", "=", RowBox[{"deltaIndex", "[", RowBox[{"i", ",", "n"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"djm", "=", RowBox[{"deltaIndex", "[", RowBox[{"j", ",", "m"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dda", "=", RowBox[{"deltaIndex", "[", RowBox[{"d", ",", "a"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"deb", "=", RowBox[{"deltaIndex", "[", RowBox[{"e", ",", "b"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dea", "=", RowBox[{"deltaIndex", "[", RowBox[{"e", ",", "a"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ddb", "=", RowBox[{"deltaIndex", "[", RowBox[{"d", ",", "b"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dFB", "=", RowBox[{"deltaIndex", "[", RowBox[{"FF", ",", "BB"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dJO", "=", RowBox[{"deltaIndex", "[", RowBox[{"JJ", ",", "OO"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dFC", "=", RowBox[{"deltaIndex", "[", RowBox[{"FF", ",", "CC"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dKO", "=", RowBox[{"deltaIndex", "[", RowBox[{"KK", ",", "OO"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "Input"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "Aia"}], ",", RowBox[{"{", RowBox[{ "i", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-Esingles H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox["0", TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-Esingles H1:\n", TraditionalForm[0]], Editable->False]], "Print", CellChangeTimes->{3.4051086668505163`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H2", "**", "Aia"}], ",", RowBox[{"{", RowBox[{ "i", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", "dim"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", "din"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-Esingles H2:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "dF", "mO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "eF", "mO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "eF", "nO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "dF", "nO"], Editable->False], TraditionalForm]}]}], ")"}]}]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-Esingles H2:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]])]], Editable->False]], "Print", CellChangeTimes->{3.405108669235981*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "AIA"}], ",", RowBox[{"{", RowBox[{ "II", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "AA", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-Nsingles H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox["0", TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-Nsingles H1:\n", TraditionalForm[0]], Editable->False]], "Print", CellChangeTimes->{3.4051086715755863`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H2", "**", "AIA"}], ",", RowBox[{"{", RowBox[{ "II", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "AA", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-Nsingles H2:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"A\"\>"], Subsuperscript["\[Delta]", "F", "A"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"I\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "I", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"de\"\>", "\<\"mn\"\>"], Subsuperscript["\!\(\*OverscriptBox[\(g\), \(_\)]\)", "de", "mn"], Editable->False], TraditionalForm]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-Nsingles H2:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["A", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["I", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] $CellContext`SQM[ $CellContext`OHead["\!\(\*OverscriptBox[\(g\), \(_\)]\)", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]]]], Editable->False]], "Print", CellChangeTimes->{3.4051086736952477`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "Aijab"}], ",", RowBox[{"{", RowBox[{ "i", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dim", "*", "djn"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"din", "*", "djm"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-EEdoubles H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"O\"\>"], Subsuperscript["T", "F", "O"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"O\"\>"], Subsuperscript["T", "F", "O"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"O\"\>"], Subsuperscript["T", "F", "O"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"O\"\>"], Subsuperscript["T", "F", "O"], Editable->False], TraditionalForm]}]}], ")"}]}]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-EEdoubles H1:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]])]], Editable->False]], "Print", CellChangeTimes->{3.4051086759550257`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H2", "**", "Aijab"}], ",", RowBox[{"{", RowBox[{ "i", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dim", "*", "djn"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"din", "*", "djm"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dda", "*", "deb"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dea", "*", "ddb"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-EEdoubles H2:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm]}], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"aO\"\>"], Subsuperscript["g", "dF", "aO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"bO\"\>"], Subsuperscript["g", "dF", "bO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"aO\"\>"], Subsuperscript["g", "eF", "aO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"bO\"\>"], Subsuperscript["g", "eF", "bO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"aO\"\>"], Subsuperscript["g", "dF", "aO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"bO\"\>"], Subsuperscript["g", "dF", "bO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"aO\"\>"], Subsuperscript["g", "eF", "aO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"bO\"\>"], Subsuperscript["g", "eF", "bO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm]}], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"iF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "iF", "mO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"iF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "iF", "nO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"jF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "jF", "mO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"jF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "jF", "nO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"iF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "iF", "mO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"iF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "iF", "nO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"jF\"\>", "\<\"mO\"\>"], Subsuperscript["g", "jF", "mO"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"jF\"\>", "\<\"nO\"\>"], Subsuperscript["g", "jF", "nO"], Editable->False], TraditionalForm]}]}], ")"}]}]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-EEdoubles H2:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["b", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["j", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]])]], Editable->False]], "Print", CellChangeTimes->{3.405108680050715*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "AiJaB"}], ",", RowBox[{"{", RowBox[{ "i", ",", "JJ", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "BB", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-NEdoubles H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"d\"\>", "\<\"m\"\>"], Subsuperscript["T", "d", "m"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"d\"\>", "\<\"n\"\>"], Subsuperscript["T", "d", "n"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"e\"\>", "\<\"m\"\>"], Subsuperscript["T", "e", "m"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"e\"\>", "\<\"n\"\>"], Subsuperscript["T", "e", "n"], Editable->False], TraditionalForm]}]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-NEdoubles H1:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] - $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] - $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] + $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] $CellContext`SQM[ $CellContext`OHead["T", $CellContext`indexSymm[-1]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]]]], Editable->False]], "Print", CellChangeTimes->{3.405108688669647*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H2", "**", "AiJaB"}], ",", RowBox[{"{", RowBox[{ "i", ",", "JJ", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "BB", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dJO", "*", "dFB"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dJO", "*", "dea"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dJO", "*", "dda"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dFB", "*", "dea"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dFB", "*", "dda"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-NEdoubles H2:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"mB\"\>"], Subsuperscript["g", "dF", "mB"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dF\"\>", "\<\"nB\"\>"], Subsuperscript["g", "dF", "nB"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dJ\"\>", "\<\"nO\"\>"], Subsuperscript["g", "dJ", "nO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"dJ\"\>", "\<\"mO\"\>"], Subsuperscript["g", "dJ", "mO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"nB\"\>"], Subsuperscript["g", "eF", "nB"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eF\"\>", "\<\"mB\"\>"], Subsuperscript["g", "eF", "mB"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eJ\"\>", "\<\"mO\"\>"], Subsuperscript["g", "eJ", "mO"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"g\"\>", "\<\"eJ\"\>", "\<\"nO\"\>"], Subsuperscript["g", "eJ", "nO"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"B\"\>"], Subsuperscript["\[Delta]", "F", "B"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"J\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "J", "O"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm]}], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"de\"\>", "\<\"am\"\>"], Subsuperscript["\!\(\*OverscriptBox[\(g\), \(_\)]\)", "de", "am"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"de\"\>", "\<\"an\"\>"], Subsuperscript["\!\(\*OverscriptBox[\(g\), \(_\)]\)", "de", "an"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"di\"\>", "\<\"mn\"\>"], Subsuperscript["\!\(\*OverscriptBox[\(g\), \(_\)]\)", "di", "mn"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\\!\\(\\*OverscriptBox[\\(g\\), \ \\(_\\)]\\)\"\>", "\<\"ei\"\>", "\<\"mn\"\>"], Subsuperscript["\!\(\*OverscriptBox[\(g\), \(_\)]\)", "ei", "mn"], Editable->False], TraditionalForm]}]}], ")"}]}]}], TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-NEdoubles H2:\n", TraditionalForm[$CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ ($CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["g", $CellContext`indexSymm[0]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]], $CellContext`indexType[$CellContext`ket]]]) + \ $CellContext`deltaIndex[ $CellContext`particleIndex["F", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["B", $CellContext`particleSpace[$CellContext`virtual, $CellContext`particleType[$CellContext`nucleus]]]] \ $CellContext`deltaIndex[ $CellContext`particleIndex["J", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]], $CellContext`particleIndex["O", $CellContext`particleSpace[$CellContext`occupied, $CellContext`particleType[$CellContext`nucleus]]]] \ (-$CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["\!\(\*OverscriptBox[\(g\), \(_\)]\)", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied]]] \ $CellContext`SQM[ $CellContext`OHead["\!\(\*OverscriptBox[\(g\), \(_\)]\)", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] - \ $CellContext`deltaIndex[ $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] $CellContext`SQM[ $CellContext`OHead["\!\(\*OverscriptBox[\(g\), \(_\)]\)", $CellContext`indexSymm[-1]], $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]] + \ $CellContext`deltaIndex[ $CellContext`particleIndex["d", $CellContext`particleSpace[$CellContext`virtual]], $CellContext`particleIndex["a", $CellContext`particleSpace[$CellContext`virtual]]] \ $CellContext`SQM[ $CellContext`OHead["\!\(\*OverscriptBox[\(g\), \(_\)]\)", $CellContext`indexSymm[-1]], $CellContext`particleIndex["e", $CellContext`particleSpace[$CellContext`virtual], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["i", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`bra]], $CellContext`particleIndex["m", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]], $CellContext`particleIndex["n", $CellContext`particleSpace[$CellContext`occupied], $CellContext`indexType[$CellContext`ket]]])]], Editable->False]], "Print", CellChangeTimes->{3.405108691938303*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "AIJAB"}], ",", RowBox[{"{", RowBox[{ "II", ",", "JJ", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "AA", ",", "BB", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-NNdoubles H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox["0", TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-NNdoubles H1:\n", TraditionalForm[0]], Editable->False]], "Print", CellChangeTimes->{3.405108694214168*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H2", "**", "AIJAB"}], ",", RowBox[{"{", RowBox[{ "II", ",", "JJ", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "AA", ",", "BB", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-NNdoubles H2:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox["0", TraditionalForm], TraditionalForm, Editable->True]}], SequenceForm["NEEtriples-NNdoubles H2:\n", TraditionalForm[0]], Editable->False]], "Print", CellChangeTimes->{3.405108696327463*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"tmp", " ", "=", " ", RowBox[{"wick", "[", RowBox[{ RowBox[{"AdeFmnO", "**", "H1", "**", "AijKabC"}], ",", RowBox[{"{", RowBox[{ "i", ",", "j", ",", "KK", ",", "m", ",", "n", ",", "OO", ",", "a", ",", "b", ",", "CC", ",", "d", ",", "e", ",", "FF"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"factorizedKOdFC", "[", "f_", "]"}], ":=", RowBox[{"Collect", "[", RowBox[{"f", ",", RowBox[{"dKO", "*", "dFC"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"dim", "*", "djn"}], ",", "factorizedKOdFC"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tmp", "=", RowBox[{"Collect", "[", RowBox[{"tmp", ",", RowBox[{"din", "*", "djm"}], ",", "factorizedKOdFC"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"tmp", "//", "TraditionalForm"}]}], "]"}], ";"}]}], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"NEEtriples-NEEtriples H1:\\n\"\>", "\[InvisibleSpace]", TagBox[ FormBox[ RowBox[{ RowBox[{ RowBox[{"-", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm]}], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"K\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "K", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"C\"\>"], Subsuperscript["T", "F", "C"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"K\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "K", "O"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"F\"\>", "\<\"C\"\>"], Subsuperscript["T", "F", "C"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"C\"\>"], Subsuperscript["\[Delta]", "F", "C"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"K\"\>", "\<\"O\"\>"], Subsuperscript["\[Delta]", "K", "O"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", RowBox[{"(", RowBox[{ RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", \ "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["T", "d", "a"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", \ "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["T", "d", "b"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", \ "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["T", "e", "a"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", \ "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["T", "e", "b"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["T", "i", "m"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "j", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["T", "i", "m"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["T", "i", "n"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "j", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["T", "i", "n"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["T", "j", "m"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"n\"\>"], Subsuperscript["\[Delta]", "i", "n"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"j\"\>", "\<\"m\"\>"], Subsuperscript["T", "j", "m"], Editable->False], TraditionalForm]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["T", "j", "n"], Editable->False], TraditionalForm]}], "-", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "d", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "e", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"i\"\>", "\<\"m\"\>"], Subsuperscript["\[Delta]", "i", "m"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"T\"\>", "\<\"j\"\>", "\<\"n\"\>"], Subsuperscript["T", "j", "n"], Editable->False], TraditionalForm]}]}], ")"}]}], "+", RowBox[{ FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"d\"\>", "\<\"b\"\>"], Subsuperscript["\[Delta]", "d", "b"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"e\"\>", "\<\"a\"\>"], Subsuperscript["\[Delta]", "e", "a"], Editable->False], TraditionalForm], " ", FormBox[ InterpretationBox[ SubsuperscriptBox["\<\"\[Delta]\"\>", "\<\"F\"\>", "\<\"C\"\>"], Subsuperscript["\[Delta]", "F", "C"], Editable->False], TraditionalForm], " ",