The van der Waals phase transition

F[t_,v_] := -t*(1+Log[(v-1/3)*t^(3/2)])-9/(8 v)

Manipulate[Plot[F[t,v],{v,0,10},PlotRange->{{0,8},{-3.5,-0.8}},PlotStyle->{Thickness[0.007]},AxesLabel->{"V","F(T,V)"},BaseStyle->{FontSize->22},ImageSize->700],{{t,1.2},0.75,1.2},LabelStyle->{Large}]

p[t_,v_] := Evaluate[-D[F[t,v],{v,1}]]

Manipulate[ParametricPlot[{{v,p[t,v]},{v,3/8},{v,0.105954}},{v,0.4,8},PlotRange->{{0,8},{-0.35,0.5}},AspectRatio->1,PlotStyle->{Thickness[0.007],Thickness[0.003],Thickness[0.003]},AxesLabel->{"V","P(T,V)"},BaseStyle->{FontSize->22},ImageSize->700],{{t,1},0.75,1.2},LabelStyle->{Large}]

HV[t_,v_,p_] := F[t,v]+p*v

Manipulate[Plot[HV[t,v,p],{v,0.334,10},PlotRange->{{0,10},{-1.6,-0.5}},PlotStyle->{Thickness[0.007]},AxesLabel->{"V","F(T,V)+PV"},BaseStyle->{FontSize->22},ImageSize->700],{{t,1},0.75,1.2},{{p,3/8},0.05,0.8},LabelStyle->{Large}]

Manipulate[ParametricPlot[{{p[t,v],HV[t,v,p[t,v]]}},{v,0.1,20},PlotRange->{{-0.35,0.5},{-2,-1}},AspectRatio->1,PlotStyle->{Thickness[0.007]},AxesLabel->{"P","G(T,P)"},BaseStyle->{FontSize->22},ImageSize->700],{{t,1},0.75,1.2},LabelStyle->{Large}]