Voyager 1 data from beyond the heliopause provide the first direct measurements of the interstellar cosmic ray spectra below 1 GeVnuc. In this paper we combine these Voyager measurements of H and He nuclei from 3-600 MeVnuc with higher energy measurements at 1 AU from the BESS and PAMELA experiments up to 100 GeVnuc. Using a Weighted Leaky Box Model for propagation in the galaxy, we obtain an excellent fit to these new Voyager observations and the much higher energy spectra up to 100 GeVnuc by using source spectra which are P-2.28, with the exponent independent of rigidity from low to high rigidities; along with a rigidity dependence of the diffusion path length which is P-0.5 at rigidities 1.00 GV, and possibly changing to P1.0 at lower rigidities.

The new Voyager measurements of cosmic ray electrons between 6-60 MeV beyond 122 AU are very sensitive indicators of cosmic ray propagation and acceleration in the galaxy at a very low modulation level. Using a Monte Carlo diffusion model with a source spectrum with a single spectral index of -2.2 at all energies we are able to fit this observed Voyager spectrum and the contemporary PAMELA electron spectrum over an energy range from 6 MeV to ~200 GeV. This spectrum has a break in it but this break is due to propagation effects, not changes in the primary spectrum. This break is gradual, starting at > 2 GeV where the spectrum is ~E^-3.2 and continuing down to ~100 MeV or below where the spectrum becomes ~E^-1.5. At the higher energies the loss terms due to synchrotron radiation and inverse Compton effects which are ~E^2.0 steepen the exponent of the source spectrum by 1.0. At lower energies, these terms become unimportant and the loss is governed by diffusion and escape from the galaxy. A diffusion term which is proportional to beta^-1 below ~0.32 GV (which also fit the H and He spectra measured at Voyager) and has a value = 3x1028 cm^2 x s^-1 at 1 GV and a boundary at +/-1 Kpc will fit the Voyager or other similar spectra at low energies.

After the disappearance of lower energy heliospheric particles at Voyager 1 starting on August 25th, 2012, spectra of H, He and C/O nuclei were revealed that resembled those to be expected for galactic cosmic rays. These spectra had intensity peaks in the range of 30-60 MeV, decreasing at both lower energies down to a few MeV and at higher energies up to several hundred MeV. We have modeled the propagation of these particles in the galaxy using an updated Leaky Box Diffusion model which determines the spectra of these components from ~2 MeV to >200 GeV. The key parameters used in the model are a galactic input spectrum ~P^-2.24, the same for all components and independent of rigidity, and a diffusion coefficient that is ~P^0.5 above a lower rigidity and increases ~beta^-1.0 below a lower rigidity ~0.56 GV. These same parameters also fit the high energy H and He data from ~10-200 GeV/nuc from the PAMELA and BESS experiments. The new Voyager spectra for all three nuclei are thus consistent with rigidity spectra ~P^-2.24 from the lowest energies to at least 100 GeV. Deviations from this spectrum can reasonably be attributed to propagation effects. Some deviations between the calculated and newly observed spectra are noted, however, below ~30 MeV/nuc, particularly for C/O nuclei, that could be significant regarding the propagation and sources of these particles.

This paper describes the unfolding of the solar modulated galactic cosmic ray H and He nuclei spectra beyond ~105 AU in the heliosheath. Between 2008.0 and 2012.3 when Voyager 1 went from about 105 to 120.5 AU the spectral intensities of these two components between about 30 and 500 MeV/nuc unfolded (increased) in a manner consistent with an average modulation potential decrease ~5 MV per AU as described by a Parker like cosmic ray transport in the heliosphere where the overall modulation is described by a modulation potential in MV. Between 120.5 and 121.7 AU, however, as a result of two sudden intensity increases starting on May 8th and August 25th, 2012, this modulation potential decreased by ~80 MV and spectra resembling possible local interstellar spectra for H and He were revealed. Considering these spectra to be the local interstellar spectra would imply that almost 1/3 of the total modulation potential of about 270 MV required to explain the spectra of these components observed at the Earth must occur in just a 1 AU radial interval in the outer heliosheath. As a result about ~80% of the total modulation potential observed at the Earth at this time occurs in the heliosheath itself. The remaining 20% of the total modulation occurs inside the heliospheric termination shock. The details of these intensity changes and their description by a simple modulation model are discussed.