So going off Glen_b's suggestion, I found a logspline package for r. I can't pretend to understand how the splines are calculated (beyond the wikipedia entry), but they do seem to provide a better fit in the tail regions of a distribution.
Here's an example of a fit to 100,000 t distributed points with 3 degrees of freedom:
The red line is the generating function, the black line the density plot, and the blue line is the logspline fit.
A couple of observations:
-The range of the logspline fit needs to be slightly larger than the range of the data, otherwise you'll get a weird inflection at the endpoints. I use a factor of 1.2 to 1.4.
-At the extremes, the lines will always become straight. This is a feature of the splines, not of the data. It's important to keep this in mind because otherwise it may fool you into thinking the tails are decreasing exponentially when they aren't (see right side of picture above).
-As such, while logsplines do increase the range at which you can visually fit the data, they also break down at extreme ranges and cannot be relied on.
-It seems to have trouble with highly concentrated data sets. I have one with ~180000 data points between 0 and 0.5 (most of which are at single particular values) and another ~500 between -6 and 6, and the package just can't deal with the data adequately. I can get the logspline function to run when a little bit of gaussian jitter is added, but the fit is weak. Otherwise the package errors out.