all random code should use default_rng

[drammock]

our minimum supported NumPy is 1.23, which has np.random.default_rng(). We should change all instances in our code to use that instead of legacy calls to randn() or RandomState()