A nonlinear BOLD model accounting for refractory effect by applying the longitudinal relaxation in NMR to the linear BOLD model

A mathematical model to regress the nonlinear blood oxygen level-dependent (BOLD) fMRI signal has been developed by incorporating the refractory effect into the linear BOLD model of the biphasic gamma variate function. The refractory effect was modeled as a relaxation of two separate BOLD capacities corresponding to the biphasic components of the BOLD signal in analogy with longitudinal relaxation of magnetization in NMR. When tested with the published fMRI data of finger tapping, the nonlinear BOLD model with the refractory effect reproduced the nonlinear BOLD effects such as reduced poststimulus undershoot and saddle pattern in a prolonged stimulation as well as the reduced BOLD signal for repetitive stimulation.