With the information that you've provided, it's difficult to really understand what you did to set up your analysis to produce your results. When I ran your model as-is, I got some very odd dCP on the horizontal tail that caused a very poor solution.
For starters:
* Did you include the fuselage and engines in this analysis or just the lifting surfaces? It's almost always best to start with a single wing and find some reasonable results before moving on to a full aircraft. You're using Sets, which is good. Start simple and build up from there.
* You're running fairly large alpha sweeps at transonic conditions. VSPAERO is essentially a potential flow solver with compressibility corrections that can solve for supersonic flow in VLM (Thin) mode but will not give very reasonable results with Thick components at supersonic speeds (for now). Transonic conditions are problematic for many solvers and VSPAERO is no exception. Try to recreate the results at M 0.5 and see what you get.
* I suggest increasing your wing resolution (Num U, Num W, clustering) to try and capture more of the pressure changes along the wing. I suspect that your wing is a bit too coarse to get a good load distribution.
* VSPAERO has some parasitic drag (CD0) included in the results but this is essentially a corrected flat-plate model based on the NACA 0012 for wings and a regular flat plate for bodies. You may have better luck adding in some drag from the parasite drag tool.
* VSPAERO does not model separation but does include some experimental stall modeling. You should find that the stall effects will significantly alter your results at higher alpha.
Work with these suggestions and then report back with your results and some more detail about how you set up your analysis. Then we can help you more.
Good luck!