We analyze the non-gaussian density perturbations generated in ekpyrotic/cyclic models based on heterotic M-theory. In this picture, two scalar fields produce nearly scale-invariant entropic perturbations during an ekpyrotic phase that are converted into curvature modes after the ekpyrotic phase is complete and just before the big bang. Both intrinsic non-linearity in the entropy perturbation and the conversion process contribute to non-gaussianity. The range of the non-gaussianity parameter fNL depends on the details of the scalar field potential during the ekpyrotic phase, and on how gradual the conversion process is. Although a wider range is possible, in principle, natural values of the parameters of the potential combined with a gradual conversion process lead to values of −60 fNL +80, typically much greater than slow-roll inflation but within the current observational bounds.Ekpyrotic [1] and cyclic [2] models of the universe use quantum fluctuations of scalar fields produced during a slowly contracting phase with equation of state w > 1 to generate the observed nearly scale-invariant spectrum of curvature (energy density) fluctuations after the big bang. One mechanism considered for converting the fluctuations of a scalar field into cosmological curvature perturbations relies on higher-dimensional effects at the collision between orbifold planes (branes) along an extra dimension [3]. Recently, however, a new "entropic" mechanism [4,5] has been proposed that relies on two (or more) scalar fields and ordinary 4d physics, stimulating new approaches to ekpyrotic and cyclic cosmology that may not require branes or extra dimensions at all [6][7][8][9].In this paper, we wish to consider an important byproduct of the entropic mechanism: a non-gaussian contribution to the density fluctuation spectrum that is more than an order of magnitude greater than for conventional inflationary models and that can satisfy current observational bounds, including the recently claimed detection of non-gaussianity [10]. Our results differ significantly from the cases considered by Koyama et al. [11] and Buchbinder et al. [12] in which they assumed an ekpyrotic (w ≫ 1) phase that continues all the way to the conversion of entropic to curvature fluctuations, as in the "new ekpyrotic" model [6]; for these cases, the non-gaussianity is amplified by the non-linear evolution on super-horizon scales to the point where f N L , the parameter characterizing the non-linear curvature perturbation [13], reaches magnitude O(c 2 1 ), where c 1 ≈ 2 √ w + 1 parameterizes the steepness of the scalar field potential during the ekpyrotic phase. A potential problem is that a minimum of c 1 ≥ 10 is required just to satisfy the current upper bound constraints on n s , the spectral tilt of the scalar (energy density) perturbation spectrum [5]; and c 1 ≥ 30 is needed to reach the best-fit value n s ≈ 0.97. Yet, excluding finelytuned cancellations, the resulting value of f N L obtained in Ref. [6][7][8][9] is marginally consistent with current observationa...