Hmm... Last estimate of V_rms applied to EOM can't be true.
If laser frequency noise is
Total frequency noise down to 1 kHz should then be 316 Hz_rms. As you've noted in your jupyter notebook the EOM frequency domain response to actuation signal goes as f, i.e.
where ms is EOM phase slope (15 mrad/V) and V_EOM is the applied actuation voltage. As you wrote, that leads to
If we assume that the EOM is taking 100% of the load for canceling laser frequency noise at a given frequency then it follows that the applied voltage to exactly cancel laser frequency noise is
This would indicate that the burden on the EOM becomes 1/f heavier in frequency because of the laser noise roll up and 1/f heaver because the EOM responds as f to signals in the Fourier domain. The above puts an upper bound on the ASD curve of load that the EOM is absorbing to cancel laser frequency noise.
Integrating the above V_EOM PSD down from high frequency will give the total V_rms load that the loop would need to apply to suppress laser frequency noise:
or
I've plotted this below and attached the notebook used for the calculation. It kind of looks like the maximum load at 1 kHz hits about 200 Vrms. Maybe I've gotten some factors wrong here but you can sort of see the scaling of how the maximum load on the EOM will look.
One thing to note is that the above estimates are an upper bound assuming that the EOM is taking all of the load down to that frequency point and that the PZT path isn't fighting or out of phase with EOM. To correctly compute the load on the EOM you are going to have to break down the EOM only portion of the loop from the laser frequency to the point of voltage injected into the EOM. This can be done by effectively nesting the PZT loop into the round trip gain in a way similar to that described in Josh Smith's Thesis section 2.6.2. Finding the actuation signal should be similar to finding the PLL actuation signal, at this point in the loop it is the G/(1G)/A copy of the sensor noise. In the high gain regime the applied EOM control signal should just be the laser frequency divided by the EOM frequency slope. Of course you can compute for G_EOM OLG to get a true value with a bunch of algebra.
Quote: 
I added some more graphs to our FSS analysis.
 I modeled AD602 variable gain amplifier with a circuit similar to how it actually works internally. The resultant circuit now has a right input resistance of 100 Ohm and capacitance of 2 pF. It is slightly noisier than AD602.
 I calculated the total suppression function and suppressed laser frequency noise.
 I did a parametric analysis by varying the PZT gain and obtained a range of cross over frequencies and phase margin we can attain.
 Everything looked awesome so far, but the real problem is EOM railing.
 So I calculated that if the said suppression is achieved, what is the actuation signal that is sent to EOM.
 The last plot is the spectral density of this actuation signal. Since only, we clearly are going to rail the EOM and the said suppression will never actually happen.
 We need to shift actuation region of EOM to even higher frequencies so that it doesn't rail. But at the same time, we need to keep the crossover frequency and phase margin of crossover decent. And to top it off, we need to do all this making sure we get at least suppression at 100 Hz.
 So from here, I guess we need to work on shifting the poles and zeros and tune in the gain values right to get point.
 I'm not 100% sure about this analysis, please let me know if I am doing something wrong here. The latest notebook and pdf are on git.
Quote: 
I'm attaching the first results of the transfer function from liso model of complete TTFSS box. I've also attached the jupyter notebook with some important formulas used.
Both files are present in git/cit_ctnlab/ctn_electronics/TTFSS_lisomodel/ git repo.
I'm just posting results for the analysis I did so far. I'll be able to make better inferences with some more work I intend to do tomorrow.
Edit: Wed Feb 6 16:59:59 2019
The updated plot is at git. Added analysis of variation of crossover frequency and phase margin with changing Fast gain.


