How to analyse data on metric space? I just ran an experiment and I'm not sure how best to analyse the data.  My data are distance values between objects in a metric space bounded by [0,1].  I have drawn up a probability density estimate as follows:
This looks to me like a kind of hybrid between an exponential and normal distribution, what tests can I run to make better sense of this?  

More background:
This metric space is a finite corpus of Documents, the distance between which indicates their similarity: if the distance is 0 they are identical, if 1 they have no commonality.  This sample represents all the distances between 1,000 randomly selected documents.
 A: I can't answer the question of what is the distribution of the distances, but I can shed some light on why you are seeing that very sharp narrow peak.
This is an aspect of the dimensionality curse, known as distance concentration. To the best of my knowledge the first paper on this phenomenon was Beyer et al "When is nearest neighbor meaningful" where the authors demonstrate that in high dimensional probability spaces one can expect the distances between points to converge to a common value with high probability. The phenomenon is very general and holds for any metric (for example), although some converge slower than others.
Some relevant links are in the post by @Denis (and my answer there).
A: One simple way would be to model your data as beta distributed. The beta is by definition between 0 and 1:
xx <- seq(.01,.99,by=.01)
plot(xx,dbeta(xx,shape1=12,shape2=2),type="l",xlab="",ylab="")


For an added twist, you could model the little bump at zero by mixing the beta with a point mass there. This is commonly done to add "additional" zeros to the Poisson distribution ("zero-inflated Poisson"), but it looks like it could be helpful in your case, too.
A: A circular normal distribution is of course unbounded and so it would not be an example where all the distances would fall in a unit circle but if we make the variance small enough the tail outside the circle is small.  The density for the distance of a point from the center of the circular normal distribution is called a Rayleigh distribution and it has the skewness behavior similar to your plot.  Remember $D=\sqrt{X^2+Y^2}$.  So the shape of the distribution is primarily dictated by the function of $X$ and $Y$ that it is.  Note that if $X$ and $Y$ were independent standard normal distributions $D$ would be the square root of a $\chi^2$ with 2 degrees of freedom.
