We don't want to advertise these parameters yet, since they will
likely be changed to support different loss functions.

Adding Hong et al. calculations for Coating Brownian Noise

See merge request !52

Anchal Gupta authored 4 years ago and Jameson Rollins committed 4 years ago

Modified getCoatBrownian() to accept different Bulk and Shear Loss
Angles and calculate brownian noise using Hong et a. PRD 87, 082001
(2013).

New Function description:
This function calculates Coating ThermoOptic noise using
Hong et al . PRD 87, 082001 (2013).
All references to 'the paper', 'Eq' adn 'Sec' are to this paper.

***Important Note***
Inside this function phi is used for denoting the phase shift suffered
by light in one way propagation through a layer. This is in conflict
with present nomenclature everywhere else where it is used as loss
angle.

The layers are assumed to be alernating low-n high-n layers, with
low-n first.
Inputs:
         f = frequency vector in Hz
       ifo = parameter struct from IFOmodel.m
     wBeam = beam radius (at 1 / e^2 power)
      dOpt = the optical thickness, normalized by lambda, of each
             coating layer.
New required arguments:
       mTi = Mirror Transmittance
        Ic = Circulating Laser Power falling on the Mirror (W)
The following new arguments should be made available in the Materials
object and lines 108-110 and 116-118 should be uncommented:
 lossBlown = Coating Bulk Loss Angle of Low Refractive Index layer
 lossSlown = Coating Shear Loss Angle of Low Refractive Index layer
lossBhighn = Coating Bulk Loss Angle of High Refractive Index layer
lossShighn = Coating Shear Loss Angle of High Refractive Index layer
   PETlown = Relevant component of Photoelastic Tensor of High n layer*
  PEThighn = Relevant component of Photoelastic Tensor of Low n layer*
Returns:
SbrZ = Brownian noise spectra for one mirror in m^2 / Hz

*
Default values of PETlown and PEThighn are chosen from sec. A.1. to get
the longitudnal coefficent of photoelasticity as -0.5 as been asserted
by the paper there for Tantala and -0.27 for Silica.
These values also need to be added in Materials object.

API change: load_budget function return separate Budget and IFO

Closes #69

See merge request !97

This addresses issues about confusion over the IFO being set as a class
attribute of the Budget, and is required to close #69.

newlines '\n' in the string provided to ipshell.ex cause it to choke.

Resolve "Mistake in suspension thermal noise calculation"

Closes #68

See merge request !96

Christopher Wipf authored 4 years ago and Jameson Rollins committed 4 years ago

The joint order in hForce/vForce must remain consistent between
this code and the thermal noise calculation

Also remove a couple of unused variables

Fixes #68

Remove carrier density noise from pygwinc

Closes #50

See merge request !90

Christopher Wipf authored 4 years ago and Jameson Rollins committed 4 years ago

Based on the new calculation with Debye screening from
arXiv:2003.05345, there appears to be no plausible scenario where
this noise will affect the design of GW interferometers.  Koji
estimated the noise is below 1e-28/rtHz for Voyager, as did Florian
Bruns (lead author of arXiv:2003.05345).  It's also negligible for
ET, as discussed in the paper.

Fixes #50

Update coating parameters for aLIGO and Aplus

See merge request !93

Nicholas Demos authored 4 years ago and Jameson Rollins committed 4 years ago

From MIT coating thermal noise experiment results:
* https://dcc.ligo.org/LIGO-G2000880
and LMA 2020 paper:
* https://iopscience.iop.org/article/10.1088/1361-6382/ab77e9

convert more types to float for easier serialization

was not matching some derived classes

also fix some associated error and log messages

convert IFO Structs to yaml for HDF5 storage