We have measured cross sections for forward neutron production from a variety of targets using proton beams from the Fermilab Main Injector. Measurements were performed for proton beam momenta of 58 GeV/c, 84 GeV/c, and 120 GeV/c. The cross section dependence on the atomic weight (A) of the targets was found to vary as $A^(alpha)$ where $\alpha$ is $0.46\pm0.06$ for a beam momentum of 58 GeV/c and 0.54$\pm$0.05 for 120 GeV/c. The cross sections show reasonable agreement with FLUKA and DPMJET Monte Carlos. Comparisons have also been made with the LAQGSM Monte Carlo.
We report a high statistics measurement of Upsilon production with an 800 GeV/c proton beam on hydrogen and deuterium targets. The dominance of the gluon-gluon fusion process for Upsilon production at this energy implies that the cross section ratio, $\sigma (p + d \to \Upsilon) / 2\sigma (p + p\to \Upsilon)$, is sensitive to the gluon content in the neutron relative to that in the proton. Over the kinematic region 0 < x_F < 0.6, this ratio is found to be consistent with unity, in striking contrast to the behavior of the Drell-Yan cross section ratio $\sigma(p+d)_{DY}/2\sigma(p+p)_{DY}$. This result shows that the gluon distributions in the proton and neutron are very similar. The Upsilon production cross sections are also compared with the p+d and p+Cu cross sections from earlier measurements.