Here are the largest common divisors known for small matrices. The larger non-square ones were sent by Philippe Fondanaiche (and are not known to be optimal). The 6×6 annd 7×7 were sent by Ulrich Schimke.
Claudio Baiocchi also considered the cases where each entry can/must be two digits. The matrix has a common divisor of 89, and the matrix has a common divisor of 288.
Philippe Fondanaiche extended this to larger matrices by finding with common divisor 1085 and with common divisor 6869.
Antonio Ianiero suggested using cyclic numbers to form large matrices of this form, like the one below. Unfortunately, each number appears twice, once in a row and once in a column.
Johan de Ruiter considered small bases. He proved there are no such matrices in base 2, but there are square matrices in all larger bases. He found the largest common divisor for some small matrices in base 3. He also considered factorial bases.
Ulrich Schimke pointed out that the common divisor can obviously never be a multiple of 10, and asked whether all other numbers could be. Also, he found that if the common divisor divides 10n–1 and a solution has no "9" in the first row or column then exchanging all digits "k" with "9–k" gives another solution.