We present two- and three-dimensional simulations of the growth of quantum fluctuations of the scalar fields in supersymmetric hybrid inflation models. For a natural range of couplings, sub-horizon quantum fluctuations undergo rapid growth due to scalar field dynamics, resulting in the formation of quasi-stable non-topological solitons (inflaton condensate lumps) which dominate the post-inflation era.

We demonstrate the existence of two-field Q-ball solutions of the scalar field equations of supersymmetric D- and F-term hybrid inflation. The solutions consist of a complex inflaton field together with a real symmetry breaking field. Such inflatonic Q-balls may play a fundamental role in reheating and the post-inflation era of supersymmetric hybrid inflation models.

It has been shown that hybrid inflation may end with the formation of non-topological solitons of inflaton field. As a first step towards a fully realistic picture of the post-inflation era and reheating in supersymmetric hybrid inflation models, we study the classical scalar field equations of a supersymmetric hybrid inflation model using a semi-analytical ansatz for the spatial dependence of the fields. Using the minimal D-term inflation model as an example, the inflaton field is evolved using the full 1-loop effective potential from the slow-rolling era to the U(1)_{FI} symmetry-breaking phase transition. Spatial perturbations of the inflaton corresponding to quantum fluctuations are introduced for the case where there is spatially coherent U(1)_{FI} symmetry breaking. The maximal growth of the dominant perturbation is found to depend only on the ratio of superpotential coupling \lambda to the gauge coupling g. The inflaton condensate fragments to non-topological solitons for \lambda/g > 0.09. Possible consequences of non-topological soliton formation in fully realistic SUSY hybrid inflation models are discussed.