Recent mutagenesis studies using the hydrophobic segment of Aβ suggest that aromatic Σ-stacking interactions may not be critical for fibril formation. We have tested this conjecture by probing the effect of Leu, Ile, and Ala mutation of the aromatic Phe residues at positions 19 and 20, on the double-layer hexametric chains of Aβ fragment Aβ16–22 using explicit solvent all-atom molecular dynamics. As these simulations rely on the accuracy of the utilized force fields, we first evaluated the dynamic and stability dependence on various force fields of small amyloid aggregates. These initial investigations led us to choose AMBER99SB-ILDN as force field in multiple long molecular dynamics simulations of 100 ns that probe the stability of the wild-type and mutants oligomers. Single-point and double-point mutants confirm that size and hydrophobicity are key for the aggregation and stability of the hydrophobic core region (Aβ16–22). This suggests as a venue for designing Aβ aggregation inhibitors the substitution of residues (especially, Phe 19 and 20) in the hydrophobic region (Aβ16–22) with natural and non-natural amino acids of similar size and hydrophobicity.