The complete Joule-Thomson (JT) inversion curve for carbon dioxide is calculated using molecular simulations. A two center Lennard-Jones model with an embedded point quadrupole is used to model the fluid-fluid interactions. The simulation results agree quantitatively with all available experimental data. Comparison with commonly used equations of state provides only a modest agreement, with the highest discrepancies being observed at the high temperature branch of the inversion curve.