We study numerically the decay of the inflaton by solving the full non-linear equations of motion on the lattice. We confirm that parametric resonance is effective in transferring energy from the inflaton to a scalar field as long as the self-interactions of the second field are very small. However, in the very broad resonance case (q>>1) the decay rate is limited by scatterings, which significantly slows down the decay. We also find that the inflaton cannot decay via parametric resonance into a scalar field with moderate self-interactions. This means that the preheating stage may be completely absent in many natural inflationary models.

We study the decay of the inflaton in a renormalizable two scalar theory. Since the dynamics of the system is dominated by states with large occupation numbers which admit a semiclassical description, the decay can be studied by solving the classical equations of motion on the lattice. Of particular interest is the case when the cross-coupling between the inflaton and the second scalar field is negative, which is naturally allowed in many realistic models. While the inflaton decays via parametric resonance in the positive coupling case we find that for negative coupling there is a new mechanism of particle production which we call negative coupling instability. Due to this new mechanism the variances of the fields grow significantly larger before the production is shut off by the backreaction of the created particles. We also find that heavy particles are produced much more efficiently with negative coupling, which is of prime importance for GUT baryogenesis. Using a simple toy model for baryogenesis and the results of our lattice simulations we show that for natural values of the cross-coupling enough 10^{14}GeV bosons are created to produce a baryon to entropy ratio consistent with observation. For positive coupling the value of the cross-coupling required to produce such massive particles is unnaturally large. In addition to our numerical results we obtain analytical estimates for the maximum variances of the fields in an expanding universe for all cases of interest in our model.

Neutral hydrogen mapping of 70 dwarf irregular (Sm, Im and BCD) galaxies is reported, with position-velocity contour maps presented for those resolved by the Arecibo beam. The galaxies were selected either from the Virgo Cluster Catalog, from similarly identified field galaxies, or from a distance-limited sample within the Arecibo declination range. We do not find any isolated dwarfs with a larger HI to optical radius ratio than DDO 154; the ``protogalaxy'' HI 1225+01 (Giovanelli \& Haynes 1989) continues to be a unique object among dwarfs that have been mapped in HI. For all dwarfs with significant rotation we are able to determine the sense of the spin. For a number of better-resolved dwarfs we are able to determine rotation curves, in most cases extending well beyond the last measured point in available synthesis array maps. Correlations among the several measures of galaxy size and mass are studied; in Paper II (Salpeter \& Hoffman 1996) we combine these data with those for the set of all available mapped dwarf irregular galaxies and for mapped spirals spanning a similar range of redshifts to investigate variations in Tully-Fisher relations and in surface densities as functions of galaxy size and luminosity or mass.