computational physics - Why can't we set the lattice spacing 'a' in lattice QCD?

My question is to do with lattice QCD. In the lattice action there is a parameter, 'a', the lattice spacing in physical units.

However, if we want to generate a configuration with a certain lattice spacing, we don't just set a=some number. We take the roundabout route of setting some dimensionless parameters (with some dependence on 'a') equal to some value and then extracting what 'a' must have been afterwards.

My question is why do we have to do this? Why can't I just go into the code and set a=something and simulate? I know that for other computational problems it makes sense to only simulate dimensionless quantities, but the reasons for this are to do with stability and reducing the computational time.

Answer

The reason is that changes in the scale a and changes in the coupling g can compensate for each other. Two simulations, one with a small lattice spacing a and gauge couling g and another with an even smaller lattice spacing a' and coupling g' give the same results at long distances when g' is adjusted properly. This is the statemen that the theory is renormalizble, so that you can take the limit a goes to zero, g' changes correspondingly, and extract a good limit. Further, as the lattice spacing a' gets smaller, to keep the physics the same, g' gets weaker. This is the statement that QCD is asymptotically free (free at short distances).

But the dependence of g on a for not-so-small lattices is annoying to calculate, because it only is simple in covariant regulators, and the log-running means it is never that small at any reasonable scale. So instead of fixing a and calculating what g should be, you use the known existence of the scaling continuum limit to fix your physics. So you just set your length scale to make a=1 and you adjust g to be approximately .5 (this makes ${g^2\over 2\pi} = .04$, and this is the perturbative parameter), and then you look to see if the gauge field randomizes over your box with this choice.

If you make g too small, the gauge field will be nearly constant in the box, if you make g too big, the typical gauge field configuration will be random from point to point, with a large SU(3) matrix for plaquettes. You want to make sure that your g is in the good spot so that the box is not too small to see the long-distance randomness, and not too big to make the lattice coarse to see the interior structure of a hadron bag.

Because the choice of g and a are intertwined in a nontrivial way, it is best to fix the simulation parameters by using the output masses. The dependence of g and a cannot be extracted from traditional dimensional analysis, because it is logarithmic in a. Classically, g is independent of a.