c1=Sqrt[(1/2)^2+1/2*((1)^2-(1/2)^2)*(1+(2/Pi)*ArcTan[(x+(8*Tan[(Pi/((1)^2-(1/2)^2))*((3/4)^2-1/2*((1/2)^2+1^2))]))/(8)])];

kappa=N[D[c1,x]/.x->0, 1200];

nuRZeroPotential=1/(Exp[2*Pi*omega/kappa]-1);
nuLZeroPotential=1/(4*omega^2)*Abs[Exp[Pi*omega/kappa](0)/(2*I*omega+0)+(2*I*omega)]^2*1/(Exp[2*Pi*omega/kappa]-1)+0^2/(4*omega^2+0^2);
nvLZeroPotential=1/(4*omega^2)*Abs[0(2*I*omega-0)/(2*I*omega-0)-0*Exp[Pi*omega/kappa]]^2*1/(Exp[2*Pi*omega/kappa]-1);
LogLinearPlot[omega*nuRZeroPotential,{omega,1/1000000,2},PlotStyle->Automatic, ImageSize->Large, WorkingPrecision->30, PlotLabel->" No potential ", AxesLabel->{ Style["ω", 20],Style[Superscript["ω n","u,R"], 20] }]

LogLinearPlot[omega*nuLZeroPotential,{omega,1/1000000,2},PlotStyle->Automatic, ImageSize->Large, WorkingPrecision->30, PlotLabel->" No potential ", AxesLabel->{ ω,Superscript["ω n","u,L"] }]

LogLinearPlot[omega*nvLZeroPotential,{omega,1/1000000,2},PlotStyle->Automatic, ImageSize->Large, WorkingPrecision->30, PlotLabel->" No potential ", AxesLabel->{ ω,Superscript["ω n","v,L"] }]

2 Dirac Delta Function Potential Case

nuRGeneric=(Abs[TH]^2)/(Exp[2*Pi*omega/kappa]-1);
nuLGeneric=(Abs[Conjugate[alpha]+RH*beta*Exp[Pi*omega/kappa]]^2)/(Exp[2*Pi*omega/kappa]-1)+Abs[TH*beta]^2;
nvLGeneric=(Abs[RH*alpha+Conjugate[beta]*Exp[Pi*omega/kappa]]^2)/(Exp[2*Pi*omega/kappa]-1);

nuR2DeltaGeneric=nuRGeneric/.{TH->2*I*omega/(2*I*omega-VR),RH->VR/(2*I*omega-VR),alpha->(2*I*omega+VL)/(2*I*omega),beta->(-VL)/(2*I*omega)};
nuL2DeltaGeneric=Simplify[ComplexExpand[nuLGeneric/.{TH->2*I*omega/(2*I*omega-VR),RH->VR/(2*I*omega-VR),alpha->(2*I*omega+VL)/(2*I*omega),beta->(-VL)/(2*I*omega)}]];
nvL2DeltaGeneric=Simplify[ComplexExpand[nvLGeneric/.{TH->2*I*omega/(2*I*omega-VR),RH->VR/(2*I*omega-VR),alpha->(2*I*omega+VL)/(2*I*omega),beta->(-VL)/(2*I*omega)}]];
Simplify[nvL2DeltaGeneric-ComplexExpand[1/(4*omega^2)*Abs[VR(2*I*omega+VL)/(2*I*omega-VR)+VL*Exp[Pi*omega/kappa]]^2*1/(Exp[2*Pi*omega/kappa]-1)]]

Simplify[nuL2DeltaGeneric-ComplexExpand[1/(4*omega^2)*Abs[-Exp[Pi*omega/kappa](VR*VL)/(2*I*omega+VR)+(2*I*omega+VL)]^2*1/(Exp[2*Pi*omega/kappa]-1)+VL^2/(4*omega^2+VR^2)]]

nuR2DeltaGenericCheck=4*omega^2/(4*omega^2+VR^2)*1/(Exp[2*Pi*omega/kappa]-1);
nuL2DeltaGenericCheck=1/(4*omega^2)*Abs[Exp[Pi*omega/kappa](VR*VL)/(2*I*omega+VR)+(2*I*omega-VL)]^2*1/(Exp[2*Pi*omega/kappa]-1)+VL^2/(4*omega^2+VR^2);
nvL2DeltaGenericCheck=Simplify[ComplexExpand[1/(4*omega^2)*Abs[VR(2*I*omega-VL)/(2*I*omega-VR)-VL*Exp[Pi*omega/kappa]]^2*1/(Exp[2*Pi*omega/kappa]-1)]];
nvLDelta          =     1/(4*omega^2)*Abs[VR(2*I*omega-VL)/(2*I*omega-VR)-VL*Exp[Pi*omega/kappa]]^2*1/(Exp[2*Pi*omega/kappa]-1);

nuRmassive=nuR2DeltaGeneric/.{VR-> 0,VL-> 0};
nuLmassive=nuL2DeltaGeneric/.{VR-> 0,VL->0};
nvLmassive=nvL2DeltaGeneric/.{VR-> 0,VL->0};
nuLmassiveCheck=nuL2DeltaGenericCheck/.{VR-> 10*kappa,VL->100*kappa};
nvLmassive=nvL2DeltaGeneric/.{VR-> 10*kappa,VL->-100*kappa};
nvLmassiveCheck=nvL2DeltaGenericCheck/.{VR-> 10*kappa,VL->100*kappa};
Simplify[nuLmassive-nuLmassiveCheck];
Simplify[nvLmassive-nvLmassiveCheck]

uRuRMassive=LogLinearPlot[omega*nuRmassive,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] }]
uLuLMassive=LogLinearPlot[omega*nvLmassive,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,L"] }]
vLvLMassive=LogLinearPlot[omega*nvLmassive,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","v,L"] }]

uRuRMassiveRatio=LogLinearPlot[nuRmassive/nuRZeroPotential,{omega,2/100,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case Ratio" , PlotRange->All, AxesLabel->{ ω,Ratio  Superscript[" n","u,R"] }]
uLuLMassiveRatio=LogLinearPlot[nuLmassive/nuLZeroPotential,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case Ratio" , PlotRange->All, AxesLabel->{ ω,Ratio Superscript["ω n","u,L"] }]

uRuRMassiveCompare=LogLinearPlot[{nuRmassive,nuRZeroPotential},{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" 2 Dirac Delta Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] } ]

Massive Case 2

nuRmassive2=nuR2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
nuLmassive2=nuL2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
nvLmassive2=nvL2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
uRuRMassive2=LogLinearPlot[omega*nuRmassive2,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] }]
uLuLMassive2=LogLinearPlot[omega*nuLmassive2,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,L"] }]
vLvLMassive2=LogLinearPlot[omega*nvLmassive2,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","v,L"] }]

Massive Case 3

nuRmassive3=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};
nuLmassive3=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};
nvLmassive3=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};
uRuRMassive3=LogLinearPlot[omega*nuRmassive3,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] }]
uLuLMassive3=LogLinearPlot[omega*nuLmassive3,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,L"] }]
vLvLMassive3=LogLinearPlot[omega*nvLmassive3,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","v,L"] }]

nuRmassive4=nuR2DeltaGeneric/.{VR-> kappa*50,VL-> -20*kappa};
nuLmassive4=nuL2DeltaGeneric/.{VR-> kappa*50,VL-> -20*kappa};
nvLmassive4=nuL2DeltaGeneric/.{VR-> kappa*50,VL-> -20*kappa};
uRuRMassive4=LogLinearPlot[omega*nuRmassive4,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] }]
uLuLMassive4=LogLinearPlot[omega*nuLmassive4,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,L"] }]
vLvLMassive4=LogLinearPlot[omega*nvLmassive4,{omega,1/100000,1/10},PlotStyle->{Orange, Thick}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","v,L"] }]

Massive Comparison

nuLmassiveMax=(omega*nuLmassive)/.omega-> 10^-6;
nuLmassive2Max=(omega*nuLmassive2)/.omega-> 10^-6;
nuLmassive3Max=(omega*nuLmassive3)/.omega-> 10^-6;
nvLmassiveMax=(omega*nvLmassive)/.omega-> 10^-6;
nvLmassive2Max=(omega*nvLmassive2)/.omega-> 10^-6;
nvLmassive3Max=(omega*nvLmassive3)/.omega-> 10^-6;
nuRmassiveMax=FindMaximum[omega*nuRmassive, {omega,10^-6,9/200}];
nuRmassive2Max=FindMaximum[omega*nuRmassive2, {omega,10^-6,9/200}];
nuRmassive3Max=FindMaximum[omega*nuRmassive3, {omega,10^-6,9/200}];

nuLZeroPotentialMax2=(omega*nuLZeroPotential)/.omega-> 10^-6;
nvLmassiveMax=(omega*nvLmassive)/.omega-> 10^-6;
nuRmassiveMax=FindMaximum[omega*nuRmassive, {omega,10^-6,9/200}];

uRuRMassiveCompare123=LogLinearPlot[{(omega*nuRmassive)/nuRmassiveMax[[1]],(omega*nuRmassive2)/nuRmassive2Max[[1]],(omega*nuRmassive3)/nuRmassive3Max[[1]]},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Orange, Thick}  },Ticks->{{10^-6, 10^-4, 10^-2, 1}, Automatic}, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{Style[ ω, 20],Style[Superscript["C*ω n","u,r"], 20] },TicksStyle->Directive["Label", 16]]

uLuLMassiveCompare123=LogLinearPlot[{((omega*nuLZeroPotential)/nuLZeroPotentialMax2)^2,((omega*nuLmassive3)/nuLmassive3Max)^2},{omega,10^-6,1},PlotStyle->{{Orange, Thick, Dashed} ,{Blue, Thick, Dotted} },FrameLabel->{Style[ ω/κ, 14],Style[ HoldForm[ Superscript["|C ω (Subscript[n, up]^int)|","2"]], 16] },FrameTicksStyle->Directive["Label", 12],FrameTicks->{{Table[{k,k },{k,0,8}], None},{ Table[{10^(2*k),Superscript[10,2*k]},{k,-8,0}] ,None}},Frame->True, ImageSize->Full]

vLvLMassiveCompare123=LogLinearPlot[{0,((omega*nvLmassive3)/nvLmassive3Max)^2},{omega,10^-6,1},PlotStyle->{{Orange, Dashed},{Blue, Thick, Dotted} },FrameLabel->{Style[ ω/κ, 14],Style[ HoldForm[ Superscript["|C ω (Subscript[n, ds]^int)|","2"]], 16] },FrameTicksStyle->Directive["Label", 12],FrameTicks->{{Table[{k,k },{k,0,8}], None},{ Table[{10^(2*k),Superscript[10,2*k]},{k,-8,0}] ,None}},Frame->True, PlotRange->All, ImageSize->Full]

Export["/Users/richarddudley/Dropbox/Correlator/PlotProducers/uRuR2DeltaCompare123.pdf",uRuRMassiveCompare123]
Export["/Users/richarddudley/Dropbox/Correlator/PlotProducers/uLuL2DeltaCompare123.pdf",uLuLMassiveCompare123]
Export["/Users/richarddudley/Dropbox/Correlator/PlotProducers/vLvL2DeltaCompare123.pdf",vLvLMassiveCompare123]

nuRmassive20=nuR2DeltaGeneric/.{VR-> -kappa,VL-> 100*kappa};
nuLmassive20=nuL2DeltaGeneric/.{VR-> -kappa,VL-> 100*kappa};
nvLmassive20=nvL2DeltaGeneric/.{VR-> -kappa,VL-> 100*kappa};
uRuRMassive20=LogLinearPlot[{omega*nuRmassive20,omega*nuRmassive2},{omega,1/100000,1/10}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] },PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Orange, Thick} }]
uLuLMassive20=LogLogPlot[{omega*nuLmassive20,omega*nuLZeroPotential},{omega,1/100000,1/10}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,L"] },PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Orange, Thick} }]
vLvLMassive20=LogLinearPlot[{omega*nvLmassive20,omega*nvLmassive2},{omega,1/100000,1/10}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","v,L"] },PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Orange, Thick} }]

uRuRMassive20=LogLinearPlot[{omega*nuRmassive20,omega*nuLZeroPotential},{omega,1/100000,1/10}, ImageSize-> Full,  PlotLabel->" Massive Case" , PlotRange->All, AxesLabel->{ ω,Superscript["ω n","u,R"] },PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Orange, Thick} }]

nuLZeroPotentialMax=(omega*nuLZeroPotential)/.omega-> 10^-1;
nuLmassiveMax2=(omega*nuLmassive)/.omega-> 10^-1;

uLuLMassiveCompareZerotoFull=LogLinearPlot[{omega*(nuLmassive)/nuLmassiveMax2,(omega*nuLZeroPotential)/nuLZeroPotentialMax},{omega,10^-1,10},PlotStyle->{{Blue, Thick, Dashed},{Red, Thick }  }, ImageSize-> Full,AxesLabel->{ Style["ω", 20],Style[Superscript["C*ω n","u,l"], 20] }, PlotRange->All, WorkingPrecision->100,TicksStyle->Directive["Label", 16] ]

Export["D:\\Dropbox\\ParticleProductionDraft\\Figures\\uLuLDiracDeltaCompareZero.pdf",uLuLMassiveCompareZerotoFull]

nuRmassive=nuR2DeltaGeneric/.{VR-> 10*kappa,VL-> -100*kappa};
nuLmassive=nuL2DeltaGeneric/.{VR-> 10*kappa,VL->-100*kappa};
nvLmassive=nvL2DeltaGeneric/.{VR-> 10*kappa,VL->-100*kappa};
nuRmassive2v2=nuR2DeltaGeneric/.{VR-> 2*kappa,VL-> -100*kappa};
nuLmassive2v2=nuL2DeltaGeneric/.{VR->2*kappa,VL-> -100*kappa};
nvLmassive2v2=nvL2DeltaGeneric/.{VR-> 2*kappa,VL-> -100*kappa};
nuRmassive2=nuR2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
nuLmassive2=nuL2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
nvLmassive2=nvL2DeltaGeneric/.{VR-> kappa,VL-> -100*kappa};
nuRmassive3=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};
nuLmassive3=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};
nvLmassive3=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -100*kappa};

nuRmassive3v2=nuR2DeltaGeneric/.{VR-> kappa*1/3,VL-> -100*kappa};
nuLmassive3v2=nuL2DeltaGeneric/.{VR-> kappa*1/3,VL-> -100*kappa};
nvLmassive3v2=nuL2DeltaGeneric/.{VR-> kappa*1/3,VL-> -100*kappa};

nuRmassive3v3=nuR2DeltaGeneric/.{VR-> kappa*1/6,VL-> -100*kappa};
nuLmassive3v3=nuL2DeltaGeneric/.{VR-> kappa*1/6,VL-> -100*kappa};
nvLmassive3v3=nuL2DeltaGeneric/.{VR-> kappa*1/6,VL-> -100*kappa};

nuLmassiveMax=(omega*nuLmassive)/.omega-> 10^-6;
nuLmassive2Max=(omega*nuLmassive2)/.omega-> 10^-6;
nuLmassive3Max=(omega*nuLmassive3)/.omega-> 10^-6;
nvLmassiveMax=(omega*nvLmassive)/.omega-> 10^-6;
nvLmassive2Max=(omega*nvLmassive2)/.omega-> 10^-6;
nvLmassive3Max=(omega*nvLmassive3)/.omega-> 10^-6;


nuLmassive2v2Max=(omega*nuLmassive2v2)/.omega-> 10^-6;
nvLmassive2v2Max=(omega*nvLmassive2v2)/.omega-> 10^-6;
nuLmassive3v2Max=(omega*nuLmassive3v2)/.omega-> 10^-6;
nvLmassive3v2Max=(omega*nvLmassive3v2)/.omega-> 10^-6;
nuLmassive3v3Max=(omega*nuLmassive3v3)/.omega-> 10^-6;
nvLmassive3v3Max=(omega*nvLmassive3v3)/.omega-> 10^-6;


nuRmassiveMax=FindMaximum[omega*nuRmassive, {omega,10^-6,9/200}];
nuRmassive2Max=FindMaximum[omega*nuRmassive2, {omega,10^-6,9/200}];
nuRmassive2v2Max=FindMaximum[omega*nuRmassive2v2, {omega,10^-6,9/200}];
nuRmassive3Max=FindMaximum[omega*nuRmassive3, {omega,10^-6,9/200}];
nuRmassive3v2Max=FindMaximum[omega*nuRmassive3v2, {omega,10^-6,9/200}];
nuRmassive3v3Max=FindMaximum[omega*nuRmassive3v3, {omega,10^-6,9/200}];


uRuRMassiveCompare123V2=LogLinearPlot[{(omega*nuRmassive)/nuRmassiveMax[[1]],(omega*nuRmassive2v2)/nuRmassive2v2Max[[1]],(omega*nuRmassive2)/nuRmassive2Max[[1]],(omega*nuRmassive3)/nuRmassive3Max[[1]],(omega*nuRmassive3v2)/nuRmassive3v2Max[[1]],(omega*nuRmassive3v3)/nuRmassive3v3Max[[1]]},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick} ,{Purple, Thick}}, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,R"] }]
uLuLMassiveCompare123V2=LogLinearPlot[{(omega*nuLmassive)/nuLmassiveMax,(omega*nuLmassive2v2)/nuLmassive2v2Max,(omega*nuLmassive2)/nuLmassive2Max,(omega*nuLmassive3)/nuLmassive3Max,(omega*nuLmassive3v2)/nuLmassive3v2Max,(omega*nuLmassive3v3)/nuLmassive3v3Max},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick} }, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,L"] }]
vLvLMassiveCompare123V2=LogLinearPlot[{(omega*nvLmassive)/nvLmassiveMax,(omega*nvLmassive2v2)/nvLmassive2v2Max,(omega*nvLmassive2)/nvLmassive2Max,(omega*nvLmassive3)/nvLmassive3Max,(omega*nvLmassive3v2)/nvLmassive3v2Max,(omega*nvLmassive3v3)/nvLmassive3v3Max},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick}}, ImageSize-> Full, PlotRange->{0, 1.1},  AxesLabel->{ ω,Superscript["C*ω n","v,L"] }]

nuRmassive3a=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1000*kappa};
nuLmassive3a=nuL2DeltaGeneric/.{VR->kappa*2/3,VL->-1000*kappa};
nvLmassive3a=nvL2DeltaGeneric/.{VR->kappa*2/3,VL->-1000*kappa};
nuRmassive3b=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -500*kappa};
nuLmassive3b=nuL2DeltaGeneric/.{VR->kappa*2/3,VL-> -500*kappa};
nvLmassive3b=nvL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -500*kappa};
nuRmassive3c=nuR2DeltaGeneric/.{VR->kappa*2/3,VL-> -150*kappa};
nuLmassive3c=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -150*kappa};
nvLmassive3c=nvL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -150*kappa};
nuRmassive3d=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -50*kappa};
nuLmassive3d=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -50*kappa};
nvLmassive3d=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -50*kappa};

nuRmassive3e=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -10*kappa};
nuLmassive3e=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -10*kappa};
nvLmassive3e=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -10*kappa};


nuLmassive3aMax=(omega*nuLmassive3a)/.omega-> 10^-6;
nuLmassive3bMax=(omega*nuLmassive3b)/.omega-> 10^-6;
nuLmassive3cMax=(omega*nuLmassive3c)/.omega-> 10^-6;
nvLmassive3aMax=(omega*nvLmassive3a)/.omega-> 10^-6;
nvLmassive3bMax=(omega*nvLmassive3b)/.omega-> 10^-6;
nvLmassive3cMax=(omega*nvLmassive3c)/.omega-> 10^-6;


nuLmassive3dMax=(omega*nuLmassive3d)/.omega-> 10^-6;
nvLmassive3dMax=(omega*nvLmassive3d)/.omega-> 10^-6;
nuLmassive3eMax=(omega*nuLmassive3e)/.omega-> 10^-6;
nvLmassive3eMax=(omega*nvLmassive3e)/.omega-> 10^-6;



nuRmassive3aMax=FindMaximum[omega*nuRmassive3a, {omega,10^-6,9/200}];
nuRmassive3bMax=FindMaximum[omega*nuRmassive3b, {omega,10^-6,9/200}];
nuRmassive3cMax=FindMaximum[omega*nuRmassive3c, {omega,10^-6,9/200}];
nuRmassive3dMax=FindMaximum[omega*nuRmassive3d, {omega,10^-6,9/200}];
nuRmassive3eMax=FindMaximum[omega*nuRmassive3e, {omega,10^-6,9/200}];



uRuRMassiveCompare123V2=LogLinearPlot[{(omega*nuRmassive3a)/nuRmassive3aMax[[1]],(omega*nuRmassive3b)/nuRmassive3bMax[[1]],(omega*nuRmassive3c)/nuRmassive3cMax[[1]],(omega*nuRmassive3)/nuRmassive3Max[[1]],(omega*nuRmassive3d)/nuRmassive3dMax[[1]],(omega*nuRmassive3e)/nuRmassive3eMax[[1]]},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick} ,{Purple, Thick}}, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,R"] }]
uLuLMassiveCompare123V2=LogLinearPlot[{(omega*nuLmassive3a)/nuLmassive3aMax,(omega*nuLmassive3b)/nuLmassive3bMax,(omega*nuLmassive3c)/nuLmassive3cMax,(omega*nuLmassive3)/nuLmassive3Max,(omega*nuLmassive3d)/nuLmassive3dMax,(omega*nuLmassive3e)/nuLmassive3eMax},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick} }, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,L"] }]
vLvLMassiveCompare123V2=LogLinearPlot[{(omega*nvLmassive3a),(omega*nvLmassive3b),(omega*nvLmassive3c),(omega*nvLmassive3),(omega*nvLmassive3d),(omega*nvLmassive3e)},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick}}, ImageSize-> Full, PlotRange->All,  AxesLabel->{ ω,Superscript["ω n","v,L"] }]
uLuLMassiveCompare123V2=LogLinearPlot[{(omega*nuLmassive3a)/nuLmassive3aMax,(omega*nuLmassive3b)/nuLmassive3bMax,(omega*nuLmassive3c)/nuLmassive3cMax,(omega*nuLmassive3)/nuLmassive3Max,(omega*nuLmassive3d)/nuLmassive3dMax,(omega*nuLmassive3e)/nuLmassive3eMax},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick} }, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,L"] }]

nuRmassive3ea=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1*kappa};
nuLmassive3ea=nuL2DeltaGeneric/.{VR->kappa*2/3,VL->-1*kappa};
nvLmassive3ea=nvL2DeltaGeneric/.{VR->kappa*2/3,VL->-1*kappa};
nuRmassive3eb=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/10*kappa};
nuLmassive3eb=nuL2DeltaGeneric/.{VR->kappa*2/3,VL-> -1/10*kappa};
nvLmassive3eb=nvL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/10*kappa};
nuRmassive3ec=nuR2DeltaGeneric/.{VR->kappa*2/3,VL-> -1/100*kappa};
nuLmassive3ec=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/100*kappa};
nvLmassive3ec=nvL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/100*kappa};
nuRmassive3ed=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/1000*kappa};
nuLmassive3ed=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/1000*kappa};
nvLmassive3ed=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/1000*kappa};

nuRmassive3ee=nuR2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/10000*kappa};
nuLmassive3ee=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/10000*kappa};
nvLmassive3ee=nuL2DeltaGeneric/.{VR-> kappa*2/3,VL-> -1/10000*kappa};


nuLmassive3eaMax=(omega*nuLmassive3ea)/.omega-> 10^-6;
nuLmassive3ebMax=(omega*nuLmassive3eb)/.omega-> 10^-6;
nuLmassive3ecMax=(omega*nuLmassive3ec)/.omega-> 10^-6;
nvLmassive3eaMax=(omega*nvLmassive3ea)/.omega-> 10^-6;
nvLmassive3ebMax=(omega*nvLmassive3eb)/.omega-> 10^-6;
nvLmassive3ecMax=(omega*nvLmassive3ec)/.omega-> 10^-6;


nuLmassive3edMax=(omega*nuLmassive3ed)/.omega-> 10^-6;
nvLmassive3edMax=(omega*nvLmassive3ed)/.omega-> 10^-6;
nuLmassive3eeMax=(omega*nuLmassive3ee)/.omega-> 10^-6;
nvLmassive3eeMax=(omega*nvLmassive3ee)/.omega-> 10^-6;



nuRmassive3eaMax=FindMaximum[omega*nuRmassive3ea, {omega,10^-6,9/200}];
nuRmassive3ebMax=FindMaximum[omega*nuRmassive3eb, {omega,10^-6,9/200}];
nuRmassive3ecMax=FindMaximum[omega*nuRmassive3ec, {omega,10^-6,9/200}];
nuRmassive3edMax=FindMaximum[omega*nuRmassive3ed, {omega,10^-6,9/200}];
nuRmassive3eeMax=FindMaximum[omega*nuRmassive3ee, {omega,10^-6,9/200}];



uRuRMassiveCompare123V2=LogLinearPlot[{(omega*nuRmassive3e)/nuRmassive3eMax[[1]],(omega*nuRmassive3ea)/nuRmassive3eaMax[[1]],(omega*nuRmassive3eb)/nuRmassive3ebMax[[1]],(omega*nuRmassive3ec)/nuRmassive3ecMax[[1]],(omega*nuRmassive3ed)/nuRmassive3edMax[[1]],(omega*nuRmassive3ee)/nuRmassive3eeMax[[1]]},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick} ,{Purple, Thick}}, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,R"] }]
uLuLMassiveCompare123V2=LogLinearPlot[{(omega*nuLmassive3e)/nuLmassive3eMax,(omega*nuLmassive3ea)/nuLmassive3eaMax,(omega*nuLmassive3eb)/nuLmassive3ebMax,(omega*nuLmassive3ec)/nuLmassive3ecMax,(omega*nuLmassive3ed)/nuLmassive3edMax,(omega*nuLmassive3ee)/nuLmassive3eeMax},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick} }, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,L"] }]
vLvLMassiveCompare123V2=LogLinearPlot[{(omega*nvLmassive3e),(omega*nvLmassive3ea),(omega*nvLmassive3eb),(omega*nvLmassive3ec),(omega*nvLmassive3ed),(omega*nvLmassive3ee)},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick}}, ImageSize-> Full, PlotRange->All,  AxesLabel->{ ω,Superscript["ω n","v,L"] }]
uLuLMassiveCompare123V2=LogLinearPlot[{(omega*nuLmassive3e)/nuLmassive3eMax,(omega*nuLmassive3ea)/nuLmassive3eaMax,(omega*nuLmassive3eb)/nuLmassive3ebMax,(omega*nuLmassive3ec)/nuLmassive3ecMax,(omega*nuLmassive3ed)/nuLmassive3edMax,(omega*nuLmassive3ee)/nuLmassive3eeMax},{omega,10^-6,1},PlotStyle->{{Blue, Thick, Dashed},{Green, Thick , Dotted},{Red, Thick},{Orange, Thick} ,{Cyan, Thick},{Purple, Thick} }, ImageSize-> Full, PlotRange->{0, 1.1},AxesLabel->{ ω,Superscript["C*ω n","u,L"] }]