Using Leaky Box Model propagation calculations for H nuclei and a Monte Carlo diffusion propagation model for electrons, starting from specific source spectra, we have matched the observed LIS spectra of these cosmic rays measured by Voyager at lower energies and AMS-2 at higher energies, a range from ~10 MeV to ~1 TeV. The source spectra required are very similar rigidity spectra. Below ~6-10 GV the source spectra for both particles are ~P-2.25 and above 10 GV the spectra are ~P-2.36-2.40. This break in the source spectral index is not seen for He and C nuclei in a match of Voyager and AMS-2 intensities both of which have source rigidity spectra with an index ~-2.24 throughout the entire range of measured energies from ~10 MeV//nuc to ~1 TeV/nuc. The absolute source intensities of electrons and H nuclei are derived and the source ratio of accelerated electrons to H nuclei is between 2-5%. The total number of accelerated electrons is much greater than that for protons, however, because the accelerated electron spectrum extends down to ~1-2 MV rigidity whereas the H nuclei spectrum cannot be observed below ~50-100 MV because of ionization energy loss. Most of these low energy electrons escape from the galaxy forming an intergalactic background.

We have used new measurements of the B/C ratio in galactic cosmic rays at both low and high energies by the Voyager and AMS-2 spacecraft, respectively, along with propagation calculations using a truncated LBM to examine the implications of these new measurements over an extended energy range from a few MeV/nuc to 1 TeV/nuc. We find that the predictions from both the truncated LBM and the Diffusive Reacceleration model for GALPROP both agree with the Voyager and AMS-2 measurements of the B/C ratio to within +/- 10 percent throughout the entire energy range from 50 MeV/nuc to 1 TeV/nuc. The two propagation approaches also agree with each other to within +/-10 percent or less throughout this energy range. In effect a diffusion model, without significant additional acceleration, provides a match within +/-10 percent to the combined data from Voyager 1 and AMS-2 on the B/C ratio from 50 MeV/nuc to 1 TeV/nuc. The B/C ratio below 50 MeV/nuc measured at V1 exceeds the predictions of both propagation models by as much as 3 sigma in the data measurement errors.

Voyager observations for over 3 years beyond the heliopause have started to define features of the low energy, less than 100 MeV per nuc, cosmic ray secondary nuclei that have a zero or negligibly small source component. As an example, the abundance of B between about 7 and 15 MeV per nuc is unexpectedly large, greater than the prediction of a LBM by 2 or 3 in the measurement and cross section error. On the other hand, for several other heavier secondary nuclei with a low source abundance such as F and Z equal 17 to 19 and 21 to 23 nuclei, in the corresponding energy channels between about 10-20 MeV/nuc, zero nuclei have been observed. The same LBM calculations would predict about 6 to 7 events for the sum of these three groups of nuclei. The B observed intensities could be more closely matched by considering a nested LBM with 0.3-1.0 g/cm2 of matter near the cosmic ray sources, essentially a source component of B. This nested LBM calculation, if extended to the production of the above groups of secondaries, would however predict a total 8 to 12 events thus leading to an even larger discrepancy with the zero events that are observed. The measurements of heavier secondaries at low energies therefore make it very problematical that a nested LBM with more than a few 0.1 g/cm2 or more of matter near the sources could be the source of the large B abundance seen by Voyager. The large abundance of B at the lowest energies is best understood, so far, as the result of production mostly from a matter traversal of 10 g/cm2 in the galaxy, possibly in combination with uncertainties in the data and cross sections at these low energies.

This paper seeks to find an explanation of the galactic cosmic ray B/C ratio newly measured in cosmic rays between ~10 MeV/nuc and 1 TeV/nuc. Voyager measurements of this ratio at low energies and AMS-2 measurements at high energies are used in this study. These measurements both considerably exceed at both low and high energies the ratio predicted using a simple Leaky Box Model for propagation of cosmic rays in the galaxy. Between 1-70 GeV/nuc, however, this same model provides an excellent fit (within +2-3%) to the new AMS-2 measurements using an escape length ~P-0.45. This would imply a diffusion coefficient ~P0.45, very close to the Kraichian cascade value of 0.50 for the exponent. Extending this same diffusion dependence to high energies, still in a LBM, along with a truncation of short path lengths in the galaxy will predict a B/C ratio of 4.5% at ~800 GeV/nuc which is very close to the AMS-2 measurement at that energy. This would indicate that the amount of material traversed near the source in a Nested LBM, for example, is less than about 0.5 g/cm2 at these energies since this process, with additional matter near the sources, would increase the B/C ratio. At low energies, however, the B/C ratio of 14% +3% at ~10 MeV/nuc measured by Voyager is more difficult to explain. The same parameters used in a LBM that fit the high energy B/C measurements predicts a B/C ratio of only ~4% at 10 MeV/nuc where the path length is about 1.5 g/cm2. If the cosmic rays have traversed ~10 g/cm2 of material in the galaxy at low energies the prediction for B/C is still only 8%, which is between 2.0-3.0 sigma below the measurement and its errors.

The intensities of the low energy part of the spectra of primary cosmic ray nuclei including C, O, Ne, Mg, Si and Fe measured by Voyager 1 beyond the heliopause are deficient relative to the spectra measured at energies above ~100 MeV/nuc as calculated using a standard Leaky Box Model with the path length a function of rigidity. Modifications to the normal exponential distribution of path lengths at a single rigidity as is used in a simple LBM will provide a good fit to this new Voyager 1 data at low energies. These modifications, sometimes called a truncation, lead to a deficiency of short path lengths relative to an exponential distribution. This deficiency in the intensities can be described by a truncation parameter = 0.12. This modification can be produced in several ways including a non-uniform local distribution of cosmic ray sources. A uniform source distribution in the galactic plane that is deficient in sources within 0.2-0.4 Kpc of the Sun is indicated by the data. Further studies of these low energy spectra in more detail will improve these estimates and help define other features of the local galactic distribution of these cosmic rays.

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.

We discuss two extraordinary increases of cosmic ray intensity observed by Voyager 1 in the last 1.1 AU before heliopause crossing, Aug 2012 at 121.7 AU. The two increases are roughly similar in amplitude and result in a total increase of cosmic ray nuclei around 1 GV of over 50 percent and of 0.01 GV electrons of a factor about 2. During the first increase, the changes in the magnetic, B, field are small. After the first increase, the B field variations and cosmic ray changes become large and during the second increase the B field variations and the cosmic ray changes are correlated to within a day. The intensity variations of H and He nuclei during these time intervals are measured from 0.1 to over 1 GV. The total GCR increse over the two events resemble those expected from a simple force-field 'like' solar modulation model with a modulation potential of about 80MV. This is nearly one third of the total modulation potential of about 250 GV required to produce the modulation of these particles observed at the earth at the 2009 sunspot minimum and adds a new aspect to ideas about heliospheric modulation.

Voyager 1 has explored the solar wind-interstellar medium interaction region between the terminal shock and heliopause following the intensity distribution of galactic cosmic ray protons above 200 MeV energy. Before this component reached the galactic level at 121.7 AU, 4 episodes of rapid intensity change occured similar to the Forbush Decreases found near the sun, rather than the expected result of models related to those describing Long Term Modulation in the inner solar system. Because the mean solar wind flow is both expected and observed to be perpendicular to the radial direction close to the heliopause, explanation is given in terms of transient radial flows related to possible heliopause boundary flapping. It is necessary that radial flows are at the sound speed found for conditions downstream of the teminal shock and that the relevant perpendicular cosmic ray diffusion is controlled by 'slab' field fluctuations accounting for 20 percent or less of the total power in field variance. However, additional radial drift motion related to possible north to south gradients in the magnetic field may allow the inclusion of some diffusion according to 2-D turbulence theory. The required field gradients may arise due to variation in the field carried by the solar plasma deflected away from the solar equatorial plane. Modulation amounting to a total 30 percent drop in galactic intensity requires explanation by a combination of several transient episodes.

Voyager 1 has explored the solar wind-interstellar medium interaction region between the Terminal Shock and the Heliopause, following the intensity distribution of the shock accelerated anomalous component of cosmic rays in the MeV energy range. The sudden disappearance of this component at 121.7 AU from the sun is discussed in terms of three models for the transition into the interstellar plasma flow. Particles trapped flowing parallel to the boundary may penetrate up to one Larmour radius beyond. If the boundary is stationary, Voyager 1 directly samples this distance. The boundary could flap, depending on Heliosheath pressure changes and Voyager 1 then samples the extent of this motion. Finally, a turbulent boundary layer is considered in which the MeV particle distribution falls off with distance, thus measuring diffusion within the layer.

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.

Voyager 1 has entered regions of different propagation conditions for energetic cosmic rays in the outer heliosheath beginning at a distance of about 111 AU from the Sun. This conclusion is based on the fact that the low energy 6-14 MeV galactic electron intensity suddenly increased by ~20% over a time period \leg 10 days and the electron radial intensity gradient abruptly decreased from ~19%/AU to ~8%/AU at 2009.7 at a radial distance of 111.2 AU. A sudden radial gradient change was also observed at this time for >200 MeV protons. The gradients were constant during the time period before and after the electron increase. At about 2011.2 at a distance of 116.6 AU a second abrupt intensity increase was observed, this time for both electrons and protons. The increase for electrons was ~25% and occurred over a time period ~15 days or less. For >200 MeV protons the increase at this time was ~5% (unusually large) and occurred over a longer time period ~50 days. Between about 2011.2 and 2011.6, radial intensity gradients ~18%/AU and 3%/AU were observed for electrons and protons, respectively. These gradients were very similar to those observed for these particles before the 1st sudden increase at 2009.7. These large positive gradients observed after 2011.2 indicate that V1, although it has entered a different propagation region, is still within the overall heliospheric modulating region at least up to a time ~2011.6 (118.0 AU). In this paper we will discuss these events in more detail and consider possibilities for their explanation that have recently been suggested.

We have extended our earlier calculations of the distance to the Heliospheric Termination Shock (HTS) - which covered the period from the launch of V1 and V2 in 1977 to 2005 - to the period from 2006 to 2011. During this latter period the solar wind speed, ram pressure and magnetic field decreased to the lowest levels in recent history, related to the sunspot minimum in 2008-2009. The HTS distance has decreased correspondingly so that V1, which was crossed by the HTS at 94 AU in late 2004, would now, in early 2011, be expected to reach the HTS at a distance ~80 AU, when the HTS distance would be expected to be at its minimum. Similarly V2, which was crossed by the HTS at 84 AU in mid 2007, would, in early 2011, reach the HTS at a distance of only 74 AU. These distances, in early 2011, are ~15% less than those at which V1 and V2 initially reached the HTS. The distance to the Heliopause (HP) is more uncertain but recent calculations place its equilibrium distance at between 1.4-1.6 times the HTS distance. Allowing for an additional 1 year for the HP to reach its equilibrium minimum distance relative to the HTS would mean that, assuming this distance remains a constant fraction larger than the HTS distance, the HP distance would be at its minimum distance of (1.4-1.6) x 80 AU = 112-128 AU at V1 in early 2012. At this time V1 will be at a distance of ~120 AU so that there is a possibility that V1 could cross the HP and enter interstellar space at the time 2012.0 \pm 1 year. If the crossing does not happen during this time period, then it is unlikely that V1 will reach this defining boundary before about 2016 because of the expected outward motion of the HTS and the HP towards their more normal distances of 85-96 and ~120 -140 AU coincident with the maximum of the new sunspot cycle.

In an attempt to understand the source and propagation of galactic cosmic rays we have employed the Modified Weighted Slab technique along with recent values of the relevant cross sections to compute primary to secondary ratios including B/C and Sub-Fe/Fe for different galactic propagation models. The models that we have considered are the disk-halo diffusion model, the dynamical halo wind model, the turbulent diffusion model and a model with minimal reacceleration. The modified weighted slab technique will be briefly discussed and a more detailed description of the models will be given. We will also discuss the impact that the various models have on the problem of anisotropy at high energy and discuss what properties of a particular model bear on this issue.