Synthetic, self-oscillating models of the human vocal folds are used to study the complex and inter-related flow, structure, and acoustical aspects of voice production. The vocal folds typically collide during each cycle, thereby creating a brief period of glottal closure that has important implications for flow, acoustic, and motion-related outcomes. Many previous synthetic models, however, have been limited by incomplete glottal closure during vibration. In this study, a low-fidelity, two-dimensional, multilayer finite element model of vocal fold flow-induced vibration was coupled with a custom genetic algorithm optimization code to determine geometric and material characteristics that would be expected to yield physiologically-realistic frequency and closed quotient values. The optimization process yielded computational models that vibrated with favorable frequency and closed quotient characteristics. A tradeoff was observed between frequency and closed quotient. A synthetic, self-oscillating vocal fold model with geometric and material properties informed by the simulation outcomes was fabricated and tested for onset pressure, oscillation frequency, and closed quotient. The synthetic model successfully vibrated at a realistic frequency and exhibited a nonzero closed quotient. The methodology described in this study provides potential direction for fabricating synthetic models using isotropic silicone materials that can be designed to vibrate with physiologically-realistic frequencies and closed quotient values. The results also show the potential for a low-fidelity model optimization approach to be used to tune synthetic vocal fold model characteristics for specific vibratory outcomes.