I believe this is the "Nyquist Frequency".
Similar, but not the same.
The problem addressed by Nyquist is whether a certain frequency can still be represented in a discrete signal. This is not really interesting in this case: vibrations faster than a certain frequency aren't visible anyway, so it doesn't matter if the simulation can't handle them. Stability of the simulation is more important; see for instance this page. It has a nice picture of the problem:
I haven't found a mathematical justification, but many people seem to follow rules of thumb like "the step size needs to be at least 5 to 10 times smaller than the smallest time scale in the model". For quantitative engineering purposes (e.g. testing the structural integrity of an airplane design in a computer model), I've even heard they use 100 time steps for one period of the fastest oscillation mode in the model. In the picture above, the Euler method is similar to what I use. Orbiter uses Runge Kutta in many cases, which has the same problem, but less dramatically. The green line in the picture has a step size about 2 times larger than the time scale of the system. With smaller step sizes, the Euler method would also approach the exact solution.