  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.
