plot1 = ListPlot[MNvsAmp, PlotRange -> All, PlotStyle -> PointSize[Large], ImagePadding -> 85, Frame -> {True, True, True, False}, FrameLabel -> {"Time, s", "Number of atoms, 1000"}, PlotMarkers -> {marker1, .035}, PlotLegends -> "gamma = 1.903 beta = 2.173*10^(-20)"];
plot2 = ListPlot[TvsTXTvalue, PlotRange -> All, PlotStyle -> PointSize[Large], ImagePadding -> 85, Frame -> {False, False, False, True}, FrameTicks ->{None, None, None, All}, FrameLabel -> {{"","Temperature, mK"},{"",""}}, PlotMarkers -> {marker2, .035}, PlotLegends -> "gamma = 1.903 beta = 2.173*10^(-20)"];
plot3 = Plot[Normal[bettafit], {tt, 3.7, 4.4}, PlotStyle -> Directive [Thick], Axes -> {False, False}]
Overlay[{plot1, plot2, plot3}]