Observable consequences of the hypothesis that the observed universe is a numerical simulation performed on a cubic space-time lattice or grid are explored. The simulation scenario is first motivated by extrapolating current trends in computational resource requirements for lattice QCD into the future. Using the historical development of lattice gauge theory technology as a guide, we assume that our universe is an early numerical simulation with unimproved Wilson fermion discretization and investigate potentially-observable consequences. Among the observables that are considered are the muon g-2 and the current differences between determinations of alpha, but the most stringent bound on the inverse lattice spacing of the universe, b^(-1) >~ 10^(11) GeV, is derived from the high-energy cut off of the cosmic ray spectrum. The numerical simulation scenario could reveal itself in the distributions of the highest energy cosmic rays exhibiting a degree of rotational symmetry breaking that reflects the structure of the underlying lattice.

The dependence of the nuclear force on standard model parameters plays an important role in bounding time and space variations of fundamental couplings over cosmological time scales. We discuss the quark-mass dependence of deuteron and di-neutron binding in a systematic chiral expansion. The leading quark-mass dependence of the nuclear force arises from one-pion exchange and from local quark-mass dependent four-nucleon operators with coefficients that are presently unknown. By varying these coefficients while leaving nuclear observables at the physical values of the quark masses invariant, we find scenarios where two-nucleon physics depends both weakly and strongly on the quark masses. While the determination of these coefficients is an exciting future opportunity for lattice QCD, we conclude that, at present, bounds on time and space variations of fundamental parameters from the two-nucleon sector are much weaker than previously claimed. This brings into question the reliability of coupling-constant bounds derived from more complex nuclei and nuclear processes.