The possibility of antihelium production in interaction of cosmic rays with the interstellar gas is studied using large-scale Monte Carlo simulations. For this purpose, an energy-dependent coalescence mechanism developed previously is extended to estimate the production of light antinuclei (${}^3\overline{\text{He}}$ and ${}^4\overline{\text{He}}$). The uncertainty in the coalescence parameter and its effect on the expected antiparticle flux is also investigated. The simulated background antihelium fluxes are found to be lower than the fluxes predicted by simplified models using numerical scaling techniques.

The study of the cosmic-ray deuteron and antideuteron flux receives an increasing interest in current astrophysics investigations. For both cases an important contribution is expected from the nuclear interactions of primary cosmic rays with intergalactic matter. In this work, deuteron and antideuteron production from 20 to 2.6$\times$10$^{7}$ GeV beam energy in p+p and p+A collisions were simulated using EPOS-LHC and Geant4's FTFP-BERT Monte Carlo models by adding an event-by-event coalescence model afterburner. These estimates depend on a single parameter ($p_0$) obtained from a fit to the data. The $p_0$ for deuterons in this wide energy range was evaluated for the first time. It was found that $p_0$ for antideuterons is not a constant at all energies as previous works suggested and as a consequence the antideuteron production cross section can be at least 20 times smaller in the low collision energy region, than earlier estimations.