Quantum dynamical theories have progressed to the stage in which state-to-state differential cross sections can now be routinely computed with high accuracy for three-atom systems since the first such calculation was carried out more than 30 years ago for the H + H(2) system. For reactions beyond three atoms, however, highly accurate quantum dynamical calculations of differential cross sections have not been feasible. We have recently developed a quantum wave packet method to compute full-dimensional differential cross sections for four-atom reactions. Here, we report benchmark calculations carried out for the prototypical HD + OH → H(2)O + D reaction on an accurate potential energy surface that yield differential cross sections in excellent agreement with those from a high-resolution, crossed-molecular beam experiment.