We present a large suite of cosmological simulations, the FORGE (F-of-R Gravity Emulator) simulation suite, which is designed to build accurate emulators for cosmological observables in galaxy clustering, weak gravitational lensing and galaxy clusters, for the 𝑓 (𝑅) gravity model. The total of 200 simulations explore the cosmological parameter space around the Planck Collaboration et al. ( 2020) cosmology with a Latin hypercube, for 50 combinations of fR0 , Ω m , 𝜎 8 and ℎ with all other parameters fixed. For each parameter combination, or node, we ran four independent simulations, one pair using 1024 3 particles in 500 ℎ −1 Mpc simulation boxes to cover small scales, and another pair using 512 3 simulation particles in 1.5 ℎ −1 Gpc boxes for larger scales. Each pair of initial conditions are selected such that sample variance on large scales is minimised on average. In this work we present an accurate emulator for the matter power spectrum in 𝑓 (𝑅)-gravity trained on FORGE. We have verified, using the crossvalidation technique, that the emulator accuracy is better than 2.5% for the majority of nodes, particularly around the center of the explored parameter space, up to scales of 𝑘 = 10 ℎMpc −1 . We have also checked the power spectrum emulator against simulations which are not part of our training set and found excellent agreement. Due to its high accuracy on small scales, the FORGE matter power spectrum emulator is well suited for weak lensing analysis and can play a key tool in constraining 𝑓 (𝑅)-gravity using current and future observational data.