I did TF, input impedance, opamp stability and noise analyses of the complete circuit with feedback part, modeled AOM (by phase shifter) and modeled Cavity (by LPF+NOT). Attached are the results. In the process, I completed wrapping some more LISO functions in python. We can do calculate input referred noise, input impedance, opamp stability and plot only dominant noise sources. The circuit schematic is attached and for now uses OP27. Next steps are to figure out what adjustments (components and values) can lower the noise. Opamps are in the stable region only.

r R14 510 OUT n1
c C7 1u n1 n2
c C4 470p n1 n4
r R15 560k n2 n4
r R16 560k n3 gnd
c C6 1u n3 gnd
op ISS2 op27 n3 n2 n4

r R18 1k n4 n5
c C8 4.7n n5 n7