Fibrillary protein aggregates rich in β-sheet structure have been implicated in the pathology of several neurodegenerative diseases. In this work, we investigate the formation of fibrils by performing discontinuous molecular dynamics simulations on systems containing 12 to 96 model Ac-KA 14 K-NH 2 peptides using our newly-developed off-lattice, implicit-solvent, intermediateresolution model, PRIME. We find that at a low concentration, random-coil peptides assemble into α-helices at low temperatures. At intermediate concentrations, random-coil peptides assemble into α-helices at low temperatures and large β-sheet structures at high temperatures. At high concentrations, the system forms β-sheets over a wide range of temperatures. These assemble into fibrils above a critical temperature which decreases with concentration and exceeds the isolated peptide's folding temperature. At very high temperatures and all concentrations, the system is in a random-coil state. All of these results are in good qualitative agreement with those by Blondelle and coworkers on Ac-KA 14 K-NH 2 peptides. The fibrils observed in our simulations mimic the structural characteristics observed in experiments in terms of the number of sheets formed, the values of the intra-and inter-sheet separations, and the parallel peptide arrangement within each β-sheet. Finally, we find that when the strength of the hydrophobic interaction between non-polar sidechains is high compared to the strength of hydrogen bonding, amorphous aggregates, rather than fibrillar aggregates, are formed.