Clear[z, sig];  
    sig = 0.29;
T = 
      9/(4*P^2) + 9*4/(4*Q^2) + 2/R^2; 
Vne = -3*z/P - 3*z/Q - z/R; 
Vee = (2/Q)*(3 - (6/10)*(P/Q)^2) - 6*sig*P^2/Q^3  +  
        2*(1 - (3/20)*(P/R)^2 + (1/32)*(P/R)^3)/R + 
        2*(1 - (3/20)*(Q/R)^2 + (1/32)*(Q/R)^3)/R + 
        3*sig/P; 
func = T + Vne + Vee;  
    Vir  =  (Vne + Vee)/T; 
  
t  =  
      Table[FindMinimum[
          func, {P, 0.3}, {Q, 1.8}, {R, 
            3.0}, {MaxIterations -> 100}], {z, 5, 18}]  
Table[
      N[Vir  /.  t[[{z - 4}, 2]], 10], {z , 5, 
        18}]
Table[
      N[Vne  /.  t[[{z - 4}, 2]], 10], {z , 5, 
        18}]
Table[
      N[Vee  /.  t[[{z - 4}, 2]], 10], {z , 5, 
        18}]
  
Table[
      N[t[[{z - 4}, 1]]*27.21138, 10], {z, 5, 
        18}]; 
PA = 
      Table[t[[{z - 4}, 2, 1]], {z, 5, 
          18}]; 
QA = 
      Table[t[[{z - 4}, 2, 2]], {z, 5, 
          18}]; 
RA = 
      Table[t[[{z - 4}, 2, 3]], {z, 5, 
          18}]; 

Do[Show[
        Graphics[{Circle[{0, 0}, P  /.  PA[[i]]], 
            Circle[{0, 0}, Q  /.  QA[[i]]], 
            Circle[{0, -R /.  RA[[i]]}, 
              R  /.  RA[[i]]], 
            Circle[{0, R  /.  RA[[i]]}, 
              R  /.  RA[[i]]]}], {AspectRatio  ->  Automatic, 
          Axes  ->  True, GridLines  ->  Automatic,  
          PlotRange  ->  {{-6.04, 6.04}, {-6.04, 6.04}},  
          Frame  ->  True}],  {i, 1, 14}]