A couple of coating stacks with better tolerance (transmission +- 10%). Attachments #1-2 show the spectral reflectivities for ETM/ITM respectively, while Attachments #3-4 show the corner plots. I think the tolerances are inflated by the fact that all the stack indices and thicknesses are varying, while in reality these two effects are degenerate because what matters is the optical thickness. I will try to reflect this in the MCMC code next. Finally, attachments # 5-6 are the hdf5 files with the optimization results.

Sat Amp 

- Changes to sat amp 15.8 k ohm resistors instead of 16k The change has been made on Sat Amp - S1103733 & S1103732 ONLY Channel 4 and 2 have been changed on both boards.

OSEM

- I developed a test bed for our OSEM to measure force 

I will attach images of the setup and some of the results from 3 different OSEMs.

Future Work

- For the current test bed, we are using a clear plastic bin although not ideal it manages to get the job done and works for now there could be a better solution for this,
- Next step for OSEM we want to use 40 m single pendulum to test OSEM and measure the transfer function.

In the following test, a single Sat Amp chassis that holds Sat Amp Board S1106078 and S1106077

Verification of Sat Amp

First, as the test of the LED driver circuits in the chassis, 8 of various color LEDs were inserted to the appropriate output pins of the chassis. This resulted in all the LED lit and the LED mon TP was confirmed to have voltage outputs of 5V. (See my previous ELOG)

OSEM tests

Connected OSEMs to the sat amp to test the OSEM LED/PD pairs. With the first test, the PD out gave us 15V. We wondered if this was just the problem of the OSEM or circuit, or just there are too much light for the transimpedance gain of 121K Ohm.

By blocking the OSEM light by a random heat shrink tube found on the table, we saw the number got reduced. This indicates that the OSEM/Satamp outputs are working and there are just too much light.

We decided to reduce the gain: The transimpedance gain R18 was reduced to 16k, which gave us a voltage range from 5V~7V  with some outlier OSEMS at 1V (See the attached table)

There are 24 total OSEMs:

• one apparently not functional
• 20 fell in the range of 5~7V
• 3 fell in the range of  ~1V

(These numbers given after the change of R18 to 16k Ohm)

Note: We originally aimed for 8~9V. However, from a previous study of OSEM at cryogenic temperature, we learned that there was about an about 30% increase in the response.
Therefore, we decided to leave a sufficient margin from 10V considering this expected increase in the response.

Yesterday, I did two rounds of slowly heating 4 samples to the maximum hot plate temperature. This was to formally test if my success with a single sample earlier in the week was a fluke. Note that the hot plate takes about 10-15 minutes to reach a stable temperature when it is turned up one notch.

First round:
I bonded 4 samples by putting methanol in the gap between the glass slides and letting it dry to create a gap.
Starting at room temperature, I heated the slides on each setting for roughly 15 minutes, then let them cool down naturally over the course of an hour. 3 broke broke at medium heat, and 1 survived the whole process. I belive these broke because the bonds were weaker and I heated them slightly too quickly. In previous tests, I would manually switch the hot plate on and off, but I wanted to see if the hot plate could heat up slow enough on its own.

Second round:
I bonded 4 samples by scrubbing the slides with methanol, using a compressed air duster to blow off the fibers, rubbing them together with the pressure of my fingers, and repeating this whole procedure until they stuck (it took 2-4 repeats).
Starting at room temperature, I heated the slides on each setting for exactly 20 minutes, then let them cool down naturally over the course of an hour. All of them survived to the maximum temperature (the pictures show them at the start and end of the heating, note the temperature). I credit this to the stronger bonding proceedure and the extra 5 minutes for them to adjust to the temperature. I did not turn the hot plate on or off at any point, I just let it heat up at its own rate.

I cannot tell if the bonds are stronger. The size and shape of the Newtons rings did not change.

Sat amp seems to be working just fine. There does seem to be a saturation issue with one of the outputs we may need to change a resistor on the board.

I have found that, after cleaning the glass with methanol (or even sometimes with just a dry lense-cleaning cloth), I can get glass slides to bond by rappidly rubbing them together until something sticks. This was inspired by watching "Wizard of Vaz" perform bonds on YoutTube. While cleaning, I now use enough strength to make the glass squeak, as advised by him.

Upon heating, I encountered the same issue as when I bonded them by putting a liquid (water, methanol, etc.) in the gap, which leads me to now believe that the broken bond is not due to the expansion of a liquid. Further, even at the low temperature of 60°C, placing the liquid-less sample on the hotplate breaks the bond in seconds, which I caught on video. In the attached video*, you can see that, before the heat, the bond is strong enough that I cannot push it appart with my fingers, but after the heat, it slides easily. Note that, outside of taking the video, I always lay the entire slide on the center of the metal so the sample is evenly heated.

*This is my first time attaching a video. If it didn't attach properly, I'll add it on to a later log. I also want to record myself performing the rubbing bonding technique.

we have 23 OSEMS they look all full built and I will try and test them this week and or next week.

Just a general update of what I have been up to deriving Lagrange for double pendulum system and also been looking at code that koji gave to me I've add comment to some of the code also working on my report.

Here are some corner plots to analyze the sensitivity of the designs in the previous elog to a 1% gaussian distributed perturbation using MCMC.

Attachment #1 shows the ETM corner plot 

Attachment #2 shows the ITM corner plot.

I let the indices of both high and low index materials vary, as well as the physical thicknesses and project their covariances to the transmission for PSL and AUX wavelengths.

The result shows that for our designs it is better to undershoot in the optimization stage rather than meet the exact number. Nevertheless, 1% level perturbations in the optical thickness of the stack result in 30% deviations in our target transmission specifications. It would be nice to have a better constraint on how much each parameter is actually varying by, e.g. I don't believe we can't fix the index of refraction to better than 1%., but exactly what its value is I don't know, and what are the layer deposition tolerances? These numbers will make our perturbation analysis more precise.

I've been running the HR coating optimization for mariner TMs. Relative to the specifications found here we now are aiming for

• ITM HR coating of 2000 ppm @ 2050.15 nm, and 1000 ppm @ 1550 nm
• ETM HR coating of 25 ppm @ 2050.15 nm, and 1000 ppm @ 1550 nm.

Both the PSL and AUX cavity finesses range the few couple of thousands, and the goal is not to optimize the coating stack for noise, but more importantly for the transmission values and tolerances. This way we ensure the average finesse and differential finesse requirements are met. Anyways, Attachment #1-2 shows the transmission plots for the optimized coating stacks (so far). Attachments #3-4 show the dielectric stacks. The code still lives in this repository.

I'm on the process of assessing the tolerance of this design stacks against perturbations in the layer thicknesses; to be posted in a follow-up elog.

Note that I am just testing out different techniques, so I have not set up the thermocouples to precisely measure the temperatue.
On Tuesday, I developed a new method of putting water, isopropanol, or methanol on one slide then squishing the other slide on top of it to fill the gap with the afformentioned liquid. The slides are slippery at first, but as they dried, which took about 15 minutes, the bond forms. The bonds were strong enough that I could just barely push the slides appart by applying pressure to the side using my thumbs. I prepared 4 samples this way, 2 with iso and 2 with meth. I took one of each and heated them on Medium for 30 minutes under the brass hunk with the aluminum square on the bottom and copper foil on both sides of the samples. Earlier in the day, I tried heating them without the weight on top, but the heat just broke the bond. I took the remain two and set them aside as controls.
On Thursday, I returned to check the bonds. The heated samples had broken. I intented to check on Wednesday, but I was sick from food poisoning, so I do not know whether the bonds broke immediately after heating or due to sitting for an extra day. For the control samples, one also had a broken bond, but the other had become even stronger.
I noticed that, when the slides are successfully bonded, the shape and appearance of the Newton's rings change, which can be seen in the pictures. I speculate that the circles on the unbroken control are the bonded regions. Ideally, we want to see no Newton's rings.

I have finished all coil driver and sat amp chassis they all seem to be functioning properly.
 

Finished all 3 Coil Drover chassis and power lines still need to install the rear cables will do that after I finish Sat Amp chassis tomorrow. 

Note that the slides have "GLOBE" printed on one side. I always bond the opposite using the opposite side without the text.

On Monday (7/11), I began experimenting with bonding, starting with "air-bonding," which is trying to make dry, gently cleaned slides stick. I achieved my first succesful optical contact with what I call "acidental water-assisted direct bonding" or "water-bonding," where I accidentally clasped two wet slides together while washing my dirty finger prints off them. After the accidental discovery, I repeated it by running water over the slides while there were clasped together and achieved the same result. After a few hours, I attempted partially sliding apart the second water-bonded sample. I could slowly push them apart by pressing my thumbs against the long edge, but it took quite a bit of force. I decided to let 4 samples sit overnight: 1 air-bonded, 1 air-bonded with the brass hunk on top of it, and 2 water-bonded. Neither time nor pressure improved the air-bonded samples as they still slid apart very easily. The first water-bonded sample slid apart easier, but one part remained stubornly attached until I began shaking it violently. The second water-bonded sample was much harder to slide apart than the last time I tested it. With all the force of my fingers, I could barely make it budge.

This was performed last Friday (7/8).

I secured a thermocouple perpendicular against the hotplate and recorded the maximum temperature the hotplate reached at Low, Medium, and High. It took about 5 minutes to reach a stable temperature, where stable means that the temperature stayed within +/- 0.5°C for a minute. The hotplate maintains a certain temperature by turning itself on and off, so the temperature would drop slightly (at most, a few °C) while the hotplate was off. The maximums were:
Low: 51°C
Medium: 185°C
High: 263°C
At the max temperature, I moved the perpendicular thermocouple around to roughly find the variation in tempearture at different locations on the hotplate. Facing the nob, the top right quadrant is about 10-20°C cooler than the other quadrants, which are within 5°C of eachother. Excluding the cooler quandrant, the center and the outer edge are within 5°C of eachother. The temperature increases as one approaches half the radius, with it being about 20-40°C greater than the center and outer edge. The highest recorded temparture was 289°C at half the radius in the bottom left quandrant. This was only meant to be a rough test to see how even the heating is.

All three coil driver boards are complete and have been tested. Modification for all 4 sat amp have been completed. Ideally, I would like to finish all the chassis on Monday I have one just about done. 
 

For our optical contacting, Jennifer and I are starting out with glass (microscope slides), with the setup in the EE shop next to the drill press (photos from Jennifer to follow).

Some interesting links:

• https://www.laserfocusworld.com/optics/article/16546805/optical-fabrication-optical-contacting-grows-more-robust is a write up on contacting, and the link to Dan Shaddock's paper is also useful (need to sign up to get the acutal TSP writeup)
• Thesis on Silicon Bonding (https://escholarship.org/uc/item/5bm8g42k)
• https://youtu.be/qvBoGoh_-AE

Almost done with coil driver boards 

Update of my current work I have finished one coil driver board and started on the last two that I need here is the progress and Ideally, I'll finish by tomorrow. 

I've completed one coil driver board. 
Hopefully next week I can finish the other 2 boards and make the modifications to the sat amp baords. 
 

I've managed to cut and crimp wires for the power board for coil driver. I will begin adding components to the coil driver board.

- Add Components to Coil Driver board 

 - Replace some Sat Amp Componetns

- Still working on moving optical table to CAML

- Unsure if cryochamber has been cleaned and moved

Table moving effort in the OMC lab: See https://nodus.ligo.caltech.edu:8081/OMC_Lab/412

- B246/QIL Skyhook

• Find the base of Skyhook. It should be in the storage room (B246). Stephen contacted Chub for lab access. Done
• Assemble Skyhook with the base and check the stability/safety/capacity/height/etc

- DOPO

• Ask Paco to move the delicate instruments from the table. Done
• Bring Skyhook to DOPO. The chamber seems already vented.
• Find the way to place the cap on the floor safely and cleanly. => Stephen
   
• Open the cap and then remove the crackle interferometer. Wrap it with something and place it somewhere in the room. How? => Stephen
   
• Move the base to a dolly or something. Then put a cap on the base. => It'd be better to ask Caltech Transp for the chamber transportation.
• Do we have to temporarily remove the laser safety curtain?

- OMC Lab

• We probably need to separate the optical table and the base. Ask Caltech Transp to check how the work should be done.
• Do we have to temporarily move anything on the way?
• The table can be rolled out to the corridor and then rolled in to the CAML.

- CAML

• Remove the grey rack and push the desk to the East.
• Place the optical table.
• Place the rack close to the table.

Today we looked at possible locations for where we will be setting up Mariner Suspension and Cryo chamber. The first option was the far left table in the CAML lab but it seems that there is going to be an issue with height clearance, so we have come up with another solution which takes a table from Koji's lab which is 3'x4' ft and moving it into CAML lab in the back right of the lab. To move the table we may need to call facilities to help us because we will most likely need to take the table apart to get it out of the lab. The aisle space in Koji's lab is about 43 inches, but the doorway, which is the tightest space, is 35 inches.

After we have set up the table in CAML we are planning on moving the Chamber in DOPO-lab to CAML. We plan to use skyhook with has a load limit of 500lbs/227kg this should be more than enough to move the chamber. We still need to get the wheeled base for skyhook we are in the works in doing so. 

Also, We want to remove the previous setup from the chamber and leave it at DOPO-lab. Stephen is going to figure out how to keep it clean (sort of). Besides these transportation logistics, I am also working on the electronics as an immediate task and the electrical arrangement in the chamber.

to do list
        - Check the table height
        - Check the chamber height (base/cap)
        - Check how much the chamber cap needs to be lifted (so that we can remove it)
        - Is the weight capacity sufficient?

cool

Does this work? Is this insane?

I've checked the validity of my state space model in a couple of ways so that we have confidence in the results that it gives. I've checked the DC gain of the transfer functions where it is non-zero. I did this by solving the static balance of forces problem in the extended body model by hand to get the DC CoM position as well as the pitch angle of both masses. In the previous ELOG entry I didn't quite do this for all transfer functions so here I completed the check. My values agree with the model's values to within 10% at the worst end and to within 0.1% at the best end. I performed a second check to see if the frequencies occur in the correct places by considering the case of very low coupling between the different resonant modes. It's difficult to check this in the case where the modes are strongly coupled (for example length-pitch is strong or the two pitch modes are close together) but if I sufficiently separate them, I get very good agreement between my analytic approximation and the state space model.

The model can easily be converted from one that gives motion in X and RY into one that gives motion in Y and RX. Running the model for both directions gives the following list of resonances (note pendulum modes in X and Y direction are identical):

Given that I think the model seems to give sensible values, I've pushed the updated model to the GitLab repository. It is now possible to quickly change the parameters of the suspension and very quickly see the corresponding shift in the resonances. To change the parameters, open the plain text file called 'params' and change the values to the new ones. Afterwards, run the file 'ss_extended.py', which will solve the state space model, save the resulting ABCD matrices to a folder and print out the values of the resonances to terminal.

I've been having a look at the transfer functions for the translation and pitch of both masses. I'm attaching the plot of all input-to-output transfer functions of interest so far. Here I've identified the pitch resonances of the two masses (one each) as well as the two pendulum modes. I need to now investigate if they occur in the correct places. I have confirmed the DC response by directly solving the statics problem on paper.

[Paco]

There was an error in the last plot of the previous log. This was correctly pointed out by rana's pointing out that the broadening from air should be independent of the CO2 concentration, so nominally both curves should coincide with each other. Nevertheless, this doesn't affect the earlier conclusions -->

The PMC loss by background, pressure broadened absorption lines at 2049.9 nm by CO2 is < 1 ppm.

The results posted here are reflected in the latest notebook commit here.

[Paco]

I used the HITRAN database to download the set of ro-vibrational absorption lines of CO2 (carbon dioxide) near 2.05 um. The lines are plotted for reference vs wavenumber in inverse cm in Attachment #1.

Then, in Attachment #2, I estimate the broadened spectrum around 2.05 um and compare it against one produced by an online tool using the 2004 HITRAN catalog.

For the broadened spectrum, I assumed 1 atm pressure, 296 K temperature (standard conditions) and a nominal CO2 density of 1.96 kg/m^3 under this conditions. Then, the line profile was Lorentzian with a HWHM width determined by self and air broadening coefficients also from HITRAN. The difference between 2050 nm and 2040 nm absorption is approximately 2 orders of magnitude; so 2040 nm would be better suited to avoid in-air absorption. Nevertheless, the estimate implies an absorption coefficient at 2050 nm of ~ 20 ppm / m, with a nearby absorption line peaking at ~ 100 ppm / m. 

For the PMC, (length = 50 cm), the roundtrip loss contribution by in-air absorption at 2050 nm would amount to ~ 40 ppm. BUT, this is nevery going to happen unless we pump out everything and pump in 1 atm of pure CO2. So ignore this part.

Tue Nov 9 08:23:56 2021 UPDATE

Taking a partial pressure of 0.05 % (~ 500 ppm concentration in air), the broadening and total absorption decrease linearly with respect to the estimate above. Attachment #3 shows the new estimate.

For the PMC, (length = 50 cm), the roundtrip loss contribution by in-air absorption at 2050 nm would amount to ~ 1 ppm.

Some more progress:

- Shaved the height of the top clamp blocks. We can extend the suspension height a bit more, but this has not been done.

- The IM OSEM arrangement was fixed.

- Some EQ stops were implemented. Not complete yet.

All parameters are temporary:

Test mass size: D150mm x L140mm
Intermediate mass size W152.4mm x D152.4mm x H101.6mm
TM Magnets: 70mm from the center

Height from the bottom of the base plate
- Test mass: 5.0" (127mm) ==> 0.5" margin for the thermal insulation etc (for optical height of 5.5")
- Suspension Top: 488.95mm
- Top suspension block bottom: 17.75" (450.85mm)
- Intermediate Mass: 287.0mm (Upper pendulum length 163.85mm / Lower pendulum length 160mm)

OSEMs
- IM OSEMs: Top x2 (V/P)<-This is a mistake (Nov 3 fixed), Face x3 (L/Y/P), Side x 1 (S)
- TM OSEMs: Face x4
- OSEM insertion can be adjusted with 4-40 screws

To Do:
- EQ Stops / Cradle (Nov 3 50% done)
- Space Consideration: Is it too tight?
- Top Clamp: We are supposed to have just two wires (Nov 3 50% done)
- Lower / Middle / Upper Clamps & Consider installation procedure
- Fine alignment adjustment
- Pendulum resonant frequencies & tuning of the parameters
- Utility holes: other sensors / RTDs / Cabling / etc

- Top clamp options: rigid mount vs blade springs
- Top plate utility holes
- IM EQ stops

Discussion with Rana

- Hoe do we decide the clear aperture size for the TM faces?
- OSEM cable stays
- Thread holes for baffles

- Light Machinery can do Si machining
- Thermal conductivity/expansion
- The bottom base should be SUS... maybe others Al except for the clamps

- Suspension eigenmodes separation and temperature dependence

# Deleted the images because they are obsolete.

I've been testing out the extended body lagrangian models and I'm trying to understand the ground motion and force coupling to the test mass displacement. I've compared the two point-mass model to the extended model and, as expected, I get very similar results for the ground coupling. Attachment 1 shows the comparison and asside from more agressive damping of the point-mass model making a small difference at high frequency, the two models look the same. If I look at the force coupling, I get a significantly different result (see attachment 2). I think this makes sense because in the point-mass model I am driving purely horizontal displacement as there is no moment of inertia. However, for the extended body I drive the horizontal position of the centre of mass, which then results in an induced rotation as the change propagates through the dynamics of the system. To obtain a consistent result with the point-mass model, I would need to apply a force through the CoM as well as a counteracting torque to maintain a purely horizontal displacement of the mass. What I am wondering now is, what's the correct/more convenient way to consider the system? Do I want my lagrangian model to (a) couple in pure forces through the CoM and torques around the CoM and then find the correct actuation matrix for driving each degree of freedom in isolation or (b) incorporate the actuation matrix into the lagrangian model so that the inputs to the plant model are a pure drive of the test mass position or tilt?

I reran the cooldown model, setting the emissivity of the inner surface of the inner shield to 0.7 (coating), and the emissivity of the outer surface to 0.03 (polished Al). Previously, the value for both surfaces was set to 0.3 (rough aluminum). 

Attachment 1: TM cooldown, varying area of the inner shield. Now, the marginal improvement in cooldown once the IS area reaches 0.22 m² is negligible. Cooldown time to 123K is ~100 hrs, just over 4 days. I've kept IS area set to 0.22 m² moving forward.

Attachment 2: TM/IS cooldown, considering 2 lengths for the test mass. Choosing l=100m instead of 150mm increases cooldown time from ~100 hrs to ~145 hrs, or 6 days.

I used the same model in [37] to consider how test mass length affects the cooldown. Attachment 1 plots the curves for TM length=100mm and 150mm. The coupling between the test mass and inner shield is proportional to the area of the test mass, and therefore increases with increasing length. Choosing l=100mm (compared to 150mm) thus reduces the radiative cooling of the test mass. The cooldown time to 123K is ~125 hrs or over 5 days for TM length=150mm (unchanged from [37]), but choosing TM length=100m increases this time to ~170 hrs or ~7 days. (Note that these times/curves are derived from choosing an arbitrary inner shield area of 0.22 m², but the relative times should stay roughly consistent with different IS area choices.)

Here are the DAC and residual displacement spectra for different suspension heights ranging from 450 mm to 600 mm. I aimed to get the Q of the lower resonance close to 5 and the DAC output RMS close to 0.5 V but as this was just tweaking values by hand I didn't get to exactly these values so I'm adding the actual values for reference. The parameters are as follows:

Now that I have correct phase and amplitude behaviour for my MIMO state space model of the suspension and the system is being correctly evaluated as stable, I'm uploading the useful plots from my analysis. File names should be fairly self-explanatory. The noise plots are for a total height of 550 mm, or wire lengths of 100 mm per stage. I've also attached a model showing the ground motion for different lengths of the suspension.

I've ironed out the issues with my MATLAB model so that it now shows correct phase behaviour. The problem seems to arise from infinite Q poles where there is an ambiguity in choosing a shift of +/- 180 deg in phase. I've changed my state space model to include finite but very high Q poles to aid with the phase behaviour. The model has been uploaded to the GitLab project under mariner40 -> mariner_sus -> models -> lagrangian.

Following our discussion at the Friday JC meeting, I gathered several resources and made a small simulation to show how frequency combs might be generated on platforms other than microcombs or mode-locked lasers.

Indeed, frequency combs generated directly from a mode-locked laser are expensive as they require ultra-broadband operation (emitting few fs pulses) to allow for f-2f interferometry.

Microcombs are a fancy way of generating combs. They are low-power-consuming, chip-scale, have a high repetition rate, and are highly compatible with Silicon technology. While these are huge advantages for industry, they might be disadvantageous for our purpose. Low-power means that the output comb will be weak (on the order of uW of average power). Microscopic/chip-scale means that they suffer from thermal fluctuations. High rep-rate means we will have to worry about tuning our lasers/comb to get beat notes with frequencies smaller than 1GHz.

Alternatively, and this is what companies like Menlo are selling as full-solution frequency combs, we could use much less fancy mode-locked lasers emitting 50fs - 1ps pulses and broaden their spectrum in a highly nonlinear waveguide, either on a chip or a fiber, either in a cavity or linear topologies. This has all the advantages:

1. High-power (typically 100mW)

2. Low rep-rate (typically 100MHz)

3. Relatively cheap

4. "Narrowband" mode-locked lasers are diverse and can come as a fiber laser which offers high stability.

As a proof of concept, I used this generalized Schrodinger equation solver python package to simulate 1d light propagation in a nonlinear waveguide. I simulated pulses coming out of this "pocket" laser (specs in attachment 1) using 50mW average power out of the available 180mW propagating in a 20cm long piece of this highly nonlinear fiber (specs in attachment 2).

The results are shown in attachments 3-4:

Attachment 3 shows the spectrum of the pulse as a function of propagation distance.

Attachment 4 shows the spectrum and the temporal shape of the pulse at the input and output of the fiber.

It can be seen that the spectrum is octave-spanning and reaches 2um at moderate powers.

One important thing to consider in choosing the parameters of the laser and fiber is the coherence of the generated supercontinuum. According to this paper and others, >100fs pulses and/or too much power (100mW average is roughly the limit for 50fs pulses) result in incoherent spectra which is useless in laser locking or 1f-2f interferometry. These limitations apply only when pumping in the anomalous dispersion regime as traditionally have been done. Pumping in an all-normal (but low) dispersion (like in this fiber) can generate coherent spectra even for 1ps pulses according to this paper and others. So even cheaper lasers can be used. ps pulses will require few meter-long fibers though.

Building on [32], I added a copper cold finger to conductively cool the inner shield, instead of holding the inner shield fixed at 77K. The cold finger draws cooling power from a cyro cooler or "cold bath" held at 60K, for simplicity. I added an outer shield and set its temperature to 100K. The outer shield supplies some radiative heating to the inner shield, but blocks out 295K heating, which is what we want. The expanded diagram can be seen in Attachment 1. 

I wanted to find the optimal choice of inner shield area (A_IS) to maximize the radiative cooling to the test mass. I chose 5 values for A_IS (from A_TM to A_OS) and plotted the test mass cooldown for each in Attachment 2. The radiative coupling between the inner shield and test mass is maximized when the ratio of the areas, A_TM/A_IS, is minimized. Therefore, the larger A_IS, the colder the test mass can be cooled. Even though choosing A_IS close to A_OS increases the coupling between the 2 shields, the resulting heating from the outer shield is negligible compared to the enhancement in cooling.

I chose A_IS = 0.22 m² to model the inner shield and test mass cooldown in Attachment 3. The test mass reaches 123 K at ~ 125 hours, or a little over 5 days. I have pushed the updated script which can be found under mariner40/CryoEngineering/MarinerCooldownEstimation.ipynb.

Agree with this. Quickly running tmm on the same "stacks" gave the Attachment #1-3. (Ignore the vertical axis units... will post corrected plots) and extend the wavelength range to 100 um.

I have to question whether this passes a sanity test. Surely in the case of Stack 2, the 10um thick Ta2O5 will absorb the majority of the incident radiation before it reaches the SiO2 layer beneath. It should at least be similar to just absorption in Ta2O5 with some Fresnel reflection from the AIr-Ta2O5 interface.

For example, at around 18um, K~2, so the amplitude attenuation factor in a 10um thick layer is 160,000x or a gain of 6E-6. So whatever is under the Ta2O5 layer should be irrelevant - there is negligible reflection.

[Paco, Nina, Aidan]

We ran our stack emissivity calculation on different AR stacks to try and make a decision for the TM barrel coatings. This code has yet to be validated by cross checking against tmm as suggested by Chris. The proposed layer structures by Aidan and Nina are:

1. *| Air || SiO2 x 800 nm || Ta2O5 x 5 um || Silicon |*
2. *| Air || Ta2O5 x 10 um || Sio2 x 20 nm || Silicon |*
3. *| Air || SiO2 x 100 nm || TiO2 x 1 um || Silicon |*

Attachments # 1-3 show the emissivity curves for these simple dielectric stacks. Attachment #4 shows the extinction coefficient data used for the three different materials. The next step is to validate these results with tmm, but so far it looks like TiO2 might be a good absorbing film option.

[Paco, Nina, Aidan]

Updated the stack emissivity code to use the Kitamura paper fused silica dispersion which has a prominent 20 um absorption peak which wasn't there before... (data was up to 15 um, and extrapolated smoothly beyond). The updated HR stack emissivities are in Attachments #1 - #2. A weird feature I don't quite understand is the discontinous jump at ~ 59 um ... 

Attachment 1 is a geometric diagram that reflects the current state of the ITM cooldown model, introduced in [30]. The inner shield is assumed to be held at 77K for simplicity, and 2 heat sources are considered: laser heating, and radiative heating from the room-temperature snout opening. The view factor F_ij between the snout opening and test mass (modeled as 2 coaxial parallel discs separated by length L - equation found in Cengel Heat Transfer) is calculated to be 0.022. The parameters used in the model are noted in the figure.

Attachment 2 is a simplified diagram that includes the heating/cooling links to the test mass. At 123K, the radiative cooling power from the inner shield (at 77K) is 161 mW. The radiative heating from the snout opening is 35 mW, and the laser heating (constant) is 101.5 mW. Due to the tiny view factor betwen the snout opening and the test mass, most of the heat emitted by the opening does not get absorbed. 

The magnitudes of heating and cooling power can be seen in Attachment 3. Lastly, Attachment 4 plots the final cooldown curve given this model. 

My next step is to add the outer shield and fix its temperature, and then determine the optimal size/location of the inner shield to maximize cooling of the test mass. This is question was posed by Koji in order to inform inner shield/outer shield geometric specs. Then, I will add a cold finger and cryo cooler (conductive cooling). Diagrams will be updated/posted accordingly.

How about a diagram so that we can understand what this model includes?

*Note: the current modeling script can be found at: CryoEngineering/MarinerCooldownEstimation.ipynb

Nina pointed me to the current mariner cooldown estimation script (path above) and we have since met a few times to discuss upgrades/changes. Nina's hand calculations were mostly consistent with the existing model, so minimal changes were necessary. The material properties and geometric parameters of the TM and snout were updated to the values recently verified by Nina. To summarize, the model considers the following heat sources onto the testmass (P_in):

- laser absorption by ITM bulk (function of incident laser power, PR gain, and bulk absorption)

- laser absorption by ITM HR coating (function of incident laser power and HR coating absorption)

- radiative heating from room-temp tube snout (function of snout radius and length, and TM radius)

The heat transfer out of the testmass (P_out) is simply the sum of the radiative heat emitted by the HR and AR faces and the barrel. Note that the script currently assumes an inner shield T of 77K, and the inner/outer shield geometric parameters need to be obtained/verified.

Nina and Paco have been working towards obtaining tabulated emissivity data as a function of temperature and wavelength. In the meantime, I created the framework to import this tabulated data, use cubic spline interpolation, and return temperature-dependent emissivities. It should be straightforward to incorporate the emissivity data once it is available. Currently, the script uses room-temperature values for the emissivities of various materials. 

Future steps:

- Incorporate tabulated emissivity data

- Verify and update inner/outer shield dimensions

We had a meeting with the code open in ZOOM. Here are some points we discussed:

• The code requires another file ground.m. It is attached here.
   
• The phase of the bode plots were not wrapped. This can be fixed by applying the "PhaseWrapping" options as
  opts=bodeoptions('cstprefs');
opts.PhaseWrapping = 'on';
bode(A,opts)
   
• We evaluated the open-loop transfer function of the system. For this purpose, we added the monitor point ('F') at the actuator and cut the loop there like:
  sys = connect(P, S, W, ADC, Winv, A2, DWinv, Dinv, DAC, DW, D, R, C, {'xg' 'nADC', 'nDAC', 'nd', 'nth'}, 'xt', {'F','VDAC'});
OLTF=getLoopTransfer(sys(1),'F');
figure(2)
clf
bode(OLTF,opts);
   
• We played with the loopgain (Ga2). When Ga2 is a positive number, the loop was stable. We had to shift the low pass cut-off frequency from 10Hz to 12Hz to make the damping of the 2nd peak stable.

[Paco, Nina]

We have been working on an estimate of the wavelength dependent emissivity for the mariner test mass HR coatings. Here is a brief summary.

We first tried extending the thin film optimization code to include extinction coefficient (so using the complex index of refraction rather than the real part only). We used cubic interpolations of the silica and tantala thin film dispersions found here for wavelengths in the 1 to 100 um range. This allowed us to recompute the field amplitude reflectivity and transmissivity over a broader range. Then, we used the imaginary part of the index of refraction and the thin film thicknesses to estimate the absorbed fraction of power from the interface. The power loss for a given layer is exponential in the product of the thickness and the extinction coefficient (see eq 2.6.16 here) . Then, the total absorption is the product of all the individual layer losses times the transmitted field at the interface. This is true when energy conservation distributes power among absorption (=emission), reflection, and transmission:

The resulting emissivity estimate using this reasoning is plotted as an example in Attachment #1 for the ETM design from April. Two things to note from this; (1) the emissivity is vanishignly small around 1419 and 2128 nm, as most of the power is reflected which kind of makes sense, and (2) the emissivity doesn't quite follow the major absorption features in the thin film interpolated data at lower wavelengths (see Attachment #2), which is dominated by Tantala... which is not naively expected?

Maybe not the best proxy for emissivity? Code used to generate this estimates is hosted here.

Here's the DAC voltage spectrum with its associated RMS.

Also, for clarity, this model is for a lossless point-mass double pendulum system with equal masses and equal lengths of 20 cm.