Turns out what was causing the instability in the aLIGO plots were the lock commands which I forgot to remove before running the simulation. Removing these also made the simulation much faster.

Other than that I improved other stuff in the simulations:

• The LO phase in the aPlus simulation is now optimized for the lowest noise at 100Hz.
• Added RF PDs diagnostics (see attachments 8 for aPlus and 9 for aLIGO). The thresholds (red dashed lines in attachments 8,9) for cutting marginal simulations are set such that roughly 30% of the simulations are discarded.
• Removed DHARD because it jacks up the RF PD readings in aPlus for some reason.
• Fixed the sign of laser frequency shift in response to CARM offset.

Still need to do:

• Incorporate Jon’s noise curves.
• Add phase noise for LO beam.
• Add arm reflectivity imbalance.
• Add mirror curvature imbalance.
• Include feedback loops using Pytickle.

Feel free to add to the todo list.