I've been testing out a few circuit simulation applications and ran into an unusual issue with Qucs and the 2N3904 transistor. For a simple NPN transistor circuit with load resistor \$R_L\$, base resistor \$R_B\$, and pull-down resistor \$R_P\$, I'm seeing an inflection in the collector-to-emitter voltage \$V_{CE}\$ as I sweep \$R_B\$ through a range of values. It also reads much higher than it should, starting out around 0.4 V and increasing to nearly 5 V. In the screenshot below, you can see how \$V_{CE}\$ rises sharply and \$I_C\$ begins to decline when \$R_B\$ is around 4 kOhm.
When I simulate the same circuit in LTSpice, the results are more typical: \$V_{CE}\$ is less than 0.2 V while in saturation and there's a nice flat \$I_C\$ curve around 120 mA or so, which correlates to the supply voltage and load resistor.
I'm thinking there must be an incorrect parameter in the 2N3904 component in Qucs that is throwing things off. I've tried to compare the parameters between the different software packages, but haven't had much luck (some of the data sheets don't even have all of those parameters listed). Interestingly, if I try to substitute a 2N2222, the inflection is gone, but I think \$V_{CE}\$ is still too high since the transistor should be in saturation.
Am I doing something wrong, or is my Qucs configured incorrectly? I'm trying to evaluate a variety of simulation packages and I really like Qucs, but I want to be sure the simulation results are going to be accurate going forward. Thanks for the help.