data=Dataset[With[{
labels=ToExpression/@Rest[StringSplit[FileBaseName[#],"_"]],
data=Import[#]
},
<|
"n"->labels[[1]],"β"->labels[[2]],
"E"->data[[1,2]]/data[[1,1]],"E2"->data[[1,3]]/data[[1,1]],
"C"->data[[2]]/data[[1,1]],
"M"->Partition[data[[3]],labels[[1]]]|>]&/@FileNames["~/doc/research/hadamard/code/monte_carlo/kauzmann_*.dat"]];

With[{d=KeySort[data[GroupBy[#n&],All,{#β,#E/((#n)^2)}&]]},
ListPlot[Sort/@Values[d],Joined->True,PlotLegends->Normal@Keys[d]]
]

With[{d=KeySort[data[GroupBy[#n&],All,{#β,((#β)^2)(#E2-((#E)^2))/((#n)^2)}&]]},
ListPlot[Sort/@Values[d],Joined->True,PlotLegends->Normal@Keys[d]]
]

Show[ListLogLogPlot[data[Select[#n==28&],Thread[{Range[Length[#C]]-1,(#C-((#E)^2))/(#C[[1]]-((#E)^2))}]&],Joined->True,PlotRange->{{0,50000},{0,1.1}}],LogPlot[0.4Exp[-t/240],{t,0,5000},PlotStyle->Black]]

100 (28^2)

Normal@data[SelectFirst[#β==9.5&],#M.Transpose[#M]&]//MatrixForm

V

Range[8]/(8.+1)

D[Sqrt[x],x]