Greetings, I'm going through the fun task of optimising an RF PA, in the UHF range. I've only recently started with RF PA design, having done bunch of LNAs and other low power / 50 Ohm stuff in the past.
For a PA, your bigger issue is you have a non-linear design. Harmonic balance simulators, like those available in ADS and Microwave Office are ideal, provided you have good models. Both programs are about the same, though people have their individual preferences. But if you have load pull data or even if you just estimate your optimum load using the load line method, you can get pretty good results with a linear simulation, at least close enough to be tunable.
With both of those programs, everything is an extra cost add-on. Harmonic balance, 2.5D EM, etc. are all options.
Sonnet software is also fine for 2.5D EM. But it can't do other types of simulation.
Remember if you do non-linear or harmonic balance simulation, all your models need to cover frequencies from DC out to several harmonics of your carrier frequency.