Neutralization and focusing of intense charged particle beam pulses by electrons forms the basis for a wide range of applications to high energy accelerators and colliders, heavy ion fusion, and astrophysics. For example, for ballistic propagation of intense ion beam pulses, background plasma can be used to effectively neutralize the beam charge and current, so that the self-electric and self-magnetic fields do not affect the ballistic propagation of the beam. From the practical perspective of designing advanced plasma sources for beam neutralization, a robust theory should be able to predict the self-electric and self-magnetic fields during beam propagation through the background plasma. The major scaling relations for the self-electric and self-magnetic fields of intense ion charge bunches propagating through background plasma have been determined taking into account the effects of transients during beam entry into the plasma, the excitation of collective plasma waves, the effects of gas ionization, finite electron temperature, and applied solenoidal and dipole magnetic fields. Accounting for plasma production by gas ionization yields a larger self-magnetic field of the ion beam compared to the case without ionization, and a wake of current density and self-magnetic field perturbations is generated behind the beam pulse. A solenoidal magnetic field can be applied for controlling the beam propagation. Making use of theoretical models and advanced numerical simulations, it is shown that even a small applied magnetic field of about 100G can strongly affect the beam neutralization. It has also been demonstrated that in the presence of an applied magnetic field the ion beam pulse can excite large-amplitude whistler waves, thereby producing a complex structure of self-electric and self-magnetic fields. The presence of an applied solenoidal magnetic field may also cause a strong enhancement of the radial self-electric field of the beam pulse propagating through the background plasma. If controlled, this physical effect can be used for optimized beam transport over long distances.