A buffer gas beam source for short, intense and slow molecular pulses

Abstract

Experiments with cold molecules usually begin with a molecular source. We describe the construction and characteristics of a cryogenic buffer gas source of CaF molecules. The source emits pulses with a typical duration of 240 $\mathrm{\mu }$s, a mean speed of about 150 m/s, and a flux of $5\times {10}^{10}$ molecules per steradian per pulse in a single rotational state.

Keywords:

• Buffer gas
• molecular beam
• radicals
• laser cooling molecules

View correction statement:

Publisher’s Note

1. Introduction

A variety of applications in physics and chemistry call for molecules cooled to low temperatures [1–6]. The starting point for many experiments with cold molecules is the pulsed molecular beam. For example, beams of cold molecules are used to test fundamental physics [7–14] and to study collisions and reactions [16–18].

There are two main methods for producing pulsed beams of cold molecules, supersonic expansion and buffer gas cooling. In the first, the molecules of interest are entrained in a carrier gas which cools to low temperature via an adiabatic expansion from high to low pressure. The initial thermal energy of the carrier gas is converted to forward kinetic energy, resulting in fast-moving molecular beams. These beams have been used for a vast array of spectroscopic and collisional studies, and have been slowed using time-varying electric [19] or magnetic fields [20] and subsequently trapped [21]. In the second method, the molecules are produced inside a cryogenically-cooled cell where they thermalize with a cold buffer gas of helium or neon and then exit through a hole in the cell to make a beam [22–25]. This method can produce slower beams that contain more molecules, though the number density may be lower because the pulse duration is typically much longer. Molecules in the low velocity tail of such a buffer gas beam can be selected [26,27], trapped [28,29] and cooled [30]. Recently, buffer-gas-cooled beams of a few molecular species have been slowed to low velocity using radiation pressure [31–35], and these slow molecules have been captured and cooled to low temperatures in magneto-optical traps [36–38]. The high intensity and low velocity of the beams produced by the buffer gas method was crucial to the success of these experiments. Electric and magnetic deceleration techniques would benefit from these advantages too. However, these techniques usually accept only a small slice of the molecular pulse, preserving its density, and so the long pulses produced by a buffer gas source are unwelcome. In this paper, we describe a buffer gas source of CaF molecules that produces pulses with short duration, high intensity and low velocity. Our aim is to provide information that is useful to others who wish to build a similar source. We present the design of the source, give practical details about how to build it, and then characterize its performance.

2. Design

The working principles of buffer-gas-cooled beams have been reviewed extensively before [25]. We aimed to design a buffer-gas cell that would ensure efficient and rapid extraction of molecules produced inside by laser ablation. To achieve that, an objective was to minimize regions of low helium flow, since molecules entering such regions are likely to be lost, and to prevent vortex formation where molecules can be trapped [39]. Figure 1 illustrates the basic design of our source, and shows the flow of helium through the cell predicted using finite-element modelling software¹ [40]. The source has some featuresin common with those developed by Smalley and co-workers to produce supersonic metal cluster beams [41], and others in common with buffer gas sources, so can be thought of as a hybrid of these two source types. We produce CaF inside the cell by laser ablation of a calcium rod in the presence of SF${}_6$. Cold helium enters the cell through an angled tube, flows along the angled base of the cell and is directed towards the ablation target. This creates a region of high flow and high density near the target. The high density confines the ablation products, preventing them from reaching the cell walls, and also ensures rapid thermalization with the cold helium. The high flow results in rapid extraction of molecules from the cell, yielding a short pulse and leaving little time for diffusion to the walls. The simulation predicts that molecules will exit the cell with a mean velocity of 120 m/s in a pulse that has a full width at half maximum of 340 $\mathrm{\mu }$s. We will see in section 4 that the measured parameters are not too different from these.

Figure 1. Design of buffer gas source showing the simulated flow of helium gas through the cell [39].

3. Construction details

After experimenting with a few different versions of this basic design, we settled on the one shown in Figure 2. This is the design we used for our recent experiments on laser slowing and magneto-optical trapping of CaF molecules [35,37]. The length of the cell is extended by 10 mm, and we have added a conical aperture with an inner diameter of 3.5 mm at the exit. The cell is discharge machined from an oxygen-free-copper block. It has a 10 mm diameter bore and a base angled at ${22}^{\circ }$ to the horizontal where the cold helium gas enters. Warm SF${}_6$ gas enters from a capillary attached to one edge of the cell. The ablation target is a Ca rod with a diameter of 5.9 mm and a length of 30 mm. We roughly cut the target from a Ca ingot and then turn it into a cylinder on a lathe. Once formed, the rod is cleaned with solvent and then stored under vacuum or mineral oil until needed. Just before the main vacuum chamber (see below) is closed, the target is polished with fine grit sand paper, attached to a finely threaded aluminium screw, and inserted into the cell through a 6 mm diameter bore. A thin layer of Apiezon N vacuum grease is applied to the thread to prevent seizing. Rotating the screw rotates and translates the target, increasing the surface area that can be ablated and thus the lifetime of the target. For ablation, we use a Nd:YAG laser producing pulses of 5 ns duration and 5 mJ energy at a wavelength of 1064 nm. We focus the output of the laser onto the target using a 75 cm focal length lens positioned about 55 cm from the target. The beam enters through a window offset from the cell by a narrow tube, which we call the snorkel. This helps to prevent ablation products coating the window.

Figure 2. Design details of the buffer gas cell.

Figure 3 shows a cut-away section through the centre of the source chamber. The cell is attached to a copper plate which is then mounted onto the cold head of a two stage Gifford–McMahon cryocooler.² The first and second stages of this cryocooler have cooling powers of 40 W at 50 K and 1.5 W at 4.2 K, respectively. Under our normal operating conditions, the second stage has a temperature of about 4 K. A vacuum chamber houses the cryocooler and provides ports for optical access and for the required gas, rotary and electrical feedthroughs. A turbomolecular pump with a pumping speed of 600 l/s for helium evacuates the chamber, reaching ${10}^{-7}$ mbar when the cryocooler is off.

An aluminium cylinder attached to the first stage of the cryocooler is cooled to about 30 K and serves as a radiation shield. This reduces the radiative heat load onto the second stage from approximately 50 mW/cm${}^2$ to below 0.05 mW/cm${}^2$. Aluminized mylar applied to the radiation shield can reduce the heat load further. Holes in the radiation shield let the ablation laser beam and gas pipes enter, and the molecular beam exit. Inside the radiation shield, surrounding the cell and attached to the 4 K cold head, is a copper cylinder. Coconut charcoal is glued to the inside of this cylinder using thermally conductive epoxy.³ When cooled below 8 K, the charcoal acts as an efficient cryopump for helium gas [42,43]. The preparation and condition of the charcoal is important. We degrease the copper and paint on a thin layer of epoxy with a brush. The layer should be thin enough that most of the charcoal surface is exposed, but thick enough to provide good thermal conduction. We sprinkle on the charcoal, leave it to dry for a few hours, and finally bake the cylinder in air at 400 K for 24 h. Once installed and under vacuum, we heat the cryopump to $\approx 330$ K for 48 h. We have observed an improved performance of the pump after a few cycles of cooling down to 4 K and warming back up to 330 K.

Room temperature helium gas, with 99.999 % purity, is fed into the vacuum chamber through 1/4-inch stainless-steel tubing. Inside the vacuum chamber, the stainless-steel tube is connected to a 1/4-inch copper tube which winds around a bobbin that is thermally anchored to the first stage of the cryocooler, cooling the gas to about 35 K. The tubing is soft-soldered to the bobbin along the length of the winding for good thermal contact. A second copper bobbin is attached to the second stage of the cryocooler, and copper tubing is wound and soft-soldered around this bobbin to cool the gas to about 4 K before it enters the cell. The end of this copper tube is soldered into the cell. The two bobbins are separated by a 20 cm long section of thin-walled 1/4-inch stainless steel tubing to minimize the thermal conductivity between them. To make each connection between copper and stainless-steel, we braze the copper tube into a stainless-steel Swagelok union and connect the stainless-steel tube to the union using a standard Swagelok connector. The brazed joints and Swagelok connectors are helium leak checked before closing the vacuum chamber. The helium flow is regulated with a flow controller.⁴ When the flow rate is 1 sccm, the mean helium density in the cell is about ${10}^{15}$ cm${}^{-3}$.

A novel feature of our cryogenic source is the formation of molecules by ablation of a pure metal in the presence of a reactant gas, SF${}_6$ in this case [40,44]. We have found this to be better than using an ablation target formed from a mixture of powders. Lower energy ablation pulses can be used, the signal is more stable, and we do not observe much powdered residue building up in the cell, as reported by others [45,46], even after long periods of operation. It has been suggested that such particulates, or ‘dust’, might be responsible for reducing the lifetime of molecules in buffer gas cells [47]. The room temperature SF${}_6$ gas is fed into the vacuum chamber through a 1/4-inch stainless steel tube which connects to a copper capillary inside the vacuum chamber. This capillary is connected to the cell at the position shown in Figure 2. At atmospheric pressure, the melting point of SF${}_6$ is 222 K. To prevent the gas from freezing, the capillary is thermally insulated from the cell by a polyimide (Vespel) spacer. The end of the capillary should be flush with the internal edge of the cell wall. If it is not inserted far enough SF${}_6$ can freeze on the walls of the insulator or cell and clog the capillary. If it is inserted too far into the cell the hot tube heats the helium resulting in a hotter and faster molecular beam. A thermistor measures the temperature of the capillary close to the point where it enters the cell, and a heater wrapped around the capillary can be used to warm it up if the SF${}_6$ freezes. In practice, the thermal isolation is good enough that the heater is never used.

At a distance of 30 mm from the exit aperture of the cell, the molecular beam passes through a 6 mm aperture in a charcoal-covered copper plate attached to the 4 K cold head. This reduces the helium gas load into the rest of the system, and is similar to the design in reference [46]. Installing it reduced the pressure measured outside the radiation shields by a factor of 2.

We operate the source at a repetition rate of 2 Hz. The molecule flux gradually decays, so we rotate the target to expose fresh calcium after about ${10}^4$ shots. The variation in molecule number for various ablation spots is up to 50%, but a spot with similar ablation yield can be found relatively quickly. The target is turned by the aluminium screw (see Figure 2) which is attached to a rotary vacuum feedthrough by a universal joint. The universal joint is thermally connected to the radiation shield via a bronze brush bushing and is thermally disconnected from the room temperature feedthrough and the cold aluminium screw using polyether ether ketone (PEEK) spacers.

After about 24 h of continuous operation, the cryopumps start to saturate. We regenerate them by warming them to a temperature exceeding 20 K. The helium then desorbs and is pumped away by the mechanical pumps. We also observe a significant reduction in the flux of molecules after about 20 h of operation which cannot be recovered either by turning the target or regenerating the charcoal pumps. It may be that in this time enough SF${}_6$ freezes to the cell walls to obstruct either the SF${}_6$ input or the exit aperture of the cell. We find that we can recover the original flux by going through a heat up and cool down cycle overnight. We heat the source to at least 240 K in 6-7 h using 50 W power resistors attached to the first and second stages of the cryocooler. We then immediately start cooling the source back down to 4 K. This takes approximately 6 h so that the source is ready to use again the next morning.

Early in the development of the source we noticed a gradual increase of the velocity of the molecular beam over weeks of operation. The thermal cycling described above did not fix this. Only by opening the source and cleaning the cell could we recover the original low beam velocity. We later discovered a small leak in one of the gas lines which was allowing air to leak into the cell. We hypothesize that this was oxidizing the ablation products deposited on the cell walls. These oxidized products then form a thermally insulating layer which prevents the helium from fully thermalizing (or re-thermalizing [48]) with the cell walls, leading to a hotter and faster molecular beam. Since fixing this small leak in the gas line we have been able to operate the source for many months without it requiring any maintenance.

Figure 3. Cut-away section through the centre of the source. The cryocooler head (blue) sits inside a vacuum chamber (grey). An aluminium radiation shield (pink) is attached to the first stage of the cryocooler. The copper inner cylinder (orange), attached to the second stage of the cryocooler, is coated with coconut charcoal and acts a sorb for helium. The buffer gas cell is attached to the second stage.

4. Characteristics

The molecules are detected by cw laser-induced fluorescence detection on the $A^2{\mathrm{\Pi }}_{1/2}(v=0,J=1/2)\leftarrow X^2{\mathrm{\Sigma }}^{+}(v=0,N=1)$ transition near 494431 GHz. This transition has four resolved components due to hyperfine and spin-rotation interactions. We apply radio-frequency sidebands to the probe laser in order to address all four components simultaneously [32,35]. The probe laser beam is circular, has a total power of 0.5 mW and an intensity distribution with an rms radius of 1.5 mm. This probe beam intersects the molecular beam at ${90}^{\circ }$.

We first searched for a signal within the cryogenic region of the source, between 1.5 and 4 cm from the cell aperture. Laser-induced fluorescence from this region was imaged onto a photomultiplier tube (PMT) through a 1 inch window in the lid of the vacuum chamber and corresponding holes in the radiation shield and cryo-sorb (see Figure 3). After optimizing the signal at this position we probed the molecular beam further downstream, first at 53 cm and then at 130 cm. By comparing the signal inside and outside the cryogenic region, we can determine whether the molecular beam is attenuated by collisions with the helium buffer gas. This can be done either by calibrating the two detectors and then comparing the flux at the two positions, or by recording the flux as a function of the He flow rate inside and outside the cryogenic region. If the maximum molecular flux occurs at a lower flow rate at the downstream position than the upstream position, it indicates that the pumping speed for helium is insufficient. The source then has to be operated at a flow rate which does not extract the most molecules from the cell. In that case, more charcoal, and smaller apertures in the radiation shield and charcoal-covered plate are likely to improve the signal and stability of the source.

Figure 4. Typical time of flight profiles recorded at (a) 2.5 cm and (b) 53 cm from the cell aperture.

Figure 4(a) shows a typical time of flight profile recorded 2.5 cm from the cell aperture. The molecular flux peaks 300 $\mathrm{\mu }$s after the ablation laser fires. Accounting for the time taken to travel the 2.5 cm from the exit, we conclude that the molecules exit the cell about 100 $\mathrm{\mu }$s after they are produced, showing that they are extracted rapidly. Modelling the interior volume of our cell as a cube of side length 1 cm, and using a typical diffusion cross section of ${10}^{-14}$ cm${}^2$ and a typical helium density of ${10}^{15}$ cm${}^{-3}$, we find that the diffusion time to the walls is about 100 $\mathrm{\mu }$s. So we see that the extraction time and the diffusion time are similar. The pulse of molecules measured 2.5 cm from the aperture has a duration (full width at half maximum) of 275 $\mathrm{\mu }$s. This is far shorter than other buffer-gas sources which typically emit pulses with durations from a few ms up to $\sim 100$ ms [46]. The short duration is a consequence of the small cell volume and rapid extraction of the molecules. Figure 4(b) shows a typical time of flight profile recorded 53 cm from the cell aperture. Here, the detection efficiency is much higher, which is why the signal rate is similar to that measured at the upstream detector. The signal peaks at 3.14 ms, indicating that the mean speed is about 170 m/s. At 4 K, an effusive beam of CaF would have a mean speed of 45 m/s, while a supersonic beam of helium has a speed of 204 m/s. In common with other buffer gas sources, the mean speeds we measure lie between these two extremes. It is possible to reach lower speeds by operating in a more effusive regime [22] or by using a two-stage buffer gas source [49], but the flux obtained is considerably lower.

Figure 5 shows velocity distributions for four different helium flows. These distributions are determined from the time-of-flight profiles measured a distance $L=130$ cm from the source, assuming the relation $v=L/t$ between the velocity v and arrival time t. This relation is accurate when the arrival time is far greater than the duration of the pulse exiting the source, which is the case for these measurements. The method has been verified previously [35]. We see from Figure 5 that the signal increases with helium flow rate in the range 0.1-0.8 sccm. In this range, the number of molecules scales approximately linearly with the flow rate. At flows exceeding 1 sccm this number saturates. We expect this saturation to occur once the extraction time from the cell is shorter than the diffusion time to the walls, because then most of the molecules produced in the cell are extracted into the beam [25]. This occurs once the flow is greater than $\approx 1$ sccm. The dashed lines in Figure 5 are fits to the model $f(v)=A{v}^{2}exp[-M{(v-v_0)}^2/(2k_{\mathrm{B}}T)]$, where M is the mass of a CaF molecule, A is an amplitude, T is the temperature and $v_0$ is a velocity boost due to entrainment with the helium. The measured distributions tend to have more molecules at high velocity than this model, but otherwise the model describes the data well. The best-fit value of $v_0$ increases with flow, from 114 m/s at 0.1 sccm to 153 m/s at 0.8 sccm. The best-fit value of T decreases from 5.8 K at 0.1 sccm to 4.8 K at 0.8 sccm. We see that higher flows produce higher velocities and lower temperatures, as a result of the flow regime shifting towards the supersonic limit and the molecules becoming more fully thermalized. Thus, while higher flows give more molecules in total, lower flows yield more molecules at the lowest velocities.

We note that the translational temperature and the shape of the velocity distribution depend on the ablation energy and alignment of the ablation laser. We sometimes see distributions characterized by $T\approx 10$ K, having long tails extending to high velocities, indicating incomplete thermalization. We can understand this with a simple model. Taking a typical mean density in the cell of ${10}^{15}$ cm${}^{-3}$ (for flows around 0.5 sccm) and a typical collision cross section of ${10}^{-14}$ cm${}^2$, the mean free path is about 0.7 mm. Therefore, molecules have roughly 40 collisions before escaping from the cell, perhaps somewhat more since the density is higher in the vicinity of the target. After n hard-sphere collisions with stationary buffer gas atoms of mass m, a molecule of initial energy $E_0$ has an energy $E_n=E_0{(1-f)}^n$ where $f=2Mm/{(m+M)}^2=0.12$. This gives $E_{40}=0.006E_0$. Given that the temperatures involved in laser ablation are often a few thousand kelvin [50], this number of collisions is barely sufficient to cool the molecules to the helium temperature. For other species, especially heavier ones or ones that require higher ablation energies for their production, more collisions may be needed to cool the molecules, requiring a longer cell or a smaller exit aperture. For those constructing a similar source for a different molecule, we recommend experimenting with the dimensions of the cell.

Figure 5. Solid lines: velocity distributions measured 130 cm from the source for four different helium flows. Dashed lines: fits to the model $f(v)=A{v}^{2}exp[-M{(v-v_0)}^2/(2k_{\mathrm{B}}T)]$.

We turn now to the determination of the flux of molecules. Let $N_{\mathrm{mol}}$ be the number of molecules in the ${\text{X}}^2{\mathrm{\Sigma }}^{+}(v=0,N=1)$ state per unit solid angle per pulse, and $N_{\mathrm{ph}}$ be the number of photons detected per pulse at a detector a distance L from the source. They are related by$$\begin{array}{c}\hfill N_{\mathrm{ph}}=\frac{N_{\mathrm{mol}}}{L^2}\int \int \mathit{ϵ}(x,y)p(y)\mathrm{d}x\mathrm{d}y.\end{array}$$

Here, x is the displacement from the centre of the detection system in the direction of the probe beam, y is the displacement along the line to the PMT, $\mathit{ϵ}(x,y)$ is the efficiency of detecting photons emitted at position (x, y), and p(y) is the number of photons emitted per molecule as it travels through the probe beam with displacement y. The fluorescence is collimated by an aspheric lens of diameter 50 mm and focal length 39 mm, and then imaged onto the PMT with a magnification of 2 by a second lens. A spherical mirror retro-reflects light emitted in the opposite direction to the PMT, doubling the collection efficiency. An aperture in the imaging plane ensures that only light from a small rectangular region can be detected. This region is 3 mm wide in the direction of the molecular beam and $w=5$ mm wide along x. The detection efficiency, $\mathit{ϵ}$, is uniform across this small area, and the depth of field is large enough that $\mathit{ϵ}$ is also independent of y for all values of y where photons are emitted. Taking into account the collection efficiency, PMT quantum efficiency and the transmission of lenses and windows, we estimate $\mathit{ϵ}=0.01$. With this simplification, we can write(1) $$\begin{array}{c}\hfill N_{\mathrm{mol}}=\frac{N_{\mathrm{ph}}{L}^2}{\mathit{ϵ}w\int p(y)\mathrm{d}y}.\end{array}$$(1)

At $L=53$ cm we detect $N_{\mathrm{ph}}=1.15\times {10}^5$ photons per pulse. We use a rate equation model [51] to calculate how many photons a molecule emits within the detection area as it travels through the probe beam, averaged over the distribution of velocities. Because the probe beam size is small, the integral rapidly converges as its limits are extended outwards. The calculation gives $\int p(y)\mathrm{d}y=19.6$ mm. To verify this rate model calculation, we measure how $N_{\mathrm{ph}}$ changes with probe laser power, and how it changes between a single frequency probe and a probe with rf sidebands added (see above), then compare these results to the predictions of the model. They are in good qualitative agreement, but the comparisons indicate that, for the parameters used in our estimate of $N_{\mathrm{mol}}$, the model overestimates $\int p(y)\mathrm{d}y$ by about 50%. We think the scattering rate may have been limited by optical pumping into dark states which were not effectively de-stabilized by the small magnetic field in the probe region, an effect not captured by the rate model. Therefore, we take $\int p(y)\mathrm{d}y=13$ mm. Using these values in equation (1(1) $$\begin{array}{c}\hfill N_{\mathrm{mol}}=\frac{N_{\mathrm{ph}}{L}^2}{\mathit{ϵ}w\int p(y)\mathrm{d}y}.\end{array}$$(1) ) we obtain $N_{\mathrm{mol}}=5\times {10}^{10}$ per steradian per shot, accurate to within a factor of 2. This flux is similar to that obtained from other buffer gas sources where molecules are produced by laser ablation [46,52].

The source is very stable when operated at a repetition rate of 2 Hz. The standard deviation of the flux measured over ${10}^4$ shots is 12% of the mean flux. This reduces to 6 % when the shots are phase-locked to the compression cycle of the cryocooler. Increasing the repetition rate to 5 Hz reduces the mean number of molecules per shot by a factor of 2 and results in a very erratic flux with far higher shot-to-shot fluctuations. At 10 Hz the signal decays very rapidly, reaching about 2% of its initial value after 500 shots. It is possible that, at higher repetition rates, local heating of the target produces a molten surface that is not amenable to ablation. We are currently investigating the cause with a view to running sources of this kind at higher repetition rates.

5. Conclusions

We have described a cryogenic buffer gas source that produces pulses of CaF molecules with a typical duration of 240 $\mathrm{\mu }$s, a mean speed of about 150 m/s, and a flux of $5\times {10}^{10}$ molecules per steradian per pulse in a single rotational state. The flux diminishes after running the source continuously for about 24 h, but is easily recovered by going through a warm-up and cool-down cycle that takes about 12 h. The source has run for many months without requiring maintenance. It has been an invaluable starting point for making cold, slow, velocity-controlled beams by frequency-chirped laser slowing [35], and for loading a magneto-optical trap of CaF [37]. As well as these applications, the source may be useful for spectroscopic studies and for loading Stark and Zeeman decelerators. The low velocity and high density of the pulses produced will help to extend the applicability of these deceleration devices to heavier molecules [53,54]. By providing slower, more intense beams, the source could be used to improve experiments that test fundamental physics, such as the measurement of the electric dipole moment of the electron or proton [7,55], or studies of nuclear parity violation using molecules [13].

6. Dedication

We dedicate this paper to Danny Segal who was a good friend and brilliant colleague to many of us. Danny helped us to get started on buffer gas cooling and worked with us on early versions of our sources [48,56]. As often happens with new projects we had plenty of enthusiasm but no equipment, no funding, and little expertise. To get going, we needed a helium cryostat. Danny went for a walk around the department and returned an hour or two later with two perfectly functional cryostats. He was endlessly generous, and he brought out the generosity in others. Unlike the rest of us, Danny knew a bit about cryogenics, so he helped us modify one of the cryostats and taught us how to use it. We used it for all our early work on buffer gas cooling. Scientific progress was a source of joy for Danny. He would have been delighted to see how those first endeavours cobbled together from borrowed parts, requiring constant maintenance, tinkering and botching, led to the trouble-free source described in this paper and to the new science with ultracold molecules that has now become feasible.

Acknowledgements

We are grateful to John Dyne, Giovanni Marinaro and Valerijus Gerulis for technical assistance.

Additional information

Funding

This work was supported by the EPSRC [grant number EP/I012044], [grant number EP/M027716] and [grant number EP/P01058X/1], and by the European Research Council under the European Unions Seventh Framework Programme (FP7/2007-2013)/ ERC [grant agreement 320789].

Notes

No potential conflict of interest was reported by the authors.

The underlying research materials for this article can be accessed at Zenodo (https://doi.org/10.5281/zenodo.995663) and may be used under the Creative Commons CCZero license.

1 Star-CCM+.

2 Sumitomo RDK-415D.

3 Stycast 2850 FT with 23 LV catalyst.

4 Alicat Scientific.