Residual gas sub-budgets and updates

Adds residual gas damping and additional molecular species to the residual gas calculations. The gas scattering and gas damping are now tracked separately for each species in an excess gas sub-budget. The damping calculation is done in the infinite volume limit by default but adds corrections for squeezed film damping if specified in the yaml file. The aLIGO excess gas sub-budget is shown below as an example.

The new yaml file fields are explained in this comment. The yaml files in this merge request are simpler and none of the interferometers include squeezed film damping.

For aLIGO, Aplus, and Voyager, H2, N2, and H2O are considered. The numbers for aLIGO are taken from the link in this comment. The numbers for Aplus and Voyager are copied from aLIGO and should be checked. The Cosmic Explorer interferometers additionally consider CO2 and Ar.

Closes #63 (closed).

gwinc/ifo/noises.py

@@ -833,6 +833,42 @@ class ExcessGasScatteringH2O(nb.Noise):
         return n * 2 / sinc_sqr
 
 
+class ExcessGasScatteringCO2(nb.Noise):
+    """Excess gas scattering for CO2
+
+    """
+    style = dict(
+        label='CO2 scattering',
+        color='xkcd:slate blue'
+    )
+
+    def calc(self):
+        cavity = arm_cavity(self.ifo)
+        species = self.ifo.Infrastructure.ResidualGas.CO2
+        n = noise.residualgas.residual_gas_scattering(
+            self.freq, self.ifo, cavity, species)
+        dhdl_sqr, sinc_sqr = dhdl(self.freq, self.ifo.Infrastructure.Length)
+        return n * 2 / sinc_sqr
+
+
+class ExcessGasScatteringAr(nb.Noise):
+    """Excess gas scattering for Ar
+
+    """
+    style = dict(
+        label='Ar scattering',
+        color='xkcd:violet'
+    )
+
+    def calc(self):
+        cavity = arm_cavity(self.ifo)
+        species = self.ifo.Infrastructure.ResidualGas.Ar
+        n = noise.residualgas.residual_gas_scattering(
+            self.freq, self.ifo, cavity, species)
+        dhdl_sqr, sinc_sqr = dhdl(self.freq, self.ifo.Infrastructure.Length)
+        return n * 2 / sinc_sqr
+
+
 class ExcessGasDampingH2(nb.Noise):
     """Residual gas damping for H2
 
@@ -865,9 +901,6 @@ class ExcessGasDampingN2(nb.Noise):
     @nb.precomp(sustf=precomp_suspension)
     def calc(self, sustf):
         species = self.ifo.Infrastructure.ResidualGas.N2
-        # n = noise.residualgas.residual_gas_damping(
-        #     self.freq, self.ifo, species, sustf)
-        # return 4 * n
         return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
 
 class ExcessGasDampingH2O(nb.Noise):
@@ -883,12 +916,40 @@ class ExcessGasDampingH2O(nb.Noise):
     @nb.precomp(sustf=precomp_suspension)
     def calc(self, sustf):
         species = self.ifo.Infrastructure.ResidualGas.H2O
-        # n = noise.residualgas.residual_gas_damping(
-        #     self.freq, self.ifo, species, sustf)
-        # return 4 * n
         return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
 
 
+class ExcessGasDampingCO2(nb.Noise):
+    """Excess gas damping for CO2
+
+    """
+    style = dict(
+        label='CO2 damping',
+        color='xkcd:slate blue',
+        linestyle='--',
+    )
+
+    @nb.precomp(sustf=precomp_suspension)
+    def calc(self, sustf):
+        species = self.ifo.Infrastructure.ResidualGas.CO2
+        return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
+
+
+class ExcessGasDampingAr(nb.Noise):
+    """Excess gas damping for Ar
+
+    """
+    style = dict(
+        label='Ar damping',
+        color='xkcd:violet',
+        linestyle='--',
+    )
+
+    @nb.precomp(sustf=precomp_suspension)
+    def calc(self, sustf):
+        species = self.ifo.Infrastructure.ResidualGas.Ar
+        return _calc_residual_gas_damping(self.freq, self.ifo, species, sustf)
+
 def _calc_residual_gas_damping(f, ifo, species, sustf):
     """Redundant calculation of res gas damping for each species