Hi, Taek-
Looking at the code, it appears that that FFT implementation is optimized for speed, so the input time series are zero-padded up to "powers of 2 combined with powers of 3 and 5 (no more than 3^3 and 5^5, though)." Note that hte output file's "TR" info has the appropriate physical frequency value (in Hz) corresponding to the upsampled time series. Do you need the spectral values at specific frequencies?
Note also that the type of the numbers in the output volumes is "complex", so if you wanted amplitudes, for example, you would perhaps want to use 3dcalc:
3dcalc -a FILE -expr 'abs(a)' -prefix AMPS.nii.gz
--pt