do you only want to show the effect of thermal noise in output spectral? if so, i think that you can add a time-domain random number scaled to the output of adc and perform the fft.
yes, you can calculate the noise of opamp firstly, and then generate the random number with the variance from the opamp's noise. the thermal noise is Gauss distribution, so it can be calculated by 3-σ rule.