Thank you guys for your hints.
I now extend the frequency range to 40G. Still, I would say the trend looks not right.
Previously I set the sweep type to interpolation. But later I found it sometimes gives me very strange result (for example, non-monotonic behavior). So now I change the type to discrete. Are there any guidelines as to when to use which? I checked the help, but didn't get it. If I want to simulate from 0.1G to 100G, which shall I use? Or do I have to break the range to a few?
Another thing I noticed is still the initial drop at low frequencies. I have been reading Howard Johnson's black magic book, and it talked about the transition from the RC region to the LC region. However, in both regions the rate should follow sqrt(freq). Yet in my simulation results, I didn't see such behavior in the dB-log plot. I am running more simulation to see if what I see is the transition region between the RC and LC regions.