Attachment 'paramC1DRMI.m'
Download 1 %---------------- p a r a m C 1 DRMI . m ---------------------
2 % Struct of optical parameters in Caltech 40m interferomter
3 %
4 % created by Kiwamu Izumi (July 20th 2011)
5 %
6 %--------------------------------------------------------
7 %
8 %[Description]
9 % This function defines all the optical parameters such as
10 % * cavity length
11 % * modulation frequencies and their depth
12 % * Mirror properties (reflectiviries, RoC, losses and etc.)
13 % * Microscopic offsets in the positions of the mirrors
14 % * Mechanical responses of the suspensions
15 % * Laser wavelength
16 % * Demodulation phases
17 %
18 % example usage :
19 % par = paramC1;
20 % opt = optC1(par);
21 % where 'opt' is an optic model instance which is created based on
22 % this parameter structure function.
23 %
24 %--------------------------------------------------------
25 %
26 % [Notes]
27 % In the current setting PR3, SR2 and SR3 are omitted for simplicity.
28 % However PR2 is included as a high reflective beam splitter so that
29 % the POP2 (light coming from BS to PRM) signal can be obtained.
30 %
31 % - remainging refinements
32 % 1. implementation of PR3, SR2 and SR3
33 % 2. Mismatch in the input beams, PRCL and SRCL
34 %
35 % Modified by KK 25 Jan 2012
36 %
37 %
38
39
40 function par = paramC1(par)
41
42 % basic constants
43 lambda = 1064e-9; % laser wavelength
44 c = 299792458; % speed of light, of course !
45
46 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
47 % Detector Geometry (distances in meters)
48 % Lengths (designed) see the 40m wiki page for details
49 lPRC = 6.75380; % PRCL: lPRC = lPR + (lIX + lIY) / 2
50 lSRC = 5.39915; % SRCL: lSRC = lSR + (lIX + lIY) / 2
51 lasy = 0.0342; % Schnupp Asy: lasy = lIX - lIY
52 lmean = 4.0; % (lIX + lIY) / 2
53
54 % Mirror curvatures (all dimensions in meters)
55 Ri = 1e9; % input mirrors (IX and IY)
56 Re = 57.37; % end mirrors (EX and EY)
57 Rpr = 115.5; % power recycling mirror
58 Rsr = 142.0; % signal recycling mirror
59
60
61 % Put together all the length parameters into the 'par' variable
62 par.Length.IX = lmean + lasy / 2; % distance [m] from BS to IX
63 par.Length.IY = lmean - lasy / 2; % distance [m] from BS to IY
64 par.Length.EX = 37.7974; % length [m] of the X arm
65 par.Length.EY = 37.7974; % length [m] of the Y arm
66 par.Length.PR = lPRC - lmean; % distance from PR to BS
67 par.Length.SR = lSRC - lmean; % distance from SR to BS
68 par.Length.PR_PR2 = 1.0; % distance from PR to PR2
69 par.Length.PR2_BS = par.Length.PR - par.Length.PR_PR2; % distance from PR2 to BS
70
71 % Put together all the Radius of Curvature [1/m]
72 par.IX.ROC = 1 / Ri;
73 par.IY.ROC = 1 / Ri;
74 par.EX.ROC = 1 / Re;
75 par.EY.ROC = 1 / Re;
76 par.BS.ROC = 0;
77 par.PR.ROC = 1 / Rpr;
78 par.SR.ROC = 1 / Rsr;
79 par.PR2.ROC = 0; % 40m doesn't use curved mirrors for PRC folding
80
81 % Microscopic length offsets
82 dETM = 0; % DARM offset, for DC readout - leave this as zero
83 par.IX.pos = 0;
84 par.IY.pos = 0;
85 par.EX.pos = 0; % Set DARMoffset in your own scripts, not here.
86 par.EY.pos = 0;
87 par.BS.pos = 0;
88 par.PR.pos = 0;
89 par.SR.pos = lambda/4; % pos = lambda/4 for signal recycling. pos = 0 for broadband signal extraction
90 par.PR2.pos = 0;
91
92
93 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
94 % Mirror Parameters
95
96 % HR Transmissivities
97 %par.IX.T = 0.014; % T = 1.4% for ITMX
98 %par.IY.T = 0.014; % T = 1.4% for ITMY
99 par.IX.T = 0; % T = 1.4% for ITMX
100 par.IY.T = 0; % T = 1.4% for ITMY
101
102
103 par.BS.T = 0.5; % T = 50% for BS
104 par.EX.T = 1; % T = 1 for DRMI (T = 15ppm)
105 par.EY.T = 1; % T = 1 for DRMI (T = 15ppm)
106
107 par.PR.T = 0.0575; % T = 5.75% for PRM
108 %par.PR.T = 1; % T = 5.75% for PRM
109 par.SR.T = 0.10; % T = 10% for SRM
110 %par.SR.T = 1; % T = 10% for SRM
111 par.PR2.T = 25e-6; % 25 ppm (assumption)
112
113 % Power reflectivity on AR Surfaces
114 par.IX.Rar = 500e-6; % designed value is 500 ppm
115 par.IY.Rar = 500e-6; % designed value is 500 ppm
116 par.EX.Rar = 200e-6; % designed value is less than 300 ppm
117 par.EY.Rar = 200e-6; % designed value is less than 300 ppm
118 par.BS.Rar = 0; % designed value is less than 600 ppm
119 par.PR.Rar = 0; % designed value is less than 300 ppm
120 par.SR.Rar = 0; % designed value is less than 300 ppm
121 par.PR2.Rar = 0;
122
123 % HR Losses (50 ppm is assumed)
124 par.IX.L = 50e-6;
125 par.IY.L = 50e-6;
126 par.EX.L = 50e-6;
127 par.EY.L = 50e-6;
128 par.BS.L = 50e-6;
129 par.PR.L = 50e-6;
130 par.SR.L = 50e-6;
131 par.PR2.L = 50e-6;
132
133
134 %par.IX.L = 0;
135 %par.IY.L = 0;
136 %par.EX.L = 0;
137 %par.EY.L = 0;
138 %par.BS.L = 0;
139 %par.PR.L = 0;
140 %par.SR.L = 0;
141 %par.PR2.L = 0;
142
143 % mechanical parameters
144 par.w = 2 * pi * 1.0; % resonance frequency of the mirror (rad/s)
145 par.mass = 0.25; % mass of the mirror (kg)
146 par.w_pit = 2 * pi * 0.6; % pitch mode resonance frequency
147
148
149 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
150 % Input Beam Parameters
151 par.Pin = 1; % input power (W)
152 f1 = 11065399; % first modulation frequency
153 f2 = 5 * f1; % second modulation frequency
154
155 Nmod1 = 2; % first modulation order
156 Nmod2 = 2; % second modulation order
157
158 % construct modulation vectors
159 n1 = (-Nmod1:Nmod1)';
160 n2 = (-Nmod2:Nmod2)';
161 vMod1 = n1*f1;
162 vMod2 = n2*f2;
163
164 % make sidebands of sidebands
165 for i=1:length(n1) % run through all the f1 components
166 for j=1:length(n2) % run through all the f2 components
167 index = j + (i-1)*length(n1); % index for new arrays
168 vFrf(index) = vMod1(i)+vMod2(j); % frequency of sidebands and carrier
169 end
170 end
171 vFrf = sortrows(vFrf');
172
173 %vFrf = sortrows(f1*(-18:18)');
174
175 % input amplitude is only carrier and zero ampliutde in sidebands.
176 % because two RF modulators are placed after the laser source
177 nCarrier = find(vFrf == 0, 1);
178 vArf = zeros(size(vFrf));
179 vArf(nCarrier) = sqrt(par.Pin);
180
181 par.Laser.vFrf = vFrf;
182 par.Laser.vArf = vArf;
183 par.Laser.Power = par.Pin;
184 par.Laser.Wavelength = lambda;
185
186 par.Mod.f1 = f1;
187 par.Mod.f2 = f2;
188 par.Mod.g1 = 0.1; % first modulation depth (radians)
189 par.Mod.g2 = 0.1; % second modulation depth (radians)
190
191 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
192 % Adjustment of demodulation phase
193 % Demodulation Phases -- tuned with newSensMat.m
194 % All the units are in 'degree'
195
196 par.phi.phREFL1 = -176.9563; % f1 : adjusted for CARM, I-phase
197 par.phi.phREFL2 = -86; % f2 : adjusted for CARM, I-phase
198 par.phi.phREFL31 = 9.1369+1.1006; % 3*f1: adjusted for PRCL, I-phase
199 par.phi.phREFL32 = -138.5568; % 3*f2: adjusted for SRCL, I-phase
200
201 par.phi.phAS1 = -87.6498 - 90.; % f1 : adjusted for DARM, Q-pjase
202 par.phi.phAS2 = -165.9168 - 90.-2; % f2 : adjusted for DARM, Q-phase
203 par.phi.phAS31 = 2.3937 - 90.; % 3f1: adjusted for MICH, Q-phase
204 par.phi.phAS32 = -15.9515 - 90;
205
206 par.phi.phPOP1 = -9.9161; % f1 : adjusted for PRCL, I-phase
207 par.phi.phPOP2 = 128.4403; % f2 : adjusted for SRCL, I-phase
208 par.phi.phPOP31 = -10.7091; % 3f1: adjusted for PRCL, I-phase
209 par.phi.phPOP32 = 24.2958; % 3f2: adjusted for SRCL, I-phase
210
211 par.phi.phPOX1 = -86.3637; % f1 : adjusted for PRCL, I-phase
212 par.phi.phPOX2 = 32.3234 - 90; % f2 : adjusted for MICH, Q-phase
213 par.phi.phPOX31 = 120.3357; % 3f1: adjusted for PRCL, I-phase
214 par.phi.phPOX32 = -85.8551 - 90; % 3f2: adjusted for MICH, Q-phase
215
216 par.phi.phPOY1 = 0;
217 par.phi.phPOY2 = 0;
218 par.phi.phPOY31 = 0;
219 par.phi.phPOY32 = 0;
220
221
222 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Attached Files
To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.You are not allowed to attach a file to this page.