This is what i did and what Stata is doing:

ppml EXP2 lnyy lnpilprocapite lndistcity lnsupsup BORDER UE Euro Lang Landlo_dest Landlo_orig lnmount Autonomy Presence Riconoscimento cris
> is _OPR* _DDES* _PPRO*, cluster(pairid)
note: checking the existence of the estimates
WARNING: EXP2 has very large values, consider rescaling
WARNING: lnyy has very large values, consider rescaling or recentering
WARNING: lnpilprocapite has very large values, consider rescaling or recentering
WARNING: lnsupsup has very large values, consider rescaling or recentering
note: starting ppml estimation

Iteration 1: deviance = 3.74e+12
Iteration 2: deviance = 2.69e+12
Iteration 3: deviance = 2.52e+12
Iteration 4: deviance = 2.49e+12
Iteration 5: deviance = 2.49e+12
Iteration 6: deviance = 2.49e+12
Iteration 7: deviance = 2.49e+12
Iteration 8: deviance = 2.49e+12
Iteration 9: deviance = 2.49e+12
Iteration 10: deviance = 2.49e+12
Iteration 11: deviance = 2.49e+12
Iteration 12: deviance = 2.49e+12
Iteration 13: deviance = 2.49e+12
Iteration 14: deviance = 2.49e+12
Iteration 15: deviance = 2.49e+12

Starting from that deviance has the same value.