We present new kinetic-gasdynamic model of the solar wind interaction with the local interstellar medium. The model incorporates several processes suggested by McComas et al. (2009) for the origin of the heliospheric ENA ribbon -- the most prominent feature seen in the all sky maps of heliospheric ENAs discovered by the Interstellar Boundary Explorer (IBEX). The ribbon is a region of enhanced fluxes of ENAs crossing almost the entire sky. Soon after the ribbon's discovery it was realized (McComas et al., 2009) that the enhancement of the fluxes could be in the directions where the radial component of the interstellar magnetic field around the heliopause is close to zero (Schwadron et al., 2009). Our model includes secondary charge exchange of the interstellar H atoms with the interstellar pickup protons outside the heliopause and is a further advancement of the kinetic-gasdynamic model by Malama et al. (2006) where pickup protons were treated as a separate kinetic component. Izmodenov et al. (2009) have shown in the frame of Malama's model that the interstellar pickup protons outside the heliopause maybe a significant source of ENAs at energies above 1 keV. The difference between the current work and that of Izmodenov et al. (2009) is in the assumption of no-scattering for newly created pickup protons outside the heliopause. In this limit the model produces a feature qualitatively similar to the ribbon observed by IBEX.

We present a new model of the heliospheric interface - the region of the solar wind interaction with the local interstellar medium. This new model performs a multi-component treatment of charged particles in the heliosphere. All charged particles are divided into several co-moving types. The coldest type, with parameters typical of original solar wind protons, is considered in the framework of fluid approximation. The hot pickup proton components created from interstellar H atoms and heliospheric ENAs by charge exchange, electron impact ionization and photoionization are treated kinetically. The charged components are considered self-consistently with interstellar H atoms, which are described kinetically as well. To solve the kinetic equation for H atoms we use the Monte Carlo method with splitting of trajectories, which allows us 1) to reduce statistical uncertainties allowing correct interpretation of observational data, 2) to separate all H atoms in the heliosphere into several populations depending on the place of their birth and on the type of parent protons.