Precise time-dependent measurements of the Z = 2 component in the cosmic radiation provide crucial information about the propagation of charged particles through the heliosphere. The PAMELA experiment, with its long flight duration (15th June 2006 - 23rd January 2016) and the low energy threshold (80 MeV/n) is an ideal detector for cosmic ray solar modulation studies. In this paper, the helium nuclei spectra measured by the PAMELA instrument from July 2006 to December 2009 over a Carrington rotation time basis are presented. A state-of-the-art three-dimensional model for cosmic-ray propagation inside the heliosphere was used to interpret the time-dependent measured fluxes. Proton-to-helium flux ratio time profiles at various rigidities are also presented in order to study any features which could result from the different masses and local interstellar spectra shapes.

The local interstellar spectra (LIS's) for galactic cosmic rays (CR's) cannot be directly observed at the Earth below certain energies, because of solar modulation in the heliosphere. With Voyager 1 crossing the heliopause in 2012, in situ experimental LIS data below 100 MeV/nuc can now constrain computed galactic CR spectra. Using galactic propagation models, galactic electron, proton and light nuclei spectra can be computed, now more reliably as very LIS's. Using the Voyager 1 observations made beyond the heliopause, and the observations made by the PAMELA experiment in Earth orbit for the 2009 solar minimum, as experimental constraints, we simultaneously reproduced the CR electron, proton, Helium and Carbon observations by implementing the GALPROP code. Below about 30 GeV/nuc solar modulation has a significant effect and a comprehensive three-dimensional (3D) numerical modulation model is used to compare the computed spectra with the observed PAMELA spectra possible at these energies. Subsequently the computed LIS's can be compared over as wide a range of energies as possible. The simultaneous calculation CR spectra with a single propagation model allows the LIS's for positrons, Boron and Oxygen to also be inferred. This implementation of the most comprehensive galactic propagation model (GALPROP), alongside a sophisticated solar modulation model to compute CR spectra for comparison with both Voyager 1 and PAMELA observations over a wide energy range, allows us to present new self-consistent very LIS's (and expressions) for electrons, positrons, protons, Helium, Carbon, Boron and Oxygen for the energy range of 3 MeV/nuc to 100 GeV/nuc.

Galactic cosmic rays (GCRs) are affected by solar modulation while they propagate through the heliosphere. The study of the time variation of GCR spectra observed at Earth can shed light on the underlying physical processes, specifically diffusion and particle drifts. Recently, the AMS-02 experiment measured with very high accuracy the time variation of the cosmic ray proton and helium flux between May 2011 and May 2017 in the rigidity range from 1 to 60 GV. In this work, a comprehensive three-dimensional (3D) steady-state numerical model is used to solve Parker's transport equation and is used to reproduce the monthly proton fluxes observed by AMS-02. We find that the rigidity slope of the perpendicular mean free path above 4 GV remains constant, while below 4 GV it increases during solar maximum. Assuming the same mean free paths for helium and protons, the models are able to reproduce the time behavior of the p/He ratio observed by AMS-02. The dependence of the diffusion tensor on the particle mass-to-charge ratio, A/Z, is found to be the main cause of the time dependence of p/He below 3 GV.

Heliospheric modulation conditions were unusually quite during the last solar minimum activity between Solar Cycles 23/24. Fortunately, the PAMELA space-experiment measured six-month averaged Galactic positron spectra for the period July 2006 to December 2009, over an energy range of 80 MeV to 30 GeV, which is important for solar modulation. The highest level of Galactic positrons was observed at Earth during the July-December 2009 period. A well-established, comprehensive three-dimensional (3D) numerical model is applied to study the modulation of the observed positron spectra. This model had been used previously to understand the modulation of Galactic protons and electrons also measured by PAMELA for the same period. First, a new very local interstellar spectrum for positrons is constructed, using the well-known GALPROP code together with the mentioned PAMELA observations. The 3D model is used to distinguish between the dominant mechanisms responsible for the heliospheric modulation of Galactic positrons, and to understand the effect of particle drift during this unusual minimum in particular, which is considered diffusion dominant, even though particle drift still had a significant role in modulating positrons. Lastly, the expected intensity of Galactic positrons during an A>0 polarity minimum, with similar heliospheric conditions than for 2006-2009, is predicted to be higher than what was observed by PAMELA for the 2006-2009 unusual minimum.

New results on the short-term galactic cosmic ray (GCR) intensity variation (Forbush decrease) in December 2006 measured by the PAMELA instrument are presented. Forbush decreases are sudden suppressions of the GCR intensities which are associated with the passage of interplanetary transients such as shocks and interplanetary coronal mass ejections (ICMEs). Most of the past measurements of this phenomenon were carried out with ground-based detectors such as neutron monitors or muon telescopes. These techniques allow only the indirect detection of the overall GCR intensity over an integrated energy range. For the first time, thanks to the unique features of the PAMELA magnetic spectrometer, the Forbush decrease commencing on 2006 December 14, following a CME at the Sun on 2006 December 13 was studied in a wide rigidity range (0.4 - 20 GV) and for different species of GCRs detected directly in space. The daily averaged GCR proton intensity was used to investigate the rigidity dependence of the amplitude and the recovery time of the Forbush decrease. Additionally, for the first time, the temporal variations in the helium and electron intensities during a Forbush decrease were studied. Interestingly, the temporal evolutions of the helium and proton intensities during the Forbush decrease were found in good agreement, while the low rigidity electrons (< 2 GV) displayed a faster recovery. This difference in the electron recovery is interpreted as a charge-sign dependence introduced by drift motions experienced by the GCRs during their propagation through the heliosphere.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

The PAMELA space experiment, in orbit since 2006, has measured cosmic rays through the most recent A < 0 solar minimum activity period. During this entire time, galactic electrons and protons have been detected down to 70 MV and 400 MV, respectively, and their differential intensity variation in time has been monitored with unprecedented accuracy. These observations are used to show how differently electrons and protons responded to the quiet modulation conditions that prevailed from 2006 to 2009. It is well known that particle drifts, as one of four major mechanisms for the solar modulation of cosmic rays, cause charge-sign dependent solar modulation. Solar minimum activity periods provide optimal conditions to study these drift effects. The observed behaviour is compared to the solutions of a three-dimensional model for cosmic rays in the heliosphere, including drifts. The numerical results confirm that the difference in the evolution of electron and proton spectra during the last prolonged solar minimum is attributed to a large extent to particle drifts. We therefore present new evidence of charge-sign dependent solar modulation, with perspective on its peculiarities for the observed period from 2006 to 2009.

Cosmic-ray electrons and positrons are a unique probe of the propagation of cosmic rays as well as of the nature and distribution of particle sources in our Galaxy. Recent measurements of these particles are challenging our basic understanding of the mechanisms of production, acceleration and propagation of cosmic rays. Particularly striking are the differences between the low energy results collected by the space-borne PAMELA and AMS-02 experiments and older measurements pointing to sign-charge dependence of the solar modulation of cosmic-ray spectra. The PAMELA experiment has been measuring the time variation of the positron and electron intensity at Earth from July 2006 to December 2015 covering the period for the minimum of solar cycle 23 (2006-2009) till the middle of the maximum of solar cycle 24, through the polarity reversal of the heliospheric magnetic field which took place between 2013 and 2014. The positron to electron ratio measured in this time period clearly shows a sign-charge dependence of the solar modulation introduced by particle drifts. These results provide the first clear and continuous observation of how drift effects on solar modulation have unfolded with time from solar minimum to solar maximum and their dependence on the particle rigidity and the cyclic polarity of the solar magnetic field.

With the cosmic ray observations made by the Voyager 1 spacecraft outside the dominant modulating influence of the heliosphere, the comparison of computed galactic spectra with experimental data at lower energies is finally possible. Spectra for specifically protons, Helium and Carbon nuclei, computed by galactic propagation models, can now be compared with observations at low energies from Voyager 1 and at high energies from the PAMELA space detector at Earth. We set out to reproduce the Voyager 1 observations in the energy range of 6 MeV/nuc to 60 MeV/nuc, and the PAMELA spectrum above 50 GeV/nuc, using the GALPROP code, similarly to our previous study for Voyager 1 electrons. By varying the galactic diffusion parameters in the GALPROP plain diffusion model, specifically the rigidity dependence of spatial diffusion, and then including reacceleration, we compute spectra simultaneously for galactic protons, Helium and Carbon.We present new local interstellar spectra, with expressions for the energy range of 3 MeV/nuc to 100 GeV/nuc, which should be of value for solar modulation modeling.

Precision measurements of the electron component in the cosmic radiation provide important information about the origin and propagation of cosmic rays in the Galaxy not accessible from the study of the cosmic-ray nuclear components due to their differing diffusion and energy-loss processes. However, when measured near Earth, the effects of propagation and modulation of galactic cosmic rays in the heliosphere, particularly significant for energies up to at least 30 GeV, must be properly taken into account. In this paper the electron (e^-) spectra measured by PAMELA down to 70 MeV from July 2006 to December 2009 over six-months time intervals are presented. Fluxes are compared with a state-of-the-art three-dimensional model of solar modulation that reproduces the observations remarkably well.

We argue that the degree of understanding the causes and mechanisms of the long-term variations (11-year and 22-year) in the galactic cosmic rays (GCR) characteristics is still insufficient and to improve it we need new approaches and methods. For the time being there is a long-lasting controversy on how these long-term variations, observed for more than 50 years in the inner heliosphere, are formed. It is widely believed that the 11-year variation is due entirely to the toroidal branch of solar activity (the area and number of sunspots, the strength of the heliospheric magnetic field etc) because of the diffusion, convection and adiabatic energy loss, while the much smaller 22-year variation is caused by the particle drifts connected with the poloidal branch of solar activity (the high-latitude solar magnetic fields). At the same time, both past and more recent numerical simulations indicate that the contribution of particle drifts could be significant for both 22- and 11-year variations in the GCR intensities. However, even those who agree on the significant influence of drifts appear to have different perceptions on the mechanisms of this influence. In this paper, we present an analysis of the possible causes of the first point of view (small role of drifts in the 11-year GCR variation) and the reasons why one can expect the significant contri- bution of the processes connected with the poloidal branch of solar activity in both types of the long-term variations of the GCR characteristics. Then we briefly discuss some numerical methods suggested in the past and recently and the approaches and perspectives for the sought-for methods are considered.

The local interstellar spectra (LIS) for cosmic rays at energies below about 30 GeV/nuc are increasingly obscured from view at Earth by solar modulation, the lower the energy. These charged particles encounter significant changes in the heliosphere, over an 11-year cycle, which include processes such as convection, diffusion, adiabatic energy losses and gradient, curvature and current sheet drifts. Particle drifts cause chargesign dependent modulation and a 22-year cycle, adding complexity to determining the respective very LIS from observations only at Earth. However, with measurements now made by the Voyager 1 spacecraft in the vicinity of the heliopause, it is possible to determine a very LIS for galactic electrons between about 5 MeV and about 120 MeV. At these low energies also galactic protons observed in the outer heliosphere had been completely obscured by the so-call anomalous component which is accelerated inside the heliosheath. Since August 2012, these anomalous cosmic rays are substantially depleted at Voyager 1 so that also for cosmic ray ions it is now possible to obtain a lower limit to their very LIS. Combining numerical modelling of solar modulation with the accurate measurements by the PAMELA mission and with Voyager observations, the lower limit of the very LIS for electrons, protons and helium, and other ions, can be determined from about 5 MeV and above. These spectra are called heliopause spectra, considered to be the lowest possible very LIS. Also from an astrophysics point of view, the determination of what can be called a very LIS, not just an averaged galactic spectrum, is encouraging. The mentioned aspects are discussed, focussing on a comparison of recent heliospheric observations and corresponding solar modulation modelling.

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a time-dependent modulation model. By introducing recent theoretical advances in the transport coefficients in the model, computed intensities are compared with Voyager 1, International Monitoring Platform (IMP) 8, and Ulysses proton observations in search of compatibility. The effect of different modulation parameters on computed intensities is also illustrated. It is shown that this approach produces, on a global scale, realistic cosmic-ray proton intensities along the Voyager 1 spacecraft trajectory and at Earth upto ~2004, whereafter the computed intensities recovers much slower towards solar minimum than observed in the inner heliosphere. A modified time dependence in the diffusion coefficients is proposed to improve compatibility with the observations at Earth after ~2004. This modified time dependence led to an improved compatibility between computed intensities and the observations along the Voyager 1 trajectory and at Earth even after ~2004. An interesting result is that the cosmic-ray modulation during the present polarity cycle is not determined only by changes in the drift coefficient and tilt angle of the wavy current sheet, but is also largely dependent on changes in the diffusion coefficients.

A two-dimensional, time-dependent numerical model is used to calculate the modulation of cosmic rays in the heliosphere. Computations are compared to spacecraft observations in the inner and outer heliosphere. It is shown that the model produces cosmic ray proton intensities compatible to different spacecraft observations on a global scale, at Earth and along both Voyager spacecraft trajectories. The study reveals that when the same modulation parameters, which resulted in compatible intensities along Voyager 1, were assumed along the Voyager 2 trajectory, the model failed to reproduce the observations. The study also found that any change in diffusion parameters alone could not reproduce the cosmic ray observations along Voyager 2 so that changes to the heliospheric geometry were necessary i.e the computed intensities along both Voyager trajectories suggest that the heliosphere is asymmetric. Furthermore, $E >$ 70 MeV and 133-242 MeV proton intensities along Voyager 1 and 2 trajectories are predicted from end of 2012 onwards. It is shown that the computed intensities along Voyager 1 increase with an almost constant rate up to the heliopause. However, the model shows that Voyager 2 is still under the influence of temporal solar activity changes because of its relatively large distance to the heliopause. Along the Voyager 2 trajectory, the intensities remained generally constant for some time and should soon start to increase steadily.

Given that the bulk of the galactic cosmic rays (CRs) originates from transient point-like sources, such as supernova remnants, the flux of the CR primary component measured at Earth depends on the local source history. Whereas the secondary omponent shows little to no variation due to nearby sources, and instead, should depend on the global distribution of CRs. Steady-state, rotational symmetric models (2D) of CR propagation assume smeared-out CR sources in the Galaxy and cannot properly take into account the influence of nearby point sources. In this paper, we infer evidence of nearby sources by illustrating that a 2D propagation model will not describe CR primaries and secondaries equally well. We adapted the 2D version of the GALPROP code to a compute-cluster environment and perform parameter studies comparing CR spectra with primary and secondary CR data separately. Doing a parameter study, one may expect different best fit values looking at the primary and secondary CR components separately, as it is unlikely that the source history mimicked by the 2D models coincides with the real local source history. We find that the primaries and secondaries are fitted with differing best fit parameters, with the results indicating that the 2D model being more suited to modelling secondary CRs. The insufficient modelling of primary CRs can be contributed to the possible presence of local CR sources.

The transport environment for particles in the heliosphere, e.g. galactic cosmic rays (GCRs) and MeV electrons (including those originating from Jupiters magnetosphere), is defined by the solar wind flow and the structure of the embedded heliospheric magnetic field. Solar wind structures, such as co-rotating interaction regions (CIR), can result in periodically modulation of both particles species. A detailed analysis of this recurrent Jovian electron events and galactic cosmic ray decreases measured by SOHO EPHIN is presented here, showing clearly a change of phase between both phenomena during the cause of the years 2007 and 2008. This effect can be explained by the change of difference in heliolongitude between the Earth and Jupiter, which is of central importance for the propagation of Jovian electrons. Furthermore, the data can be ordered such that the 27-day Jovian electron variation vanishes in the sector which does not connect the Earth with Jupiter magnetically using observed solar wind speeds.

The recent solar minimum activity period and the consequent minimum modulation conditions for cosmic rays were unusual compared to previous solar minimum periods. The highest spectra of galactic protons (and other cosmic rays) were recorded by the PAMELA instrument at Earth in late 2009, in contrast to expectations. The spectrum, between 100 MeV and 50 GeV, for December 2009 is compared to proton spectra observed during previous solar minimum periods, back to 1965. Corresponding numerical modeling is presented which predicts that the next solar minimum spectra could even be higher if similar modulation conditions then would occur as in 2008-2009. The reason is that incorporating gradient and curvature drifts in modulation models causes proton spectra for A > 0 solar magnetic cycles (e.g., around 1976, 1997) to always be higher than during A < 0 cycles (e.g. around 1965, 1987, 2009) at energies below a few GeV, if the same modulation conditions would prevail.

We present a heliopause spectrum for galactic electrons over an energy range from 1 MeV to 50 GeV at 122 AU from the Sun. This spectrum can be considered the lowest possible local interstellar spectrum (LIS). This is obtained by using a comprehensive numerical model for solar modulation in comparison with Voyager 1 observations during 2010 and PAMELA data at Earth. Below about 1 GeV, this LIS exhibits a power law with E^(-(1.5 plus minus 0.1), where E is the kinetic energy of these electrons. Reproducing the PAMELA electron spectrum for 2009 requires a LIS with a different power law above about 5 GeV of the form E^(-3.15 plus minus 0.05). The heliopause spectrum thus consists of two power laws with a break between about 800 MeV and about 2 GeV.

Observationally, it is difficult to establish at what rigidity the modulation of galactic cosmic rays (CRs) actually begins in the heliosphere. It should be possible to do if the relevant local interstellar CR spectra were known and reliable measurements were made between 10 GV and about 200 GV, inside and outside the heliosphere. Numerical models for solar modulation studies are based on simply assuming that CR modulation begins at a given spatial boundary and at rigidities between 30 to 50 GV, usually handled as an initial condition. The Stochastic Differential Equations approach to numerical modelling presents the opportunity to determine the level of modulation from high to low rigidities because an initial condition at a prescribed high rigidity is not required. We present the results of such an approach and show the percentage modulation of CR protons as a function of kinetic energy between 100 MeV and 250 GeV.

We present a new local interstellar spectrum for cosmic ray electrons over an energy range from 1 MeV to 70 GeV. Below (0.8 plus minus 0.2) GeV it has a power law form, E^-(1.55 plus minus 0.05), with E kinetic energy, which is consistent with previous studies. This is derived from comparing Voyager 1 electron data observed during 2010 with a comprehensive modulation model. However, reproducing the PAMELA electron spectrum observed at Earth during late 2006 and to address an unexpected increase in the electron spectrum between about 2 GeV and up to about 20 GeV, a spectral index of -(3.3 plus minus 0.1) instead of the reported -(3.18 plus minus 0.05) was found for this energy range. This feature cannot be caused by solar modulation or any other process inside the heliosphere.

This is an overview of the solar modulation of cosmic rays in the heliosphere. It is a broad topic with numerous intriguing aspects so that a research framework has to be chosen to concentrate on. The review focuses on the basic paradigms and departure points without presenting advanced theoretical or observational details for which there exists a large number of comprehensive reviews. Instead, emphasis is placed on numerical modeling which has played an increasingly significant role as computational resources have become more abundant. A main theme is the progress that has been made over the years. The emphasis is on the global features of CR modulation and on the causes of the observed 11-year and 22-year cycles and charge-sign dependent modulation. Illustrative examples of some of the theoretical and observational milestones are presented, without attempting to review all details or every contribution made in this field of research. Controversial aspects are discussed where appro- priate, with accompanying challenges and future prospects. The year 2012 was the centennial celebration of the discovery of cosmic rays so that several general reviews were dedicated to historical aspects so that such developments are brie y presented only in a few cases.

A heliopause spectrum at 122 AU from the Sun is presented for galactic electrons over an energy range from 1 MeV to 50 GeV that can be considered the lowest possible local interstellar spectrum (LIS). The focus is on the spectral shape of the LIS below 1.0 GeV. The study is done by using a comprehensive numerical model for solar modulation in comparison with Voyager 1 observations at 110 AU from the Sun and PAMELA data at Earth. Below 1.0 GeV, this LIS exhibits a power law,E to the power -(1.55+-0.05), where E is the kinetic energy. However, reproducing the PAMELA electron spectrum averaged for 2009, requires a LIS with a different power law of the form E to the power -(3.15+-0.05) above about 5 GeV. Combining the two power laws with a smooth transition from low to high energies yields a LIS over the full energy range that is relevant and applicable to the modulation of cosmic ray electrons in the heliosphere. The break occurs between 800 MeV and 2 GeV as a characteristic feature of this LIS.

The last solar minimum activity period, and the consequent minimum modulation conditions for cosmic rays, was unusual. The highest levels of galactic protons were recorded at Earth in late 2009 in contrast to expectations. Proton spectra observed for 2006 to 2009 from the PAMELA cosmic ray detector on-board the Resurs-DK1 satellite are presented together with the solutions of a comprehensive numerical model for the solar modulation of cosmic rays. The model is used to determine what mechanisms were mainly responsible for the modulation of protons during this period, and why the observed spectrum for 2009 was the highest ever recorded. From mid-2006 until December 2009 we find that the spectra became significantly softer because increasingly more low energy protons had reached Earth. To simulate this effect, the rigidity dependence of the diffusion coefficients had to decrease significantly below ~3 GeV. The modulation minimum period of 2009 can thus be described as relatively more "diffusion dominated" than previous solar minima. However, we illustrate that drifts still had played a significant role but that the observable modulation effects were not as well correlated with the waviness of the heliospheric current sheet as before. Protons still experienced global gradient and curvature drifts as the heliospheric magnetic field had decreased significantly until the end of 2009, in contrast to the moderate decreases observed during previous minimum periods. We conclude that all modulation processes contributed to the observed increases in the proton spectra for this period, exhibiting an intriguing interplay of these major mechanisms.

The energy spectra of galactic cosmic rays carry fundamental information regarding their origin and propagation. These spectra, when measured near Earth, are significantly affected by the solar magnetic field. A comprehensive description of the cosmic radiation must therefore include the transport and modulation of cosmic rays inside the heliosphere. During the end of the last decade the Sun underwent a peculiarly long quiet phase well suited to study modulation processes. In this paper we present proton spectra measured from July 2006 to December 2009 by PAMELA. The large collected statistics of protons allowed the time variation to be followed on a nearly monthly basis down to 400 MV. Data are compared with a state-of-the-art three-dimensional model of solar modulation.

Geminga and B0656+14 are the closest pulsars with characteristic ages in the ran ge of 100 kyr to 1 Myr. They both have spindown powers of the order 3e34 erg/s at present. The winds of these pulsars had most probably powered pulsar wind nebulae (PWNe) that broke up less than about 100 kyr after the birth of the pulsars. Assuming that leptonic particles accelerated by the pulsars were confined in th e PWNe and were released into the interstellar medium (ISM) on breakup of the PW Ne, we show that, depending on the pulsar parameters, both pulsars make a non-ne gligible contribution to the local cosmic ray (CR) positron spectrum, and they m ay be the main contributors above several GeV. The relatively small angular dist ance between Geminga and B0656+14 thus implies an anisotropy in the local CR po sitron flux at these energies. We calculate the contribution of these pulsars to the locally observed CR electr on and positron spectra depending on the pulsar birth period and the magnitude o f the local CR diffusion coefficient. We further give an estimate of the expecte d anisotropy in the local CR positron flux. Our calculations show that within the framework of our model, the local CR posit ron spectrum imposes constraints on pulsar parameters for Geminga and B0656+14, notably the pulsar period at birth, and also the local interstellar diffusion co efficient for CR leptons.

High-energy collisions of cosmic-ray nuclei with interstellar gas are believed to be the mechanism producing the majority of cosmic ray antiprotons. Due to the kinematics of the process they are created with a nonzero momentum; the characteristic spectral shape with a maximum at ~2 GeV and a sharp decrease towards lower energies makes antiprotons a unique probe of models for particle propagation in the Galaxy and modulation in the heliosphere. On the other hand, accurate calculation of the secondary antiproton flux provides a ``background'' for searches for exotic signals from the annihilation of supersymmetric particles and primordial black hole evaporation. Recently new data with large statistics on both low and high energy antiproton fluxes have become available which allow such tests to be performed. We use our propagation code GALPROP to calculate interstellar cosmic-ray propagation for a variety of models. We show that there is no simple model capable of accurately describing the whole variety of data: boron/carbon and sub-iron/iron ratios, spectra of protons, helium, antiprotons, positrons, electrons, and diffuse gamma rays. We find that only a model with a break in the diffusion coefficient plus convection can reproduce measurements of cosmic-ray species, and the reproduction of primaries (p, He) can be further improved by introducing a break in the primary injection spectra. For our best-fit model we make predictions of proton and antiproton fluxes near the Earth for different modulation levels and magnetic polarity using a steady-state drift model of propagation in the heliosphere.

High energy collisions of cosmic ray (CR) nuclei with interstellar gas are believed to be the mechanism producing the majority of CR antiprotons. The distinguishing spectral shape with a maximum at 2 GeV and a sharp decrease towards lower energies makes antiprotons a unique probe of the models of particle propagation in the Galaxy and modulation in the heliosphere. Besides, accurate calculation of the secondary antiproton flux provides a ``background'' for searches for exotic signals from the annihilation of supersymmetric particles and primordial black hole evaporation. Recently new data with large statistics on the antiproton flux have become available which allow for such tests to be performed. We use our 3D Galactic cosmic ray propagation code GALPROP to calculate interstellar propagation in several models. For our best model we make predictions of proton and antiproton fluxes near the Earth for different modulation levels and polarity using a steady-state drift model for heliospheric modulation.