Computational simulations reproduce in vivo population receptive field mapping fMRI results

Abstract

Population receptive field (pRF) mapping is widely used to characterize retinotopic organization based on functional magnetic resonance imaging (fMRI) data. Despite its broad adoption, the factors governing intra- and inter-subject variability in pRF estimates remain incompletely understood, limiting the ability to evaluate and optimize visual stimulation paradigms prior to data collection. Here, we investigate whether large-scale simulations can reproduce in vivo run-to-run variability patterns observed in pRF mapping and provide mechanistic insight into their origins. With GEMSim-pRF, our newly proposed computational framework for large-scale simulation and estimation of pRF responses, we generated millions of synthetic fMRI time courses across a wide range of receptive field parameters and noise conditions. We analyzed the variability of pRF estimation results derived from simulations and compared them with in vivo data from the publicly available NYU Retinotopy Dataset. Here we show that our simulation results matched the characteristic eccentricity-dependent variability observed in empirical pRF estimates. These findings show that key variability patterns observed in empirical pRF mapping can be successfully reproduced in large-scale simulations, establishing simulation-based analysis as a practical approach for understanding, predicting, evaluating and ultimately improving the behaviour of retinotopic mapping paradigms before empirical data collection.