kB=1.3806503`*^-23/(1.62*10^-19);
elementsName={"lithium","sodium","potassium","rubidium","cesium","copper","silver","gold","beryllium","magnesium","calcium","strontium","barium","niobium","iron","manganese","zinc","cadmium","mercury","aluminum","gallium","indium","thallium","tin","lead","bismuth","antimony"};

w[Elements_]:=Switch[Elements,
"lithium",4.74,"sodium",3.24,"potassium",2.12,"rubidium",1.85,"cesium",1.59,"copper",7.,"silver",5.49,"gold",5.53,"beryllium",14.3,"magnesium",7.08,"calcium",4.69,"strontium",3.93,"barium",3.64,"niobium",5.32,"iron",11.1,"manganese",10.9,"zinc",9.47,"cadmium",7.47,"mercury",7.13,"aluminum",11.7,"gallium",10.4,"indium",8.63,"thallium",8.15,"tin",10.2,"lead",9.47,"bismuth",9.9,"antimony",10.9]

f[e_,  T_]:= 1/(Exp[(e-4.74)/(kB *T)]+1)

Manipulate[(Plot[{1/(Exp[(ϵ-#)/(kB *80)]+1),1/(Exp[(ϵ-#)/(kB temperature)]+1)},{ϵ,0,2*#},PlotRange->{{0,2#},{-0.01,1.05}},Frame->True,Axes->False,PlotStyle->{{RGBColor[.48,.11,.56],Thick},{RGBColor[.88,.63,.23],Dashed,Thick}},FrameLabel->{{"Fermi-Dirac function",""},{Row[{"energy (", eV,")"}],Row[{" Fermi energy of ",elements== # eV}] }},Filling->{1->{{2},{Red,Blue}}},Frame->True,Axes->False,ImageSize->{700,450}, Epilog->{{Red, Text["Fermi Energy",{#-2,0.1},{-1,-1}]},{Black, Arrowheads[{0.04}],Arrow[{{#-1,0.1},{#,-0.01}}]}}])&[w[elements]],{{elements,"lithium"},elementsName,ControlType->PopupMenu},
{{temperature,.01,"temperature (°K)"}, 0.01,4.0*10^5,Appearance -> "Labeled"},SaveDefinitions->True]

Plot[{f[e,0.01],f[e,1000], f[e,10000], f[e,30000],f[e,100000]},{e, 0, 2*4.47}, PlotLabels->Placed[{"0K","10^3K", "10^4K", "3*10^4"K, "10^5K"},Above],Frame->True,FrameLabel->{{ "Fermi-Dirac Function",""} ,{"Energy eV","Fermi-Dirac Distribution of Lithium"}},Epilog->{{Red, Text["Fermi Energy 4.47eV",{1.74,0.1},{-1,-1}]},{Blue, Arrowheads[{0.04}],Arrow[{{3.74,0.1},{4.74,0}}]}},ImageSize->{450,350}]

Plot[{f[e,80],f[e,10000]},{e, 0, 2*4.47}, PlotLabels->Placed[{"T=10^4K"},Above],Filling->{2->{Bottom, LightBrown},1->{{2},{Red,Blue}}}, PlotRange->{{0,2*4.47},{0,1}},Frame->True,FrameLabel->{{ "Fermi-Dirac Function",""} ,{"Energy eV","Fermi-Dirac Distribution of Lithium"}},Epilog->{{Red, Text["Fermi Energy 4.47eV",{2,0.1},{-1,-1}]},{Black, Arrowheads[{0.04}],Arrow[{{3.74,0.1},{4.74,0}}]}, {Black, Text["U_{E<E_F}",{2,0.5},{-1,-1}]},{Black, Text["U_{E>E_F}",{5,0.07},{-1,-1}]}},ImageSize->{450,350}]