Second Order Butterworth LPF Cut-off frequency determining from graph
Whew.  Several problems here.
First, you accidentally labeled your C2 reference designator as 16n.  I don't think it affects the simulation (since you also set the value to 16n), but it's better to fix it so people looking at your schematic can better understand it and make clear comments on it.
Second, you're doing a linear sweep on the frequencies and trying to compare it to a Bode plot which is a decade sweep.  Simply change the lin to dec such that you get: .ac dec 10000 1 5k
Last, your schematic doesn't represent a 2nd order Butterworth, which has complex conjugate poles to acheive the &quot;maximally flat&quot; amplitude characteristic.  What you have is two 1st order RC sections cascaded and then buffered with the opamp at the end.  Not only do each of these have a single real pole, but they &quot;load&quot; each other down so their transfer functions can't be multiplied together without introducing a buffer in between.  If you were trying to match the topology of the circuit shown in the website you referenced (duplicated below), you didn't quite get it 100% correct.  Your C2 should be in feedback from the opamp output and not tied to ground.

For more information regarding the last point, see here:

Is this cascaded RC filter transfer function explanation wrong?
and here:

Transfer function of three cascaded RC filters?