Water in biological structures often displays non-Gaussian diffusion behavior. The objective of this study was to test the feasibility of non-Gaussian fitting by using the kurtosis model of the signal intensity decay curves obtained from DWI by using an extended range of b-values in studies of phantoms and HNSCC. Seventeen patients with HNSCC underwent DWI by using 6 b-factors (0, 50-1500 s/mm(2)) at 1.5T. Monoexponential (yielding ADC(mono)) and non-Gaussian kurtosis (yielding apparent diffusion coefficient D(app) and apparent kurtosis coefficient K(app)) fits were performed on a voxel-by-voxel basis in selected regions of interest (primary tumors, metastatic lymph nodes, and spinal cord). DWI studies were also performed on phantoms containing either water or homogenized asparagus. To determine whether the kurtosis model provided a significantly better fit than did the monoexponential model, an F test was performed. Spearman correlation coefficients were calculated to assess correlations between K(app) and D(app). The kurtosis model fit the experimental data points significantly better than did the monoexponential model (P < .05). D(app) was approximately twice the value of ADC(mono) (eg, in neck nodal metastases D(app) was 1.54 and ADC(mono) was 0.84). K(app) showed a weak Spearman correlation with D(app) in a homogenized asparagus phantom and for 44% of tumor lesions. The use of kurtosis modeling to fit DWI data acquired by using an extended b-value range in HNSCC is feasible and yields a significantly better fit of the data than does monoexponential modeling. It also provides an additional parameter, K(app), potentially with added value.