Abstract
The Lieb-Robinson correlation function captures propagation of quantum correlations in a many-body system. We calculate the value of the leading order of the correlation function, not its bound, for a system of interacting qubits at early times. The general analytical result is compared to numerical calculations and is applied to regular qubit lattices in one, two, and three dimensions. The Lieb-Robinson velocity and the approximately exponential leading edge of correlations emerge in the limit of large arrays.