Perhaps you are bound by a regulatory agency of some sort. In that case, you must yield to what they say to do. My argument is a philosphical one so here goes.
The correlation coefficient tells you some things but it's not the end-all. I took your data and I calculated the 1-parameter (y = mx) least-squares slope for it. It came out 42439. Your 2-parameter fit (y = mx + b) yielded 43446. Those are within 2% of each other, which is very good agreement in my opinion. To me, this is the first indicator that the intercept is not terribly large compared to the instrument responses used to calculate it. Know that the slope and intercept in the 2-parameter fit are not independent of each other AND the intercept is always going to be the dumping ground for the majority of the error in your calibration data. The intercept is actually dictated by the midpoint of your calibration data (average of x and average of y).
I had to do some fooling around to get your data to work in the US version of Excel so I could have some roundoff error going on. Please try to ignore the fact that I didn't get the exact values for your slope and intercept.
I also created some data to overlay with yours. Series 1 in the graphic is your data. Series 2 is:
y = 42439*x^0.8 (this data slightly concave down at higher concentrations)
and Series 3 is:
y = 42439*x^1.2 (this data is slightly concave up at higher concentrations)
In effect, these treatments add nonlinearity to the data. Here's the plot:
https://s20.postimg.org/c2oc3r2xp/Linea ... mple_I.jpg
The R^2 values are still "good" yet there's a fair amount of curvature in the data.
If your samples are running in the middle of your calibrated range, you are far from zero where the intercept would could/would have a big effect on your calculated concentration. Far from the origin:
(y - b)/m ~ y/m
Take an instrument response of 1x10^7 cts:
1x10^7/43446 = 230
Accounting for the nonzero intercept you get (my 2 parameter fit says the intercept is -231607):
(1x10^7 + 231607)/43446 = 236
1x10^7/42439 = 236 (1 parameter least-squares response factor)
How big of a deal is 230 vs. 236 in your application? Obviously, as you get closer to the origin, the intercept becomes a bigger factor in your calculation. Your relative standard deviations on your triplicate analyses of your standards are wonderful (0.5% or less). Perhaps being off by 2% on the "sample" I described above is a problem?
I fit your data to y = a1*x + a2*x^2 and the a2 term is small but positive. This is probably why your intercept is negative (your data is slightly concaved upward).