Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found

Target

Select target project
  • wenxuan.jia/pygwinc
  • sean-leavey/pygwinc
  • sebastian.steinlechner/pygwinc
  • nicholas.demos/pygwinc
  • chris.whittle/pygwinc
  • raymond.robie/pygwinc
  • mateusz.bawaj/pygwinc
  • anchal.gupta/pygwinc
  • 40m/pygwinc
  • evan.hall/pygwinc
  • kevin.kuns/pygwinc
  • geoffrey-lovelace/pygwinc
  • brittany.kamai/pygwinc
  • daniel-brown/pygwinc
  • lee-mcculler/pygwinc
  • jameson.rollins/pygwinc
  • gwinc/pygwinc
17 results
Show changes
Commits on Source (3)
  • Anchal Gupta's avatar
    Modified getCoatBrownian() · 15b14b43
    Anchal Gupta authored and Jameson Rollins's avatar Jameson Rollins committed
    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.
    15b14b43
  • Jameson Rollins's avatar
    Merge branch 'master' into 'master' · 8d4dbe85
    Jameson Rollins authored
    Adding Hong et al. calculations for Coating Brownian Noise
    
    See merge request gwinc/pygwinc!52
    8d4dbe85
  • Jameson Rollins's avatar
    remove comments about bulk/shear loss parameters in aLIGO · 6348af30
    Jameson Rollins authored
    We don't want to advertise these parameters yet, since they will
    likely be changed to support different loss functions.
    6348af30
......@@ -188,7 +188,7 @@ Suspension:
## Optic Material -------------------------------------------------------
Materials:
MassRadius: 0.17 # m;
MassRadius: 0.17 # m;
MassThickness: 0.200 # m; Peter F 8/11/2005
## Dielectric coating material parameters----------------------------------
......
......@@ -296,10 +296,21 @@ class CoatingBrownian(nb.Noise):
w0, wBeam_ITM, wBeam_ETM = arm_cavity(self.ifo)
dOpt_ITM = coating_thickness(self.ifo, 'ITM')
dOpt_ETM = coating_thickness(self.ifo, 'ETM')
mTi_ITM = None
mTi_ETM = None
Ic = None
if 'IncCoatBrAmpNoise' in materials.Coating:
if materials.Coating.IncCoatBrAmpNoise.lower() == 'yes':
mTi_ITM = self.ifo.Optics.ITM.Transmittance
mTi_ETM = self.ifo.Optics.ETM.Transmittance
Ic = ifo_power(self.ifo)[1] # Circulating Arm Power
precomp_mirror(self.freq, self.ifo)
nITM = noise.coatingthermal.coating_brownian(
self.freq, materials, wavelength, wBeam_ITM, dOpt_ITM)
self.freq, materials, wavelength, wBeam_ITM,
dOpt_ITM, Ic, mTi_ITM)
nETM = noise.coatingthermal.coating_brownian(
self.freq, materials, wavelength, wBeam_ETM, dOpt_ETM)
self.freq, materials, wavelength, wBeam_ETM,
dOpt_ETM, Ic, mTi_ETM)
return (nITM + nETM) * 2
......
This diff is collapsed.