Certainly nothing keeps you from avergaging bin $0$ and bin $1$; the question is what you're trying to calculate. The FFT calculates a discrete Fourier transform, and you're trying to interpret it as something continuous. That's possible, but it requires a careful interpretation.
One interpretation goes like this: The data you transform represent samples taken over a finite interval with equal spacing. The Fourier coefficients that the FFT computes are the discrete Fourier transform of the data, which you can interpret as the Fourier sum for the unique periodic function that interpolates the sampled values and is band-restricted by the Nyquist frequency. You can view these Fourier coefficients as the coefficients of delta peaks at the corresponding frequencies in a continuous Fourier integral. Now if you restrict the periodic function to a single period, you effectively multiply it by a rectangular window, and in Fourier space that corresponds to convoluting it with the $\operatorname{sinc}$ function. Thus each of the delta peaks gets replaced by a copy of the $\operatorname{sinc}$ function centred at the corresponding frequency. Since these $\operatorname{sinc}$ functions are $1$ at their central frequency and $0$ at the frequencies corresponding to all other coefficients, this gives rise to a continuous function that interpolates between the discrete Fourier coefficients. It's this function that has a well-defined interpretation, namely as the Fourier transform of a single period of the band-limited interpolation of your data.
Thus, to get something meaningful, you'd have to interpolate with $\operatorname{sinc}$ functions, so you'd need a full sum over all the coefficients to get a single interpolated value, not just a combination of the two adjacent coefficients. However, if your data are sufficienly well-behaved, it will often be a good approximation to do linear interpolation instead as you suggested.